Identification of TAF1, SAT1, and ARHGEF9 as DNA methylation biomarkers for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a major cause of cancer-related deaths worldwide. More than 90% of primary HCC is HCC. Hepatitis C virus (HCV) infection and alcohol consumption have been widely accepted as two major risk factors for developing HCC. Herein, we aimed to identify DNA methylation genes related to both HCV infection and alcohol consumption. In this study, we identified methylation genes that were associated with the risk of HCV infection and alcohol consumption, respectively, by a large-scale bioinformatic analysis. Through PPI network analysis, we revealed the associations between the two types of genes and found six hub genes—TAF1, SAT1, Phospholipase C-beta 2, FGD1, ARHGAP4, and ARHGEF9—that may be associated with both HCV infection and alcohol consumption. Gene Ontology enrichment analysis was used to analyze the function which these genes in the network enriched. Among them, TAF1, SAT1, and ARHGEF9 were methylated genes that have been found to be related to tumor progression in HCC patients. Through independent data sets, we verified the methylation pattern of these six genes in HCC samples that had both HCV infection and alcohol consumption risks. Furthermore, we found that three of the six methylated genes were also associated with the prognosis of HCC patients. To summarize, we identified six hub genes that were associated with both HCV infection and alcohol consumption in the progress of HCC. The six methylation genes that might play an important role in both HCV infection and alcohol consumption would be potential therapy targets for HCC.     

Identification of IMPA2 as the hub gene associated with colorectal cancer and liver metastasis by integrated bioinformatics analysis

Background and ObjectivesColorectal cancer (CRC) is one of the most common malignant tumors worldwide with high incidence and mortality rate, while colorectal liver metastasis (CRLM) is one of the major causes of cancer-related deaths. Therefore, the present study aims to identify the hub gene associated with CRC carcinogenesis and liver metastasis, and then explore its diagnostic and prognostic value as well as the potential regulation mechanism.MethodsThe overlapping differential co-expression genes among CRC, CRLM, and normal tissues were explored on the GSE49355 and GSE81582 datasets from the Gene Expression Omnibus (GEO) database by integrated bioinformatics analysis. Then, the hub prognostic genes were selected from the overlapping genes by univariate Cox proportional hazard analysis and online database Gene Expression Profiling Interactive Analysis 2 (GEPIA2). Subsequently, the clinical value of the hub genes was evaluated in the TCGA and GSE39582 cohorts. Finally, the underlying mechanisms of the hub gene regulating CRC carcinogenesis and metastasis were explored by Gene function annotation and DNA methylation analysis.ResultsInositol mono-phosphatase 2 (IMPA2) was identified as the hub gene associated with CRC carcinogenesis and liver metastasis. IMPA2 had an excellent diagnostic efficiency, and its expression was significantly decreased in CRC and liver metastasis samples, being positively correlated with poor prognosis. Moreover, its low expression was associated with AJCC stage III+IV, T4, N1+2, and M1. In addition, our results revealed that the potential mechanisms used by IMPA2 to mediate CRC carcinogenesis and metastasis could be associated with lipid metabolism and epithelial mesenchymal transition (EMT). Finally, IMPA2 expression could be regulated by DNA methylation.ConclusionsIMPA2 was identified and reported for the first time as a hub gene biomarker in the diagnosis and prognosis of CRC, which could regulate CRC carcinogenesis and liver metastasis through the regulation of lipid metabolism, EMT, and DNA methylation.     

DNA methylation-based diagnostic and prognostic biomarkers of nonsmoking lung adenocarcinoma patients

Currently, there are few studies on patients with nonsmoking lung adenocarcinoma, and the pathogenesis is still unclear. The role of DNA methylation in the pathogenesis of cancer is gradually being recognized. The purpose of this study was to determine the abnormal methylation genes and pathways involved in nonsmoking lung adenocarcinoma patients. Gene expression microarray data (GSE10072, GSE43458) and gene methylation microarray data (GSE62948) were downloaded from the Gene Expression Omnibus (GEO) database and differentially expressed genes were obtained through GEO2R. Next, we analyzed the function and enrichment of the selected genes using Database for Annotation, Visualization, and Integrated Discovery. The protein-protein interaction (PPI) networks were constructed using the Search Tool for the Retrieval of Interacting Genes database and visualized in Cytoscape. Finally, we performed module analysis of the PPI network using Molecular Complex Detection. And we obtained 10 hub genes by Cytoscape Centiscape. We analyzed the independent prognostic value of each hub gene in nonsmoking nonsmall cell lung cancer patients through Kaplan-Meier plotter. Seven hub genes (CXCL12, CDH1, CASP3, CREB1, COL1A1, ERBB2, and ENO2) were closely related to the overall survival time. This study provides an effective bioinformatics basis for further understanding the pathogenesis and prognosis of nonsmoking lung adenocarcinoma patients. Hub genes with prognostic value could be selected as effective biomarkers for timely diagnosis and prognostic of nonsmoking lung adenocarcinoma patients.     

Transcriptome and DNA methylation analysis reveals molecular mechanisms underlying intrahepatic cholangiocarcinoma progression

Intrahepatic cholangiocarcinoma (iCCA) is an aggressive malignancy with increasing incidence. It has been suggested that DNA methylation drives cancer development. However, the molecular mechanisms underlying iCCA progression and the roles of DNA methylation still remain elusive. In this study, weighted correlation networks were constructed to identify gene modules and hub genes associated with the tumour stage. We identified 12 gene modules, two of which were significantly positively or negatively related to the tumour stage, respectively. Key hub genes SLC2A1, CDH3 and EFHD2 showed increased expression across the tumour stage and were correlated with poor survival, whereas decrease of FAM171A1, ONECUT1 and PHYHIPL was correlated with better survival. Pathway analysis revealed hedgehog pathway was activated in CDH3 up-regulated tumours, and chromosome separation was elevated in tumours expressing high EFHD2. JAK-STAT pathway was overrepresented in ONECUT1 down-regulated tumours, whereas Rho GTPases-formins signalling was activated in PHYHIPL down-regulated tumours. Finally, significant negative associations between expression of EFHD2, PHYHIPL and promoter DNA methylation were detected, and alterations of DNA methylation were correlated with tumour survival. In summary, we identified key genes and pathways that may participate in progression of iCCA and proposed putative roles of DNA methylation in iCCA.     

Integrative analysis of DNA methylation and gene expression through machine learning identifies stomach cancer diagnostic and prognostic biomarkers

DNA methylation is an early event in tumorigenesis. Here, by integrative analysis of DNA methylation and gene expression and utilizing machine learning approaches, we introduced potential diagnostic and prognostic methylation signatures for stomach cancer. Differentially-methylated positions (DMPs) and differentially-expressed genes (DEGs) were identified using The Cancer Genome Atlas (TCGA) stomach adenocarcinoma (STAD) data. A total of 256 DMPs consisting of 140 and 116 hyper- and hypomethylated positions were identified between 443 tumour and 27 nontumour STAD samples. Gene expression analysis revealed a total of 2821 DEGs with 1247 upregulated and 1574 downregulated genes. By analysing the impact of cis and trans regulation of methylation on gene expression, a dominant negative correlation between methylation and expression was observed, while for trans regulation, in hypermethylated and hypomethylated genes, there was mainly a negative and positive correlation with gene expression, respectively. To find diagnostic biomarkers, we used 28 hypermethylated probes locating in the promoter of 27 downregulated genes. By implementing a feature selection approach, eight probes were selected and then used to build a support vector machine diagnostic model, which had an area under the curve of 0.99 and 0.97 in the training and validation (GSE30601 with 203 tumour and 94 nontumour samples) cohorts, respectively. Using 412 TCGA-STAD samples with both methylation and clinical data, we also identified four prognostic probes by implementing univariate and multivariate Cox regression analysis. In summary, our study introduced potential diagnostic and prognostic biomarkers for STAD, which demands further validation.     

Identification of EPHX2 and RMI2 as two novel key genes in cervical squamous cell carcinoma by an integrated bioinformatic analysis

Cervical cancer is the fourth most common malignancy in women worldwide and cervical squamous cell carcinoma (CESC) is the most common histological type of cervical cancer. The dysregulation of genes plays a significant role in cancer. In the present study, we screened out differentially expressed genes (DEGs) of CESC in the GSE63514 data set from the Gene Expression Omnibus database. An integrated bioinformatics analysis was used to select hub genes, as well as to investigate their related prognostic signature, functional annotation, methylation mechanism, and candidate molecular drugs. As a result, a total of 1907 DEGs were identified (944 were upregulated and 963 were downregulated). In the protein–protein interaction network, three hub modules and 30 hub genes were identified. And two hub modules and 116 hub genes were screened out from four CESC-related modules by the weighted gene coexpression network analysis. The gene ontology term enrichment analysis and Kyoto encyclopedia of genes and genomes pathway analysis were performed to better understand functions and pathways. Genes with a significant prognostic value were found by prognostic signature analysis. And there were five genes (EPHX2, CHAF1B, KIAA1524, CDC45, and RMI2) identified as significant CESC-associated genes after expression validation and survival analysis. Among them, EPHX2 and RMI2 were noted as two novel key genes for the CESC-associated methylation and expression. In addition, four candidate small molecule drugs for CESC (camptothecin, resveratrol, vorinostat, and trichostatin A) were defined. Further studies are required to explore these significant CESC-associated genes for their potentiality in diagnosis, prognosis, and targeted therapy.     

Identification of prognostic signature in cancer based on DNA methylation interaction network

Background

The identification of prognostic biomarkers for cancer patients is essential for cancer research. These days, DNA methylation has been proved to be associated with cancer prognosis. However, there are few methods which identify the prognostic markers based on DNA methylation data systematically, especially considering the interaction among DNA methylation sites.

Methods

In this paper, we first evaluated the stabilities of microRNA, mRNA, and DNA methylation data in prognosis of cancer. After that, a rank-based method was applied to construct a DNA methylation interaction network. In this network, nodes with the largest degrees (10% of all the nodes) were selected as hubs. Cox regression was applied to select the hubs as prognostic signature. In this prognostic signature, DNA methylation levels of each DNA methylation site are correlated with the outcomes of cancer patients. After obtaining these prognostic genes, we performed the survival analysis in the training group and the test group to verify the reliability of these genes.

Results

We applied our method in three cancers (ovarian cancer, breast cancer and Glioblastoma Multiforme).In all the three cancers, there are more common ones of prognostic genes selected from different samples in DNA methylation data, compared with gene expression data and miRNA expression data, which indicates the DNA methylation data may be more stable in cancer prognosis. Power-law distribution fitting suggests that the DNA methylation interaction networks are scale-free. And the hubs selected from the three networks are all enriched by cancer related pathways. The gene signatures were obtained for the three cancers respectively, and survival analysis shows they can distinguish the outcomes of tumor patients in both the training data sets and test data sets, which outperformed the control signatures.

Conclusions

A computational method was proposed to construct DNA methylation interaction network and this network could be used to select prognostic signatures in cancer.

     

DNA Methylation Patterns Can Estimate Nonequivalent Outcomes of Breast Cancer with the Same Receptor Subtypes

Breast cancer has various molecular subtypes and displays high heterogeneity. Aberrant DNA methylation is involved in tumor origin, development and progression. Moreover, distinct DNA methylation patterns are associated with specific breast cancer subtypes. We explored DNA methylation patterns in association with gene expression to assess their impact on the prognosis of breast cancer based on Infinium 450K arrays (training set) from The Cancer Genome Atlas (TCGA). The DNA methylation patterns of 12 featured genes that had a high correlation with gene expression were identified through univariate and multivariable Cox proportional hazards models and used to define the methylation risk score (MRS). An improved ability to distinguish the power of the DNA methylation pattern from the 12 featured genes (p = 0.00103) was observed compared with the average methylation levels (p = 0.956) or gene expression (p = 0.909). Furthermore, MRS provided a good prognostic value for breast cancers even when the patients had the same receptor status. We found that ER-, PR- or Her2- samples with high-MRS had the worst 5-year survival rate and overall survival time. An independent test set including 28 patients with death as an outcome was used to test the validity of the MRS of the 12 featured genes; this analysis obtained a prognostic value equivalent to the training set. The predict power was validated through two independent datasets from the GEO database. The DNA methylation pattern is a powerful predictor of breast cancer survival, and can predict outcomes of the same breast cancer molecular subtypes.     

NDRG2 and TLR7 as novel DNA methylation prognostic signatures for acute myelocytic leukemia

Acute myelocytic leukemia (AML) is an aggressive malignant tumor and typically fatal without treatment. Identification and development of novel biomarkers could be beneficial for the diagnosis and prognosis of AML patients. Here, we aimed to identify the accurate DNA methylation prognostic signatures for AML patients. The DNA methylation data of AML patients and corresponding clinical information were retrieved from The Cancer Genome Atlas database. CPG sites that correlates closely with the survival of the AML patients were identified and further combined into CPG sites pairs to screen the survival-related pairs. The prognostic signatures were identified by the C-index and forward search algorithms and validated by the verification group. Finally, the functional enrichment analysis was performed on these CPG sites. As a result, a total of 498 CPG sites associated with the overall survival of AML patients was obtained. A prognostic signature composed of 10 CPG sites pairs was obtained and validated. The functional enrichment analysis showed prognostic genes were mainly enriched in tumor protein processing, cell differentiation, blood leukocyte immunity, and platelet growth factor pathways. In summary, we identified two accurate prognostic methylation signatures (NDRG2 and TLR7), which would be served as a novel therapy target for AML.     

Hub gene identification and prognostic model construction for isocitrate dehydrogenase mutation in glioma

Our study attempted to identify hub genes related to isocitrate dehydrogenase (IDH) mutation in glioma and develop a prognostic model for IDH-mutant glioma patients. In a first step, ten hub genes significantly associated with the IDH status were identified by weighted gene coexpression analysis (WGCNA). The functional enrichment analysis demonstrated that the most enriched terms of these hub genes were cadherin binding and glutathione metabolism. Three of these hub genes were significantly linked with the survival of glioma patients. 328 samples of IDH-mutant glioma were separated into two datasets: a training set (N = 228) and a test set (N = 100). Based on the training set, we identified two IDH-mutant subtypes with significantly different pathological features by using consensus clustering. A 31 gene-signature was identified by the least absolute shrinkage and selection operator (LASSO) algorithm and used for establishing a differential prognostic model for IDH-mutant patients. In addition, the test set was employed for validating the prognostic model, and the model was proven to be of high value in classifying prognostic information of samples. The functional annotation revealed that the genes related to the model were mainly enriched in nuclear division, DNA replication, and cell cycle. Collectively, this study provided novel insights into the molecular mechanism of IDH mutation in glioma, and constructed a prognostic model which can be effective for predicting prognosis of glioma patients with IDH-mutation, which might promote the development of IDH target agents in glioma therapies and contribute to accurate prognostication and management in IDH-mutant glioma patients.     

Systematic prediction of key genes for ovarian cancer by co-expression network analysis

Ovarian cancer (OC) is the most lethal gynaecological malignancy, characterized by high recurrence and mortality. However, the mechanisms of its pathogenesis remain largely unknown, hindering the investigation of the functional roles. This study sought to identify key hub genes that may serve as biomarkers correlated with prognosis. Here, we conduct an integrated analysis using the weighted gene co-expression network analysis (WGCNA) to explore the clinically significant gene sets and identify candidate hub genes associated with OC clinical phenotypes. The gene expression profiles were obtained from the MERAV database. Validations of candidate hub genes were performed with RNASeqV2 data and the corresponding clinical information available from The Cancer Genome Atlas (TCGA) database. In addition, we examined the candidate genes in ovarian cancer cells. Totally, 19 modules were identified and 26 hub genes were extracted from the most significant module (R² = .53) in clinical stages. Through the validation of TCGA data, we found that five hub genes (COL1A1, DCN, LUM, POSTN and THBS2) predicted poor prognosis. Receiver operating characteristic (ROC) curves demonstrated that these five genes exhibited diagnostic efficiency for early-stage and advanced-stage cancer. The protein expression of these five genes in tumour tissues was significantly higher than that in normal tissues. Besides, the expression of COL1A1 was associated with the TAX resistance of tumours and could be affected by the autophagy level in OC cell line. In conclusion, our findings identified five genes could serve as biomarkers related to the prognosis of OC and may be helpful for revealing pathogenic mechanism and developing further research.     

An 11-gene-based prognostic signature for uveal melanoma metastasis based on gene expression and DNA methylation profile

Uveal melanoma (UM) is the most common intraocular tumor worldwide. We proposed to identify a vital gene signature that has prognostic value for UM metastasis. For this purpose, we obtained a published DNA methylation and gene expression data set associated with UM from the Gene Expression Omnibus. The genes whose aberrant expression significantly associated with UM patients’ metastasis-free survival (MFS) were identified by applying a univariate Cox proportional hazards model to the gene expression data set followed by a robust likelihood-based survival analysis to screen the optimal prognostic gene signatures (PGS). A formula for calculating the risk score that represents UM metastasis risk was constructed by including the PGSs’ expression values weighted by their regression coefficients, which were obtained by a multivariate Cox regression analysis. As a result, aberrant expression of 2884 genes were found to be significantly associated with UM patients’ MFS, which were referred to as MFSGs, and 11 out of those MFSGs, GJC1, TCEA1, MFSD3, FAF2, TLCD1, GPAA1, CYC1, ASAP1, JPH1, LDB3, and KDELR3, were identified as PGSs through which we could accurately separate UM samples with shorter MFS from those with longer MFS. By combining the DNA methylation data set and MFSGs, we further identified 265 MFSGs, which contained CpG sites that significantly hyper- or hypo-methylated in UM samples compared with control samples. Functional enrichment analysis and pathway crosstalk analysis of those genes indicated significant enrichment of cancer-related pathways. In conclusion, we identified an 11-gene-based prognostic signature and several gene biomarkers for UM metastasis, which should be helpful for selecting an appropriate treatment method for specific patients with UM.     

Analysis of methylation-driven genes for predicting the prognosis of patients with head and neck squamous cell carcinoma

To help provide evidence for prognosis prediction and personalized targeted therapy for patients with head and neck squamous cell carcinoma (HNSCC), we investigated prognosis-specific methylation-driven genes in HNSCC. Survival time data, RNA sequencing data, and methylation data for HNSCC patients were downloaded from The Cancer Genome Atlas. The MethylMix R package based on the β mixture model was utilized to screen genes with different methylation statuses in tumor tissues and adjacent normal tissues, and a total of 182 HNSCC-related methylation-driven genes were then identified. A survival prediction scoring model based on multivariate Cox analysis was developed to screen the genes related to the prognosis of HNSCC, and a linear risk model of the methylation status of six genes (INA, LINC01354, TSPYL4, MAGEB2, EPHX3, and ZNF134) was constructed. The prognostic values of the six genes were further independently explored by survival analysis combined with methylation and gene expression analyses. The 5-year survival rate in the high-risk group of patients in the test set was 30.4% (95% CI: 22.7%-40.8%) and that in the low-risk group of patients was 65.5% (95% CI: 56.1%-76.5%). The area under the receiver operating characteristic curve for the model was 0.723, which further verified the specificity and sensitivity of the model. In addition, subsequent combined survival analysis revealed that all six genes could be used as independent prognostic markers and thus might be potential drug targets. The innovative method provides new insight into the molecular mechanism and prognosis of HNSCC.     

ABCC9, NKAPL,and TMEM132C are potential diagnostic and prognostic markers in triple‐negative breast cancer

Triple‐negative breast cancer (TNBC) is a highly heterogeneous disease. The aim of this study is to identify the diagnostic and poor prognostic signatures in TNBC by exploring the aberrant DNA methylation and gene expression. Differential expression and methylation analysis of the TNBC and paracancer samples from The Cancer Genome Atlas were performed. Gene set enrichment and protein–protein interaction (PPI) network analysis was used to explore the mechanisms of TNBC. Methylation‐gene expression correlation analysis was performed, and multivariate Cox analysis and receiver operating characteristics analysis were used to further screen the hub genes for TNBC. We identified 1,525 differentially expressed genes and 150 differentially methylated genes between TNBC and paracancer samples. About 96.64% of the methylation sites were located on the CpG island. A total of 17 Gene Ontology biological process terms and 18 signal pathways were significantly enriched. GNG4, GNG11, PENK, MAOA, and AOX1 were identified as the core genes of the PPI network. Methylation‐expression correlations revealed that ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) showed promise as diagnostic and prognostic markers in TNBC. ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) were potential diagnostic and prognostic markers in TNBC.     

Circulating ctDNA methylation quantification of two DNA methyl transferases in papillary thyroid carcinoma

Papillary thyroid cancer (PTC) is the most common type of cancer among thyroid malignancies. Tumor-related methylation of circulating tumor DNA (ctDNA) in plasma could represent tumor specific alterations can be considered as good biomarkers in circulating tumor cells. In this study, we studied the methylation status of seven promoter regions of two DNA methyl Transferases (MGMT and DNMT1) genes as the methylated ctDNA in plasma and tissue samples of patients with PTC and goiter patients as noncancerous controls. Methods: Both ctDNA and tissue genomic DNA of 57 PTC and 45 Goiter samples were isolated. After bisulfite modification, the methylation status was studied by Methylation-Sensitive High Resolution Melting (MS-HRM) assay technique. Four promoter regions of O6-methylguanine-DNA methyltransferase (MGMT) and three promoter regions of DNA methyltransferase 1 (DNMT1) were assessed. Results: From seven candidate promoter regions of two methyltrasferase coding genes, the methylation status of ctDNA within MGMT (a), MGMT (c), MGMT (d), and DNMT1 (b) were meaningfully different between PTC cases and controls. However, the most significant differences were seen in circulating ctDNA MGMT (c) which was hypermethylated in 25 (43.9 %) of patients with PTC vs 2 (4. 4 %) of goiter samples. Between two selected DNA methyl transferase, the methylation of MGMT as the maintenance methyltransferase was significantly higher in PTC cases than goiter controls (P-value < .001). The resulting areas under the receiver operating characteristic (ROC) curve were 0.78 for MGMT (d) for PTC versus goiter samples that can represent the overall ability of MGMT (d) methylation status to discriminate between PTC and goiter patients. Conclusion: Among seven candidate regions of ctDNA the MGMT (c) and MGMT (d) showed higher sensitivity and specificity for PTC as a suitable candidates as biomarkers of PTC.     

Genome‐wide methylomic analyses identify prognostic epigenetic signature in lower grade glioma

Glioma is the most malignant and aggressive type of brain tumour with high heterogeneity and mortality. Although some clinicopathological factors have been identified as prognostic biomarkers, the individual variants and risk stratification in patients with lower grade glioma (LGG) have not been fully elucidated. The primary aim of this study was to identify an efficient DNA methylation combination biomarker for risk stratification and prognosis in LGG. We conducted a retrospective cohort study by analysing whole genome DNA methylation data of 646 patients with LGG from the TCGA and GEO database. Cox proportional hazard analysis was carried out to screen and construct biomarker model that predicted overall survival (OS). The Kaplan‐Meier survival curves and time‐dependent ROC were constructed to prove the efficiency of the signature. Then, another independent cohort was used to further validate the finding. A two‐CpG site DNA methylation signature was identified by multivariate Cox proportional hazard analysis. Further analysis indicated that the signature was an independent survival predictor from other clinical factors and exhibited higher predictive accuracy compared with known biomarkers. This signature was significantly correlated with immune‐checkpoint blockade, immunotherapy‐related signatures and ferroptosis regulator genes. The expression pattern and functional analysis showed that these two genes corresponding with two methylation sites contained in the model were correlated with immune infiltration level, and involved in MAPK and Rap1 signalling pathway. The signature may contribute to improve the risk stratification of patients and provide a more accurate assessment for precision medicine in the clinic.