scholarly journals Prognostic Immune-Related Analysis Based on Differentially Expressed Genes in Left- and Right-Sided Colon Adenocarcinoma

2021 ◽  
Vol 11 ◽  
Author(s):  
Jun-Nan Guo ◽  
Ming-Qi Li ◽  
Shen-Hui Deng ◽  
Chen Chen ◽  
Yin Ni ◽  
...  

BackgroundColon adenocarcinoma (COAD) can be divided into left-sided and right-sided COAD (LCCs and RCCs, respectively). They have unique characteristics in various biological aspects, particularly immune invasion and prognosis. The purpose of our study was to develop a prognostic risk scoring model (PRSM) based on differentially expressed immune-related genes (IRGs) between LCCs and RCCs, therefore the prognostic key IRGs could be identified.MethodsThe gene sets and clinical information of COAD patients were derived from TCGA and GEO databases. The comparison of differentially expressed genes (DEGs) of LCCs and RCCs were conducted with appliance of “Limma” analysis. The establishment about co-expression modules of DEGs related with immune score was conducted by weighted gene co-expression network analysis (WGCNA). Furthermore, we screened the module genes and completed construction of gene pairs. The analysis of the prognosis and the establishment of PRSM were performed with univariate- and lasso-Cox regression. We employed the PRSM in the model group and verification group for the purpose of risk group assignment and PRSM accuracy verification. Finally, the identification of the prognostic key IRGs was guaranteed by the adoption of functional enrichment, “DisNor” and protein-protein interaction (PPI).ResultsA total of 215 genes were screened out by differential expression analysis and WGCNA. A PRSM with 16 immune-related gene pairs (IRGPs) was established upon the genes pairing. Furthermore, we confirmed that the risk score was an independent factor for survival by univariate- and multivariate-Cox regression. The prognosis of high-risk group in model group (P < 0.001) and validation group (P = 0.014) was significantly worse than that in low-risk group. Treg cells (P < 0.001) and macrophage M0 (P = 0.015) were highly expressed in the high-risk group. The functional analysis indicated that there was significant up-regulation with regard of lymphocyte and cytokine related terms in low-risk group. Finally, we identified five prognostic key IRGs associated with better prognosis through PPI and prognostic analysis, including IL2RB, TRIM22, CIITA, CXCL13, and CXCR6.ConclusionThrough the analysis and screening of the DEGs between LCCs and RCCs, we constructed a PRSM which could predicate prognosis of LCCs and RCCs, and five prognostic key IRGs were identified as well. Therefore, the basis for identifying the benefits of immunotherapy and immunomodulatory was built.

2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Dakui Luo ◽  
Zezhi Shan ◽  
Qi Liu ◽  
Sanjun Cai ◽  
Qingguo Li ◽  
...  

A metabolic disorder is considered one of the hallmarks of cancer. Multiple differentially expressed metabolic genes have been identified in colon cancer (CC), and their biological functions and prognostic values have been well explored. The purpose of the present study was to establish a metabolic signature to optimize the prognostic prediction in CC. The related data were downloaded from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx) database, and Gene Expression Omnibus (GEO) combined with GSE39582 set, GSE17538 set, GSE33113 set, and GSE37892 set. The differentially expressed metabolic genes were selected for univariate Cox regression and lasso Cox regression analysis using TCGA and GTEx datasets. Finally, a seventeen-gene metabolic signature was developed to divide patients into a high-risk group and a low-risk group. Patients in the high-risk group presented poorer prognosis compared to the low-risk group in both TCGA and GEO datasets. Moreover, gene set enrichment analyses demonstrated multiple significantly enriched metabolism-related pathways. To sum up, our study described a novel seventeen-gene metabolic signature for prognostic prediction of colon cancer.


Author(s):  
Yue Li ◽  
Ruoyi Shen ◽  
Anqi Wang ◽  
Jian Zhao ◽  
Jieqi Zhou ◽  
...  

BackgroundLung adenocarcinoma (LUAD) originates mainly from the mucous epithelium and glandular epithelium of the bronchi. It is the most common pathologic subtype of non-small cell lung cancer (NSCLC). At present, there is still a lack of clear criteria to predict the efficacy of immunotherapy. The 5-year survival rate for LUAD patients remains low.MethodsAll data were downloaded from The Cancer Genome Atlas (TCGA) database. We used Gene Set Enrichment Analysis (GSEA) database to obtain immune-related mRNAs. Immune-related lncRNAs were acquired by using the correlation test of the immune-related genes with R version 3.6.3 (Pearson correlation coefficient cor = 0.5, P < 0.05). The TCGA-LUAD dataset was divided into the testing set and the training set randomly. Based on the training set to perform univariate and multivariate Cox regression analyses, we screened prognostic immune-related lncRNAs and given a risk score to each sample. Samples were divided into the high-risk group and the low-risk group according to the median risk score. By the combination of Kaplan–Meier (KM) survival curve, the receiver operating characteristic (ROC) (AUC) curve, the independent risk factor analysis, and the clinical data of the samples, we assessed the accuracy of the risk model. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the differentially expressed mRNAs between the high-risk group and the low-risk group. The differentially expressed genes related to immune response between two risk groups were analyzed to evaluate the role of the model in predicting the efficacy and effects of immunotherapy. In order to explain the internal mechanism of the risk model in predicting the efficacy of immunotherapy, we analyzed the differentially expressed genes related to epithelial-mesenchymal transition (EMT) between two risk groups. We extracted RNA from normal bronchial epithelial cell and LUAD cells and verified the expression level of lncRNAs in the risk model by a quantitative real-time polymerase chain reaction (qRT-PCR) test. We compared our risk model with other published prognostic signatures with data from an independent cohort. We transfected LUAD cell with siRNA-LINC0253. Western blot analysis was performed to observed change of EMT-related marker in protein level.ResultsThrough univariate Cox regression analysis, 24 immune-related lncRNAs were found to be strongly associated with the survival of the TCGA-LUAD dataset. Utilizing multivariate Cox regression analysis, 10 lncRNAs were selected to establish the risk model. The K-M survival curves and the ROC (AUC) curves proved that the risk model has a fine predictive effect. The GO enrichment analysis indicated that the effect of the differentially expressed genes between high-risk and low-risk groups is mainly involved in immune response and intercellular interaction. The KEGG enrichment analysis indicated that the differentially expressed genes between high-risk and low-risk groups are mainly involved in endocytosis and the MAPK signaling pathway. The expression of genes related to the efficacy of immunotherapy was significantly different between the two groups. A qRT-PCR test verified the expression level of lncRNAs in LUAD cells in the risk model. The AUC of ROC of 5 years in the independent validation dataset showed that this model had superior accuracy. Western blot analysis verified the change of EMT-related marker in protein level.ConclusionThe immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.


2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Ying Ye ◽  
Qinjin Dai ◽  
Hongbo Qi

AbstractOvarian cancer (OC) is a highly malignant gynaecological tumour that has a very poor prognosis. Pyroptosis has been demonstrated in recent years to be an inflammatory form of programmed cell death. However, the expression of pyroptosis-related genes in OC and their correlations with prognosis remain unclear. In this study, we identified 31 pyroptosis regulators that were differentially expressed between OC and normal ovarian tissues. Based on these differentially expressed genes (DEGs), all OC cases could be divided into two subtypes. The prognostic value of each pyroptosis-related gene for survival was evaluated to construct a multigene signature using The Cancer Genome Atlas (TCGA) cohort. By applying the least absolute shrinkage and selection operator (LASSO) Cox regression method, a 7-gene signature was built and classified all OC patients in the TCGA cohort into a low- or high-risk group. OC patients in the low-risk group showed significantly higher survival possibilities than those in the high-risk group (P < 0.001). Utilizing the median risk score from the TCGA cohort, OC patients from a Gene Expression Omnibus (GEO) cohort were divided into two risk subgroups, and the low-risk group had increased overall survival (OS) time (P = 0.014). Combined with the clinical characteristics, the risk score was found to be an independent factor for predicting the OS of OC patients. Gene ontology (GO) and Kyoto Encylopedia of Genes and Genomes (KEGG) analyses indicated that immune-related genes were enriched and that the immune status was decreased in the high-risk group. In conclusion, pyroptosis-related genes play important roles in tumour immunity and can be used to predict the prognosis of OCs.


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.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e8128 ◽  
Author(s):  
Cheng Yue ◽  
Hongtao Ma ◽  
Yubai Zhou

Background Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient’s outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.


2020 ◽  
Author(s):  
Jihang Luo ◽  
Puyu Liu ◽  
Leibo Wang ◽  
Yi Huang ◽  
Yuanyan Wang ◽  
...  

Abstract Background Colon cancer is the most common type of gastrointestinal cancer and has high morbidity and mortality. Colon adenocarcinoma(COAD) is the main pathological type of colon cancer. There is a lot of evidence describing the correlation between the prognosis of COAD and the immune system. The objective of the current study was the development of a robust prognostic immune-related gene pairs (IRGPs) model for estimating overall survival of COAD. Methods The gene expression profiles and clinical information of patients with colon adenocarcinoma come from TCGA and GEO databases and are divided into training and validation cohorts. Immune genes were selected which show significantly association with prognosis. Results Among 1647 immune genes, a 17 IRGPs model was built which was significantly associated with OS in the training cohort. In the training and validation data set, the IRGPs model divided patients into high-risk groups and low-risk groups, and the prognosis of the high-risk group was significantly worse( P <0.001). Univariate and multivariate Cox proportional hazard analysis confirmed the feasibility of this model. Functional analysis confirmed that multiple tumor progression and stem cell growth-related pathways in high-risk groups were up-regulated. T cells regulatory and Macrophage M0 were significantly highly expressed in the high-risk group. Conclusion We successfully constructed an IRGPs model that can predict the prognosis of COAD, which provides new insights into the treatment strategy of COAD.


2020 ◽  
Vol 13 (1) ◽  
Author(s):  
Yinglian Pan ◽  
Li Ping Jia ◽  
Yuzhu Liu ◽  
Yiyu Han ◽  
Qian Li ◽  
...  

Abstract Background In this study we aimed to identify a prognostic signature in BRCA1/2 mutations to predict disease progression and the efficiency of chemotherapy ovarian cancer (OV), the second most common cause of death from gynecologic cancer in women worldwide. Methods Univariate Cox proportional-hazards and multivariate Cox regression analyses were used to identifying prognostic factors from data obtained from The Cancer Genome Atlas (TCGA) database. The area under the curve of the receiver operating characteristic curve was assessed, and the sensitivity and specificity of the prediction model were determined. Results A signature consisting of two long noncoding RNAs(lncRNAs), Z98885.2 and AC011601.1, was selected as the basis for classifying patients into high and low-risk groups (median survival: 7.2 years vs. 2.3 years). The three-year overall survival (OS) rates for the high- and low-risk group were approximately 38 and 100%, respectively. Chemotherapy treatment survival rates indicated that the high-risk group had significantly lower OS rates with adjuvant chemotherapy than the low-risk group. The one-, three-, and five-year OS were 100, 40, and 15% respectively in the high-risk group. The survival rate of the high-risk group declined rapidly after 2 years of OV chemotherapy treatment. Multivariate Cox regression associated with other traditional clinical factors showed that the 2-lncRNA model could be used as an independent OV prognostic factor. Analyses of data from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) indicated that these signatures are pivotal to cancer development. Conclusion In conclusion, Z98885.2 and AC011601.1 comprise a novel prognostic signature for OV patients with BRCA1/2 mutations, and can be used to predict prognosis and the efficiency of chemotherapy.


2020 ◽  
Author(s):  
Li Liu ◽  
She Tian ◽  
Zhu Li ◽  
Yongjun Gong ◽  
Hao Zhang

Abstract Background : Hepatocellular carcinoma (HCC) is one of the most common clinical malignant tumors, resulting in high mortality and poor prognosis. Studies have found that LncRNA plays an important role in the onset, metastasis and recurrence of hepatocellular carcinoma. The immune system plays a vital role in the development, progression, metastasis and recurrence of cancer. Therefore, immune-related lncRNA can be used as a novel biomarker to predict the prognosis of hepatocellular carcinoma. Methods : The transcriptome data and clinical data of HCC patients were obtained by using The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA‑LIHC), and immune-related genes were extracted from the Molecular Signatures Database (IMMUNE RESPONSE M19817 and IMMUNE SYSTEM PROCESS M13664). By constructing the co-expression network and Cox regression analysis, 13 immune-lncRNAs was identified to predict the prognosis of HCC patients. Patients were divided into high risk group and low risk group by using the risk score formula, and the difference in overall survival (OS) between the two groups was reflected by Kaplan-Meier survival curve. The time - dependent receiver operating characteristics (ROC) analysis and principal component analysis (PCA) were used to evaluate 13 immune -lncRNAs signature. Results : Through TCGA - LIHC extracted from 343 cases of patients with hepatocellular carcinoma RNA - Seq data and clinical data, 331 immune-related genes were extracted from the Molecular Signatures Database , co-expression networks and Cox regression analysis were constructed, 13 immune-lncRNAs signature was identified as biomarkers to predict the prognosis of patients. At the same time using the risk score median divided the patients into high risk and low risk groups, and through the Kaplan-Meier survival curve analysis found that high-risk group of patients' overall survival (OS) less low risk group of patients. The AUC value of the ROC curve is 0.828, and principal component analysis (PCA) results showed that patients could be clearly divided into two parts by immune-lncRNAs, which provided evidence for the use of 13 immune-lncRNAs signature as prognostic markers. Conclusion : Our study identified 13 immune-lncRNAs signature that can effectively predict the prognosis of HCC patients, which may be a new prognostic indicator for predicting clinical outcomes.


2021 ◽  
Vol 12 ◽  
Author(s):  
Qian Zhang ◽  
Liping Lv ◽  
Ping Ma ◽  
Yangyang Zhang ◽  
Jiang Deng ◽  
...  

BackgroundPancreatic adenocarcinoma (PAAD) spreads quickly and has a poor prognosis. Autophagy research on PAAD could reveal new biomarkers and targets for diagnosis and treatment.MethodsAutophagy-related genes were translated into autophagy-related gene pairs, and univariate Cox regression was performed to obtain overall survival (OS)-related IRGPs (P&lt;0.001). LASSO Cox regression analyses were performed to construct an autophagy-related gene pair (ARGP) model for predicting OS. The Cancer Genome Atlas (TCGA)-PAAD cohort was set as the training group for model construction. The model predictive value was validated in multiple external datasets. Receiver operating characteristic (ROC) curves were used to evaluate model performance. Tumor microenvironments and immune infiltration were compared between low- and high-risk groups with ESTIMATE and CIBERSORT. Differentially expressed genes (DEGs) between the groups were further analyzed by Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses and used to identify potential small-molecule compounds in L1000FWD.ResultsRisk scores were calculated as follows: ATG4B|CHMP4C×(-0.31) + CHMP2B|MAP1LC3B×(0.30) + CHMP6|RIPK2 ×(-0.33) + LRSAM1|TRIM5×(-0.26) + MAP1LC3A|PAFAH1B2×(-0.15) + MAP1LC3A|TRIM21×(-0.08) + MET|MFN2×(0.38) + MET|MTDH×(0.47) + RASIP1|TRIM5×(-0.23) + RB1CC1|TPCN1×(0.22). OS was significantly shorter in the high-risk group than the low-risk group in each PAAD cohort. The ESTIMATE analysis showed no difference in stromal scores but a significant difference in immune scores (p=0.0045) and ESTIMATE scores (p=0.014) between the groups. CIBERSORT analysis showed higher naive B cell, Treg cell, CD8 T cell, and plasma cell levels in the low-risk group and higher M1 and M2 macrophage levels in the high-risk group. In addition, the results showed that naive B cells (r=-0.32, p&lt;0.001), Treg cells (r=-0.31, p&lt;0.001), CD8 T cells (r=-0.24, p=0.0092), and plasma cells (r=-0.2, p&lt;0.026) were statistically correlated with the ARGP risk score. The top 3 enriched GO-BPs were signal release, regulation of transsynaptic signaling, and modulation of chemical synaptic transmission, and the top 3 enriched KEGG pathways were the insulin secretion, dopaminergic synapse, and NF-kappa B signaling pathways. Several potential small-molecule compounds targeting ARGs were also identified.ConclusionOur results demonstrate that the ARGP-based model may be a promising prognostic indicator for identifying drug targets in patients with PAAD.


Author(s):  
Peng Gu ◽  
Lei Zhang ◽  
Ruitao Wang ◽  
Wentao Ding ◽  
Wei Wang ◽  
...  

Background: Female breast cancer is currently the most frequently diagnosed cancer in the world. This study aimed to develop and validate a novel hypoxia-related long noncoding RNA (HRL) prognostic model for predicting the overall survival (OS) of patients with breast cancer.Methods: The gene expression profiles were downloaded from The Cancer Genome Atlas (TCGA) database. A total of 200 hypoxia-related mRNAs were obtained from the Molecular Signatures Database. The co-expression analysis between differentially expressed hypoxia-related mRNAs and lncRNAs based on Spearman’s rank correlation was performed to screen out 166 HRLs. Based on univariate Cox regression and least absolute shrinkage and selection operator Cox regression analysis in the training set, we filtered out 12 optimal prognostic hypoxia-related lncRNAs (PHRLs) to develop a prognostic model. Kaplan–Meier survival analysis, receiver operating characteristic curves, area under the curve, and univariate and multivariate Cox regression analyses were used to test the predictive ability of the risk model in the training, testing, and total sets.Results: A 12-HRL prognostic model was developed to predict the survival outcome of patients with breast cancer. Patients in the high-risk group had significantly shorter median OS, DFS (disease-free survival), and predicted lower chemosensitivity (paclitaxel, docetaxel) compared with those in the low-risk group. Also, the risk score based on the expression of the 12 HRLs acted as an independent prognostic factor. The immune cell infiltration analysis revealed that the immune scores of patients in the high-risk group were lower than those of the patients in the low-risk group. RT-qPCR assays were conducted to verify the expression of the 12 PHRLs in breast cancer tissues and cell lines.Conclusion: Our study uncovered dozens of potential prognostic biomarkers and therapeutic targets related to the hypoxia signaling pathway in breast cancer.


Sign in / Sign up

Export Citation Format

Share Document