scholarly journals MYD88 Is a Potential Prognostic Gene and Immune Signature of Tumor Microenvironment for Gliomas

2021 ◽  
Vol 11 ◽  
Author(s):  
Qinglong Guo ◽  
Xing Xiao ◽  
Jinsen Zhang

PurposeTo explore the profiles of immune and stromal components of the tumor microenvironment (TME) and their related key genes in gliomas.MethodsWe applied bioinformatic techniques to identify the core gene that participated in the regulation of the TME of the gliomas. And immunohistochemistry staining was used to calculate the gene expressions in clinical cases.ResultsThe CIBERSORT and ESTIMATE were used to figure out the composition of TME in 698 glioma cases from The Cancer Genome Atlas (TCGA) database. Differential expression analysis identified 2103 genes between the high and the low-score group. Then the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, univariate Cox regression analysis, and protein–protein interaction (PPI) network construction were conducted based on these genes. MYD88 was identified as the key gene by the combination univariate Cox and PPI analysis. Furthermore, MYD88 expression was significantly associated with the overall survival and WHO grade of glioma patients. The genes in the high-expression MYD88 group were mainly in immune-related pathways in the Gene Set Enrichment Analysis (GSEA). We found that macrophage M2 accounted for the largest portion with an average of 27.6% in the glioma TIICs and was associated with high expression of MYD88. The results were verified in CGGA database and clinical cases in our hospital. Furthermore, we also found the MYD88 expression was higher in IDH1 wild types. The methylation rate was lower in high grade gliomas.ConclusionMYD88 had predictive prognostic value in glioma patients by influencing TIICs dysregulation especially the M2-type macrophages.

2021 ◽  
Author(s):  
Yucheng Qian ◽  
Lina Zhang ◽  
Jihang Wen ◽  
Yanxia Mei ◽  
Jingsun Wei ◽  
...  

Abstract ColorColorectal cancer is one of the most common cancer worldwide. Recently, tumor microenvironment (TME), especially its remoulding , is thought to control the colorectal cancer genesis and progression. In this study, we use ESTIMATEscore to make out the proportion of immune and stromal components in colorectal adenocarcinoma (CRA) samples from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were found by COX regression analysis and protein-protein interaction (PPI) network, among which TGFβ1 was supposed to be a prognosis factor and tumor environment indicator. Continuous analysis showed that TGFβ1 expression is positively correlated with lymph node metastasis (N stage) but negatively correlated with survival. Gene Set Enrichment Analysis (GSEA) revealed that the genes of the high-expression TGFβ1 group were mainly enriched in immune-related activities. Cluster analysis divided the samples into 2 subgroups. 24 HLA-related genes and 8 immune checkpoint genes were found upregulated in the high immunity group as well as TGFβ1, which suggests the possibility of novel therapies targeting immune checkpoints combined with TGFβ1. Tumor-infiltrating immune cell (TIC) profile of CRA patients was described by CIBERSORT analysis. Further analysis showed that the infiltration of Tregs and Neutrophils were positively correlated with TGFβ1 high expression. Then 3 TGFβ1-related genes were picked out to construct a prognostic model, which matches the survival data well.


2021 ◽  
Author(s):  
Chengran Xu ◽  
Jinhai Huang ◽  
Yi Yang ◽  
Lun Li ◽  
Guangyu Li

Abstract Background: The homeobox gene 5 (HOXB5) encodes a transcription factor that regulates the central nervous system embryonic development. Of note, its expression pattern and prognostic role in glioma remain unelucidated. This study aimed to identify the relationship between HOXB5 and glioma by investigating the HOXB5 expression data from the The Cancer Genome Atlas (TCGA) and The Genotype Tissue Expression (GTEx) databases and validating the obtained data using the Chinese Glioma Genome Atlas (CGGA) database. Kaplan-Meier and univariate cox regression analyses were performed to assess the prognostic value of HOXB5. The key functions and signaling pathways of HOXB5 were analyzed using GSEA and GSVA. Immune infiltration was calculated using Microenvironment Cell Populations-counter (MCP-counter), single-sample Gene Set Enrichment Analysis (ssGSEA), and ESTIMATE algorithms.Result: HOXB5 expression was elevated in glioma tissues. The increased levels of HOXB5 were significantly correlated with a higher WHO grade and aggressive cancer phenotypes. HOXB5 overexpression represented a risk factor that was associated with shorter overall survival (OS) while exhibiting a moderate forecast efficiency in most clinical subgroups. These results were validated using the CGGA and Rembrandt datasets. Furthermore, the functional analysis showed enrichment of angiogenesis, the IL6/JAK-STAT3 pathway, and inflammatory response in the tissues that showed high expression of HOXB5. Lastly, the high expression of HOXB5 was associated with enrichment of Tregs and MDSCs, and HOXB5 expression was shown to play a role in several immune checkpoint genes.Conclusions: HOXB5 may serve as a predictive factor of glioma malignancy and prognostic status and represents potential as a molecular treatment candidate.


2021 ◽  
Vol 11 ◽  
Author(s):  
Long-hao Chen ◽  
Jin-Fu Liu ◽  
Yan- Lu ◽  
Xin-yu He ◽  
Chi- Zhang ◽  
...  

The tumor microenvironment (TME) has important effects on the tumorigenesis and development of osteosarcoma (OS). However, the dynamic mechanism regulating TME immune and matrix components remains unclear. In this study, we collected quantitative data on the gene expression of 88 OS samples from The Cancer Genome Atlas (TCGA) database and downloaded relevant clinical cases of OS from the TARGET database. The proportions of tumor-infiltrating immune cells (TICs) and the numbers of immune and matrix components were determined by CIBERSORT and ESTIMATE calculation methods. Protein-protein interaction (PPI) network construction and Cox regression analysis were conducted to analyze differentially expressed genes (DEGs). The complement components C1qA, C1qB and C1qC were then determined to be predictive factors through univariate Cox analysis and PPI cross analysis. Further analysis found that the levels of C1qA, C1qB and C1qC expression were positively linked to OS patient survival time and negatively correlated with the clinicopathological feature percent necrosis at definitive surgery. The results of gene set enrichment analysis (GSEA) demonstrated that genes related to immune functions were significantly enriched in the high C1qA, C1qB and C1qC expression groups. Proportion analysis of TICs by CIBERSORT showed that the levels of C1qA, C1qB and C1qC expression were positively related to M1 and M2 macrophages and CD8+ cells and negatively correlated with M0 macrophages. These results further support the influence of the levels of C1qA, C1qB and C1qC expression on the immune activity of the TME. Therefore, C1qA, C1qB and C1qC may be potential indicators of remodeling in the OS TME, which is helpful to predict the prognosis of patients with OS and provide new ideas for immunotherapy for OS.


2020 ◽  
Vol 40 (8) ◽  
Author(s):  
Sihan Chen ◽  
Guodong Cao ◽  
Wei Wu ◽  
Yida Lu ◽  
Xiaobo He ◽  
...  

Abstract Colon adenocarcinoma (COAD) is a malignant gastrointestinal tumor, often occurring in the left colon, which is regulated by glycolysis-related processes. In past studies, multiple genes that influence the prognosis for survival have been discovered through bioinformatics analysis. However, the prediction of disease prognosis using a single gene is not an accurate method. In the present study, a mechanistic model was established to achieve better prediction for the prognosis of COAD. COAD-related data downloaded from The Cancer Genome Atlas (TCGA) were correlated with the glycolysis process using gene set enrichment analysis (GSEA) to determine the glycolysis-related genes that regulate COAD. Using COX regression analysis, glycolysis-related genes associated with the prognosis of COAD were identified, and the genes screened to establish a predictive model. The risk scores of this model were correlated with relevant clinical data to obtain a connection diagram between the model and survival rate, tumor characteristic data, etc. Finally, genes in the model were correlated with cells in the tumor microenvironment, finding that they affected specific immune cells in the model. Seven genes related to glycolysis were identified (PPARGC1A, DLAT, 6PC2, P4HA1, STC2, ANKZF1, and GPC1), which affect the prognosis of patients with COAD and constitute the model for prediction of survival of COAD patients.


2021 ◽  
Vol 12 ◽  
Author(s):  
Guomin Wu ◽  
Qihao Wang ◽  
Ting Zhu ◽  
Linhai Fu ◽  
Zhupeng Li ◽  
...  

This study aimed to establish a prognostic risk model for lung adenocarcinoma (LUAD). We firstly divided 535 LUAD samples in TCGA-LUAD into high-, medium-, and low-immune infiltration groups by consensus clustering analysis according to immunological competence assessment by single-sample gene set enrichment analysis (ssGSEA). Profile of long non-coding RNAs (lncRNAs) in normal samples and LUAD samples in TCGA was used for a differential expression analysis in the high- and low-immune infiltration groups. A total of 1,570 immune-related differential lncRNAs in LUAD were obtained by intersecting the above results. Afterward, univariate COX regression analysis and multivariate stepwise COX regression analysis were conducted to screen prognosis-related lncRNAs, and an eight-immune-related-lncRNA prognostic signature was finally acquired (AL365181.2, AC012213.4, DRAIC, MRGPRG-AS1, AP002478.1, AC092168.2, FAM30A, and LINC02412). Kaplan–Meier analysis and ROC analysis indicated that the eight-lncRNA-based model was accurate to predict the prognosis of LUAD patients. Simultaneously, univariate COX regression analysis and multivariate COX regression analysis were undertaken on clinical features and risk scores. It was illustrated that the risk score was a prognostic factor independent from clinical features. Moreover, immune data of LUAD in the TIMER database were analyzed. The eight-immune-related-lncRNA prognostic signature was related to the infiltration of B cells, CD4+ T cells, and dendritic cells. GSEA enrichment analysis revealed significant differences in high- and low-risk groups in pathways like pentose phosphate pathway, ubiquitin mediated proteolysis, and P53 signaling pathway. This study helps to treat LUAD patients and explore molecules related to LUAD immune infiltration to deeply understand the specific mechanism.


2021 ◽  
Vol 2021 ◽  
pp. 1-8
Author(s):  
Honglan Guo ◽  
Qinqiao Fan

Background. We aimed to investigate the expression of the hyaluronan-mediated motility receptor (HMMR) gene in hepatocellular carcinoma (HCC) and nonneoplastic tissues and to investigate the diagnostic and prognostic value of HMMR. Method. With the reuse of the publicly available The Cancer Genome Atlas (TCGA) data, 374 HCC patients and 50 nonneoplastic tissues were used to investigate the diagnostic and prognostic values of HMMR genes by receiver operating characteristic (ROC) curve analysis and survival analysis. All patients were divided into low- and high-expression groups based on the median value of HMMR expression level. Univariate and multivariate Cox regression analysis were used to identify prognostic factors. Gene set enrichment analysis (GSEA) was performed to explore the potential mechanism of the HMMR genes involved in HCC. The diagnostic and prognostic values were further validated in an external cohort from the International Cancer Genome Consortium (ICGC). Results. HMMR mRNA expression was significantly elevated in HCC tissues compared with that in normal tissues from both TCGA and the ICGC cohorts (all P values <0.001). Increased HMMR expression was significantly associated with histologic grade, pathological stage, and survival status (all P values <0.05). The area under the ROC curve for HMMR expression in HCC and normal tissues was 0.969 (95% CI: 0.948–0.983) in the TCGA cohort and 0.956 (95% CI: 0.932–0.973) in the ICGC cohort. Patients with high HMMR expression had a poor prognosis than patients with low expression group in both cohorts (all P < 0.001 ). Univariate and multivariate analysis also showed that HMMR is an independent predictor factor associated with overall survival in both cohorts (all P values <0.001). GSEA showed that genes upregulated in the high-HMMR HCC subgroup were mainly significantly enriched in the cell cycle pathway, pathways in cancer, and P53 signaling pathway. Conclusion. HMMR is expressed at high levels in HCC. HMMR overexpression may be an unfavorable prognostic factor for HCC.


2021 ◽  
Vol 19 (1) ◽  
pp. 169-190
Author(s):  
Peiyuan Li ◽  
◽  
Gangjie Qiao ◽  
Jian Lu ◽  
Wenbin Ji ◽  
...  

<abstract> <p>Plasmacytoma variant translocation 1 (PVT1) is involved in multiple signaling pathways and plays an important regulatory role in a variety of malignant tumors. However, its role in the prognosis and immune invasion of bladder urothelial carcinoma (BLCA) remains unclear. This study investigated the expression of PVT1 in tumor tissue and its relationship with immune invasion, and determined its prognostic role in patients with BLCA. Patients were identified from the cancer genome atlas (TCGA). The enrichment pathway and function of PVT1 were explained by gene ontology (GO) term analysis, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA), and the degree of immune cell infiltration was quantified. Kaplan–Meier analysis and Cox regression were used to analyze the correlation between PVT1 and survival rate. PVT1-high BLCA patients had a lower 10-year disease-specific survival (DSS P &lt; 0.05) and overall survival (OS P &lt; 0.05). Multivariate Cox regression analysis showed that PVT1 (high vs. low) (P = 0.004) was an independent prognostic factor. A nomogram was used to predict the effect of PVT1 on the prognosis. PVT1 plays an important role in the progression and prognosis of BLCA and can be used as a medium biomarker to predict survival after cystectomy.</p> </abstract>


2021 ◽  
Vol 2021 ◽  
pp. 1-9
Author(s):  
Xiangfeng Jin ◽  
Haiqing Zhou ◽  
Jianfang Song ◽  
Hong Cui ◽  
Yiren Luo ◽  
...  

Background. The present study is aimed at evaluating the functional and clinical values of P3H4 in lung adenocarcinoma. Moreover, we also investigated the downstream pathways that P3H4 might participate in. Methods. The differential expression analysis was used to identify genes differentially expressed in lung adenocarcinoma tissues as compared with normal tissues. Survival analysis was used to test the association between P3H4 and survival time. Gene set enrichment analysis was conducted to explore the downstream pathways. CCK8 and transwell were employed to examine the impact of P3H4 on cell phenotypes. Results. P3H4 was highly upregulated in LUAD tissues at both RNA and protein levels. Moreover, the LUAD patients, who had high expression of P3H4, were also observed to have shorter disease-free survival and overall survival. These results demonstrated that P3H4 could be used as a prognostic biomarker for LUAD. Moreover, we also found that it was the copy number alterations (CNAs), not DNA methylation, that regulated the RNA expression of P3H4, indicating that its upregulation might be partially resulted from the CNAs. Furthermore, functional experiments revealed that the A549 and H1299 cells with siRNA treatment (siP3H4) exhibited significantly decreased cell proliferation after 24 hours, migratory ability, and invasiveness. Functionally, the upregulated proteins in the P3H4 high expression group were mainly enriched in tumor microenvironment-related pathways such as phagosome, focal adhesion, and ECM-receptor interaction and cancer-related pathways such as bladder cancer pathway, proteoglycans in cancer, and hippo signaling pathway. Conclusion. The present study systematically evaluated the functional and clinical values of P3H4 in LUAD, and explored the related biological pathways. P3H4 might promote LUAD progression through regulating tumor microenvironment-related pathways.


2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Jin Zhou ◽  
Zheming Liu ◽  
Huibo Zhang ◽  
Tianyu Lei ◽  
Jiahui Liu ◽  
...  

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.


2021 ◽  
Author(s):  
Lijun Ning ◽  
Yuqing Yan ◽  
Tianying Tong ◽  
Ziyun Gao ◽  
Zhe Cui ◽  
...  

Abstract Background: As tumor microenvironment (TME) play an indispensable role in tumorigenesis of colorectal cancer, this study performs a bunch of bioinformatics analysis to identify the indicator of the status of TME in Colorectal cancer (CRC). Results: In the presented study, we applied CIBERSORT and ESTIMATE computational methods to calculate the proportion of tumor-infiltrating immune cells (TICs) and the amount of immune and stromal components in 444 COAD-READ cases from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were analyzed by COX regression analysis and protein–protein interaction (PPI) network construction. Then, fatty acid-binding protein four ( FABP4 ) was determined as a predictive factor by the intersection analysis of univariate COX and PPI. Further analysis revealed that FABP4 expression was positively correlated with the clinical pathologic characteristics (clinical stage, distant metastasis) and negatively correlated with the survival of CRC patients. Gene Set Enrichment Analysis (GSEA) showed that the genes in the high-expression FABP4 group were mainly enriched in immune-related activities. In the low-expression FABP4 group, the genes were enriched in metabolic pathways. CIBERSORT analysis for the proportion of TICs revealed that NK cell, CD4 + T cells and CD8 + T cells were negatively correlated with FABP4 expression, suggesting that FABP4 might be a potential prognostic factor of CRC patients. Conclusion: Our study has developed a new biomarker (FABP4) that can predict the status of tumor microenvironment in Colorectal cancer. Keywords: FABP4, tumor microenvironment, ESTIMATE, CIBERSORT, colorectal cancer


Sign in / Sign up

Export Citation Format

Share Document