2.4. Gene Ontology and Pathway Enrichment Analysis

After identifying the different modules in the co-expression network, gene set enrichment and pathway enrichment analysis were performed to explore functional categories associated with modules. Gene ontology (GO) and pathway enrichment analyses were performed using ShinyGO [36], Metascape [37], and Parametric Gene Set Enrichment Analysis (PGSEA) [38] with a p-value cutoff <0.05. KEGG and REACTOME pathway databases were consulted during pathway enrichment analysis. Fisher’s exact test followed by Benjamini–Hochberg multiple testing correction was employed for selecting the GO category and the relevant pathway.