To calculate single-sample gene set enrichment, we used the GSVA program to derive the absolute enrichment scores of gene sets from several publications and previously experimentally validated gene signatures from MsigDB [74 (link)] as follows: tumor proliferation signature [27 (link)], tumor inflammation signature [75 (link)], cellular response to hypoxia, MYC targets, DNA replication, G2M checkpoint, P53 pathway, PI3K/AKT/mTOR pathway, IFNG signaling, genes up-regulated by ROS, DNA repair, degradation of ECM, collagen formation, angiogenesis, EMT markers, apoptosis, TGFB signaling, hallmark glycolysis, hypoxia signature and metabolism related pathways. To derive the GSVA score of each signature in each tumor sample, the normalized log2 (TPM+1) values were passed on as input for GSVA in the RNA-seq mode. Differentially enriched gene sets between the glycolysis high and low tumor groups were defined by GSVA adj.p < 0.05 (we used the limma R pachage because the GSVA scores were normally distributed around zero) [76 (link)].
Full text: Click here