Development and validation of a kidney renal clear cell carcinoma prognostic model relying on pyroptosis-related LncRNAs-A multidimensional comprehensive bioinformatics exploration
European Journal of Medical Research volume 28, Article number: 341 (2023)
Renal cell carcinoma (RCC) is a malignant tumour that may develop in the kidney. RCC is one of the most common kinds of tumours of this sort, and its most common pathological subtype is kidney renal clear cell carcinoma (KIRC). However, the aetiology and pathogenesis of RCC still need to be clarified. Exploring the internal mechanism of RCC contributes to diagnosing and treating this disease. Pyroptosis is a critical process related to cell death. Recent research has shown that pyroptosis is a critical factor in the initiation and progression of tumour formation. Thus far, researchers have progressively uncovered evidence of the regulatory influence that long noncoding RNAs (lncRNAs) have on pyroptosis.
In this work, a comprehensive bioinformatics approach was used to produce a predictive model according to pyroptosis-interrelated lncRNAs for the purpose of predicting the overall survival and molecular immune specialties of patients diagnosed with KIRC. This model was verified from multiple perspectives.
First, we discovered pyroptosis-associated lncRNAs in KIRC patients using the TCGA database and a Sankey diagram. Then, we developed and validated a KIRC patient risk model based on pyroptosis-related lncRNAs. We demonstrated the grouping power of PLnRM through PCA and used PLnRM to assess the tumour immune microenvironment and response to immunotherapy. Immunological and molecular traits of diverse PLnRM subgroups were evaluated, as were clinical KIRC patient characteristics and predictive risk models. On this basis, a predictive nomogram was developed and analyzed, and novel PLnRM candidate compounds were identified. Finally, we investigated possible medications used by KIRC patients.
The results demonstrate that the model generated has significant value for KIRC in clinical practice.
There are approximately 210,000 new instances of renal cell carcinoma (RCC) diagnosed each year worldwide, accounting for between 2 and 3% of all tumour cases. The most prevalent kind of RCC is kidney renal clear cell carcinoma, which it accounts for approximately 85% of all RCC cases. Individuals diagnosed with KIRC almost always have poor prognosis, affecting both their health and life . Patients with early-stage kidney cancer can be treated effectively with surgery. Nonetheless, after major surgery, approximately 30% of patients experience recurrence or metastasis, with poor prognosis for overall survival (OS) . Typically, metastatic KIRC, a more advanced form of kidney cancer, is not entirely curable and has a short median survival. In recent years, immune checkpoint inhibitors, a new type of tumour therapy, have benefited some kidney cancer patients, especially PD-L1 and PD-1 inhibitors . In reality, immunotherapy has a poor overall response rate of approximately 33%, and a significant minority of kidney cancer patients do not respond to this treatment [4, 5]. One explanation for the ineffectiveness of therapy, as shown by several study findings, is that some kidney cancer patients have limited sensitivity to immunosuppressive drugs; another cause is drug resistance and not the tumour per se. Consequently, to increase the chances of survival among patients with kidney cancer, it is crucial to investigate the molecular mechanisms responsible for the onset and development (as referred to as O and D) of this disease.
A growing number of studies have shed light in recent years on the impact of pyroptosis on the O and D of malignancies. According to conventional wisdom, pyroptosis describes the cell necrosis caused by certain bacterial invasions carried out by the proteinase caspase-1, which is specific for the amino acid cysteinyl aspartate. After the identification of inflammatory compounds in 2002 and discovery of nonclassical inflammasomes in 2011, pyroptosis has been regarded as a mechanism related to cell death and inflammatory bodies [6,7,8]. It has never been thought that caspase cleavage causes pyroptosis. Nevertheless, some research shows that the ability to cause pyroptosis may be achieved by expression of the N-terminal domain of Gasdermin D (GSDMD) or another Gasdermin [9,10,11,12]. In addition, activation of caspase is not necessary for pyroptosis. For instance, granzyme A and B act as upstream molecules of GSDMB and GSDME to cleave them, which is irrelevant to the action of caspase [13, 14]. Currently, pyroptosis has been redefined as apoptosis mediated by GSDM proteins. Pore-forming effector proteins comprise the GSDM family. For pyroptosis to occur, these proteins must first be cleaved and then have the capacity to create holes in the cell membrane [9, 15, 16].
An RNA that is 200 bp in length or more is referred to as long noncoding RNA (lncRNA). While lncRNAs are abundant in the cytoplasm and nucleus, they do not encode proteins . The process of gene expression involves "noise", which has previously been thought to be lncRNAs . The reality that synthesis of lncRNAs is comparable to that of coding genes was shown by DERRIEN et al. . Splicing patterns, exon/intron organization, and histone modifications have all been discovered to be similar. Coding genes produce lncRNAs . Presently, research suggests that lncRNAs play a critical role in the O and D stages of kidney cancer. OSRC-2 kidney cancer cells were found to overexpress the oncogene BCL-W protein, and WANG  discovered that the more aggressive the cancer cells were, the more lncRNA-RP11-436H11.5 was overexpressed. After we treated these cells with the inhibitor ATB-737, the tumour cells became less aggressive, and the inhibitory effect became more distinct when the concentration of ATB-737 was increased. The content of the lncRNA GIHCG was clearly greater in 46 RCC patient tissue and plasma samples than in normal tissues (P < 0.01), according to an analytical investigation by HE et al. .In addition, there is increasing evidence for the pre-target value of lncrnas in aging and aging-related diseases .
LncRNAs regulate proteins related to pyroptosis signalling through downstream pathways. MALA T1 may considerably raise NLRP3 levels by upregulating production of ELA VL1 proteins, which would further cause pyroptosis to develop in renal cells and damage in diabetic nephropathy model mice. Patients with uric acid kidney disease have increased expression of ANRIL, an opposing lncRNA at the INK4 locus. One study found that ANRIL might activate BRCC3, leading to production of NLRP3 and IL-1β/18, which is linked to renal disorders [24, 25]. Yi et al. shed further light on the processes involved in the elevation of NLRP3 expression triggered by lncRNA Gm4419 and subsequent release of IL-1β. Their research showed that in mouse mesangial cells cultivated with high glucose, NF-κB (p50) acts as a promoter of the NLRP3 inflammasome . LncRNAs might be considered possible biomarkers of kidney diseases because they often have a favourable correlation with pyroptosis in kidney diseases.
The prognosis of KIRC was examined in this work using bioinformatics tools to assess the impact of pyroptosis-related lncRNAs. In addition, the application potential of a constructed prognostic model in tumour immunity and drug response was explored. Moreover, the accuracy of this model was verified in a separate cohort.
Materials and methods
We used API v3.0.0 to download KIRC patient TCGA database mutation information, corresponding clinical data, and transcriptome RNA-seq data available at https://portal.gdc.cancer.gov (Release date: October 29, 2021). We used expression profile information from 91 RCC samples in the ICGC database (https://dcc.icgc.org/) as our validation cohort.
pyroptosis-related prognostic LncRNA identification
According to Additional file 1: Table S1, 52 genes relevant to pyroptosis were gathered from earlier research and papers. Using a screening process based on gene annotation, 2876 lncRNAs were ultimately found in the TCGA cohort. The association study between 52 genes and lncRNAs associated with pyroptosis was conducted using the Pearson correlation coefficient. We identified lncRNAs associated with pyroptosis based on their absolute correlation coefficient (> 0.4) and P-value (< 0.001). As a consequence, 576 lncRNAs associated with pyroptosis were examined. After that, we performed univariate Cox regression analysis for OS, choosing lncRNAs with and without value in the process. P-values < 0.05 were used to determine whether lncRNAs are related to prognosis. Overall, 295 lncRNAs of prognostic relevance were found to be linked to pyroptosis.
Construction and evaluation of the prognostic model
The training set served as the basis for construction of the prognostic model, and the test set and the whole TCGA set were used to evaluate the model's prediction performance. The training set was analyzed using multivariate Cox regression, revealing six pyroptosis-related lncRNAs. The model, named "PLnRM", uses the following formula to calculate the risk score: risk score = coef (lncRNA1) × expr (lncRNA1) + coef (lncRNA2) × expr (lncRNA2) + …… + coef (lncRNAn) × expr (lncRNAn). In this equation, "coef" represents coefficients. Specifically, "coef (lncRNAn)" refers to the survival-related coefficient of lncRNAs, and "expr (lncRNAn)" refers to expression of lncRNAs.
To verify differentially expressed genes (DEGs), we carried out Gene Ontology (GO) analysis using R-pack cluster profilers. Since P < 0.05 shows that functional annotations are enriched, this threshold is defined by the p-value.
Analysis of the model for immunotherapeutic therapy
We evaluated and calculated abrupt change data using maftools of the R package, and we estimated the tumour mutation burden (TMB) by tumour-specific mutation genes. In addition, we used tumour immune dysfunction and exclusion (TIDE) computation to estimate immunotherapy response probability .
PCA with Kaplan‒Meier survival analysis
Six pyrogenic lncRNAs, as well as genes and lncRNAs associated with pyroptosis, were studied by examining expression patterns . In addition, differences in OS between the two models were assessed using Kaplan‒Meier survival analysis.
Investigation of possible PLnRM-targeting compounds for therapeutic use
Using the GDSC website, we were able to identify the semimaximum inhibitory concentration (IC50) of medications to locate therapeutic chemicals for KIRC patients and enhance clinical effectiveness. To calculate the IC50 for KIRC patients, we employed a pRophetic formula found in the R package.
Cox regression models, including single- and multiple-variable models, were utilized to show that individuals with KIRC do not have a prognostic pattern that is distinct from or unrelated to other clinical characteristics .
Establishing and verifying the predictive nomogram
We created a graph throughout the procedure to predict OS at 1, 3, and 5 years for KIRC patients to improve practicability and make the findings clear. We use the Hosmer‒Lemeshow test to demonstrate the consistency of the corrective curve.
Molecular and immune characteristics and ICI therapy in PLnRM risk groups
During signalling pathway analysis, we first analyzed differential expression of gene samples with R's limma package to assess samples with high (n = 285) and low (n = 245) PLnRM scores. Enrichment analysis was conducted to identify signalling pathways in which differentially expressed genes are involved using the gene set enrichment analysis (GSEA) approach. The KEGG and HALLMARK gene sets were used, along with the cluster profile package of R. A significance threshold of P < 0.05 and FDR < 0.25 was applied. For a few common gene sets, we next conducted single-sample GSEA (ssGSEA). Differences in patient survival were then explored using Kaplan‒Meier survival curves. The cBioPortal database provides information on gene changes for gene mutation research. Using the R package Maftools, we examined the quantity and quality of mutations in two PLnRM risk groups. Correlation analysis was performed on expression of PD1, PDL1, CTLA4, BTLA, and CD24 between the two PLnRM risk groups. We imported expression data from 530 KIRC samples into CIBERSORT (https://cibersort.stanford.edu/). After 1000 iterations, we estimated the relative percentage of 22 categories of immune cells to identify their immunological features. Additionally, the matching proportions of 22 different kinds of immune cells as well as clinicopathological variables were compared between the two PLnRM risk groups, and the findings are shown in a landscape map. ssGSEA was conducted on specific gene signatures to better elucidate the immunological and molecular functions between the two PLnRM risk groups, and the resulting scores were also compared [30,31,32,33].
Discovery of pyroptosis-associated LncRNAs in KIRC patients
Figure 1 illustrates the precise procedure for building the risk model and subsequent analyses. The TCGA database was used to derive matrix expression of 52 pyroptosis genes and 2,876 lncRNAs. We classified pyroptosis-related lncRNAs as those strongly linked to at least one of the 52 pyroptosis genes (|Pearson R|> 0.4 and P < 0.001). The pyroptosis-lncRNA coexpression network is depicted using a Sankey diagram (Fig. 2A). The results revealed 576 lncRNAs to be associated with pyroptosis. Figure 2B shows the relationship between pyroptosis genes and pyroptosis-related lncRNAs in the whole TCGA dataset.
KIRC patient risk model development and validation based on pyroptosis-related LncRNAs
We conducted univariate Cox regression analysis on 2876 pyroptosis-related lncRNAs in the KIRC dataset included in the TCGA database to identify prognostic lncRNAs. The TCGA database contains 295 pyroptosis-related lncRNAs with a strong correlation with the OS of KIRC patients (Fig. 3A). Multiple regression analysis often uses LASSO-penalized Cox analysis. It may concurrently perform selection and regularization of variables in addition to improving the statistical model's capacity for and accuracy at making predictions. This technique has been used widely to choose the best features in data in high dimensions with little connection and a strong projected value to prevent overfitting. As a result, this technique can efficiently identify the best predictive markers and provide a prognostic indicator to predict clinical outcomes. The value of log l that ranks first is shown on the dashed perpendicular line with the least bias in segment probability. Hence, Fig. 3B and C shows that 10 lncRNAs associated with pyroptosis were chosen for further multivariate analysis. Based on this, autocephalous prognostic proteins were identified using multivariate Cox ratio hazard regression analysis. To create a risk model for assessing the prognosis of KIRC patients, we utilized 6 pyroptosis-related lncRNAs (listed in Additional file 2: Table S2) that were independently associated with OS in the training set (as shown in Additional file 3: Table S3).
We conducted Kaplan‒Meier survival analysis on KIRC samples, dividing them into low-risk and high-risk groups using median prognostic risk scores. Figure 3D–F displays the survival status of patients in both groups across the full TCGA set, training set, and test set. The results showed a significant (P < 0.001) difference in outcomes between the high-risk and low-risk groups. Figure 4 A1 and 4A2 displays the patient survival status, survival time, and risk level distribution for both groups. The six lncRNAs associated with pyroptosis and their relative expression standards in each patient are shown in Fig. 4A3.
We used the same formula to calculate each patient's risk score in both the training and test sets. This allowed us to assess how effective our model is at predicting prognosis. Figures 4B, C depict expression of pyroptosis-related lncRNAs in the training set (Fig. 4B1–B3) and test set (Fig. 4 C1–C3). Additionally, they illustrate the distribution of risk grades, survival status patterns, and survival times (Fig. 4C1–C3).
To confirm the precision and applicability of the model, we verified expression profile data for 91 RCCs in the ICGC database. The results demonstrate the model's continued effectiveness in predicting survival time, particularly long-term survival (Fig. 5).
We compared differences in OS between low- and high-risk groups across the entire TCGA dataset, stratifying by common clinicopathologic characteristics. The low-risk group's OS remained superior to that of the high-risk group across subgroups categorized by age, sex, stage, and grade (see Additional file 4: Figure S1).
Validation of the PLnRM grouping capability through PCA
Comprehensive gene expression profiles, 52 pyroptosis genes, lncRNAs related to pyroptosis, and a risk model separated by lncRNA expression profiles all contributed to this model, which was then subjected to principal component analysis (PCA) to validate the distinction between low-risk and high-risk groups (Fig. 6A–D). As shown in Fig. 6A–C, there was absolutely no uniformity in the distribution of the high-risk and low-risk groups. Nonetheless, the findings from this model showed that the two groups’ distributions were significantly different (Fig. 6D). Based on these findings, it seems that there is a difference in the prognostic signature between the two groups.
PLnRM to assess the tumour immune microenvironment and response to immunotherapy
PLnRM was utilized to analyze the enrichment and activity levels of various immune cells, pathways, and activities in 530 KIRC patients. The findings showed that most immunological indicator expression levels between the two groups were significantly different (Fig. 7B). A GO enrichment study was performed to explore the putative molecular mechanisms related to PLnRM. Regarding BP, the model genes are related to a variety of biological functions relevant to the immune system (Fig. 7A). The relationship between PLnRM and the effectiveness of immunotherapy was, therefore, investigated. PLnRM can predict TIDE by indicating that the high-risk group will be more responsive to immunotherapies than the low-risk group, as hypothesized (Fig. 7H). In addition, maftools of the R package were used to investigate and assemble the data associated with mutations, and the consequences of mutation were used for stratification. Figure 7C and D displays the top 20 driver genes that show the greatest frequency change across all categories. TMB scores were then computed on the basis of the TGCA somatic mutation data. The findings indicate no significant difference in tumour mutation between the high-risk and low-risk groups, as shown in Fig. 7E. Kaplan‒Meier survival analysis of TMB was conducted using tumour tissues. The results (Fig. 7F) indicated that the low-mutation group had a higher chance of surviving than the high-mutation group. Furthermore, it was found that the prognosis for those in the high-mutation and high-risk groups was worse than that for those in the low-mutation and low-risk groups. Even when there were two groups with either a high or low mutation risk (Fig. 7G), individuals in the high-risk group still had a poorer prognosis than those in the low-risk group. This risk model seems to be reliable and stable based on these findings and our earlier findings.
Immunological and molecular traits of diverse PLnRM subgroups
The Wilcoxon test was used to examine the immune cell distribution in various PLnRM subgroups with the intention of examining the immune cell composition in various PLnRM subgroups. According to the results, the group with PLnRM-high risk had a higher prevalence of regulatory T cells (Tregs), activated memory CD4 T cells, CD8 T cells, memory B cells (MBCs), and macrophages (M0). Resting dendritic cells, macrophages (M2), resting NK cells, gamma delta (γδ) T cells, T follicular helper (Tfh) cells, resting memory CD4 T cells, plasma cells, naïve B cells, and activated dendritic cells were more prevalent in the PLnRM-low risk group, as shown in Fig. 8A, B. Figure 8C displays traits connected to the immunological landscape.
The gene sets enriched in various PLnRM subgroups were determined using GSEA. The P53 signalling pathway, ribosome, cytokine receptor interaction, systemic lupus erythematosus, and taste transduction pathways were enriched in the gene sets for samples with significant levels of PLnRM (Fig. 8D). However, according to Fig. 8E, the PLnRM-low samples showed gene sets that were enriched in pathways related to tumour metastasis (P < 0.05, FDR < 0.25). To further investigate how the two groups' responses to immunotherapy differed, variations in immune checkpoint expression and correlation were evaluated. As shown in Fig. 8F–J, CD47, B- and T-lymphocyte attenuator (BTLA), cytotoxic T lymphocyte-associated antigen (CTLA), PD ligand 1 (PD-L1), and programmed cell death (PD-1) were all higher in the high-risk group.
Clinical KIRC patient features and evaluation of the predictive risk model
This investigation utilized univariate and multivariate Cox regression analyses to determine whether the risk model incorporates independent prognostic factors for KIRC. Multivariate Cox regression analysis showed a hazard ratio (HR) of 1.067 and a 95% confidence interval (CI) of 1.050–1.085, with a significance level of P < 0.001 (Fig. 9A). In univariate Cox regression analysis, the odds ratios (HRs) were found to be 1.034, with a corresponding CI of 1.014–1.054, and significant results at P < 0.001 (Fig. 9B). The findings support the notion that the risk model can predict prognosis regardless of other clinical variables. We evaluated the distinctiveness and sensitivity of risk scores in predicting KIRC patient outcomes using the concordance index and area under the receiver operating characteristic (ROC) curve (AUC), as shown in Fig. 9D, E. Longer follow-up periods resulted in a greater concordance index of risk ratings, which ultimately trumped other clinical criteria and increased risk severity. Thus, the prognosis of KIRC patients may be accurately predicted using this model's risk score (Fig. 9C). Moreover, the AUC representing the risk level rose to a point where it exceeded that of the vast majority of other clinicopathological characteristics. Our results tended to support the idea that PLnRM is a viable option for incorporation into the predictive risk model that is currently being used for KIRC patients (Fig. 9C).
Development and analysis of the predictive nomogram
The 1-, 2-, and 3-year OS of KIRC patients can be predicted using a nomogram that takes into account both risk levels and clinical risk factors. The nomogram revealed that the prediction model's risk level has high predictive power when compared with clinical parameters (Fig. 10A). The 1-, 2-, and 3-year OS observation and prediction rates were all shown to be in good agreement by means of appropriate diagrams (Fig. 10B).
Identification of novel PLnRM candidate compounds
We used the prediction approach to evaluate treatment response for KIRC patients by analyzing the IC50 of each sample in Genomics of Drug Sensitivity in Cancer (GDSC). Our goal was to identify potential medications for PLnRM. There were 110 discovered compounds in all, and the predicted IC50 values between the two groups varied significantly. Partially sensitive chemicals are shown in Additional file 5: Figure S2.
Investigation of possible medicines for patients with KIRC
Several possible therapeutic medications for KIRC were evaluated using the connection map (CMAP) database to elucidate a new therapy plan for the disease. The top 10 closely related molecules of medications were chosen to depict the three-dimensional structure (Additional file 6: Figure S3).
The third most common tumour in the urinary system by incidence, RCC, is a kind of cancer that develops from the epithelium that lines the renal tubules, and its prevalence is growing . Although surgical resection is the most effective treatment for RCC, many patients receive their diagnosis in the middle or late stage. Furthermore, these tumours are unresponsive to chemotherapies, immunotherapies, and radiotherapies. Targeted treatments may also lead to temporary drug resistance. RCC patients often have poor prognosis as a result [35, 36]. RCC is influenced by a number of different variables that may impact its onset and progression, as well as by a number of different genes. Another significant contributor to RCC is abnormal alterations in the network that controls gene expression . Control of gene expression is influenced by various factors, including the level of genes and their regulation at different stages, such as transcription, translation, and protein degradation. Understanding lncRNA roles and identification allows for new perspectives on how gene expression is regulated.
Pyroptosis is a type of cell death that involves rapid rupture of cell membranes, resulting in release of cellular contents and proinflammatory substances such as interleukin (IL), IL-1β, and IL-18. This process is characterized by cell swelling and the formation of large bubbles in the cytoplasm [38, 39]. The consequences of pyroptosis on several inflammation-related disorders, such as heart disease, sepsis, diabetes, nephropathy, and atherosclerosis, are mostly due to release of these distinctive inflammatory factors [40,41,42].
In recent years, there have been many breakthroughs in understanding the mechanisms, molecules, and pathways related to pyroptosis. Pyroptosis has been shown to be related to the incidence, growth, prognosis, and treatment of a number of cancers [43, 44]. In addition, pyroptosis has been verified to participate in chimeric antigen receptor T-cell immunotherapy (CAR-T therapy), cytokine release syndrome (CRS), and chemotherapy [45,46,47]. In macrophages, certain factors induce the release of DNA to trigger the cGAS-STING induced IFN response, thereby regulating pyroptosis and inflammation, which provides further evidence for the regulation of pyroptosis [48,49,50]. Exploring the mechanism of pyroptosis and its correlation with tumours from a comprehensive perspective is conducive to broadening our understanding of tumours, which might indicate a new direction for cancer therapy.
The process of pyroptosis may be directly regulated by lncRNAs, as previously indicated, in addition to having an indirect effect. LncRNAs have been shown to directly control pyroptosis, according to recent studies. For example, it was found that the lncRNA Neat1 stabilizes mature caspase-1 tetramers (p20: p10)2 and (p33: p10)2 in mouse bone marrow-derived macrophages (BMEMs) treated with flagellin and poly(I:C) after induction with LPS. This study found that lncRNA directly binds to pro-caspase-1, promoting assembly of NLRP3 and AIM2 inflammasomes. Additionally, it was discovered that lncRNAs play a regulatory role in the pyroptosis signalling pathway affecting the inflammasome. Hence, a KIRC predictive model produced using lncRNAs associated with pyroptosis would be valuable.
In this work, we built an independent prognostic model based on lncRNAs associated with pyroptosis in KIRC, considering the functions of both, and whether there are any medications that may be useful for treating KIRC is thus being explored. It is possible to investigate the prognostic role of pyroptosis-related lncRNAs using the 576 pyroptosis-related lncRNAs in the TCGA database. The TCGA database findings show that 10 pyroptosis-related lncRNAs have predictive significance and that a model including these lncRNAs may be useful to predict the overall survival rate in KIRC patients after utilizing these 6 lncRNAs. The median prognostic risk score was used to classify KIRC patients as having high or low risk. The high-risk group had poorer prognosis, as indicated by outcomes. Multivariate Cox regression analysis revealed that the lncRNA model associated with pyroptosis is an autologous risk factor for OS. According to the findings of ROC analysis, this model is superior to the majority of common clinical characteristics in predicting the overall survival of individuals with KIRC. In addition, a nomogram was shown to illustrate the degree to which the observed and predicted OS rates correspond over the first, third, and fifth years of follow-up. The prediction rates for the first, third, and fifth years all agreed quite well. The risk model built using 10 pyroptosis-associated lncRNAs that are independently related to KIRC OS has greater accuracy. This prediction model can help to identify new biomarkers for future research.
The TIDE method is also used to predict the probability of an immunotherapeutic response, revealing that the high-risk group had a higher immune response rate than the low-risk group. This suggests that immune-related medications may be more effective in predicting outcomes for the high-risk group. This discovery may be applied to guidelines for immune-related medications.
Furthermore, the model was used to analyze the immunological and biochemical traits of various subgroups. According to the findings, there are some variations between the high-risk and low-risk groups in terms of immune cell enrichment and infiltration. While analyzing expression and connection of popular immune-related genes such as PD1 and PD-L1, significant differences in their expression between the high-risk and low-risk groups were found. Moreover, there was an inverse relationship between the risk score and expression of these genes.
This model was validated using datasets found in the TCGA database. The ICGC RCC dataset is integrated as an external cohort so that the accuracy of this model and its application may be evaluated. The findings of the survival research indicate that there is a significant gap between the high-risk and low-risk groups. As a result, our model is capable of accurately predicting the chances of survival for KIRC patients. With the help of AI and deep learning models, we can well utilize radio-genomics combined with the prediction model of KIRC reported in this paper .
The prognosis of KIRC patients is largely determined by the pathological stage and grade, yet tumours with the same clinical stage and grade may not always have the same prognosis . Exploring more detailed and focused prediction signs or biomarkers is thus very important. The pyroptosis-related lncRNA model that was developed is intended to provide a novel strategy for predicting the prognosis of KIRC patients. These results offer a new approach to studying lncRNAs associated with pyroptosis and modification processes. Several approaches are used in this work to validate the new model, making it possible to choose and implement the best model. Due to the lack of external data verification, it may be presumed that this prediction model is appropriate.
However, this research still has several shortcomings. For example, some questions about the biological mechanism of lncRNAs associated with pyroptosis remain. Investigating lncRNA function and how lncRNAs interact with genes involved in pyroptosis is thus crucial. As a whole, the conclusions provide fresh perspectives for predicting KIRC patient survival and prognosis, which may help to shed light on the mechanism behind pyroptosis-related lncRNAs. As a result of the development of this immunotherapy-sensitive model, several preliminary medication candidates were also found that may be useful, impacting the therapeutic approach used for KIRC patients.
Availability of data and materials
The datasets presented in this study can be found in online repositories. The names of the repositories/repositories and accession number(s) can be found in the article/Supplementary Material.
Renal cell carcinoma
Kidney renal clear cell carcinoma
Long noncoding RNAs
Tumour mutation burden
Tumour immune dysfunction and exclusion
Gene set enrichment analysis
Memory B cells
T follicular helper
Principal component analysis
B- and T-lymphocyte attenuator
Genomics of Drug Sensitivity in Cancer
Cytokine release syndrome
Bone marrow-derived macrophages
Ljungberg B, Albiges L, Abu-Ghanem Y, Bensalah K, Dabestani S, Fernández-Pello S, et al. European association of urology guidelines on renal cell carcinoma: the 2019 update. Eur Urol. 2019;75(5):799–810. https://doi.org/10.1016/j.eururo.2019.02.011.
Finelli A, Ismaila N, Bro B, Durack J, Eggener S, Evans A, et al. Management of small renal masses: american society of clinical oncology clinical practice guideline. J Clin Oncol. 2017;35(6):668. https://doi.org/10.1200/JCO.2016.69.9645.
Ge W, Song S, Qi X, Chen F, Che X, Sun Y, et al. Review and prospect of immune checkpoint blockade therapy represented by PD-1/PD-L1 in the treatment of clear cell renal cell carcinoma. Oncol Res. 2023;31(3):255–70. https://doi.org/10.32604/or.2023.027942.
Motzer RJ, Escudier B, McDermott DF, George S, Hammers HJ, Srinivas S, et al. Nivolumab versus Everolimus in advanced renal-cell carcinoma. N Engl J Med. 2015;373(19):1803–13. https://doi.org/10.1056/NEJMoa1510665.
Zeng SY, Liu YF, Liu JH, Zeng ZL, Xie H, Liu JH. Potential effects of Akkermansia Muciniphila in aging and aging-related diseases: current evidence and perspectives. Aging Dis. 2023. https://doi.org/10.14336/AD.2023.0325.
Martinon F, Burns K, Tschopp J. The inflammasome: a molecular platform triggering activation of inflammatory caspases and processing of proIL-beta. Mol Cell. 2002;10(2):417–26. https://doi.org/10.1016/s1097-2765(02)00599-3.
Hersh D, Monack DM, Smith MR, Ghori N, Falkow S, Zychlinsky A. The Salmonella invasin SipB induces macrophage apoptosis by binding to caspase-1. Proc Natl Acad Sci U S A. 1999;96(5):2396–401. https://doi.org/10.1073/pnas.96.5.2396.
Brennan MA, Cookson BT. Salmonella induces macrophage death by caspase-1-dependent necrosis. Mol Microbiol. 2000;38(1):31–40. https://doi.org/10.1046/j.1365-2958.2000.02103.x.
Cookson BT, Brennan MA. Pro-inflammatory programmed cell death. Trends Microbiol. 2001;9(3):113–4. https://doi.org/10.1016/s0966-842x(00)01936-3.
Tsuchiya K. Switching from apoptosis to pyroptosis: gasdermin-elicited inflammation and antitumor immunity. Int J Mol Sci. 2021;22(1):426. https://doi.org/10.3390/ijms22010426.
Zhang Z, Lieberman J. Lighting a fire on the reef. Sci Immunol. 2020;5(54):eabf0905. https://doi.org/10.1126/sciimmunol.abf0905.
Xu YJ, Zheng L, Hu YW, Wang Q. Pyroptosis and its relationship to atherosclerosis. Clin Chim Acta. 2018;476:28–37. https://doi.org/10.1016/j.cca.2017.11.005.
Pirzada RH, Javaid N, Choi S. The roles of the NLRP3 inflammasome in neurodegenerative and metabolic diseases and in relevant advanced therapeutic interventions. Genes (Basel). 2020;11(2):131. https://doi.org/10.3390/genes11020131.
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. https://doi.org/10.1038/s41586-020-2071-9.
Bergsbaken T, Fink SL, Cookson BT. Pyroptosis: host cell death and inflammation. Nat Rev Microbiol. 2009;7(2):99–109. https://doi.org/10.1038/nrmicro2070.
Slack FJ, Chinnaiyan AM. The role of non-coding RNAs in Oncology. Cell. 2019;179(5):1033–55. https://doi.org/10.1016/j.cell.2019.10.017.
Gandhi M, Caudron-Herger M, Diederichs S. RNA motifs and combinatorial prediction of interactions, stability and localization of noncoding RNAs. Nat Struct Mol Biol. 2018;25(12):1070–6. https://doi.org/10.1038/s41594-018-0155-0.
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22(9):1775–89. https://doi.org/10.1101/gr.132159.111.
Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, et al. The transcriptional landscape of the mammalian genome. Science. 2005;309(5740):1559–63. https://doi.org/10.1126/science.1112014.
Wang K, Jin W, Song Y, Fei X. LncRNA RP11-436H11.5, competitive endogenous RNA, upregulates BCL-W expression by sponging miR-335–5p and promotes proliferation and invasion in renal cell carcinoma. Mol Cancer. 2017;16(1):166. https://doi.org/10.1186/s12943-017-0735-3.
Liu Y, Chen X, Che Y, Li H, Zhang Z, Peng W, Yang J. LncRNAs as the regulators of brain function and therapeutic targets for Alzheimer’s disease. Aging Dis. 2022;13(3):837–51. https://doi.org/10.14336/AD.2021.1119.
Li X, Zeng L, Cao C, Lu C, Lian W, Han J, et al. Long noncoding RNA MALAT1 regulates renal tubular epithelial pyroptosis by modulated miR-23c targeting of ELAVL1 in diabetic nephropathy. Exp Cell Res. 2017;350(2):327–35. https://doi.org/10.1016/j.yexcr.2016.12.006.
Hu J, Wu H, Wang D, Yang Z, Dong J. LncRNA ANRIL promotes NLRP3 inflammasome activation in uric acid nephropathy through miR-122-5p/BRCC3 axis. Biochimie. 2019;157:102–10. https://doi.org/10.1016/j.biochi.2018.10.011.
Py BF, Kim MS, Vakifahmetoglu-Norberg H, Yuan J. Deubiquitination of NLRP3 by BRCC3 critically regulates inflammasome activity. Mol Cell. 2013;49(2):331–8. https://doi.org/10.1016/j.molcel.2012.11.009.
Yi H, Peng R, Zhang LY, Sun Y, Peng HM, Liu HD, et al. LincRNA-Gm4419 knockdown ameliorates NF-κB/NLRP3 inflammasome-mediated inflammation in diabetic nephropathy. Cell Death Dis. 2017;8(2):e2583. https://doi.org/10.1038/cddis.2016.451.
Xu F, Zhan X, Zheng X, Xu H, Li Y, Huang X, et al. A signature of immune-related gene pairs predicts oncologic outcomes and response to immunotherapy in lung adenocarcinoma. Genomics. 2020;112(6):4675–83. https://doi.org/10.1016/j.ygeno.2020.08.014.
Li X, Li Y, Yu X, Jin F. Identification and validation of stemness-related lncRNA prognostic signature for breast cancer. J Transl Med. 2020;18(1):331. https://doi.org/10.1186/s12967-020-02497-4.
Xu F, Huang X, Li Y, Chen Y, Lin L. m6A-related lncRNAs are potential biomarkers for predicting prognoses and immune responses in patients with LUAD. Mol Ther Nucleic Acids. 2021;24:780–91. https://doi.org/10.1016/j.omtn.2021.04.003.
He Y, Jiang Z, Chen C, Wang X. Classification of triple-negative breast cancers based on Immunogenomic profiling. J Exp Clin Cancer Res. 2018;37(1):327. https://doi.org/10.1186/s13046-018-1002-1.
Zeng D, Ye Z, Wu J, Zhou R, Fan X, Wang G, et al. Macrophage correlates with immunophenotype and predicts anti-PD-L1 response of urothelial cancer. Theranostics. 2020;10(15):7002–14. https://doi.org/10.7150/thno.46176.
Hwang S, Kwon AY, Jeong JY, Kim S, Kang H, Park J, et al. Immune gene signatures for predicting durable clinical benefit of anti-PD-1 immunotherapy in patients with non-small cell lung cancer. Sci Rep. 2020;10(1):643. https://doi.org/10.1038/s41598-019-57218-9.
Tirosh I, Izar B, Prakadan SM, Wadsworth MH, Treacy D, Trombetta JJ, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science (New York, NY). 2016;352(6282):189–96. https://doi.org/10.1126/science.aad0501.
Pang C, Guan Y, Li H, Chen W, Zhu G. Urologic cancer in China. Jpn J Clin Oncol. 2016;46(6):497–501. https://doi.org/10.1093/jjco/hyw034.
Yamamoto K, Ioroi T, Shinomiya K, Yoshida A, Harada K, Fujisawa M, et al. STAT3 polymorphism associates With mTOR inhibitor-induced interstitial lung disease in patients with renal cell carcinoma. Oncol Res. 2022;29(1):11–23. https://doi.org/10.3727/096504022X16418911579334.
Funakoshi T, Lee CH, Hsieh JJ. A systematic review of predictive and prognostic biomarkers for VEGF-targeted therapy in renal cell carcinoma. Cancer Treat Rev. 2014;40(4):533–47. https://doi.org/10.1016/j.ctrv.2013.11.008.
Jonasch E, Walker CL, Rathmell WK. Clear cell renal cell carcinoma ontogeny and mechanisms of lethality. Nat Rev Nephrol. 2021;17(4):245–61. https://doi.org/10.1038/s41581-020-00359-2.
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:111–6. https://doi.org/10.1038/nature18590.
Broz P, Pelegrin P, Shao F. The gasdermins, a protein family executing cell death and inflammation. Nat Rev Immunol. 2020;20:143–57. https://doi.org/10.1038/s41577-019-0228-2.
Yin Y, Li X, Sha X, Xi H, Li YF, Shao Y, et al. Early hyperlipidemia promotes endothelial activation via a caspase-1-sirtuin 1 pathway. Arterioscler Thromb Vasc Biol. 2015;35:804–16. https://doi.org/10.1161/ATVBAHA.115.305282.
Mamun AA, Wu Y, Nasrin F, Akter A, Taniya MA, Munir F, et al. Role of pyroptosis in diabetes and its therapeutic implications. J Inflamm Res. 2021;14:2187–206. https://doi.org/10.2147/JIR.S291453.
Jia C, Chen H, Zhang J, Zhou K, Zhuge Y, Niu C, et al. Role of pyroptosis in cardiovascular diseases. Int Immunopharmacol. 2019;67:311–8. https://doi.org/10.1016/j.intimp.2018.12.028.
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. https://doi.org/10.1016/j.intimp.2019.105713.
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:415–20. https://doi.org/10.1038/s41586-020-2071-9.
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:99–103. https://doi.org/10.1038/nature22393.
Kumagai S, Togashi Y, Kamada T, Sugiyama E, Nishinakamura H, Takeuchi Y, et al. The PD-1 expression balance between effector and regulatory T cells predicts the clinical efficacy of PD-1 blockade therapies. Nat Immunol. 2020;21:1346–58. https://doi.org/10.1038/s41590-020-0769-3.
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:1264–75. https://doi.org/10.1038/s41556-020-0575-z.
Zhang P, Cao L, Zhou R, Yang X, Wu M. The lncRNA Neat1 promotes activation of inflammasomes in macrophages. Nat Commun. 2019;10(1):1495. https://doi.org/10.1038/s41467-019-09482-6.
Xu Y, Chen C, Liao Z, Xu P. cGAS-STING signaling in cell death: mechanisms of action and implications in pathologies. Eur J Immunol. 2023. https://doi.org/10.1002/eji.202350386.
Zheng W, Feng D, Xiong X, Liao X, Wang S, Xu H, Le W, Wei Q, Yang L. The role of cGAS-STING in age-related diseases from mechanisms to therapies. Aging Dis. 2023. https://doi.org/10.14336/AD.2023.0117.
Huang A, Zhou W. Mn-based cGAS-STING activation for tumor therapy. Chin J Cancer Res. 2023;35(1):19–43. https://doi.org/10.21147/j.issn.1000-9604.2023.01.04.
Ferro M, Musi G, Marchioni M, Maggi M, Veccia A, Del Giudice F, et al. Radiogenomics in renal cancer management-current evidence and future prospects. Int J Mol Sci. 2023;24(5):4615. https://doi.org/10.3390/ijms24054615.
Feng S. Gasdermins: making pores for pyroptosis. Nat Rev Immunol. 2021. https://doi.org/10.1038/s41577-021-00602-2.
The related research work in this paper was funded by the Natural Science Foundation of Anhui Province (No. 2108085MH297), Research Foundation of Anhui Medical University: Effect and mechanism of exosome transport circ55033 on the enhancement of bladder cancer stem cell dryness and EMT by cigarette exposure (2022xkj025), Study on the mechanism of LncRNAAGAP2-AS1 regulation of TAMs polarization and exosomal miR-296 secretion on RCC Sunitinib resistance (2021xkj040) and Anhui Provincial Education Department Scientific Research Project Funding (YJS20210311).
Ethics approval and consent to participate
Consent for publication
All patients signed the relevant paper informed consent form.
All participants were in the same subject group. They are in a cooperative relationship with each other. There is no conflict of interest between them.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Gene names of the 52 pyroptosis-related genes.
Names of the six lncRNAs in the model that are connected to pyroptosis.
In the risk model, the clinical characteristics are shared by the complete set, the training set, and the test set.
Kaplan‒Meier curves show differences in overall survival between the two groups in the entire TCGA set, stratified by age, sex, tumour grade, and stage.
Compounds with partial sensitivity.
Three-dimensional structure of the top ten related drugs.
About this article
Cite this article
Liu, C., Dai, S., Geng, H. et al. Development and validation of a kidney renal clear cell carcinoma prognostic model relying on pyroptosis-related LncRNAs-A multidimensional comprehensive bioinformatics exploration. Eur J Med Res 28, 341 (2023). https://doi.org/10.1186/s40001-023-01277-2