- 1The Fifth Affiliated Hospital of Zhengzhou University, Zhengzhou University, Zhengzhou, China
- 2Department of Orthopaedics, Affiliated National Hospital of Guangxi Medical University, Nanning, China
- 3College of Pharmacy, Shantou University School of Medicine, Shantou, China
- 4Department of Biology, Chemistry, Pharmacy, Free University of Berlin, Berlin, Germany
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.
Introduction
Bladder urothelial carcinoma (BLCA) is a malignancy of the bladder uroepithelium with approximately 500,000 newly diagnosed cases and 200,000 deaths worldwide each year (Siegel et al., 2022). A quarter of BLCA cases are muscle-invasive bladder cancer (MIBC) by statistics (Robertson et al., 2017), while non-muscle invasive bladder cancer (NMIBC) has a fast progression and a high recurrence rate (Sylvester et al., 2016). The recurrence probability of NMIBC patients undergoing transurethral resection and combination therapy is as high as 50–70% (Robertson et al., 2018). However, MIBC often has a poor prognosis due to early occult metastatic spread, and its 5-years survival rate does not exceed 50% (Hurst and Knowles, 2022). Despite multiple therapies, overall survival and prognosis for BLCA remain poor severely. Therefore, emerging biomarkers are urgently identified for the clinic.
Increased eosinophil counts are commonly seen in various malignancies (Samoszuk, 1997), it was reported that eosinophils in the peripheral blood of cancer patients increased compared with normal people (Olsen and Spry, 1985). Eosinophils are generally considered granulocytes involved in allergic disease responses and parasitic infection responses. Eosinophilic tumor infiltration occurs during the proliferation of certain types of tumors, such as breast cancer and gastric adenocarcinoma (Costello et al., 2005). Eosinophils readily respond to a variety of stimuli by synthesizing and secreting a large number of molecules, including unique granulin proteins that may kill tumor cells (Grisaru-Tal et al., 2020). Studies have shown that some cytokines (IL-4, IL-5) can increase the anti-tumor activity of eosinophils (Furbert-Harris et al., 2003).
Immunotherapy made tremendous progress and has been an indispensable component in the treatment of advanced solid tumors in the past decade. Cancer immunotherapy is not limited to understanding the relationship between cancer and the immune system, and immunocytes have also been used to predict the prognosis of cancer patients (Havel et al., 2019; Kruger et al., 2019). Immune checkpoint inhibitors (ICIs) have been formally used in BLCA in 2016 (Özdemir et al., 2018), moreover, immunotherapy targeting immune checkpoints is still based on PD-L1 inhibitors (Rhea et al., 2021). The safety and efficacy of immunotherapy for BLCA have long been proven relative to other second-line treatments in clinical trials (Bellmunt et al., 2017). T cells, macrophages, and tumor-associated fibroblasts, myeloid-derived suppressor cells in the microenvironment, almost significantly inhibit effective immunotherapy (Hessmann et al., 2020). Among them, eosinophil and eosinophil-associated signaling pathways play a key role in antitumor immunity and induction of epithelial mesenchymal transition (De Palma et al., 2017; Korbecki et al., 2020). Therefore, the strategy of autoimmune profiling may be the most reliable and promising strategy to comprehensively assess tumor sensitivity to clinical therapy. However, the survival value of eosinophils in BLCA patients and the biological role of the tumor immune microenvironment (TIME) are still not evident.
We investigated the potential role of eosinophils in BLCA TIME with the TCGA-BLCA and GSE31684 cohorts. Eosinophil analysis was obtained using the CIBERSORT algorithm, and significant modules related to eosinophils were identified by weighted gene co-expression network analysis (WGCNA). We established a nomogram combining risk signature and clinical features according to eosinophil-associated genes. In addition, we also verified the synergistic effect of risk score (RS) and tumor mutational burden (TMB), the relationship of immunotherapy responsiveness and drug sensitivity.
Methods
Data Collection
Transcriptome data and clinical annotations were downloaded from the TCGA-BLCA (The Cancer Genome Atlas-Bladder urothelial carcinoma) cohort consisting of 433 RNA-sequencing samples, including 19 normal samples and 414 tumor samples. After removing the sample without clinical follow-up information from the TCGA-BLCA cohort, 404 cases were applied for the following analysis. The four types of somatic mutation data in BLCA patients were also obtained from the TCGA database. Furthermore, we took the GSE31684 cohort data from the GEO database for external validation.
The Landscape of Immune Cell Infiltration
Sequencing data from computational samples were analyzed via the CIBERSORT algorithm to obtain the abundance of 22 tumor-infiltrating immune cell (TIC) subtypes. The cellularity of the tumor immune microenvironment were characterized by these subtypes (Langfelder and Horvath, 2008).
Weighted Gene Co-Expression Network Analysis
The 16,394 gene matrix from the TCGA-BLCA cohort was imported into the CIBERSORT. The soft threshold power (β) from 1 to 20 was selected as a candidate, and the corresponding power values were calculated using the pickSoftThreshold function, the best power value was selected to build the proximity matrix, and our gene distribution was made to conform to the scale-free network according to the connectivity. Using the TOM matrix obtained from gene expression, the genes were again continued to be clustered, the minimum number of module genes was set, and the gene clustering results were cut to obtain different gene modules. The “dynamic tree cutting” algorithm was used to introduce similar genes into the same candidate module. In addition, to realize the correlation analysis between the module eigengenes and the sample traits, we performed the Pearson correlation test was performed. Finally, we focused on the “eosinophils” population and extracted the modules most significantly associated with eosinophils for subsequent analysis.
Construction of Eosinophil-Related Prognostic Features
The 313 genes from the most important module were taken for further screening to discuss the predictive effects of eosinophil-related genes. First, candidate genes significantly associated with overall survival were screened out by univariate Cox regression analysis (p < 0.05). Then the range of candidate genes is further narrowed by LASSO algorithm. Finally, 9 genes associated with overall survival were finally screened out by multivariate Cox regression. A prognostic risk model including these 9 eosinophil-related genes was established.
RS was calculated according to the expression of differential genes and regression analysis coefficient values. The formula is shown below:
Here, β is the regression coefficient in the multivariate Cox regression analysis (Lossos et al., 2004).
Validation of Prognostic Eosinophil-Related Features
According to the above formula, there is one RS for each BLCA sample. All BLCA samples were divided into two subgroups, high-risk group (HRG) and low-risk group (LRG), using the median of RS as the cutoff value. A K-M survival curve was drawn to compare the prognostic differences between the two subgroups, and this process was performed in the “survival” package in R. In addition, transient receiver operating characteristic (ROC) curves were plotted to assess performance. Then, univariate and multivariate Cox regression analysis was performed to confirm the RS as a meaningful independent influencing factor for BLCA patients. The correlations of RS with clinical variables were visualized by the “pheatmap” package.
Gene Set Enrichment Analysis
The c2.cp.kegg.v7.4.symbols collection were used to explore the function annotation in HRG and LRG by GSEA software. Results with p value < 0.05 were being considered statistically significant. The first eight results were selected for visualization.
Establishment and Verification of the Nomogram
ROC analysis of gender, age, RS, tumor grade, and clinicopathological stage was performed to determine the best prognostic indicators for 1, 3, and 5-years overall survival (OS) in patients (Blanche et al., 2013). To provide a quantifiable tool to predict 1-, 3-, and 5-years OS, the nomogram integrating risk signature and other clinicopathological features were constructed. Finally, the calibration curves were drawn to verify the performance of the nomogram.
Correlation Between Tumor Mutational Burden and Risk Score
Somatic mutation data from the TCGA-BLCA cohort were obtained from the TCGA database. TMB refers to the number of somatic non-synonymous mutations in a specific genomic region, which can indirectly reflect the ability and degree of tumor production of neoantigens, and predict the immunotherapy efficacy of various tumors. The “maftools” R package (Yoshihara et al., 2013) was used to detect the number of somatic non-synonymous point mutations in each sample.
Correlation of Risk Score With Tumor Immune Microenvironment Characterization
Seven methods were applied to assess immune cell infiltration in the tumor microenvironment, including XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT, and CIBERSORT-ABS. ESTIMATE algorithm, which can be based on gene expression data, estimates the stromal score and immune score of a tumor sample for representing the presence of stromal and immune cells. The two scores are summed to obtain the ESTIMATE score, which can be used to estimate tumor purity. Correlation between RS and TICs was performed using Spearman correlation analysis.
Gene Set Variation Analysis
We downloaded hallmark gene sets from the MSigDB database (Chan et al., 2019). To assess relative pathway activity in individual samples, the applied gene set variation analysis (GSVA) (Rizvi et al., 2015) in the GSVA package was used to assign pathway activity estimates.
Predict Patient Response to Immunotherapy
Immune checkpoints have been defined as key targets for the inhibition of immune cell function. In this study, we analyzed the expression levels of 47 immune checkpoint blockage-related genes in HRG and LRG. Immunophenoscore (IPS) determines the immunogenicity of a tumor and predicts the response to immune checkpoint inhibitor therapy. IPS calculates scores for each of the four different immunophenotypes (antigen-presenting, effector, suppressor, and checkpoint), and the IPS z-score is an integration of all four, and the higher the IPS z-score, the more immunogenic the sample (Charoentong et al., 2017).
Predicting Chemotherapy Effect
To investigate the correlation of RS with chemotherapeutic drug sensitivity, we used the R package “pRRophetic” to appraise the half-maximal inhibitory concentration (IC50) of BLCA samples between HRG and LRG. The package “pRRophetic” can estimate the IC50 of chemotherapeutic drugs by constructing a ridge regression model based on Genomics of Cancer Drug Sensitivity (GDSC) (www.cancerrxgene.org/) cell line expression profiles and TCGA gene expression profiles (Rosenberg et al., 2016).
Statistical Analysis
The Wilcoxon test was used to compare two groups, and the Kruskal–Wallis test was used to compare more than two groups. The Kaplan-Meier log-rank test was used to analyze the survival curves. Associations between RS subgroups and somatic mutational burden were analyzed with the chi-square test, and Spearman analysis was used to calculate correlation coefficients. Results were further analyzed using the CIBERSORT algorithm (p < 0.05). Two-sided p < 0.05 was considered statistically significant. All statistical analyses were performed using R software (version 4.1.1).
Results
Landscape of Tumor Immune Microenvironment in Bladder Urothelial Carcinoma
We employed the CIBERSORT algorithm (Supplementary Table S1) to clarify the synthesis of TIME. The abundance of the 22 TIC types is shown in Figure 1A. To explore and characterize the relationship between TIME patterns and clinical phenotypes, we created a comprehensive heatmap (Figure 1B). We visualized the integrated landscape of TIME to show the correlation between immune cells (Figure 1C), and found that T cells CD4 memory activated and T cells CD8 have the largest positive correlation (r = 0.6, p < 0.05), while T cells CD8 and T cells CD4 memory resting with the largest negative correlation (r = −0.55, p < 0.05). The largest positive correlation with Eosinophils was T cells CD4 memory resting and Monocytes (r = 0.13, p < 0.05), and the largest negative correlation was T cells CD8 (r = −0.13, p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g001.jpg)
FIGURE 1. (A) Immune cell infiltration in the BLCA tumor immune environment. Subpopulations of 22 immune cell subtypes (B) and scale heatmap of 22 TICs in each BLCA sample (C). Intrinsic correlation of 22 infiltrating immune cells in EBLCA.
Establishment of the Weighted Gene Co-Expression Network Analysis Network
Based on the 16,394 genes as well as immune-infiltrating subsets, we construct the scale-free network by defining the first power value when the scale-free topology index reaches 0.90, which is 9 (Figure 2A) and setting it as the optimal soft threshold power (β). The similarly expressed pattern genes were introduced into the same module through a dynamic tree-cutting algorithm (module size = 60) to form a hierarchical clustering tree with different modules. Hierarchical cluster analysis was performed based on weighted correlations and the clustering results were segmented according to the set criteria to obtain 13 gene modules (Figure 2B). Each column of Figure 2C represents a TIC type, and each row represents candidate modules with trait carrier genes. The green-yellow module is highly correlated with eosinophils (r = 0.2, p < 0.01) in Figure 2C, the genes in this module (Supplementary Table S2) were applied for further studies.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g002.jpg)
FIGURE 2. Choosing an appropriate soft threshold (power) and building a hierarchical clustering tree (A) The choice of the soft threshold enables the scale-free topology to achieve an exponent of 0.90, and the average connectivity for 1–20 soft threshold powers is analyzed (B) M2 macrophage-related genes with similar expression patterns were merged into the same module using a dynamic tree-cutting algorithm, creating a hierarchical clustering tree. Heatmap of correlations between (C) modules and immune-infiltrating cells (traits).
Development of Risk Signatures
Twenty-seven eosinophil-related genes were screened by univariate Cox regression analysis (p < 0.05, Supplementary Table S3), and a forestogram was drawn (Figure 3A). To avoid overfitting, we performed Lasso regression on these 27 genes and identified 18 eosinophil-related genes (Figure 3B), furthermore determining the optimal value of the penalty parameter by cross-validation (Figure 3C). We hence identified 9 eosinophil-related genes (AGXT, B3GALT2, CCDC62, CLEC1B, CLEC2D, CYP19A1, DNM3, SLC5A9, SLC26A8) as hub genes by multivariate Cox regression. All nine genes were considered beneficial prognostic indicators (Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g003.jpg)
FIGURE 3. (A) Nomogram between 27 eosinophil-related genes and overall survival screened by univariate Cox regression analysis (B) LASSO coefficient profiles of 22 candidate genes. A vertical line is drawn at the value chosen by 10-fold cross-validation (C) Ten-fold cross-validation for tuning parameter selection in lasso regression. Vertical lines are drawn from the best data according to the minimum criterion and 1 standard error criterion. The vertical lines on the left represent the final 9 genes identified (D,G) Kaplan-Meier curve analysis showing the difference in overall survival between high-risk and low-risk groups (E) Univariate Cox regression results for overall survival (F) Multivariate Cox regression results for overall survival.
More importantly, nine hub genes were incorporated into the risk profile of BLCA patients. The RS was computed:
Finally, BLCA samples were divided into HRG and LRG according to the Bit RS.
Validation of Risk Prognostic Features
We found the OS time of HRG patients in the TCGA cohort is significantly shorter than that of LRG patients through the K-M survival curve (p < 0.001, Figure 3D), which proved that our RS has a good performance. There was also a difference in OS between the two groups of patients with high and low hub gene expression (Figure 3G). In addition, we showed RS in the TCGA cohort (Figure 4A), polygenic model RS distribution (Figure 4B), survival status, and longevity of BLCA patients (Figure 4C). By univariate Cox regression analysis, we obtained a hazard ratio (HR) of 1.096 for RS (95% CI 1.066–1.127, Figure 3E). The results of multivariate Cox regression analysis (HR = 1.087, 95% CI 1.057–1.117, Figure 3F) supported RS as an independent prognostic indicator in BLCA. The above results clarified these 9 hub genes have a good prediction of clinical outcomes. More importantly, this feature was validated by an external database GSE31684. The outcomes showed the distribution of nine gene expression patterns, sample prognostic information, and RS in the external validation cohort (Figures 4D–F). These results suggested that this feature has a stable and robust prognostic value.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g004.jpg)
FIGURE 4. (A) Confirmation of prognostic risk scores in the TCGA cohort (B) Polygenic model risk score distribution in the TCGA cohort (C) Survival status and duration of BLCA patients in the TCGA cohort (D) Confirmation of prognostic risk scores in the GEO cohort (E) Polygenic model risk score distribution in the GEO cohort (F) Survival status and duration of BLCA patients in the GEO cohort.
Functional Analysis of Eosinophil-Related Genes
All samples were divided into low and high expression groups by the median expression of the hub genes that served as the cut-off point. Functional enrichment of high and low hub gene expression was identified by GSEA subsequently. The KEGG enrichment analysis showed that the high expression of DNM3 was related to signaling pathways such as arrhythmogenic right ventricular cardiomyopathy and neuroactive ligand-receptor interaction; the high expression of SLC5A9 was related to signaling pathways such as drug metabolism—cytochrome p450 and neuroactive ligand-receptor interaction; the high expression of SLC26A8 was related to olfactory transduction and drug metabolism—cytochrome p450 and other signaling pathways (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g005.jpg)
FIGURE 5. GSEA for samples with high and low expression of 9 central genes (A) Enriched gene set collected in KEGG for samples with high B3GALT2 expression (B) Enriched gene set collected in KEGG for samples with high CCDC62 expression (C) Enriched gene set collected in KEGG for samples with high CLEC1B expression (D) Enriched gene set collected in KEGG for samples with high CLEC2D expression (E) Enriched gene set collected in KEGG for samples with high CYP19A1 expression (F) High DNM3 expression sample enriched gene set collected in KEGG (G) Enriched gene set collected in KEGG for samples with high SLC5A9 expression (H) Enriched gene set in the KEGG collection for samples with high SLC26A8 expression.
Association of Risk Characteristics With Clinicopathological Variables
To visualize the distribution of clinical variables in the LRG/HRG subgroup, we plotted Figure 6A. Clinical subtype scores between HRG and LRG based on sex, tumor grade, clinical stage, T, N, and M stage are shown in Figures 6B–G.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g006.jpg)
FIGURE 6. Clinical significance of prognostic risk characteristics (A) heatmap showing the distribution of clinical characteristics and corresponding risk scores in each sample. Incidence of clinical variable subtypes of LRG/HRG (B) Gender (C) Grade (D) Stage (E) Stage T (F) Stage N, and (G) Stage M.
Construction of Prognostic Nomogram
The ROC curves indicated a high prognostic validity, with AUC values of 0.692, 0.755, and 0.772 for 1-, 3-, and 5-years OS, respectively (Figure 7C). Moreover, we combined gender, tumor grade, T, N, and M stages for AUC analysis of 1-year (Figure 7D), 3-years (Figure 7E), and 5-years (Figure 7F) OS to further validate that RS had the best prognostic value among multiple clinicopathological variables. The results elucidated that RS obtained the highest AUC value. Then, we developed a prognostic nomogram composed of RS, gender, tumor grade, clinical stage, T, N, and M stage for quantitative prognostic prediction (Figure 7A). The survival rate of BLCA patients can be accurately predicted based on nomogram. For example, an 80-year-old male was diagnosed with bladder cancer. The patient’s clinical information was T3N1M0 with stage III and had a low RS. We could then calculate the patient’s total score as 495, corresponding to survival rates of 0.911, 0.713, and 0.609 for 1, 3, and 5 years, respectively. Finally, the calibration curve was established to evaluate the predictive performance of the nomogram, which showed good prognostic prediction performance (Figure 7B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g007.jpg)
FIGURE 7. Validation of prognostic efficiency of risk signatures (A) Nomogram was used to predict survival in BLCA patients (B) 1-, 3-, and 5-years nomogram calibration curves (C) ROC analysis was used to estimate the predictive value of prognostic features. The area under the (D–F) curve (AUC) of the risk score for predicting overall survival at 1, 3, and 5 years and other clinical characteristics.
Association of Risk Signature With Tumor Mutational Burden
To reveal the relationship between RS and TMB and to reveal the genetic variation of RS subtypes, we plotted K-M survival curves based on TBM levels, indicating that high TMB values led to longer OS (p < 0.001, Figure 8A). We divided the patients into four subgroups according to the median of RS and TMB. As shown in Figure 8B, the RS subgroup exhibited significant prognostic differences with the low and high TMB subgroup, and TMB levels did not interfere with the RS prognostic prediction performance (p < 0.001). These results further confirm that RS may serve as an independent prognostic predictor and has the potential to evaluate the clinical outcomes of antitumor immunotherapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g008.jpg)
FIGURE 8. Correlation between risk score and TMB (A) Kaplan-Meier curves of high TMB and low TMB groups (B) Kaplan-Meier curve stratification of patients according to TMB and risk signature. The oncoPrint was constructed using high risk score (C) and low risk score (D).
To explore the relationship between RS and gene mutations, we analyzed the distribution of gene mutations between HRG and LRG and mapped the top 20 gene waterfalls with the most frequent somatic mutations between HRG and LRG (Figures 8C,D). Significantly, mutated gene (SMG) mutation profiles indicated that TP53 (49 vs. 45%) experienced a higher rate of somatic mutation in HRG core subtypes, while FGFR3 (21 vs. 7%) had a higher rate of LRG somatic mutation rate. These findings may contribute to new insights into the intrinsic link between eosinophil infiltration and somatic variation in BLCA immunotherapy.
Risk Signature in Tumor Immune Microenvironment Context
To explore the role of RS and TIME, we analyzed the association of RS with immune infiltration by Spearman correlation (Figure 9A, Supplementary Table S5). The ESTIMATE analysis showed a significant upward trend in HRG matrix scores. Furthermore, ESTIMATE scores were significantly upregulated in higher-risk samples (Figure 9B). Validation of the correlation predicted by the two methods CIBERSORT (Figure 9C) and XCELL (Figure 9D) showed that our analysis is dependable.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g009.jpg)
FIGURE 9. Estimated abundance of tumor-infiltrating cells. Patients in the (A) high-risk group had a stronger correlation with tumor-infiltrating immune cells, as shown by the Spearman correlation analysis (B) Association between prognostic risk signatures and central immune checkpoint genes. The correlations predicted by the two methods CIBERSORT (C) and XCELL (D) were validated.
Enrichment of Signaling Pathways in High-Risk Group/Low-Risk Group
As shown in Figure 10A, we performed gene set variation analysis (GSVA) to further reveal the biological roles of different risk groups in tumorigenicity and progression. Subjects in the LRG showed enhanced KEGG/PPAR pathway activity. Most genes with high expression levels in the HRG were enriched in KEGG/WNT, KEGG/TGF/BETA, KEGG/MAPK, and KEGG/CHEMOKINE signaling pathway.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963455/fmolb-09-963455-HTML/image_m/fmolb-09-963455-g010.jpg)
FIGURE 10. Enrichment pathways for GSVA (A) heatmap showing the correlation of representative pathway items on KEGG with risk scores. Predicting immunotherapy response (B) Association of immune checkpoint blockade gene expression levels with risk scores (C–F) IPS score distribution map. Estimates of chemotherapy effect risk scores. Sensitivity analysis of (G) Veliparib in patients with high and low risk scores. Sensitivity analysis of (H) cisplatin in patients with high and low risk scores (I) Sensitivity analysis of paclitaxel in patients with high and low risk scores (J) Sensitivity analysis of methotrexate in patients with high and low risk scores.
Immunotherapy Prediction
Since there is no data information for immunotherapy in the TCGA-BLCA dataset, we indirectly analyzed the response to immunotherapy. As shown in Figure 10B, most immune checkpoint blockade-related genes (TNFSF9, TNFSF4, TNFRSF8) were significantly positively associated with RS.
In this risk scoring system, patients in HRG have a higher IPS score at the time of immunotherapy (Figure 10C: PD1-positive and CTLA4-negative; Figure 10D: PD1-negative and CTLA4-positive). It is easy to draw from Figure 10F that LRG is more suitable for novel immunotherapy (Figure 10F: PD1-negative and CTLA4-negative). These results confirmed that RS is associated with immunotherapy.
Prediction of Response to Chemotherapy
We estimated the IC50s of 4 chemotherapeutic agents (Veliparib, paclitaxel, cisplatin, and methotrexate) in BLCA patients. Veliparib and methotrexate exhibited higher IC50s in HRG signature patients (all p < 0.001; Figures 10G,J). In contrast, paclitaxel and cisplatin exhibited higher IC50s in LRG patients (all p < 0.001; Figures 10H,I). These results suggested that RS effects the drug sensitivity of patients.
Discussion
As the ninth most common cancer in the world, BLCA with high recurrence and high invasiveness poses a great threat to the health of all human beings (Cumberbatch et al., 2018). Cisplatin-based chemotherapy is the universal treatment for locally advanced or metastatic BLCA in the clinic currently (von der Maase et al., 2005), however, the efficacy is usually unsatisfactory due to the physical condition of the patients and various adverse reactions (Chen et al., 2021). In recent years, immunotherapy has developed rapidly, with a variety of ICIs appearing and being used to treat advanced BLCA (Chism, 2017). Although the safety and efficacy of these treatments have been confirmed by clinical studies, little is known about the effectiveness of ICI treatment. Although investigators have identified various markers associated with the immunotherapy response, such as TMB, TIC, and immune gene signatures, the application of immunotherapy has been limited due to limitations and differences between studies (Gibney et al., 2016; Chan et al., 2019; Samstein et al., 2019).
In studies of BLCA, increasing attention has been paid to infiltrating immune cells (Jing et al., 2022), especially eosinophils. It is well known that eosinophils have been shown to infiltrate multiple tumors and are able to regulate tumor progression either directly by interacting with tumor cells or indirectly by shaping TIME. Previous reports elucidate the complex interactions between different immune cell types in TIME and the effects on the prognosis of cancer patients (Varn et al., 2017). The pharmacological blocker of the eosinophilic immune checkpoint is an approach for cancer immunotherapy. These molecules primarily control T-cell activity, including cytotoxic T-lymphocyte-associated protein4, PD1, and PDL1 (Pardoll, 2012; Webster, 2014).
In the current work, we collected two separated cohorts from TCGA-BLCA and GSE31684 to explore the potential roles of eosinophil-related genes in BLCA patients. 414 tumor samples and 16,394 associated genes were used for further studies. Twenty-two tumor-infiltrating immune cell subsets were obtained using the CIBERSORT algorithm. Next, we identified a module (greenyellow) that was positively associated with eosinophil-related genes using the WGCNA approach, which contained 313 candidate genes. After a series of screening, nine eosinophil-related genes (AGXT, B3GALT2, CCDC62, CLEC1B, CLEC2D, CYP19A1, DNM3, SLC5A9, SLC26A8) were finally identified. RS are then calculated and prognostic features are constructed. Risk characteristics were considered independent prognostic predictors of BLCA and performed well in univariate and multivariate regression analyses.
Using known samples, we constructed a nomogram based on RS and other clinical data, including gender, tumor grade, clinical stage, T, N, and M stage. It reduces statistical prediction models to a single numerical estimate of 1-, 3-, and 5-years OS tailored to individual patient conditions. These predictions of individual clinical outcomes can provide valuable information for clinical decision-making. The nomogram transforms the complex regression formula into a visual graph, making the results of the prediction model more readable and more valuable (Iasonos et al., 2008; Wang et al., 2022). Therefore, it is suggested that nomogram can be more widely concerned and applied in future medical research and clinical practice.
Studies have found a correlation between genetic alternation and immunotherapy response (Burr et al., 2017; George et al., 2017). Our results showed that TMB, a predictor of immunotherapy sensitivity, increased significantly with RS. By plotting the stratified survival curves, we found that RS had prognostic predictive power independent of TMB, which also confirmed that TMB and RS represent different aspects of immunobiology. In addition, our study also found significant differences between HRG and LRG gene variant frequencies and transcriptome levels.
The PPAR signaling pathway of LRG was activated in GSVA, while HRG was associated with WNT, TGF-BETA, MAPK, and CHEMOKINE signaling pathways. These results suggest that there exist significant differences in the molecular mechanisms between the high and low risk groups. In addition, risk scoring schemes showed that chemotherapeutic drug sensitivity was associated with RS. Based on the above results, we speculated BLCA patients may be better suited for different combinations of molecularly targeted and chemotherapeutic agents based on risk stratification.
In addition, RS was significantly and positively correlated with ICI-related gene CTLA4, indicating that high-risk samples might be more affected by immune checkpoint blockade. In the IPS risk scoring system, there were high IPS scores in the HRG patients (except PD-1/PD-L1/PD-L2 positivity and CTLA-4 positivity), suggesting low-risk patients are more suitable for novel ICB-targeted therapy, whereas not PD1/CTLA4 combination immunotherapy.
Through the previous RS formula, it is not difficult to find that the expression of SLC5A9 and CLEC2D genes is negatively correlated with the RS. Interestingly, however, we found that SLC5A9 and CLEC2D genes had opposite associations with immune checkpoint-related genes (Figure 10B). CLEC2D gene was positively correlated with most immune checkpoint-related genes, while SLC5A9 gene was negatively correlated with most immune checkpoint-related genes. This may seem contradictory, but it is not. The impact of eosinophil-associated genes on the risk score is comprehensive and comprehensive, and its impact on the prognosis of BLAC includes tumor occurrence and progression, as well as the therapeutic effects of chemotherapy, surgery, and immunotherapy. Therefore, it is understandable that there are differences in the effects of SLC5A9 and CLEC2D genes on the effect of immunotherapy. Nowadays, human research on tumor progression and therapeutic effects has only taken a small step, and there are too many unknowns waiting to be explored and discovered by human beings.
Among the nine eosinophil-related genes we screened, DNM3 gene plays an important role in tumor biological progression and clinical prognosis. Ma et al. found that the DNM3 gene can inhibit the proliferation and metastasis of colon cancer cells through the AKT signaling pathway (Ma et al., 2019). Cheng and his colleagues found that the expression level of the DNM3 gene in colorectal cancer was significantly correlated with tumor size and clinical case type (Cheng et al., 2022). Zhang et al. found that DNM3 can attenuate the growth of hepatocellular carcinoma by activating p53 (Zhang et al., 2016). Jing Fa found that DNM3 overexpression could inhibit the proliferation, migration, and invasion of cervical cancer cells (Fa, 2021). The specific mechanism of DNM3 gene in BLCA in tumor biological progression and clinical prognosis is still unclear. This study is expected to have a certain reference value for DNM3 gene tumor research.
In addition to DNM3 gene, the research of other 8 eosinophil-related genes in the tumor also made gratifying progress in recent years. The CCDC62 gene plays a role in tumor progression of gastric (Domae et al., 2009) and prostate cancers (Chen et al., 2009). CLEC1B gene is associated with pancreatic cancer cachexia (Narasimhan et al., 2020). Overexpression of the CLEC2D gene in prostate cancer affects the biological process of the tumor (Mathew et al., 2016; Buller et al., 2020). The CYP19A1 gene is associated with breast (Wu et al., 2017; Lv et al., 2021) and lung cancer (Shi et al., 2008). This paper corroborates the biological roles of these genes in tumors in the above studies, and deepens the understanding of the biological roles of eosinophil-related genes in BCLA.
In short, this study identifies potential directions for research related to BLCA tumor progression and clinical prognosis. Under the general trend of early diagnosis and precision medicine in the future, the results of this study are expected to provide certain guidance for individualized treatment.
Conclusion
Overall, the distinction between eosinophil-based risk scoring schemes helps predict clinical outcomes, genetic mutations, TIME heterogeneity, and treatment response. Nonetheless, these findings require further experimental and clinical validation in different centers and larger cohorts.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
CX, LS, ZK, CL, and DP designed this work. PH, YY, YL, JG, XL, and NL analyzed the data. JL and LS wrote this manuscript. ZK and PH edited and revised the manuscript. All authors approved this manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.963455/full#supplementary-material
Abbreviations
BLCA, bladder urothelial carcinoma; TIME, tumor immune microenvironment; WGCNA, weighted gene co-expression network analysis; RS, risk score; HRG, high-risk group; LRG, low-risk group; ROC, receiver operating characteristic; TMB, tumor mutational burden; IPS, Immune Phenomenon Score; IC50s, Half-maximal inhibitory concentrations; ICI, immune checkpoint inhibitors; MIBC, muscle-invasive bladder cancer; NMIBC, non-muscle invasive bladder cancer; ORR, objective response rate; CRR, complete response rate; TCGA-BLCA, The Cancer Genome Atlas-Bladder urothelial carcinoma; TCGA, The Cancer Genome Atlas; TIC, tumor-infiltrating immune cell; OS, overall survival; GSVA, Applied Gene Set Variation Analysis; IPS, Immune Phenomenon Score; GDSC, Genomics of Cancer Drug Sensitivity; SMG, Significantly mutated gene.
References
Bellmunt, J., de Wit, R., Vaughn, D., Fradet, Y., Lee, J., Fong, L., et al. (2017). Pembrolizumab as Second-Line Therapy for Advanced Urothelial Carcinoma. N. Engl. J. Med. 376 (11), 1015–1026. doi:10.1056/NEJMoa1613683
Blanche, P., Dartigues, J., and Jacqmin-Gadda, H. (2013). Estimating and Comparing Time-dependent Areas under Receiver Operating Characteristic Curves for Censored Event Times with Competing Risks. Stat. Med. 32 (30), 5381–5397. doi:10.1002/sim.5958
Buller, C. W., Mathew, P. A., and Mathew, S. O. (2020). Roles of NK Cell Receptors 2B4 (CD244), CS1 (CD319), and LLT1 (CLEC2D) in Cancer. Cancers 12 (7), 1755. doi:10.3390/cancers12071755
Burr, M., Sparbier, C., Chan, Y., Williamson, J., Woods, K., Beavis, P., et al. (2017). CMTM6 Maintains the Expression of PD-L1 and Regulates Anti-tumour Immunity. Nature 549 (7670), 101–105. doi:10.1038/nature23643
Chan, T., Yarchoan, M., Jaffee, E., Swanton, C., Quezada, S., Stenzinger, A., et al. (2019). Development of Tumor Mutation Burden as an Immunotherapy Biomarker: Utility for the Oncology Clinic. Ann. Oncol. 30 (1), 44–56. doi:10.1093/annonc/mdy495
Charoentong, P., Finotello, F., Angelova, M., Mayer, C., Efremova, M., Rieder, D., et al. (2017). Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 18 (1), 248–262. doi:10.1016/j.celrep.2016.12.019
Chen, M., Ni, J., Chang, H. C., Lin, C. Y., Muyan, M., and Yeh, S. (2009). CCDC62/ERAP75 Functions as a Coactivator to Enhance Estrogen Receptor Beta-Mediated Transactivation and Target Gene Expression in Prostate Cancer Cells. Carcinogenesis 30 (5), 841–850. doi:10.1093/carcin/bgn288
Chen, S., Wang, C., Su, X., Dai, X., Li, S., and Mo, Z. (2021). KCNN4 Is a Potential Prognostic Marker and Critical Factor Affecting the Immune Status of the Tumor Microenvironment in Kidney Renal Clear Cell Carcinoma. Transl. Androl. Urol. 10 (6), 2454–2470. doi:10.21037/tau-21-332
Cheng, S. A., Huang, X., Jiang, L., Liang, Q. L., Hong, X. C., Yang, H. X., et al. (2022). Expression of DNM3 Is Associated with Good Outcome in Colorectal Cancer. Open Med. Wars. Pol. 17 (1), 221–226. doi:10.1515/med-2022-0420
Chism, D. D. (2017). Urothelial Carcinoma of the Bladder and the Rise of Immunotherapy. J. Natl. Compr. Canc Netw. 15 (10), 1277–1284. doi:10.6004/jnccn.2017.7036
Costello, R., O'Callaghan, T., and Sébahoun, G. (2005). Eosinophils and Antitumour Response. Rev. Med. Interne 26 (6), 479–484. doi:10.1016/j.revmed.2005.02.013
Cumberbatch, M., Jubber, I., Black, P., Esperto, F., Figueroa, J., Kamat, A., et al. (2018). Epidemiology of Bladder Cancer: A Systematic Review and Contemporary Update of Risk Factors in 2018. Eur. Urol. 74 (6), 784–795. doi:10.1016/j.eururo.2018.09.001
De Palma, M., Biziato, D., and Petrova, T. J. NrC. (2017). Microenvironmental Regulation of Tumour Angiogenesis. Nat. Rev. Cancer 17 (8), 457–474. doi:10.1038/nrc.2017.51
Domae, S., Nakamura, Y., Nakamura, Y., Uenaka, A., Wada, H., Nakata, M., et al. (2009). Identification of CCDC62-2 as a Novel Cancer/testis Antigen and its Immunogenicity. Int. J. Cancer 124 (10), 2347–2352. doi:10.1002/ijc.24220
Fa, J. (2021). Dynamin 3 Overexpression Suppresses the Proliferation, Migration and Invasion of Cervical Cancer Cells. Oncol. Lett. 22 (1), 524. doi:10.3892/ol.2021.12785
Furbert-Harris, P., Laniyan, I., Harris, D., Dunston, G., Vaughn, T., Abdelnaby, A., et al. (2003). Activated Eosinophils Infiltrate MCF-7 Breast Multicellular Tumor Spheroids. Anticancer Res. 23, 71–78.
George, S., Miao, D., Demetri, G., Adeegbe, D., Rodig, S., Shukla, S., et al. (2017). Loss of PTEN Is Associated with Resistance to Anti-PD-1 Checkpoint Blockade Therapy in Metastatic Uterine Leiomyosarcoma. Immunity 46 (2), 197–204. doi:10.1016/j.immuni.2017.02.001
Gibney, G. T., Weiner, L. M., and Atkins, M. B. (2016). Predictive Biomarkers for Checkpoint Inhibitor-Based Immunotherapy. Lancet Oncol. 17 (12), e542–e551. doi:10.1016/s1470-2045(16)30406-5
Grisaru-Tal, S., Itan, M., Klion, A. D., and Munitz, A. (2020). A New Dawn for Eosinophils in the Tumour Microenvironment. Nat. Rev. Cancer 20 (10), 594–607. doi:10.1038/s41568-020-0283-9
Havel, J., Chowell, D., and Chan, T. J. NrC. (2019). The Evolving Landscape of Biomarkers for Checkpoint Inhibitor Immunotherapy. Nat. Rev. Cancer 19 (3), 133–150. doi:10.1038/s41568-019-0116-x
Hessmann, E., Buchholz, S., Demir, I., Singh, S., Gress, T., Ellenrieder, V., et al. (2020). Microenvironmental Determinants of Pancreatic Cancer. Physiol. Rev. 100 (4), 1707–1751. doi:10.1152/physrev.00042.2019
Hurst, C., and Knowles, M. A. (2022). Mutational Landscape of Non-muscle-invasive Bladder Cancer. Urologic Oncol. Seminars Orig. Investigations 40, 295–303. doi:10.1016/j.urolonc.2018.10.015
Iasonos, A., Schrag, D., Raj, G. V., and Panageas, K. S. (2008). How to Build and Interpret a Nomogram for Cancer Prognosis. Jco 26 (8), 1364–1370. doi:10.1200/jco.2007.12.9791
Jing, W., Wang, G., Cui, Z., Xiong, G., Jiang, X., Li, Y., et al. (2022). FGFR3 Destabilizes PD-L1 via NEDD4 to Control T-Cell-Mediated Bladder Cancer Immune Surveillance. Cancer Res. 82 (1), 114–129. doi:10.1158/0008-5472.Can-21-2362
Korbecki, J., Kojder, K., Simińska, D., Bohatyrewicz, R., Gutowska, I., Chlubek, D., et al. (2020). CC Chemokines in a Tumor: A Review of Pro-cancer and Anti-cancer Properties of the Ligands of Receptors CCR1, CCR2, CCR3, and CCR4. Int J Mol Sci. 21 (21), 8412. doi:10.3390/ijms21218412
Kruger, S., Ilmer, M., Kobold, S., Cadilha, B., Endres, S., Ormanns, S., et al. (2019). Advances in Cancer Immunotherapy 2019 - Latest Trends. J. Exp. Clin. Cancer Res. 38 (1), 268. doi:10.1186/s13046-019-1266-0
Langfelder, P., and Horvath, S. (2008). WGCNA: an R Package for Weighted Correlation Network Analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Lossos, I., Czerwinski, D., Alizadeh, A., Wechser, M., Tibshirani, R., Botstein, D., et al. (2004). Prediction of Survival in Diffuse Large-B-Cell Lymphoma Based on the Expression of Six Genes. N. Engl. J. Med. 350 (18), 1828–1837. doi:10.1056/NEJMoa032520
Lv, F., Huang, W., and Wang, Y. (2021). The CYP19A1 Rs700519 Polymorphism and Breast Cancer Susceptibility in China: A Case-Control Study and Updated Meta-Analysis. Genet. Test. Mol. biomarkers 25 (7), 486–495. doi:10.1089/gtmb.2021.0032
Ma, Y., Guan, L., Han, Y., Zhou, Y., Li, X., Liu, Y., et al. (2019). siPRDX2-elevated DNM3 Inhibits the Proliferation and Metastasis of Colon Cancer Cells via AKT Signaling Pathway. Cancer Manag. Res. 11, 5799–5811. doi:10.2147/cmar.S193805
Mathew, S. O., Chaudhary, P., Powers, S. B., Vishwanatha, J. K., and Mathew, P. A. (2016). Overexpression of LLT1 (OCIL, CLEC2D) on Prostate Cancer Cells Inhibits NK Cell-Mediated Killing through LLT1-Nkrp1a (CD161) Interaction. Oncotarget 7 (42), 68650–68661. doi:10.18632/oncotarget.11896
Narasimhan, A., Shahda, S., Kays, J. K., Perkins, S. M., Cheng, L., Schloss, K. N. H., et al. (2020). Identification of Potential Serum Protein Biomarkers and Pathways for Pancreatic Cancer Cachexia Using an Aptamer-Based Discovery Platform. Cancers 12 (12), 3787. doi:10.3390/cancers12123787
Olsen, E. G. J., and Spry, C. J. F. (1985). Relation between Eosinophilia and Endomyocardial Disease. Prog. Cardiovasc. Dis. 27, 241–254. doi:10.1016/0033-0620(85)90008-8
Özdemir, B. C., Siefker-Radtke, A. O., Campbell, M. T., and Subudhi, S. K. (2018). Current and Future Applications of Novel Immunotherapies in Urological Oncology: A Critical Review of the Literature. Eur. Urol. Focus 4 (3), 442–454. doi:10.1016/j.euf.2017.10.001
Pardoll, D. M. (2012). The Blockade of Immune Checkpoints in Cancer Immunotherapy. Nat. Rev. Cancer 12 (4), 252–264. doi:10.1038/nrc3239
Rhea, L. P., Mendez-Marti, S., Kim, D., and Aragon-Ching, J. B. (2021). Role of Immunotherapy in Bladder Cancer. Cancer Treat. Res. Commun. 26, 100296. doi:10.1016/j.ctarc.2020.100296
Rizvi, N., Hellmann, M., Snyder, A., Kvistborg, P., Makarov, V., Havel, J., et al. (2015). Mutational Landscape Determines Sensitivity to PD-1 Blockade in Non-small Cell Lung Cancer. Science 348 (6230), 124–128. doi:10.1126/science.aaa1348
Robertson, A., Kim, J., Al-Ahmadie, H., Bellmunt, J., Guo, G., Cherniack, A., et al. (2017). Compr. Mol. Charact. Muscle-Invasive Bladder Cancer 171 (3), 540–556. e25. doi:10.1016/j.cell.2017.09.007
Robertson, A., Kim, J., Al-Ahmadie, H., Bellmunt, J., Guo, G., Cherniack, A., et al. (2018). Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer. Cell 174 (4), 1033. doi:10.1016/j.cell.2018.07.036
Rosenberg, J. E., Hoffman-Censits, J., Powles, T., van der Heijden, M. S., Balar, A. V., Necchi, A., et al. (2016). Atezolizumab in Patients with Locally Advanced and Metastatic Urothelial Carcinoma Who Have Progressed Following Treatment with Platinum-Based Chemotherapy: a Single-Arm, Multicentre, Phase 2 Trial. Lancet 387 (10031), 1909–1920. doi:10.1016/s0140-6736(16)00561-4
Samstein, R., Lee, C., Shoushtari, A., Hellmann, M., Shen, R., Janjigian, Y., et al. (2019). Tumor Mutational Load Predicts Survival after Immunotherapy across Multiple Cancer Types. Nat. Genet. 51 (2), 202–206. doi:10.1038/s41588-018-0312-8
Shi, X., Zhou, S., Wang, Z., Zhou, Z., and Wang, Z. (2008). CYP1A1 and GSTM1 Polymorphisms and Lung Cancer Risk in Chinese Populations: a Meta-Analysis. Lung Cancer 59 (2), 155–163. doi:10.1016/j.lungcan.2007.08.004
Siegel, R., Miller, K., Fuchs, H., and Jemal, A. (2022). Cancer Statistics, 2022. CA A Cancer J. Clin. 72 (1), 7–33. doi:10.3322/caac.21708
Sylvester, R., Gontero, P., and Oddens, L. (2016). Reply to Stephen B. Williams and Ashish M. Kamat's Letter to the Editor re: Samantha Cambier, Richard J. Sylvester, Laurence Collette, et al. EORTC Nomograms and Risk Groups for Predicting Recurrence, Progression, and Disease-specific and Overall Survival in Non-Muscle-invasive Stage Ta-T1 Urothelial Bladder Cancer Patients Treated with 1-3 Years of Maintenance Bacillus Calmette-Guérin. Eur Urol 2016;69:60-9. Eur. Urol. 69 (6):e123-4. doi:10.1016/j.eururo.2016.01.055
Varn, F. S., Wang, Y., Mullins, D. W., Fiering, S., and Cheng, C. (2017). Systematic Pan-Cancer Analysis Reveals Immune Cell Interactions in the Tumor Microenvironment. Cancer Res. 77 (6), 1271–1282. doi:10.1158/0008-5472.Can-16-2490
von der Maase, H., Sengelov, L., Roberts, J., Ricci, S., Dogliotti, L., Oliver, T., et al. (2005). Long-term Survival Results of a Randomized Trial Comparing Gemcitabine Plus Cisplatin, with Methotrexate, Vinblastine, Doxorubicin, Plus Cisplatin in Patients with Bladder Cancer. J. Clin. Oncol. 23 (21), 4602–4608. doi:10.1200/jco.2005.07.757
Wang, C., Chen, S., Li, S., and Mi, H. (2022). A Prognostic Model for Predicting Tumor Mutation Burden and Tumor-Infiltrating Immune Cells in Bladder Urothelial Carcinoma. Front. Genet. 13, 708003. doi:10.3389/fgene.2022.708003
Webster, R. M. (2014). The Immune Checkpoint Inhibitors: where Are We Now? Nat. Rev. Drug Discov. 13 (12), 883–884. doi:10.1038/nrd4476
Wu, H., Ouyang, Q., Tian, C., Xie, N., Cao, M., and Shui, Z. R. (2017). Cytochrome P450 1A1 (CYP1A1) Gene Polymorphisms and Susceptibility to Breast Cancer: a Meta-Analysis in the Chinese Population. Clin. Lab. 63 (1), 67–72. doi:10.7754/Clin.Lab.2016.160535
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612
Keywords: eosinophil, bladder urothelial carcinoma, WGCNA, risk score, immunotherapy
Citation: Xu C, Song L, Peng H, Yang Y, Liu Y, Pei D, Guo J, Liu N, Liu J, Li X, Li C and Kang Z (2022) Clinical Eosinophil-Associated Genes can Serve as a Reliable Predictor of Bladder Urothelial Cancer. Front. Mol. Biosci. 9:963455. doi: 10.3389/fmolb.2022.963455
Received: 07 June 2022; Accepted: 23 June 2022;
Published: 22 July 2022.
Edited by:
Tao Xu, Anhui Medical University, ChinaReviewed by:
Di Shen, NovaRock Biotherapeutics Ltd., United StatesJin Zhang, University of Mississippi Medical Center, United States
Copyright © 2022 Xu, Song, Peng, Yang, Liu, Pei, Guo, Liu, Liu, Li, Li and Kang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhengjun Kang, kzj7153@sina.com; Chen Li, chen.li.scholar@gmail.com; Xiaoyong Li, LixyGDYXWK64@163.com; Jiabang Liu, jiabangliu6@gmail.com
†These authors have contributed equally to this work