- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
Three allele combinations associated with Multiple Sclerosis
© Favorova et al; licensee BioMed Central Ltd. 2006
Received: 17 March 2006
Accepted: 26 July 2006
Published: 26 July 2006
Multiple sclerosis (MS) is an immune-mediated disease of polygenic etiology. Dissection of its genetic background is a complex problem, because of the combinatorial possibilities of gene-gene interactions. As genotyping methods improve throughput, approaches that can explore multigene interactions appropriately should lead to improved understanding of MS.
286 unrelated patients with definite MS and 362 unrelated healthy controls of Russian descent were genotyped at polymorphic loci (including SNPs, repeat polymorphisms, and an insertion/deletion) of the DRB1, TNF, LT, TGFβ1, CCR5 and CTLA4 genes and TNFa and TNFb microsatellites. Each allele carriership in patients and controls was compared by Fisher's exact test, and disease-associated combinations of alleles in the data set were sought using a Bayesian Markov chain Monte Carlo-based method recently developed by our group.
We identified two previously unknown MS-associated tri-allelic combinations:
-509TGFβ1*C, DRB1*18(3), CTLA4*G and -238TNF*B1,-308TNF*A2, CTLA4*G, which perfectly separate MS cases from controls, at least in the present sample. The previously described DRB1*15(2) allele, the microsatellite TNFa9 allele and the biallelic combination CCR5Δ32, DRB1*04 were also reidentified as MS-associated.
These results represent an independent validation of MS association with DRB1*15(2) and TNFa9 in Russians and are the first to find the interplay of three loci in conferring susceptibility to MS. They demonstrate the efficacy of our approach for the identification of complex-disease-associated combinations of alleles.
Investigation of polygenic human diseases, which arise from the combined contribution of multiple independently acting and/or interacting polymorphic genes, remains a great challenge [1–3]. A common constituent of the complexity that underlies the susceptibility to polygenic diseases is heterogeneity [3, 4]. MS [MIM 126200] is an immune-mediated hereditary disease [5, 6], and can be considered as a prototype for polygenic human diseases . The results of linkage-based whole genome screen studies [7, 8] and a global meta-analysis [9, 10] document the concept that MS is the result of the interaction of several genes. The effects of individual genes are small or modest, making association studies more informative than others because of their greater statistical power [1, 11]. Association testing is extensively employed in candidate-gene studies, which are usually conducted in population-based case-control studies. To date, the HLA class II DRB1*1501/DQA1*0102/DQB1*0602 (DR2) haplotype is the only region repeatedly confirmed as being associated with MS in most Caucasians [5, 10]. Other candidate genes for MS predisposition studies have been selected mainly because their encoded proteins are involved in autoimmune pathogenesis. These include genes for immunorelevant molecules such as cytokines, cytokine receptors, immunoglobulins, T-cell receptors and specific adapter protein, potential autoantigens of the myelin sheath, ICAM1, and others .
Some studies searched for candidate gene combinations as MS risk factors; however, these studies have not yet extended beyond MS associations with alleles of each candidate gene coupled with HLA DRB1 alleles [13–19]. To examine the possibility that the combined effect of definite genes is a risk factor for a polygenic disease, it is necessary to explore a massive number of potential combinations of allelic variants identified at candidate-gene polymorphic loci. This number grows exponentially with the number of candidate variants that may interact, causing computational and statistical restrictions on the use of standard enumerational methodologies. In a recent paper , we described a novel algorithm based on Markov chain Monte Carlo exploration using a Bayesian statistical basis, APSampler, which allows the exploration of genotypes tied to phenotypic trait levels to identify possible combinations of allelic variants at multiple loci that could affect disease development. The aim of this study is to simultaneously examine multiple candidate genes in single groups of unrelated MS patients and healthy unrelated controls, all of Russian descent, and to search for disease-associated combinations of allelic variants at multiple loci using our novel nonparametric methodology .
Subjects and DNA samples
Two hundred eighty six unrelated patients (110 men and 176 women, mean age 33 ± 12 years) had a diagnosis of MS . Of these 187 had a relapsing-remitting MS course, 39 a primary progressive MS course, and 60 were secondary progressive. The mean age at onset was 23 ± 9 years. Three hundred sixty-two unrelated controls (203 men and 159 women, mean age 30 ± 11 years) were studied. All controls were free of acute or chronic internal and neurological diseases as determined by physical examinations. All subjects were living in the Moscow area; both their parents were ethnic Russians. Informed consent to the study was obtained from all participants and was approved by the local Ethical Committee. Blood was adjusted to 25 mM EDTA. Genomic DNA was isolated from 5 mL of peripheral blood by phenol-chloroform extraction using standard procedures.
Polymorphous loci at or near genes of immune response included in the database
Polymorphism type* [refSNP ID]**
Names of alleles considered
Method of analysis (the restriction endonuclease used)
Numbers of genotyped MS patients/controls
Allele groups corresponding to serological specificities DR1-DR18(3)
SNP -376A→G [rs1800750]
SNP -308G→A [rs1800629]
SNP -238A→G [rs361525]
SNP -509C→T [rs17551290]
SNP +72 wild type→C insertion
SNP +869T→C (10Leu→Pro)
SNP +915G→C (25Arg→Pro)
SNP +1632C→T (263Thr→Ile)
Wild type→32 base pair deletion
wt, Δ 32
SNP +49A→G (17Thr→Ala) [rs231775]
HLA DRB1 gene
For genomic typing of the DRB1 gene, a sequence-specific primer (PCR-SSP) technique was used. The two-step PCR allowed amplification of all known DRB1 alleles and their separation into groups corresponding to the specificities from DR1 to DR18 .
TNFa and TNFb microsatellites
For a length polymorphism analysis of (AC)n and (TC)n microsatellites, designated as TNFa and TNFb, which are located 3.5 kb upstream of the LT gene, nested PCR was used. The second PCR was carried out in the presence of [α-32P]dATP, then PCR products were treated by the Klenow fragment of DNA polymerase I and electrophoresed in an 8% polyacrylamide sequencing gel .
SNPs in TNF gene
Genotyping of the -238A→G polymorphism was performed by analysis of restriction fragment length polymorphism of PCR products (PCR-RFLP method). PCR amplification was carried out using the exactly homologous forward primer, while the reverse primer contained two sequence mismatches, which made it possible to evaluate mutation status using the restriction enzyme BamHI . Analysis of the -308G→A polymorphism was performed by PCR-SSP . For genotyping of the -376A→G polymorphism, the PCR-SSP method was also used. The forward primers are 5'-CTT TTT CCT GCA TCC TGT CTG GAA A-3' for -376A and 5'-CTT TTT CCT GCA TCC TGT CTG GAA G-3' for -376G, the common reverse primer was 5'-TTC TGT CTC GGT TTC TTC TCC ATC G-3'. These primers were constructed using free on-line software GeneFisher-Interactive PCR Primer Design . PCR was performed in 10 μL volumes containing 0.1 units of Taq polymerase (from SileksM, Moscow, Russia), 200 ng of genomic DNA, 5 pmol of the reverse PCR primer, 5 pmol of the sequence-specific PCR primer for allele A or 2.5 pmol of the sequence-specific PCR primer for allele G, four dNTPs (each at 5 mM), 70 mM Tris HCl pH 9.0, 20 mM (NH4)2SO4, 0.025% Tween 20, 0.025% NP-40, 1.0 mM MgCl2. Thermocycling consisted of 35 cycles of 92°C for 60 s, 60°C for 90 s, 72°C for 90 s. The resulting fragment length was 238 bp.
SNPs in LT gene
Genotyping of the +252G→A and +319C→G polymorphisms was performed by PCR-RFLP method. For analysis of the SNP +252G→A, PCR amplification was performed and the product aliquots were digested using the restriction enzyme NcoI. In parallel, other aliquots of the same PCR products were digested using the restriction enzyme Alw21I (AspHI) for analysis of SNP +319C→G .
SNPs in TGFβ1 gene
For genotyping of five SNPs in the TGFβ1 gene, four PCR fragments were obtained (one PCR fragment was common for the SNP +869T→C (10Leu→Pro) and SNP +915G→C (25Arg→Pro)). Then resulting fragments were immobilized on Hybond-N+ membrane (Amersham Pharmacia Biotech) followed by hybridization with sequence-specific oligonucleotide probes (PCR-SSO method). PCR conditions, primers and probes were as those described in , hybridization procedure was modified according to  as follows. Prehybridization was performed for at least 2 h at 52°C in 3.0 M tetramethylammonium chloride (TMAC), 50 mM Tris HCl pH 8.0, 2 mM EDTA, 5×Denhardt's solution, 0.1% SDS, 100 ng/mL heat-denatured herring sperm DNA and 10–20 pmol/mL of unlabeled oligonucleotide probe specific for the other allele. Then 10–20 pmol/mL of [γ-32P]ATP-labeled allele-specific oligonucleotide was added and hybridization was performed for at least 2 h at 52°C. The membranes were washed at room temperature in 2×SSC for 30 min followed by 2 × 15 min in 2×SSC, 0.1% SDS. Then the membranes were washed twice at 58°C in 3.0 M TMAC, 50 mM Tris HCl pH 8.0, 2 mM EDTA, 0.1% SDS, rinsed in 2×SSC and autoradiographed.
32 base pair deletion in CCR5 gene
CCR5 genomic typing was performed using one-step PCR with primers flanking the region of the 32-nucleotide deletion .
SNP in CTLA4 gene
Genotyping of SNP +49A→G in CTLA4 gene was performed by PCR-RFLP analysis using the restriction enzyme BstEII .
The genotypes and personal data for all patients and controls were entered into a database, together with clinical characteristics for patients. Microsoft Visual FoxPro was used to develop a standalone database management system that was used for input and analysis of the data.
We used the APSampler algorithm  that identifies combinations (patterns) of alleles at different loci that are potentially associated with a phenotypic trait. Here, we provide a brief description of the algorithm. The a posteriori probability for support of each pattern given the data is evaluated using a likelihood obtained from a battery of conditional rank sum statistics that insures that each pattern is evaluated after removing the effect of other patterns, as in a statistical adjustment for multiple regression. Each step of the algorithm is an update of a current set of allelic patterns by a variation of the Metropolis-Hastings algorithm [32, 33]. The sets of patterns that receive high support from the data are stored as potential results. For each pattern, only those individuals that have all the pattern's alleles genotyped were considered. Thus the algorithm proceeded from incomplete genotypic information, temporarily omitting those individuals that cannot be unambiguously classified into the pattern carriers and non-carriers at a given step. The executable files for Win32 console and for FreeBSD, together with documentation, are available free to academic users (contact AVF or MFO).
In the current application to the MS case-control data, the prior probability that a locus has no effect on the phenotype was set to 0.99. This favored a priori the case in which each locus has no effect on phenotype, so that the search would not explore patterns with a very large number of alleles, which are unlikely to be reliably investigated in this data set. We looked for sets of 2 or 3 patterns.
To compare allelic distributions in MS patients and controls and to verify the APSampler's results, we used Fisher's exact test. The analysis was performed using the GraphPad InStat software package and original computer scripts under Microsoft Visual FoxPro. As a statistical control for systematic genotyping error and population stratification, the expected genotype proportions according to the Hardy-Weinberg equilibrium were calculated and compared to observed genotypes. For polymorphic DRB1 and TNFa loci, the p-values were corrected for the number of comparisons (13 in case of DRB1 allele groups and 12 in case of TNFa alleles) according to the Bonferroni method. The p-values and corrected p-values (p corr ) were considered to be significant at a level smaller than 0.01. The odds ratio (OR) was calculated with a 95% confidence interval (CI).
For the haplotype frequency estimation, for the linkage disequilibrium test, and for the Hardy-Weinberg disequilibrium test, the population genetics data-analyzing software Arlequin (version 2.0)  was used. We considered a linkage disequilibrium or a Hardy-Weinberg disequilibrium to be significant when the p-value from the Arlequin output was less than or equal to 0.05. For estimation of haplotype frequencies, the default parameters were used. To assess the significance of the APSampler-discovered patterns, we computed the Empirical Bayes false discovery rate, which is the estimated fraction of patterns unrelated to the phenotype among those that have a significance score as high or higher than the pattern reported . The vector of disease labels was permuted in a balanced way by assigning the same number of diseased and healthy individuals to two groups generated by permutation. We repeated the permutation 100 times and examined all possible combinations with the same number of alleles as the pattern considered. For each permutation and each combination, we computed a test of association and counted the fraction of these tests that were larger than the observed test for the pattern.
Phenotypic frequencies (carriership) of the three-allele combinations (trios) and their two-element subsets in Russian MS patients and healthy individuals
MS patients, N (%)
Controls, N (%)
OR (95% CI)##
-509TGFβ1*C, DRB1*18(3), CTLA4*G (trio 1)*
-238TNF*B1,-308TNF*A2, CTLA4*G (trio 2)**
Based on the results presented, several (partially overlapping) subgroups may be identified in the common group of Russian MS patients, depending on carriership of distinct minimal patterns of susceptibility including one (DRB1*15(2) or TNFa9), two (CCR5Δ32, DRB1*04) and three (-509TGFβ1*C, DRB1*18(3), CTLA4*G or -238TNF*B1,-308TNF*A2, CTLA4*G) alleles of candidate genes. These results are evidence of the genetic heterogeneity of MS.
MS associations with DRB1*15(2) and TNFa9 were previously identified in our studies of independent groups of ethnic Russians [36, 37] and were replicated here for a new dataset of patients and control subjects. Thus, the HLA class II DRB1*15(2) is validated in this study as being associated with MS in Russians as in most other Caucasians. As follows from our data, DRB1*15(2) and TNFa9 represent two independent predisposing markers, in agreement with the proposal that two MS susceptibility loci exist within the MHC [5, 38].
In previous studies, identification of MS-predisposing combinations of allelic variants at multiple loci consisted in stratification of affected and unaffected individuals mostly according to the carriership of MS-predisposing DR2 haplotype or its constituents and a subsequent pair comparison of phenotypic frequencies of distinct alleles of another gene of interest in subgroups of MS patients and controls. This kind of analysis found reliable associations of MS in subgroups of DRB1*15(2) (or DQB1*0602)-positive or negative individuals with the carriership of some alleles of genes coding TCRβ receptor , TGFβ1 , CTLA4 , ICAM-1  and interleukin 4 receptor . In our recent studies, we have extended this approach and stratified individuals according to any DRB1 phenotypes. For DRB1*04-positive individuals, associations of MS with CCR5Δ 32 mutation  and with alleles of myelin basic protein (MBP) gene  were found. The latter association was found also for DRB1*05 positive individuals .
A novel nonparametric methodology used in this paper provides the capability to explore potential combinations of more than two allelic variants of polymorphic candidate genes. Two previously unknown tri-allelic combinations: -509TGFβ1*C, DRB1*18(3), CTLA4*G and -238TNF*B1,-308TNF*A2, CTLA4*G, which perfectly separate cases from controls in the present sample, were identified. Importantly, at least within our data set, non-overlapping subgroups of individuals bearing predisposing trios 1 and 2 constitute about 5% and 9% of considered MS patients, respectively, and 0% of considered healthy controls (see Table 2). This suggests that carriership of these patterns might predict MS development. Due to the linkage disequilibrium of chromosomal loci, a genetic epidemiological approach cannot prove unambiguously that a disease-associated gene is causal. However, a biological role for the DRB1, CCR5, TGFβ1, TNF and CTLA4 gene products in the pathogenesis of MS is plausible, and supports the idea that the genes are actual MS susceptibility genes.
Trios 1 and 2 have striking similarities that may determine their MS-predisposing properties due to dysregulation of inflammatory pathways by protein products encoded by their genes. First, both trios include the allele G of the gene for co-stimulatory molecule CTLA4, which is an important inhibitor of T-cell activation . Carriers of the +49CTLA4*G allele in exon 1, coding for the peptide leader sequence, are characterized by a reduced CTLA4 inhibitory function , i.e. by impaired negative regulation of the immune response. Second, both trios include alleles of cytokine genes which promote inflammatory immune response owing to decreased level of antiinflammatory cytokine TGFβ1 (trio 1) or increased level of proinflammatory cytokine TNF (trio 2). Indeed, it is known that the circulating concentration of TGFβ1 is under genetic control, being lower in carriers of the promoter -509TGFβ1*C allele, in comparison with -509TGFβ1*T allele . The -308TNF*A2 allele included in trio 2 is more actively transcribed than -308TNF*A1 . The sequence of the promoter region of -238TNF*B1, another trio 2 TNF allele, suggests that it is also associated with increased TNF expression . In trio 1 (-509TGFβ1*C, DRB1*18(3), CTLA4*G) carriers, a two-hit repression of negative regulation of the immune response is coupled with a contribution of the HLA class II DRB1*18(3), which is not associated alone with MS (see Figure 2A). Characteristics of autoantigen presentation in MS by the products encoded by DRB1*18(3) have not been studied; however, it is clear that presence of this allele in trio 1 may be an important factor of the individual immune response regulation, supplementing effects provided by two other alleles of the trio.
Our results support the notion that genetic susceptibility to MS arises as a result of contribution of several predisposing alleles involved in the autoimmune inflammatory response. The nature of the interplay between the alleles in the trios remains unclear. If it arises as a result of a cumulative effect of multiple hits, these patterns provide a picture of small, individually insufficient changes combining to provide an overall significant response leading to predisposition to MS. However, this scenario does not exclude possible epistatic interactions of genes involved in the trios.
The results presented here provide the first identification of a combination of more than two alleles conferring a genetic predisposition to MS as a complex polygenic disease. The two newly identified trios perfectly separate cases from controls, at least in the present sample, as would occur in a classic monogenic dominant disease, where all observed carriers are patients. Because of sample size limitations, the result can only be considered as preliminary. However, the biological properties of the genes included in the identified trios suggest a coherent picture of dysregulation of inflammatory pathways, implying the validity of the MS-associated patterns and encouraging efforts required to replicate associations in independent data. The APSampler algorithm used to identify the predisposing patterns is highly efficient, as it requires only hours of computation time on a laptop computer for this data set, and is flexible, as it can handle incomplete genotypic information. The algorithm provides a valuable resource for the growing volume of polygenic disease-related genomic data, allowing efficient exploration of such data to identify genetic predisposition and potential therapeutic targets.
We thank Erica Golemis and Kenneth Morgan for very useful discussion and for advice. All authors from Russia thank Russian Foundation for Basic Research (05-04-48982-a and 04-04-49601-a). AVF thanks the Howard Hughes Medical Institute (55005610), the Ludwig Institute for Cancer Research (CRDF 12039), Russian Ministry of Science and Education (State Contract # 02.434.11008) and Russian Academy of Sciences Presidium Program in Molecular and Cellular Biology. GP thanks NSF (DMS034211) and NCI (CA06973). MFO thanks NCI (support under CCSG grant CA06927), the National Library of Medicine (LM008932), and the Pennsylvania Department of Health for support.
- Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005, 6 (2): 95-108. 10.1038/nrg1521.View ArticlePubMedGoogle Scholar
- Thornton-Wells TA, Moore JH, Haines JL: Genetics, statistics and human disease: analytical retooling for complexity. Trends Genet. 2004, 20 (12): 640-647. 10.1016/j.tig.2004.09.007.View ArticlePubMedGoogle Scholar
- Threadgill DW, Hunter KW, Williams RW: Genetic dissection of complex and quantitative traits: from fantasy to reality via a community effort. Mamm Genome. 2002, 13 (4): 175-178. 10.1007/s00335-001-4001-y.View ArticlePubMedGoogle Scholar
- Bomprezzi R, Kovanen PE, Martin R: New approaches to investigating heterogeneity in complex traits. J Med Genet. 2003, 40 (8): 553-559. 10.1136/jmg.40.8.553.View ArticlePubMedPubMed CentralGoogle Scholar
- Herrera BM, Ebers GC: Progress in deciphering the genetics of multiple sclerosis. Curr Opin Neurol. 2003, 16 (3): 253-258. 10.1097/00019052-200306000-00002.View ArticlePubMedGoogle Scholar
- Kalman B, Albert RH, Leist TP: Genetics of multiple sclerosis: determinants of autoimmunity and neurodegeneration. Autoimmunity. 2002, 35 (4): 225-234. 10.1080/08916930290005611.View ArticlePubMedGoogle Scholar
- Kenealy SJ, Babron MC, Bradford Y, Schnetz-Boutaud N, Haines JL, Rimmler JB, Schmidt S, Pericak-Vance MA, Barcellos LF, Lincoln RR, et al: A second-generation genomic screen for multiple sclerosis. Am J Hum Genet. 2004, 75 (6): 1070-1078. 10.1086/426459.View ArticlePubMedPubMed CentralGoogle Scholar
- Sawcer S, Ban M, Maranian M, Yeo TW, Compston A, Kirby A, Daly MJ, De Jager PL, Walsh E, Lander ES, et al: A high-density screen for linkage in multiple sclerosis. Am J Hum Genet. 2005, 77 (3): 454-467. 10.1086/444547.View ArticlePubMedGoogle Scholar
- GAMES, Transatlantic Multiple Sclerosis Genetics Cooperative: A meta-analysis of whole genome linkage screens in multiple sclerosis. J Neuroimmunol. 2003, 143 (1–2): 39-46.Google Scholar
- Kenealy SJ, Pericak-Vance MA, Haines JL: The genetic epidemiology of multiple sclerosis. J Neuroimmunol. 2003, 143 (1–2): 7-12. 10.1016/j.jneuroim.2003.08.005.View ArticlePubMedGoogle Scholar
- Risch NJ: Searching for genetic determinants in the new millennium. Nature. 2000, 405 (6788): 847-856. 10.1038/35015718.View ArticlePubMedGoogle Scholar
- Dyment DA, Ebers GC, Sadovnick AD: Genetics of multiple sclerosis. Lancet Neurol. 2004, 3 (2): 104-110. 10.1016/S1474-4422(03)00663-X.View ArticlePubMedGoogle Scholar
- Alizadeh M, Babron MC, Birebent B, Matsuda F, Quelvennec E, Liblau R, Cournu-Rebeix I, Momigliano-Richiardi P, Sequeiros J, Yaouanq J, et al: Genetic interaction of CTLA-4 with HLA-DR15 in multiple sclerosis patients. Ann Neurol. 2003, 54 (1): 119-122. 10.1002/ana.10617.View ArticlePubMedGoogle Scholar
- Dyment DA, Steckley JL, Morrison K, Willer CJ, Cader MZ, DeLuca GC, Sadovnick AD, Risch N, Ebers GC: TCR beta polymorphisms and multiple sclerosis. Genes Immun. 2004, 5 (5): 337-342. 10.1038/sj.gene.6364091.View ArticlePubMedGoogle Scholar
- Favorova OO, Andreewski TV, Boiko AN, Sudomoina MA, Alekseenkov AD, Kulakova OG, Slanova AV, Gusev EI: The chemokine receptor CCR5 deletion mutation is associated with MS in HLA-DR4-positive Russians. Neurology. 2002, 59 (10): 1652-1655.View ArticlePubMedGoogle Scholar
- Green AJ, Barcellos LF, Rimmler JB, Garcia ME, Caillier S, Lincoln RR, Bucher P, Pericak-Vance MA, Haines JL, Hauser SL, et al: Sequence variation in the transforming growth factor-beta1 (TGFB1) gene and multiple sclerosis susceptibility. J Neuroimmunol. 2001, 116 (1): 116-124. 10.1016/S0165-5728(01)00283-1.View ArticlePubMedGoogle Scholar
- Guerini FR, Ferrante P, Losciale L, Caputo D, Lombardi ML, Pirozzi G, Luongo V, Sudomoina MA, Andreewski TV, Alekseenkov AD, et al: Myelin basic protein gene is associated with MS in DR4- and DR5-positive Italians and Russians. Neurology. 2003, 61 (4): 520-526.View ArticlePubMedGoogle Scholar
- Mirel DB, Barcellos LF, Wang J, Hauser SL, Oksenberg JR, Erlich HA: Analysis of IL4R haplotypes in predisposition to multiple sclerosis. Genes Immun. 2004, 5 (2): 138-141. 10.1038/sj.gene.6364048.View ArticlePubMedGoogle Scholar
- Nejentsev S, Laaksonen M, Tienari PJ, Fernandez O, Cordell H, Ruutiainen J, Wikstrom J, Pastinen T, Kuokkanen S, Hillert J, et al: Intercellular adhesion molecule-1 K469E polymorphism: study of association with multiple sclerosis. Hum Immunol. 2003, 64 (3): 345-349. 10.1016/S0198-8859(02)00825-X.View ArticlePubMedGoogle Scholar
- Favorov AV, Andreewski TV, Sudomoina MA, Favorova OO, Parmigiani G, Ochs MF: A Markov Chain Monte Carlo Technique for Identification of Combinations of Allelic Variants Underlying Complex Diseases. Genetics. 2005, 171 (4): 2113-2121. 10.1534/genetics.105.048090.View ArticlePubMedPubMed CentralGoogle Scholar
- McDonald WI, Compston A, Edan G, Goodkin D, Hartung HP, Lublin FD, McFarland HF, Paty DW, Polman CH, Reingold SC, et al: Recommended diagnostic criteria for multiple sclerosis: guidelines from the International Panel on the diagnosis of multiple sclerosis. Ann Neurol. 2001, 50 (1): 121-127. 10.1002/ana.1032.View ArticlePubMedGoogle Scholar
- Sudomoina MA, Boiko AN, Demina TL, Gusev EI, Boldyreva MN, Trofimov D, Alekseev LP, Favorova OO: Connection of multiple sclerosis in the Russian population with alleles of the major histocompatibility complex DRB1 gene. Mol Biol (Mosk). 1998, 32 (2): 291-296.Google Scholar
- Nedospasov SA, Udalova IA, Kuprash DV, Turetskaya RL: DNA sequence polymorphism at the human tumor necrosis factor (TNF) locus. Numerous TNF/lymphotoxin alleles tagged by two closely linked microsatellites in the upstream region of the lymphotoxin (TNF-beta) gene. J Immunol. 1991, 147 (3): 1053-1059.PubMedGoogle Scholar
- Gallagher G, Eskdale J, Oh HH, Richards SD, Campbell DA, Field M: Polymorphisms in the TNF gene cluster and MHC serotypes in the West of Scotland. Immunogenetics. 1997, 45 (3): 188-194. 10.1007/s002510050188.View ArticlePubMedGoogle Scholar
- He B, Navikas V, Lundahl J, Soderstrom M, Hillert J: Tumor necrosis factor alpha-308 alleles in multiple sclerosis and optic neuritis. J Neuroimmunol. 1995, 63 (2): 143-147. 10.1016/0165-5728(95)00138-7.View ArticlePubMedGoogle Scholar
- GeneFisher. [http://bibiserv.techfak.uni-bielefeld.de/genefisher]
- Roth MP, Nogueira L, Coppin H, Clanet M, Clayton J, Cambon-Thomsen A: Tumor necrosis factor polymorphism in multiple sclerosis: no additional association independent of HLA. J Neuroimmunol. 1994, 51 (1): 93-99. 10.1016/0165-5728(94)90133-3.View ArticlePubMedGoogle Scholar
- Cambien F, Ricard S, Troesch A, Mallet C, Generenaz L, Evans A, Arveiler D, Luc G, Ruidavets JB, Poirier O: Polymorphisms of the transforming growth factor-beta 1 gene in relation to myocardial infarction and blood pressure. The Etude Cas-Temoin de l'Infarctus du Myocarde (ECTIM) Study. Hypertension. 1996, 28 (5): 881-887.View ArticlePubMedGoogle Scholar
- Bignon JD, Fernandez-Vina MF: Technical Handbook: Twelfth International Histocombatibility Workshop and Conference. 1996, Paris: HLA et Me'decineGoogle Scholar
- Sandford AJ, Pare PD: Direct PCR of small genomic DNA fragments from serum. Biotechniques. 1997, 23 (5): 890-892.PubMedGoogle Scholar
- Marron MP, Raffel LJ, Garchon HJ, Jacob CO, Serrano-Rios M, Martinez Larrad MT, Teng WP, Park Y, Zhang ZX, Goldstein DR, et al: Insulin-dependent diabetes mellitus (IDDM) is associated with CTLA4 polymorphisms in multiple ethnic groups. Hum Mol Genet. 1997, 6 (8): 1275-1282. 10.1093/hmg/6.8.1275.View ArticlePubMedGoogle Scholar
- Robert C, Casella G: Monte Carlo Statistical Methods. 1999, New York: Springer VerlagView ArticleGoogle Scholar
- Metropolis N, Rosenbluth A, Rosenbluth M, Teller A, Teller E: Equations of state calculations by fast computing machines. J Chem Phys. 1953, 21: 1087-1091. 10.1063/1.1699114.View ArticleGoogle Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin: A software for population genetics data analysis. Ver 2.000. 2000, Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva; Arlequin, [http://cmpg.unibe.ch/software/arlequin3/]Google Scholar
- Efron B, Tibshirani R, Storey JD, Tusher V: Empirical Bayes Analysis of a Microarray Experiment. J Am Statist Assoc. 2001, 96 (456): 1151-1160. 10.1198/016214501753382129.View ArticleGoogle Scholar
- Boiko AN, Gusev EI, Sudomoina MA, Alekseenkov AD, Kulakova OG, Bikova OV, Maslova OI, Guseva MR, Boiko SY, Guseva ME, et al: Association and linkage of juvenile MS with HLA-DR2(15) in Russians. Neurology. 2002, 58 (4): 658-660.View ArticlePubMedGoogle Scholar
- Gusev E, Sudomoina M, Boiko A, Deomina T, Favorova O: TNF gene polymorphisms: associations with multiple sclerosis susceptibility and severity. Frontiers in multiple sclerosis. Edited by: O A, H O. 1997, London: Martin Dunitz, 35-41.Google Scholar
- de Jong BA, Huizinga TW, Zanelli E, Giphart MJ, Bollen EL, Uitdehaag BM, Polman CH, Westendorp RG: Evidence for additional genetic risk indicators of relapse-onset MS within the HLA region. Neurology. 2002, 59 (4): 549-555.View ArticlePubMedGoogle Scholar
- Alegre ML, Frauwirth KA, Thompson CB: T-cell regulation by CD28 and CTLA-4. Nat Rev Immunol. 2001, 1 (3): 220-228. 10.1038/35105024.View ArticlePubMedGoogle Scholar
- Maurer M, Loserth S, Kolb-Maurer A, Ponath A, Wiese S, Kruse N, Rieckmann P: A polymorphism in the human cytotoxic T-lymphocyte antigen 4 (CTLA4) gene (exon 1 +49) alters T-cell activation. Immunogenetics. 2002, 54 (1): 1-8. 10.1007/s00251-002-0429-9.View ArticlePubMedGoogle Scholar
- Grainger DJ, Heathcote K, Chiano M, Snieder H, Kemp PR, Metcalfe JC, Carter ND, Spector TD: Genetic control of the circulating concentration of transforming growth factor type beta1. Hum Mol Genet. 1999, 8 (1): 93-97. 10.1093/hmg/8.1.93.View ArticlePubMedGoogle Scholar
- Uboldi de Capei MU, Dametto E, Fasano ME, Rendine S, Curtoni ES: Genotyping for cytokine polymorphisms: allele frequencies in the Italian population. Eur J Immunogenet. 2003, 30 (1): 5-10. 10.1046/j.1365-2370.2003.00361.x.View ArticlePubMedGoogle Scholar
- Grove J, Daly AK, Bassendine MF, Day CP: Association of a tumor necrosis factor promoter polymorphism with susceptibility to alcoholic steatohepatitis. Hepatology. 1997, 26 (1): 143-146.View ArticlePubMedGoogle Scholar
- Marsh SG, Albert ED, Bodmer WF, Bontrop RE, Dupont B, Erlich HA, Geraghty DE, Hansen JA, Mach B, Mayr WR, et al: Nomenclature for factors of the HLA system, 2002. Hum Immunol. 2002, 63 (12): 1213-1268. 10.1016/S0198-8859(02)00769-3.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2350/7/63/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.