Fig. 3

The whole analysis procedures: 1) 176 DEGs are obtained by paired-t test of 19777 genes expression of patients before and after treatment; 2) 219 gene co-expression modules are obtained by Weighted Gene Co-expression Network Analysis (WGCNA) [23], 3) 22 gene co-expression modules are significant correlated with treatment response traits; 4) 79 DEGs are refined by the 22 treatment response correlated modules as M_DEGs; 5) 287 GO pathways are significantly enriched (pathway size < 100 and p-value < 0.1) from 79 M_DEGs; 6) 30 GO pathways gene expression activity are significantly correlated with treatment response traits (|correlation|> 0.7; correlation p-value < 0.01); 7) 19 M_DEGs are refined by the 30 treatment response correlated GO pathways as GO_DEGs