*2.5. GO and KEGG Analysis of Di*ff*erentially Expressed Genes (DEGs)*

Based on the HISAT2 results, HTSeq software was used to calculate the expression level of each mRNA gene in samples. The differentially expressed genes (DEGs, *p* < 0.05) were identified using the DESeq R package and are displayed in the Supplementary Materials File S1. DEGs with a false discovery rate (FDR) < 0.05 and absolute value of fold change ≥ 1.5 were used to perform the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The GO functional enrichment of DEGs was analyzed using Metascape software [18]. The signaling pathway annotations was mapped in the KEGG database. The key KEGG signaling pathway was draught by the Pathview [19].
