The occurrences of transcription factor binding site motifs in the promoters of differential expression genes were scanned with the FIMO software [18], which computes a log-likelihood ratio score for each motif corresponding to each sequence position and converts this score to a p-value.