Skip to main content

Identification of 4-genes model in papillary renal cell tumor microenvironment based on comprehensive analysis



The tumor microenvironment acts a pivotal part in the occurrence and development of tumor. However, there are few studies on the microenvironment of papillary renal cell carcinoma (PRCC). Our study aims to explore prognostic genes related to tumor microenvironment in PRCC.


PRCC expression profiles and clinical data were extracted from The Cancer Gene Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Immune/stromal scores were performed utilizing the ESTIMATE algorithm. Three hundred fifty-seven samples were split into two groups on the basis of median immune/stromal score, and comparison of gene expression was conducted. Intersect genes were obtained by Venn diagrams. Hub genes were selected through protein-protein interaction (PPI) network construction, and relevant functional analysis was conducted by DAVID. We used Kaplan–Meier analysis to identify the correlations between genes and overall survival (OS) and progression-free survival (PFS). Univariate and multivariate cox regression analysis were employed to construct survival model. Cibersort was used to predict the immune cell composition of high and low risk group. Combined nomograms were built to predict PRCC prognosis. Immune properties of PRCC were validated by The Cancer Immunome Atlas (TCIA).


We found immune/stromal score was correlated with T pathological stages and PRCC subtypes. Nine hundred eighty-nine differentially expressed genes (DEGs) and 1169 DEGs were identified respectively on the basis of immune and stromal score. Venn diagrams indicated that 763 co-upregulated genes and 4 co-downregulated genes were identified. Kaplan-Meier analysis revealed that 120 genes were involved in tumor prognosis. Then PPI network analysis identified 22 hub genes, and four of which were significantly related to OS in patients with PRCC confirmed by cox regression analysis. Finally, we constructed a prognostic nomogram which combined with influence factors.


Four tumor microenvironment-related genes (CD79A, CXCL13, IL6 and CCL19) were identified as biomarkers for PRCC prognosis.

Peer Review reports


The incidence of renal cell carcinoma (RCC) is approximately 431,288 in 2020, accounting for 2.2% of common malignant neoplasm on a global scale [1]. Clear cell, papillary, and chromophobe renal cell carcinoma are three main histologic subtypes of RCC. Study reported papillary renal cell carcinoma (PRCC) is the second most common type [2]. Clinicopathological features are the main diagnostic criteria for PRCC, but their predictive outcome is not accurate because of the lack of consistent standards [3]. Surgery is the primary choice of treatment for early stage PRCC, and comprehensive therapeutic approaches which include surgery, and immunotherapy are applied for advanced stage patients [4]. However, the prognosis for PRCC is various on account of tumor metastasis and complications. The World Health Organization (WHO) classification (2016), according to different clinicopathological and immunohistochemical features, split PRCC into two subtypes: type 1 and type 2 [5]. In general, the prognosis of type 2 PRCC was poorer than type 1 [6].

Tumor microenvironment (TME) is an intricate system comprised of tumor cells, surrounding immune and inflammatory cells, tumor-related fibroblasts and nearby stromal tissues, and varieties of cytokines and chemokines [7]. Relevant studies have indicated macrophages play a crucial part in tumor initiation and invasion of adjacent tissues [8]. Tumor-related fibroblasts, a kind of stromal cell, affect cancer progression partly by interacting with tumor cells and immune regulation [9]. Studies have demonstrated the highly enriched CD8+ T cells in PRCC are significantly connected with tumor development, progression, and mortality [10]. Yoshihara et al. put forward the ESTIMATE algorithm that utilizes gene expression patterns to calculate the immune/stromal scores in different tumor tissues [11]. Accumulating evidence showed ESTIMATE helped to clarify the importance of TME in numerous cancers [12,13,14].

In this study, we extracted PRCC datasets from TCGA and GEO and obtained relevant immune and stromal score calculated by ESTIMATE. After a series of bioinformatics analyses, several tumor microenvironment-related genes in connection with the prognosis of PRCC patients were selected and validated, and were used to establish nomograms as potential tools to evaluate the outcomes of PRCC.

Materials and methods

Data preparation

Level_3 gene expression patterns for PRCC was downloaded from TCGA via UCSC Xena ( Another dataset (GSE2748) was obtained from the GEO database. And corresponding clinical information including overall survival (OS) time and status was extracted from the ONCOMINE database ( Immune and stromal scores were calculated by ESTIMATE packages of R.

Screening of differentially expressed genes (DEGs)

PRCC samples were split into high and low groups on the basis of median immune/stromal score, limma package of R was utilized to compare DEGs, |Log2FC| > 1 and FDR < 0.05 were deemed significant. Volcano plot was drawn by the ggplot2 package of R.

Survival analysis

We have adopted Kaplan–Meier analysis to evaluate the relations between DEGs and OS. And Tarone-Ware test was applied for comparing the differences. P < 0.05 was considered statistically significant.

PPI construction and hub genes selection

The String database ( was used to analyze the molecular interactions of DEGs and construct a PPI network visualized by Cytoscape (v3.7.2). Then we used the MCODE tool to find densely connected hub genes based on topology.

Functional enrichment analysis

Gene ontology (GO) analysis was performed by DAVID (, comprising three aspects: biological process, cellular component, and molecular function. And Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was applied for searching functional pathways of hub genes. P < 0.05 was deemed statistically significant.

Abundance analysis of tumor-infltrating immune cells and immunophenoscore (IPS) analysis

CIBERSORT algorithm ( was used to estimate the composition of 22 immune cell types of PRCC. The percentages of these types were visualized in a bar graph by using ggplot2 package of R. Immunophenogram of TCIA (, influenced by four factors: MHC molecules, immunomodulators, effector cells and suppressor cells, was applied to evaluating the immune properties of PRCC. And IPS scores are positively associated with the immunogenicity.

Nomogram building

To construct a nomogram to predict patients’ OS and PFS of PRCC, the C index was used to evaluate the discriminating ability, and the calibration chart was drawn to evaluate the accuracy of the nomogram. The rms package of R was used to build the nomogram.

Statistical analysis

Statistical differences between 2 groups were identified by unpaired Student’s t test or Mann Whitney U test. 1-way ANOVA test and Kruskal-Wallis test were applied when more than 2 groups were compared. Based on accessible clinical data, we used univariate and multivariate cox regression analysis to identify the prognostic values of certain hub genes. Kaplan–Meier analysis was used to estimate the OS and PFS. All statistical analysis was conducted by software packages of R (v3.6.3).


Immune/stromal score were significantly associated with T pathological stage and PRCC subtypes

We downloaded 323 and 34 PRCC samples from the TCGA and GEO database. Gene expression patterns and clinical characteristics were included. Based on the ESTIMATE algorithm, immune score (− 1952.04 to 3372.51) and stromal score (− 1821.04 to 1469.55) of 291 samples of TCGA and 34 samples of GEO were calculated. We attempted to analyze the relations between immune/stromal score and TNM, stage, subtype and grade (Fig. 1a-e). Results showed the differences in the stromal score of T pathological stages were significant (T3–4 > T1–2, P < 0.05). And the average immune scores of type 2 were higher than type 1 (P < 0.05). Similarly, type 2 had higher stromal scores compared to type 1 (P < 0.05, Fig. 1f, g). Finally, we divided the samples into two groups on the basis of median immune/stromal score, and revealed stromal score was significantly associated with PFS (P = 0.021, Fig. 1h). However, no significant differences were found in OS via Kaplan–Meier analysis.

Fig. 1

The correlations between the stromal score and clinical characteristics, and OS analysis. ae. The relationship between stromal score and T stage, N stage, M stage, clinical stage and grade. f, g. The relationship between immune/stromal score and PRCC subtypes. h. Kaplan-Meier analysis of stromal score with PRCC patients

Identification of DEGs with immune/stromal score

To determine whether all genes were associated with immune/stromal score, gene expression profile of 291 patients were analyzed. We mapped volcano plots with the cutoffs of |Log2FC| > 1 and FDR < 0.05 (Fig. S1A, B). Nine hundred eighty-three upregulated and 6 downregulated genes were screened by comparing the high and low immune scores. Meanwhile, we compared 2 groups of cases dichotomized by high and low stromal scores, and extracted 1142 upregulated and 27 downregulated genes. Venn diagrams identified 763 upregulated common genes and 4 downregulated common genes (Fig. S1C, D, Supplementary Table S1). GO analysis showed 763 common upregulated genes were mostly associated with immune and inflammatory response, chemokine activity, IgG binding and cell adhesion molecule binding (Fig. S1 E-G).

Association of upregulated common genes with OS

In order to obtain the relationships between each gene and OS, Kaplan–Meier was applied for survival analysis of 763 upregulated genes, and the impact of 120 genes (Supplementary Table S2) on OS was statistically significant.

Selection of molecular complexes from protein interaction networks

To study protein interactions of 120 genes, String tool was performed to structure a PPI network comprised of 84 nodes and 297 edges. Then we used MCODE to select prominent modules: module 1, composed of 17 hub genes (CXCR5, CD19, PDCD1, IL21R, GZMB, TNFRSF9, CD38, CXCL10, SELL, LAG3, CD44, CD80, CCL21, ITGA4, CCL19, IL6, and CXCL13, Fig. S2A) and module 2, including 5 hub genes (BLK, CD79A, FCRLA, MS4A1 and POU2AF, Fig. S2B), which were strongly associated with other genes, suggesting that they might play an important role in PRCC.

Functional enrichment analysis of valuable genes

To understand whether 22 genes have an impact on the TME, we performed KEGG analysis by DAVID, and the results showed that these genes mainly participated in the following pathways: cytokine-cytokine receptor interaction, chemokine signaling pathway and primary immunodeficiency (Fig. S2C).

Construction and validation of a survival model to predict prognosis of PRCC patients

Three hundred twenty-three PRCC samples from TCGA and 34 samples from GEO were used as train and test group. Univariate cox regression analysis has approved all of 22 genes were significantly with OS (Supplementary Table S3). Then, we applied multivariate cox regression analysis to confirm 4 genes (CD79A, CXCL13, IL6 and CCL19) significantly related with patients’ prognosis (Supplementary Table S3 and S4) and constructed a formula for risk score calculation after extracting the coefficients from the results: expression level of CD79A * (− 0.24020) + CXCL13 * 0.22031+ IL6 * 0.12025+ CCL19 * 0.17997). According to median risk score, 323 PRCC samples from TCGA were divided into high and low risk group, Kaplan–Meier survival curves both indicated patients’ OS and PFS in high risk group were worse than that in low risk group (p < 0.05, Fig. 2a, c). And area under curve (AUC) related to OS and PFS was 0.76 and 0.745 (Fig. 2b, d). Similarly, 34 samples from GEO were analyzed, and AUC was 0.708 (Fig. 2e, f). We depicted ROC curves to discuss the predictive ability of this gene model in different subgroups of PRCC, such as subtype, grade, pathologic stage, gender and age. And the results demonstrated the immune signature well performed in age > =65, female and type 1 PRCC (Fig. S3A-D). Compared with clinical characteristics (age, gender, TNM, grade and stage), the 4-genes model showed a great performance in PRCC prognosis (Fig. S3E, F). Then, the relationship between the TMB-signature and clinicopathological factors were analyzed, and the results indicated significant differences were observed in gender, TNM, stage and type (Fig. S4).

Fig. 2

Kaplan-Meier survival and ROC analysis.. OS and PFS were compared respectively between low and high risk group of TCGA samples (a, c), OS was compared between low and high risk group of GEO samples (e), ROC of 4-gene model of OS (c) and PFS (d) analysis of TCGA samples and OS analysis of GEO samples (f)

Immune infiltration and IPS analysis

In order to predict the composition of immune cells in TCGA samples, we divided patients into low and high risk group based on median risk score. The composition of 22 types of immune cells in 2 groups was depicted in the bar graph (Fig. 3a) which indicated low risk group was connected with recruitment of plasma cell, resting memory CD4+ T cell, macrophage M0/2, and resting mast cell (P < 0.05). We used IPS to predict the potential response of ICI (immune checkpoint inhibitors) for PRCC patients and investigated the association between IPS and 4-genes model of PRCC. Statistical analysis indicated that the IPS of low risk group was higher than high risk group (Fig. 3b). Conclusively, low risk group of immune-associated signature might have stronger immunogenicity.

Fig. 3

Immune cell composition analysis and the difference analysis of IPS. a Box plot of 22 immune cells proportion between low and high risk group of 4-gene model. b the difference analysis of IPS between low and high risk group of 4-gene model

Construction of predictive nomogram

By comparing the C-index and removing invalid variables, we constructed a nomogram to predict patients’ OS and PFS of PRCC which combined with two independent factors including stage and risk score (Fig. 4). The calibration curve showed the nomogram had the best performance with the C-index of 0.855 and 0.834. Conclusively, this nomogram could be a potential tool to evaluate the outcomes of PRCC.

Fig. 4

Nomogram for predicting OS (a) and PFS (c) for PRCC patients of TCGA based on risk score and stage. The calibration plot for internal validation of the nomogram (b, d)


Recently, TME has been identified to have profound significance and effect on the occurrence, progression, and treatment of tumors, and it is extremely necessary to study its related mechanism. Study has shown changes in mesenchymal stromal cell differentiation promoted the progression of multiple tumors by altering the tumor microenvironment [15]. Maeda-Otsuka et al. demonstrated hypoxia could accelerate the proliferation and migration of angiosarcoma by regulating tumor microenvironment [16]. Bader et al. reported that TME was closely linked to immunotherapy, and the understanding of TME could better guide treatment and achieve precision therapy [17]. Meanwhile, study indicated that immunotherapy was one of the promising options for advanced stage PRCC patients [18]. This study aimed to identify genes in connection with tumor microenvironment and significantly with OS in PRCC patients.

Firstly, we discovered that the differences in the stromal score of T pathological stages were statistically significant, indicating that stromal cells of TME were meaningful in the correlation of T pathological stages. Relevant study has demonstrated that changes in stromal components, such as the presence of myofibroblasts, initiated tumor invasion, and metastasis [19]. Delahunt reported type 1 PRCC inclined to be localized in kidney whereas type 2 generally tended to invade the surrounding organs and had a poor prognosis [20]. Results showed the mean immune/stromal scores of type 2 PRCC cases were significantly higher than type 1, which indicated higher immune/stromal scores were associated with type 2 PRCC, and might be connected with more aggressive biological behaviors. GO analysis of genes co-upregulated in the immune/stromal groups indicated that they were involved in immune and inflammatory responses, suggesting that most of the genes were related to TME. Then, Kaplan–Meier analysis of 763 upregulation genes showed that 120 genes were associated with PRCC prognosis. Combined with PPI network construction and significant modules extraction, we obtained 22 hub genes. Then KEGG enrichment analysis showed that they were active in immune-related pathways, such as cytokine-cytokine receptor interaction, chemokine signaling pathway and primary immunodeficiency. Univariate and multivariate regression analysis put forward a four-genes (CD79A, CXCL13, IL6 and CCL19) survival model based on 323 PRCC samples from the TCGA database, of which AUC was 0.76 of OS and 0.745, and confirmed by a set of GEO data, with the AUC of 0.708. Based on risk score and pathologic stage, we constructed two nomograms separately associated with patients’ OS and PFS. According to IPS program, results showed IPS was significantly higher in low risk group, which indicated lower expression level of 4 immune-related genes were connected with weaker immunogenicity.

Among the 4 genes, Certain studies have demonstrated that CCL19 could act as an immunomodulator by activating dendritic cells, T and B cells in secondary lymphoid tissue to modulate primary (or secondary) adaptive immune responses [21]. Besides, overexpression of CCL19 was discovered to be implicated in tumor progression in cervical cancer, but might be contribute to anti-vascular treatment in colorectal cancer through inhibiting angiogenesis [22, 23]. Numerous studies have shown that CXCL13, as a chemokine secreted by stromal cells was implicated in the occurrence, invasion, and lymph node metastasis of tumors by binding to its receptor CXCR5 [24,25,26,27]. Rachana et al. reported it could be a novel and compelling target for prostate cancer therapy [28]. Several studies indicated that CD79A was connected with aggressive hematological malignancy, could be a popular marker in the detection and treatment of hematopathy such as Hodgkin’s lymphoma and B-cell lymphoma [29,30,31]. IL6, a proinflammatory cytokine, has been reported to be influenced by cancer-associated fibroblasts (CAFs) in TME, and participated in the process of human tumor immunity improvement, tumor metastasis and colonization [32,33,34].

The interaction of PRCC and its TME affected the whole process from tumor occurrence and progress to metastasis and recurrence, which created plenty of opportunities for diagnosis, treatment and prognosis of PRCC patients. In present study, we attempted to find tumor microenvironment-genes associated with PRCC patients’ survival. Our results may provide some related data to future research on the relations of PRCC and TME. However, the mechanism of PRCC and its microenvironment could be extremely complex, our analysis based on the TCGA and GEO database and bioinformatic tools is only one part of it. Further research and analysis based on large samples will be essential.


By analyzing the TCGA and GEO database and applying ESTIMATE algorithm and a series of bioinformatic methods, we obtained tumor microenvironment associated genes (CD79A, CXCL13, IL6 and CCL19), which were related to the clinical outcome of PRCC patients. These genes could be used as biomarkers for predicting PRCC prognosis.

Availability of data and materials

The datasets analyzed during the current study are available in the TCGA and GEO repository. [,].



Papillary renal cell carcinoma


The Cancer Gene Atlas


Gene Expression Omnibus


Protein-protein interaction


Overall survival


Progression-free survival


The Cancer Immunome Atlas


Differentially expressed genes


Renal cell carcinoma


World Health Organization


Tumor microenvironment


Gene ontology


Kyoto Encyclopedia of Genes and Genomes




Immune checkpoint inhibitors


Receiver operating characteristic


Area under curve


Cancer-associated fibroblasts


  1. 1.

    Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.

  2. 2.

    Znaor A, Lortet-Tieulent J, Laversanne M, Jemal A, Bray F. International variations and trends in renal cell carcinoma incidence and mortality. Eur Urol. 2015;67(3):519–30.

    Article  PubMed  Google Scholar 

  3. 3.

    Sukov WR, Lohse CM, Leibovich BC, Thompson RH, Cheville JC. Clinical and pathological features associated with prognosis in patients with papillary renal cell carcinoma. J Urol. 2012;187(1):54–9.

    Article  PubMed  Google Scholar 

  4. 4.

    Vrdoljak E, Ciuleanu T, Kharkevich G, Mardiak J, Mego M, Padrik P, et al. Optimizing treatment for patients with metastatic renal cell carcinoma in the central and eastern European region. Expert Opin Pharmacother. 2012;13(2):159–74.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Moch H, Cubilla AL, Humphrey PA, Reuter VE, Ulbright TM. The 2016 WHO classification of Tumours of the urinary system and male genital organs-part a: renal, penile, and testicular Tumours. Eur Urol. 2016;70(1):93–105.

    Article  PubMed  Google Scholar 

  6. 6.

    Alomari AK, Nettey OS, Singh D, Kluger H, Adeniran AJ. Clinicopathological and immunohistochemical characteristics of papillary renal cell carcinoma with emphasis on subtyping. Hum Pathol. 2015;46(10):1418–26.

    Article  PubMed  Google Scholar 

  7. 7.

    Arneth B. Tumor Microenvironment. Medicina-Lithuania. 2020;56(1).

  8. 8.

    Noy R, Pollard JW. Tumor-associated macrophages: from mechanisms to therapy (vol 41, pg 49, 2014). Immunity. 2014;41(5):866.

    CAS  Article  Google Scholar 

  9. 9.

    Miyai Y, Esaki N, Takahashi M, Enomoto A. Cancer-associated fibroblasts that restrain cancer progression: hypotheses and perspectives. Cancer Sci. 2020;111(4):1047–57.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Eich M-L, Chaux A, Rodriguez MAM, Guner G, Taheri D, Pena MDCR, et al. Tumour immune microenvironment in primary and metastatic papillary renal cell carcinoma. Histopathology. 2020;76(3):423–32.

    Article  PubMed  Google Scholar 

  11. 11.

    Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4(1):2612.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Li X, Gao Y, Xu Z, Zhang Z, Zheng Y, Qi F. Identification of prognostic genes in adrenocortical carcinoma microenvironment based on bioinformatic methods. Cancer Med. 2020;9(3):1161–72.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Luo YX, Zeng GH, Wu S. Identification of microenvironment-related prognostic genes in bladder cancer based on gene expression profile. Front Genet. 2019;10.

  14. 14.

    Zhao X, Hu D, Li J, Zhao G, Tang W, Cheng H. Database Mining of Genes of Prognostic Value for the Prostate Adenocarcinoma Microenvironment Using the Cancer Gene Atlas. BioMed Res Int. 2020;2020:5019793.

  15. 15.

    Hughes RM, Simons BW, Khan H, Miller R, Kugler V, Torquato S, et al. Asporin restricts Mesenchymal stromal cell differentiation, alters the tumor microenvironment, and drives metastatic progression. Cancer Res. 2019;79(14):3636–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Maeda-Otsuka S, Kajihara I, Tasaki Y, Yamada-Kanazawa S, Sakamoto R, Sawamura S, et al. Hypoxia accelerates the progression of angiosarcoma through the regulation of angiosarcoma cells and tumor microenvironment. J Dermatol Sci. 2019;93(2):123–32.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Bader JE, Voss K, Rathmell JC. Targeting metabolism to improve the tumor microenvironment for cancer immunotherapy. Mol Cell. 2020;78(6):1019–33.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Courthod G, Tucci M, Di Maio M, Scagliotti GV. Papillary renal cell carcinoma: a review of the current therapeutic landscape. Crit Rev Oncol Hematol. 2015;96(1):100–12.

    Article  PubMed  Google Scholar 

  19. 19.

    De Wever O, Mareel M. Role of tissue stroma in cancer cell invasion. J Pathol. 2003;200(4):429–47.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Delahunt B, Eble JN, McCredie MRE, Bethwaite PB, Stewart JH, Bilous AM. Morphologic typing of papillary renal cell carcinoma: comparison of growth kinetics and patient survival in 66 cases. Hum Pathol. 2001;32(6):590–5.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Hwang H, Shin C, Park J, Kang E, Choi B, Han J-A, et al. Human breast cancer-derived soluble factors facilitate CCL19-induced chemotaxis of human dendritic cells. Sci Rep. 2016;6(1).

  22. 22.

    Zhang X, Wang Y, Cao Y, Zhang X, Zhao H. Increased CCL19 expression is associated with progression in cervical cancer. Oncotarget. 2017;8(43):73817–25.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Xu ZQ, Zhu CC, Chen C, Zong YP, Feng H, Hui D, et al. CCL19 suppresses angiogenesis through promoting miR-206 and inhibiting met/ERK/Elk-1/HIF-1 alpha/VEGF-A pathway in colorectal cancer. Cell Death Dis. 2018;9(10):974.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Hussain M, Adah D, Tariq M, Lu Y, Zhang J, Liu J. CXCL13/CXCR5 signaling axis in cancer. Life Sci. 2019;227:175–86.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Zhu Z, Zhang X, Guo H, Fu L, Pan G, Sun Y. CXCL13-CXCR5 axis promotes the growth and invasion of colon cancer cells via PI3K/AKT pathway. Mol Cell Biochem. 2015;400(1–2):287–95.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Legler DF, Loetscher M, Roos RS, Clark-Lewis I, Baggiolini M, Moser B. B cell-attracting chemokine 1, a human CXC chemokine expressed in lymphoid tissues, selectively attracts B lymphocytes via BLR1/CXCR5. J Exp Med. 1998;187(4):655–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Biswas S, Sengupta S, Chowdhury SR, Jana S, Mandal G, Mandal PK, et al. CXCL13-CXCR5 co-expression regulates epithelial to mesenchymal transition of breast cancer cells during lymph node metastasis (vol 143, pg 265, 2014). Breast Cancer Res Treat. 2016;155(3):615–6.

    Article  PubMed  Google Scholar 

  28. 28.

    Garg R, Blando JM, Perez CJ, Abba MC, Benavides F, Kazanietz MG. Protein kinase C epsilon cooperates with PTEN loss for prostate tumorigenesis through the CXCL13-CXCR5 pathway. Cell Rep. 2017;19(2):375–88.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Hoeller S, Zihler D, Zlobec I, Obermann EC, Pileri SA, Dirnhofer S, et al. BOB.1, CD79a and cyclin E are the most appropriate markers to discriminate classical Hodgkin's lymphoma from primary mediastinal large B-cell lymphoma. Histopathology. 2010;56(2):217–28.

    Article  PubMed  Google Scholar 

  30. 30.

    Mason DY, Cordell JL, Brown MH, Borst J, Jones M, Pulford K, et al. CD79a: a novel marker for B-cell neoplasms in routinely processed tissue samples. Blood. 1995;86(4):1453–9.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Naylor TL, Tang H, Ratsch BA, Enns A, Loo A, Chen L, et al. Protein kinase C inhibitor Sotrastaurin selectively inhibits the growth of CD79 mutant diffuse large B-cell lymphomas. Cancer Res. 2011;71(7):2643–53.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Kato T, Noma K, Ohara T, Kashima H, Katsura Y, Sato H, et al. Cancer-associated fibroblasts affect Intratumoral CD8(+) and FoxP3(+) T cells via IL6 in the tumor microenvironment. Clin Cancer Res. 2018;24(19):4820–33.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Cho H, Seo Y, Loke KM, Kim S-W, Oh S-M, Kim J-H, et al. Cancer-stimulated CAFs enhance monocyte differentiation and Protumoral TAM activation via IL6 and GM-CSF secretion. Clin Cancer Res. 2018;24(21):5407–21.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Toyoshima Y, Kitamura H, Xiang H, Ohno Y, Homma S, Kawamura H, et al. IL6 modulates the immune status of the tumor microenvironment to facilitate metastatic colonization of colorectal cancer cells. Cancer Immunol Res. 2019;7(12):1944–57.

    CAS  Article  PubMed  Google Scholar 

Download references


Not applicable.


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information




LL and HZ analyzed and interpreted the data regarding the papillary renal cell tumor. LL and HS prepared the manuscript and revised important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Liang Luo.

Ethics declarations

Ethics approval and consent to participate

All data of the study were acquired from public database so that no more ethical approval was needed.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1

Comparison of differentially expressed genes. A. Volcano plots of differentially expressed genes based on immune score. B. Volcano plots of differentially expressed genes based on stromal score. C. Venn plots of co-downregulated genes. D. Venn plots of co-upregulated genes. GO analysis of 763 co-upregulated genes. E. Biological process. F. Cellular component. G. Molecular function.

Additional file 2: Supplementary Table S1

763 upregulated common genes and 4 downregulated common genes between TCGA and GEO samples.

Additional file 3: Supplementary Table S2

120 genes associated with PRCC overall survival.

Additional file 4: Fig. S2

PPI networks of significant module. A. module1, B. module2. The node color changes gradually from pink to red indicating the ascending order of the degree of the genes. C KEGG pathway analysis for 17 hub genes.

Additional file 5: Supplementary Table S3

Univariate cox proportional hazard regression analysis of 22 genes based on TCGA database.

Additional file 6: Supplementary Table S4

Multivariate cox regression results of 4-gene prognostic model.

Additional file 7: Fig. S3

ROC curve analysis. The predictive ability of 4-gene model in different ages (A), genders (B), stages (C) and subtypes (D) of PRCC. AUCs of 4-gene model were compared with age, gender, TNM, stage, tumor grade of TCGA (E) and GEO (F).

Additional file 8: Fig. S4

The association between 4-gene model and clinicopathologic features. A. age, B. gender, C. T, D. M, F. N, G. stage, H. subtype.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Luo, L., Zhou, H. & Su, H. Identification of 4-genes model in papillary renal cell tumor microenvironment based on comprehensive analysis. BMC Cancer 21, 553 (2021).

Download citation


  • Papillary renal cell carcinoma
  • Tumor microenvironment
  • Hub genes
  • Prognosis