Genome-wide association analysis suggests novel loci for Hashimoto's thyroiditis

Associations between vitamin D levels and dietary patterns in patients with Hashimoto's thyroiditis

Introduction

Vitamin D insufficiency is a global health problem affecting healthy and diseased individuals, including patients with Hashimoto's thyroiditis (HT). Identifying dietary factors that may affect vitamin D levels and providing dietary guidelines accordingly can alleviate this problem. We therefore aimed to identify still unknown associations of dietary patterns, assessed through the Food Frequency Questionnaire (FFQ) with vitamin D blood levels.

Materials and methods

FFQ was collected from 459 patients from Croatian Biobank of Patients with Hashimoto's thyroiditis (CROHT), while total 25(OH)D was measured from their stored serum samples. We performed linear regression analysis between vitamin D levels and weekly intake of 24 food groups in 459 patients with HT (ALL), and in two disease-severity groups (MILD and OVERT).

Results

The main results of our study are observations of: (1) an inverse association between vitamin D levels and coffee consumption (ALL: β = -0.433, p = 0.005; OVERT: β = -0.62, p = 0.008); (2) an inverse association between vitamin D levels and sweets consumption (ALL: β = -0.195, p = 0.034; OVERT: β = -0.431, p = 0.006); (3) positive association between vitamin D levels and vegetable consumption (ALL: β = 0.182, p = 0.019; OVERT, β = 0.311, p = 0.009). Importantly, effect sizes of all three associations were more prominent in HT patients with prolonged and more severe disease (OVERT).

Conclusion

Further research into the functional and causal relationships of the observed associations is important to provide guidance regarding coffee/sugar intake on vitamin D status. A well-balanced diet can help prevent vitamin D deficiency and improve the quality of life of patients with HT, especially those in later stages of disease characterized by greater metabolic imbalance.

Presence or severity of Hashimoto's thyroiditis does not influence basal calcitonin levels: observations from CROHT biobank

PURPOSE

The influence of Hashimoto's thyroiditis (HT) on calcitonin (Ct) production is unresolved question. The aim of this study was to explore if basal Ct levels are influenced by the presence/severity of HT or correlated with clinical phenotypes of HT patients.

METHODS

We included 467 HT patients and 184 control participants, from Croatian Biobank of HT patients (CROHT), in this retrospective study. Calcitonin levels between HT patients and controls were compared using Mann-Whitney test. Ct levels between two subgroups of HT patients, divided by intake of levothyroxine (LT4) therapy, were additionally tested to take into account the illness severity. Spearman rank correlation test was used to analyze correlations between Ct levels and 14 relevant phenotypes.

RESULTS

We have not detected significant differences in median Ct levels between HT patients and controls (2.2 vs 2.35 pg/mL, respectively, P = 0.717) nor in-between two subgroups of HT patients (P = 0.347). We have not detected statistically significant correlations between Ct levels and clinical phenotypes, although we identified three weak nominal correlations: negative correlation of Ct with TgAb in all HT patients (r = - 0.1, P = 0.04); negative correlation of Ct with age in subgroup of HT patients without LT4 therapy (r = - 0.13, P = 0.04); positive correlation of Ct with BSA in subgroup of HT patients on LT4 therapy (r = 0.16, P = 0.042).

CONCLUSION

Our results suggest that HT patients of all disease stages preserve Ct production as healthy individuals and there is no need for Ct measurements in the absence of a nodule. Additional confirmation and clarification of observed nominal correlations are needed due to potential clinical relevance of TgAb and age-dependent Ct decrease in HT women.

Clinical effect of methimazole combined with selenium in the treatment of toxic diffuse goiter in children

BACKGROUND

The incidence of toxic diffuse goiter (Graves’ disease) is higher in adolescents and preschool-aged children, with an upward trend. The incidence at 6–13 years of age is approximately 11.0%, and the incidences in men and women are 7.8% and 14.3%, respectively.

AIM

To explore the clinical effect of methimazole combined with selenium in the treatment of toxic diffuse goiter (Graves’ disease) in children and its effect on serum anti-thyroglobulin antibody (TRAb) and anti-thyroid peroxidase antibody (TPOAb).

METHODS

A total of 103 children with Graves’ disease treated in our hospital from January 2018 to June 2021 were divided into a traditional group and a combined group (15-20 mg methimazole orally given to children) and a combined group (50 µg selenium added on the basis of traditional treatment) according to different treatment methods to explore the therapeutic effects of the two methods and to observe the changes in thyroid volume and serum TRAb, TPOAb, free thyroxine (FT4) and inflammatory factor levels before and after treatment. The time taken for FT4 to return to normal was compared between the two groups.

RESULTS

Treatment was significantly more effective in the combined group than in the traditional group (P < 0.05). The thyroid volumes of the children in the two groups was measured before and after treatment. Thyroid volume decreased significantly after treatment in both groups, and the thyroid volume was significantly lower in the combined group than in the traditional group (P < 0.05). The serum levels of interleukin-6 (IL-6), IL-8, TRAb, TPOAb and FT4 in the two groups were detected before and after treatment. The levels of IL-6, IL-8, TRAb, TPOAb and FT4 were significantly lower in the combined group than in the traditional group (P < 0.05). Follow-up of the children in the two groups showed that compared with the traditional group, it took less time for children in the combined group to return to the normal level (P < 0.05).

CONCLUSION

Methimazole combined with selenium can effectively treat Graves’ disease in children, reduce the expression of TRAb, TPOAb, FT4 and inflammatory factors, and improve the curative effect. Thus, the combined treatment warrants further clinical research.

DNA hypermethylation of PTPN22 gene promoter in children and adolescents with Hashimoto thyroiditis

PURPOSE

Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is an inhibitor of T-cell activation, regulating intracellular signal transduction and thereby being implicated in the pathogenesis of autoimmune thyroid disease (AITD). The exact molecular mechanisms have not been fully elucidated. The aim of the present study was to quantitate DNA methylation within the PTPN22 gene promoter in children and adolescents with AITD and healthy controls.

METHODS

60 Patients with Hashimoto thyroiditis (HT), 25 patients with HT and type 1 diabetes (HT + T1D), 9 patients with Graves' disease (GD) and 55 healthy controls without any individual or family history of autoimmune disease were enrolled. Whole blood DNA extraction, DNA modification using sodium bisulfate and quantification of DNA methylation in the PTPN22 gene promoter, based on melting curve analysis of the selected DNA fragment using a Real-Time PCR assay, were implemented.

RESULTS

DNA methylation in the PTPN22 gene promoter was found to be significantly higher in HT patients (39.9 ± 3.1%) in comparison with other study groups (20.3 ± 2.4% for HT + T1D, 32.6 ± 7.8% for GD, 27.1 ± 2.4% for controls, p < 0.001). PTPN22 gene promoter DNA methylation was also associated marginally with thyroid autoimmunity in general (p = 0.059), as well as considerably with thyroid volume (p = 0.004) and the presence of goiter (p = 0.001) but not thyroid function tests.

CONCLUSIONS

This study demonstrates for the first time that a relationship between autoimmune thyroiditis and PTPN22 gene promoter DNA methylation state is present, thus proposing another possible etiological association between thyroiditis and abnormalities of PTPN22 function. Further expression studies are required to confirm these findings.

Effects of Thyroid Function on Hemostasis, Coagulation, and Fibrinolysis: A Mendelian Randomization Study

Background: Untreated hypothyroidism is associated with acquired von Willebrand syndrome, and hyperthyroidism is associated with increased thrombosis risk. However, the causal effects of thyroid function on hemostasis, coagulation, and fibrinolysis are unknown. Methods: In a two-sample Mendelian randomization (MR) study with genome-wide association variants, we assessed causality of genetically predicted hypothyroidism (N = 134,641), normal-range thyrotropin (TSH; N = 54,288) and free thyroxine (fT4) (N = 49,269), hyperthyroidism (N = 51,823), and thyroid peroxidase antibody positivity (N = 25,821) on coagulation (activated partial thromboplastin time, von Willebrand factor [VWF], factor VIII [FVIII], prothrombin time, factor VII, fibrinogen) and fibrinolysis (D-dimer, tissue plasminogen activator [TPA], plasminogen activator inhibitor-1) from the CHARGE Hemostasis Consortium (N = 2583-120,246). Inverse-variance-weighted random effects were the main MR analysis followed by sensitivity analyses. Two-sided p < 0.05 was nominally significant, and p < 0.0011[ = 0.05/(5 exposures × 9 outcomes)] was Bonferroni significant for the main MR analysis. Results: Genetically increased TSH was associated with decreased VWF [β(SE) = -0.020(0.006), p = 0.001] and with decreased fibrinogen [β(SE) = -0.008(0.002), p = 0.001]. Genetically increased fT4 was associated with increased VWF [β(SE) = 0.028(0.011), p = 0.012]. Genetically predicted hyperthyroidism was associated with increased VWF [β(SE) = 0.012(0.004), p = 0.006] and increased FVIII [β(SE) = 0.013(0.005), p = 0.007]. Genetically predicted hypothyroidism and hyperthyroidism were associated with decreased TPA [β(SE) = -0.009(0.024), p = 0.024] and increased TPA [β(SE) = 0.022(0.008), p = 0.008], respectively. MR sensitivity analyses showed similar direction but lower precision. Other coagulation and fibrinolytic factors were inconclusive. Conclusions: In the largest genetic studies currently available, genetically increased TSH and fT4 may be associated with decreased and increased synthesis of VWF, respectively. Since Bonferroni correction may be too conservative given the correlation between the analyzed traits, we cannot reject nominal associations of thyroid traits with coagulation or fibrinolytic factors.

Vitamin D and Hashimoto's Thyroiditis: Observations from CROHT Biobank

The aims of this study were to evaluate: (1) associations of vitamin D with the presence/severity of Hashimoto’s thyroiditis (HT) and (2) correlations of vitamin D with thyroid-related phenotypes. Total 25(OH)D (vitamin D in the text) was measured from stored serum samples of 461 HT patients and 176 controls from a Croatian Biobank of HT patients (CROHT). (1) Vitamin D levels, and proportions of vitamin D deficiency, were compared between HT cases and controls. HT patients were additionally divided into two groups (MILD and OVERT) to take into account HT severity. (2) Correlations between vitamin D and 10 clinical phenotypes in all HT patients and two subgroups of HT patients were tested using the Spearman correlation test. Our analyses were adjusted for age, gender, BMI, smoking status and seasonality of blood sampling. (1) No significant differences in vitamin D levels, or proportions of vitamin D deficiency, were detected between HT patients of all disease stages and controls. However, a nominally significant difference in vitamin D levels between MILD and OVERT subgroups (OR = 1.038, p = 0.023) was observed. Proportions of individuals with vitamin D deficiency during winter–spring were high: all HT cases (64.69%), MILD (60.64%), OVERT (68.7%), controls (60.79%). (2) A nominally significant negative correlation between vitamin D and TSH in all HT patients (r = −0.113, p = 0.029) and a positive correlation between vitamin D and systolic blood pressure in OVERT HT patients (r = 0.205, p = 0.025) were identified. Our study indicates that there is no association between vitamin D and HT; however, there may be a subtle decrease in vitamin D levels associated with overt hypothyroidism.

An update on the pathogenesis of Hashimoto's thyroiditis

It is 70 years since Noel Rose embarked on his pioneering studies that lead to the discovery of autoimmune thyroiditis and the elucidation of Hashimoto's thyroiditis. This short review to honour his passing focuses on the developments in our understanding of the causes and pathogenesis of HT over the last five years. Recent genetic studies have reported heritability estimates for HT and associated diseases for the first time, and emphasised the complexity of the genetic factors involved, including monogenic forms of HT. Environmental factors continue to be elucidated, especially as a side effect of drugs which modulate the immune system therapeutically. Regarding pathogenetic mechanisms, multiple cytokine networks have been identified which involve the thyroid cells in a circuit of escalating proinflammatory effects, such as the expression of inflammasome components, and an array of different defects in T regulatory cells may underlie the loss of self-tolerance to thyroid autoantigens. Finally, a number of studies have revealed fresh insights into disease associations with HT which may have both pathological and clinical significance, the most intriguing of which is a possible direct role of the autoimmune process itself in causing some of the persistent symptoms reported by a minority of patients with levothyroxine-treated HT.

Genetic Susceptibility to Joint Occurrence of Polycystic Ovary Syndrome and Hashimoto's Thyroiditis: How Far Is Our Understanding?

Polycystic ovary syndrome (PCOS) and Hashimoto’s thyroiditis (HT) are endocrine disorders that commonly occur among young women. A higher prevalence of HT in women with PCOS, relative to healthy individuals, is observed consistently. Combined occurrence of both diseases is associated with a higher risk of severe metabolic and reproductive complications. Genetic factors strongly impact the pathogenesis of both PCOS and HT and several susceptibility loci associated with a higher risk of both disorders have been identified. Furthermore, some candidate gene polymorphisms are thought to be functionally relevant; however, few genetic variants are proposed to be causally associated with the incidence of both disorders together.

Thyroid function, sex hormones and sexual function: a Mendelian randomization study

Hypothyroidism and hyperthyroidism are observationally associated with sex hormone concentrations and sexual dysfunction, but causality is unclear. We investigated whether TSH, fT4, hypo- and hyperthyroidism are causally associated with sex hormones and sexual function. We used publicly available summary statistics from genome-wide association studies on TSH and fT4 and hypo- and hyperthyroidism from the ThyroidOmics Consortium (N ≤ 54,288). Outcomes from UK Biobank (women ≤ 194,174/men ≤ 167,020) and ReproGen (women ≤ 252,514) were sex hormones (sex hormone binding globulin [SHBG], testosterone, estradiol, free androgen index [FAI]) and sexual function (ovulatory function in women: duration of menstrual period, age at menarche and menopause, reproductive lifespan, and erectile dysfunction in men). We performed two-sample Mendelian randomization (MR) analyses on summary level, and unweighted genetic risk score (GRS) analysis on individual level data. One SD increase in TSH was associated with a 1.332 nmol/L lower (95% CI: - 0.717,- 1.946; p = 2 × 10^-5) SHBG and a 0.103 nmol/l lower (- 0.051,V0.154; p = 9 × 10^-5) testosterone in two-sample MR, supported by the GRS approach. Genetic predisposition to hypothyroidism was associated with decreased and genetic predisposition to hyperthyroidism with increased SHBG and testosterone in both approaches. The GRS for fT4 was associated with increased testosterone and estradiol in women only. The GRS for TSH and hypothyroidism were associated with increased and the GRS for hyperthyroidism with decreased FAI in men only. While genetically predicted thyroid function was associated with sex hormones, we found no association with sexual function.

Genotype association of IP6K3 gene with Hashimoto’s thyroiditis in Algerian population (Aures region)

Background

Hashimoto’s thyroiditis (HT) is a chronic autoimmune disease of the thyroid gland and is also the main cause of hypothyroidism. A recent genome-wide association study (GWAS) suggested an association of three novel genetic variants with HT in a population of Caucasian origin (Croatian). A case-control study was performed to investigate the association of these three newly suggested genetic variants with HT in a non-Caucasian ethnic group, an Arab-Berber from Algeria.

Three variants (rs12944194 located 206 kb from SDK2, rs791903 inside IP6K3, and rs75201096 inside GNA14) were genotyped using real-time PCR.

Results

There were no significant differences in allele frequencies of the three genetic variants between HT cases and controls. However, the present study showed nominal significance in the genotype distribution of rs791903 (IP6K3 gene) between HT patients and the control group (P = 0.024); we observed a decrease in the frequency of rs791903 recessive homozygotes (CC) in HT cases versus controls (OR = 0.476, P = 0.025).

Conclusion

This is the first study that showed the genotypic association of IP6K3 intronic variant with decreased risk for HT in non-Caucasian, Algerian, population, whereas we did not confirm the association of SDK2 and GNA14 genetic variants with HT. The IP6K3 gene (inositol hexaphosphate kinase 3), located near major histocompatibility complex (MHC), has previously been associated with other common autoimmune diseases beside HT, such as Graves’s disease and rheumatoid arthritis, which is providing more evidence of a good candidacy for the genetic contribution to the development of HT and autoimmunity.

Urinary Iodine and Genetic Predisposition to Hashimoto's Thyroiditis in a Chinese Han Population: A Case-Control Study

Background: We aimed to examine the association of urinary iodine concentration with Hashimoto's thyroiditis (HT) risk, and particularly, to investigate whether the HT-related genetic variations might modify the effects of urinary iodine on HT in the Chinese Han population. Methods: We conducted a case-control study with 1723 Chinese (731 cases, 992 controls). The associations between urinary iodine concentration and HT risk were analyzed using logistic regression models. The effects of interactions between the genetic risk scores (GRSs) and urinary iodine on HT risk were assessed by including the respective interaction terms in the models. We also applied restricted cubic spline regression to estimate the possible nonlinear relationship. The multinomial logistic regression models were performed to determine the associations of urinary iodine with euthyroid-HT and hypothyroidism-HT. Results: After controlling for potential confounders, the odds of HT increased with increasing quartiles of urinary iodine concentration: adjusted odds ratios (ORs) and 95% confidence intervals [CIs] were 1.45 [1.06-1.99], 1.66 [1.17-2.34], and 2.07 [1.38-3.10] for the quartiles 2, 3, and 4, respectively, compared with the first quartile (p for trend <0.001). Multivariable restricted cubic spline regression analysis further demonstrated that there was a near-linear association between urinary iodine concentration and HT risk (p-overall <0.001; p-nonlinear = 0.074). However, we did not find significant interactions between urinary iodine and GRSs on the risk of HT (all p for interaction >0.05). Interestingly, we found that each increment of urinary iodine was associated with a more than twofold increase in the odds of hypothyroidism-HT (adjusted OR = 2.64 [CI = 1.73-4.05]), but not with euthyroid-HT (p > 0.05). Conclusions: Higher urinary iodine concentration was associated with increased risk of HT, and this association was near linear, indicating that increased urinary iodine has a continuous and graded impact on HT risk. Moreover, the iodine-HT association was not modified by genetic predisposition to HT. Interestingly, urinary iodine concentration was significantly associated with increased risk of hypothyroidism.

Differences in food consumption between patients with Hashimoto's thyroiditis and healthy individuals

Food is considered as important environmental factor that plays a role in development of Hashimoto's thyroiditis (HT). The goal of our study was to identify food groups, assessed by food frequency questionnaire, that differ in consumption frequency between 491 patients with HT and 433 controls. We also analysed association of food groups with the wealth of HT-related clinical traits and symptoms. We found significantly increased consumption of animal fat (OR 1.55, p < 0.0001) and processed meat (OR 1.16, p = 0.0012) in HT cases, whereas controls consumed significantly more frequently red meat (OR 0.80, p < 0.0001), non-alcoholic beverages (OR 0.82, p < 0.0001), whole grains (OR 0.82, p < 0.0001) and plant oil (OR 0.87, p < 0.0001). We also observed association of plant oil consumption with increased triiodothyronine levels in HT patients (β = 0.07, p < 0.0001), and, association of olive oil consumption with decreased systolic blood pressure (β = − 0.16, p = 0.001) in HT patients on levothyroxine (LT4) therapy. Analysis of food consumption between HT patients with and without LT4 therapy suggest that patients do not tend to modify their diet upon HT diagnosis in our population. Our study may be of relevance to nutritionists, nutritional therapists and clinicians involved in developing dietary recommendations for HT patients.

AATF and SMARCA2 are associated with thyroid volume in Hashimoto's thyroiditis patients

Thyroid volume of Hashimoto’s thyroiditis (HT) patients varies in size over the course of disease and it may reflect changes in biological function of thyroid gland. Patients with subclinical hypothyroidism predominantly have increased thyroid volume whereas patients with more pronounced hypothyroidism have smaller thyroid volumes. Suggested mechanism for thyroid atrophy is thyrocyte death due to apoptosis. We performed the first genome-wide association study (GWAS) of thyroid volume in two groups of HT patients, depending on levothyroxine (LT4) therapy, and then meta-analysed across. Study included 345 HT patients in total and 6 007 322 common autosomal genetic variants. Underlying hypothesis was that genetic components that are involved in regulation of thyroid volume display their effect in specific pathophysiologic conditions of thyroid gland of HT patients. We additionally performed immunohistochemical analysis using thyroid tissues and analysed differences in expression levels of identified proteins and apoptotic marker between HT patients and controls. We found genome-wide significant association of two loci, both involved in apoptosis, with thyroid volume of HT patients: rs7212416 inside apoptosis-antagonizing transcription factor AATF (P = 8.95 × 10⁻⁹) and rs10738556 near chromatin-remodeling SMARCA2 (P = 2.83 × 10⁻⁸). In immunohistochemical analysis we observed that HT patients with homozygous AATF risk genotypes have decreased AATF expression (0.46-fold, P < 0.0001) and increased apoptosis (3.99-fold, P = 0.0001) in comparison to controls. HT patients with heterozygous SMARCA2 genotypes have decreased SMARCA2 expression, albeit without reaching statistical significance (1.07-fold, P = 0.5876), and significantly increased apoptosis (4.11-fold, P < 0.0001). By two lines of evidence we show that two highly plausible genetic loci, AATF and SMARCA2, may be involved in determining the thyroid volume of HT patients. The results of our study significantly add to the current knowledge of disturbed biological mechanisms in thyroid gland of HT patients.

Interaction of Genetic Variations in NFE2L2 and SELENOS Modulates the Risk of Hashimoto's Thyroiditis

Background: Several single-nucleotide polymorphisms (SNPs) are known to increase the risk of Hashimoto's thyroiditis (HT); such SNPs reside in thyroid-specific genes or in genes related to autoimmunity, inflammation, and/or cellular defense to stress. The transcription factor Nrf2, encoded by NFE2L2, is a master regulator of the cellular antioxidant response. This study aimed to evaluate the impact of genetic variation in NFE2L2 on the risk of developing HT. Methods: In a case-control candidate gene association study, functional SNPs in the NFE2L2 promoter (rs35652124, rs6706649, and rs6721961) were examined either as independent risk factors or in combination with a previously characterized HT risk allele (rs28665122) in the gene SELENOS, encoding selenoprotein S (SelS). A total of 997 individuals from the north of Portugal (Porto) were enrolled, comprising 481 HT patients and 516 unrelated healthy controls. SELENOS and NFE2L2 SNPs were genotyped using TaqMan^® assays and Sanger sequencing, respectively. Odds ratios (ORs) were calculated using logistic regression, with adjustment for sex and age. Expression of SelS was analyzed by immunohistochemistry in thyroid tissue from HT patients and control subjects. Molecular interactions between the Nrf2 and SelS pathways were investigated in thyroid tissues from mice and in rat PCCL3 thyroid follicular cells. Results: When all three NFE2L2 SNPs were considered together, the presence of one or more minor alleles was associated with a near-significant increased risk (OR = 1.43, p = 0.072). Among subjects harboring only major NFE2L2 alleles, there was no increased HT risk associated with heterozygosity or homozygosity for the SELENOS minor allele. Conversely, in subjects heterozygous or homozygous for the SELENOS risk allele, the presence of an NFE2L2 minor allele significantly increased HT risk by 2.8-fold (p = 0.003). Immunohistochemistry showed reduced expression of SelS in thyroid follicular cells of HT patients. In Nrf2 knockout mice, there was reduced expression of SelS in thyroid follicular cells; conversely, in PCCL3 cells, reducing SelS expression caused reduced activity of Nrf2 signaling. Conclusions: The NFE2L2 promoter genotype interacts with the SELENOS promoter genotype to modulate the risk of HT in a Portuguese population. This interaction may be due to a bidirectional positive feedback between the Nrf2 and SelS pathways.

Correlation of TSHR and CTLA-4 Single Nucleotide Polymorphisms with Graves Disease

This study was designed to explore the association between Graves disease (GD) and thyroid-stimulating hormone receptor (TSHR) and cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) single nucleotide polymorphisms (SNPs). We studied a total of 1217 subjects from a Han population in northern Anhui province in China. Six SNPs within TSHR (rs179247, rs12101261, rs2284722, rs4903964, rs2300525, and rs17111394) and four SNPs within CTLA-4 (rs10197319, rs231726, rs231804, and rs1024161) were genotyped via a Taqman probe technique using a Fluidigm EP1 platform. The TSHR alleles rs179247-G, rs12101261-C, and rs4903964-G were negatively correlated with GD, whereas the rs2284722-A and rs17111394-C alleles were positively correlated with GD. Analyzing TSHR SNPs at rs179247, rs2284722, rs12101261, and rs4903964 yielded 8 different haplotypes. There were positive correlations between GD risk and the haplotypes AGTA and AATA (OR = 1.27, 95%CI = 1.07-1.50, P = 0.005; OR = 1.45, 95%CI = 1.21-1.75, P < 0.001, respectively). There were negative correlations between GD risk and the haplotype GGCG (OR = 0.56, 95%CI = 0.46-0.67, P < 0.001). With respect to haplotypes based on SNPs at the TSHR rs2300525 and rs17111394 loci, the CC haplotype was positively correlated with GD risk (OR = 1.32, 95%CI = 1.08-1.60, P = 0.006). Analyzing CTLA-4 SNPs at rs231804, rs1024161, and rs231726 yielded four haplotypes, of which AAA was positively correlated with GD risk (OR = 1.21, 95%CI = 1.02-1.43, P = 0.029). Polymorphisms at rs179247, rs12101261, rs2284722, rs4903964, and rs17111394 were associated with GD susceptibility. Haplotypes of both TSHR and CTLA-4 were additionally related to GD risk.