Group differences were evaluated by unpaired Student's t-test or one-way analysis of variance (ANOVA), where appropriate. Fisher Exact test was used to measure the gene-enrichment in annotation terms. E-value comparing common motif and associated transcriptional factors was calculated based on the Wasserman and Sandelin method (Wasserman and Sandelin, 2004). P-values within WGCNA were calculated based on a regression-based p-value for assessing the statistical significance between the matrix X and the sample trait T, where matrix X is ‘network nodes times column indices’ of all genes (X = n × m).