International Journal of Genomics

International Journal of Genomics / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6982623 | 12 pages | https://doi.org/10.1155/2019/6982623

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

Academic Editor: Byung-Hoon Jeong
Received27 Feb 2019
Revised16 May 2019
Accepted17 Jun 2019
Published03 Sep 2019

Abstract

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 (, , ; , , , respectively). There were negative correlations between GD risk and the haplotype GGCG (, , ). With respect to haplotypes based on SNPs at the TSHR rs2300525 and rs17111394 loci, the CC haplotype was positively correlated with GD risk (, , ). Analyzing CTLA-4 SNPs at rs231804, rs1024161, and rs231726 yielded four haplotypes, of which AAA was positively correlated with GD risk (, , ). 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.

1. Introduction

GD is a common organ-specific autoimmune disease and the most common cause of thyrotoxicosis. At present, however, the molecular mechanisms underlying GD have not been elucidated. GD susceptibility stems from a confluence of genetic, environmental, and immunological factors [1]. In an individual with a genetic predisposition for GD, the disease may develop as a consequence of environmental influences that induce or exacerbate immune dysfunction, ultimately leading to the onset of autoimmunity, with clear evidence for the existence of genetic factors predisposing individuals to autoimmune thyroid disease [2, 3]. A recent study from the National Health and Nutritional Examination Survey further sought to identify the role of genetic susceptibility in the etiology of GD [4]. GD is a complex polygenetic disease, with multiple risk genes influencing its onset. GWAS studies have been very popular in the identification of such susceptibility loci for thyroid autoimmune diseases [513].

TSHR is a primary candidate gene believed to be related to GD susceptibility. TSHR is a specific protein expressed in thyroid cells in the thyroid follicular membrane. TSH regulates both thyroid growth and functionality via TSHR signaling. A study by Zhan et al. [14] identified a novel susceptibility loci for serum TSH levels in Chinese populations using a GWAS approach.

TSHR is a member of the G-protein-coupled receptor superfamily encoded on chromosome 14. The protein is a single 764 amino acid peptide chain encoded for across 10 exons with a molecular weight of 84000 daltons. TSH binding to TSHR promotes G-protein signaling, leading to activation of the cAMP and/or phosphoinositide Ca2+ signal transduction pathways. There are many antithyroid autoantibodies present in the serum of patients with GD, including thyrotrophin receptor antibody (TRAb), thyroglobulin antibody (TGAb), and thyroid peroxidase antibody (TPOAb). TRAb is an antibody which is specific to TSHR, and it is believed to be the autoantibody most important for the development of hyperthyroidism. Most patients with GD exhibit TRAb autoantibodies in the peripheral blood. TRAb is an immunoglobulin (IgG) that, when present in the serum, competes with TSH to bind to TSHR, activating the receptor and inducing biological effects similar to those of TSH. Research has shown that a targeted immunotherapy strategy in mice aimed at disrupting antigen processing and presentation in HLA-DR3 transgenic mice blocks the immune response to TSHR, thus offering a potential avenue for GD treatment [15]. These findings emphasize the importance of immunological factors in driving GD and highlight the potential value of immunotherapy in its treatment. Further work has confirmed that SNPs in the TSHR gene are associated with defects in central immune tolerance that can lead to the onset of autoimmunity [16].

CTLA-4 is a member of the immunoglobulin gene superfamily and a negative regulator of T cell responses that is associated with immune tolerance. CTLA-4 is expressed on the surface of T cells mainly in the form of a dimer, and when it interacts with its cognate ligands, this induces inhibitory signals which terminate T cell activation and proliferation. Polymorphisms in CTLA-4 may alter its functionality such that the activation of T cells cannot be inhibited, resulting in a loss of immune tolerance and the occurrence of autoimmunity, making it vital that normal CTLA-4 activity be maintained. CTLA-4 is a major susceptibility gene associated with autoimmune thyroid disease (AITD). An association study aimed at identifying SNPs in the CTLA-4 gene present in GD patients and control subjects has confirmed that CTLA-4 is indeed a susceptibility gene for GD in the Chinese Han population [17], and studies in children have confirmed this result [18]. Interactions among SNPs at rs231775, rs231779, and rs3087243 significantly increase an individual’s susceptibility to GD [19].

A number of GD susceptibility genes have been identified to date including human leukocyte antigen (HLA) I and II, cluster of differentiation 40 (CD40), TSHR, protein tyrosine phosphatase nonreceptor 22 (PTPN22), interferon-inducible helicase domain 1 (IFIH1), CTLA-4, forkhead Box P3 (FoxP3), Ikaros family of zinc finger3 (IKZF3), FC-receptor-like 3 (FCRL3), and thyroglobulin (TG) [2027]. Of these, TSHR is a thyroid-specific gene and CTLA-4 is an immunoregulatory gene affecting GD development, both of which are related to the pathogenesis of hyperthyroidism. Individual patient disease phenotypes may be derived from interactions between genetic and environmental factors. Therefore, we carried out a study to assess the relationship between TSHR and CTLA-4 SNPs and GD, and we analyzed the interaction between these two genes in this context. Most previous studies have focused on the TSHR coding region; however, in recent years the relationship between the TSHR noncoding region and GD has been increasingly studied. SNPs in noncoding promoter regions can alter gene expression or localization. Similarly, certain introns contain regulatory elements or regulate transcript splicing, thereby allowing intronic SNPs to control gene expression and function. Chu et al. [5] found the susceptible sites of GD through GWAS screening. Referring to the research results and retrieving the database NCBI (https://www.ncbi.nlm.nih.gov/snp/), we therefore focused on six SNPs in intron 1 of the TSHR gene (rs179247, rs12101261, rs2284722, rs4903964, rs2300525, and rs17111394) and four SNPs in the CTLA-4 gene (rs10197319, rs231726, rs231804, and rs1024161) in order to carry out a case-control study in a Han population from northern Anhui province in China.

2. Materials and Methods

We studied a total of 1217 subjects, divided into a case group and a control group. The control group was composed of 620 healthy subjects including 113 males and 507 females. The GD case group was composed of 597 individuals including 127 male patients and 470 female patients. All patients and control subjects were from a Chinese Han population from northern Anhui province and were unrelated to each other. Patients with other autoimmune diseases or a family history thereof were excluded.

Subjects were diagnosed with GD based upon clinical and laboratory examinations that confirmed hyperthyroidism, which was accompanied by symptoms of a high metabolism, diffuse goiter, thyroid ophthalmopathy, and pretibial myxedema, as well as high serum levels of free thyroxine (FT4) and free T3 (FT3), very low levels of circulating thyroid-stimulating hormone (TSH), and positive TRAb circulation. This study was approved by our local ethics committee.

2.1. Genotyping of SNPs

From each individual, a 5 mL sample of peripheral blood was collected in an EDTA-treated tube and use for DNA extraction with a DNA purification kit (Fujifilm Company) based on the provided instructions. All DNA samples were genotyped using Illumina Human660-Quad BeadChips. Illumina BeadStudio 3.3 software was used for genotype clustering. All samples had a mean call rate of 99.8%. The genotypes of the TSHR and CTLA-4 gene SNPs were determined using a Taqman probe technique with a Fluidigm EP1 platform. Polymerase chain reaction (PCR) was employed to amplify each target gene sequence as previously described [5]. Our targets were six SNPs in intron 1 of TSHR (rs179247, rs12101261, rs2284722, rs4903964, rs2300525, and rs17111394) and four in CTLA-4 (rs10197319, rs231726, rs231804, and rs1024161).

2.2. Statistical Analyses

Quantitative data are given as , and differences between groups were compared via -test. Qualitative data were described as percentages or proportions, and differences between groups were compared via the test. The genotype and allele frequencies for TSHR and CTLA-4 gene in cases and controls were assessed by chi-square test. The correlation between SNPs and GD was analyzed via a binary classification single-factor and multifactor logistic regression model, and the OR and 95% CI were calculated to assess the relationship between the genotype and incidence of GD. Rank sum test was used for rank data. SPSS v19.0 was used for all above statistical analyses. The Haploview 4.2 software was used to determine whether the distributions of the 10 assessed loci conformed to a Hardy-Weinberg equilibrium in the control group and analyze the relationship between the haplotype of TSHR and CTLA-4 and GD. Dimensionality reduction (MDR) was used to analyze the relationship between gene-gene interactions and the incidence of GD. All the tests were two-tailed, with a test level of and with as the threshold of significance.

3. Results

Clinical data is shown in Table 1. We used 1.5 U/L as a cutoff value for TRAb levels. A Hardy-Weinberg equilibrium test was performed on the control group, revealing all 10 SNPS in TSHR and CTLA-4 to conform to a Hardy-Weinberg equilibrium () (Table 2).


Basic informationCase group ()Control group ()/

SexMale1271131.780.18
Female470507
Age6.76<0.001
Diffuse goiter015
I186
II386
III10
ExophthalmosYes252
No345
TRAb≥1.5 U/L456
<1.5 U/L96


SNPsCase group ()Control group ()HWE ()
(%) (%)

TSHRrs179247 (A>G)0.80
A895 (74.96)782 (63.06)40.16<0.001
G299 (25.04)458 (36.94)
AA342 (57.29)248 (40.00)39.44<0.001
AG211 (35.34)286 (46.13)
GG44 (7.37)86 (13.87)
rs2284722 (G>A)0.06
G833 (69.77)937 (75.56)10.310.001
A361 (30.23)303 (24.44)
GG291 (48.74)345 (55.65)12.970.002
AG251 (42.04)247 (39.84)
AA55 (9.22)28 (4.51)
rs12101261 (T>C)0.56
T872 (73.03)756 (60.97)39.97<0.001
C322 (26.97)484 (39.03)
TT320 (53.60)234 (37.74)38.60<0.001
CT232 (38.86)288 (46.45)
CC45 (7.54)98 (15.81)
rs4903964 (A>G)0.29
A824 (69.01)701 (56.53)40.49<0.001
G370 (30.99)539 (43.47)
AA279 (46.73)205 (33.06)41.47<0.001
AG266 (44.56)291 (46.94)
GG52 (8.71)124 (20.00)
rs2300525 (T>C)0.68
T836 (70.02)909 (73.31)3.240.072
C358 (29.98)331 (26.69)
TT297 (49.75)335 (54.03)3.260.196
CT242 (40.54)239 (38.55)
CC58 (9.71)46 (7.42)
rs17111394 (T>C)0.61
T913 (76.47)1001 (80.73)6.570.010
C281 (23.53)239 (19.27)
TT355 (59.46)406 (65.48)6.550.038
CT203 (34.00)189 (30.48)
CC39 (6.54)25 (4.04)
CTLA-4rs231804 (A>G)1.00
A1031 (86.35)1036 (83.55)3.720.054
G163 (13.65)204 (16.45)
AA447 (74.87)432 (69.68)4.100.129
AG137 (22.95)172 (27.74)
GG13 (2.18)16 (2.58)
rs1024161 (A>G)0.44
A876 (73.37)868 (70.00)3.400.065
G318 (26.63)372 (30.00)
AA326 (54.61)308 (49.68)3.300.192
AG224 (37.52)252 (40.65)
GG47 (7.87)60 (9.67)
rs231726 (A>G)0.79
A813 (68.09)794 (64.03)4.470.035
G381 (31.91)446 (35.97)
AA283 (47.40)256 (41.29)4.750.093
AG247 (41.37)282 (45.48)
GG67 (11.23)82 (13.23)
rs10197319 (C>T)0.45
C937 (78.48)938 (75.65)2.760.097
T257 (21.52)302 (24.35)
CC367 (61.48)351 (56.61)3.000.223
CT203 (34.00)236 (38.07)
TT27 (4.52)33 (5.32)

HWE: Hardy-Weinberg equilibrium; .

In the TSHR gene, there was no difference in the distribution of alleles or genotypes at rs2300525 between the case group and the control group (). For rs179247, rs2284722, rs12101261, rs4903964, and rs17111394, the distribution of the corresponding alleles or genotypes between the two groups was statistically significant ().

In the CTLA-4 gene, there was no difference in the distribution of alleles or genotypes at sites rs231804, rs1024161, or rs10197319 between the case group and the control group (). There was no significant difference in the distribution of the genotype at rs231726 between the two groups (), while the distribution difference of alleles at this site between the two groups was statistically significant ().

After adjusting for age, a correlation analysis of the association between these ten SNPs and GD revealed that rs179247-G, rs12101261-C, and rs4903964-G were negatively correlated with the incidence of GD in both the male and female populations. In addition, rs2284722-A and rs17111394-C were positively correlated with the incidence of GD in the overall and female populations, while rs2300525-C was positively correlated with the incidence of GD only in the female population. CTLA-4 rs231726-G was negatively correlated with the risk of GD in the overall population. In the dominant model of rs10197319, TT+CT carriers had a higher risk of GD in the overall and male populations than those with the homozygous CC genotype (Table 3).


SNPsAllelesDominant modelRecessive modelHomozygous modelHeterozygous model
OR (95% CI)OR (95% CI)OR (95% CI)OR (95% CI)OR (95% CI)

TSHRrs179247G/AGG+AG/AAGG/AG+AAGG/AAAG/AA
TP0.57 (0.48, 0.68)<0.0010.50 (0.40, 0.63)<0.0010.50 (0.34, 0.73)<0.0010.37 (0.25, 0.55)<0.0010.54 (0.42, 0.68)<0.001
SexM0.56 (0.38, 0.84)0.0050.57 (0.34, 0.95)0.0300.30 (0.11, 0.80)0.0160.25 (0.09, 0.69)0.0070.67 (0.39, 1.15)0.148
F0.57 (0.47, 0.70)<0.0010.48 (0.37, 0.62)<0.0010.55 (0.36, 0.83)0.0050.40 (0.26, 0.62)<0.0010.51 (0.39, 0.66)<0.001
rs2284722A/GAA+AG/GGAA/AG+GGAA/GGAG/GG
TP1.34 (1.12, 1.60)0.0011.33 (1.06, 1.67)0.0132.14 (1.34, 3.42)0.0022.33 (1.44, 3.78)0.0011.22 (0.96, 1.55)0.097
SexM0.91 (0.61, 1.38)0.6690.90 (0.54, 1.51)0.6940.88 (0.34, 2.30)0.7970.85 (0.32, 2.28)0.7480.91 (0.53, 1.57)0.744
F1.47 (1.21, 1.80)<0.0011.46 (1.14, 1.88)0.0032.79 (1.61, 4.83)<0.0013.15 (1.79, 5.53)<0.0011.31 (1.01, 1.70)0.046
rs12101261C/TCC+CT/TTCC/CT+TTCC/TTCT/TT
TP0.58 (0.49, 0.69)<0.0010.53 (0.42, 0.66)<0.0010.44 (0.30, 0.63)<0.0010.34 (0.23, 0.50)<0.0010.59 (0.46, 0.75)<0.001
SexM0.58 (0.40, 0.86)0.0060.54 (0.32, 0.90)0.0180.43 (0.18, 1.01)0.0530.34 (0.14, 0.82)0.0170.60 (0.35, 1.04)0.066
F0.58 (0.48, 0.70)<0.0010.52 (0.41, 0.68)<0.0010.44 (0.29, 0.66)<0.0010.34 (0.22, 0.52)<0.0010.59 (0.45, 0.77)<0.001
rs4903964G/AGG+AG/AAGG/AG+AAGG/AAAG/AA
TP0.58 (0.50, 0.69)<0.0010.57 (0.45, 0.71)<0.0010.38 (0.27, 0.54)<0.0010.31 (0.21, 0.45)<0.0010.67 (0.53, 0.86)0.020
SexM0.60 (0.41, 0.87)0.0080.56 (0.33, 0.94)0.0270.46 (0.21, 0.98)0.0440.36 (0.16, 0.82)0.0140.64 (0.37, 1.10)0.107
F0.58 (0.48, 0.70)<0.0010.57 (0.44, 0.74)<0.0010.37 (0.25, 0.54)<0.0010.30 (0.20, 0.45)<0.0010.68 (0.52, 0.90)0.007
rs2300525C/TCC+CT/TTCC/CT+TTCC/TTCT/TT
TP1.18 (0.99, 1.40)0.0721.18 (0.95, 1.48)0.1431.34 (0.90, 2.01)0.1541.42 (0.93, 2.15)0.1011.14 (0.90, 1.44)0.288
SexM0.92 (0.62, 1.36)0.6680.87 (0.53, 1.45)0.6000.98 (0.40, 2.39)0.9590.91 (0.36, 2.32)0.8490.87 (0.51, 1.47)0.593
F1.25 (1.03, 1.52)0.0281.28 (0.99, 1.64)0.0581.45 (0.92, 2.29)0.1061.58 (0.99, 2.52)0.0551.22 (0.93, 1.59)0.146
rs17111394C/TCC+CT/TTCC/CT+TTCC/TTCT/TT
TP1.29 (1.06, 1.57)0.0101.29 (1.02, 1.63)0.0311.66 (0.99, 2.78)0.0551.78 (1.05, 3.00)0.0311.22 (0.96, 1.56)0.111
SexM1.02 (0.66, 1.57)0.9341.06 (0.63, 1.78)0.8340.88 (0.30, 2.60)0.8220.91 (0.30, 2.72)0.8651.06 (0.61, 1.82)0.844
F1.37 (1.10, 1.70)0.0051.36 (1.05, 1.76)0.0211.99 (1.10, 3.59)0.0232.15 (1.18, 3.91)0.0121.26 (0.96, 1.66)0.092
CTLA-4rs231804G/AGG+AG/AAGG/AG+AAGG/AAAG/AA
TP0.80 (0.64, 1.00)0.0540.77 (0.60, 1.00)0.0450.86 (0.41, 1.80)0.6800.80 (0.38, 1.86)0.5560.77 (0.59, 1.00)0.050
SexM0.66 (0.39, 1.11)0.1160.69 (0.39, 1.23)0.2030.75 (0.42, 1.36)0.346
F0.84 (0.66, 1.08)0.1760.80 (0.60, 1.05)0.1081.08 (0.50, 2.36)0.8451.01 (0.46, 2.21)0.9770.78 (0.58, 1.04)0.086
rs1024161G/AGG+AG/AAGG/AG+AAGG/AAAG/AA
TP0.85 (0.71, 1.01)0.0650.82 (0.66, 1.03)0.0850.81 (0.54, 1.20)0.2930.75 (0.49, 1.13)0.1620.84 (0.66, 1.06)0.145
SexM0.93 (0.62, 1.39)0.7120.94 (0.57, 1.57)0.8250.77 (0.27, 2.18)0.6170.76 (0.26, 2.22)0.6120.98 (0.58, 1.65)0.929
F0.83 (0.68, 1.01)0.0650.79 (0.62, 1.02)0.0700.81 (0.53, 1.26)0.3510.74 (0.48, 1.16)0.1920.81 (0.62, 1.05)0.112
rs231726G/AGG+AG/AAGG/AG+AAGG/AAAG/AA
TP0.83 (0.71, 0.99)0.0350.78 (0.62, 0.98)0.0310.83 (0.59, 1.18)0.3000.74 (0.52, 1.07)0.1080.79 (0.62, 1.01)0.055
SexM0.76 (0.52, 1.12)0.1630.74 (0.44, 1.23)0.2480.62 (0.27, 1.41)0.2540.55 (0.23, 1.31)0.1770.79 (0.46, 1.36)0.398
F0.85 (0.71, 1.03)0.0970.79 (0.61, 1.02)0.0670.89 (0.61, 1.30)0.5410.79 (0.53, 1.18)0.2520.79 (0.60, 1.03)0.084
rs10197319T/CTT+CT/CCTT/CT+CCTT/CCCT/CC
TP0.85 (0.70, 1.03)0.0970.68 (0.54, 0.86)0.0010.85 (0.51, 1.43)0.5430.75 (0.44, 1.28)0.2910.67 (0.53, 0.86)0.001
SexM0.94 (0.62, 1.45)0.7930.32 (0.18, 0.57)<0.0011.12 (0.29, 4.27)0.8720.78 (0.20, 3.03)0.7250.28 (0.15, 0.52)<0.001
F0.83 (0.67, 1.03)0.0850.80 (0.62, 1.03)0.0800.81 (0.46, 1.43)0.4660.75 (0.42, 1.33)0.3190.80 (0.62, 1.05)0.107

TP: total population; M: male; F: female; .

The linkage disequilibrium analysis of the tested polymorphisms in TSHR and CTLA-4 revealed that four TSHR SNPs rs179247, rs2284722, rs12101261, and rs4903964, two SNPs (rs2300525 and rs17111394), and three SNPs in CTLA-4 (rs231804, rs1024161, and rs231726), respectively, formed three haplotype regions, each with a linkage disequilibrium coefficient (Figures 1 and 2 and Table 4; ).


rs179247rs2284722rs12101261rs4903964

rs1792470.8450.9550.862
rs22847220.1210.850.793
rs121012610.8320.1340.881
rs49039640.5620.1410.645
rs2300525rs17111394
rs23005250.977
rs171113940.657
rs231804rs1024161rs231726
rs2318040.9870.99
rs10241610.4370.993
rs2317260.3380.758

The italic part represents , and the non-italic part represents .

Eight different haplotypes were identified based on the rs179247, rs2284722, rs12101261, and rs4903964 sites in the TSHR gene. The TSHR sites rs2300525 and rs17111394 composed three different haplotypes, while four different haplotypes were identified based on the rs231804, rs1024161, and rs231726 alleles of CTLA-4. The relationship between these haplotypes and GD revealed that haplotypes AGTA, GGCG, and AATA, haplotype CC, and haplotype AAA, respectively, were associated with the risk of GD (Table 5). No correlation was observed between TSHR and CTLA-4 genotypes and clinical phenotypes (Table 6).


SNPHaplotypeCase groupControl groupOR (95% CI)
%%

rs179247AGTA22437.620032.21.27 (1.07, 1.50)0.005
rs2284722GGCG12721.320332.70.56 (0.46, 0.67)<0.001
rs12101261AATA16828.213221.31.45 (1.21, 1.75)<0.001
rs4903964AGTG325.3335.31.00 (0.70, 1.43)0.997
AGCG142.4162.60.94 (0.57, 1.57)0.820
GGCA111.9101.61.20 (0.66, 2.20)0.548
AATG61.081.30.76 (0.36, 1.64)0.488
GACG40.691.50.40 (0.18, 0.92)0.031
rs2300525TT41769.845172.80.86 (0.73, 1.03)0.105
rs17111394CC13923.311718.81.32 (1.08, 1.60)0.006
CT396.6497.90.83 (0.61, 1.12)0.223
rs231804AAA40567.939563.71.21 (1.02, 1.43)0.029
rs1024161GGG8113.510216.40.80 (0.64, 1.00)0.046
rs231726AGG7913.28313.40.98 (0.78, 1.24)0.863
AAG325.3386.20.84 (0.60, 1.19)0.323

Four of the eight haplotypes (composed of rs179247, rs2284722, rs12101261, and rs4903964) had frequencies greater than 0.03. There are four haplotypes (composed of rs231804, rs1024161, and rs231726) that had frequencies greater than 0.03. .

GoiterExophthalmosTRAb (U/L)
SNPGenotype0IIIIIIYesNo≥1.5<1.5

TSHRrs179247AA79923243.670.301501921.110.57274474.230.12
AG77112768312815340
GG1162701925299
rs2284722GG98119741.100.781191720.410.82215532.720.26
AG687153510914219437
AA01836552431476
rs12101261TT710320463.140.371361840.170.92260434.840.09
CT66815449613616445
CC2152802025328
rs4903964AA78617972.610.461211580.320.85226365.600.06
AG583175310915719253
GG3173202230387
rs2300525TT89019543.480.321291681.620.45227490.110.95
CT576155610313918437
CC22036020384510
rs17111394TT1011123042.870.411532022.240.33267634.420.11
CT45813568711616124
CC1172101227289
CTLA-4rs231804AA1213828891.540.671872601.030.60345720.500.78
AG345881617610221
GG031004993
rs1024161AA1110120772.560.461401860.780.68250531.510.47
AG26915039512917433
GG21629017303210
rs231726AA98618440.700.871231600.840.66218460.840.66
AG481157510414319037
GG21945125424813
rs10197319CC1210823984.780.191572100.710.70282580.140.93
CT36912928212115334
TT091801314214

Using MDR analysis of TSHR and CTLA-4 multiple site interaction yielded a significance of , indicating that there is no interaction between TSHR and CTLA-4 polymorphisms (Table 7). False-positive report probabilities (FPRPs) for the identified SNPs significantly associated with GD were <0.2 for rs179247, rs12101261, and rs4903964 (Table 8).


ModelTraining set balance precisionTest set balance precisionCross-validation consistencyOR (95% CI)

rs1792470.58640.586410/100.05652.01 (0.98, 4.14)
rs179247, rs49039640.59190.57346/100.10531.81 (0.88, 3.71)
rs179247, rs4903964, rs101973190.60740.57065/100.11931.77 (0.86, 3.62)


SNPOR (95% CI)Prior probability
0.250.10.010.001

rs1792470.57 (0.48, 0.68)<0.001<0.001<0.001<0.001
rs22847221.34 (1.12, 1.60)0.0040.0120.1190.576
rs121012610.58 (0.49, 0.69)<0.001<0.001<0.001<0.001
rs49039640.58 (0.50, 0.69)<0.001<0.001<0.001<0.001
rs171113941.29 (1.06, 1.57)0.0340.0960.5400.922
rs2317260.83 (0.71, 0.99)0.1040.2580.7920.975

Preset the critical ; .

4. Discussion

GD is a common autoimmune thyroid disease. The prevalence of clinical hyperthyroidism in China is about 0.8%, and 80% of these cases are the result of GD. GD develops as a consequence of complex interactions between genetic, environmental, and immunological factors. TRAb is the most frequently encountered autoantibodies in those with GD (present in >90% of patients), and it can compete with TSH for TSHR binding. TSAb is a pathogenic antibody associated with GD, and recent work has shown that TSAb is linked with oxidative stress present in those with GD [28]. Oxidative stress is associated with GD inflammation. These TSHR agonist antibodies that trigger TSHR signaling are characteristic of GD. Genetic polymorphisms both in TSHR and in immune genes that regulate central and peripheral tolerance, such as CTLA-4 which constrains T cell activation, are linked to the pathogenesis of GD. As such, we assessed SNPs in both TSHR and CTLA-4 in the present study.

We confirmed that all 6 TSHR SNPS and all 4 CTLA-4 SNPs conformed to the Hardy-Weinberg equilibrium, indicating that our research subjects were a good representative population ().

We next analyzed the correlation of the six studied TSHR SNPs with GD, and through different analysis models, we identified the risk alleles and genotypes for GD associated with these sites.

For TSHR site rs179247, GD patients had a higher frequency of the A allele than healthy controls, and this allele was found to be the primary risk factor. In contrast, the G allele was negatively correlated with GD in the total population and in the individual male and female populations. Bufalo et al. demonstrated that TSHR intronic polymorphisms are associated with GD and Graves’ ophthalmopathy susceptibility in a Brazilian population, with the AA genotype for rs179247 increasing GD risk [29]. rs179247 SNP AA or AG individuals have significantly lower TSHR mRNA expression levels in the thymus in nonautoimmune donors relative to those with the GG phenotype [30]. An analysis of polymorphisms in TSHR rs179247 pertaining to the pathogenesis of autoimmune thyroid diseases in children found that A alleles were more frequent in patients with GD in comparison to healthy subjects and that polymorphisms at this site contributed to the development of AITDs in children [21].

Several meta-analyses of this site have reported that it is associated with GD [3133], and our study results are consistent with these reports.

In our study, the main GD risk factor at site rs2284722 was the A allele; A alleles were positively correlated with GD incidence in both the total population and the female population, and the AA genotype was associated with a higher risk for GD, especially in females.

For TSHR site rs12101261, the C alleles were negatively correlated with the incidence of GD in the total population, the male population, and the female population. The GD main risk factor for rs12101261 was the T alleles. There have been reports identifying this site as a possible causal SNP for GD susceptibility in the TSHR gene, potentially serving as a genetic marker to predict the outcome of persistent TSHR autoantibody positivity in GD patients [34]. Research has shown that the rs12101261 disease-associated T alleles are associated with lower TSHR expression in the thymus [16].

The main GD risk factors for rs4903964 were the A alleles; the G alleles were negatively correlated with the incidence of GD in the total population, the male population, and the female population.

There were no difference in the distribution of alleles or genotypes at rs2300525 between the case group and the control group (). We further analyzed the correlation of TSHR rs2300525 alleles with GD, and we found that the C allele was positively correlated with GD incidence only in the female population and that the risk of GD was 1.25x for C alleles relative to T alleles.

For TSHR site rs17111394, the C alleles were positively correlated with the incidence of GD in the total population and the female population. The risk of GD in mutant homozygous carriers was 1.78x and 2.15x, respectively, relative to homozygous TT carriers.

Linkage disequilibrium is a measure of the correlation between alleles at different loci. The linkage disequilibrium analysis of the four TSHR SNP sites rs179247, rs2284722, rs12101261, and rs4903964 and the linkage disequilibrium coefficient were , indicating that these four SNP sites were located in the same haplotype region. Two polymorphism sites, rs2300525 and rs17111394, have an imbalance coefficient of , indicating that these two SNP sites are located in the same haplotype region.

The TSHR SNPs rs179247, rs2284722, rs12101261, and rs4903964 together formed 8 distinct haplotypes. We analyzed the relationship between these haplotypes and GD. Haplotypes AGTA, GGCG, and AATA were the most frequent of all haplotypes, accounting for more than 80 percent of the entire population. There was a positive correlation between GD risk and haplotypes AGTA and AATA. The risk was negatively correlated with the haplotypes GGCG and GACG; however, the GACG haplotype was a very small percentage of the study population. No association with GD was found for the other haplotypes (). For the rs2300525 and rs17111394 haplotypes, there was a positive correlation between haplotype CC and GD risk.

As for CTLA-4, there was no difference in the distribution of genotypes at four sites on the CTLA-4 gene between the case group and the control group (), while the distribution of alleles at rs231726 between the two groups was statistically significant (). CTLA-4 rs231726-G alleles were negatively correlated with the risk of GD in the total population, and the risk of GD was 0.73x for G allele carriers relative to A allele carriers. In the dominant model, the risk of GD in the total population was 0.78x for GG+AG carriers relative to homozygous AA carriers. In the dominant model of the rs10197319 site, TT+CT carriers had a 0.68x risk of GD in the total population relative to homozygous CC carriers, and the risk of GD in the male population was 0.32x relative to homozygous CC carriers. No correlation was found between the rs231804 and rs1024161 sites and GD risk. The linkage disequilibrium coefficients of three polymorphism sites, namely, rs231804, rs1024161, and rs231726, , indicated that these 3 SNP sites were located in the same haplotype region.

We further analyzed the relationship between GD and the haplotypes formed by CTLA-4 SNP sites rs231804, rs1024161, and rs231726. A positive correlation was identified between haplotype AAA and GD risk (). No association with GD was found for the other haplotypes ().

We analyzed comparisons between TSHR and CTLA-4 genotypes and clinical GD characteristics such as diffuse goiter, exophthalmos, and different levels of TRAb, but we did not identify any correlations between these variables (). It is important to note that many patients were being treated using antithyroid drugs which have the potential to alter TRAb/TSH levels and to impact these clinical characteristics.

GD is a complex disease triggered by multiple genes and other factors. The information provided by a single SNP site is thus very limited in the study of such complex diseases. The variation and the interaction between multiple SNP sites are what ultimately lead to the development of such complex diseases. The pathogenesis of GD is related to the interaction between genes and the environment [35, 36]. We therefore analyzed the interactions between TSHR and CTLA-4 polymorphisms to assess their impact on GD. The best cross-validation consistency was for the single-factor model rs179247, whose training balance precision was 0.5864, the test set balance precision is 0.5864, and the cross-validation consistency rate is 10/10, but . Therefore, no interaction between TSHR gene and CTLA-4 gene was found.

We calculated the FPRP for the TSHR and CTLA-4 SNPS significantly associated with GD, yielding for rs179247, rs12101261, and rs4903964, indicating that the association between these genotypes and GD was significant.

In conclusion, genetic factors play an important role in the development of GD. The five studied SNPs in TSHR intron 1 and the linkage disequilibrium haplotype composed of those related loci were all associated with GD. Only one linkage disequilibrium haplotype (AAA, composed of CTLA-4 sites rs231804, rs1024161, and rs231726) was found to be related to GD.

Previously, functional SNPs associated with diseases have been found to be more concentrated in the regulatory or coding regions of the genome. In contrast, most SNPs in the noncoding regions of the genome have been ignored by researchers. Most detected SNPs are normal variants and do not have biological functions. In spite of the fact that these SNPs function weakly or are even nonfunctional, their combinations often show a good correlation with disease or related phenotypes. To detect different specific SNP combinations associated with disease groups and to analyze the combination of these SNPs, it is possible to determine the indicators of disease sensitivity and predict therapeutic efficacy with the goal of improving early treatment and prevention. There have been many reports regarding the relationship between TSHR and CTLA-4 and GD in different populations [17, 19, 29, 37], including the Chinese Han population, revealing polymorphisms in these two genes to be related to GD. While these studies assessed different sites and populations, they all employed case-control methodology, as we did in the present study. We also had the advantage of assessing more SNPs than previous studies and of using different statistical models to achieve a more in-depth statistical analysis. Even so, our results were limited by the fact that our study included medicated patients, potentially impacting the association between identified genotypes and clinical presentation. While our results are promising, any statistical tests have the potential for false-positive results or correlations due to multiple confounding variables such as sample size or population stratification. Therefore, these results need to be further verified in a larger sample and among different populations.

Data Availability

The genotype data used to support the findings of this study are restricted by the Medical Ethics Committee of The First Affiliated Hospital of Bengbu Medical College in order to protect patient privacy. Data are available from the author sunwh2007@163.com for researchers who meet the criteria for access to confidential data.

Conflicts of Interest

No competing interests exist.

Authors’ Contributions

Weihua Sun and Xiaomei Zhang contributed equally to this work.

Acknowledgments

We are grateful to our patients and healthy volunteers for donating blood for our study as well as for their agreement for using their data for this study. This work was supported by the Anhui Provincial Department of Education Key Projects (KJ2013A187).

References

  1. S. K. Shukla, G. Singh, S. Ahmad, and P. Pant, “Infections, genetic and environmental factors in pathogenesis of autoimmune thyroid diseases,” Microbial Pathogenesis, vol. 116, pp. 279–288, 2018. View at: Publisher Site | Google Scholar
  2. M. Dittmar, C. Libich, T. Brenzel, and G. J. Kahaly, “Increased familial clustering of autoimmune thyroid diseases,” Hormone and Metabolic Research, vol. 43, no. 3, pp. 200–204, 2011. View at: Publisher Site | Google Scholar
  3. K. Hemminki, X. Li, J. Sundquist, and K. Sundquist, “The epidemiology of Graves disease: evidence of a genetic and an environmental contribution,” Journal of Autoimmunity, vol. 34, no. 3, pp. J307–J313, 2010. View at: Publisher Site | Google Scholar
  4. D. S. A. McLeod, D. S. Cooper, P. W. Ladenson, D. C. Whiteman, and S. J. Jordan, “Race/ethnicity and the prevalence of thyrotoxicosis in young Americans,” Thyroid, vol. 25, no. 6, pp. 621–628, 2015. View at: Publisher Site | Google Scholar
  5. X. Chu, C. M. Pan, S. X. Zhao et al., “A genome-wide association study identifies two new risk loci for Graves' disease,” Nature Genetics, vol. 43, no. 9, pp. 897–901, 2011. View at: Publisher Site | Google Scholar
  6. J. D. Cooper, M. J. Simmonds, N. M. Walker et al., “Seven newly identified loci for autoimmune thyroid disease,” Human Molecular Genetics, vol. 21, no. 23, pp. 5202–5208, 2012. View at: Publisher Site | Google Scholar
  7. N. Eriksson, J. Y. Tung, A. K. Kiefer et al., “Novel associations for hypothyroidism include known autoimmune risk loci,” PLoS One, vol. 7, no. 4, article e34442, 2012. View at: Publisher Site | Google Scholar
  8. S. X. Zhao, L. Q. Xue, W. Liu et al., “Robust evidence for five new Graves disease risk loci from a staged genome-wide association analysis,” Human Molecular Genetics, vol. 22, no. 16, pp. 3347–3362, 2013. View at: Publisher Site | Google Scholar
  9. X. Chu, M. Shen, F. Xie et al., “An X chromosome-wide association analysis identifies variants in GPR174 as a risk factor for Graves disease,” Journal of Medical Genetics, vol. 50, no. 7, pp. 479–485, 2013. View at: Publisher Site | Google Scholar
  10. S. H. Kwak, Y. J. Park, M. J. Go et al., “A genome-wide association study on thyroid function and anti-thyroid peroxidase antibodies in Koreans,” Human Molecular Genetics, vol. 23, no. 16, pp. 4433–4442, 2014. View at: Publisher Site | Google Scholar
  11. M. Medici, E. Porcu, G. Pistis et al., “Identification of novel genetic loci associated with thyroid peroxidase antibodies and clinical thyroid disease,” PLoS Genetics, vol. 10, no. 2, article e1004123, 2014. View at: Publisher Site | Google Scholar
  12. D. Oryoji, S. Ueda, K. Yamamoto et al., “Identification of a Hashimoto thyroiditis susceptibility locus via a genome-wide comparison with Graves disease,” The Journal of Clinical Endocrinology and Metabolism, vol. 100, no. 2, pp. E319–E324, 2015. View at: Publisher Site | Google Scholar
  13. L. Brčić, A. Barić, S. Gračan et al., “Genome-wide association analysis suggests novel loci for Hashimoto’s thyroiditis,” Journal of Endocrinological Investigation, vol. 42, no. 5, pp. 567–576, 2019. View at: Publisher Site | Google Scholar
  14. M. Zhan, G. Chen, C. M. Pan et al., “Genome-wide association study identifies a novel susceptibility gene for serum TSH levels in Chinese populations,” Human Molecular Genetics, vol. 23, no. 20, pp. 5505–5517, 2014. View at: Publisher Site | Google Scholar
  15. L. Jansson, K. Vrolix, A. Jahraus, K. F. Martin, and D. C. Wraith, “Immunotherapy with apitopes blocks the immune response to TSH receptor in HLA-DR transgenic mice,” Endocrinology, vol. 159, no. 9, pp. 3446–3457, 2018. View at: Publisher Site | Google Scholar
  16. M. Stefan, C. Wei, A. Lombardi et al., “Genetic-epigenetic dysregulation of thymic TSH receptor gene expression triggers thyroid autoimmunity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 111, no. 34, pp. 12562–12567, 2014. View at: Publisher Site | Google Scholar
  17. S. X. Zhao, C. M. Pan, H. M. Cao et al., “Association of the CTLA-4 gene with Graves disease in the Chinese Han population,” PLoS One, vol. 5, no. 3, article e9821, 2010. View at: Publisher Site | Google Scholar
  18. W. H. Ting, M. N. Chien, F. S. Lo et al., “Association of cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) gene polymorphisms with autoimmune thyroid disease in children and adults: case-control study,” PLoS One, vol. 11, no. 4, article e0154394, 2016. View at: Publisher Site | Google Scholar
  19. X. Chen, Z. Hu, M. Liu et al., “Correlation between CTLA-4 and CD40 gene polymorphisms and their interaction in Graves disease in a Chinese Han population,” BMC Medical Genetics, vol. 19, no. 1, p. 171, 2018. View at: Publisher Site | Google Scholar
  20. S. Ramgopal, C. Rathika, R. Padma Malini, V. Murali, K. Arun, and K. Balakrishnan, “Critical amino acid variations in HLA-DQB1 molecules confers susceptibility to autoimmune thyroid disease in South India,” Genes & Immunity, vol. 20, no. 1, pp. 32–38, 2019. View at: Publisher Site | Google Scholar
  21. M. Rydzewska, A. Góralczyk, J. Gościk 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, vol. 51, no. 4, pp. 183–190, 2018. View at: Publisher Site | Google Scholar
  22. F. Shehjar, D. Afroze, R. A. Misgar, S. A. Malik, and B. A. Laway, “Association of FoxP3 promoter polymorphisms with the risk of Graves' disease in ethnic Kashmiri population,” Gene, vol. 672, pp. 88–92, 2018. View at: Publisher Site | Google Scholar
  23. L. Li, X. Ding, X. Wang et al., “Polymorphisms of IKZF3 gene and autoimmune thyroid diseases: associated with Graves disease but not with Hashimoto’s thyroiditis,” Cellular Physiology and Biochemistry, vol. 45, no. 5, pp. 1787–1796, 2018. View at: Publisher Site | Google Scholar
  24. D. Wang, J. Chen, H. Zhang, F. Zhang, L. Yang, and Y. Mou, “Role of different CD40 polymorphisms in Graves disease and Hashimoto’s thyroiditis,” Immunological Investigations, vol. 46, no. 6, pp. 544–551, 2017. View at: Publisher Site | Google Scholar
  25. J. J. Khong, K. P. Burdon, Y. Lu et al., “Pooled genome wide association detects association upstream of FCRL3 with Graves’ disease,” BMC Genomics, vol. 17, no. 1, p. 939, 2016. View at: Publisher Site | Google Scholar
  26. A. Lombardi, F. Menconi, D. Greenberg et al., “Dissecting the genetic susceptibility to Graves’ disease in a cohort of patients of Italian origin,” Frontiers in Endocrinology, vol. 7, p. 21, 2016. View at: Publisher Site | Google Scholar
  27. L. Q. Gu, W. Zhu, S. X. Zhao et al., “Clinical associations of the genetic variants of CTLA-4, Tg, TSHR, PTPN22, PTPN12 and FCRL3 in patients with Graves’ disease,” Clinical Endocrinology, vol. 72, no. 2, pp. 248–255, 2010. View at: Publisher Site | Google Scholar
  28. T. Diana, A. Daiber, M. Oelze et al., “Stimulatory TSH-receptor antibodies and oxidative stress in Graves disease,” The Journal of Clinical Endocrinology & Metabolism, vol. 103, no. 10, pp. 3668–3677, 2018. View at: Publisher Site | Google Scholar
  29. N. E. Bufalo, R. B. dos Santos, M. A. Marcello et al., “TSHR intronic polymorphisms (rs179247 and rs12885526) and their role in the susceptibility of the Brazilian population to Graves’ disease and Graves’ ophthalmopathy,” Journal of Endocrinological Investigation, vol. 38, no. 5, pp. 555–561, 2015. View at: Publisher Site | Google Scholar
  30. R. Colobran, M. . P. Armengol, R. Faner et al., “Association of an SNP with intrathymic transcription of TSHR and Graves' disease: a role for defective thymic tolerance,” Human Molecular Genetics, vol. 20, no. 17, pp. 3415–3423, 2011. View at: Publisher Site | Google Scholar
  31. H. Xiong, M. Wu, H. Yi et al., “Genetic associations of the thyroid stimulating hormone receptor gene with Graves diseases and Graves ophthalmopathy: a meta-analysis,” Scientific Reports, vol. 6, no. 1, article 30356, 2016. View at: Publisher Site | Google Scholar
  32. W. Qian, K. Xu, W. Jia et al., “Association between TSHR gene polymorphism and the risk of Graves' disease: a meta-analysis,” Journal of Biomedical Research, vol. 30, no. 6, pp. 466–475, 2016. View at: Publisher Site | Google Scholar
  33. J. Gong, S. J. Jiang, D. K. Wang et al., “Association of polymorphisms of rs179247 and rs12101255 in thyroid stimulating hormone receptor intron 1 with an increased risk of Graves disease: a meta-analysis,” Journal of Huazhong University of Science and Technology. Medical Sciences, vol. 36, no. 4, pp. 473–479, 2016. View at: Publisher Site | Google Scholar
  34. B. L. Liu, S. Y. Yang, W. Liu et al., “Refined association of TSH receptor susceptibility locus to Graves disease in the Chinese Han population,” European Journal of Endocrinology, vol. 170, no. 1, pp. 109–119, 2014. View at: Publisher Site | Google Scholar
  35. M. J. Simmonds, “GWAS in autoimmune thyroid disease: redefining our understanding of pathogenesis,” Nature Reviews Endocrinology, vol. 9, no. 5, pp. 277–287, 2013. View at: Publisher Site | Google Scholar
  36. Y. H. Dong and D. G. Fu, “Autoimmune thyroid disease: mechanism, genetics and current knowledge,” European Review for Medical and Pharmacological Sciences, vol. 18, no. 23, pp. 3611–3618, 2014. View at: Google Scholar
  37. L. Liu, H.-Q. Wu, Q. Wang et al., “Association between thyroid stimulating hormone receptor gene intron polymorphisms and autoimmune thyroid disease in a Chinese Han population,” Endocrine Journal, vol. 59, no. 8, pp. 717–723, 2012. View at: Publisher Site | Google Scholar

Copyright © 2019 Weihua Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

478 Views | 281 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.