Clustering analysis and prognostic signature of lung adenocarcinoma based on the tumor microenvironment
Ontology highlight
ABSTRACT:
SUBMITTER: Shan Q
PROVIDER: S-EPMC9283441 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Lung adenocarcinoma (LUAD) represents the major histological type of lung cancer with high mortality globally. Due to the heterogeneous nature, the same treatment strategy to various patients may result in different therapeutic responses. Hence, we aimed to elaborate an effective signature for predicting patient survival outcomes. The TCGA-LUAD cohort from the TCGA portal was used as a training dataset. The GSE26939 and GSE68465 cohorts from the GEO database were taken as validation datasets. All immunologically relevant genes were extracted from the ImmPort. The ESTIMATE algorithm was employed to explore LUAD microenvironment in the training dataset. Further, the DEGs were picked out based on the immune-associated genes reflecting different statuses in the immune context of TME. Univariate/multivariate Cox regression was performed to determine six prognosis- specific genes (PIK3CG, BTK, VEGFD, INHA, INSL4, and PTPRC) and established a risk predictive signature. The time-dependent ROC indicated that AUC values were all greater than 0.70 at 1-, 3-, and 5- year intervals. Corresponding RiskScore of each LUAD patient was calculated from the signature, and they were stratified into the high- and low-risk groups by the median value of RiskScore. K-M curves and Log-rank test demonstrated significant survival differences between the two groups (P < 0.05). Similar results were exhibited in the validation datasets. The RiskScore was incredibly relevant to clinicopathological factors like gender, AJCC stage, and T stage. Also, it can mirror the distribution state of 15 kinds of TIICs and have some predictive value for the sensitivity of therapeutic drugs.
| S-EPMC8806418 | biostudies-literature
Project description:Background: Immune and stromal cells in the tumor microenvironment (TME) significantly contribute to the prognosis of lung adenocarcinoma; however, the TME-related immune prognostic signature is unknown. The aim of this study was to develop a novel immune prognostic model of the TME in lung adenocarcinoma. Methods: First, the immune and stromal scores among lung adenocarcinoma patients were determined using the ESTIMATE algorithm in accordance with The Cancer Genome Atlas (TCGA) database. Differentially expressed immune-related genes (IRGs) between high and low immune/stromal score groups were analyzed, and a univariate Cox regression analysis was performed to identify IRGs significantly correlated with overall survival (OS) among patients with lung adenocarcinoma. Furthermore, a least absolute shrinkage and selection operator (LASSO) regression analysis was performed to generate TME-related immune prognostic signatures. Gene set enrichment analysis was performed to analyze the mechanisms underlying these immune prognostic signatures. Finally, the functions of hub IRGs were further analyzed to delineate the potential prognostic mechanisms in comprehensive TCGA datasets. Results: In total, 702 intersecting differentially expressed IRGs (589 upregulated and 113 downregulated) were screened. Univariate Cox regression analysis revealed that 58 significant differentially expressed IRGs were correlated with patient prognosis in the training cohort, of which three IRGs (CLEC17A, INHA, and XIRP1) were identified through LASSO regression analysis. A robust prognostic model was generated on the basis of this three-IRG signature. Furthermore, functional enrichment analysis of the high-risk-score group was performed primarily on the basis of metabolic pathways, whereas analysis of the low-risk-score group was performed primarily on the basis of immunoregulation and immune cell activation. Finally, hub IRGs CLEC17A, INHA, and XIRP1 were considered novel prognostic biomarkers for lung adenocarcinoma. These hub genes had different mutation frequencies and forms in lung adenocarcinoma and participated in different signaling pathways. More importantly, these hub genes were significantly correlated with the infiltration of CD4+ T cells, CD8+ T cells, macrophages, B cells, and neutrophils. Conclusions: The robust novel TME-related immune prognostic signature effectively predicted the prognosis of patients with lung adenocarcinoma. Further studies are required to further elucidate the regulatory mechanisms of these hub IRGs in the TME and to develop new treatment strategies.
| S-EPMC7538811 | biostudies-literature
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.
| S-EPMC7382940 | biostudies-literature
Project description:Ferroptosis is closely linked to various cancers, including lung adenocarcinoma (LUAD); however, the factors involved in the regulation of ferroptosis-related genes are not well established. In this study, we identified and characterized ferroptosis-related long noncoding RNAs (lncRNAs) in LUAD. In particular, a coexpression network of ferroptosis-related mRNAs and lncRNAs from The Cancer Genome Atlas (TCGA) was constructed. Univariate and multivariate Cox proportional hazards analyses were performed to establish a prognostic ferroptosis-related lncRNA signature (FerRLSig). We obtained a prognostic risk model consisting of 10 ferroptosis-related lncRNAs: AL606489.1, AC106047.1, LINC02081, AC090559.1, AC026355.1, FAM83A-AS1, AL034397.3, AC092171.5, AC010980.2, and AC123595.1. High risk scores according to the FerRLSig were significantly associated with poor overall survival (hazard ratio (HR) = 1.412, 95% CI = 1.271-1.568; P < 0.001). Receiver operating characteristic (ROC) curves and a principal component analysis further supported the accuracy of the model. Next, a prognostic nomogram combining FerRLSig with clinical features was established and showed favorable predictive efficacy for survival risk stratification. In addition, gene set enrichment analysis (GSEA) revealed that FerRLSig is involved in many malignancy-associated immunoregulatory pathways. Based on the risk model, we found that the immune status and response to immunotherapy, chemotherapy, and targeted therapy differed significantly between the high-risk and low-risk groups. These results offer novel insights into the pathogenesis of LUAD, including the contribution of ferroptosis-related lncRNAs, and reveal a prognostic indicator with the potential to inform immunological research and treatment.
| S-EPMC8313561 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common histological type of lung cancers, which is the primary cause of cancer-related mortality worldwide. Growing evidence has suggested that tumor microenvironment (TME) plays a pivotal role in tumorigenesis and progression. Hence, we investigate the correlation of TME related genes with LUAD prognosis.MethodThe information of LUAD gene expression data was obtained from The Cancer Genome Atlas (TCGA). According to their immune/stromal scores calculated by the ESTIMATE algorithm, differentially expressed genes (DEGs) were identified. Then, we performed univariate Cox regression analysis on DEGs to obtain genes that are apparently bound up with LUAD survival (SurGenes). Functional annotation and protein-protein interaction (PPI) was also conducted on SurGenes. By validating the SurGenes with data sets of lung cancer from the Gene Expression Omnibus (GEO), 106 TME related SurGenes were generated. Further, intersection analysis was executed between the 106 TME related SurGenes and hub genes from PPI network, PTPRC and CD19 were obtained. Gene Set Enrichment Analysis and CIBERSORT analysis were performed on PTPRC and CD19. Based on the TCGA LUAD dataset, we conducted factor analysis and Step-wise multivariate Cox regression analysis for 106 TME related SurGenes to construct the prognostic model for LUAD survival prediction. The LUAD dataset in GEO (GSE68465) was used as the testing dataset to confirm the prognostic model. Multivariate Cox regression analysis was used between risk score from the prognostic model and clinical parameters.ResultA total of 106 TME related genes were collected in our research totally, which were markedly correlated with the overall survival (OS) of LUAD patient. Bioinformatics analysis suggest them mainly concentrated on immune response, cell adhesion, and extracellular matrix. More importantly, among 106 TME related SurGenes, PTPRC and CD19 were highly interconnected nodes among PPI network and correlated with immune activity, exhibiting significant prognostic potential. The prognostic model was a weighted linear combination of the 106 genes, by which the low-OS LUAD samples could be separated from the high-OS samples with success. This model was also able to rebustly predict the situation of survival (training set: p-value < 0.0001, area under the curve (AUC) = 0.649; testing set: p-value = 0.0009, AUC = 0.617). By combining with clinical parameters, the prognostic model was optimized. The AUC achieved 0.716 for 3 year and 0.699 for 5 year.ConclusionA series of TME-related prognostic genes were acquired in this research, which could reflect immune disorders within TME, and PTPRC and CD19 show the potential to be an indicator for LUAD prognosis and tumor microenvironment modulation. The prognostic model constructed base on those prognostic genes presented a high predictive ability, and may have clinical implications in the overall survival prediction of LUAD.
| S-EPMC7811787 | biostudies-literature
Project description:BackgroundClinically, only a minority of patients benefit from immunotherapy and few efficient biomarkers have been identified to distinguish patients who would respond to immunotherapy. The tumor microenvironment (TME) is reported to contribute to immunotherapy response, but details remain unknown. We aimed to construct a prognostic model based on the TME of lung adenocarcinoma (LUAD) to predict the prognosis and immunotherapy efficacy.MethodsWe integrated computational algorithms to describe the immune infiltrative landscape of LUAD patients. With the least absolute shrinkage and selection operator (LASSO) and Cox regression analyses, we developed a LUAD tumor microenvironment prognostic signature (LATPS). Subsequently, the immune characteristics and the benefit of immunotherapy in LATPS-defined subgroups were analyzed. RNA sequencing of tumor samples from 28 lung cancer patients treated with anti-PD-1 therapy was conducted to verify the predictive value of the LATPS.ResultsWe constructed the LATPS grounded on four genes, including UBE2T, KRT6A, IRX2, and CD3D. The LATPS-low subgroup had a better overall survival (OS) and tended to have a hot immune phenotype, which was characterized by an elevated abundance of immune cell infiltration and increased activity of immune-related pathways. Additionally, tumor immune dysfunction and exclusion (TIDE) score was markedly decreased in the LATPS-low subgroup, indicating an enhanced opportunity to benefit from immunotherapy. Survival analysis in 28 advanced lung cancer patients treated with an anti-PD-1 regimen at Nanfang hospital revealed that the LATPS-low subgroup had better immunotherapy benefit.ConclusionLATPS is an effective predictor to distinguish survival, immune characteristics, and immunotherapy benefit in LUAD patients.
| S-EPMC9729252 | biostudies-literature
Project description:BackgroundAlternative splicing (AS), a pivotal post-transcriptional process across more than 95% of human transcripts, is involved in transcript structural variations and protein complexity. Clinical implications of AS events and their interaction with tumor immunity were systematically analyzed in lung adenocarcinoma (LUAD).MethodsTranscriptome profiling as well as AS data of LUAD were retrospectively curated. Then, the network of the overall survival (OS)-relevant AS events with splicing factors was established. After screening OS-relevant AS events, a LASSO prognostic model was conducted and evaluated with ROC curves. A nomogram that integrated independent prognostic indicators was created. Immune response and immune cell infiltration were estimated with ESTIMATE, CIBERSORT, and ssGSEA algorithms. Drug sensitivity was inferred with pRRophetic package.ResultsIn total, 2415 OS-relevant AS events were identified across LUAD patients. The interaction network of splicing factors with OS-relevant AS events uncovered the underlying regulatory mechanisms of AS events in LUAD. Thereafter, a prognostic model containing 12 AS events was developed, which acted as a reliable and independent prognostic indicator following verification. A nomogram that constituted stage and risk score displayed great effectiveness in evaluating the survival likelihood. Moreover, the AS-based prognostic model was in relation to immune response and immune cell infiltration. Patients with a high-risk score displayed therapeutic superiority to cisplatin, erlotinib, gefitinib, and gemcitabine. Finally, three AS-relevant genes (CDKN2A, TTC39C, and PKIB) were identified as prognostic markers.ConclusionCollectively, our findings developed an AS event signature with powerful prognostic predictive efficacy in LUAD.
| S-EPMC9271776 | biostudies-literature
Project description: Not available
| S-EPMC11336613 | biostudies-literature
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.
| S-EPMC6886493 | biostudies-literature
Project description:Abnormal DNA methylation persists throughout carcinogenesis and cancer development. Hence, gene promoter methylation may act as a prognostic tool and provide new potential therapeutic targets for patients with lung adenocarcinoma (LUAD). In this study, to explore prognostic methylation signature, data regarding DNA methylation and RNA-seq, and clinical data of patients with LUAD from the Cancer Genome Atlas database (TCGA) were downloaded. After data preprocessing, the methylation data were divided into training (N = 405) and test sets (N = 62). Then, patients in the training set were assigned to five subgroups based on their different methylation levels using the consensus clustering method. We comprehensively analyzed the survival information, methylation levels, and clinical variables, including American Joint Committee on Cancer (AJCC) stage, tumor-node-metastasis (TNM) staging, age, smoking history, and gender of these five groups. Subsequently, we identified a 16-CpG prognostic signature and constructed a prognostic model, which was verified in the test set. Further analyses showed stable prognostic performance in the stratified cohorts. In conclusion, the new predictive DNA methylation signature proposed in this study may be used as an independent biomarker to assess the overall survival of LUAD patients and provide bioinformatics information for development of targeted therapy.
| S-EPMC7571836 | biostudies-literature