Data normalization, scaling and 2 -way Anova was used to identify differentially expressed genes, using Partek Geno mics Suite (v6.5). Average linkage gene clustering was performed with a Euclidean distance using Hierarchical clustering. The genes with statistically significant changes between the treatments and strains were identified using Anova streamlined (Partek Genomic Suite (v6.5)). Only genes with a fold change higher than 1.3 and p-value lower than 0.05 were considered.