scholarly journals Identification of Prognostic Risk Model Based on DNA Methylation-Driven Genes in Esophageal Adenocarcinoma

2021 ◽  
Vol 2021 ◽  
pp. 1-20
Author(s):  
Yuhua Chen ◽  
Jinjie Wang ◽  
Hao Zhou ◽  
Zhanghao Huang ◽  
Li Qian ◽  
...  

Background. DNA methylation is an important part of epigenetic modification, and its abnormality is closely related to esophageal adenocarcinoma (EAC). This study was aimed at using bioinformatics analysis to identify methylation-driven genes (MDGs) in EAC patients and establish a risk model as a biological indicator of EAC prognosis. Method. Downloaded EAC DNA methylation, transcriptome, and related clinical data from TCGA database. MethylMix was used to identify MDGs. R package clusterProfiler and the ConsensusPathDB online database were used to analyze the rich functions and pathways of these MDGs. The prognostic risk model was established by univariate Cox regression, Lasso regression, and multivariate Cox regression analysis. Finally each MDG in the model were carried out through the survival R package. Results. A total of 273 MDGs were identified, which were enriched in transcriptional regulation and embryonic organ morphogenesis. Cox regression analysis established a risk model consisting of GPBAR1, OLFM4, FOXI2, and CASP10. In addition, further survival analysis revealed that OLFM4 and its two related sites were significantly related to the EAC patients’ survival. Conclusion. In summary, this study used bioinformatics methods to identify EAC MDGs and established a reliable risk prognosis model. It provided potential biomarkers for the early treatment and prognosis evaluation of EAC.

2021 ◽  
Author(s):  
Yuan Li ◽  
Hao Huang ◽  
Jun Feng ◽  
Yulan Zhu ◽  
Tianwei Jiang ◽  
...  

Abstract BackgroundAlthough some advanced colorectal cancer (CRC) patients could select immunotherapy, but still most microsatellite stability (MSS) CRC patients did not respond. Our present study aims to set up a novel system for prognostic prediction and immunotherapeutic responsiveness for MSS CRC patients.MethodsUnivariable Cox regression survival analysis and least absolute shrinkage and selector operation (LASSO) regression analysis were performed to identify prognostic genes and establish immune risk signatures. Multivariate Cox regression analysis was performed to verify whether these clinical features could predict prognosis. R package was used to analyze the relationship between the immune-related risk model and these immune cells, effector molecules, and immune checkpoints.ResultsWe constructed an immune-related signature and verified its predictive capability. Immune-related signature included 12 differentially expressed IRGs (12 DE IR MSSGs), including CXCL1, CD36, FABP4, MS4A2, NRG1, VGF, GRP, HDC, XCL1, NGF, MAGEA1, and IL13. The signature consisting of 12 DE IR MSSGs was an independent and effective prognostic factor for the overall survival of CRC patients. In addition, the signature consisting of 12 DE IR MSSGs reflected the infiltration characteristics of different immunocytes in tumor immune microenvironment. The signature consisting of 12 DE IR MSSGs also had a significant correlation with immune checkpoint molecules.


2021 ◽  
Author(s):  
Cheng Yan ◽  
Qingling Liu ◽  
Ruoling Jia

Abstract Background: Autophagy plays an important role in triple negative breast cancer (TNBC). However, the prognostic value of autophagy-related genes (ARGs) in TNBC remains unknown. In this study, we established a survival model to evaluate the prognosis of TNBC patients using ARGs signature.Methods: A total of 222 autophagy-related genes were downloaded from The Human Autophagy Database. The RNA-sequencing data and corresponding clinical data of TNBC were obtained from the TCGA database. Differential gene expression of ARGs (DE-ARGs) between normal samples and TNBC samples was determined by the EdgeR software package. Then, univariate Cox, Lasso, and multivariate Cox regression analyses were performed. According to the Lasso regression results based on univariate Cox, we identified a prognostic signature for overall-survival (OS), which was further validated by using GEO cohort. We also found an independent prognostic marker that can predict the clinicopathological features of TNBC. Furthermore, a nomogram was drawn to predict the survival probability of TNBC patients, which could help in clinical decision for TNBC treatment. Finally, we validated the requirement of a ARG in our model for TNBC cell survival and metastasis.Results: There are 43 differentially expressed ARGs (DE-ARGs) were identified between normal and tumor samples. A risk model for OS using CDKN1A, CTSD, CTSL, EIF4EBP1, TMEM74 and VAMP3 by Lasso regression analysis was established based on univariate Cox regression analysis. Overall survival of TNBC patients was significantly shorter in the high-risk group than in the low-risk group for both the training and validation cohorts. Using the Kaplan-Meier curves and ROC curves, we demonstrated the accuracy of the prognostic model. Multivariate Cox regression analysis was used to verify risk score as independent predictor. Then a nomogram was proposed to predict 1-, 3-, and 5-year survival for TNBC patients. The calibration curves showed great accuracy of the model for survival prediction. Finally, we found that depletion of EIF4EBP1, one of ARGs in our model, significantly reduced cell proliferation and metastasis of TNBC cells. Conclusion: An autophagy-related prognosis model in TNBCs was constructed using ARGs signature containing CDKN1A, CTSD, CTSL, EIF4EBP1, TMEM74 and VAMP3. It could serve as an independent prognostic biomarker in TNBC.


2021 ◽  
Author(s):  
Xin-Yu Li ◽  
Lei Hou ◽  
Lu-yu Zhang ◽  
Xue-yuan Li ◽  
xi-tao Yang

Abstract Aim: A glioblastoma (GBM) prognostic model was developed with GBM -related alternative splicing (AS) data and prognostic markers were identified. Methods: AS data and clinical data of GBM patients were retrieved from The Cancer Genome Atlas (TCGA) SpliceSeq database and TCGA database, respectively. The data from these two databases were intersected to screen the prognosis-associated AS events, which was subsequently examined in Univariate Cox regression models. To avoid model overfitting, LASSO regression analysis was conducted. On the basis of these AS events, we established a prognostic model of GBM with the use of multivariate Cox regression analysis. On the strength of this model, the patients were assigned into high-risk and low-risk groups with a median risk score as the threshold. Kaplan-Meier survival, receiver operating characteristic (ROC), and calibration curves were applied to evaluate the performance of this model. Finally, combined with the risk model and clinicopathological characteristics, Cox regression analysis was utilized to identify the independent prognostic markers of GBM, and a nomogram was constructed. Results: The AS and clinical data of 169 GBM patients from the TCGA SpliceSeq and TCGA databases were collected. Univariate Cox regression analysis identified 1000 prognosis-related AS events in GBM, and then Lasso regression analysis identified 16 AS events. A GBM prognostic risk model was constructed based on AS events of 7 genes (FAM86B1, ZNF302, C19orf57, RPL39L, CBLL1, RWDD1, IGF2BP2). Through this model, we found lower overall survival (OS) rates of the high-risk population versus the low-risk population (p < 0.05). ROC and calibration curve analyses demonstrated the good ability of this model to predict the OS of GBM patients. Cox regression analysis suggested risk score as an independent prognostic factor for GBM. We also found that IGF2BP2 is associated with patient prognosis and have a strong relationship with immunotherapy response. Conclusion: The prognostic model based on AS events can significantly distinguish the survival rate of high-risk and low-risk GBM patients and IGF2BP2 were identified as a novel prognostic biomarker and immunotherapeutic target.


2021 ◽  
Author(s):  
Rui Feng ◽  
Jian Li ◽  
Weiling Xuan ◽  
Hanbo Liu ◽  
Dexin Cheng ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) is a prevalent primary liver cancer and the main cause of cancer mortality. Its high complexity and dismal prognosis bring dramatic difficulty to treatment. Due to the disclosed dual functions of autophagy in cancer development, understanding autophagy-related genes devotes into seeking novel biomarkers for HCC. Methods Differential expression of genes in normal and tumor groups was analyzed to acquire autophagy-related genes in HCC. GO and KEGG pathway analyses were conducted on these genes. Genes were then screened by univariate regression analysis. The screened genes were subjected to multivariate Cox regression analysis to build a prognostic model. The model was validated by ICGC validation set. Results Altogether, 42 autophagy-related differential genes were screened by differential expression analysis. Enrichment analysis showed that they were mainly enriched in pathways including regulation of autophagy and cell apoptosis. Genes were screened by univariate analysis and multivariate Cox regression analysis to build a prognostic model. The model was constituted by 6 feature genes: EIF2S1, BIRC5, SQSTM1, ATG7, HDAC1, FKBP1A. Validation confirmed the accuracy and independence of this model in predicting HCC patient’s prognosis. Conclusion A total of 6 feature genes were identified to build a prognostic risk model. This model is conducive to investigating interplay between autophagy-related genes and HCC prognosis.


2020 ◽  
Author(s):  
Xiang Zhou ◽  
Keying Zhang ◽  
Fa Yang ◽  
Chao Xu ◽  
Jianhua Jiao ◽  
...  

Abstract Background: Hepatocellular carcinoma (HCC) is a disease with higher morbidity, mortality, and poor prognosis in the whole world. Understanding the crosslink between HCC and the immune system is essential for people to uncover a few potential and valuable therapeutic strategies. This study aimed to reveal the correlation between HCC and immune-related genes and establish a clinical evaluation model. Methods: We had analyzed the clinical information consisted of 373 HCC and 49 normal samples from the cancer genome atlas (TCGA). The differentially expressed genes (DEGs) were selected by the Wilcoxon test and the immune-related differentially expressed genes (IRDEGs) in DEGs were identified by matching DEGs with immune-related genes downloaded from the ImmPort database. Furthermore, the univariate Cox regression analysis and multivariate Cox regression analysis were performed to construct a prognostic risk model. Then, twenty-two types of tumor immune-infiltrating cells (TIICs) were downloaded from Tumor Immune Estimation Resource (TIMER) and were used to construct the correlational graphs between the TIICs and risk score by the CIBERSORT. Subsequently, the transcription factors (TFs) were gained in the Cistrome website and the differentially expressed TFs (DETFs) were achieved. Finally, the KEGG pathway analysis and GO analysis were performed to further understand the molecular mechanisms between DETFs and PDIRGs.Results: In our study, 5839 DEGs, 326 IRDEGs, and 31 prognosis-related IRDEGs (PIRDEGs) were identified. And 8 optimal PIRDEGs were employed to construct a prognostic risk model by multivariate Cox regression analysis. The correlation between risk genes and clinical characterizations and TIICs has verified that the prognostic model was effective in predicting the prognosis of HCC patients. Finally, several important immune-related pathways and molecular functions of the eight PIRDEGs were significantly enriched and there was a distinct association between the risk IRDEGs and TFs. Conclusion: The prognostic risk model showed a more valuable predicting role for HCC patients, and produced many novel therapeutic targets and strategies for HCC.


PLoS ONE ◽  
2021 ◽  
Vol 16 (12) ◽  
pp. e0260876
Author(s):  
Jun Yang ◽  
Jiaying Zhou ◽  
Cuili Li ◽  
Shaohua Wang

Background Neuroblastoma (NB) is the most common solid tumor in children. NB treatment has made significant progress; however, given the high degree of heterogeneity, basic research findings and their clinical application to NB still face challenges. Herein, we identify novel prognostic models for NB. Methods We obtained RNA expression data of NB and normal nervous tissue from TARGET and GTEx databases and determined the differential expression patterns of RNA binding protein (RBP) genes between normal and cancerous tissues. Lasso regression and Cox regression analyses identified the five most important differentially expressed genes and were used to construct a new prognostic model. The function and prognostic value of these RBPs were systematically studied and the predictive accuracy verified in an independent dataset. Results In total, 348 differentially expressed RBPs were identified. Of these, 166 were up-regulated and 182 down-regulated RBPs. Two hubs RBPs (CPEB3 and CTU1) were identified as prognostic-related genes and were chosen to build the prognostic risk score models. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analysis using proportional hazards regression model. A five gene prognostic model: Risk score = (-0.60901*expCPEB3)+(0.851637*expCTU1) was built. Based on this model, the overall survival of patients in the high-risk subgroup was lower (P = 2.152e-04). The area under the curve (AUC) of the receiver-operator characteristic curve of the prognostic model was 0.720 in the TARGET cohort. There were significant differences in the survival rate of patients in the high and low-risk subgroups in the validation data set GSE85047 (P = 0.1237e-08), with the AUC 0.730. The risk model was also regarded as an independent predictor of prognosis (HR = 1.535, 95% CI = 1.368–1.722, P = 2.69E-13). Conclusions This study identified a potential risk model for prognosis in NB using Cox regression analysis. RNA binding proteins (CPEB3 and CTU1) can be used as molecular markers of NB.


2021 ◽  
Author(s):  
Boxuan Liu ◽  
Yun Zhao ◽  
Shuanying Yang

Abstract Background: Lung adenocarcinoma is the most occurred pathological type among non-small cell lung cancer. Although huge progress has been made in terms of early diagnosis, precision treatment in recent years, the overall 5-year survival rate of a patient remains low. In our study, we try to construct an autophagy-related lncRNA prognostic signature that may guide clinical practice.Methods: The mRNA and lncRNA expression matrix of lung adenocarcinoma patients were retrieved from TCGA database. Next, we constructed a co-expression network of lncRNAs and autophagy-related genes. Lasso regression and multivariate Cox regression were then applied to establish a prognostic risk model. Subsequently, a risk score was generated to differentiate high and low risk group and a ROC curve and Nomogram to visualize the predictive ability of current signature. Finally, gene ontology and pathway enrichment analysis were executed via GSEA.Results: A total of 1,703 autophagy-related lncRNAs were screened and five autophagy-related lncRNAs (LINC01137, AL691432.2, LINC01116, AL606489.1 and HLA-DQB1-AS1) were finally included in our signature. Judging from univariate(HR=1.075, 95% CI: 1.046–1.104) and multivariate(HR =1.088, 95%CI = 1.057 − 1.120) Cox regression analysis, the risk score is an independent factor for LUAD patients. Further, the AUC value based on the risk score for 1-year, 3-year, 5-year, was 0.735, 0.672 and 0.662 respectively. Finally, the lncRNAs included in our signature were primarily enriched in autophagy process, metabolism, p53 pathway and JAK/STAT pathway. Conclusions: Overall, our study indicated that the prognostic model we generated had certain predictability for LUAD patients’ prognosis.


2021 ◽  
Vol 12 ◽  
Author(s):  
Xiaoping Li ◽  
Jishang Chen ◽  
Qihe Yu ◽  
Hui Huang ◽  
Zhuangsheng Liu ◽  
...  

Background: A surge in newly diagnosed breast cancer has overwhelmed the public health system worldwide. Joint effort had beed made to discover the genetic mechanism of these disease globally. Accumulated research has revealed autophagy may act as a vital part in the pathogenesis of breast cancer.Objective: Aim to construct a prognostic model based on autophagy-related lncRNAs and investigate their potential mechanisms in breast cancer.Methods: The transcriptome data and clinical information of patients with breast cancer were obtained from The Cancer Genome Atlas (TCGA) database. Autophagy-related genes were obtained from the Human Autophagy Database (HADb). Long non-coding RNAs (lncRNAs) related to autophagy were acquired through the Pearson correlation analysis. Univariate Cox regression analysis as well as the least absolute shrinkage and selection operator (LASSO) regression analysis were used to identify autophagy-related lncRNAs with prognostic value. We constructed a risk scoring model to assess the prognostic significance of the autophagy-related lncRNAs signatures. The nomogram was then established based on the risk score and clinical indicators. Through the calibration curve, the concordance index (C-index) and receiver operating characteristic (ROC) curve analysis were evaluated to obtain the model's predictive performance. Subgroup analysis was performed to evaluate the differential ability of the model. Subsequently, gene set enrichment analysis was conducted to investigate the potential functions of these lncRNAs.Results: We attained 1,164 breast cancer samples from the TCGA database and 231 autophagy-related genes from the HAD database. Through correlation analysis, 179 autophagy-related lncRNAs were finally identified. Univariate Cox regression analysis and LASSO regression analysis further screened 18 prognosis-associated lncRNAs. The risk scoring model was constructed to divide patients into high-risk and low-risk groups. It was found that the low-risk group had better overall survival (OS) than those of the high-risk group. Then, the nomogram model including age, tumor stage, TNM stage and risk score was established. The evaluation index (C-index: 0.78, 3-year OS AUC: 0.813 and 5-year OS AUC: 0.785) showed that the nomogram had excellent predictive power. Subgroup analysis showed there were difference in OS between high-risk and low-risk patients in different subgroups (stage I-II, ER positive, Her-2 negative and non-TNBC subgroups; all P &lt; 0.05). According to the results of gene set enrichment analysis, these lncRNAs were involved in the regulation of multicellular organismal macromolecule metabolic process in multicellular organisms, nucleotide excision repair, oxidative phosphorylation, and TGF-β signaling pathway.Conclusions: We identified 18 autophagy-related lncRNAs with prognostic value in breast cancer, which may regulate tumor growth and progression in multiple ways.


2022 ◽  
Author(s):  
Yuying Tan ◽  
Liqing Lu ◽  
Xujun Liang ◽  
Yongheng Chen

Abstract Background: Colon adenocarcinoma (COAD) is one of the most common malignant tumors and diagnosed at an advanced stage with poor prognosis in the world. Pyroptosis is involved in the initiation and progression of tumors. This research focused on constructing a pyroptosis-related ceRNA network to generate a reliable risk model for risk prediction and immune infiltration analysis of COAD.Methods: Transcriptome data, miRNA-sequencing data and clinical information were downloaded from the TCGA database. Firstly, differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and lncRNAs (DElncRNAs) were identified to construct a pyroptosis-related ceRNA network. Secondly, a pyroptosis-related lncRNA risk model was developed applying univariate Cox regression analysis and least absolute shrinkage and selection operator method (LASSO) regression analysis. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses were utilized to functionally annotate RNAs contained in the ceRNA network. In addition, Kaplan-Meier analysis, receiver operating characteristic (ROC) curves, univariate and multivariate Cox regression, and nomogram were applied to validate this risk model. Finally, the relationship of this risk model with immune cells and immune checkpoint blockade (ICB) related genes were analyzed.Results: Totally 5373 DEmRNAs, 1159 DElncRNAs and 355 DEmiRNAs were identified. A pyroptosis-related ceRNA regulatory network containing 132 lncRNAs, 7miRNAs and 5 mRNAs was constructed and a ceRNA-based pyroptosis-related risk model including 11 lncRNAs was built. Tumor tissues were classified into high- and low- risk groups according to the median risk score. Kaplan-Meier analysis showed that the high-risk group had a shorter survival time; ROC analysis, independent prognostic analysis and nomogram further indicated the risk model was a significant independent prognostic factor which had excellent ability to predict patients’ risk. Moreover, immune infiltration analysis indicated that the risk model was related to immune infiltration cells (i.e., B cells naïve, T cells follicular helper, Macrophages M1, etc.) and ICB-related genes (i.e., PD-1, CTLA4, HAVCR2, etc).Conclusions: This pyroptosis-related lncRNA risk model possessed good prognostic value and the ability to predict the outcome of ICB immunotherapy in COAD.


2021 ◽  
Author(s):  
Liusheng Wu ◽  
Xiaoqiang Li ◽  
Jixian Liu ◽  
Da Wu ◽  
Dingwang Wu ◽  
...  

Abstract Objective: Autophagy-related LncRNA genes play a vital role in the development of esophageal adenocarcinoma.Our study try to construct a prognostic model of autophagy-related LncRNA esophageal adenocarcinoma, and use this model to calculate patients with esophageal adenocarcinoma. The survival risk value of esophageal adenocarcinoma can be used to evaluate its survival prognosis. At the same time, to explore the sites of potential targeted therapy genes to provide valuable guidance for the clinical diagnosis and treatment of esophageal adenocarcinoma.Methods: Our study have downloaded 261 samples of LncRNA-related transcription and clinical data of 87 patients with esophageal adenocarcinoma from the TCGA database, and 307 autophagy-related gene data from www.autuphagy.com. We applied R software (Version 4.0.2) for data analysis, merged the transcriptome LncRNA genes, autophagy-related genes and clinical data, and screened autophagy LncRNA genes related to the prognosis of esophageal adenocarcinoma. We also performed KEGG and GO enrichment analysis and GSEA enrichment analysis in these LncRNA genes to analysis the risk characteristics and bioinformatics functions of signal transduction pathways. Univariate and multivariate Cox regression analysis were used to determine the correlation between autophagy-related LncRNA and independent risk factors. The establishment of ROC curve facilitates the evaluation of the feasibility of predicting prognostic models, and further studies the correlation between autophagy-related LncRNA and the clinical characteristics of patients with esophageal adenocarcinoma. Finally, we also used survival analysis, risk analysis and independent prognostic analysis to verify the prognosis model of esophageal adenocarcinoma.Results: We screened and identified 22 autophagic LncRNA genes that are highly correlated with the overall survival (OS) of patients with esophageal adenocarcinoma. The area under the ROC curve(AUC=0.941)and the calibration curve have a good lineup, which has statistical analysis value. In addition, univariate and multivariate Cox regression analysis showed that the autophagy LncRNA feature of this esophageal adenocarcinoma is an independent predictor of esophageal adenocarcinoma.Conclusion: These LncRNA screened and identified may participate in the regulation of cellular autophagy pathways, and at the same time affect the tumor development and prognosis of patients with esophageal adenocarcinoma. These results indicate that risk signature and nomogram are important indicators related to the prognosis of patients with esophageal adenocarcinoma.


Sign in / Sign up

Export Citation Format

Share Document