Fig. 3

Gene co-expression network analysis looks for co-expressed gene modules and explores the association between gene networks and disease phenotypes. A: Heat map of all gene clustering. Based on the selected power values, a weighted co-expression network model was built to finally classify the 862 genes into 9 modules. B: Trait module association heat map (absolute value of correlation coefficient greater than or equal to 0.3 and pValue less than 0.05 is the threshold to filter the modules associated with each trait). C: Blue module characteristic gene histogram. D: KEGG enrichment analysis of genes in the blue module top20 (upregulated differences) bubble plots. Where the x-axis Enrichment Score is the enrichment score and the y-axis is the pathway information for top20. The larger the bubble the more entries contain the number of differential proteins, the bubble color changes from red–green–blue-purple, and the smaller its enrichment p-value, the more significant it is