PurposeTherapeutic targets of tumor-associated macrophages have been discovered and used clinically as immunotherapy. M2 macrophages are tumor-associated macrophages that promote cancer progression. This article explores the related factors and the effects of type M2 macrophages.MethodWe obtained bladder cancer (BC) sequencing data from TCGA and GSE31189. We used the CIBERSORT algorithm calculate M2 macrophage proportions among 22 type immune cells. The Estimate package was used to measure BC purity. M2 macrophage-related genes were selected using WGCNA. Receiver operating characteristic curves and Kaplan–Meier analyses were performed to determine the risk score, conducted for M2 macrophage-related factors. The Pearson test was used to determine the correlation among M2 macrophage-related genes, clinical phenotype, immune phenotype and tumor mutation burden (TMB). The TIMER database was used to calculate correlations among M2 macrophages and other cancers.ResultsExpression of four M2 macrophages co-expressed genes (CD163, CD209, CSF1, MMD) positively correlated with infiltration of M2 macrophages, which were enriched in the negative regulation of immune system process and the positive regulation of tumor necrosis factor production. M2 macrophage-related factors are robust biomarkers for predicting the BC and immune phenotypes. The Cox regression model built on these four co-expression factors showed a close correlation with outcome (AUC = 0.64). The four co-expression factors negatively correlated outcome and TMB.ConclusionFour co-expressed genes promote high levels of infiltration of type M2 macrophages in the negative regulation of immune system processes and the positive regulation of tumor necrosis factor production processes. These co-expressed genes and the biological process they involve might suggest new strategies for regulation of chemotaxis in M2 macrophages.

PurposeTo identify immune-related co-expressed genes that promote CD8+ T cell infiltration in bladder cancer, and to explore the interactions among relevant genes in the tumor microenvironment.MethodWe obtained bladder cancer gene matrix and clinical information data from TCGA, GSE32894 and GSE48075. The “estimate” package was used to calculate tumor purity and immune score. The CIBERSORT algorithm was used to assess CD8+ T cell proportions. Weighted gene co-expression network analysis was used to identify the co-expression modules with CD8+ T cell proportions and bladder tumor purity. Subsequently, we performed correlation analysis among angiogenesis factors, angiogenesis inhibitors, immune inflammatory responses, and CD8+ T cell related genes in tumor microenvironment.ResultsA CD8+ T cell related co-expression network was identified. Eight co-expressed genes (PSMB8, PSMB9, PSMB10, PSME2, TAP1, IRF1, FBOX6, ETV7) were identified as CD8+ T cell-related genes that promoted infiltration of CD8+ T cells, and were enriched in the MHC class I tumor antigen presentation process. The proteins level encoded by these genes (PSMB10, PSMB9, PSMB8, TAP1, IRF1, and FBXO6) were lower in the high clinical grade patients, which suggested the clinical phenotype correlation both in mRNA and protein levels. These factors negatively correlated with angiogenesis factors and positively correlated with angiogenesis inhibitors. PD-1 and PD-L1 positively correlated with these genes which suggested PD-1 expression level positively correlated with the biological process composed by these co-expression genes. In the high expression group of these genes, inflammation and immune response were more intense, and the tumor purity was lower, suggesting that these genes were immune protective factors that improved the prognosis in patients with bladder cancer.ConclusionThese co-expressed genes promote high levels of infiltration of CD8+ T cells in an immunoproteasome process involved in MHC class I molecules. The mechanism might provide new pathways for treatment of patients who are insensitive to PD-1 immunotherapy due to low degrees of CD8+ T cell infiltration.

Purpose: In the tumor microenvironment, the functional differences among various tumor-associated macrophages (TAM) are not completely clear. Tumor-associated macrophages are thought to promote the progression of cancer. This article focuses on exploring M2 macrophage-related factors and behaviors of renal clear cell carcinoma.Method: We obtained renal clear cell carcinoma data from TCGA-KIRC-FPKM, GSE8050, GSE12606, GSE14762, and GSE3689. We used the “Cibersort” algorithm to calculate type M2 macrophage proportions among 22 types of immune cells. M2 macrophage-related co-expression module genes were selected using weighted gene co-expression network analysis (WGCNA). A renal clear cell carcinoma prognosis risk score was built based on M2 macrophage-related factors. The ROC curve and Kaplan–Meier analysis were performed to evacuate the risk score in various subgroups. The Pearson test was used to calculate correlations among M2 macrophage-related genes, clinical phenotype, immune phenotype, and tumor mutation burden (TMB). We measured differences in co-expression of genes at the protein level in clear renal cell carcinoma tissues.Results: There were six M2 macrophage co-expressed genes (F13A1, FUCA1, SDCBP, VSIG4, HLA-E, TAP2) related to infiltration of M2 macrophages; these were enriched in neutrophil activation and involved in immune responses, antigen processing, and presentation of exogenous peptide antigen via MHC class I. M2-related factor frequencies were robust biomarkers for predicting the renal clear cell carcinoma patient clinical phenotype and immune microenvironment. The Cox regression model, built based on M2 macrophage-related factors, showed a close prognostic correlation (AUC = 0.78). The M2 macrophage-related prognosis model also performed well in various subgroups. Using western blotting, we found that VSIG4 protein expression levels were higher in clear renal cell carcinoma tissues than in normal tissues.Conclusion: These co-expressed genes were most related to the M2 macrophage phenotype. They correlated with the immune microenvironment and predicted outcomes of renal clear cell carcinoma. These co-expressed genes and the biological processes associated with them might provide the basis for new strategies to intervene via chemotaxis of M2 macrophages.

The morbidity and mortality of prostate carcinoma has increased in recent years and has become the second most common ale malignant carcinoma worldwide. The interaction mechanisms between different genes and signaling pathways, however, are still unclear. Methods Variation analysis of GSE38241, GSE69223, GSE46602 and GSE104749 were realized by GEO2R in Gene Expression Omnibus database. Function enrichment was analyzed by DAVID.6.8. Furthermore, the PPI network and the significant module were analyzed by Cytoscape, STRING and MCODE.GO. Pathway analysis showed that the 20 candidate genes were closely related to mitosis, cell division, cell cycle phases and the p53 signaling pathway. A total of six independent prognostic factors were identified in GSE21032 and TCGA PRAD. Oncomine database and The Human Protein Atlas were applied to explicit that six core genes were over expression in prostate cancer compared to normal prostate tissue in the process of transcriptional and translational. Finally, gene set enrichment were performed to identified the related pathway of core genes involved in prostate cancer. Result Hierarchical clustering analysis revealed that these 20 core genes were mostly related to carcinogenesis and development. CKS2, TK1, MKI67, TOP2A, CCNB1 and RRM2 directly related to the recurrence and prognosis of prostate cancer. This result was verified by TCGA database and GSE21032. Conclusion These core genes play a crucial role in tumor carcinogenesis, development, recurrence, metastasis and progression. Identifying these genes could help us to understand the molecular mechanisms and provide potential biomarkers for the diagnosis and treatment of prostate cancer.

No abstract

PurposeTo identify CD8+ T cell-related factors and the co-expression network in melanoma and illustrate the interactions among CD8+ T cell-related genes in the melanoma tumor microenvironment.MethodWe obtained melanoma and paracancerous tissue mRNA matrices from TCGA-SKCM and GSE65904. The CIBERSORT algorithm was used to assess CD8+ T cell proportions, and the “estimate” package was used to assess melanoma tumor microenvironment purity. Weighted gene co-expression network analysis was used to identify the most related co-expression modules in TCGA-SKCM and GSE65904. Subsequently, a co-expression network was built based on the joint results in the two cohorts. Subsequently, we identified the core genes of the two most relevant modules of CD8+T lymphocytes according to the module correlation, and constructed the signature using ssGSEA. Later, we compared the signature with the existing classical pathways and gene sets, and confirmed the important prognostic significance of the signature in this paper.ResultsNine co-expressed genes were identified as CD8+ T cell-related genes enriched in the cellular response to interferon−gamma process and antigen processing and presentation of peptide antigen. In the low expression level group, inflammation and immune responses were weaker. Single-cell sequencing and immunohistochemistry indicated that these nine genes were highly expressed in CD8+ T cells group.ConclusionWe identified nine-gene signature, and the signature is considered as the biomarker for T lymphocyte response and clinical response to immune checkpoint inhibitors for melanoma

Background There is evidence that long non-coding RNA (lncRNA) is related to genetic stability. However, the complex biological functions of these lncRNAs are unclear. Method TCGA - KIRC lncRNAs expression matrix and somatic mutation information data were obtained from TCGA database. “GSVA” package was applied to evaluate the genomic related pathway in each samples. GO and KEGG analysis were performed to show the biological function of lncRNAs-mRNAs. “Survival” package was applied to determine the prognostic significance of lncRNAs. Multivariate Cox proportional hazard regression analysis was applied to conduct lncRNA prognosis model. Results In the present study, we applied computational biology to identify genome-related long noncoding RNA and identified 26 novel genomic instability-associated lncRNAs in clear cell renal cell carcinoma. We identified a genome instability-derived six lncRNA-based gene signature that significantly divided clear renal cell samples into high- and low-risk groups. We validated it in test cohorts. To further elucidate the role of the six lncRNAs in the model’s genome stability, we performed a gene set variation analysis (GSVA) on the matrix. We performed Pearson correlation analysis between the GSVA scores of genomic stability-related pathways and lncRNA. It was determined that LINC00460 and LINC01234 could be used as critical factors in this study. They may influence the genome stability of clear cell carcinoma by participating in mediating critical targets in the base excision repair pathway, the DNA replication pathway, homologous recombination, mismatch repair pathway, and the P53 signaling pathway. Conclusion subsections These data suggest that LINC00460 and LINC01234 are crucial for the stability of the clear cell renal cell carcinoma genome.

Aim. In this paper, we aimed to develop and validate a risk prediction method using independent prognosis genes selected robustly in prostate cancer. Method. We considered 723 samples obtained from TCGA (the Cancer Genome Atlas), GSE46602, and GSE21032. Prostate cancer prognosis-related genes with P<0.05 were selected using Univariable Cox regression analysis. We then built the lowest AIC (Akaike information criterion score) optimal gene model using the “Rbsurv” package in TCGA train set. The coefficients were obtained by Multivariable Cox regression analysis. We named the new prognosis method CMU5. The CMU5 risk score was verified in TCGA test set, GSE46602, and GSE21032. Results. FAM72D, ARHGAP33, TACR2, PLEK2, and FA2H were identified as independent prognosis factors in prostate cancer patients. We built the computing model as follows: CMU5 risk score = 1.158∗FAM72D + 1.737∗ARHGAP33 − 0.737∗TACR2 − 0.651∗PLEK2 − 0.793∗FA2H. The AUC of DFS was 0.809 in the train set (274 samples), 0.710 in the test set (273 samples), and 0.768 in the complete set (547 samples). The benign prediction capacity of CMU5 was verified by GSE46602 (36 samples; AUC=0.6039) and GSE21032 GPL5188 (140 samples; AUC=0.7083). Using the cut-off point of 2.056, a significant difference was shown between high- and low-risk groups. Conclusion. A prognosis-related risk score formula named CMU5 was built and verified, providing reliable prediction of prostate cancer outcome. This signature might provide a basis for individualized treatment of prostate cancer.