3.1 IL20RB is a potential key prognosis-related gene in pRCC. We used the transcriptome RNA-seq data from the TCGA database to obtain the OS differential genes, DSS differential genes, and PFI differential genes of pRCC, respectively, and Venn diagrams were used to visualize their intersection. As shown in Fig.1a, there are 47 genes that are common intersections of the OS differential gene, the DSS differential gene, and the PFI differential gene. Lasso regression can narrow down the gene screening scope and improve the accuracy and interpretability of the analysis results through some mechanisms such as automatic feature selection and prevention of overfitting. Therefore, by lasso regression analysis, we further screened OS differential genes, DSS differential genes and PFI differential genes to obtain four key prognostic genes of pRCC, which were IGF2BP3, SLC7A11, FBXO43 and IL20RB, and visualized them in a Venn diagram (Fig.1b). As shown in Fig.1c and Fig.1d, four genes were significantly more expressed in tumor tissues than in normal tissues in both non-paired and paired samples of pRCC. To determine whether these four genes have diagnostic value for pRCC, we performed ROC correlation analysis in combination with TCGA data and visualized the results. As shown in Fig.1e, the AUC (Area Under Curve) of IGF2BP3, SLC7A11, FBXO43, and IL20RB were all greater than 0.5, but the AUC of IL20RB was significantly larger, which indicates that IL20RB has a higher diagnostic value in pRCC and is worth focusing on.
3.2 Specific Expression of IL20RB and IL20 Subfamily Genes in pRCC. In order to have a visual panorama of IL20RB expression in pRCC, we used the GEPIA2 website to obtain the average expression of IL20RB in tumor tissues as well as normal tissues of pRCC. Also considering that IL20RB is a key receptor subunit of IL20, we also visualized the average expression of IL20 subfamily genes in tumor tissues as well as normal tissues of pRCC. As shown in Fig.1f, IL20RB was highly expressed in pRCC tumor tissues relative to normal tissues. In addition, except for IL24, other IL20 subfamily genes were not expressed in tumor tissues or were more prominently expressed in normal tissues. It has been demonstrated that IL24 is an immunomodulatory cytokine and is mainly characterized by the induction of tumor cell apoptosis to inhibit tumor growth in cancer[29, 30]. The above suggests that IL20RB may not play a biological role in pRCC through direct binding with IL20 subfamily genes. To further understand the correlation between IL20RB expression in pRCC and clinical features, we compared IL20RB expression in pRCC patients with different T stages as well as clinical stages based on TCGA data. As illustrated in Figure 2, the expression of IL20RB was different and the difference was statistically significant at different pathologic stages as well as clinical stages. The expression of IL20RB gradually increased with the elevation of pathologic T stage, N stage, and M stage. Similarly, with the increase of clinical T stage, N stage and M stage, the expression of IL20RB also gradually increased. This suggests that IL20RB expression is positively correlated with clinical features.
3.3 Immunohistochemistry of IL20RB in pRCC. We collected a total of 42 postoperative tissue samples from patients with pRCC at Zhongda Hospital, Southeast University. All participants provided informed consent, and the study was approved by the Ethics Committee of Zhongda Hospital, Southeast University (2022ZDSYLL392-P01). Fig.3A and Fig.3B demonstrates the carcinoma tissue of pRCC, while Fig.3C and Fig.3D show paracarcinoma renal tissue of pRCC with IL20RB staining in brown. We observed the staining of carcinoma and paracancerous tissues in each pRCC specimen separately and compared the staining intensity. Staining intensity of carcinoma tissues was greater than that of paracancerous tissues implying upregulation of IL20RB expression, while staining intensity of carcinoma tissues was lower than that of paracancerous tissues implying decreased IL20RB expression. We found that IL20RB was mainly expressed in the cytoplasm of cancer cells in papillary renal carcinoma tissues, and some of it was also expressed in the cell membrane, and most of the expression was positive and strongly positive. Meanwhile, IL20RB was also widely expressed in the cytoplasm of renal tubular epithelial cells adjacent to cancer cells. Among the 42 patients with pRCC in our center, 33 patients presented high expression of IL20RB in immunohistochemically stained sections (Fig.3E and Fig.3F) and 9 patients presented low expression of IL20RB in immunohistochemically stained sections (Fig.3G and Fig.3H), which shows that IL20RB is overall upregulated in expression in tumor cells of pRCC. Meanwhile, in the IL20RB high expression group, we compared IL20RB staining results in different T stages (Fig.4A-4D). We found that IL20RB expression gradually increased with increasing T stage, which could also be reflected by the positive cell index of IL20RB (Fig.4E).
3.4 Survival Analysis of IL20RB. We obtained the Kaplan-Meier curve for IL20RB by analyzing the data of pRCC patients in TCGA. As shown in Figure 5a, pRCC patients with high expression of IL20RB had a significantly worse prognosis than pRCC patients with low expression of IL20RB. To further validate this result, we combined the data from our center and divided the patients into IL20RB high-expression and low-expression groups based on the results of immunohistochemical staining, thus obtaining the Kaplan-Meier curve of IL20RB in pRCC patients in our center (Figure 5b). We suggest that IL20RB may be an oncogene and its high expression may be associated with a worse prognosis in pRCC.
3.5 Immune Infiltration Analysis. To analyze the potential correlation between the expression of IL20RB and the immune infiltration status in the tumor microenvironment of pRCC, we utilized the Quantiscore method to explore the relationship between IL20RB and a range of major immune cells. A significant negative correlation was found between IL20RB expression and the infiltration of CD4+ T cells (Fig.6g), which indicated that IL20RB may promote the malignant progression of pRCC through downregulation of adaptive anticancer immunity. In addition, although the CD8+ T cells also exhibit positive correlation with the expression of IL20RB, the overall expressive extent is relatively low (Fig.6h). More importantly, as the major immune suppressive cells[31], regulatory T cells were found to be significantly positively-associated with the up-regulation of IL20RB (Fig.6i), which demonstrated that IL20RB may also induce immune tolerance in pRCC through induction of polarization or recruitment of regulatory T cells. In addition to T cells, IL20RB was also found to be correlated with a range of myeloid cells, such as neutrophils (Fig.6e), M1 macrophages (Fig.6b), and monocytes (Fig.6d), highlighting its involvement in shaping the myeloid component of the pRCC immune microenvironment.
3.6 Functional Analysis of IL20RB in pRCC. To validate the potential function of IL20RB, we performed KEGG and GO enrichment analyses based on pRCC expression profiling data from TCGA. Figure 7a illustrates that, in KEGG enrichment analysis, differentially expressed genes (DEGs) associated with IL20RB upregulation were mainly enriched in the PI3K−Akt signaling pathway and Cytokine−cytokine receptor interaction. PI3K-Akt signaling has been found to be involved in regulating tumor cell metastasis and invasion in a variety of tumors[32-34]. Figure 7b depicts the results of GO enrichment analysis, indicating that DEGs associated with IL20RB upregulation were primarily enriched in the T cell activation, extracellular structure organization, extracellular matrix organization and positive regulation of cell adhesion. In contrast, DEGs associated with IL20RB down-regulation were predominantly enriched in the Carbon metabolism and small molecule catabolic process (Fig.7c and Fig.7d). The heatmap of IL20RB-related differential genes in pRCC is shown in Supplementary Figure 1. Significant differential genes associated with IL20RB up-regulation include TPM2, CIR, BST2, CIS, ANGPTL4 and so on. Significant differential genes associated with IL20RB down-regulation mainly include OGDHL, MIOX, AKR7A3 and SLC13A1. Moreover, we utilized the String method to explore the protein-protein interactions and potential functional molecular partners centered on IL20RB. As is demonstrated in the PPI map (Figure.8), significant interactions are found between IL20RB and a range of crucial cytokines, such as IFN-γ, IL-6, IL-10, and TNF, etc., which indicated that IL20RB may be involved in the regulation of the inflammatory microenvironment of pRCC. Moreover, multiple key molecular adaptors and signaling pathways in oncologic processes are also found to be in interaction with IL20RB, such as JAK1 and TYK2, which indicated that IL20RB may also mediates the oncogenic pathways of pRCC through classical signaling pathways.
3.7 Drug Sensitivity Analysis of IL20RB in pRCC. By using the GDSC database, we performed drug sensitivity analysis of pRCC patients in TCGA. As shown in Figure 9, individuals with high IL20RB expression were sensitized to common targeted drugs including axitinib, pazotinib, sunitinib, tivozanib, and the chemotherapeutic drug cisplatin. Notably, tivozanib had the lowest half-maximal inhibitory concentration (Fig.5e), suggesting its potential as a highly effective therapeutic option for pRCC patients with elevated IL20RB expression. Although there are no prospective studies evaluating the efficacy of tivozanib in pRCC, a phase III clinical trial by Rini, Brian I et al. also found that tivozanib as third-line or fourth-line therapy improved progression-free survival and was better tolerated compared with sorafenib in patients with metastatic renal cell carcinoma[35], which may provide some inspiration for future clinical trials of tivozanib in pRCC. Conversely, our analysis revealed that pRCC patients with high expression of IL20RB display low sensitivity to erlotinib (Fig.9b), indicating that erlotinib may have limited therapeutic effect in this subgroup of pRCC patients in clinical practice.