Fig. 4

Results of weighted gene coexpression network analysis. A Clustering dendrogram of identified coexpression gene modules in the radiogenomics analysis set represented in the color-coded row. B Bubble charts show the results of GSVA values and radscore correlation of the 9 modules, FDR less than 0.01 is considered as significant correlation of modules with radscore (modules on the right side of the red line). C A heatmap of the distribution of clinical features and gene set variation analysis (GSVA) scores of modules enrichment-related pathways in 132 patients in the radiogenomics analysis set. After performing a Pearson correlation, 759 pathways were selected as prognostically relevant (FDR < 0.01) and clustered according to their specific function. Patients were ranked according to their radscore, and pathway activity (GSVA values) demonstrated differences between high and low risk groups. D Bar plot shows the FDR values for the top four pathways enriched by each gene set. E Bubble diagram showing enriched pathways count and generation for 16 representative pathways (according to FDR) from 6 gene sets