This study is an extension to previously published work that has linked variation in the human leukocyte antigen (HLA) class I region with susceptibility to multiple sclerosis (MS) in Australians from the Island State of Tasmania. Single nucleotide polymorphism (SNP) mapping was performed on an 865-kb candidate region (D6S1683-D6S265) in 166 Tasmanian MS families, and seven candidate genes [ubiquitin D (UBD), olfactory receptor 2H3 (OR2H3), gamma-aminobutyric acid B receptor 1 (GABBR1), myelin oligodendrocyte glycoprotein (MOG), HLA-F, HLA complex group 4 (HCG4) and HLA-G] were resequenced. SNPs tagging the extended MS susceptibility haplotype were genotyped in an independent sample of 356 Australian MS trios and SNPs in the MOG gene were significantly over-transmitted to MS cases. We identified significant effects on MS susceptibility of HLA-A*2 (OR: 0.51; P = 0.05) and A*3 (OR: 2.85; P = 0.005), and two coding polymorphisms in the MOG gene (V145I: P = 0.01, OR: 2.2; V142L: P = 0.04, OR: 0.45) after full conditioning on HLA-DRB1. We have therefore identified plausible candidates for the causal MS susceptibility allele, and although not conclusive at this stage, our data provide suggestive evidence for multiple class I MS susceptibility genes.

A recent genome-wide association study (GWAS) conducted by the International Multiple Sclerosis Genetics Consortium (IMSGC) identified a number of putative MS susceptibility genes. Here we have performed a replication study in 1134 Australian MS cases and 1265 controls for 17 risk-associated single nucleotide polymorphisms (SNPs) reported by the IMSGC. Of 16 SNPs that passed quality control filters, four, each corresponding to a different non-human leukocyte antigen (HLA) gene, were associated with disease susceptibility: KIAA0350 (rs6498169) P =0.001, IL2RA (rs2104286) P =0.033, RPL5 (rs6604026) P =0.041 and CD58 (rs12044852) P =0.042. There was no association ( P =0.58) between rs6897932 in the IL7R gene and the risk of MS. No interactions were detected between the replicated IMSGC SNPs and HLA-DRB1*15 , gender, disease course, disease progression or age-at-onset. We used a novel Bayesian approach to estimate the extent to which our data increased or decreased evidence for association with the six most-associated IMSGC loci. These analyses indicated that even modest P -values, such as those reported here, can contribute markedly to the posterior probability of ‘true’ association in replication studies. In conclusion, these data provide support for the involvement of four non-HLA genes in the pathogenesis of MS, and combined with previous data, increase to genome-wide significance ( P =3 × 10 −8 ) evidence of an association between KIAA0350 and risk of disease.

Association of multiple sclerosis (MS) with the human leukocyte antigen (HLA) class II haplotype DRB1*1501-DQB1*0602 is the most consistently replicated finding of genetic studies of the disease. However, the high level of linkage disequilibrium (LD) in the HLA region has hindered the identification of other loci that single-marker tests for association are unlikely to resolve. In order to address this issue, we generated haplotypes spanning 14.754 Mb (5 cM) across the entire HLA region. The haplotypes, which were inferred by genotyping relatives of 152 patients with MS and 105 unaffected control subjects of Tasmanian ancestry, define a genomic segment from D6S276 to D6S291, including 13 microsatellite markers integrated with allele-typing data for DRB1 and DQB1. Association to the DRB1*1501-DQB1*0602 haplotype was replicated. In addition, we found that the class I/extended class I region, defined by a genomic segment of approximately 400 kb between MOGCA and D6S265, harbors genes that independently increase risk of, or provide protection from, MS. Log-linear modeling analysis of constituent haplotypes that represent genomic regions containing class I (MOGCA-D6S265), class III (TNFa-TNFd-D6S273), and class II (DRB1-DQB1) genes indicated that having class I and class II susceptibility variants on the same haplotype provides an additive effect on risk. Moreover, we found no evidence for a disease locus in the class III region defined by a 150-kb genomic segment containing the TNF locus and 14 other genes. A global overview of LD performed using GOLD identified two discrete blocks of LD in the HLA region that correspond well with previous findings. We propose that the analysis of haplotypes, by use of the types of approaches outlined in the present article, should make it possible to more accurately define the contribution of the HLA to MS.