TCGA Database-Based Screening of Tumor Microenvironment Immunomodulators Related to Bladder Cancer Prognosis.

Bladder cancer (BC), as the most common malignant tumor of the urinary tract, has a complex biological behavior. Currently, there are still some limitations in the diagnosis and treatment of BC. Despite the great progress made in immunotherapy, there is still a lack of key genes for the diagnosis of BC. Therefore, it is particularly important to explore the differentially expressed genes (DEGs) and their effectiveness on prognosis of BC with different tumor microenvironment scores. The gene expression dataset of BC was downloaded from the Cancer Genome Atlas (TCGA) database. The correlation between clinicopathological characteristics of patients and scores of immune and stromal components was analyzed. Patients were divided into high and low score groups according to their tumor microenvironment score (Immune score, Stromal score, ESTIMATE score). DEGs between high and low score groups were identified using R software and then subjected to enrichment analyses to assess their potential biological functions and signaling pathways. The protein-protein interaction (PPI) network was constructed using the STRING database to further identify hub genes. The expression levels of hub genes in BC were verified by TCGA database. Subsequently, the hub genes were evaluated for overall survival (OS), disease-free survival (DFS), progression-free survival (PFS), and disease-specific survival (DSS), and corresponding forest plots were created. A total of 2346 DEGs were obtained, including 1120 up-regulated genes and 1226 down-regulated genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses found DEGs were mainly enriched in cell migration and immune-related pathways. Meanwhile, The PPI network finally yielded top 10 hub genes with predictive value, which included actin beta (ACTB), interleukin 6 (IL-6), Jun proto-oncogene (JUN), CD4 molecule (CD4), heat shock protein 90 alpha family class A member 1 (HSP90AA1), protein tyrosine phosphatase receptor type C (PTPRC), tumor protein p53 (TP53), SRC proto-oncogene (SRC), fibronectin 1 (FN1), and tumor necrosis factor (TNF). Among them, CD4, PTPRC, and SRC were potential protective factors for BC. The top 10 hub genes (ACTB, IL-6, JUN, CD4, HSP90AA1, PTPRC, TP53, SRC, FN1, TNF) obtained based on tumor microenvironment scores all had potential predictive value. Elevated expression of protective factors (CD4, PTPRC, and SRC) indicates better survival outcome of BC subjects. Further exploration of the molecular developmental mechanisms of these hub genes will help to develop novel personalized therapies and improve BC prognosis.

Ye JN ，Gan MP ，Wu XL ，Dong Z ，Mai HH ... -  《-》

CX3CR1 Acts as a Protective Biomarker in the Tumor Microenvironment of Colorectal Cancer.

The tumor microenvironment (TME) plays an important role in the pathogenesis of many cancers. We aimed to screen the TME-related hub genes of colorectal adenoma (CRAD) and identify possible prognostic biomarkers. The gene expression profiles and clinical data of 464 CRAD patients in The Cancer Genome Atlas (TCGA) database were downloaded. The Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm was performed to calculate the ImmuneScore, StromalScore, and EstimateScore. Thereafter, differentially expressed genes (DEGs) were screened. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analysis were performed to explore the roles of DEGs. Furthermore, univariate and multivariate Cox analyses were accomplished to identify independent prognostic factors of CRAD. CX3CR1 was selected as a hub gene, and the expression was confirmed in colorectal cancer (CRC) patients and cell lines. The correlations between CX3CR1 and tumor-infiltrating immune cells were estimated by Tumor IMmune Estimation Resource database (TIMER) and CIBERSORT analysis. Besides, we investigated the effects of coculture with THP-1-derived macrophages with HCT8 cells with low CX3CR1 expression on immune marker expression, cell viability, and migration. There were significant differences in the ImmuneScore and EstimateScore among different stages. Patients with low scores presented significantly lower lifetimes than those in the high-score group. Moreover, we recognized 1,578 intersection genes in ImmuneScore and StromalScore, and these genes were mainly enriched in numerous immune-related biological processes. CX3CR1 was found to be associated with immune cell infiltration levels, immune marker expression, and macrophage polarization. Simultaneous silencing of CX3CR1 and coculture with THP-1 cells further regulated macrophage polarization and promoted the cell proliferation and migration of CRC cells. CX3CR1 was decreased in CRAD tissues and cell lines and was related to T and N stages, tumor differentiation, and prognosis. Our results suggest that CX3CR1 contributes to the recruitment and regulation of immune-infiltrating cells and macrophage polarization in CRC and TAM-induced CRC progression. CX3CR1 may act as a prognostic biomarker in CRC.

Yue Y ，Zhang Q ，Sun Z 《Frontiers in Immunology》

Identification of Prognosis-Related Genes in Bladder Cancer Microenvironment across TCGA Database.

Bladder cancer (BCa) is a common urothelial malignancy. The Cancer Genome Atlas (TCGA) database allows for an opportunity to analyze the relationship between gene expression and clinical outcomes in bladder cancer patients. This study is aimed at identifying prognosis-related genes in the bladder cancer microenvironment. Immune scores and stromal scores were calculated by applying the ESTIMATE algorithm. We divided bladder cancer patients into high and low groups based on their immune/stromal scores. Then, differentially expressed genes (DEGs) were identified in bladder cancer patients based on the TCGA database. We evaluated the correlation between immune/stromal scores and clinical characteristics as well as prognosis. Finally, we validated identified genes associated with bladder cancer prognosis through a cohort study in the Gene Expression Omnibus (GEO) database. A higher stromal score was associated with female (vs. malep = 0.037), age > 65 (vs.age ≤ 65 p = 0.015), T3/4 (vs. T1/2,p < 0.001), N status(p = 0.016), and pathological high grade (vs. low gradeP < 0.001). By analyzing DEGs, there were 1125 genes commonly upregulated, and 209 genes were commonly downregulated. Protein-protein interaction networks further showed the important protein that may be involved in the biological behavior and prognosis of BCa, such as FN1, CXCL12, CD3E, LCK, and ZAP70. A total of 14 DEGs were found to be associated with overall survival of bladder cancer. After validation by a cohort of 165 BCa cases with detailed follow-up information from GSE13507, 10 immune-associated DEGs were demonstrated to be predictive of prognosis in BCa. Among them, 5 genes have not been reported previously associated with the prognosis of BCa, including BTBD16, OLFML2B, PRRX1, SPINK4, and SPON2. Our study elucidated tight associations between stromal score and clinical characteristics as well as prognosis in BCa. Moreover, we obtained a group of genes closely related to the prognosis of BCa in the tumor microenvironment.

Zhao X ，Tang Y ，Ren H ，Lei Y ... -  《-》

Database Mining of Genes of Prognostic Value for the Prostate Adenocarcinoma Microenvironment Using the Cancer Gene Atlas.

Prostate adenocarcinoma (PRAD) is a common malignant tumor in elderly men. Our research uses The Cancer Gene Atlas (TCGA) database to find potential related genes for predicting the prognosis of patients with PRAD. We downloaded gene expression profiles and clinical sample information from TCGA for 490 patients with PRAD (patient age: 41-78 years). We calculated stromal and immune scores using the ESTIMATE algorithm to predict the level of stromal and immune cell infiltration. We categorized patients with PRAD in TCGA into high and low score arrays according to their median immune/stromal scores and identified differentially expressed genes (DEGs) that were significantly correlated with the prognosis of PRAD. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The association between DEGs and overall survival was investigated by weighted Kaplan-Meier survival analysis and multivariate analysis. Furthermore, the protein-protein interaction network (PPI) of DEGs was constructed using the STRING tool. Finally, the hub genes were identified by analyzing the degree of association of PPI networks. We found that 8 individual DEGs, C6, S100A12, MLC1, PAX5, C7, FAM162B, CAMK1G, and TCEAL5, were significantly predictive of favorable overall survival and one DEG, EPYC, was associated with poor overall survival. GO and KEGG pathway analyses revealed that the DEGs were associated with immune responses. Moreover, 30 hub genes were obtained using the PPI network of DEGs: ITGAM, CD4, CD3E, IL-10, LCP2, ITGB2, ZAP-70, C3, CCL19, CXCL13, CXCL9, BTK, CCL21, CD247, CD28, CD3D, FCER1G, PTPRC, TYROBP, CCR5, ITK, CCL13, CCR1, CCR2, CD79B, CYBB, IL2RG, JAK3, PLCG2, and CD19. These prominent nodes had the most associations with other genes, indicating that they might play crucial roles in the prognosis of PRAD. We extracted a list of genes associated with the prostate adenocarcinoma microenvironment, which might contribute to the prediction and interpretation of PRAD prognosis.

Zhao X ，Hu D ，Li J ，Zhao G ，Tang W ，Cheng H ... -  《-》

Identification of Genes with Prognostic Value in the Breast Cancer Microenvironment Using Bioinformatics Analysis.

Ren H ，Hu D ，Mao Y ，Su X ... -  《MEDICAL SCIENCE MONITOR》