The unigene expression level was calculated using the fragments per kilo bases per million reads (FPKM) method [18]. The threshold for significantly differentially expressed genes (DEGs) was set at a fold change ≥2.0 and a FDR value <0.001. DEG functions were explored through GO and KEGG ...