[Skip to Navigation]
Sign In
Figure. 
Carriership of HLA-DRB1*1501 and focal spinal cord lesions. Scatterplot of focal spinal cord lesions for carriers and noncarriers of HLA-DRB1*1501 (measured as the presence of the A allele of rs3135388). Carriers of HLA-DRB1*1501 had a median of 4 spinal cord lesions compared with 2 lesions for noncarriers (P < .001, Mann-Whitney test). Line reflects the median number of focal abnormalities.

Carriership of HLA-DRB1*1501 and focal spinal cord lesions. Scatterplot of focal spinal cord lesions for carriers and noncarriers of HLA-DRB1*1501 (measured as the presence of the A allele of rs3135388). Carriers of HLA-DRB1*1501 had a median of 4 spinal cord lesions compared with 2 lesions for noncarriers (P < .001, Mann-Whitney test). Line reflects the median number of focal abnormalities.

Table 1. 
Patient Demographic, Clinical, and Magnetic Resonance (MR) Imaging Characteristics According to Multiple Sclerosis Subtype
Patient Demographic, Clinical, and Magnetic Resonance (MR) Imaging Characteristics According to Multiple Sclerosis Subtype
Table 2. 
Correlation by Kruskal-Wallis Test of Single-Nucleotide Polymorphisms (SNPs) With Magnetic Resonance Imaging Variables
Correlation by Kruskal-Wallis Test of Single-Nucleotide Polymorphisms (SNPs) With Magnetic Resonance Imaging Variables
1.
Barcellos  LFOksenberg  JRBegovich  AB  et al. Multiple Sclerosis Genetics Group, HLA-DR2 dose effect on susceptibility to multiple sclerosis and influence on disease course.  Am J Hum Genet 2003;72 (3) 710- 716PubMedGoogle ScholarCrossref
2.
Hoffmann  SCepok  SGrummel  V  et al.  HLA-DRB1*0401 and HLA-DRB1*0408 are strongly associated with the development of antibodies against interferon-beta therapy in multiple sclerosis [published correction appears in Am J Hum Genet. 2008;83(4):541].  Am J Hum Genet 2008;83 (2) 219- 227PubMedGoogle ScholarCrossref
3.
Yamasaki  KHoriuchi  IMinohara  M  et al.  HLA-DPB1*0501-associated opticospinal multiple sclerosis: clinical, neuroimaging and immunogenetic studies.  Brain 1999;122 (pt 9) 1689- 1696PubMedGoogle ScholarCrossref
4.
DeLuca  GCRamagopalan  SVHerrera  BM  et al.  An extremes of outcome strategy provides evidence that multiple sclerosis severity is determined by alleles at the HLA-DRB1 locus.  Proc Natl Acad Sci U S A 2007;104 (52) 20896- 20901PubMedGoogle ScholarCrossref
5.
Kantarci  OHde Andrade  MWeinshenker  BG Identifying disease modifying genes in multiple sclerosis.  J Neuroimmunol 2002;123 (1-2) 144- 159PubMedGoogle ScholarCrossref
6.
Barcellos  LFSawcer  SRamsay  PP  et al.  Heterogeneity at the HLA-DRB1 locus and risk for multiple sclerosis.  Hum Mol Genet 2006;15 (18) 2813- 2824PubMedGoogle ScholarCrossref
7.
Fazekas  FEnzinger  CRopele  SSchmidt  HSchmidt  RStrasser-Fuchs  S The impact of our genes: consequences of the apolipoprotein E polymorphism in Alzheimer disease and multiple sclerosis.  J Neurol Sci 2006;245 (1-2) 35- 39PubMedGoogle ScholarCrossref
8.
van Veen  Tvan Winsen  LCrusius  JB  et al.  αB-crystallin genotype has impact on the multiple sclerosis phenotype.  Neurology 2003;61 (9) 1245- 1249PubMedGoogle ScholarCrossref
9.
Lucchinetti  CBrück  WParisi  JScheithauer  BRodriguez  MLassmann  H Heterogeneity of multiple sclerosis lesions: implications for the pathogenesis of demyelination.  Ann Neurol 2000;47 (6) 707- 717PubMedGoogle ScholarCrossref
10.
Lee  MASmith  SPalace  J  et al.  Spatial mapping of T2 and gadolinium-enhancing T1 lesion volumes in multiple sclerosis: evidence for distinct mechanisms of lesion genesis?  Brain 1999;122 (pt 7) 1261- 1270PubMedGoogle ScholarCrossref
11.
Nijeholt  GJBergers  EKamphorst  W  et al.  Post-mortem high-resolution MRI of the spinal cord in multiple sclerosis: a correlative study with conventional MRI, histopathology and clinical phenotype.  Brain 2001;124 (pt 1) 154- 166PubMedGoogle ScholarCrossref
12.
Bakshi  RThompson  AJRocca  MA  et al.  MRI in multiple sclerosis: current status and future prospects.  Lancet Neurol 2008;7 (7) 615- 625PubMedGoogle ScholarCrossref
13.
Fukazawa  TYamasaki  KIto  H  et al.  Both the HLA-DPB1 and -DRB1 alleles correlate with risk for multiple sclerosis in Japanese: clinical phenotypes and gender as important factors.  Tissue Antigens 2000;55 (3) 199- 205PubMedGoogle ScholarCrossref
14.
Fukazawa  TKikuchi  SMiyagishi  R  et al.  HLA-DPB1*0501 is not uniquely associated with opticospinal multiple sclerosis in Japanese patients: important role of DPB1*0301.  Mult Scler 2006;12 (1) 19- 23PubMedGoogle ScholarCrossref
15.
Shibasaki  HMcDonald  WIKuroiwa  Y Racial modification of clinical picture of multiple sclerosis: comparison between British and Japanese patients.  J Neurol Sci 1981;49 (2) 253- 271PubMedGoogle ScholarCrossref
16.
Matsuoka  TMatsushita  TOsoegawa  M  et al.  Association of the HLA-DRB1 alleles with characteristic MRI features of Asian multiple sclerosis.  Mult Scler 2008;14 (9) 1181- 1190PubMedGoogle ScholarCrossref
17.
Deen  SBacchetti  PHigh  AWaubant  E Predictors of the location of multiple sclerosis relapse.  J Neurol Neurosurg Psychiatry 2008;79 (10) 1190- 1193PubMedGoogle ScholarCrossref
18.
Barcellos  LFOksenberg  JRGreen  AJ  et al. Multiple Sclerosis Genetics Group, Genetic basis for clinical expression in multiple sclerosis.  Brain 2002;125 (pt 1) 150- 158PubMedGoogle ScholarCrossref
19.
Butterfield  RJBlankenhorn  EPRoper  RJZachary  JFDoerge  RWTeuscher  C Identification of genetic loci controlling the characteristics and severity of brain and spinal cord lesions in experimental allergic encephalomyelitis.  Am J Pathol 2000;157 (2) 637- 645PubMedGoogle ScholarCrossref
20.
Poser  CMPaty  DWScheinberg  L  et al.  New diagnostic criteria for multiple sclerosis: guidelines for research protocols.  Ann Neurol 1983;13 (3) 227- 231PubMedGoogle ScholarCrossref
21.
Polman  CHReingold  SCEdan  G  et al.  Diagnostic criteria for multiple sclerosis: 2005 revisions to the “McDonald criteria.”  Ann Neurol 2005;58 (6) 840- 846PubMedGoogle ScholarCrossref
22.
Kurtzke  JF Rating neurologic impairment in multiple sclerosis: an Expanded Disability Status Scale (EDSS).  Neurology 1983;33 (11) 1444- 1452PubMedGoogle ScholarCrossref
23.
Cutter  GRBaier  MLRudick  RA  et al.  Development of multiple sclerosis functional composite as a clinical trial outcome measure.  Brain 1999;122 (pt 5) 871- 882PubMedGoogle ScholarCrossref
24.
Tejedor  DCastillo  SMozas  P  et al. Spanish FH Group, Reliable low-density DNA array based on allele-specific probes for detection of 118 mutations causing familial hypercholesterolemia.  Clin Chem 2005;51 (7) 1137- 1144PubMedGoogle ScholarCrossref
25.
van Walderveen  MABarkhof  FHommes  OR  et al.  Correlating MRI and clinical disease activity in multiple sclerosis: relevance of hypointense lesions on short-TR/short-TE (T1-weighted) spin-echo images.  Neurology 1995;45 (9) 1684- 1690PubMedGoogle ScholarCrossref
26.
Bot  JCBarkhof  FPolman  CH  et al.  Spinal cord abnormalities in recently diagnosed MS patients: added value of spinal MRI examination.  Neurology 2004;62 (2) 226- 233PubMedGoogle ScholarCrossref
27.
Benjamini  YHochberg  Y Controlling the false discovery rate: a practical and powerful approach to multiple testing.  J R Stat Soc Ser B 1995;57289- 300Google Scholar
28.
de Bakker  PIMcVean  GSabeti  PC  et al.  A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC.  Nat Genet 2006;38 (10) 1166- 1172PubMedGoogle ScholarCrossref
29.
Hafler  DACompston  ASawcer  S  et al. International Multiple Sclerosis Genetics Consortium, Risk alleles for multiple sclerosis identified by a genomewide study.  N Engl J Med 2007;357 (9) 851- 862PubMedGoogle ScholarCrossref
30.
Haines  JLTer-Minassian  MBazyk  A  et al. Multiple Sclerosis Genetics Group, A complete genomic screen for multiple sclerosis underscores a role for the major histocompatibility complex.  Nat Genet 1996;13 (4) 469- 471PubMedGoogle ScholarCrossref
31.
Zivadinov  RUxa  LZacchi  T  et al.  HLA genotypes and disease severity assessed by magnetic resonance imaging findings in patients with multiple sclerosis.  J Neurol 2003;250 (9) 1099- 1106PubMedGoogle ScholarCrossref
32.
Honig  LSSheremata  WA Magnetic resonance imaging of spinal cord lesions in multiple sclerosis.  J Neurol Neurosurg Psychiatry 1989;52 (4) 459- 466PubMedGoogle ScholarCrossref
33.
Lycklama à Nijeholt  GJBarkhof  FScheltens  P  et al.  MR of the spinal cord in multiple sclerosis: relation to clinical subtype and disability.  AJNR Am J Neuroradiol 1997;18 (6) 1041- 1048PubMedGoogle Scholar
34.
Furby  JHayton  TAnderson  V  et al.  Magnetic resonance imaging measures of brain and spinal cord atrophy correlate with clinical impairment in secondary progressive multiple sclerosis.  Mult Scler 2008;14 (8) 1068- 1075PubMedGoogle ScholarCrossref
35.
Okuda  DTSrinivasan  ROksenberg  JR  et al.  Genotype-phenotype correlations in multiple sclerosis: HLA genes influence disease severity inferred by 1HMR spectroscopy and MRI measures.  Brain 2009;132 (pt 1) 250- 259PubMedGoogle ScholarCrossref
36.
van Veen  TNielsen  JBerkhof  J  et al.  CCL5 and CCR5 genotypes modify clinical, radiological and pathological features of multiple sclerosis.  J Neuroimmunol 2007;190 (1-2) 157- 164PubMedGoogle ScholarCrossref
37.
Zivadinov  RUxa  LBratina  A  et al.  HLA-DRB1*1501, -DQB1*0301, -DQB1*0302, -DQB1*0602, and -DQB1*0603 alleles are associated with more severe disease outcome on MRI in patients with multiple sclerosis.  Int Rev Neurobiol 2007;79521- 535PubMedGoogle Scholar
38.
Sepulcre  JMasdeu  JCPalacios  R  et al.  HLA-DR2 and white matter lesion distribution in MS.  J Neuroimaging 2008;18 (3) 328- 331PubMedGoogle ScholarCrossref
39.
Fernando  MMStevens  CRWalsh  EC  et al.  Defining the role of the MHC in autoimmunity: a review and pooled analysis.  PLoS Genet 2008;4 (4) e1000024PubMedGoogle ScholarCrossref
40.
Storch  MKStefferl  ABrehm  U  et al.  Autoimmunity to myelin oligodendrocyte glycoprotein in rats mimics the spectrum of multiple sclerosis pathology.  Brain Pathol 1998;8 (4) 681- 694PubMedGoogle ScholarCrossref
41.
Stromnes  IMCerretti  LMLiggitt  DHarris  RAGoverman  JM Differential regulation of central nervous system autoimmunity by TH1 and TH17 cells.  Nat Med 2008;14 (3) 337- 342PubMedGoogle ScholarCrossref
Original Contribution
December 2009

HLA-DRB1*1501 and Spinal Cord Magnetic Resonance Imaging Lesions in Multiple Sclerosis

Author Affiliations

Author Affiliations: Departments of Neurology (Drs Sombekke, Killestein, Uitdehaag, and Polman), Radiology (Drs Lukas, Geurts, Barkhof, and Vrenken), Pathology (Drs Crusius, Peña, and Geurts), Epidemiology and Biostatistics (Drs Uitdehaag and Knol), and Physics and Medical Technology (Dr Vrenken), VU University Medical Center, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands; Department of Diagnostic and Interventional Radiology, St Josef Hospital, Ruhr University, Bochum, Germany (Dr Lukas); Progenika Biopharma, SA, Derio, Spain (Drs Tejedor, Arteta, and Martínez); and Center for Neurologic Diseases, Department of Neurology, Brigham and Women's Hospital and Harvard Medical School, and Partners Center for Personalized Genetic Medicine, Boston, Massachusetts (Dr De Jager).

Arch Neurol. 2009;66(12):1531-1536. doi:10.1001/archneurol.2009.278
Abstract

Background  Multiple sclerosis (MS) is a heterogeneous neurologic disease with extensive variation with respect to the most affected central nervous system region (brain vs spinal cord).

Objective  To test the hypothesis that this variation in lesion location (brain vs spinal cord) might be (partially) genetically determined.

Design  Candidate gene study.

Setting  Academic research.

Patients  Patients were selected for the availability of DNA material, clinical variables, and brain and spinal cord magnetic resonance images (evaluating T2-weighted lesion load in the brain and the number of spinal cord lesions).

Main Outcome Measures  For genotyping, we used a DNA chip containing a set of genes mentioned in previous publications noting their relation to different phenotypes of MS. We assessed the association between brain and spinal cord abnormalities and the genotypes of the patients.

Results  One hundred fifty patients were included in the analysis. Five single-nucleotide polymorphisms within the major histocompatibility complex region were associated with the number of focal abnormalities in the spinal cord. The most significant was rs3135388 (surrogate marker for the HLA-DRB1*1501 allele). Carriers of HLA-DRB1*1501 had a median of 4 spinal cord lesions compared with 2 lesions for noncarriers (P < .001). No significant association was noted between the single-nucleotide polymorphisms and T2-weighted lesion load in the brain.

Conclusions  Carriership of HLA-DRB1*1501 (via rs3135388) was associated with the extent of focal abnormalities in the spinal cord. Spinal cord lesions might be an explanation for increased MS disease severity in patients carrying HLA-DRB1*1501.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system characterized by inflammation, demyelination, and axonal loss in the brain and spinal cord. The current hypothesis is that MS is caused by a complex interplay of genetic and environmental factors. The genetic influence is characterized by the interaction of multiple genes that exert modest effects. The most striking among these is the association of the major histocompatibility complex (MHC) with MS susceptibility, which has consistently been reported over the past decades.1

Clinically, MS is a heterogeneous disease with a diverse spectrum of neurologic deficits and variable outcome. Some studies1-8 have focused on genetic predictors of disease phenotypes (such as disease severity, disease subtypes, magnetic resonance [MR] imaging characteristics, and response to treatment), showing some effect of different genes.

Besides clinical heterogeneity, considerable variation exists between patients regarding type and anatomic location of the lesions.9,10 Varying degrees of involvement of the cerebrum, brainstem, cerebellum, and spinal cord have been described in postmortem and MR imaging studies.11,12 This variability is unexplained, and much effort is being put in unraveling it. Several arguments point toward a genetic role in the regional distribution of lesions in the brain and spinal cord in MS. First, differences in pathologic manifestation of MS have been observed between Asian and Western populations. Asian-type MS is characterized by predominant involvement of the optic nerve and spinal cord, whereas Western-type MS predominantly involves the brain. Genetic differences within the HLA region between these populations could underlie this variability.3,13-16 Second, recent studies17,18 confirm a tendency of patients to develop relapses in the same locations, including spinal cord lesions. This tendency for localized exacerbations could be genetically predetermined. Third, involvement of a genetic factor has already been demonstrated in lesion distribution in experimental autoimmune encephalomyelitis (an animal model that shows many similarities to MS); Butterfield et al19 showed that in mouse experimental autoimmune encephalomyelitis the lesions in the brain and spinal cord were controlled by different quantitative trait loci.

Many researchers have demonstrated the clinical relevance of spinal cord lesions in diagnosing MS and disability accumulation. However, we are unaware of any studies that included spinal cord MR imaging variables in genetic analyses searching for genotype-phenotype correlations.

In this study, we used a DNA chip containing a set of single-nucleotide polymorphisms (SNPs) in candidate genes to assess the genetic effect on regional lesion distribution in the brain and spinal cord as seen on MR imaging. The SNPs on this chip were selected based on published associations with MS pathogenesis, prognosis, or response to treatment.

Methods
Study participants

Unrelated Dutch Caucasian patients were selected retrospectively from natural history studies conducted in the Department of Neurology and Radiology, Vrije Universiteit Medical Center, Amsterdam, the Netherlands. Patients were selected for the availability of DNA material, as well as brain and spinal cord MR images that fulfilled certain standardized requirements (described herein) and were obtained less than 2 years apart. The study was performed with the approval of the medical ethics committee of the VU University Medical Center, and informed consent was obtained from all participants. Patients were diagnosed as having MS as ascertained by Poser et al20 or revised McDonald criteria.21 For patients included in the analyses, clinical data were collected, including age, sex, type of disease, age at onset, disease course, and duration of disease. Disability status was determined for all subjects using Kurtzke's Expanded Disability Status Scale22 and, whenever available, the Multiple Sclerosis Functional Composite Scale.23

SELECTION OF SNPs

Single-nucleotide polymorphisms were selected based on involvement in MS pathogenesis, prognosis, or response to treatment according to the literature published before July 2007. The polymorphisms were confirmed and associated with an identifier using dbSNP database (http://www.ncbi.nlm.nih.gov/SNP). Nucleotide sequences for the design of allele-specific probes and polymerase chain reaction primers were retrieved in the SNPper database (http://snpper.chip.org/bio). Sequence-specific probes and primers were designed using freely available Primer3 software (http://frodo.wi.mit.edu/).

If a polymorphism was not present in the database, position and sequences were established by performing a BLAST search (http://blast.ncbi.nlm.nih.gov/Blast.cgi) using the data available in the literature.

Genotyping

Genomic DNA was isolated from anticoagulated blood using DNAzol reagent (Molecular Research Center, Inc, Cincinnati, Ohio). Genotyping was performed using a newly developed low-density DNA microarray based on allele-specific probes. The design, fabrication, validation, and analysis of the arrays were performed using the procedure described by Tejedor et al24 with minor modifications.

Brain mr imaging

Scans were acquired on 1.0-T or 1.5-T scanners (Siemens AG, Erlangen, Germany) with standard head coils using 2-dimensional conventional or fast spin-echo proton density–weighted and T2-weighted images (repetition time, 2200-3000 milliseconds; echo times, 20-30 and 80-100 milliseconds with section thicknesses of 3 to 5 mm, a maximum gap between sections of 0.5 mm, and an in-plane resolution of 1 × 1 mm2. Lesions were identified by an expert reader and then outlined on the corresponding proton density–weighted image using home-developed semiautomated seed-growing software based on a local thresholding technique (Show_images).25 Lesion areas were multiplied by the intersection distance to obtain total T2-weighted brain lesion volume for each patient.

Spinal cord mr imaging

Spinal cord scanning included a cardiac-triggered sagittal proton density–weighted and T2-weighted dual-echo spin-echo sequence with a section thickness of 3 mm covering the whole spinal cord (repetition time, 2500-3000 milliseconds; echo times, 20-30 and 80-100 milliseconds), with a gap between sections of 0.3 mm, and an in-plane resolution of 1 × 1 mm2. From this sequence, the number of focal lesions and the presence of diffuse abnormalities were scored by an experienced reader (C.L.). Diffuse abnormalities were defined as poorly delineated areas with increased signal intensity compared with the signal intensity of spinal cerebrospinal fluid best seen on proton density–weighted images.26 Lesion volume was assessed in the spinal cord as the total number of focal pathologic segments involved (total extension of all spinal cord lesions over several corresponding vertebral segments).

Statistical analysis

First, the associations between the brain variable (T2-weighted lesion load) and the spinal cord variables (the number of focal lesions and the presence of diffuse abnormalities) were tested per SNP and per clinical variable. We used the nonparametric Spearman ρ rank correlation test, Kruskal-Wallis test, Mann-Whitney test, and χ2 test as appropriate, applying the false discovery rate according to Benjamini and Hochberg27 to correct for multiple testing. The corrected number represents the expected proportion of false discoveries for a given P value cutoff. We used the cutoff point of P < .05 after false discovery rate correction. Spearman rank correlation coefficient was used to test the correlations between 2 scaled variables.

Second, multivariate analysis by general linear models was performed on the associated SNPs (P < .10 significance after correction for multiple testing). We used log-transformed brain T2-weighted lesion loads and the square root–transformed number of focal abnormalities in the spinal cord to correct for significant clinical variables and the type of MR imaging system (1.0 T vs 1.5 T).

All analyses were performed using commercially available software packages. These included SPSS (version 15; SPSS Inc, Chicago, Illinois), Excel 2003 (Microsoft Corporation, Redmond, Washington), and HelixTree (Golden Helix, Inc, Bozeman, Montana).

Results
Patient characteristics

One hundred fifty patients were included in the analysis. Our patient group reflects a representative population with MS, with 36.7% being male and 20.0% having primary progressive MS (Table 1). Most patients (132 of 150) demonstrated abnormalities (mostly focal lesions) on spinal cord MR imaging, while all patients had abnormalities on brain MR imaging.

Genotyping

In total, 80 SNPs in 44 genes were studied on the DNA chip. Twelve SNPs were excluded from further analysis (5 SNPs were monomorphic, and 7 SNPs had a minor allele frequency <5% [eTable]). Hardy-Weinberg equilibrium was calculated for all SNPs.

Correlation between lesion load in the brain and genotypes

In the univariate analyses of T2-weighted lesion load in the brain vs all SNPs on the DNA chip, the only suggestive correlation was rs2107538 (CCL5) (OMIM 187011) (Table 2). Two clinical covariates were significant and were included in the general linear model for brain T2-weighted lesion load, namely, disease duration and MS subtype. The type of MR imaging system (1.0 T vs 1.5 T) was not associated with T2-weighted lesion load in the brain. After inclusion of rs2107538 (CCL5) in the model that contained the clinical covariates, this SNP showed significant association with brain T2-weighted lesion load (P = .03).

Correlation between spinal cord abnormalities and genotypes

Several MHC SNPs were found to be related to the number of focal spinal cord abnormalities (Table 2). The most significant is SNP rs3135388. Carrier status of the A allele (surrogate marker for HLA-DRB1*1501) was associated with significantly more lesions in the spinal cord (Figure). The median number of focal abnormalities in carriers of HLA-DRB1*1501 (n = 80) was 4 lesions and in noncarriers (n = 70) was 2 lesions (P < .001, Mann-Whitney test).

When corrected for multiple testing, 5 SNPs within the MHC class II region (rs3135388, rs2395182, rs2239802, rs2227139, and rs2213584) remained significant, and 1 SNP within the C2TA gene (class II transactivator) showed a suggestive result. The 5 HLA SNPs are in strong linkage disequilibrium.

In addition, the aforementioned MHC class II SNPs were also found to be correlated with a higher lesion volume in the spinal cord. Specifically, carriership of HLA-DRB1*1501 was associated with more segments of the spinal cord affected by focal lesions (P = .01, Mann-Whitney test).

No clinical covariate (disease duration, age of patients, or MS subtype) was significantly correlated with the number of focal lesions in the spinal cord. No general linear model that included the 5 significantly associated MHC class II SNPs could be produced because of the high level of collinearity owing to the high linkage disequilibrium (r2 range, 0.33-0.99). Only rs3087456 (C2TA) (OMIM 600005) showed a marginally significant P value (P = .05) in the general linear model that included rs3135388 (P = .002). No association was observed between the presence of diffuse abnormalities and the evaluated SNPs.

Comment

In the present candidate gene study, we observed an association between several SNPs within the MHC class II region and the number of focal abnormalities in the spinal cord. The most significant results were found for rs3135388. This SNP is a surrogate marker for the HLA-DRB1*1501 allele (r2 = 0.97).28

After correction for multiple testing, no genes were significantly related to T2-weighted lesion load in the brain. The MHC class II SNPs that were associated with more lesions in the spinal cord and with a higher lesion volume in the spinal cord do not seem to affect the presence of T2-weighted lesions in the brain.

Studies6,29,30 have unambiguously shown that HLA-DRB1*1501 strongly influences MS susceptibility. There is some evidence that this haplotype might also be associated with a more severe disease course.1,4,6,31 Recent findings show a correlation between carriership of HLA-DRB1*1501 and Expanded Disability Status Scale scores using an extremes of outcome strategy (comparing patients with severe and mild disability only).4 In our study, we could not demonstrate a correlation between carriership of HLA-DRB1*1501 and Expanded Disability Status Scale scores. This might be because of the impossibility of performing the extremes of outcome analysis owing to a small sample size. However, findings from this study suggest that spinal cord lesions might be an additional explanation for the described relationship between HLA-DRB1*1501 and MS disease severity.

Abnormalities in the spinal cord correlate with the degree of disability and with the date of diagnosis.32,33 In our study, we also noticed an association between Expanded Disability Status Scale scores and the number of focal abnormalities in the spinal cord (P = .02, ρ = 0.185). Future studies should include additional MR imaging variables of the spinal cord such as atrophy (cross-sectional area), as this variable might correlate better with clinical disability.34

In a recent article, Okuda et al35 showed in their patients that HLA-DRB1*1501 affects disease severity as measured by clinical variables and by brain MR imaging variables. In our patients, we could not demonstrate this effect on T2-weighted lesion load in the brain or on cognitive function (using PASAT [Paced Auditory Serial Addition Task]23 scores). Post hoc power analysis detected a power of less than 25% to detect a T2-weighted lesion load difference of 1 mL, suggesting that a type II error cannot be ruled out. This also warrants careful interpretation of the suggestive positive finding of rs2107538 (CCL5) in association with brain T2-weighted lesion load, although previous findings have shown an influence of this gene on other MR imaging variables.36

Previous data were sometimes conflicting about the effect of HLA haplotypes on brain MR imaging features. Some studies31,35 found a relationship between HLA genotypes and brain MR imaging quantitative markers, while other studies31,37 did not observe this effect of HLA-DRB1*1501. In addition, no effect of HLA-DRB1*1501 was shown on lesion distribution within the brain using T1-weighted and T2-weighted lesion maps.38 To date, we are unaware of any study that specifically relates HLA-DRB1*1501 to spinal cord abnormalities on MR imaging in a Caucasian population.

The molecular basis of the association with greater involvement of the spinal cord in patients carrying the HLA-DRB1*1501 allele remains unknown. Major histocompatibility complex class II genes are involved in self vs nonself immune recognition.39 These genes encode for polymorphic surface glycoproteins. Variability in this region may determine individual differences in T-cell responses. It is postulated that MHC class II may present variable central nervous system antigens to T cells, possibly producing different lesion distribution in animals and in humans.40

Moreover, Stromnes et al41 reported different lesion distribution in the central nervous system (spinal cord vs brain parenchyma) in 2 mouse experimental autoimmune encephalomyelitis models with different MHC strains. This was found to be mediated by variable preferential MOG (myelin oligodendrocyte glycoprotein) epitope presentation and ultimately by an alternative ratio of T-helper type 17 to T-helper type 1. The present study indicates a different mechanism of lesion formation in the brain vs the spinal cord, with a possible indirect role of the MHC class II genes. The role of the genes on this mechanism is unclear. The translation of these findings from animal models to the human situation warrants further studies.

Because of the complexity of the MHC class II region (epistatic effects and high linkage disequilibrium), future studies using high-density HLA mapping are warranted to unravel the genetic influence on lesion distribution and to determine the responsible HLA allele.

In conclusion, this study demonstrates more focal lesions in the spinal cord among carriers of the HLA-DRB1*1501 allele. If confirmed in independent samples, these observations may provide important insight in MS disease heterogeneity and its underlying mechanisms.

Correspondence: Madeleine H. Sombekke, MD, Department of Neurology, VU University Medical Center, Vrije Universiteit Amsterdam, De Boelelaan 1117, 1081 HV Amsterdam, the Netherlands (m.sombekke@vumc.nl).

Accepted for Publication: June 12, 2009.

Author Contributions:Study concept and design: Sombekke, Lukas, Crusius, Tejedor, Killestein, Arteta, Martínez, Uitdehaag, Peña, Vrenken, and Polman. Acquisition of data: Sombekke, Lukas, Crusius, and Barkhof. Analysis and interpretation of data: Sombekke, Lukas, Crusius, Tejedor, Killestein, Arteta, Martínez, Uitdehaag, Knol, Geurts, De Jager, Barkhof, Vrenken, and Polman. Drafting of the manuscript: Sombekke, Crusius, Tejedor, Arteta, Martínez, and Vrenken. Critical revision of the manuscript for important intellectual content: Crusius, Tejedor, Killestein, Arteta, Martínez, Uitdehaag, Knol, Peña, Geurts, De Jager, Barkhof, Vrenken, and Polman. Statistical analysis: Sombekke, Uitdehaag, Knol, and De Jager. Obtained funding: Tejedor, Arteta, Martínez, and Polman. Administrative, technical, and material support: Crusius, Tejedor, Arteta, Martínez, and Vrenken. Study supervision: Crusius, Uitdehaag, Barkhof, Vrenken, and Polman.

Financial Disclosure: Dr Lukas received a research grant from Bayer Schering Germany. Drs Tejedor, Arteta, and Martínez are employees of the corporate sponsor. Dr Killestein has been involved in clinical trials of companies that market drugs for multiple sclerosis (Bayer Schering Germany, AG, Biogen Idec, Merck Serono, and Teva Neuroscience) and some companies that have development programs for future drugs to treat multiple sclerosis. Dr Uitdehaag received honoraria for consultancy from Novartis and Merck Serono. Dr De Jager received speaking and advisory board fees from Biogen Idec, honoraria for serving on the advisory board of Teva Neuroscience, and research support from Affymetrix Inc and Source MDx Inc. Dr Vrenken received program grant 05-358c and project grants 05-566MS and 08-633MS from the Dutch Multiple Sclerosis Research Foundation and project grant EADNI-06-10004 from the Alzheimer's Association. Dr Polman received consulting fees or speaking fees from Actelion, Biogen Idec, Bayer Schering Germany, Teva Neuroscience, Merck Serono, Novartis, GlaxoSmithKline, UCB, Roche, and Antisense Therapeutics, as well as grant support from Biogen Idec, Bayer Schering Germany, GlaxoSmithKline, Novartis, UCB, Merck Serono, and Teva Neuroscience.

Funding/Support: This study was supported by Progenika Biopharma, SA. The MS Center Amsterdam is supported by program grant 05-358 from the Dutch Multiple Sclerosis Research Foundation. Dr De Jager is a Harry Weaver Neuroscience Scholar of the National Multiple Sclerosis Society.

Role of the Sponsor: The corporate sponsor participated in the design and conduct of the study (including selection of single-nucleotide polymorphisms for the DNA chip, which is intellectually protected by the corporate sponsor), the genetic analysis, and the review and approval of the final manuscript.

Additional Contributions: M. van de Wiel, PhD, provided statistical advice. We thank all the patients for their participation.

References
1.
Barcellos  LFOksenberg  JRBegovich  AB  et al. Multiple Sclerosis Genetics Group, HLA-DR2 dose effect on susceptibility to multiple sclerosis and influence on disease course.  Am J Hum Genet 2003;72 (3) 710- 716PubMedGoogle ScholarCrossref
2.
Hoffmann  SCepok  SGrummel  V  et al.  HLA-DRB1*0401 and HLA-DRB1*0408 are strongly associated with the development of antibodies against interferon-beta therapy in multiple sclerosis [published correction appears in Am J Hum Genet. 2008;83(4):541].  Am J Hum Genet 2008;83 (2) 219- 227PubMedGoogle ScholarCrossref
3.
Yamasaki  KHoriuchi  IMinohara  M  et al.  HLA-DPB1*0501-associated opticospinal multiple sclerosis: clinical, neuroimaging and immunogenetic studies.  Brain 1999;122 (pt 9) 1689- 1696PubMedGoogle ScholarCrossref
4.
DeLuca  GCRamagopalan  SVHerrera  BM  et al.  An extremes of outcome strategy provides evidence that multiple sclerosis severity is determined by alleles at the HLA-DRB1 locus.  Proc Natl Acad Sci U S A 2007;104 (52) 20896- 20901PubMedGoogle ScholarCrossref
5.
Kantarci  OHde Andrade  MWeinshenker  BG Identifying disease modifying genes in multiple sclerosis.  J Neuroimmunol 2002;123 (1-2) 144- 159PubMedGoogle ScholarCrossref
6.
Barcellos  LFSawcer  SRamsay  PP  et al.  Heterogeneity at the HLA-DRB1 locus and risk for multiple sclerosis.  Hum Mol Genet 2006;15 (18) 2813- 2824PubMedGoogle ScholarCrossref
7.
Fazekas  FEnzinger  CRopele  SSchmidt  HSchmidt  RStrasser-Fuchs  S The impact of our genes: consequences of the apolipoprotein E polymorphism in Alzheimer disease and multiple sclerosis.  J Neurol Sci 2006;245 (1-2) 35- 39PubMedGoogle ScholarCrossref
8.
van Veen  Tvan Winsen  LCrusius  JB  et al.  αB-crystallin genotype has impact on the multiple sclerosis phenotype.  Neurology 2003;61 (9) 1245- 1249PubMedGoogle ScholarCrossref
9.
Lucchinetti  CBrück  WParisi  JScheithauer  BRodriguez  MLassmann  H Heterogeneity of multiple sclerosis lesions: implications for the pathogenesis of demyelination.  Ann Neurol 2000;47 (6) 707- 717PubMedGoogle ScholarCrossref
10.
Lee  MASmith  SPalace  J  et al.  Spatial mapping of T2 and gadolinium-enhancing T1 lesion volumes in multiple sclerosis: evidence for distinct mechanisms of lesion genesis?  Brain 1999;122 (pt 7) 1261- 1270PubMedGoogle ScholarCrossref
11.
Nijeholt  GJBergers  EKamphorst  W  et al.  Post-mortem high-resolution MRI of the spinal cord in multiple sclerosis: a correlative study with conventional MRI, histopathology and clinical phenotype.  Brain 2001;124 (pt 1) 154- 166PubMedGoogle ScholarCrossref
12.
Bakshi  RThompson  AJRocca  MA  et al.  MRI in multiple sclerosis: current status and future prospects.  Lancet Neurol 2008;7 (7) 615- 625PubMedGoogle ScholarCrossref
13.
Fukazawa  TYamasaki  KIto  H  et al.  Both the HLA-DPB1 and -DRB1 alleles correlate with risk for multiple sclerosis in Japanese: clinical phenotypes and gender as important factors.  Tissue Antigens 2000;55 (3) 199- 205PubMedGoogle ScholarCrossref
14.
Fukazawa  TKikuchi  SMiyagishi  R  et al.  HLA-DPB1*0501 is not uniquely associated with opticospinal multiple sclerosis in Japanese patients: important role of DPB1*0301.  Mult Scler 2006;12 (1) 19- 23PubMedGoogle ScholarCrossref
15.
Shibasaki  HMcDonald  WIKuroiwa  Y Racial modification of clinical picture of multiple sclerosis: comparison between British and Japanese patients.  J Neurol Sci 1981;49 (2) 253- 271PubMedGoogle ScholarCrossref
16.
Matsuoka  TMatsushita  TOsoegawa  M  et al.  Association of the HLA-DRB1 alleles with characteristic MRI features of Asian multiple sclerosis.  Mult Scler 2008;14 (9) 1181- 1190PubMedGoogle ScholarCrossref
17.
Deen  SBacchetti  PHigh  AWaubant  E Predictors of the location of multiple sclerosis relapse.  J Neurol Neurosurg Psychiatry 2008;79 (10) 1190- 1193PubMedGoogle ScholarCrossref
18.
Barcellos  LFOksenberg  JRGreen  AJ  et al. Multiple Sclerosis Genetics Group, Genetic basis for clinical expression in multiple sclerosis.  Brain 2002;125 (pt 1) 150- 158PubMedGoogle ScholarCrossref
19.
Butterfield  RJBlankenhorn  EPRoper  RJZachary  JFDoerge  RWTeuscher  C Identification of genetic loci controlling the characteristics and severity of brain and spinal cord lesions in experimental allergic encephalomyelitis.  Am J Pathol 2000;157 (2) 637- 645PubMedGoogle ScholarCrossref
20.
Poser  CMPaty  DWScheinberg  L  et al.  New diagnostic criteria for multiple sclerosis: guidelines for research protocols.  Ann Neurol 1983;13 (3) 227- 231PubMedGoogle ScholarCrossref
21.
Polman  CHReingold  SCEdan  G  et al.  Diagnostic criteria for multiple sclerosis: 2005 revisions to the “McDonald criteria.”  Ann Neurol 2005;58 (6) 840- 846PubMedGoogle ScholarCrossref
22.
Kurtzke  JF Rating neurologic impairment in multiple sclerosis: an Expanded Disability Status Scale (EDSS).  Neurology 1983;33 (11) 1444- 1452PubMedGoogle ScholarCrossref
23.
Cutter  GRBaier  MLRudick  RA  et al.  Development of multiple sclerosis functional composite as a clinical trial outcome measure.  Brain 1999;122 (pt 5) 871- 882PubMedGoogle ScholarCrossref
24.
Tejedor  DCastillo  SMozas  P  et al. Spanish FH Group, Reliable low-density DNA array based on allele-specific probes for detection of 118 mutations causing familial hypercholesterolemia.  Clin Chem 2005;51 (7) 1137- 1144PubMedGoogle ScholarCrossref
25.
van Walderveen  MABarkhof  FHommes  OR  et al.  Correlating MRI and clinical disease activity in multiple sclerosis: relevance of hypointense lesions on short-TR/short-TE (T1-weighted) spin-echo images.  Neurology 1995;45 (9) 1684- 1690PubMedGoogle ScholarCrossref
26.
Bot  JCBarkhof  FPolman  CH  et al.  Spinal cord abnormalities in recently diagnosed MS patients: added value of spinal MRI examination.  Neurology 2004;62 (2) 226- 233PubMedGoogle ScholarCrossref
27.
Benjamini  YHochberg  Y Controlling the false discovery rate: a practical and powerful approach to multiple testing.  J R Stat Soc Ser B 1995;57289- 300Google Scholar
28.
de Bakker  PIMcVean  GSabeti  PC  et al.  A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC.  Nat Genet 2006;38 (10) 1166- 1172PubMedGoogle ScholarCrossref
29.
Hafler  DACompston  ASawcer  S  et al. International Multiple Sclerosis Genetics Consortium, Risk alleles for multiple sclerosis identified by a genomewide study.  N Engl J Med 2007;357 (9) 851- 862PubMedGoogle ScholarCrossref
30.
Haines  JLTer-Minassian  MBazyk  A  et al. Multiple Sclerosis Genetics Group, A complete genomic screen for multiple sclerosis underscores a role for the major histocompatibility complex.  Nat Genet 1996;13 (4) 469- 471PubMedGoogle ScholarCrossref
31.
Zivadinov  RUxa  LZacchi  T  et al.  HLA genotypes and disease severity assessed by magnetic resonance imaging findings in patients with multiple sclerosis.  J Neurol 2003;250 (9) 1099- 1106PubMedGoogle ScholarCrossref
32.
Honig  LSSheremata  WA Magnetic resonance imaging of spinal cord lesions in multiple sclerosis.  J Neurol Neurosurg Psychiatry 1989;52 (4) 459- 466PubMedGoogle ScholarCrossref
33.
Lycklama à Nijeholt  GJBarkhof  FScheltens  P  et al.  MR of the spinal cord in multiple sclerosis: relation to clinical subtype and disability.  AJNR Am J Neuroradiol 1997;18 (6) 1041- 1048PubMedGoogle Scholar
34.
Furby  JHayton  TAnderson  V  et al.  Magnetic resonance imaging measures of brain and spinal cord atrophy correlate with clinical impairment in secondary progressive multiple sclerosis.  Mult Scler 2008;14 (8) 1068- 1075PubMedGoogle ScholarCrossref
35.
Okuda  DTSrinivasan  ROksenberg  JR  et al.  Genotype-phenotype correlations in multiple sclerosis: HLA genes influence disease severity inferred by 1HMR spectroscopy and MRI measures.  Brain 2009;132 (pt 1) 250- 259PubMedGoogle ScholarCrossref
36.
van Veen  TNielsen  JBerkhof  J  et al.  CCL5 and CCR5 genotypes modify clinical, radiological and pathological features of multiple sclerosis.  J Neuroimmunol 2007;190 (1-2) 157- 164PubMedGoogle ScholarCrossref
37.
Zivadinov  RUxa  LBratina  A  et al.  HLA-DRB1*1501, -DQB1*0301, -DQB1*0302, -DQB1*0602, and -DQB1*0603 alleles are associated with more severe disease outcome on MRI in patients with multiple sclerosis.  Int Rev Neurobiol 2007;79521- 535PubMedGoogle Scholar
38.
Sepulcre  JMasdeu  JCPalacios  R  et al.  HLA-DR2 and white matter lesion distribution in MS.  J Neuroimaging 2008;18 (3) 328- 331PubMedGoogle ScholarCrossref
39.
Fernando  MMStevens  CRWalsh  EC  et al.  Defining the role of the MHC in autoimmunity: a review and pooled analysis.  PLoS Genet 2008;4 (4) e1000024PubMedGoogle ScholarCrossref
40.
Storch  MKStefferl  ABrehm  U  et al.  Autoimmunity to myelin oligodendrocyte glycoprotein in rats mimics the spectrum of multiple sclerosis pathology.  Brain Pathol 1998;8 (4) 681- 694PubMedGoogle ScholarCrossref
41.
Stromnes  IMCerretti  LMLiggitt  DHarris  RAGoverman  JM Differential regulation of central nervous system autoimmunity by TH1 and TH17 cells.  Nat Med 2008;14 (3) 337- 342PubMedGoogle ScholarCrossref
×