scholarly journals Methylation and transcriptome analysis reveal lung adenocarcinoma-specific diagnostic biomarkers

2019 ◽  
Vol 17 (1) ◽  
Author(s):  
Rui Li ◽  
Yi-E Yang ◽  
Yun-Hong Yin ◽  
Meng-Yu Zhang ◽  
Hao Li ◽  
...  
Keyword(s):  
Gene Expression ◽  
Dna Methylation ◽  
Lung Adenocarcinoma ◽  
Cox Regression ◽  
R Package ◽  
Functional Enrichment ◽  
Independent Effect ◽  
Cox Regression Analysis ◽  
Kaplan Meier ◽  
Differentially Methylated Genes

Abstract Background DNA methylation can regulate the role of long noncoding RNAs (lncRNAs) in the development of lung adenocarcinoma (LUAD). The present study aimed to identify methylation-driven lncRNAs and mRNAs as biomarkers in the prognosis of LUAD using bioinformatics analysis. Methods Differentially expressed RNAs were obtained using the edge R package from 535 LUAD tissues and 59 adjacent non-LUAD tissues. Differentially methylated genes were obtained using the limma R package from 475 LUAD tissues and 32 adjacent non-LUAD tissues. Methylation-driven mRNA and lncRNA were obtained using the MethylMix R package from 465 LUAD tissues with matched DNA methylation and RNA expression and 32 non-LUAD tissues with DNA methylation. Gene ontology and ConsensusPathDB pathway analysis were performed to identify functional enrichment of methylation-driven mRNAs. Univariate and multivariate Cox regression analyses were performed to identify the independent effect of each variable for predicting the prognosis of LUAD. Kaplan–Meier curve analysis of DNA methylation and gene expression might provide potential prognostic biomarkers for LUAD patients. Results A total of 99 methylation-driven mRNAs and 17 methylation-driven lncRNAs were obtained. Univariate and multivariate Cox regression analysis showed that 6 lncRNAs (FOXE1, HOXB13-AS1_2, VMO1, HIST1H3F, AJ003147.8, ASXL3) were retrieved to construct a predictive model associated with overall survival in LUAD patients. Combined DNA methylation and gene expression survival analysis revealed that 4 lncRNAs (AC023824.1, AF186192.1, LINC01354 and WASIR2) and 8 mRNAs (S1PR1, CCDC181, F2RL1, EFS, KLHDC9, MPV17L, GKN2, ITPRIPL1) might act as independent biomarkers for the prognosis of LUAD. Conclusions Methylation-driven lncRNA and mRNA contribute to the survival of LUAD, and 4 lncRNAs and 8 mRNAs might be potential biomarkers for the prognosis of LUAD.

scholarly journals Identification of Methylation Driven Biomarkers for Diagnosis and Prognosis in Colon Cancer by Integrative Analysis of TCGA, GTEx, and GEO Data

2021 ◽  
Author(s):  
Lichao Cao ◽  
Erfei Chen ◽  
Jin Yang
Keyword(s):  
Gene Expression ◽  
Dna Methylation ◽  
Mrna Expression ◽  
Cox Regression ◽  
Integrative Analysis ◽  
The Cancer Genome Atlas ◽  
Expression Data ◽  
Cox Regression Analysis ◽  
Diagnosis And Prognosis ◽  
Differentially Methylated Genes

Abstract Background: The intention of the present work was to investigate methylation driven biomarkers for diagnosis and prognosis in colorectal cancer (CRC) by integrative analysis of DNA methylation and gene expression data from The Cancer Genome Atlas (TCGA), The Genotype-Tissue Expression (GTEx), and Gene Expression Omnibus (GEO). Methods and Results: The differentially expression genes (DEGs) and differentially methylated genes (DMGs) were screened using mRNA expression and DNA methylation data from TCGA, respectively. And the methylation driven genes (MDGs) of CRC were further identified using MethylMix R package. Subsequently, the MDGs were underwent Random Forest (RF) analyses to establish diagnosis prediction model using mRNA expression data from TCGA and GTEx, which were then validated by GSE39582 from GEO dataset. In addition, prognostic biomarkers that were used to establish the methylation related risk score model was generated by the univariate and multivariate Cox regression analysis. 9 out of 10 DMGs revealed excellent performance as independent diagnostic predictors, including CLDN1, EPHX4, TCN1, ARHGAP20, LY6G6D, FAM150A, KCNJ12, KRT7 and STK33. Furthermore, STK33 and EPHX4 were found to be associated with Overall survival (OS). Conclusions: Our findings suggest that the identified MDGs could be potential biomarkers for diagnosis and prognosis of CRC.

scholarly journals Establishment of a 5 DNA Methylation Site Prognostic Signature Associated with N6-Methyladenosine in Lung Adenocarcinoma

2021 ◽  
Author(s):  
Liqiang Yuan ◽  
Wei Jiang ◽  
Zhanyu Xu ◽  
Kung Deng ◽  
Yu Sun ◽  
...  
Keyword(s):  
Dna Methylation ◽  
Lung Adenocarcinoma ◽  
Cox Regression ◽  
Clinicopathological Parameters ◽  
Main Body ◽  
Methylation Site ◽  
Cox Regression Analysis ◽  
Site Model ◽  
Kaplan Meier ◽  
Cox Analysis

Abstract Background: There is a high incidence of lung adenocarcinoma (LUAD). Even with surgery, targeted therapy and immunotherapy, the survival rate of LUAD patients is still low. N6-methyladenosine (m6A) and DNA methylation markers can help with the diagnosis and treatment of LUAD patients. Therefore, it is necessary to identify a novel m6A-related DNA methylation sites signature to predict the survival of patients with LUAD. Methods: In this study, we screened 15 m6A-related genes and their 217 methylation sites. RNA sequencing data of 15 genes and the clinicopathological parameters of TCGA-LUAD were obtained from the TCGA database (http://cancergenome.nih.gov/). The LUAD-DNA CpG site information was obtained from the Illumina Human Methylation 450 BeadChip (Illumina, San Diego, CA, United States). The methylation sites related to prognosis were screened using univariate COX analysis, and the independent predictors of LUAD patients were identified using multivariate COX analysis of least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Finally, a model with 5 methylation sites as the main body to predict the prognosis of OS in patients with LUAD was obtained. According to the risk grouping of the prediction model, Kaplan-Meier curve and the receiver operating characteristic (ROC) curve were performed in the test and training sets to assess the predicted capacity of the model. In addition, a nomogram constructed by combining the risk score of methylation group and other related clinicopathological factors to verify the reliability of our model.Results: We constructed a m6A-related 5-DNA methylation site model to predict OS in LUAD patients. According to the results of the Kaplan-Meier curve, both the test set and the training set, the high-risk group showed a worse prognosis. The AUCs of the 5 DNA methylation signature at 1, 5 and 10 years in test datasets were 0.730, 0.649 and 0.726, respectively, and 0.679, 0.656 and 0.732 in training datasets. Finally, we constructed a nomogram to further verify the reliability of the model.Conclusion: In this study, we analyzed the methylation sites of m6A-related genes and established a m6A-related 5-DNA methylation site model to predict OS in LUAD patients.

A Six CT83-Related Genes Based Prognostic Signature for Lung Adenocarcinoma

2021 ◽  
Vol 24 ◽  
Author(s):  
Yongmei Wang ◽  
Guimin Zhang ◽  
Ruixian Wang
Keyword(s):  
Regression Analysis ◽  
Lung Adenocarcinoma ◽  
Differential Expression ◽  
Risk Score ◽  
Prognostic Model ◽  
Cox Regression ◽  
Differential Expression Analysis ◽  
Functional Enrichment ◽  
Risk Scores ◽  
Cox Regression Analysis

Background: This study aims to explore the prognostic values of CT83 and CT83-related genes in lung adenocarcinoma (LUAD). Methods: We downloaded the mRNA profiles of 513 LUAD patients (RNA sequencing data) and 246 NSCLC patients (Affymetrix Human Genome U133 Plus 2.0 Array) from TCGA and GEO databases. According to the median expression of CT83, the TCGA samples were divided into high and low expression groups, and differential expression analysis between them was performed. Functional enrichment analysis of differential expression genes (DEGs) was conducted. Univariate Cox regression analysis and LASSO Cox regression analysis were performed to screen the optimal prognostic DEGs. Then we established the prognostic model. A Nomogram model was constructed to predict the overall survival (OS) probability of LUAD patients. Results: CT83 expression was significantly correlated to the prognosis of LUAD patients. A total of 59 DEGs were identified, and a predictive model was constructed based on six optimal CT83-related DEGs, including CPS1, RHOV, TNNT1, FAM83A, IGF2BP1, and GRIN2A, could effectively predict the prognosis of LUAD patients. The nomogram could reliably predict the OS of LUAD patients. Moreover, the six important immune checkpoints (CTLA4, PD1, IDO1, TDO2, LAG3, and TIGIT) were closely correlated with the Risk Score, which was also differentially expressed between the LUAD samples with high and low-Risk Scores, suggesting that the poor prognosis of LUAD patients with high-Risk Score might be due to the immunosuppressive microenvironments. Conclusion: A prognostic model based on six optimal CT83 related genes could effectively predict the prognosis of LUAD patients.

scholarly journals Identification of a Novel Gene Model-Based Homologous Recombination Deficiency Score to Improve Survival Prediction of TNBC.

2021 ◽  
Author(s):  
Wenxiang Zhang ◽  
Bolun Ai ◽  
Xiangyi Kong ◽  
Xiangyu Wang ◽  
Jie Zhai ◽  
...  
Keyword(s):  
Breast Cancer ◽  
Gene Expression ◽  
Homologous Recombination ◽  
Risk Score ◽  
Cox Regression ◽  
Roc Curves ◽  
Cox Regression Analysis ◽  
Homologous Recombination Deficiency ◽  
Ajcc Stage ◽  
Kaplan Meier

Abstract Background Triple-negative breast cancer (TNBC) is a specific histological type of breast cancer with a poor prognosis, early recurrence, which lacks durable chemotherapy responses and effective targeted therapies. We aimed to construct an accurate prognostic risk model based on homologous recombination deficiency (HRD) - gene expression profiles for improving prognosis prediction of TNBC. Methods Triple-negative breast cancer RNA sequencing data and sample clinical information were downloaded from the breast invasive carcinoma (BRCA) cohort in the Cancer Genome Atlas (TCGA) database. Combined with the HRD database, tumor samples were divided into two sets. We screened differentially expressed genes (DEGs) and then identified HRD-related prognostic genes using weighted gene co-expression network analysis (WGCNA) and Cox regression analysis. The least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis were used to identifying key prognostic genes. Risk scores were calculated and compared with HRD score, Kaplan–Meier (KM) survival analysis were used to assess its prognostic power. GSE103091 dataset from GEO (Gene Expression Omnibus) database was used to validate the signature. Univariate and multivariate Cox regression were performed to independently verify the prognosis of the risk score. A nomogram was constructed and revealed by time-dependent ROC curves to guide clinical practice. Results We found that HRD tumor samples (HRD score > = 42) in TNBC patients were associated with poor overall survival (p = 0.027). We identified a total of 147 differential genes including 203 up-regulated and 213 down-regulated genes, among which 29 were prognosis-related genes. Through the LASSO method, 6 key prognostic genes ((MUCL1, IVL, FAM46C, CHI3L1, PRR15L, and CLEC3A) were selected and a 6-gene risk score was constructed. We found risk score was negatively associated with homologous recombination deficiency (HRD) scores (r = -0.22, p = 0.019). Compared with the low-risk group, Kaplan-Meier survival analysis shows that the high-risk group has an obvious poorer prognosis (P < 0.0001). Finally, we integrated the risk score model and clinical factors of TNBC (AJCC-stage, HRD score, T stage, and N stage) to construct a compound nomogram. Time-dependent ROC curves showed the risk score performed better in 1-, 3- and 5-year survival predictions compared with AJCC-stage. Conclusions Based on HRD gene expression data, our six HRD-related gene signature and nomogram could be practical and reliable tools for predicting OS in patients with TNBC.

scholarly journals Identification and Validation of a Potential Prognostic 7-lncRNA Signature for Predicting Survival in Patients with Multiple Myeloma

2020 ◽  
Vol 2020 ◽  
pp. 1-16 ◽  
Author(s):  
Yun Zhong ◽  
Zhe Liu ◽  
Dangchi Li ◽  
Qinyuan Liao ◽  
Jingao Li
Keyword(s):  
Gene Expression ◽  
Risk Score ◽  
Cox Regression ◽  
External Validation ◽  
Survival Rates ◽  
Enrichment Analysis ◽  
Functional Enrichment Analysis ◽  
Functional Enrichment ◽  
Test Set ◽  
Cox Regression Analysis

Background. An increasing number of studies have indicated that the abnormal expression of certain long noncoding RNAs (lncRNAs) is linked to the overall survival (OS) of patients with myeloma. Methods. Gene expression data of myeloma patients were downloaded from the Gene Expression Omnibus (GEO) database (GSE4581 and GSE57317). Cox regression analysis, Kaplan-Meier, and receiver operating characteristic (ROC) analysis were performed to construct and validate the prediction model. Single sample gene set enrichment (ssGSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to predict the function of a specified lncRNA. Results. In this study, a seven-lncRNA signature was identified and used to construct a risk score system for myeloma prognosis. This system was used to stratify patients with different survival rates in the training set into high-risk and low-risk groups. Test set, the entire test set, the external validation set, and the myeloma subtype achieved the authentication of the results. In addition, functional enrichment analysis indicated that 7 prognostic lncRNAs may be involved in the tumorigenesis of myeloma through cancer-related pathways and biological processes. The results of the immune score showed that IF_I was negatively correlated with the risk score. Compared with the published gene signature, the 7-lncRNA model has a higher C-index (above 0.8). Conclusion. In summary, our data provide evidence that seven lncRNAs could be used as independent biomarkers to predict the prognosis of myeloma, which also indicated that these 7 lncRNAs may be involved in the progression of myeloma.

scholarly journals Identification of a Ferroptosis-Related LncRNA Signature as a Novel Prognosis Model for Lung Adenocarcinoma

2021 ◽  
Vol 11 ◽  
Author(s):  
Lu Lu ◽  
Le-Ping Liu ◽  
Qiang-Qiang Zhao ◽  
Rong Gui ◽  
Qin-Yu Zhao
Keyword(s):  
Regression Analysis ◽  
Prediction Model ◽  
Lung Adenocarcinoma ◽  
Tumor Cells ◽  
Risk Prediction ◽  
Prognostic Value ◽  
Cox Regression ◽  
Risk Prediction Model ◽  
Functional Enrichment ◽  
Cox Regression Analysis

Lung adenocarcinoma (LUAD) is a highly heterogeneous malignancy, which makes prognosis prediction of LUAD very challenging. Ferroptosis is an iron-dependent cell death mechanism that is important in the survival of tumor cells. Long non-coding RNAs (lncRNAs) are considered to be key regulators of LUAD development and are involved in ferroptosis of tumor cells, and ferroptosis-related lncRNAs have gradually emerged as new targets for LUAD treatment and prognosis. It is essential to determine the prognostic value of ferroptosis-related lncRNAs in LUAD. In this study, we obtained RNA sequencing (RNA-seq) data and corresponding clinical information of LUAD patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and ferroptosis-related lncRNAs by co-expression analysis. The best predictors associated with LUAD prognosis, including C5orf64, LINC01800, LINC00968, LINC01352, PGM5-AS1, LINC02097, DEPDC1-AS1, WWC2-AS2, SATB2-AS1, LINC00628, LINC01537, LMO7DN, were identified by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis, and the LUAD risk prediction model was successfully constructed. Kaplan-Meier analysis, receiver operating characteristic (ROC) time curve analysis and univariate and multivariate Cox regression analysis and further demonstrated that the model has excellent robustness and predictive ability. Further, based on the risk prediction model, functional enrichment analysis revealed that 12 prognostic indicators involved a variety of cellular functions and signaling pathways, and the immune status was different in the high-risk and low-risk groups. In conclusion, a risk model of 12 ferroptosis related lncRNAs has important prognostic value for LUAD and may be ferroptosis-related therapeutic targets in the clinic.

scholarly journals Construction and Analysis of Survival-Associated Competing Endogenous RNA Network in Lung Adenocarcinoma

2021 ◽  
Vol 2021 ◽  
pp. 1-17
Author(s):  
Lixian Chen ◽  
Zhonglu Ren ◽  
Yongming Cai
Keyword(s):  
Prognostic Factors ◽  
Lung Adenocarcinoma ◽  
Molecular Mechanisms ◽  
Cox Regression ◽  
Enrichment Analysis ◽  
Functional Enrichment Analysis ◽  
Functional Enrichment ◽  
Competing Endogenous Rna ◽  
Cox Regression Analysis ◽  
Cerna Network

Increasing evidence has shown that noncoding RNAs play significant roles in the initiation, progression, and metastasis of tumours via participating in competing endogenous RNA (ceRNA) networks. However, the survival-associated ceRNA in lung adenocarcinoma (LUAD) remains poorly understood. In this study, we aimed to investigate the regulatory mechanisms underlying ceRNA in LUAD to identify novel prognostic factors. mRNA, lncRNA, and miRNA sequencing data obtained from the GDC data portal were utilized to identify differentially expressed (DE) RNAs. Survival-related RNAs were recognized using univariate Kaplan-Meier survival analysis. We performed functional enrichment analysis of survival-related mRNAs using the clusterProfiler package of R and STRING. lncRNA-miRNA and miRNA-mRNA interactions were predicted based on miRcode, Starbase, and miRanda. Subsequently, the survival-associated ceRNA network was constructed for LUAD. Multivariate Cox regression analysis was used to identify prognostic factors. Finally, we acquired 15 DE miRNAs, 49 DE lncRNAs, and 843 DE mRNAs associated with significant overall survival. Functional enrichment analysis indicated that survival-related DE mRNAs were enriched in cell cycle. The survival-associated lncRNA-miRNA-mRNA ceRNA network was constructed using five miRNAs, 49 mRNAs, and 21 lncRNAs. Furthermore, seven hub RNAs (LINC01936, miR-20a-5p, miR-31-5p, TNS1, TGFBR2, SMAD7, and NEDD4L) were identified based on the ceRNA network. LINC01936 and miR-31-5p were found to be significant using the multifactorial Cox regression model. In conclusion, we successfully constructed a survival-related lncRNA-miRNA-mRNA ceRNA regulatory network in LUAD and identified seven hub RNAs, which provide novel insights into the regulatory molecular mechanisms associated with survival of LUAD, and identified two independent prognostic predictors for LUAD.

scholarly journals Six-long Non-coding RNA Signature to Predict the Prognosis of Lung Adenocarcinoma

2020 ◽  
Author(s):  
Yang Wang ◽  
Chengping Hu
Keyword(s):  
Lung Adenocarcinoma ◽  
Cox Regression ◽  
Survival Rates ◽  
Enrichment Analysis ◽  
Functional Enrichment Analysis ◽  
The Cancer Genome Atlas ◽  
Functional Enrichment ◽  
Protein Coding ◽  
Cox Regression Analysis ◽  
Protein Coding Genes

Abstract Background: Long non-coding RNAs (lncRNAs) have been reported to play essential roles in tumorigenesis and cancers prognosis, and they can be a potential cancer prognostic markers. However, in lung adenocarcinoma(LUAD), how lncRNA signatures predict the survival of patients is poorly understood. Our study aims to explore lncRNA signatures and prognostic function in LUAD.Methods: The expression and prognosis data of lncRNAs in LUAD patients was collected from the Cancer Genome Atlas (TCGA) data. All analyses were performed using the R package (version 3.6.2). Metascape, STRING and Cytoscape were used for enrichment analysis and function prediction of the lncRNA co-expressed protein-coding genes.Results: We have collected lncRNA expression data in 466 LUAD tumors, and a six-lncRNA signature(RP11-79H23.3, RP11-309M7.1, CTD-2357A8.3, RP11-108P20.4, U47924.29, LHFPL3-AS2) has been shown to be significantly related to LUAD patients’ overall survival. According to the lncRNA signatures, the high-risk and low-risk groups were divided in LUAD patients with different survival rates. Further multivariable cox regression analysis showed that the prognostic value of this signature was independent of clinical factors. The potential functional roles and hub co-expressed protein-coding genes in the six prognostic lncRNAs are shown in the functional enrichment analysis.Conclusions: These results showed that these six lncRNAs could be independent predicted prognostic biomarkers in LUAD patients.

scholarly journals BACH1: A Potential Predictor of Survival in Early-Stage Lung Adenocarcinoma

2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Jin Zhou ◽  
Zheming Liu ◽  
Huibo Zhang ◽  
Tianyu Lei ◽  
Jiahui Liu ◽  
...  
Keyword(s):  
Lung Adenocarcinoma ◽  
Prognostic Model ◽  
Cox Regression ◽  
Early Stage ◽  
External Validation ◽  
Enrichment Analysis ◽  
Gene Set Enrichment Analysis ◽  
The Cancer Genome Atlas ◽  
Functional Enrichment ◽  
Cox Regression Analysis

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.

scholarly journals Integrated DNA Methylation and Gene Expression Analysis Identified S100A8 and S100A9 in the Pathogenesis of Obesity

2021 ◽  
Vol 8 ◽  
Author(s):  
Ningyuan Chen ◽  
Liu Miao ◽  
Wei Lin ◽  
Donghua Zou ◽  
Ling Huang ◽  
...  
Keyword(s):  
Gene Expression ◽  
Dna Methylation ◽  
Mrna Expression ◽  
Regulatory Networks ◽  
Enrichment Analysis ◽  
Microarray Dataset ◽  
Functional Enrichment Analysis ◽  
Normal Subjects ◽  
Functional Enrichment ◽  
Differentially Methylated Genes

Background: To explore the association of DNA methylation and gene expression in the pathology of obesity.Methods: (1) Genomic DNA methylation and mRNA expression profile of visceral adipose tissue (VAT) were performed in a comprehensive database of gene expression in obese and normal subjects. (2) Functional enrichment analysis and construction of differential methylation gene regulatory networks were performed. (3) Validation of the two different methylation sites and corresponding gene expression was done in a separate microarray dataset. (4) Correlation analysis was performed on DNA methylation and mRNA expression data.Results: A total of 77 differentially expressed mRNAs matched with differentially methylated genes. Analysis revealed two different methylation sites corresponding to two unique genes—s100a8-cg09174555 and s100a9-cg03165378. Through the verification test of two interesting different expression positions [differentially methylated positions (DMPs)] and their corresponding gene expression, we found that methylation in these genes was negatively correlated to gene expression in the obesity group. Higher S100A8 and S100A9 expressions in obese subjects were validated in a separate microarray dataset.Conclusion: This study confirmed the relationship between DNA methylation and gene expression and emphasized the important role of S100A8 and S100A9 in the pathogenesis of obesity.

Sign in / Sign up

Export Citation Format

Share Document