Bioinformatics and Connectivity Map Analysis Suggest Viral Infection as a Critical Causative Factor of Hashimoto’s Thyroiditis

Abstract

Hashimoto’s thyroiditis (HT) is a common autoimmune disease, with its prevalence rapidly increasing. Both genetic and environmental risk factors contribute to the development of HT. Recently, viral infection has been suggested to act as a trigger of HT by eliciting the host immune response and subsequent autoreactivity. We analyzed features of HT through bioinformatics analysis so as to identify markers of HT development. We accessed public microarray data of HT patients from the Gene Expression Omnibus (GEO) and obtained differentially expressed genes (DEGs) under HT. Gene Ontology (GO) and KEGG pathway enrichment analyses were performed for functional clustering of our protein-protein interaction (PPI) network. Utilizing ranked gene lists, we performed Gene Set Enrichment Analysis (GSEA) using the clusterprofiler R package. By comparing the expression signatures of the huge perturbation database with the queried rank-ordered gene list, connectivity map (CMap) analysis was performed to screen potential therapeutic targets and agents. The gene expression profile of the HT group was in line with general characteristics of HT. Biological processes related to the immune response and viral infection pathways were obtained for the upregulated DEGs. GSEA results revealed activation of autoimmune disease-related pathways and several viral infection pathways. Autoimmune disease and viral infection pathways were highly interconnected by common genes, while the HLA genes which are shared by both were significantly upregulated. CMap analysis suggested that perturbagens, including SRRM1, NLK, and CCDC92, have the potential to reverse the HT expression profile. Several lines of evidence suggested that viral infection and the host immune response are activated during HT. Viral infection is suspected to act as a key trigger of HT by causing autoimmunity. SRRM1, an alternative splicing factor, which responds to viral activity, might serve as potential marker of HT.