Abstract
BACKGROUND: Autoimmune thyroid disease (AITD) is induced by various factors, including inheritability, which regulates gene expression. Multiple loci correlated with AITD have been discovered utilizing genome-wide association studies (GWASs). Nevertheless, demonstrating the biological relevance and function of these genetic loci is difficult. METHODS: The FUSION software was utilized to define genes that were expressed differentially in AITD using a transcriptome-wide association study (TWAS) method in accordance with GWAS summary statistics from the largest genome-wide association study of 755,406 AITD individuals (30,234 cases and 725,172 controls) and levels of gene expression from two tissue datasets (blood and thyroid). Further analyses were performed such as colocalization, conditional, and fine-mapping analyses to extensively characterize the identified associations, using functional mapping and annotation (FUMA) to conduct functional annotation of the summary statistics of 23329 significant risk SNPs (P < 5 × 10(-8)) recognized by GWAS, together with summary-data-based mendelian randomization (SMR) for identifying functionally related genes at the loci in GWAS. RESULTS: There were 330 genes with transcriptome-wide significant differences between cases and controls, and the majority of these genes were new. 9 of the 94 unique significant genes had strong, colocalized, and potentially causal correlations with AITD. Such strong associations included CD247, TPO, KIAA1524, PDE8B, BACH2, FYN, FOXK1, NKX2-3, and SPATA13. Subsequently, applying the FUMA approach, novel putative AITD susceptibility genes and involved gene sets were detected. Furthermore, we detected 95 probes that showed strong pleiotropic association with AITD through SMR analysis, such as CYP21A2, TPO, BRD7, and FCRL3. Lastly, we selected 26 genes by integrating the result of TWAS, FUMA, and SMR analysis. A phenome-wide association study (pheWAS) was then carried out to determine the risk of other related or co-morbid phenotypes for AITD-related genes. CONCLUSIONS: The current work provides further insight into widespread changes in AITD at the transcriptomic level, as well as characterized the genetic component of gene expression in AITD by validating identified genes, establishing new correlations, and uncovering novel susceptibility genes. Our findings indicate that the genetic component of gene expression plays a significant part in AITD.