Figures
Abstract
Multiple sclerosis (MS) is a complex disease of the central nervous system of unknown etiology. The human leukocyte antigen (HLA) locus on chromosome 6 confers a considerable part of the susceptibility to MS, and the most important factor is the class II allele HLA-DRB1*15:01. In addition, we and others have previously established a protective effect of HLA-A*02. Here, we genotyped 1,784 patients and 1,660 healthy controls from Scandinavia for the HLA-A, HLA-B, HLA-C and HLA-DRB1 genes and investigated their effects on MS risk by logistic regression. Several allele groups were found to exert effects independently of DRB1*15 and A*02, in particular DRB1*01 (OR = 0.82, p = 0.034) and B*12 (including B*44/45, OR = 0.76, p = 0.0028), confirming previous reports. Furthermore, we observed interaction between allele groups: DRB1*15 and DRB1*01 (multiplicative: OR = 0.54, p = 0.0041; additive: AP = 0.47, p = 4×10−06), DRB1*15 and C*12 (multiplicative: OR = 0.37, p = 0.00035; additive: AP = 0.58, p = 2.6×10−05), indicating that the effect size of these allele groups varies when taking DRB1*15 into account. Analysis of inferred haplotypes showed that almost all DRB1*15 bearing haplotypes were risk haplotypes, and that all A*02 bearing haplotypes were protective as long as they did not carry DRB1*15. In contrast, we found one class I haplotype, carrying A*02-C*05-B*12, which abolished the risk of DRB1*15. In conclusion, these results confirms a complex role of HLA class I and II genes that goes beyond DRB1*15 and A*02, in particular by including all three classical HLA class I genes as well as functional interactions between DRB1*15 and several alleles of DRB1 and class I genes.
Citation: Link J, Kockum I, Lorentzen ÅR, Lie BA, Celius EG, Westerlind H, et al. (2012) Importance of Human Leukocyte Antigen (HLA) Class I and II Alleles on the Risk of Multiple Sclerosis. PLoS ONE 7(5): e36779. https://doi.org/10.1371/journal.pone.0036779
Editor: Jeong-Sun Seo, Seoul National University College of Medicine, Republic of Korea
Received: August 1, 2011; Accepted: April 6, 2012; Published: May 7, 2012
Copyright: © 2012 Link et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was funded by the Bibbi & Nils Jensen Foundation, the Swedish Research Council, the Swedish Council for Working Life and Social Research, Karolinska Institutet, and the 6th Framework Program of the European Union, NeuroproMiSe (LSHM-CT-2005-018637). The Norwegian work was supported by the University of Oslo, the Research Council of Norway (166005/V5), South and Eastern Regional Health Authority, Ullevål University Hospital Scientific Advisory Council (UUHSAC) and Odd Fellow MS Society. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
In the past few years, the number of genes known to influence the risk of multiple sclerosis (MS, OMIM:126200), a chronic inflammatory disease of the central nervous system, has increased dramatically. Even so, the by far most strongly acting genetic signal originates from the human leukocyte antigen (HLA) gene complex on the short arm of chromosome 6. The discovery of the class II risk haplotype, later shown to be best represented by the HLA-DRB1*15:01 allele, was made almost 40 years ago and was only recently shown to be accompanied by an independent protective allele group in class I, HLA-A*02 [1], [2], [3], [4], [5], [6]. In addition to these two established alleles a series of reports have suggested additional MS risk alleles as well as interactions between alleles [7]. Within the HLA-DRB1 gene, the allele groups DRB1*03, DRB1*14, DRB1*07 and DRB1*11 have all been suggested to influence risk in MS [8], [9], [10], [11]. Furthermore, interactions have been implicated among DRB1 allele groups: DRB1*08, DRB1*10 and DRB1*01 were reported to interact with DRB1*15 in modifying susceptibility to MS [8], [10].
An increasing number of studies has recently confirmed the importance of genes in the class I region on MS susceptibility [12], [13]. The HLA-A*02 association was first identified in a limited Swedish cohort [2] and was later confirmed in a larger Scandinavian cohort [1]. In addition, studies including not only HLA-DRB1 and HLA-A, but also HLA-B and/or HLA-C have indicated additionally associated alleles [14], [15], [16], [17], including HLA-C*08, HLA-C*05 and HLA-B*44. A single nucleotide polymorphism (SNP) genotyping approach, where HLA alleles were imputed based on linkage disequilibrium with SNP markers, found a signal of association from HLA-B*44:02 when conditioning on DRB1*15 [18]. HLA-B*44 was also found associated independently from DRB1*15 in a meta-analysis of available SNP genome scans [19]. The MOG gene, encoding the myelin oligodendrocyte glycoprotein and positioned telomeric of HLA-A, has been repeatedly studied [20], [21] but its signal of association seems most likely secondary to HLA-A*02 [15].
Several studies have tried to distinguish which of the class I alleles or genes that are responsible for the class I region signal. When carried together, HLA-C*05 was reported to enhance the protective effect of HLA-A*02 [15], but this could not be replicated in the Scandinavian population [17]. HLA-C*05 has been seen to be secondary to both HLA-A*02 and HLA-B*44 whereas B*44 was not secondary to HLA-A*02 [22]. In addition to this, a SNP screen revealed a SNP associated to MS in the untranslated region of HLA-G that could not be ruled out to be independent from either HLA-A, -B or -C alleles [23].
In this study, we aimed at assessing signals from the three classical class I genes together with HLA-DRB1 in a large cohort consisting of Swedish and Norwegian MS patients and controls. All three major genes in class I as well as HLA-DRB1 have been genotyped with robust direct HLA genotyping techniques and the association with MS and potential interactions were assessed with logistic regression.
Results
The aim of this study was to investigate the association to MS of all three classical HLA class I genes, HLA-A, HLA-B, HLA-C as well as the HLA class II gene, HLA-DRB1, in a Scandinavian cohort of 1,784 patients and 1,660 controls. The nature of the previously reported effects within class I in MS was also evaluated.
HLA class I and II allele groups in MS risk
We started by assessing the mode of inheritance of the two already established genetic factors within the HLA region to determine how to proceed with the analysis. HLA-DRB1*15 was seen to behave multiplicatively as the odds ratio for carrying one copy of DRB1*15 (95% confidence interval) was 3.03 (2.62–3.52) and 11.47 (7.74–16.99) for carrying two copies. In contrast HLA-A*02 behaved dominantly, the odds ratio for carrying one copy of A*02 was 0.70 (0.61–0.81) and 0.57 (0.45–0.71) for carrying two copies. Therefore these particular allele groups were encoded accordingly throughout this study (see material and methods section for details).
There was no difference between coding the additional allele groups as gene dose (0, 1 or2) or carriage (0 or 1) when assessing which allele groups should be in the final regression model (data not shown), but the overall fit of the model (deviance) was superior in the dominant model. Therefore only carriage coded results will be presented. In the final model DRB1*15 and B*14 were positively associated with MS while A*02, B*12, DRB1*01 and DRB1*07 were negatively associated (Table 1). Thus, HLA-B seemed to influence MS risk in addition to HLA-A and HLA-DRB1 but HLA-C was not independently associated with MS at this stage.
Next, we wanted to assess whether groups stratified by DRB1*15 carriage display different genetic etiology. 1,051 cases and 476 controls were DRB1*15 positive while 733 cases and 1,184 controls were negative for DRB1*15 in this cohort. Allele group frequencies are shown in Table S1. Six allele groups reached significance in the HLA-DRB1*15 negative group whereas four allele groups reached significance in the DRB1*15 positive group (Table 2). Surprisingly, apart from A*02, which was negatively associated in both strata, the other allele groups were different between the groups. In the negative group we found associations to B*14, A*02, DRB1*07, B*12, A*10 and B*18 and in the positive group both DRB1*15 and A*02 as well as C*12 and DRB1*01 were associated (Table 2). An importance of DRB1*01 in DRB1*15 positives has previously been reported [8] and was thus confirmed. In addition, even allele groups of HLA-C may influence the risk of MS, at least in the DRB1*15 positive individuals.
The difference in association pattern between the DRB1*15 positive and negative subpopulation prompted us to formally assess the presence of interaction. The allele groups DRB1*07, B*12, B*18, A*10 and B*14, failed to show a significant interaction with DRB1*15. However two allele groups (C*12 and DRB1*01) did reach significant statistical interaction when using either departure from multiplicativity as definition, or when analyzed by attributable proportion due to interaction (AP) (i.e. departure from additivity) (Table 3). In addition, A*02 showed statistically significant interaction with DRB1*15 in AP analysis only.
In conclusion, we observed significant support for interaction for MS risk between HLA class I and class II genes as well as between allele groups of HLA-DRB1 itself. This confirms and extends previous reports of such interactions [8].
Evaluation of protective effects of class I allele groups in MS
Previous reports have indicated HLA-C*05 and HLA-B*44 as important allele groups in MS, and Healy et al used logistic regression to study these allele groups with the previously established HLA-DRB1*15 and HLA-A*02 as covariates [22]. It has also been reported that the protective effect of A*02 is enhanced when also carrying C*05 in the Italian population [14]. In a previous study, we failed to confirm this role of C*05 [17]. Now, when being able to study also B*12, the protective effect of HLA-C*05 was seen to be secondary to both B*12 and HLA-A*02 due to linkage disequilibrium (Table S2). Again, we failed to confirm an augmented protective effect of A*02 when carrying C*05 (A*02 alone OR = 0.66, p = 5.95×10−09; carrying A*02 and C*05 OR = 0.64, p = 2.65×10−04).
In recent years, SNP genotyping has widely been applied in attempts to map genetic effects within the HLA gene region. Early attempts failed to find much residual effects after stratifying for DRB1*15 [13]. More recently, Cree et al reported an association to the SNP rs4959039 (A/G), located in the downstream untranslated region of HLA-G [23]. This SNP was found to be associated independently of a DRB1*15 tagging SNP but it could not be determined if it tagged HLA-A*02 or was independent of this allele group. We genotyped this marker in 1,978 patients and 1,782 controls previously genotyped for HLA-A and HLA-DRB1 and observed a minor allele frequency (G) of 0.26 in cases and 0.33 in controls. A haplotype analysis was performed using UNPHASED to elucidate the correlation between rs4959039 and HLA-A. The correlation between minor allele rs4959039G and HLA-A*02 was clearly higher compared to the other top five HLA-A allele groups with a D′ of 0.90 and r2 of 0.78. There was no residual association with the SNP when correcting for A*02 (p = 0.90) in logistic regression, convincingly showing that the rs4959039 effect was secondary to A*02.
Haplotype analysis
Allele group association analysis assesses the combined risk of all the haplotypes carrying certain allele groups. The effects of associated HLA allele groups in logistic regression are independent of each other but a shown risk or protection can, alternatively, indicate a stronger effect at the haplotype level. The property of the HLA gene complex to harbor these extended conserved haplotypes often spanning several megabases, may be seen as an argument against case-control design where phase cannot be deduced. To compensate for this void, we attempted to infer haplotypes using UNPHASED 3.0.13, hypothesizing that the greater available numbers in this case control study compared with available family cohorts could outweigh the disadvantages. Table 4 contains the 20 most frequent inferred haplotypes from this analysis. Only haplotypes with a probability of at least 80% were used in the association analysis. Similar results were obtained when all haplotypes, irrespectively of probability, were included (data not shown).
The haplotype analysis revealed several nominally significant associated haplotypes in spite of low frequencies. Interestingly, almost all of the most common haplotypes contain one or more allele groups found to be associated in the final full logistic regression model (Table 1). A striking finding was that almost all DRB1*15 carrying haplotypes in Table 4 were risk haplotypes and all A*02 carrying haplotypes were protective as long as they did not carry DRB1*15. There was however an interesting exception (haplotype number 12 in Table 4) where A*02 together with C*05 and B*12 abolished the risk of DRB1*15. To further investigate these allele groups, logistic regressions with the 20 most common haplotypes carrying each of these allele groups was created (Table S3, S4, S5, S6). Indeed almost all DRB1*15 haplotypes were risk haplotypes even though not all reached significance for association with MS (Table S3). The exception was again the haplotype carrying A*02, C*05 and B*12, OR (95% C.I.) = 1.13 (0.62–2.07) and p-value = 0.71 thus somewhat lowering the risk of DRB1*15. The same analysis was done for A*02 bearing haplotypes (Table S4), and also here the pattern from Table 4 was repeated. Most of the B*12 haplotypes were protective unless carrying DRB1*15 or DRB1*06 (Table S5). When looking at the C*05 bearing haplotypes, C*05 was only seen in haplotypes with either B*12, B*18 or B*05 and varied in risk for MS (Table S6). One should bear in mind that these haplotypes were indeed estimated and the haplotype frequencies were small.
The protective haplotype carrying A*02, C*05 and B*12 (hereafter referred to as “protective class I haplotype” or abbreviated as “H”) was further studied for the nature of its protective properties in relation to DRB1*15 and if any of its alleles were sole responsible for the protective effect. The protective class I haplotype was independent of DRB1*15 and was protective on both DRB1*15 positive and DRB1*15 negative haplotypes (Table 5).
Chromosomes positive for both A*02 and B*12 but lacking C*05 was the only combination of the protective class I haplotype that did not have a significantly higher risk than chromosomes carrying all three allele groups (H+)(Figure 1). Thereby indicating that a haplotype containing at least A*02 and B*12 seems responsible for this protective effect in our cohort.
Odds ratios for haplotypes containing one or two of the allele groups of the protective class I haplotype compared to the full protective class I haplotype disregarding the DRB1 locus. Chromosomes carrying both A*02 and B*12 was the only combination within this setting that had the same risk as the protective haplotype, the risk of MS was significantly higher for all other combinations.“+” = Chromosomes carrying the allele group specified to the left, “H−” = chromosomes carrying other allele group than the one specified to the left (i.e. not the complete protective haplotype). OR = odds ratio with 95% confidence interval, nominal p-value and Bonferroni corrected p-values derived from chi square test with full protective haplotype as reference.
When studying the effect of the protective class I haplotype in the context of DRB1*15, it was seen that compared to doubly negative chromosomes (“H− DRB1*15-” neither carrying the protective class I haplotype nor DRB1*15), chromosomes carrying DRB1*15 but negative for protective class I haplotype had an OR of 3 (p-value 1.08×10−63) (Figure 2). Chromosomes carrying the protective class I haplotype but not DRB1*15 had a negative association compared to doubly negative chromosomes (OR 0.42, p-value 3.42×10−06). Chromosomes carrying both A*02, C*05, B*12 and DRB1*15 did not differ in risk from doubly negatives (OR 1.25, p-value 1). Thus, this haplotype diminishes the risk of DRB1*15 also in this analysis.
Chromosomes carrying both DRB1*15 and the protective class I haplotype do not differ in risk of MS from chromosomes not carrying either of them (reference). DRB1*15 carriage increase risk of MS significantly and carriage of the class I haplotype significantly decrease risk of MS. Counts are measured in chromosomes and only individuals with estimated haplotypes with a probability of more than 80% were included. “+” = Chromosomes positive for the protective haplotype or DRB1*15, “−” and “H−” = chromosomes negative for DRB1*15 or the protective haplotype. OR = odds ratio with 95% confidence interval, nominal p-value and Bonferroni corrected p-values derived from chi square test with doubly negative chromosomes as reference.
Analysis of the haplotype was also performed with logistic regression to be able to account for carriage of A*02, C*05 and B*12 as well as baseline covariates as country of origin, gender and carriage of DRB1*15. After including the three class I allele groups and adding carriage of the haplotype into the model both A*02 (OR 0.73 p = 0.00015), B*12 (OR 0.69 p = 0.0032) and the protective haplotype (OR 0.54 p = 0.0080) were significantly associated to MS. Carriage of C*05 was not significant and DRB1*15 was independently associated to MS in the model (OR 3.28 p = <2×10−16). The overall fit (deviance) of the model was significantly better than the model with only baseline variables and carriage of the three class I allele groups (p = 0.0076).
Discussion
The aim of this study was to investigate the class I and class II allele groups and their effect on MS susceptibility in the Scandinavian population. All 3,444 individuals were genotyped for HLA-A, -B, -C, and -DRB1.
We have confirmed and extended previous reports of a complex genetic influence of the HLA region on MS susceptibility. The two strongest effects are clearly conferred by DRB1*15 and A*02 but we could also show that additional DRB1 allele groups contribute to risk as well as allele groups of HLA-B which seem more important than those of HLA-C. Further, several alleles were interacting with DRB1*15, exerting their effects either in DRB1*15 positive or negative individuals. However, it was clear that no other class I allele group equaled the effect of A*02.
A novelty of this study is the application of formal interaction analysis. We could see that several of the associations from the full analysis persisted in only one of the two DRB1*15 strata. Two HLA allele groups did interact with DRB1*15, namely C*12 and DRB1*01. A protective effect of DRB1*01 was seen both in the full sample set and in the DRB1*15 positives. This trans interaction has been discussed previously where an under transmission of DRB1*01 was detected in DRB1*15 positive affected children in family trios [8]. How such an interaction could be explained molecularly is potentially important as it may give clues on how HLA genotype exerts its role on MS risk. Explanations could potentially involve competition for peptides or trans-coded HLA class II heterodimers [24]. The latter mechanism would be more likely to occur at the HLA-DQ locus where both α- and β-chains are polymorphic (in contrast to the almost non-polymorphic DRα-chain).
The interaction between DRB1*15 and A*02 and DRB1*15 and C*12 suggests that they are involved in the same pathway of developing disease but not necessarily exerting their effect at the same time or interact physically since they act mainly in different immune pathways. This data supports the hypothesis that viral infections play an important role in susceptibility to MS.
Traditionally, it has been argued that the existence of extended haplotypes spanning HLA class I and class II indicates a functional property being evolutionarily advantageous and potentially important for the risk of autoimmune disease. As shown in Table 4, all haplotypes were of limited frequency and none of the most common ones showed stronger risk than the most clearly associated allele groups. Haplotypes carrying DRB1*15 were almost always risk haplotypes regardless of which other allele groups that were present, indicating that it is probably the DRB1*15 allele itself that confers the risk of MS. The protective class I haplotype suggests a protection against MS that must be further studied to evaluate two things: 1) is this haplotype tagging some other marker in the HLA or 2) is the protection inflated by combining several independent protective markers that when carried together overcome the effect of DRB1*15.
Recently, the International MS Genetics Consortium (IMSGC) published an analysis of SNP-derived inferred HLA genotypes as part of a genome wide association study on close to 10.000 patients [25]. Also here, a dominance of a DRB1*15 effect was followed by an influence of A*02. In addition, residual effects were explained by other DRB1 alleles and little further signal in class I was detected which would be in line with our suggestion of the protective class I haplotype. Besides from harboring the classic HLA genes, the HLA gene complex is rich in other genes of importance for immune function. It is plausible that such genes are of importance also in MS but the proximity of the two most strongly acting MS genes is likely to obscure such signals making them difficult to identify. In essence the main result of the IMSGC screen of the HLA gene complex is that the strongest effects indeed map to the classic HLA genes themselves, both for the class II and class I subregions.
The confirmation of an importance for one or even several HLA class I genes in MS, paralleled by findings in other autoimmune diseases [18], points toward a possible importance of class I restricted CD8+ T-cells. Our group has previously studied this in the MS-like model experimental autoimmune encephalomyelitis (EAE) after immunization with myelin basic protein (MBP). Here, protection against disease was found in rats carrying the MHC class I u and MHC class II a haplotype, an effect efficiently abrogated by antibodies against CD8+ cells [26]. A series of studies using double and triple transgenic mice has further illustrated how class I restricted lymphocytes may contribute to MS-like autoimmune events [27], [28]. The explanation could be that a combination of class I allele groups, such as A*02-C*05-B*44, immunologically counteract the negative effect of carrying DRB1*15 in class II when encountering some pathogen involved in MS susceptibility. This mechanism could vary in effect size between populations since the environmental challenge may be different in different areas.
In conclusion, evidence is suggesting that HLA class I and class II effects contribute and interact in determining the risk of MS, as in autoimmune diseases in general. These genes are the most strongly associated with MS and the functional background to these effects should be given increased attention.
Materials and Methods
Subjects
All 1784 patients were diagnosed with MS according to the McDonald or Poser criteria [29], [30]. The 1,229 Swedish patients were gathered through the neurology clinics at Karolinska University Hospital, Stockholm, Sweden. The ethical board of Karolinska Institutet approved of the study and oral or written informed consent was obtained from all participants according to the ethical permission. From Norway 555 MS patients were collected from the Oslo MS registry at the Department of Neurology, Ullevål, Oslo University Hospital. The study was approved by the Regional Committee for Medical and Health Research Ethics, and written informed consent was obtained from all individuals.
Of the 1,660 healthy controls used in this study, 1,062 were Swedish blood donors of Nordic origin giving informed consent. From Norway 598 healthy controls were recruited from the Norwegian Bone Marrow Donor Registry, (http://www.nordonor.org/). Age and gender distribution for cases and controls are shown in Table 6.
Methods
HLA genotyping was performed using several techniques. The Swedish cohort was genotyped for HLA-DRB1, HLA-A and HLA-C using Olerup SSP™ HLA Low resolution Kit [31] and for HLA-B using a Luminex based reverse PCR-SSO (LABType® SSO from One Lambda, Inc., Canoga Park, CA, USA). The Norwegian patients were genotyped using PCR with sequencing-based approach [32] and the controls were genotyped with various techniques used in the Norwegian Bone Marrow Donor registry laboratory, (http://www.nordonor.org/).
SNP genotyping and analysis
Swedish MS cases (N = 2,576) and controls (N = 2,031) were genotyped in the study for correlation analysis between a SNP (rs4959039) in the HLA region and HLA-A*02 status. This group of individuals consisted of the above described patients and controls plus newly diagnosed patients and matched controls from the EIMS study described elsewhere [33]. SNP genotyping was performed using the TaqMan® SNP Genotyping Assays (Applied Biosystems, USA) according to the manufacturer's manual. HLA-A and HLA-DRB1 genotyping was formerly performed using Olerup SSP™ HLA Low resolution Kit [31]. In total, 1,978 patients and 1,782 controls had full information on gender, the SNP, HLA-A and HLA-DRB1 genotypes and were used in the correlation analysis.
Serological allele groups
As some Norwegian controls were not genotyped down to two digit level, major allele groups were created where there was an information shortage. The allele grouping is shown in Table 7. It should be noted that the B*44 allele group which has previously been reported to be associated to MS is included in the B*12 major allele group together with B*12 and B*45. However, B*44 is the dominating allele within this group (90%).
Statistical methods
The software used for statistical calculations in this paper was R version 2.6.2 [34] and UNPHASED 3.0.13 [35]. Nested logistic regression models were compared using analysis of deviance and chi square test in R. In the logistic regressions all allele groups were coded 0, 1, 2 (dose) or 0, 1 (carrier) depending on question and those allele groups with a total count of less than 150 chromosomes (2.2% allele frequency) in the full cohort were grouped within each gene into a group called X (see Table 7). This cutoff was set to decrease complexity of the cohort and increase power. DRB1*08, A*11, B*35 and C*02 were used as intercepts for the full analysis as these were the most evenly distributed between cases and controls and had a frequency of at least 5% in both groups. First, a baseline model including the already known risk allele groups for MS was created. This model included country of origin, gender, HLA-DRB1*15 and HLA-A*02. Allele groups within each gene were added to this baseline model to determine if any allele group had an additional independent effect. Within each model only those allele groups with a statistically significant odds ratio were kept in the next step in which a new model was created containing the baseline plus the significant allele groups from all genes. Allele groups were removed using stepwise regression until only significant ones were left. The model with the least number of included allele groups was compared to baseline. This analysis was performed twice with all allele groups coded as dose (0, 1, 2) or as carriage (0, 1) (in previous literature commonly referred to as “phenotype” referring to the co-dominant expression of HLA genes), with the exceptions being DRB1*15 which was previously seen to act multiplicatively and was therefore always coded as dose and A*02 which had a dominant effect in our cohort, hence always coded for carriage as described above. When stratifying the cohort into DRB1*15 positive or negative strata, all allele groups with a total count of less than 75 (2.5%) or 100 (2.6%) chromosomes, respectively, were grouped as an X allele within each gene (see Table S1). When stratifying for DRB1*15 status, A*19, B*15, C*03 and DRB1*04 were used as intercepts in the positive group and A*01, B*07, C*04 and DRB1*04 were used as intercept in the DRB1*15 negative group . DRB1*15 was coded as dose (i.e. 2, 1 or 0 alleles) in the DRB1*15 positive group, all other alleles were coded as carriage (i.e. 1 or 0). The significant allele groups did not change in either stratum when using the same intercept as in the full analysis (data not shown).
Allele group associations were also analyzed for interaction using deviation from multiplicativity of effects (i.e. significant interaction term in the logistic regression model) as well as using deviation from additivity of effects (estimated by attributable proportion (AP) due to interaction) as criteria for interaction [36], [37]. The deviation from multiplicativity of effects analysis was performed with 5 variables: country of origin, gender, DRB1*15 allele group carriage status, allele group carriage status and an interaction term between DRB1*15 and the allele group investigated. The individuals were coded as carriers or non carriers of the particular risk allele group as well as for DRB1*15. In the AP analyses, risk was defined as “not carrying the protective allele group” for those allele groups whose nominal effects were protective.
Haplotypes for individuals were estimated using UNPHASED 3.0.13 since this version of the program could handle larger number of alleles and the tolerance for convergence was set to 10−6 due to complexity of the matrix, only printing most likely haplotypes to file. Only individuals with a probability over 80% of carrying the assigned haplotypes were used in the subsequent analyses leaving2,886 individuals of 3,444 in the analysis and excluding 48 haplotypes out of 1,279.
There was no difference in allele group distribution in different age groups or between males and females (data not shown). Gender was included as a covariate anyway since gender is a risk factor for MS in itself. Country of origin was included as a covariate to exclude possible hidden effects from population stratification even though both populations were of Scandinavian origin.
The correlation between the genotyped SNP (rs4959039) and HLA-A allele groups was estimated with UNPHASED 3.0.13 as for the haplotype analysis, but the complexity did not require changing the tolerance for convergence in this case. To study the effect of the SNP on MS susceptibility, logistic regression was performed in R.
Supporting Information
Table S1.
Allele group frequencies in the DRB1*15 positive and DRB1*15 negative strata.
https://doi.org/10.1371/journal.pone.0036779.s001
(DOC)
Table S2.
Associations of each locus in regression analysis of all subjects (n = 3444). There is no difference in MS risk between carriage of only A*02 or both A*02 and C*05. All investigated allele groups except C*05 show robust association with MS. Country of origin and gender were included in the models as covariates
https://doi.org/10.1371/journal.pone.0036779.s002
(DOC)
Table S3.
Frequencies of estimated haplotypes, odds ratios and p-values from logistic regression with the 20 most common DRB1*15 carrying haplotypes.
https://doi.org/10.1371/journal.pone.0036779.s003
(DOC)
Table S4.
Frequencies of estimated haplotypes, odds ratios and p-values from logistic regression with the 20 most common A*02 carrying haplotypes.
https://doi.org/10.1371/journal.pone.0036779.s004
(DOC)
Table S5.
Frequencies of estimated haplotypes, odds ratios and p-values from logistic regression with the 20 most common B*12 carrying haplotypes.
https://doi.org/10.1371/journal.pone.0036779.s005
(DOC)
Table S6.
Frequencies of estimated haplotypes, odds ratios and p-values from logistic regression with the 20 most common C*05 carrying haplotypes.
https://doi.org/10.1371/journal.pone.0036779.s006
(DOC)
Acknowledgments
The authors would like to thank Emilie Sundqvist and Magdalena Lindén for help with HLA genotyping of the EIMS cohort. The Norwegian Bone Marrow donor registry is acknowledged for contribution of the Norwegian controls included in the study, and the Oslo MS genetics group is thanked for contributions to the Norwegian MS data set. The Swedish neurology clinics and all patients as well as the Swedish blood donors are thanked for their contributions to this study.
Author Contributions
Conceived and designed the experiments: JL IK BB HH JH. Performed the experiments: JL IK ÅL BL EC HW MS LA TO BB HH JH. Analyzed the data: JL IK BB JH. Contributed reagents/materials/analysis tools: JL IK ÅL BL EC HW MS LA TO BB HH JH. Wrote the paper: JL IK ÅL BL EC HW MS LA TO BB HH JH.
References
- 1. Brynedal B, Duvefelt K, Jonasdottir G, Roos IM, Akesson E, et al. (2007) HLA-A confers an HLA-DRB1 independent influence on the risk of multiple sclerosis. PloS one 2: e664.
- 2. Fogdell-Hahn A, Ligers A, Gronning M, Hillert J, Olerup O (2000) Multiple sclerosis: a modifying influence of HLA class I genes in an HLA class II associated autoimmune disease. Tissue antigens 55: 140–148.
- 3. Harbo HF, Lie BA, Sawcer S, Celius EG, Dai KZ, et al. (2004) Genes in the HLA class I region may contribute to the HLA class II-associated genetic susceptibility to multiple sclerosis. Tissue Antigens 63: 237–247.
- 4. Jersild C, Fog T (1972) Histocompatibility (HL-A) antigens associated with multiple sclerosis. Acta Neurol Scand Suppl 51: 377.
- 5. Jersild C, Fog T, Hansen GS, Thomsen M, Svejgaard A, et al. (1973) Histocompatibility determinants in multiple sclerosis, with special reference to clinical course. Lancet 2: 1221–1225.
- 6. Naito S, Namerow N, Mickey MR, Terasaki PI (1972) Multiple sclerosis: association with HL-A3. Tissue Antigens 2: 1–4.
- 7. Hoppenbrouwers IA, Hintzen RQ (2011) Genetics of multiple sclerosis. Biochim Biophys Acta 1812: 194–201.
- 8. Dyment DA, Herrera BM, Cader MZ, Willer CJ, Lincoln MR, et al. (2005) Complex interactions among MHC haplotypes in multiple sclerosis: susceptibility and resistance. Hum Mol Genet 14: 2019–2026.
- 9. Hillert J, Gronning M, Nyland H, Link H, Olerup O (1992) An immunogenetic heterogeneity in multiple sclerosis. J Neurol Neurosurg Psychiatry 55: 887–890.
- 10. Ramagopalan SV, Morris AP, Dyment DA, Herrera BM, DeLuca GC, et al. (2007) The inheritance of resistance alleles in multiple sclerosis. PLoS Genet 3: 1607–1613.
- 11. Wu JS, James I, Wei Q, Castley A, Christiansen FT, et al. (2010) Influence of HLA-DRB1 allele heterogeneity on disease risk and clinical course in a West Australian MS cohort: a high-resolution genotyping study. Mult Scler 16: 526–532.
- 12. Chao MJ, Barnardo MC, Lui GZ, Lincoln MR, Ramagopalan SV, et al. (2007) Transmission of class I/II multi-locus MHC haplotypes and multiple sclerosis susceptibility: accounting for linkage disequilibrium. Hum Mol Genet 16: 1951–1958.
- 13. Lincoln MR, Montpetit A, Cader MZ, Saarela J, Dyment DA, et al. (2005) A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis. Nat Genet 37: 1108–1112.
- 14. Bergamaschi L, Ban M, Barizzone N, Leone M, Ferrante D, et al. (2011) Association of HLA class I markers with multiple sclerosis in the Italian and UK population: evidence of two independent protective effects. J Med Genet.
- 15. Bergamaschi L, Leone MA, Fasano ME, Guerini FR, Ferrante D, et al. (2010) HLA-class I markers and multiple sclerosis susceptibility in the Italian population. Genes Immun 11: 173–180.
- 16. IMSGC IMSGC (2007) A second major histocompatibility complex susceptibility locus for multiple sclerosis. Ann Neurol 61: 228–236.
- 17. Link J, Lorentzen AR, Kockum I, Duvefelt K, Lie BA, et al. (2010) Two HLA class I genes independently associated with multiple sclerosis. J Neuroimmunol 226: 172–176.
- 18. Rioux JD, Goyette P, Vyse TJ, Hammarstrom L, Fernando MM, et al. (2009) Mapping of multiple susceptibility variants within the MHC region for 7 immune-mediated diseases. Proc Natl Acad Sci U S A 106: 18680–18685.
- 19. De Jager PL, Jia X, Wang J, de Bakker PI, Ottoboni L, et al. (2009) Meta-analysis of genome scans and replication identify CD6, IRF8 and TNFRSF1A as new multiple sclerosis susceptibility loci. Nat Genet 41: 776–782.
- 20. Burfoot RK, Jensen CJ, Field J, Stankovich J, Varney MD, et al. (2008) SNP mapping and candidate gene sequencing in the class I region of the HLA complex: searching for multiple sclerosis susceptibility genes in Tasmanians. Tissue Antigens 71: 42–50.
- 21. D'Alfonso S, Bolognesi E, Guerini FR, Barizzone N, Bocca S, et al. (2008) A sequence variation in the MOG gene is involved in multiple sclerosis susceptibility in Italy. Genes Immun 9: 7–15.
- 22. Healy BC, Liguori M, Tran D, Chitnis T, Glanz B, et al. (2010) HLA B*44: protective effects in MS susceptibility and MRI outcome measures. Neurology 75: 634–640.
- 23. Cree BA, Rioux JD, McCauley JL, Gourraud PA, Goyette P, et al. (2010) A major histocompatibility Class I locus contributes to multiple sclerosis susceptibility independently from HLA-DRB1*15:01. PLoS ONE 5: e11296.
- 24. Nepom BS, Schwarz D, Palmer JP, Nepom GT (1987) Transcomplementation of HLA genes in IDDM. HLA-DQ alpha- and beta-chains produce hybrid molecules in DR3/4 heterozygotes. Diabetes 36: 114–117.
- 25. Sawcer S, Hellenthal G, Pirinen M, Spencer CC, Patsopoulos NA, et al. (2011) Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature 476: 214–219.
- 26. Mustafa M, Vingsbo C, Olsson T, Issazadeh S, Ljungdahl A, et al. (1994) Protective influences on experimental autoimmune encephalomyelitis by MHC class I and class II alleles. J Immunol 153: 3337–3344.
- 27. Friese MA, Jensen LT, Willcox N, Fugger L (2006) Humanized mouse models for organ-specific autoimmune diseases. Curr Opin Immunol 18: 704–709.
- 28. McMahon RM, Friis L, Siebold C, Friese MA, Fugger L, et al. (2011) Structure of HLA-A*0301 in complex with a peptide of proteolipid protein: insights into the role of HLA-A alleles in susceptibility to multiple sclerosis. Acta Crystallogr D Biol Crystallogr 67: 447–454.
- 29. McDonald WI, Compston A, Edan G, Goodkin D, Hartung HP, et al. (2001) Recommended diagnostic criteria for multiple sclerosis: guidelines from the International Panel on the diagnosis of multiple sclerosis. Ann Neurol 50: 121–127.
- 30. Poser CM, Paty DW, Scheinberg L, McDonald WI, Davis FA, et al. (1983) New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann Neurol 13: 227–231.
- 31. Olerup O, Zetterquist H (1992) HLA-DR typing by PCR amplification with sequence-specific primers (PCR-SSP) in 2 hours: an alternative to serological DR typing in clinical practice including donor-recipient matching in cadaveric transplantation. Tissue antigens 39: 225–235.
- 32. 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.
- 33. Hedstrom AK, Sundqvist E, Baarnhielm M, Nordin N, Hillert J, et al. (2011) Smoking and two human leukocyte antigen genes interact to increase the risk for multiple sclerosis. Brain 134: 653–664.
- 34. R_Development_Core_Team (2008) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing.
- 35. Dudbridge F (2008) Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data. Hum Hered 66: 87–98.
- 36. Hosmer DW, Lemeshow S (1992) Confidence interval estimation of interaction. Epidemiology 3: 452–456.
- 37.
Rothman K (2002) Epidemiology, an introduction. New York: Oxford University Press.