AIM: To uncover the underlying pathogenesis of thyroid associated ophthalmopathy (TAO) and explore potential biomarkers of this disease. METHODS: The expression profile GSE9340, which was downloaded from Gene Expression Omnibus database, included 18 specimens from 10 TAO patients and 8 hyperthyroidism patients without ophthalmopathy. The platform was HumanRef-8 v2 Expression BeadChip. Raw data were normalized using preprocess. Core package and the differentially expressed genes (DEGs) were identified based on t-test with limma package of R. Functional enrichment analyses were performed recruiting the DAVID tool. Based on STRING database, a protein-protein interaction (PPI) network was constructed, from which a module was extracted. The functional enrichment for genes in the module was performed by the BinGO plugin. RESULTS: In total, 861 DEGs (433 up-regulated and 428 down-regulated) between TAO patients and hyperthyroidism patients without ophthalmopathy were identified. Crucial nodes in the PPI network included TPX2, CDCA5, PRC1, KIF23 and MKI67, which were also remarkable in the module and all enriched in cell cycle process. Additionally, MKI67 was highly correlated with TAO. Besides, the DEGs of GTF2F1, SMC3, USF1 and ZNF263 were predicted as transcription factors (TFs). CONCLUSION: Several crucial genes are identified such as TPX2, CDCA5, PRC1 and KIF23, which all might play significant roles in TAO via the regulation of cell cycle process. Regulatory relationships between TPX2 and CDCA5 as well as between PRC1 and KIF23 may exist. Additionally, MKI67 may be a potent biomarker of TAO, and SMC3 and ZNF263 may exert their roles as TFs in TAO progression.