scholarly journals Development and Validation of a Prognostic Gene Signature Correlated With M2 Macrophage Infiltration in Esophageal Squamous Cell Carcinoma

2021 ◽  
Vol 11 ◽  
Author(s):  
Jiannan Yao ◽  
Ling Duan ◽  
Xuying Huang ◽  
Jian Liu ◽  
Xiaona Fan ◽  
...  

BackgroundEsophageal squamous cell carcinoma (ESCC) is the most common type of esophageal cancer and the seventh most prevalent cause of cancer-related death worldwide. Tumor microenvironment (TME) has been confirmed to play an crucial role in ESCC progression, prognosis, and the response to immunotherapy. There is a need for predictive biomarkers of TME-related processes to better prognosticate ESCC outcomes.AimTo identify a novel gene signature linked with the TME to predict the prognosis of ESCC.MethodsWe calculated the immune/stromal scores of 95 ESCC samples from The Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm, and identified differentially expressed genes (DEGs) between high and low immune/stromal score patients. The key prognostic genes were further analyzed by the intersection of protein–protein interaction (PPI) networks and univariate Cox regression analysis. Finally, a risk score model was constructed using multivariate Cox regression analysis. We evaluated the associations between the risk score model and immune infiltration via the CIBERSORT algorithm. Moreover, we validated the signature using the Gene Expression Omnibus (GEO) database. Within the ten gene signature, five rarely reported genes were further validated with quantitative real time polymerase chain reaction (qRT-PCR) using an ESCC tissue cDNA microarray.ResultsA total of 133 up-regulated genes were identified as DEGs. Ten prognostic genes were selected based on intersection analysis of univariate COX regression analysis and PPI, and consisted of C1QA, C1QB, C1QC, CD86, C3AR1, CSF1R, ITGB2, LCP2, SPI1, and TYROBP (HR>1, p<0.05). The expression of 9 of these genes in the tumor samples were significantly higher compared to matched adjacent normal tissue based on the GEO database (p<0.05). Next, we assessed the ability of the ten-gene signature to predict the overall survival of ESCC patients, and found that the high-risk group had significantly poorer outcomes compared to the low-risk group using univariate and multivariate analyses in the TCGA and GEO cohorts (HR=2.104, 95% confidence interval:1.343-3.295, p=0.001; HR=1.6915, 95% confidence interval:1.053-2.717, p=0.0297). Additionally, receiver operating characteristic (ROC) curve analysis demonstrated a relatively sensitive and specific profile for the signature (1-, 2-, 3-year AUC=0.672, 0.854, 0.81). To identify the basis for these differences in the TME, we performed correlation analyses and found a significant positive correlation with M1 and M2 macrophages and CD8+ T cells, as well as a strong correlation to M2 macrophage surface markers. A nomogram based on the risk score and select clinicopathologic characteristics was constructed to predict overall survival of ESCC patients. For validation, qRT-PCR of an ESCC patient cDNA microarray was performed, and demonstrated that C1QA, C3AR1, LCP2, SPI1, and TYROBP were up-regulated in tumor samples and predict poor prognosis.ConclusionThis study established and validated a novel 10-gene signature linked with M2 macrophages and poor prognosis in ESCC patients. Importantly, we identified C1QA, C3AR1, LCP2, SPI1, and TYROBP as novel M2 macrophage-correlated survival biomarkers. These findings may identify potential targets for therapy in ESCC patients.

2021 ◽  
Author(s):  
Shaopei Ye ◽  
Wenbin Tang ◽  
Ke Huang

Abstract Background: Autophagy is a biological process to eliminate dysfunctional organelles, aggregates or even long-lived proteins. . Nevertheless, the potential function and prognostic values of autophagy in Wilms Tumor (WT) are complex and remain to be clarifed. Therefore, we proposed to systematically examine the roles of autophagy-associated genes (ARGs) in WT.Methods: Here, we obtained differentially expressed autophagy-related genes (ARGs) between healthy and Wilms tumor from Therapeutically Applicable Research To Generate Effective Treatments(TARGET) and The Cancer Genome Atlas (TCGA) database. The functionalities of the differentially expressed ARGs were analyzed using Gene Ontology. Then univariate COX regression analysis and multivariate COX regression analysis were performed to acquire nine autophagy genes related to WT patients’ survival. According to the risk score, the patients were divided into high-risk and low-risk groups. The Kaplan-Meier curve demonstrated that patients with a high-risk score tend to have a poor prognosis.Results: Eighteen DEARGs were identifed, and nine ARGs were fnally utilized to establish the FAGs based signature in the TCGA cohort. we found that patients in the high-risk group were associated with mutations in TP53. We further conducted CIBERSORT analysis, and found that the infiltration of Macrophage M1 was increased in the high-risk group. Finally, the expression levels of crucial ARGs were verifed by the experiment, which were consistent with our bioinformatics analysis.Conclusions: we emphasized the clinical significance of autophagy in WT, established a prediction system based on autophagy, and identified a promising therapeutic target of autophagy for WT.


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.


2021 ◽  
Author(s):  
Sijia Li ◽  
Hongyang Zhang ◽  
Wei Li

Abstract Background: The purpose of our study is establishing a model based on ferroptosis-related genes predicting the prognosis of patients with head and neck squamous cell carcinoma (HNSCC).Methods: In our study, transcriptome and clinical data of HNSCC patients were from The Cancer Genome Atlas, ferroptosis-related genes and pathways were from Ferroptosis Signatures Database. Differentially expressed genes (DEGs) were screened by comparing tumor and adjacent normal tissues. Functional enrichment analysis of DEGs, protein-protein interaction network and gene mutation examination were applied. Univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression were used to identified DEGs. The model was constructed by multivariate Cox regression analysis and verified by Kaplan-Meier analysis. The relationship between risk scores and other clinical features was also analyzed. Univariate and multivariate Cox analysis was used to verified the independence of our model. The model was evaluated by receiver operating characteristic analysis and calculation of the area under the curve (AUC). A nomogram model based on risk score, age, gender and TNM stages was constructed.Results: We analyzed data including 500 tumor tissues and 44 adjacent normal tissues and 259 ferroptosis-related genes, then obtained 73 DEGs. Univariate Cox regression analysis screened out 16 genes related to overall survival, and LASSO analysis fingered out 12 of them with prognostic value. A risk score model based on these 12 genes was constructed by multivariate Cox regression analysis. According to the median risk score, patients were divided into high-risk group and low-risk group. The survival rate of high-risk group was significantly lower than that of low-risk group in Kaplan-Meier curve. Risk scores were related to T and grade. Univariate and multivariate Cox analysis showed our model was an independent prognostic factor. The AUC was 0.669. The nomogram showed high accuracy predicting the prognosis of HNSCC patients.Conclusion: Our model based on 12 ferroptosis-related genes performed excellently in predicting the prognosis of HNSCC patients. Ferroptosis-related genes may be promising biomarkers for HNSCC treatment and prognosis.


2021 ◽  
Vol 16 ◽  
Author(s):  
Dongqing Su ◽  
Qianzi Lu ◽  
Yi Pan ◽  
Yao Yu ◽  
Shiyuan Wang ◽  
...  

Background: Breast cancer has plagued women for many years and caused many deaths around the world. Method: In this study, based on the weighted correlation network analysis, univariate Cox regression analysis and least absolute shrinkage and selection operator, 12 immune-related genes were selected to construct the risk score for breast cancer patients. The multivariable Cox regression analysis, gene set enrichment analysis and nomogram were also conducted in this study. Results: Good results were obtained in the survival analysis, enrichment analysis, multivariable Cox regression analysis and immune-related feature analysis. When the risk score model was applied in 22 breast cancer cohorts, the univariate Cox regression analysis demonstrated that the risk score model was significantly associated with overall survival in most of the breast cancer cohorts. Conclusion: Based on these results, we could conclude that the proposed risk score model may be a promising method, and may improve the treatment stratification of breast cancer patients in the future work.


2021 ◽  
Author(s):  
Menglin He ◽  
Cheng Hu ◽  
Jian Deng ◽  
Hui Ji ◽  
Weiqian Tian

Abstract Background: Breast cancer (BC) is a kind of cancer with high incidence and mortality in female. Conventional clinical characteristics are far from accurate to predict individual outcomes. Therefore, we aimed to develop a novel signature to predict the survival of patients with BC. Methods: We analyzed the data of a training cohort from the TCGA database and a validation cohort from GEO database. After the applications of GSEA and Cox regression analyses, a glycolysis-related signature for predicting the survival of patients with BC was developed. The signature contains AK3, CACNA1H, IL13RA1, NUP43, PGK1, and SDC1. Then, we constructed a risk score formula to classify the patients into high and low-risk groups based on the expression levels of six-gene in patients. The receiver operating characteristic (ROC) curve and the Kaplan-Meier curve were used to assess the predicted capacity of the model. Next, a nomogram was developed to predict the outcomes of patients with risk score and clinical features in 1, 3, and 5 years. We further used Human Protein Atlas (HPA) database to validate the expressions of the six biomarkers in tumor and sample tissues.Results: We constructed a six-gene signature to predict the outcomes of patients with BC. The patients in high-risk group showed poor prognosis than that in low-risk group. The AUC values were 0.719 and 0.702, showing that the prediction performance of the signature is acceptable. Additionally, Cox regression analysis revealed that these biomarkers could independently predict the prognosis of BC patients without being affected by clinical factors. The expression levels of all six biomarkers in BC tissues were higher than that in normal tissues except AK3. Conclusion: We developed a six-gene signature to predict the prognosis of patients with BC. Our signature has been proved to have the ability to make an accurate and obvious prediction and might be used to expand the prediction methods in clinical.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Gongmin Zhu ◽  
Hongwei Xia ◽  
Qiulin Tang ◽  
Feng Bi

Abstract Background Tumor metastasis is one of the leading reasons of the dismal prognosis of hepatocellular carcinoma (HCC). Epithelial-mesenchymal transition (EMT) is closely associated with tumor metastasis including HCC. The purpose of this study is to construct and validate an EMT-related gene signature for predicting the prognosis of HCC patients. Methods Gene expression data of HCC patients was downloaded from The Cancer Genome Atlas (TCGA) database. Gene set enrichment analysis (GSEA) was performed to found the EMT-related gene sets which were obviously distinct between normal samples and paired HCC samples. Cox regression analysis was used to develop an EMT-related prognostic signature, and the performance of the signature was evaluated by Kaplan–Meier curves and time-dependent receiver operating characteristic (ROC) curves. A nomogram incorporating the independent predictors was established. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect the expression levels of the hub genes in HCC cell lines, and the role of PDCD6 in the metastasis of HCC was determined by functional experiments. Results An EMT-related 5-gene signature (PDCD6, TCOF1, TRIM28, EZH2 and FAM83D) was constructed using univariate and multivariate Cox regression analysis. Based on the signature, the HCC patients were classified into high- and low-risk groups, and patients in high-risk group had a poor prognosis. Time-dependent ROC and Cox regression analyses suggested that the signature could predict HCC prognosis exactly and independently. The predictive capacity of the signature was also validated in two external cohorts. GSEA results showed that many cancer-related signaling pathways such as PI3K/Akt/mTOR pathway and TGF-β/SMAD pathway were enriched in high-risk group. The result of qRT-PCR revealed that PDCD6, TCOF1 and FAM83D were highly expressed in HCC cancer cells. Among them, PDCD6 were found to promote cell migration and invasion. Conclusion The EMT-related 5-gene signature can serve as a promising prognostic biomarker for HCC patients and may provide a novel mechanism of HCC metastasis.


2021 ◽  
Author(s):  
Fei Li ◽  
Dongcen Ge ◽  
Shu-lan Sun

Abstract Background. Ferroptosis is a newly discovered form of cell death characterized by iron-dependent lipid peroxidation. The aim of this study is to investigate the relationship between ferroptosis and the prognosis of lung adenocarcinoma (LUAD).Methods. RNA-seq data was collected from the LUAD dataset of The Cancer Genome Altas (TCGA) database. We used ferroptosis-related genes as the basis, and identify the differential expression genes (DEGs) between cancer and paracancer. The univariate Cox regression analysis were used to screen the prognostic-related genes. We divided the patients into training and validation sets. Then, we screened out key genes and built a 5 genes prognostic prediction model by the applications of the least absolute shrinkage and selection operator (LASSO) 10-fold cross-validation and the multi-variate Cox regression analysis. We divided the cases by the median value of risk score and validated this model in the validation set. Meanwhile, we analyzed the somatic mutations, and estimated the score of immune infiltration in the high- and low-risk groups, as well as performed functional enrichment analysis of DEGs.Results. The result revealed that the high-risk score triggered the worse prognosis. The maximum area under curve (AUC) of the training set and the validation set of in this study was 0.7 and 0.69. Moreover, we integrated the age, gender, and tumor stage to construct the composite nomogram. The charts indicated that the AUC of cases with survival time of 1, 3 and 5 years are 0.698, 0.71 and 0.73. In addition, the mutation frequency of patients in the high-risk group was higher than that in the low-risk group. Simultaneously, DEGs were mainly enriched in ferroptosis-related pathways by analyzing the functional results.Conclusion. This study constructed a novel LUAD prognosis prediction model base on 5 ferroptosis-related genes, which can provide a prognostic evaluation tool for the clinical therapeutic decision.


2021 ◽  
Vol 11 ◽  
Author(s):  
Chao Zhu ◽  
Liqun Gu ◽  
Mianfeng Yao ◽  
Jiang Li ◽  
Changyun Fang

The prognosis and immunotherapy response rates are unfavorable in patients with oral squamous cell carcinoma (OSCC). The tumor microenvironment is associated with tumor prognosis and progression, and the underlying mechanisms remain unclear. We obtained differentially expressed immune-related genes from OSCC mRNA data in The Cancer Genome Atlas (TCGA) database. Overall survival-related risk signature was constructed by univariate Cox regression analysis and LASSO Cox regression analysis. The prognostic performance was validated with receiver operating characteristic (ROC) analysis and Kaplan–Meier survival curves in the TCGA and Gene Expression Omnibus (GEO) datasets. The risk score was confirmed to be an independent prognostic factor and a nomogram was built to quantify the risk of outcome for each patient. Furthermore, a negative correlation was observed between the risk score and the infiltration rate of immune cells, as well as the expression of immunostimulatory and immunosuppressive molecules. Functional enrichment analysis between different risk score subtypes detected multiple immune-related biological processes, metabolic pathways, and cancer-related pathways. Thus, the immune-related gene signature can predict overall survival and contribute to the personalized management of OSCC patients.


2022 ◽  
Vol 20 (1) ◽  
Author(s):  
Junyu Huo ◽  
Ge Guan ◽  
Jinzhen Cai ◽  
Liqun Wu

Abstract Background Stromal cells in tumor microenvironment could promote immune escape through a variety of mechanisms, but there are lacking research in the field of gastric cancer (GC). Methods We identified differential expressed immune-related genes (DEIRGs) between the high- and low-stromal cell abundance GC samples in The Cancer Genome Atlas and GSE84437 datasets. A risk score was constructed basing on univariate cox regression analysis, LASSO regression analysis, and multivariate cox regression analysis in the training cohort (n=772). The median value of the risk score was used to classify patients into groups with high and low risk. We conducted external validation of the prognostic signature in four independent cohorts (GSE26253, n=432; GSE62254, n=300; GSE15459, n=191; GSE26901, n=109) from the Gene Expression Omnibus (GEO) database. The immune cell infiltration was quantified by the CIBERSORT method. Results The risk score contained 6 genes (AKT3, APOD, FAM19A5, LTBP3, NOV, and NOX4) showed good performance in predicting 5-year overall survival (OS) rate and 5-year recurrence-free survival (RFS) rate of GC patients. The risk death and recurrence of GC patients growing with the increasing risk score. The patients were clustered into three subtypes according to the infiltration of 22 kinds of immune cells quantified by the CIBERSORT method. The proportion of cluster A with the worst prognosis in the high-risk group was significantly higher than that in the low-risk group; the risk score of cluster C subtype with the best prognosis was significantly lower than that of the other two subtypes. Conclusion This study established and validated a robust prognostic model for gastric cancer by integrated analysis 1804 samples of six centers, and its mechanism was explored in combination with immune cell infiltration characterization.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Guangyao Li ◽  
Xiyi Wei ◽  
Shifeng Su ◽  
Shangqian Wang ◽  
Wei Wang ◽  
...  

Abstract Background Considerable evidence has indicated an association between the immune microenvironment and clinical outcome in ccRCC. The purpose of this study is to extensively figure out the influence of immune-related genes of tumors on the prognosis of patients with ccRCC. Methods Files containing 2498 immune-related genes were obtained from the Immunology Database and Analysis Portal (ImmPort), and the transcriptome data and clinical information relevant to patients with ccRCC were identified and downloaded from the TCGA data-base. Univariate and multivariate Cox regression analyses were used to screen out prognostic immune genes. The immune risk score model was established in light of the regression coefficient between survival and hub immune-related genes. We eventually set up a nomogram for the prediction of the overall survival for ccRCC. Kaplan-Meier (K-M) and ROC curve was used in evaluating the value of the predictive risk model. A P value of < 0.05 indicated statistically significant differences throughout data analysis. Results Via differential analysis, we found that 556 immune-related genes were expressed differentially between tumor and normal tissues (p < 0. 05). The analysis of univariate Cox regression exhibited that there was a statistical correlation between 43 immune genes and survival risk in patients with ccRCC (p < 0.05). Through Lasso-Cox regression analysis, we established an immune genetic risk scoring model based on 18 immune-related genes. The high-risk group showed a bad prognosis in K-M analysis. (p < 0.001). ROC curve showed that it was reliable of the immune risk score model to predict survival risk (5 year over survival, AUC = 0.802). The model indicated satisfactory AUC and survival correlation in the validation data set (5 year OS, Area Under Curve = 0.705, p < 0.05). From Multivariate regression analysis, the immune-risk score model plays an isolated role in the prediction of the prognosis of ccRCC. Under multivariate-Cox regression analysis, we set up a nomogram for comprehensive prediction of ccRCC patients’ survival rate. At last, it was identified that 18 immune-related genes and risk scores were not only tremendously related to clinical prognosis but also contained in a variety of carcinogenic pathways. Conclusion In general, tumor immune-related genes play essential roles in ccRCC development and progression. Our research established an unequal 18-immune gene risk index to predict the prognosis of ccRCC visually. This index was found to be an independent predictive factor for ccRCC.


Sign in / Sign up

Export Citation Format

Share Document