Microarray Data of Lacrimal Gland Implicates Dysregulated Protein Processing in Endoplasmic Reticulum in Graves’ Ophthalmopathy

Graves’ ophthalmopathy (GO) has become one of the most common orbital diseases. Although some evidences announced the potential mechanism of pathological changes in extraocular muscle and orbital adipose tissue, little is known about that in lacrimal enlargement of GO patients. Thus, gene expression profiles of lacrimal gland derived from GO patients and normal controls were investigated using the microarray datasets of GSE105149 and GSE58331. The raw data and annotation files of GSE105149 and GSE58331 were downloaded from Gene Expression Omnibus (GEO) database. Bioinformatics including differentially expressed genes (DEGs), Gene Ontology, Kyoto Encyclopedia of Gene and Genome (KEGG) pathway, protein-protein interaction (PPI) network construction, hub gene identification, and gene set variation analysis (GSVA) were successively performed. A total of 173 overlapping DEGs in GSE105149 and GSE58331 were screened out, including 20 up-regulated and 153 down-regulated genes. Gene Ontology, KEGG and GSVA analyses of these DEGs showed that the most significant mechanism was closely associated with endoplasmic reticulum (ER). Moreover, we identified 40 module genes and 13 hub genes which were also enriched in the ER-associated terms and pathways. Among the hub genes, five genes including HSP90AA1, HSP90B1, DNAJC10, HSPA5, and CANX may be involved in the dysfunction of protein processing in ER. Taken together, our observations revealed a dysregulated gene network which is essential for protein processing in ER in GO patients. These findings provided a potential mechanism in the progression of lacrimal enlargement in GO patients, as a new insight into GO pathogenesis.