Abstract
Hashimoto's thyroiditis (HT) is an autoimmune disruption manifested by immune cell infiltration in thyroid tissue and the production of antibodies against thyroid-specific antigens, such as the thyroid peroxidase antibody (TPOAb) and thyroglobulin antibody (TGAb). TPOAb and TGAb are commonly used in clinical tests; however, handy indicators of the diagnosis and progression of HT are still scarce. Extracellular proteins are glycosylated and are likely to enter body fluids and become readily available and detectable biomarkers. Our research aimed to discover extracellular biomarkers and potential treatment targets associated with HT through integrated bioinformatics analysis and clinical sample validations. A total of 19 extracellular protein-differentially expressed genes (EP-DEGs) were screened by the GSE138198 dataset from the Gene Expression Omnibus (GEO) database and protein annotation databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to analyze the function and pathway of EP-DEGs. STRING, Cytoscape, MCODE, and Cytohubba were used to construct a protein-protein interaction (PPI) network and screen key EP-DEGs. Six key EP-DEGs (CCL5, GZMK, CXCL9, CXCL10, CXCL11, and CXCL13) were further validated in the GSE29315 dataset and the diagnostic curves were evaluated, which all showed high diagnostic accuracy (AUC > 0.95) for HT. Immune profiling revealed the correlation of the six key EP-DEGs and the pivotal immune cells in HT, such as CD8+ T cells, dendritic cells, and Th2 cells. Further, we also confirmed the key EP-DEGs in clinical thyroid samples. Our study may provide bioinformatics and clinical evidence for revealing the pathogenesis of HT and improving the potential diagnosis biomarkers and therapeutic strategies for HT.