Skip to main content

Polymorphisms in the ADRB2 gene and Graves disease: a case-control study and a meta-analysis of available evidence



The beta-2-Adrenergic receptor (ADRB2) gene on chromosome 5q33.1 is an important immunoregulatory factor. We and others have previously implicated chromosomal region 5q31-33 for contribution to the genetic susceptibility to Graves disease (GD) in East-Asian populations. Two recent studies showed associations between the single nucleotide polymorphism (SNP) rs1042714 in the ADRB2 gene and GD. In this study, we aimed to fully investigate whether the ADRB2 gene conferred susceptibility to GD in Chinese population, and to perform a meta-analysis of association between ADRB2 and GD.


Approximately 1 kb upstream the transcription start site and the entire coding regions of the ADRB2 gene were resequenced in 48 Han Chinese individuals to determine the linkage disequilibrium (LD) patterns. Tag SNPs were selected and genotyped in a case-control collection of 1,118 South Han Chinese subjects, which included 428 GD patients and 690 control subjects. A meta-analysis was performed with the data obtained in the present samples and those available from prior studies.


Fifteen SNPs in the ADRB2 gene were identified by resequencing and one SNP was novel. Ten tag SNPs were investigated further to assess association of ADRB2 in the case-control collection. Neither individual tag SNP nor haplotypes showed association with GD in Han Chinese population (P > 0.05). Our meta-analysis of the ADRB2 SNP rs1042714 measured heterogeneity between the ethnic groups (I2 = 53.1%) and no association to GD was observed in the overall three studies with a random effects model (OR = 1.13, 95% CI, 0.95 to 1.36; P = 0.18). However, significant association was found from the combined data of Caucasian population with a fixed effects model (OR = 1.18, 95% CI, 1.06 to 1.32; P = 0.002; I2 = 5.9%).


Our study indicated that the ADRB2 gene did not exert a substantial influence on GD susceptibility in Han Chinese population, but contributed to a detectable GD risk in Caucasian population. This inconsistency resulted largely from between-ethnicity heterogeneity.

Peer Review reports


Graves disease (GD) affects 0.5–1% of the general population[1], and results from the presence of autoantibodies to the thyroid-stimulating hormone receptor (TSHR), leading to over-activity of the thyroid gland. Genome-wide screens conducted in GD have identified a number of regions of linkage, although independent replication and novel candidate loci are to be awaited. Our previous study[2], along with others[3, 4], has suggested that chromosomal region 5q31-33 may contain a locus that contributes to the genetic susceptibility to GD in Eastern-Asian populations. ADRB2, a 1242-base pair intronless gene, is located on the long arm of chromosome 5q33.1. At least thirteen SNPs have been identified in the gene region[5, 6]. Two activity-related polymorphisms (rs1042713 and rs1042714) with high allelic frequency in the general population are single nucleic acid substitutions at positions 46 (A-G) and 79 (C-G), corresponding to substitutions of glycine for arginine at amino acid position 16 and glutamate for glutamine at amino acid position 27[5, 6].

Recently, these two SNPs (rs1042713 and rs1042714) and another SNP located at position -367 (rs11959427) were reported association with GD in a Polish case-control collection[7]. The result indicated that Gln27 carriers (rs1042714CC or rs1042714GC genotypes) had increased risk of GD. Additionally, the genome-wide association study (GWAS) of the Wellcome Trust Case Control Consortium (WTCCC) et al also showed positive signal in the ADRB2 gene[8]. Although three SNPs (rs12654778, rs1042713 and rs1042714) in the ADRB2 gene genotyped in the study of WTCCC et al. did not meet a point-wise significance level of P < 10-3 for the Cochran-Armitage test for trend, the original analyses showed significant association of rs1042713 and rs1042714 to GD (P = 0.04 for rs1042713 and P = 0.02 for rs1042714).

ADRB2 is a G protein coupled receptor that is expressed by all lymphoid cells, with the exception of T helper (Th) 2 cells[9]. Norepinephrine and epinephrine, through stimulation of the ADRB2-cAMP-protein kinase A pathway, cause a selective suppression of Th1 responses and cellular immunity, and a Th2 shift toward dominance of humoral immunity[9]. Moreover, various polymorphisms in the ADRB2 gene were demonstrated as leading to a variation of ADRB2 expression on the cell surface[10, 11]. Therefore, the polymorphisms may influence the activity of ADRB2, and disturb the Th1/Th2 cytokine balance, which is an important factor in the pathogenesis of GD. In addition, ADRB2 had been reported association with another autoimmune disease, rheumatoid arthritis [12]. This increased the possibility of ADRB2 association with GD, as different autoimmune diseases may share similar genetic etiologies. Given its chromosomal location within a region of genetic linkage to GD and the role of ADRB2 on the immune response, the ADRB2 gene makes a plausible candidate for GD and autoimmunity in general.

In the present study, we set to determine whether ADRB2 was a susceptibility locus for GD in Han Chinese population and conduct a global meta-analysis of the ADRB2 variation association with GD. Firstly, we identified the variants of ADRB2 by resequencing the 5' flanking region and the entire coding region; then we defined linkage disequilibrium (LD) among SNPs across the ADRB2 genomic interval, and selected tag SNPs. Secondly, we conducted a case-control association study (428 GD cases and 690 controls) using the tag SNPs, which included three aforementioned polymorphisms. Finally, we performed a meta-analysis combining the result of our present study and the data sets of the two previously published studies.



Unrelated southern Han Chinese patients with GD were recruited from hospitals in Shanghai. The diagnosis of GD was based on documented clinical and biochemical evidence of hyperthyroidism, diffused goiter, and the presence of at least one of the following: positive TSH receptor antibody tests, diffusely increased 131I (iodine-131) uptake in the thyroid gland, or presence of exophthalmos. The severity of Graves' ophthalmopathy (GO) was assessed according to the NOSPECS classification[13]. Information was also collected from examination of medical records and all other available sources. All control subjects were given tests for thyroid function and autoantibody status, and anyone showing evidence of subclinical autoimmune thyroid disease (AITD) was removed from the study. Control subjects with known family history of autoimmune disease were also excluded.

In total, 428 GD patients and 690 unrelated control subjects were recruited for the case-control study. Of the GD patients, 132 showed ophthalmopathy defined as ATA class III or worse and were classified as GO (+). A sample set of 24 GD cases and 24 controls was used for resequencing. With informed consent, blood samples were collected from all the participants for DNA preparation as well as biochemical measurements. Genomic DNA was isolated from peripheral blood leukocytes with FlexiGene DNA Kit (Qiagenn Hilden, Germany). This project was approved by the Ethical Committee of the Chinese National Human Genome Center at Shanghai for the involvement of human subjects.

SNP identification by resequencing

We carried out SNP identification by resequencing a 1050 bp stretch located 5'-upstream of transcription start codon and the entire coding regions of ADRB2 in 24 patients and 24 control subjects. The sequencing reactions were performed using Applied Biosystems BigDye (version 3.1) chemistry (Applied Biosystem, Foster City, CA, USA), and the sequences were resolved using an ABI 3730 Genetic Analyzer (Applied Biosystem, Foster City, CA, USA). Analyses of the sequence traces were performed using the Staden package[14] and were double scored by a second operator.


Genotyping for the tag SNPs was performed using TaqMan chemistry (Applied Biosystem, Foster City, CA, USA). The TaqMans reaction was carried out following Applied Biosystems standard protocol. Fluorescent emissions were captured using the ABI Prisms 7900HT Sequence Detection System (Applied Biosystem, Foster City, CA, USA) and converted into genotypes using the ABI Prisms SDS 2.0 software package. The data completion rate was 98.5%.

Statistical analysis

We used χ2 analysis to evaluate the significance of differences in genotype and allele frequencies between cases and controls. Allele frequencies for cases and controls were used to calculate the odds ratio (OR) and the 95% confidence interval (95% CI). The software used for statistical calculations was the SPSS 15.0 (SPSS Inc., Chicago, IL, USA) unless specified. Genotype distributions in patients and controls were evaluated for departure from Hardy-Weinberg equilibrium by the Haploview v4.0 software[15]. From the SNPs detected by resequencing, we selected tag SNPs with a minimum threshold value of 0.8 for the r2 parameter. The tag SNPs were further genotyped in the collections of 428 GD patients and 690 unrelated controls. Blocks were determined using the default method of Gabriel et al[16] and haplotypes frequencies were estimated using an implementation of the expectation maximization algorithm by Haploview v4.0. We tested association with χ2 analysis for each possible combination of haplotypes derived from rs1042713 and rs1042714 using Haploview v4.0. P-values less than 0.05 were considered statistically significant.

Heterogeneity between studies was evaluated using the I2 statistic for inconsistency[17] and the χ2 distributed Cochran Q-statistic[18]. I2 is given by 100% × (Q - df)/Q, where df denotes degrees of freedom. It describes the proportion of variation that is unlikely due to chance and is considered large for values > 50%[19]. The Q-test is the most widely used test for heterogeneity, but is recognized to have poor power when there are few studies. Hence Q is considered statistically significant for P < 0.10[17, 18]. However I2 is unaffected by the number of studies and consequently is useful when comparing subgroups within the overall study. The I2 estimate and the Q-statistic were both computed for comparison. A fixed-effect model using Mantel-Haenszel method and a random-effects model using DerSimonian and Laird method were used to pool the results [20]. Random effects models are more conservative than fixed effects models and generate a wider confidence interval. Meta-analysis was performed using Review Manager software (version 4.2)


Fifteen polymorphisms were identified through the resequencing of ADRB2 for 48 individuals (Table 1), all of which were SNPs; one SNP (-332) located at position -332 upstream from the ADRB2 transcription initiation site was novel when compared with dbSNP Build 128. The nonsynonymous SNP rs1800888 at nt 491 (+491 C→T, 164 Thr→Ile) in the coding region was found not to be polymorphic in these 48 Chinese individuals. Among the fifteen SNPs, twelve SNPs had minor allele frequencies (MAF) great than 5% (Table 1), with frequencies ranging from 5.2% (rs1042711, rs1801704 and rs1042714) to 38.9% (rs1042719). The frequencies of the fifteen SNPs in 48 individuals were in Hardy-Weinberg equilibrium (P > 0.05, Table 1). The linkage disequilibrium (LD) structure among the SNPs was examined by program Haploview 4.0 (Figure 1). The LD between rs11959427 and rs1042714 (D' = 1, r2 = 0.82) is similar to a previous report by Japanese researchers [6], but much stronger than that in the Polish study [7].

Table 1 Fifteen SNPs in the ADRB2 gene identified from 24 GD patients and 24 control subjects
Figure 1

Linkage disequilibrium (LD) plots of the fifteen SNPs in ADRB2 identified from 48 subjects. We constructed the plots with the Haploview program [15], and r2 (×100) values were depicted in the diamonds. Blocks were determined using the default method of Gabriel et al [16].

All the fifteen SNPs were subsequently included in the tag SNP selection, and ten tag SNPs were selected and genotyped in the case-control collection. Allele and genotype frequencies for these tag SNPs are shown in Table 2. All tag SNPs genotypes in cases and controls conformed to Hardy-Weinberg equilibrium (P > 0.001). Genotype distributions and allele frequencies of the ten tag SNPs in ADRB2 were not significantly different between the GD group and the control group (Table 2). When under the assumption of dominant and recessive model, all the tag SNPs failed to show associations with GD either (data not shown here). Since the association of rs1042714 with GD was stronger in the patients without symptoms of GO in the Polish study[7], we also compared the allele and genotype frequencies of rs1042714 in GD subjects with or without GO to controls, and observed no difference existed (Table 3). No blocks were detected based on the ten tag SNPs. To explore the associations of haplotype rs11959427/rs1042713/rs1042714 (-367T/+46A/+79C) with GD as previously reported[7], we inferred the haplotypes based on rs1042713 and rs1042714, since rs11959427 captured by rs1042714 at r2 ≥ 0.8 were not directly genotyped in our case-control collection. Three common haplotypes were observed and none of them showed statistically significant association with GD (Table 4).

Table 2 Association analysis of the ten tag SNPs in GD patients and control subjects
Table 3 Association analysis of the rs1042714 polymorphism in control subjects and GD patients with or without GO.
Table 4 Estimated frequencies of the haplotypes inferred from rs1042713 and rs1042714 in GD patients and controls

A meta-analysis for rs1042714 was conducted combining the initial findings from the Polish populations, the data set from the study of WTCCC et al. and the result of present study, including a total of 1,728 GD patients and 2,491 control samples. The forest plots for the allele model are shown in Figure 2. The p-value for Q-test is just above the significance level (P = 0.12, Figure 2A), whereas the I2 statistic suggested there was significant heterogeneity among the studies (I2 = 53.1%, Figure 2A). Since only three studies were included in our meta-analysis, we chose the result of I2 statistic as measurement of heterogeneity. Then the pooled OR was calculated by random effects approaches using DerSimonian and Laird methods[20]. The pooled OR was 1.13 but the confidence interval spanned unity (95% CI, 0.95 to 1.36), which indicated the combined data did not show significant association for rs1042714 with GD. A sensitivity analysis showed that the present study in Chinese population was the main cause of the heterogeneity. When only the two European studies (including 1300 cases and 1801 controls) were combined, the heterogeneity was no longer significant (I2 = 5.9%, Figure 2B), and the result demonstrated significant association by fixed effect (OR = 1.18, 95% CI, 1.06 to 1.32; P = 0.002, Figure 2B).

Figure 2

Meta-analysis of the association between rs1042714 and GD measured by allele frequency data. OR (black squares) and 95% CI (bar) are shown for each study. The pooled ORs and their 95% CIs are represented by the shaded diamonds. Summary ORs are given for each ethnic group as well as all groups combined. The symbol n indicates the total number of C alleles, and N indicates the total number of C plus G alleles. A) DerSimonian-Laird OR meta-analysis (random effects) of the three studies; B) Mantel-Haenszel OR meta-analysis (fixed effects) of the two European studies.


We have performed a comprehensive study investigating association of GD susceptibility with the ADRB2 variants and failed to replicate association in Chinese population. However, meta-analysis revealed that rs1042714 was significantly associated with GD in Caucasian population. Genetic association studies of GD applied to other non-MHC candidate genes have also led to encouraging yet somewhat inconsistent findings in recent years. Several reasons may explain these inconsistent findings, including allelic heterogeneity, true variation in disease association between populations, modifying genetic and/or environmental factors, statistically underpowered small sample sizes, and a low prior probability of obtaining a true result[21]. The most parsimonious explanation for the lack of replication of initial ADRB2 association results in our study would be variation in disease association between populations and relatively small sample size.

It has long been noted that the MAF of rs1042714 locus varies significantly among the major populations around the world[5, 22]. The minor allele of rs1042714 loci is much more common in Caucasians than in Africans and Asians according to the previous reports and Hapmap data[5, 22]. Xie et al [22]and Chang et al [23] observed that the MAF of rs1042714 locus was 7.2% among 104 healthy Chinese and 7.7% among 260 nondiabetic Chinese respectively. The Hapmap data revealed the MAF of rs1042714 locus was 12.2% among 45 Chinese. In our study, the MAF of rs1042714 locus among the 690 healthy controls was 7.2%, consistent with the previous report. In the previous reports of Jazdzewski K et al in Polish population[7] and WTCCC et al in UK population[8], carriers of the minor allele G were present less frequently in the patient groups (38% and 43% respectively) than in the control groups (45% and 46% respectively), whereas in our data of Chinese population, the proportion of minor allele G carriers was larger in cases (8%) than in controls (7%). The trend of rs1042714 allele effect in cases and controls was substantially different in the two populations, which was also reflected by the OR value in the meta-analysis (less or greater than 1, Figure 2A). The allelic heterogeneity and opposite effect directions might indicate variations in disease associations between populations.

Genetic heterogeneity across populations had previously been reported in GD and other polygenic autoimmune diseases, such as in the case of SUMO4 gene association with type 1 diabetes, GD, or rheumatoid arthritis (RA) [2433]. SUMO4 gene was suggested as a general autoimmunity locus in the Asian population, associations with type 1 diabetes, GD and RA have been reported[27, 3133]. In contrast to the consistent association observed in the Asian populations, no evidence of SUMO4 association with susceptibility to type 1 diabetes, GD and RA was found from several studies in European descent[26, 2830]. One study did show SUMO4 was associated with T1D in association with high risk HLADR3 and DR4 genotypes in Swedish patients. But SUMO4 by itself was not associated with the disease[34]. Interestingly, a tendency for an opposite association was reported in type 1 diabetes of British subjects[25].

The major limitation of our case-control study was that the sample size was not very compelling. Sample size in thousands of subjects might be more sufficient to detect small genetic effects. However, positive association was reported from analysis of 300 GD patients and 301 control samples in Polish population, and positive signal was also found from the GWA screen of 1,000 GD patients and 1,500 control samples in UK population[7, 8]. This prompted us to investigate the potential association between ADRB2 and GD in fair number of well-defined cases and controls (428 GD cases and 690 controls) in Chinese population. Assuming a frequency of the risk allele of 0.075, our study in 1,118 subjects had only 70% power to detect the OR of 1.5 for rs1042714 with a significance level of 0.05. We sought to increase the effective sample size and statistical power by meta-analysis and to draw a more compelling result, but heterogeneity was found to exist between European studies and our study, and the data involved in ADRB2 association with GD in Chinese and Caucasian populations must be analyzed separately. The possibility that SNPs in ADRB2 may influence the disease risk modestly in Chinese population cannot be excluded completely. Further accumulation of samples would be needed for ascertainment. The overall sample size included in the meta-analysis was relative large (3101 Caucasian individuals) and this meta-analysis provided an up-to-date summary of the association analysis between ADRB2 and GD. However, the number of published studies available for meta-analysis was small, which made our analysis more prone to publication bias. Thus, it seems a little premature to draw conclusions that ADRB2 is associated with GD susceptibility in Caucasian populations.

A GD linkage locus had been mapped to chromosome 5q31-33 in Asian population in two GD genome screen studies [2, 4], but not detected in a Caucasian population[35]. The genome wide association study in Caucasians revealed some positive signal in this region[8], albeit not as strong as other loci. Our meta-analysis showed association of the ADRB2 variation with GD in Caucasians, which provided new evidence for the presence of GD susceptibility genes in 5q31-33. Much effort has been spent in exploring the associated genes and/or the causal allele(s) of Crohn disease[36, 37] and rheumatoid arthritis[38] located on the chromosome region of 5q31-33, and great progress has been made. However, no confirmed GD association genes have been identified in this region so far. Of note, different autoimmune diseases may share similar genetic etiologies and susceptibility alleles; therefore, the data sets from Crohn disease and rheumatoid arthritis studies will be helpful in further clarifying the variant accounted for susceptibility to GD in this region. A more comprehensive analysis of chromosome 5q31-33, including a wider range of SNPs in ADRB2 and neighboring genes in a large, well matched dataset, is required to determine whether a causal variant within any candidates in this linkage region confers susceptibility to GD.


We systematically studied the association of the ADRB2 gene by a case-control analysis using tag SNPs and a meta-analysis. The case-control study did not support ADRB2 as a GD susceptibility gene in Chinese population, but our meta-analysis provided strong evidence of a positive association of the ADRB2 SNP rs1042714 with GD in Caucasians. More studies with large sample size in Chinese and other populations are needed before solid conclusions about the association could be made.


  1. 1.

    Vanderpump PJ, Tunbridge WMG: The Epidemiology of Autoimmune Thyroid Disease. 1999, Humana Press: New Jerse

    Chapter  Google Scholar 

  2. 2.

    Jin Y, Teng W, Ben S, Xiong X, Zhang J, Xu S, Shugart YY, Jin L, Chen J, Huang W: Genome-wide scan of Graves' disease: evidence for linkage on chromosome 5q31 in Chinese Han pedigrees. J Clin Endocrinol Metab. 2003, 88 (4): 1798-1803. 10.1210/jc.2001-011980.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Akamizu T, Hiratani H, Ikegami S, Rich SS, Bowden DW: Association study of autoimmune thyroid disease at 5q23-q33 in Japanese patients. J Hum Genet. 2003, 48 (5): 236-242. 10.1007/s10038-003-0017-3.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Sakai K, Shirasawa S, Ishikawa N, Ito K, Tamai H, Kuma K, Akamizu T, Tanimura M, Furugaki K, Yamamoto K, et al: Identification of susceptibility loci for autoimmune thyroid disease to 5q31-q33 and Hashimoto's thyroiditis to 8q23-q24 by multipoint affected sib-pair linkage analysis in Japanese. Hum Mol Genet. 2001, 10 (13): 1379-1386. 10.1093/hmg/10.13.1379.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Drysdale CM, McGraw DW, Stack CB, Stephens JC, Judson RS, Nandabalan K, Arnold K, Ruano G, Liggett SB: Complex promoter and coding region beta 2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc Natl Acad Sci USA. 2000, 97 (19): 10483-10488. 10.1073/pnas.97.19.10483.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Munakata M, Harada Y, Ishida T, Saito J, Nagabukuro A, Matsushita H, Koga N, Ohsaki M, Imagawa K, Shiratsuchi T: Molecular-based haplotype analysis of the beta 2-adrenergic receptor gene (ADRB2) in Japanese asthmatic and non-asthmatic subjects. Allergol Int. 2006, 55 (2): 191-198. 10.2332/allergolint.55.191.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Jazdzewski K, Bednarczuk T, Stepnowska M, Liyanarachchi S, Suchecka-Rachon K, Limon J, Narkiewicz K: beta-2-adrenergic receptor gene polymorphism confers susceptibility to Graves disease. Int J Mol Med. 2007, 19 (1): 181-186.

    PubMed  PubMed Central  Google Scholar 

  8. 8.

    Burton PR, Clayton DG, Cardon LR, Craddock N, Deloukas P, Duncanson A, Kwiatkowski DP, McCarthy MI, Ouwehand WH, Samani NJ, et al: Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007, 39 (11): 1329-1337. 10.1038/ng.2007.17.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Elenkov IJ, Wilder RL, Chrousos GP, Vizi ES: The sympathetic nerve – an integrative interface between two supersystems: the brain and the immune system. Pharmacol Rev. 2000, 52 (4): 595-638.

    CAS  PubMed  Google Scholar 

  10. 10.

    McGraw DW, Forbes SL, Kramer LA, Liggett SB: Polymorphisms of the 5' leader cistron of the human beta2-adrenergic receptor regulate receptor expression. J Clin Invest. 1998, 102 (11): 1927-1932. 10.1172/JCI4862.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Chong LK, Chowdry J, Ghahramani P, Peachell PT: Influence of genetic polymorphisms in the beta2-adrenoceptor on desensitization in human lung mast cells. Pharmacogenetics. 2000, 10 (2): 153-162. 10.1097/00008571-200003000-00007.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Malysheva O, Pierer M, Wagner U, Wahle M, Wagner U, Baerwald CG: Association between beta2 adrenergic receptor polymorphisms and rheumatoid arthritis in conjunction with human leukocyte antigen (HLA)-DRB1 shared epitope. Ann Rheum Dis. 2008, 67 (12): 1759-1764. 10.1136/ard.2007.083782.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Werner SC: Modification of the classification of the eye changes of Graves' disease: recommendations of the Ad Hoc Committee of the American Thyroid Association. J Clin Endocrinol Metab. 1977, 44 (1): 203-204.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Bonfield JK, Rada C, Staden R: Automated detection of point mutations using fluorescent sequence trace subtraction. Nucleic Acids Res. 1998, 26 (14): 3404-3409. 10.1093/nar/26.14.3404.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Higgins JP, Thompson SG, Deeks JJ, Altman DG: Measuring inconsistency in meta-analyses. Bmj. 2003, 327 (7414): 557-560. 10.1136/bmj.327.7414.557.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Lau J, Ioannidis JP, Schmid CH: Quantitative synthesis in systematic reviews. Ann Intern Med. 1997, 127 (9): 820-826.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Higgins JP, Thompson SG: Quantifying heterogeneity in a meta-analysis. Stat Med. 2002, 21 (11): 1539-1558. 10.1002/sim.1186.

    Article  PubMed  Google Scholar 

  20. 20.

    DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7 (3): 177-188. 10.1016/0197-2456(86)90046-2.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Tello-Ruiz MK, Curley C, DelMonte T, Giallourakis C, Kirby A, Miller K, Wild G, Cohen A, Langelier D, Latiano A, et al: Haplotype-based association analysis of 56 functional candidate genes in the IBD6 locus on chromosome 19. Eur J Hum Genet. 2006, 14 (6): 780-790. 10.1038/sj.ejhg.5201612.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Xie HG, Stein CM, Kim RB, Xiao ZS, He N, Zhou HH, Gainer JV, Brown NJ, Haines JL, Wood AJ: Frequency of functionally important beta-2 adrenoceptor polymorphisms varies markedly among African-American, Caucasian and Chinese individuals. Pharmacogenetics. 1999, 9 (4): 511-516.

    CAS  PubMed  Google Scholar 

  23. 23.

    Chang TJ, Tsai MH, Jiang YD, Lee B, Lee KC, Lin JY, Chiu KC, Tai TY, Chuang LM: The Arg16Gly polymorphism of human beta2-adrenoreceptor is associated with type 2 diabetes in Taiwanese people. Clin Endocrinol (Oxf). 2002, 57 (5): 685-690. 10.1046/j.1365-2265.2002.01661.x.

    CAS  Article  Google Scholar 

  24. 24.

    Guo D, Li M, Zhang Y, Yang P, Eckenrode S, Hopkins D, Zheng W, Purohit S, Podolsky RH, Muir A, et al: A functional variant of SUMO4, a new I kappa B alpha modifier, is associated with type 1 diabetes. Nat Genet. 2004, 36 (8): 837-841. 10.1038/ng1391.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Bohren KM, Nadkarni V, Song JH, Gabbay KH, Owerbach D: A M55V polymorphism in a novel SUMO gene (SUMO-4) differentially activates heat shock transcription factors and is associated with susceptibility to type I diabetes mellitus. J Biol Chem. 2004, 279 (26): 27233-27238. 10.1074/jbc.M402273200.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Park Y, Park S, Kang J, Yang S, Kim D: Assessing the validity of the association between the SUMO4 M55V variant and risk of type 1 diabetes. Nat Genet. 2005, 37 (2): 112-10.1038/ng0205-112a. author reply 112–113

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Qu H, Bharaj B, Liu XQ, Curtis JA, Newhook LA, Paterson AD, Hudson TJ, Polychronakos C: Assessing the validity of the association between the SUMO4 M55V variant and risk of type 1 diabetes. Nat Genet. 2005, 37 (2): 111-112. 10.1038/ng0205-111. author reply 112–113

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Smyth DJ, Howson JM, Lowe CE, Walker NM, Lam AC, Nutland S, Hutchings J, Tuomilehto-Wolf E, Tuomilehto J, Guja C, et al: Assessing the validity of the association between the SUMO4 M55V variant and risk of type 1 diabetes. Nat Genet. 2005, 37 (2): 110-111. 10.1038/ng0205-110. author reply 112–113

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Gibbons LJ, Thomson W, Zeggini E, Worthington J, Barton A, Eyre S, Donn R, Hinks A: The type 1 diabetes susceptibility gene SUMO4 at IDDM5 is not associated with susceptibility to rheumatoid arthritis or juvenile idiopathic arthritis. Rheumatology (Oxford). 2005, 44 (11): 1390-1393. 10.1093/rheumatology/kei041.

    CAS  Article  Google Scholar 

  30. 30.

    Jennings CE, Owen CJ, Wilson V, Pearce SH: No association of the codon 55 methionine to valine polymorphism in the SUMO4 gene with Graves' disease. Clin Endocrinol (Oxf). 2005, 62 (3): 362-365. 10.1111/j.1365-2265.2005.02224.x.

    CAS  Article  Google Scholar 

  31. 31.

    Noso S, Ikegami H, Fujisawa T, Kawabata Y, Asano K, Hiromine Y, Tsurumaru M, Sugihara S, Lee I, Kawasaki E, et al: Genetic Heterogeneity in Association of the SUMO4 M55V Variant With Susceptibility to Type 1 Diabetes. Diabetes. 2005, 54 (12): 3582-3586. 10.2337/diabetes.54.12.3582.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Noso S, Ikegami H, Fujisawa T, Kawabata Y, Asano K, Hiromine Y, Sugihara S, Lee I, Kawasaki E, Awata T, et al: Association of SUMO4, as a candidate gene for IDDM5, with susceptibility to type 1 diabetes in Asian populations. Ann N Y Acad Sci. 2006, 1079: 41-46. 10.1196/annals.1375.006.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Tsurumaru M, Kawasaki E, Ida H, Migita K, Moriuchi A, Fukushima K, Fukushima T, Abiru N, Yamasaki H, Noso S, et al: Evidence for the Role of Small Ubiquitin-Like Modifier 4 as a General Autoimmunity Locus in the Japanese Population. J Clin Endocrinol Metab. 2006, 91 (8): 3138-3143. 10.1210/jc.2006-0206.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Sedimbi SK, Luo XR, Sanjeevi CB: SUMO4 M55V polymorphism affects susceptibility to type I diabetes in HLA DR3- and DR4-positive Swedish patients. Genes Immun. 2007, 8 (6): 518-521. 10.1038/sj.gene.6364406.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Tomer Y, Ban Y, Concepcion E, Barbesino G, Villanueva R, Greenberg DA, Davies TF: Common and unique susceptibility loci in Graves and Hashimoto diseases: results of whole-genome screening in a data set of 102 multiplex families. Am J Hum Genet. 2003, 73 (4): 736-747. 10.1086/378588.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Rioux JD, Daly MJ, Silverberg MS, Lindblad K, Steinhart H, Cohen Z, Delmonte T, Kocher K, Miller K, Guschwan S, et al: Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease. Nat Genet. 2001, 29 (2): 223-228. 10.1038/ng1001-223.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Silverberg MS, Duerr RH, Brant SR, Bromfield G, Datta LW, Jani N, Kane SV, Rotter JI, Philip Schumm L, Hillary Steinhart A, et al: Refined genomic localization and ethnic differences observed for the IBD5 association with Crohn's disease. Eur J Hum Genet. 2007, 15 (3): 328-335. 10.1038/sj.ejhg.5201756.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Tokuhiro S, Yamada R, Chang X, Suzuki A, Kochi Y, Sawada T, Suzuki M, Nagasaki M, Ohtsuki M, Ono M, et al: An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat Genet. 2003, 35 (4): 341-348. 10.1038/ng1267.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


This work was supported by the grants from Chinese High-Tech Program (2006AA020706, 2006AA02A406), National Key Project for Basic Research (2004CB518605), Chinese National Natural Science Fund for Distinguished Young Scholars (30625019), Shanghai Science and Technology Committee (06XD14015) and Chinese National Natural Science Fund (30771029).

Author information



Corresponding author

Correspondence to Wei Huang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

XC designed the study, carried out the statistical analysis, and wrote the manuscript. MS, BW, KZ and QH conducted the molecular genotyping of the variants. LS and YD were involved in the recruitment of the patients and controls. CD, YW and SX assisted in manuscript revising. HW conceptualized the study, supervised the study and obtained the funding. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Chu, X., Dong, Y., Shen, M. et al. Polymorphisms in the ADRB2 gene and Graves disease: a case-control study and a meta-analysis of available evidence. BMC Med Genet 10, 26 (2009).

Download citation


  • Grave Disease
  • Grave Disease Patient
  • Wellcome Trust Case Control Consortium
  • Rs1042714 Locus
  • ADRB2 Gene