Identification of the immune cell infiltration landscape in pancreatic cancer to assist immunotherapy

2021 ◽  
Author(s):  
Zizheng Wang ◽  
Wenbo Zou ◽  
Fei Wang ◽  
Gong Zhang ◽  
Kuang Chen ◽  
...  

Background: A malignant tumor's immune environment, including infiltrating immune cell status, can be critical to patient outcomes. Recent studies have shown that immune cell infiltration (ICI) in pancreatic cancer (PC) is highly correlated with the response to immunotherapy and patient prognosis. Therefore, we aimed to create an ICI score that accurately predicts patient outcomes and immunotherapeutic efficacy. Methods: The ICI statuses of patients with PC were estimated from the publicly available The Cancer Genome Atlas (TCGA) pancreatic ductal adenocarcinoma and GSE57495 gene expression datasets using two computational algorithms (CIBERSORT and ESTIMATE). ICI and transcriptome subsets were defined using a clustering algorithm, and survival analysis was also performed. Principal component analysis was used to calculate the novel ICI score, and gene set enrichment analysis was performed to identify the pathways underlying the defined clusters. The tumor mutational burden (TMB) was further explored in TCGA cohort, and survival analysis was used to assess the capability of the ICI and TMB scores to predict overall survival. Additionally, common driver gene mutations and their differential expression in the different ICI score group were investigated. Results: The ICI landscapes of 240 patients were generated using the devised algorithm, revealing three ICI and three gene clusters whose use improved the prediction of overall survival (p = 0.019 and p < 0.001, respectively). Crucial immune checkpoint genes were differentially expressed among these subtypes; the RIG-I-LIKE and NOD-LIKE receptor signaling pathways were enriched in samples with low ICI scores (p < 0.05). We also found that the TMB scores could predict survival outcomes, whereas the ICI scores also could predict prognoses independent of TMB. Notably, ICI scores could effectively predict responses to immunotherapy. KRAS, TP53, CDKN2A, SMAD4 and TTN remained the most commonly mutated genes in PC; moreover, KRAS and TP53 mutation rates were significantly different between the two ICI score groups. Conclusions: We developed a novel ICI score that could independently predict the response to immunotherapy and survival of patients with PC. Evaluation of the ICI landscape in a larger cohort could clarify the interactions between these infiltrating cells, the tumor microenvironment and response to immunotherapy.

2021 ◽  
Vol 12 ◽  
Author(s):  
Zhuomao Mo ◽  
Daiyuan Liu ◽  
Dade Rong ◽  
Shijun Zhang

Background: Generally, hepatocellular carcinoma (HCC) exists in an immunosuppressive microenvironment that promotes tumor evasion. Hypoxia can impact intercellular crosstalk in the tumor microenvironment. This study aimed to explore and elucidate the underlying relationship between hypoxia and immunotherapy in patients with HCC.Methods: HCC genomic and clinicopathological datasets were obtained from The Cancer Genome Atlas (TCGA-LIHC), Gene Expression Omnibus databases (GSE14520) and International Cancer Genome Consortium (ICGC-LIRI). The TCGA-LIHC cases were divided into clusters based on single sample gene set enrichment analysis and hierarchical clustering. After identifying patients with immunosuppressive microenvironment with different hypoxic conditions, correlations between immunological characteristics and hypoxia clusters were investigated. Subsequently, a hypoxia-associated score was established by differential expression, univariable Cox regression, and lasso regression analyses. The score was verified by survival and receiver operating characteristic curve analyses. The GSE14520 cohort was used to validate the findings of immune cell infiltration and immune checkpoints expression, while the ICGC-LIRI cohort was employed to verify the hypoxia-associated score.Results: We identified hypoxic patients with immunosuppressive HCC. This cluster exhibited higher immune cell infiltration and immune checkpoint expression in the TCGA cohort, while similar significant differences were observed in the GEO cohort. The hypoxia-associated score was composed of five genes (ephrin A3, dihydropyrimidinase like 4, solute carrier family 2 member 5, stanniocalcin 2, and lysyl oxidase). In both two cohorts, survival analysis revealed significant differences between the high-risk and low-risk groups. In addition, compared to other clinical parameters, the established score had the highest predictive performance at both 3 and 5 years in two cohorts.Conclusion: This study provides further evidence of the link between hypoxic signals in patients and immunosuppression in HCC. Defining hypoxia-associated HCC subtypes may help reveal potential regulatory mechanisms between hypoxia and the immunosuppressive microenvironment, and our hypoxia-associated score could exhibit potential implications for future predictive models.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Shaokun Wang ◽  
Li Pang ◽  
Zuolong Liu ◽  
Xiangwei Meng

Abstract Background The change of immune cell infiltration essentially influences the process of colorectal cancer development. The infiltration of immune cells can be regulated by a variety of genes. Thus, modeling the immune microenvironment of colorectal cancer by analyzing the genes involved can be more conducive to the in-depth understanding of carcinogenesis and the progression thereof. Methods In this study, the number of stromal and immune cells in malignant tumor tissues were first estimated by using expression data (ESTIMATE) and cell-type identification with relative subsets of known RNA transcripts (CIBERSORT) to calculate the proportion of infiltrating immune cell and stromal components of colon cancer samples from the Cancer Genome Atlas database. Then the relationship between the TMN Classification and prognosis of malignant tumors was evaluated. Results By investigating differentially expressed genes using COX regression and protein-protein interaction network (PPI), the candidate hub gene serine protease inhibitor family E member 1 (SERPINE1) was found to be associated with immune cell infiltration. Gene Set Enrichment Analysis (GSEA) further projected the potential pathways with elevated SERPINE1 expression to carcinogenesis and immunity. CIBERSORT was subsequently utilized to investigate the relationship between the expression differences of SERPINE1 and immune cell infiltration and to identify eight immune cells associated with SERPINE1 expression. Conclusion We found that SERPINE1 plays a role in the remodeling of the colon cancer microenvironment and the infiltration of immune cells.


2022 ◽  
Vol 2022 ◽  
pp. 1-14
Author(s):  
Yuanyuan Feng ◽  
Xinfang Tang ◽  
Changcheng Li ◽  
Ying Su ◽  
Xiaoyu Wang ◽  
...  

Objective. ARID1A has been discovered as a potential cancer biomarker. But its role in hepatocellular carcinoma (HCC) is subject to considerable dispute. Methods. The relationship between ARID1A and clinical factors was investigated. Clinicopathological variables related to overall survival in HCC subjects were identified using Cox and Kaplan–Meier studies. The connection between immune infiltrating cells and ARID1A expression was investigated using the tumor Genome Atlas (TCGA) dataset for gene set enrichment analysis (GSEA). Finally, a cell experiment was used to confirm it. Results. The gender and cancer topography (T) categorization of HCC were linked to increased ARID1A expression. Participants with advanced levels of ARID1A expression had a worse prognosis than someone with lower levels. ARID1A was shown to be a risk indicator of overall survival on its own. ARID1A expression is inversely proportional to immune cell infiltration. In vitro, decreasing ARID1A expression substantially slowed the cell cycle and decreased HCC cell proliferation, migration, and invasion. Conclusion. The expression of ARID1A could be used to predict the outcome of HCC. It is closely related to tumor immune cell infiltration.


2022 ◽  
Vol 2022 ◽  
pp. 1-24
Author(s):  
Bin-Bin Da ◽  
Shuai Luo ◽  
Ming Huang ◽  
Fei Song ◽  
Rong Ding ◽  
...  

It has been demonstrated that the inflammatory response influences cancer development and can be used as a prognostic biomarker in various tumors. However, the relevance of genes associated with inflammatory responses in hepatocellular carcinoma (HCC) remains unknown. The Cancer Genome Atlas (TCGA) database was analyzed using weighted gene coexpression network analysis (WGCNA) and differential analysis to discover essential inflammatory response-related genes (IFRGs). Cox regression studies, both univariate and multivariate, were employed to develop a prognostic IFRGs signature. Additionally, Gene Set Enrichment Analysis (GSEA) was used to deduce the biological function of the IFRGs signature. Finally, we estimated immune cell infiltration using a single sample GSEA (ssGSEA) and x-cell. Our results revealed that, among the major HCC IFRGs, two (DNASE1L3 and KLKB1) were employed to create a predictive IFRG signature. The IFRG signature could correctly predict overall survival (O.S) as per Kaplan-Meier time-dependent roc curves analysis. It was also linked to pathological tumor stage and T stage and might be used as a prognostic predictor in HCC. GSEA analysis concluded that the IFRG signature might influence the immune response in HCC. Immunological cell infiltration and immune checkpoint molecule expression differed in the high-risk and low-risk groups. As a result of our findings, DNASILE may play a role in the tumor microenvironment. However, more research is necessary to confirm the role of DNASE1L3 and KLKB1.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xiaoxi Shi ◽  
Yuanlin Liu ◽  
Shuai Cheng ◽  
Haidi Hu ◽  
Jian Zhang ◽  
...  

BackgroundCancer stem cells (CSCs) have been proven to influence drug resistance, recurrence, and metastasis in tumors. Our study aimed to identify stemness-related prognostic biomarkers for new therapeutic strategies in adrenocortical carcinoma.MethodsRNA-seq data and clinical characteristics were downloaded from The Cancer Genome Atlas (TCGA). The stemness indexes, mDNAsi and mRNAsi, were calculated to classify all samples into low-score and high-score groups. Two algorithms, based on the R language, ESTIMATE and single-sample Gene Set Enrichment Analysis (ssGSEA) were used to assess the immune cell infiltration states of adrenocortical carcinoma patients. Weighted Gene Co-expression Network Analysis (WGCNA) was used to find genes that were related to the stemness of cancer. By bioinformatics methods, the correlations between biomarkers capable of predicting immune checkpoint inhibitors (ICIs) responses and stemness of cancer were explored.ResultsHigh-mRNAsi predicted shorter overall survival (OS) and a higher metastatic trend in adrenocortical carcinoma (ACC) patients. Compared with the low-mRNAsi group, the high-mRNAsi group had a lower ImmuneScore and StromalScroe. Twenty-two stemness-related prognostic genes were obtained by WGCNA, which focused on the function of the cell cycle and cell mitosis. Immune cell infiltration, especially CD8+T cell, increased in the low-mRNAsi group compared with the high-mRNAsi group. Lower expression of PD-L1, CTLA-4, and TIGHT was evaluated in the high-mRNAsi group.ConclusionsACC patients with high-mRNAsi have poor prognosis and less immune cell infiltration. Combined with the finding of lower expression of CTLA-4, TIGHT, and PD-L1 in the high-mRNAsi group, we came to the conclusion that stemness index is a potential biomarker to predict the effectiveness of ICIs.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Qiuxian Zheng ◽  
Qin Yang ◽  
Jiaming Zhou ◽  
Xinyu Gu ◽  
Haibo Zhou ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) has a poor prognosis and has become the sixth most common malignancy worldwide due to its high incidence. Advanced approaches to therapy, including immunotherapeutic strategies, have played crucial roles in decreasing recurrence rates and improving clinical outcomes. The HCC microenvironment is important for both tumour carcinogenesis and immunogenicity, but a classification system based on immune signatures has not yet been comprehensively described. Methods HCC datasets from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC) were used in this study. Gene set enrichment analysis (GSEA) and the ConsensusClusterPlus algorithm were used for clustering assessments. We scored immune cell infiltration and used linear discriminant analysis (LDA) to improve HCC classification accuracy. Pearson's correlation analyses were performed to assess relationships between immune signature indices and immunotherapies. In addition, weighted gene co-expression network analysis (WGCNA) was applied to identify candidate modules closely associated with immune signature indices. Results Based on 152 immune signatures from HCC samples, we identified four distinct immune subtypes (IS1, IS2, IS3, and IS4). Subtypes IS1 and IS4 had more favourable prognoses than subtypes IS2 and IS3. These four subtypes also had different immune system characteristics. The IS1 subtype had the highest scores for IFNγ, cytolysis, angiogenesis, and immune cell infiltration among all subtypes. We also identified 11 potential genes, namely, TSPAN15, TSPO, METTL9, CD276, TP53I11, SPINT1, TSPO, TRABD2B, WARS2, C9ORF116, and LBH, that may represent potential immunological biomarkers for HCC. Furthermore, real-time PCR revealed that SPINT1, CD276, TSPO, TSPAN15, METTL9, and WARS2 expression was increased in HCC cells. Conclusions The present gene-based immune signature classification and indexing may provide novel perspectives for both HCC immunotherapy management and prognosis prediction.


2020 ◽  
Vol 2020 ◽  
pp. 1-13
Author(s):  
Qian Zhang ◽  
Yuping Cai ◽  
Vadim Kurbatov ◽  
Sajid A. Khan ◽  
Lingeng Lu ◽  
...  

N6-methyladenosine (m6A) plays an important role in many cancers. However, few studies have examined the role of m6A in colorectal CRC. To examine the effect of m6A on CRC, we studied the genome of 591 CRC cases from The Cancer Genome Atlas (TCGA). The relationship between the messenger RNA (mRNA) expression, copy number variation (CNVs), and mutations of m6A “Writers,” “Readers,” and “Erasers,” prognosis, immune cell infiltration, and genetic mutations in CRC cases were analyzed. CNVs and mutations were found in thirteen m6A regulators. As expected, gain and amplification of m6A regulators increased the mRNA expression of these regulators, while deletion led to reduction in the mRNA expression. Moreover, CNVs and mutation of these regulators were significantly associated with APC, TP53, and microsatellite instability (MSI) status ( p < 0.001 , p < 0.001 , and p = 0.029 , respectively). CNVs of m6A regulators also correlated with inferred immune cell infiltration in CRC tissues, especially in colon tissues. Additionally, alterations of RBM15, YTHDF2, YTHDC1, YTHDC2, and METTL14 genes were related to the worse overall survival and disease-free survival (DFS) of CRC patients. Specifically, the deletion status of “Writers” was also correlated to the DFS of CRC patients ( p = 0.02 ). Gene set enrichment analysis found that FTO was involved in mRNA 3 ′ end processing, polyubiquitin binding, and RNA polymerase promoter elongation, while YTHDC1 was related to interferon-alpha and gamma response. In conclusion, a novel relationship was identified between CNVs and mutations of m6A regulators with prognosis and inferred immune function of CRC. These findings will improve the understanding of the relationship of m6A in CRC.


2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Xingwang Zhao ◽  
Longlong Zhang ◽  
Juan Wang ◽  
Min Zhang ◽  
Zhiqiang Song ◽  
...  

Abstract Background Systemic lupus erythematosus (SLE) is a multisystemic, chronic inflammatory disease characterized by destructive systemic organ involvement, which could cause the decreased functional capacity, increased morbidity and mortality. Previous studies show that SLE is characterized by autoimmune, inflammatory processes, and tissue destruction. Some seriously-ill patients could develop into lupus nephritis. However, the cause and underlying molecular events of SLE needs to be further resolved. Methods The expression profiles of GSE144390, GSE4588, GSE50772 and GSE81622 were downloaded from the Gene Expression Omnibus (GEO) database to obtain differentially expressed genes (DEGs) between SLE and healthy samples. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were performed by metascape etc. online analyses. The protein–protein interaction (PPI) networks of the DEGs were constructed by GENEMANIA software. We performed Gene Set Enrichment Analysis (GSEA) to further understand the functions of the hub gene, Weighted gene co‐expression network analysis (WGCNA) would be utilized to build a gene co‐expression network, and the most significant module and hub genes was identified. CIBERSORT tools have facilitated the analysis of immune cell infiltration patterns of diseases. The receiver operating characteristic (ROC) analyses were conducted to explore the value of DEGs for SLE diagnosis. Results In total, 6 DEGs (IFI27, IFI44, IFI44L, IFI6, EPSTI1 and OAS1) were screened, Biological functions analysis identified key related pathways, gene modules and co‐expression networks in SLE. IFI27 may be closely correlated with the occurrence of SLE. We found that an increased infiltration of moncytes, while NK cells resting infiltrated less may be related to the occurrence of SLE. Conclusion IFI27 may be closely related pathogenesis of SLE, and represents a new candidate molecular marker of the occurrence and progression of SLE. Moreover immune cell infiltration plays important role in the progession of SLE.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12304
Author(s):  
Zhengyuan Wu ◽  
Leilei Chen ◽  
Chaojie Jin ◽  
Jing Xu ◽  
Xingqun Zhang ◽  
...  

Background Cutaneous melanoma (CM) is a life-threatening destructive malignancy. Pyroptosis significantly correlates with programmed tumor cell death and its microenvironment through active host-tumor crosstalk. However, the prognostic value of pyroptosis-associated gene signatures in CM remains unclear. Methods Gene profiles and clinical data of patients with CM were downloaded from The Cancer Genome Atlas (TCGA) to identify differentially expressed genes associated with pyroptosis and overall survival (OS). We constructed a prognostic gene signature using LASSO analysis, then applied immune cell infiltration scores and Kaplan-Meier, Cox, and pathway enrichment analyses to determine the roles of the gene signature in CM. A validation cohort was collected from the Gene Expression Omnibus (GEO) database. Results Four pyroptosis-associated genes were identified and incorporated into a prognostic gene signature. Integrated bioinformatics findings showed that the signature correlated with patient survival and was associated with tumor growth and metastasis. The results of Gene Set Enrichment Analysis of a risk signature indicated that several enriched pathways are associated with cancer and immunity. The risk signature for immune status significantly correlated with tumor stem cells, the immune microenvironment, immune cell infiltration and immune subtypes. The expression of four pyroptosis genes significantly correlated with the OS of patients with CM and was related to the sensitivity of cancer cells to several antitumor drugs. A signature comprising four genes associated with pyroptosis offers a novel approach to the prognosis and survival of patients with CM and will facilitate the development of individualized therapy.


Sign in / Sign up

Export Citation Format

Share Document