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

BMC Cancer. 2021 May 17;21(1):553. doi: 10.1186/s12885-021-08319-0.

Abstract

Background: 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.

Methods: 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).

Results: 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.

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

Keywords: Hub genes; Papillary renal cell carcinoma; Prognosis; Tumor microenvironment.

MeSH terms

  • Age Factors
  • Aged
  • Biomarkers, Tumor / genetics*
  • CD79 Antigens / genetics
  • Carcinoma, Renal Cell / diagnosis
  • Carcinoma, Renal Cell / genetics
  • Carcinoma, Renal Cell / immunology
  • Carcinoma, Renal Cell / mortality*
  • Chemokine CCL19 / genetics
  • Chemokine CXCL13 / genetics
  • Datasets as Topic
  • Feasibility Studies
  • Female
  • Gene Expression Regulation, Neoplastic / immunology
  • Humans
  • Interleukin-6 / genetics
  • Kaplan-Meier Estimate
  • Kidney / immunology
  • Kidney / pathology
  • Kidney Neoplasms / diagnosis
  • Kidney Neoplasms / genetics
  • Kidney Neoplasms / immunology
  • Kidney Neoplasms / mortality*
  • Male
  • Middle Aged
  • Models, Genetic
  • Neoplasm Grading
  • Nomograms*
  • Protein Interaction Maps / genetics
  • Protein Interaction Maps / immunology
  • ROC Curve
  • Risk Assessment / methods
  • Risk Assessment / statistics & numerical data
  • Sex Factors
  • Tumor Microenvironment / genetics*
  • Tumor Microenvironment / immunology
  • Up-Regulation / immunology

Substances

  • Biomarkers, Tumor
  • CCL19 protein, human
  • CD79 Antigens
  • CD79A protein, human
  • CXCL13 protein, human
  • Chemokine CCL19
  • Chemokine CXCL13
  • IL6 protein, human
  • Interleukin-6