Immune profile of the tumor microenvironment and the identification of a four-gene signature for lung adenocarcinoma.
Ontology highlight
ABSTRACT:
SUBMITTER: Fan T
PROVIDER: S-EPMC7880407 | biostudies-literature | 2020 Dec
REPOSITORIES: biostudies-literature
Publications
Immune profile of the tumor microenvironment and the identification of a four-gene signature for lung adenocarcinoma.
Fan Tao T Zhu Mingchuang M Wang Liyu L Liu Yu Y Tian He H Zheng Yujia Y Tan Fengwei F Sun Nan N Li Chunxiang C He Jie J
Aging 20201209 2
The composition and relative abundances of immune cells in the tumor microenvironment are key factors affecting the progression of lung adenocarcinomas (LUADs) and the efficacy of immunotherapy. Using the cancer gene expression dataset from The Cancer Genome Atlas (TCGA) program, we scored stromal and immune cells for tumor purity prediction by CIBERSORT and ESTMATE. Differential expression analysis was employed to identify 374 genes between the high-score group and the low-score group, which we ...[more]
Similar Datasets
Project description:Tumor microenvironment (TME) changes are related to the occurrence and development of colon adenocarcinoma (COAD). This study aimed to analyze the characteristics of the immune microenvironment in CC, as well as the microenvironment's relationship with the clinical features of CC. Based on The Cancer Genome Atlas (TCGA) and GSE39582 cohorts, the scores of 22 tumor infiltrating lymphocytes (TILs) were calculated using CIBERSORT. ConsensusClusterPlus was used for unsupervised clustering. Three TME subtypes (TMEC1, TMEC2, and TME3) were identified based on TIL scores. TMEC2 was associated with the worst prognosis. Random forest, k-means clustering, and principal component analysis were used to construct the TME score risk signature. The median TME score was used to divide the samples into high- and low-risk groups. The prognoses of the patients with high TME scores were worse than those of the patients with low TME scores. A high TME score was an independent prognostic risk factor for patients with colon cancer. The Gene Set Enrichment Analysis (GSEA) results showed that those with high TME scores were enriched in FOCAL_ADHESION, ECM_RECEPTOR_INTERACTION, and PATHWAYS_IN_CANCER. Our findings will provide a new strategy for immunotherapy in patients with CC.
| S-EPMC8020110 | 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:BackgroundLung cancer has become the most common cancer type and caused the most cancer deaths. Lung adenocarcinoma (LUAD) is one of the major types of lung cancer. Accumulating evidence suggests the tumor microenvironment is correlated with the tumor progress and the patient's outcome. This study aimed to establish a gene signature based on tumor microenvironment that can predict patients' outcomes for LUAD.MethodsDataset TCGA-LUAD, downloaded from the TCGA portal, were taken as training cohort, and dataset GSE72094, obtained from the GEO database, was set as validation cohort. In the training cohort, ESTIMATE algorithm was applied to find intersection differentially expressed genes (DEGs) among tumor microenvironment. Kaplan-Meier analysis and univariate Cox regression model were performed on intersection DEGs to preliminarily screen prognostic genes. Besides, the LASSO Cox regression model was implemented to build a multi-gene signature, which was then validated in the validation cohorts through Kaplan-Meier, Cox, and receiver operating characteristic curve (ROC) analyses. In addition, the correlation between tumor mutational burden (TMB) and risk score was evaluated by Spearman test. GSEA and immune infiltrating analyses were conducted for understanding function annotation and the role of the signature in the tumor microenvironment.ResultsAn eight-gene signature was built, and it was examined by Kaplan-Meier analysis, revealing that a significant overall survival difference was seen. The eight-gene signature was further proven to be independent of other clinico-pathologic parameters via the Cox regression analyses. Moreover, the ROC analysis demonstrated that this signature owned a better predictive power of LUAD prognosis. The eight-gene signature was correlated with TMB. Furthermore, GSEA and immune infiltrating analyses showed that the exact pathways related to the characteristics of eight-genes signature, and identified the vital roles of Mast cells resting and B cells naive in the prognosis of the eight-gene signature.ConclusionIdentifying the eight-gene signature (INSL4, SCN7A, STAP1, P2RX1, IKZF3, MS4A1, KLRB1, and ACSM5) could accurately identify patients' prognosis and had close interactions with Mast cells resting and B cells naive, which may provide insight into personalized prognosis prediction and new therapies for LUAD patients.
| S-EPMC7546815 | 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:Lung adenocarcinoma (LUAD) remains the most common deadly disease and has a poor prognosis. More and more studies have reported that mitochondrial-related genes (MTRGs) were associated with the clinical outcomes of multiple tumors solely. In this study, we aimed to develop a novel prognostic model based on MTRGs. Differentially expressed MTRGs were identified from TCGA-LUAD and GSE31210 cohorts. Univariate Cox regression analysis was utilized to screen differentially expressed MTRGs that were related to prognosis of LUAD. Then, LASSO Cox regression analysis was used to develop a prognostic signature. ESTIMATE was used for estimating the fractions of immune cell types. In this study, we identified 44 overlapping differentially expressed MTRGs in TCGA-LUAD and GSE31210 cohorts. Among 44 overlapping differentially expressed MTRGs, nine genes were associated with prognosis of LUAD. When the penalty parameter lambda was the minimum, there were six genes meeting the conditions of constructing the signature, including SERPINB5, CCNB1, FGR MAOB, SH3BP5, and CYP24A1. The survival analysis suggested that prognosis of patients in the high-risk group was significantly worse than that in the low-risk group. Cox regression analyses showed that the risk score was an independent predictor of LUAD prognosis. As with the results of ESTIMATE score, the degree of immune cell infiltration in the low-risk group was higher than that in the high-risk group, such as TIL, Treg, and B cells. In addition, TMB and cancer stem cell infiltration were higher in the low-risk group than the high-risk group. In conclusion, we developed a novel MTRG signature acting as a negative independent prognostic factor. In the future, individualized treatments and medical decision-making may benefit from using the predicted model.
| S-EPMC9159837 | 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:Background: TEA domain transcription factor 4 (TEAD4) is a member of the transcriptional enhancer factor (TEF) family of transcription factors, which is studied to be linked to the tumorigenesis and progression of various forms of cancers, including lung adenocarcinoma (LUAD). However, the specific function of this gene in the progression of LUAD remains to be explored. Method: A total of 19 genes related to the Hippo pathway were analyzed to identify the significant genes involved in LUAD progression. The TCGA-LUAD data (n = 585) from public databases were mined, and the differentially expressed genes (DEGs) in patients with the differential level of TEAD4 were identified. The univariate Cox regression, zero LASSO regression coefficients, and multivariate Cox regression were performed to identify the independent prognostic signatures. The immune microenvironment estimation in the two subgroups, including immune cell infiltration, HLA family genes, and immune checkpoint genes, was assessed. The Gene Set Enrichment Analysis (GSEA) and GO were conducted to analyze the functional enrichment of DEGs between the two risk groups. The potential drugs for the high-risk subtypes were forecasted via the mode of action (moa) module of the connectivity map (CMap) database. Results: TEAD4 was found to be significantly correlated with poor prognosis in LUAD-patients. A total of 102 DEGs in TEAD4-high vs. TEAD4-low groups were identified. Among these DEGs, four genes (CPS1, ANLN, RHOV, and KRT6A) were identified as the independent prognostic signature to conduct the Cox risk model. The immune microenvironment estimation indicated a strong relationship between the high TEAD4 expression and immunotherapeutic resistance. The GSEA and GO showed that pathways, including cell cycle regulation, were enriched in the high-risk group, while immune response-related and metabolism biological processes were enriched in the low-risk group. Several small molecular perturbagens targeting CFTR or PLA2G1B, by the mode of action (moa) modules of the glucocorticoid receptor agonist, cyclooxygenase inhibitor, and NFkB pathway inhibitor, were predicted to be suited for the high-risk subtypes based on the high TEAD4 expression. Conclusion: The current study revealed TEAD4 is an immune regulation-related predictor of prognosis and a novel therapeutic target for LUAD.
| S-EPMC9114646 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) is the most common malignancy, leading to more than 1 million related deaths each year. Due to low long-term survival rates, the exploration of molecular mechanisms underlying LUAD progression and novel prognostic predictors is urgently needed to improve LUAD treatment. In our study, cancer-specific differentially expressed genes (DEGs) were identified using the robust rank aggregation (RRA) method between tumor and normal tissues from six Gene Expression Omnibus databases (GSE43458, GSE62949, GSE68465, GSE115002, GSE116959, and GSE118370), followed by a selection of prognostic modules using weighted gene co-expression network analysis. Univariate Cox regression, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were applied to identify nine hub genes (CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6) that constructed a prognostic risk model. The RNA expressions of nine hub genes were validated in tumor and normal tissues by RNA-sequencing and single-cell RNA-sequencing, while immunohistochemistry staining from the Human Protein Atlas database showed consistent results in the protein levels. The risk model revealed that high-risk patients were associated with poor prognoses, including advanced stages and low survival rates. Furthermore, a multivariate Cox regression analysis suggested that the prognostic risk model could be an independent prognostic factor for LUAD patients. A nomogram that incorporated the signature and clinical features was additionally built for prognostic prediction. Moreover, the levels of hub genes were related to immune cell infiltration in LUAD microenvironments. A CMap analysis identified 13 small molecule drugs as potential agents based on the risk model for LUAD treatment. Thus, we identified a prognostic risk model including CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6 as novel biomarkers and validated their prognostic and predicted values for LUAD.
| S-EPMC8440811 | 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:Background: Lung cancer is the leading cause of cancer-related death globally. Hypoxia can suppress the activation of the tumor microenvironment (TME), which contributes to distant metastasis. However, the role of hypoxia-mediated TME in predicting the diagnosis and prognosis of lung adenocarcinoma (LUAD) patients remains unclear. Methods: Both RNA and clinical data from the LUAD cohort were downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Both univariate and multivariate Cox regression analyses were used to further screen prognosis-related hypoxia gene clusters. Time-dependent receiver operation characteristic (ROC) curves were established to evaluate the predictive sensitivity and specificity of the hypoxia-related risk signature. The characterization of gene set enrichment analysis (GSEA) and TME immune cell infiltration were further explored to identify hypoxia-related immune infiltration. Results: Eight hypoxia-related genes (LDHA, DCN, PGK1, PFKP, FBP1, LOX, ENO3, and CXCR4) were identified and established to construct a hypoxia-related risk signature. The high-risk group showed a poor overall survival compared to that of the low-risk group in the TCGA and GSE68465 cohorts (p < 0.0001). The AUCs for 1-, 3-, and 5-year overall survival were 0.736 vs. 0.741, 0.656 vs. 0.737, and 0.628 vs. 0.649, respectively. The high-risk group was associated with immunosuppression in the TME. Conclusion: The hypoxia-related risk signature may represent an independent biomarker that can differentiate the characteristics of TME immune cell infiltration and predict the prognosis of LUAD.
| S-EPMC8634728 | biostudies-literature