Human Molecular Genetics Advance Access originally published online on June 27, 2007
Human Molecular Genetics 2007 16(18):2149-2153; doi:10.1093/hmg/ddm165
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
A novel and major association of HLA-C in Graves' disease that eclipses the classical HLA-DRB1 effect
1 Division of Medical Sciences, Institute of Biomedical Research, University of Birmingham, Birmingham B15 2TT, UK and 2 Juvenile Diabetes Research Foundation/Wellcome Trust Diabetes and Inflammation Laboratory, Cambridge Institute for Medical Research, University of Cambridge, Cambridge CB2 0XY, UK
* To whom correspondence should be addressed. Tel: +44 1214158819; Fax: +44 1214158712; Email: s.c.gough{at}bham.ac.uk
Received April 27, 2007; Accepted June 22, 2007
| ABSTRACT |
|---|
|
|
|---|
Association of the major histocompatibility complex (MHC) class II-encoded HLA-DRB1-DQA1-DQB1 haplotype with Graves' disease (GD) has been known for several years. Recent evidence from other autoimmune diseases has suggested that the HLA class I encoded HLA-B/-C molecules could be conferring HLA-DRB1-DQA1-DQB1 independent effects on disease. The aim of this study was to determine the effect of HLA-B and HLA-C in GD in a white ethnic group of 806 patients with GD and 487 control subjects from the UK. Of the five loci (HLA-B, -C, -DRB1, -DQA1, -DQB1), HLA-C demonstrated the strongest association (P = 1.20 x 10–20) with HLA-C*07 predisposing [OR = 1.63, 95% CI (1.23–2.17)] and both HLA-C*03 [OR = 0.54, 95% CI (0.38–0.77)], HLA-C*16 [OR = 0.36, 95% CI (0.21–0.61)] protective. The other loci were then tested for HLA-C-independent associations. HLA-B was found to be associated independently of HLA-C (P = 1.54 x 10–6) with the other three loci, HLA-DRB1, HLA-DQB1 and HLA-DQA1, also improving the model but with less confidence (P > 10–5). This study has for the first time provided evidence of a primary association of HLA-C, and to a lesser extent HLA-B, with GD. Class II loci could still have effects on GD, but they appear smaller than the HLA-C association. A full investigation of the MHC region, including all class I and II loci is now required. Our results point to a primary role for class I-mediated responses in GD, a condition classically assumed to be a straightforward HLA-class II-restricted autoantibody response to the thyroid stimulating hormone receptor.
| INTRODUCTION |
|---|
|
|
|---|
Components of the major histocompatibility complex (MHC) region are attractive candidate loci for autoimmune disease (AID), with the DRB1–DQA1–DQB1 molecules extensively examined for association with a series of AIDs including Graves' disease (GD). The HLA class II region, including the DRB1*03 and DQA1*0501 alleles and the DRB1*03–DQA1*0501–DQB1*02 (DR3) haplotype, have been consistently associated with GD (1). We have previously reported that the individual contributions of HLA-DQA1 and HLA-DRB1, to the association of the HLA class II-encoded HLA-DRB1_-DQA1_-DQB1 susceptibility haplotypes with GD were indistinguishable (2). We, and others (2,3), have proposed that this association maybe attributable, in part, to HLA-DRB1 exon 2-encoded position ß74.
Recently, in the AIDs, type 1 diabetes and multiple sclerosis, evidence has arisen suggesting the existence of susceptibility loci within the HLA class I region, independent of known class II effects (4–7) (J.M.M. Howson and J.A. Todd, manuscript in preparation). Interestingly, upregulation of HLA class I molecules on thyroid cells, including HLA-B and HLA-C molecules, has been demonstrated in response to immune cell infiltration of the thyroid gland, which may be one of the earliest features of autoimmune attack in GD (8). Methimazole used in the treatment of GD has been shown to function, in part, by reducing HLA class I and class II expression on thyrocytes (8,9).
Owing to the strong evidence for a role of HLA class I-encoded molecules in GD susceptibility, the aim of this study was to extend the typing from our original MHC class II study (2) to investigate the contribution of the MHC class I encoded HLA-B and HLA-C genes to GD.
| RESULTS |
|---|
|
|
|---|
The newly typed HLA-B and HLA-C loci were combined with the previously typed HLA-DRB1, HLA-DQB1 and HLA-DQA1 data. The strongest association signal came from HLA-C (P = 1.20 x 10–20; Table 1) with the most common allele, HLA-C*07 being the most predisposing (OR = 1.63, 95% CI = 1.23–2.17; Table 2). Protective effects were observed for HLA-C*03 and HLA-C*16: OR = 0.54, 95% CI = 0.38–0.77 and OR = 0.36, 95% CI = 0.21–0.61, respectively, using the approximately neutral HLA-C*05 allele as a reference (Table 2). HLA-B was also strongly associated (P = 2.31 x 10–7; Table 1). The allele, HLA-B*08 (OR = 1.62, 95% CI = 1.19–2.19; Table 3) was the most predisposing (using the neutral B*07 allele as reference), with a protective association for B*44 (OR = 0.64, 95% CI = 0.47–0.88; Table 3).
|
|
|
Next, we tested if HLA-B, HLA-DRB1, HLA-DQA1 or HLA-DQB1, had an association with GD that was independent of HLA-C. Each locus in turn was added to the alleles or genotypes of HLA-C in a logistic model and a likelihood ratio test used to evaluate whether a given locus improved the model (10). HLA-B improved a model containing HLA-C whether the alleles or the genotypes of HLA-C were used as the null model (P = 1.54 x 10–6 and P = 1.51 x 10–5, respectively; Table 4). All three MHC class II loci, HLA-DRB1, HLA-DQB1 and HLA-DQA1, also improved HLA-C. Conversely, HLA-C significantly improved all models when added to HLA-B, HLA-DRB1, HLA-DQA1 or HLA-DQB1 (P = 8.68 x 10–13–3.26 x 10–19).
|
As HLA-B may be contributing to the association of the MHC with GD, we examined the allelic ORs of both HLA-C and HLA-B having conditioned on the effects of the other MHC class I locus. Conditioning on HLA-B had little effect on the HLA-C allelic ORs (except the C*07 allele whose 95% CI crossed 1; Table 2). However, conditioning on HLA-C had a greater effect on the HLA-B ORs, with B*44 no longer showing a protective association (Table 3).
We found there was linkage disequilibrium (LD) between the previously associated HLA-DRB1*03, *07 and the HLA-C*07, *03 alleles. HLA-C*03 has D' = 0.88 with HLA-DRB1*03 and D' = 0.75 with HLA-DRB1*07, and HLA-C*07 has D' = 0.66 and 0.65 with HLA-DRB1*03 and HLA-DRB1*07, respectively. The most common HLA-C_-B_-DRB1 haplotype was HLA-C*07_HLA-B*08_HLA-DRB1*03 in both cases and controls. When LD was investigated between HLA-C*07 and the HLA-B*08-HLA-DRB1*03 haplotype, it was shown that the HLA-B*08-HLA-DRB1*03 haplotype did not add to HLA-C*07 (P = 0.167), suggesting that the results seen at HLA-C*07 were not attributed to LD with the HLA-B*08-HLA-DRB1*03 haplotype.
| DISCUSSION |
|---|
|
|
|---|
Strong LD confounds determination of the exact location(s) of aetiological variant(s) in the MHC region. This study has for the first time provided evidence in support of a primary association of HLA-C with GD, which appears to eclipse the classical association of HLA-DRB1 in GD. Furthermore, a second class I locus, HLA-B, may also be associated with GD. This work therefore calls into doubt a primary contribution of position ß74 to GD. However, we cannot rule out independent effects for components of the class II-encoded HLA-DRB1/HLA-DQA1/HLA-DQB1 haplotypes in disease susceptibility in the British population but they do appear, unexpectedly, to be smaller than the MHC class I gene, HLA-C.
Early studies investigating the HLA region and GD over 30 years ago reported association of the class I-encoded HLA-B*08 (11–15) and HLA-A*08 (16) alleles. However, owing to strong LD and an inability of statistical methodology at the time to split these effects, association between class I loci and disease was presumed to be secondary to the class II region which appeared to be exerting a larger effect (17). Consequently, most subsequent studies focused on the class II loci. Interestingly, a more recent study investigated the HLA class I region, reported an association of the class I encoded HLA-A*02 in a small Taiwanese population, although conditional analyses to exclude a primary class II effect was not performed (18).
HLA class I molecules could be playing a greater role than class II molecules in GD for a number of reasons. HLA class I molecules bind and present internally derived peptides including viral or bacterial antigens which may be linked to disease initiation through molecular mimicry (19). The most postulated bacterial trigger for GD is Yersinia enterocolitica. Antibodies raised against Y. enterocolitica have been shown to displace the thyroid stimulating hormone (TSH) from binding to its receptor (TSHR) (20,21), suggesting immune cross-reactivity between Y. enterocolitica and the TSHR. Several other viral triggers including human intracisternal type A retroviral particle (HIAP) (22) and gag protein from the human foamy virus (HPV) (23,24) have also been proposed, although further work is needed to determine a pathogenic role for these in GD. It has been suggested that viral or bacterial peptides could act as superantigens, causing a series of non-specific T cells to bind to the bacterial peptide/class I complex and become activated, some of which could cross-react, causing autoimmune attack of the thyroid gland (25).
Cytotoxic T cells, including natural killer (NK) cells, specific for TSHR peptides could have an early regulatory or pathogenic role in GD, as has been proposed in type 1 diabetes. Killer immunoglobulin-like receptors (KIR) on the surface of NK cells interact with specific HLA molecules and play a key role in inhibiting NK cell-mediated destruction. For HLA-C there are two sets of KIRs, 2DL1 and 2DL2/2DL3 which recognize specific sets of HLA-C molecules, group C1 molecules (containing C*02, *04, *05 and *06) and group C2 molecules (containing C*01, *03, *07 and *08), respectively (26). Groups C1 and C2 are distinguished by the presence of serine at position 77 and asparagine at position 80 of the
1 helix or asparagine at position 77 and a lysine at position 80, respectively. KIR/HLA-C interactions can be altered by peptide loading and presentation by HLA-C (27,28), which could suggest that interaction of the associated HLA-C molecules with a given autoantigen/s such as TSHR could be affecting KIR binding and that this could play a role in the onset of autoimmunity. Finally, specific class I molecules may be prone to misfolding or when they are themselves presented to the immune system by HLA class II molecules an autoimmune response could be triggered (29). These mechanisms are neither mutually exclusive nor exhaustive but represent hypotheses warranting further study.
Given these findings, a full investigation of the MHC region is now justified and required in larger, more statistically powerful and independent datasets, including those of differing ethnic origin, using the dense SNP maps now available. This should include, for example, the third class I locus, HLA-A, as well as MICA and MICB and the HLA-DPB1 locus, in order to localize all the primary GD variants and establish whether HLA-C is causative or simply in LD with the aetiological variant(s).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Subjects
Unrelated white patients of UK origin with GD were recruited from thyroid clinics in Birmingham, Walsall, Bournemouth and Exeter, UK, as previously described (2,30). Patients were defined as having GD by the presence of biochemical hyperthyroidism together with either the presence of dysthyroid eye disease (NOPPSECS score > 2) or two of the following criteria: diffuse goiter, a significant titer of microsomal, thyroglobulin or thyrotropin (TSH) receptor autoantibodies as described in detail previously (2,30). Control subjects were recruited from Birmingham, UK. All patients and subjects gave informed written consent and the project approved by the Local Research Ethics Committee.
In total, DNA was obtained from 950 unrelated white British GD patients and 621 British, age- and gender-matched controls. The cases were collected from three locations in Great Britain: 773 cases were from Birmingham in the Midlands; 100 cases were from Exeter in the Southwest; and 77 cases were from Bournemouth in the South. All 621 controls were collected from the Birmingham area. We found that HLA-C allele frequencies were geographically variable (P = 0.0008) and therefore, we only report analyses for the 1394 subjects from the Midlands. Geographical information and analyses for the full collection are reported in the Supplementary information (Supplementary Material, Tables S1–S8).
Materials and methods
DNA was prepared from whole blood using the Nucleon Bacc II kit (Tepnel Life Sciences PLC, UK). Genotyping was performed at the University of Birmingham, UK using the primer sequences and methods published by Bunce et al. (31) with primers obtained from Sigma Genosys (Poole, UK).
The HLA-DRB1, -DQB1 and -DQA1 regions were typed previously (2), however, 79 additional cases were typed for this study using the same method.
Statistical analysis
We were well powered to find the primary effects of the MHC loci in GD. With 800 cases and 500 controls assuming a multiplicative model and a minor allele frequency of 0.10, we had 98% power to find an effect size of 2.0 at a type 1 error rate,
= 0.0001. Logistic regression in STATA (http://www.stata.com) with software written by David Clayton for use within that package (http://www.gene.cimr.cam.ac.uk/clayton/software/stata) was used to test for association and calculate odds ratios (ORs) at all loci (10). No evidence was obtained for non-multiplicative inheritance at any locus studied (P > 0.3); hence, only allelic models, which assume a multiplicative model are reported. Analyses stratified by sex were performed; however, they were not qualitatively different to the unstratified analyses; hence, unstratified results are reported. Haplotypes were reconstructed in PHASE version 2.1.1 (32,33). All loci were in Hardy–Weinberg equilibrium in the controls (P > 0.06).
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at HMG Online.
| ACKNOWLEDGEMENTS |
|---|
We acknowledge the Wellcome Trust and the Juvenile Diabetes Research Foundation for funding.
Conflict of Interest statement. None of the authors have any conflicts of interest.
| REFERENCES |
|---|
|
|
|---|
- Simmonds M.J., Gough S.C. Genetic insights into disease mechanisms of autoimmunity. Br. Med. Bull. (2004) 71:93–113.[CrossRef][Web of Science][Medline]
- Simmonds M.J., Howson J.M., Heward J.M., Cordell H.J., Foxall H., Carr-Smith J., Gibson S.M., Walker N., Tomer Y., Franklyn J.A., et al. Regression mapping of association between the human leukocyte antigen region and Graves disease. Am. J. Hum. Genet. (2005) 76:157–163.[CrossRef][Web of Science][Medline]
- Ban Y., Davies T.F., Greenberg D.A., Concepcion E.S., Osman R., Oashi T., Tomer Y. Arginine at position 74 of the HLA-DR beta1 chain is associated with Graves' disease. Genes Immunol. (2004) 5:203–208.[CrossRef]
- Yeo T.W., De Jager P.L., Gregory S.G., Barcellos L.F., Walton A., Goris A., Fenoglio C., Ban M., Taylor C.J., Goodman R.S., et al. A second major histocompatibility complex susceptibility locus for multiple sclerosis. Ann. Neurol. (2007) 61:228–236.[CrossRef][Web of Science][Medline]
- Fogdell-Hahn A., Ligers A., Gronning M., Hillert J., Olerup O. Multiple sclerosis: a modifying influence of HLA class I genes in an HLA class II associated autoimmune disease. Tissue Antigens (2000) 55:140–148.[CrossRef][Web of Science][Medline]
- Lie B.A., Thorsby E. Several genes in the extended human MHC contribute to predisposition to autoimmune diseases. Curr. Opin. Immunol. (2005) 17:526–531.[Web of Science][Medline]
- Johansson S., Lie B.A., Todd J.A., Pociot F., Nerup J., Cambon-Thomsen A., Kockum I., Akselsen H.E., Thorsby E., Undlien D.E. Evidence of at least two type 1 diabetes susceptibility genes in the HLA complex distinct from HLA-DQB1, -DQA1 and -DRB1. Genes Immunol. (2003) 4:46–53.[CrossRef]
- Singer D.S., Mozes E., Kirshner S., Kohn L.D. Role of MHC class I molecules in autoimmune disease. Crit. Rev. Immunol. (1997) 17:463–468.[Web of Science][Medline]
- Saji M., Moriarty J., Ban T., Singer D.S., Kohn L.D. Major histocompatibility complex class I gene expression in rat thyroid cells is regulated by hormones, methimazole, and iodide as well as interferon. J. Clin. Endocrinol. Metab. (1992) 75:871–878.[Abstract]
- Cordell H.J., Clayton D.G. A unified stepwise regression procedure for evaluating the relative effects of polymorphisms within a gene using case/control or family data: application to HLA in type 1 diabetes. Am. J. Hum. Genet. (2002) 70:124–141.[CrossRef][Web of Science][Medline]
- Farid N.R., Barnard J.M., Marshall W.H. The association of HLA with autoimmune thyroid disease in Newfoundland. The influence of HLA homozygosity in Graves' disease. Tissue Antigens (1976) 8:181–189.[Web of Science][Medline]
- Mather B.A., Roberts D.F., Scanlon M.F., Mukhtar E.D., Davies T.F., Smith B.R., Hall R. HLA antigens and thyroid autoantibodies in patients with Graves' disease and their first degree relatives. Clin. Endocrinol. (Oxf.) (1980) 12:155–163.[Medline]
- Stenszky V., Kozma L., Balazs C., Bear J.C., Farid N.R. The role of HLA antigens in the manifestation and course of Graves' disease. Mol. Biol. Med. (1986) 3:53–62.[Web of Science][Medline]
- Farid N.R., Stone E., Johnson G. Graves' disease and HLA: clinical and epidemiologic associations. Clin. Endocrinol. (Oxf.) (1980) 13:535–544.[Medline]
-
Dahlberg P.A., Holmlund G., Karlsson F.A., Safwenberg J. HLA-A, -B, -C and -DR antigens in patients with Graves' disease and their correlation with signs and clinical course. Acta. Endocrinol. (Copenh.) (1981) 97:42–47.
[Abstract/Free Full Text] -
Grumet F.C., Payne R.O., Konishi J., Kriss J.P. HL-A antigens as markers for disease susceptibility and autoimmunity in Graves' disease. J. Clin. Endocrinol. Metab. (1974) 39:1115–1119.
[Abstract/Free Full Text] -
Bech K., Lumholtz B., Nerup J., Thomsen M., Platz P., Ryder L.P., Svejgaard A., Siersbaek-Nielsen K., Hansen J.M., Larsen J.H. HLA antigens in Graves' disease. Acta. Endocrinol. (Copenh.) (1977) 86:510–516.
[Abstract/Free Full Text] - Huang S.M., Wu T.J., Lee T.D., Yang E.K., Shaw C.K., Yeh C.C. The association of HLA -A, -B, and -DRB1 genotypes with Graves' disease in Taiwanese people. Tissue Antigens (2003) 61:154–158.[Web of Science][Medline]
-
Albert L.J., Inman R.D. Molecular mimicry and autoimmunity. N. Engl. J. Med. (1999) 341:2068–2074.
[Free Full Text] - Wenzel B.E., Peters A., Zubaschev I. Bacterial virulence antigens and the pathogenesis of autoimmune thyroid diseases (AITD). Exp. Clin. Endocrinol. Diabetes (1996) 104(Suppl. 4):75–78.[Web of Science][Medline]
- Wenzel B.E., Franke T.F., Heufelder A.E., Heesemann J. Autoimmune thyroid diseases and enteropathogenic Yersinia enterocolitica. Autoimmunity (1990) 7:295–303.[Web of Science][Medline]
- Jaspan J.B., Sullivan K., Garry R.F., Lopez M., Wolfe M., Clejan S., Yan C., Tenenbaum S., Sander D.M., Ahmed B., et al. The interaction of a type A retroviral particle and class II human leukocyte antigen susceptibility genes in the pathogenesis of Graves' disease. J. Clin. Endocrinol. Metab. (1996) 81:2271–2279.[Abstract]
- Wick G., Trieb K., Aguzzi A., Recheis H., Anderl H., Grubeck-Loebenstein B. Possible role of human foamy virus in Graves' disease. Intervirology (1993) 35:101–107.[Web of Science][Medline]
- Lee H., Kim S., Kang M., Kim W., Cho B. Prevalence of human foamy virus-related sequences in the Korean population. J. Biomed. Sci. (1998) 5:267–273.[CrossRef][Web of Science][Medline]
- Moller E. Mechanisms for induction of autoimmunity in humans. Acta. Paediatr. Suppl. (1998) 424:16–20.[Medline]
- Bashirova A.A., Martin M.P., McVicar D.W., Carrington M. The killer immunoglobulin-like receptor gene cluster: tuning the genome for defense. Annu. Rev. Genom. Hum. Genet. (2006) 7:277–300.[CrossRef][Web of Science][Medline]
-
Malnati M.S., Peruzzi M., Parker K.C., Biddison W.E., Ciccone E., Moretta A., Long E.O. Peptide specificity in the recognition of MHC class I by natural killer cell clones. Science (1995) 267:1016–1018.
[Abstract/Free Full Text] -
Mandelboim O., Wilson S.B., Vales-Gomez M., Reyburn H.T., Strominger J.L. Self and viral peptides can initiate lysis by autologous natural killer cells. Proc. Natl Acad. Sci. USA (1997) 94:4604–4609.
[Abstract/Free Full Text] - Marker-Hermann E., Meyer zum Buschenfelde K.H., Wildner G. HLA-B27-derived peptides as autoantigens for T lymphocytes in ankylosing spondylitis. Arthritis Rheum. (1997) 40:2047–2054.[Web of Science][Medline]
-
Simmonds M.J., Heward J.M., Carr-Smith J., Foxall H., Franklyn J.A., Gough S.C. Contribution of single nucleotide polymorphisms within FCRL3 and MAP3K7IP2 to the pathogenesis of Graves' disease. J. Clin. Endocrinol. Metab. (2006) 91:1056–1061.
[Abstract/Free Full Text] - Bunce M., O'Neill C.M., Barnardo M.C., Krausa P., Browning M.J., Morris P.J., Welsh K.I. Phototyping: comprehensive DNA typing for HLA-A, B, C, DRB1, DRB3, DRB4, DRB5 and DQB1 by PCR with 144 primer mixes utilizing sequence-specific primers (PCR-SSP). Tissue Antigens (1995) 46:355–367.[Web of Science][Medline]
- Stephens M., Smith N.J., Donnelly P. A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. (2001) 68:978–989.[CrossRef][Web of Science][Medline]
-
Stephens M., Scheet P. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am. J. Hum. Genet. (2005) 76:449–462.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
O. J. Brand, J. C. Barrett, M. J. Simmonds, P. R. Newby, C. J. McCabe, C. K. Bruce, B. Kysela, J. D. Carr-Smith, T. Brix, P. J. Hunt, et al. Association of the thyroid stimulating hormone receptor gene (TSHR) with Graves' disease Hum. Mol. Genet., May 1, 2009; 18(9): 1704 - 1713. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Anderson Update in Endocrine Autoimmunity J. Clin. Endocrinol. Metab., October 1, 2008; 93(10): 3663 - 3670. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Huber, E. M. Jacobson, K. Jazdzewski, E. S. Concepcion, and Y. Tomer Interleukin (IL)-23 Receptor Is a Major Susceptibility Gene for Graves' Ophthalmopathy: The IL-23/T-helper 17 Axis Extends to Thyroid Autoimmunity J. Clin. Endocrinol. Metab., March 1, 2008; 93(3): 1077 - 1081. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

