The GO function enrichment analysis was performed by DAVID database and GOplot
The results of GO analyses of 627 common DEGs were shown in barplot (A) and bubbleplot (B), consisting of BP, CC and MF. Circle plot (C) and Gocluster plot (D) visualized the distribution of high-expression and low-expression genes contained in each term, as well as the Z-score.