scholarly journals Signatures of Multi-Omics Reveal Distinct Tumor Immune Microenvironment Contributing to Immunotherapy in Lung Adenocarcinoma

2021 ◽  
Vol 12 ◽  
Author(s):  
Ziqi Huang ◽  
Baihui Li ◽  
Yan Guo ◽  
Lei Wu ◽  
Fan Kou ◽  
...  

BackgroundLung adenocarcinoma (LUAD) contains a variety of genomic and epigenomic abnormalities; the effective tumor markers related to these abnormalities need to be further explored.MethodsClustering analysis was performed based on DNA methylation (MET), DNA copy number variation (CNV), and mRNA expression data, and the differences in survival and tumor immune microenvironment (TIME) between subtypes were compared. Further, we evaluated the signatures in terms of both prognostic value and immunological characteristics.ResultsThere was a positive correlation between MET and CNV in LUAD. Integrative analysis of multi-omics data from 443 samples determined molecular subtypes, iC1 and iC2. The fractions of CD8+ T cells and activated CD4+ T cells were higher, the fraction of Tregs was lower, and the expression level of programmed death-ligand 1 (PD-L1) was higher in iC2 with a poor prognosis showing a higher TIDE score. We selected PTTG1, SLC2A1, and FAM83A as signatures of molecular subtypes to build a prognostic risk model and divided patients into high-risk group and low-risk group representing poor prognosis and good prognosis, respectively, which were validated in 180 patients with LUAD. Further, the low-risk group with lower TIDE score had more infiltrating immune cells. In 100 patients with LUAD, the high-risk group with an immunosuppressive state had a higher expression of PD-L1 and lower counts of CD8+ T cells and dendritic cells.ConclusionsThese findings demonstrated that combined multi-omics data could determine molecular subtypes with significant differences of prognosis and TIME in LUAD and suggested potent utility of the signatures to guide immunotherapy.

2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Xi Chen ◽  
Lijun Yan ◽  
Yu Lu ◽  
Feng Jiang ◽  
Ni Zeng ◽  
...  

Adrenocortical carcinoma (ACC) is a rare malignancy with dismal prognosis. Hypoxia is one of characteristics of cancer leading to tumor progression. For ACC, however, no reliable prognostic signature on the basis of hypoxia genes has been built. Our study aimed to develop a hypoxia-associated gene signature in ACC. Data of ACC patients were obtained from TCGA and GEO databases. The genes included in hypoxia risk signature were identified using the Cox regression analysis as well as LASSO regression analysis. GSEA was applied to discover the enriched gene sets. To detect a possible connection between the gene signature and immune cells, the CIBERSORT technique was applied. In ACC, the hypoxia signature including three genes (CCNA2, COL5A1, and EFNA3) was built to predict prognosis and reflect the immune microenvironment. Patients with high-risk scores tended to have a poor prognosis. According to the multivariate regression analysis, the hypoxia signature could be served as an independent indicator in ACC patients. GSEA demonstrated that gene sets linked to cancer proliferation and cell cycle were differentially enriched in high-risk classes. Additionally, we found that PDL1 and CTLA4 expression were significantly lower in the high-risk group than in the low-risk group, and resting NK cells displayed a significant increase in the high-risk group. In summary, the hypoxia risk signature created in our study might predict prognosis and evaluate the tumor immune microenvironment for ACC.


Author(s):  
Qingshan Huang ◽  
Yilin Lin ◽  
Chenglong Chen ◽  
Jingbing Lou ◽  
Tingting Ren ◽  
...  

Background: Abnormal expression of lncRNA is closely related to the occurrence and metastasis of osteosarcoma. The tumor immune microenvironment (TIM) is considered to be an important factor affecting the prognosis and treatment of osteosarcoma. This study aims to explore the effect of immune-related lncRNAs (IRLs) on the prognosis of osteosarcoma and its relationship with the TIM.Methods: Ninety-five osteosarcoma samples from the TARGET database were included. Iterative LASSO regression and multivariate Cox regression analysis were used to screen the IRLs signature with the optimal AUC. The predict function was used to calculate the risk score and divide osteosarcoma into a high-risk group and low-risk group based on the optimal cut-off value of the risk score. The lncRNAs in IRLs signature that affect metastasis were screened for in vitro validation. Single sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithms were used to evaluate the role of TIM in the influence of IRLs on osteosarcoma prognosis.Results: Ten IRLs constituted the IRLs signature, with an AUC of 0.96. The recurrence and metastasis rates of osteosarcoma in the high-risk group were higher than those in the low-risk group. In vitro experiments showed that knockdown of lncRNA (AC006033.2) could increase the proliferation, migration, and invasion of osteosarcoma. ssGSEA and ESTIMATE results showed that the immune cell content and immune score in the low-risk group were generally higher than those in the high-risk group. In addition, the expression levels of immune escape-related genes were higher in the high-risk group.Conclusion: The IRLs signature is a reliable biomarker for the prognosis of osteosarcoma, and they alter the prognosis of osteosarcoma. In addition, IRLs signature and patient prognosis may be related to TIM in osteosarcoma. The higher the content of immune cells in the TIM of osteosarcoma, the lower the risk score of patients and the better the prognosis. The higher the expression of immune escape-related genes, the lower the risk score of patients and the better the prognosis.


2021 ◽  
Author(s):  
Congli Jia ◽  
Fu Yang ◽  
Ruining Li

Abstract Background: Breast cancer (BC) is the most common cancer among women, with high rates of metastasis and recurrence. Some studies have confirmed that pyroptosis is an immune-related programmed cell death. However, the correlation between the expression of pyroptosis-related genes in BC and its prognosis remains unclear. Methods: In this study, we identified 38 pyroptosis-related genes that were differentially expressed between BC and normal tissues. The prognostic value of each pyroptosis-related gene was evaluated using patient data from The Cancer Genome Atlas (TCGA). The Cox regression method was performed to establish a prognostic model for 16-gene signature, classifying all BC patients in the TCGA database into a low-or high-risk group. Results: The survival rate of BC patients in the high-risk group was significantly lower than that in the low-risk group (P<0.01). Prognostic model is independent prognostic factor for BC patients compared to clinical features. Single sample gene set enrichment analysis (ssGSEA) showed a decrease for immune cells and immune function in the high-risk group. Conclusions: Pyroptosis-related genes influence the tumor immune microenvironment and can predict the prognosis of BC.


QJM ◽  
2021 ◽  
Vol 114 (Supplement_1) ◽  
Author(s):  
Aisha Yasin ◽  
Hala Talkhan ◽  
Afaf Abd El-Aleem ◽  
Rania Abo Shady ◽  
Rania El-Kabareety ◽  
...  

Abstract Background Hepatitis C virus (HCV) infection is a major cause of chronic liver disease, with approximately 71 million chronically infected individuals worldwide. After acute infection, a majority of healthy adults will develop persistent viremia. An infection with HCV activates the immune system to defend the host with a broad range of innate and adaptive immune responses. Cellular immune response shows antiviral immunity by means of virus specific CD8+ T cytotoxic lymphocytes and CD4+ T helper cells. After resolution of infection, a population of the virus-specific memory CD8+ T cells protects the host against reinfection. These subsets of memory cells differ with respect to their anatomical localization, their phenotype and their protective potential. Objectives The aim of the present study is to find out the phenotype pattern of HCV specific memory CD8+ T cells in chronic HCV patients in comparison to HCV exposed seronegative high-risk individuals. Patients and Methods This case control study included 60 individuals categorized into two main groups, (group I) included 20 chronic HCV patients and (group II= control group) included 40 HCV seronegative healthcare workers. Both groups were subjected to carboxyfluorescein diacetate succinimidyl ester (CFSE) proliferation assay in response to HCV specific peptides. Group II was divided into 2 subgroups according to proliferation index (cut-off<2) into, proliferation index +ve high-risk group (+PIHRG) (20 healthcare workers) and proliferation index -ve high-risk group (-PIHRG) (20 healthcare workers). Then flow cytometric analysis of HCV specific CD8+ T cells was done for the all subjects to determine the expression of CD127, CCR7, CD45RO and CD45RA on both proliferating and non-proliferating fractions of cells. Results We found that expression of CD45RA was significantly increased and CCR7 was significantly decreased on HCV CD8+ T cells in the chronic HCV group than in the high risk groups in response to different HCV peptides. Conclusion During persistent viremia in chronic HCV infection subthreshold signal levels established by HCV different proteins provide anergy of virus specific memory CD8+ T cells that are unable to differentiate into full effector cells upon a possible reencounter with the antigen. So, chronic HCV patients in our study showed predominance of less differentiated effector memory phenotype (CD45RA+, CCR7-) in peripheral blood in comparison to HCV exposed seronegative high-risk individuals who don’t show overt infection.


2021 ◽  
Author(s):  
Yi Wang ◽  
Gui-Qi Zhu ◽  
Di Tian ◽  
Chang-Wu Zhou ◽  
Na Li ◽  
...  

Abstract Background N6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal role in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approach of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between them enrolled in the shaping of TIME remains unclear. Methods The RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering was implemented to classify GC patients into 2 subtypes. Survival analysis, the infiltration of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal genes and risk score model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the discrepancies of response to immune checkpoints inhibitor (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect the expression pattern in 35 tumor tissues and paired adjacent normal tissue, and validated the prognostic value of risk model in the our cohort (N=35). Results The expression profiles of 23 lncRNAs were included to cluster patients into different subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and mutation rate of genes. Different immune cell infiltration pattern were also displayed between different clusters. GSEA showed that cluster1 was preferentially enriched with tumor hallmarks and tumor correlated biological pathways. Next, 9 lncRNAs were selected by LASSO regression model to construct risk model. Patients in the high risk group had poor prognosis. The prognosis value of this model was also validated in our cohort. As for predicting responses to the ICIs therapy, we found that patients from high risk group had lower TMB score and lower proportion of MSI-H subtype. Moreover, patients had distinct immunophenoscores in different risk groups. Conclusion Our study revealed that the potential interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to immune checkpoints inhibitors therapy.


2021 ◽  
Author(s):  
Wei Dalong ◽  
Xiaoling Lan ◽  
Qiang Tang ◽  
Zhiqun Huang ◽  
Qianli Tang

Abstract Background: Cutaneous melanoma is cancer that is both malignant and aggressive, with a poor prognosis. Pyroptosis can affect the prognosis of cancer patients by controlling tumor cell growth, migration, and metastasis, as well as is closely related to the tumor immune microenvironment. The significance of pyroptosis-related genes (PRGs) in cutaneous melanoma, however, is unknown. Methods: The training set and external validation sets were cutaneous melanoma samples from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO), respectively. By using univariate Cox regression analysis and selection operator (Lasso) regression model, prognostic genes for overall survival (OS) were found. Candidate genes that were screened were used to calculate risk scores and construct a PRG risk model. The Kaplan Meier curve, time-dependent receiver operating characteristic (ROC) curve, and area under the curve (AUC)were used to assess the functional and prognostic usefulness of gene signatures in the risk model. Furthermore, to speculate on the activity of immune cell infiltration and immune-related pathways in the tumor immune microenvironment and calculate corresponding scores, a single sample gene set enrichment analysis (ssGSEA) was used.Results: An eight PRGs risk signature (AIM2, CASP4, CASP5, CASP8, IL18, NLRC4, NLRP6, PRKACA) were conducted and divided all cutaneous melanoma patients in the TCGA cohort into two groups: Low-risk and High-risk. Both the training and external validation sets showed that patients in the low-risk group showed a significantly higher likelihood of survival than those in the high-risk group (p < 0.001). Except for PRKACA, all the other eight PRGs in our study appeared to be longer survival times for patients. The results of ssGSEA in terms of 16 types of immune cells and the activity of 13 immune-related pathways showed that the High-risk group had lower immune pathway activity and lower levels of immune cell infiltration. In conclusion, the PRG-signature may be a significant predictor of prognosis and may play an essential role in UM patients' tumor immunity.


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.


2021 ◽  
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 ◽  
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.


Sign in / Sign up

Export Citation Format

Share Document