- 1Department of Oncology, Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu University of Traditional Chinese Medicine, Chengdu, China
- 2Department of Medical Oncology, Cancer Center, West China Hospital, Sichuan University, Chengdu, China
- 3Department of Respiratory and Critical Care Medicine, Chongqing Traditional Chinese Medicine Hospital, Chongqing, China
- 4Department of Oncology, Taizhou People’s Hospital, Taizhou, China
Background: Bladder carcinoma (BLCA) is a common malignant tumor with high morbidity and mortality in the urinary system. Pyroptosis is a pattern of programmed cell death that is closely associated with progression of tumors. Therefore, it is significant to probe the expression of pyroptosis-related genes (PRGs) in BLCA.
Methods: The differentially expressed genes in normal and BLCA tissues were first obtained from the Cancer Genome Atlas (TCGA) database analysis, as well as PRGs from the National Center for Biotechnology Information (NCBI) database, intersecting to obtain differentially expressed pyroptosis-related genes (DEPRGs) in BLCA. With the construction of a prognostic model of pyroptosis by regression analysis, we derived and validated key genes, which were ascertained as a separate prognostic marker by individual prognostic and clinical relevance analysis. In addition, we gained six immune cells from the Tumor Immune Evaluation Resource (TIMER) website and analyzed the relationship between pyroptosis prognostic genes and immune infiltration.
Result: Our results revealed that 31 DEPRGs were available by comparing normal and BLCA tissues with |log2 (fold change, FC)| > 0.5 and FDR <0.05. Four key genes (CRTAC1, GSDMB, AIM2, and FOXO3) derived from the pyroptosis prognostic model were experimentally validated for consistent expression in BLCA patients. Following risk scoring, the low-risk group of BLCA patients had noticeably higher overall survival (OS) than the high-risk group (p < 0.001). Risk score was still an independent prognostic factor (HR = 1.728, 95% CI =1.289–2.315, p < 0.001). In addition, we found remarkable correlations among the expression of pyroptosis-related prognostic genes and the immune infiltration of CD4+ T cells, CD8+ T cells, B cells, dendritic cells, macrophages, and neutrophils.
Conclusion: Genes (CRTAC1, GSDMB, AIM2, and FOXO3) associated with pyroptosis are potential BLCA prognostic biomarkers that act as an essential part in the predictive prognosis of survival and immunotherapy of BLCA.
Introduction
Bladder carcinoma (BLCA) is the prevalent malignancy of the urinary tract and the fourth most common malignancy in men, with a high mortality and morbidity rate (1, 2). In 2020, the number of new BLCA cases was predicted to be about 81,400, and the number of new deaths was calculated to be about 17,980 in the United States (3). Of these, non-muscle invasive bladder carcinoma (NMIBC), which is prone to recurrence, accounts for approximately 75%, and muscle invasive bladder carcinoma (MIBC) and metastatic disease, which are more aggressive, account for 25% (4, 5). Currently, cisplatin-based chemotherapy remains the first-line treatment option for advanced BLCA, but the response period to chemotherapy is limited and the prognosis of patients remains dismal, with a 5-year relative overall survival rate of only 15% for metastatic disease (5, 6). Moreover, studies have indicated that the wild-type TP53 gene expression profile was resistant to neoadjuvant chemotherapy in MIBC tumors (7). The highly heterogeneous and genomically unstable nature of BLCA have also been identified by an increasing number of studies (8, 9). The genetic prognostic biomarkers of BLCA have been emerging in recent years. An example was minichromosome maintenance protein 5 (MCM5), a biomarker of cell proliferation, which is expressed in normal urinary epithelium only in cells within the basal proliferative lumen and can be used to diagnose BLCA (10). A clinical trial demonstrated that the detection of MCM5 in urine sediment was a biomarker for predicting BLCA recurrence and prognosis. However, it was not widely used in clinical practice because its false-positive rate was higher than that of cystoscopy (11). Thus, searching for gene expression characteristics of BLCA patients and identifying newly effective prognostic markers are essential procedures for personalized treatment strategies to improve clinical outcomes.
Pyroptosis is a novel non-apoptotic mode of programmed cell death (PCD) mediated by the gasdermin family, characterized by pore formation on the cytoplasmic membrane, cell swelling, and rupture, and accompanied by the efflux of cell contents. Currently, the gasdermin family comprises six human homologs, gasdermin A–D (GSDMA–D), gasdermin E (GSDME, alternatively known as DFNA5), and Pejvakin (PJVK, alternatively known as DFNB59) (12). Pyroptosis can be induced by both the classical inflammasome pathway, which is activation of caspase-1 cutting the GSDMD protein, and the non-classical inflammasome pathway, which is bacterial lipopolysaccharide binding to caspase-4/5/11 cutting the GSDMD protein (13, 14). In addition, it was recently found that the serine protease granzyme A (GZMA) in cytotoxic lymphocytes could access target cells via perforin and evoke cancer cells to undergo pyroptosis by cleaving activated GSDMB protein, as well as activate immune response (15). Meanwhile, emerging lines of evidence also suggested that selective delivery of active gasdermin protein to cancer cells or targeted cleavage of GSDME by granzymes can inspire potent antitumor immunity (16, 17). Produced during pyroptosis, IL-1β and IL-18 can impair the recruitment of neutrophils, whereas GSDMD acted as an anti-inflammatory agent to stimulate neutrophil death (18). Moreover, GSDMD also increased cleavage in activated CD8+ T cells and enhanced their cytolytic ability (19). Indeed, immune infiltration of natural killer (NK), M1 macrophages, CD8+ T, and type 1 T helper (Th1) cells typically represented a good prognosis, while high-level infiltration of neutrophils, M2 macrophages, MDSCs, Treg cells, and Th2 cells was frequently correlated with a poor prognosis (20, 21). Interestingly, Hou J et al. demonstrated that GSDMC expression was correlated with breast cancer survival prognosis (22).
These studies have demonstrated that pyroptosis was inextricably bound to tumorigenesis, invasion, and tumor immune microenvironment, implicating that induction of tumor cell pyroptosis would be a potential strategy to strengthen the efficacy of anti-cancer therapy. Moreover, it has been determined that pyroptosis-related genes (PRGs) have been developed as markers of tumor prognosis in a variety of carcinomas, such as lung cancer (23), hepatocellular carcinoma (24), and breast cancer (25), by establishing prognostic models of pyroptosis genes. However, there was relatively less research on its specific functions in BLCA. Hence, we systematically researched and validated the expression levels of genes associated with pyroptosis in normal bladder and BLCA tissues to study the prognostic value of these genes. We also analyzed the relationship between PRGs and tumor immune microenvironment by the TIMER website.
Materials and Methods
Data Collection
As shown in the flowchart (Figure 1), search for the word “pyroptosis” AND “Homo sapiens” on the NCBI database (https://www.ncbi.nlm.nih.gov/gene/) to retrieve PRGs. About 80 pyroptosis genes were screened by this method. We obtained transcriptome data and relevant clinical profile of BLCA and normal samples from the TCGA database (https://portal.gdc.cancer.gov/).
Figure 1 Workflow chart. Detailed data analysis workflow chart. Create BioRender.com.
Enrichment Analysis of Differentially Expressed Pyroptosis-Related Genes
The condition |log2 (fold change, FC)| > 0.5 and FDR < 0.05 was established with R software to calculate the DEPRGs of BLCA vs. normal samples. Biological properties were analyzed by means of functional analysis, including Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). We conducted biological process (BP), cellular component (CC), and molecular function (MF) enrichment analysis of DEPRGs using the R package “clusterprofler”. The same approach was also used for the enrichment analysis of KEGG. The p-value of less than 0.05 was statistically valid.
Construction of Pyroptosis Prognosis Model Based on Pyroptosis-Related Genes
The DEPRGs expression signatures were intersected from the TCGA and NCBI databases. Univariate Cox regression analysis was applied to probe and identify the DEPRGs associated to survival of BLCA patients with a threshold value of p < 0.05 by the R package “survival”. Subsequently, the pyroptosis prognostic mode was structured and optimized to filter the key genes. Next, the risk score of individual patients was computed by the following equation: risk score = ∑X J * coef J, where coef J means coefficient, and X J means correlative expression of each DEPRGs normalized by Z-score. Finally, the median risk score categorized BLCA cases into two groups: high risk and low risk.
Validation of the Key PRGs by Immunohistochemistry
To validate the pyroptosis prognosis model, we collected BLCA samples and adjacent non-cancerous bladder tissue from BLCA patients with primary surgical resection at Taizhou People’s Hospital. The Clinical Research Ethics Committee of Taizhou People’s Hospital has confirmed the tissues used in this study. Paraffin-embedded sections were stained for IHC analysis at 4 µm thickness. IHC was conducted with antibodies for CRTAC1 (1:200, bs-12948R, bioss), GSDMB (1:200, 12885-1-AP, Proteintech), AIM2 (1:200, PB9683, Boster), and FOXO3 (1:200, PB9196, Boster). Immunohistochemistry analysis was used following the manufacturer’s protocol. Then, two experienced pathologists, blinded to detailed information, would randomly select five regions and evaluate the immunostaining score.
Prognostic Mode Evaluation and Clinical Relevance Assessment
Kaplan–Meier (K-M) survival and risk profile were depicted through the R package comparing the survival of the two groups. Next, we screened for separate prognostic factors by assessing the association between clinicopathological features and risk scores with the R package “Survival”. Finally, the ROC curve was calculated with the R package “survivalROC”, where the ROC area under the curve (AUC) value was analyzed to ascertain the prediction accuracy of the model. For clinical relevance, the scatter plot was pictured with the R package “beeswarm”. According to clinicopathological characteristics, key pyroptosis genes and risk value were analyzed to forecast the biological markers for the prognosis of BLCA.
Analysis of Immune Cells Infiltration in PRGs of BLCA Patients
Systematic assessment of various immune cell infiltrations and their clinical influences can be conducted by TIMER (https://cistrome.shinyapps.io/timer/) (26). The relevance between PRG levels and immune cells infiltration was valued using the “Gene module”. For clinical outcomes, relevance to immune cells infiltration and PRG expression would be estimated by “Survival module”.
Results
Identification of DEPRGs Between Normal and BLCA Tissues
Gene expression profiles of 19 normal and 414 BLCA tissue samples and consequent clinical information were derived from the TCGA. Based on the search terms “pyroptosis” AND “Homo sapiens” in the Gene database, we identified 80 PRGs. Then, we found that 31 DEPRGs (Additional file 1) met the criteria of FDR < 0.05 and |log2(FC)| > 0.5 by R language screening, including 17 upregulated genes (CTSV, CAPN1, GSDMB, MIR25, PYCARD, TFAM, GSDMD, AIM2, NLRP7, IL36G, AGER, CDKN2B-AS1, METTL3, CASP3, TREM2, CASP8, and BSG) and 14 downregulated genes (VIM, PRDM1, TUBB6, MEG3, PECAM1, TXNIP, CRTAC1, KLF3-AS1, NEK7, NLRP3, EEF2K, FOXO3, NFE2L2, and NLRP1) (Figures 2A, B).
Functional Enrichment Analysis of DEPRGs
We scanned and analyzed the above 31 DEPRGs in BLCA with GO and KEGG to explore their potential signaling pathways. As shown in Figure 3A, BP indicated that DEPRGs were mainly enriched in the positive regulation of IL-1 and IL-1β production, response to LPS, etc. Among the CC, the most significantly altered pathway was the inflammasome complex. As for MF, cysteine-type endopeptidase activity, scaffold protein binding, and death receptor binding were significantly enhanced. The top 10 enriched signaling pathways with significant relevance in KEGG analysis were in the circle diagram (Figure 3B), including the NLR signaling pathway, apoptosis, necroptosis, and the C-type lectin receptor signaling pathway. Briefly, the analysis indicated a significant relevance between DEPRGs and the development of BLCA.
Figure 3 GO and KEGG enrichment analysis. (A) GO analysis of 31 DEPRGs. (B) The top 10 most enriched KEGG pathways analysis of 31 DEPRGs.
Construction and Evaluation of Pyroptosis Prognosis Model in BLCA
Based on univariate Cox regression analysis, we gained 12 genes associated with prognosis and structured a pyroptosis prognostic model (Figure 4A). The differential expression of 12 genes was displayed in the box plot of Figure 4B. Furthermore, the model was optimized by multivariate Cox regression analysis, and four key genes (CRTAC1, GSDMB, AIM2, and FOXO3) were obtained (Table 1). A heat map indicated the gene expression profiles of high- and low-risk groups in normal and BLCA (Figure 5A). Using K-M survival curve analysis, the 3- and 5-year survival rates for patients in the high- and low-risk group were 34.8% and 65.2%, and 31.1% and 54.8%, respectively. In a word, patients in group L had significantly higher survival rates relative to group H (p < 0.001) (Figure 5B). Risk score was computed for each patient based on the expression levels and risk factors of PRGs, and as shown in the risk curve results, the higher the risk score was, the shorter the survival time would be for BLCA patients (Figures 5C, D). It was suggested that these four DEPRGs may be potential prognostic indicators for BLCA.
Figure 4 Regression analysis to select pyroptosis genes related to prognosis of BLCA. (A) Forest map of pyroptosis genes analyzed by univariate cox regression. (B) Boxplot of pyroptosis genes analyzed by LASSO regression. “N” stands for “Normal” and “T” stands for “Tumor”.
Figure 5 Structure of pyroptosis prognosis model in BLCA. (A) Heatmap displaying the expression patterns of those 4 key genes (B) K-M survival curve analysis. (C) Risk score distribution for BLCA cases. (D) OS for BLCA patients.
Expression Level of the Key PRGs in BLCA Tissues
IHC analysis confirmed that CRTAC1 and FOXO3 proteins were highly expressed in adjacent non-cancerous bladder tissues, 59.26% and 51.85%, respectively, while GSDMB and AIM2 proteins were highly expressed in BLCA tissues, 59.26% and 85.19%, respectively (Figure 6). Interestingly, it was consistent with the results of bioinformatics analysis, except that FOXO3 protein was also highly expressed in BLCA.
Figure 6 The expression protein of CRTAC1 FOXO3, GSDMB, and AIM2 in BLCA tissues and non-cancerous tissues (IHC).
Correlation Between Pyroptosis-Related Prognostic Genes and Clinicopathological Characteristics
The relationship between clinicopathological features and hazard scores was investigated by single and multifactorial Cox independent prognostic analyses. Univariate Cox regression analysis revealed that age, stage, T stage, N stage, and risk score were remarkably associated with BLCA OS (Figure 7A). The risk score was retained after multifactorial Cox regression analysis (Figure 7B). Furthermore, the accuracy of the model was assessed by combining risk scores with clinicopathological characteristics to picture ROC curves (Figure 7C). The larger the area under the ROC curve was, the more precise the prognostic model was structured.
Figure 7 Evaluation of pyroptosis prognosis model in BLCA. (A) Univariate cox regression analysis. (B) Multiple cox independent prognosis analysis. (C) ROC curve of clinicopathological features.
Furthermore, we investigated the association of four prognostic genes (CRTAC1, GSDMB, AIM2, and FOXO3) and risk scores with clinical characteristics by independent samples t-test. As demonstrated in the scatter plot of Figure 8, decreased expression of CRTAC1 and GSDMB was significantly associated with advanced pathological T3–4, N1–3, and M1 stages, stages III and IV, and high grade (Figures 8A–J). Instead, high expression of FOXO3 and AIM2 was notably correlated with high grade, and overexpression of FOXO3 was also strongly correlated with stages III and IV and T3–4 (Figures 8K–N). Finally, we observed that higher risk values occurred in advanced pathological T3–4, late N1–3 stage, and higher stage (III and IV) (Figures 8O–R).
Figure 8 The clinical relevance analysis in BLCA. (A–E) CRTAC1 expression in ages, grade, pathological stages, pathological T and N stages. (F–J) GSDMB expression in grade, pathological stages, pathological T, N, and M stages. (K–M) FOXO3 expression in grade, pathological stages, pathological T stages. (N) AIM2 expression in grade. (O–R) Risk scores expressed in grades of pathological T and N stages.
Immune Cells Infiltration of Pyroptosis-Related Prognostic Genes in BLCA
It was demonstrated that pyroptosis participated in the modulation of the inflammatory response and tumor immune microenvironment. Therefore, we proceeded to comprehensively probe the association between DEPRGs in BLCA and immune cells infiltration. The results indicated that CRTAC1 expression was negatively relevant with the infiltration of CD8+ T cells (p = 2.51e−05), neutrophils (p = 1.89e−05), and dendritic cells (p = 2.70e−10) and positively relevant with the infiltration of B cells (p = 2.01e−02) (Figure 9A). Similarly, the expression of GSDMB was negatively associated with the infiltration of CD8+ T cells (p =7.25e−10), macrophages (p = 8.86e−05), and dendritic cells (p =3.47e−05) and positively associated with the infiltration of B cells (p = 6.05e−04) and CD4+ T cells (p = 3.05e−02) (Figure 9B). There was a positive relevance between FOXO3 expression and CD4+ T cells, CD8+ T cells, B cells, macrophages, dendritic cells, and neutrophils, all p < 0.05 (Figure 9C). In addition, there was a negative relevance between AIM2 expression and the infiltration of macrophages (p = 5.90e−04) and a positive association between AIM2 expression and the infiltration of CD8+ T cells (p = 1.06e−05), CD4+ T cells (p = 6.06e−09), neutrophils (p = 5.58e−24), and dendritic cells (p = 1.05e−23) (Figure 9D). Interestingly, FOXO3 and AIM2 showed opposite correlations with macrophages, which may be related to the lack of detailed typing of macrophages (M1 macrophages behaved as tumor suppressors, whereas M2 macrophages behaved as tumor promoters). In short, the infiltration of six types of immune cells (CD4+ T cells, CD8+ T cells, B cells, macrophages, dendritic cells, and neutrophils) was remarkably relevant with the expression of pyroptosis-related prognostic genes.
Figure 9 The relevance between pyroptosis-related prognostic genes and immune cell infiltration. The correlation between the abundance of immune cells and the expression of (A) CRTAC1, (B) GSDMB, (C) FOXO3, (D) AIM2 in BLCA.
Discussion
Pyroptosis is a lysis form of programmed death triggered by inflammasome and can rapidly release cellular contents such as pro-inflammatory molecules (IL-1β and IL-18) and antigens into the tumor microenvironment to activate the immune response of the body (27). A recent study has identified GSDM-mediated pyroptosis as a mechanism of cytotoxic lymphocyte killing to potentiate antitumor immunity (15). It was also confirmed that GSDME-mediated pyroptosis could promote T-cell recruitment through the release of relevant cytokines and trigger the activation of T cells to stimulate T-cell infiltration into tumor tissues (28). However, GSDME was often silenced in tumor tissues with promoter methylation, thus avoiding the occurrence of chemotherapeutic drug-induced pyroptosis and resisting antitumor effects (27, 29, 30). Hence, it is worthwhile to clarify the expression of PRGs in tumor tissues, and inducing pyroptosis in tumor cells would be possible and deemed as a robust method to promote the effect of oncotherapy.
With the development of bioinformatics and high-throughput sequencing technologies, it has become a reality to predict the prognosis of BLCA by assessing their risk level. For example, the risk models constructed based on m6A-related non-coding RNAs (31), pyroptosis-related non-coding RNAs (32), and ferroptosis-related genes (33), respectively, have evaluated the prognosis of BLCA patients or the response to immunotherapy. These models were constructed from different mechanisms to predict the prognosis of BLCA, yet unfortunately, the prognostic relationship and immune infiltration characteristics of pyroptosis-related mRNAs with BLCA patients were not considered.
In our study, we first explored the DEPRGs in normal and BLCA tissues, in which 31 genes were differentially expressed compared to normal tissue. These data suggested that differentially expressed genes may exert an underestimated role in the tumorigenesis and progression of BLCA. Then, we highlighted the function of DEPRGs using GO and KEGG pathway enrichment analysis. Predictably, we discovered that these genes were primarily relevant to the NLR signaling pathway, necroptosis, and regulation of IL-1β production. NLRs played an important regulatory role in inflammation-associated tumorigenesis. NLRs can specifically and sensitively recognize pathogen-associated molecular patterns, and activate caspase-1 or caspase-11, enhancing the production and release of IL-1β and IL-18 to recruit immune cells and promote the pyroptosis process (34, 35). Moreover, a study indicated that activation of the NF-κB/NLRP3/Caspase-1 pyroptosis signaling pathway exerted antitumor effects by increasing GSDMD levels and inhibiting cell proliferation and epithelial–mesenchymal transition progression (36). Therefore, DEPRGs may be potential targets for therapeutic and prognostic markers of BLCA.
Next, the expression of four key genes (CRTAC1, GSDMB, AIM2, and FOXO3) has been confirmed to be consistent with the bioinformatics analysis by IHC. However, the absence of difference in FOXO3 expression in cancerous and non-cancerous tissue may be related to both small sample size and the possibility that the FOXO3 gene was not transcribed in some adjacent non-cancerous bladder tissue. After that, they were screened to be associated with prognostic survival in BLCA by univariate and multifactorial Cox regression analysis. In addition, patient risk scores were calculated from mRNA expression levels and risk factors of each patient to obtain risk curves and K-M survival curves. Single- and multi-factor independent prognostic and clinical relevance analyses were conducted using the pyroptosis prognostic model to identify the relationship between gene expression levels and clinical prognosis. It was indicated that these genes were correlated with age, sex, grade, stage, and TNM subtype in the clinicopathological features (p < 0.05) and could be treated as a predictor of BLCA prognosis.
Cartilage acidic protein 1 (CRTAC1) is a glycosylated extracellular matrix molecule of human articular cartilage secreted by chondrocytes (37). They used gene expression data analysis from public databases to find that CRTAC1 expression was lower in lung adenocarcinoma (LUAD) than in normal tissues, demonstrating that the CRTAC1 gene was a protective gene in LUAD, and the high expression of CRTAC1 gene had a better prognosis (38). Equally, our study indicated that high expression of CRTAC1 was associated with earlier T (T1–2), N (N0) subtype, earlier stage (I and II), and low grades in BLCA patients. It suggested that the CRTAC1 gene could also be regarded as a protective gene in BLCA. Presently, there were no studies of CRTAC1 in BLCA or the relationship between CRTAC1 and pyroptosis, but our studies may provide new findings on prognostic markers in BLCA. The mechanism of CRTAC1 also needs further validation in BLCA.
GSDMB, as a member of the GSDM family, has also been observed to resemble the pyroptosis feature (39). However, the mechanism of GSDMB-induced pyroptosis was not quite clear and consistent. It was found that GSDMB-induced pyroptosis may be achieved by caspase-1 cleavage of GSDMB protein in 293T cells (40, 41). In contrast, Chen et al. demonstrated that GSDMB induced a conformational change in caspase-4 protein, which, in turn, increased caspase-4 enzyme activity and accelerated the cleavage of GSDMD, consequently causing non-canonical pyroptosis (42). In addition, the association between the high or low expression of GSDMB in normal and tumor tissues and prognosis was also controversial. Based on several studies (39, 43–45), hyperexpression of GSDMB in hepatocellular carcinoma, gastric cancer, and cervical cancer was considered as an oncogene that would appear to be involved in cancer invasion and metastasis. However, Shao Feng’s team indicated that CTLs could mediate target cell pyroptosis through GZMA cleavage of GSDMB in target cells. Respectively, they engrafted CT26 cells and B16-F10 cells with recombinant human GSDMB into mice subcutaneously, and found that GZMA/GSDMB pathway-mediated pyroptosis strengthened tumor clearance in the mouse model by enhancing the therapeutic efficacy of antitumor immunity. Meanwhile, they revealed that high GSDMB expression was relevant with OS benefit in BLCA patients (15). This is consistent with our findings. Our research also indicated that GSDMB was highly expressed in BLCA and hyper-expression of GSDMB was associated with low risk, low-grade infiltration, early stage (I and II), no lymph node infiltration (N0), and no distant metastasis (M0) of BLCA patients, all p < 0.05. Thus, GSDMB could serve as a protective prognostic biomarker for BLCA.
The forkhead box O class (FOXO) family is a widely expressed transcription factor, mainly including FOXO1, FOXO3, FOXO4, and FOXO6. The FOXO family is involved in several biological processes that were related to cell metabolism, cell proliferation, and apoptosis (46). Among them, FOXO3, also known as FOXO3a, plays a significant role in the regulation of cellular processes by targeting the expression and activity of effector genes. Meanwhile, the FOXO3 activity can be inhibited by an imbalance between kinases and phosphatases, significantly affecting cellular processes and even inducing carcinogenesis (47). The expression of FOXO3 in various tumors exhibited variable pro-carcinogenic and anti-carcinogenic effects. Qian et al. demonstrated that high expression of FOXO3 was relevant with the progression of glioblastoma, suggesting a poor prognosis for patients with glioblastoma (48), the same result for triple-negative breast cancer (49) and hepatocellular carcinoma (50). However, downregulation of FOXO3 promoted DNMT3B overexpression, leading to tumor growth in lung cancer, while clinical findings also showed that there was a high DNMT3B, low FOXO3a, and high MDM2 expression associated with low survival rates in lung cancer patients (51). Moreover, transcriptional inhibition of the FOXO3 gene directly and downregulation of FOXO3 protein phosphorylation promoted proliferation and migration of uroepithelial cancer cells in vivo and in vitro (52). Overexpression of circular RNA FOXO3 suppressed the proliferation, migration, and invasion of BLCA cell lines as well as promoted their apoptosis (53, 54). In contrast to our findings, FOXO3 was an indicator of poor prognosis in BLCA because it was enriched in high-risk groups, and the mechanism remains to be investigated. The study showed that downregulation of FOXO3 could inhibit NLRP3 inflammasome-mediated endothelial cell pyroptosis (55). However, the relationship between FOXO3 and pyroptosis in tumors has been rarely reported. Therefore, it is essential to further explore the regulation of FOXO3 on tumor cell pyroptosis.
Absent in melanoma 2 (AIM2), a typical receptor for cytoplasmic DNA contained a positively charged HIN-200 structural domain that recognized the sugar-phosphate backbone of host ectopic negatively charged double-stranded DNA (dsDNA) forming an electrostatic attraction (56). Host ectopic dsDNA activated AIM2 inflammasomes that drove the secretion of proinflammatory cytokines IL-18 and IL-1β, inducing cellular pyroptosis (57). AIM2 exhibited both growth inhibitory effects in colon, hepatoma, nasopharyngeal, and lung cancers, and proliferation-promoting effects in cutaneous squamous cell carcinomas (SCCs) and endometrial cancers (57–60). Compared with normal tissues, components of NLRP3 and AIM2 inflammasomes were highly expressed in nasopharyngeal carcinoma tissues, which correlated with an increased chance of survival in nasopharyngeal carcinoma (58). On the other hand, the absence of AIM2 inhibited the growth of cutaneous SCC (59). Therefore, the role of AIM2 in tumorigenesis was bipartite, and its outcome depended on the type of cancer. According to our findings, AIM2 appeared to be an oncogenic gene, as it was highly expressed in BLCA tissues and significantly upregulated in high-grade BLCA.
Briefly, all four genes (CRTAC1, GSDMB, FOXO3, and AIM2) were prognostic indicators of BLCA pyroptosis, with CRTAC1 and GSDMB suggesting a better prognosis, and FOXO3 and AIM2 suggesting a poor prognosis.
PRGs worked through cell membrane pore formation to activate pyroptosis, enhance antitumor immunity, and inhibit tumor growth. Expression of pyroptosis genes strengthened phagocytosis of tumor cells by tumor-associated macrophages (TAM), as well as the count and function of tumor-infiltrating NK cells (17). However, TAM infiltration in BLCA disrupted the adaptive immune response and promoted tumor vascularization, which, in turn, led to tumor progression (61). Meanwhile, one study described that BLCA can be stratified as high or low risk on the basis of immune cells status (62). It was evident that the infiltration of immune cells was of great importance to the survival prognosis of BLCA. It was shown that BLCA recruited more T cells than normal bladder cells, and that infiltrating T cells increased the proliferation and invasion of BLCA cells by modulating the ERβ/c-MET or ERβ/IL-1/c-MET signaling pathways, promoting metastasis (63). Neutrophils may promote tumor aggressiveness and progression in BLCA owing to their stimulation of excessive inflammation, release of growth factors to malignant cells, and secretion of human neutrophil polypeptides (64). Within the present research, we indicated a remarkable relevance between the expression of PRGs (CRTAC1, GSDMB, FOXO3, and AIM2) and the infiltration of six types of immune cells, including CD4+ T cells, CD8+ T cells, B cells, macrophages, dendritic cells, and neutrophils, implying that these genes not only could serve as prognostic indicators, but also may reflect the immune status.
Chen et al. (65) have previously constructed a pyroptosis-related prognostic model to predict survival in BLCA patients, although their study contained only 53 pyroptosis genes compared with 80 pyroptosis genes in our study. Moreover, there was a lack of experimental validation in their study. Unfortunately, our study also has some limitations. It is unclear how these DEPRGs in the prognostic model we constructed modulate the underlying mechanisms of the BLCA process. A further exploration of their biological functions is required by elaborate experiments. Furthermore, this was a retrospective study, and we would like to conduct further well-designed prospective studies to confirm our findings.
Conclusion
In summary, our study has constructed a prognostic model based on pyroptosis in BLCA and identified four key prognostic genes (CRTAC1, GSDMB, FOXO3, and AIM2), which also reflected immune status. We expect that our findings will provide new insights into emerging immunotherapeutic targets, identify biomarkers to more precisely predict survival in BLCA patients, and provide data to support clinical formulation of personalized treatment plans.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the Clinical Research Ethics Committee of Taizhou People’s Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
Conception: CY, QG, and XY. Equally contributing to the paper, study design, data analysis, experiment execution, and drafting: JY and HL. Experiment execution: YW. Study supervision: PF and YZ. Final approval: all authors. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the Taizhou Science and Technology Support Plan (Social Development) project (TS201903), the China Postdoctoral Science Foundation (2019M663505), and the Postdoctoral Interdisciplinary Innovation Foundation, Sichuan University (No. 0040204153243).
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.
Acknowledgments
The authors thank Ping Lin, Jie Zhang, and Qin Li of the lab of experimental oncology for modification of the manuscript. The authors thank TCGA for their research network. The authors thank Jiao Liu of Taizhou People’s Hospital for providing sectioning techniques and Chengdu Xavier for providing immunohistochemical technical support. We gratefully acknowledge contributions from TIMER platforms and BioRender for modification of Figure 1.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.881860/full#supplementary-material
References
1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68(6):394–424. doi: 10.3322/caac.21492
2. Lenis AT, Lec PM, Chamie K, Mshs MD. Bladder Cancer: A Review. JAMA (2020) 324(19):1980–91. doi: 10.1001/jama.2020.17598
3. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590
4. Dy GW, Gore JL, Forouzanfar MH, Naghavi M, Fitzmaurice C. Global Burden of Urologic Cancers, 1990-2013. Eur Urol. (2017) 71(3):437–46. doi: 10.1016/j.eururo.2016.10.008
5. Nadal R, Bellmunt J. Management of Metastatic Bladder Cancer. Cancer Treat Rev (2019) 76:10–21. doi: 10.1016/j.ctrv.2019.04.002
6. Alifrangis C, McGovern U, Freeman A, Powles T, Linch M. Molecular and Histopathology Directed Therapy for Advanced Bladder Cancer. Nat Rev Urol. (2019) 16(8):465–83. doi: 10.1038/s41585-019-0208-0
7. Choi W, Porten S, Kim S, Willis D, Plimack ER, Hoffman-Censits J, et al. Identification of Distinct Basal and Luminal Subtypes of Muscle-Invasive Bladder Cancer With Different Sensitivities to Frontline Chemotherapy. Cancer Cell (2014) 25(2):152–65. doi: 10.1016/j.ccr.2014.01.009
8. Robertson AG, Kim J, Al-Ahmadie H, Bellmunt J, Guo G, Cherniack AD, et al. Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer. Cell (2018) 174(4):1033. doi: 10.1016/j.cell.2018.07.036
9. Kamoun A, de Reynies A, Allory Y, Sjodahl G, Robertson AG, Seiler R, et al. A Consensus Molecular Classification of Muscle-Invasive Bladder Cancer. Eur Urol. (2020) 77(4):420–33. doi: 10.1016/j.eururo.2019.11.011
10. Ferro M, La Civita E, Liotti A, Cennamo M, Tortora F, Buonerba C, et al. Liquid Biopsy Biomarkers in Urine: A Route Towards Molecular Diagnosis and Personalized Medicine of Bladder Cancer. J Pers Med (2021) 11(3):237. doi: 10.3390/jpm11030237
11. Roupret M, Gontero P, McCracken SRC, Dudderidge T, Stockley J, Kennedy A, et al. Diagnostic Accuracy of MCM5 for the Detection of Recurrence in Nonmuscle Invasive Bladder Cancer Followup: A Blinded, Prospective Cohort, Multicenter European Study. J Urol. (2020) 204(4):685–90. doi: 10.1097/JU.0000000000001084
12. Broz P, Pelegrin P, Shao F. The Gasdermins, a Protein Family Executing Cell Death and Inflammation. Nat Rev Immunol (2020) 20(3):143–57. doi: 10.1038/s41577-019-0228-2
13. Liu X, Lieberman J. A Mechanistic Understanding of Pyroptosis: The Fiery Death Triggered by Invasive Infection. Adv Immunol (2017) 135:81–117. doi: 10.1016/bs.ai.2017.02.002
14. Shi J, Zhao Y, Wang K, Shi X, Wang Y, Huang H, et al. Cleavage of GSDMD by Inflammatory Caspases Determines Pyroptotic Cell Death. Nature (2015) 526(7575):660–5. doi: 10.1038/nature15514
15. Zhou Z, He H, Wang K, Shi X, Wang Y, Su Y, et al. Granzyme A from cytotoxic lymphocytes cleaves GSDMB to trigger pyroptosis in target cells. Science (2020) 368(6494):eaaz7548. doi: 10.1126/science.aaz7548
16. Wang Q, Wang Y, Ding J, Wang C, Zhou X, Gao W, et al. A Bioorthogonal System Reveals Antitumour Immune Function of Pyroptosis. Nature (2020) 579(7799):421–6. doi: 10.1038/s41586-020-2079-1
17. Zhang Z, Zhang Y, Xia S, Kong Q, Li S, Liu X, et al. Gasdermin E Suppresses Tumour Growth by Activating Anti-Tumour Immunity. Nature (2020) 579(7799):415–20. doi: 10.1038/s41586-020-2071-9
18. Kambara H, Liu F, Zhang X, Liu P, Bajrami B, Teng Y, et al. Gasdermin D Exerts Anti-Inflammatory Effects by Promoting Neutrophil Death. Cell Rep (2018) 22(11):2924–36. doi: 10.1016/j.celrep.2018.02.067
19. Xi G, Gao J, Wan B, Zhan P, Xu W, Lv T, et al. GSDMD is Required for Effector CD8(+) T Cell Responses to Lung Cancer Cells. Int Immunopharmacol. (2019) 74:105713. doi: 10.1016/j.intimp.2019.105713
20. Fridman WH, Zitvogel L, Sautes-Fridman C, Kroemer G. The Immune Contexture in Cancer Prognosis and Treatment. Nat Rev Clin Oncol (2017) 14(12):717–34. doi: 10.1038/nrclinonc.2017.101
21. Becht E, Giraldo NA, Germain C, de Reynies A, Laurent-Puig P, Zucman-Rossi J, et al. Immune Contexture, Immunoscore, and Malignant Cell Molecular Subgroups for Prognostic and Theranostic Classifications of Cancers. Adv Immunol (2016) 130:95–190. doi: 10.1016/bs.ai.2015.12.002
22. Hou J, Zhao R, Xia W, Chang CW, You Y, Hsu JM, et al. PD-L1-Mediated Gasdermin C Expression Switches Apoptosis to Pyroptosis in Cancer Cells and Facilitates Tumour Necrosis. Nat Cell Biol (2020) 22(10):1264–75. doi: 10.1038/s41556-020-0575-z
23. Lin W, Chen Y, Wu B, Chen Y, Li Z. Identification of the Pyroptosisrelated Prognostic Gene Signature and the Associated Regulation Axis in Lung Adenocarcinoma. Cell Death Discovery (2021) 7(1):161. doi: 10.1038/s41420-021-00557-2
24. Deng M, Sun S, Zhao R, Guan R, Zhang Z, Li S, et al. The Pyroptosis-Related Gene Signature Predicts Prognosis and Indicates Immune Activity in Hepatocellular Carcinoma. Mol Med (2022) 28(1):16. doi: 10.1186/s10020-022-00445-0
25. Wu J, Zhu Y, Luo M, Li L. Comprehensive Analysis of Pyroptosis-Related Genes and Tumor Microenvironment Infiltration Characterization in Breast Cancer. Front Immunol (2021) 12:748221. doi: 10.3389/fimmu.2021.748221
26. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77(21):e108–e10. doi: 10.1158/0008-5472.CAN-17-0307
27. Fan JX, Deng RH, Wang H, Liu XH, Wang XN, Qin R, et al. Epigenetics-Based Tumor Cells Pyroptosis for Enhancing the Immunological Effect of Chemotherapeutic Nanocarriers. Nano Lett (2019) 19(11):8049–58. doi: 10.1021/acs.nanolett.9b03245
28. Peng Z, Wang P, Song W, Yao Q, Li Y, Liu L, et al. GSDME Enhances Cisplatin Sensitivity to Regress non-Small Cell Lung Carcinoma by Mediating Pyroptosis to Trigger Antitumor Immunocyte Infiltration. Signal Transduct Target Ther (2020) 5(1):159. doi: 10.1038/s41392-020-00274-9
29. Wang Y, Gao W, Shi X, Ding J, Liu W, He H, et al. Chemotherapy Drugs Induce Pyroptosis Through Caspase-3 Cleavage of a Gasdermin. Nature (2017) 547(7661):99–103. doi: 10.1038/nature22393
30. Wang CJ, Tang L, Shen DW, Wang C, Yuan QY, Gao W, et al. The Expression and Regulation of DFNA5 in Human Hepatocellular Carcinoma DFNA5 in Hepatocellular Carcinoma. Mol Biol Rep (2013) 40(12):6525–31. doi: 10.1007/s11033-013-2581-8
31. Lu M, Zhan H, Liu B, Li D, Li W, Chen X, et al. N6-Methyladenosine-Related non-Coding RNAs are Potential Prognostic and Immunotherapeutic Responsiveness Biomarkers for Bladder Cancer. EPMA J (2021) 12(4):589–604. doi: 10.1007/s13167-021-00259-w
32. Lu Z, Tang F, Li Z, Lai Y, Lu Z, Zhang J, et al. Prognosis Risk Model Based on Pyroptosis-Related lncRNAs for Bladder Cancer. Dis Markers. (2022) 2022:7931393. doi: 10.1155/2022/7931393
33. Sun J, Yue W, You J, Wei X, Huang Y, Ling Z, et al. Identification of a Novel Ferroptosis-Related Gene Prognostic Signature in Bladder Cancer. Front Oncol (2021) 11:730716. doi: 10.3389/fonc.2021.730716
34. Lee C, Do HTT, Her J, Kim Y, Seo D, Rhee I. Inflammasome as a Promising Therapeutic Target for Cancer. Life Sci (2019) 231:116593. doi: 10.1016/j.lfs.2019.116593
35. Liu P, Lu Z, Liu L, Li R, Liang Z, Shen M, et al. NOD-Like Receptor Signaling in Inflammation-Associated Cancers: From Functions to Targeted Therapies. Phytomedicine (2019) 64:152925. doi: 10.1016/j.phymed.2019.152925
36. Tan YF, Wang M, Chen ZY, Wang L, Liu XH. Inhibition of BRD4 Prevents Proliferation and Epithelial-Mesenchymal Transition in Renal Cell Carcinoma via NLRP3 Inflammasome-Induced Pyroptosis. Cell Death Dis (2020) 11(4):239. doi: 10.1038/s41419-020-2431-2
37. Steck E, Braun J, Pelttari K, Kadel S, Kalbacher H, Richter W. Chondrocyte Secreted CRTAC1: A Glycosylated Extracellular Matrix Molecule of Human Articular Cartilage. Matrix Biol (2007) 26(1):30–41. doi: 10.1016/j.matbio.2006.09.006
38. Yu P, Tong L, Song Y, Qu H, Chen Y. Systematic Profiling of Invasion-Related Gene Signature Predicts Prognostic Features of Lung Adenocarcinoma. J Cell Mol Med (2021) 25(13):6388–402. doi: 10.1111/jcmm.16619
39. Hergueta-Redondo M, Sarrio D, Molina-Crespo A, Vicario R, Bernado-Morales C, Martinez L, et al. Gasdermin B Expression Predicts Poor Clinical Outcome in HER2-Positive Breast Cancer. Oncotarget (2016) 7(35):56295–308. doi: 10.18632/oncotarget.10787
40. Ding J, Wang K, Liu W, She Y, Sun Q, Shi J, et al. Pore-Forming Activity and Structural Autoinhibition of the Gasdermin Family. Nature (2016) 535(7610):111–6. doi: 10.1038/nature18590
41. Panganiban RA, Sun M, Dahlin A, Park HR, Kan M, Himes BE, et al. A Functional Splice Variant Associated With Decreased Asthma Risk Abolishes the Ability of Gasdermin B to Induce Epithelial Cell Pyroptosis. J Allergy Clin Immunol (2018) 142(5):1469–78.e2. doi: 10.1016/j.jaci.2017.11.040
42. Chen Q, Shi P, Wang Y, Zou D, Wu X, Wang D, et al. GSDMB Promotes non-Canonical Pyroptosis by Enhancing Caspase-4 Activity. J Mol Cell Biol (2019) 11(6):496–508. doi: 10.1093/jmcb/mjy056
43. Hergueta-Redondo M, Sarrio D, Molina-Crespo A, Megias D, Mota A, Rojo-Sebastian A, et al. Gasdermin-B Promotes Invasion and Metastasis in Breast Cancer Cells. PLoS One (2014) 9(3):e90099. doi: 10.1371/journal.pone.0090099
44. Carl-McGrath S, Schneider-Stock R, Ebert M, Rocken C. Differential Expression and Localisation of Gasdermin-Like (GSDML), a Novel Member of the Cancer-Associated GSDMDC Protein Family, in Neoplastic and non-Neoplastic Gastric, Hepatic, and Colon Tissues. Pathology (2008) 40(1):13–24. doi: 10.1080/00313020701716250
45. Sun Q, Yang J, Xing G, Sun Q, Zhang L, He F. Expression of GSDML Associates With Tumor Progression in Uterine Cervix Cancer. Transl Oncol (2008) 1(2):73–83. doi: 10.1593/tlo.08112
46. Yang T, Li Y, Zhao F, Zhou L, Jia R. Circular RNA Foxo3: A Promising Cancer-Associated Biomarker. Front Genet (2021) 12:652995. doi: 10.3389/fgene.2021.652995
47. Liu Y, Ao X, Ding W, Ponnusamy M, Wu W, Hao X, et al. Critical Role of FOXO3a in Carcinogenesis. Mol Cancer. (2018) 17(1):104. doi: 10.1186/s12943-018-0856-3
48. Qian Z, Ren L, Wu D, Yang X, Zhou Z, Nie Q, et al. Overexpression of FoxO3a Is Associated With Glioblastoma Progression and Predicts Poor Patient Prognosis. Int J Cancer. (2017) 140(12):2792–804. doi: 10.1002/ijc.30690
49. Rehman A, Kim Y, Kim H, Sim J, Ahn H, Chung MS, et al. FOXO3a Expression Is Associated With Lymph Node Metastasis and Poor Disease-Free Survival in Triple-Negative Breast Cancer. J Clin Pathol (2018) 71(9):806–13. doi: 10.1136/jclinpath-2018-205052
50. Ahn H, Kim H, Abdul R, Kim Y, Sim J, Choi D, et al. Overexpression of Forkhead Box O3a and Its Association With Aggressive Phenotypes and Poor Prognosis in Human Hepatocellular Carcinoma. Am J Clin Pathol (2018) 149(2):117–27. doi: 10.1093/ajcp/aqx132
51. Yang YC, Tang YA, Shieh JM, Lin RK, Hsu HS, Wang YC. DNMT3B Overexpression by Deregulation of FOXO3a-Mediated Transcription Repression and MDM2 Overexpression in Lung Cancer. J Thorac Oncol (2014) 9(9):1305–15. doi: 10.1097/JTO.0000000000000240
52. Wu WR, Lin JT, Pan CT, Chan TC, Liu CL, Wu WJ, et al. Amplification-Driven BCL6-Suppressed Cytostasis is Mediated by Transrepression of FOXO3 and Post-Translational Modifications of FOXO3 in Urinary Bladder Urothelial Carcinoma. Theranostics (2020) 10(2):707–24. doi: 10.7150/thno.39018
53. Wang C, Tao W, Ni S, Chen Q. Circular RNA Circ-Foxo3 Induced Cell Apoptosis in Urothelial Carcinoma via Interaction With miR-191-5p. Onco Targets Ther (2019) 12:8085–94. doi: 10.2147/OTT.S215823
54. Li Y, Qiao L, Zang Y, Ni W, Xu Z. Circular RNA FOXO3 Suppresses Bladder Cancer Progression and Metastasis by Regulating MiR-9-5p/TGFBR2. Cancer Manag Res (2020) 12:5049–56. doi: 10.2147/CMAR.S253412
55. Li P, Zhong X, Li J, Liu H, Ma X, He R, et al. MicroRNA-30c-5p Inhibits NLRP3 Inflammasome-Mediated Endothelial Cell Pyroptosis Through FOXO3 Down-Regulation in Atherosclerosis. Biochem Biophys Res Commun (2018) 503(4):2833–40. doi: 10.1016/j.bbrc.2018.08.049
56. Jin T, Perry A, Jiang J, Smith P, Curry JA, Unterholzner L, et al. Structures of the HIN Domain:DNA Complexes Reveal Ligand Binding and Activation Mechanisms of the AIM2 Inflammasome and IFI16 Receptor. Immunity (2012) 36(4):561–71. doi: 10.1016/j.immuni.2012.02.014
57. Sharma BR, Karki R, Kanneganti TD. Role of AIM2 Inflammasome in Inflammatory Diseases, Cancer and Infection. Eur J Immunol (2019) 49(11):1998–2011. doi: 10.1002/eji.201848070
58. Karki R, Kanneganti TD. Diverging Inflammasome Signals in Tumorigenesis and Potential Targeting. Nat Rev Cancer. (2019) 19(4):197–214. doi: 10.1038/s41568-019-0123-y
59. Farshchian M, Nissinen L, Siljamaki E, Riihila P, Piipponen M, Kivisaari A, et al. Tumor Cell-Specific AIM2 Regulates Growth and Invasion of Cutaneous Squamous Cell Carcinoma. Oncotarget (2017) 8(28):45825–36. doi: 10.18632/oncotarget.17573
60. Xue Y, Enosi Tuipulotu D, Tan WH, Kay C, Man SM. Emerging Activators and Regulators of Inflammasomes and Pyroptosis. Trends Immunol (2019) 40(11):1035–52. doi: 10.1016/j.it.2019.09.005
61. Miyake M, Hori S, Morizawa Y, Tatsumi Y, Nakai Y, Anai S, et al. CXCL1-Mediated Interaction of Cancer Cells With Tumor-Associated Macrophages and Cancer-Associated Fibroblasts Promotes Tumor Progression in Human Bladder Cancer. Neoplasia (2016) 18(10):636–46. doi: 10.1016/j.neo.2016.08.002
62. Qiu H, Hu X, He C, Yu B, Li Y, Li J. Identification and Validation of an Individualized Prognostic Signature of Bladder Cancer Based on Seven Immune Related Genes. Front Genet (2020) 11:12. doi: 10.3389/fgene.2020.00012
63. Tao L, Qiu J, Slavin S, Ou Z, Liu Z, Ge J, et al. Recruited T Cells Promote the Bladder Cancer Metastasis via Up-Regulation of the Estrogen Receptor Beta/IL-1/C-MET Signals. Cancer Lett (2018) 430:215–23. doi: 10.1016/j.canlet.2018.03.045
64. Thompson DB, Siref LE, Feloney MP, Hauke RJ, Agrawal DK. Immunological Basis in the Pathogenesis and Treatment of Bladder Cancer. Expert Rev Clin Immunol (2015) 11(2):265–79. doi: 10.1586/1744666X.2015.983082
Keywords: bladder carcinoma, pyroptosis, prognostic biomarkers, immune infiltration, Lasso-Cox regression
Citation: You J, Li H, Wei Y, Fan P, Zhao Y, Yi C, Guo Q and Yang X (2022) Novel Pyroptosis-Related Gene Signatures Identified as the Prognostic Biomarkers for Bladder Carcinoma. Front. Oncol. 12:881860. doi: 10.3389/fonc.2022.881860
Received: 23 February 2022; Accepted: 30 May 2022;
Published: 30 June 2022.
Edited by:
Yair Lotan, University of Texas Southwestern Medical Center, United StatesReviewed by:
Jian Zhou, Chinese Academy of Sciences (CAS), ChinaMatteo Ferro, European Institute of Oncology (IEO), Italy
Raj Bhanvadia, University of Texas Southwestern Medical Center, United States
Copyright © 2022 You, Li, Wei, Fan, Zhao, Yi, Guo and Yang. 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: Cheng Yi, eWljaGVuZzY4MzRAMTI2LmNvbQ==; Qing Guo, YmlnbGVtb25vbmVAMTYzLmNvbQ==; Xi Yang, eWFuZ3hpMDA3M0AxNjMuY29t
†These authors have contributed equally to this work