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.

The ongoing COVID-19 pandemic is one of the biggest health and societal challenges of the recent decades. Among the causes of mortality triggered by SARS-CoV-2 infection, the presence of an inflammatory "cytokine storm" (CS) at later stages of the disease has been found to play a determinant role. Here, we used available transcriptomic data from the bronchoalveolar lavage fluid (BALF) of COVID-19 patients suffering from a CS to obtain gene-signatures associated to this pathological process. Using these signatures, we interrogated the Connectivity MAP (CMap) dataset that contains the effects of over 5,000 small molecules on the transcriptome of human cancer cell lines, and looked for molecules which effects on transcription mimic or oppose those associated to the CS. Consistent with their medical use, this analysis found a significant enrichment of glucocorticoids or inhibitors of the Janus Kinases (JAK) as drugs that could revert the CS. On the other hand, molecules that potentiate the immune response such as PKC activators are predicted to worsen the CS. Besides these expected findings, our analysis also reveals a potential effect of the antibiotic doxycycline or MAPK/RAF/MEK inhibitors in reverting the CS, or of topoisomerase inhibitors and the anti-alcohol abuse drug disulfiram in potentiating its effects. Finally, our analyses support that the gender-related differences in the severity of COVID-19 are related to the anti-inflammatory properties of female hormones. While acknowledging that this is an analysis based on limited available data, we decided to share it as a resource that might help others in the selection of drugs that could be tested in the context of experimental models of CS triggered by viral infections.