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.

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.

Background: Growing evidence suggests that infiltrating neutrophils are key players in hepatocellular carcinoma (HCC) tumor progression. However, a comprehensive analysis of the biological roles of neutrophil infiltration and related genes in clinical outcomes and immunotherapy is lacking.Methods: HCC samples were obtained from the TCGA and GEO databases. The CIBERSORT algorithm was used to reveal the TIME landscape. Gene modules significantly associated with neutrophils were found using weighted gene co-expression network analysis (WGCNA), a “dynamic tree-cut” algorithm, and Pearson correlation analysis. Genes were screened using Cox regression analysis and LASSO and prognostic value validation was performed using Kaplan-Meier curves and receiver operating characteristic (ROC) curves. Risk scores (RS) were calculated and nomograms were constructed incorporating clinical variables. Gene set variation analysis (GSVA) was used to calculate signaling pathway activity. Immunophenoscore (IPS) was used to analyze differences in immunotherapy among samples with different risk scores. Finally, the relationship between RS and drug sensitivity was explored using the pRRophetic algorithm.Results: 10530 genes in 424 samples (50 normal samples, 374 tumor samples) were obtained from the TCGA database. Using WGCNA, the “MEbrown” gene module was most associated with neutrophils. Nine genes with prognostic value in HCC (PDLIM3, KLF2, ROR2, PGF, EFNB1, PDZD4, PLN, PCDH17, DOK5) were finally screened. Prognostic nomograms based on RS, gender, tumor grade, clinical stage, T, N, and M stages were constructed. The nomogram performed well after calibration curve validation. There is an intrinsic link between risk score and TMB and TIME. Samples with different risk scores differed in different signaling pathway activity, immunopharmaceutical treatment and chemotherapy sensitivity.Conclusion: In conclusion, a comprehensive analysis of neutrophil-related prognostic features will help in prognostic prediction and advance individualized treatment.

BackgroundWith advances in early diagnosis and treatment, the number of cancer survivors continues to grow, and more and more cancer survivors face the threat of second primary cancer (SPM). Second primary pancreatic ductal adenocarcinoma (spPDAC) is an important subclass of SPM, but its prognostic characteristics are poorly understood.MethodsA total of 5,439 spPDAC samples and 67,262 primary pancreatic ductal adenocarcinoma (pPDAC) samples were extracted from the SEER database for this study. Survival differences between spPDAC and pPDAC samples were compared using Kaplan–Meier curves and log-rank tests. The Fine and Gray proportional subdistributed hazard method was used to analyze potential associations between clinical variables and pancreatic ductal adenocarcinoma-specific death (PDACSD) and death from other causes. After that, the clinical variables significantly related to PDACSD were screened out to construct a competing risk nomogram, which was used to evaluate the probability of the occurrence of PDACSD. The C-index was used to evaluate the discriminative ability of the model. The area under the curve (AUC) was used to verify the discrimination of the model. The calibration curve was used to verify the calibration of the model. Decision curve analysis (DCA) was used to validate the clinical utility of the model.ResultsCompared with patients with spPDAC, the pPDAC sample had a better prognosis (p = 0.0017). Across all spPDAC samples, the three most common sites of first-present cancer were the prostate, breast, and digestive system. Age (p < 0.001), race (p = 0.006), interval (p = 0.016), location (p < 0.001), T stage (p = 0.003), M stage (p < 0.001), chemotherapy (p < 0.001), and radiotherapy (p = 0.006) were the clinical variables associated with PDACSD screened by multivariate competing risks analysis. The concordance index values for the training and validation sets were 0.665 (95% CI, 0.655, 0.675) and 0.666 (95% CI, 0.650, 0.682), respectively. AUC, calibration curve, and DCA indicated that the model we constructed had good discrimination, calibration, and clinical utility.ConclusionsIn conclusion, we first analyzed the impact of previous cancer history on prognosis. We then constructed a competing risk model that can predict the probability of developing PDACSD in spPDAC. This model has good discriminative ability, calibration, and clinical practicability and has certain guiding value for clinical decision-making.