Fig. 4

WGCNA analysis results. A Clustering results. Cluster dendrogram and the co-expression network modules identified in average linkage hierarchical clustering of DEGs based on topological overlaps. Each branch within the dendrogram represents a single gene. Height represents the Euclidean distance. Each color indicates a single module containing weighted co-expressed genes; B Weighted network heat map of all genes; C Eigengene adjacency heatmap. The heatmap illustrates the relationship among the distinctive co-expression modules; D Heat map of correlation between the modules and the phenotypes. Each row represents a color module, and each column indicates a clinical trait. Each cell contains the R2 values of Pearson’s correlations between the modules and the clinical features, and the corresponding P-values are inside the parentheses. The gradient color of each cell indicates the R2 values of Pearson’s correlations (red = 1, blue =  − 1). DEGs differentially expressed genes; ME module eigengene