Identification of Novel Genetic Loci Associated with Thyroid Peroxidase Antibodies and Clinical Thyroid Disease: e1004123

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.

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.

Author Summary

Individuals with thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune thyroid diseases (AITD), which are common in the general population and associated with increased cardiovascular, metabolic and psychiatric morbidity and mortality. As the causative genes of TPOAbs and AITD remain largely unknown, we performed a genome-wide scan for TPOAbs in 18,297 individuals, with replication in 8,990 individuals. Significant associations were detected with variants at TPO, ATXN2, BACH2, MAGI3, and KALRN. Individuals carrying multiple risk variants also had a higher risk of increased thyroid-stimulating hormone levels (including subclinical and overt hypothyroidism), and a decreased risk of goiter. The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, and the MAGI3 variant was also associated with an increased risk of hypothyroidism. This first genome-wide scan 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. These results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which individuals are particularly at risk of developing clinical thyroid dysfunction.

Citation: Medici M, Porcu E, Pistis G, Teumer A, Brown SJ, et al. (2014) Identification of Novel Genetic Loci Associated with Thyroid Peroxidase Antibodies and Clinical Thyroid Disease. PLoS Genet 10(2): e1004123. doi:10.1371/journal.pgen.1004123

Editor: Chris Cotsapas, Yale School of Medicine, United States of America

Received: August 22, 2013; Accepted: December 3, 2013; Published: February 27, 2014

This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Funding: The Asklepios Study was supported by a Fonds voor Wetenschappelijk Onderzoek-Vlaanderen FWO research grant G.0427.03 and G.0838.10N (Asklepios Study). The 1994-5 Busselton Health Survey was funded by Healthway, Western Australia. The Busselton Health Studies are supported by the National Health and Medical Research Council of Australia and the Great Wine Estates Auctions. The CHS research reported in this article was supported by NHLBI contracts HHSN268201200036C, N01HC85239, N01HC55222, N01HC85079, N01HC85080, N01HC85081, N01HC85082, N01HC85083, N01HC85086; and NHLBI grants HL080295, HL087652, HL105756 with additional contribution from the National Institute of Neurological Disorders and Stroke (NINDS). Additional support was provided through AG023629 from the National Institute on Aging (NIA). DNA handling and genotyping at Cedars-Sinai Medical Center was supported in part by the National Center for Research Resources, grant UL1RR033176, and is now at the National Center for Advancing Translational Sciences, CTSI grant UL1TR000124; in addition to the National Institute of Diabetes and Digestive and Kidney Disease grant DK063491 to the Southern California Diabetes Endocrinology Research Center. Additional funding was provided by the Cedars-Sinai Board of Governors' Chair in Medical Genetics (JIR). The CARLA Study was founded by a grant from the Deutsche Forschungsgemeinschaft as part of the Collaborative Research Center 598 "Heart failure in the elderly - cellular mechanisms and therapy" at the Medical Faculty of the Martin-Luther-University Halle-Wittenberg, by a grant of the Wilhelm-Roux Programme of the Martin-Luther-University Halle-Wittenberg; by the Ministry of Education and Cultural Affairs of Saxony-Anhalt, and by the Federal Employment Office. The Exeter Family Study of Childhood Health (EFSOCH) was supported by South West NHS Research and Development, Exeter NHS Research and Development, the Darlington Trust, and the Peninsula NIHR Clinical Research Facility at the University of Exeter. Genotyping of EFSOCH DNA samples was supported by the Endocrine Research Fund. ATH and BMS are employed as core members of the Peninsula NIHR Clinical Research Facility. RMF is funded by a Sir Henry Wellcome Postdoctoral Fellowship (Wellcome Trust grant: 085541/Z/08/Z). The Health2006 Study is funded by grants from The Velux Foundation; The Danish Medical Research Council, Danish Agency for Science, Technology and Innovation; The Aase and Ejner Danielsens Foundation; ALK-Abelló A/S (Hørsholm, Denmark), Timber Merchant Vilhelm Bangs Foundation, MEKOS Laboratories (Denmark), the Health Insurance Foundation, and Research Centre for Prevention and Health, the Capital Region of Denmark. Helsinki Birth Cohort Study has been supported by grants from the Academy of Finland, the Finnish Diabetes Research Society, Finnish Society for Cardiovascular Research, Folkhälsan Research Foundation, Novo Nordisk Foundation, Finska Läkaresällskapet, Signe and Ane Gyllenberg Foundation, University of Helsinki, European Science Foundation (EUROSTRESS), Ministry of Education, Ahokas Foundation, Emil Aaltonen Foundation, Juho Vainio Foundation, and Wellcome Trust (grant number WT089062). This work was supported by KORA, which is a research platform initiated and financed by the Helmholtz Center Munich, German Research Center for Environmental Health, by the German Federal Ministry of Education and Research and by the State of Bavaria. The work of KORA is supported by the German Federal Ministry of Education and Research (BMBF) in the context of the German National Genome Research Network (NGFN-2 and NGFN-plus). The present research was supported within the Munich Center of Health Sciences (MC Health) as part of LMUinnovativ. Thyroid examinations in KORA-F4 were supported by Sanofi-Aventis in the framework of the Papillon Initiative. Collection and genotyping of the NBS samples was funded in part by the European Commission (POLYGENE: LSHC-CT-2005-018827) and a research investment grant of the Radboud University Nijmegen Medical Centre. This work was sponsored by the National Computing Facilities Foundation (NCF) for the use of supercomputer facilities, with financial support from the NWO. The Thyroid Cancer Program (P.I. Matthew Ringel) at the Ohio State University is supported by grants P30 CA16058 and P01 CA124570 from the National Cancer Institute, USA. The generation and management of GWAS genotype data for the Rotterdam Study is supported by the Netherlands Organisation of Scientific Research NWO Investments (no. 175.010.2005.011, 911-03-012). This study is funded by the Research Institute for Diseases in the Elderly (014-93-015; RIDE2), the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO) project no. 050-060-810. The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. The SardiNIA study is supported by the Intramural Research Program of the National Institute on Aging (NIA), National Institutes of Health (NIH). The SardiNIA ("Progenia") team was supported by Contract NO1-AG-1-2109 from the NIA; the efforts of SS were supported in part by contract 263-MA-410953 from the NIA to the University of Michigan and by research grant HG002651. SHIP is part of the Community Medicine Research net of the University of Greifswald, Germany, which is funded by the Federal Ministry of Education and Research (grants no. 01ZZ9603, 01ZZ0103, and 01ZZ0403), the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg-West Pomerania, and the network 'Greifswald Approach to Individualized Medicine (GANI_MED)' funded by the Federal Ministry of Education and Research (grant 03IS2061A). Genome-wide data have been supported by the Federal Ministry of Education and Research (grant no. 03ZIK012) and a joint grant from Siemens Healthcare, Erlangen, Germany and the Federal State of Mecklenburg-West Pomerania. Data analyses were further supported by the German Research Foundation (DFG Vo 955/10-2; SPP 1629: THYROID TRANS ACT WA 1328/5-1) and the Federal Ministry of Nutrition, Agriculture and Consumer's Safety (BMELV 07 HS 003). SHIP-Trend is part of the Community Medicine Research net of the University of Greifswald, Germany, which is funded by the Federal Ministry of Education and Research (grants no. 01ZZ9603, 01ZZ0103, and 01ZZ0403), the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg - West Pomerania. Thyroid-related examinations have been funded by the Federal Ministry of Nutrition, Agriculture and Consumer's Safety (BMELV 07 HS 003) and the German Research Foundation (DFG Vo 955/10-1; SPP 1629: THYROID TRANS ACT WA 1328/5-1). Genome-wide data have been supported by the Federal Ministry of Education and Research (grant no. 03ZIK012) and a joint grant from Siemens Healthcare, Erlangen, Germany and the Federal State of Mecklenburg, West Pomerania. Whole-body MR imaging was supported by a joint grant from Siemens Healthcare, Erlangen, Germany and the Federal State of Mecklenburg West Pomerania. TwinsUK received funding from the Wellcome Trust; the Chronic Disease Research Foundation; the European Community's Seventh Framework Program grant agreement (FP7/2007-2013); ENGAGE project grant agreement (HEALTH-F4-2007-201413); the Department of Health via the National Institute for Health Research (NIHR) Comprehensive Biomedical Research Centre award to Guy's & St Thomas' NHS Foundation Trust in partnership with King's College London; the Canadian Institutes of Health Research, Canadian Foundation for Innovation, Fonds de la Recherche en Santé Québec, Ministère du Développement Économique, de l'Innovation et de l'Exportation Québec and the Lady Davis Institute of the Jewish General Hospital; the Australian National Health and Medical Research Council (Project Grants 1010494, 1031422) and the Sir Charles Gairdner Hospital Research Fund. Val Borbera was supported by funds from Compagnia di San Paolo, Torino, Italy; Fondazione Cariplo, Italy and Ministry of Health, Ricerca Finalizzata 2008. The UK Graves' disease cohort was funded by the Wellcome Trust grant 068181. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: I have read the journal's policy and have the following conflicts: Dr. Bruce M Psaty reported serving on a DSMB for a clinical trial of a device funded by the manufacturer (Zoll LifeCor) and on the Yale Open Data Access Project funded by Medtronic. All other authors have declared that no competing interests exist.

¶ SS, SN and RPP also contributed equally to this work.

Introduction

Autoimmune thyroid disease (AITD), including Hashimoto's thyroiditis and Graves' disease, is one of the most common autoimmune diseases, affecting 2-5% of the general population [1], [2], [3]. Thyroid dysfunction has been associated with osteoporosis, depression, atrial fibrillation, heart failure, metabolic syndrome, and mortality [4], [5], [6], [7], [8], [9], [10], [11]. High serum antibodies against the enzyme thyroid peroxidase (TPO), which is located in the thyroid and plays a key role in thyroid hormone synthesis, are present in 90% of patients with Hashimoto's thyroiditis [12], [13], the most frequent cause of hypothyroidism and goiter. Although TPO antibodies (TPOAbs) are a useful clinical marker for the detection of early AITD, it remains controversial if these antibodies play a causative role in the pathogenesis of Hashimoto's thyroiditis [14], [15], [16].

Interestingly, TPOAb-positive persons also have an increased risk of developing autoimmune hyperthyroidism (Graves' disease) [17], [18], which is caused by stimulating antibodies against the thyroid stimulating hormone (TSH) receptor [19]. Numerous studies have shown that Graves' hyperthyroidism and Hashimoto's thyroiditis show co-inheritance [17], [20], [21]. Finally, thyroid autoimmunity is the most common autoimmune disorder in women of childbearing age, and TPOAb-positive women have an increased risk of developing pregnancy complications such as miscarriage and pre-term delivery [17], [18], [22], [23], [24], [25], [26].

The prevalence of TPOAb-positivity in the general population ranges from 5-24%, but it is currently unknown why these people develop TPOAbs, nor is it known why not all individuals with thyroid autoimmunity develop clinical thyroid disease [27], [28]. It is estimated that around 70% of the susceptibility to develop thyroid autoantibodies is due to genetic factors [29]. In this context it is remarkable to note that little is known about the genetic factors that determine TPOAb-positivity and the risk of AITD.

We therefore performed a genome wide association study (GWAS) meta-analysis for TPOAbs in the general population in 18,297 individuals from 11 populations. Newly identified genetic variants were studied in relation to subclinical and overt hypo- and hyperthyroidism, goiter, thyroid autoimmunity during pregnancy and thyroid cancer risk.

Results

Characteristics of the studied populations are shown in Table 1 and the Supplementary Material S1. Heritability estimates in the family-based cohorts SardiNIA, TwinsUK and Val Borbera were, respectively, 0.65, 0.66, and 0.54 for TPOAb-positivity, and 0.43, 0.66, and 0.30 for TPOAb levels.

Table 1. Population characteristics and serum TPOAb, TSH, and FT4 level measurements specifications.

Loci associated with TPOAb-positivity and TPOAb levels

See Table 1 and Supplementary Figure S1 for TPOAb measurements and Supplementary Table S1 for genotyping procedures. In most autoimmune diseases, both the presence and the level of autoantibodies are relevant for the disease onset [18], [30], [31]. Furthermore, different pathophysiological processes may be involved in the initiation and severity of the autoimmune response. We therefore performed a GWAS on TPOAb-positivity (including 1769 TPOAb-positives and 16,528 TPOAb-negatives), as well as a GWAS on continuous TPOAb levels (including 12,353 individuals) in stage 1. See Supplementary Figures S2 and S3 for QQ (quantile-quantile) and Manhattan plots.

In stage 2, we followed-up 20 stage 1 SNPs (P<5×10-6; 13 TPOAb-positivity and 10 TPOAb level SNPs, with 3 SNPs overlapping) in 5 populations, including up to 8,990 individuals for TPOAb-positivity (922 TPOAb-positives and 8068 TPOAb-negatives) and 8,159 individuals for TPOAb level analyses (see Supplementary Material S1). Results of the combined stage 1 and 2 meta-analyses, including heterogeneity analyses, are shown in Supplementary Tables S2 and S3. Regional association plots are shown in Supplementary Figures S4 and S5. In the combined stage 1 and 2 meta-analyses GWAS significant associations (P<5×10-8) were observed near TPO (Chr 2p25; rs11675434), at ATXN2 (Chr 12q24.1; rs653178), and BACH2 (Chr 6q15; rs10944479) for TPOAb-positivity, and near TPO (rs11675434), at MAGI3 (Chr 6q15; rs1230666), and KALRN (Chr 3q21; rs2010099) for TPOAb levels (Table 2 and Figure 1). The TPOAb level meta-analysis P-values for the 3 GWAS significant TPOAb-positivity loci were: TPO-rs11675434: P = 7.4×10-13, ATXN2-rs653178: P = 1.3×10-7, and BACH2-rs10944479: P = 2.0×10-4.

[Figure omitted, see PDF]

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 - log10 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.

Table 2. Newly identified loci associated with TPOAb-positivity and/or serum TPOAb levels reaching genome wide significance.

As the 3 GWAS significant loci for TPOAb levels also showed associations with TPOAb-positivity (TPO-rs11675434: OR, 1.21 [95% CI, 1.15-1.28)], P = 1.5×10-16; MAGI3-rs1230666: OR, 1.23 [95% CI, 1.14-1.33], P = 1.5×10-6; KALRN-rs2010099: OR, 1.24 [95% CI, 1.12-1.37], P = 7.4×10-5), we subsequently studied the (combined) effects of these 5 SNPs on clinical thyroid disease. Genetic risk scores were calculated as described in the Supplementary Material. The variance explained by these 5 SNPs was 3.1% for TPOAb-positivity and 3.2% for TPOAb levels. Subjects with a high genetic risk score had a 2.2 times increased risk of TPOAb-positivity compared to subjects with a low genetic risk score (P = 8.1×10-8) (Table 3).

Table 3. Genetic risk score and the risk of TPOAb-positivity.

Table S4 shows the stage 1 TPOAb-positivity and TPOAb level meta-analyses results for GWAS significant SNPs reported in previous GWAS on thyroid related phenotypes.

Associations with hypo- and hyperthyroidism

The associations between the 5 GWAS significant SNPs and the risk of abnormal thyroid function tests are shown in Table 4. MAGI3- rs1230666 was associated with an increased risk of overt hypothyroidism and increased TSH levels below the Bonferroni threshold (i.e., P = 0.05/5 = 0.01). Borderline significant signals were observed at BACH2- rs10944479 with a higher risk of increased TSH levels as well as overt hyperthyroidism (P = 0.011 and P = 0.012), and at the KALRN-rs2010099 SNP with a lower risk of decreased TSH levels (P = 0.010).

Table 4. Newly identified TPOAb associated loci and the risk of thyroid disease in stage 1 and 2 populations.

Furthermore, a higher genetic risk score was associated with a higher risk of increased TSH levels (Supplementary Table S5). No effects of the genetic risk score on the risk of overt hypothyroidism, hyperthyroidism or decreased TSH levels were observed.

Associations with goiter

Individuals with a high genetic risk score had a 30.4% risk of sonographically-proven goiter, compared to 35.2% in subjects with a low score (P = 6.5×10-4) (Table 5). None of the individual SNPs was significantly associated with goiter risk.

Table 5. Newly identified TPOAb associated loci, genetic risk scores and the risk of goiter.

Thyroid autoimmunity during pregnancy

As autoimmunity significantly changes during pregnancy [25], we additionally studied these effects in an independent pregnant population. Pregnant women with a high genetic risk score had a 2.4 times increased risk of TPOAb-positivity compared to women with a low score (10.3% vs 4.8%, P = 0.03). These women did not have a higher risk of increased TSH levels. However, a borderline significant signal with a lower risk of increased TSH levels was observed at ATXN2- rs653178 (OR, 0.54 [95% CI, 0.34-0.87], P = 0.012).

Associations with thyroid disease in independent populations

a) Graves' disease.

As MAGI3- rs1230666 and BACH2- rs10944479 showed promising associations (i.e., P≤0.05) with hyperthyroidism in our meta-analyses, we tested these SNPs in an independent population of 2478 patients with Graves' disease and 2682 controls (see Supplementary Material for further details). Both were associated with an increased risk of Graves' disease (MAGI3- rs1230666: OR, 1.37 [95% CI, 1.22-1.54]; P = 1.2×10-7; BACH2- rs10944479: OR, 1.25 [1.12-1.39]; P = 6.2×10-5).

b) Thyroid cancer.

Supplementary Table S6 shows the associations of the 5 GWAS significant SNPs with thyroid cancer. No statistically significant associations were detected, but a borderline significant signal with an increased risk of thyroid cancer was observed at ATXN2- rs653178 (OR, 1.32 [95% CI, 1.02-1.70], P = 0.03).

Pathway analyses

Ingenuity Pathway Analyses (IPA; Ingenuity Systems, Ca, USA) and GRAIL analyses [32] were performed to identify potential pathways involved in AITD, the results of which are shown in Supplementary Tables S7 and S8, and Figure S6. The identified top pathways involved cell death, survival, movement, and OX40 signalling.

Discussion

This is the first GWAS meta-analysis investigating the genetics of TPOAbs in the normal population in up to 18,297 individuals from 11 populations with replication in up to 8,990 individuals from 5 populations. We identified 5 GWAS significant loci associated with TPOAb-positivity and/or levels.

The most significant hit for both TPOAb-positivity and TPOAb levels was located near the TPO gene itself. TPO is a membrane-bound protein located on the apical membranes of the thyroid follicular cell, catalyzing key reactions in thyroid hormone synthesis [33]. Mutations in TPO have been found in patients with congenital hypothyroidism [34], [35]. Although TPOAbs are valid clinical biomarkers of AITD, they are generally considered to be secondary to the thyroid damage inflicted by T-cells.

The FOXE1 gene has been previously associated with hypothyroidism [36], [37] and is known to regulate transcription of TPO [38]. In this context it is interesting to note that we did not find any associations of the variant near TPO with hypothyroidism. Most genes that have been associated with AITD (predominantly Graves' disease) by candidate gene and GWAS studies so far are located in the HLA class I and II regions, or in genes involved in T-cell (i.e., CTLA-4, PTPN22) or other autoimmune responses [28], [39]. Until now, the TPO gene itself had not been associated with AITD, except in one recent candidate gene analysis in a small cohort (n = 188) without replication [40]. A variant near TPO (rs11694732), which is in LD with rs11675434 (r2 = 0.97 in HapMap2), has previously been associated with TSH levels by Gudmundsson et al [41]. However, various other GWAS on serum TSH and FT4 levels have not found any significant associations in or near this locus, including a recent similar sized GWAS by Porcu et al [42].

Three of the other four loci identified here are located in or are in linkage disequilibrium (LD) with genes previously associated with other autoimmune diseases. Rs1230666 is located in intron 9 of MAGI3, encoding a protein that modulates activity of AKT/PKB. AKT/PKB is expressed in the thyroid and regulates apoptosis [43], which seems to play an important role in the development of AITD [44], [45]. In addition, rs1230666 is in LD with rs2476601 (r2 = 0.70 in HapMap2), a variant causing a R620W substitution in PTPN22. PTPN22 is a lymphoid-specific intracellular phosphatase involved in the T-cell receptor signaling pathway. Variations in PTPN22, and specifically R620W, are associated with various autoimmune disorders including type 1 diabetes, rheumatoid arthritis, systemic lupus erythematosus and Graves' disease [46], [47], [48], [49]. The associations of the MAGI3 locus with TPOAb-positivity and Graves' disease may therefore also be explained by linkage with disease-associated variants in PTPN22 [50]. Of note, the association signal at rs2476601 is one order weaker than that of the top variant rs1230666.

The BACH2 locus has been implicated in the susceptibility to several autoimmune diseases, including celiac disease, type 1 diabetes, vitiligo, Crohn's disease, and multiple sclerosis [46], [51], [52], [53], [54]. A recent candidate gene analysis associated the BACH2 locus with an increased risk of AITD, including Hashimoto's thyroiditis and Graves' disease [55]. However, the associations were not significant when Hashimoto's thyroiditis and Graves' disease were studied separately. BACH2 is specifically expressed in early stages of B-cell differentiation and represses different immunoglobulin genes [56]. Interestingly, BACH2 can bind to the co-repressor SMRT (silencing mediator of retinoid and thyroid receptor), which may suggest a more direct effect on thyroid hormone secretion and action as well.

Polymorphisms in ATXN2 have been associated with multiple neurodegenerative diseases, including spinocerebellar ataxia and Parkinson's disease [57], [58], [59]. Different epidemiological studies have associated thyroid dysfunction with cerebellar ataxia [60], [61]. Furthermore, the identified SNP in ATXN2 has been previously associated with renal function, serum urate levels and blood pressure [62], [63], [64]. However, this SNP is in high LD with rs3184504 (r2 = 0.873), a variant causing a Trp262Arg substitution of SH2B adaptor protein 3 (SH2B3). SH2B3 encodes the adaptor protein LNK, a key negative regulator of cytokine signaling playing a critical role in hematopoiesis. This variant is associated with susceptibility to several autoimmune diseases, including celiac disease, type 1 diabetes, vitiligo, and rheumatoid arthritis [46], [51], [53], [65], suggesting more relevance for TPOAb levels than ATXN2. This is supported by a recent study which showed that variants in LD with SH2B3, BACH2, and PTPN22 are associated with TPOAb levels in patients with type 1 diabetes [66].

Whereas the above four loci are located in genes involved in the immune response or the autoantigen, the KALRN (Kalirin) gene encodes a multi-domain guanine nucleotide exchange factor for GTP-binding proteins of the Rho family. The relation of KALRN with levels of TPOAbs is unclear. This gene has recently been found to be associated with megakaryopoiesis and platelet formation [67], which may suggest a function in the immune system [68]. We furthermore performed pathway analyses on the stage 1 TPOAb-positivity and TPOAb level lead SNPs, and identified the cell death, survival and movement pathway as an important pathway for TPOAbs. This finding is supported by previous studies, which show an important role for apoptosis in the development of AITD [44], [45]. Another top pathway involved was the OX40 signalling pathway, and it is of interest to note that OX40 is a T-cell activator promoting the survival of CD4+ T-cells at sites of inflammation [69].

Our results have potential clinical relevance for several reasons. Genetic risk scores based on these novel common (risk allele frequencies: 9-40%) TPOAb-associated SNPs enabled us to identify a large subgroup in the general population with a two-fold increased risk of TPOAb-positivity (10.4% vs 5.4%). These individuals also have a higher risk of increased TSH levels and a lower risk of goiter, suggesting an advanced stage of destruction of the thyroid due to autoimmune processes. Furthermore, pregnant women with high genetic risk scores had a 2.4 times increased risk of TPOAb-positivity during pregnancy. In this context it is interesting to note that TPOAb-positive pregnant women have an increased risk of miscarriages and preterm births independent of thyroid function [70].

Associations with thyroid disease were also found on an individual SNP level. The MAGI3 SNP was associated with a substantially increased risk of hypothyroidism, and the BACH2 SNP showed a borderline significant association (P = 0.011) with a higher risk of increased TSH levels, which includes subjects with subclinical and overt hypothyroidism. Furthermore, both loci were significantly associated with an increased risk of Graves' hyperthyroidism in an independent population. To predict which patients with first or second degree relatives with documented Hashimoto's or Graves' disease will develop clinical thyroid disease, a clinical algorithm has been developed (i.e., the THEA score) [18]. Future studies should analyze if these genetic markers increase the sensitivity of the THEA score. Graves' hyperthyroidism and Hashimoto's thyroiditis co-segregate in families and subjects with TPOAbs have an increased risk of both diseases [17], [18], [20], [21], [22], [26]. The current study provides insight into this phenomenon by showing that specific loci associated with TPOAbs and (subclinical) hypothyroidism, i.e. MAGI3 and BACH2, are also associated with Graves' hyperthyroidism in an independent case-control study.

The prevalence of TPOAb-positivity in the general population is high (5-24%), but it is currently unknown why part of the individuals with thyroid autoimmunity develop clinical thyroid disease whereas others do not [27], [28]. In this context it is interesting to note that the TPOAb-associated SNPs located in TPO and ATXN2 were not associated with clinical thyroid disease. This suggests that the TPOAbs in these individuals may be of less clinical relevance, providing insight into why TPOAb-positive individuals do or do not eventually develop clinical thyroid disease.

Our study has some limitations. The validity of the results is restricted to individuals from populations of European ancestry. Future GWASs in populations from non-European descent will be required to determine to which extent our results can be generalized to other ethnic groups. Secondly, we did not perform conditional analyses to further identify secondary association signals within the identified loci, nor did we perform functional studies for the identified variants. Further research is therefore needed to unravel the exact biological mechanism behind the observed associations. The fact that various TPOAb assays were used across the participating cohorts could lead to bias. We therefore used TPOAb-positivity cut-off values as provided by the respective assay manufacturer, instead of using one fixed cut-off value. This is also of clinical importance as in clinical practice most institutions rely on the TPOAb-positivity cut-off as provided by the assay manufacturer. Furthermore, we did not detect heterogeneity in our results, supporting the fact that results obtained with different assays can be combined across cohorts using the z-score based meta-analysis. Finally, as AITD coincides with other autoimmune diseases, our results could be driven by indirect associations with other autoimmune diseases. However, AITD is the most common autoimmune disease in the general population. We furthermore show that carriage of multiple risk alleles is associated with an increased risk of thyroid dysfunction, which underlines the clinical importance of our findings.

In conclusion, this first GWAS for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. Furthermore, we show that carriage of multiple risk variants is not only associated with a substantial increased risk of TPOAb-positivity, but also with a higher risk of increased TSH levels (including subclinical and overt hypothyroidism) and a lower risk of goiter. These genetic markers not only help to identify large groups in the general population with an increased risk of TPOAb-positivity, but may also predict which TPOAb-positive persons are particularly at risk of developing clinical thyroid disease.

Materials and Methods

Study cohorts

For the TPOAb GWAS stage 1 and 2 analyses, and the hypothyroidism, hyperthyroidism and goiter analyses, individuals were recruited from 16 independent community-based and family studies. For the Graves' disease analyses, cases were recruited from the United Kingdom Graves' disease cohort and controls from the British 1958 Birth Cohort. Thyroid cancer cases and controls were recruited from the Nijmegen and Ohio thyroid cancer cohorts. A detailed description of the original cohorts contributing samples is provided in Table 1 and in the Supplementary Material. All participants provided written informed consent and protocols were approved by the institutional review boards or research ethics committees at the respective institutions, and conducted according to the Declaration of Helsinki.

Phenotype definitions

Serum TPOAb levels were determined with a range of assays. TPOAb-positives were defined as subjects with TPOAb levels above the assay-specific TPOAb-positivity cut-off, as defined by the manufacturer (Table 1). Serum TSH and free thyroxine (FT4) levels were determined using a range of assays (Table 1). Assay-specific TSH and FT4 reference ranges were used, as provided by the manufacturer (Table 1). Overt hypothyroidism was defined as a high TSH (i.e., a TSH level above the TSH reference range) and a low FT4. Increased TSH was defined as a high TSH, including persons with overt hypothyroidism or subclinical hypothyroidism (i.e., high TSH with a normal FT4). Overt hyperthyroidism was defined as a low TSH and a high FT4. Decreased TSH was defined as a low TSH, including persons with subclinical or overt hyperthyroidism.

The diagnosis of goiter is described in the Supplementary Material, and the diagnosis of Graves' disease and thyroid cancer in the respective cohorts have been described previously [41].

Genotyping

Samples were genotyped with a range of GWAS genotyping arrays (Supplementary Table S1). Sample and SNP quality control procedures were undertaken within each study. For each GWAS, over 2.5 million SNPs were imputed using CEU samples from Phase 2 of the International HapMap project (www.hapmap.org). Genotyping procedures in the stage 2, Graves' disease and thyroid cancer populations are described in the Supplementary Material.

Association analyses

The heritabilities of TPOAb-positivity and serum TPOAb levels were estimated, as described in the Supplementary Material.

In stage 1, we performed a GWAS on TPOAb-positivity as well as a GWAS on continuous TPOAb levels. Persons taking thyroid medication were excluded. Each SNP was tested for association with TPOAb-positivity using logistic regression analyses, adjusting for age and sex. For cohorts with family structure, we approximated the probability of being affected with a linear mixed model adjusting for age and sex. The produced model was used to predict the expected proportion of "risk" (effective) alleles in cases and controls, hence giving the means to estimate odds ratios. Only unrelated individuals were considered for the SardiNIA cohort. For the GWAS of continuous TPOAb levels, samples with a TPOAb level lower than the minimum TPOAb assay detection limit (Table 1) were excluded. TPOAb levels were natural log-transformed, and sex-specific, age adjusted standardized residuals were calculated. Each SNP was tested for association with these TPOAb level residuals using linear regression analyses (additive model), correcting for relatedness in studies with family structure. See Supplementary Table S1 for the software used for these analyses.

Before meta-analysis, SNPs with a minor allele frequency (MAF) <1% or a low imputation quality were excluded (Supplementary Material), after which the results of each GWAS were combined in a population size weighted z-score based meta-analysis using METAL [71]. Genomic control was applied to individual studies if λ>1.0.

In stage 2, we followed-up stage 1 GWAS significant SNPs, as well as promising SNPs not reaching GWAS significance, in an attempt to reach GWAS significant associations by increasing sample size (Supplementary Material). Results from stage 1 and 2 were combined in a population size weighted z-score based meta-analysis using METAL [71]. A z-score based meta-analysis was used to reduce bias that might be induced by different assays. As this method does not provide betas, and we wanted to provide a rough estimate of the actual effect sizes for convenience, we calculated betas using the fixed effects (inverse variance based) meta-analysis method. Heterogeneity was tested, applying bonferroni based P-value thresholds of P = 0.004 for the TPOAb-positivity analyses and P = 0.005 for the TPOAb level analyses.

All studies assessed and, if present, corrected for population stratification using principal-component analysis (PCA) and/or multidimensional-scaling (MDS), with the exception of SardiNIA and ValBorbera where the high isolation substantiates a lack of stratification (Table S1) [72], [73]. Lambda values were all ~1, indicating that population stratification was overall properly accounted for (Table S1). To fully remove residual effects, we applied genomic correction to studies were lambda was >1. The final meta-analyses reported a lambda of 1.01 for both the TPOAb-positivity and the TPOAb level GWAS, thus no genomic correction was applied.

The variances explained by the GWAS significant SNPs were calculated. We subsequently studied the individual as well as the combined effects of the GWAS significant SNPs on the risk of clinical thyroid disease, as specified in the Supplementary Material. In short, to study combined effects, a genetic risk score was calculated for every person as the weighted sum of TPOAb risk alleles. The associations between the individual SNPs, genetic risk scores and the risk of abnormal thyroid function tests were studied using logistic regression analyses. Logistic regression analyses were used to study the associations with goiter, Graves' disease and thyroid cancer (Supplementary Material). The results of each study were combined in a population size weighted z-score based meta-analysis using METAL [71].

Various bioinformatic tools were searched for evidence for functional relevance of the GWAS significant SNPs and pathway analyses were performed on the Stage 1 lead SNPs (see Supplementary Material).

Supporting Information

Figure S1.

TPOAb level distributions in persons with detectable TPOAb levels in stage 1 and 2 populations.

(PPTX)

Figure S2.

Quantile-quantile (QQ) plots for the TPOAb-positivity and TPOAb level stage 1 meta-analyses.

(PPTX)

Figure S3.

Manhattan plots for stage 1 meta-analyses for TPOAb-positivity (a) and TPOAb levels (b). SNPs are plotted on the x-axis according to their chromosomal position against TPOAb-positivity (a) or TPOAb levels (b) (shown as - log10 P value) on the y-axis. The horizontal grey line indicates the threshold for genome-wide statistical significance (P<5×10-8). Genome-wide significant associations were observed near TPO (Chr 2p25; P = 1.5×10-12), at ATXN2 (Chr 12q24.1; P = 1.6×10-9) and near HCP5 (Chr 6p21.3; P = 4.1×10-8) for TPOAb-positivity, and near TPO (Chr 2p25; P = 5.4×10-13) and at ATXN2 (Chr 12q24.1; P = 1.1×10-8) for TPOAb levels.

(PPTX)

Figure S4.

Regional association plots of stage 1 lead loci for TPOAb-positivity (panels a-m). The y-axis on the left indicates the - log10 P value for the association with TPOAb -positivity. SNPs are plotted on the x-axis according to their chromosomal position. 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.

(PPTX)

Figure S5.

Regional association plots of stage 1 lead loci for TPOAb levels (panels a-j). The y-axis on the left indicates the - log10 P value for the association with TPOAb levels. SNPs are plotted on the x-axis according to their chromosomal position. 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.

(PPTX)

Figure S6.

GRAIL results for the stage 1 TPOAb-positivity and TPOAb level lead SNPs. GRAIL circle plot of locus connectivity where each locus is plotted in a circle, where significant connections (P<0.05) based on PubMed abstracts are drawn spanning the circle. Analyses were based on the 20 stage 1 TPOAb-positivity and TPOAb level lead SNPs.

(PPTX)

Table S1.

Study sample genotyping, quality control and association analyses for stage 1 populations.

(DOCX)

Table S2.

Associations of stage 1 lead SNPs with TPOAb-positivity in stage 1 and 2.

(DOCX)

Table S3.

Associations of stage 1 lead SNPs with serum TPOAb levels in stage 1 and 2.

(DOCX)

Table S4.

Stage 1 TPOAb-positivity and TPOAb level meta-analyses results for GWAS significant SNPs reported in previous GWAS on thyroid related phenotypes.

(XLSX)

Table S5.

Genetic risk score and the risk of increased TSH levels.

(DOCX)

Table S6.

Newly identified TPOAb associated loci and the risk of thyroid cancer.

(DOCX)

Table S7.

Top IPA associated networks for the Stage 1 TPOAb-positivity and TPOAb level lead SNPs.

(DOCX)

Table S8.

Top IPA associated canonical pathways for the Stage 1 TPOAb-positivity and TPOAb level lead SNPs.

(DOCX)

Text S1.

Supplementary methods.

(DOCX)

Acknowledgments

We thank all study participants, volunteers and study personnel that made this work possible.

The Asklepios study is indebted to Femke van Hoeke, Bianca Leydens, and Caroline van Daele, and the residents and general practitioners of Erpe-Mere and Nieuwerkerken for their help in completing the study.

The Busselton Health Study thanks the Busselton Population Medical Research Foundation for approving the study. We thank Siemens Ltd. Australia and New Zealand Healthcare Sector for donating assay reagents.

The Rotterdam Study thanks Pascal Arp, Mila Jhamai, Marijn Verkerk, Lizbeth Herrera and Marjolein Peters for their help in creating the GWAS database, and Karol Estrada and Maksim V. Struchalin for their support in creation and analysis of imputed data. The authors are grateful to the study participants, the staff from the Rotterdam Study and the participating general practitioners and pharmacists. We would like to thank Karol Estrada, Dr. Fernando Rivadeneira, Dr. Tobias A. Knoch, Anis Abuseiris, Luc V. de Zeeuw, and Rob de Graaf (Erasmus MC Rotterdam, The Netherlands), for their help in creating GRIMP, and BigGRID, MediGRID, and Services@MediGRID/D-Grid for access to their grid computing resources. We would like to thank Symen Ligthart for his help with the IPA and GRIMP pathway analyses.

The SardiNIA study thanks the many individuals who generously participated in this study, Monsignore Piseddu, Bishop of Ogliastra, the mayors and citizens of the Sardinian towns (Lanusei, Ilbono, Arzana, and Elini), and the head of the Public Health Unit ASL4 for their volunteerism and cooperation; the team also thanks the physicians, Marco Orrù, Maria Grazia Pilia, Liana Ferreli, Francesco Loi, Stefano Angius, nurses Paola Loi, Monica Lai and Anna Cau who carried out participant physical exams, and the recruitment personnel Susanna Murino. We thank Francesco Cucca, PI of the SardiNIA study.

The SHIP study is grateful to the contribution of Florian Ernst, Anja Wiechert and Astrid Petersmann in generating the SNP data.

The SHIP-Trend study is grateful to Mario Stanke for the opportunity to use his Server Cluster for the SNP imputation as well as to Holger Prokisch and Thomas Meitinger (Helmholtz Zentrum München) for the genotyping of the SHIP-TREND cohort.

TwinsUK thanks the staff from the Genotyping Facilities at the Wellcome Trust Sanger Institute, UK, for sample preparation, quality control, and genotyping; Le Centre National de Génotypage, France, for genotyping; Duke University, NC, USA, for genotyping; and the Finnish Institute of Molecular Medicine, Finnish Genome Center, University of Helsinki. We thank the volunteer twins who made available their time.

The United Kingdom (UK) Graves' disease cohort would like to thank all principle investigators (Amit Allahabadia, Northern General Hospital; Sheffield, UK, Mary Armitage Royal Bournemouth Hospital, Bournemouth, UK; Krishna V. Chatterjee, University of Cambridge, Addenbrookes Hospital, Cambridge, UK; John H. Lazarus Centre for Endocrine and Diabetes Sciences, Cardiff University, Cardiff, UK; Simon H. Pearce, Institute of Human Genetics, Newcastle University, Newcastle-upon-Tyne, Newcastle, UK and Bijay Viadya, Royal Devon and Exeter Hospital, Exeter, UK), doctors and nurses for recruiting AITD subjects into the AITD National Collection.

Val Borbera thanks the inhabitants of the Val Borbera for participating in the study, the local administrations and the ASL-Novi Ligure for support and Fiammetta Viganò for technical help. We also thank Prof. Clara Camaschella, Prof Federico Caligaris-Cappio and the MDs of the Medicine Dept. of the San Raffaele Hospital for help with clinical data collection.

Author Contributions

Conceived and designed the experiments: MM SJB RAJ RR AA HJG ER JIR HH LC DTi BV TdM TJ JGE BMP AHo DS HW AdlC TMF AL KR LAK AGU JPW KS EWic CMe MdH TJV TDS SGW HV AC DTo SS SN RPP. Performed the experiments: MM EP GP AT LC SJB RAJ RR GLR TSP SHV JL MJS LLNH RMF BMS CG YSA AL TJV SS SN RPP. Analyzed the data: MM EP GP AT SJB RAJ RR GLR TSP SHV JL MJS LLNH RMF SLi BMS DP LC LB CG TC EK BT YET AA MvdB CMa TEG MT NP YSA AdlC RTNM SCLG JMK AL JWAS FR MdH SS RPP. Contributed reagents/materials/analysis tools: MM RR GLR TSP JL MJS LLNH BMS RN MGP CSa UV JBR FCS TIMK WEV ATH JK LC AHa WL GH ML SM NS MC MN CSp AR MH EML ER PJL SLa MV GA EWid AP AD APB DIWP JPB AM TF AJ JH HP EER PF SJF JIR AK DR GLS EB HH JAF BV TdM TJ JGE PCO ARH BMP TI AHo HW AdlC RTNM SCLG HMzS TMF AL FR AGU JPW CMe TJV TDS SGW HV AC DTo RPP. Wrote the paper: MM AT LC TJV SGW AC SS SN RPP.

References

Gough SC (2000) The genetics of Graves' disease. Endocrinol Metab Clin North Am 29: 255-266. doi: 10.1016/s0889-8529(05)70130-4

Simmonds MJ, Gough SC (2004) Unravelling the genetic complexity of autoimmune thyroid disease: HLA, CTLA-4 and beyond. Clin Exp Immunol 136: 1-10. doi: 10.1111/j.1365-2249.2004.02424.x

Tunbridge WM, Evered DC, Hall R, Appleton D, Brewis M, et al. (1977) The spectrum of thyroid disease in a community: the Whickham survey. Clin Endocrinol (Oxf) 7: 481-493. doi: 10.1111/j.1365-2265.1977.tb01340.x

Biondi B (2012) Mechanisms in endocrinology: Heart failure and thyroid dysfunction. Eur J Endocrinol 167: 609-618. doi: 10.1530/eje-12-0627

Collet TH, Gussekloo J, Bauer DC, den Elzen WP, Cappola AR, et al. (2012) Subclinical hyperthyroidism and the risk of coronary heart disease and mortality. Arch Intern Med 172: 799-809. doi: 10.1001/archinternmed.2012.402

Davis JD, Tremont G (2007) Neuropsychiatric aspects of hypothyroidism and treatment reversibility. Minerva Endocrinol 32: 49-65

Gencer B, Collet TH, Virgini V, Bauer DC, Gussekloo J, et al. (2012) Subclinical thyroid dysfunction and the risk of heart failure events: an individual participant data analysis from 6 prospective cohorts. Circulation 126: 1040-1049. doi: 10.1161/circulationaha.112.096024

Nicholls JJ, Brassill MJ, Williams GR, Bassett JH (2012) The skeletal consequences of thyrotoxicosis. J Endocrinol 213: 209-221. doi: 10.1530/joe-12-0059

Rodondi N, den Elzen WP, Bauer DC, Cappola AR, Razvi S, et al. ((201)) Subclinical hypothyroidism and the risk of coronary heart disease and mortality. Jama 304: 1365-1374. doi: 10.1001/jama.2010.1361

Ruhla S, Weickert MO, Arafat AM, Osterhoff M, Isken F, et al. (2010) A high normal TSH is associated with the metabolic syndrome. Clin Endocrinol (Oxf) 72: 696-701. doi: 10.1111/j.1365-2265.2009.03698.x

Selmer C, Olesen JB, Hansen ML, Lindhardsen J, Olsen AM, et al. (2012) The spectrum of thyroid disease and risk of new onset atrial fibrillation: a large population cohort study. Bmj 345: e7895. doi: 10.1136/bmj.e7895

Pearce EN, Farwell AP, Braverman LE (2003) Thyroiditis. N Engl J Med 348: 2646-2655. doi: 10.1056/nejmra021194

Schweizer U, Chiu J, Kohrle J (2008) Peroxides and peroxide-degrading enzymes in the thyroid. Antioxid Redox Signal 10: 1577-1592. doi: 10.1089/ars.2008.2054

Brix TH, Hegedus L, Gardas A, Banga JP, Nielsen CH (2011) Monozygotic twin pairs discordant for Hashimoto's thyroiditis share a high proportion of thyroid peroxidase autoantibodies to the immunodominant region A. Further evidence for genetic transmission of epitopic "fingerprints". Autoimmunity 44: 188-194. doi: 10.3109/08916934.2010.518575

Huber G, Staub JJ, Meier C, Mitrache C, Guglielmetti M, et al. (2002) Prospective study of the spontaneous course of subclinical hypothyroidism: prognostic value of thyrotropin, thyroid reserve, and thyroid antibodies. J Clin Endocrinol Metab 87: 3221-3226. doi: 10.1210/jcem.87.7.8678

Nielsen CH, Brix TH, Leslie RG, Hegedus L (2009) A role for autoantibodies in enhancement of pro-inflammatory cytokine responses to a self-antigen, thyroid peroxidase. Clin Immunol 133: 218-227. doi: 10.1016/j.clim.2009.07.014

Strieder TG, Prummel MF, Tijssen JG, Endert E, Wiersinga WM (2003) Risk factors for and prevalence of thyroid disorders in a cross-sectional study among healthy female relatives of patients with autoimmune thyroid disease. Clin Endocrinol (Oxf) 59: 396-401. doi: 10.1046/j.1365-2265.2003.01862.x

Strieder TG, Tijssen JG, Wenzel BE, Endert E, Wiersinga WM (2008) Prediction of progression to overt hypothyroidism or hyperthyroidism in female relatives of patients with autoimmune thyroid disease using the Thyroid Events Amsterdam (THEA) score. Arch Intern Med 168: 1657-1663. doi: 10.1001/archinte.168.15.1657

Weetman AP (2000) Graves' disease. N Engl J Med 343: 1236-1248. doi: 10.1056/nejm200010263431707

Brix TH, Hegedus L (2011) Twins as a tool for evaluating the influence of genetic susceptibility in thyroid autoimmunity. Ann Endocrinol (Paris) 72: 103-107. doi: 10.1016/j.ando.2011.03.013

Tomer Y, Davies TF (2003) Searching for the autoimmune thyroid disease susceptibility genes: from gene mapping to gene function. Endocr Rev 24: 694-717. doi: 10.1210/er.2002-0030

Kordonouri O, Deiss D, Danne T, Dorow A, Bassir C, et al. (2002) Predictivity of thyroid autoantibodies for the development of thyroid disorders in children and adolescents with Type 1 diabetes. Diabet Med 19: 518-521. doi: 10.1046/j.1464-5491.2002.00699.x

Medici M, de Rijke YB, Peeters RP, Visser W, de Muinck Keizer-Schrama SM, et al. (2012) Maternal early pregnancy and newborn thyroid hormone parameters: the Generation R study. J Clin Endocrinol Metab 97: 646-652. doi: 10.1210/jc.2011-2398

Negro R, Mestman JH (2011) Thyroid disease in pregnancy. Best Pract Res Clin Endocrinol Metab 25: 927-943. doi: 10.1016/j.beem.2011.07.010

Poppe K, Velkeniers B, Glinoer D (2008) The role of thyroid autoimmunity in fertility and pregnancy. Nat Clin Pract Endocrinol Metab 4: 394-405. doi: 10.1038/ncpendmet0846

Vanderpump MP, Tunbridge WM, French JM, Appleton D, Bates D, et al. (1995) The incidence of thyroid disorders in the community: a twenty-year follow-up of the Whickham Survey. Clin Endocrinol (Oxf) 43: 55-68. doi: 10.1111/j.1365-2265.1995.tb01894.x

Hollowell JG, Staehling NW, Flanders WD, Hannon WH, Gunter EW, et al. (2002) Serum TSH, T(4), and thyroid antibodies in the United States population (1988 to 1994): National Health and Nutrition Examination Survey (NHANES III). J Clin Endocrinol Metab 87: 489-499. doi: 10.1210/jcem.87.2.8182

Weetman AP (2011) Diseases associated with thyroid autoimmunity: explanations for the expanding spectrum. Clin Endocrinol (Oxf) 74: 411-418. doi: 10.1111/j.1365-2265.2010.03855.x

Hansen PS, Brix TH, Iachine I, Kyvik KO, Hegedus L (2006) The relative importance of genetic and environmental effects for the early stages of thyroid autoimmunity: a study of healthy Danish twins. Eur J Endocrinol 154: 29-38. doi: 10.1530/eje.1.02060

Linnik MD, Hu JZ, Heilbrunn KR, Strand V, Hurley FL, et al. (2005) Relationship between anti-double-stranded DNA antibodies and exacerbation of renal disease in patients with systemic lupus erythematosus. Arthritis Rheum 52: 1129-1137. doi: 10.1002/art.20980

Nielsen SF, Bojesen SE, Schnohr P, Nordestgaard BG (2012) Elevated rheumatoid factor and long term risk of rheumatoid arthritis: a prospective cohort study. Bmj 345: e5244. doi: 10.1136/bmj.e5244

Raychaudhuri S, Plenge RM, Rossin EJ, Ng AC, Purcell SM, et al. (2009) Identifying relationships among genomic disease regions: predicting genes at pathogenic SNP associations and rare deletions. PLoS Genet 5: e1000534. doi: 10.1371/journal.pgen.1000534

Ruf J, Carayon P (2006) Structural and functional aspects of thyroid peroxidase. Arch Biochem Biophys 445: 269-277. doi: 10.1016/j.abb.2005.06.023

Bakker B, Bikker H, Vulsma T, de Randamie JS, Wiedijk BM, et al. (2000) Two decades of screening for congenital hypothyroidism in The Netherlands: TPO gene mutations in total iodide organification defects (an update). J Clin Endocrinol Metab 85: 3708-3712. doi: 10.1210/jcem.85.10.6878

Bikker H, Baas F, De Vijlder JJ (1997) Molecular analysis of mutated thyroid peroxidase detected in patients with total iodide organification defects. J Clin Endocrinol Metab 82: 649-653. doi: 10.1210/jcem.82.2.3729

Denny JC, Crawford DC, Ritchie MD, Bielinski SJ, Basford MA, et al. (2011) Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies. Am J Hum Genet 89: 529-542. doi: 10.1016/j.ajhg.2011.09.008

Eriksson N, Tung JY, Kiefer AK, Hinds DA, Francke U, et al. (2012) Novel associations for hypothyroidism include known autoimmune risk loci. PLoS One 7: e34442. doi: 10.1371/journal.pone.0034442

Ortiz L, Aza-Blanc P, Zannini M, Cato AC, Santisteban P (1999) The interaction between the forkhead thyroid transcription factor TTF-2 and the constitutive factor CTF/NF-1 is required for efficient hormonal regulation of the thyroperoxidase gene transcription. J Biol Chem 274: 15213-15221. doi: 10.1074/jbc.274.21.15213

Simmonds MJ (2013) GWAS in autoimmune thyroid disease: redefining our understanding of pathogenesis. Nat Rev Endocrinol 9: 277-287. doi: 10.1038/nrendo.2013.56

Faam B, Daneshpour MS, Azizi F, Salehi M, Hedayati M (2012) Association between TPO gene polymorphisms and Anti-TPO level in Tehranian population: TLGS. Gene 498: 116-119. doi: 10.1016/j.gene.2012.01.077

Gudmundsson J, Sulem P, Gudbjartsson DF, Jonasson JG, Masson G, et al. (2012) Discovery of common variants associated with low TSH levels and thyroid cancer risk. Nat Genet 44: 319-322. doi: 10.1038/ng.1046

Porcu E, Medici M, Pistis G, Volpato CB, Wilson SG, et al. (2013) A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function. PLoS Genet 9: e1003266. doi: 10.1371/journal.pgen.1003266

Wu Y, Dowbenko D, Spencer S, Laura R, Lee J, et al. (2000) Interaction of the tumor suppressor PTEN/MMAC with a PDZ domain of MAGI3, a novel membrane-associated guanylate kinase. J Biol Chem 275: 21477-21485. doi: 10.1074/jbc.m909741199

Bossowski A, Czarnocka B, Bardadin K, Stasiak-Barmuta A, Urban M, et al. (2008) Identification of apoptotic proteins in thyroid gland from patients with Graves' disease and Hashimoto's thyroiditis. Autoimmunity 41: 163-173. doi: 10.1080/08916930701727749

Wang SH, Baker JR (2007) The role of apoptosis in thyroid autoimmunity. Thyroid 17: 975-979. doi: 10.1089/thy.2007.0208

Barrett JC, Clayton DG, Concannon P, Akolkar B, Cooper JD, et al. (2009) Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet 41: 703-707. doi: 10.1038/ng.381

Orru V, Tsai SJ, Rueda B, Fiorillo E, Stanford SM, et al. (2009) A loss-of-function variant of PTPN22 is associated with reduced risk of systemic lupus erythematosus. Hum Mol Genet 18: 569-579. doi: 10.1093/hmg/ddn363

Raychaudhuri S, Remmers EF, Lee AT, Hackett R, Guiducci C, et al. (2008) Common variants at CD40 and other loci confer risk of rheumatoid arthritis. Nat Genet 40: 1216-1223. doi: 10.1038/ng.233

Velaga MR, Wilson V, Jennings CE, Owen CJ, Herington S, et al. (2004) The codon 620 tryptophan allele of the lymphoid tyrosine phosphatase (LYP) gene is a major determinant of Graves' disease. J Clin Endocrinol Metab 89: 5862-5865. doi: 10.1210/jc.2004-1108

Steer S, Abkevich V, Gutin A, Cordell HJ, Gendall KL, et al. (2007) Genomic DNA pooling for whole-genome association scans in complex disease: empirical demonstration of efficacy in rheumatoid arthritis. Genes Immun 8: 57-68. doi: 10.1038/sj.gene.6364359

Dubois PC, Trynka G, Franke L, Hunt KA, Romanos J, et al. (2010) Multiple common variants for celiac disease influencing immune gene expression. Nat Genet 42: 295-302. doi: 10.1038/ng.543

Franke A, McGovern DP, Barrett JC, Wang K, Radford-Smith GL, et al. (2010) Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci. Nat Genet 42: 1118-1125. doi: 10.1038/ng.717

Jin Y, Birlea SA, Fain PR, Ferrara TM, Ben S, et al. (2012) Genome-wide association analyses identify 13 new susceptibility loci for generalized vitiligo. Nat Genet 44: 676-680. doi: 10.1038/ng.2272

Sawcer S, Hellenthal G, Pirinen M, Spencer CC, Patsopoulos NA, et al. (2011) Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature 476: 214-219

Cooper JD, Simmonds MJ, Walker NM, Burren O, Brand OJ, et al. (2012) Seven newly identified loci for autoimmune thyroid disease. Hum Mol Genet 21: 5202-5208. doi: 10.1093/hmg/dds357

Muto A, Hoshino H, Madisen L, Yanai N, Obinata M, et al. (1998) Identification of Bach2 as a B-cell-specific partner for small maf proteins that negatively regulate the immunoglobulin heavy chain gene 3' enhancer. Embo J 17: 5734-5743. doi: 10.1093/emboj/17.19.5734

Houlden H, Singleton AB (2012) The genetics and neuropathology of Parkinson's disease. Acta Neuropathol 124: 325-338. doi: 10.1007/s00401-012-1013-5

Liu X, Lu M, Tang L, Zhang N, Chui D, et al. (2013) ATXN2 CAG repeat expansions increase the risk for Chinese patients with amyotrophic lateral sclerosis. Neurobiol Aging 34: 2236.e5-8. doi: 10.1016/j.neurobiolaging.2013.04.009

Magana JJ, Velazquez-Perez L, Cisneros B (2013) Spinocerebellar ataxia type 2: clinical presentation, molecular mechanisms, and therapeutic perspectives. Mol Neurobiol 47: 90-104. doi: 10.1007/s12035-012-8348-8

Bonuccelli U, Nuti A, Monzani F, De Negri F, Muratorio A (1991) Familial occurrence of hypothyroidism and cerebellar ataxia. Funct Neurol 6: 171-175

Tandeter H, Levy A, Gutman G, Shvartzman P (2001) Subclinical thyroid disease in patients with Parkinson's disease. Arch Gerontol Geriatr 33: 295-300. doi: 10.1016/s0167-4943(01)00196-0

Kottgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, et al. (2013) Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet 45: 145-154

Kottgen A, Pattaro C, Boger CA, Fuchsberger C, Olden M, et al. (2010) New loci associated with kidney function and chronic kidney disease. Nat Genet 42: 376-384. doi: 10.1055/s-0030-1266452

Wain LV, Verwoert GC, O'Reilly PF, Shi G, Johnson T, et al. (2011) Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure. Nat Genet 43: 1005-1011

Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, et al. (2010) Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat Genet 42: 508-514. doi: 10.1038/ng.582

Plagnol V, Howson JM, Smyth DJ, Walker N, Hafler JP, et al. (2011) Genome-wide association analysis of autoantibody positivity in type 1 diabetes cases. PLoS Genet 7: e1002216. doi: 10.1371/journal.pgen.1002216

Gieger C, Radhakrishnan A, Cvejic A, Tang W, Porcu E, et al. (2011) New gene functions in megakaryopoiesis and platelet formation. Nature 480: 201-208

Knaus UG (2000) Rho GTPase signaling in inflammation and transformation. Immunol Res 21: 103-109. doi: 10.1385/ir:21:2-3:103

Gramaglia I, Weinberg AD, Lemon M, Croft M (1998) Ox-40 ligand: a potent costimulatory molecule for sustaining primary CD4 T cell responses. J Immunol 161: 6510-6517

Thangaratinam S, Tan A, Knox E, Kilby MD, Franklyn J, et al. (2011) Association between thyroid autoantibodies and miscarriage and preterm birth: meta-analysis of evidence. Bmj 342: d2616. doi: 10.1136/bmj.d2616

Willer CJ, Li Y, Abecasis GR (2010) METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26: 2190-2191. doi: 10.1093/bioinformatics/btq340

Milani G, Masciullo C, Sala C, Bellazzi R, Buetti I, et al. (2011) Computer-based genealogy reconstruction in founder populations. J Biomed Inform 44: 997-1003. doi: 10.1016/j.jbi.2011.08.004

Pilia G, Chen WM, Scuteri A, Orru M, Albai G, et al. (2006) Heritability of cardiovascular and personality traits in 6,148 Sardinians. PLoS Genet 2: e132. doi: 10.1371/journal.pgen.0020132