Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease

Abstract

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

Figure 1. Genome wide association studies meta-analyses: Loci associated with TPOAb-positivity (a–c) and TPOAb levels (d–f) on a genome-wide level of significance.

Regional association plots of the genome-wide significant loci associated with TPOAb positivity (a–c) and TPOAb levels (d–f). The y-axis on the left indicates the – log₁₀ P value for the association with TPOAb –positivity (a–c) or TPOAb levels (d–f). SNPs are plotted on the x-axis according to their chromosomal position against the association with the phenotype on the y-axis. The most significant stage 1 SNP is indicated in purple. The combined stage 1 and 2 result of this SNP is indicated in yellow. The SNPs surrounding the most significant SNP are color-coded to reflect their LD with this SNP. Symbols reflect functional genomic annotation, as indicated in the legend. The blue y-axes on the right of each plot indicate the estimated recombination rates (based on HapMap Phase II); the bottom of each panel shows the respective annotated genes at the locus and their transcriptional direction. Mb, megabases.