Exploration of Tumor Mutation Burden Combined with Immune Infiltrates in the Prognosis of Lung Adenocarcinoma

Author(s):  
Zhenyu Zhao ◽  
Boxue He ◽  
Qidong Cai ◽  
Pengfei Zhang ◽  
Xiong Peng ◽  
...  

Abstract Background: Lung adenocarcinoma (LUAD) accounts for a majority of cancer-related deaths worldwide annually. A recent study shows that immunotherapy is an effective method of LUAD treatment, and tumor mutation burden (TMB) was associated with the immune microenvironment and affected the immunotherapy. Exploration of the gene signature associated with tumor mutation burden and immune infiltrates in predicting prognosis in lung adenocarcinoma in this study, we explored the correlation of TMB with immune infiltration and prognosis in LUAD.Materials and Methods: In this study, we firstly got mutation data and LUAD RNA-Seq data of the LUAD from The Cancer Genome Atlas (TCGA), and according to the TMB we divided the patients into high/low-TMB levels groups. The gene ontology (GO) pathway enrichment analysis and KOBAS-Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis were utilized to explore the molecular function of the differentially expressed genes (DEGs) between the two groups. The function enrichment analyses of DEGs were related to the immune pathways. Then, the ESTIMATE algorithm, CIBERSORT, and ssGSEA analysis were utilized to identify the relationship between TMB subgroups and immune infiltration. According to the results, Venn analysis was utilized to select the immune-related genes in DEGs. Univariate and Lasso Cox proportional hazards regression analyses were performed to construct the signature which positively associated with the immune infiltration and affected the survival. Finally, we verified the correlation between the signature and immune infiltration. Result: The exploration of the immune infiltration suggested that high-TMB subgroups positively associated with the high level of immune infiltration in LUAD patients. According to the TMB-related immune signature, the patients were divided into High/Low-risk groups, and the high-risk group was positively associated with poor prognostic. The results of the PCA analysis confirmed the validity of the signature. We also verified the effectiveness of the signature in GSE30219 and GSE72094 datasets. The ROC curves and C-index suggested the good clinical application of the TMB-related immune signature in LUAD prognosis. Another result suggested that the patients of the high-risk group were positively associated with higher TMB levels, PD-L1expression, and immune infiltration levels.Conclusion: In conclusion, our signature provides potential biomarkers for studying aspects of the TMB in LUAD such as TMB affected immune microenvironment and prognosis. This signature may provide some biomarkers which could improve the biomarkers of PD-L1 immunotherapy response and were inverted for the clinical application of the TMB in LUAD. LUAD male patients with higher TMB-levels and risk scores may benefit from immunotherapy. The high-risk patients along with higher PD-L1 expression of the signature may suitable for immunotherapy and improve their survival by detecting the TMB of LUAD.

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Songhua Cai ◽  
Xiaotong Guo ◽  
Chujian Huang ◽  
Youjun Deng ◽  
Longde Du ◽  
...  

AbstractAngiogenesis is the process of capillary sprouting from pre-existing vessels and it plays a critical role in the carcinogenic process of lung adenocarcinoma (LUAD). However, the association of angiogenesis regulators with the prognosis and progression of LUAD needs to be further elucidated. In this study, we adopted differential expression analysis, Cox proportional hazards (PH) regression analysis and experimental validation to identify angiogenesis regulators correlated with a poor prognosis, immune infiltration and cancer progression in LUAD. These results showed that the diagnostic and prognostic models based on COL5A2 and EPHB2 served as independent biomarkers with superior predictive ability. The patients in the high-risk group exhibited a worse prognosis in the TCGA cohort (P < 0.001, HR = 1.72, 95% CI 1.28–2.30), GSE310210 cohort (P = 0.005, HR = 2.87, 95% CI 1.46–5.61), and GSE31019 cohort (P = 0.01, HR = 2.14, 95% CI 1.19–3.86) than patients in the low-risk group. The high prognostic risk patients had a higher TMB (P < 0.001); higher fractions of M0 macrophages, neutrophils, NK cells resting, and T cells CD4 memory activated (P < 0.05); and higher expression of immune checkpoints PD-1, PDL-1, PDL-2, and B7H3 (P < 0.001). Patients in the high-risk group were more sensitive to chemotherapeutic drugs and molecular targeted drugs such as cisplatin, doxorubicin, gefitinib, and bosutinib (P < 0.0001). In addition, inhibition of COL5A2 and EPHB2 effectively suppressed the proliferation and migration of LUAD cells. The current study identified angiogenesis regulators as potential biomarkers and therapeutic targets for LUAD and may help to further optimize cancer therapy.


2022 ◽  
Vol 2022 ◽  
pp. 1-17
Author(s):  
Lei Zhang ◽  
Dahai Hu ◽  
Shuchen Huangfu ◽  
Jiaxin Zhou ◽  
Wei Wang ◽  
...  

The genomic variant features (mutations, deletions, structural variants, etc.) within gastric cancer impact its evolution and immunogenicity. The tumor has developed several coping strategies to respond to these changes by DNA repair and replication (DRR). However, the intrinsic relationship between the associated DRR-related genes and gastric cancer progression remained unknown. This study selected DRR-related genes with tumor mutation burden based on the TCGA (The Cancer Genome Atlas) database of gastric cancer transcriptome and mutation data. The prognosis model of seven genes (LAMA2, CREB3L3, SELP, ABCC9, CYP1B1, CDH2, and GAMT) was constructed by a univariate and LASSO regression analysis and divided into high-risk and low-risk groups with the median risk score. Survival analysis showed that overall survival (OS) was lower in the high-risk group than that in the low-risk group. Moreover, patients with gastric cancer in the high-risk group have worse survival in different subgroups, including age, gender, histological grade, and TNM stage. The nomogram that included risk scores for DRR-related genes could accurately foresee OS of patients with gastric cancer. Interestingly, the tumor mutation burden score was higher in the low-risk group than that in the high-risk group, and the risk score for DRR-related genes was negatively correlated with tumor mutation burden in gastric cancer. Next, we further combined the risk score and tumor mutation burden to evaluate the prognosis of gastric cancer patients. The low-risk cohort had a better prognosis than the high-risk cohort in the high tumor mutation burden subgroup. The number of mutation types in the high-risk group was lower than that in the low-risk group. In the immune microenvironment of gastric cancer, more naïve B cells, memory resting CD4+ T cells, Treg cells, monocytes cells, and resting mast cells were infiltrated in the high-risk group. At last, PD-L1 and IAP expressions were negatively correlated with the risk scores; patients with gastric cancer in the low-risk group showed better immunotherapy outcomes than those in the high-risk group. Overall, the DRR-related gene signature based on tumor mutation burden is a novel biomarker for prognostic and immunotherapy response in patients with gastric cancer.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zhanghao Huang ◽  
Muqi Shi ◽  
Hao Zhou ◽  
Jinjie Wang ◽  
Hai-Jian Zhang ◽  
...  

AbstractLung adenocarcinoma (LUAD) is characterized by high infiltration and rapid growth. The function of the stem cell population is to control and maintain cell regeneration. Therefore, it is necessary to study the prognostic value of stem cell-related genes in LUAD. Signature genes were screened out from 166 stem cell-related genes according to the least absolute shrinkage operator (LASSO) and subsequently multivariate Cox regression analysis, and then established risk model. Immune infiltration and nomogram model were used to evaluate the clinical efficacy of signature. A signature consisting of 10 genes was used to dichotomize the LUAD patients into two groups (cutoff, 1.314), and then validated in GSE20319 and GSE42127. There was a significant correlation between signature and clinical characteristics. Patients with high-risk had a shorter overall survival. Furthermore, significant differences were found in multiple immune cells between the high-risk group and low-risk group. A high correlation was also reflected between signature and immune infiltration. What’s more, the signature could effectively predict the efficacy of chemotherapy in patients with LUAD, and a nomogram based on signature might accurately predict the prognosis of patients with LUAD. The signature-based of stem cell-related genes might be contributed to predicting prognosis of patients with LUAD.


Cancers ◽  
2021 ◽  
Vol 13 (1) ◽  
pp. 155
Author(s):  
Pankaj Ahluwalia ◽  
Meenakshi Ahluwalia ◽  
Ashis K. Mondal ◽  
Nikhil Sahajpal ◽  
Vamsi Kota ◽  
...  

Lung cancer is one of the leading causes of death worldwide. Cell death pathways such as autophagy, apoptosis, and necrosis can provide useful clinical and immunological insights that can assist in the design of personalized therapeutics. In this study, variations in the expression of genes involved in cell death pathways and resulting infiltration of immune cells were explored in lung adenocarcinoma (The Cancer Genome Atlas: TCGA, lung adenocarcinoma (LUAD), 510 patients). Firstly, genes involved in autophagy (n = 34 genes), apoptosis (n = 66 genes), and necrosis (n = 32 genes) were analyzed to assess the prognostic significance in lung cancer. The significant genes were used to develop the cell death index (CDI) of 21 genes which clustered patients based on high risk (high CDI) and low risk (low CDI). The survival analysis using the Kaplan–Meier curve differentiated patients based on overall survival (40.4 months vs. 76.2 months), progression-free survival (26.2 months vs. 48.6 months), and disease-free survival (62.2 months vs. 158.2 months) (Log-rank test, p < 0.01). Cox proportional hazard model significantly associated patients in high CDI group with a higher risk of mortality (Hazard Ratio: H.R 1.75, 95% CI: 1.28–2.45, p < 0.001). Differential gene expression analysis using principal component analysis (PCA) identified genes with the highest fold change forming distinct clusters. To analyze the immune parameters in two risk groups, cytokines expression (n = 265 genes) analysis revealed the highest association of IL-15RA and IL 15 (> 1.5-fold, p < 0.01) with the high-risk group. The microenvironment cell-population (MCP)-counter algorithm identified the higher infiltration of CD8+ T cells, macrophages, and lower infiltration of neutrophils with the high-risk group. Interestingly, this group also showed a higher expression of immune checkpoint molecules CD-274 (PD-L1), CTLA-4, and T cell exhaustion genes (HAVCR2, TIGIT, LAG3, PDCD1, CXCL13, and LYN) (p < 0.01). Furthermore, functional enrichment analysis identified significant perturbations in immune pathways in the higher risk group. This study highlights the presence of an immunocompromised microenvironment indicated by the higher infiltration of cytotoxic T cells along with the presence of checkpoint molecules and T cell exhaustion genes. These patients at higher risk might be more suitable to benefit from PD-L1 blockade or other checkpoint blockade immunotherapies.


2020 ◽  
Vol 29 ◽  
pp. 096368972097713
Author(s):  
Xueping Jiang ◽  
Yanping Gao ◽  
Nannan Zhang ◽  
Cheng Yuan ◽  
Yuan Luo ◽  
...  

Tumor microenvironment (TME) has critical impacts on the pathogenesis of lung adenocarcinoma (LUAD). However, the molecular mechanism of TME effects on the prognosis of LUAD patients remains unclear. Our study aimed to establish an immune-related gene pair (IRGP) model for prognosis prediction and internal mechanism investigation. Based on 702 TME-related differentially expressed genes (DEGs) extracted from The Cancer Genome Atlas (TCGA) training cohort using the ESTIMATE algorithm, a 10-IRGP signature was established to predict LUAD patient prognosis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs were significantly associated with tumor immune response. In both TCGA training and Gene Expression Omnibus validation datasets, the risk score was an independent prognostic factor for LUAD patients using Lasso-Cox analysis, and patients in the high-risk group had poorer prognosis than those in the low-risk one. In the high-risk group, M2 macrophage and neutrophil infiltrations were higher, while the levels of T cell follicular helpers were significantly lower. The gene set enrichment analysis results showed that DNA repair signaling pathways were involved. In summary, we established an IRGP signature as a potential biomarker to predict the prognosis of LUAD patients.


2021 ◽  
Author(s):  
Jing Liu ◽  
Ting Ye ◽  
Xue fang Zhang ◽  
Yong jian Dong ◽  
Wen feng Zhang ◽  
...  

Abstract Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs.Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes (DEGs) were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic (ROC) curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT, Xcell and ssGSEA in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 were significantly different in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes (ALOX5AP, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) from the nine-IRG prognostic model, of which the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, we analyzed the prognostic ability and expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 in metastatic melanoma. Overall, a prognostic model for metastatic melanoma based on the characteristics of the tumor immune microenvironment was established, which was helpful for further studies.It could function well in helping people to understand the characteristics of the immune microenvironment in metastatic melanoma and to find possible therapeutic targets.


2021 ◽  
Vol 2021 ◽  
pp. 1-18
Author(s):  
Yifang Hu ◽  
Jiahang Song ◽  
Zhen Wang ◽  
Jingbao Kan ◽  
Yaoqi Ge ◽  
...  

Background. Glioma is the most common central nervous system (CNS) cancer with a short survival period and a poor prognosis. The S100 family gene, comprising 25 members, relates to diverse biological processes of human malignancies. Nonetheless, the significance of S100 genes in predicting the prognosis of glioma remains largely unclear. We aimed to build an S100 family-based signature for glioma prognosis. Methods. We downloaded 665 and 313 glioma patients, respectively, from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database with RNAseq data and clinical information. This study established a prognostic signature based on the S100 family genes through multivariate COX and LASSO regression. The Kaplan–Meier curve was plotted to compare overall survival (OS) among groups, whereas Receiver Operating Characteristic (ROC) analysis was performed to evaluate model accuracy. A representative gene S100B was further verified by in vitro experiments. Results. An S100 family-based signature comprising 5 genes was constructed to predict the glioma that stratified TCGA-derived cases as a low- or high-risk group, whereas the significance of prognosis was verified based on CGGA-derived cases. Kaplan–Meier analysis revealed that the high-risk group was associated with the dismal prognosis. Furthermore, the S100 family-based signature was proved to be closely related to immune microenvironment. In vitro analysis showed S100B gene in the signature promoted glioblastoma (GBM) cell proliferation and migration. Conclusions. We constructed and verified a novel S100 family-based signature associated with tumor immune microenvironment (TIME), which may shed novel light on the glioma diagnosis and treatment.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11911
Author(s):  
Lei Liu ◽  
Huayu He ◽  
Yue Peng ◽  
Zhenlin Yang ◽  
Shugeng Gao

Background The prognosis of patients for lung adenocarcinoma (LUAD) is known to vary widely; the 5-year overall survival rate is just 63% even for the pathological IA stage. Thus, in order to identify high-risk patients and facilitate clinical decision making, it is vital that we identify new prognostic markers that can be used alongside TNM staging to facilitate risk stratification. Methods We used mRNA expression from The Cancer Genome Atlas (TCGA) cohort to identify a prognostic gene signature and combined this with clinical data to develop a predictive model for the prognosis of patients for lung adenocarcinoma. Kaplan-Meier curves, Lasso regression, and Cox regression, were used to identify specific prognostic genes. The model was assessed via the area under the receiver operating characteristic curve (AUC-ROC) and validated in an independent dataset (GSE50081) from the Gene Expression Omnibus (GEO). Results Our analyses identified a four-gene prognostic signature (CENPH, MYLIP, PITX3, and TRAF3IP3) that was associated with the overall survival of patients with T1-4N0-2M0 in the TCGA dataset. Multivariate regression suggested that the total risk score for the four genes represented an independent prognostic factor for the TCGA and GEO cohorts; the hazard ratio (HR) (high risk group vs low risk group) were 2.34 (p < 0.001) and 2.10 (p = 0.017). Immune infiltration estimations, as determined by an online tool (TIMER2.0) showed that CD4+ T cells were in relative abundance in the high risk group compared to the low risk group in both of the two cohorts (both p < 0.001). We established a composite prognostic model for predicting OS, combined with risk-grouping and clinical factors. The AUCs for 1-, 3-, 5- year OS in the training set were 0.750, 0.737, and 0.719; and were 0.645, 0.766, and 0.725 in the validation set. The calibration curves showed a good match between the predicted probabilities and the actual probabilities. Conclusions We identified a four-gene predictive signature which represents an independent prognostic factor and can be used to identify high-risk patients from different TNM stages of LUAD. A new prognostic model that combines a prognostic gene signature with clinical features exhibited better discriminatory ability for OS than traditional TNM staging.


2021 ◽  
Vol 2021 ◽  
pp. 1-6
Author(s):  
Yuwang Bao ◽  
Jianxiong Luo ◽  
Tianxing Yu ◽  
Yang Liu ◽  
Xiaohua Li ◽  
...  

We constructed a prognostic-related risk prediction for patients with lung adenocarcinoma by integrating multiple omics information of lung adenocarcinoma clinical information group and genome and transcriptome. Blood samples and cancer and paracancerous lung tissue samples were collected from 480 patients with lung adenocarcinoma. DNA and RNA sequencing was performed on DNA samples and RNA samples. The first follow-up was carried out 3 months after discharge. Clinical information of patients including age, gender, smoking history, and TNM stage was collected. The Cox proportional hazard model evaluated more than 600 potential SNPs related to the prognosis of lung adenocarcinoma. After LASSO analysis, we obtained 4 SNPs related to the prognosis of lung adenocarcinoma (including rs1059292, rs995343, rs2013335, and rs8078328). Through the Cox proportional hazard model, 260 candidate genes related to the prognosis of lung adenocarcinoma were evaluated. After subsequent analysis, 3 genes related to the prognosis of lung adenocarcinoma (LDHA, SDHC, and TYMS) were obtained. All survived patients were spilt into a high-risk group ( n = 170 ) and a low-risk group ( n = 170 ) according to 4 SNPs and 3 genes related to the prognosis of lung adenocarcinoma. The overall survival rate of patients in the high-risk group was lower than that in the low-risk group. The prognostic risk prediction index constructed by combining clinical information group and genomic and transcriptome characteristics of multiomics information can effectively distinguish the prognosis of patients with lung adenocarcinoma, which will provide effective support for the precise treatment of patients with lung adenocarcinoma.


Author(s):  
Junyu Huo ◽  
Jinzhen Cai ◽  
Ge Guan ◽  
Huan Liu ◽  
Liqun Wu

Background: Due to the heterogeneity of tumors and the complexity of the immune microenvironment, the specific role of ferroptosis and pyroptosis in hepatocellular carcinoma (HCC) is not fully understood, especially its impact on prognosis.Methods: The training set (n = 609, merged by TCGA and GSE14520) was clustered into three subtypes (C1, C2, and C3) based on the prognosis-related genes associated with ferroptosis and pyroptosis. The intersecting differentially expressed genes (DEGs) among C1, C2, and C3 were used in univariate Cox and LASSO penalized Cox regression analysis for the construction of the risk score. The median risk score served as the unified cutoff to divide patients into high- and low-risk groups.Results: Internal (TCGA, n = 370; GSE14520, n = 239) and external validation (ICGC, n = 231) suggested that the 12-gene risk score had high accuracy in predicting the OS, DSS, DFS, PFS, and RFS of HCC. As an independent prognostic indicator, the risk score could be applicable for patients with different clinical features tested by subgroup (n = 26) survival analysis. In the high-risk patients with a lower infiltration abundance of activated B cells, activated CD8 T cells, eosinophils, and type I T helper cells and a higher infiltration abundance of immature dendritic cells, the cytolytic activity, HLA, inflammation promotion, and type I IFN response in the high-risk group were weaker. The TP53 mutation rate, TMB, and CSC characteristics in the high-risk group were significantly higher than those in the low-risk group. Low-risk patients have active metabolic activity and a more robust immune response. The high- and low-risk groups differed significantly in histology grade, vascular tumor cell type, AFP, new tumor event after initial treatment, main tumor size, cirrhosis, TNM stage, BCLC stage, and CLIP score.Conclusion: The ferroptosis and pyroptosis molecular subtype-related signature identified and validated in this work is applicable for prognosis prediction, immune microenvironment estimation, stem cell characteristics, and clinical feature assessment in HCC.


Sign in / Sign up

Export Citation Format

Share Document