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

A Polymorphism in the HLA-DPB1 Gene Is Associated with Susceptibility to Multiple Sclerosis

  • Judith Field,

    Affiliations Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia, Centre for Neuroscience, University of Melbourne, Melbourne, Victoria, Australia

  • Sharon R. Browning,

    Affiliation Department of Statistics, The University of Auckland, Auckland, New Zealand

  • Laura J. Johnson,

    Affiliations Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia, Centre for Neuroscience, University of Melbourne, Melbourne, Victoria, Australia

  • Patrick Danoy,

    Affiliation The University of Queensland Diamantina Institute, Princess Alexandra Hospital, University of Queensland, Brisbane, Queensland, Australia

  • Michael D. Varney,

    Affiliation The Australian Red Cross Blood Service, Melbourne, Victoria, Australia

  • Brian D. Tait,

    Affiliation The Australian Red Cross Blood Service, Melbourne, Victoria, Australia

  • Kaushal S. Gandhi,

    Affiliation The Westmead Millenium Institute, University of Sydney, Sydney, New South Wales, Australia

  • Jac C. Charlesworth,

    Affiliation Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas, United States of America

  • Robert N. Heard,

    Affiliation The Westmead Millenium Institute, University of Sydney, Sydney, New South Wales, Australia

  • The Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) ,

    Members of the ANZgene Consortium are listed in the Acknowledgments.

  • Graeme J. Stewart,

    Affiliation The Westmead Millenium Institute, University of Sydney, Sydney, New South Wales, Australia

  • Trevor J. Kilpatrick,

    Affiliations Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia, Centre for Neuroscience, University of Melbourne, Melbourne, Victoria, Australia, Royal Melbourne Hospital, Melbourne, Victoria, Australia

  • Simon J. Foote,

    Affiliation Menzies Research Institute, University of Tasmania, Hobart, Tasmania, Australia

  • Melanie Bahlo,

    Affiliation The Walter and Eliza Hall Institute of Medical Research, Melbourne, Victoria, Australia

  • Helmut Butzkueven,

    Affiliations Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia, Department of Medicine, University of Melbourne, Melbourne, Victoria, Australia, Department of Neurology, Box Hill Hospital, Box Hill, Victoria, Australia

  • James Wiley,

    Affiliation Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia

  • David R. Booth,

    Affiliation The Westmead Millenium Institute, University of Sydney, Sydney, New South Wales, Australia

  • Bruce V. Taylor,

    Affiliation Menzies Research Institute, University of Tasmania, Hobart, Tasmania, Australia

  • Matthew A. Brown,

    Affiliations The University of Queensland Diamantina Institute, Princess Alexandra Hospital, University of Queensland, Brisbane, Queensland, Australia, Nuffield Department of Orthopaedic Surgery, Botnar Research Centre, University of Oxford, Oxford, United Kingdom

  • Justin P. Rubio,

    Affiliations Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia, GlaxoSmithKline, Harlow, Essex, United Kingdom

  •  [ ... ],
  • Jim Stankovich

    J.Stankovich@utas.edu.au

    Affiliation Menzies Research Institute, University of Tasmania, Hobart, Tasmania, Australia

  • [ view all ]
  • [ view less ]

Abstract

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each ). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls () and were highly significant in the combined dataset (). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set , replication set , combined ). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.

Introduction

Multiple sclerosis (MS, OMIM 126200) is an autoimmune disease of the central nervous system which is most prevalent in young adults of European ancestry. It is hypothesized that immune dysregulation leads to autoimmune attack on central nervous system myelin, and that both demyelination and axonal injury play key roles in disability progession. MS aggregates to some extent in families, and is clearly triggered by a complex mix of genetic and environmental factors. In Europeans, the HLA class II haplotype DRB1*15:01–DQB1*06:02 (DR15) is by far the strongest genetic risk factor (odds ratio approximately 3). There is evidence of interaction between this haplotype and environmental risk factors such as vitamin D [1] and exposure to infant siblings [2]. Recently genome-wide association (GWA) studies and follow-up studies have identified many other non-HLA loci associated with MS susceptibility [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14]. Most of these loci have been associated with other autoimmune diseases or encompass genes with known immune functions; the notable exceptions are associations with two genes involved in axonal transport [9], [13].

There is now also clear evidence of multiple, independent associations with MS across the major histocompatibility complex (MHC), as is the case for other autoimmune diseases [15], [16]. There have been associations with other alleles of HLA-DRB1 [17], [18], [19], [20], [21], independent associations in the class I region [4], [16], [22], [23], [24], [25], [26], [27], and cis [28] and trans [17], [29] interactions between alleles of class II genes. The outstanding candidate genes across the MHC are the highly polymorphic class II and class I genes which present antigens to CD4+ and CD8+ T cells respectively, but there are many other compelling candidates in the HLA complex, which is rich in genes of immune function. Elucidating causative genes and polymorphisms is harder than elsewhere in the genome due to the extensive linkage disequilibrium (LD) across the region.

To find evidence of independent associations, it is necessary either to perform stratified analyses or to condition on other associations. In a recent study, dense SNP genotyping data across the MHC was explored by conditioning on the predominant DRB1*15:01 association [16]. In the current study we extend this approach, conditioning on multiple MS-associated SNPs and imputed MHC alleles (in addition to DRB1*15:01) to identify further independent associations and build up multivariate models of MS-susceptibility. SNPs and alleles in these models were then validated by genotyping them (or proxies for them) in an independent replication case-control dataset.

Results

A multivariate SNP model of MS-susceptibility

Our discovery dataset consisted of the genotypes of 1618 MS cases and 3413 controls of European ancestry at 1927 SNPs across the MHC and flanking regions, samples which were previously analysed in a genome-wide association study [6]. We began by testing all SNPs individually for association with MS using a trend test. The most strongly associated SNP was rs9271366, a proxy for the HLA-DRB1*15:01 allele (genotypic squared correlation in HapMap CEU data [30], with DR15 in two-digit typing of 740 cases from the discovery dataset). In the second step, we tested all other SNPs using a trend test, conditioning on genotype at rs9271366 as a three-level factor. The most significantly associated SNP was rs2394160 (), located in an intron of LOC285830 approximately 15 kb centromeric of the HLA-F gene in the class I region. This SNP is in strong LD () with MS-associated SNPs rs2743951 and rs2523393 reported in previous studies [4], [16]. There were six other SNPs with P-values less than : five SNPs correlated with rs2394160 (rs2517912, rs1262126, rs2975033, rs6904029, rs4711207, see Figure 1) and rs2854050 in an intron of NOTCH4.

thumbnail
Figure 1. LD between the seven SNPs in the fitted model (Table 1) and other MHC polymorphisms.

The plot shows LD between the 7 model SNPs and other SNPs and classical HLA alleles genotyped in the HapMap CEU population. Points are only plotted for SNP pairs where . Points are plotted at the genomic positions of the HLA alleles/other SNPs; named points at the top of the plot () show the locations of the SNPs in the model. The various colours are used to distinguish LD with different SNPs in the model (rs2394160 pink, rs2854050 khaki, rs3830041 green, rs3129939 blue, rs9271366 purple, rs2187668 brown, rs9277535 red). SNPs are shown as small dots and HLA alleles are shown as triangles: refer to Table 2 for the names of HLA alleles correlated with . Pink stars denote SNPs in LD with rs2394160 that were strongly associated with MS in the second stage of the model-building (just with adjustment for the DRB1*15:01 proxy SNP rs9271366). The black triangles and crossess show the locations of some classical HLA genes and other genes respectively that are mentioned in the text. The red bars at the bottom of the graph show recombination rates averaged over 25kb windows. The vertical dashed lines mark the boundaries of the class I, class III and class II regions of the MHC.

https://doi.org/10.1371/journal.pone.0013454.g001

In the third step, all SNPs were tested conditioning on genotypes at rs9271366 and rs2394160. The strongest association was with rs2854050. In subsequent steps we identified four more independently-associated SNPs (rs3830041, rs3129939, rs2187668, rs9277535), including a proxy (rs2187668) for the HLA-DRB1*03:01 allele ( in HapMap CEU) which has been associated with MS previously [17], [18], [19], [20], [29]. No further SNPs were significant after applying a Bonferroni correction for 1927 tests. Table 1 shows a fitted model containing additive terms for these seven SNPs.

thumbnail
Table 1. Seven SNPs across the MHC independently associated with MS: terms in the fitted logistic regression model.

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

One might expect the form of this model to be heavily influenced by the terms included in the model-building procedure (the 1927 SNPs in the discovery dataset). However when we repeated the model-building procedure, augmenting the 1927 SNPs with 10,260 haplotype clusters selected using a variable-length Hidden Markov Model of LD [31], [32], the fitted model was remarkably similar (data not shown). This SNP-plus-haplotype model included 6 terms: 4 SNPs from the SNP-only model (rs2394160, rs2854050, rs9271366, rs9277535), and two 2-SNP haplotypes which each included a SNP from the SNP-only model (rs2854050, rs2187668). Hence we focused on the SNP-only model for replication and interpretation.

We attempted to replicate the 7-SNP model (Table 1) by genotyping these SNPs, and others across the MHC, in an independent dataset comprising 2212 cases and 2251 controls. In the replication data all terms were significant after correcting for multiple testing (one-sided , Table 1); hence they were all highly significant in the combined dataset (). In the combined dataset we tested for departures from multiplicativity (dominance effects) and interactions with the DRB1*15:01 proxy rs9271366. There were significant dominance effects for the DRB1*15:01 and DRB1*03:01 proxy SNPs, and there was modest evidence of interaction () between rs2854050 and rs9271366 (Table S1).

Correlations with classical MHC alleles

The seven SNPs in the SNP-only model (Table 1) are in LD with many other polymorphisms across the MHC. Figure 1 shows LD () between the SNPs in the model and other SNPs, as well as alleles of the classical MHC genes HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1 and HLA-DQB1, in 90 unrelated individuals from the HapMap CEU population [30]. Table 2 shows classical alleles correlated with .

thumbnail
Table 2. Alleles of the classical MHC genes HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DQA1 and HLA-DQB1 correlated with SNPs in the fitted logistic regression model ().

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

We investigated the effects of some of these HLA alleles, by imputing them in the discovery dataset and genotyping proxy SNPs in the replication dataset. We added imputed allele dosages and proxy SNPs to the logistic regression model one-by-one, both with and without the corresponding correlated SNP in the original model (see Table S2 for imputed dosages and Table S3 for proxy SNPs).

These analyses provided evidence that associations with the three most telomeric model SNPs (rs2394160, rs2854050 and rs3830041) may be driven by LD with MHC alleles that have been associated with MS previously: HLA-A*02:01 (protective) [22], [23], [33], DRB1*04:01 (protective) [20], [34] and DRB1*13:03 (susceptibility) [21], [35] respectively. In the discovery dataset, imputed dosages for these three HLA alleles show strong associations when they are substituted for the relevant SNP (HLA-A*02:01, ; DRB1*04:01, ; DRB1*13:03, , Table S2). In the discovery dataset rs6904029, which was strongly associated in the second step of the model building (pink stars in Figure 1), is strongly correlated with HLA-A*02:01 (). In the replication dataset the association with rs2394160 completely disappears () after adjustment for a good proxy for HLA-A*02:01 (rs2523822, , Table S3). This SNP is highly significant in the model, both with and without adjustment for rs2394160 ( and respectively). After adjustment for rs2394160 in the replication dataset, there is little evidence of association with candidate polymorphisms rs2857766 and rs1233334 in the nearby genes MOG [22], [36] and HLA-G [37] (Table S3).

These analyses (Table S2) did not provide any evidence that the association with rs3129939 is caused by LD with correlated alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1 or HLA-DQB1 (Table 2), and rs9277535 shows low correlation with all alleles of these genes. Table 3 shows an alternative logistic regression model fitted to the discovery dataset, with the other five SNPs replaced by imputed dosages for correlated HLA alleles.

thumbnail
Table 3. Alternative logistic regression model fitted to the discovery dataset, with five SNPs from the model in Table 1 replaced by imputed dosages for correlated MHC alleles.

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

Association with a SNP in the gene HLA-DPB1

rs9277535 lies in the 3¿ untranslated region of another class II gene, HLA-DPB1. All moderately correlated SNPs (, Figure 1) lie within a 131-kb region (33,155kb–33,286kb) which also encompasses the genes HLA-DPB2, COL11A2, RXRB, SLC39A7, HSD17B8, MIR219-1 and RING1. HLA-DPB1 is centromeric of the class II genes HLA-DRB1, HLA-DQA1 and HLA-DQB1, and is separated by a region with a high recombination rate (Figure 1). Consequently, LD is low between rs9277535 and all DRB1, DQA1 and DQB1 alleles in the HapMap CEU dataset (). In the discovery dataset, the association with rs9277535 remains strong after conditioning on imputed allele dosages at all alleles of DRB1, DQA1 and DQB1 (), and with full conditioning on maximum likelihood genotypes at DRB1 (); hence the association is unlikely to be caused by complex interactions between DRB1 alleles [17], [29]. There was no evidence that rs9277535G exerts different effects in DR15-positive and -negative individuals ( testing for interaction with rs9271366 in the combined dataset).

The association with rs9277535 is somewhat masked due to negative confounding by DRB1*15:01, even though LD between the DRB1 and DPB1 loci is weak. Table S4 shows correlations between the risk-associated rs9277535 allele and imputed allele dosages for various classical alleles in the discovery dataset. In particular, rs9277535G is negatively-correlated with the DRB1*15:01 allele (r = ¿0.09). In an unadjusted analysis of this SNP in the discovery dataset, the estimated odds ratio for rs9277535G is 1.19 (). The estimated odds ratio increases to 1.34 with adjustment for DRB1*15:01 and to 1.36 in the full model (Table 3) after further adjustment for DRB1*03:01 (negative confounder) and DRB1*04:01 (positive confounder). In the replication dataset, negative confounding means that the association with rs9277535 is significant after adjustment for other terms, even though the overall frequency of the minor allele is very similar in cases and controls (0.265 and 0.261 respectively).

A recent study showed that there is a highly-conserved vitamin D response element (VDRE) GGGTGGAGGGGGTTCA immediately upstream of the HLA-DRB1 gene on DRB1*15:01 haplotypes [1]. We searched for Vitamin D Response Elements (VDREs) associated with the HLA-DPB1*03:01 allele [38], using the homozygous COX cell line which carries this allele. This in silico analysis revealed one predicted VDRE. However this sequence differed from the sequence proximal to HLA-DRB1*15:01, was not located in the proximal promoter region of HLA-DPB1*03:01 and was found in all HLA-DPB1 alleles assessed.

Association with a copy number variant

We observed a copy number variant in the course of genotyping the HLA-B*44:02 proxy SNP rs2256583CT in the replication dataset. This SNP lies between HLA-C and HLA-B in a region known to contain copy number variants. Plotting areas under the C- and T-allele peaks revealed three groups of heterozygotes, which appear to correspond to C∶T ratios of 2∶1, 1∶1 and 1∶2 (Figure S1). There were different proportions of cases in the three groups (Table S5, using Fisher's exact test). The heterozygote groupings are not well-correlated with any of the seven SNPs in the model of Table 1, and these differences remain almost significant after adjustment for all terms in the SNP model ( using a likelihood ratio test with two degrees of freedom). There was too much inter-sample variation in peak heights to detect copy number differences among homozygotes.

Discussion

Using multivariate analysis we identified seven SNPs that were independently associated with MS, and confirmed all seven associations in a replication dataset. All associations were highly significant in the combined dataset (, Table 1).

Two of these SNPs are proxies for the known associations with DRB1*15:01 and DRB1*03:01. As in previous analyses [4], [16], the most prominent associations after conditioning on DRB1*15:01 were with SNPs in the vicinity of HLA-A. Hence our data add to previous, very strong evidence of a functional effect correlated with rs2394160, rs2743951 [16], rs2523393 [4], rs6904029 and HLA-A*02:01 [33]. We also found evidence to support previously-reported associations with the alleles DRB1*04:01 [20], [34] and DRB1*13:03 [21], [35].

In MS, strong protective associations have been reported previously for the class I alleles HLA-C*05:01 [26] and HLA-B*44:02 [16], which are highly correlated with one another ( in HapMap CEU data). In an unadjusted analysis, HLA-B*44:02 is strongly associated with MS in our discovery dataset (OR 0.67, ), however this association is positively confounded by the DRB1*15:01 association (r = ¿0.05), and by the (putative) protective effects of HLA-A*02:01 () and DRB1*04:01 (). Hence, HLA-B*44:02 shows little evidence of association when added to the model in Table 3 (OR = 0.92, ). Our data provide stronger evidence of an association with DRB1*04:01 than with HLA-C*05:01 or HLA-B*44:02, but this evidence is far from definitive. The multivariate analyses suggest that not all three of these alleles are exerting causative effects, but further analyses in larger datasets and functional studies will be required to determine which (if any) are causative.

The most novel association in the model is with the SNP rs9277535 in the class II gene HLA-DPB1, centromeric of the DR/DQ region and separated by a recombination hotspot. The DP region has not been investigated as thoroughly as the DR/DQ region in genetic association studies. In normal conditions DP gene products are expressed at much lower levels than DR gene products [39]; therefore DP typing was difficult with serological techniques prior to the advent of molecular typing. Associations with DR and DQ alleles were detected early on for many diseases including MS [40], so there has been a natural focus on these genes.

Nevertheless there have been several studies reporting associations between MS and HLA-DPB1. A SNP rs3135021 in an intron of HLA-DPB1 was recently associated with MS in African Americans [41]. The SNP is in LD with some HLA-DRB1 alleles in African Americans [41], and it is in linkage equilibrium () with rs9275535 in the HapMap CEU population. The DPB1*05:01 allele has been associated with opticospinal MS in Japan [42] and with MS in Southern Han Chinese [43]. The DPB1*03:01 allele has been implicated in epitope spreading in MS [44]. It has also been associated with opticospinal MS in Japan [45], and with MS in a small Australian study [46] and in Sardinia [35]. In the Sardinian study, comparing 835 relapsing MS patients with 592 controls, carriage of the DPB1*03:01 allele was associated with increased risk of MS (OR 1.30, ) after conditioning on an associated microsatellite in the class I region and on presence or absence of five DRB1–DQB1 risk haplotypes (three of which include DRB1*15:01, DRB1*03:01 and DRB1*13:03). We tested for presence/absence of the DPB1*03:01 allele in 422 Australian MS cases using a combination of sequence specific priming and sequence based typing, and found a correlation with rs9277535 (, Table S6): all but two individuals carrying the DPB1*03:01 allele also carry the rs9277535G allele.

Mechanisms other than peptide affinities for HLA molecules may be underpinning some of these associations. The HLA complex has a particularly high density of copy number variants, and these represent one possible source of functional variation. In the course of replication genotyping we identified a copy number variant between the genes HLA-B and HLA-C which showed some evidence of independent association with MS susceptibility.

Disease-associated variants may also be acting via effects on gene expression [47]. A recent study suggests that a vitamin D response element (VDRE) upstream of HLA-DRB1 plays an important role in expression of this gene [1]. VDREs are bound by the vitamin D receptor (VDR) which acts as a ligand activated transcription factor, and the VDRE upstream of HLA-DRB1 is highly conserved on DRB1*15:01 haplotypes. It was found that the DRB1*15:01 haplotype binds the VDR more efficiently than other haplotypes, in particular the 98% of DRB1*04, DRB1*07 and DRB1*09 haplotypes where this sequence differs at two nucleotides (GGGTGGAGAGGGGTCA) [1]. We observed that the homozygous COX and QBL lymphoblastoid cell lines, which both carry the risk-associated HLA-DRB1*03:01 allele, have an identical proximal VDRE to that occurring on the DRB1*15:01 haplotype. This is noteworthy given that the amino acid sequences of DRB1*15:01 and DRB1*03:01 are not particularly similar compared to other DRB1 alleles.

We also searched for VDREs associated with the HLA-DPB1*0301 allele, but did not find convincing evidence of an active response element. One of the genes in the region of LD around HLA-DPB1 is retinoid X receptor beta (RXRB). This gene encodes a receptor which forms a dimer with VDR and increases its DNA binding and transcriptional function [48].

In summary, we have detected a highly significant association with MS in the vicinity of the HLA-DPB1 gene that is independent of other MS associations in the HLA complex, confirming suggestive reports from previous smaller studies. Our data indicate that DPB1*03:01 may be the causative HLA allele underpinning this association, however positive replication studies followed by fine mapping will be required to confirm and resolve this association in other populations. Our study also highlights the utility of multivariate analysis to disentangle HLA associations for autoimmune diseases, such as MS, which have multiple genetic risk factors in the HLA complex. This approach will be useful more generally at other loci and for other diseases, as more clusters of associations are identified [14], [49].

Materials and Methods

Approval for this research was granted by the Melbourne Health Human Research Ethics Committee and other institutional ethics committees. Written consent was given by the patients for their information to be stored in the study database and used for research.

To identify SNPs across the MHC independently associated with MS, we analysed genotype data assembled for a GWA study [6]. After data cleaning followed by principal components analysis to remove population outliers [6], this dataset comprised 1618 MS cases of European ancestry from Australia and New Zealand, and 3413 European-ancestry controls from Australia, the UK and US. All samples were genotyped with Illumina arrays.

The cleaned discovery dataset included genotypes at 1927 SNPs across the MHC and flanking regions (chr6: 24–36 Mb, NCBI dbSNP genome build 128). We used a stepwise logistic regression procedure to identify SNPs independently associated with disease status, building a model with increasing numbers of independently associated SNPs. At each step we used a trend test to identify the most significantly associated SNP after conditioning on SNPs already in the model, then added this SNP to the model. Conditioning on SNPs already in the model was performed by coding their genotypes as three-level factors. The procedure was repeated until no further SNPs were significant after applying a Bonferroni correction to correct for 1927 tests. In the final model fitted to the discovery dataset (Table 1), all SNPs in the model were coded as quantitative variables (trend test) rather than factors. The final model was also fitted to the discovery dataset with adjustment for position along the first three eigenvectors of the final principal components analysis [6]; with this additional adjustment no coefficients of SNPs changed by more than 7% and all remained significant at (full data not shown).

Twenty-five SNPs across the HLA region were genotyped in an independent replication dataset of MS cases and controls from Australia and New Zealand described previously [6], using the Sequenom MassARRAY system and iPlex Gold chemistry. These SNPs included the seven SNPs in the fitted model (Table 1), candidate polymorphisms in the genes MOG and HLA-G, proxies for several classical MHC alleles, and backups in case of assay failure. After excluding 147 samples with call rates less than 90%, the filtered replication dataset comprised 2212 cases and 2251 controls. In this filtered dataset call rates were above 99.7% for all 25 SNPs.

We fitted the same 7-SNP logistic regression model to the replication dataset and to the combined dataset. Two-sided values were used throughout except in the analysis of the replication dataset. Analysis of the combined dataset was adjusted for sample group (discovery/replication). In the replication dataset we also fitted a 7-SNP model using an alternative proxy rs3135388 for HLA-DRB1*15:01 [30]; in this alternative model the magnitudes of the coefficients of rs2854050 and rs3830041 decreased by 32% and 13% respectively but remained significant (), and no other coefficients changed by more than 6% (full data not shown).

To investigate LD with polymorphisms not genotyped in the discovery dataset we downloaded phased data for individuals from the HapMap CEU population (Utah residents with ancestry from northern and western Europe), including classical HLA typing data for the genes HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1 and HLA-DQB1 (http://www.inflammgen.org) [30]. In the discovery dataset we imputed genotypes at classical HLA alleles using the software BEAGLE (http://www.stat.auckland.ac.nz/bbrowning/beagle/beagle.html) [50]. Imputation was done by dividing the samples randomly into batches of 50, and running each batch with the phased HapMap CEU data. The default settings in BEAGLE were changed to allow consideration of long haplotypes spanning the entire region, which has proven beneficial for imputation of HLA alleles [51]. For each classical allele we calculated posterior allele dosage, which is twice the posterior probability of carrying two copies of the allele plus the posterior probability of carrying exactly one copy of the allele. These posterior allele dosages were used as predictive variables in logistic regression models (Table 3, Table S2).

MHC sequences for HLA-homozygous cell lines were downloaded from the website of the MHC Haplotype project (https://www.sanger.ac.uk/HGP/Chr6/MHC/). We searched for VDREs associated with the HLA-DPB1*03:01 allele using the sequence for this allele (carried by the homozygous cell line COX). The sequence was scanned for VDREs using JASPAR (http://jaspar.genereg.net) [38] with a profile threshold of 80%.

Supporting Information

Table S1.

The 7-SNP model in the combined dataset, with two significant dominance effects and one significant SNP-SNP interaction added

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

(0.02 MB PDF)

Table S2.

The effects of adding correlated MHC alleles to the fitted model in the discovery dataset (Table 1), both with and without the correlated SNP. For these analyses imputed MHC allele dosages were calculated using the software BEAGLE (50).

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

(0.03 MB PDF)

Table S3.

The effects of adding other SNPs genotyped in the replication dataset to the fitted model (Table 1), both with and without the correlated SNP from the model.

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

(0.03 MB PDF)

Table S4.

Correlations between various SNPs and imputed classical alleles in the discovery dataset

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

(0.01 MB PDF)

Table S5.

A copy number variant at SNP rs2256583: numbers of heterozygous cases and controls with various C∶T allele ratios in the replication dataset

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

(0.01 MB PDF)

Table S6.

Correlation between HLA-DPB1*0301 and rs9277535 in 422 MS cases

https://doi.org/10.1371/journal.pone.0013454.s006

(0.01 MB PDF)

Figure S1.

A copy number variant at SNP rs2256583. A plot of area under the C- allele peak (y-axis) versus area under the T-allele peak (x-axis) for heterozygotes, showing three groups of individuals with C∶T ratios of 1∶2 (red), 1∶1 (green) and 2∶1 (blue).

https://doi.org/10.1371/journal.pone.0013454.s007

(0.16 MB PDF)

Acknowledgments

We thank individuals with MS in Australia and New Zealand for supporting this research. We are grateful to J. Wright and C. Remediakis from Multiple Sclerosis Research Australia for expediting this research. Replication genotyping was conducted at the Murdoch Children's Research Institute Sequenom Platform Facility. Healthy control DNAs were provided by the Australian Red Cross Blood Service. We acknowledge use of genotype data from the British 1958 Birth Cohort DNA collection.

The ANZgene Consortium members are: Melanie Bahlo, 1 David R Booth, 2 Simon A Broadley, 3,4 Matthew A Brown, 5,6 Brian L Browning, 7 Sharon R Browning, 7 Helmut Butzkueven, 8–10 William M Carroll, 11–13 Patrick Danoy, 5 Judith Field, 9,14 Simon J Foote, 15 Lyn R Griffiths, 16 Robert N Heard, 2 Allan G Kermode, 11–13 Trevor J Kilpatrick, 9,14,17 Jeanette Lechner-Scott, 18,19 Pablo Moscato, 19,20 Victoria M Perreau, 14 Rodney J Scott, 18–20 Mark Slee, 21 Jim Stankovich, 15 Graeme J Stewart, 2 Bruce V Taylor, 15 James Wiley9 (Chair).

(1) The Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria Australia.

(2) The Westmead Millenuium Institute, Westmead, New South Wales, Australia.

(3) School of Medicine, Griffith University, Queensland, Australia.

(4) Department of Neurology, Gold Coast Hospital, Queensland, Australia.

(5) The University of Queensland Diamantina Institute, Princess Alexandra Hospital, University of Queensland, Brisbane, Queensland, Australia.

(6) Botnar Research Centre, Nuffield Department of Orthopaedic Surgery, University of Oxford, Oxford, UK.

(7) Department of Statistics, The University of Auckland, Auckland, New Zealand.

(8) Department of Medicine, University of Melbourne, Melbourne, Victoria, Australia.

(9) Florey Neuroscience Institutes, University of Melbourne, Melbourne, Victoria, Australia.

(10) Department of Neurology, Box Hill Hospital, Victoria, Australia.

(11) Sir Charles Gairdner Hospital, Nedlands, Western Australia, Australia.

(12) Australian Neuromuscular Research Institute, Nedlands, Western Australia, Australia.

(13) Centre for Neuromuscular and Neurological Disorders, University of Western Australia, Nedlands, Western Australia, Australia.

(14) Centre for Neuroscience, University of Melbourne, Victoria, Australia.

(15) Menzies Research Institute, University of Tasmania, Hobart, Tasmania, Australia.

(16) Genomics Research Centre, Griffith University, Queensland, Australia.

(17) Royal Melbourne Hospital, Parkville, Victoria, Australia.

(18) John Hunter Hospital, Hunter New England Health Service, Newcastle, New South Wales, Australia.

(19) Hunter Medical Research Institute, Newcastle, New South Wales, Australia.

(20) Centre for Bioinformatics, Biomarker Discovery and Information-based Medicine, University of Newcastle, New South Wales, Australia.

(21) School of Medicine, Department of Neurology, Flinders University, Bedford Park, Adelaide, South Australia, Australia.

Author Contributions

Conceived and designed the experiments: JF SRB JW JR JS. Performed the experiments: JF LJJ PD MDV BDT JR. Analyzed the data: JF JS. Contributed reagents/materials/analysis tools: SRB MDV BDT KG JCC RNH GS TJK SJF MB HB JW DB BT MAB JR. Wrote the paper: JF SRB HB BT JR JS.

References

  1. 1. Ramagopalan SV, Maugeri NJ, Handunnetthi L, Lincoln MR, Orton SM, et al. (2009) Expression of the multiple sclerosis-associated MHC class II Allele HLA-DRB1*1501 is regulated by vitamin D. PLoS Genet 5: e1000369.
  2. 2. van der Mei IA, Ponsonby AL, Taylor BV, Stankovich J, Dickinson JL, et al. (2009) Human leukocyte antigen-DR15, low infant sibling exposure and multiple sclerosis: Gene-environment interaction. Ann Neurol 67: 261–265.
  3. 3. Hafler JP, Maier LM, Cooper JD, Plagnol V, Hinks A, et al. (2009) CD226 Gly307Ser association with multiple autoimmune diseases. Genes Immun 10: 5–10.
  4. 4. 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.
  5. 5. Ban M, Goris A, Lorentzen AR, Baker A, Mihalova T, et al. (2009) Replication analysis identifies TYK2 as a multiple sclerosis susceptibility factor. Eur J Hum Genet 17: 1309–1313.
  6. 6. Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) (2009) Genome-wide association study identifies new multiple sclerosis susceptibility loci on chromosomes 12 and 20. Nat Genet 41: 824–828.
  7. 7. International Multiple Sclerosis Genetics Consortium (IMSGC) (2008) Refining genetic associations in multiple sclerosis. Lancet Neurol 7: 567–569.
  8. 8. Hoppenbrouwers IA, Aulchenko YS, Ebers GC, Ramagopalan SV, Oostra BA, et al. (2008) EVI5 is a risk gene for multiple sclerosis. Genes Immun 9: 334–337.
  9. 9. Aulchenko YS, Hoppenbrouwers IA, Ramagopalan SV, Broer L, Jafari N, et al. (2008) Genetic variation in the KIF1B locus influences susceptibility to multiple sclerosis. Nat Genet 40: 1402–1403.
  10. 10. Hafler DA, Compston A, Sawcer S, Lander ES, Daly MJ, et al. (2007) Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med 357: 851–862.
  11. 11. Rubio JP, Stankovich J, Field J, Tubridy N, Marriott M, et al. (2008) Replication of KIAA0350, IL2RA, RPL5 and CD58 as multiple sclerosis susceptibility genes in Australians. Genes Immun 9: 624–630.
  12. 12. Jakkula E, Leppa V, Sulonen AM, Varilo T, Kallio S, et al. (2010) Genome-wide association study in a high-risk isolate for multiple sclerosis reveals associated variants in STAT3 gene. Am J Hum Genet 86: 285–291.
  13. 13. International Multiple Sclerosis Genetics Consortium (IMSGC) (2010) Comprehensive follow-up of the first genome-wide association study of multiple sclerosis identifies KIF21B and TMEM39A as susceptibility loci. Hum Mol Genet 19: 953–962.
  14. 14. 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.
  15. 15. Nejentsev S, Howson JM, Walker NM, Szeszko J, Field SF, et al. (2007) Localization of type 1 diabetes susceptibility to the MHC class I genes HLA-B and HLA-A. Nature 450: 887–892.
  16. 16. 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.
  17. 17. Barcellos LF, Sawcer S, Ramsay PP, Baranzini SE, Thomson G, et al. (2006) Heterogeneity at the HLA-DRB1 locus and risk for multiple sclerosis. Hum Mol Genet 15: 2813–2824.
  18. 18. Fernando MM, Stevens CR, Walsh EC, De Jager PL, Goyette P, et al. (2008) Defining the role of the MHC in autoimmunity: a review and pooled analysis. PLoS Genet 4: e1000024.
  19. 19. Marrosu MG, Murru MR, Costa G, Murru R, Muntoni F, et al. (1998) DRB1-DQA1-DQB1 loci and multiple sclerosis predisposition in the Sardinian population. Hum Mol Genet 7: 1235–1237.
  20. 20. Stankovich J, Butzkueven H, Marriott M, Chapman C, Tubridy N, et al. (2009) HLA-DRB1 associations with disease susceptibility and clinical course in Australians with multiple sclerosis. Tissue Antigens 74: 17–21.
  21. 21. Wu JS, James I, Qiu W, 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.
  22. 22. 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.
  23. 23. 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.
  24. 24. 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.
  25. 25. Marrosu MG, Murru R, Murru MR, Costa G, Zavattari P, et al. (2001) Dissection of the HLA association with multiple sclerosis in the founder isolated population of Sardinia. Hum Mol Genet 10: 2907–2916.
  26. 26. Yeo TW, De Jager PL, Gregory SG, Barcellos LF, Walton A, et al. (2007) A second major histocompatibility complex susceptibility locus for multiple sclerosis. Ann Neurol 61: 228–236.
  27. 27. Cree BAC, Rioux JD, McCauley JL, Gourraud P-AFD, 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.
  28. 28. Lincoln MR, Ramagopalan SV, Chao MJ, Herrera BM, Deluca GC, et al. (2009) Epistasis among HLA-DRB1, HLA-DQA1, and HLA-DQB1 loci determines multiple sclerosis susceptibility. Proc Natl Acad Sci U S A 106: 7542–7547.
  29. 29. 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.
  30. 30. de Bakker PI, McVean G, Sabeti PC, Miretti MM, Green T, et al. (2006) A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nat Genet 38: 1166–1172.
  31. 31. Browning SR (2006) Multilocus association mapping using variable-length Markov chains. Am J Hum Genet 78: 903–913.
  32. 32. Browning BL, Browning SR (2007) Efficient multilocus association testing for whole genome association studies using localized haplotype clustering. Genet Epidemiol 31: 365–375.
  33. 33. Bergamaschi L, Leone MA, Fasano ME, Guerini FR, Ferrante D, et al. HLA-class I markers and multiple sclerosis susceptibility in the Italian population. Genes Immun 11: 173–180.
  34. 34. McDonnell GV, Mawhinney H, Graham CA, Hawkins SA, Middleton D (1999) A study of the HLA-DR region in clinical subgroups of multiple sclerosis and its influence on prognosis. J Neurol Sci 165: 77–83.
  35. 35. Marrosu MG, Cocco E, Costa G, Murru MR, Mancosu C, et al. (2006) Interaction of loci within the HLA region influences multiple sclerosis course in the Sardinian population. J Neurol 253: 208–213.
  36. 36. 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.
  37. 37. Wiendl H, Feger U, Mittelbronn M, Jack C, Schreiner B, et al. (2005) Expression of the immune-tolerogenic major histocompatibility molecule HLA-G in multiple sclerosis: implications for CNS immunity. Brain 128: 2689–2704.
  38. 38. Portales-Casamar E, Thongjuea S, Kwon AT, Arenillas D, Zhao X, et al. (2010) JASPAR 2010: the greatly expanded open-access database of transcription factor binding profiles. Nucleic Acids Res 38: D105–110.
  39. 39. Pawelec G, Buhring HJ (1990) Expression of MHC class II epitopes on human T lymphocyte clones. Cell Immunol 127: 520–526.
  40. 40. Tiwari JL, Terasaki PI (1981) HLA-DR and disease associations. Prog Clin Biol Res 58: 151–163.
  41. 41. McElroy JP, Cree BAC, Caillier SJ, Gregersen PK, Herbert J, et al. (2010) Refining the association of MHC with multiple sclerosis in African Americans. Hum Mol Genet 19: 3080–3088.
  42. 42. Yamasaki K, Horiuchi I, Minohara M, Kawano Y, Ohyagi Y, et al. (1999) HLA-DPB1*0501-associated opticospinal multiple sclerosis: clinical, neuroimaging and immunogenetic studies. Brain 122(Pt 9): 1689–1696.
  43. 43. Wu XM, Wang C, Zhang KN, Lin AY, Kira J, et al. (2009) Association of susceptibility to multiple sclerosis in Southern Han Chinese with HLA-DRB1, -DPB1 alleles and DRB1-DPB1 haplotypes: distinct from other populations. Mult Scler 15: 1422–1430.
  44. 44. Yu M, Kinkel RP, Weinstock-Guttman B, Cook DJ, Tuohy VK (1998) HLA-DP: a class II restriction molecule involved in epitope spreading during the development of multiple sclerosis. Hum Immunol 59: 15–24.
  45. 45. Fukazawa T, Kikuchi S, Miyagishi R, Miyazaki Y, Yabe I, et al. (2006) HLA-dPB1*0501 is not uniquely associated with opticospinal multiple sclerosis in Japanese patients. Important role of DPB1*0301. Mult Scler 12: 19–23.
  46. 46. Dekker JW, Easteal S, Jakobsen IB, Gao X, Stewart GJ, et al. (1993) HLA-DPB1 alleles correlate with risk for multiple sclerosis in Caucasoid and Cantonese patients lacking the high-risk DQB1*0602 allele. Tissue Antigens 41: 31–36.
  47. 47. Handel AE, Handunnetthi L, Berlanga AJ, Watson CT, Morahan JM, et al. (2010) The effect of single nucleotide polymorphisms from genome wide association studies in multiple sclerosis on gene expression. PLoS ONE 5: e10142.
  48. 48. Yu VC, Delsert C, Andersen B, Holloway JM, Devary OV, et al. (1991) RXR beta: a coregulator that enhances binding of retinoic acid, thyroid hormone, and vitamin D receptors to their cognate response elements. Cell 67: 1251–1266.
  49. 49. Al Olama AA, Kote-Jarai Z, Giles GG, Guy M, Morrison J, et al. (2009) Multiple loci on 8q24 associated with prostate cancer susceptibility. Nat Genet 41: 1058–1060.
  50. 50. Browning BL, Browning SR (2009) A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet 84: 210–223.
  51. 51. Leslie S, Donnelly P, McVean G (2008) A statistical method for predicting classical HLA alleles from SNP data. Am J Hum Genet 82: 48–56.