scholarly journals Profiles of Immune Infiltration and Prognostic Immunoscore in Lung Adenocarcinoma

2020 ◽  
Vol 2020 ◽  
pp. 1-15
Author(s):  
Yanyan Li ◽  
Liping Tao ◽  
Weiyang Cai

Lung tissue is abundant with immune cells that form a powerful first defense against exotic particles and microbes. The malignant phenotype of lung adenocarcinoma (LUAD) is defined not only by intrinsic tumor cells but also by tumor-infiltrating immune cells (TIICs) recruited to the immune microenvironment. Understanding more about the immune microenvironment of LUAD could function in sorting out patients more likely with high risk and benefit from immunotherapy. Twenty-two types of TIICs were estimated based on large public LUAD cohorts from the TCGA and GEO datasets using the CIBERSORT algorithm. Then principal component analysis (PCA), meta-analysis, and single-sample gene set enrichment analysis (ssGSEA) were used to measure and evaluate the specific immune responses and relative mechanisms. Moreover, an immunoscore model based on the percent of immune cells was constructed via the univariate and multivariate Cox regression models, which provided an in-depth overview of the LUAD immune microenvironment and shed light on the immune regulatory mechanism. The differential expression genes (DEGs) were acquired based on the immunoscore model, and prognostic immune-related genes were further identified. GSEA and the protein–protein interaction network (PPI) further revealed that these genes were mostly enriched in many immune-related biological processes. It is hoped that this immune landscape could provide a more accurate understanding for LUAD development and tumor immune therapy.

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11029
Author(s):  
Miaomiao Zhang ◽  
Peiyan Zheng ◽  
Yuan Wang ◽  
Baoqing Sun

Background It is well accepted that both competitive endogenous RNAs (ceRNAs) and immune microenvironment exert crucial roles in the tumor prognosis. The present study aimed to find prognostic ceRNAs and immune cells in lung adenocarcinoma (LUAD). Materials and Methods More specifically, we explored the associations of crucial ceRNAs with the immune microenvironment. The Cancer Genome Atlas (TCGA) database was employed to obtain expression profiles of ceRNAs and clinical data. CIBERSORT was utilized to quantify the proportion of 22 immune cells in LUAD. Results We constructed two cox regression models based on crucial ceRNAs and immune cells to predict prognosis in LUAD. Subsequently, seven ceRNAs and seven immune cells were involved in prognostic models. We validated both predicted models via an independent cohort GSE72094. Interestingly, both predicted models proved that the longer patients were smoking, the higher risk scores would be obtained. We further investigated the relationships between seven genes and immune/stromal scores via the ESTIMATE algorithm. The results indicated that CDC14A and H1F0 expression were significantly related to stromal scores/immune scores in LUAD. Moreover, based on the result of the ceRNA model, single-sample gene set enrichment analysis (ssGSEA) suggested that differences in immune status were evident between high- and low-risk groups.


2021 ◽  
Vol 12 ◽  
Author(s):  
Guomin Wu ◽  
Qihao Wang ◽  
Ting Zhu ◽  
Linhai Fu ◽  
Zhupeng Li ◽  
...  

This study aimed to establish a prognostic risk model for lung adenocarcinoma (LUAD). We firstly divided 535 LUAD samples in TCGA-LUAD into high-, medium-, and low-immune infiltration groups by consensus clustering analysis according to immunological competence assessment by single-sample gene set enrichment analysis (ssGSEA). Profile of long non-coding RNAs (lncRNAs) in normal samples and LUAD samples in TCGA was used for a differential expression analysis in the high- and low-immune infiltration groups. A total of 1,570 immune-related differential lncRNAs in LUAD were obtained by intersecting the above results. Afterward, univariate COX regression analysis and multivariate stepwise COX regression analysis were conducted to screen prognosis-related lncRNAs, and an eight-immune-related-lncRNA prognostic signature was finally acquired (AL365181.2, AC012213.4, DRAIC, MRGPRG-AS1, AP002478.1, AC092168.2, FAM30A, and LINC02412). Kaplan–Meier analysis and ROC analysis indicated that the eight-lncRNA-based model was accurate to predict the prognosis of LUAD patients. Simultaneously, univariate COX regression analysis and multivariate COX regression analysis were undertaken on clinical features and risk scores. It was illustrated that the risk score was a prognostic factor independent from clinical features. Moreover, immune data of LUAD in the TIMER database were analyzed. The eight-immune-related-lncRNA prognostic signature was related to the infiltration of B cells, CD4+ T cells, and dendritic cells. GSEA enrichment analysis revealed significant differences in high- and low-risk groups in pathways like pentose phosphate pathway, ubiquitin mediated proteolysis, and P53 signaling pathway. This study helps to treat LUAD patients and explore molecules related to LUAD immune infiltration to deeply understand the specific mechanism.


2021 ◽  
Vol 8 ◽  
Author(s):  
Li Zhang ◽  
Xianzhe Tang ◽  
Jia Wan ◽  
Xianghong Zhang ◽  
Tao Zheng ◽  
...  

Background: N6-methylandenosine-related long non-coding RNAs (m6A-related lncRNAs) are critically involved in cancer development. However, the roles and clinical significance of m6A-related lncRNAs in soft tissue sarcomas (STS) are inconclusive, thereby warranting further investigations.Methods: Transcriptome profiling data were extracted from The Cancer Genome Atlas (TCGA) database and Genotype-Tissue Expression (GTEx). Consensus clustering was employed to divide patients into clusters and Kaplan–Meier analysis was used to explore the prognostic differences between the subgroups. Gene set enrichment analysis (GSEA) was conducted to identify the biological processes and signaling pathways associated with m6A-Related lncRNAs. Finally, patients were randomly divided into training and validation cohorts, and least absolute shrinkage and selection operator (LASSO) Cox regression was conducted to establish the m6A-related lncRNA-based risk signature.Results: A total of 259 STS patients from TCGA-SARC dataset were enrolled in our study. Thirteen m6A-Related lncRNAs were identified to be closely related to the prognosis of STS patients. Patients were divided into two clusters, and patients in cluster 2 had a better overall survival (OS) than those in cluster 1. Patients in different clusters also showed differences in immune scores, infiltrating immune cells, and immune checkpoint expression. Patients were further classified into high-risk and low-risk subgroups according to risk scores, and high-risk patients were found to have a worse prognosis. The receiver operating characteristic (ROC) curve indicated that the risk signature displayed excellent performance at predicting the prognosis of patients with STS. Further, the risk signature was remarkably connected with the immune microenvironment and chemosensitivity in STS.Conclusion: Our study demonstrated that m6A-related lncRNAs were significantly associated with prognosis and tumor immune microenvironment and could function as independent prognosis-specific predictors in STS, thereby providing novel insights into the roles of m6A-related lncRNAs in STS.


2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Jin Zhou ◽  
Zheming Liu ◽  
Huibo Zhang ◽  
Tianyu Lei ◽  
Jiahui Liu ◽  
...  

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.


2020 ◽  
Vol 2020 ◽  
pp. 1-21
Author(s):  
Jungang Zhao ◽  
Wenming Bao ◽  
Weiyang Cai

Intrinsic cancer cells and the tumor-infiltrating immune cells (TIICs) recruited to the immune microenvironment define the malignant phenotype of lung squamous cell carcinoma (LUSC). Understanding more about the immune microenvironment of LUSC enables the selection of high-risk patients who would derive benefit from immunotherapy. Based on large public LUSC cohorts obtained from TCGA and GEO datasets, 22 types of infiltrating immune cell subgroups were evaluated by CIBERSORT. Meta-analysis, principal component analysis (PCA), single-sample gene set enrichment analysis (ssGSEA), and hierarchical clustering analysis were used to evaluate specific immune responses of LUSC. The distribution of TIICs of LUSC was entirely different from normal. TIIC subpopulations were also found to be closely associated with clinical features and molecular subtypes. Unsupervised clustering analysis revealed that three distinct TIIC subgroups existed with different survival patterns. TIICs are extensively implicated in the pathogenesis and development of LUSC. Characterizing the composition of TIICs influences the metabolism, pathological stage, and survival of tumor patients. It is hoped that this immune landscape could provide a more accurate understanding of the development and immunotherapy of LUSC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yumei Fan ◽  
Jiajie Hou ◽  
Xiaopeng Liu ◽  
Bihui Han ◽  
Yanxiu Meng ◽  
...  

Hepatocellular carcinoma (HCC) is one of the most common malignancies and ranks as the second leading cause of cancer-related mortality worldwide. Heat shock factor 2 (HSF2) is a transcription factor that plays a critical role in development, particularly corticogenesis and spermatogenesis. However, studies examining the expression and prognostic value of HSF2 and its association with tumor-infiltrating immune cells in HCC are still rare. In the present study, we found that HSF2 expression was significantly upregulated in HCC tissues compared with normal liver tissues using the TCGA, ICGC, GEO, UALCAN, HCCDB and HPA databases. High HSF2 expression was associated with shorter survival of patients with HCC. Cox regression analyses and nomogram were used to evaluate the association of HSF2 expression with the prognosis of patients with HCC. Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and gene set enrichment analysis (GSEA) revealed that HSF2 was associated with various signaling pathways, including the immune response. Notably, HSF2 expression was significantly correlated with the infiltration levels of different immune cells using the TIMER database and CIBERSORT algorithm. HSF2 expression also displayed a significant correlation with multiple immune marker sets in HCC tissues. Knockdown of HSF2 significantly inhibited the proliferation, migration, invasion and colony formation ability of HCC cells. In summary, we explored the clinical significance of HSF2 and provided a therapeutic basis for the early diagnosis, prognostic judgment, and immunotherapy of HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Qianhui Xu ◽  
Shaohuai Chen ◽  
Yuanbo Hu ◽  
Wen Huang

BackgroundIncreasing evdence supports the suggestion that the immune cell infiltration (ICI) patterns play a pivotal role in tumor progression in breast cancer (BRCA). Nonetheless, there has been no comprehensive analysis of the ICI patterns effects on the clinical outcomes and immunotherapy.MethodsMultiomic data for BRCA samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method, and CIBERSORT analysis were used to uncover the landscape of the tumor immune microenvironment (TIME). BRCA subtypes based on the ICI pattern were identified by consensus clustering and principal-component analysis was performed to obtain the ICI scores to quantify the ICI patterns in individual tumors. Their prognostic value was validated by the Kaplan-Meier survival curves. Gene set enrichment analysis (GSEA) was applied for functional annotation. Immunophenoscore (IPS) was employed to explore the immunotherapeutic role of the ICI scores. Finally, the mutation data was analyzed by using the “maftools” R package.ResultsThree different immune infiltration patterns with a distinct prognosis and biological signature were recognized among 1,198 BRCA samples. The characteristics of TIME under these three patterns were highly consistent with three known immune profiles: immune- excluded, immune-desert, and immune-inflamed phenotypes, respectively. The identification of the ICI patterns within individual tumors based on the ICI score, developed under the ICI-related signature genes, contributed into dissecting biological processes, clinical outcome, immune cells infiltration, immunotherapeutic effect, and genetic variation. High ICI score subtype, characterized with a suppression of immunity, suggested an immune-exhausted phenotype. Abundant effective immune cells were discovered in the low ICI score patients, which corresponded to an immune-activated phenotype and might present an immunotherapeutic advantage. Immunophenoscore was implemented as a surrogate of immunotherapeutic outcome, low-ICI scores samples obtained a significantly higher immunophenoscore. Enrichment of the JAK/STAT and VEGF signal pathways were activated in the ICI low-score subgroup. Finally, the synergistic effect between the ICI score and the tumor mutation burden (TMB) was confirmed.ConclusionThis work comprehensively elucidated that the ICI patterns served as an indispensable player in complexity and diversity of TIME. Quantitative identification of the ICI patterns in individual tumor will contribute into mapping the landscape of TIME further optimizing precision immunotherapy.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Lei Zhu ◽  
Fugui Yang ◽  
Lingwei Wang ◽  
Lin Dong ◽  
Zhiyuan Huang ◽  
...  

Abstract Background Ferroptosis is a recently recognized non-apoptotic cell death that is distinct from the apoptosis, necroptosis and pyroptosis. Considerable studies have demonstrated ferroptosis is involved in the biological process of various cancers. However, the role of ferroptosis in esophageal adenocarcinoma (EAC) remains unclear. This study aims to explore the ferroptosis-related genes (FRG) expression profiles and their prognostic values in EAC. Methods The FRG data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Univariate and multivariate cox regressions were used to identify the prognostic FRG, and the predictive ROC model was established using the independent risk factors. GO and KEGG enrichment analyses were performed to investigate the bioinformatics functions of significantly different genes (SDG) of ferroptosis. Additionally, the correlations of ferroptosis and immune cells were assessed through the single-sample gene set enrichment analysis (ssGSEA) and TIMER database. Finally, SDG were verified in clinical EAC specimens and normal esophageal mucosal tissues. Results Twenty-eight significantly different FRG were screened from 78 EAC and 9 normal tissues. Enrichment analyses showed these SDG were mainly related to the iron-related pathways and metabolisms of ferroptosis. Gene network demonstrated the TP53, G6PD, NFE2L2 and PTGS2 were the hub genes in the biology of ferroptosis. Cox regression analyses demonstrated four FRG (CARS1, GCLM, GLS2 and EMC2) had prognostic values for overall survival (OS) (all P < 0.05). ROC curve showed better predictive ability using the risk score (AUC = 0.744). Immune cell enrichment analysis demonstrated that the types of immune cells and their expression levels in the high-risk group were significant different with those in the low-risk group (all P < 0.05). The experimental results confirmed the ALOX5, NOX1 were upregulated and the MT1G was downregulated in the EAC tissues compared with the normal esophageal mucosal tissues (all P < 0.05). Conclusions We identified differently expressed ferroptosis-related genes that may involve in EAC. These genes have significant values in predicting the patients’ OS and targeting ferroptosis may be an alternative for therapy. Further studies are necessary to verify these results of our study.


2021 ◽  
Author(s):  
Yueren Yan ◽  
Zhendong Gao ◽  
Han Han ◽  
Yue Zhao ◽  
Yang Zhang ◽  
...  

Abstract Purpose: NRAS plays a pivotal role in progression of various kinds of somatic malignancies; however, the correlation between NRAS and lung adenocarcinoma is less known. We aim to analyze the prognostic value of NRAS expression in lung adenocarcinoma, and explore the relationship between NRAS and tumor immune microenvironment. Methods: We obtained the transcriptome pofiles and clinical data of LUAD from The Cancer Genome Atlas database and three Genome Expression Omnibus datasets. Specimens from 325 patients with completely resected lung adenocarcinoma were collected for immunohistochemical assays of NRAS, PD-L1, PD-1 and TIM-3. Then we performed gene set enrichment analysis to investigate cancer-related and immune-related signaling pathways. TIMER algorithms were performed to evaluate tumor immune infiltrating cells and immune-related biomarkers.Results: Compared with adjacent non-tumor tissue, NRAS expression was significantly upregulated in LUAD tissue. NRAS expression was significantly correlated with more advanced stage and positive lymph nodes. Kaplan-Meier curves and Cox analysis suggested that high NRAS expression led to a poor prognosis, and could be an independent prognostic factor in LUAD patients. Besides, NRAS expression was positively correlated with CD8+ T cells, macrophages, and neutrophils, and negatively correlated with B cells and CD4+ T cells. The expression level of NRAS was positively correlated with PD-L1, PD-1, and TIM-3 both at RNA and protein level. Conclusions: To conclude, we found NRAS a novel prognostic biomarker in LUAD. Besides, the expression level of NRAS may influence the prognosis of LUAD via various kinds of cancer-related pathways and remodeling TIM.


2020 ◽  
Author(s):  
Lei Zhu ◽  
Fugui Yang ◽  
Lingwei Wang ◽  
Lin Dong ◽  
Zhiyuan Huang ◽  
...  

Abstract Background Ferroptosis is a recently recognized non-apoptotic cell death that is distinct from the apoptosis, necroptosis and pyroptosis. Considerable studies have demonstrated ferroptosis is involved in the biological process of various cancers. However, the role of ferroptosis in esophageal adenocarcinoma (EAC) remains unclear. The aim of this study was to explore the ferroptosis-related genes (FRG) expression profiles and their prognostic values in EAC.Methods The FRG data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Univariate and multivariate cox regressions were used to identify the prognostic FRG, and the predictive ROC model was established using the independent risk factors. GO and KEGG enrichment analyses were performed to investigate the bioinformatics functions of significantly different genes (SDG) of ferroptosis. Additionally, the correlations of ferroptosis and immune cells were assessed through the single-sample gene set enrichment analysis (ssGSEA). Finally, significantly different genes were verified in our clinical EAC specimens and normal esophageal mucosal tissues.Results: Twenty-eight significantly different FRG were screened from 78 EAC and 9 normal tissues. GO and KEGG enrichments showed these SDG were mainly related to the iron-related pathways and metabolisms of ferroptosis. Gene network demonstrated the TP53, G6PD, NFE2L2 and PTGS2 were the hub genes in the biology of ferroptosis. Cox regression analyses demonstrated four FRG (CARS1, GCLM, GLS2 and EMC2) had prognostic values for overall survival (OS) (all P<0.001). ROC curves showed better efficacy to predict survival using the risk score (AUC=0.744). Immune cell enrichment analysis demonstrated that the types of immune cells and their expression levels in the high-risk group were significant different with those in the low-risk group (all P<0.05). The experimental results confirmed the ALOX5, NOX1 were upregulated and the MT1G was downregulated in the EAC tissues compared with the normal esophageal mucosal tissues (all P<0.05).Conclusions: We identified differently expressed ferroptosis-related genes that may involve in the process in EAC. These genes have significant values in predicting the patients’ OS and targeting ferroptosis may be an alternative for therapy. Further studies are necessary to verify these results of our study.


Sign in / Sign up

Export Citation Format

Share Document