Computation and visualization were achieved using heatmap.plus package in R. Differentially expressed gene identification By using Cufflinks and cuffdiff, we identified pairwisely differentially expressed genes (DEG) based on the following criteria: fold change ≥2 (or log2 fold ≥ 1) and FDR ≤ 0.05.