Identification of potential biomarkers for papillary thyroid carcinoma by comprehensive bioinformatics analysis

To perform bioinformatics analysis on the papillary thyroid carcinoma (PTC) gene chip dataset to explore new biological markers for PTC. The gene expression profiles of GSE3467 and GSE6004 chip data were collected by GEO2R, and the differentially expressed genes (DEGs) were selected for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Protein-protein interaction (PPI) relationship analysis was achieved using STRING, and the hub genes were obtained using the Cytoscape software. GEPIA was used to validate the expressions of the hub genes in the normal and tumor tissues and to conduct survival analyses. Pertinent genetic pathology results were fetched using the HPA database. Finally, the key genes were clinically verified by reverse transcription-polymerase chain reaction. 97 genes were jointly up-regulated and 107 genes were jointly down-regulated in GSE3467 and GSE6004. GO function enrichment analysis revealed that the DEGs were involved in the regulation of calcium ion transport into cytosol, integrin binding, and cell adhesion molecule binding. KEGG pathway enrichment analysis indicated that the DEGs were chiefly associated with thyroid cancer and non-small cell lung cancer. According to the PPI network, 30 key target genes were identified. Only the expressions of ANK2, TLE1, and TCF4 matched between the normal and tumor tissues, and were associated with disease prognosis. When compared with the normal thyroid tissues, the protein and mRNA expressions of ANK2, TLE1, and TCF4 were down-regulated in PTC. Significant differences exist in overall gene expression between the thyroid tissues of patients with PTC and those of healthy people. Furthermore, the differential genes ANK2, TLE1, and TCF4 are expected to be reliable molecular markers for the mechanism study and diagnosis of PTC.