Background: Numerous studies have shown that infiltrating eosinophils play a key role in the tumor progression of bladder urothelial carcinoma (BLCA). However, the roles of eosinophils and associated hub genes in clinical outcomes and immunotherapy are not well known.Methods: BLCA patient data were extracted from the TCGA database. The tumor immune microenvironment (TIME) was revealed by the CIBERSORT algorithm. Candidate modules and hub genes associated with eosinophils were identified by weighted gene co-expression network analysis (WGCNA). The external GEO database was applied to validate the above results. TIME-related genes with prognostic significance were screened by univariate Cox regression analysis, lasso regression, and multivariate Cox regression analysis. The patient’s risk score (RS) was calculated and divided subjects into high-risk group (HRG) and low-risk group (LRG). The nomogram was developed based on the risk signature. Models were validated via receiver operating characteristic (ROC) curves and calibration curves. Differences between HRG and LRG in clinical features and tumor mutational burden (TMB) were compared. The Immune Phenomenon Score (IPS) was calculated to estimate the immunotherapeutic significance of RS. Half-maximal inhibitory concentrations (IC50s) of chemotherapeutic drugs were predicted by the pRRophetic algorithm.Results: 313 eosinophil-related genes were identified by WGCNA. Subsequently, a risk signature containing 9 eosinophil-related genes (AGXT, B3GALT2, CCDC62, CLEC1B, CLEC2D, CYP19A1, DNM3, SLC5A9, SLC26A8) was finally developed via multiplex analysis and screening. Age (p < 0.001), grade (p < 0.001), and RS (p < 0.001) were independent predictors of survival in BLCA patients. Based on the calibration curve, our risk signature nomogram was confirmed as a good predictor of BLCA patients’ prognosis at 1, 3, and 5 years. The association analysis of RS and immunotherapy indicated that low-risk patients were more credible for novel immune checkpoint inhibitors (ICI) immunotherapy. The chemotherapeutic drug model suggests that RS has an effect on the drug sensitivity of patients.Conclusions: In conclusion, the eosinophil-based RS can be used as a reliable clinical predictor and provide insights into the precise treatment of BLCA.

BackgroundThe tumor microenvironment (TME) regulates the proliferation and metastasis of solid tumors and the effectiveness of immunotherapy against them. We investigated the prognostic role of TME-related genes based on transcriptomic data of bladder urothelial carcinoma (BLCA) and formulated a prediction model of TME-related signatures.MethodsMolecular subtypes were identified using the non-negative matrix factorization (NMF) algorithm based on TME-related genes from the TCGA database. TME-related genes with prognostic significance were screened with univariate Cox regression analysis and lasso regression. Nomogram was developed based on risk genes. Receiver operating characteristic (ROC) curve and decision curve analysis (DCA) were used for inner and outer validation of the model. Risk scores (RS) of patients were calculated and divided into high-risk group (HRG) and low-risk group (LRG) to compare the differences in clinical characteristics and PD-L1 treatment responsiveness between HRG and LRG.ResultsWe identified two molecular subtypes (C1 and C2) according to the NMF algorithm. There were significant differences in overall survival (OS) (p<0.05), progression-free survival (PFS) (p<0.05), and immune cell infiltration between the two subtypes. A total of eight TME-associated genes (CABP4, ZNF432, BLOC1S3, CXCL11, ANO9, OAS1, FBN2, CEMIP) with independent prognostic significance were screened to build prognostic risk models. Age (p<0.001), grade (p<0.001), and RS (p<0.001) were independent predictors of survival in BLCA patients. The developed RS nomogram was able to predict the prognosis of BLCA patients at 1, 3, and 5 years more potentially than the models of other investigators according to ROC and DCA. RS showed significantly higher values (p = 0.047) in patients with stable disease (SD)/progressive disease (PD) compared to patients with complete response (CR)/partial response (PR).ConclusionsWe successfully clustered and constructed predictive models for TME-associated genes and helped guide immunotherapy strategies.

BackgroundNumerous studies have found that infiltrating M2 macrophages play an important role in the tumor progression of lung adenocarcinoma (LUAD). However, the roles of M2 macrophage infiltration and M2 macrophage-related genes in immunotherapy and clinical outcomes remain obscure.MethodsSample information was extracted from TCGA and GEO databases. The TIME landscape was revealed using the CIBERSORT algorithm. Weighted gene co-expression network analysis (WGCNA) was used to find M2 macrophage-related gene modules. Through univariate Cox regression, lasso regression analysis, and multivariate Cox regression, the genes strongly associated with the prognosis of LUAD were screened out. Risk score (RS) was calculated, and all samples were divided into high-risk group (HRG) and low-risk group (LRG) according to the median RS. External validation of RS was performed using GSE68571 data information. Prognostic nomogram based on risk signatures and other clinical information were constructed and validated with calibration curves. Potential associations of tumor mutational burden (TMB) and risk signatures were analyzed. Finally, the potential association of risk signatures with chemotherapy efficacy was investigated using the pRRophetic algorithm.ResultsBased on 504 samples extracted from TCGA database, 183 core genes were identified using WGCNA. Through a series of screening, two M2 macrophage-related genes (GRIA1 and CLEC3B) strongly correlated with LUAD prognosis were finally selected. RS was calculated, and prognostic risk nomogram including gender, age, T, N, M stage, clinical stage, and RS were constructed. The calibration curve shows that our constructed model has good performance. HRG patients were suitable for new ICI immunotherapy, while LRG was more suitable for CTLA4-immunosuppressive therapy alone. The half-maximal inhibitory concentrations (IC50) of the four chemotherapeutic drugs (metformin, cisplatin, paclitaxel, and gemcitabine) showed significant differences in HRG/LRG.ConclusionsIn conclusion, a comprehensive analysis of the role of M2 macrophages in tumor progression will help predict prognosis and facilitate the advancement of therapeutic techniques.

BackgroundThe tumor microenvironment (TME) is a complex and evolving environment, and the tumor immune microenvironment in kidney renal clear cell carcinoma (KIRC) has a strong suppressive profile. This study investigates the potential prognostic role and value of genes of the tumor microenvironment in KIRC.MethodsThe transcriptome sequencing data of 530 cases and 39 cases of KIRC and the corresponding clinical prognosis information were downloaded from TCGA data and GEO data, respectively, and TME-related gene expression profiles were extracted. A prognostic signature was constructed and evaluated using univariate Cox regression analysis and LASSO regression analysis. Gene set enrichment analysis (GSEA) was used to obtain the biological process of gene enrichment in patients with high and low-risk groups.ResultsA prognostic signature consisting of eight TME-related genes (LRFN1, CSF1, UCN, TUBB2B, SERPINF1, ADAM8, ABCB4, CCL22) was constructed. Kaplan-Meier survival analysis yielded significantly lower survival times for patients in the high-risk group than in the low-risk group, and the AUC values for the ROC curves of this prognostic signature were essentially greater than 0.7, and univariate and multifactorial Cox regression analyses indicated that the risk score was independent risk factors for KIRC prognosis. GSEA analysis showed that immune-related biological processes were enriched in the high-risk group and that risk values were strongly associated with multiple immune cell scores and immune checkpoint-related genes (PDCD1, CTLA4).ConclusionsThe prognostic signature can accurately predict the prognosis of KIRC patients, which may provide new ideas for future precision immunotherapy of KIRC.

Background. Previous studies have shown that RNA N6-methyladenosine (m6A) plays an important role in the construction of the tumor microenvironment (TME). However, how m6A plays a role in the TME of clear cell renal cell carcinoma remains unclear. Methods. Based on 23 m6A modulators, we applied consensus cluster analysis to explore the different m6A modification profiles of ccRCC. The CIBERSORT method was employed to reveal the correlation between TME immune cell infiltration and different m6A modification patterns. A m6A score was constructed using a principal component analysis algorithm to assess and quantify the m6A modification patterns of individual tumors. Results. Three distinct m6A modification patterns of ccRCC were identified. The characteristics of TME cell infiltration in these three patterns were consistent with immune rejection phenotype, immune inflammation phenotype, and immune desert phenotype. In particular, when m6A scores were high, TME was characterized by immune cell infiltration and patient survival was higher ( p < 0.05 ). When m6A scores were low, TME was characterized by immunosuppression and patient survival was lower ( p < 0.05 ). The immunotherapy cohort confirmed that patients with higher m6A scores had significant therapeutic advantages and clinical benefits. Conclusions. The m6A modification plays an important role in the formation of TME. The m6A scoring system allows the identification of m6A modification patterns in individual tumors, discriminates the immune infiltrative features of TME, and provides more effective prognostic indicators and treatment strategies for immunotherapy.