Gene expression-based immune infiltration analyses of renal cancer and their associations with survival outcome
BMC Cancer volume 21, Article number: 595 (2021)
Renal cancer is a common malignant tumor with an increasing incidence rate.
In this study, based on the gene expression profiles, we analyzed the compositions of tumor-infiltrating immune cells (TIICs) in renal cancer and paracancerous samples using CIBERSORT. The proportions of 22 TIICs subsets in 122 paired renal carcinoma and paracancerous samples, and 224 Wilms tumor (WT) samples varied between intragroup and intergroup.
After analyzed the difference of TIICs composition between renal cancer and paired paracancerous samples, we found that M0 macrophages and CD8 T cells were significantly elevated, while naive B cells were significantly decreased in renal cancer samples compared with paracancerous samples. Survival analysis showed that high overall TIICs proportion, the low proportion of resting mast cells and the high proportion of activated memory CD4 T cells were associated with poor prognosis of renal cancer patients. In addition, 3 clusters were identified by hierarchical clustering analysis, and they presented a distinct prognosis. Cluster 1 had superior survival outcomes, while cluster 2 had an inferior survival outcome.
Our study indicated that overall TIICs proportion, certain TIICs subset proportion, including resting mast cells and activated memory CD4 T cells, and distinct cluster patterns were associated with the prognosis of renal cancer, which was significant for the clinical surveillance and treatment of renal cancer.
Renal cancer is a common malignant tumor with an increasing incidence rate over time [1, 2]. Wilms tumor (WT), also known as nephroblastoma, is the second most common solid abdominal organ tumor in children and the most frequent primary malignant renal cancer . It accounts for 6% of pediatric cancers, and an estimated 95% of renal cancer is WT among children . Clinically, WT or nephroblastoma is a common childhood tumor that is intimately linked to early kidney development and is often associated with persistent embryonic renal tissue and other kidney abnormalities, which is different from other types of renal cancer . As an embryonic tumor of the kidney, WT arises from metanephric mesenchyme and presents high a incidence rate among 2–3-year-old children . The treatment of WT is determined by several factors, such as genetic, pathologic, and demographic factors, and current therapies for WT consist of surgery, chemotherapy, and radiation therapy [3, 4]. Although the survival of WT has improved with the development of treatment strategies, the prognosis-related factors, including hematogenous and pulmonary metastasis, relatively high recurrence rate, and late effects, make the prospects of treatment noteworthy [7, 8]. Therefore, further investigation of the biological processes and underlying molecular mechanisms of WT may facilitate the prognosis improvement of patients.
Currently, increasing evidence has proved that the occurrence and development of tumors are associated with tumor cells as well as tumor microenvironment . The tumor microenvironment is composed of complex components, such as mesenchymal stem cells, fibroblasts, and immune cells . Among them, tumor-infiltrating immune cells (TIICs) are the pivotal components and are considered the leading players of the tumor microenvironment . Accumulating studies have shown that TIICs are closely associated with the development and survival outcome of cancers. For example, infiltrated macrophages were found to promote the progression of prostate cancer, and NK cells showed anti-tumor effects against tumor development . In colorectal cancer, TIICs were associated with patients’ clinical outcomes and considered key signatures for prognosis . The immunohistologic features and immunotype of TIICs played a crucial role in metastatic melanoma and were related to the survival outcome . In addition, dysfunctional infiltrating lymphocytes were observed in renal cell carcinoma (RCC), immunogenic renal cancer . It was indicated that tumors might exert an impairment effect on the immune system, and more studies on the associations between TIICs and tumor occurrence and prognosis should be helpful for the application of immunotherapeutic strategies on renal cancer.
Immunohistochemistry is a commonly used method to evaluate the TIICs in tumors. However, the results are often inaccurate due to the extensive-expression of markers in non-immune cells . In our research, the CIBERSORT algorithm was adopted to assess the TIICs subsets in renal cancer based on gene expression profiles retrieved from The Cancer Genome Atlas (TCGA). Meanwhile, we also analyzed the relationship of TIICs fraction and immune patterns with the survival outcome of renal cancer to explore the prognostic values of TIICs fraction and immune clusters in renal cancer.
We downloaded the gene expression profiles of renal cancer from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/, Affymetrix HG-U133A platform) and The Cancer Genome Atlas (TCGA, www.cancergenome.nih.gov, Illumina HiSeq platform). The TCGA dataset was composed of 889 renal cancer samples and 128 paracancerous samples, including chromophobe renal cell carcinoma, kidney renal clear cell carcinoma, and kidney renal papillary cell carcinoma samples. The GEO dataset (access no. GSE31403) included 224 Wilms tumor (WT) samples.
Tumor-infiltrating immune cells calculation
CIBERSORT (http://cibersort.stanford.edu), a deconvolution algorithm based on gene expression profiles , was used to calculate the relative proportions of 22 tumor-infiltrating immune cells (TIICs) subsets. CIBERSORT is able to evaluate the composition of TIICs with 547 barcode gene expression values, using P-value as a parameter for measurement of the confidence in results.
In addition, In the complex cancer immune microenvironment, cytotoxic T cells (Tc) and NK cells are two main effector cell types that can attack tumor cells directly. Upon exposure to transformed cells, cytotoxic T cells and NK cells secrete granzymes (a family of serine proteases) and perforin (a pore-forming protein) that will ultimately lead to target cell death. Thus, the local immune cytolytic activity can be quantified based on the transcript levels of perforin (PRF1) and granzyme A (GZMA). The first cytotoxin polymerizes and creates a channel in the membrane of the target cell. Through these pores, granzymes will then enter the cytoplasm and trigger a caspase cascade, composed of cysteine proteases that will ultimately lead to apoptosis . GZMA is a tryptase that induces caspase-independent programmed cell death, and PRF1 serves as a pore-forming enzyme that regulates entry of granzymes into target cells . Accordingly, it has been well-identified that the mean expression levels of genes GZMA and PRF1 represent the immune cytolytic activity of immune cells, which also reflect the fraction of TIICs [19,20,21,22,23,24]. Therefore, the expression levels of GZMA and PRF1 were determined to comprehensively assess the immune cytolytic activity and evaluate TIICs composition.
All 1241 samples were grouped according to results of CIBERSORT with P-values greater than or equal to 0.05 and less than 0.05, then the proportion of samples in each group and the average expression levels of genes GZMA and PRF1 were calculated. Seven samples of patiens without survival information were excluded, then based on kaplan-Meier method, survival analyses of the 1010 samples from TCGA, which were also stratified by P-value of 0.05 were performed by survival package (https://cran.r-project.org/web/packages/survival/) and survminer package (https://cran.r-project.org/web/packages/survminer/). With the relative proportions of 22 TIICs subsets as continuous variables, we calculated the hazard ratios (HR) using Cox regression analysis. We performed survival analyses to explore the effects of significant factors on survival outcomes.
Based on K-mean method, cluster the samples with the relative content of immune infiltrated cells by R software v3.5.2.
We used R language to calculate the proportions of 22 TIICs subpopulations in the 1241 samples from TCGA and GEO datasets. The correlation of the results from different datasets was analyzed by Pearson correlation analysis. The compositions of TIICs subpopulations in 122 paired cancer and paracancerous samples from TCGA were computed, and the difference in compositions was further analyzed. Besides, the fraction of TIICs subsets in WT samples from GEO was also calculated.
Then we stratified the 1241 samples by P-value of 0.05, and calculated the compositions of TIICs subsets as well as mean expression values of GZMA and PRF1 in the groups with P ≥ 0.05 and P < 0.05 respectively.
Performance of CIBERSORT for TIICs evaluation in renal cancer
The TIICs composition was analyzed by CIBERSORT. As shown in Fig. 1a, the samples from GEO and TCGA presented distinct proportions of 22 TIICs subsets. As the renal cancer samples in GEO were all WT patients, we speculated that the TIICs in WT samples were significantly different from those in other renal cancer samples. Then we analyzed the correlation of immune cell proportions in renal cancer samples from TCGA and GEO. It was revealed that the renal cancer samples from the two datasets had highly consistent immune cell proportions (Fig. 1b), illustrating that CIBERSORT could evaluate the fraction of TIICs independently of data sources and platforms.
Landscape of TIICs in renal cancer
The compositions of TIICs subsets in 122 paired paracancerous and renal cancer samples from the TCGA and 224 WT samples from GEO were calculated by CIBERSORT, respectively. As shown in Fig. 2a-c and Table S1 the intragroup and intergroup differences in TIICs fractions were manifest. Thus, we inferred that the proportion of TIICs subsets was an inherent characteristic which varied significantly among different individuals. In addition, we analyzed the difference in TIICs composition between renal cancer and paired paracancerous samples. As shown in Fig. 2d, M0 macrophages and CD8 T cells were significantly elevated, while naive B cells were significantly decreased in renal cancer samples compared with those in paired paracancerous samples, and the results were basically consistent with previous researches [25, 26]. Moreover, the data in TCGA was performed by CIBERSORT (Fig. S1), the results were consistent with previous research .
P-value of CIBERSORT represents the overall proportion of TIICs
It should be noted that instead of determining the actual values, CIBERSORT only calculates the relative ratios of TIICs subsets, which contributes to the dependency of results on each other. Therefore, we further analyzed the association between the P-value provided by CIBERSORT and TIICs composition. CIBERSORT P-value < 0.05 correlate with higher immune cell infiltrates, while no significantly different of the P-value ≥0.05 [28, 29]. Figure 3a showed that the proportions of samples with P-value < 0.05 and P-value ≥0.05 in TCGA and GEO were obviously different.
It has been well proved that the mean expression values of GZMA and PRF1 represent the immune cytolytic activity and are positively associated with TIICs proportions [19,20,21,22,23,24]. Accordingly, we futher evaluated the mean expression values, which represented the immune cytolytic activity in the system. After analyzing their mean expression values, we found that the samples with P-value < 0.05 had higher immune cytolytic activity in both TCGA and GEO cohorts (p < 2.22e-16 and p = 0.043, respectively, Fig. 3b and c). These results indicated that the proportion of TIICs in samples with P-value < 0.05 was higher in comparison to samples with P-value ≥0.05.
Overall TIICs proportion, resting mast cells and activated memory CD4 T cells are associated with prognosis of renal cancer
To investigate the effects of overall TIICs proportion or 22 individual TIICs subset on renal cancer prognosis, the survival analyses of 1010 renal cancer samples with survival information were performed. The survival curve of renal cancer samples stratified by P-value of 0.05 showed that samples with P-value < 0.05 presented inferior survival outcome in comparison to samples with P-value ≥0.05 (p < 0.0001, Fig. 3d). It was suggested that high TIICs proportions might be associated with poor prognosis of renal cancer patients.
Subsequently, to further explore the effects of 22 individual TIICs subset on renal cancer prognosis, the univariate Cox regression analysis was conducted with 22 TIICs subsets as continuous variables. The relevant 95% confidence intervals and Hazard Ratios (HRs) were shown in Fig. 4a. It was found that activated memory CD4 T cells and resting mast cells were significantly associated with the prognosis of renal cancer patients (HR = 9.4e + 05, p = 0.004 and HR = 3.7e – 03, p = 0.034, respectively). The survival curves showed a low proportion of resting mast cells and a high proportion of activated memory CD4 T cells were related to poor prognosis of renal cancer patients (both p < 0.0001).
Different immune clusters are associated with the prognosis of renal cancer
The above results indicated TIICs alteration might affect the prognosis of renal cancer patients. Therefore, we speculated whether different immune clusters could be identified with the TIICs data. Firstly, the optimal number of clusters was determined as 3 using the within cluster sum of square errors (WSS) method (Fig. S1). Subsequently, the hierarchical clustering of the samples was conducted by the Euclidean distance model (Fig. 5a). After analyzed the association of different immune patterns with prognosis, cluster 1 exhibited superior survival outcomes, while cluster 2 exhibited inferior survival outcomes (p = 0.0057, Fig. 5b). Moreover, the compositions of 22 TIICs subsets is significantly different among the cluster 1, cluster 2, and cluster 3 (Fig. 5c). In additional, the immunes cell type abundances differ between clusters was revealed in violin plot Fig. S2 and clinical differences of the three clusters was showed in the Fig. S3. The results suggested that the relative proportion of 20 types of immune cells was significantly different in 3 clusters (P < 0.05, ANOVA). In Cluster1, with the best prognosis, the relative proportion of these 4 types of immune cells, including T cells CD4 memory resting, Macrophages M1, Mast cells resting and Monocytes, and Monocyte, were significantly higher than that of the other two clusters.
The tumor microenvironment is considered as a complex “society”, and exerts a regulatory effect on tumor progression with the participation of multiple cell types and extracellular matrix . The various cell types in tumor microenvironment include mesenchymal stem cells, fibroblasts, endothelial cells, and immune cells . It is believed that the TIICs are indispensable members of the tumor microenvironment and reflect the host immune reaction to tumors . Accumulating evidence have proved that TIICs are associated with tumor invasion, metastasis, prognosis, and response to therapy [32, 33]. Given their important role, comprehensive researches on the TIICs in the tumor microenvironment may provide novel therapeutic approaches for tumors.
Numerous efforts have been devoted to exploring the specific role of TIICs in renal cancer by multiple approaches. Jensen HK et al. evaluated the TIICs in localized renal cell carcinoma by immunohistochemistry and found that neutrophils were independent prognostic signatures of survival outcome for localized renal cell carcinoma . Donskov F et al. identified the positive correlation of CD57+ NK cells and negative correlation of neutrophils with the prognosis of metastatic renal cell carcinoma by immunohistochemistry . Webster WS et al. proved that the infiltrated mononuclear cells were able to predict the survival outcome of renal cell carcinoma patients independently through hematoxylin and eosin staining and flow cytometric analysis . Unlike the above researches, the TIICs in renal cancer were analyzed using a distinct approach in our study. We assessed the 22 TIICs subpopulations composition of renal cancer and paracancerous samples from TCGA and GEO using CIBERSORT algorithm based on their gene expression profiles. The TIICs subsets proportions of renal cancer samples from TCGA and GEO showed consistent results after correlation analysis. It was suggested that CIBERSORT was able to assess the composition of the TIICs subset, which was independent of data platforms and sources. Besides, CIBERSORT could be performed to characterize cell heterogeneity using RNA mixtures from nearly any source . Many researchers used CIBERSORT for estimating the infiltration of immune cells, and there is no differences between using CIBERSORT on microarray datasets [38,39,40] and RNA sequencing datasets [41,42,43], even on both microarray datasets and RNA sequencing datasets , and chip platform . Meanwhile, it has been reported that CIBERSORT has been used in several previous studies to analyzed the TIICs in the renal cell carcinoma, which have identified that CD8+ T cells were associated with prolonged overall survival and the potential biomarker relaed to CD8+ T cells, respectively [46, 47]. However, several researches reported that elevated CD8+ T cells were negatively related to prognosis in the patients with glioma  and hepatocellular carcinoma . Besides, a different method of ESTIMATE algorithm has also been applied to analyze prognostic microenvironment-related genes and stromal and immune scores in Clear Cell Renal Cell Carcinoma . In this study, our CIBERSORT analysis identified the landscape of TIICs in renal cancer and demonstated high TIICs proportions might be associated with poor prognosis of renal cancer patients. These data providing new evidence of the application of CIBERSORT in renal cancer and provide valuable information of TIICs in renal cancer patients.
Then we investigated the association of 22 individual TIICs subpopulation with the survival outcome of renal cancer patients by Cox regression analysis and survival analysis. The results showed that elevated proportion of activated memory CD4 T cells and decreased proportion of resting mast cells were associated with poor prognosis of renal cancer. It was known that the activated memory CD4 T cells could secret interleukin (IL) 17, a proinflammatory cytokine that promotes the proliferation and growth of cervical cancer . It has been reported that IL-17 is able to promote the development of colorectal cancer and associate with poor prognosis . In the present study, a high proportion of activated memory CD4 T cells were related to poor prognosis of renal cancer patients. Hence, we speculated that the increased proportion of activated memory CD4 T cells contributed to inferior survival outcomes in renal cancer through enhanced secretion of IL-17, which might be one of the possible mechanisms and needed to be validated further in the future investigations. Mast cells, which reside in vascularized tissues, have two states, including resting condition and activated condition . They are involved in multiple tumor-related processes. In thyroid cancer, thyroid cancer cells activated the mast cells, which could secrete extensive proinflammatory, angiogenic, and growth factors and exert protumorigenic effects . We inferred that the decreased proportion of resting mast cells, which reflected the enhancement of activated mast cells activities, was associated with a poor prognosis due to the protumorigenic effects of activated mast cells. In addition, a previous study showed a four immune-related genes signature based on CXCL2, SEMA3G, PDGFD, and UCN is closely associated with the prognosis of renal clear cell carcinoma . In this study, we identified that TIICs proportion is closely correlated with the prognosis of renal cancer, which was diffent from the previous study.
In addition, a total of 3 immune clusters was identified by hierarchical clustering analysis. After analyzing their fractions, we found that the proportions of TIICs subpopulations presented remarkable difference in the different clusters. The survival analysis revealed that the survival outcomes were significantly diffenrt among these 3 clusters, in which cluster 1 exhibited superior survival outcomes, while cluster 2 exhibited inferior survival outcomes. These findings indicated that immune infiltrate is heterogeneously different in the renal cancer patients and the difference of immune infiltrate is closely associated with the prognosis of of renal cancer patients.
In summary, the analysis of 22 TIICs proportions in renal cancer samples showed elevated activated memory CD4 T cells proportion, and decreased resting mast cells proportion predicted poor prognosis in renal cancer. Different immune clusters also presented distinct survival outcomes. The results might unveil novel prognosis prediction and immunotherapeutic strategies on renal cancer.
Corgna E, Betti M, Gatta G, Roila F, Mulder PHMD. Renal cancer. Crit Rev Oncol Hematol. 2007;64(3):247–62. https://doi.org/10.1016/j.critrevonc.2007.04.007.
Weir HK, Johnson CJ, Ward KC, Coleman MP. The effect of multiple primary rules on cancer incidence rates and trends. Cancer Causes Control. 2016;27(3):377–90. https://doi.org/10.1007/s10552-016-0714-9.
Ehrlich PF (2017) Wilms tumor. Fundamentals of Pediatric Surgery. Springer International Publishing.
Barber TD, Derinkuyu BE, Wickiser J, Joglar J, Koral K, Baker LA. Wilms tumor: preoperative risk factors identified for intraoperative tumor spill. J Urol. 2011;185(4):1414–8. https://doi.org/10.1016/j.juro.2010.11.047.
Rivera MN, Haber DA. Wilms' tumour: connecting tumorigenesis and organ development in the kidney. Nat Rev Cancer. 2005;5(9):699–712. https://doi.org/10.1038/nrc1696.
Perotti D, Hohenstein P, Bongarzone I, Maschietto M, Weeks M, Radice P. Is wilms tumor a candidate neoplasia for treatment with wnt/β-catenin pathway modulators?—a report from the renal tumors biology-driven drug development workshop. Mol Cancer Ther. 2013;12(12):2619–27. https://doi.org/10.1158/1535-7163.MCT-13-0335.
Davidoff AM. Wilms' tumor. Adv Pediatr. 2012;59(1):247–67. https://doi.org/10.1016/j.yapd.2012.04.001.
Njuguna F, Martijn HA, Kuremu RT, Saula P, Kirtika P, Olbara G. Wilms tumor treatment outcomes: perspectives from a low-income setting. J Glob Oncol. 2017;3(5):555–62. https://doi.org/10.1200/JGO.2016.005389.
Liu X, Wu S, Yang Y, Zhao M, Zhu G, Hou Z. The prognostic landscape of tumor-infiltrating immune cell and immunomodulators in lung cancer. Biomed Pharmacother. 2017;95:55–61. https://doi.org/10.1016/j.biopha.2017.08.003.
Isaac A, Sharma B. Modulation of the tumor microenvironment for cancer treatment: a biomaterials approach. J Funct Biomat. 2015;6:81–103.
Ino Y, Yamazaki-Itoh R, Shimada K, Iwasaki M, Kosuge T, Kanai Y. Immune cell infiltration as an indicator of the immune microenvironment of pancreatic cancer. Br J Cancer. 2013;108(4):914–23. https://doi.org/10.1038/bjc.2013.32.
Gannon PO, Poisson AO, Delvoye N, Lapointe R, Mes-Masson AM, Saad F. Characterization of the intra-prostatic immune cell infiltration in androgen-deprived prostate cancer patients. J Immunol Methods. 2009;348(1-2):9–17. https://doi.org/10.1016/j.jim.2009.06.004.
Caro GD, Bergomas F, Grizzi F, Doni A, Bianchi P, Malesci A. Occurrence of tertiary lymphoid tissue is associated with t-cell infiltration and predicts better prognosis in early-stage colorectal cancers. Clin Cancer Res. 2014;20(8):2147–58. https://doi.org/10.1158/1078-0432.CCR-13-2590.
Erdag G, Schaefer JT, Smolkin ME, Deacon DH, Shea SM, Dengel LT. Immunotype and immunohistologic characteristics of tumor-infiltrating immune cells are associated with clinical outcome in metastatic melanoma. Cancer Res. 2012;72(5):1070–80. https://doi.org/10.1158/0008-5472.CAN-11-3218.
Thompson RH, Dong H, Lohse CM, Leibovich BC, Blute ML, Cheville JC. Pd-1 is expressed by tumor-infiltrating immune cells and is associated with poor outcome for patients with renal cell carcinoma. Clin Cancer Res. 2007;13(6):1757–61. https://doi.org/10.1158/1078-0432.CCR-06-2599.
Xiong Y, Wang K, Zhou H, Peng L, You W, Fu Z. Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 2018;7(9):4496–508. https://doi.org/10.1002/cam4.1745.
Trapani JA, Smyth MJ. Functional significance of the perforin/granzyme cell death pathway. Nat Rev Immunol. 2002;2(10):735–47. https://doi.org/10.1038/nri911.
Chowdhury D, Lieberman J. Death by a thousand cuts: granzyme pathways of programmed cell death. Annu Rev Immunol. 2008;26(1):389–420. https://doi.org/10.1146/annurev.immunol.26.021607.090404.
Constantinos R, Dimitrios C, Anestis M, Efstathiades C, Dimopoulos C, Zaravinos A. The expression and prognostic impact of immune Cytolytic activity-related markers in human malignancies: a comprehensive meta-analysis. Front Oncol. 2018;8:27.
Gong Z, Jia Q, Chen J, Diao X, Gao J, Wang X. Impaired cytolytic activity and loss of clonal neoantigens in elderly patients with lung adenocarcinoma. J Thorac Oncol. 2019;14(5):857–66. https://doi.org/10.1016/j.jtho.2019.01.024.
Gu XL, Boldrup LD, Coates PJ, Fahraeus R, Wang LX, Wilms TB. High immune cytolytic activity in tumor-free tongue tissue confers better prognosis in patients with squamous cell carcinoma of the oral tongue. J Pathol. 2019;5(4):240–7. https://doi.org/10.1002/cjp2.138.
Narayanan S, Kawaguchi T, Yan L, Peng X, Qi Q, Takabe K. Cytolytic activity score to assess anticancer immunity in colorectal cancer. Ann Surg Oncol. 2018;25(8):2323–31. https://doi.org/10.1245/s10434-018-6506-6.
Rooney M, Shukla S, Wu C, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1-2):48–61. https://doi.org/10.1016/j.cell.2014.12.033.
Tian T, Ji W, Yan H, Haoran L, Feng G, Youdong P. A comprehensive survey of immune cytolytic activity-associated gene co-expression networks across 17 tumor and normal tissue types. Cancers. 2018;10(9):307. https://doi.org/10.3390/cancers10090307.
Rohr-Udilova N, Klinglmüller F, Schulte-Hermann R, Stift J, Herac M, Salzmann M, et al. Deviations of the immune cell landscape between healthy liver and hepatocellular carcinoma. Sci Rep. 2018;8(1):6220. https://doi.org/10.1038/s41598-018-24437-5.
Zhong QY, Fan EX, Feng GY, Chen QY, Zhang GH. A gene expression-based study on immune cell subtypes and glioma prognosis. BMC Cancer. 2019;19(1):1116. https://doi.org/10.1186/s12885-019-6324-7.
Wang Y, Yin C, Geng L, Cai W. Immune infiltration landscape in clear cell renal cell carcinoma implications. Front Oncol. 2021;10:491621.
Raza AH, Leon C, Pharoah PDP, Florian M, Carlos C, Marc L. Patterns of immune infiltration in breast cancer and their clinical implications: a gene-expression-based retrospective study. PLoS Med. 2016;13:e1002194.
Yongfu X, Kang W, He Z, Linglong P, Wenxian Y, Zhongxue F. Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 2018;7:4496–508.
Albini A, Magnani E, Noonan DM. The tumor microenvironment: biology of a complex cellular and tissue society. Q J Nucl Med Mol Imaging. 2010;54(3):244–8.
Ungefroren H, Sebens S, Seidl D, Lehnert H, Hass R. Interaction of tumor cells with the microenvironment. Cell Commun Signal. 2011;9(1):18. https://doi.org/10.1186/1478-811X-9-18.
Jochems C, Schlom J. Tumor-infiltrating immune cells and prognosis: the potential link between conventional cancer therapy and immunity. Exp Biol Med. 2011;236(5):567–79. https://doi.org/10.1258/ebm.2011.011007.
Jiang B, Mason J, Jewett A, Liu ML, Man YG. Tumor-infiltrating immune cells: triggers for tumor capsule disruption and tumor progression? Int J Med Sci. 2013;10(5):475–97. https://doi.org/10.7150/ijms.5798.
Jensen HK, Donskov F, Marcussen N, Nordsmark M, Lundbeck F, Von DMH. Presence of intratumoral neutrophils is an independent prognostic factor in localized renal cell carcinoma. J Clin Oncol. 2009;27(28):4709–17. https://doi.org/10.1200/JCO.2008.18.9498.
Donskov F. Impact of immune parameters on long-term survival in metastatic renal cell carcinoma. J Clin Oncol. 2006;24(13):1997–2005. https://doi.org/10.1200/JCO.2005.03.9594.
Webster WS, Lohse CM, Thompson RH, Dong H, Frigola X, Dicks DL. Mononuclear cell infiltration in clear-cell renal cell carcinoma independently predicts patient survival. Cancer. 2006;107(1):46–53. https://doi.org/10.1002/cncr.21951.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7. https://doi.org/10.1038/nmeth.3337.
Stahl D, Gentles AJ, Thiele R, Gütgemann I. Prognostic profiling of the immune cell microenvironment in ewing’s sarcoma family of tumors. OncoImmunology. 2019;8(12):e1674113. https://doi.org/10.1080/2162402X.2019.1674113.
Xu F, Zhang H, Chen J, Lin L, Chen Y. Immune signature of t follicular helper cells predicts clinical prognostic and therapeutic impact in lung squamous cell carcinoma. Int Immunopharmacol. 2019;81:105932.
Zhou R, Zhang JW, Zeng DQ, Sun HY, Rong XX, Shi M. Immune cell infiltration as a biomarker for the diagnosis and prognosis of stage i-iii colon cancer. Cancer Immunol Immunother. 2018;68:433–42.
Chaumette B, Kebir O, Dion PA, Rouleau GA, Krebs MO. Reliability and correlation of mixture cell correction in methylomic and transcriptomic blood data. BMC Res Notes. 2020;13:74.
Diaz-Mejia JJ, Meng EC, Pico AR, MacParland SA, Ketela T, Pugh TJ. Evaluation of methods to assign cell type labels to cell clusters from single-cell RNA-sequencing data. F1000Res. 2019;8:ISCB Comm J-296.
Yazdani S, Callemeyn J, Gazut S, Lerut E, De Loor H, Wevers M. Natural killer cell infiltration is discriminative for antibody-mediated rejection and predicts outcome after kidney transplantation. Kidney Int. 2018;95:188–98.
Wang X, Wu Y, Wen D, Wu LY, Yang H. An individualized immune prognostic index is a superior predictor of survival of hepatocellular carcinoma. Med Sci Monit. 2020;26:e921786.
Tang C, Ma J, Liu X, Liu Z. Development and validation of a novel stem cell subtype for bladder cancer based on stem genomic profiling. Stem Cell Res Ther. 2020;11(1):457. https://doi.org/10.1186/s13287-020-01973-4.
Lin JX, Yu M, Xu X, Wang YT, Xing HT, An J. Identification of biomarkers related to CD8 + T cell infiltration with gene co-expression network in clear cell renal cell carcinoma. Aging (Albany NY). 2020;12(4):3694–712. https://doi.org/10.18632/aging.102841.
Zhang S, Zhang E, Long J, Hu Z, Peng J, Liu L. Immune infiltration in renal cell carcinoma. Cancer Sci. 2019;110(5):1564–72. https://doi.org/10.1111/cas.13996.
Chen B, Chen W, Jin J, Wang X, He Y. Data mining of prognostic microenvironment-related genes in clear cell renal cell carcinoma: a study with tcga database. Dis Markers. 2019:1–11.
Tartour E, Fossiez F, Joyeux I, Galinha A, Gey A, Claret E. Interleukin 17, a t-cell-derived cytokine, promotes tumorigenicity of human cervical tumors in nude mice. Cancer Res. 1999;59:3698.
Liu J, Duan Y, Cheng X, Chen X, Xie W, Long H. Il-17 is associated with poor prognosis and promotes angiogenesis via stimulating vegf production of cancer cells in colorectal carcinoma. Biochem Biophys Res Community. 2011;407(2):348–54. https://doi.org/10.1016/j.bbrc.2011.03.021.
Carla V, Nella P, Federica L, Gianni M. Tumor-associated mast cells in thyroid cancer. Int J Endocrinol. 2015;2015:1–8.
Gao X, Yang J, Chen Y. Identification of a four immune-related genes signature based on an immunogenomic landscape analysis of clear cell renal cell carcinoma. J Cell Physiol. 2020;235(12):9834–50. https://doi.org/10.1002/jcp.29796.
Ethics approval and consent to participate
Consent for publication
The authors declare that there are no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Relative proportions of immune cells in TCGA dataset analyzed by CIBERSORT. a The colum plot of 22 types of immune cells in TCGA renal cancer and paracancerous samples. b The box plot of 6 types of immune cells with significant differences between renal cancer and paracancerous samples.
. Selection of the optimal number of clusters.
. The violin plot of immune cell type abundances differ between clusters. Violin colors represent different clusters, and the vertical axis represents the relative proportion of immune cells.
. The violin plot of clinical characters in three clusters. Violin colors represent different clusters, and the vertical axis represents samples. a Age. b Sex. c Stage.
About this article
Cite this article
Chen, L., Yin, L., Qi, Z. et al. Gene expression-based immune infiltration analyses of renal cancer and their associations with survival outcome. BMC Cancer 21, 595 (2021). https://doi.org/10.1186/s12885-021-08244-2