Differentially expressed genes (DEGs) were identified through limma package (log2 fold change >0.585, false discovery rate (FDR) < 0.05), and protein-protein interaction (PPI) network was constructed using STRING. 
limma
Differentially expressed genes (DEGs) were identified through limma package (log2 fold change >0.585, false discovery rate (FDR) < 0.05), and protein-protein interaction (PPI) network was constructed using STRING. 
STRING
Moreover, perturbation and pathway enrichment analyses were performed through PerturbationAnalyzer in Cytoscape (iterative criteria <1×e-10) and clusterProfiler package (FDR <0.05), respectively.



PerturbationAnalyzer
Moreover, perturbation and pathway enrichment analyses were performed through PerturbationAnalyzer in Cytoscape (iterative criteria <1×e-10) and clusterProfiler package (FDR <0.05), respectively.



Cytoscape
Moreover, perturbation and pathway enrichment analyses were performed through PerturbationAnalyzer in Cytoscape (iterative criteria <1×e-10) and clusterProfiler package (FDR <0.05), respectively.



clusterProfiler