ALKBH1-8 and FTO: Potential Therapeutic Targets and Prognostic Biomarkers in Lung Adenocarcinoma Pathogenesis
Geting Wu1† 
Yuanliang Yan2† 
Yuan Cai1 
Bi Peng1 Juanni Li1 Jinzhou Huang3 
Zhijie Xu1,4* 
Jianhua Zhou1*- 1Department of Pathology, Xiangya Hospital, Central South University, Changsha, China
- 2Department of Pharmacy, Xiangya Hospital, Central South University, Changsha, China
- 3Department of Oncology, Mayo Clinic, Rochester, MN, United States
- 4National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
The AlkB family consists of Fe(II)- and α-ketoglutarate-dependent dioxygenases that can catalyze demethylation on a variety of substrates, such as RNA and DNA, subsequently affecting tumor progression and prognosis. However, their detailed functional roles in lung adenocarcinoma (LUAD) have not been clarified in a comprehensive manner. In this study, several bioinformatics databases, such as ONCOMINE, TIMER, and DiseaseMeth, were used to evaluate the expression profiles and prognostic significance of the AlkB family (ALKBH1-8 and FTO) in LUAD. The expression levels of ALKBH1/2/4/5/7/8 were significantly increased in LUAD tissues, while the expression levels of ALKBH3/6 and FTO were decreased. The main functions of differentially expressed AlkB homologs are related to the hematopoietic system and cell adhesion molecules. We also found that the expression profiles of the AlkB family are highly correlated with infiltrating immune cells (i.e., B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils and dendritic cells). In addition, DNA methylation analysis indicated that the global methylation levels of ALKBH1/2/4/5/6/8 and FTO were decreased, while the global methylation levels of ALKBH3/7 were increased. In addition, the patients with upregulated ALKBH2 have significantly poor overall survival (OS) and post-progressive survival (PPS). Taken together, our work could provide insightful information about aberrant AlkB family members as potential biomarkers for the diagnostic and prognostic evaluation of LUAD. Especially, ALKBH2 could be served as a therapeutic candidate for treating LUAD.
Introduction
Lung cancer is one of the most common malignant tumors in the world (Yan et al., 2019). Thousands of patients die from this malignant disease every year. In addition, the incidence of lung cancer continues to increase (Lin et al., 2019; Subedi et al., 2019; Arbyn et al., 2020). Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer (Wei et al., 2019). At present, great progress has been made in the diagnosis and treatment of LUAD (Lim and Ma, 2019; Zhang et al., 2019b), but there is still a lack of advanced diagnosis and treatment programs (Oberndorfer and Mullauer, 2018). Therefore, it is urgent to identify more therapeutic targets and prognostic biomarkers.
The AlkB family of Fe(II) and α-ketoglutarate-dependent dioxygenases is a universal class of direct reversal DNA repair enzymes (Pilzys et al., 2019; Xu et al., 2020). This family of enzymes can remove alkyl adducts from nucleobases by oxidative dealkylation, protecting the bacterial genome from alkylation damage (Trewick et al., 2002; Liu et al., 2018). Research has found that there are currently 9 homologs of the AlkB protein, including ALKBH1-8 and FTO (Du et al., 2019; Xiao et al., 2020). They have multiple biological functions, such as regulation of RNA metabolism, involvement in the DNA damage response or participation in fatty acid metabolism (Wu et al., 2016; Bian et al., 2019; Rajecka et al., 2019).
Many studies have found that the AlkB family plays a key role in the occurrence and development of tumors either directly or indirectly, including breast cancer, ovarian cancer, and bladder cancer (Fujii et al., 2013; Tao et al., 2016; Zhu et al., 2019). However, the detailed role of the AlkB family in LUAD remains to be further elucidated. With the rapid development of second-generation gene sequencing technology and the establishment of various databases, a comprehensive analysis of the AlkB family is beneficial to the clarification of the AlkB family in LUAD pathogenesis and treatment. In this study, we conducted a thorough and comprehensive bioinformatics analysis of the AlkB family (including ALKBH1-8 and FTO) in LUAD. Moreover, we evaluated their potential as therapeutic targets and prognostic biomarkers based on multiple public bioinformatics databases. The purpose of this study is to help clinicians choose appropriate therapeutic drugs and more accurately predict the long-term prognosis of patients with LUAD.
Materials and Methods
Oncomine3.0
Oncomine3.01 is a cancer microarray database and web-based data mining platform with 40 microarray data sets and approximately 100 differential expression analyses. It provides users with powerful, comprehensive genome-wide expression analysis (Rhodes et al., 2004, 2007). In this research, p < 0.05 and genes ranked in the top 10% were taken as the significance thresholds. Student’s t test was used to analyze the difference in AlkB family expression in LUAD. Specific information is summarized in Supplementary Table 1.
GEPIA2
GEPIA22, Gene Expression Profiling Interactive Analysis, is a web-based that provides key interactive and customizable features, including differential expression analysis, correlation analysis, and patient survival analysis (Tang et al., 2017, 2019). In this study, we used the “single gene analysis” module of GEPIA to analyze mRNA expression differences between tumors and normal tissues. At the same time, multiple gene comparison analysis of the AlkB family was performed using the “Multiple Gene Comparison” module and the “KIRC” GEPIA dataset. Student’s t test was used to generate p values for expression, and p < 0.05 was considered statistically significant.
UALCAN
UALCAN3 is an interactive network resource based on TCGA datasets that can compare the primary tumor and normal tissue samples based on pathological stage, tumor grade and other clinicopathological characteristics (Chandrashekar et al., 2017). In our research, the expression data of the AlkB family were obtained through the “stage analysis” module and the “KIRC” data set of UALCAN. Differences in transcriptional expression were compared by Student’s t test, and p < 0.05 was considered statistically significant.
The Human Protein Atlas
The Human Protein Atlas is a database of tools that can be used to identify clinically useful biomarkers using produced antibody and protein expression data (Asplund et al., 2012). Researchers can study the expression patterns of different proteins expressed in specific tumors. In this study, we directly compared the protein expression of Alkb family members in normal and LUAD tissues by immunohistochemistry.
Kaplan-Meier Plotter
Kaplan-Meier plotter4 was used to analyze the prognostic value of the AlkB family in LUAD (Sun et al., 2019; Wang et al., 2019). To analyze OS and PPS in patients with LUAD, the patient samples were divided into two groups by median expression (high and low expression) and were evaluated by the K-M survival chart. Information on the number of high-risk cases, median mRNA expression levels, HR, 95% CI, and p values can be found on the K-M plotter website. A p < 0.05 was considered statistically significant.
cBioPortal
c-BioPortal5, a comprehensive network of resources, can be used to explore, visualize and analyze multidimensional cancer genomics and clinical data. It contains more than 200 cancer genomics studies, including all the data on TCGA (Cerami et al., 2012; Wu et al., 2019). In this study, we analyzed the genome map of the AlkB family, which includes mutation and mRNA expression data. The mRNA expression z scores (RNA Seq V2 RSEM) were obtained using a z score threshold of ± 0.75.
GeneMANIA
Given a list of genes to query, GeneMANIA6 can use vast amounts of genomics and proteomics data to find genes with similar functions. In addition, GeneMANIA can predict gene function. Given a query gene, GeneMANIA finds genes that might share the target gene’s function based on their interactions (Zuberi et al., 2013; Franz et al., 2018).
Cytoscape
Cytoscape can integrate biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework (Shannon et al., 2003). In this study, we performed functional integration on 284 co-expressed molecules of AlkB family members screened from cBioPortal (the molecular names are provided in Supplementary Table 2). According to the degree values between the interacting proteins decided the nodes size. The higher the degree, the larger the circles.
WebGestalt
WebGestalt7 is a comprehensive, powerful, flexible, and interactive web-based analysis toolkit (Wang et al., 2017). Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed in this study using WebGestalt.
TIMER2.0
TIMER2.08 can assess immune cell infiltration and the clinical impact of 10,897 tumors from 32 cancer types (Li et al., 2017). In our study, the “Gene Module” was used to assess the association between the AlkB family and immune cell infiltration. The “survival module” was used to assess the correlation of clinical outcomes with immune cell infiltration and the AlkB family. The multivariable cox proportional hazard model is used as the statistical method.
DiseaseMeth2.0
DiseaseMeth2.09 aims to provide information on abnormal DNA methylation in human diseases, especially all sorts of cancer, in the most complete collection and comments (Lv et al., 2012; Xiong et al., 2017). We used Wanderer to screen the possible methylation value of the AlkB family. A p < 0.05 was considered as statically significant.
Results
Abnormal Expression of the AlkB Family in Patients With LUAD
We first searched the expression levels of the AlkB family (ALKBH1-ALKBH8, FTO) in LUAD and normal lung tissue using the Oncomine database. The results are shown in Table 1. The results demonstrated that the expression levels of ALKBH1/2/4/5/7/8 were significantly elevated in LUAD. Five datasets showed that ALKBH1 expression in LUAD was higher than that in normal lung tissue (Garber et al., 2001; Su et al., 2007; Landi et al., 2008; Hou et al., 2010; Okayama et al., 2012). In the Selamat and Hou datasets, the expression level of ALKBH2 in LUAD was significantly upregulated (Hou et al., 2010; Selamat et al., 2012). The results of Selamat (Selamat et al., 2012), Okayama (Okayama et al., 2012), Hou (Hou et al., 2010), and Su (Su et al., 2007) all suggest that ALKBH4 expression was significantly higher in LUAD compared with normal samples. The expression levels of ALKBH5 and ALKBH7 in LUAD were markedly higher than those in normal lung tissues in the Okayama (Okayama et al., 2012) and Hou (Hou et al., 2010) datasets. However, Selamat et al. (2012) found a decreased level of ALKBH8 in LUAD. Next, mRNA expression levels of the AlkB family in LUAD and normal lung tissue were verified in the GEPIA database. As shown in Figure 1A, the results showed clear support that the expression levels of ALKBH1/2/4/5/7/8 in LUAD were significantly elevated, while the expression levels of ALKBH3/6 and FTO were reduced compared with normal lung tissues. GEPIA was also used to compare the relative expression levels of the AlkB family in LUAD, and this analysis revealed that ALKBH5/7 had the highest relative expression levels among all AlkB family molecules (Figure 1B). After a comprehensive analysis of the mRNA expression levels of Alkb family in LUAD, we used the Human Protein Atlas to explore the protein expression levels of Alkb family in LUAD. The results showed that, ALKB1/2/4/5/7/8 was medium or high expression in LUAD. However, ALKB3/6 and FTO had low or no detect expression in LUAD. The results are shown in Figure 1C. This result is consistent with our previous findings on mRNA levels of expression. Moreover, UALCAN was used to analyze the relationship between the mRNA expression of AlkB family members and the clinicopathological staging of LUAD. From the results, it is clear that the mRNA expression level of ALKBH1/2/4/6 was positively correlated with tumor stage. In contrast, the mRNA expression levels of ALKBH7 and FTO were negatively correlated with tumor stage (Figure 2A). Similarly, the mRNA expression level of AlkB family members was also significantly correlated with lymph node metastasis. The highest mRNA expression levels of ALKBH2 and ALKBH3 were found in tumor N3 and were significantly correlated. At the same time, the mRNA expression levels of ALKBH2, ALKBH3, and ALKBH6 had a trend to upper expression in tumors with lymph node metastasis, although this was not significant. That may be caused by insufficient sample size and other reasons. Conversely, the expression levels of ALKBH5, ALKBH8 and FTO were lower than those in normal tissues at any stage of lymph node metastasis (Figure 2B).
Table 1. The expression profiles of the AlkB family in patients with LUAD using the Oncomine database.
Figure 1. (A) mRNA expression levels of different AlkB family members in LUAD and normal lung samples (GEPIA). T and N indicated the LUAD tissues and normal tissues, respectively. (B) The relative expression of the AlkB family in LUAD. (C) Immunohistochemical expression of ALKB family in LUAD tissue and normal lung tissue (the Human Protein Atlas).
Figure 2. (A) The relationship between mRNA expression and pathological stage of LUAD patients with different members of the AlkB family (UALCAN). (B) The relationship between mRNA expression of distinct AlkB family members with lymph node metastasis of LUAD patients (UALCAN). ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001.
Prognostic Value of mRNA Expression of the AlkB Family in LUAD Patients
Next, we used Kaplan-Meier plotter to analyze the prognostic values of the mRNA expression of AlkB family members in LUAD patients, including overall survival (OS) and post-progressive survival (PPS). The curves of OS are presented in Figure 3A. It is interesting to note that patients with higher transcription levels of ALKBH2, ALKBH4, and ALKBH7 displayed shorter OS times. The downregulations of ALKBH1, ALKBH3, ALKBH8, and FTO were significantly correlated with longer OS. However, ALKBH5 and ALKBH6 mRNA transcription showed no obvious correlation with the OS of LUAD patients. Moreover, the higher the transcription level of ALKBH2, the shorter the PPS in patients, while the higher the transcription levels of ALKBH3, ALKBH8 and FTO, the longer the PPS in patients. These results are summarized in Figure 3B.
Figure 3. (A) The overall survival curve of AlkB molecules in LUAD patients (Kaplan-Meier plotter). (B) The post-progressive survival curve of the AlkB molecules in LUAD patients (Kaplan-Meier plotter).
Genetic Alteration and Functional Analysis of the AlkB Family in LUAD Patients
We conducted a comprehensive biological function analysis to further grasp the molecular characteristics of the differentially expressed AlkB family. The genetic alterations of the differentially expressed AlkB family were evaluated in the temporary TCGA dataset. The results are shown in Figure 4A. ALKBH1, ALKBH2, ALKBH3, ALKBH4, ALKBH5, ALKBH6, ALKBH7, ALKBH8, and FTO were all altered, with 12, 6, 2, 2.8, 2.2, 3, 1, 2, and 1.8% alterations in the LUAD samples, respectively. We calculated the correlation by mRNA expression in the GEPIA online database, and Pearson’s correction was included. There was a strongly negative correlation between ALKBH2 and FTO, ALKBH7 and ALKBH8. ALKBH5 was positively correlated with FTO, and ALKBH1 was positively correlated with ALKBH8 (Figure 4B).
Figure 4. Genetic alternatives and correlation analysis of the AlkB family in LUAD. (A) Summary of the mutation rates in each AlkB family member in LUAD. (B) Correlation between different AlkB family in LUAD (GEPIA).
In fact, multiple of cellular proteins can work in a synergistic manner by forming multi-protein complexes. Moreover, protein-protein interactions (PPIs) play important roles among the different functions of a cell. Therefore, exploring PPI patterns is essential for understanding the structure and function of protein complexes (Planas-Iglesias et al., 2013; Furmanova et al., 2018). Afterward, to discover the PPI patterns between the differentially expressed AlkB family, we downloaded 248 co-expressed molecules with the highest correlation with AlkB family members using the cBioPortal database and modified them with Cytoscape (Figure 5A). These data suggest that TNF, ITGAM, ITGB2, ITGAX, DOCK2, IRF4, IRF8, CXCR5, CSF2RB, ZAP70, and TBX21 were primarily associated with the modulation and function of AlkB family in LUAD. In addition, GO analysis could be used to annotate the biological meaning of candidate biomarkers and their functions in various organisms. It is composed of molecular functions, biological processes and cellular components (Gene Ontology Consortium, 2019). And KEGG, a genome encyclopedia, can be used to perform biological interpretation of fully sequenced genomes and infer the systematic behaviors of targeted gene sets (Kanehisa et al., 2017). To probe the biological functions of these co-expressed molecules, WebGestalt was used to conduct GO annotation and KEGG pathway analysis for this study. GO annotation recognized that the AlkB family is mainly located on the cell membrane and is involved in biological regulation processes. In terms of molecular function, the family mainly participates in protein binding (Figure 5B). The corresponding pathways, as shown Figure 5C, were hematopoietic cell lineage, viral myocarditis and cell adhesion molecules (CAMs).
Figure 5. Coexpression network analysis of the AlkB family in LUAD. (A) The PPI network of AlkB family interaction partners generated by the frequently used c-BioPortal and Cytoscape. (B) Primary molecular functions, biological processes and cell components related to the biology of the AlkB family were identified using WebGestalt. (C) Bar plot of KEGG enriched terms analyzed by WebGestalt.
Immune Cell Infiltration of the AlkB Family in LUAD Patients
Members of the AlkB family are involved in the inflammatory response and the infiltration of immune cells in the tumor (Zhang et al., 2019a), thereby affecting the clinical outcomes of LUAD patients. The TIMER database was used to discuss the correlation between the differential expression of the AlkB family and immune cell infiltration. ALKBH1 and FTO were positively correlated with B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils and dendritic cells (Figures 6A,I), but ALKBH2 was negatively correlated with these immune cells (Figure 6B). There were positive correlations between ALKBH3 expression and the infiltration of B cells, CD4 + cells, neutrophils and dendritic cells (Figure 6C) and negative correlations between ALKBH4 expression and the infiltration of CD8 + T cells and macrophages (Figure 6D). Similarly, ALKBH5 was positively correlated with CD4 + T cells (Figure 6E). In contrast, ALKBH6 and ALKBH7 were negatively correlated with CD8 + T cells and macrophages (Figures 6F,G). We also found that ALKBH8 was positively correlated with CD4 + T cells, CD8 + T cells, macrophages, neutrophils and dendritic cells (Figure 6H). We also evaluated the correlation between the expression levels of the AlkB family and immune cell infiltration. After adjusting for some confounding factors, including B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils and dendritic cells, we used the Cox proportional hazard model to find that B cells were significantly associated with the clinical outcome of LUAD patients (Table 2).
Table 2. The Cox proportional hazard model of the ALKB family and six tumor-infiltrating immune cells in LUAD.
Figure 6. The correlations between differentially expressed AlkB family and immune cell infiltration. (A–I) The effect of ALKBH1-8 and FTO on the immune cell infiltration analyzed by TIMER2.0.
Methylation Expression Levels of the AlkB Family in LUAD Patients
DNA methylation is negatively correlated with gene expression (Gensous et al., 2020; Kronfol et al., 2020). And we want to evaluate the effect of methylation on the expression levels of Alkb family members. We downloaded the methylation expression data of AlkB family members from the DiseaseMeth database. The results reflect that the methylation expression levels of ALKBH1, ALKBH2, and ALKBH6 are lower in LUAD samples than in normal samples. Furthermore, the methylation expression levels of ALKBH3 were higher in LUAD samples (Supplementary Figure 1). Previous findings indicated the upregulated ALKBH1/2 and downregulated ALKBH3 in LUAD samples (Figure 1). Thus, the abnormal expression of ALKBH1/2/3 might be due to the changes in their methylation values.
Discussion
Emerging reports have revealed that DNA damage induces autophagy, apoptosis and senescence of cells, thus playing pivotal roles in the development of malignant tumors (Lord and Ashworth, 2012; Roos et al., 2016). As a classic reversal DNA repair enzyme, the AlkB family cannot be ignored in the development of multiple cancers (Xie et al., 2018; Zhu et al., 2019). At present, some scholars have discovered that ALKBH5 regulates epithelial-mesenchymal transition and angiogenesis to promote cancer progression (Panneerdoss et al., 2018). Chen et al. found that ALKBH3 can promote the proliferation, migration and invasion of cancer cells (Panneerdoss et al., 2018). Although a few members of the AlkB family have been discovered to play critical roles in lung cancer (Tasaki et al., 2011; Wu et al., 2011), the specific roles of other members of the AlkB family in LUAD have not yet been elucidated.
In this research, the AlkB family was studied comprehensively in terms of its expression and clinical prognostic value in LUAD. For the first time, we summarized the mRNA expression levels of nine members of the ALKBH family in LUAD compared with normal tissues. Next, a further novel finding is that the AlkB family is closely related to individual clinicopathological stages and lymph node metastasis. In addition, higher expression levels of ALKBH2/4/7 were significantly associated with shorter OS times. Lower expression levels of ALKBH1/3/8 and FTO were significantly associated with longer OS times. The patients with elevated ALKBH3/8 and FTO levels had long PPS, while patients with decreased ALKBH2 levels had shortened PPS. The process of tumor occurrence is intricate and multifaceted, and genetic alterations play a significant role in this process (Martincorena et al., 2017; Sondka et al., 2018; Coelho et al., 2019). The AlkB family shows evidence of frequent gene alterations in LUAD. The expression alteration of mRNA is one of the most typical mutations. These results comprehensively demonstrated that differential expression of the AlkB family probably plays a critical role in LUAD.
After discovering that the AlkB family is differentially expressed in LUAD and normal tissues, we next explored their molecular biological functions. We focused on the functional analysis of the AlkB family and their 247 frequently altered neighboring genes in LUAD. DNA repair- and protein methyltransferase activity-related genes, including ASCC3 and TRMT112, were highly linked to AlkB family alterations. Protein network interactions showed that TNF, ITGAM, ITGB2, ITGAX, DOCK2, IRF4, and IRF8 were primarily associated with the modulation and function of the differentially expressed AlkB family in LUAD. KEGG pathway analysis showed that AlkB family members were closely related to the CAMs. Multiple of CAM molecules have been proved to involve in the synaptic connection formation, maturation, function, and plasticity among cells. In a variety of human malignancies, including LUAD, tumor progression has been observed to be associated with altered CAM signaling pathways (Johnson, 1991). Similarly, Fujii and colleague reported the stimulative roles of ALKBH2 on the MUC1, an anti-adhesion molecule, on tumor cell growth (Fujii et al., 2013).
Increasing evidence shows that tumor progression and recurrence will be affected by immune cell infiltration and become an important determinant of immunotherapy response and clinical outcome (Lopes et al., 2019; Rotte, 2019). In our study, we found a significant correlation between the expression of the AlkB family and the infiltration of six immune cells, including B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils and dendritic cells. Previous studies have shown that the relationship between tumor cells and immune microenvironment plays a key role in the occurrence and development of tumors (Varn et al., 2017). Many researchers have found that the infiltration of immune cells is closely related to the OS or disease-specific survival rate (DSS) of many cancers, including LUAD (Kadara et al., 2017; Hofman, 2020; Luo et al., 2020). As a kind of immune cells, B cells can fight against tumors through the production of antibodies, antigen royalty or the secretion of immune regulatory factors (Stoycheva et al., 2021). Han et al. (2020) found that B cells can be used as an independent indicator of OS prediction in LUAD. It is well known that methylation plays an important role in the development of cancer. In this study, we have found that Alkb family has abnormal expression of DNA methylation in LUAD. In fact, many researchers have also made some breakthroughs in the modification of RNA methylation in Alkb family. N6-methyladenosine is the most common RNA modification. ALKBH5, A member of the Alkb family, is the most widespread m6A demethylase. It has been shown that ALKBH5 can inhibit tumor growth and metastasis (Ji et al., 2020; Jin et al., 2020). Zhu et al. (2020) found that m6A-related gene ALKBH5 had a good predictive value in pathological staging and prognosis of LUAD. The results of our study reflect that the AlkB family not only can be used as prognostic indicators but also reflect immune status.
In this study, we comprehensively analyzed the expression and prognosis of the AlkB family in LUAD from the perspective of bioinformatics. Through a comprehensive analysis, we found that high ALKBH2 expression was significantly related to poor OS and PPS in LUAD patients. In addition, abnormally expressed ALKBH2 could also alter the immune cell infiltration in LUAD. These results provide novel biomarkers for the diagnosis of LUAD in the future and will help clinicians design potential therapeutic strategies for LUAD patients to improve the survival rate.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
GW, JL, and ZX: acquisition of the data. YY: analysis and interpretation of the data. JZ: conception and design. YC and BP: data curation. JH: development of the methodology. GW and ZX: writing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Natural Science Foundation of China (No. 81803035), the Natural Science Foundation of Hunan Province (2020JJ5934 and 2019JJ50932), and the China Postdoctoral Science Foundation (2020M672521).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.633927/full#supplementary-material
Supplementary Figure 1 | The methylation expression of AlkB family members in LUAD patients (DiseaseMeth).
Footnotes
- ^ http://www.oncomine.org
- ^ http://gepia.cancer-pku.cn/index.html
- ^ http://ualcan.path.uab.edu/index.html
- ^ http://kmplot.com/analysis/
- ^ http://cbioportal.org
- ^ http://genemania.org
- ^ http://www.webgestalt.org
- ^ https://cistrome.shinyapps.io/timer/
- ^ http://bioinfo.hrbmu.edu.cn/diseasemeth/
References
Arbyn, M., Weiderpass, E., Bruni, L., de Sanjose, S., Saraiya, M., Ferlay, J., et al. (2020). Estimates of incidence and mortality of cervical cancer in 2018: a worldwide analysis. Lancet Glob. Health 8, e191–e203. doi: 10.1016/S2214-109X(19)30482-6
Asplund, A., Edqvist, P. H., Schwenk, J. M., and Ponten, F. (2012). Antibodies for profiling the human proteome-The Human Protein Atlas as a resource for cancer research. Proteomics 12, 2067–2077. doi: 10.1002/pmic.201100504
Bian, K., Lenz, S. A. P., Tang, Q., Chen, F., Qi, R., Jost, M., 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. doi: 10.1093/nar/gkz395
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2, 401–404. doi: 10.1158/2159-8290.CD-12-0095
Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B., et al. (2017). UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 19, 649–658. doi: 10.1016/j.neo.2017.05.002
Coelho, M. C., Pinto, R. M., and Murray, A. W. (2019). Heterozygous mutations cause genetic instability in a yeast model of cancer evolution. Nature 566, 275–278. doi: 10.1038/s41586-019-0887-y
Du, K., Zhang, L., Lee, T., and Sun, T. (2019). m(6)A RNA methylation controls neural development and is involved in human diseases. Mol. Neurobiol. 56, 1596–1606. doi: 10.1007/s12035-018-1138-1
Franz, M., Rodriguez, H., Lopes, C., Zuberi, K., Montojo, J., Bader, G. D., et al. (2018). GeneMANIA update 2018. Nucleic Acids Res. 46, W60–W64. doi: 10.1093/nar/gky311
Fujii, T., Shimada, K., Anai, S., Fujimoto, K., and Konishi, N. (2013). ALKBH2, a novel AlkB homologue, contributes to human bladder cancer progression by regulating MUC1 expression. Cancer Sci. 104, 321–327. doi: 10.1111/cas.12089
Furmanova, K., Byska, J., Groller, E. M., Viola, I., Palecek, J. J., and Kozlikova, B. (2018). COZOID: contact zone identifier for visual analysis of protein-protein interactions. BMC Bioinformatics 19:125. doi: 10.1186/s12859-018-2113-6
Garber, M. E., Troyanskaya, O. G., Schluens, K., Petersen, S., Thaesler, Z., Pacyna-Gengelbach, M., et al. (2001). Diversity of gene expression in adenocarcinoma of the lung. Proc. Natl. Acad. Sci. U.S.A. 98, 13784–13789. doi: 10.1073/pnas.241500798
Gene Ontology Consortium (2019). The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 47, D330–D338. doi: 10.1093/nar/gky1055
Gensous, N., Garagnani, P., Santoro, A., Giuliani, C., Ostan, R., Fabbri, C., et al. (2020). One-year Mediterranean diet promotes epigenetic rejuvenation with country- and sex-specific effects: a pilot study from the NU-AGE project. Geroscience 42, 687–701. doi: 10.1007/s11357-019-00149-0
Han, L., Shi, H., Luo, Y., Sun, W., Li, S., Zhang, N., et al. (2020). Gene signature based on B cell predicts clinical outcome of radiotherapy and immunotherapy for patients with lung adenocarcinoma. Cancer Med. 9, 9581–9594. doi: 10.1002/cam4.3561
Hofman, P. (2020). New insights into the interaction of the immune system with non-small cell lung carcinomas. Transl. Lung Cancer Res. 9, 2199–2213. doi: 10.21037/tlcr-20-178
Hou, J., Aerts, J., den Hamer, B., van Ijcken, W., den Bakker, M., Riegman, P., et al. (2010). Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One 5:e10312. doi: 10.1371/journal.pone.0010312
Ji, G., Huang, C., He, S., Gong, Y., Song, G., Li, X., et al. (2020). Comprehensive analysis of m6A regulators prognostic value in prostate cancer. Aging 12, 14863–14884. doi: 10.18632/aging.103549
Jin, D., Guo, J., Wu, Y., Yang, L., Wang, X., Du, J., et al. (2020). m(6)A demethylase ALKBH5 inhibits tumor growth and metastasis by reducing YTHDFs-mediated YAP expression and inhibiting miR-107/LATS2-mediated YAP activity in NSCLC. Mol. Cancer 19:40. doi: 10.1186/s12943-020-01161-1
Johnson, J. P. (1991). Cell adhesion molecules of the immunoglobulin supergene family and their role in malignant transformation and progression to metastatic disease. Cancer Metastasis Rev. 10, 11–22. doi: 10.1007/BF00046840
Kadara, H., Choi, M., Zhang, J., Parra, E. R., Rodriguez-Canales, J., Gaffney, S. G., et al. (2017). Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma with fully annotated clinical follow-up. Ann. Oncol. 28, 75–82. doi: 10.1093/annonc/mdw436
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 45, D353–D361. doi: 10.1093/nar/gkw1092
Kronfol, M. M., Jahr, F. M., Dozmorov, M. G., Phansalkar, P. S., Xie, L. Y., Aberg, K. A., et al. (2020). DNA methylation and histone acetylation changes to cytochrome P450 2E1 regulation in normal aging and impact on rates of drug metabolism in the liver. Geroscience 42, 819–832. doi: 10.1007/s11357-020-00181-5
Landi, M. T., Dracheva, T., Rotunno, M., Figueroa, J. D., Liu, H., Dasgupta, A., et al. (2008). Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS One 3:e1651. doi: 10.1371/journal.pone.0001651
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77, e108–e110. doi: 10.1158/0008-5472.CAN-17-0307
Lim, Z. F., and Ma, P. C. (2019). Emerging insights of tumor heterogeneity and drug resistance mechanisms in lung cancer targeted therapy. J. Hematol. Oncol. 12:134. doi: 10.1186/s13045-019-0818-2
Lin, L., Yan, L., Liu, Y., Yuan, F., Li, H., and Ni, J. (2019). Incidence and death in 29 cancer groups in 2017 and trend analysis from 1990 to 2017 from the Global Burden of Disease Study. J. Hematol. Oncol. 12:96. doi: 10.1186/s13045-019-0783-9
Liu, Y., Yuan, Q., and Xie, L. (2018). The AlkB family of Fe (II)/alpha-ketoglutarate-dependent dioxygenases modulates embryogenesis through epigenetic regulation. Curr. Stem Cell Res. Ther. 13, 136–143. doi: 10.2174/1574888X12666171027105532
Lopes, A., Vandermeulen, G., and Preat, V. (2019). Cancer DNA vaccines: current preclinical and clinical developments and future perspectives. J. Exp. Clin. Cancer Res. 38:146. doi: 10.1186/s13046-019-1154-7
Lord, C. J., and Ashworth, A. (2012). The DNA damage response and cancer therapy. Nature 481, 287–294. doi: 10.1038/nature10760
Luo, C., Lei, M., Zhang, Y., Zhang, Q., Li, L., Lian, J., et al. (2020). Systematic construction and validation of an immune prognostic model for lung adenocarcinoma. J. Cell. Mol. Med. 24, 1233–1244. doi: 10.1111/jcmm.14719
Lv, J., Liu, H., Su, J., Wu, X., Liu, H., Li, B., et al. (2012). DiseaseMeth: a human disease methylation database. Nucleic Acids Res. 40, D1030–D1035. doi: 10.1093/nar/gkr1169
Martincorena, I., Raine, K. M., Gerstung, M., Dawson, K. J., Haase, K., Van Loo, P., et al. (2017). Universal patterns of selection in cancer and somatic tissues. Cell 171, 1029–1041.e21. doi: 10.1016/j.cell.2017.09.042
Oberndorfer, F., and Mullauer, L. (2018). Molecular pathology of lung cancer: current status and perspectives. Curr. Opin. Oncol. 30, 69–76. doi: 10.1097/CCO.0000000000000429
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2012). Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 72, 100–111. doi: 10.1158/0008-5472.CAN-11-1403
Panneerdoss, S., Eedunuri, V. K., Yadav, P., Timilsina, S., Rajamanickam, S., Viswanadhapalli, S., et al. (2018). Cross-talk among writers, readers, and erasers of m(6)A regulates cancer growth and progression. Sci. Adv. 4:eaar8263. doi: 10.1126/sciadv.aar8263
Pilzys, T., Marcinkowski, M., Kukwa, W., Garbicz, D., Dylewska, M., Ferenc, K., et al. (2019). ALKBH overexpression in head and neck cancer: potential target for novel anticancer therapy. Sci. Rep. 9:13249. doi: 10.1038/s41598-019-49550-x
Planas-Iglesias, J., Bonet, J., Garcia-Garcia, J., Marin-Lopez, M. A., Feliu, E., and Oliva, B. (2013). Understanding protein-protein interactions using local structural features. J. Mol. Biol. 425, 1210–1224. doi: 10.1016/j.jmb.2013.01.014
Rajecka, V., Skalicky, T., and Vanacova, S. (2019). The role of RNA adenosine demethylases in the control of gene expression. Biochim. Biophys. Acta Gene Regul. Mech. 1862, 343–355. doi: 10.1016/j.bbagrm.2018.12.001
Rhodes, D. R., Kalyana-Sundaram, S., Mahavisno, V., Varambally, R., Yu, J., Briggs, B. B., et al. (2007). Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia 9, 166–180. doi: 10.1593/neo.07112
Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia 6, 1–6. doi: 10.1016/s1476-5586(04)80047-2
Roos, W. P., Thomas, A. D., and Kaina, B. (2016). DNA damage and the balance between survival and death in cancer biology. Nat. Rev. Cancer 16, 20–33. doi: 10.1038/nrc.2015.2
Rotte, A. (2019). Combination of CTLA-4 and PD-1 blockers for treatment of cancer. J. Exp. Clin. Cancer Res. 38:255. doi: 10.1186/s13046-019-1259-z
Selamat, S. A., Chung, B. S., Girard, L., Zhang, W., Zhang, Y., Campan, M., et al. (2012). Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 22, 1197–1211. doi: 10.1101/gr.132662.111
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Sondka, Z., Bamford, S., Cole, C. G., Ward, S. A., Dunham, I., and Forbes, S. A. (2018). The COSMIC cancer gene census: describing genetic dysfunction across all human cancers. Nat. Rev. Cancer 18, 696–705. doi: 10.1038/s41568-018-0060-1
Stoycheva, D., Simsek, H., Weber, W., Hauser, A. E., and Klotzsch, E. (2021). External cues to drive B cell function towards immunotherapy. Acta Biomater. 2021:S1742-7061(21)00121-5. doi: 10.1016/j.actbio.2021.02.026
Su, L. J., Chang, C. W., Wu, Y. C., Chen, K. C., Lin, C. J., Liang, S. C., et al. (2007). Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC Genomics 8:140. doi: 10.1186/1471-2164-8-140
Subedi, P., Nembrini, S., An, Q., Zhu, Y., Peng, H., Yeh, F., et al. (2019). Telomere length and cancer mortality in American Indians: the Strong Heart Study. Geroscience 41, 351–361. doi: 10.1007/s11357-019-00080-4
Sun, C. C., Li, S. J., Hu, W., Zhang, J., Zhou, Q., Liu, C., et al. (2019). Comprehensive analysis of the expression and prognosis for E2Fs in human breast cancer. Mol. Ther. 27, 1153–1165. doi: 10.1016/j.ymthe.2019.03.019
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 47, W556–W560. doi: 10.1093/nar/gkz430
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45, W98–W102. doi: 10.1093/nar/gkx247
Tao, Y., Hu, K., Tan, F., Zhang, S., Zhou, M., Luo, J., et al. (2016). SH3-domain binding protein 1 in the tumor microenvironment promotes hepatocellular carcinoma metastasis through WAVE2 pathway. Oncotarget 7, 18356–18370. doi: 10.18632/oncotarget.7786
Tasaki, M., Shimada, K., Kimura, H., Tsujikawa, K., and Konishi, N. (2011). ALKBH3, a human AlkB homologue, contributes to cell survival in human non-small-cell lung cancer. Br. J. Cancer 104, 700–706. doi: 10.1038/sj.bjc.6606012
Trewick, S. C., Henshaw, T. F., Hausinger, R. P., Lindahl, T., and Sedgwick, B. (2002). Oxidative demethylation by Escherichia coli AlkB directly reverts DNA base damage. Nature 419, 174–178. doi: 10.1038/nature00908
Varn, F. S., Wang, Y., Mullins, D. W., Fiering, S., and Cheng, C. (2017). Systematic pan-cancer analysis reveals immune cell interactions in the tumor microenvironment. Cancer Res. 77, 1271–1282. doi: 10.1158/0008-5472.CAN-16-2490
Wang, B., Ran, Z., Liu, M., and Ou, Y. (2019). Prognostic significance of potential immune checkpoint member HHLA2 in human tumors: a comprehensive analysis. Front. Immunol. 10:1573. doi: 10.3389/fimmu.2019.01573
Wang, J., Vasaikar, S., Shi, Z., Greer, M., and Zhang, B. (2017). WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 45, W130–W137. doi: 10.1093/nar/gkx356
Wei, C., Dong, X., Lu, H., Tong, F., Chen, L., Zhang, R., et al. (2019). LPCAT1 promotes brain metastasis of lung adenocarcinoma by up-regulating PI3K/AKT/MYC pathway. J. Exp. Clin. Cancer Res. 38:95. doi: 10.1186/s13046-019-1092-4
Wu, P., Heins, Z. J., Muller, J. T., Katsnelson, L., de Bruijn, I., Abeshouse, A. A., et al. (2019). Integration and analysis of CPTAC proteomics data in the context of cancer genomics in the cBioPortal. Mol. Cell. Proteomics 18, 1893–1898. doi: 10.1074/mcp.TIR119.001673
Wu, S. S., Xu, W., Liu, S., Chen, B., Wang, X. L., Wang, Y., et al. (2011). Down-regulation of ALKBH2 increases cisplatin sensitivity in H1299 lung cancer cells. Acta Pharmacol. Sin. 32, 393–398. doi: 10.1038/aps.2010.216
Wu, T. P., Wang, T., Seetin, M. G., Lai, Y., Zhu, S., Lin, K., et al. (2016). DNA methylation on N(6)-adenine in mammalian embryonic stem cells. Nature 532, 329–333. doi: 10.1038/nature17640
Xiao, M. Z., Liu, J. M., Xian, C. L., Chen, K. Y., Liu, Z. Q., and Cheng, Y. Y. (2020). Therapeutic potential of ALKB homologs for cardiovascular disease. Biomed. Pharmacother. 131:110645. doi: 10.1016/j.biopha.2020.110645
Xie, Q., Wu, T. P., Gimple, R. C., Li, Z., Prager, B. C., Wu, Q., et al. (2018). N(6)-methyladenine DNA modification in glioblastoma. Cell 175, 1228–1243.e20. doi: 10.1016/j.cell.2018.10.006
Xiong, Y., Wei, Y., Gu, Y., Zhang, S., Lyu, J., Zhang, B., et al. (2017). DiseaseMeth version 2.0: a major expansion and update of the human disease methylation database. Nucleic Acids Res. 45, D888–D895. doi: 10.1093/nar/gkw1123
Xu, Z., Peng, B., Cai, Y., Wu, G., Huang, J., Gao, M., et al. (2020). N6-methyladenosine RNA modification in cancer therapeutic resistance: current status and perspectives. Biochem. Pharmacol. 182:114258. doi: 10.1016/j.bcp.2020.114258
Yan, Y., Chen, X., Wang, X., Zhao, Z., Hu, W., Zeng, S., et al. (2019). The effects and the mechanisms of autophagy on the cancer-associated fibroblasts in cancer. J. Exp. Clin. Cancer Res. 38:171. doi: 10.1186/s13046-019-1172-5
Zhang, C., Fu, J., and Zhou, Y. A. (2019a). Review in research progress concerning m6A methylation and immunoregulation. Front. Immunol. 10:922. doi: 10.3389/fimmu.2019.00922
Zhang, C., Leighl, N. B., Wu, Y. L., and Zhong, W. Z. (2019b). Emerging therapies for non-small cell lung cancer. J. Hematol. Oncol. 12:45. doi: 10.1186/s13045-019-0731-8
Zhu, H., Gan, X., Jiang, X., Diao, S., Wu, H., and Hu, J. (2019). ALKBH5 inhibited autophagy of epithelial ovarian cancer through miR-7 and BCL-2. J. Exp. Clin. Cancer Res. 38:163. doi: 10.1186/s13046-019-1159-2
Zhu, J., Wang, M., and Hu, D. (2020). Deciphering N(6)-methyladenosine-related genes signature to predict survival in lung adenocarcinoma. Biomed Res. Int. 2020:2514230. doi: 10.1155/2020/2514230
Keywords: AlkB family, lung adenocarcinoma, expression profiles, prognosis, methylation, immune cell infiltration
Citation: Wu G, Yan Y, Cai Y, Peng B, Li J, Huang J, Xu Z and Zhou J (2021) ALKBH1-8 and FTO: Potential Therapeutic Targets and Prognostic Biomarkers in Lung Adenocarcinoma Pathogenesis. Front. Cell Dev. Biol. 9:633927. doi: 10.3389/fcell.2021.633927
Received: 26 November 2020; Accepted: 03 May 2021;
Published: 03 June 2021.
Edited by:
Rodrigo Nalio Ramos, INSERM U1138 Centre de Recherche des Cordeliers (CRC), FranceReviewed by:
Jimena Tosello, Institut Curie, FranceKewa Gao, University of California, Davis, United States
Copyright © 2021 Wu, Yan, Cai, Peng, Li, Huang, Xu and Zhou. 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: Zhijie Xu, xzj1322007@csu.edu.cn; Jianhua Zhou, zhoujh15@163.com
†These authors have contributed equally to this work