PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Genes Immun. Author manuscript; available in PMC 2010 April 1.
Published in final edited form as:
PMCID: PMC2765552
NIHMSID: NIHMS118482

Follow-up examination of linkage and association to chromosome 1q43 in multiple sclerosis

Abstract

Multiple sclerosis is a debilitating neuroimmunological and neurodegenerative disease affecting more than 400,000 individuals in the United States. Population and family-based studies have suggested that there is a strong genetic component. Numerous genomic linkage screens have identified regions of interest for MS loci. Our own second-generation genome-wide linkage study identified a handful of non-MHC regions with suggestive linkage. Several of these regions were further examined using single-nucleotide polymorphisms (SNPs) with average spacing between SNPs of approximately 1.0 Mb in a dataset of 173 multiplex families. The results of that study provided further evidence for the involvement of the chromosome 1q43 region. This region is of particular interest given linkage evidence in studies of other autoimmune and inflammatory diseases including rheumatoid arthritis and systemic lupus erythematosus. In this follow-up study, we saturated the region with ~700 SNPs (average spacing of 10kb per SNP) in search of disease associated variation within this region. We found preliminary evidence to suggest that common variation within the RGS7 locus may be involved in disease susceptibility.

Keywords: multiple sclerosis, linkage, association, 1q43, RGS7

Introduction

Multiple sclerosis (MS [MIM 126200]) is a neurodegenerative autoimmune disease characterized by demyelination within the central nervous system (CNS. Demyelination and the resulting formation of scar tissue in the CNS impair the saltatory conduction along axons that is necessary for normal functioning of nerve impulses. Though little is known about the underlying etiology of the disease, MS is a heterogeneous disorder with several characteristics common to autoimmune disorders —including polygenic inheritance, evidence of environmental exposure, and partial susceptibility conferred by a human leukocyte antigen (HLA)-associated gene (Barcellos, et al. 2002). There is also a greater prevalence of the disease in women than in men. This inflammatory disorder results from an autoimmune response directed against CNS antigens—particularly myelin proteins. Class II major histocompatibility complex (MHC) molecules, such as HLA (human leukocyte antigen), normally function to bind and present peptide antigens to antigen-specific T cells. It is thought that the dysregulation of this process in MS results in damage to the myelin sheath, producing the pathophysiological phenotype seen in the disease. Despite the degenerative nature of the disease, an affected individual’s lifespan is rarely shortened.

The clinical heterogeneity and complex etiology of MS have been confounding factors for genetic studies of the disease. Yet despite these complexities, it is clear that genes play a vital role in the susceptibility to MS. Numerous research groups, including our own, have conducted genomic linkage screens for MS in an attempt to identify regions that harbor disease loci (Akesson, et al. 2002, Ban, et al. 2002, Broadley, et al. 2001, Coraddu, et al. 2001, Dyment, et al. 2004, Ebers, et al. 1996, Eraksoy, et al. 2003, Haines, et al. 1996, Hensiek, et al. 2003, Kenealy, et al. 2004, Kuokkanen, et al. 1997, Sawcer, et al. 1996, Sawcer, et al. 2005). Over 70 genomic regions have been investigated, revealing varying levels of support. However, the lack of replication of results from these studies has been problematic for using linkage to detect loci involved in MS.

The strongest and most consistent finding for linkage to MS is chromosome 6p21.3, the location of the MHC containing HLA. Until just recently, the MHC has been the only region clearly and consistently demonstrating linkage and association with MS (Haines, et al. 1996, Sawcer, et al. 1996, Sawcer, et al. 2005, Yaouanq, et al. 1997). However, the MHC has been estimated to account for only 10–50% of the genetic component of MS susceptibility, at least in the Caucasian populations of northern European descent (Dyment, et al. 2004, Haines, et al. 1998). It appears that the association with the HLA-DRB1*1501 allele explains this linkage signal (Barcellos, et al. 2002, Haines, et al. 1998), although this issue has been debated (Ligers, et al. 2001). The exact mechanism by which a gene or genes in the MHC increase disease risk has yet to be determined, although recent work has sought to provide the tools needed to explore this region in greater detail (de Bakker, et al. 2006, Horton, et al. 2008). Recent work by our group has identified a strongly associated single-nucleotide polymorphism (SNP) (rs6897932) in the interleukin-7 receptor α gene (IL7RA), a gene located on chromosome 5p13. This SNP introduces a coding change (T244I) that alters the ratio of soluble to bound protein (Gregory, et al. 2007). This result was further replicated in a genome-wide association study (GWAS) performed by a broader collaborative group of investigators (The International Multiple Sclerosis Genetics Consortium (IMSGC)). The GWAS also identified additional replicated associations with variations in other genes (including CD25/IL2RA and CD58) (International Multiple Sclerosis Genetics Consortium, et al. 2007). The involvement of these genes plus several others (e.g. CLEC16A) are now being tested and confirmed in multiple international studies (International Multiple Sclerosis Genetics Consortium (IMSGC). 2008, Rubio, et al. 2008, Weber, et al. 2008).

Although genome-wide association studies have now begun, dense SNP follow-up studies of narrowed linkage regions are still underway because of their potential for identifying additional disease susceptibility genes. The American-French Multiple Sclerosis Genetics Group published one of the largest MS genomic linkage screens conducted (Kenealy, et al. 2004). Follow-up of candidate regions identified in this study entailed genotyping SNPs at approximately 1.0 Mb intervals flanking 10 Mb on each side of peak screen SNPs in an expanded U.S. dataset (Kenealy, et al. 2006). The 1q43 region not only continued to demonstrate evidence for linkage in the follow-up analysis (non-parametric multipoint LOD = 2.99), but also generated a narrowed linkage interval when using the ordered subset analysis (OSA) method to condition on the other loci (Ghosh, et al. 2000, Hauser and Boehnke. 1998, Hauser, et al. 2001, Kenealy, et al. 2006). This narrowed interval is ~7.0 Mb for a LOD score cut-off of 1.8 (corresponding to a –2.0 LOD score confidence interval from the 3.8 LOD score HLA-based peak) (see Figure 1). Chromosome 1q43 was selected for further detailed investigation based on several additional lines of evidence. SNPs in the 1q43 region have also demonstrated suggestive linkage and/or association in several other MS screens conducted in a variety of study populations (Ban, et al. 2002, Broadley, et al. 2001, Goedde, et al. 2002, Laaksonen, et al. 2003, Sawcer, et al. 2002). None of these earlier studies, which lack the power and marker density of the current generation of genome-wide association studies, demonstrate significant evidence for either linkage or association. However, they do provide support from independent sample populations of an effect within a common region of interest. Another compelling piece of evidence for 1q is linkage to this region in studies for the autoimmune diseases rheumatoid arthritis and systemic lupus erythematosus, suggesting the potential presence of a gene common to autoimmune processes in this region (Gaffney, et al. 1998, Gaffney, et al. 2000, Jawaheer, et al. 2001, Jawaheer, et al. 2003, MacKay, et al. 2002, Shai, et al. 1999). Following these lines of evidence, we set out to test whether common variation within this region is involved in susceptibility to develop multiple sclerosis.

Figure 1
Region of Interest on Chromosome 1q

Results

After quality control procedures, 578 SNPs were included in our initial analysis. Using the Tagger implementation in the Haploview program we selected 268/578 SNPs, with r2<0.2, for inclusion in a multipoint linkage analysis of the narrow interval. These results continue to demonstrate considerable evidence for linkage to this region (peak Lod Score 3.42 at ~238Mb) (Figure 2a). As was analyzed in the initial follow-up (Kenealy, et al. 2006), we examined this narrowed interval using the ordered subset analysis (OSA) method to condition on the same loci (namely HLA-DRB1*1501 and Chr2) in an attempt to further resolve this effect (Figure 2a). The goal of this analysis was to determine if other previously suggested linkage regions, within our initial dataset, have a significant effect on the 1q43 region such that we might be able to identify specific subsets of our sample that may account for the effect within this region. Neither of these conditional analyses demonstrated significant differences in lod scores from the overall multipoint analysis, suggesting that neither of these two loci can be used to predict families which will demonstrate strong evidence for a genetic effect within the 1q43 region.

Figure 2Figure 2
Multipoint Linkage Analyses across 1q43

Finally, as recent genome-wide association studies have found genetic overlap (by way of identifying common loci, which are definitively associated with multiple autoimmune-related disorders) within autoimmune disorders, it is feasible that a gene (or genes) within this region of chromosome 1q43 is affecting multiple autoimmune phenotypes (Maier and Hafler. 2008). Given the evidence for linkage in other autoimmune diseases, we split our dataset into a) families with evidence of additional autoimmune disease (n=100); b) families with no evidence of additional autoimmune diseases (n=53); note there were 20 families for which a decision could not be made given our available information. While the lod scores are generally higher for the families with evidence for multiple autoimmune diseases, this difference is not significant (Figure 2b).

Subsequent to the multipoint linkage analysis, we conducted tests of association across the 578 SNPs. Twelve SNPs demonstrated nominal p-values (p ≤0.05) (Table 1). We then chose an additional set of 51 SNPs to be genotyped (using Taqman) in our extended dataset of 831 multiplex and trio families. These SNPs provide a higher density across regions demonstrating either evidence for linkage (based on 2pt lod scores) or association. Analysis of these additional SNPs pointed to a region of interest within the RGS7 gene. Two adjacent SNPs within the follow-up demonstrated both the highest 2pt lod score (rs4660010, lod=3.24) and the most significant association (rs261809, PDT p = 0.01). While none of these results survive any sort of multiple comparisons correction, they suggest the need for further examination. We then examined the data from the recent IMSGC genome-wide association study to determine if there was any evidence for association across this region (International Multiple Sclerosis Genetics Consortium, et al. 2007). Several of these SNPs demonstrate moderately significant p-values within this same region (Table 2). It should be noted that 301 trios in our follow-up dataset (n= 831 families) overlap with the 931 trios within the GWAS. While the p-values listed in Table 2 represent independent results, the single exception is rs1380304 for which data is represented for our follow-up, the non-overlapping families in the GWAS, and the GWAS result itself. These preliminary findings further suggest that there may be a common variant within this region that is involved in MS.

Table 1
Association Results
Table 2
Follow-up and GWAS Association Results within 50kb region of RGS7

Discussion

Despite the success of previous studies in narrowing the linkage peak on chromosome 1q43, the search for a susceptibility locus in this region remains a challenging task. Even in this narrowed interval there are over 20 known or predicted genes currently listed in the public databases, several of which could serve as candidate genes due to their proposed function (e.g. involvement in autoimmunity, viral susceptibility, oxidative stress/mitochondrial function, or neuronal processes) or due to their tissue expression patterns (Table 3). In addition, many of the remaining genes in this region have undergone little or no functional characterization, likely resulting in failure of these genes to be selected for investigation in a candidate gene study despite their potential involvement in disease. These problems are not unique to studies of MS—studies of many complex diseases are hindered by flawed assumptions and lack of information concerning genes in a given candidate region. New methods that address these confounding factors are warranted. With this in mind we devised (and previously described) a process that prioritizes annotated SNPs for genotyping studies based on their location within Multi-species Conserved Sequences (MCSs) (Margulies, et al. 2003) and used this process to select SNPs across this ~7.0 Mb region of interest (McCauley, et al. 2007).

Table 3
Potential Functional Candidate Genes in the ~7 Mb region on chromosome 1q43

The linkage signal in this region remains strong, although our results indicate that this appears to be independent of HLA status, a putative chromosome 2 locus, or the presence of multiple autoimmune disorders within a family. One possible explanation is that there may be multiple rare variants present in this region, leading to the consistent detection of linkage, but lending difficulties to detect strong association to a common variant.

Interestingly, the previous reports demonstrating linkage in the MS, rheumatoid arthritis, and systemic lupus erythematosus screens are within a few Mb of the peak SNP from the MSGG genomic screen and follow-up. However, genome-wide association studies (GWAS) across multiple autoimmune disorders have failed to convincingly identify common variation within this region predisposing to any or all of these autoimmune diseases. This is based on the current literature detailing the findings of these studies in numerous other autoimmune diseases. This does not exclude the possibility that the effects of common variation within this region are too small within these other diseases to be detected with the current available datasets. When examining the IMSGC GWAS results (International Multiple Sclerosis Genetics Consortium, et al. 2007) across this region, we find further potential evidence for association to SNPs within RGS7. Future examination of this region will undoubtedly require examination of copy-number variation (CNVs) and possible resequencing to identify a potential functional variant predisposing to MS.

Although the results of this study demonstrate some modest evidence for involvement of a common variant within the RGS7 gene to be associated with MS, additional well-powered studies need to be conducted in order to confirm these findings. RGS7 is a member of the family of proteins called regulators of G protein signaling (RGS) and it has been shown that altering RGS7 function could impair the normal dampening of the inflammatory response (Hausmann, et al. 2002). Hausmann and colleagues found that expression of RGS7 in microglia and/or invading peripheral macrophages was induced by the inflammatory response created by experimental spinal cord injury in rats. It may follow that impaired regulation of RGS7 expression may have detrimental effects on the normal inflammatory cascade and that these effects may play a role in multiple sclerosis.

Patients and Methods

Subjects and phenotypes

The dataset used for the first stage of this follow-up study consisted of 173 multiplex families with 405 total affected individuals (115 male, 290 female). Of these initial multiplex families, 91 families had previous evidence for positive linkage to 1q43 (Kenealy, et al. 2006). The second stage (or extended) follow-up dataset included 831 multiplex and trio families with 1077 total affected individuals (266 male, 811 female). All families were ascertained at the University of California at San Francisco (UCSF) according to well established diagnostic and inclusion criteria.

All protocols were approved by the appropriate Institutional Review Boards and all individuals provided informed consent before participating in the study. Positive family histories were investigated by direct contact with other family members, request for medical records, and by clinical examination, laboratory testing, or paraclincial studies (MRI scanning and evoked-response testing). Individuals were placed into one of four categories: definite MS, probable MS, possible MS, and no evidence of MS. Consistent and stringent clinical criteria were applied as described elsewhere (Goodkin, et al. 1991, Haines, et al. 1996) and all clinically definite MS cases met the Poser criteria (Poser, et al. 1983). Only definite MS cases were considered as affected in the analyses.

Molecular analysis

For the first stage of follow-up, 768 SNPs were genotyped across this region in 173 multiplex families. The Illumina BeadArray™ platform was used for rapid and accurate SNP genotyping in our dataset. All genotyping was performed by the Duke Genomics Resource Laboratory Core using the Illumina BeadArray™ platform. This system provides automated outputs that ease transfer between data generation and the PEDIGENE® database system used in our statistical analyses (Haynes, et al. 1995). Due to our own SNP density requirements and the nature of the Illumina system, we chose 768 SNPs which completed two Illumina oligo pool assays (or OPAs). A detailed description of the SNP selection process for this project has been previously described ((McCauley, et al. 2007)). In general SNPs were selected by considering informativeness, validation status, location/density, and putative function. Average spacing of the 768 SNPs in the ~ 7.0 Mb region on 1q43–1q44 region was < 10 kb— allowing for coverage of the region that is appropriate for observed patterns of linkage disequilibrium in Caucasian populations (Gabriel, et al. 2002). Most SNPs fell within intronic and intergenic areas.

We examined the Illumina panel of SNPs using standard measures of quality control, dropping 190 SNPs from downstream analysis. Of the 190 SNPs that failed: 167 SNPs were monomorphic in our dataset; 22 SNPs demonstrated poor assay performance as measured by call rate (<95%), cluster separation, or other assay failure; and 1 SNP mapped to multiple locations in the genome. For stage one of the follow-up, eleven samples were dropped (including 1 affected male, 3 affected females) due to low call rates (<95%) and an excess of Mendelian errors.

Subsequent to the initial linkage and association analysis of the 173 multiplex family dataset, more detailed association analysis was performed within the extended dataset of 831 total families. SNPs, within the Illumina panel, demonstrating evidence for association as well as additional novel SNPs (n=51) within high priority regions based on linkage and association results were selected for further genotyping and analysis using Taqman assays from Applied Biosystems (Foster City, CA). The final set of the additional 51 SNPs met the design criteria for Taqman, provided an additional level of marker density in key regions, and had relatively high minor allele frequencies.

Statistical analysis

Following completion of genotyping, quality control (QC) procedures were conducted to ensure the accuracy of all genotype data. These included SNP checks for call rate and assay reliability, as well as sample checks for call rate and Mendelian inconsistencies. Thresholds used to exclude data are discussed in the results.

After assessing the quality of both SNPs and samples, we conducted a battery of analyses to localize and identify common susceptibility variants underlying our narrowed linkage peak. Multipoint linkage analysis was performed using MERLIN (Abecasis, et al. 2002). Two-point LOD scores were calculated in FASTLINK (Cottingham, et al. 1993, Schaffer, et al. 1994), with two-point heterogeneity LOD (HLOD) scores being calculated through HOMOG (Ott. 1986). Ordered subset analyses (OSA) were performed using the FLOSS software (Browning. 2006). SNPs included in both multipoint linkage analysis and FLOSS were chosen using the implementation of Tagger within the Haploview software package (Barrett, et al. 2005, de Bakker, et al. 2005). Association analyses for this study were conducted using the PDT statistic (Martin, et al. 2000).

Acknowledgements

This work was funded through NIH grants NS051695 and NS32830 and a post-doctoral fellowship (FG 1718-A-1) to JLM from the National Multiple Sclerosis Society (NMSS). This project was performed through the use of the Vanderbilt Center for Human Genetics Research Core (CHGR) facilities (DNA Resources Core and Bioinformatics Core), as well as the Duke Center for Human Genetics Molecular Genetics Core. We also recognize the International Multiple Sclerosis Genetics Consortium (IMSGC) for the use of data generated by the efforts and support of this broader collaboration. Most notably, we thank the patients and families that contributed to this study.

References

  • Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30:97–101. [PubMed]
  • Akesson E, Oturai A, Berg J, Fredrikson S, Andersen O, Harbo HF, et al. A genome-wide screen for linkage in nordic sib-pairs with multiple sclerosis. Genes Immun. 2002;3:279–285. [PubMed]
  • Ban M, Stewart GJ, Bennetts BH, Heard R, Simmons R, Maranian M, et al. A genome screen for linkage in australian sibling-pairs with multiple sclerosis. Genes Immun. 2002;3:464–469. [PubMed]
  • Barcellos LF, Oksenberg JR, Green AJ, Bucher P, Rimmler JB, Schmidt S, et al. Genetic basis for clinical expression in multiple sclerosis. Brain. 2002;125:150–158. [PubMed]
  • Barrett JC, Fry B, Maller J, Daly MJ. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]
  • Broadley S, Sawcer S, D'Alfonso S, Hensiek A, Coraddu F, Gray J, et al. A genome screen for multiple sclerosis in italian families. Genes Immun. 2001;2:205–210. [PubMed]
  • Browning BL. FLOSS: Flexible ordered subset analysis for linkage mapping of complex traits. Bioinformatics. 2006;22:512–513. [PubMed]
  • Coraddu F, Sawcer S, D'Alfonso S, Lai M, Hensiek A, Solla E, et al. A genome screen for multiple sclerosis in sardinian multiplex families. Eur J Hum Genet. 2001;9:621–626. [PubMed]
  • Cottingham RW, Jr, Idury RM, Schaffer AA. Faster sequential genetic linkage computations. Am J Hum Genet. 1993;53:252–263. [PubMed]
  • de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, Altshuler D. Efficiency and power in genetic association studies. Nat Genet. 2005;37:1217–1223. [PubMed]
  • de Bakker PI, McVean G, Sabeti PC, Miretti MM, Green T, Marchini J, et al. A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nat Genet. 2006;38:1166–1172. [PMC free article] [PubMed]
  • Dyment DA, Sadovnick AD, Willer CJ, Armstrong H, Cader ZM, Wiltshire S, et al. An extended genome scan in 442 canadian multiple sclerosis-affected sibships: A report from the canadian collaborative study group. Hum Mol Genet. 2004;13:1005–1015. [PubMed]
  • Ebers GC, Kukay K, Bulman DE, Sadovnick AD, Rice G, Anderson C, et al. A full genome search in multiple sclerosis. Nat Genet. 1996;13:472–476. [PubMed]
  • Eraksoy M, Kurtuncu M, Akman-Demir G, Kilinc M, Gedizlioglu M, Mirza M, et al. A whole genome screen for linkage in turkish multiple sclerosis. J Neuroimmunol. 2003;143:17–24. [PubMed]
  • Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–2229. [PubMed]
  • Gaffney PM, Kearns GM, Shark KB, Ortmann WA, Selby SA, Malmgren ML, et al. A genome-wide search for susceptibility genes in human systemic lupus erythematosus sib-pair families. Proc Natl Acad Sci U S A. 1998;95:14875–14879. [PubMed]
  • Gaffney PM, Ortmann WA, Selby SA, Shark KB, Ockenden TC, Rohlf KE, et al. Genome screening in human systemic lupus erythematosus: Results from a second minnesota cohort and combined analyses of 187 sib-pair families. Am J Hum Genet. 2000;66:547–556. [PubMed]
  • Ghosh S, Watanabe RM, Valle TT, Hauser ER, Magnuson VL, Langefeld CD, et al. The finland-united states investigation of non-insulin-dependent diabetes mellitus genetics (FUSION) study. I. an autosomal genome scan for genes that predispose to type 2 diabetes. Am J Hum Genet. 2000;67:1174–1185. [PubMed]
  • Goedde R, Sawcer S, Boehringer S, Miterski B, Sindern E, Haupts M, et al. A genome screen for linkage disequilibrium in HLA-DRB1*15-positive germans with multiple sclerosis based on 4666 microsatellite markers. Hum Genet. 2002;111:270–277. [PubMed]
  • Goodkin DE, Doolittle TH, Hauser SL, Ransohoff RM, Roses AD, Rudick RA. Diagnostic criteria for multiple sclerosis research involving multiply affected families. Arch Neurol. 1991;48:805–807. [PubMed]
  • Gregory SG, Schmidt S, Seth P, Oksenberg JR, Hart J, Prokop A, et al. Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis. Nat Genet. 2007;39:1083–1091. [PubMed]
  • Haines JL, Terwedow HA, Burgess K, Pericak-Vance MA, Rimmler JB, Martin ER, et al. Linkage of the MHC to familial multiple sclerosis suggests genetic heterogeneity. the multiple sclerosis genetics group. Hum Mol Genet. 1998;7:1229–1234. [PubMed]
  • Haines JL, Ter Minassian M, Bazyk A, Gusella JF, Kim DJ, Terwedow H, et al. A complete genomic screen for multiple sclerosis underscores a role for the major histocompatability complex. the multiple sclerosis genetics group. Nat Genet. 1996;13:469–471. [PubMed]
  • Hauser ER, Boehnke M. Genetic linkage analysis of complex genetic traits by using affected sibling pairs. Biometrics. 1998;54:1238–1246. [PubMed]
  • Hauser ER, Bass MP, Martin ER, Watanabe RM, Duren WL, Boehnke M. Power of the ordered subset method for detection and localization of genes in linkage analysis of complex traits. Am J Hum Genet. 2001;69:529.
  • Hausmann ON, Hu WH, Keren-Raifman T, Witherow DS, Wang Q, Levay K, et al. Spinal cord injury induces expression of RGS7 in microglia/macrophages in rats. Eur J Neurosci. 2002;15:602–612. [PubMed]
  • Haynes C, Speer MC, Peedin M, Roses AD, Haines JL, Vance JM, et al. PEDIGENE: A comprehensive data management system to facilitate efficient and rapid disease gene mapping. Am J Hum Genet. 1995;57:A193.
  • Hensiek AE, Roxburgh R, Smilie B, Coraddu F, Akesson E, Holmans P, et al. Updated results of the united kingdom linkage-based genome screen in multiple sclerosis. J Neuroimmunol. 2003;143:25–30. [PubMed]
  • Horton R, Gibson R, Coggill P, Miretti M, Allcock RJ, Almeida J, et al. Variation analysis and gene annotation of eight MHC haplotypes: The MHC haplotype project. Immunogenetics. 2008;60:1–18. [PMC free article] [PubMed]
  • International Multiple Sclerosis Genetics Consortium. Hafler DA, Compston A, Sawcer S, Lander ES, Daly MJ, et al. Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med. 2007;357:851–862. [PubMed]
  • International Multiple Sclerosis Genetics Consortium (IMSGC) Refining genetic associations in multiple sclerosis. Lancet Neurol. 2008;7:567–569. [PMC free article] [PubMed]
  • Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Monteiro J, et al. A genomewide screen in multiplex rheumatoid arthritis families suggests genetic overlap with other autoimmune diseases. Am J Hum Genet. 2001;68:927–936. [PubMed]
  • Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Etzel C, et al. Screening the genome for rheumatoid arthritis susceptibility genes: A replication study and combined analysis of 512 multicase families. Arthritis Rheum. 2003;48:906–916. [PubMed]
  • Kenealy SJ, Herrel LA, Bradford Y, Schnetz-Boutaud N, Oksenberg JR, Hauser SL, et al. Examination of seven candidate regions for multiple sclerosis: Strong evidence of linkage to chromosome 1q44. Genes Immun. 2006 [PubMed]
  • Kenealy SJ, Babron MC, Bradford Y, Schnetz-Boutaud N, Haines JL, Rimmler JB, et al. A second-generation genomic screen for multiple sclerosis. Am J Hum Genet. 2004;75:1070–1078. [PubMed]
  • Kuokkanen S, Gschwend M, Rioux JD, Daly MJ, Terwilliger JD, Tienari PJ, et al. Genomewide scan of multiple sclerosis in finnish multiplex families. Am J Hum Genet. 1997;61:1379–1387. [PubMed]
  • Laaksonen M, Jonasdottir A, Fossdal R, Ruutiainen J, Sawcer S, Compston A, et al. A whole genome association study in finnish multiple sclerosis patients with 3669 markers. J Neuroimmunol. 2003;143:70–73. [PubMed]
  • Ligers A, Dyment DA, Willer CJ, Sadovnick AD, Ebers G, Risch N, et al. Evidence of linkage with HLA-DR in DRB1*15-negative families with multiple sclerosis. Am J Hum Genet. 2001;69:900–903. [PubMed]
  • MacKay K, Eyre S, Myerscough A, Milicic A, Barton A, Laval S, et al. Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the united kingdom. Arthritis Rheum. 2002;46:632–639. [PubMed]
  • Maier LM, Hafler DA. The developing mosaic of autoimmune disease risk. Nat Genet. 2008;40:131–132. [PubMed]
  • Margulies EH, Blanchette M, Haussler D, Green ED. NISC Comparative Sequencing Program. Identification and characterization of multi-species conserved sequences. Genome Res. 2003;13:2507–2518. [PubMed]
  • Martin ER, Monks SA, Warren LL, Kaplan NL. A test for linkage and association in general pedigrees: The pedigree disequilibrium test. Am J Hum Genet. 2000;67:146–154. [PubMed]
  • McCauley JL, Kenealy SJ, Margulies EH, Schnetz-Boutaud N, Gregory SG, Hauser SL, et al. SNPs in multi-species conserved sequences (MCS) as useful markers in association studies: A practical approach. BMC Genomics. 2007;8:266. [PMC free article] [PubMed]
  • Ott J. The number of families required to detect or exclude linkage heterogeneity. Am J Hum Genet. 1986;39:159–165. [PubMed]
  • Poser CM, paty DW, Scheinberg L, McDonald WI, Davis FA, Ebers GC, et al. New diagnostic criteria for multiple sclerosis: Guidelines for research protocols. Ann Neurol. 1983;13:227–231. [PubMed]
  • Rubio JP, Stankovich J, Field J, Tubridy N, Marriott M, Chapman C, et al. Replication of KIAA0350, IL2RA, RPL5 and CD58 as multiple sclerosis susceptibility genes in australians. Genes Immun. 2008 [PubMed]
  • Sawcer S, Jones HB, Feakes R, Gray J, Smaldon N, Chataway J, et al. A genome screen in multiple sclerosis reveals susceptibility loci on chromosome 6p21 and 17q22. Nat Genet. 1996;13:464–476. [PubMed]
  • Sawcer S, Maranian M, Setakis E, Curwen V, Akesson E, Hensiek A, et al. A whole genome screen for linkage disequilibrium in multiple sclerosis confirms disease associations with regions previously linked to susceptibility. Brain. 2002;125:1337–1347. [PubMed]
  • Sawcer S, Ban M, Maranian M, Yeo TW, Compston A, Kirby A, et al. A high-density screen for linkage in multiple sclerosis - international multiple sclerosis genetics consortium. Am J Hum Genet. 2005;77:454–467. [PubMed]
  • Schaffer AA, Gupta SK, Shriram K, Cottingham RW. Avoiding recomputation in linkage analysis. Hum Hered. 1994;44:225–237. [PubMed]
  • Shai R, Quismorio FP, Jr, Li L, Kwon OJ, Morrison J, Wallace DJ, et al. Genome-wide screen for systemic lupus erythematosus susceptibility genes in multiplex families. Hum Mol Genet. 1999;8:639–644. [PubMed]
  • Weber F, Fontaine B, Cournu-Rebeix I, Kroner A, Knop M, Lutz S, et al. IL2RA and IL7RA genes confer susceptibility for multiple sclerosis in two independent european populations. Genes Immun. 2008;9:259–263. [PubMed]
  • Yaouanq J, Semana G, Eichenbaum S, Quelvennec E, Roth MP, Clanet M, et al. Evidence for linkage disequilibrium between HLA-DRB1 gene and multiple sclerosis. the french research group on genetic susceptibility to MS. Science. 1997;276:664–665. [PubMed]