*4.4. Identification of Di*ff*erentially Expressed Genes*

To identify DEGs between two samples, the gene expression levels were quantified with the FPKM method. The read counts were adjusted with the edgeR program package, with one scaling-normalized factor for each sequenced library. The DEGs between two samples were analyzed with the DEGSeq R package (version 1.20.0). The *p*-values were adjusted according to the Benjamini and Hochberg method. A corrected *p*-value of 0.005 and a log<sup>2</sup> (fold-change) of 1 were set as the threshold for identifying significant DEGs. Significantly enriched GO terms and KEGG pathways were determined based on a corrected *p*-value ≤ 0.05. The GO functional enrichment and KEGG pathway enrichment analyses of the DEGs were completed with GOseq R packages and KOBAS (version 2.0) (http://kobas.cbi.pku.edu.cn/home.do), respectively.
