In terms of the molecular function (MF), the enriched functional gene set associated with hyperthyroidism in children was most significant for function “hormone activity” with the p.adjust value of 1.63×10-10 and the corresponding gene set with mutations included 14 genes: IGF1, TG, CALCA, POMC, CGA, PTH, GHRL, FBN1, TRH, PRL, LEP, ADIPOQ, INS, GH1. Here, PRL is linked to hyperthyroidism.