Variation of the Myelin Oligodendrocyte Glycoprotein gene is not primarily associated with multiple sclerosis in the Sardinian population
BMC Genetics volume 8, Article number: 25 (2007)
Multiple sclerosis (MS) is consistently associated with particular HLA-DRB1-DQB1 haplotypes. However, existing evidence suggests that variation at these loci does not entirely explain association of the HLA region with the disease. The MOG locus is a prime positional and functional candidate for such additional predisposing effects but the analysis is complicated by the strong, albeit labyrinthine pattern of linkage disequilibrium in the region. Here we have assessed the association of MOG variation with MS in the Sardinian population to see if it represents an independent contributor to MS predisposition.
After re-sequencing the MOG gene in 21 healthy parents of MS patients we detected 134 variants, 33 of which were novel. A set of 40 informative SNPs was then selected and assessed for disease association together with 1 intragenic microsatellite in an initial data set of 239 MS families. This microsatellite and 11 SNPs were found to be positively associated with MS, using the transmission disequilibrium test, and were followed up in an additional 158 families (total families analysed = 397). While in these 397 families, 8 markers showed significant association with MS, through conditional tests we determined that these MOG variants were not associated with MS independently of the main DRB1-DQB1 disease associations.
These results indicate that variation within the MOG gene is not an important independent determinant of MS-inherited risk in the Sardinian population.
Multiple sclerosis (MS) is a serious chronic inflammatory and demyelinizating disorder of the central nervous system which results from an autoimmune attack on components of the oligodendrocyte cell . The disease is more common in European and European-derived populations. Within Europe it shows a north-south gradient, with the notable exception of the Mediterranean island of Sardinia, which has one of the highest prevalences worldwide . The disease risk for a monozygotic twin of an affected patient is about 30%, showing a quick fall-off rate with decreased genetic relatedness to affected individuals. These data, and the increasing incidence of disease reported in some populations over the last few decades [3, 4], suggest that the chance of this inflammatory process occurring depends on the complex interplay between a polygenic trait and unknown environmental factors influencing the penetrance of susceptibility genes [5, 6]. MS has been found to be consistently associated with specific HLA class II variants and notably with the HLA-DRB1*1501-DQB1*0602 haplotype, which represents the main risk factor for disease occurrence in different ethnic backgrounds [7–11]. The relative contribution of variation at the DRB1 and DQB1 loci to disease predisposition is still not completely clear although some studies based on cross-comparing rare HLA DRB1-DQB1 haplotype splits have indicated that the main contribution comes from variation at the DRB1 locus . Analysis of large data sets from Finland and Canada have also suggested that in those populations, any important additional modifiers of MS susceptibility were likely to be contained in the regions close to DRB1 . In Sardinia, the HLA-DRB1*1501-DQB1*0602 haplotype is rare but is still significantly positively associated with MS together with an additional 4 haplotypes including DRB1*1303-DQB1*0301, DRB1*0405-DQB1*0301, DRB1*0301-DQB1*0201 and DRB1*0405-DQB1*0302 . Some of these haplotypes, such as DRB1*0301-DQB1*0201 are also associated with MS in some non-Sardinian populations [11, 14, 15] but not in others [8, 9, 16–20]. The presence of HLA non-DRB1-DQB1 predisposing effects related to a different distribution of extended HLA-DRB1*0301-DQB1*0201 haplotypes in different populations can explain these findings. Indeed, some studies have suggested that within the HLA region there are further independent predisposing effects determined by as yet unidentified non-DRB1-DQB1 variants [13, 21–24]. However, strong linkage disequilibrium (LD) between the variants contained in the HLA region makes it difficult to detect which polymorphisms, outside the exon 2 sequences of the DR/DQ loci but within the HLA region, further influence disease risk. Within the HLA region, the MOG gene is a prime candidate for additional MS associations. This gene is located, 2.9 Mb telomeric of the DRB1 locus, in a chromosome interval which has shown some evidence of association with MS independent of DRB1-DQB1 . Furthermore, the rodent ortholog of this gene encodes for an autoantigen which triggers autoimmune responses in experimental models of disease [25, 26]. Moreover, in MS patients both T-cell and antibody responses against this protein have been detected [27, 28] and an aetiologic role of anti-MOG antibodies has been suggested in acute lesions of MS patients , although this role is controversial [30, 31]. These findings could be consistent with a model in which specific polymorphisms in the MOG gene could determine amino acid variation or differences in the level of expression of this protein in the central nervous system and affect immune responses against it, thus acting as primary aetiologic determinants of disease pathogenesis.
A few studies have tested the association between MOG polymorphisms and disease, with conflicting results. Some of these studies indicated a potential role of MOG [32, 33] while others failed to find any evidence of primary association with MS [34, 35, 12]. However, none of these studies analysed the MOG gene comprehensively, nor assessed the association of genetic variation of this gene in the context of extended DRB1-DQB1-MOG haplotypes.
The aim of this study was to establish the SNP content of the MOG gene in the Sardinian population and use this information together with appropriate conditional association tests to assess the individual contribution of the MOG gene to MS inherited risk to determine if it is an important independent risk locus.
Results and discussion
In order to assess the association of MOG variation with MS we first established polymorphism content of the MOG gene in the Sardinian population by resequencing the gene in 21 healthy parents of MS patients. The primers used are listed in Table 1, and technical details are reported in the Methods section. We detected 134 SNPs, 101 of which had already been reported  and 33 were novel. A set of 40 informative SNPs was then selected to be tested for disease association (see Methods). These were initially genotyped in a data set of 239 MS families. To identify the haplotype-tagging patterns, we computed the intermarker D' and r2 values of these 40 informative SNPs, computed for 478 parents of the MS patients, shown in Figure 1. The physical positions, with base 1165 corresponding to the translation start site of MOG, and other relevant information are indicated in Table 2. D' values indicate that these SNPs are part of the same linkage disequilibrium block (average 0.97 between contiguous marker pairs) that, as we have previously shown, encompasses HLA class I and III sub-regions and also includes the HLA class II DR and DQ loci [37, 13]. As expected, considering selected SNPs were chosen in order to avoid genetically redundant markers, the r2 values showed relatively low values with an average of 0.13 between contiguous marker pairs (Figure 1).
We then assessed disease association of the 41 variants of interest, including these 40 SNPs of interest and the MOG51 microsatellite, in a sample set of 239 families using the extended transmission disequilibrium test (ETDT). Without taking linkage disequilibrium with variation at the HLA-DRB1-DQB1 loci into account, 11 SNPs -SNP07 (novel), SNP06 (novel), SNP05 (novel), rs2071652, rs2857767, SNP04 (novel), rs2535260, rs2535261, SNP03 (novel), rs2252734, rs29228 and the MOG51 microsatellite showed evidence of association at a significance level of 5% with MS (data not shown). We then evaluated MS association of these markers in an additional 158 MS families (total n = 397 families). In this larger sample set, 8 markers showed some evidence of association at 5% significance level (Figure 2) with one haplotype defined by variants at rs2071652 (allele C), rs2857767 (allele C) and rs2535260 (allele C) showing the strongest disease association (P = 3.29 × 10-6). This level of significance was slightly lower than that shown by association of HLA-DRB1 (P = 7.9 × 10-7) and -DQB1 in the same sample set (P = 1.3 × 10-7). The transmission data of these variants individually and in the context of the three locus DRB1-DQB1-MOG haplotypes are shown in the online appendix along with the D' values among MOG variants and the various DRB1-DQB1 haplotypes. The eight disease associated variants are essentially part of two main LD blocks formed, on the one hand by rs2071652, rs2857767 and rs2535260 and on the other hand by the other five variants. From these data it is evident that while SNP07, SNP06; SNP04, SNP03 are in positive LD only with the DRB1*0301-DQB1*0201 haplotype, SNPs rs2071652, rs2857767 and rs2535260 occur frequently (albeit at different degrees of positive LD from strong to modest) and show a significant association with more than one HLA-DRB1-DQB1 haplotype. The co-occurrence of these MOG variants with different DRB1-DQB1 haplotypes could just be a chance event or could be due to any of the various forces which impact LD . Importantly, we also assessed the disease association of MOG alleles which individually showed significant evidence of association with MS together with the other DRB1-DQB1 haplotypes grouped (that is, all DRB1-DQB1 haplotypes which did not show any evidence of significant association with MS). These extended DRB1-DQB1-MOG haplotypes showed a consistent negative transmission from parents to affected patients (see Additional file 1). This suggested that MOG disease associations were most likely driven by DRB1-DQB1 variation.
Next, we wanted to directly assess whether these MOG associations were independent from DRB1-DQB1 associations. To this end, we analysed the overall effect of allelic and haplotype variation at the MOG gene using the conditional ETDT (CETDT) . Using this conditional analysis, which takes into account linkage disequilibrium between overall MOG variation and overall variation at the HLA-DQB1 and -DRB1 loci, MOG associations were no longer significant. Conversely, the overall associations of DRB1-DQB1 haplotypes, taking into account the most associated MOG haplotype, gave a P value of 9.7 × 10-3.
We also evaluated the relative transmission of the MOG alleles and haplotypes according to which specific DRB1-DQA1-DQB1 haplotypes they were on (Figure 3). In order to carry out these analyses we used the haplotype method (HM)-TDT (Methods) . Using this test we initially evaluated association of specific variants of MOG, taking into account linkage disequilibrium with the HLA-DRB1-DQB1 haplotypes for which there is prior evidence of association with MS in the Sardinian population (i.e. DRB1*1303-DQB1*0301, DRB1*0405-DQB1*0301, DRB1*0301-DQB1*0201, DRB1*1501-DQB1*0602 and DRB1*0405-DQB1*0302 haplotypes) . Again, allelic variation at the MOG locus did not significantly affect the transmission of these predisposing DRB1-DQB1 haplotypes (Figure 3). Also the MOG haplotype defined by rs2071652, rs2857767 and rs2535260 did not produce significant evidence of primary association when analysed with the HM-TDT conditional on DRB1-DQB1 positively associated haplotypes (data not shown). Moreover, we also conducted the HM-TDT to assess whether these MOG variants were influencing association of the non-predisposing DRB1-DQB1 haplotypes grouped. No significant heterogeneity in the association of DRB1-DQB1 haplotype conditional on MOG was seen, again suggesting that the initial observations of association with MS risk were likely to be only a consequence of LD with the DR-DQ mediated MS risk. This is shown in Figure 3, Graph F.
Considering our results globally, they provide evidence that haplotypes identical at the DRB1-DQB1 loci but different at the MOG locus do not have different associations with MS. Thus the significant association of MOG alleles and haplotypes appear to be secondary to hitchhiking effects due to variation at the DRB1-DQB1 loci. In interpreting these results it should be considered that the HLA region as a whole confers only a relatively modest contribution to the familial clustering and to the individual inherited risk of MS. Thus, the statistical power required to rule out a rare disease variant not detected in our resequencing, or even more dramatically, one that only slightly modulates the disease risk conferred by HLA-DRB1-DQB1 variants would require unrealistic sample sizes which exceed those used in this study. Still, our results and conclusions are consistent with those obtained in other populations , although there is still the suggestion that other non-DR-DQ genes may convey an independent risk to MS [41, 42, 12, 24].
The genetic results presented here have relevant mechanistic biological implications because they suggest that the protein product of the candidate gene MOG is unlikely to represent the primary autoantigen in human MS and thus cannot explain the target organ specificity observed in this disease. Contrast these data with the primary and consistent association detected with variation in the insulin gene (INS) region in type 1 diabetes and the epistatic deviations from a multiplicative model detected in the joint analysis of HLA class II and INS .
From our results we can conclude that variation at the MOG gene does not provide a significant independent contribution to the inherited risk for MS in the Sardinian population. The data suggests that the protein product of MOG is unlikely to represent the primary autoantigen in human MS.
In the study we considered a sample set of 397 independent MS families, which included 391 trios families (one affected child and both healthy parents), 4 multiplex families (more than one affected sibling and their healthy parents) and 2 vertical families (one parent and one affected offspring), all being at least third generation Sardinian. This data set of 397 MS families is a proper subset of the 490 Sardinian MS families described previously . All patients met the criteria of definite MS . The average age of onset for MS was 28.9 years (range 6 – 51 years), and the female/male ratio was 2.3. The disease course was relapsing-remitting in 65%, secondary progressive in 27% and primary progressive in 8% of patients. All families had previously been analysed at the DRB1-DQA1-DQB1 loci .
To examine the feasibility of verifying whether variants in the MOG gene, or in strong LD with it could explain the association peak that was previously observed in Sardinia , we focused on the parent-affected child trios that we had already studied during the initial examination of the role of the HLA in MS susceptibility. We postulated that our search would be for a fairly common variant (in line with the frequency of the predisposing DRB1-DQB1 haplotypes) with a maximal relative risk of 1.9 and using a disease prevalence of 0.0014 (corresponding to the prevalence in Sardinia). The strength of LD between the DRB1-DQB1 locus and MOG was estimated at a mean value of r2 = 0.3 as previously shown . Purcell and Sham power calculations for our 391 trios gave for an alpha value of 0.10 a power of 0.83 . It should be noted that for the conditional analyses, in the case of conditioning on the absence of high-risk HLA haplotypes, the trio size is 121 and the corresponding power is reduced to 0.42 for the same alpha criteria.
MOG polymorphism analysis
We resequenced the MOG gene in 21 parents of MS patients heterozygous for the 226 bp allele at the MOG51 microsatellite locus, which was found to be significantly associated with MS . This selection of informative parents increases the chance of detecting relevant polymorphisms. 37 overlapping fragments covering the whole gene (8 exons and 7 introns) and -384 and 545 bp 5' and 3' from the transcriptional unit were amplified and capillary sequenced (Megabace 1000). Overall, 16,561 bp were sequenced. Primers were designed according to the published human genomic DNA sequence  and are listed in Table 1). PCR was carried out in a final volume of 25 μl using a specific combination of primers and 5 μl of 4 ng/μl genomic DNA. After an initial denaturation of 4 minutes at 95°C, there were 37 cycles at 95°C for 30 seconds, at the specific TM primer for 30 s, 72°C for 30 s and a final extension at 72°C for 7 minutes.
A set of 40 informative SNPs was selected and dot-blot genotyped in the family sample set, using the set of primers listed in Table 1 and the SSO probes reported in Table 2. Selection of these SNPs was based on the patterns of LD (in order to provide good map coverage and avoid genotyping and testing genetically redundant variants), on the position in the gene (in order to maximise the chance of testing functionally relevant variants) and on the minor allele frequency (MAF) (in order to exclude particularly rare variants for which we did not have sufficient statistical power in the subsequent association stage).
Single-point disease association of individual MOG variants as well as HLA-DRB1, DQB1 alleles was evaluated using ETDT . This test takes into account the transmission or non-transmission of alleles of a locus relative to those present on the other parental chromosome. In the case of multiallelic markers, the ETDT takes multiple alleles into account and obtains a global P value indicative of the degree of significance of the association with the disease at each individual locus. In order to distinguish primary associations from those due to linkage disequilibrium at the various assessed loci, we used a variant of the ETDT, called CETDT . This test allows us to analyse the overall effect of one locus, in our case MOG, taking the overall association of variation at other loci into account; in this study DRB1 and DQB1. In order to study the transmission of specific MOG alleles, conditioned on alleles or haplotypes at the DRB1-DQB1 disease loci, we used a TDT variant of the haplotype method ; the HM-TDT . This method tests the null hypothesis of equality of transmission of haplotypes identical at one variant or combination of variants, but different at another closely-linked test locus . If there is heterogeneity in the transmission of two marker haplotypes identical at a predisposing marker (variant A) but different at a putative predisposing marker at another site (variant B), then this is evidence that variant A does not entirely explain disease predisposition and that variant B itself, or another marker in linkage disequilibrium with variant B, is influencing the transmission of variant A and, thus, disease susceptibility .
LD between individual SNPs at the MOG gene was assessed computing the D' and r2 values. Both measures are built on the basic pairwise-disequilibrium coefficient, D, described by Lewontin , representing the difference between the probability of observing two marker alleles on the same haplotype and observing them independently in the population. The D and normalised disequilibrium D', or |D/Dmax|, values were computed on parental haplotypes using the software package . A value of D' equal to 0.0 implies complete independence, whereas 1.0 implies that all copies of the rare allele occur exclusively with one of the two possible alleles at the other assessed marker. R2, computed using a programme written by Jason Cooper, assesses not only LD between two markers but also the degree of matching between all allelic frequencies of both SNPs , thus encompassing the effects of both D and all marker allele frequencies. As for D', an r2 value of 0.0 also implies independence, but r2 = 1.0 occurs only when the marker loci have identical allele frequencies and every occurrence of an allele at each of the markers perfectly predicts the allele at the other locus (i.e. in the presence of genetically indistinguishable makers). The two LD measures have different features: D' levels are useful to define the extent of LD, for instance to establish the presence of the so-called LD blocks along chromosome intervals. Instead, r2 are useful to establish marker informativity and the degree of genetic redundancy between selected SNPs and thus determine map coverage .
For all analyses performed in this manuscript, haplotypes were established following the co-segregation of alleles within families and using computer programmes written by F. Dudbridge and B. Koeleman. Only haplotypes showing definitive parental genotype data are considered in the analyses shown in this paper. In those families with more than one affected patient, only the probands were evaluated.
extended transmission disequilibrium test
haplotype method transmission disequilibrium test
- MOG :
Myelin Oligodendrocyte Glycoprotein
Withaker JM: Clinical feature of multiple sclerosis. Multiple Sclerosis: Clinical and Pathogenetic Basis. Edited by: In Raine CMFHTW. 1997, London, 3-19.
Granieri E, Casetta I, Govoni V, Tola MR, Marchi D, Murgia SB, Ticca A, Pugliatti M, Murgia B, Rosati G: The increasing incidence and prevalence of MS in a Sardinian province. Neurology. 2000, 55 (6): 842-848.
Pugliatti M, Riise T, Sotgiu MA, Sotgiu S, Satta WM, Mannu L, Sanna G, Rosati G: Increasing incidence of multiple sclerosis in the province of Sassari, northern Sardinia. Neuroepidemiology. 2005, 25 (3): 129-134. 10.1159/000086677.
Pugliatti M, Sotgiu S, Solinas G, Castiglia P, Pirastru MI, Murgia B, Mannu L, Sanna G, Rosati G: Multiple sclerosis epidemiology in Sardinia: evidence for a true increasing risk. Acta Neurol Scand. 2001, 103 (1): 20-26. 10.1034/j.1600-0404.2001.00207.x.
Ebers GC, Sadovnick AD, Risch NJ: A genetic basis for familial aggregation in multiple sclerosis. Canadian Collaborative Study Group. Nature. 1995, 377 (6545): 150-151. 10.1038/377150a0.
Compston A, Coles A: Multiple sclerosis. Lancet. 2002, 359 (9313): 1221-1231. 10.1016/S0140-6736(02)08220-X.
Barcellos LF, Oksenberg JR, Begovich AB, Martin ER, Schmidt S, Vittinghoff E, Goodin DS, Pelletier D, Lincoln RR, Bucher P, Swerdlin A, Pericak-Vance MA, Haines JL, Hauser SL: HLA-DR2 dose effect on susceptibility to multiple sclerosis and influence on disease course. Am J Hum Genet. 2003, 72 (3): 710-716. 10.1086/367781.
Stewart GJ, Teutsch SM, Castle M, Heard RN, Bennetts BH: HLA-DR, -DQA1 and -DQB1 associations in Australian multiple sclerosis patients. Eur J Immunogenet. 1997, 24 (2): 81-92. 10.1046/j.1365-2370.1997.00252.x.
Laaksonen M, Pastinen T, Sjoroos M, Kuokkanen S, Ruutiainen J, Sumelahti ML, Reijonen H, Salonen R, Wikstrom J, Panelius M, Partanen J, Tienari PJ, Ilonen J: HLA class II associated risk and protection against multiple sclerosis-a Finnish family study. J Neuroimmunol. 2002, 122 (1-2): 140-145. 10.1016/S0165-5728(01)00456-8.
Chataway J, Feakes R, Coraddu F, Gray J, Deans J, Fraser M, Robertson N, Broadley S, Jones H, Clayton D, Goodfellow P, Sawcer S, Compston A: The genetics of multiple sclerosis: principles, background and updated results of the United Kingdom systematic genome screen. Brain. 1998, 121 ( Pt 10): 1869-1887. 10.1093/brain/121.10.1869.
Oksenberg JR, Barcellos LF, Cree BA, Baranzini SE, Bugawan TL, Khan O, Lincoln RR, Swerdlin A, Mignot E, Lin L, Goodin D, Erlich HA, Schmidt S, Thomson G, Reich DE, Pericak-Vance MA, Haines JL, Hauser SL: Mapping multiple sclerosis susceptibility to the HLA-DR locus in African Americans. Am J Hum Genet. 2004, 74 (1): 160-167. 10.1086/380997.
Lincoln MR, Montpetit A, Cader MZ, Saarela J, Dyment DA, Tiislar M, Ferretti V, Tienari PJ, Sadovnick AD, Peltonen L, Ebers GC, Hudson TJ: A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis. Nat Genet. 2005, 37 (10): 1108-1112. 10.1038/ng1647.
Marrosu MG, Murru R, Murru MR, Costa G, Zavattari P, Whalen M, Cocco E, Mancosu C, Schirru L, Solla E, Fadda E, Melis C, Porru I, Rolesu M, Cucca F: Dissection of the HLA association with multiple sclerosis in the founder isolated population of Sardinia. Hum Mol Genet. 2001, 10 (25): 2907-2916. 10.1093/hmg/10.25.2907.
Masterman T, Ligers A, Olsson T, Andersson M, Olerup O, Hillert J: HLA-DR15 is associated with lower age at onset in multiple sclerosis. Ann Neurol. 2000, 48 (2): 211-219. 10.1002/1531-8249(200008)48:2<211::AID-ANA11>3.0.CO;2-R.
C.Alvarado-de la Barrera, J. Zùniga-Ramos, J.A. Ruìz-Morales, B. Estanol, J. Granados, L. Llorente: HLA class II genotypes in Mexican Mestizos with familial and nonfamilial multiple sclerosis. Neurology. 2000, 55: 1897-1900.
Kwon OJ, Karni A, Israel S, Brautbar C, Amar A, Meiner Z, Abramsky O, Karussis D: HLA class II susceptibility to multiple sclerosis among Ashkenazi and non-Ashkenazi Jews. Arch Neurol. 1999, 56 (5): 555-560. 10.1001/archneur.56.5.555.
Saruhan-Direskeneli G, Esin S, Baykan-Kurt B, Ornek I, Vaughan R, Eraksoy M: HLA-DR and -DQ associations with multiple sclerosis in Turkey. Hum Immunol. 1997, 55 (1): 59-65. 10.1016/S0198-8859(97)00086-4.
Luomala M, Elovaara I, Ukkonen M, Koivula T, Lehtimaki T: The combination of HLA-DR1 and HLA-DR53 protects against MS. Neurology. 2001, 56 (3): 383-385.
Coraddu F, Reyes-Yanez MP, Parra A, Gray J, Smith SI, Taylor CJ, Compston DA: HLA associations with multiple sclerosis in the Canary Islands. J Neuroimmunol. 1998, 87 (1-2): 130-135. 10.1016/S0165-5728(98)00074-5.
Caballero A, Alves-Leon S, Papais-Alvarenga R, Fernandez O, Navarro G, Alonso A: DQB1*0602 confers genetic susceptibility to multiple sclerosis in Afro-Brazilians. Tissue Antigens. 1999, 54 (5): 524-526. 10.1034/j.1399-0039.1999.540511.x.
Fogdell-Hahn A, Ligers A, Gronning M, Hillert J, Olerup O: Multiple sclerosis: a modifying influence of HLA class I genes in an HLA class II associated autoimmune disease. Tissue Antigens. 2000, 55 (2): 140-148. 10.1034/j.1399-0039.2000.550205.x.
Ebers GC, Kukay K, Bulman DE, Sadovnick AD, Rice G, Anderson C, Armstrong H, Cousin K, Bell RB, Hader W, Paty DW, Hashimoto S, Oger J, Duquette P, Warren S, Gray T, O'Connor P, Nath A, Auty A, Metz L, Francis G, Paulseth JE, Murray TJ, Pryse-Phillips W, Nelson R, Freedman M, Brunet D, Bouchard JP, Hinds D, Risch N: A full genome search in multiple sclerosis. Nat Genet. 1996, 13 (4): 472-476. 10.1038/ng0896-472.
Bitti PP, Murgia BS, Ticca A, Ferrai R, Musu L, Piras ML, Puledda E, Campo S, Durando S, Montomoli C, Clayton DG, Mander AP, Bernardinelli L: Association between the ancestral haplotype HLA A30B18DR3 and multiple sclerosis in central Sardinia. Genet Epidemiol. 2001, 20 (2): 271-283. 10.1002/1098-2272(200102)20:2<271::AID-GEPI9>3.0.CO;2-L.
Rubio JP, Bahlo M, Stankovich J, Burfoot RK, Johnson LJ, Huxtable S, Butzkueven H, Lin L, Taylor BV, Speed TP, Kilpatrick TJ, Mignot E, Foote SJ: Analysis of extended HLA haplotypes in multiple sclerosis and narcolepsy families confirms a predisposing effect for the class I region in Tasmanian MS patients. Immunogenetics. 2007, 59 (3): 177-186. 10.1007/s00251-006-0183-5.
Stefferl A, Brehm U, Storch M, Lambracht-Washington D, Bourquin C, Wonigeit K, Lassmann H, Linington C: Myelin oligodendrocyte glycoprotein induces experimental autoimmune encephalomyelitis in the "resistant" Brown Norway rat: disease susceptibility is determined by MHC and MHC-linked effects on the B cell response. J Immunol. 1999, 163 (1): 40-49.
Mendel I, Gur H, Kerlero de Rosbo N, Ben-Nun A: Experimental autoimmune encephalomyelitis induced in B6.C-H-2bm12 mice by myelin oligodendrocyte glycoprotein: effect of MHC class II mutation on immunodominant epitope selection and fine epitope specificity of encephalitogenic T cells. J Neuroimmunol. 1999, 96 (1): 9-20. 10.1016/S0165-5728(99)00009-0.
Kerlero de Rosbo N, Hoffmann M, Mendel I, Yust I, Kaye J, Bakimer R, Flechter S, Abramsky O, Milo R, Karni A, A BN: Predominance of the autoimmune response to myelin oligodendrocyte glycoprotein (MOG) in multiple sclerosis: reactivity to the extracellular domain of MOG is directed against three main regions. Eur J Immunol. 1997, 27: 3059-3069. 10.1002/eji.1830271144.
Wallstrom E, Khademi M, Andersson M, Weissert R, Linington C, Olsson T: Increased reactivity to myelin oligodendrocyte glycoprotein peptides and epitope mapping in HLA DR2(15)+ multiple sclerosis. Eur J Immunol. 1998, 28 (10): 3329-3335. 10.1002/(SICI)1521-4141(199810)28:10<3329::AID-IMMU3329>3.0.CO;2-B.
Raine CS, Cannella B, Hauser SL, Genain CP: Demyelination in primate autoimmune encephalomyelitis and acute multiple sclerosis lesions: a case for antigen-specific antibody mediation. Ann Neurol. 1999, 46 (2): 144-160. 10.1002/1531-8249(199908)46:2<144::AID-ANA3>3.0.CO;2-K.
O'Connor. KC, McLaughlin KA, De Jager PL, Chitnis T, Bettelli E, Xu C, Robinson WH, Cherry SV, Bar-Or A, Banwell B, Fukaura H, Fukazawa T, Tenembaum S, Wong SJ, Tavakoli NP, Idrissova Z, Viglietta V, Rostasy K, Pohl D, Dale RC, Freedman M, Steinman L, Buckle GJ, Kuchroo VK, Hafler DA, Wucherpfennig KW: Self-antigen tetramers discriminate between myelin autoantibodies to native or denatured protein. Nat Med. 2007, 13 (2): 211-217. 10.1038/nm1488.
Kuhle J, Pohl C, Mehling M, Edan G, Freedman MS, Hartung HP, Polman CH, Miller DH, Montalban X, Barkhof F, Bauer L, Dahms S, Lindberg R, Kappos L, Sandbrink R: Lack of association between antimyelin antibodies and progression to multiple sclerosis. N Engl J Med. 2007, 356 (4): 371-378. 10.1056/NEJMoa063602.
Seboun E, Oksenberg JR, Rombos A, Usuku K, Goodkin DE, Lincoln RR, Wong M, Pham-Dinh D, Boesplug-Tanguy O, Carsique R, Fitoussi R, Gartioux C, Reyes C, Ribierre F, Faure S, Fizames C, Gyapay G, Weissenbach J, Dautigny A, Rimmler JB, Garcia ME, Pericak-Vance MA, Haines JL, Hauser SL: Linkage analysis of candidate myelin genes in familial multiple sclerosis. Neurogenetics. 1999, 2 (3): 155-162. 10.1007/s100480050076.
Gomez-Lira M, Moretto G, Bonamini D, Benedetti MD, Pignatti PF, Rizzuto N, Salviati A: Myelin oligodendrocyte glycoprotein polymorphisms and multiple sclerosis. J Neuroimmunol. 2002, 133 (1-2): 241-243. 10.1016/S0165-5728(02)00361-2.
Roth MP, Malfroy L, Offer C, Sevin J, Enault G, Borot N, Pontarotti P, Coppin H: The human myelin oligodendrocyte glycoprotein (MOG) gene: complete nucleotide sequence and structural characterization. Genomics. 1995, 28 (2): 241-250. 10.1006/geno.1995.1137.
Ohlenbusch A, Pohl D, Hanefeld F: Myelin oligodendrocyte gene polymorphisms and childhood multiple sclerosis. Pediatr Res. 2002, 52 (2): 175-179.
National Center for Biotechnology Information, http://www.ncbi.nlm.nih.gov/: .
Zavattari P, Deidda E, Whalen M, Lampis R, Mulargia A, Loddo M, Eaves I, Mastio G, Todd JA, Cucca F: Major factors influencing linkage disequilibrium by analysis of different chromosome regions in distinct populations: demography, chromosome recombination frequency and selection. Hum Mol Genet. 2000, 9 (20): 2947-2957. 10.1093/hmg/9.20.2947.
Koeleman BP, Dudbridge F, Cordell HJ, Todd JA: Adaptation of the extended transmission/disequilibrium test to distinguish disease associations of multiple loci: the Conditional Extended Transmission/Disequilibrium Test. Ann Hum Genet. 2000, 64 (Pt 3): 207-213. 10.1046/j.1469-1809.2000.6430207.x.
Cucca F, Dudbridge F, Loddo M, Mulargia AP, Lampis R, Angius E, De Virgiliis S, Koeleman BP, Bain SC, Barnett AH, Gilchrist F, Cordell H, Welsh K, Todd JA: The HLA-DPB1--associated component of the IDDM1 and its relationship to the major loci HLA-DQB1, -DQA1, and -DRB1. Diabetes. 2001, 50 (5): 1200-1205. 10.2337/diabetes.50.5.1200.
Rubio JP, Bahlo M, Tubridy N, Stankovich J, Burfoot R, Butzkueven H, Chapman C, Johnson L, Marriott M, Mraz G, Tait B, Wilkinson C, Taylor B, Speed TP, Foote SJ, Kilpatrick TJ: Extended haplotype analysis in the HLA complex reveals an increased frequency of the HFE-C282Y mutation in individuals with multiple sclerosis. Hum Genet. 2004, 114 (6): 573-580. 10.1007/s00439-004-1095-9.
Yeo TW, De Jager PL, Gregory SG, Barcellos LF, Walton A, Goris A, Fenoglio C, Ban M, Taylor CJ, Goodman RS, Walsh E, Wolfish CS, Horton R, Traherne J, Beck S, Trowsdale J, Caillier SJ, Ivinson AJ, Green T, Pobywajlo S, Lander ES, Pericak-Vance MA, Haines JL, Daly MJ, Oksenberg JR, Hauser SL, Compston A, Hafler DA, Rioux JD, Sawcer S: A second major histocompatibility complex susceptibility locus for multiple sclerosis. Ann Neurol. 2007, 61 (3): 228-236. 10.1002/ana.21063.
Rubio JP, Bahlo M, Butzkueven H, van Der Mei IA, Sale MM, Dickinson JL, Groom P, Johnson LJ, Simmons RD, Tait B, Varney M, Taylor B, Dwyer T, Williamson R, Gough NM, Kilpatrick TJ, Speed TP, Foote SJ: Genetic dissection of the human leukocyte antigen region by use of haplotypes of Tasmanians with multiple sclerosis. Am J Hum Genet. 2002, 70 (5): 1125-1137. 10.1086/339932.
Motzo C, Contu D, Cordell HJ, Lampis R, Congia M, Marrosu MG, Todd JA, Devoto M, Cucca F: Heterogeneity in the magnitude of the insulin gene effect on HLA risk in type 1 diabetes. Diabetes. 2004, 53 (12): 3286-3291. 10.2337/diabetes.53.12.3286.
Purcell S, Cherny SS, Sham PC: Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19 (1): 149-150. 10.1093/bioinformatics/19.1.149.
Sham PC, Curtis D: An extended transmission/disequilibrium test (TDT) for multi-allele marker loci. Ann Hum Genet. 1995, 59 ( Pt 3): 323-336.
Valdes AM, Thomson G: Detecting disease-predisposing variants: the haplotype method. Am J Hum Genet. 1997, 60 (3): 703-716.
Lewontin RC: On measures of gametic disequilibrium. Genetics. 1988, 120 (3): 849-852.
A software for population genetics data analysis, http://lgb.unige.ch/arlequin/: .
Zondervan KT, Cardon LR: The complex interplay among factors that influence allelic association. Nat Rev Genet. 2004, 5 (2): 89-100. 10.1038/nrg1270.
We wish to thank Rebecca Lewis for editorial assistance in the preparation of the manuscript and Jason Cooper for providing us with a programme to compute r2 square values. This work was supported by FISM – Fondazione Italiana Sclerosi Multipla – Cod. 2004/R/5.
The author(s) declare that they have no competing interests.
MGM1participated in the study design, discussed results with FC3and reviewed the manuscript.RM1checked all data files for genotyping problems, performed statistical analysis on the data and designed all figures and tables. GC1set up and supervised the different protocols and methods performed.MCM1, MR1, LS1and ES1performed amplification and re-sequencing of all exons and introns. SC1completed all DNA extraction from peripheral blood. MAS1optimised DNA concentration for all samples and prepared all aliquots and dilutions. EC1organised family selection and maintained contacts with them for sample collection. MBW3assisted in manuscript revision and statistical analysis. MP2, SS2and GR2provided useful suggestions and collected DNA for part of the patients studied together with the related clinical data. FC3supervised the study, analysed the results and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: TDT results for Alleles of the MOG haplotypes treated separately or with DRB1-DQB1. Transmission test for linkage disequilibrium was conducted with each of the alleles for the MOG haplotypes. The column on the left shows the results for each allele singly, while the column on the right shows the transmission results for HLA DRB1-DQB1 when the MOG allele is fixed. (DOC 128 KB)
About this article
Cite this article
Marrosu, M.G., Murru, R., Costa, G. et al. Variation of the Myelin Oligodendrocyte Glycoprotein gene is not primarily associated with multiple sclerosis in the Sardinian population. BMC Genet 8, 25 (2007). https://doi.org/10.1186/1471-2156-8-25