For hypergeometric analysis, we used enrichGO and enrichDGN from the clusterProfiler and DOSE packages, respectively. Like the GSEA analysis, we defined the gene set ‘universe’ as all unique genes tested for differential expression.