Development of a prognostic model based on an immunogenomic landscape analysis of colorectal cancer

2020 ◽  
Author(s):  
HJ Li ◽  
YL Wang ◽  
L Ming ◽  
XQ Guo ◽  
YL Li ◽  
...  

Screening and therapeutic programs for colorectal cancer (CRC) are invasive or not effective and unable to meet patient needs. Major advances in immunogenomics may change this status but need more exploration. Differentially expressed genes and immune-related genes (IRGs) were identified by computational methods. A prognostic model was established and validated based on survival-related IRGs via stepwise multivariate Cox regression analysis. Nine IRGs were selected and identified as survival-related genes. A 7-gene prognostic model could offer a preliminary and valid determination of risk in CRC patients. The area under the curve of the receiver operating characteristic was 0.672. The 7-gene prognostic model might be used as a novel prognostic tool in CRC patients.

2020 ◽  
Vol 19 ◽  
pp. 153303382098417
Author(s):  
Ting-ting Liu ◽  
Shu-min Liu

Objective: The incidence of colorectal cancer is increasing every year, and autophagy may be related closely to the pathogenesis of colorectal cancer. Autophagy is a natural catabolic mechanism that allows the degradation of cellular components in eukaryotic cells. However, autophagy plays a dual role in tumorigenesis. It not only promotes normal cell survival and tumor growth but also induces cell death and suppresses tumors survival. In addition, the pathogenesis of various conditions, including inflammation, neurodegenerative diseases, or tumors, is associated with abnormal autophagy. The present work aimed to examine the significance of autophagy-related genes (ARGs) in prognosis prediction, to construct an autophagy prognostic model, and to identify independent prognostic factors for colorectal cancer (CRC). Methods: This study discovered a total of 36 ARGs in CRC cases using The Cancer Genome Atlas (TCGA) and Human Autophagy-dedicated (HADd) databases along with functional enrichment analysis. Then, an autophagy prognostic model was constructed using univariate Cox regression analysis, and the key prognostic genes were screened. Finally, independent prognostic markers were determined through independent prognostic analysis and clinical correlation analysis of key genes. Results: Of the 36 differentially expressed ARGs, 13 were related to prognosis, as determined by univariate Cox regression analysis. A total of 6 key genes were obtained by a multivariate Cox regression analysis. Independent prognostic values were shown by 3 genes, namely, microtubule-associated protein 1 light chain 3 (MAP1LC3C), small GTPase superfamily and Rab family (RAB7A), and WD-repeat domain phosphoinositide-interacting protein 2 (WIPI2) by independent prognostic analysis and clinical correlation. Conclusions: In this study, molecular bioinformatics technology was employed to determine and construct a prognostic model of autophagy for colon cancer patients, which revealed 3 autophagy-related features, namely, MAP1LC3C, WIPI2, and RAB7A.


2021 ◽  
Vol 11 ◽  
Author(s):  
Zhipeng Zhu ◽  
Mengyu Song ◽  
Wenhao Li ◽  
Mengying Li ◽  
Sihan Chen ◽  
...  

Hepatocellular carcinoma is a common malignant tumor with poor prognosis, poor treatment effect, and lack of effective biomarkers. In this study, bioinformatics analysis of immune-related genes of hepatocellular carcinoma was used to construct a multi-gene combined marker that can predict the prognosis of patients. The RNA expression data of hepatocellular carcinoma were downloaded from The Cancer Genome Atlas (TCGA) database, and immune-related genes were obtained from the IMMPORT database. Differential analysis was performed by Wilcox test to obtain differentially expressed genes. Univariate Cox regression analysis, lasso regression analysis and multivariate Cox regression analysis were performed to establish a prognostic model of immune genes, a total of 5 genes (HDAC1, BIRC5, SPP1, STC2, NR6A1) were identified to construct the models. The expression levels of 5 genes in HCC tissues were significantly different from those in paracancerous tissues. The Kaplan-Meier survival curve showed that the risk score calculated according to the prognostic model was significantly related to the overall survival (OS) of HCC. The receiver operating characteristic (ROC) curve confirmed that the prognostic model had high accuracy. Independent prognostic analysis was performed to prove that the risk value can be used as an independent prognostic factor. Then, the gene expression data of hepatocellular carcinoma in the ICGC database was used as a validation data set for the verification of the above steps. In addition, we used the CIBERSORT software and TIMER database to conduct immune infiltration research, and the results showed that the five genes of the model and the risk score have a certain correlation with the content of immune cells. Moreover, through Gene Set Enrichment Analysis (GSEA) and the construction of protein interaction networks, we found that the p53-mediated signal transduction pathway is a potentially important signal pathway for hepatocellular carcinoma and is positively regulated by certain genes in the prognostic model. In conclusion, this study provides potential targets for predicting the prognosis and treatment of hepatocellular carcinoma patients, and also provides new ideas about the correlation between immune genes and potential pathways of hepatocellular carcinoma.


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):  
Zhen Kang ◽  
Wei Li ◽  
Yan-Hong Yu ◽  
Meng Che ◽  
Mao-Lin Yang ◽  
...  

Background:To identify the immune-related genes of bladder cancer (BLCA) based on immunological characteristics and explore their correlation with the prognosis. Methods:We downloaded the gene and clinical data of BLCA from the Cancer Genome Atlas (TCGA) as the training group, and obtained immune-related genes from the Immport database. We downloaded GSE31684 and GSE39281 from the Gene Expression Omnibus (GEO) as the external validation group. R (version 4.0.5) and Perl were used to analyze all data. Result:Univariate Cox regression analysis and Lasso regression analysis revealed that 9 prognosis-related immunity genes (PIMGs) of differentially expressed immune genes (DEIGs) were significantly associated with the survival of BLCA patients (p < 0.01), of which 5 genes, including NPR2, PDGFRA, VIM, RBP1, RBP1 and TNC, increased the risk of the prognosis, while the rest, including CD3D, GNLY, LCK, and ZAP70, decreased the risk of the prognosis. Then, we used these genes to establish a prognostic model. We drew receiver operator characteristic (ROC) curves in the training group, and estimated the area under the curve (AUC) of 1-, 3- and 5-year survival for this model, which were 0.688, 0.719, and 0.706, respectively. The accuracy of the prognostic model was verified by the calibration chart. Combining clinical factors, we established a nomogram. The ROC curve in the external validation group showed that the nomogram had a good predictive ability for the survival rate, with a high accuracy, and the AUC values of 1-, 3-, and 5-year survival were 0.744, 0.770, and 0.782, respectively. The calibration chart indicated that the nomogram performed similarly with the ideal model. Conclusion:We had identified nine genes, including PDGFRA, VIM, RBP1, RBP1, TNC, CD3D, GNLY, LCK, and ZAP70, which played important roles in the occurrence and development of BLCA. The prognostic model based on these genes had good accuracy in predicting the OS of patients and might be promising candidates of therapeutic targets. This study may provide a new insight for the diagnosis, treatment and prognosis of BLCA from the perspective of immunology. However, further experimental studies are necessary to reveal the underlying mechanisms by which these genes mediate the progression of BLCA.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Jie Zhu ◽  
Han Wang ◽  
Ting Ma ◽  
Yan He ◽  
Meng Shen ◽  
...  

AbstractBladder cancer is one of the most common cancers worldwide. The immune response and immune cell infiltration play crucial roles in tumour progression. Immunotherapy has delivered breakthrough achievements in the past decade in bladder cancer. Differentially expressed genes and immune-related genes (DEIRGs) were identified by using the edgeR package. Gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed for functional enrichment analysis of DEIRGs. Survival-associated IRGs were identified by univariate Cox regression analysis. A prognostic model was established by univariate COX regression analysis, and verified by a validation prognostic model based on the GEO database. Patients were divided into high-risk and low-risk groups based on the median risk score value for immune cell infiltration and clinicopathological analyses. A regulatory network of survival-associated IRGs and potential transcription factors was constructed to investigate the potential regulatory mechanisms of survival-associated IRGs. Nomogram and ROC curve to verify the accuracy of the model. Quantitative real-time PCR was performed to validate the expression of relevant key genes in the prognostic model. A total of 259 differentially expressed IRGs were identified in the present study. KEGG pathway analysis of IRGs showed that the “cytokine-cytokine receptor interaction” pathway was the most significantly enriched pathway. Thirteen survival-associated IRGs were selected to establish a prognostic index for bladder cancer. In both TCGA prognostic model and GEO validation model, patients with high riskscore had worse prognosis compared to low riskscore group. A high infiltration level of macrophages was observed in high-risk patients. OGN, ELN, ANXA6, ILK and TGFB3 were identified as hub survival-associated IRGs in the network. EBF1, WWTR1, GATA6, MYH11, and MEF2C were involved in the transcriptional regulation of these survival-associated hub IRGs. The present study identified several survival-associated IRGs of clinical significance and established a prognostic index for bladder cancer outcome evaluation for the first time.


2020 ◽  
Vol 7 ◽  
Author(s):  
Shuting Wen ◽  
Long He ◽  
Zhuotai Zhong ◽  
Hong Mi ◽  
Fengbin Liu

BackgroundColorectal cancer (CRC) is a common malignant tumor of the digestive tract with a high mortality rate. Growing evidence demonstrates that immune-related genes play a prominent role in the occurrence and development of CRC. The aim of this study was to investigate the prognostic value of immune-related genes in CRC.MethodsGene expression profiles and clinical data of 568 CRC and 44 non-tumorous tissues were obtained from The Cancer Genome Atlas (TCGA) database. First, we performed a differentially expressed gene (DEG) analysis and univariate Cox regression analysis to determine the DEGs associated with overall survival. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were subsequently performed for prognostic immune-related genes. Then, a multivariate Cox regression analysis was performed to establish the immune prognostic model and identify the independent prognostic factors of CRC. Next, in vitro experiments were done to further validate the model. Finally, we analyzed the correlation among immune-related genes, clinical traits, and immune cell infiltration.ResultsIn total, 3,702 DEGs were obtained, and 338 prognostic immune-related genes were identified. Among them, 45 genes were significantly correlated with the prognosis of CRC patients. A TF-mediated network was set up to explore its internal mechanism. GO and KEGG analyses further illustrated that these genes were enriched in immune-and inflammatory-related pathways. Then, a prognostic prediction model composed of eight immune-related genes (SLC10A2, UTS2, FGF2, UCN, IL1RL2, ESM1, ADIPOQ, and VIP) was constructed. The AUC of the ROC curve for 1, 3, 5, and 10 years overall survival (OS) was 0.751, 0.707, 0.680, and 0.729, respectively. The survival analysis suggested that the OS of the high-risk group was significantly poorer than that of the low-risk group. Meanwhile, in vitro assays revealed that ESM1 and SLC10A2 exert opposing roles in colon cancer cell proliferation, validating the accuracy of the model. The correlation analysis indicated that immune cell infiltration was positively related to the model.ConclusionThis study screened prognosis-related immune genes and developed a prognostic prediction model of CRC. These findings may help provide potential novel prognostic biomarkers and therapeutic targets for CRC. At the same time, the understanding of the CRC immune microenvironment status was deepened.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Lin Chen ◽  
Yuxiang Dong ◽  
Yitong Pan ◽  
Yuhan Zhang ◽  
Ping Liu ◽  
...  

Abstract Background Breast cancer is one of the main malignant tumors that threaten the lives of women, which has received more and more clinical attention worldwide. There are increasing evidences showing that the immune micro-environment of breast cancer (BC) seriously affects the clinical outcome. This study aims to explore the role of tumor immune genes in the prognosis of BC patients and construct an immune-related genes prognostic index. Methods The list of 2498 immune genes was obtained from ImmPort database. In addition, gene expression data and clinical characteristics data of BC patients were also obtained from the TCGA database. The prognostic correlation of the differential genes was analyzed through Survival package. Cox regression analysis was performed to analyze the prognostic effect of immune genes. According to the regression coefficients of prognostic immune genes in regression analysis, an immune risk scores model was established. Gene set enrichment analysis (GSEA) was performed to probe the biological correlation of immune gene scores. P < 0.05 was considered to be statistically significant. Results In total, 556 immune genes were differentially expressed between normal tissues and BC tissues (p < 0. 05). According to the univariate cox regression analysis, a total of 66 immune genes were statistically significant for survival risk, of which 30 were associated with overall survival (P < 0.05). Finally, a 15 immune genes risk scores model was established. All patients were divided into high- and low-groups. KM survival analysis revealed that high immune risk scores represented worse survival (p < 0.001). ROC curve indicated that the immune genes risk scores model had a good reliability in predicting prognosis (5-year OS, AUC = 0.752). The established risk model showed splendid AUC value in the validation dataset (3-year over survival (OS) AUC = 0.685, 5-year OS AUC = 0.717, P = 0.00048). Moreover, the immune risk signature was proved to be an independent prognostic factor for BC patients. Finally, it was found that 15 immune genes and risk scores had significant clinical correlations, and were involved in a variety of carcinogenic pathways. Conclusion In conclusion, our study provides a new perspective for the expression of immune genes in BC. The constructed model has potential value for the prognostic prediction of BC patients and may provide some references for the clinical precision immunotherapy of patients.


2021 ◽  
Vol 15 ◽  
pp. 117955492110241
Author(s):  
Hongkai Zhuang ◽  
Zixuan Zhou ◽  
Zuyi Ma ◽  
Shanzhou Huang ◽  
Yuanfeng Gong ◽  
...  

Background: The prognosis of patients with pancreatic ductal adenocarcinoma (PDAC) of pancreatic head remains poor, even after potentially curative R0 resection. The aim of this study was to develop an accurate model to predict patients’ prognosis for PDAC of pancreatic head following pancreaticoduodenectomy. Methods: We retrospectively reviewed 112 patients with PDAC of pancreatic head after pancreaticoduodenectomy in Guangdong Provincial People’s Hospital between 2014 and 2018. Results: Five prognostic factors were identified using univariate Cox regression analysis, including age, histologic grade, American Joint Committee on Cancer (AJCC) Stage 8th, total bilirubin (TBIL), CA19-9. Using all subset analysis and multivariate Cox regression analysis, we developed a nomogram consisted of age, AJCC Stage 8th, perineural invasion, TBIL, and CA19-9, which had higher C-indexes for OS (0.73) and RFS (0.69) compared with AJCC Stage 8th alone (OS: 0.66; RFS: 0.67). The area under the curve (AUC) values of the receiver operating characteristic (ROC) curve for the nomogram for OS and RFS were significantly higher than other single parameter, which are AJCC Stage 8th, age, perineural invasion, TBIL, and CA19-9. Importantly, our nomogram displayed higher C-index for OS than previous reported models, indicating a better predictive value of our model. Conclusions: A simple and practical nomogram for patient prognosis in PDAC of pancreatic head following pancreaticoduodenectomy was established, which shows satisfactory predictive efficacy and deserves further evaluation in the future.


Author(s):  
Yongmei Wang ◽  
Guimin Zhang ◽  
Ruixian Wang

Background: This study aims to explore the prognostic values of CT83 and CT83-related genes in lung adenocarcinoma (LUAD). Methods: We downloaded the mRNA profiles of 513 LUAD patients (RNA sequencing data) and 246 NSCLC patients (Affymetrix Human Genome U133 Plus 2.0 Array) from TCGA and GEO databases. According to the median expression of CT83, the TCGA samples were divided into high and low expression groups, and differential expression analysis between them was performed. Functional enrichment analysis of differential expression genes (DEGs) was conducted. Univariate Cox regression analysis and LASSO Cox regression analysis were performed to screen the optimal prognostic DEGs. Then we established the prognostic model. A Nomogram model was constructed to predict the overall survival (OS) probability of LUAD patients. Results: CT83 expression was significantly correlated to the prognosis of LUAD patients. A total of 59 DEGs were identified, and a predictive model was constructed based on six optimal CT83-related DEGs, including CPS1, RHOV, TNNT1, FAM83A, IGF2BP1, and GRIN2A, could effectively predict the prognosis of LUAD patients. The nomogram could reliably predict the OS of LUAD patients. Moreover, the six important immune checkpoints (CTLA4, PD1, IDO1, TDO2, LAG3, and TIGIT) were closely correlated with the Risk Score, which was also differentially expressed between the LUAD samples with high and low-Risk Scores, suggesting that the poor prognosis of LUAD patients with high-Risk Score might be due to the immunosuppressive microenvironments. Conclusion: A prognostic model based on six optimal CT83 related genes could effectively predict the prognosis of LUAD patients.


2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Dan Chen ◽  
Xiaoting Li ◽  
Hui Li ◽  
Kai Wang ◽  
Xianghua Tian

Background. As the most common hepatic malignancy, hepatocellular carcinoma (HCC) has a high incidence; therefore, in this paper, the immune-related genes were sought as biomarkers in liver cancer. Methods. In this study, a differential expression analysis of lncRNA and mRNA in The Cancer Genome Atlas (TCGA) dataset between the HCC group and the normal control group was performed. Enrichment analysis was used to screen immune-related differentially expressed genes. Cox regression analysis and survival analysis were used to determine prognostic genes of HCC, whose expression was detected by molecular experiments. Finally, important immune cells were identified by immune cell infiltration and detected by flow cytometry. Results. Compared with the normal group, 1613 differentially expressed mRNAs (DEmRs) and 1237 differentially expressed lncRNAs (DElncRs) were found in HCC. Among them, 143 immune-related DEmRs and 39 immune-related DElncRs were screened out. These genes were mainly related to MAPK cascade, PI3K-AKT signaling pathway, and TGF-beta. Through Cox regression analysis and survival analysis, MMP9, SPP1, HAGLR, LINC02202, and RP11-598F7.3 were finally determined as the potential diagnostic biomarkers for HCC. The gene expression was verified by RT-qPCR and western blot. In addition, CD4 + memory resting T cells and CD8 + T cells were identified as protective factors for overall survival of HCC, and they were found highly expressed in HCC through flow cytometry. Conclusion. The study explored the dysregulation mechanism and potential biomarkers of immune-related genes and further identified the influence of immune cells on the prognosis of HCC, providing a theoretical basis for the prognosis prediction and immunotherapy in HCC patients.


Sign in / Sign up

Export Citation Format

Share Document