Project description:BackgroundLung adenocarcinoma is the most common lung cancer subtype and accounts for the highest proportion of cancer-related deaths. The tumor microenvironment influences prognostic outcomes in lung adenocarcinoma (LUAD).Materials and methodsWe used the ESTIMATE algorithm (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) to investigate the role of microenvironment-related genes and stromal cells in lung adenocarcinoma prognosis. This analysis was done on lung adenocarcinoma cases from The Cancer Genome Atlas (TCGA). The cases were divided into high and low groups on the basis of immune and stromal scores, respectively.ResultsThere were close correlations between immune scores with prognosis and disease stage. There were 367 differentially expressed genes. Combining the Gene Expression Omnibus (GEO) database, we found 14 prognosis-related genes.ResultsThere were close correlations between immune scores with prognosis and disease stage. There were 367 differentially expressed genes. Combining the Gene Expression Omnibus (GEO) database, we found 14 prognosis-related genes. Results. Based on the enrichment levels of the immune cell types, we clustered LUAD into Immunity_H and Immunity_L subtypes. Most of these genes were upregulated in Immunity_H subtype. Finally, using the Human Protein Atlas (HPA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC) databases, most of the proteins corresponding to prognostic genes were verified to be differentially expressed between the tumor and normal groups.ConclusionsThe key genes identified in this study are involved in molecular mechanisms of LUAD.
Project description:Lung adenocarcinoma (LUAD) is a major subtype of non-small cell lung cancer. Despite significant progress in its diagnosis and treatment, the mortality and morbidity rate of LUAD remains high worldwide. The aim of the present study was to perform a systematic investigation of the tumor microenvironment (TME) and identify TME-related genes of prognostic value in patients with LUAD. Firstly, the immune scores and stromal scores of patients with LUAD from The Cancer Genome Atlas were calculated using the Estimation of STromal and Immune cells in MAlignant Tumors using Expression data algorithm, and a total of 281 prognostic TME-related genes were identified. Subsequently, functional analysis and protein-protein interaction network analysis revealed that these genes were mainly related to immune response, inflammatory response and chemotaxis. Finally, two independent LUAD cohorts from the Gene Expression Omnibus database were used to validate these genes, and 4 genes (GTPase IMAP family member 1, T-cell surface glycoprotein CD1b, integrin alpha-L and leukocyte surface antigen CD53) were identified, and downregulation of these genes was indicated to be associated with poor overall survival rate in patients with LUAD. In conclusion, a comprehensive analysis of TME was performed and 4 prognostic TME-related genes in patients with LUAD were identified.
Project description:The tumor microenvironment (TME) has been shown to be involved in angiogenesis, tumor metastasis, and immune response, thereby affecting the treatment and prognosis of patients. This study aims to identify genes that are dysregulated in the TME of patients with colon adenocarcinoma (COAD) and to evaluate their prognostic value based on RNA omics data. We obtained 512 COAD samples from the Cancer Genome Atlas (TCGA) database and 579 COAD patients from the independent dataset (GSE39582) in the Gene Expression Omnibus (GEO) database. The immune/stromal/ESTIMATE score of each patient based on their gene expression was calculated using the ESTIMATE algorithm. Kaplan-Meier survival analysis, Cox regression analysis, gene functional enrichment analysis, and protein-protein interaction (PPI) network analysis were performed. We found that immune and stromal scores were significantly correlated with COAD patients' overall survival (log rank p < 0.05). By comparing the high immune/stromal score group with the low score group, we identified 688 intersection differentially expressed genes (DEGs) from the TCGA dataset (663 upregulated and 25 downregulated). The functional enrichment analysis of intersection DEGs showed that they were mainly enriched in the immune process, cell migration, cell motility, Toll-like receptor signaling pathway, and PI3K-Akt signaling pathway. The hub genes were revealed by PPI network analysis. Through Kaplan-Meier and Cox analysis, four TME-related genes that were significantly related to the prognosis of COAD patients were verified in GSE39582. In addition, we uncovered the relationship between the four prognostic genes and immune cells in COAD. In conclusion, based on the RNA expression profiles of 1091 COAD patients, we screened four genes that can predict prognosis from the TME, which may serve as candidate prognostic biomarkers for COAD.
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common type of lung cancer and is a severe threat to human health. Although many therapies have been applied to LUAD, the long-term survival rate of patients remains unsatisfactory. We aim to find reliable immune microenvironment-related lncRNA biomarkers to improve LUAD prognosis.MethodsESTIMATE analysis was performed to evaluate the degree of immune infiltration of each patient in TAGA LUAD cohort. Correlation analysis was used to identify the immune microenvironment-related lncRNAs. Univariate cox regression analysis, LASSO analysis, and Kaplan Meier analysis were used to construct and validate the prognostic model based on microenvironment-related lncRNAs.ResultsWe obtained 1,178 immune microenvironment-related lncRNAs after correlation analysis. One hundred and eighty of them are independent prognostic lncRNAs. Sixteen key lncRNAs were selected by LASSO method. This lncRNA-based model successfully predicted patients' prognosis in validation cohort, and the risk score was related to pathological stage. Besides, we also found that TP53 had the highest frequency mutation in LUAD, and the mutation of TP53 in the high-risk group, which was identified by our survival model, has a poor prognosis. lncRNA-mRNA co-expression network further suggested that these lncRNAs play a vital role in the prognosis of LUAD.ConclusionHere, we filtered 16 key lncRNAs, which could predict the survival of LUAD and may be potential biomarkers and therapeutic targets.
Project description:Non-small cell lung cancer (NSCLC), which consists mainly of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), are the leading cause of cancer deaths worldwide. In this study, we performed a comprehensive analysis of the tumor microenvironmental and genetic factors to identify prognostic biomarkers for NSCLC. We evaluated the immune and stromal scores of patients with LUAD and LUSC using data from The Cancer Genome Atlas database with the ESTIMATE algorithm. Based on these scores, the differentially expressed genes were obtained and immune-related prognostic genes were identified. Functional analysis and protein-protein interaction network further revealed the immune-related biological processes in which these genes participated. Additionally, 22 subsets of tumor-infiltrating immune cells (TIICs) in the tumor microenvironment were analyzed with the CIBERSORT algorithm. Finally, we validated these valuable genes using an independent cohort from the Gene Expression Omnibus database. The associations of the immune and stromal scores with patients' clinical characteristics and prognosis were positive in LUAD but negative in LUSC and the correlations of TIICs with clinical characteristics were clarified. Several differentially expressed genes were identified to be potential immune-related prognostic genes. This study comprehensively analyzed the tumor microenvironment and presented immune-related prognostic biomarkers for NSCLC.
Project description:ObjectiveCurrent advances in immunotherapy requires accurate tumor sub-classification due to the heterogeneity of lung adenocarcinoma (LUAD). This study aimed to develop a LUAD sub-classification system based on immune cell signatures and identified prognostic gene markers.MethodsSignatures related to the prognosis of TCGA-LUAD and 4 GSE cohorts were screened and intersected from 184 previously published immune cell signatures. The LUAD samples in the TCGA were clustered by ConsensusClusterPlus. Molecular characteristics, immune characteristics and sensitivity to immunotherapies/chemotherapies were compared. LDA score was established through Linear Discriminant Analysis (LDA). Co-expression module was constructed by Weighted Gene Co-Expression Network Analysis (WGCNA).ResultsFour LUAD subtypes with different molecular and immune characteristics were identified. Significant differences in prognosis among the four subtypes were observed. The IS1 subtype with the worst prognosis showed the highest number of TMB, mutant genes, IFN γ score, angiogenesis score and immune score. Twenty co-expression modules were generated by WGCNA. Blue module, sky blue module and light yellow module were significantly correlated with LUAD prognosis. The hub genes (CCDC90B, ARNTL2, RIPK2, SMCO2 and ADA and NBN) showing great prognostic significance were identified from the blue module. A total of 8 hub genes (NLRC3, CLEC2D, GIMAP5, CXorf65, PARP15, AKNA, ZC3H12D, and ARRDC5) were found in the light yellow module. Except for CXorf65, the expression of the other seven genes were significantly correlated with LUAD prognosis.ConclusionThis study determined four LUAD subtypes with different molecular and immune characteristics and 13 genes closely related to the prognosis of LUAD. The current findings could help understand the heterogeneity of LUAD immune classes.
Project description:Background:Prognostic genes in the tumor microenvironment play an important role in immune biological processes and the response of cancer to immunotherapy. Thus, we aimed to assess new biomarkers that are associated with immune/stromal cells in lung adenocarcinomas (LUAD) using the ESTIMATE algorithm, which also significantly affects the prognosis of cancer. Methods:The RNA sequencing (RNA-Seq) and clinical data of LUAD were downloaded from the the Cancer Genome Atlas (TCGA ). The immune and stromal scores were calculated for each sample using the ESTIMATE algorithm. The LUAD gene chip expression profile data and the clinical data (GSE37745, GSE11969, and GSE50081) were downloaded from the Gene Expression Omnibus (GEO) for subsequent validation analysis. Differentially expressed genes were calculated between high and low score groups. Univariate Cox regression analysis was performed on differentially expressed genes (DEGs) between the two groups to obtain initial prognosis genes. These were verified by three independent LUAD cohorts from the GEO database. Multivariate Cox regression was used to identify overall survival-related DEGs. UALCAN and the Human Protein Atlas were used to analyze the mRNA /protein expression levels of the target genes. Immune cell infiltration was evaluated using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT methods, and stromal cell infiltration was assessed using xCell. Results:In this study, immune scores and stromal scores are significantly associated with the clinical characteristics of LUAD, including T stage, M stage, pathological stage, and overall survival time. 530 DEGs (18 upregulated and 512 downregulated) were found to coexist in the difference analysis with the immune scores and stromal scores subgroup. Univariate Cox regression analysis showed that 286 of the 530 DEGs were survival-related genes (p < 0.05). Of the 286 genes initially identified, nine prognosis-related genes (CSF2RB, ITK, FLT3, CD79A, CCR4, CCR6, DOK2, AMPD1, and IGJ) were validated from three separate LUAD cohorts. In addition, functional analysis of DEGs also showed that various immunoregulatory molecular pathways, including regulation of immune response and the chemokine signaling pathways, were involved. Five genes (CCR6, ITK, CCR4, DOK2, and AMPD1) were identified as independent prognostic indicators of LUAD in specific data sets. The relationship between the expression levels of these genes and immune genes was assessed. We found that CCR6 mRNA and protein expression levels of LUAD were greater than in normal tissues. We evaluated the infiltration of immune cells and stromal cells in groups with high and low levels of expression of CCR6 in the TCGA LUAD cohort. In summary, we found a series of prognosis-related genes that were associated with the LUAD tumor microenvironment.
Project description:Background:Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient's outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods:The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results:A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.
Project description:Background: The incidence of lung adenocarcinoma (LUAD) increased substantially in recent years. A systematic investigation of the metabolic genomics pattern is critical to improve the treatment and prognosis of LUAD. This study aimed to analyze the relationship between tumor microenvironment (TME) and metabolism-related genes of LUAD. Methods: The data was extracted from TCGA and GEO datasets. The metabolism-related gene expression profile and the corresponding clinical data of LUAD patients were then integrated. The survival-related genes were screened out using univariate COX regression and lasso regression analysis. The latent properties and molecular mechanisms of these LUAD-specific metabolism-related genes were investigated by computational biology. Results: A novel prognostic model was established based on 8 metabolism-related genes, including TYMS, ALDH2, PKM, GNPNAT1, LDHA, ENTPD2, NT5E, and MAOB. The immune infiltration of LUAD was also analyzed using CIBERSORT algorithms and TIMER database. In addition, the high- and low-risk groups exhibited distinct layout modes in the principal component analysis. Conclusions: In summary, our studies identified clinically significant metabolism-related genes, which were potential signature for LUAD diagnosis, monitoring, and prognosis.
Project description:Background:There is plenty of evidence showing that immune-related genes (IRGs) and epigenetic modifications play important roles in the biological process of cancer. The purpose of this study is to establish novel IRG prognostic markers by integrating mRNA expression and methylation in lung adenocarcinoma (LUAD). Methods and Results:The transcriptome profiling data and the RNA-seq data of LUAD with the corresponding clinical information of 543 LUAD cases were downloaded from The Cancer Genome Atlas (TCGA) database, which were analyzed by univariate Cox proportional regression and multivariate Cox proportional regression to develop an independent prognostic signature. On the basis of this signature, we could divide LUAD patients into the high-risk, medium-risk, and low-risk groups. Further survival analyses demonstrated that high-risk patients had significantly shorter overall survival (OS) than low-risk patients. The signature, which contains 8 IRGs (S100A16, FGF2, IGKV4-1, CX3CR1, INHA, ANGPTL4, TNFRSF11A, and VIPR1), was also validated by data from the Gene Expression Omnibus (GEO) database. We also conducted analyses of methylation levels of the relevant IRGs and their CpG sites. Meanwhile, their associations with prognosis were examined and validated by the GEO database, revealing that the methylation levels of INHA, S100A16, the CpG site cg23851011, and the CpG site cg06552037 may be used as the potential regulators for the treatment of LUAD. Conclusion:Collectively, INHA, S100A16, the CpG site cg23851011, and the CpG site cg06552037 are promising biomarkers for monitoring the outcomes of LUAD.