Abstract
Recent studies have identified a role for ALKBH7 in the occurrence and progression of cancer, and this protein is related to cellular immunity and immune cell infiltration. However, the prognostic and immunotherapeutic value of ALKBH7 in different cancers have not been explored. In this study, we observed high ALKBH7 expression in 17 cancers and low expression in 5 cancers compared to paired normal tissues. Although ALKBH7 expression did not correlate relatively significantly with the clinical parameters of age (6/33), sex (3/33) and stage (3/27) in the cancers studied, the results of the survival analysis reflect the pan-cancer prognostic value of ALKBH7. In addition, ALKBH7 expression was significantly correlated with the TMB (7/33), MSI (13/33), mDNAsi (12/33) and mRNAsi (13/33) in human cancers. Moreover, ALKBH7 expression was associated and predominantly negatively correlated with the expression of immune checkpoint (ICP) genes in many cancers. Furthermore, ALKBH7 correlated with infiltrating immune cells and ESTIMATE scores, especially in PAAD, PRAD and THCA. Finally, the ALKBH7 gene coexpression network is involved in the regulation of cellular immune, oxidative, phosphorylation, and metabolic pathways. In conclusion, ALKBH7 may serve as a potential prognostic pan-cancer biomarker and is involved in the immune response. Our pan-cancer analysis provides insight into the role of ALKBH7 in different cancers.
Introduction
The AlkB family consists of Fe (II) and α-ketoglutarate-dependent dioxygenases. Nine AlkB homologues have been identified, including ALKBH1-8 and FTO. Previous experimental studies have found that these proteins are involved in biological processes such as RNA modification and fatty acid metabolism and the DNA damage response (Wu et al., 2016; Bian et al., 2019; Rajecka et al., 2019). In addition, recent studies have also discovered the potential of the AlkB family to participate in immune responses. ALKBH5 regulates the immune response by controlling CD4+ T cells (Zhou et al., 2021), regulating lactic acid levels (Li et al., 2020), and regulating HMGB1 expression (Chen et al., 2021). Moreover, many studies have reported a role for the AlkB family in the development of BLCA, HNSC, LUAD and OV, and this family is involved in regulating the immune response (Fujii et al., 2013; Pilžys et al., 2019; Cai et al., 2021; Wu et al., 2021), suggesting that AlkB homologues may be promising therapeutic targets.
ALKBH7 is a member of the AlkB family. Multiple studies have shown that ALKBH7 participates in biological processes such as lipid metabolism and programmed necrosis (Wang et al., 2014). ALKBH7 deficiency increases body weight and body fat in mice (Solberg et al., 2013) and protects mouse hearts from ischaemia-reperfusion (IR) injury (Kulkarni et al., 2020). ALKBH7 plays a key role in the process of alkylation and oxidation-induced programmed necrosis (Fu et al., 2013) and drives tissue- and sex-specific necrotic cell death responses (Jordan et al., 2017). In addition, recent studies have identified a role for ALKBH7 in the progression of several cancers and its relationship with immune cell infiltration. ALKBH7 expression is significantly elevated in hepatocellular carcinoma and negatively correlates with CD4+ cells, macrophages and neutrophils (Peng et al., 2021). ALKBH7 is associated with overall survival in individuals with lung adenocarcinoma and negatively correlates with CD8+ T cells and macrophages (Wu et al., 2021). ALKBH7 correlates with the pathological stage of ovarian serous carcinoma and positively correlates with the infiltration of CD8+ T cells, dendritic cells and neutrophils (Cai et al., 2021). ALKBH7 is involved in cellular immunity and the proliferation of HeLa cervical cancer cell lines (Meng et al., 2019).
However, the prognostic value and immunological role of ALKBH7 in cancer have not been systematically investigated. In the present study, we explored changes in the expression and prognostic value of ALKBH7 in 33 cancers. Then, we investigated ALKBH7 expression in different cancer immune and molecular subtypes. In addition, we performed an in-depth study of the immune mechanism of ALKBH7 in different cancers to explore its potential immunotherapeutic value. Overall, this work provides evidence to elucidate the immunotherapeutic role of ALKBH7 in cancer, which may be helpful for further functional experiments.
Materials and Methods
Data Acquisition and Software Availability
The genomic and clinicopathological information, somatic mutation and stemness score data of 33 cancers were obtained from TCGA (https://cancergenome.nih.gov/) and UCSC Xena (https://xena.ucsc.edu/) database (Goldman et al., 2020). In addition, to obtain more normal tissue genomic data, we downloaded tumor and normal tissue gene expression data combined with TCGA and GTEx database on the UCSCXenaShiny (https://shiny.hiplot.com.cn/ucsc-xena-shiny/) website (Wang S. et al., 2021). R 4.1.0 was used to integrate, analyse the original data and visualize the results.
Differential ALKBH7 Expression Analysis in the Normal, Tumor, Various Age, Gender, and Pathological Stage Tissues
The discrepancy of the gene expression between various types of cancer and paired normal tissues was investigated to explore whether ALKBH7 is associated with cancer development. Differential expression analysis of ALKBH7 has been investigated in a variety of cancers with patients’ age, gender and pathological stage by using wilcox test.
Immunohistochemical Staining
The Human Protein Atlas (https://www.proteinatlas.org/) contains over 25,000 antibodies and a collection of over 10 million immunohistochemical (IHC) images (Thul and Lindskog, 2018). To further compare the expression of ALKBH7 gene in tumors and corresponding normal tissues, antibody-based ALKBH7 protein profiles using immunohistochemistry were obtained from the HPA database.
Prognostic Analysis of ALKBH7 in Human Cancers
Univariate Cox regression analysis and Kaplan-Meier curve were used to analyse the relationship between ALKBH7 expression and clinical survival data including overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI) and progression-free interval (PFI) in 33 cancers.
Analysis of ALKBH7 Expression in Different Subtypes of Human Cancers
The TISIDB database (http://cis.hku.hk/TISIDB/) is an online integrated repository portal integrating multiple types of data resources in oncoimmunology (Ru et al., 2019). The relationship between ALKBH7 expression and immune or molecular subtypes of different cancer types was explored through the TISIDB database.
Correlation Analysis of ALKBH7 Expression With Immune Checkpoint Genes, Tumor Mutational Burden, Microsatellite Instability and Tumor Stemness Index in Human Cancers
The correlation between ALKBH7 expression and the expression of immune checkpoint (ICP) genes, was explored via the SangerBox website (http://sangerbox.com/). The tumor mutational burden (TMB), microsatellite instability (MSI) score and tumor stemness index of each TCGA tumor case were obtained from somatic mutation data and previously published studies respectively (Topalian et al., 2015; Bonneville et al., 2017). Tumor stemness indices are indicators for assessing the degree of oncogenic dedifferentiation. Among them, mRNAsi is a gene expression-based stemness index while mDNAsi is a DNA methylation-based stemness index. Correlations between ALKBH7 expression and TMB, MSI, mRNAsi and mDNAsi were analyzed using Spearman’s method.
Analysis of Immune Infiltration-Related Factors and Pathways
The TIMER database (https://cistrome.shinyapps.io/timer/), which collected 10,897 samples across 32 cancer types from TCGA, was created to analyze the level of tumor-associated immune cell infiltration in the TME (Li et al., 2016; Li et al., 2017). The correlation between ALKBH7 expression and six immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages and dendritic cells) and tumor-infiltrating lymphocyte (TIL) marker genes in human cancers was investigated using the TIMER database. ESTIMATE is an algorithm that predicts the presence of immune and stromal cells in tumor tissue which based on gene expression profiles (Yoshihara et al., 2013). We calculated the stromal score, immune score and estimate socre of each case by using the ESTIMATE package. xCell is a powerful web tool for inferring the proportion of immune cell subtypes in tumor tissue (Aran et al., 2017). A spearman correlation heat map of ALKBH7 expression with 36 immunoinfiltrating subtypes of cells in human cancers was established. Finally, to further investigate the relevant signalling pathways, gene set enrichment analysis (GSEA) was performed to explore pathways of ALKBH7 coexpression gene network.
Results
Clinical Landscape of ALKBH7 Expression in 33 Cancers
The details of the analysis are summarized and presented in Figure 1 for a more comprehensive perspective. As illustrated in Figure 2A, significantly higher ALKBH7 expression was detected in most human cancers than in adjacent normal tissues, such as ACC, BRCA, COAD, DLBC, GBM, KICH, KIRP, LGG, LIHC, OV, PAAD, PRAD, READ, SKCM, STAD, THYM and UCEC. In contrast, significantly lower ALKBH7 expression was observed in a few human cancers (ESCA, HNSC, KIRC, LAML and TGCT). ALKBH7 was highly differentially expressed among elderly patients in the THCA, BRCA, KIRP, READ and COAD groups, whereas it was weakly expressed in patients with THYM (Figure 2B). Meanwhile, the results indicated significant sex-based differences in ALKBH7 expression in HNSC, KIRP and LUAD (Figure 2C). In addition, ALKBH7 expression was significantly correlated with the pathological stage of some cancers, including BLCA, KIRC and UCS (Figure 2D). Finally, we used immunohistochemistry to validate ALKBH7 expression. Compared with normal tissues, ALKBH7 was highly expressed in BRCA, LUAD, LUSC, OV, PRAD and UCEC (Figure 3).
FIGURE 1
FIGURE 2
FIGURE 3
Pan-Cancer Analysis of the Multifaceted Prognostic Value of ALKBH7
The association between ALKBH7 expression and patient prognosis was estimated in the pan-cancer dataset. The survival metrics included OS, DSS, DFI, and PFI. Univariate Cox regression analysis of the results from 33 types of cancer suggested that ALKBH7 expression significantly correlated with OS in 4 types of cancer, including BLCA, HNSC, KIRP, and PAAD. Kaplan–Meier survival curves indicated that downregulated ALKBH7 expression was remarkably associated with shorter OS of patients with KIRP, LAML, MESO, SARC, and UCEC (Figure 4). The relationship between ALKBH7 expression and DSS in patients with cancer was examined. ALKBH7 expression affected DSS in six types of cancer, including BLCA, KIRP, LIHC, LUSC, PAAD, and PCPG. The Kaplan–Meier analysis indicated that decreased ALKBH7 expression indicated shorter DSS of patients with BLCA, KIRP, MESO, and UCEC, while increased ALKBH7 expression corresponded with shorter DSS of patients with KIRC (Figure 5). Cox regression analysis of the DFI revealed that ALKBH7 expression significantly correlated with DFI in 4 types of cancer, including LUSC, OV, PAAD, and THCA. The results from the Kaplan–Meier analysis suggested that increased ALKBH7 expression was associated with a poor prognosis for patients with PRAD, while decreased ALKBH7 expression was associated with a poor prognosis for patients with THCA (Supplementary Figure S1). We also assessed the association between ALKBH7 expression and PFI and identified that ALKBH7 expression influenced PFI in patients with BLCA, KIRC, LUSC and PAAD. Kaplan–Meier PFI curves revealed that decreased ALKBH7 mRNA expression correlated with an unfavourable PFI in patients with BLCA and PAAD. In contrast, increased ALKBH7 mRNA expression correlated with an unfavourable PFI in patients with KIRC (Supplementary Figure S2).
FIGURE 4
FIGURE 5
ALKBH7 Expression Is Related to Immune and Molecular Subtypes in Human Cancers
Based on accumulating evidence, immunophenotyping reflects the comprehensive immune status of a tumor, which is closely related to immunotherapy and the tumor microenvironment (Ma et al., 2021). Different molecular subtypes correspond to the unique molecular biology of cancer and may facilitate the selection of molecular targeted therapies and immunotherapy strategies (Kim et al., 2019; Bai et al., 2021). Next, ALKBH7 expression in immune and molecular subtypes of human cancer was explored using the TISIDB website. ALKBH7 expression was significantly different in different immune subtypes of BLCA, BRCA, KIRC, LIHC, PRAD, SKCM, TGCT, and UCEC (Figure 6). In addition, the trends for the up- and downregulation of ALKBH7 expression were also different in different immune subtypes of a specific cancer type. Taking SKCM as an example, low ALKBH7 expression was detected in C2 and C4 types and high expression was observed in the C3 type. Regarding different molecular subtypes of cancers, a significant correlation with ALKBH7 expression was observed in BRCA, COAD, HNSC, KIRP, LGG, LUSC, OV, PRAD, STAD and UCEC (Figure 7). Based on the results described above, we suggest that ALKBH7 may play an important role in the tumor immune microenvironment and modulate the effect of immunotherapy.
FIGURE 6
FIGURE 7
ALKBH7 Expression Is Related to Immune Checkpoint Genes in Human Cancers
Studies have shown that immune checkpoint genes have important implications for immunotherapy in many cancers (Topalian et al., 2015; Li et al., 2019). Here, we collected expression patterns of 47 common immune checkpoint genes and analysed the relationship between ALKBH7 expression and immune checkpoint gene expression to explore the potential role of ALKBH7 in immunotherapy. As shown in Figure 8, ALKBH7 expression significantly correlated with the expression of most ICP genes in many cancers, such as BRCA COAD, HNSC, KIRC, LUAD, OV, PAAD, PRAD, READ, SKCM, THCA, THYM, and UVM. Among them, a negative correlation was the main trend; for example, in PRAD, ALKBH7 expression was negatively correlated with the expression of 30 ICP genes and positively correlated with the expression of 5 ICP genes. Thus, high levels of ALKBH7 expression may predict unsatisfactory immunotherapy outcomes when targeting ICP genes. On the other hand, ALKBH7 inhibitors may be potential alternative therapeutic approaches. Therefore, we hypothesized that ALKBH7, a potential pan-cancer biomarker or a novel immunotherapeutic target, may predict the response to immunotherapy or achieve promising therapeutic outcomes, respectively.
FIGURE 8
ALKBH7 Expression Is Related to the Tumor Mutational Burden, Microsatellite Instability, and Tumor Stemness Index
We analysed the correlations between ALKBH7 expression and the TMB, MSI, and tumor stemness index to explore the role of ALKBH7 in the immune mechanism and immune response of the tumor microenvironment (TME). The TMB, MSI, and tumor stemness index in the tumor microenvironment are related to antitumor immunity and might predict the therapeutic efficacy of tumor immunotherapy (Lee et al., 2016; Yarchoan et al., 2017; Malta et al., 2018). As presented in Figure 9, ALKBH7 was associated with the TMB in 7 cancers and MSI in 13 cancers. In addition, ALKBH7 was related to mDNAsi in 12 cancers and mRNAsi in 13 cancers. Among them, ALKBH7 expression was negatively correlated with the TMB and MSI in COAD and READ, while it was positively correlated with the TMB and MSI in UCEC. Based on this finding, ALKBH7 might exert an indirect effect on the immunotherapeutic response of COAD, READ and UCEC. ALKBH7 was positively correlated with mRNAsi and mDNAsi in TGCT and HNSC, but negatively correlated with mRNAsi and mDNAsi in BRCA. High ALKBH7 expression in TGCT and HNSC may be related to the low sensitivity to immune checkpoint blockade therapy; in contrast, high ALKBH7 expression in BRCA may be related to the high sensitivity to immune checkpoint blockade therapy. Interestingly, ALKBH7 was positively correlated with mRNAsi but negatively correlated with mDNAsi in THCA and THYM. This result might arise from the discrepancies between mRNAsi and mDNAsi caused by DNA hypermethylation (Malta et al., 2018).
FIGURE 9
Correlation Analysis Between ALKBH7 Expression and Infiltrating Immune Cells and the ESTIMATE Score
The tumor microenvironment contains immune cells and fibroblasts, which affect the effect of immunotherapy (Aran et al., 2015). We analysed the correlation between ALKBH7 expression and six types of infiltrating immune cells, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells. The results revealed a significant correlation in 31 cancer types. ALKBH7 expression displayed a strong relationship with dendritic cells in 8 cancer types, macrophages in 9 cancer types, neutrophils in 11 cancer types, CD8+ T cells in 14 cancer types, B cells in 8 cancer types and CD4+ T cells in 3 cancer types (Figure 10A). Results from the TIMER database included these results, and all details are shown in Supplementary Table S1. Subsequently, the correlation between ALKBH7 expression and stromal, immune and ESTIMATE scores was analysed (Figure 10B). All results are presented in Supplementary Table S2. Interestingly, the most significant correlation between ALKBH7 expression and the two parameters described above was observed in PAAD, PRAD and THCA. As shown in Figures 10C,D, ALKBH7 expression was negatively correlated with TILs and stromal, immune, and ESTIMATE scores. Therefore, ALKBH7 may be involved in inhibiting immune cell infiltration in PAAD, PRAD and THCA.
FIGURE 10
We also used the xCell web tool to explore the association between ALKBH7 gene expression and the infiltration of various subtypes of immune cells. The Spearman correlation heat map is shown in Figure 11. NK T cells and CD4+ Th1 T cells were positively correlated with ALKBH7 gene expression in most cancers. In contrast, CD4+ Th2 T cells, memory CD4+ T cells, monocytes and mast cells were negatively correlated with ALKBH7 gene expression in most cancers. In PAAD, PRAD and THCA, ALKBH7 expression was associated with most subtypes of immune cells and generally exhibited negative correlations.
FIGURE 11
Correlation Between ALKBH7 Expression and Various Immune Markers
We validated the correlations between ALKBH7 expression and diverse immune signatures in PAAD, PRAD and THCA using the TIMER database to obtain a better understanding of ALKBH7 crosstalk with the immune response. The genes listed in Table 1 were used to characterize immune cells, including CD8+ T cells, T cells, B cells, monocytes, tumor-associated macrophages (TAMs), M1 macrophages, M2 macrophages, neutrophils and dendritic cells. Tumor purity is an important aspect affecting the number of infiltrating immune cells in clinical cancer biopsies. After adjusting for tumor purity, ALKBH7 expression was significantly negatively correlated with most of the immune markers of divergent types of immune cells in PAAD, PRAD and THCA (Table 1).
TABLE 1
| Description | Gene markers | THCA | PRAD | PAAD | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| None | Purity | None | Purity | None | Purity | ||||||||
| Cor | Pvalue | Cor | Pvalue | Cor | Pvalue | Cor | Pvalue | Cor | Pvalue | Cor | Pvalue | ||
| CD8+ T cell | CD8A | −0.202 | a | −0.192 | a | −0.314 | a | −0.254 | a | −0.292 | a | −0.233 | b |
| CD8B | −0.152 | a | −0.143 | b | −0.157 | a | −0.124 | c | −0.279 | a | −0.217 | b | |
| T cell (general) | CD3D | −0.26 | a | −0.25 | a | −0.153 | a | −0.071 | 0.146 | −0.19 | c | −0.122 | 0.112 |
| CD3E | −0.305 | a | −0.295 | a | −0.259 | a | −0.2 | a | −0.238 | b | −0.17 | c | |
| CD2 | −0.316 | a | −0.303 | a | −0.253 | a | −0.16 | b | −0.292 | a | −0.224 | b | |
| B cell | CD19 | −0.183 | a | −0.167 | a | −0.059 | 0.187 | −0.015 | 0.759 | −0.186 | c | −0.118 | 0.123 |
| CD79A | −0.209 | a | −0.194 | a | −0.153 | a | −0.089 | 0.068 | −0.194 | b | −0.123 | 0.109 | |
| Monocyte | CD86 | −0.404 | a | −0.393 | a | −0.414 | a | −0.336 | a | −0.399 | a | −0.332 | a |
| CSF1R | −0.346 | a | −0.343 | a | −0.37 | a | −0.299 | a | −0.347 | a | −0.289 | a | |
| TAM | CCL2 | −0.26 | a | −0.244 | a | −0.223 | a | −0.165 | a | −0.282 | a | −0.253 | a |
| CD68 | −0.358 | a | −0.34 | a | −0.35 | a | −0.289 | a | −0.228 | b | −0.157 | c | |
| IL10 | −0.292 | a | −0.279 | a | −0.338 | a | −0.251 | a | −0.258 | a | −0.208 | b | |
| M1 Macrophage | IRF5 | −0.324 | a | −0.317 | a | −0.234 | a | −0.263 | a | 0.027 | 0.719 | 0.058 | 0.454 |
| PTGS2 | −0.319 | a | −0.31 | a | −0.268 | a | −0.179 | a | −0.237 | b | −0.253 | a | |
| M2 Macrophage | CD163 | −0.335 | a | −0.32 | a | −0.414 | a | −0.343 | a | −0.438 | a | −0.377 | a |
| VSIG4 | −0.353 | a | −0.345 | a | −0.379 | a | −0.303 | a | −0.347 | a | −0.278 | a | |
| MS4A4A | −0.334 | a | −0.321 | a | −0.363 | a | −0.297 | a | −0.379 | a | −0.309 | a | |
| Neutrophils | CEACAM8 | −0.196 | a | −0.197 | a | 0.004 | 0.922 | 0.017 | 0.724 | −0.059 | 0.433 | −0.005 | 0.95 |
| ITGAM | −0.411 | a | −0.401 | a | −0.361 | a | −0.309 | a | −0.285 | a | −0.191 | c | |
| CCR7 | −0.293 | a | −0.276 | a | −0.274 | a | −0.217 | a | −0.178 | c | −0.124 | 0.106 | |
| Dendritic cell | HLA-DPB1 | −0.305 | a | −0.295 | a | −0.119 | b | −0.057 | 0.246 | −0.232 | b | −0.159 | c |
| HLA-DQB1 | −0.291 | a | −0.295 | a | −0.198 | a | −0.147 | b | −0.243 | b | −0.189 | c | |
| HLA-DRA | −0.377 | a | −0.365 | a | −0.351 | a | −0.295 | a | −0.35 | a | −0.289 | a | |
| HLA-DPA1 | −0.364 | a | −0.351 | a | −0.339 | a | −0.257 | a | −0.356 | a | −0.3 | a | |
| CD1C | −0.336 | a | −0.319 | a | −0.319 | a | −0.242 | a | −0.199 | b | −0.14 | 0.067 | |
| NRP1 | −0.258 | a | −0.245 | a | −0.29 | a | −0.271 | a | −0.471 | a | −0.439 | a | |
| ITGAX | −0.363 | a | −0.348 | a | −0.299 | a | −0.277 | a | −0.185 | c | −0.091 | 0.236 | |
Correlation analysis between ALKBH7 and gene markers of immune cells in TIMER.
p < 0.05.
p < 0.01.
p < 0.001.
We also examined the correlations between ALKBH7 expression and various functional T cells, including Th1, Th1-like, Th2, Th17, Tfh, Treg, resting Tregs, effector Tregs, effector T cells, naïve T cells, effector memory T cells, resistant memory T cells, and exhausted T cells (Table 2). Using the TIMER database, the ALKBH7 expression level was also significantly negatively correlated with 44 of 50 T cell markers in PRAD and THCA and with 31 of 50 T cell markers in PAAD after adjusting for tumor purity (Table 2). These findings further support the hypothesis that ALKBH7 may be involved in inhibiting immune cell infiltration in PAAD, PRAD and THCA.
TABLE 2
| Description | Gene markers | THCA | PRAD | PAAD | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| None | Purity | None | Purity | None | Purity | ||||||||
| Cor | P | Cor | P | Cor | P | Cor | P | Cor | P | Cor | P | ||
| Th1 | TBX21 | −0.199 | a | −0.192 | a | −0.2 | a | −0.17 | a | −0.194 | b | −0.143 | 0.062 |
| STAT4 | −0.317 | a | −0.317 | a | −0.281 | a | −0.216 | a | −0.207 | b | −0.179 | c | |
| STAT1 | −0.451 | a | −0.436 | a | −0.421 | a | −0.339 | a | −0.355 | a | −0.302 | a | |
| IFNG | −0.231 | a | −0.221 | a | −0.206 | a | −0.144 | b | −0.277 | a | −0.238 | b | |
| TNF | −0.236 | a | −0.229 | a | −0.261 | a | −0.168 | a | −0.125 | 0.096 | −0.076 | 0.32 | |
| IL12A | −0.064 | 0.152 | −0.057 | 0.206 | −0.21 | a | −0.155 | b | −0.127 | 0.09 | −0.108 | 0.159 | |
| IL12B | −0.254 | a | −0.245 | a | −0.211 | a | −0.145 | b | −0.134 | 0.075 | −0.101 | 0.187 | |
| Th1-like | HAVCR2 | −0.395 | a | −0.382 | a | −0.364 | a | −0.3 | a | −0.352 | a | −0.28 | a |
| IFNG | −0.231 | a | −0.221 | a | −0.206 | a | −0.144 | b | −0.277 | a | −0.238 | b | |
| CXCR3 | −0.15 | a | −0.136 | a | −0.226 | a | −0.185 | a | −0.017 | 0.825 | 0.046 | 0.549 | |
| BHLHE40 | −0.375 | a | −0.367 | a | −0.345 | a | −0.311 | a | −0.097 | 0.196 | −0.091 | 0.238 | |
| CD4 | −0.394 | a | −0.384 | a | −0.4 | a | −0.322 | a | −0.337 | a | −0.268 | a | |
| Th2 | GATA3 | −0.077 | 0.082 | −0.058 | 0.199 | −0.127 | b | −0.035 | 0.479 | −0.157 | c | −0.119 | 0.122 |
| STAT6 | −0.301 | a | −0.285 | a | −0.308 | a | −0.269 | a | −0.034 | 0.65 | −0.019 | 0.806 | |
| STAT5A | −0.273 | a | −0.266 | a | −0.24 | a | −0.164 | a | −0.072 | 0.339 | −0.006 | 0.942 | |
| IL13 | −0.074 | 0.094 | −0.071 | 0.116 | −0.061 | 0.177 | −0.088 | 0.073 | −0.058 | 0.441 | −0.049 | 0.521 | |
| Th17 | STAT3 | −0.356 | a | −0.336 | a | −0.449 | a | −0.385 | a | −0.338 | a | −0.298 | a |
| IL17A | −0.135 | a | −0.13 | b | −0.139 | b | −0.028 | 0.565 | −0.235 | b | −0.228 | b | |
| Tfh | BCL6 | −0.231 | a | −0.202 | a | −0.316 | a | −0.306 | a | −0.26 | a | −0.238 | b |
| IL21 | −0.144 | b | −0.133 | b | −0.139 | b | −0.116 | c | −0.089 | 0.234 | −0.048 | 0.532 | |
| Treg | FOXP3 | −0.377 | a | −0.363 | a | −0.332 | a | −0.337 | a | −0.264 | a | −0.199 | b |
| CCR8 | −0.421 | a | −0.402 | a | −0.451 | a | −0.396 | a | −0.392 | a | −0.342 | a | |
| STAT5B | −0.209 | a | −0.189 | a | −0.465 | a | −0.409 | a | −0.126 | 0.093 | −0.127 | 0.098 | |
| TGFB1 | −0.05 | 0.264 | −0.042 | 0.359 | −0.179 | a | −0.193 | a | 0.07 | 0.353 | 0.133 | 0.083 | |
| Resting Treg | FOXP3 | −0.377 | a | −0.363 | a | −0.332 | a | −0.337 | a | −0.264 | a | −0.199 | b |
| IL2RA | −0.414 | a | −0.404 | a | −0.414 | a | −0.356 | a | −0.397 | a | −0.334 | a | |
| Effector Treg | FOXP3 | −0.377 | a | −0.363 | a | −0.332 | a | −0.337 | a | −0.264 | a | −0.199 | b |
| T-cell | CCR8 | −0.421 | a | −0.402 | a | −0.451 | a | −0.396 | a | −0.392 | a | −0.342 | a |
| TNFRSF9 | −0.38 | a | −0.361 | a | −0.454 | a | −0.377 | a | −0.362 | a | −0.309 | a | |
| Effector | CX3CR1 | −0.17 | a | −0.162 | a | −0.383 | a | −0.262 | a | −0.262 | a | −0.236 | b |
| T-cell | FGFBP2 | 0.055 | 0.217 | 0.051 | 0.259 | −0.118 | b | −0.094 | 0.056 | −0.22 | b | −0.196 | c |
| FCGR3A | −0.336 | a | −0.329 | a | −0.38 | a | −0.317 | a | −0.388 | a | −0.326 | a | |
| Naive T-cell | CCR7 | −0.293 | a | −0.276 | a | −0.274 | a | −0.217 | a | −0.178 | c | −0.124 | 0.106 |
| SELL | −0.355 | a | −0.352 | a | −0.4 | a | −0.331 | a | −0.249 | a | −0.183 | c | |
| Effector memory | DUSP4 | −0.29 | a | −0.278 | a | −0.162 | a | −0.153 | a | −0.111 | 0.137 | 0.09 | 0.242 |
| T-cell | GZMK | −0.253 | a | −0.242 | a | −0.258 | a | −0.19 | a | −0.216 | b | −0.149 | 0.052 |
| GZMA | −0.242 | a | −0.235 | a | −0.228 | a | −0.151 | b | −0.222 | b | −0.173 | c | |
| Resident memory | CD69 | −0.324 | a | −0.311 | a | −0.389 | a | −0.306 | a | −0.335 | a | −0.297 | a |
| T-cell | CXCR6 | −0.282 | a | −0.269 | a | −0.292 | a | −0.188 | a | −0.339 | a | −0.283 | a |
| MYADM | −0.144 | a | −0.123 | b | −0.335 | a | −0.306 | a | −0.186 | c | −0.159 | c | |
| General memory | CCR7 | −0.293 | a | −0.276 | a | −0.274 | a | −0.217 | a | −0.178 | c | −0.124 | 0.106 |
| T-cell | SELL | −0.355 | a | −0.352 | a | −0.4 | a | −0.331 | a | −0.249 | a | −0.183 | c |
| IL7R | −0.415 | a | −0.4 | a | −0.458 | a | −0.394 | a | −0.425 | a | −0.38 | a | |
| Exhausted T cell | PDCD1 | −0.152 | a | −0.154 | a | −0.088 | c | −0.072 | 0.142 | −0.117 | 0.119 | −0.044 | 0.568 |
| CTLA4 | −0.317 | a | −0.304 | a | −0.151 | a | −0.115 | c | −0.227 | b | −0.159 | c | |
| LAG3 | −0.211 | a | −0.207 | a | −0.031 | 0.486 | 0.001 | 0.98 | −0.074 | 0.323 | −0.047 | 0.543 | |
| HAVCR2 | −0.395 | a | −0.382 | a | −0.364 | a | −0.3 | a | −0.352 | a | −0.28 | a | |
| GZMB | −0.231 | a | −0.23 | a | −0.162 | a | −0.114 | c | −0.33 | a | −0.273 | a | |
| CXCL13 | −0.264 | a | −0.262 | a | −0.192 | a | −0.147 | b | −0.218 | b | −0.161 | c | |
| LAYN | −0.073 | 0.098 | −0.074 | 0.101 | −0.237 | a | −0.189 | a | −0.217 | b | −0.162 | c | |
Correlation analysis between ALKBH7 and gene markers of different types of T cells in TIMER.
p < 0.05.
p < 0.01.
p < 0.001.
The ALKBH7 Coexpression Network Relevant Signalling Pathways
The aforementioned results identified significant associations between ALKBH7 expression and the prognosis and immunity of cancers. Considering the robust correlation between ALKBH7 expression and PAAD, PRAD and THCA, GSEA was performed to investigate the potential signalling pathways of ALKBH7 in these cancers. The results presented in Figure 12 indicate that genes coexpressed with ALKBH7 are enriched in the regulation of immune and inflammatory responses and are negatively associated with these pathways, such as the JAK/STAT signalling pathway and TGF-β signalling pathway. These results suggested that ALKBH7 expression might play an essential role in human cancers by suppressing the immune response of the TME.
FIGURE 12
Discussion
ALKBH7 is a mitochondrial protein involved in programmed necrosis, fatty acid metabolism and obesity development. In this paper, a comprehensive pan-cancer study of ALKBH7 revealed the potential prognostic and immunotherapeutic value of ALKBH7 in human cancers. First, significantly higher ALKBH7 expression was detected in most cancers compared to paired normal tissues, consistent with previous studies. Cai et al. observed high ALKBH7 expression in ovarian plasmacytoma (Cai et al., 2021), and Peng et al. found high ALKBH7 expression in hepatocellular carcinoma (Peng et al., 2021). However, ALKBH7 expression correlates with clinical parameters (age, sex and pathological stage) in only a few patients with cancer. For example, in BLCA, ALKBH7 expression correlated with the pathological stage of the tumor. Interestingly, ALKBH7 expression has some prognostic value for some cancers. For example, in a univariate survival analysis, ALKBH7 expression was significantly associated with four clinical survival datasets (OS, DSS, DFI and PFI) in patients with PAAD; in Kaplan–Meier survival estimates, downregulated ALKBH7 expression was significantly associated with shorter OS and DSS of patients with UCEC. These results suggest that ALKBH7 is a potential prognostic biomarker.
Next, ALKBH7 expression in different immune subtypes and molecular subtypes of human cancers was explored to determine its potential mechanism of action. ALKBH7 expression was significantly different in different immune subtypes and molecular subtypes in many cancer types, suggesting that ALKBH7 is a promising diagnostic pan-cancer biomarker and participates in immune regulation. Moreover, we documented significant differences in ALKBH7 expression in different immune and molecular subtypes of BRCA, PRAD and UCEC. In fact, differential ALKBH7 expression was detected in the cancers listed above and their normal tissue, indicating that ALKBH7 might play a role in the growth and progression of cancers.
Tumor cells use the immune checkpoint pathway to suppress immune cells and achieve immune escape (Topalian et al., 2015). Based on this principle, immune checkpoint inhibitors (ICIs) have emerged as new therapeutic approaches for cancer treatment (Muenst et al., 2016) and have been successfully applied in the clinic. The most commonly used ICI predictive biomarkers are programmed cell death ligand-1 (PD-L1), microsatellite instability (MSI) and tumor mutational burden (TMB) (Wang Y. et al., 2021). In addition, a study by Malta et al. found that a high tumor stemness index was associated with reduced PD-L1 expression in most cancers (Malta et al., 2018). In the present study, immunotherapy biomarkers (TMB and MSI) and the tumor stemness index showed significant associations with ALKBH7 in some cancers. Moreover, a strong relationship between the expression of ALKBH7 and ICP genes was identified. These results indicate that ALKBH7 has a strong association with ICIs.
Based on accumulating evidence, the tumor microenvironment (TME) is involved in tumor progression and significantly affects the treatment response and clinical outcome (Wu and Dai, 2017; Hinshaw and Shevde, 2019). Tumor-infiltrating lymphocytes (TILs) in the TME have been proven to be an independent predictor of the prognosis of patients with cancer and immunotherapeutic efficacy (Azimi et al., 2012). Our study found that ALKBH7 was related to the immune, stromal, and ESTIMATE scores and immune cell infiltration in the TME of most human cancer types, especially in PAAD, PRAD and THCA. Then, we explored the function of ALKBH7 in PAAD, PRAD and THCA by performing a KEGG analysis. ALKBH7 and its coexpression network were indeed involved in the regulation of the immune response and inflammatory response. In summary, these results strongly indicated the potential of ALKBH7 as a target of anticancer immunotherapy.
Overall, our pan-cancer analysis of ALKBH7 is the first to explore the relationship between ALKBH7 expression in human cancers and clinical prognostic factors, immune subtypes, molecular subtypes, immune checkpoints (ICPs), tumor mutational burden (TMB), microsatellite instability (MSI), tumor stemness index, tumor microenvironment (TME) and tumor-infiltrating lymphocytes (TILs). This information contributes to the understanding of the function of ALKBH7 in cancer development and its role in immunology. However, more experimental studies are required to explore the specific mechanisms of ALKBH7 action in cancer.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repositories and accession numbers can be found in the article/Supplementary Material.
Author contributions
XG, CL, and YG conceived and designed the study. KC and DS were responsible for the collection and assembly of data, data analysis, and interpretation. KC, DL, and YH were involved in writing of the manuscript. LT provided help in revising the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (81772829 and 81830052), the Special Program for Collaborative Innovation, the Construction project of Shanghai Key Laboratory of Molecular Imaging (18DZ2260400) and Hunan Province National Science Foundation of China (2018JJ3895), “Top-100 Talent Cultivation Plan” of Shanghai University of Medicine and Health Sciences.
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, orclaim 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/fgene.2022.822261/full#supplementary-material
References
1
AranD.HuZ.ButteA. J. (2017). xCell: Digitally Portraying the Tissue Cellular Heterogeneity Landscape. Genome Biol.18, 220. 10.1186/s13059-017-1349-1
2
AranD.SirotaM.ButteA. J. (2015). Systematic Pan-Cancer Analysis of Tumour Purity. Nat. Commun.6, 8971. 10.1038/ncomms9971
3
AzimiF.ScolyerR. A.RumchevaP.MoncrieffM.MuraliR.McCarthyS. W.et al (2012). Tumor-infiltrating Lymphocyte Grade Is an Independent Predictor of sentinel Lymph Node Status and Survival in Patients with Cutaneous Melanoma. Jco30, 2678–2683. 10.1200/JCO.2011.37.8539
4
BaiR.LiL.LiL.ChenX.ZhaoY.SongW.et al (2021). Advances in Novel Molecular Typing and Precise Treatment Strategies for Small Cell Lung Cancer. Chin. J. Cancer Res.33, 522–534. 10.21147/j.issn.1000-9604.2021.04.09
5
BianK.LenzS. A. P.TangQ.ChenF.QiR.JostM.et al (2019). DNA Repair Enzymes ALKBH2, ALKBH3, and AlkB Oxidize 5-methylcytosine to 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine In Vitro. Nucleic Acids Res.47, 5522–5529. 10.1093/nar/gkz395
6
BonnevilleR.KrookM. A.KauttoE. A.MiyaJ.WingM. R.ChenH.-Z.et al (2017). Landscape of Microsatellite Instability across 39 Cancer Types. JCO Precision Oncol.PO.17, 1–15. 10.1200/PO.17.00073
7
CaiY.WuG.PengB.LiJ.ZengS.YanY.et al (2021). Expression and Molecular Profiles of the AlkB Family in Ovarian Serous Carcinoma. Aging13, 9679–9692. 10.18632/aging.202716
8
ChenG.ZhaoQ.YuanB.WangB.ZhangY.LiZ.et al (2021). ALKBH5-Modified HMGB1-STING Activation Contributes to Radiation Induced Liver Disease via Innate Immune Response. Int. J. Radiat. Oncology*Biology*Physics111, 491–501. 10.1016/j.ijrobp.2021.05.115
9
FuD.JordanJ. J.SamsonL. D. (2013). Human ALKBH7 Is Required for Alkylation and Oxidation-Induced Programmed Necrosis. Genes Dev.27, 1089–1100. 10.1101/gad.215533.113
10
FujiiT.ShimadaK.AnaiS.FujimotoK.KonishiN. (2013). ALKBH2, a Novel AlkB Homologue, Contributes to Human Bladder Cancer Progression by Regulating MUC1 Expression. Cancer Sci.104, 321–327. 10.1111/cas.12089
11
GoldmanM. J.CraftB.HastieM.RepečkaK.McDadeF.KamathA.et al (2020). Visualizing and Interpreting Cancer Genomics Data via the Xena Platform. Nat. Biotechnol.38, 675–678. 10.1038/s41587-020-0546-8
12
HinshawD. C.ShevdeL. A. (2019). The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res.79, 4557–4566. 10.1158/0008-5472.CAN-18-3962
13
JordanJ. J.ChhimS.MarguliesC. M.AlloccaM.BronsonR. T.KlunglandA.et al (2017). ALKBH7 Drives a Tissue and Sex-specific Necrotic Cell Death Response Following Alkylation-Induced Damage. Cell Death Dis8, e2947. 10.1038/cddis.2017.343
14
KimT. S.da SilvaE.CoitD. G.TangL. H. (2019). Intratumoral Immune Response to Gastric Cancer Varies by Molecular and Histologic Subtype. Am. J. Surg. Pathol.43, 851–860. 10.1097/PAS.0000000000001253
15
KulkarniC. A.NadtochiyS. M.KennedyL.ZhangJ.ChhimS.AlwaseemH.et al (2020). ALKBH7 Mediates Necrosis via Rewiring of Glyoxal Metabolism. Elife9, e58573. 10.7554/eLife.58573
16
LeeV.MurphyA.LeD. T.DiazL. A.Jr (2016). Mismatch Repair Deficiency and Response to Immune Checkpoint Blockade. Oncologist21, 1200–1211. 10.1634/theoncologist.2016-0046
17
LiB.ChanH. L.ChenP. (2019). Immune Checkpoint Inhibitors: Basics and Challenges. Cmc26, 3009–3025. 10.2174/0929867324666170804143706
18
LiB.SeversonE.PignonJ.-C.ZhaoH.LiT.NovakJ.et al (2016). Comprehensive Analyses of Tumor Immunity: Implications for Cancer Immunotherapy. Genome Biol.17, 174. 10.1186/s13059-016-1028-7
19
LiN.KangY.WangL.HuffS.TangR.HuiH.et al (2020). ALKBH5 Regulates Anti-PD-1 Therapy Response by Modulating Lactate and Suppressive Immune Cell Accumulation in Tumor Microenvironment. Proc. Natl. Acad. Sci. USA117, 20159–20170. 10.1073/pnas.1918986117
20
LiT.FanJ.WangB.TraughN.ChenQ.LiuJ. S.et al (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res.77, e108–e110. 10.1158/0008-5472.CAN-17-0307
21
MaS.BaY.JiH.WangF.DuJ.HuS. (2021). Recognition of Tumor-Associated Antigens and Immune Subtypes in Glioma for mRNA Vaccine Development. Front. Immunol.12, 738435. 10.3389/fimmu.2021.738435
22
MaltaT. M.SokolovA.GentlesA. J.BurzykowskiT.PoissonL.WeinsteinJ. N.et al (2018). Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation. Cell173, 338–e15. e15. 10.1016/j.cell.2018.03.034
23
MengS.ZhanS.DouW.GeW. (2019). The Interactome and Proteomic Responses of ALKBH7 in Cell Lines by In-Depth Proteomics Analysis. Proteome Sci.17, 8. 10.1186/s12953-019-0156-x
24
MuenstS.LäubliH.SoysalS. D.ZippeliusA.TzankovA.HoellerS. (2016). The Immune System and Cancer Evasion Strategies: Therapeutic Concepts. J. Intern. Med.279, 541–562. 10.1111/joim.12470
25
PengB.YanY.XuZ. (2021). The Bioinformatics and Experimental Analysis of AlkB Family for Prognosis and Immune Cell Infiltration in Hepatocellular Carcinoma. PeerJ9, e12123. 10.7717/peerj.12123
26
PilžysT.MarcinkowskiM.KukwaW.GarbiczD.DylewskaM.FerencK.et al (2019). ALKBH Overexpression in Head and Neck Cancer: Potential Target for Novel Anticancer Therapy. Sci. Rep.9, 13249. 10.1038/s41598-019-49550-x
27
RajeckaV.SkalickyT.VanacovaS. (2019). The Role of RNA Adenosine Demethylases in the Control of Gene Expression. Biochim. Biophys. Acta (Bba) - Gene Regul. Mech.1862, 343–355. 10.1016/j.bbagrm.2018.12.001
28
RuB.WongC. N.TongY.ZhongJ. Y.ZhongS. S. W.WuW. C.et al (2019). TISIDB: an Integrated Repository portal for Tumor-Immune System Interactions. Bioinformatics35, 4200–4202. 10.1093/bioinformatics/btz210
29
SolbergA.RobertsonA. B.AronsenJ. M.RognmoO.SjaastadI.WisloffU.et al (2013). Deletion of Mouse Alkbh7 Leads to Obesity. J. Mol. Cel Biol.5, 194–203. 10.1093/jmcb/mjt012
30
ThulP. J.LindskogC. (2018). The Human Protein Atlas: A Spatial Map of the Human Proteome. Protein Sci.27, 233–244. 10.1002/pro.3307
31
TopalianS. L.DrakeC. G.PardollD. M. (2015). Immune Checkpoint Blockade: a Common Denominator Approach to Cancer Therapy. Cancer Cell27, 450–461. 10.1016/j.ccell.2015.03.001
32
WangG.HeQ.FengC.LiuY.DengZ.QiX.et al (2014). The Atomic Resolution Structure of Human AlkB Homolog 7 (ALKBH7), a Key Protein for Programmed Necrosis and Fat Metabolism. J. Biol. Chem.289, 27924–27936. 10.1074/jbc.M114.590505
33
WangS.XiongY.ZhaoL.GuK.LiY.ZhaoF.et al (2021). UCSCXenaShiny: An R/CRAN Package for Interactive Analysis of UCSC Xena Data. Bioinformatics38, 527–529. 10.1093/bioinformatics/btab561
34
WangY.TongZ.ZhangW.ZhangW.BuzdinA.MuX.et al (2021). FDA-approved and Emerging Next Generation Predictive Biomarkers for Immune Checkpoint Inhibitors in Cancer Patients. Front. Oncol.11, 683419. 10.3389/fonc.2021.683419
35
WuG.YanY.CaiY.PengB.LiJ.HuangJ.et al (2021). ALKBH1-8 and FTO: Potential Therapeutic Targets and Prognostic Biomarkers in Lung Adenocarcinoma Pathogenesis. Front. Cel Dev. Biol.9, 633927. 10.3389/fcell.2021.633927
36
WuT.DaiY. (2017). Tumor Microenvironment and Therapeutic Response. Cancer Lett.387, 61–68. 10.1016/j.canlet.2016.01.043
37
WuT. P.WangT.SeetinM. G.LaiY.ZhuS.LinK.et al (2016). DNA Methylation on N6-Adenine in Mammalian Embryonic Stem Cells. Nature532, 329–333. 10.1038/nature17640
38
YarchoanM.HopkinsA.JaffeeE. M. (2017). Tumor Mutational Burden and Response Rate to PD-1 Inhibition. N. Engl. J. Med.377, 2500–2501. 10.1056/NEJMc1713444
39
YoshiharaK.ShahmoradgoliM.MartínezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun.4, 2612. 10.1038/ncomms3612
40
ZhouJ.ZhangX.HuJ.QuR.YuZ.XuH.et al (2021). m 6 A Demethylase ALKBH5 Controls CD4 + T Cell Pathogenicity and Promotes Autoimmunity. Sci. Adv.7, eabg0470. 10.1126/sciadv.abg0470
Summary
Keywords
ALKBH7, pan-cancer, prognosis, immunotherapy, immune infiltration
Citation
Chen K, Shen D, Tan L, Lai D, Han Y, Gu Y, Lu C and Gu X (2022) A Pan-Cancer Analysis Reveals the Prognostic and Immunotherapeutic Value of ALKBH7. Front. Genet. 13:822261. doi: 10.3389/fgene.2022.822261
Received
02 November 2021
Accepted
14 January 2022
Published
11 February 2022
Volume
13 - 2022
Edited by
Guiyou Liu, Tianjin Institute of Industrial Biotechnology (CAS), China
Reviewed by
Li Hongdong, Gannan Medical University, China
Hailin Zhao, Imperial College London, United Kingdom
Updates
Copyright
© 2022 Chen, Shen, Tan, Lai, Han, Gu, Lu and Gu.
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: Xuefeng Gu, guxf@sumhs.edu.cn; Changlian Lu, lvcl@sumhs.edu.cn; Yonggang Gu, gyg0217@qq.com
†These authors have contributed equally to this work and share first authorship
This article was submitted to Statistical Genetics and Methodology, a section of the journal Frontiers in Genetics
Disclaimer
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.