To calculate the significance of differential expression of any given gene set M of size N, empirical distribution was determined. Specifically, differential expression of the set M was calculated as the average of the differential expression of genes of the set M, where the differential expression of each gene was quantified as: , which accounts for both fold change and P value of the differential expression. Then, a random gene set R of size N was sampled, and (i.e., average differential expression of genes in the set R) was calculated. Repeating this procedure 1E7 times, we constructed the empirical distribution of . P value of the significance of upregulation of M was calculated as , and significance of downregulation of M was calculated as
We compared cells of lymphoid lineage (i.e., CD4+ T, CD8+ T, B, and NK cells) between ESRRAi treatment and control (vehicle). The significance of the differential expression of a marker set is determined by the empirical distribution procedure described above.