Heatmap and hierarchical clustering of log2 transformed RPKM was performed and visualized in R using the gplots package
. A gene expression pattern transition matrix was created between Hanoch and 53 expression clusters. For comparing between expression clusters, each cluster of Hanoch was compared with cluster of 53. GO enrichment tests were performed by using the GO enrichment tool of the Blast2GO .
To obtain the functional classification of DE genes in the development stages as well as between the cultivars Fisher exact tests were used for the determination of significant biological, molecular or cellular components (Qvalue < 0.05).