Construction of the optimization prognostic model based on differentially expressed immune genes of lung adenocarcinoma.
Ontology highlight
ABSTRACT: Background
Lung adenocarcinoma (LUAD) is the most common pathology subtype of lung cancer. In recent years, immunotherapy, targeted therapy and chemotherapeutics conferred a certain curative effects. However, the effect and prognosis of LUAD patients are different, and the efficacy of existing LUAD risk prediction models is unsatisfactory.Methods
The Cancer Genome Atlas (TCGA) LUAD dataset was downloaded. The differentially expressed immune genes (DEIGs) were analyzed with edgeR and DESeq2. The prognostic DEIGs were identified by COX regression. Protein-protein interaction (PPI) network was inferred by STRING using prognostic DEIGs with p valueResultsIn total,1654 DEIGs were identified, of which 436 genes were prognostic. Gene functional enrichment analysis indicated that the DEIGs were involved in inflammatory pathways. We constructed 4 models using DEIGs. Finally, model 4, which was constructed using the 436 DEIGs performed the best in prognostic predictions, the receiver operating characteristic curve (ROC) was 0.824 for 3?years, 0.838 for 5?years, 0.834 for 10?years. High levels of FERMT2, FKBP3 and low levels of SMAD9, GATA2, ITIH4 expression are related to the poor overall survival in LUAD (p?ConclusionsIn our study, we built an optimal prognostic signature for LUAD using DEIGs and verified the expression of selected genes in LUAD. Our result suggests immune signature can be harnessed to obtain prognostic insights.
SUBMITTER: Zhai Y
PROVIDER: S-EPMC7923649 | biostudies-literature | 2021 Mar
REPOSITORIES: biostudies-literature
Publications
Construction of the optimization prognostic model based on differentially expressed immune genes of lung adenocarcinoma.
Zhai Yang Y Zhao Bin B Wang Yuzhen Y Li Lina L Li Jingjin J Li Xu X Chang Linhan L Chen Qian Q Liao Zijun Z
BMC cancer 20210301 1
<h4>Background</h4>Lung adenocarcinoma (LUAD) is the most common pathology subtype of lung cancer. In recent years, immunotherapy, targeted therapy and chemotherapeutics conferred a certain curative effects. However, the effect and prognosis of LUAD patients are different, and the efficacy of existing LUAD risk prediction models is unsatisfactory.<h4>Methods</h4>The Cancer Genome Atlas (TCGA) LUAD dataset was downloaded. The differentially expressed immune genes (DEIGs) were analyzed with edgeR ...[more]
Similar Datasets
Project description:BackgroundLong noncoding RNAs (lncRNAs) could function as competitive endogenous RNAs (ceRNAs) to competitively adsorb microRNAs (miRNAs), thereby regulating the expression of their target protein-coding mRNAs. In this study, we aim to identify more effective diagnostic and prognostic markers for lung adenocarcinoma (LUAD).MethodsWe obtained differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) for LUAD by using The Cancer Genomes Atlas (TCGA) portal. Weighted gene coexpression network analysis (WGCNA) was performed to unveil core gene modules associated with LUAD. The Cox proportional hazards model was performed to determine the prognostic significance of DElncRNAs. The diagnostic and prognostic significance of DElncRNAs was further verified based on the receiver operating characteristic curve (ROC). Cytoscape was used to construct the ceRNA networks comprising the lncRNAs-miRNAs-mRNAs axis based on the correlation obtained from the miRcode, miRDB, and TargetScan.ResultsCompared with normal lung tissues, 2355 DElncRNAs, 820 DEmiRNAs, and 17289 DEmRNAs were identified in LUAD tissues. We generated 8 WGCNA core modules in the lncRNAs coexpression network, 5 modules in the miRNAs, and 12 modules in the mRNAs coexpression network, respectively. One lncRNA module (blue) consisting of 441 lncRNAs, two miRNA modules (blue and turquoise) containing 563 miRNAs, and one mRNA module (turquoise), which consisted of 15162 mRNAs, were mostly significantly related to LUAD status. Furthermore, 67 DEmRNAs were found to be tumor-associated as well as the target genes of the DElncRNAs-DEmiRNAs axis. Survival analyses showed that 6 lncRNAs (LINC01447, WWC2-AS2, OGFRP1, LINC00942, LINC01168, and AC005863.1) were significantly correlated with the prognosis of LUAD patients. Ultimately, the potential ceRNA networks including 6 DElncRNAs, 4 DEmiRNAs, and 22 DEmRNAs were constructed.ConclusionOur study indicated that 6 DElncRNAs had the possibilities as diagnostic and prognostic biomarkers for LUAD. The lncRNA-mediated ceRNA networks might provide novel insights into the molecular mechanisms of LUAD progression.
| S-EPMC8723861 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) originates mainly from the mucous epithelium and glandular epithelium of the bronchi. It is the most common pathologic subtype of non-small cell lung cancer (NSCLC). At present, there is still a lack of clear criteria to predict the efficacy of immunotherapy. The 5-year survival rate for LUAD patients remains low.MethodsAll data were downloaded from The Cancer Genome Atlas (TCGA) database. We used Gene Set Enrichment Analysis (GSEA) database to obtain immune-related mRNAs. Immune-related lncRNAs were acquired by using the correlation test of the immune-related genes with R version 3.6.3 (Pearson correlation coefficient cor = 0.5, P < 0.05). The TCGA-LUAD dataset was divided into the testing set and the training set randomly. Based on the training set to perform univariate and multivariate Cox regression analyses, we screened prognostic immune-related lncRNAs and given a risk score to each sample. Samples were divided into the high-risk group and the low-risk group according to the median risk score. By the combination of Kaplan-Meier (KM) survival curve, the receiver operating characteristic (ROC) (AUC) curve, the independent risk factor analysis, and the clinical data of the samples, we assessed the accuracy of the risk model. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the differentially expressed mRNAs between the high-risk group and the low-risk group. The differentially expressed genes related to immune response between two risk groups were analyzed to evaluate the role of the model in predicting the efficacy and effects of immunotherapy. In order to explain the internal mechanism of the risk model in predicting the efficacy of immunotherapy, we analyzed the differentially expressed genes related to epithelial-mesenchymal transition (EMT) between two risk groups. We extracted RNA from normal bronchial epithelial cell and LUAD cells and verified the expression level of lncRNAs in the risk model by a quantitative real-time polymerase chain reaction (qRT-PCR) test. We compared our risk model with other published prognostic signatures with data from an independent cohort. We transfected LUAD cell with siRNA-LINC0253. Western blot analysis was performed to observed change of EMT-related marker in protein level.ResultsThrough univariate Cox regression analysis, 24 immune-related lncRNAs were found to be strongly associated with the survival of the TCGA-LUAD dataset. Utilizing multivariate Cox regression analysis, 10 lncRNAs were selected to establish the risk model. The K-M survival curves and the ROC (AUC) curves proved that the risk model has a fine predictive effect. The GO enrichment analysis indicated that the effect of the differentially expressed genes between high-risk and low-risk groups is mainly involved in immune response and intercellular interaction. The KEGG enrichment analysis indicated that the differentially expressed genes between high-risk and low-risk groups are mainly involved in endocytosis and the MAPK signaling pathway. The expression of genes related to the efficacy of immunotherapy was significantly different between the two groups. A qRT-PCR test verified the expression level of lncRNAs in LUAD cells in the risk model. The AUC of ROC of 5 years in the independent validation dataset showed that this model had superior accuracy. Western blot analysis verified the change of EMT-related marker in protein level.ConclusionThe immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.
| S-EPMC8044985 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD), the most common non-small-cell lung cancer, is characterized by a dense lymphocytic infiltrate, which indicates that the immune system plays an active role in the development and growth of this cancer. However, no investigations to date have proposed robust models for predicting survival outcome for patients with LUAD in terms of tumour immunology. A total of 761 LUAD patients were included in this study, in which the database of The Cancer Genome Atlas (TCGA) was utilized for discovery, and the Gene Expression Omnibus (GEO) database was utilized for validation. Bioinformatics analysis and R language tools were utilized to construct an immune prognostic model and annotate biological functions. Lung adenocarcinoma showed a weakened immune phenotype compared with adjacent normal tissues. Immune-related gene sets were profiled, an immune prognostic model based on 2 immune genes (ANLN and F2) was developed with the TCGA database to distinguish cases as having a low or high risk of unfavourable prognosis, and the model was verified with the GEO database. The model was prognostically significant in stratified cohorts, including stage I-II, stage III-IV and epidermal growth factor receptor (EGFR) mutant subsets, and was considered to be an independent prognostic factor for LUAD. Furthermore, the low- and high-risk groups showed marked differences in tumour-infiltrating leucocytes, tumour mutation burden, aneuploidy and PD-L1 expression. In conclusion, an immune prognostic model was proposed for LUAD that is capable of independently identifying patients at high risk for poor survival, suggesting a relationship between local immune status and prognosis.
| S-EPMC6991688 | biostudies-literature
Project description:In many cancers, significantly down- or upregulated genes are found within chromosomal regions with DNA copy number alteration opposite to the expression changes. Generally, this paradox has been overlooked as noise, but can potentially be a consequence of interference of epigenetic regulatory mechanisms, including microRNA-mediated control of mRNA levels. To explore potential associations between microRNAs and paradoxes in non-small-cell lung cancer (NSCLC) we curated and analyzed lung adenocarcinoma (LUAD) data, comprising gene expressions, copy number aberrations (CNAs) and microRNA expressions. We integrated data from 1,062 tumor samples and 241 normal lung samples, including newly-generated array comparative genomic hybridization (aCGH) data from 63 LUAD samples. We identified 85 "paradoxical" genes whose differential expression consistently contrasted with aberrations of their copy numbers. Paradoxical status of 70 out of 85 genes was validated on sample-wise basis using The Cancer Genome Atlas (TCGA) LUAD data. Of these, 41 genes are prognostic and form a clinically relevant signature, which we validated on three independent datasets. By meta-analysis of results from 9 LUAD microRNA expression studies we identified 24 consistently-deregulated microRNAs. Using TCGA-LUAD data we showed that deregulation of 19 of these microRNAs explains differential expression of the paradoxical genes. Our results show that deregulation of paradoxical genes is crucial in LUAD and their expression pattern is maintained epigenetically, defying gene copy number status.
| S-EPMC5823624 | biostudies-literature
Project description:Background: Lung adenocarcinoma is one of the most common malignant tumors of the respiratory system, ranking first in morbidity and mortality among all cancers. This study aims to establish a ferroptosis-related gene-based prognostic model to investigate the potential prognosis of lung adenocarcinoma. Methods: We obtained gene expression data with matching clinical data of lung adenocarcinoma from the The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The ferroptosis-related genes (FRGs) were downloaded from three subgroups in the ferroptosis database. Using gene expression differential analysis, univariate Cox regression, and LASSO regression analysis, seven FRGs with prognostic significance were identified. The result of multivariate Cox analysis was utilized to calculate regression coefficients and establish a risk-score formula that divided patients with lung adenocarcinoma into high-risk and low-risk groups. The TCGA results were validated using GEO data sets. Then we observed that patients divided in the low-risk group lived longer than the overall survival (OS) of the other. Then we developed a novel nomogram including age, gender, clinical stage, TNM stage, and risk score. Results: The areas under the curves (AUCs) for 3- and 5-years OS predicted by the model were 0.823 and 0.852, respectively. Calibration plots and decision curve analysis also confirmed the excellent predictive performance of the model. Subsequently, gene function enrichment analysis revealed that the identified FRGs are important in DNA replication, cell cycle regulation, cell adhesion, chromosomal mutation, oxidative phosphorylation, P53 signaling pathway, and proteasome processes. Conclusions: Our results verified the prognostic significance of FRGs in patients with lung adenocarcinoma, which may regulate tumor progression in a variety of pathways.
| S-EPMC8493116 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD), which is the most common type of non-small cell lung cancer (NSCLC), is one of the most aggressive and fatal tumors. Therefore, the identification of key biomarkers affecting prognosis is important to improving the prognosis of patients with LUAD. Cell membranes have long been understood; however, few studies have focused on the role of membrane tension in LUAD. The present study aimed to construct a prognostic model associated with membrane-tension-related genes (MRGs) and explore its prognostic value in LUAD patients.MethodsRNA sequencing data and the corresponding clinical characteristics data of LUAD were obtained from The Cancer Genome Atlas (TCGA) database. Five membrane-tension prognosis-related genes (5-MRG) were screened by univariate and multifactorial COX regression and least absolute shrinkage and selection operator (LASSO) regression analyses. The data were then divided into testing, training, and all groups to build a prognostic model, and Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), copy number variations (CNV), tumor mutation burden (TMB), and tumor microenvironment (TME) analyses were performed to explore the potential mechanisms of MRGs. Finally, single-cell data from the GSE200972 dataset in the Gene Expression Omnibus (GEO) database were obtained to determine the distribution of prognostic MRGs.ResultsConstruction and validation of the prognostic risk models were conducted using 5-MRG in the trial, test, and all data sets. Patients in the low-risk group had a better prognosis than those in the high-risk group, and the Kaplan-Meier survival curve and receiver operating characteristic curve (ROC) confirmed that the model had a better predictive value for LUAD patients. GO and KEGG analyses of differential genes in the high- and low-risk groups were significantly enriched in immune-related pathways. Immune checkpoint (ICP) differential genes differed significantly in the high- and low-risk groups. By analyzing the single-cell sequencing data, the cells were divided into nine subpopulations and cell subpopulation localization through 5-MRG.ConclusionsThe results of this study suggest that a prognostic model based on prognosis-associated MRGs can be used to predict the prognosis of LUAD patients. Therefore, prognosis-related MRGs could be potential prognostic biomarkers and therapeutic targets.
| S-EPMC10183546 | biostudies-literature
Project description:BackgroundAutophagy related genes (ARGs) regulate lysosomal degradation to induce autophagy, and are involved in the occurrence and development of a variety of cancers. The expression of ARGs in tumor tissues has a great prospect in predicting the survival of patients. The aim of this study was to construct a prognostic risk score model for lung adenocarcinoma (LUAD) based on ARGs.Methods5,786 ARGs were obtained from GeneCards database. Gene expression profiles and clinical data of 395 LUAD patients were collected from The Cancer Genome Atlas (TCGA) database. All ARGs expression data were extracted, and The ARGs differentially expressed were identified by R software. Survival analysis of differentially expressed ARGs was performed to screen for ARGs with prognostic value, and functional enrichment analysis was performed. The least absolute selection operator (LASSO) regression and Cox regression model were used to construct a prognostic risk scoring model for ARGs. The receiver operating characteristic (ROC) curve was drawn to obtain the optimal cut-off value of risk score. According to the cut-off value, the patients were divided into high-risk group and low-risk group. The area under curve (AUC) and the Kaplan-Meier survival curve was plotted to evaluate the model performance, which was verified in external data sets. Finally, univariate and multivariate Cox regression analysis was applied to evaluate the independent prognostic value of the model, and its clinical relevance was analyzed.ResultsSurvival analysis, Lasso regression and Cox regression analysis were used to construct a LUAD prognostic risk score model with five ARGs (ADAM12, CAMP, DKK1, STRIP2 and TFAP2A). The survival time of patients with low-risk score in this model was significantly better than that of patients with high-risk score (P<0.001). The model showed good prediction performance for LUAD in both the training set (AUCmax=0.78) and two external validation sets (AUCmax=0.88). Risk score was significantly associated with the prognosis of LUAD patients in univariate and multivariate Cox regression analyses, suggested that risk score could be a potential independent prognostic factor for LUAD. Correlation analysis of clinical characteristic showed that high risk score was closely associated with high T stage, high tumor stage and poor prognosis.ConclusionsWe constructed a LUAD risk score model consisting of five ARGs, which can provide a reference for predicting the prognosis of LUAD patients, and may be used in combination with tumor node metastasis (TNM) staging for prognosis prediction of LUAD patients in the future.
| S-EPMC8387651 | biostudies-literature
Project description:ObjectivesThe immune cell infiltration (ICI) in the tumor microenvironment (TME) can provide a reference for prognosis after immunotherapy. We aim to establish an ICI scoring model and evaluate its predictive ability for the immunotherapy efficacy and the prognosis in lung adenocarcinoma (LUAD) patients.MethodsWe developed and analyzed the landscape of infiltrative immune cells based on the CIBERSORT and ESTIMATE algorithms. Then, three clusters of LUAD patients were discerned from TCGA-LUAD and GSE11969 data. Furthermore, two gene clusters were classified based on the PCA.ResultsLUAD patients with better prognoses tend to have higher immune checkpoint expression and immune/stromal scores. There is a correlation between TMB and ICI, and their relationship deserves further exploration. Moreover, the early-stage and male patients with high ICI scores have more prolonged survival.ConclusionsThe feasibility of the ICI score model in evaluating prognosis after immune checkpoint therapy for LUAD patients was verified, specifically reflected in the screening of sensitive immune checkpoints as a treatment reference. The scoring system can accurately predict the overall survival of LUAD patients, which has clinical value to monitor disease and evaluate prognosis.
| S-EPMC10086903 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) remains an incurable disease with a poor prognosis. This study aimed to explore neutrophil‑related genes (NRGs) and develop a prognostic signature for predicting the prognosis of LUAD. NRGs were obtained by intersecting modular genes identified by weighted gene co-expression network analysis (WGCNA) using bulk RNA-seq data and the marker genes of neutrophils identified from single-cell RNA-sequencing(scRNA-seq) data. Univariate Cox regression, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses were run to construct a prognostic signature, follow by delineation of risk groups, and external validation. Analyses of ESTIMAT, immune function, Tumor Immune Dysfunction and Exclusion (TIDE) scores, Immune cell Proportion Score (IPS), and immune checkpoint genes between high- and low-risk groups were performed, and then analyses of drug sensitivity to screen for sensitive anticancer drugs in high-risk groups. A total of 45 candidate NRGs were identified, of which PLTP, EREG, CD68, CD69, PLAUR, and CYP27A1 were considered to be significantly associated with prognosis in LUAD and were used to construct a prognostic signature. Correlation analysis showed significant differences in the immune landscape between high- and low-risk groups. In addition, our prognostic signature was important for predicting drug sensitivity in the high-risk group. Our study screened for NRGs in LUAD and constructed a novel and effective signature, revealing the immune landscape and providing more appropriate guidance protocols in LUAD treatment.
| S-EPMC10600204 | biostudies-literature
Project description:The purpose of this study was to construct a circular RNA (circRNA)-related competing endogenous RNA (ceRNA) regulatory network and risk score model for lung adenocarcinoma (LUAD). The relationship of the risk score to immune landscape and sensitivity to chemotherapy and targeted therapy of LUAD was assessed. We downloaded mRNA and miRNA expression data, along with clinical information, from The Cancer Genome Atlas (TCGA) program, and circRNA expression data from the Gene Expression Omnibus (GEO) database and identified differently expressed circRNA (DEcircRNA), miRNA (DEmiRNA), and mRNA (DEmRNA) using R software. We then constructed the circRNA-related network using bioinformatics method. The risk score model was established by LASSO Cox regression analysis based on 10 hub genes. In addition, the risk score model was an independent predictor for overall survival (OS) in both the TCGA and CPTAC datasets. Patients in the high-risk group had shorter OS and disease-free survival (DFS) than those in the low-risk group and were more sensitive to chemotherapy and targeted therapy. The types of tumor-infiltrating immune cells were different in the high- and low-risk groups. Our data revealed that the circRNA-related risk score model is closely associated with the level of immune cell infiltration in the tumor and the effects of adjuvant treatment. This network may be useful in designing personalized treatments for LUAD patients.
| S-EPMC8381365 | biostudies-literature