scholarly journals ­ Characterization of the Immune Infiltration Landscape and Identification of Prognostic Biomarkers for Esophageal Cancer

Author(s):  
Yuanmei Chen ◽  
Xinyi Huang ◽  
Haiyan Peng ◽  
Guibin Weng ◽  
Zhengrong Huang ◽  
...  

Abstract Background. Esophageal cancer (EC) is the 7th most common neoplasm and the 6th most common cause of cancer-related death worldwide. Immunotherapy is an effective treatment for EC patients. However, there are no dependable markers for predicting prognosis and immunotherapy responses in EC. Our study aims to explore prognostic models and markers in EC as well as predictors for immunotherapy. Methods. The expression profiles of EC were obtained from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO), and International Cancer Genome Consortium (ICGC) databases. Cox regression analysis was performed to construct a prognostic model. Overall survival and receiver operating characteristic curve analyses were applied to verify the accuracy of the model. The CIBERSORT algorithm was conducted to quantify the infiltration of different immune cells, and EC was grouped into three immune cell infiltration (ICI) clusters. PD-1 and PD-L1 expressions were compared between the ICI clusters. Overall survival analysis between ICI score and tumor mutation burden was conducted. The immunotherapy response of patients in different ICI score clusters was also compared. The copy number variations, somatic mutations, and single nucleotide polymorphisms were analyzed. Enrichment analyses were also performed. Results. A prognostic model was successfully constructed. Three ICI clusters were identified, and the clusters with high immune and stromal scores tended to have more immune-activated phenotypes and higher expressions of PD-1 and PDL1. The ICI score may be used as a predictor independent of tumor mutation burden. Patients with higher ICI score tended to have better immunotherapeutic responses than those with lower scores. Enrichment analyses showed that the differentially expressed genes were mostly enriched in microvillus and the KRAS and IL6/JAK/STAT3 pathways. The top eight genes with the highest mutation frequencies in EC were identified and all related to the prognosis of EC patients. Conclusions. Our study established an effective prognostic model and identified markers for predicting the prognosis and immunotherapy response of EC patients.

2021 ◽  
Author(s):  
Liusheng Wu ◽  
Xiaoqiang Li ◽  
Jixian Liu ◽  
Da Wu ◽  
Dingwang Wu ◽  
...  

Abstract Objective: Autophagy-related LncRNA genes play a vital role in the development of esophageal adenocarcinoma.Our study try to construct a prognostic model of autophagy-related LncRNA esophageal adenocarcinoma, and use this model to calculate patients with esophageal adenocarcinoma. The survival risk value of esophageal adenocarcinoma can be used to evaluate its survival prognosis. At the same time, to explore the sites of potential targeted therapy genes to provide valuable guidance for the clinical diagnosis and treatment of esophageal adenocarcinoma.Methods: Our study have downloaded 261 samples of LncRNA-related transcription and clinical data of 87 patients with esophageal adenocarcinoma from the TCGA database, and 307 autophagy-related gene data from www.autuphagy.com. We applied R software (Version 4.0.2) for data analysis, merged the transcriptome LncRNA genes, autophagy-related genes and clinical data, and screened autophagy LncRNA genes related to the prognosis of esophageal adenocarcinoma. We also performed KEGG and GO enrichment analysis and GSEA enrichment analysis in these LncRNA genes to analysis the risk characteristics and bioinformatics functions of signal transduction pathways. Univariate and multivariate Cox regression analysis were used to determine the correlation between autophagy-related LncRNA and independent risk factors. The establishment of ROC curve facilitates the evaluation of the feasibility of predicting prognostic models, and further studies the correlation between autophagy-related LncRNA and the clinical characteristics of patients with esophageal adenocarcinoma. Finally, we also used survival analysis, risk analysis and independent prognostic analysis to verify the prognosis model of esophageal adenocarcinoma.Results: We screened and identified 22 autophagic LncRNA genes that are highly correlated with the overall survival (OS) of patients with esophageal adenocarcinoma. The area under the ROC curve(AUC=0.941)and the calibration curve have a good lineup, which has statistical analysis value. In addition, univariate and multivariate Cox regression analysis showed that the autophagy LncRNA feature of this esophageal adenocarcinoma is an independent predictor of esophageal adenocarcinoma.Conclusion: These LncRNA screened and identified may participate in the regulation of cellular autophagy pathways, and at the same time affect the tumor development and prognosis of patients with esophageal adenocarcinoma. These results indicate that risk signature and nomogram are important indicators related to the prognosis of patients with esophageal adenocarcinoma.


2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Jiaqiong Lin ◽  
Yan Lin ◽  
Zena Huang ◽  
Xiaoyong Li

Background. Immunotherapy offers a novel approach for the treatment of cutaneous melanoma, but the clinical efficiency varies for individual patients. In consideration of the high cost and adverse effects of immunotherapy, it is essential to explore the predictive biomarkers of outcomes. Recently, the tumor mutation burden (TMB) has been proposed as a predictive prognosticator of the immune response. Method. RNA-seq and somatic mutation datasets of 472 cutaneous melanoma patients were downloaded from The Cancer Genome Atlas (TCGA) database to analyze mutation type and TMB. Differently expressed genes (DEGs) were identified for functional analysis. TMB-related signatures were identified via LASSO and multivariate Cox regression analysis. The association between mutants of signatures and immune cells was evaluated from the TIMER database. Furthermore, the Wilcox test was applied to assess the difference in immune infiltration calculated by the CIBERSORT algorithm in risk groupings. Results. C>T substitutions and TTN were the most common SNV and mutated gene, respectively. Patients with low TMB presented poor prognosis. DEGs were mainly implicated in skin development, cell cycle, DNA replication, and immune-associated crosstalk pathways. Furthermore, a prognostic model consisting of eight TMB-related genes was developed, which was found to be an independent risk factor for treatment outcome. The mutational status of eight TMB-related genes was associated with a low level of immune infiltration. In addition, the immune infiltrates of CD4+ and CD8+ T cells, NK cells, and M1 macrophages were higher in the low-risk group, while those of M0 and M2 macrophages were higher in the high-risk group. Conclusion. Our study demonstrated that a higher TMB was associated with favorable survival outcome in cutaneous melanoma. Moreover, a close association between prognostic model and immune infiltration was identified, providing a new potential target for immunotherapy.


2021 ◽  
Vol 10 ◽  
Author(s):  
Jun Liu ◽  
Shanqiang Zhang ◽  
Wenjie Dai ◽  
Chongwei Xie ◽  
Ji-Cheng Li

SLC41A3, as a member of the 41st family of solute carriers, participates in the transport of magnesium. The role of SLC41A3 in cancer prognosis and immune regulation has rarely been reported. This study was designed to analyze the expression status and prognostic significance of SLC41A3 in pan-cancers. The mRNA expression profiles of SLC41A3 were obtained from The Cancer Genome Atlas (TCGA), the Genotype-Tissue Expression (GTEx), the Broad Institute Cancer Cell Line Encyclopedia (CCLE), and the International Cancer Genome Consortium (ICGC). The Cox regression and Kaplan-Meier analyses were used to evaluate the prognostic value of SLC41A3 in pan-cancer. Furthermore, the correlation between SLC41A3 expression and immune cells infiltration, immune checkpoint, mismatch repair (MMR), DNA methyltransferase (DNMT), tumor mutation burden (TMB), and microsatellite instability (MSI) were calculated using data form TCGA database. The results showed that the expression of SLC41A3 was down-regulated in kidney renal clear cell carcinoma (KIRC), and was associated with poor overall survival and tumor-specific mortality. Whereas, the expression of SLC41A3 was up-regulated in liver hepatocellular carcinoma (LIHC), and the results of Cox regression analysis revealed that SLC41A3 was an independent factor for LIHC prognosis. Meanwhile, a nomogram including SLC41A3 and stage was built and exhibited good predictive power for the overall survival of LIHC patients. Additionally, correlation analysis suggested a significant correlation between SLC41A3 and TMB, MSI, MMR, DNMT, and immune cells infiltration in various cancers. The overall survival and disease-specific survival analysis revealed that the combined SLC41A3 expression and immune cell score, TMB, and MSI were significantly associated with clinical outcomes in ACC, LIHC, and UVM patients. Therefore, we proposed that SLC41A3 may serve as a potential prognostic biomarker for cancer.


2020 ◽  
Author(s):  
Shuwen Han ◽  
Kefeng Ding

Abstract Background: Colorectal cancer (CRC) is one of the most common malignancies. The purpose of this study is to construct a prognostic model for predicting the overall survival (OS) in patients with CRC. Methods: The mRNA-seq and miRNA-seq data of colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ) were downloaded from The Cancer Genome Atlas (TCGA) database. The differentially expressed RNAs (DE-RNAs) between tumor and normal tissues were screened. The Kaplan-Meier and univariate Cox regression analysis were used to screen the survival-related genes. Functional enrichment analysis of survival-related genes was conducted, followed by protein-protein interaction (PPI) analysis. Subsequently, the potential drugs targeting differentially expressed mRNAs (DE-mRNAs) were investigated. Multivariate Cox regression analysis was then conducted to screen the independent prognostic factors, and these genes were used to establish a prognostic model. A receiver operator characteristic (ROC) curve was constructed, and the area under the curve (AUC) value of ROC was calculated to evaluate the specificity and sensitivity of the model. Results: A total of 855 survival-related genes were screened. These genes were mainly enriched in Gene Ontology (GO) terms, such as methylation, synapse organization, and methyltransferase activity; and pathway analysis showed that these genes were significantly involved in N-Glycan biosynthesis and the calcium signaling pathway. PPI analysis showed that aminolevulinate dehydratase (ALAD) and cholinergic receptor muscarinic 2 (CHRM2) served vital roles in the development of CRC. Aminolevulinic acid, levulinic acid, and loxapine might be potential drugs for CRC treatment. The prognostic models were built and the patients were divided into high-risk and low-risk groups based on the median of risk score (RS) as screening threshold. The OS for patients in the high-risk group was markedly shorter than that for patients in the low-risk group. Meanwhile, kazal type serine peptidase inhibitor domain 1 (KAZALD1), hippocalcin like 4 (HPCAL4), cadherin 8 (CDH8), synaptopodin 2 (SYNPO2), cyclin D3 (CCND3), and hsa_mir_26b may be independent prognostic factors that could be considered as therapeutic targets for CRC.Conclusion: We established prognostic models that could predict the OS for CRC patients and may assist clinicians in providing personalized and precision treatment in this patient population.Highlights:1. ALAD served a vital role in the development of CRC.2. CHRM2 played a role in CRC development by affecting the calcium signaling pathway.3. Aminolevulinic acid, levulinic acid, and loxapine might be potential drugs for treating CRC.4. KAZALD1 and HPCAL4 were associated with the OS of CRC.5. CDH8, SYNPO2, CCND3, and hsa-mir-26b were closely related to the prognostic of CRC staging.


2021 ◽  
Vol 12 ◽  
Author(s):  
Junyu Huo ◽  
Liqun Wu ◽  
Yunjin Zang

Recently, growing evidence has revealed the significant effect of reprogrammed metabolism on pancreatic cancer in relation to carcinogenesis, progression, and treatment. However, the prognostic value of metabolism-related genes in pancreatic cancer has not been fully revealed. We identified 379 differentially expressed metabolic-related genes (DEMRGs) by comparing 178 pancreatic cancer tissues with 171 normal pancreatic tissues in The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression project (GTEx) databases. Then, we used univariate Cox regression analysis together with Lasso regression for constructing a prognostic model consisting of 15 metabolic genes. The unified risk score formula and cutoff value were taken into account to divide patients into two groups: high risk and low risk, with the former exhibiting a worse prognosis compared with the latter. The external validation results of the International Cancer Genome Consortium (IGCC) cohort and the Gene Expression Omnibus (GEO) cohort further confirm the effectiveness of this prognostic model. As shown in the receiver operating characteristic (ROC) curve, the area under curve (AUC) values of the risk score for overall survival (OS), disease-specific survival (DSS), and progression-free survival (PFS) were 0.871, 0.885, and 0.886, respectively. Based on the Gene Set Enrichment Analysis (GSEA), the 15-gene signature can affect some important biological processes and pathways of pancreatic cancer. In addition, the prognostic model was significantly correlated with the tumor immune microenvironment (immune cell infiltration, and immune checkpoint expression, etc.) and clinicopathological features (pathological stage, lymph node, and metastasis, etc.). We also built a nomogram based on three independent prognostic predictors (including individual neoplasm status, lymph node metastasis, and risk score) for the prediction of 1-, 3-, and 5-year OS of pancreatic cancer, which may help to further improve the treatment strategy of pancreatic cancer.


2021 ◽  
Vol 18 (6) ◽  
pp. 8045-8063
Author(s):  
Han Zhao ◽  
◽  
Yun Chen ◽  
Peijun Shen ◽  
Lan Gong ◽  
...  

<abstract> <sec><title>Background</title><p>Uveal melanoma (UM) is the most aggressive intraocular tumor worldwide. Accurate prognostic models are urgently needed. The present research aimed to construct and validate a prognostic signature is associated with overall survival (OS) for UM patients based on metabolism-related genes (MRGs).</p> </sec> <sec><title>Methods</title><p>MRGs were obtained from molecular signature database (MSigDB). The gene expression profiles and patient clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. In the training datasets, MRGs were analyzed through univariate Cox regression analyses and least absolute shrinkage and selection operator (LASSO) Cox analyses to build a prognostic model. The GSE84976 was treated as the validation cohort. In addition, time-dependent receiver operating characteristic (ROC) and Kaplan-Meier survival curve analyses the reliability of the developed model. Then, gene set enrichment analysis (GSEA) was used for gene enrichment analysis. Nomogram that combined the five-gene signature was used to evaluate the predictive OS value of UM patients.</p> </sec> <sec><title>Results</title><p>Five MRGs were identified and used to establish the prognostic model for UM patients. The model was successfully validated using the testing cohort. Moreover, ROC analysis demonstrated a strong predictive ability that our prognostic signature had for UM prognosis. Multivariable Cox regression analysis revealed that the risk model was an independent predictor of prognosis. UM patients with a high-risk score showed a higher level of immune checkpoint molecules.</p> </sec> <sec><title>Conclusion</title><p>We established a novel metabolism-related signature that could predict survival and might be therapeutic targets for the treatment of UM patients.</p> </sec> </abstract>


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zhixiang Yu ◽  
Haiyan He ◽  
Yanan Chen ◽  
Qiuhe Ji ◽  
Min Sun

AbstractOvarian cancer (OV) is a common type of carcinoma in females. Many studies have reported that ferroptosis is associated with the prognosis of OV patients. However, the mechanism by which this occurs is not well understood. We utilized Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA) to identify ferroptosis-related genes in OV. In the present study, we applied Cox regression analysis to select hub genes and used the least absolute shrinkage and selection operator to construct a prognosis prediction model with mRNA expression profiles and clinical data from TCGA. A series of analyses for this signature was performed in TCGA. We then verified the identified signature using International Cancer Genome Consortium (ICGC) data. After a series of analyses, we identified six hub genes (DNAJB6, RB1, VIMP/ SELENOS, STEAP3, BACH1, and ALOX12) that were then used to construct a model using a training data set. The model was then tested using a validation data set and was found to have high sensitivity and specificity. The identified ferroptosis-related hub genes might play a critical role in the mechanism of OV development. The gene signature we identified may be useful for future clinical applications.


2021 ◽  
Vol 20 ◽  
pp. 153303382110414
Author(s):  
Xiaoyong Li ◽  
Jiaqong Lin ◽  
Yuguo pan ◽  
Peng Cui ◽  
Jintang Xia

Background: Liver progenitor cells (LPCs) play significant roles in the development and progression of hepatocellular carcinoma (HCC). However, no studies on the value of LPC-related genes for evaluating HCC prognosis exist. We developed a gene signature of LPC-related genes for prognostication in HCC. Methods: To identify LPC-related genes, we analyzed mRNA expression arrays from a dataset (GSE57812 & GSE 37071) containing LPCs, mature hepatocytes, and embryonic stem cell samples. HCC RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to explore the differentially expressed genes (DEGs) related to prognosis through DEG analysis and univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed to construct the LPC-related gene prognostic model in the TCGA training dataset. This model was validated in the TCGA testing set and an external dataset (International Cancer Genome Consortium [ICGC] dataset). Finally, we investigated the relationship between this prognostic model with tumor-node-metastasis stage, tumor grade, and vascular invasion of HCC. Results: Overall, 1770 genes were identified as LPC-related genes, of which 92 genes were identified as DEGs in HCC tissues compared with normal tissues. Furthermore, we randomly assigned patients from the TCGA dataset to the training and testing cohorts. Twenty-six DEGs correlated with overall survival (OS) in the univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed in the TCGA training set, and a 3-gene signature was constructed to stratify patients into 2 risk groups: high-risk and low-risk. Patients in the high-risk group had significantly lower OS than those in the low-risk group. Receiver operating characteristic curve analysis confirmed the signature's predictive capacity. Moreover, the risk score was confirmed to be an independent predictor for patients with HCC. Conclusion: We demonstrated that the LPC-related gene signature can be used for prognostication in HCC. Thus, targeting LPCs may serve as a therapeutic alternative for HCC.


2020 ◽  
Author(s):  
Xiaohong - Liu ◽  
Qian - Xu ◽  
Zi-Jing - Li ◽  
Bin - Xiong

Abstract BackgroundMetabolic reprogramming is an important hallmark in the development of malignancies. Numerous metabolic genes have been demonstrated to participate in the progression of hepatocellular carcinoma (HCC). However, the prognostic significance of the metabolic genes in HCC remains elusive. MethodsWe downloaded the gene expression profiles and clinical information from the GEO, TCGA and ICGC databases. The differently expressed metabolic genes were identified by using Limma R package. Univariate Cox regression analysis and LASSO (Least absolute shrinkage and selection operator) Cox regression analysis were utilized to uncover the prognostic significance of metabolic genes. A metabolism-related prognostic model was constructed in TCGA cohort and validated in ICGC cohort. Furthermore, we constructed a nomogram to improve the accuracy of the prognostic model by using the multivariate Cox regression analysis.ResultsThe high-risk score predicted poor prognosis for HCC patients in the TCGA cohort, as confirmed in the ICGC cohort (P < 0.001). And in the multivariate Cox regression analysis, we observed that risk score could act as an independent prognostic factor for the TCGA cohort (HR (hazard ratio) 3.635, 95% CI (confidence interval)2.382-5.549) and the ICGC cohort (HR1.905, 95%CI 1.328-2.731). In addition, we constructed a nomogram for clinical use, which suggested a better prognostic model than risk score.ConclusionsOur study identified several metabolic genes with important prognostic value for HCC. These metabolic genes can influence the progression of HCC by regulating tumor biology and can also provide metabolic targets for the precise treatment of HCC.


Sign in / Sign up

Export Citation Format

Share Document