The inclusive pathways were ordered according to the normalized enrichment score (NES) calculation and the top-ten most enriched KEGG-based annotations among these modules were constructed, meanwhile, mRNAs count as well as the − log10(P-value), namely their hypothesis testing parameters, were all displayed in Fig. 4, in which the thyroid cancer and p53 signaling pathway were significantly enriched with the maximum score values. Here, TP53 is linked to thyroid cancer.