Warning: fopen(/home/virtual/enm-kes/journal/upload/ip_log/ip_log_2024-10.txt): failed to open stream: Permission denied in /home/virtual/lib/view_data.php on line 100 Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 101 <i xmlns="">TSHR</i> Gene (rs179247) Polymorphism and Susceptibility to Autoimmune Thyroid Disease: A Systematic Review and Meta-Analysis

TSHR Gene (rs179247) Polymorphism and Susceptibility to Autoimmune Thyroid Disease: A Systematic Review and Meta-Analysis

Article information

Endocrinol Metab. 2024;39(4):603-614
Publication date (electronic) : 2024 August 1
doi : https://doi.org/10.3803/EnM.2024.1987
1Divisions of Endocrinology, Metabolism, and Diabetes, Thyroid Center, Department of Internal Medicine, Faculty of Medicine, Universitas Syiah Kuala (University Syiah Kuala), Banda Aceh, Indonesia
2Divisions of Endocrinology, Metabolism, and Diabetes, Department of Internal Medicine, Dr. Zainoel Abidin Hospital, Banda Aceh, Indonesia
3Innovation and Research Center of Endocrinology, Faculty of Medicine, Universitas Syiah Kuala (University Syiah Kuala), Banda Aceh, Indonesia
4Department of Internal Medicine, Faculty of Medicine, Pelita Harapan University, Tangerang, Indonesia
Corresponding author: Hendra Zufry Divisions of Endocrinology, Metabolism, and Diabetes, Thyroid Center, Department of Internal Medicine, Faculty of Medicine, Universitas Syiah Kuala (University Syiah Kuala), Jl. Teungku Tanoh Abee, Kopelma Darussalam, Kec. Syiah Kuala, Banda Aceh, Indonesia Tel: +62-8126909131, Fax: +62-65152053, E-mail: hendra_zufry@usk.ac.id
Received 2024 March 14; Revised 2024 May 3; Accepted 2024 May 30.

Abstract

Background

Both Graves’ disease (GD) and Hashimoto’s thyroiditis (HT) are classified as autoimmune thyroid diseases (AITDs). It has been hypothesized that changes in the thyroid-stimulating hormone receptor (TSHR) gene may contribute to the development of these conditions. This study aimed to analyze the correlation between the TSHR rs179247 gene polymorphism and susceptibility to AITD.

Methods

We conducted a thorough search of the Google Scholar, Scopus, Medline, and Cochrane Library databases up until March 2, 2024, utilizing a combination of relevant keywords. This review examines data on the association between TSHR rs179247 and susceptibility to AITD. Random-effect models were employed to assess the odds ratio (OR), and the findings are presented along with their respective 95% confidence intervals (CIs).

Results

The meta-analysis included 12 studies. All genetic models of the TSHR rs179247 gene polymorphism were associated with an increased risk of developing GD. Specifically, the associations were observed in the dominant model (OR, 1.65; P<0.00001), recessive model (OR, 1.65; P<0.00001), as well as for the AA genotype (OR, 2.09; P<0.00001), AG genotype (OR, 1.39; P<0.00001), and A allele (OR, 1.44; P<0.00001). Further regression analysis revealed that these associations were consistent regardless of the country of origin, sample size, age, and sex distribution. However, no association was found between TSHR rs179247 and the risk of HT across all genetic models.

Conclusion

This study suggests that the TSHR rs179247 gene polymorphism is associated with an increased risk of GD, but not with HT, and may therefore serve as a potential biomarker.

GRAPHICAL ABSTRACT

INTRODUCTION

Autoimmune thyroid disease (AITD) encompasses a range of clinical disorders that affect the thyroid gland due to an autoimmune process in which the immune system mistakenly targets thyroid gland cells [1]. AITD includes two forms of thyroid disorders with distinct clinical presentations: Graves’ disease (GD) leads to hyperthyroidism, while Hashimoto thyroiditis (HT) results in hypothyroidism [1]. In GD, B lymphocyte cells produce thyroid-stimulating immunoglobulin, which binds to the thyroid-stimulating hormone receptor (TSHR) on the thyroid gland, leading to the secretion of thyroid hormone and the symptoms of hyperthyroidism [1,2]. GD is the primary cause of hyperthyroidism, accounting for about 60% to 80% of all cases of thyrotoxicosis worldwide [1,2]. In HT, lymphocyte cells, including T cells and B cells, infiltrate the thyroid gland, causing damage and eventually leading to hypothyroidism due to decreased thyroid hormone synthesis [1,3]. Like GD, HT is a common cause of hypothyroidism [1,3].

Globally, the estimated prevalence of GD is approximately 1% to 1.5% of the population, with an annual incidence of 20 to 30 cases per 100,000 people. This number is expected to continue rising [4]. The detrimental effects of GD should not be underestimated [5]. Research indicates that individuals with GD face a 23% higher risk of all-cause mortality than healthy individuals [5]. Meanwhile, HT has a worldwide prevalence of 10% to 12% and an annual incidence of 0.3 to 1.5 cases per 1,000 people [6]. Both GD and HT continue to impose a significant burden on society.

The mechanisms underlying why certain individuals are more susceptible to developing AITD remain elusive [7]. Like other diseases, both genetic and environmental factors can influence an individual’s vulnerability to AITD [7]. Studies have suggested that genetic factors account for approximately 70% to 80% of the susceptibility to AITD, while environmental factors account for the remaining 20% to 30% [7]. These genetic factors are divided into immune regulatory genes and thyroid-specific genes [8]. Among the thyroid-specific genes, the TSHR gene has been extensively studied for its link to AITD [8]. Genomewide association studies have identified significant associations between TSHR and major histocompatibility complex class II polymorphisms, consistently correlating with thyrotropin receptor antibody positivity and GD [9,10]. Consequently, it is hypothesized that the cleavage of the A subunit from TSHR may trigger the production of thyroid-stimulating antibodies, potentially leading to the onset of GD [9,10]. Previous research has shown that intronic polymorphisms are particularly intriguing due to their possible role in regulating TSHR expression [11]. Polymorphisms in intronic regions, such as rs179247, may influence the generation of short RNAs and provide multiple start sites for alternative mRNA transcription, which could affect gene expression or posttranslational modifications [11]. However, the link between TSHR gene polymorphism (rs179247) and AITD susceptibility remains unclear and a matter of debate. The aim of this study was to comprehensively investigate the relationship between TSHR gene polymorphism (rs179247) and susceptibility to AITD, encompassing both GD and HT.

METHODS

Eligibility criteria

The Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guideline was selected as the standard for reporting in this review (Supplemental Table S1) [12]. The protocol of this review has been registered in the International Prospective Register of Systematic Reviews (PROSPERO; CRD 42024520343). To be eligible for our study, studies were required to follow a case-control design, with participants consisting of patients diagnosed with AITD, specifically GD or HT (as the case group), and healthy individuals (as the control group). The research was also required to include data regarding TSHR gene polymorphisms, particularly rs179247, which can manifest as major (A) and minor (G) alleles, along with the distribution of genotype variants (AA, AG, and GG) in both patient groups.

Furthermore, we excluded the following studies: (1) research involving individuals with thyroid disorders other than AITD; (2) studies without a control group; (3) studies with limited information on rs179247; (4) studies not available in full-text format; and (5) studies employing methodologies other than the case-control design.

Search strategy and study selection

A thorough literature search was conducted across four international databases: Google Scholar, Scopus, Medline, and Cochrane Library, up until March 6, 2024. The search employed specific keywords: “(Autoimmune thyroid disease OR AITD OR Graves’ Disease OR GD OR Graves’ hyperthyroidism OR Basedow disease OR Hashimoto Thyroiditis OR HT OR Hashimoto’s disease OR autoimmune thyroiditis) AND (thyroid stimulating hormone receptor gene OR TSHR gene OR rs179247).” Two authors independently performed an initial screening of these databases using the aforementioned keywords. Additionally, they verified the citations exported to the reference manager to ensure their consistency and completeness. Citation tracking was utilized to identify additional relevant papers, which involved reviewing the references of the initially identified research, following citations, and exploring related articles. Two authors independently reviewed the titles and abstracts, excluding those that were not relevant to the study. Subsequently, the same two authors assessed the eligibility of the studies based on the full texts, resolving any discrepancies through discussion.

Data extraction

Two authors separately collected data from the eligible articles, such as participant characteristics (age and sex distribution) and study details (author’s last name, publication year, country, Hardy-Weinberg equation [HWE] testing, the genotyping method, and the number of participants).

Risk of bias assessment

Two authors independently evaluated each publication for potential bias using standardized tools. The Newcastle-Ottawa Scale (NOS) was employed to assess the overall quality of each observational study [13]. This scale measures the selection of study subjects, the comparability of participant groups, and the outcomes of the studies [13]. The scoring range of the instrument is from 0 to 9 [13]. Studies achieving a total score of 7 or higher were considered to be of good quality [13].

Statistical analysis

The outcomes, expressed as dichotomous variables, were calculated using odds ratios (ORs) with 95% confidence intervals (95% CIs), employing the Mantel-Haenszel technique. This approach facilitated the comparison of the TSHR (rs179247) gene polymorphism between two patient groups. Due to expected substantial heterogeneity stemming from variations in population characteristics, a random-effect model was used. The study included three common genetic models—dominant, recessive, and allele type—as well as three additional genotype models (AA vs. GG, AA vs. AG, and AG vs. GG). These models were used to thoroughly assess the association between the TSHR (rs179247) gene polymorphism and susceptibility to AITD. We used the I-squared (I2) statistic to evaluate heterogeneity among studies [14]. I2 values were categorized as ≤25%, 26%–50%, and >50%, indicating low, moderate, and high heterogeneity, respectively [14]. Meta-regression with a random-effects model was conducted using a restricted-maximum likelihood approach for pre-specified variables, including country of origin, sample size, mean age, and sex distribution. This analysis aimed to explore the interaction effects between TSHR gene polymorphism and these variables on the risk of AITD. A publication bias analysis was performed when at least 10 studies were available for each outcome of interest. Statistical analyses were carried out using Review Manager 5.4, a software developed by the Cochrane Collaboration (London, UK), and Comprehensive Meta-Analysis version 3 (Biostat, Englewood, NJ, USA).

Ethics approval

This was a systematic review and meta-analysis. The Faculty of Medicine, University of Syiah Kuala Research Ethics Committee confirmed that no ethical approval was required.

RESULTS

Study selection and characteristics

We searched four databases: Google Scholar (n=129), Scopus (n=13), Medline (n=12), and Cochrane Library (n=8). In total, we identified 162 citations. After removing 136 duplicates and ineligible citations based on title/abstract screening, we examined 26 records. Of these, 14 papers were excluded for the following reasons: eight were review articles, two had incomplete data for analysis, two lacked a control group, one had missing data on the outcomes of interest, and one was an abstract-only article. The list of excluded studies is available in Supplemental Table S2. The final selection included 12 papers [15-26] involving 8,594 patients with GD, 536 individuals with HT, and 8,543 healthy controls (Fig. 1). Among these studies, six originated from China, two from Japan, and one each from Brazil; Poland, the United Kingdom (UK), and a combined effort from the UK and Poland. The sample sizes varied, with the GD group ranging from 112 to 2,504, the HT group from 57 to 230, and the healthy control group from 56 to 2,784 across the studies reviewed. Most studies employed TaqMan real-time polymerase chain reaction for genotyping, with one exception that used matrix-assisted laser desorption/ionization coupled to time-of-flight mass spectrometry genotyping. Table 1 provides a comprehensive overview of the baseline parameters for each paper analyzed.

Fig. 1.

Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) diagram showing the detailed process of study selection for inclusion in the systematic review and meta-analysis.

Characteristics of the Included Studies

Quality of study assessment

All case-control studies included in our analysis were of high quality, scoring between 8 and 9 on the NOS, as shown in Table 2. Each study met the criteria for inclusion in the meta-analysis.

Quality Assessment of Observational Studies Using the Newcastle-Ottawa Scale

Graves’ disease

Classical model

1) Dominant (AA+AG vs. GG)

Our pooled analysis of datasets from 12 studies, encompassing 16,570 participants, demonstrated that the dominant model of TSHR gene polymorphism (rs179247) (AA+AG vs. GG) is associated with an increased risk of developing GD (OR, 1.65; 95% CI, 1.51 to 1.80; P<0.00001; I2=0%; random-effect model) (Fig. 2A).

Fig. 2.

Forest plot demonstrating the association between classical genetic models of the relationship between the thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism and the risk of Graves’ disease: (A) dominant (AA+AG vs. GG), (B) recessive (AA vs. AG+GG), and (C) allele (A vs. G). M-H, Mantel-Haenszel; CI, confidence interval.

2) Recessive (AA vs. AG+GG)

Based on our pooled analysis of datasets from 12 studies, involving a total of 16,570 participants, we found that the recessive model of TSHR gene polymorphism (rs179247) (AA vs. AG+GG) is associated with an increased risk of developing GD (OR, 1.65; 95% CI, 1.52 to 1.80; P<0.00001; I2=30%; random-effect model) (Fig. 2B).

3) Alleles (A vs. G)

A meta-analysis of datasets from 13 datasets derived from 11 studies (n=30,884), demonstrated that the A allele of the TSHR (rs179247) gene polymorphism is associated with a higher risk of developing GD compared to the G allele (OR, 1.44; 95% CI, 1.37 to 1.52; P<0.00001; I2=13%; random-effect model) (Fig. 2C).

Additional model

1) AA vs. GG

A meta-analysis including 14 datasets from 12 studies (n=9,106), demonstrated that the AA genotype of the TSHR (rs179247) gene polymorphism is associated with a higher risk of GD than the GG genotype (OR, 2.09; 95% CI, 1.89 to 2.31; P<0.00001; I2=0%; random-effect model) (Supplemental Fig. S1A).

2) AA vs. AG

A meta-analysis involving 14 datasets from 12 studies (n=14,006) showed that the genetic model consisting of the AA genotype of the TSHR (rs179247) gene polymorphism was associated with a higher risk of GD than the AG genotype (OR, 1.54; 95% CI, 1.40 to 1.68; P<0.0001, I2=31%; random-effect model) (Supplemental Fig. S1B).

3) AG vs. GG

A meta-analysis of 14 datasets from 12 studies (n=10,024) showed that the genetic model consisting of the AG genotype of TSHR (rs179247) gene polymorphism was associated with a higher risk of GD than the GG genotype (OR, 1.39; 95% CI, 1.27 to 1.53; P<0.00001; I2=0%; random-effect model) (Supplemental Fig. S1C).

Hashimoto thyroiditis

Classical model

1) Dominant (AA+AG vs. GG)

A pooled analysis of four studies (n=1,083) showed no association between the dominant model of TSHR (rs179247) gene polymorphism (AA+AG vs. GG) and the risk of developing HT (OR, 0.70; 95% CI, 0.40 to 1.24; P=0.23; I2=59%; random-effect model) (Fig. 3A).

Fig. 3.

Forest plot demonstrating the association between classical genetic models of the relationship between the thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism and the risk of Hashimoto’s thyroiditis: (A) dominant (AA+AG vs. GG), (B) recessive (AA vs. AG+GG), and (C) allele (A vs. G). M-H, Mantel-Haenszel; CI, confidence interval.

2) Recessive (AA vs. AG+GG)

A pooled analysis including four studies (n=1,083) showed no association between the recessive model of the TSHR (rs179247) gene polymorphism (AA vs. AG+GG) and the risk of developing HT (OR, 1.01; 95% CI, 0.71 to 1.43; P=0.96; I2=37%; random-effect model) (Fig. 3B).

3) Alleles (A vs. G)

A pooled analysis of four studies (n=2,156) showed no association between the A versus G allele of the TSHR (rs179247) gene polymorphism and the risk of developing HT (OR, 0.90; 95% CI, 0.70 to 1.15; P=0.41; I2=42%; random-effect model) (Fig. 3C).

Additional model

1) AA vs. GG

A pooled analysis of four studies (n=591) revealed no association between the AA genotype of the TSHR (rs179247) gene polymorphism and the risk of developing HT when compared to the GG genotype (OR, 0.80; 95% CI, 0.40 to 1.60; P=0.53; I2=63%; random-effect model) (Supplemental Fig. S2A).

2) AA vs. AG

A pooled analysis including four studies (n=889) revealed no association between the AA genotype of the TSHR (rs179247) gene polymorphism with the risk of developing HT when compared to the AG genotype (OR, 1.12; 95% CI, 0.77 to 1.62; P= 0.56; I2=36%; random-effect model) (Supplemental Fig. S2B).

3) AG vs. GG

A pooled analysis of four studies (n=686) revealed no association between AG genotype of the TSHR (rs179247) gene polymorphism with the risk of developing HT when compared to the GG genotype (OR, 0.67; 95% CI, 0.38 to 1.19; P=0.18; I2=55%; random-effect model) (Supplemental Fig. S2C).

Meta-regression

Meta-regression was conducted solely for the GD genetic models, as HT did not produce any statistically significant results. The results of the meta-regression analysis, which explored several factors that might influence the relationship between the TSHR (rs179247) gene polymorphism and GD, are presented in Table 3. Factors such as country of origin (with Asia as the reference) (P=0.7361), sample size (P=0.8064), age (P=0.9804), and sex distribution (P=0.9403) did not significantly impact the association between the TSHR (rs179247) gene polymorphism and GD risk in the dominant model (AA+AG vs. GG). Similarly, in the recessive model (AA vs. AG+GG), none of these factors, including country of origin (Asia as reference) (P=0.7408), sample size (P=0.1194), age (P=0.6127), and sex distribution (P=0.9327), significantly influenced the risk of GD. Furthermore, in the allele model (A vs. G) of the TSHR (rs179247) gene polymorphism, the analysis showed that none of the examined factors—country of origin (Asia as reference) (P=0.3731), sample size (P=0.3558), age (P=0.8114), and sex distribution (P=0.8939)—affected GD susceptibility. Lastly, the meta-regression analysis indicated that the relationships of all additional genetic models, including AA vs. GG, AA vs. GA, and AG vs. GG, were not significantly influenced by any of the confounding factors, such as country of origin, sample size, age, and sex distribution.

Results of Meta-Regression for Each Genetic Model of the Relationship between the rs179247 Gene Polymorphism and the Development of Graves’ Disease

Publication bias

Funnel plot analysis was employed to assess publication bias. The analysis revealed relatively symmetrical inverted plots across all genetic models in GD, suggesting an absence of publication bias (Supplemental Fig. S3). However, the current study did not assess publication bias for genetic models in HT, as fewer than 10 studies (only four) were included for each outcome of interest. Therefore, the evaluation of publication bias does not achieve the same level of robustness as it would with more than 10 studies available for analysis [27,28].

DISCUSSION

Our systematic review and meta-analysis indicated that all genetic models of the TSHR (rs179247) gene polymorphism significantly elevate the risk of developing GD. These models encompass the dominant model (AA+AG vs. GG), the recessive model (AA vs. AG+GG), the AA genotype, the AG genotype, and the A allele. Further regression analysis also demonstrated that these associations were not significantly influenced by confounding factors such as country of origin (Asia vs. non-Asia), sample size, age, and sex distribution. Conversely, our meta-analysis did not identify any statistically significant associations between all genetic models of the TSHR (rs179247) gene polymorphism and susceptibility to HT. This suggests that other single-nucleotide polymorphisms (SNPs) or genes besides TSHR may influence HT.

Two potential mechanisms may explain the association between TSHR intron 1 (rs179247) gene polymorphism and the increased risk of GD. The first mechanism involves alternative splicing that leads to the production of soluble TSHR, a process known as peripheral tolerance [29,30]. Additionally, post-translational intramolecular cleavage of the TSHR protein could result in the release of the TSHR A-subunit [29,30]. This process is thought to potentially initiate or exacerbate autoimmunity in GD [29,30]. Empirical evidence indicates that TSHR autoantibodies predominantly target the extracellular A-subunit [29,30]. Moreover, it has been demonstrated that administering the TSHR A-subunit via intramuscular injection in a murine model is essential for inducing TSHR autoantibody production and the subsequent development of hyperthyroidism [29,30]. Three specific types of TSHR are produced by the major transcripts of TSHR: the full-length TSHR (flTSHR), ST4, and ST5 [29,30]. Both ST4 and ST5 include the first 8 exons of flTSHR and an additional ninth exon [29,30]. Intron 8 is divided into two segments that encode this ninth exon, leading to intron retention [29,30]. The study by Brand et al. [15], included in our meta-analysis, analyzed flTSHR, ST4, and ST5 expression in 12 thyroid tissue samples. The findings indicated that ST4 and ST5 expression was higher in individuals with the AA genotype than in those with the GG genotype [15]. The risk allele (A allele) of TSHR SNP (rs179 247AA) was associated with a reduced ratio of flTSHR to ST4 and flTSHR to ST5 [15]. The increased expression of ST4 and ST5 resulted in elevated synthesis of the soluble “A” subunit of TSHR in the peripheral region [15]. This led to the production of thyroid autoantibodies and the development of GD [15,29,30].

The subsequent mechanism involves the regulation of TSHR expression in the thymus, a process known as central tolerance [30,31]. The risk allele (A) leads to reduced TSHR expression in the thymus, which results in decreased elimination of TSHR self-reactive T cells [30,31]. Consequently, the number of TSHR self-reactive T cells increases [30,31]. As a result, germinal centers in the lymph nodes draining the thyroid can produce more TSHR antibodies, triggering the development of GD [30,31]. Colobran highlighted this concept in his research [32]. Carriers of the risk allele (A) at rs179247 exhibited fewer TSHR mRNA transcripts in the thymic glands of non-autoimmune donors who were homozygous compared to carriers of the protective allele (G) [32]. In heterozygous individuals, the TSHR risk allele (A) was present at a lower level than the protective allele (G) [32].

Our meta-analysis results are consistent with the findings of the 2016 meta-analysis by Xiong et al. [33]. Their study established a correlation between the TSHR (rs179247) gene polymorphism and susceptibility to GD [33]. Although both studies reached similar conclusions, there are several important differences between the study by Xiong et al. [33] and our current investigation.

First, Xiong et al. [33] included only eight papers that investigated the rs179247 gene polymorphism. Of these, five studies focused on the rs179247 gene polymorphism in individuals with GD compared to healthy controls, while the remaining three examined the differences between GD patients with ophthalmopathy and those without [33]. Our meta-analysis, in contrast, encompassed 12 studies that explored the rs179247 gene polymorphism in both individuals with GD and healthy populations. The objective was to assess whether a specific SNP in intron 1 of the TSHR gene could influence the risk of developing GD. Including a greater number of studies in our analysis increases the robustness of the evidence produced.

Second, Xiong et al. [33] employed different approaches to determine the dominant and recessive models of the TSHR rs179247 gene polymorphism. In their study, they identified GG+GA vs. AA as the dominant model and GG vs. GA+AA as the recessive model. This classification impacts the conclusions, suggesting that the dominant and recessive models of the TSHR rs179247 gene polymorphism provide protection against GD [33]. When describing the TSHR intron 1 gene, specifically the rs179247 A/G variant with allele A listed before allele G, the correct notation for the dominant model should be AA+AG vs. GG, and for the recessive model, it should be AA vs. AG+GG. Consequently, in our meta-analysis, we adhered to this standard notation, which led to a revised conclusion that both the dominant and recessive models are associated with an increased risk of GD.

Third, our meta-analysis incorporated meta-regression to evaluate the influence of various confounding factors on the outcomes. Further regression analysis indicated that variables such as country of origin, sample size, age, and sex distribution did not significantly affect the association between the TSHR (rs179247) gene polymorphism and the likelihood of developing GD. Additionally, we examined HT alongside GD, as both conditions are categorized as AITDs. This approach differs from the previous study by Xiong et al. [33], which focused solely on GD. However, our findings did not reveal a significant association between the TSHR (rs179247) gene polymorphism and the risk of developing HT.

This investigation has several limitations. First, the research findings are based on case-control studies, which are susceptible to biases such as information bias, selection bias, and the presence of confounding factors. Consequently, our results should be interpreted with caution. Although a Mendelian randomization study is the optimal design for examining the impact of genetic factors on a specific disease, research of this nature is still lacking on the topic of TSHR gene polymorphisms and GD. Second, the studies included in our analysis did not provide sufficient information on goiter size, thyroid function test results, thyroid antibody test results, and environmental factors that could influence the outcomes. As a result, a comprehensive analysis of how these factors interact with genetic factors is not feasible. Incorporating this data would significantly improve our understanding of the complex pathophysiology of GD and HT. Finally, some of the studies we analyzed had small sample sizes—specifically, fewer than 500 individuals. To validate the findings of our study, further research with larger sample sizes that include data on hormonal and environmental factors is necessary.

In conclusion, our systematic review and meta-analysis have identified an association between the TSHR (rs179247) gene polymorphism and an increased risk of developing GD. This elevated risk was consistent across all genetic models tested, including dominant (AA+AG vs. GG), recessive (AA vs. AG+GG), AA genotype, AG genotype, and the A allele of the TSHR (rs179247) gene polymorphism. Further regression analysis indicated that this association is not affected by factors such as country of origin, sample size, age, and sex distribution. However, our analysis did not find a significant association between the TSHR (rs179247) gene polymorphism and susceptibility to HT. Therefore, our findings suggest that rs179247 may be a significant biomarker for predicting the risk of GD, but not HT. Further research is necessary to investigate the role of other SNPs in the TSHR gene or adjacent genes in the development of both GD and HT, as our current study focused solely on one SNP (rs179247). Future studies should also consider examining the gene-environment interactions that may influence the development of AITDs.

Supplementary Material

Supplemental Table S1.

Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2020 Checklist [12]

enm-2024-1987-Supplemental-Table-S1.pdf

Supplemental Table S2.

List of Excluded Studies with Their Reasons of Exclusion

enm-2024-1987-Supplemental-Table-S2.pdf

Supplemental Fig. S1.

Forest plot that demonstrates the association between additional genetic models of thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism: (A) AA vs. GG, (B) AA vs. AG, and (C) AG vs. GG with the risk of Graves’ disease. M-H, Mantel-Haenszel; CI, confidence interval.

enm-2024-1987-Supplemental-Fig-S1.pdf

Supplemental Fig. S2.

Forest plot that demonstrates the association between additional genetic models of thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism: (A) AA vs. GG, (B) AA vs. AG, and (C) AG vs. GG with the risk of Hashimoto thyroiditis. M-H, Mantel-Haenszel; CI, confidence interval.

enm-2024-1987-Supplemental-Fig-S2.pdf

Supplemental Fig. S3.

Funnel plot analysis that demonstrates relatively symmetrical inverted plot for the association between genetic models of thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism: (A) dominant, (B) recessive, (C) allele, (D) AA vs. GG, (E) AA vs. AG, and (F) AG vs. GG with the risk of Graves’ disease. SE, standard error; OR, odds ratio;

enm-2024-1987-Supplemental-Fig-S3.pdf

Notes

CONFLICTS OF INTEREST

No potential conflict of interest relevant to this article was reported.

AUTHOR CONTRIBUTIONS

Conception or design: H.Z., T.I.H. Acquisition, analysis, or interpretation of data: H.Z., T.I.H. Drafting the work or revising: H.Z., T.I.H. Final approval of the manuscript: H.Z., T.I.H.

References

1. Antonelli A, Ferrari SM, Corrado A, Di Domenicantonio A, Fallahi P. Autoimmune thyroid disorders. Autoimmun Rev 2015;14:174–80.
2. Davies TF, Andersen S, Latif R, Nagayama Y, Barbesino G, Brito M, et al. Graves’ disease. Nat Rev Dis Primers 2020;6:52.
3. Ragusa F, Fallahi P, Elia G, Gonnella D, Paparo SR, Giusti C, et al. Hashimotos’ thyroiditis: epidemiology, pathogenesis, clinic and therapy. Best Pract Res Clin Endocrinol Metab 2019;33:101367.
4. Kahaly GJ. Management of Graves thyroidal and extrathyroidal disease: an update. J Clin Endocrinol Metab 2020;105:3704–20.
5. Okosieme OE, Taylor PN, Evans C, Thayer D, Chai A, Khan I, et al. Primary therapy of Graves’ disease and cardiovascular morbidity and mortality: a linked-record cohort study. Lancet Diabetes Endocrinol 2019;7:278–87.
6. Ahmed R, Al-Shaikh S, Akhtar M. Hashimoto thyroiditis: a century later. Adv Anat Pathol 2012;19:181–6.
7. Brent GA. Environmental exposures and autoimmune thyroid disease. Thyroid 2010;20:755–61.
8. Tomer Y. Mechanisms of autoimmune thyroid diseases: from genetics to epigenetics. Annu Rev Pathol 2014;9:147–56.
9. Chu X, Pan CM, Zhao SX, Liang J, Gao GQ, Zhang XM, et al. A genome-wide association study identifies two new risk loci for Graves’ disease. Nat Genet 2011;43:897–901.
10. Chen CR, Pichurin P, Nagayama Y, Latrofa F, Rapoport B, McLachlan SM. The thyrotropin receptor autoantigen in Graves disease is the culprit as well as the victim. J Clin Invest 2003;111:1897–904.
11. Dechairo BM, Zabaneh D, Collins J, Brand O, Dawson GJ, Green AP, et al. Association of the TSHR gene with Graves’ disease: the first disease specific locus. Eur J Hum Genet 2005;13:1223–30.
12. Page MJ, McKenzie JE, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. The PRISMA 2020 statement: an updated guideline for reporting systematic reviews. BMJ 2021;372:n71.
13. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol 2010;25:603–5.
14. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003;327:557–60.
15. Brand OJ, Barrett JC, Simmonds MJ, Newby PR, McCabe CJ, Bruce CK, et al. Association of the thyroid stimulating hormone receptor gene (TSHR) with Graves’ disease. Hum Mol Genet 2009;18:1704–13.
16. Bufalo NE, Dos Santos RB, Marcello MA, Piai RP, Secolin R, Romaldini JH, et al. TSHR intronic polymorphisms (rs179247 and rs12885526) and their role in the susceptibility of the Brazilian population to Graves’ disease and Graves’ ophthalmopathy. J Endocrinol Invest 2015;38:555–61.
17. Fujii A, Inoue N, Watanabe M, Kawakami C, Hidaka Y, Hayashizaki Y, et al. TSHR gene polymorphisms in the enhancer regions are most strongly associated with the development of Graves’ disease, especially intractable disease, and of Hashimoto’s disease. Thyroid 2017;27:111–9.
18. Inoue N, Watanabe M, Katsumata Y, Hidaka Y, Iwatani Y. Different genotypes of a functional polymorphism of the TSHR gene are associated with the development and severity of Graves’ and Hashimoto’s diseases. Tissue Antigens 2013;82:288–90.
19. Liu L, Wu HQ, Wang Q, Zhu YF, Zhang W, Guan LJ, et al. Association between thyroid stimulating hormone receptor gene intron polymorphisms and autoimmune thyroid disease in a Chinese Han population. Endocr J 2012;59:717–23.
20. Ploski R, Brand OJ, Jurecka-Lubieniecka B, Franaszczyk M, Kula D, Krajewski P, et al. Thyroid stimulating hormone receptor (TSHR) intron 1 variants are major risk factors for Graves’ disease in three European Caucasian cohorts. PLoS One 2010;5e15512.
21. Rydzewska M, Goralczyk A, Goscik J, Wawrusiewicz-Kurylonek N, Bossowska A, Kretowski A, et al. Analysis of chosen polymorphisms rs2476601 a/G - PTPN22, rs1990760 C/T - IFIH1, rs179247 a/G - TSHR in pathogenesis of autoimmune thyroid diseases in children. Autoimmunity 2018;51:183–90.
22. Sun W, Zhang X, Wu J, Zhao W, Zhao S, Li M. Correlation of TSHR and CTLA-4 single nucleotide polymorphisms with Graves disease. Int J Genomics 2019;2019:6982623.
23. Wang H, Yang L, Jia Z, Wang B, Sun Y, Pan C, et al. Association of TSHR gene polymorphisms and haplotypes with Graves’ disease in Han population from coastal areas in Shandong province. Chin J Endocrinol Metab 2013;29:108–11.
24. Wang R, Zhang X, Liu B, Yang J, Teng X, Sun W, et al. TSHR gene intron 1 region polymorphism and its relationship with Graves’ disease. Chin J Endocrinol Metab 2012;28:306–10.
25. Wu C, Zhang X, Shun W, Wang R, Zhao W, Liu T. The relationship between thyroid stimulating hormone receptor gene intron 1 area polymorphism and Graves disease and its clinical manifestations in Han Chinese. J Bengbu Medl Coll 2016;41:1405–12. 1405-7, 12.
26. Yang S, Liu W, Xue L, Zhao S, Pan C, Liang J, et al. Association of TSHR gene intron 1 polymorphisms with Graves’ disease. Chin J Endocrinol Metab 2011;27:478–81.
27. Thornton A, Lee P. Publication bias in meta-analysis: its causes and consequences. J Clin Epidemiol 2000;53:207–16.
28. Terrin N, Schmid CH, Lau J, Olkin I. Adjusting for publication bias in the presence of heterogeneity. Stat Med 2003;22:2113–26.
29. Marin-Sanchez A, Alvarez-Sierra D, Gonzalez O, LucasMartin A, Selles-Sanchez A, Rudilla F, et al. Regulation of TSHR expression in the thyroid and thymus may contribute to TSHR tolerance failure in Graves’ disease patients via two distinct mechanisms. Front Immunol 2019;10:1695.
30. McLachlan SM, Rapoport B. Breaking tolerance to thyroid antigens: changing concepts in thyroid autoimmunity. Endocr Rev 2014;35:59–105.
31. Pujol-Borrell R, Alvarez-Sierra D, Jaraquemada D, Marin-Sanchez A, Colobran R. Central tolerance mechanisms to TSHR in Graves’ disease: contributions to understand the genetic association. Horm Metab Res 2018;50:863–70.
32. Colobran R, Armengol Mdel P, Faner R, Gartner M, Tykocinski LO, Lucas A, et al. Association of an SNP with intrathymic transcription of TSHR and Graves’ disease: a role for defective thymic tolerance. Hum Mol Genet 2011;20:3415–23.
33. Xiong H, Wu M, Yi H, Wang X, Wang Q, Nadirshina S, et al. Genetic associations of the thyroid stimulating hormone receptor gene with Graves diseases and Graves ophthalmopathy: a meta-analysis. Sci Rep 2016;6:30356.

Article information Continued

Fig. 1.

Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) diagram showing the detailed process of study selection for inclusion in the systematic review and meta-analysis.

Fig. 2.

Forest plot demonstrating the association between classical genetic models of the relationship between the thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism and the risk of Graves’ disease: (A) dominant (AA+AG vs. GG), (B) recessive (AA vs. AG+GG), and (C) allele (A vs. G). M-H, Mantel-Haenszel; CI, confidence interval.

Fig. 3.

Forest plot demonstrating the association between classical genetic models of the relationship between the thyroid-stimulating hormone receptor (TSHR) rs179247 gene polymorphism and the risk of Hashimoto’s thyroiditis: (A) dominant (AA+AG vs. GG), (B) recessive (AA vs. AG+GG), and (C) allele (A vs. G). M-H, Mantel-Haenszel; CI, confidence interval.

Table 1.

Characteristics of the Included Studies

Study Country HWE test Genotyping method Cases
Control
Sample size Age, yr (mean±SD) Male, % Sample size Age, yr (mean±SD) Male, %
Graves’ disease
 Brand et al. (2009) [15] UK 0.454 TaqMan PCR 768 NR NR 768 NR NR
 Bufalo et al. (2015) [16] Brazil 0.213 TaqMan PCR 279 39.8±11.7 17.2 296 36.8±12.9 20.2
 Fujii et al. (2017) [17] Japan >0.050 TaqMan PCR 180 NR NR 111 NR NR
 Inoue et al. (2013) [18] Japan 0.066 TaqMan PCR 112 34.3±14.5 15.1 56 46.3±12.5 19.6
 Liu et al. (2012) [19] China 0.504 MALDI-TOF-MS 404 34.2±13.7 28.2 242 34.8±12.8 27.7
 Ploski et al. (2010) [20]
  Cohort 1 Poland 0.443 TaqMan PCR 558 39.6±14.2 19.7 520 39.8±13.7 18.5
  Cohort 2 Poland 0.857 TaqMan PCR 196 44.2±10.8 17.9 198 45.3±12.3 18.8
  Cohort 3 UK 0.393 TaqMan PCR 2,504 43.1±11.3 17.2 2,784 44.2±11.7 18.6
 Rydzewska et al. (2018) [21] Poland 0.710 TaqMan PCR 142 16.5±2.0 24.6 160 16.3±3.0 53.1
 Sun et al. (2019) [22] China 0.800 TaqMan PCR 597 38.8±14.2 21.2 620 44.3±14.1 18.2
 Wang et al. (2013) [23] China 0.713 TaqMan PCR 471 NR NR 472 NR NR
 Wang et al. (2012) [24] China 0.859 TaqMan PCR 618 NR NR 646 NR NR
 Wu et al. (2016) [25] China >0.050 TaqMan PCR 699 NR NR 563 NR NR
 Yang et al. (2011) [26] China 0.471 TaqMan PCR 1,066 NR NR 1,107 NR NR
Hashimoto’s thyroiditis
 Fujii et al. (2017) [17] Japan >0.050 TaqMan PCR 151 NR NR 111 NR NR
 Inoue et al. (2013) [18] Japan >0.050 TaqMan PCR 98 45.9±14.9 18.3 56 46.3±12.5 19.6
 Liu et al. (2012) [19] China >0.050 MALDI-TOF-MS 230 31.5±12.9 12.1 242 34.8±12.8 27.7
 Rydzewska et al. (2018) [21] Poland 0.650 TaqMan PCR 57 15.2±2.2 17.5 160 16.3±3.0 53.1

HWE, Hardy-Weinberg equation; SD, standard deviation; PCR, polymerase chain reaction; NR, not reported; MALDI-TOF-MS, matrix-assisted laser desorption/ionization coupled to time-of-flight mass spectrometry.

Table 2.

Quality Assessment of Observational Studies Using the Newcastle-Ottawa Scale

Study Study design Selectiona Comparabilityb Outcomec Total score Result
Brand et al. (2009) [15] Case-control **** ** ** 8 Good
Bufalo et al. (2015) [16] Case-control **** ** *** 9 Good
Fujii et al. (2017) [17] Case-control **** ** *** 9 Good
Inoue et al. (2013) [18] Case-control *** ** *** 8 Good
Liu et al. (2012) [19] Case-control **** ** *** 9 Good
Ploski et al. (2010) [20] Case-control **** ** *** 9 Good
Rydzewska et al. (2018) [21] Case-control **** ** *** 9 Good
Sun et al. (2019) [22] Case-control *** ** *** 8 Good
Wang et al. (2013) [23] Case-control *** ** *** 8 Good
Wang et al. (2012) [24] Case-control *** ** *** 8 Good
Wu et al. (2016) [25] Case-control *** ** *** 8 Good
Yang et al. (2011) [26] Case-control *** ** *** 8 Good
a

(1) is the case definition adequate, (2) representativeness of the cases, (3) selection of controls, (4) definition of controls;

b

(1) comparability of cases and controls on the basis of design or analysis (maximum two stars);

c

(1) ascertainment of exposure, (2) same method of ascertainment for cases and controls, (3) non-response rate.

Table 3.

Results of Meta-Regression for Each Genetic Model of the Relationship between the rs179247 Gene Polymorphism and the Development of Graves’ Disease

Covariate Coefficient 95% CI (min) 95% CI (max) SE P value
Dominant (AA+AG vs. GG)
 Country (Asia as reference) 0.0321 –0.1545 0.2187 0.0952 0.7361
 Sample size, n –0.0000 –0.0000 0.0000 0.0000 0.8064
 Age, yr 0.0002 –0.0182 0.0187 0.0094 0.9804
 Male sex prevalence, % –0.0008 –0.0224 0.0208 0.0110 0.9403
Recessive (AA vs. AG+GG)
 Country (Asia as reference) 0.0314 –0.1544 0.2171 0.0948 0.7408
 Sample size, n –0.0000 –0.0001 0.0000 0.0000 0.1194
 Age, yr –0.0065 –0.0315 0.0185 0.0128 0.6127
 Male sex prevalence, % –0.0012 –0.0288 0.0264 0.0141 0.9327
Allele (A vs. G)
 Country (Asia as reference) 0.0539 –0.0648 0.1727 0.0606 0.3731
 Sample size, n –0.0000 –0.0001 0.0000 0.0000 0.3558
 Age, yr –0.0016 –0.0150 0.0117 0.0068 0.8114
 Male sex prevalence, % –0.0010 –0.0163 0.0142 0.0078 0.8939
AA vs. GG
 Country (Asia as reference) 0.0034 –0.2183 0.2251 0.1131 0.9759
 Sample size, n –0.0000 –0.0001 0.0000 0.0000 0.1741
 Age, yr –0.0099 –0.0351 0.0153 0.0129 0.4427
 Male sex prevalence, % 0.0086 –0.0205 0.0377 0.0149 0.5615
AA vs. AG
 Country (Asia as reference) 0.0977 –0.0816 0.2770 0.0915 0.2855
 Sample size, n –0.0000 –0.0001 0.0000 0.0000 0.0506
 Age, yr –0.0073 –0.0335 0.0189 0.0134 0.5864
 Male sex prevalence, % –0.0016 –0.0307 0.0276 0.0149 0.9170
AG vs. GG
 Country (Asia as reference) –0.0832 –0.2805 0.1141 0.1007 0.4087
 Sample size, n 0.0000 –0.0000 0.0001 0.0000 0.5411
 Age, yr 0.0021 –0.0173 0.0215 0.0099 0.8328
 Male sex prevalence, % –0.0017 –0.0243 0.0210 0.0116 0.8866

CI, confidence interval; SE, standard error.