Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Interaction Analysis between HLA-DRB1 Shared Epitope Alleles and MHC Class II Transactivator CIITA Gene with Regard to Risk of Rheumatoid Arthritis

  • Marcus Ronninger ,

    marcus.ronninger@ki.se

    Affiliation Rheumatology Unit, Department of Medicine, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden

  • Maria Seddighzadeh,

    Affiliation Rheumatology Unit, Department of Medicine, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden

  • Morten Christoph Eike,

    Affiliation Department of Medical Genetics, University of Oslo and Oslo University Hospital, Oslo, Norway

  • Darren Plant,

    Affiliation 1arc-Epidemiology Unit, University of Manchester, Manchester, United Kingdom

  • Nina A. Daha,

    Affiliation Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands

  • Beate Skinningsrud,

    Affiliation Department of Medical Genetics, University of Oslo and Oslo University Hospital, Oslo, Norway

  • Jane Worthington,

    Affiliation 1arc-Epidemiology Unit, University of Manchester, Manchester, United Kingdom

  • Tore K. Kvien,

    Affiliation Department of Rheumatology, Diakonhjemmet Hospital, Oslo, Norway

  • Rene E. M. Toes,

    Affiliation Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands

  • Benedicte A. Lie,

    Affiliation Department of Medical Genetics, University of Oslo and Oslo University Hospital, Oslo, Norway

  • Lars Alfredsson,

    Affiliation Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden

  • Leonid Padyukov

    Affiliation Rheumatology Unit, Department of Medicine, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden

Abstract

HLA-DRB1 shared epitope (SE) alleles are the strongest genetic determinants for autoantibody positive rheumatoid arthritis (RA). One of the key regulators in expression of HLA class II receptors is MHC class II transactivator (CIITA). A variant of the CIITA gene has been found to associate with inflammatory diseases.

We wanted to explore whether the risk variant rs3087456 in the CIITA gene interacts with the HLA-DRB1 SE alleles regarding the risk of developing RA. We tested this hypothesis in a case-control study with 11767 individuals from four European Caucasian populations (6649 RA cases and 5118 controls).

We found no significant additive interaction for risk alleles among Swedish Caucasians with RA (n = 3869, attributable proportion due to interaction (AP) = 0.2, 95%CI: −0.2–0.5) or when stratifying for anti-citrullinated protein antibodies (ACPA) presence (ACPA positive disease: n = 2945, AP = 0.3, 95%CI: −0.05–0.6, ACPA negative: n = 2268, AP = −0.2, 95%CI: −1.0–0.6). We further found no significant interaction between the main subgroups of SE alleles (DRB1*01, DRB1*04 or DRB1*10) and CIITA. Similar analysis of three independent RA cohorts from British, Dutch and Norwegian populations also indicated an absence of significant interaction between genetic variants in CIITA and SE alleles with regard to RA risk.

Our data suggest that risk from the CIITA locus is independent of the major risk for RA from HLA-DRB1 SE alleles, given that no significant interaction between rs3087456 and SE alleles was observed. Since a biological link between products of these genes is evident, the genetic contribution from CIITA and class II antigens in the autoimmune process may involve additional unidentified factors.

Introduction

Rheumatoid arthritis (RA) is a relatively common disease of poorly understood aetiology that affects approximately 1% of the world's population. Even though the pathophysiology of the disease is well studied only a limited number of risk factors with low to moderate effect has been described [1]. Even the strongest genetic risk factor for RA, the variants in the HLA-DRB1 gene suggested by the shared epitope (SE) hypothesis [2], confer only a moderate risk increase for RA, with an odds ratio (OR) of 4–6 in European Caucasians with regard to anti-citrullinated protein antibodies (ACPA) positive RA [3], [4]. Also, these variations are quite common in the normal population and their predictive value is very low. Therefore, the major fraction of RA risk remains unexplained by existing information and interaction between known risk factors may account for putative “missing” risk factors. Indeed, some risk factors for RA has been shown to moderate the risk for disease in the context of the SE [3], [5], [6] suggesting that interactions with SE may play an important role in the development of RA [7].

We have earlier reported on a variant of the MHC class II transactivator (CIITA) gene, the −168A/G promoter SNP (rs3087456), which associates with inflammatory diseases and also the expression of CIITA and downstream HLA expression [8]. This association was not consistently replicated in different populations [9], [10]. However, other variants in the CIITA locus than rs3087456 have been reported in association with autoimmune disease [11], [12], [13] warranting further exploration of this locus in the context of RA.

Also, due to involvement of HLA class II in RA, the study of CIITA may reveal more detailed mechanisms of disease development, since the protein is known to be a key regulator of MHC class II expression and therefore may be involved in development of RA in combination with SE alleles [14], [15]. A complete lack of expression of CIITA leads to the bare lymphocyte syndrome with a complete abolishment of classical MHC class II gene expression [16]. This is unlikely to be relevant to RA, but less severe changes in efficiency of CIITA expression might be important for autoimmunity development and statistical evaluation of genetic interaction of CIITA and shared epitope alleles may reveal “missing” risk factors.

With this as a background, we set out to define a possible gene-gene interaction between HLA-DRB1 and the CIITA locus in development of RA with a study population of 11767 individuals from four European Caucasian cohorts (6649 RA cases and 5118 controls).

Results

First, we tested the hypothesis of an interactive effect between risk alleles of HLA-DRB1 SE and rs3087456 for developing of RA in the Swedish cohort (Cohort I, Table S1). Interaction was estimated between SE positivity and the risk allele G of rs3087456 in a homozygous state (GG) [8]. The analysis demonstrated no significant evidence of interaction in this model (attributable proportion (AP) = 0.2, 95% CI: −0.2–0.5). Since SE is primarily a risk factor for ACPA positive disease we stratified data according to ACPA status of RA cases. Still, no significant evidence for interaction was found, although a tendency was apparent (AP = 0.3, 95%CI: −0.05–0.6 for ACPA positive status in RA cases, Table 1). Similarly, no significant interaction in additive and multiplicative models was found in the British and the Dutch cohorts. In the Norwegian cohort, however, a significant interaction was detected, both in the total material and in the ACPA positive RA cases (AP = 0.4, 95%CI: 0.03–0.0.7 for RA in total and AP = 0.4, 95%CI: 0.05–0.7 for ACPA positive status in RA cases, Table 1).

thumbnail
Table 1. Risk of developing RA for combinations of the HLA-DRB1 SE and rs3087456 alleles in Swedish.

https://doi.org/10.1371/journal.pone.0032861.t001

To investigate the interaction between HLA-DRB1 SE alleles and rs3087456 in depth, we used a more detailed description of the HLA-DRB1 SE alleles by introducing the allelic groups DRB1*01, DRB1*04 and DRB1*10 as separate risk factors. These analyses did not reveal an SE subgroup allele specific interaction with SNP rs3087456 (Table 2).

thumbnail
Table 2. Summary data of the interaction analysis for HLA-DRB1 SE allelic groups and SNP rs3087456 for the Swedish cohort.

https://doi.org/10.1371/journal.pone.0032861.t002

Since other variants in the CIITA locus have been reported in association with autoimmune disease we genotyped an additional 22 SNPs across the CIITA locus for the Swedish cohort (Chr16: 10842650–10931606, details for RA association tests of these SNPs can be found in Tables S2 and S3). Of these 22 SNPs, including rs3087456, only rs8048002 was significantly associated with RA after correction for multiple comparisons (ACPA negative patients, adjusted for 44 test: p = 0.013, data submitted elsewhere, Eike et al. [17]). Rs3087456 was significantly associated with RA before adjusting for multiple comparisons. To exhaust the possibility of an interaction between CIITA variants and SE, we screened all CIITA SNPs for interaction in the Swedish cohort, using the additive model and two alternative models of dominance: with the minor and major allele of each SNP. In these analyses the SNP rs4781019 showed significant interaction with SE for the ACPA positive subgroup (Table 3, see Table S4). However, the statistical significance did not hold after correction for multiple testing (nominal p = 0.02 for RA in total, nominal p = 0.03 for ACPA positive RA, significance threshold for 44 tests is p = 0.0011, Bonferroni correction). In addition the SE interaction with this SNP could not be confirmed in the independent Norwegian cohort (Table 3.)

thumbnail
Table 3. Summary data for interaction analysis between CIITA rs4781019 and HLA-DRB1 SE.

https://doi.org/10.1371/journal.pone.0032861.t003

Discussion

To assess the combinatorial risk of CIITA and HLA-DRB1 we have investigated the interaction between HLA-DRB1 SE alleles and the CIITA −168A/G polymorphism rs3087456, which was previously found to be associated with RA [8]. It may be that the risk for disease is only detectable in certain combinations and also in certain population. This may be the underlying reason why rs3087456 has been shown as a genetic risk factor for immunological disease in some cohorts [8], [18], [19], but has not been consistently replicated in other cohorts [20], [21], [22]. This was why we set out to define a more specific role of this polymorphism through genetic interaction with HLA-DRB1 SE alleles. However, a straightforward interaction model of SE and the rs3087456 G allele did not reveal significant interaction with regard to RA. In addition, we performed a detailed analysis with the specific SE alleles DRB1*01, DRB1*04, DRB1*10 for a more strict allelic interaction. From this we could conclude that the interaction trend stayed with DRB1*04 but it did not reach statistical significance. We observed that the models with DRB1*01 and DRB1*10 showed negative interaction, which led us to remove individuals with these genotypes from the dataset for a more fair measure of effect from DRB1*04. This analysis resulted in a significant interaction, but was not replicated in the Norwegian cohort (Table S5). A reason for this could be the reduction in size of the dataset and decreased statistical power.

Although the involvement of SNP rs3087456 was the main focus of our study, we also addressed genetic variability in this locus on a broader scale by scrutinizing the CIITA locus for other putative risk markers in the Swedish RA cohort. In a recent article addressing the influence of CIITA and HLA-DRB1 in multiple sclerosis [11], a complex risk relationship between these loci is presented. The polymorphism rs4774 is described as the major risk variant in the CIITA locus instead of rs3087456 and with an increased risk in individuals carrying the DRB1*1501 allele. Rs4774 is also reported to be associated with the production of donor-specific HLA antibodies in renal allograft recipients [23]. Indeed, we found some evidence for interaction with another polymorphism than rs3087456 (rs4781019), but this could not be replicated in the independent Norwegian cohort.

In general, it is difficult to estimate the statistical power for measuring low effects of interaction. Also, except for the Swedish cohort, no particular measures were made to match controls with the RA patients, which may reduce the power of our study and to increase genetic heterogeneity of the study. Therefore, to conclude that an interaction is completely absent between CIITA and HLA-DRB1 SE alleles in development of RA is not possible. However, absence of convincing results in four large, independent cohorts makes it highly unlikely that any strong interaction is present or a sizeable subgroup of the disease could be explained by this interaction.

Our study is not directly comparable to association studies of CIITA in RA, since the major aim was to discover a hypothetical interaction between CIITA and HLA-DRB1 SE. Even so, this and previous studies indicates that CIITA plays an ambiguous role for RA where association signals are difficult to replicate. In a recent article by Eike et al. (unpublished) [17], an updated meta-analysis supports the association of CIITA with RA, in particular in the Scandinavian populations. Also, it seems likely that CIITA is involved in other autoimmune disease, with multiple sclerosis being the most pronounced [11], [18], [24], [25]. According to our data, the previously found association of CIITA variations with RA that could not be replicated in different Caucasian populations is not due to a putative interaction between CIITA and SE alleles. Thus, the missing genetic risk factors for RA remain to be discovered.

To conclude, we did not observe any significant interaction between rs3087456 or 22 other SNPs in the CIITA locus and HLA-DRB1 SE alleles with regard to risk of RA.

Materials and Methods

Description of cohorts

Interaction between HLA-DRB1 SE alleles and rs3087456 in CIITA was primarily investigated in a cohort consisting of 2520 incidence RA cases and 1349 matched controls from Swedish EIRA study, which is described elsewhere [3], [26]. The analysis was repeated in three other cohorts: from UK (1916 cases and 1270 controls), the Netherlands (1260 cases and 346 controls) and Norway (953 cases and 2153 controls) with overall 11767 individuals in the study (6649 RA cases and 5118 controls). All RA patients met the American College for Rheumatology 1987 (ACR-87) revised criteria for RA [27]. For the Swedish cohort, patients were recruited to the study by practioners not responsible for the study, who after informing registered a verbal consent in the patients journal. If consent was given, an extensive questionere was filled in by the patient. Controls were invited by letter and were also asked to fill in and send back an extensive questionere and to visit the closest primary care for leaving samples. This active participation was the foundation for informed consent. All subjects can at any moment withdraw from the study. This procedure is in line with the ethical permit and regulation in Sweden. All individuals in the Dutch, Brittish and Norwegian cohort gave their written informed consent to participate. The local regional ethical review boards approved this study (Regionala etikprövningsnämnden i Stockholm, Sweden; The Leiden institutional review board, Commissie Medische Ethiek, Netherlands; The North-West Multi-Centre Research Ethics Committee and University of Manchester Committee on the Ethics of Research on Human Beings, United Kingdom; Regionale komiteer for medisinsk og helsefaglig forskningsetikk (REK) sør-øst, Oslo, Norway).

Swedish cohort

Genotyping for Swedish cohort was performed by TaqMan predesigned genotyping assay (CIITA) (Applied Biosystems, Foster City, California) and by SSP-PCR (Olerup SSP, Saltsjöbaden, Sweden). The ACPA status was previously identified for all RA patients by Immunoscan RA (Mark 2, cut-off for positivity ≥25 U/ml) enzyme-linked immunosorbent assay (Euro-Diagnostica, Malmö, Sweden).

In the EIRA study, only a minor number of controls were detected to be ACPA positive (1.8%, n = 24), while controls from other cohorts were not tested. When stratifying by ACPA status, controls were always considered as a whole group and not divided in strata.

Dutch cohort

Patient characteristics have been described previously [28]. The healthy controls were randomely selected by the immunogenetics and Transplantation Immunology section of the Leiden University Medical Center.

Genotyping was performed as described for the Swedish cohort. ACPA status was determined by Immunoscan RA (Mark 2, cut-off for positivity ≥25 U/ml) enzyme-linked immunosorbent assay (Euro-Diagnostica, Malmö, Sweden). HLA typing has been described elsewhere [29].

British cohort

All subjects were white Caucasians and all patients satisfied ACR-87 criteria modified for genetic studies [30].

Genotyping for the UK cohort was performed using the Sequenom MassARRAY iPLEX system in accordance with the manufacturer's instructions (www.Sequenom.com). ACPA was tested using the Axis-Shield DIASTAT kit according to the manufacturer's instructions (positivity: concentration >5 U/ml). The presence of HLA-DRB1 SE copy number (0, 1 or 2 copies) was detected using a semi-automated reverse hybridization method (Dynal Biotech, Wirral, UK).

Norwegian cohort

Norwegian RA patients were from the Oslo RA Registry (ORAR) and the European Research on Incapacitating Disease and Social Support (EURIDISS) cohorts [31], [32]. Healthy Norwegian control samples were collected from the Norwegian Bone Marrow Donor Registry (NBMDR), Oslo University Hospital, Rikshospitalet (Controls-1, n = 1121) and blood donors recruited at Oslo University Hospital, Ullevål (Controls-2, n = 1032). An ELISA kit assay (INOVA Diagnostics, San Diego, California, USA) was used to measure ACPA concentrations in the RA samples, with a positivity cut-off defined as levels >25 U/ml. Genotyping of CIITA SNPs in the Norwegian cohort was performed with TaqMan predesigned assays. Genotyping for HLA-DRB1 in the Norwegian RA patients and controls from NBMDR was done by sequence-based genotyping [33], whereas blood donors were genotyped by PCR-based sequence-specific oligonucleotide probe system [34].

Statistical analyses

Additive interaction was defined by departure from additivity of effects originally described by Rothman [35] and was estimated by calculating the attributable proportion due to interaction (AP) [36]. For each individual, variables were defined for having none, either or both risk factors. The statistical tool R was used for logistic regression and estimating ORs for variables and AP was calculated with 95% confidence intervals (version 2.9.0, http://www.r-project.org/). This procedure was facilitated by scripts developed by Kallberg et al., 2006 [37]. Interaction was calculated between SNP rs3087456 (homozygous for G allele) and HLA-DRB1 SE alleles defined as any of DRB1*01, DRB1*04 and DRB1*10. For the other polymorphisms in the CIITA locus both dominant and recessive models for the risk allele were used. For calculating multiplicative interaction and performing allelic association analysis we used the software PLINK (version 1.06, http://pngu.mgh.harvard.edu/purcell/plink/).

We used Bonferroni correction for multiple testing when appropriate.

Supporting Information

Table S1.

Description of the cohorts included in the study.

https://doi.org/10.1371/journal.pone.0032861.s001

(DOC)

Table S2.

Association analyses of the CIITA locus with risk of RA in the Swedish cohort. Pos refers to the genomic position in chromosome 16. MAF is short for minor allele frequency. Plink was used for statistical analysis (http://pngu.mgh.harvard.edu/purcell/plink/). P-values are unadjusted.

https://doi.org/10.1371/journal.pone.0032861.s002

(DOC)

Table S3.

Additional data from association analysis. P values were calculated with trend test in Plink (http://pngu.mgh.harvard.edu/purcell/plink/). P values are unadjusted.

https://doi.org/10.1371/journal.pone.0032861.s003

(DOC)

Table S4.

A summary of interaction analysis for HLA-DRB1 SE with the CIITA locus for the Swedish cohort. * For recessive models the complementary to the recessive risk allele is used for calcualtion due to a low allele frequency. AP = attributable proportion; SE = shared epitope; ACPA = anit citrullinated protein antibodies. Positions with missing data (−) were not possible to calculate.

https://doi.org/10.1371/journal.pone.0032861.s004

(DOC)

Table S5.

Additive interaction between rs3087456 and HLA-DRB1 SE subgroups. # Exclusion of individuals with DBR1*01 or DRB1*10 alleles. Interaction between rs3087456 and DBR1*10 could not be calculated for the Norwegian cohort. AP = attributable proportion; SE = shared epitope; ACPA = anit citrullinated protein antibodies.

https://doi.org/10.1371/journal.pone.0032861.s005

(DOC)

Acknowledgments

We would like to thank Emeli Lundström for HLA typing the Swedish cohort. Sample collection for the Swedish cohort could not have been done without members of EIRA secretariat, Eva Jemseby, Gull-Britt Almgren and Julia Boström.

The Norwegian Bone Marrow Donor Registry is acknowledged for contributing parts of the Norwegian controls. Siri Flåm is thanked for SNP genotyping of the Norwegian data set.

Author Contributions

Conceived and designed the experiments: MR JW REMT LA TK BL LP. Performed the experiments: MR MS MCE DP NAD BS. Analyzed the data: MR MS MCE DP NAD BS LP. Contributed reagents/materials/analysis tools: JW REMT LA TK BL LP. Wrote the paper: MR LP. Critical revising of manuscript: MCE.

References

  1. 1. Klareskog L, Catrina AI, Paget S (2009) Rheumatoid arthritis. Lancet 373: 659–672.
  2. 2. Gregersen PK, Silver J, Winchester RJ (1987) The shared epitope hypothesis. An approach to understanding the molecular genetics of susceptibility to rheumatoid arthritis. Arthritis Rheum 30: 1205–1213.
  3. 3. Klareskog L, Stolt P, Lundberg K, Kallberg H, Bengtsson C, et al. (2006) A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination. Arthritis Rheum 54: 38–46.
  4. 4. Huizinga TW, Amos CI, van der Helm-van Mil AH, Chen W, van Gaalen FA, et al. (2005) Refining the complex rheumatoid arthritis phenotype based on specificity of the HLA-DRB1 shared epitope for antibodies to citrullinated proteins. Arthritis Rheum 52: 3433–3438.
  5. 5. Kallberg H, Padyukov L, Plenge RM, Ronnelid J, Gregersen PK, et al. (2007) Gene-gene and gene-environment interactions involving HLA-DRB1, PTPN22, and smoking in two subsets of rheumatoid arthritis. Am J Hum Genet 80: 867–875.
  6. 6. Padyukov L, Silva C, Stolt P, Alfredsson L, Klareskog L (2004) A gene-environment interaction between smoking and shared epitope genes in HLA-DR provides a high risk of seropositive rheumatoid arthritis. Arthritis Rheum 50: 3085–3092.
  7. 7. Klareskog L, Malmstrom V, Lundberg K, Padyukov L, Alfredsson L (2011) Smoking, citrullination and genetic variability in the immunopathogenesis of rheumatoid arthritis. Semin Immunol 23: 92–98.
  8. 8. Swanberg M, Lidman O, Padyukov L, Eriksson P, Akesson E, et al. (2005) MHC2TA is associated with differential MHC molecule expression and susceptibility to rheumatoid arthritis, multiple sclerosis and myocardial infarction. Nat Genet 37: 486–494.
  9. 9. Bronson PG, Criswell LA, Barcellos LF (2008) The MHC2TA -168A/G polymorphism and risk for rheumatoid arthritis: a meta-analysis of 6861 patients and 9270 controls reveals no evidence for association. Ann Rheum Dis 67: 933–936.
  10. 10. Bronson PG, Ramsay PP, Seldin MF, Gregersen PK, Criswell LA, et al. (2011) CIITA is not associated with risk of developing rheumatoid arthritis. Genes Immun 12: 235–238.
  11. 11. Bronson PG, Caillier S, Ramsay PP, McCauley JL, Zuvich RL, et al. (2010) CIITA variation in the presence of HLA-DRB1*1501 increases risk for multiple sclerosis. Hum Mol Genet 19: 2331–2340.
  12. 12. Koizumi K, Okamoto H, Iikuni N, Nakamura T, Kawamoto M, et al. (2005) Single nucleotide polymorphisms in the gene encoding the major histocompatibility complex class II transactivator (CIITA) in systemic lupus erythematosus. Ann Rheum Dis 64: 947–950.
  13. 13. Skinningsrud B, Husebye ES, Pearce SH, McDonald DO, Brandal K, et al. (2008) Polymorphisms in CLEC16A and CIITA at 16p13 are associated with primary adrenal insufficiency. J Clin Endocrinol Metab 93: 3310–3317.
  14. 14. Reith W, Mach B (2001) The bare lymphocyte syndrome and the regulation of MHC expression. Annu Rev Immunol 19: 331–373.
  15. 15. Ting JP, Trowsdale J (2002) Genetic control of MHC class II expression. Cell 109: SupplS21–33.
  16. 16. Steimle V, Otten LA, Zufferey M, Mach B (1993) Complementation cloning of an MHC class II transactivator mutated in hereditary MHC class II deficiency (or bare lymphocyte syndrome). Cell 75: 135–146.
  17. 17. Eike MC, Skinningsrud B, Ronninger M, Stormyr A, Kvien TK, et al. (2012) CIITA gene variants are associated with rheumatoid arthritis in Scandinavian populations: an association study and meta-analysis. Manuscript submitted for publication.
  18. 18. Martinez A, Sanchez-Lopez M, Varade J, Mas A, Martin MC, et al. (2007) Role of the MHC2TA gene in autoimmune diseases. Ann Rheum Dis 66: 325–329.
  19. 19. Iikuni N, Ikari K, Momohara S, Tomatsu T, Hara M, et al. (2007) MHC2TA is associated with rheumatoid arthritis in Japanese patients. Ann Rheum Dis 66: 274–275.
  20. 20. Eyre S, Bowes J, Spreckley K, Potter C, Ring S, et al. (2006) Investigation of the MHC2TA gene, associated with rheumatoid arthritis in a Swedish population, in a UK rheumatoid arthritis cohort. Arthritis Rheum 54: 3417–3422.
  21. 21. Harrison P, Pointon JJ, Farrar C, Harin A, Wordsworth BP (2007) MHC2TA promoter polymorphism (−168*G/A, rs3087456) is not associated with susceptibility to rheumatoid arthritis in British Caucasian rheumatoid arthritis patients. Rheumatology (Oxford) 46: 409–411.
  22. 22. Orozco G, Robledo G, Linga Reddy MV, Garcia A, Pascual-Salcedo D, et al. (2006) Study of the role of a functional polymorphism of MHC2TA in rheumatoid arthritis in three ethnically different populations. Rheumatology (Oxford) 45: 1442–1444.
  23. 23. Martin J, Worthington J, Harris S, Martin S (2009) The influence of class II transactivator and interleukin-6 polymorphisms on the production of antibodies to donor human leucocyte antigen mismatches in renal allograft recipients. Int J Immunogenet 36: 235–239.
  24. 24. Martinez A, Alvarez-Lafuente R, Mas A, Bartolome M, Garcia-Montojo M, et al. (2007) Environment-gene interaction in multiple sclerosis: human herpesvirus 6 and MHC2TA. Hum Immunol 68: 685–689.
  25. 25. Zuvich RL, Bush WS, McCauley JL, Beecham AH, De Jager PL, et al. (2011) Interrogating the complex role of chromosome 16p13.13 in multiple sclerosis susceptibility: independent genetic signals in the CIITA-CLEC16A-SOCS1 gene complex. Hum Mol Genet 20: 3517–3524.
  26. 26. Stolt P, Kallberg H, Lundberg I, Sjogren B, Klareskog L, et al. (2005) Silica exposure is associated with increased risk of developing rheumatoid arthritis: results from the Swedish EIRA study. Ann Rheum Dis 64: 582–586.
  27. 27. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, et al. (1988) The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 31: 315–324.
  28. 28. Kurreeman FA, Padyukov L, Marques RB, Schrodi SJ, Seddighzadeh M, et al. (2007) A candidate gene approach identifies the TRAF1/C5 region as a risk factor for rheumatoid arthritis. PLoS Med 4: e278.
  29. 29. van der Woude D, Lie BA, Lundstrom E, Balsa A, Feitsma AL, et al. (2010) Protection against anti-citrullinated protein antibody-positive rheumatoid arthritis is predominantly associated with HLA-DRB1*1301: a meta-analysis of HLA-DRB1 associations with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in four European populations. Arthritis Rheum 62: 1236–1245.
  30. 30. MacGregor AJ, Bamber S, Silman AJ (1994) A comparison of the performance of different methods of disease classification for rheumatoid arthritis. Results of an analysis from a nationwide twin study. J Rheumatol 21: 1420–1426.
  31. 31. Kvien TK, Uhlig T (2003) The Oslo experience with arthritis registries. Clin Exp Rheumatol 21: S118–122.
  32. 32. Guillemin F, Gerard N, van Leeuwen M, Smedstad LM, Kvien TK, et al. (2003) Prognostic factors for joint destruction in rheumatoid arthritis: a prospective longitudinal study of 318 patients. J Rheumatol 30: 2585–2589.
  33. 33. Sayer DC, Whidborne R, De Santis D, Rozemuller EH, Christiansen FT, et al. (2004) A multicenter international evaluation of single-tube amplification protocols for sequencing-based typing of HLA-DRB1 and HLA-DRB3,4,5. Tissue Antigens 63: 412–423.
  34. 34. Erlich H, Valdes AM, Noble J, Carlson JA, Varney M, et al. (2008) HLA DR-DQ haplotypes and genotypes and type 1 diabetes risk: analysis of the type 1 diabetes genetics consortium families. Diabetes 57: 1084–1092.
  35. 35. Rothman KJ (2002) Epidemiology: an introduction. New York, N.Y.: Oxford University Press. viii, 223 p.
  36. 36. Andersson T, Alfredsson L, Kallberg H, Zdravkovic S, Ahlbom A (2005) Calculating measures of biological interaction. Eur J Epidemiol 20: 575–579.
  37. 37. Kallberg H, Ahlbom A, Alfredsson L (2006) Calculating measures of biological interaction using R. Eur J Epidemiol 21: 571–573.