Gene set variation analysis (GSVA)
The R/Bioconductor software package GSVA32 was used as a non-parametric, unsupervised method to estimate enrichment of immune cell and pathway gene sets in RNAseq data as previously described12,33,34. Gene sets used as input for GSVA can be found in Supplementary Table 2.