Fig. 1

Enrichment analysis and weighted correlation network analysis for differentially expressed genes (DEGs). a An UpSet plot illustrating the overlaps among DEGs identified by edgeR, Limma, and DESeq. b Bubble plots for enriched GO and KEGG terms. The x-axis represents the -log10P-value) of each term and the y-axis represents the number of genes in each term. c Dendrogram generated using the WGCNA. d PCC matrix between MEs and clinical traits. The PCC values range from − 1 to 1 depending on the strength of the relationship. A positive value indicates that the genes within a particular co-expression module increase as the clinical trait increases, whereas the opposite is true if the PCC is negative. Each PCC value is accompanied by the corresponding P-value in brackets