Fig. 3

Pathway analysis. Enrichment of the top clusters was performed using IPA GO algorithms for all significant DEGs. a Enrichment of clusters for DEGs for T1 versus T0 (top) and T2 versus T0 (bottom). The threshold of minimum significance (−log(FDR)) was set to 1.3. The z-score predicts the direction of change for the function, and the threshold was set to an absolute value of z-score > 2. The height of the bars represents the probability that the DEGs in the dataset are associated with each pathway. The blue bars indicate predicted pathway activation, and orange indicate predicted inhibition. The orange points connected by a line represent the ratio which indicates the proportion of DEGs in the datasets that map to a given canonical pathway. b Canonical pathways associated with cytokine signaling. Enrichment of clusters associated with cytokine signaling was performed using IPA analysis (FDR < 0.05). The Benjamini–Hochberg method was used as an FDR adjustment for multiple testing correction. The top 20 enriched canonical pathways (ranked by average -log(FDR) value) are presented by showing the −log(FDR) and z-scores in heatmaps