Liexi Xu1†
Zihang Zeng
Jiali Li
Conghua Xie- 1Department of Radiation and Medical Oncology, Wuhan University of Zhongnan Hospital, Wuhan, China
- 2Tumor Precision Diagnosis and Treatment Technology and Translational Medicine, Hubei Engineering Research Center, Zhongnan Hospital of Wuhan University, Wuhan, China
- 3Hubei Key Laboratory of Tumor Biological Behaviors, Zhongnan Hospital of Wuhan University, Wuhan, China
Background: Abnormal DNA methylation of gene promoters is an important feature in lung adenocarcinoma (LUAD). However, the prognostic value of DNA methylation remains to be further explored. Objectives. We sought to explore DNA methylation characteristics and develop a quantifiable criterion related to DNA methylation to improve survival prediction for LUAD patients.
Methods: Illumina Human Methylation450K array data, level 3 RNA-seq data and corresponding clinical information were obtained from TCGA. Cox regression analysis and the Akaike information criterion were used to construct the best-prognosis methylation signature. Receiver operating characteristic curve analysis was used to validate the prognostic ability of the DNA methylation-related feature score. qPCR was used to measure the transcription levels of the identified genes upon methylation.
Results: We identified a set of DNA methylation features composed of 11 genes (MYEOV, KCNU1, SLC27A6, NEUROD4, HMGB4, TACR3, GABRA5, TRPM8, NLRP13, EDN3 and SLC34A1). The feature score, calculated based on DNA methylation features, was independent of tumor recurrence and TNM stage in predicting overall survival. Of note, the combination of this feature score and TNM stage provided a better overall survival prediction than either of them individually. The transcription levels of all the hypermethylated genes were significantly increased after demethylation, and the expression levels of 3 hypomethylated proteins were significantly higher in tumor tissues than in normal tissues, as indicated by immunohistochemistry data from the Human Protein Atlas. Our results suggested that these identified genes with prognostic features were regulated by DNA methylation of their promoters.
Conclusion: Our studies demonstrated the potential application of DNA methylation markers in the prognosis of LUAD.
1 Introduction
Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer, accounting for approximately 50% of all lung cancer cases in most countries (Goldstraw et al., 2011; Bray et al., 2018). Previous studies have revealed that in addition to cigarette smoking, risk factors such as age, environmental pollution, occupational exposure, race, sex, and preexisting lung disease are also substantially involved in lung cancer. With the development and popularization of public databases in recent years, an increasing number of researchers have tried to identify prognostic biomarkers for LUAD by analyzing clinical characteristics and molecular information (Edge and Compton, 2010). The TNM staging system of the American Joint Commission on Cancer (AJCC) was reported to have great value in LUAD prognosis (Folkman, 1971). Liu et al. stated that conventional staging alone was not enough to predict prognosis and guide treatment decisions. They analyzed large cohorts from The Cancer Genome Atlas (TCGA) database and developed a 4-gene feature related to glycolysis (Kaishang et al., 2018). Su et al. identified an RNA sequencing network of 29 key lncRNAs, 72 mRNAs and 24 miRNAs as potential biomarkers to optimize the diagnosis and prognosis of LUAD patients by using the TCGA database (Li et al., 2017). The findings from these studies indicate that it is feasible to use different molecular markers and clinical features in public databases to establish practical models that have great application potential. Although the effectiveness of these prediction models has not been tested in clinical practice, it is necessary to continue to mine and improve the gene signatures related to the prognosis of LUAD.
Epigenetic disorders, especially abnormal DNA methylation in gene promoters, are a fundamental feature of human malignant tumors (Liu et al., 2019). As one of the most well-studied epigenetic modifications, DNA methylation mainly occurs at 5′-cytosine-phosphate-guanine-3' (CpG) dinucleotides and is regulated by DNA methyltransferases and DNA demethylases (Yang et al., 2022). Methylation and cancer formation are associated in 2 main ways: one is the regulation of tumor suppressor gene expression by gene hypermethylation in the promoter, and the other is genome-wide hypomethylation, which plays an important role in the stability of the heterochromatin structure (Feinberg, 2007). In virtually every step of tumor progression, there is abnormal promoter methylation regulation (Sui et al., 2016). APC, CDH13, MLH1 and IRX1 have hypermethylation in promoter CpG islands (CGIs). The hypermethylation of the APC and CDH13 genes in LUAD is associated with cancer cell adhesion, and the loss of MLH1 and IRX1 expression is associated with poor tumor survival (Goto et al., 2009; Küster et al., 2020). The hypomethylation of LINE-1 and ELF3 induces protein overexpression in LUAD. The overexpression of ELF3 can stimulate the carcinogenic phenotype of LUAD cells and reduce the survival time of patients, suggesting that the hypomethylation of LINE-1 is a prognostic marker of LUAD development and progression (Ikeda et al., 2013; Enfield et al., 2019).
TCGA has disclosed the clinical information of more than 10,000 patients and the molecular phenotype information of their tumor tissues. This information covers 33 different types of tumors and multiple data from different sources, including transcriptomic, methylomic and proteomic sources (Mazor et al., 2016). By integrating data from different sources, we can identify specific events in the carcinogenic process and identify potential biomarkers associated with patient survival.
In this study, we obtained the TCGA Illumina Human Methylation 450K microarray data, RNA-seq data, and clinical data of LUAD patients and performed an integrative analysis to identify a set of DNA methylation features for 11 genes. We performed area under the receiver operating characteristic curve (AUC-ROC) analysis to verify the ability of the identified DNA methylation feature to predict the survival of LUAD. In addition, we performed qPCR, and the results suggested that these identified genes with prognostic features were regulated by DNA methylation in their promoters.
2 Materials and methods
2.1 Data preparation
The steps of data acquisition and analysis, as well as methylation feature acquisition and verification, are shown in the flow chart (Figure 1). Illumina Human Methylation 450K array data were obtained from TCGA, and a total of 24,587 DMSs from 478 pretreated methylation arrays were screened using the camp and Minfi R software packages. After that, 478 samples (449 LUAD samples and 29 normal samples) were included after being filtered, inspected, and standardized with the ChAMP R package. Level 3 RNA-seq data from TCGA were normalized and log2 transformed by the edgeR package. For the preprocessing of clinical information corresponding to the sample, patients with nonsurvival status or survival time less than 1 month were excluded because of other disease-related deaths.
FIGURE 1. Flow chart for obtaining and verifying methylation feature.
2.2 Differential methylation analysis and differential expression analysis
In total, 449 LUAD samples and 29 normal samples were subjected to differential methylation analysis with the ChAMP R package and the Minfi R package (Li et al., 2019). Principal component analysis was used to detect the sample quality. The ChAMP DMP function and the Minfi R package defined the methylation loci with an average methylation difference >0.2 and a false discovery rate <0.05 as differentially methylated sites (DMSs), and the final DMSs were obtained through the intersection of the two (Li et al., 2018). Differentially expressed genes (DEGs) between the 526 LUAD samples and the 59 normal samples were analyzed with the Limma R package (p < 0.05 and |log2FC| ≥ 1). Metascape (https://metascape.org/gp/index.html#/main/step1) was used to analyze the pathway enrichment of the hyper-down and hypo-up methylated related differential expression genes (mrDEGs) groups.
2.3 Survival model construction process
A prognosis prediction model was established according to the DNA methylation β value of mrDEGs and matched prognostic data of patients. According to the methylation β value, univariate Cox regression analysis was used to screen mrDEGs (p < 0.01) that were significantly associated with overall survival (OS). Then, mrDEGs identified in the univariate Cox regression analysis were subjected to multivariate Cox regression analysis (Lian et al., 2019). At the same time, the Akaike information criterion (AIC) was used to screen out the genes with subtly individual but significantly synergistic effects to determine the most appropriate gene feature (Tozzi et al., 2020). A Kaplan‒Meier (K-M) curve with a log rank test was used to validate the survival difference of patients (Zhao et al., 2020). Harrell’s concordance index (C-index) and the corresponding 95% confidence intervals (CIs) were calculated to determine the prognostic model’s ability. These steps were performed by R with the survival and survcomp R packages (Schröder et al., 2011).
2.4 Consensus clustering analysis
We selected 155 LUAD samples with complete clinical information. For our standard, we considered the standard deviation of the β value in tumor samples to be greater than 0.2 and the average β value in normal tissues to be less than 0.05. With this approach, we selected 641 methylation probes according to the standard. According to the PAM algorithm and Euclidean distance, we then performed unsupervised consistent clustering on 641 probes of 155 samples. The Consensus Cluster Plus R package was used for the clustering analysis (Wilkerson and Hayes, 2010). The Kruskal‒Wallis test was used to validate the significance of clinical features among clusters.
2.5 Validation experiments in cell lines
qPCR was used to verify the changes in gene transcription levels upon methylation. A549, PC9 and H1975 cells were purchased from the Canadian Standards Association (CSA). All cell lines were cultured in RPMI 1640 with 10% fetal bovine serum. All experimental cells were treated with 5-Aza-2′- deoxycytidine (5-aza, Aladdin) for 96 h at 1 µM. qPCR analyses of all cell lines were repeated at least 3 times (Christman, 2002). All primer sequences used in qPCR are listed in Attachment 1: Table 1.
TABLE 1. All primer sequences used in qPCR.
2.6 Statistical analysis
The correlation between feature scores and clinical factors was analyzed by the chi square or Fisher exact test (Jung, 2014; Pandis, 2016). Multivariate Cox regression combined with hierarchical data analysis was used to evaluate the predictive power of the clinical features, TNM stage, and methylation feature score for prognosis. The forest map was drawn by Prism, and other statistical tests were performed by R using the corresponding R packages.
3 Results
3.1 Differential methylation and the identification of mrDMGs
Illumina Human Methylation 450K array data were obtained from TCGA. We screened a total of 24,587 DMSs from data from 478 pretreated methylation arrays using the camp and Minfi R software packages. We then divided DMSs into 15,387 hypermethylated and 9,200 hypomethylated sites and evaluated their distribution in the genome. Compared to 63% in the whole genome, hypermethylated sites increased significantly in the promoter, CGIs, and CGI promoters (71%, 96% and 99%, respectively, Figure 2A). At the same time, most DMSs on CGIs were hypermethylated (96%), and most DMSs on shelf CpG positions were hypomethylated (Figure 2B). When we detected the distribution of DMSs around the gene, we found that the hypermethylation of CpGs was higher near the transcription start site (TSS). For example, the proportions of hypermethylated CpGs in the 5′UTR, tss200 and first exon were 70%, 77%, and 76%, respectively (Figure 2C). Next, we located DMSs on the gene and obtained 5,900 differentially methylated genes (DMGs) (Figure 2D). Next, we identified 1,887 DEGs based on the RNA-seq data. Then, 406 mrDEGs were determined through the intersection of the DMGs and the identified DEGs (Figure 2E). Among them, 43 mrDEGs were in the hypermethylation downregulation group (hyper-down group), and 126 mrDEGs were in the hypomethylation upregulation group (hypo-up group) (Figure 2F).
FIGURE 2. Distribution of DMSs and obtain mrDEGs of LUAD. (A) Distribution of DMSs across various genomic regions, including CpG islands (CGI), promoters, CGI promoter, and the whole genome (all). (B) Distribution of DMSs in various areas related to CGI distance, including CpG shelves, CpG shores and CpG islands. (C) Distribution of DMSs in gene location, including 3′ UTRs, gene bodies, first exons, 5′ UTRs, TSS200 and TSS1500. (D) Distribution of DMGs across various genomic regions. (E) Scatter plot shows mean methylation difference versus log2 expression change, and each point represents a pair of methylation site and gene. (F) Venn diagrams shows the intersection between DEGs and hypermethylated genes (left) and between DEGs and hypomethylated genes (right).
3.2 mrDEGs involved in biological processes
The Metascape website was used to analyze the pathway enrichment of the hyper-down and hypo-up mrDEGs. In the hypo-up group, the genes showed a significant abundance in fatty acid degradation, cyclic adenosine monophosphate (cAMP)-mediated signaling, glycolysis/gluconeogenesis, etc. (Figure 3A). Cancer is usually accompanied by nutritional metabolic imbalances, such as abnormal glucose and lipid metabolism (Li and Liao, 2021). cAMP was the first second messenger to be discovered, and it plays key roles in physiological defects caused by metabolic disorders (Zhang et al., 2020; Chi et al., 2021). Interestingly, in the hyper-down group, there was also gene enrichment related to fatty acid degradation (Figure 3B). The effects of lipid metabolism disorders on cancer have attracted increasing attention in recent years (Karagiota et al., 2022). There is no doubt that to reprogram their metabolic state and ensure cell survival, tumor cells need epigenetic modifications to regulate gene expression. For example, brother of the regulator of imprinted sites can mediate the Warburg effect and promote breast cancer by regulating the methylation of pyruvate kinase M1/2 (PKM) exons (Singh et al., 2017; Huo et al., 2021). Our results suggest that abnormal lipid metabolism in LUAD may be closely mediated by DNA methylation (Figure 3B).
FIGURE 3. Pathway enrichment analysis of mrDEGs in LUAD. (A) The pathway enrichment analysis of the upregulated mrDEGs. Each node represents a gene group. The node size is proportional to the total number of genes in each gene set. The width of the line between nodes represents the proportion of genes shared among gene sets. (B) The pathway enrichment analysis of the downregulated mrDEGs.
3.3 Identification of prognostic mrDEGs
We established a prognosis prediction model according to the methylation data of mrDEGs and the matched prognostic data. We first analyzed 406 mrDEGs in 453 patients by univariate Cox regression and identified 32 mrDEGs related to prognosis (p < 0.05). Next, multivariate Cox regression analysis was performed on the 32 mrDEGs. The AIC, as the indicator for model fitness, determined the most suitable prognostic model. Finally, we identified 11 mrDEGs (MYEOV, KCNU1, SLC27A6, NEUROD4, HMGB4, TACR3, GABRA5, TRPM8, NLRP13, EDN3 and SLC34A1) to be included in a DNA methylation feature prognostic model. There were 5 genes (MYEOV, KCNU1, SLC27A6, NEUROD4 and HMGB4) with statistically nonsignificant p values in multivariate Cox regression analysis (Table 2). However, the AIC of this prognostic model was the lowest (AIC = 1733.8, p = 2e-05), indicating that this model was the most suitable, and the overall effect of the model was significant. The C-index of the identified DNA methylation feature model was 0.666 (95% CI = 0.641–0.690), indicating great discrimination ability. The correlation between methylation level and gene expression of these 11 genes is shown in the appendix (Supplementary Figure.S1).
TABLE 2. Eleven mrDEGs identified as a DNA methylation signature prognostic model.
3.4 DNA methylation feature model for predicting the OS of LUAD patients
According to the correlation coefficients of the eleven mrDEGs obtained by multivariate Cox regression analysis, we established a feature score formula.
The LUAD patients were ranked according to their calculated methylation-related feature scores and divided into higher-risk (n = 226) and lower-risk groups (n = 227) according to the median. The K-M curve showed that the median OS of the higher-risk group was significantly shorter than that of the lower-risk group (log rank test p < 0.0001) (Figure 4A). We also analyzed the distribution of the methylation feature scores, patient survival statuses and methylated β values in LUAD patients, as well as the methylation β value spectra of 11 DNA methylation feature genes (Figures 4B–E).
FIGURE 4. The DNA methylation feature for predicting OS prediction in LUAD patients. (A) K-M assessed OS based on the DNA methylation feature. The LUAD patients were divided into the lower-risk (n = 226) and higher-risk (n = 227) subgroups according to the median of the methylation scores. Log rank test was used between curves (p < 0.0001). (B) The distribution of feature scores for DNA methylation feature of patients. (C–D) The distribution of survival status of LUAD patients in the lower- and higher-risk groups (Chi-square test, p < 0.0001). (E) The methylation β value spectrum of 11 DNA methylation feature genes.
3.5 DNA methylation feature model with clinicopathological features
First, we attempted to validate the correlation between DNA methylation levels and clinicopathological features in LUAD. We then performed unsupervised consistent clustering of the 641 most variable DNA methylation probes in 155 samples (with complete clinical information) into 4 clusters: CGI Methylator Phenotype (CIMP) high, CIMP medium high, CIMP medium low and CIMP low (Figure 5A). The average methylation levels among the different clusters were significant (p < 2.2e-16) (Figure 5B). DNA methylation was significantly correlated with tumor subtype (X-squared = 31.457, p = 2.073e-05) (Figure 5C). Most of magnoid tumors were enriched in the CIMP low and CIMP medium low clusters, while the CIMP high group had more squamoid tumors Furthermore, DNA methylation showed a trend related to tumor recurrence; however, there was no statistical Figure 5A Significance (p = 0.2864) (Figure 5D).
FIGURE 5. The DNA methylation feature with clinicopathological features. (A) Unsupervised cluster analysis of the methylation levels in LUAD. A total of 155 samples are presented in rows, and 641 CpG loci with the largest variation (mean methylation level β < 0.05 in normal samples and standard deviation σ > 0.20 in tumor samples) are listed. The 4 identified clusters are represented as CIMP high (n = 28), CIMP medium high (n = 45), CIMP medium low (n = 52) and CIMP low (n = 30). (B) Significant differences (p < 0.0001) in the methylation levels of the 4 clusters. (C) The sample distributions in terms of tumor subtype (Chi-square test, p < 0.0001). (D) The sample distributions in terms of recurrence (Chi-square test, p = 0.2864). (E) Forest map: multivariate Cox regression analysis was used to analyze the prognostic values of age, gender, smoking, tumor expression subtype, tumor recurrence, TNM stage, feature score and other clinicopathological features in 391 cases.
Then, we analyzed the correlation between DNA methylation feature scores and clinicopathological features. The results showed that the methylation feature score was significantly correlated with tumor subtype (Table 3, p = 0.032). In addition, the DNA methylation feature score was associated with smoking history in LUAD patients (Table 3, p = 0.004). Previous studies have shown that smoking is associated with methylation levels. For example, hypomethylation at cg05575921 in the aryl hydrocarbon receptor repressor gene was strongly associated with the smoking behavior of an individual (Jamieson et al., 2020). Therefore, the hypothesis that prognostic signals are related to smoking is reasonable. To delve into the effects of DNA methylation and clinicopathological features on prognosis, 391 patients with complete clinicopathological features were analyzed in a Cox regression model. The forest map showed that the feature score (HR = 1.69, 95% CI = 1.17–2.45, p = 0.006), the TNM staging system, and tumor recurrence (HR = 2.58, 95% CI = 1.79–3.72, p < 0.0001) were independent prognostic factors for LUAD, while smoking history and tumor subtype were not (Figure 5E). DNA methylation is closely related to tumor immune microenvironment (Chiappinelli and Baylin, 2022; Li et al., 2022). We investigated the relationship between methylation feature scores and tumor immune infiltration. In the group with low methylation feature scores, we observed an increase in monocyte, dendritic cell (resting) and mast cell (resting) infiltration, as well as a decrease in macrophage (M0) cell infiltration. However, there is no changes of immune effector cells observed (Supplementary Figure S2).
TABLE 3. The correlation between DNA methylation feature scores and clinicopathological feature.
3.6 Prognostic value of the DNA methylation feature score is independent of TNM stage and cancer recurrence
Since a high feature score, tumor recurrence and a high TNM stage were independent adverse prognostic factors for LUAD (Figure 5D), we performed a combined analysis between DNA methylation features and the other 2 influencing factors. We found that the prognosis of patients in the higher-feature score group was poorer, whether in the recurrence (log rank test, p < 0.0009) or nonrecurrence subgroup (log rank test, p = 0.02) (Figure 6A). In the combined analysis of TNM stages and DNA methylation features, we found that patients in the lower TNM stage (I and II) subgroups had a notably worse prognosis when they were also in the high-feature score subgroups (p = 2e-05) but not in the higher TNM stage (III and IV) subgroups (p = 0.3) (Figure 6B). We further classified patients at low TNM stages and found that the p value of the K-M curve in the stage I subgroup (p = 9e-04) was more significant than that in the stage II subgroup (p = 0.01) (Figures 6C,D). These results suggested that the DNA methylation feature score was more valuable in patients at a lower TNM stage. AUC-ROC analysis was used to evaluate the sensitivity and specificity of the prediction model (Figure 6E). The combination of this feature score and TNM stage was significantly superior to that of TNM stage alone (0.697 vs 0.658, p = 0.0275) or feature score alone (0.697 vs 0.603, p = 0.0001). These results suggested that the combination of the DNA methylation feature score and TNM stage might help to improve OS prediction in LUAD patients.
FIGURE 6. The prognostic value of DNA methylation feature score was not associated with tumor recurrence status and TNM stage. (A) K-M analysis of OS based on the DNA methylation feature and recurrence status. The LUAD patients were divided into the lower- and higher-risk subgroups according to the median of the methylation feature scores, with or without recurrence. Log-rank test (p < 0.0001). (B) K-M analysis of OS based on the DNA methylation feature and TNM stage. The LUAD patients were divided into the lower- and higher-risk subgroups according to the median of the methylation feature scores, and divided into the low (stage I+ II, n = 350) and high (stage III+IV, n = 98) stages according to the TNM stage (Log-rank test, p < 0.0001). (C) K-M curves for patients in the TNM stage I subgroup (n = 350). (D) K-M curves for patients in the TNM stage II subgroup (n = 98). (E) ROC analysis assessed the sensitivity and specificity of DNA methylation feature score, TNM stage and the combination of the 2 factors in predicting OS.
3.7 The expression of the eleven identified genes
The prognostic methylation signature consists of 11 genes. Three of them are hypermethylated in LUAD (TACR3, EDN3 and SLC27A6). We selected the broad-spectrum demethylation drug 5-aza-2′-deoxycytidine (5-aza) to treat LUAD cells (A549, PC9 and H1975) and measured the mRNA levels of TACR3, EDN3 and SLC27A6 by qPCR 4 days after treatment. The results showed that, compared with the control group, the transcription levels of these 3 genes were significantly increased after treatment with 5-aza (Figures 7A–C), suggesting that the transcriptional regulation of TACR3, EDN3 and SLC27A6 was related to promoter methylation. The other 8 genes showed low methylation and high mRNA expression in LUAD cells (MYEOV, NLRP13, SLC34A1, NEUROD4, HMGB4, KCNU1, GABRA5 and TRPM8). Since there is no broad-spectrum drug to improve DNA methylation, we used the Human Protein Atlas (HPA) (https://www.proteinatlas.org/) to verify the expression of the proteins encoded by these genes in LUAD and normal tissues. Five of them (MYEOV, NLRP13, SLC34A1, NEUROD4 and HMGB4) had protein expression data in the database, and MYEOV, NLRP13 and SLC34A1 were highly expressed in tumor tissues (Figure 7D), while NEUROD4 and HMGB4 showed no significant difference (Figure 7E).
FIGURE 7. The expression of TACR3, EDN3 and SLC27A6 is related to promoter region methylation. (A–C) qPCR was used to detect the mRNA levels of TACR3, EDN3 and SLC27A6 in A549, PC9 and H1975 cells before and after the 5-Aza-2′-deoxycytidine treatment. (D–E) Immunohistochemistry images obtained from the HPA database demonstrated the protein expression of the 5 hypomethylated genes.
4 Discussion
DNA methylation, as one of the most studied epigenetic alterations related to tumor phenotype, is of great significance for tumor research (Mazor et al., 2016; Biswas S Rao, 2017). Previous research showed that the overall DNA methylation pattern in tumor cell genomes is hypomethylation, while many CGIs associated with promoters showed focal hypermethylation (Hansen et al., 2011). Promoter hypermethylation was associated with tumor suppressor-related gene silencing, while the hypomethylation of the tumor cell genome could increase genomic instability (Feinberg and Vogelstein, 1983; Goelz et al., 1985; Pfeifer, 2018). Abnormal DNA methylation can be used not only as a target for tumor therapy but also as a biomarker for diagnosis and prognosis (Yamashita et al., 2018; Szejniuk et al., 2019). With the public information provided in TCGA, we conducted a comprehensive analysis and identified 11 methylation-related genes (TACR3, SLC27A6, EDN3, TRPM8, MYEOV, NLRP13, KCNU1, NEUROD4, GABRA5, SLC34A1 and HMGB4) to predict the prognosis of LUAD. These genes not only were differentially methylated and expressed in LUAD tumor tissues from TCGA but also were related to the prognosis of patients. The survival curves showed that there was a significant difference in the survival curve between the higher-risk and lower-risk groups, especially in patients with early LUAD. Some of the 11 identified methylation-related genes have been shown to be abnormally expressed and important in cancer or other diseases. For example, TRPM8 is a calcium permeability channel abnormally expressed in multiple malignant tumors. There is evidence that TRPM8 plays a major role in promoting cell invasion and preventing replicative senescence (Yee, 2016). In the present study, TRPM8 showed hypermethylation and low RNA expression in LUAD samples from TCGA, and the hypermethylation and low expression of TRPM8 were associated with long survival. Previous studies have shown that decreased expression or inactivation of EDN3 can inhibit the migration of cancer cells and improve survival (Wang et al., 2013; Kim et al., 2017). Our results revealed that EDN3 was hypermethylated and expressed at low levels in LUAD, which was associated with longer OS. Another gene, MYEOV, is a region of cancer-associated genomic amplification. The amplification of this gene was reported to promote the progression of NSCLC pancreatic ductal adenocarcinoma and colorectal cancer (Lawlor et al., 2010; Fang et al., 2019). Subsequent mechanistic studies showed that the overexpression of MYEOV might be regulated by promoter hypomethylation (Liang et al., 2020). In accordance with the above results, our studies showed that MYEOV was hypermethylated and expressed at low levels in LUAD and that hypermethylation was positively correlated with survival time. As for other genes in the DNA methylation feature. TACR3 was found to be highly elevated in endometrial carcinoma. Although the role of TAC1-TACR3 axis is not clear. Haixu et al. found that highly methylated TAC1 promoted the development of endometrial carcinoma through the deregulation of TAC (Xu et al., 2018). Kyoichi Obata et al. found that TACR3 protein showed significant and significant overexpression at the onset of bone matrix invasion in oral squamous cell carcinoma (Obata et al., 2016). SLC27A6 is used as a predictor in the genetic analysis of colorectal cancer, prostate cancer, pancreatic cancer, and other tumors (Mohammed et al., 2019; Uhan et al., 2020; Verma et al., 2020; Zhong et al., 2021). However, basic research on this topic is still very limited. As a neuron differentiation factor, NEUROD4 has been reported overexpressed in neuroendocrine tumors. Studies have shown that the continuous expression of NEUROD4 in neuronal cells may be related to the regeneration of neural cells, and its expression level gradually decreases with the maturation of neurons (Masserdotti et al., 2015; Cecil et al., 2016). SLC27A1 is rarely studied in tumors, but in recent years, some articles have pointed out that SLC27A1 is highly expressed in melanoma and breast cancer and enhances tumor invasion, migration, and growth (Kwaepila et al., 2006; Zhang et al., 2018). The expression of GABRA5, which encodes the α 5-GABAA receptor, has a synthetic lethal role in MYC-driven medulloblastoma (Sengupta et al., 2014). NLRP13, KCNU1, and HMGB4, although not as studied in tumors compared to the other genes in the model, need further exploration. Many articles verify the effectiveness of their own prediction formulas by comparing them with TNM staging (Peng et al., 2020; Zhu et al., 2021; Qiu et al., 2022). In our research, the DNA methylation feature score is an independent predictor and is not associated with the TNM stage. ROC curve analysis showed that the combination of the DNA methylation feature score and TNM stage was better for prognosis than TNM stage alone, suggesting that the combination of the 2 might help to improve the prediction of OS in LUAD patients. In addition, we found that the transcription levels of TACR3, EDN3 and SLC27A6 in LUAD cells were significantly increased by treatment with broad-spectrum demethylating drugs. These results suggested that the low expression of these genes was related to promoter hypermethylation. At the same time, the expression of some hypomethylated genes (MYEOV, NLRP13 and SLC34A1)in immunohistochemical sections of LUAD was significantly stronger than that in lung tissues. These results showed that the identified genes are worthy of further study as biomarkers of methylation in LUAD.
Because of its noninvasive and fast characteristics, detecting circulating tumor DNA (ctDNA) in blood to monitor epigenetic changes in tumor DNA has become a very promising technology. Although this technology is not sufficiently mature, blood testing based on a single DNA methylation biomarker has been approved (Frankell and Jamal-Hanjani, 2022). Our study shows that the methylation signals of these 11 genes may be used as candidate markers to detect ctDNA methylation in LUAD patients. This model can predict the prognosis of patients with low cost and high efficiency.
However, our research still has several limitations. First, DNA methylation biomarkers are not effective in predicting advanced LUAD. Considering the small sample size of the advanced LUAD group, the results will have a certain deviation. Second, because it was difficult to obtain data with a sufficient sample size and consistent methylation detection platform, we did not use other datasets to verify the methylation formula. However, considering the large sample size of this study, this model was less likely to be an accidental feature of methylation noise but more likely to be a determinant of LUAD survival. Finally, our basic experiments were limited. We did not regulate the specific methylation site of genes. Further experimental studies on these genes will help to determine further their therapeutic potential.
5 Conclusion
In conclusion, we explored the characteristics of DNA methylation in LUAD. Furthermore, we confirmed a DNA methylation feature consisting of 11 genes. DNA methylation is associated with the survival of LUAD patients and can provide a better OS predictive ability when combined with TNM stage. Unfortunately, it was difficult for us to obtain a sufficient sample size and consistent methylation detection platform data to verify this methylation formula. However, our experiments indicated that the transcription of the hypermethylated genes was increased after demethylation with 5-aza, suggesting the validity of these results and indicating the potential value of these 11 genes in the study of LUAD prognosis.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
All authors have contributed significantly. LX, CX contributed to the conception of the study. LX, ZH and HX performed the experiment. LX, ZH, ZZ and CX contributed significantly to analysis and manuscript preparation. LX, ZH and JL performed the data analyses and wrote the manuscript. LX, ZZ, CX helped perform the analysis with constructive discussions. All authors approve of the version to be submitted finally.
Funding
This work was supported by National Natural Science Foundation of China (grant number 81972852), Key Research & Development Project of Hubei Province (grant number 2020BCA069), Health Commission of Hubei Province Medical Leading Talent Project, and Chinese Society of Clinical Oncology TopAlliance Tumor Immune Research Fund (grant number Y-JS2019-036).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.970507/full#supplementary-material
References
Biswas S Rao, C. M. (2017). Epigenetics in cancer: Fundamentals and beyond. Pharmacol. Ther. 173, 118–134. doi:10.1016/j.pharmthera.2017.02.011
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin. 68 (6), 394–424. doi:10.3322/caac.21492
Cecil, C. A., Walton, E., Smith, R. G., Viding, E., McCrory, E. J., Relton, C. L., et al. (2016). DNA methylation and substance-use risk: A prospective, genome-wide study spanning gestation to adolescence. Transl. Psychiatry 6 (12), e976. doi:10.1038/tp.2016.247
Chi, Z., Kong, C. C., Wang, Z. Z., Wang, Z., Liu, G. L., Hu, Z., et al. (2021). The signaling pathways involved in metabolic regulation and stress responses of the yeast-like fungi Aureobasidium spp. Biotechnol. Adv. 55, 107898. doi:10.1016/j.biotechadv.2021.107898
Chiappinelli, K. B., and Baylin, S. B. (2022). Inhibiting DNA methylation improves antitumor immunity in ovarian cancer. J. Clin. Invest. 132 (14), e160186. doi:10.1172/JCI160186
Christman, J. K. (2002). 5-Azacytidine and 5-aza-2'-deoxycytidine as inhibitors of DNA methylation: Mechanistic studies and their implications for cancer therapy. Oncogene 21 (35), 5483–5495. doi:10.1038/sj.onc.1205699
Edge, S. B., and Compton, C. C. (2010). The American Joint committee on cancer: The 7th edition of the AJCC cancer staging manual and the future of TNM. Ann. Surg. Oncol. 17 (6), 1471–1474. doi:10.1245/s10434-010-0985-4
Enfield, K. S. S., Marshall, E. A., Anderson, C., Ng, K. W., Rahmati, S., Xu, Z., et al. (2019). Epithelial tumor suppressor ELF3 is a lineage-specific amplified oncogene in lung adenocarcinoma. Nat. Commun. 10 (1), 5438. doi:10.1038/s41467-019-13295-y
Fang, L., Wu, S., Zhu, X., Cai, J., Wu, J., He, Z., et al. (2019). MYEOV functions as an amplified competing endogenous RNA in promoting metastasis by activating TGF-β pathway in NSCLC. Oncogene 38 (6), 896–912. doi:10.1038/s41388-018-0484-9
Feinberg, A. P. (2007). Phenotypic plasticity and the epigenetics of human disease. Nature 447 (7143), 433–440. doi:10.1038/nature05919
Feinberg, A. P., and Vogelstein, B. (1983). Hypomethylation distinguishes genes of some human cancers from their normal counterparts. Nature 301 (5895), 89–92. doi:10.1038/301089a0
Folkman, J. (1971). Tumor angiogenesis: Therapeutic implications. N. Engl. J. Med. 285 (21), 1182–1186. doi:10.1056/NEJM197111182852108
Frankell, A., and Jamal-Hanjani, M. (2022). When can we be confident of surgical cure with ctDNA? Nat. Rev. Clin. Oncol. doi:10.1038/s41571-022-00664-8
Goelz, S. E., Vogelstein, B., Hamilton, S. R., and Feinberg, A. P. (1985). Hypomethylation of DNA from benign and malignant human colon neoplasms. Science 228 (4696), 187–190. doi:10.1126/science.2579435
Goldstraw, P., Ball, D., Jett, J. R., Le Chevalier, T., Lim, E., Nicholson, A. G., et al. (2011). Non-small-cell lung cancer. Lancet 378 (9804), 1727–1740. doi:10.1016/S0140-6736(10)62101-0
Goto, Y., Shinjo, K., Kondo, Y., Shen, L., Toyota, M., Suzuki, H., et al. (2009). Epigenetic profiles distinguish malignant pleural mesothelioma from lung adenocarcinoma. Cancer Res. 69 (23), 9073–9082. doi:10.1158/0008-5472.Can-09-1595
Hansen, K. D., Timp, W., Bravo, H. C., Sabunciyan, S., Langmead, B., McDonald, O. G., et al. (2011). Increased methylation variation in epigenetic domains across cancer types. Nat. Genet. 43 (8), 768–775. doi:10.1038/ng.865
Huo, M., Zhang, J., Huang, W., and Wang, Y. (2021). Interplay among metabolism, epigenetic modifications, and gene expression in cancer. Front. Cell Dev. Biol. 9, 793428. doi:10.3389/fcell.2021.793428
Ikeda, K., Shiraishi, K., Eguchi, A., Shibata, H., Yoshimoto, K., Mori, T., et al. (2013). Long interspersed nucleotide element 1 hypomethylation is associated with poor prognosis of lung adenocarcinoma. Ann. Thorac. Surg. 96 (5), 1790–1794. doi:10.1016/j.athoracsur.2013.06.035
Jamieson, E., Korologou-Linden, R., Wootton, R. E., Guyatt, A. L., Battram, T., Burrows, K., et al. (2020). Smoking, DNA methylation, and lung function: A mendelian randomization analysis to investigate causal pathways. Am. J. Hum. Genet. 106 (3), 315–326. doi:10.1016/j.ajhg.2020.01.015
Jung, S. H. (2014). Stratified Fisher's exact test and its sample size calculation. Biom. J. 56 (1), 129–140. doi:10.1002/bimj.201300048
Kaishang, Z., Xue, P., Shaozhong, Z., Yingying, F., Yan, Z., Chanjun, S., et al. (2018). Elevated expression of Twinfilin-1 is correlated with inferior prognosis of lung adenocarcinoma. Life Sci. 215, 159–169. doi:10.1016/j.lfs.2018.10.067
Karagiota, A., Chachami, G., and Paraskeva, E. (2022). Lipid metabolism in cancer: The role of acylglycerolphosphate acyltransferases (AGPATs). Cancers (Basel) 14 (1), 228. doi:10.3390/cancers14010228
Kim, I. S., Heilmann, S., Kansler, E. R., Zhang, Y., Zimmer, M., Ratnakumar, K., et al. (2017). Microenvironment-derived factors driving metastatic plasticity in melanoma. Nat. Commun. 8, 14343. doi:10.1038/ncomms14343
Küster, M. M., Schneider, M. A., Richter, A. M., Richtmann, S., Winter, H., Kriegsmann, M., et al. (2020). Epigenetic inactivation of the tumor suppressor IRX1 occurs frequently in lung adenocarcinoma and its silencing is associated with impaired prognosis. Cancers 12 (12), E3528. doi:10.3390/cancers12123528
Kwaepila, N., Burns, G., and Leong, A. S. (2006). Immunohistological localisation of human FAT1 (hFAT) protein in 326 breast cancers. Does this adhesion molecule have a role in pathogenesis? Pathology 38 (2), 125–131. doi:10.1080/00313020600559975
Lawlor, G., Doran, P. P., MacMathuna, P., and Murray, D. W. (2010). MYEOV (myeloma overexpressed gene) drives colon cancer cell migration and is regulated by PGE2. J. Exp. Clin. Cancer Res. 29 (1), 81. doi:10.1186/1756-9966-29-81
Li, C. H., and Liao, C. C. (2021). The metabolism reprogramming of microRNA let-7-mediated glycolysis contributes to autophagy and tumor progression. Int. J. Mol. Sci. 23 (1), 113. doi:10.3390/ijms23010113
Li, J., Liang, Y., Fan, J., Xu, C., Guan, B., Zhang, J., et al. (2022). DNA methylation subtypes guiding prognostic assessment and linking to responses the DNA methyltransferase inhibitor SGI-110 in urothelial carcinoma. BMC Med. 20 (1), 222. doi:10.1186/s12916-022-02426-w
Li, Y., Gu, J., Xu, F., Zhu, Q., Ge, D., and Lu, C. (2018). Transcriptomic and functional network features of lung squamous cell carcinoma through integrative analysis of GEO and TCGA data. Sci. Rep. 8 (1), 15834. doi:10.1038/s41598-018-34160-w
Li, Y., Sun, N., Lu, Z., Sun, S., Huang, J., Chen, Z., et al. (2017). Prognostic alternative mRNA splicing signature in non-small cell lung cancer. Cancer Lett. 393, 40–51. doi:10.1016/j.canlet.2017.02.016
Li, Z., Zhang, R., Yang, X., Zhang, D., Li, B., Zhang, D., et al. (2019). Analysis of gene expression and methylation datasets identified ADAMTS9, FKBP5, and PFKBF3 as biomarkers for osteoarthritis. J. Cell. Physiol. 234 (6), 8908–8917. doi:10.1002/jcp.27557
Lian, H., Han, Y. P., Zhang, Y. C., Zhao, Y., Yan, S., Li, Q. F., et al. (2019). Integrative analysis of gene expression and DNA methylation through one-class logistic regression machine learning identifies stemness features in medulloblastoma. Mol. Oncol. 13 (10), 2227–2245. doi:10.1002/1878-0261.12557
Liang, E., Lu, Y., Shi, Y., Zhou, Q., and Zhi, F. (2020). MYEOV increases HES1 expression and promotes pancreatic cancer progression by enhancing SOX9 transactivity. Oncogene 39 (41), 6437–6450. doi:10.1038/s41388-020-01443-4
Liu, C., Li, Y., Wei, M., Zhao, L., and Yu Y Li, G. (2019). Identification of a novel glycolysis-related gene signature that can predict the survival of patients with lung adenocarcinoma. Cell Cycle 18 (5), 568–579. doi:10.1080/15384101.2019.1578146
Masserdotti, G., Gillotin, S., Sutor, B., Drechsel, D., Irmler, M., Jørgensen, H. F., et al. (2015). Transcriptional mechanisms of proneural factors and REST in regulating neuronal reprogramming of astrocytes. Cell Stem Cell 17 (1), 74–88. doi:10.1016/j.stem.2015.05.014
Mazor, T., Pankov, A., and Song J S Costello, J. F. (2016). Intratumoral heterogeneity of the epigenome. Cancer Cell 29 (4), 440–451. doi:10.1016/j.ccell.2016.03.009
Mohammed, A., Janakiram, N. B., Suen, C., Stratton, N., Lightfoot, S., Singh, A., et al. (2019). Targeting cholecystokinin-2 receptor for pancreatic cancer chemoprevention. Mol. Carcinog. 58 (10), 1908–1918. doi:10.1002/mc.23084
Obata, K., Shimo, T., Okui, T., Matsumoto, K., Takada, H., Takabatake, K., et al. (2016). Tachykinin receptor 3 distribution in human oral squamous cell carcinoma. Anticancer Res. 36 (12), 6335–6341. doi:10.21873/anticanres.11230
Pandis, N. (2016). The chi-square test. Am. J. Orthod. Dentofac. Orthop. 150 (5), 898–899. doi:10.1016/j.ajodo.2016.08.009
Peng, Y., Wu, Q., Wang, L., Wang, H., and Yin, F. (2020). A DNA methylation signature to improve survival prediction of gastric cancer. Clin. Epigenetics 12 (1), 15. doi:10.1186/s13148-020-0807-x
Pfeifer, G. P. (2018). Defining driver DNA methylation changes in human cancer. Int. J. Mol. Sci. 19 (4), E1166. doi:10.3390/ijms19041166
Qiu, Z., Ji, J., Xu, Y., Zhu, Y., Gao, C., Wang, G., et al. (2022). Common DNA methylation changes in biliary tract cancers identify subtypes with different immune characteristics and clinical outcomes. BMC Med. 20 (1), 64. doi:10.1186/s12916-021-02197-w
Schröder, M. S., Culhane, A. C., Quackenbush, J., and Haibe-Kains, B. (2011). survcomp: an R/Bioconductor package for performance assessment and comparison of survival models. Bioinformatics 27 (22), 3206–3208. doi:10.1093/bioinformatics/btr511
Sengupta, S., Weeraratne, S. D., Sun, H., Phallen, J., Rallapalli, S. K., Teider, N., et al. (2014). α5-GABAA receptors negatively regulate MYC-amplified medulloblastoma growth. Acta Neuropathol. 127 (4), 593–603. doi:10.1007/s00401-013-1205-7
Singh, S., Narayanan, S. P., Biswas, K., Gupta, A., Ahuja, N., Yadav, S., et al. (2017). Intragenic DNA methylation and BORIS-mediated cancer-specific splicing contribute to the Warburg effect. Proc. Natl. Acad. Sci. U. S. A. 114 (43), 11440–11445. doi:10.1073/pnas.1708447114
Sui, J., Li, Y-H., Zhang, Y-Q., Li, C-Y., Shen, X., Yao, W-Z., et al. (2016). Integrated analysis of long non-coding RNA-associated ceRNA network reveals potential lncRNA biomarkers in human lung adenocarcinoma. Int. J. Oncol. 49 (5), 2023–2036. doi:10.3892/ijo.2016.3716
Szejniuk, W. M., Robles, A. I., McCulloch, T., Falkmer U G I Røe, O. D., and Roe, O. D. (2019). Epigenetic predictive biomarkers for response or outcome to platinum-based chemotherapy in non-small cell lung cancer, current state-of-art. Pharmacogenomics J. 19 (1), 5–14. doi:10.1038/s41397-018-0029-1
Tozzi, R., Masci, F., and Pezzopane, M. (2020). A stress test to evaluate the usefulness of Akaike information criterion in short-term earthquake prediction. Sci. Rep. 10 (1), 21153. doi:10.1038/s41598-020-77834-0
Uhan, S., Zidar, N., and Tomažič A Hauptman, N. (2020). Hypermethylated promoters of genes UNC5D and KCNA1 as potential novel diagnostic biomarkers in colorectal cancer. Epigenomics 12 (19), 1677–1688. doi:10.2217/epi-2020-0118
Verma, S., Shankar, E., Chan, E. R., and Gupta, S. (2020). Metabolic reprogramming and predominance of solute carrier genes during acquired enzalutamide resistance in prostate cancer. Cells 9 (12), E2535. doi:10.3390/cells9122535
Wang, R., Löhr, C. V., Fischer, K., Dashwood, W. M., Greenwood, J. A., Ho, E., et al. (2013). Epigenetic inactivation of endothelin-2 and endothelin-3 in colon cancer. Int. J. Cancer 132 (5), 1004–1012. doi:10.1002/ijc.27762
Wilkerson, M. D., and Hayes, D. N. (2010). ConsensusClusterPlus: A class discovery tool with confidence assessments and item tracking. Bioinformatics 26 (12), 1572–1573. doi:10.1093/bioinformatics/btq170
Xu, H., Sun, Y., Ma, Z., Xu, X., Qin, L., and Luo, B. (2018). LOC134466 methylation promotes oncogenesis of endometrial carcinoma through LOC134466/hsa-miR-196a-5p/TAC1 axis. Aging (Albany NY) 10 (11), 3353–3370. doi:10.18632/aging.101644
Yamashita, K., Hosoda, K., Nishizawa, N., Katoh, H., and Watanabe, M. (2018). Epigenetic biomarkers of promoter DNA methylation in the new era of cancer treatment. Cancer Sci. 109 (12), 3695–3706. doi:10.1111/cas.13812
Yang, S., Huang, Y., and Zhao, Q. (2022). Epigenetic alterations and inflammation as emerging use for the advancement of treatment in non-small cell lung cancer. Front. Immunol. 13, 878740. doi:10.3389/fimmu.2022.878740
Yee, N. S. (2016). TRPM8 ion channels as potential cancer biomarker and target in pancreatic cancer. Adv. Protein Chem. Struct. Biol. 104, 127–155. doi:10.1016/bs.apcsb.2016.01.001
Zhang, H., Kong, Q., Wang, J., Jiang, Y., and Hua, H. (2020). Complex roles of cAMP-PKA-CREB signaling in cancer. Exp. Hematol. Oncol. 9 (1), 32. doi:10.1186/s40164-020-00191-1
Zhang, M., Di Martino, J. S., Bowman, R. L., Campbell, N. R., Baksh, S. C., Simon-Vermot, T., et al. (2018). Adipocyte-Derived lipids mediate melanoma progression via FATP proteins. Cancer Discov. 8 (8), 1006–1025. doi:10.1158/2159-8290.Cd-17-1371
Zhao, J., Guo, C., Ma, Z., Liu, H., Yang, C., and Li, S. (2020). Identification of a novel gene expression signature associated with overall survival in patients with lung adenocarcinoma: A comprehensive analysis based on TCGA and geo databases. Lung Cancer 149, 90–96. doi:10.1016/j.lungcan.2020.09.014
Zhong, X., Yang, Y., Li, B., Liang, P., Huang, Y., Zheng, Q., et al. (2021). Downregulation of SLC27A6 by DNA hypermethylation promotes proliferation but suppresses metastasis of nasopharyngeal carcinoma through modulating lipid metabolism. Front. Oncol. 11, 780410. doi:10.3389/fonc.2021.780410
Keywords: DNA methylation, integrative analysis, lung adenocarcinoma, prognosis, TCGA
Citation: Xu L, Huang Z, Zeng Z, Li J, Xie H and Xie C (2022) An integrative analysis of DNA methylation and gene expression to predict lung adenocarcinoma prognosis. Front. Genet. 13:970507. doi: 10.3389/fgene.2022.970507
Received: 16 June 2022; Accepted: 03 August 2022;
Published: 29 August 2022.
Edited by:
Zhongxiao Wang, Biogen Idec, United StatesReviewed by:
Yajuan Zhou, Hubei Cancer Hospital, ChinaWanzun Lin, Shanghai Cancer Center, Fudan University, China
Copyright © 2022 Xu, Huang, Zeng, Li, Xie and Xie. 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: Conghua Xie, chxie_65@whu.edu.cn
†These authors contributed equally to this work