BackgroundThe detection rate of thyroid cancer (TC) has been continuously improved due to the development of detection technology. Compared with other cancers, the gene profile plays a more prominent role in the diagnosis and treatment of TC. MethodsFour datasets from Gene Expression Omnibus (GEO) was used to perform transcriptomic profile analysis. The overlapping differentially expressed genes (DEGs) were analyzed by R package “limma” and “RobustRankAggreg”. Then the hub genes, which had a higher degree, were identified by the protein-protein interaction (PPI) network. Gene expression analysis between the TC and normal data, the disease-free survival (DFS) analysis of TC patients from Gene Expression Profiling Interactive Analysis (GEPIA) database, function analysis, and immunohistochemistry (IHC) were performed to verify the importance of the hub genes.ResultsA total of 80 DEGs (34 upregulated and 46 downregulated) were identified. Then FN1, TIMP1, ITGA2, and KIT were considered hub genes, which had a high degree of connectivity in the PPI network. GEPIA identified that FN1, TIMP1, and ITGA2 were upregulated, and KIT was downregulated. Upregulations of FN1 expression (P=0.024) and ITGA2 expression (P=0.029) and downregulation of KIT expression (P=0.012) increased risks of decreased DFS for patients. IHC showed that the expression of FN1, TIMP1, and ITGA2 protein were upregulated, while the expression of KIT protein was downregulated in the TC clinical specimens. Besides, five hub genes were enriched in the PI3K-Akt signaling pathway and ECM-receptor interaction.ConclusionsIn summary, these hub genes were potential biomarkers of diagnosis and prognosis of TC.