Fig. 4

Comparison of temporal gene regulation of EAE processes between NAWM and WML using spatial transcriptomics. A GO term analysis results of DE genes that have a higher (left) or lower (right) fold change in “WML” than in “NAWM” and have p < 0.05 in both subsets and Bonferroni-corrected p < 0.05 in at least one. Networks of “biological processes” GO terms are plotted by gene set overlap (width of grey lines scaled by Jaccard index). Images were created with R package enrichplot. B Log₂-fold changes from early-vs-late DGE-testing in “NAWM” (y axis) and “WML” (x axis) subset. Common DE genes between both subsets enriched in selected GO terms are plotted separately per functional modules with ellipse and regression line. A selection of the modules in C, ellipses are colored according to the 3 groups of modules. C Log₂-fold changes early-vs-late in the subsets “NAWM” and “WML”. Common DE genes between both subsets enriched in selected GO terms are plotted per functional modules of “biological processes” GO terms as box-and-whiskers in violin-plots. Genes may appear in multiple modules. The modules were grouped into glia-related, immune-related, or neuron-related processes. Under the modules the pairwise overlap of their gene sets is drawn as blue arcs (line width scaled by Jaccard index, minimum 0.15)