Abstract
BACKGROUND: Tumor immunological heterogeneity potentially influences the prognostic disparities among patients with clear cell renal cell carcinoma (ccRCC); however, there is a lack of macroscopic imaging tools that can be used to predict immune-related gene expression in ccRCC. METHODS: A novel non-invasive radiogenomics biomarker was constructed for immune-related gene expression in ccRCC. First, 520 ccRCC transcriptomic datasets from The Cancer Genome Atlas (TCGA) were analyzed using a non-negative matrix decomposition (NMF) clustering to identify immune-related molecular subtypes. Immune-related prognostic genes were analyzed through Cox regression and Gene Set Enrichment Analysis (GSEA). We then built a risk model based on an immune-related gene subset to predict prognosis in patients with ccRCC. CT images corresponding to the ccRCC patients in The Cancer Imaging Archive (TCIA) database were used to extract radiomic features. To stratify immune-related gene expression levels, extracted radiogenomics features were identified according to standard consecutive steps. A nomogram was built to combine radiogenomics and clinicopathological information through multivariate logistic regression to further enhance the radiogenomics model. Mann-Whitney U test and ROC curves were used to assess the effectiveness of the radiogenomics marker. RESULTS: NMF methods successfully clustered patients into diverse subtypes according to gene expression levels in the tumor microenvironment (TME). The relative abundance of 10 immune cell populations in each tissue was also analyzed. The immune-related genomic signature (consisting of eight genes) of the tumor was shown to be significantly associated with survival in patients with ccRCC in TCGA database. The immune-related genomic signature was delineated by grouping the signature expression as either low- or high-risk. Using TCIA database, we constructed a radiogenomics biomarker consisting of 11 radiomic features that were optimal predictors of immune-related gene signature expression levels, which demonstrated AUC (area under the ROC curve) values of 0.76 and 0.72 in the training and validation groups, respectively. The nomogram built by combining radiomics and clinical pathological information could further improve the predictive efficacy of the radiogenomics model (AUC = 0.81, 074). CONCLUSIONS: The novel prognostic radiogenomics biomarker achieved excellent correlation with the immune-related gene expression status of patients with ccRCC and could successfully stratify the survival status of patients in TCGA database. It is anticipated that this work will assist in selecting precise clinical treatment strategies. This study may also lead to precise theranostics for patients with ccRCC in the future.