Human Molecular Genetics, 2001, Vol. 10, No. 25 2907-2916
© 2001 Oxford University Press
Dissection of the HLA association with multiple sclerosis in the founder isolated population of Sardinia
1Dipartimento di Neuroscienze, University of Cagliari, Centro Sclerosi Multipla, Ospedale Binaghi, Via Is Guadazzonis 2, Cagliari 09126, Italy and 2Laboratorio di Immunogenetica, Dipartimento di Scienze Biomediche e Biotecnologie, University of Cagliari, Ospedale Microcitemico, Via Jenner, Cagliari 09121, Italy
Received August 2, 2001; Revised and Accepted October 8, 2001.
| ABSTRACT |
|---|
|
|
|---|
Several studies have indicated that multiple sclerosis (MS) is associated and linked to the major histocompatibility complex (MHC)/human leukocyte antigen (HLA) region of chromosome 6p21.3, but the exact location and nature of the primarily associated locus within the HLA complex is still controversial and largely presumptive. By linkage disequilibrium mapping, we have systematically investigated this chromosome region in the founder population of Sardinia to determine the relative associations of the various loci with MS. An overall 11.4 Mb region, which encompasses the whole HLA complex, was scanned with 19 microsatellite markers and with single nucleotide polymorphisms within 12 functional candidate genes and assessed for MS association using the extended transmission disequilibrium test (ETDT). A peak of association represented by the three adjacent DRB1, DQA1 and DQB1 loci was detected in the class II region. Two additional less significant areas of association were detected, respectively, in the centromeric side of the class II region at the DPB1 locus and, telomeric of the classically defined class I loci, at the D6S1683 microsatellite. Conditional ETDT analysis indicated that these regions of association could be independent of each other. Within the main peak of association, DRB1 and DQB1 contribute to the disease association independently of each other whereas DQA1 had no detectable primary genetic effects. We evaluated the haplotype distribution at the region showing the strongest association and found five DQB1DRB1 haplotypes positively associated with MS in Sardinia. These consistently included all the haplotypes previously found associated with MS in the various human populations, thus supporting a primary effect of the products of these loci in MS. Overall these results are consistent with a multilocus model of the MHC encoded susceptibility to MS.
| INTRODUCTION |
|---|
|
|
|---|
Multiple sclerosis (MS) is a chronic inflammatory and demyelinating disease of the central nervous system which leads to severe neurological disability due to axonal loss (1). The disease, which is more common in young adults, is strongly clustered in families: the sibling lifetime risk/population prevalence ratio,
s, ranges from 15 to 20 (2). Adopted siblings of MS patients have the same risk of MS as individuals in the general population; thus the observed familial clustering is more likely due to sharing of genes rather than environmental factors (3). The rapid reduction in risk from monozygotic twins (30.8%), to full siblings (3.46%), to half siblings (1.47%) and to cousins (0.88%) suggests the presence of a polygenic trait with multiplicative interactions between the various disease loci (2,4,5). Three independent whole genome scans were not able to locate major predisposing loci (69). However, some of these linkage studies demonstrated a modest, but definite effect of the major histocompatibility complex (MHC)/human leukocyte antigen (HLA) region. This is in agreement with a long series of case-control studies showing an association of specific HLA class II haploypes with MS (10). However, despite this prior evidence indicating that MS is linked (68,11) and associated (reviewed in 10) to the MHC region, to date there are no systematic studies addressed to fine map the primary disease associated locus (i). Moreover, the association of MS with the HLA class II haplotypes appears to be heterogeneous in different populations. The most common association is with the DRB1*1501DQA1*0102DQB1*0602 haplotype (the serologically defined DR2, DQw1 haplotype) detected in European (12) and in European-derived populations (13) as well as in Ashkenazi and non-Ashkenazi Jews (14), in Turkish (15) and in Mexican Mestizo (16) patients. In the Sardinian population, the disease has been found to be associated with the DRB1*0301DQA1*0501DQB1*0201 and DRB1*0405DQA1*0501DQB1*0301 haplotypes (17,18). Until now little effort has been made in the identification of the disease loci within the HLA/MHC. Within the MHC region, identification of the primary variant responsible for MS susceptibility is complicated by the high degree and irregular distribution of the linkage disequilibrium (LD) observed between the different loci (19,20). Proof of a primary association of any HLA locus must take into account the LD patterns and therefore requires dense maps of markers, very large clinical resources and statistical methods that are able to delineate the relative associations and the individual contributions of the various loci once LD has been taken into account. Moreover, even the choice of the population to study appears critical in LD mapping surveys. The Mediterranean island of Sardinia appears well suited to identify the MHC loci involved in MS predisposition. With approximately 140 cases per 100 000 inhabitants, it has one of the highest incidences of MS in Europe (21). Hence, large clinical resources can be collected on the island. The Sardinian population does not show large-scale genetic heterogeneity (22). This reduces the risk of artifacts owing to population admixture, which can seriously complicate case-control studies. Sardinia is a genetic isolate; the present population is the result of a fixation of alleles and haplotypes, rare or absent elsewhere (23). Since the patterns of LD around the etiologic allele(s) could vary in different populations, this genetic differentiation might be extremely useful for fine-mapping analysis in order to define which HLA alleles and motifs are shared in the patients from distantly related populations. These inter-population shared regions are most likely to contain the etiologic polymorphisms.
We mapped the main HLA disease component to the loci DRB1 and DQB1 encoding the ß subunits of the HLA class II presenting molecules. Other additional genetic effects are associated with the locus encoding the ß subunit of the third class II presenting molecule DPB1, and at the D6S1683 microsatellite.
| RESULTS |
|---|
|
|
|---|
Linkage disequilibrium mapping
We initially evaluated a sample set of 375 Sardinian MS families using the extended transmission disequilibrium test (ETDT). Nineteen microsatellite markers and a number of single nucleotide polymorphisms (SNPs) defining allelic variation at 12 expressed genes across an 11.4 Mb region which encompasses the whole HLA complex, were scored. Six markers in three regions showed some evidence of association at or below the 5% significance level after correction by 31, the number of markers tested (Fig. 1). A main area of association was located in the class II sub-region and peaked at DQB1 (P = 4.4 x 106, Pc = 1.4 x 104) including also D6S2447 (P = 1.5 x 104, Pc = 4.7 x 103), DQA1 (P = 1.1 x 105, Pc = 3.4 x 104) and DRB1 (P = 1.1 x 105, Pc = 3.4 x 104). In the centromeric portion of the class II region there was a second area of association defined by DPB1 (P = 2.5 x 104, Pc = 7.8 x 103). Evidence for a third potential association was located telomeric of the classical class I loci and peaked at D6S1683 (P = 9.1 x 104, Pc = 2.8 x 102).
|
Next, we wanted to assess whether these associations were independent of each other or whether some of them reflected hitchhiking effects owing to the LD between these loci. To address this, we used a variant of the ETDT, the conditional ETDT (CETDT) (Materials and Methods), which allowed us to analyse the overall effects of allelic variation at these loci once LD between them had been taken into account. To provide more power for this conditional analysis, an additional 115 Sardinian MS families were genotyped (total number of families analysed 490) at the peak loci in the three regions of association (DPB1, DQB1, DQA1, DRB1 and D6S1683) (Table 1). We initially evaluated, with the CETDT, the individual contribution of the DQB1, DQA1 and DRB1 loci to MS within the individual peak of association. DRB1 conditioned on DQB1DQA1 and DQB1 conditioned on DRB1DQA1 were significantly independent at or below the 5% level of significance in these 490 MS families (CETDT P = 4.9 x 103 and 5.3 x 102, respectively). The DQA1 locus, when conditioned on DQB1DRB1 had no significantly detectable genetic effect in this population (CETDT P = 0.2). The association of DQB1DRB1 was also independent of DPB1 and D6S1683 (CETDT P = 3.6 x 106 and 1.7 x 103, respectively).
|
We evaluated the individual contribution of the other areas of association. In the most centromeric region, DPB1 was analysed conditional on DQB1DRB1 and D6S1683, respectively, and found that it was associated in all cases (CETDT P = 3.5 x 104 and 2.7 x 102, respectively). In the most telomeric area of association, D6S1683 was independent of DQB1DRB1 and DPB1 (CETDT P = 2.0 x 102 and 7.1 x 103, respectively).
A critical issue in LD mapping projects is the marker spacing and intermarker LD within the map. Based on previous analyses showing that the LD is stronger in the HLA class I and III in comparison with the class II region (20,24,25), we used an uneven distribution of markers with higher marker density in the class II than in the class I and III. This distribution is also consistent with the observed hot spots for recombination detected within the class II region (19,20,26). Nevertheless, we wanted to assess the coverage of the map we used in this study by checking the intermarker LD, measured as normalized D' values and corresponding P-values (Materials and Methods), between the various markers considered in this study. As shown in Figure 2, there was good coverage, with strong LD for all the contiguous markers in the HLA regions of interest. Moreover, as found with the CETDT analysis, it is also evident that the average LD between the markers defining the peaks of association is consistent with the independence of effects.
|
Association of the individual alleles and haplotypes with MS
By analysis of the relative predisposition of the two-locus DRB1DQB1 haplotypes (Table 2), we can see which haplotypes underlie the associations. We have included in this analysis the familial sample set of 490 families (previously considered in the CETDT analysis) and an additional independent case-control sample set of 378 sporadic patients and 631 newborns. The results obtained in these two data sets were very similar and were merged to establish the relative associations of the various haplotypes with MS in the Sardinian population (Materials and Methods). Overall, there were five DRB1DQB1 positively associated haplotypes going from HLADRB1*1303DQB1*0301 [odds ratio (OR) = 2.7, P = 2.8 x 104], followed by DRB1*0405DQB1*0301 (OR = 2.2, P = 6.0 x 107), DRB1*0301DQB1*0201 (OR = 1.7, P = 6.1 x 1012), DRB1*1501DQB1*0602 (OR = 1.6, P = 4.5 x 102) and DRB1*0405DQB1*0302 (OR = 1.5, P = 1.4 x 102). There were several DRB1DQB1 haplotypes showing an apparent negative association with MS (Table 2). However, no significantly negatively associated haplotypes were present when displacement effects driven by the positively associated haplotypes were taken into account by using the relative predispositional effect (RPE) method (27) (Table 3). Additionally, it should be noted that the DRB1*1601DQB1*0502 haplotype, which shows an apparent negative association with the disease, and the DRB1*1501DQB1*0602, which is positively associated with the disease in this sample set and in several populations, share the same DQA1*0102 allele, further suggesting that the DQA1 locus is not contributing in a major way to the genetic susceptibility to MS.
|
|
We also considered the association of particular alleles at the DPB1 and D6S1683 loci in the 490 MS families (Tables 4 and 5). At DPB1, the DPB1*0301 allele was positively associated with the disease (OR = 1.6, P = 1.4 x 105; Table 4). As for the DRB1DQB1 haplotypes also at DPB1 no significantly negatively associated alleles were present when displacement effects driven by positively associated alleles were taken into account by using the RPE method (27) (data not shown).
|
At D6S1683, allele 4 (186 bp) and allele 3 (184 bp) were, respectively, significantly increased and decreased in patients (OR = 1.5, P = 2.8 x 105 and OR = 0.5, P = 2.0 x 105, respectively; Table 5). The negative associations of allele 3 at D6S1683 appear to be genuine in that significant results for this allele were obtained also using the RPE method (27) (data not shown).
|
| DISCUSSION |
|---|
|
|
|---|
We have analysed a large sample set of MS patients from the isolated founder population of Sardinia to establish the disease-associated component(s) within the MHC/HLA region. Using a LD mapping approach, we have found that within this region the DRB1DQB1 loci are associated with MS in this population. Moreover, at least two other markers, namely DPB1 and the microsatellite D6S1683, were also associated with risk of MS. Conditional analysis revealed that these multiple associations were independent of each other after taking into account LD. The overall genetic contributions of the HLA loci to MS susceptibility were relatively modest compared to other HLA-associated diseases such as type 1 diabetes (T1D) (25). For instance, contrast the strength of association and the statistical significance of the main peak in the class II region at DRB1 observed in 257 Sardinian T1D families (ETDT = 252.5, P =3.6 x 1042) (25) and in the 375 MS families reported in this study (ETDT = 73.4, P = 1.0 x 105). These results are also consistent with the previously reported relatively low contribution of the MHC/HLA region to the familial aggregation of susceptibility to MS (9,28). Despite the relatively small genetic effect of the HLA region in MS predisposition, a number of DRB1DQB1 haplotypes conferring susceptibility to MS were detected in this large Sardinian sample and consistently accounted for the main associations previously reported in the other human populations. For instance, the DRB1*1501DQA1*0201DQB1*0602 and the DRB1*1303DQA1*0501DQB1*0301 haplotypes were found to be respectively associated with MS in all the northern European populations (12,13) and in non-Ashkenazi Jews (29) as well as patients from Minnesota (USA) (30) and were among the most positively associated haplotypes also in this study of Sardinian patients. These haplotypes are too rare in the general Sardinian population (1.4 and 0.8%, respectively) to have permitted the detection of their positive association with MS in previous studies on smaller sample sets (18). Similarly, the DRB1*0405DQA1*03DQB1*0302 and the DRB1*0301DQA1*0501DQB1*0201 haplotypes, associated with MS in the Sardinian population, have been found to be associated with MS, respectively, also in Turkey (15) and in some studies of northern European (31) and non-European (16) patients. The presence of a multiallelic association at these class II loci, is also suggested by a recent report in which some evidence of linkage was obtained not only on affected sib pair families with parents carrying DRB1*15 but also in a subgroup of affected sib pair families with parents negative for DRB1*15 (28). These findings indicate that, using an adequate sample size to provide sufficient statistical power and an informative population, virtually all the worldwide HLA-susceptibility determinants for MS can be found.
These similarities between distinct populations indicate a primary association of these class II loci with MS, thus supporting a direct role of their products in the disease pathogenesis. Also, alleles of the locus HLADPB1, encoding for the ß subunit of the third HLA class II presenting molecule, might be primarily associated with MS. This is consistent with the fact that the DPB1*0301 allele found to be associated with MS in this Sardinian sample, was also associated with the disease in the distantly related population of Japan (32) and of Australia (33). Moreover, the DPB1*0301 allele was found to preferentially restrict the T-cell responses against an epitope within the myelin proteolipid protein (PLP) which together with the myelin basic protein (MBP) and the myelin oligodendrocyte glycoprotein (MOG) is one of the candidate autoantigens in MS (34).
Further sequencing and association analyses on larger sample sets are necessary to allow the exclusion mapping of MHC polymorphisms responsible for the MS association of the D6S1683 microsatellite in the telomeric portion of the HLA region. Irrespective of the nature of the DNA variation primarily responsible for the HLA telomeric association, these results are consistent with a report showing, in a Scadinavian sample set, a modifying effect of HLA class I markers on the association of the HLA class II DRB1*1501DQA1*0102DQB1*0602 haplotype (35). In another report the microsatellite D6S461, 10.7 cM telomeric of HLADRB1, showed significant evidence for LD with the disease, in three independent data sets, suggesting a gene effect in this region (6).
Taken together these observations support a complex multilocus, and in the case of DRB1DQB1, multiallelic model, of the MHC-encoded susceptibility to MS. However, in the absence of all the markers and haplotypes in the region, it is difficult to firmly establish the actual number and the relative importance of loci involved, and formally, it is still possible, although highly unlikely, that an untyped locus might explain all the associations detected in the MHC.
Overall, our results are fully consistent with the rat animal model in which specific MHC haplotypes determined the degree of disease susceptibility, clinical course and central nervous system pathology in a hierarchical and allele-specific way, with the major effect being mapped on the MHC class II subregion and with modifiers in the MHC (36). Interestingly, as we observed in our patients, as in the rat model, no single MHC class II haplotype was able to provide dominant protection (36).
These observations have some mechanistic implications. The observed lack of dominant protection might indicate that HLA class II presenting molecules are unlikely to play a major role in dominant tolerance mechanisms, such as the clonal deletion of autoreactive T cells in the thymus, if these mechanisms are operative at all in conferring a resistance from MS. In contrast, T1D HLA class II associations show a highly skewed distribution from susceptibility to dominant protection (37). Moreover, the multiple associations of MS with specific HLA class II haplotypes seen in humans might underlie the presence of distinct immunopathogenetic pathways involved in the presentation of different critical autoantigens (or different epitopes from the same autoantigen) to CD4 T cells. However, we cannot exclude the possibility that this apparent class II allelic heterogeneity is instead related to the presence of common structural motifs shared in a complex way by different HLADR, DQ, DP disease-associated molecules with one or very few primary autoantigenic epitopes. Mechanistic considerations aside, we can conclude that MS is primarily associated in humans, as it is in the rat animal model, to the MHC loci encoding the ß chains of the class II antigen presenting molecules and to additional unknown susceptibility gene polymorphism(s) in LD with the microsatellite D6S1683.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Subjects
The whole family data set consisted of 490 Sardinian MS families and included 444 simplex families (one affected child and both parents), 40 MS multiplex families (more than one affected sibling and their parents) and six vertical families with one parent and one affected offspring. 264 of these families have been previously analysed at the DRB1DQA1DQB1 loci (18). We also considered an independent sample set of 378 MS sporadic patients (total independent Sardinian patients equal to 868). 56.4% of the patients were from the province of Cagliari, 25.7% were from the province of Sassari, 11.4% were from the province of Nuoro and 6.5% were from the province of Oristano. Among the MS patients there were 593 females and 275 males (F/M ratio 2.16), mean age 35.4, range 1072 years. The diseases course was relapsing remitting in 73%, secondary progressive in 19% and primary progressive in 8% of patients. All index patients and affected first-degree relatives of index patients met the criteria of clinically or laboratory-supported definite MS (38). Patients were born and lived in Sardinia and were of Sardinian descent from at least three generations.
In the present study we have also considered the DRB1DQA1DQB1 frequencies of 631 Sardinian newborns that have been previously reported by Lampis et al. (22).
Microsatellites and sequence-specific oligonucleotide (SSO) typing
375 MS families were genotyped for 19 microsatellite markers (D6S291, D6S439, D6S1629, D6S1560, D6S2445, D6S2444, D6S2447, D3A, D6S273, 82-1, TNFc, TNFa, C1-2-A, D6S265, MOG51, D6S1683, D6S258, D6S306 and D6S2223) and for a number of SNPs defining allelic variation at 12 expressed genes (Tapasin, DPB1, DMB, TAP-1, TAP-2, DOB, DQB1, DQA1, DRB1, HSP70-1, HSP70-2 and HSP70HOM). An additional 115 MS families were genotyped at the loci peaking the regions of association (DPB1, DQB1, DQA1, DRB1 and D6S1683) (Results). The primer sequences for D6S291, D6S439, D6S1629, D6S1560, D6S2445, D6S2444, D6S2447, D6S273, C1-2-A, D6S265, D6S1683, D6S258 and D6S306, were obtained from Foissac and Cambon-Thomsen (39). The primer sequences for MOG51 were obtained from Malfroy et al. (40). Sequences for TNFa and TNFc were obtained from Udalova et al. (41). The primer sequence for D6S2223 was obtained from The Genome Database (GDB; http://www.gdb.org). Sequences for D3A and 82-1 were established by Hsieh et al. (42). Genotyping was performed by separating fluorescently tagged PCR products on a polyacrylamide gel using ABI 377 sequencer and the GeneScan 3.1, Genotyper 2.0 software (Perkin-Elmer, Applied Biosystems). PCR product standards consisting of the amplification product of two different standards for each marker were loaded on each gel for correct allele assignment. The two standards consisted of the Centre dEtude des Polymorphisms Human (CEPH) individual number 1347.02 and a pool of DNAs. The alleles at each microsatellite were given a numerical value (1, 2, 3, etc.) starting with the allele with the lowest number of base pairs. The physical map of the region with relative order, map position and distances between markers were obtained from the Sanger Centre (http://www.sanger.ac.uk/HGP/Chr6/MHC.shtml; A.Mungall, personal communication). Physical map position of the microsatellite markers are shown in Figure 1. Markers were PCR amplified and genotyped a second time when showing failures during the first round of amplifications. On average 75.9% of all parents were heterozygous for the microsatellites investigated, with outlying markers D6S2444, TNFc, D6S2223, D6S2445, showing heterozygosity lower than 60%. The average number of alleles per microsatellite was 11, but it was only three when we considered alleles with a parental frequency of at least 10%.
The expressed genes considered in this study were typed as follows. The polymorphic second exons of the HLADRB1, DQA1, DQB1, DPB1 genes were amplified and the amplified products were dot-blot analysed using primers and SSO probes described previously (4345). This includes characterization of 36, 9, 15 and 22 alleles for the DRB1, DQA1, DQB1 and DPB1 loci, respectively.
DMB (six alleles) polymorphisms were typed using primers and conditions previously described by Djilali-Saiah et al. (46). DOB (two alleles) was typed with amplification of the fourth exon and subsequent dot-blot analysis of the amplified products using 5'-GTGTCTAGTACAGATTCTG-3', 5'-CACTCCTCACAGGCTCAT-3' as PCR primers and 5'-GTGGGAATCATCATCCAG-3', 5'-GTGGGAATCGTCATCCAG-3' as (SSO) probes.
The Tapasin gene (two alleles) was typed by amplification of the fourth exon and subsequent dot-blot analysis of the amplified products using 5'-AAATGGGACCTTCTGGCTGC-3', 5'-AAGCTCCAGGGTGACCTGTC-3' as PCR primers and 5'-GGCTGCCTAGAGTTCAACCC-3', 5'-GGCTGCCTACAGTTCAACCC-3' as (SSO) probes. (J.Copeman, personal communication).
TAP-1 (four alleles) and TAP-2 (eight alleles) were typed according to Powis et al. (47).
HSP70-1 (two alleles), HSP71-2 (two alleles) and HSPHOM (two alleles) were typed according to Vinasco et al. (48).
Statistical analysis of the data
The degree of association of the various loci with MS was established using the ETDT (49). This test takes into account the transmission or non-transmission of alleles of a marker relative to the alleles of the marker present on the other parental chromosome. 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 of all allelic variation at each individual locus. Therefore, the ETDT could be considered a particularly suitable method for the analysis of a region showing a multiallelic association. The P-values obtained with the ETDT were corrected for the number of markers tested, 31.
In order to distinguish primary associations from those due to LD at the test loci, we used a variant of the ETDT, the CETDT (50). This test allows us to analyse the overall effect of one locus taking the association of other linked loci into account. The CETDT compares the transmission of haplotypes constructed from all the loci against the null hypothesis that all haplotypes identical at the conditioning loci have equal transmission probabilities. This conditional method appears particularly suitable in the presence of multiallelic association at the conditioned locus.
To establish which alleles and haplotypes were accounting for the main locus contributions detected with the ETDT and CETDT analyses, we performed an analysis of the relative association of the various markers using a case-control design. The control frequencies of the various alleles and haplotypes, referred to as affected family-based control (AFBAC), have been calculated as described by Thomson for single alleles (51). The distribution of the patient and control alleles and haplotypes, determined with a gene counting procedure, was then arranged in a 2 x 2 contingency table and tested by Fishers exact or Pearsons chi-squared test and the OR computed. This analysis was used to establish the relative association of alleles and haplotypes at loci with prior evidence of association (ETDT and CETDT results). For this reason, in this case we did not apply a correction for multiple testing. Moreover, in the case of particular DRB1DQB1 haplotypes and DPB1 alleles, the prior evidence of association was also strongly supported by specific prior evidence obtained in different studies from distinct populations (Discussion).
When more than one allele or haplotype is associated with the susceptibility to a complex trait, the calculation of ORs might be influenced by the frequency and penetrance of the other associated alleles (or haplotypes) at the same disease locus (i). In particular, it is difficult to distinguish whether the negatively associated alleles are such by default; that is, their reduced frequencies in patients are due to algebraic constraints (the sum of the frequencies must be 100%) driven by the positively associated alleles or whether they are genuinely negatively associated and hence indicative of active protection. To establish the RPE of the various alleles (or haplotypes) we have applied the RPE method (27). The RPE method calculates a P-value for each allele according to the formula
2 = (Oi Ei)2 / Ei, where Oi is the observed number of occurrence of allele i in patients and Ei is the number expected according to the control allele frequency. The total P-value for the data set is calculated by adding together the non-zero alleles. If the total P-value is less than the significance criteria (in our analysis P = 0.0001) then the alleles with the smallest P-value are sequentially removed from the data set and the expected frequencies are re-calculated as though that allele did not exist. This procedure is iterated until the total P-value becomes non-significant.
Haplotypes in MS families were established following the co-segregation of alleles within families and using computer program TDTPHASE written by F.Dudbridge and available at http://www-gene.cimr.cam.ac.uk/tdt/. Only certain haplotypes from parental genotype data, and in the absence of inter-crosses (that is when both parents are heterozygous for the same alleles), were considered in the analyses shown in this work. In all the association analyses performed, we selected the probands from all the families with more than one affected sibling. The DRB1DQB1 haplotypes in the 631 control newborns and in the 378 MS sporadic patients were assigned following the known patterns of LD in Caucasians and Sardinians (22,52). In the case of rare associations, the haplotypes were accepted only when the haplotype present on the other chromosome was well defined. Ambiguous assignments were resolved by excluding those individuals.
The LD patterns between the various marker loci measured in this study were calculated on the parental chromosomes from 375 MS families using a normalized disequilibrium (total D') multiallelic extension of Lewontins standardized measure of disequilibrium (53,54). The D' values range from 0 to 1, with 0 reflecting perfect independence between alleles at the two loci compared and 1 reflecting complete LD (53,55). The respective P-values were calculated using the Markov chain method described by Guo and Thompson (54) (available at http://anthropologie.unige.ch/arlequin). In all cases, 100 000 tables were explored. This P-value is reported with each pairwise comparison and was considered significant each time the ratio was <0.05. LD P-values were corrected for multiple testing using the step-down Holm-Sidak procedure as described in Lautenberger et al. (56). Each P-value is corrected by the formula Pc = 1 [1 P(uncorrected)]n where n is the number of uncorrected P-values for that locus less than or equal to the value to be corrected. For example if the eighth most significant P-value was 0.01 before correction, the corrected value would be Pc = 1 (1 0.01)8 = 0.0773.
| ACKNOWLEDGEMENTS |
|---|
We are grateful to the patients and their relatives for their kind co-operation in our study. We also wish to thank John Todd and Antonio Cao for continuous help, advise and support, Cesare Zavattari for writing a program that allows assignment of the allele sizes into their appropriate allele bins, Heather Cordell and Frank Dudbridge for statistical advice, James Copeman for information about tapasin, Rosanna Lampis for information about primers and PCR conditions to genotype the SNPs within the MHC region, Andrew Mungall (Sanger Center) for the establishment of the physical map. The work has been supported by Federazione Italiana Sclerosi Multipla and Dompé-Biotec grants. The DNA MS bank is supported by a grant of the Federazione Italiana Sclerosi Multipla.
| FOOTNOTES |
|---|
+ To whom correspondence should be addressed. Tel: +39 070 6095681; Fax: +39 070 6095558; Email: fcucca@mcweb.unica.it Correspondence may also be addressed to M.G.Marrosu. Tel: +39 070 6092806; Fax: +39 070 6092929; Email: gmarrosu@unica.it
| REFERENCES |
|---|
|
|
|---|
1 Withaker, J. and Mitchell, G. (1997) Clinical features of multiple sclerosis. In Raine, C., McFarland, H. and Tourtellotte, W. (eds), Multiple Sclerosis: Clinical and Pathogenetic Basis. Chapman and Hall, London, pp. 319.
2 Sadovnick, A.D. and Ebers, G.C. (1993) Epidemiology of multiple sclerosis: a critical overview. Can. J. Neurol. Sci., 20, 1729.[Web of Science][Medline]
3 Ebers, G.C., Sadovnick, A.D. and Risch, N.J. (1995) A genetic basis for familial aggregation in multiple sclerosis. Canadian Collaborative Study Group. Nature, 377, 150151.[Medline]
4 Robertson, N.P., Fraser, M., Deans, J., Clayton, D., Walker, N. and Compston, D.A. (1996) Age-adjusted recurrence risks for relatives of patients with multiple sclerosis. Brain, 119, 449455.
5 Carton, H., Vlietinck, R., Debruyne, J., De Keyser, J., DHooghe, M.B., Loos, R., Medaer, R., Truyen, L., Yee, I.M. and Sadovnick, A.D. (1997) Risks of multiple sclerosis in relatives of patients in Flanders, Belgium. J. Neurol. Neurosurg. Psychiatry, 62, 329333.
6 Ebers, G.C., Kukay, K., Bulman, D.E., Sadovnick, A.D., Rice, G., Anderson, C., Armstrong, H., Cousin, K., Bell, R.B., Hader, W. et al. (1996) A full genome search in multiple sclerosis. Nat. Genet., 13, 472476.[Web of Science][Medline]
7 Sawcer, S., Jones, H.B., Feakes, R., Gray, J., Smaldon, N., Chataway, J., Robertson, N., Clayton, D., Goodfellow, P.N. and Compston, A. (1996) A genome screen in multiple sclerosis reveals susceptibility loci on chromosome 6p21 and 17q22. Nat. Genet., 13, 464468.[Web of Science][Medline]
8 Haines, J.L., Ter-Minassian, M., Bazyk, A., Gusella, J.F., Kim, D.J., Terwedow, H., Pericak-Vance, M.A., Rimmler, J.B., Haynes, C.S., Roses, A.D. et al. (1996) A complete genomic screen for multiple sclerosis underscores a role for the major histocompatability complex. The Multiple Sclerosis Genetics Group. Nat. Genet., 13, 469471.[Web of Science][Medline]
9 The Transatlantic Multiple Sclerosis Genetics Cooperative (2001) A meta-analysis of genomic screens in multiple sclerosis. Mult. Scler., 7, 311.
10 Sawcer, S. and Goodfellow, P.N. (1998) Inheritance of susceptibility to multiple sclerosis. Curr. Opin. Immunol., 10, 697703.[Web of Science][Medline]
11 Kuokkanen, S., Gschwend, M., Rioux, J.D., Daly, M.J., Terwilliger, J.D., Tienari, P.J., Wikstrom, J., Palo, J., Stein, L.D., Hudson, T.J. et al. (1997) Genomewide scan of multiple sclerosis in Finnish multiplex families. Am. J. Hum. Genet., 61, 13791387.[Web of Science][Medline]
12 Hillert, J. and Olerup, O. (1993) Multiple sclerosis is associated with genes within or close to the HLADRDQ subregion on a normal DR15, DQ6, Dw2 haplotype. Neurology, 43, 163168.
13 Hauser, S.L., Fleischnick, E., Weiner, H.L., Marcus, D., Awdeh, Z., Yunis, E.J. and Alper, C.A. (1989) Extended major histocompatibility complex haplotypes in patients with multiple sclerosis. Neurology, 39, 275277.
14 Kwon, O.J., Karni, A., Israel, S., Brautbar, C., Amar, A., Meiner, Z., Abramsky, O. and Karussis, D. (1999) HLA class II susceptibility to multiple sclerosis among Ashkenazi and non-Ashkenazi Jews. Arch. Neurol., 56, 555560.
15 Saruhan-Direskeneli, G., Esin, S., Baykan-Kurt, B., Ornek, I., Vaughan, R. and Eraksoy, M. (1997) HLADR and DQ associations with multiple sclerosis in Turkey. Hum. Immunol., 55, 5965.[Medline]
16 Alvarado-de la Barrera, C., Zuniga-Ramos, J., Ruiz-Morales, J.A., Estanol, B., Granados, J. and Llorente, L. (2000) HLA class II genotypes in Mexican Mestizos with familial and nonfamilial multiple sclerosis. Neurology, 55, 18971900.
17 Marrosu, M.G., Murru, M.R., Costa, G., Cucca, F., Sotgiu, S., Rosati, G. and Muntoni, F. (1997) Multiple sclerosis in Sardinia is associated and in linkage disequilibrium with HLADR3 and DR4 alleles. Am. J. Hum. Genet., 61, 454457.[Web of Science][Medline]
18 Marrosu, M.G., Murru, M.R., Costa, G., Murru, R., Muntoni, F. and Cucca, F. (1998) DRB1DQA1DQB1 loci and multiple sclerosis predisposition in the Sardinian population. Hum. Mol. Genet., 7, 12351237.
19 Carrington, M. (1999) Recombination within the human MHC. Immunol. Rev., 167, 245256.[Web of Science][Medline]
20 Zavattari, P., Deidda, E., Whalen, M., Lampis, R., Mulargia, A., Loddo, M., Eaves, I., Mastio, G., Todd, J.A. and Cucca, F. (2000) Major factors influencing linkage disequilibrium by analysis of different chromosome regions in distinct populations: demography, chromosome recombination frequency and selection. Hum. Mol. Genet., 9, 29472957.
21 Granieri, E., Casetta, I., Govoni, V., Tola, M.R., Marchi, D., Murgia, S.B., Ticca, A., Pugliatti, M., Murgia, B. and Rosati, G. (2000) The increasing incidence and prevalence of MS in a Sardinian province. Neurology, 55, 842848.
22 Lampis, R., Morelli, L., Congia, M., Macis, M.D., Mulargia, A., Loddo, M., De Virgiliis, S., Marrosu, M.G., Todd, J.A. and Cucca, F. (2000) The inter-regional distribution of HLA class II haplotypes indicates the suitability of the Sardinian population for case-control association studies in complex diseases. Hum. Mol. Genet., 9, 29592965.
23 Lampis, R., Morelli, L., De Virgiliis, S., Congia, M. and Cucca, F. (2000) The distribution of HLA class II haplotypes reveals that the Sardinian population is genetically differentiated from the other Caucasian populations. Tissue Antigens, 56, 515521.[Web of Science][Medline]
24 Herr, M., Dudbridge, F., Zavattari, P., Cucca, F., Guja, C., March, R., Campbell, R.D., Barnett, A.H., Bain, S.C., Todd, J.A. et al. (2000) Evaluation of fine mapping strategies for a multifactorial disease locus: systematic linkage and association analysis of IDDM1 in the HLA region on chromosome 6p21. Hum. Mol. Genet., 9, 12911301.
25 Zavattari, P., Lampis, R., Mulargia, A., Loddo, M., Angius, E., Todd, J.A. and Cucca, F. (2000) Confirmation of the DRB1DQB1 loci as the major component of IDDM1 in the isolated founder population of Sardinia. Hum. Mol. Genet., 9, 29672972.
26 Jeffreys, A.J., Ritchie, A. and Neumann, R. (2000) High resolution analysis of haplotype diversity and meiotic crossover in the human TAP2 recombination hotspot. Hum. Mol. Genet., 9, 725733.
27 Payami, H., Joe, S., Farid, N.R., Stenszky, V., Chan, S.H., Yeo, P.P., Cheah, J.S. and Thomson, G. (1989) Relative predispositional effects (RPEs) of marker alleles with disease: HLADR alleles and Graves disease. Am. J. Hum. Genet., 45, 541546.[Web of Science][Medline]
28 Ligers, A., Dyment, D.A., Willer, C.J., Sadovnick, A.D., Ebers, G., Risch, N. and Hillert, J. (2001) Evidence of linkage with hladr in drb1*15-negative families with multiple sclerosis. Am. J. Hum. Genet., 69, 900903.[Web of Science][Medline]
29 Karni, A., Kohn, Y., Safirman, C., Abramsky, O., Barcellos, L., Oksenberg, J.R., Kahana, E., Karussis, D., Chapman, J. and Brautbar, C. (1999) Evidence for the genetic role of human leukocyte antigens in low frequency DRB1*1501 multiple sclerosis patients in Israel. Mult. Scler., 5, 410415.
30 Weinshenker, B.G., Santrach, P., Bissonet, A.S., McDonnell, S.K., Schaid, D., Moore, S.B. and Rodriguez, M. (1998) Major histocompatibility complex class II alleles and the course and outcome of MS: a population-based study. Neurology, 51, 742747.
31 Masterman, T., Ligers, A., Olsson, T., Andersson, M., Olerup, O. and Hillert, J. (2000) HLADR15 is associated with lower age at onset in multiple sclerosis. Ann. Neurol., 48, 211219.[Web of Science][Medline]
32 Fukazawa, T., Yamasaki, K., Ito, H., Kikuchi, S., Minohara, M., Horiuchi, I., Tsukishima, E., Sasaki, H., Hamada, T., Nishimura, Y. et al. (2000) Both the HLADPB1 and DRB1 alleles correlate with risk for multiple sclerosis in Japanese: clinical phenotypes and gender as important factors. Tissue Antigens, 55, 199205.[Web of Science][Medline]
33 Dekker, J.W., Easteal, S., Jakobsen, I.B., Gao, X., Stewart, G.J., Buhler, M.M., Hawkins, B.R., Higgins, D.A., Yu, Y.L. and Serjeantson, S.W. (1993) HLADPB1 alleles correlate with risk for multiple sclerosis in Caucasoid and Cantonese patients lacking the high-risk DQB1*0602 allele. Tissue Antigens, 41, 3136.[Web of Science][Medline]
34 Yu, M., Kinkel, R.P., Weinstock-Guttman, B., Cook, D.J. and Tuohy, V.K. (1998) HLADP: a class II restriction molecule involved in epitope spreading during the development of multiple sclerosis. Hum. Immunol., 59, 1524.[Web of Science][Medline]
35 Fogdell-Hahn, A., Ligers, A., Gronning, M., Hillert, J. and Olerup, O. (2000) Multiple sclerosis: a modifying influence of HLA class I genes in an HLA class II associated autoimmune disease. Tissue Antigens, 55, 140148.[Web of Science][Medline]
36 Weissert, R., Wallstrom, E., Storch, M.K., Stefferl, A., Lorentzen, J., Lassmann, H., Linington, C. and Olsson, T. (1998) MHC haplotype-dependent regulation of MOG-induced EAE in rats. J. Clin. Invest., 102, 12651273.[Web of Science][Medline]
37 Cucca, F., Lampis, R., Congia, M., Angius, E., Nutland, S., Bain, S.C., Barnett, A.H. and Todd, J.A. (2001) A correlation between the relative predisposition of MHC class II alleles to type 1 diabetes and the structure of their proteins. Hum. Mol. Genet., 10, 20252037
38 Poser, C.M., Paty, D.W., Scheinberg, L., McDonald, W.I., Davis, F.A., Ebers, G.C., Johnson, K.P., Sibley, W.A., Silberberg, D.H. and Tourtellotte, W.W. (1983) New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann. Neurol., 13, 227231.[Web of Science][Medline]
39 Foissac, A. and Cambon-Thomsen, A. (1998) Microsatellites in the HLA region: 1998 update. Tissue Antigens, 52, 318352.[Web of Science][Medline]
40 Malfroy, L., Viratelle, C., Coppin, H., Borot, N. and Roth, M.P. (1995) Polymorphic tri- and tetranucleotide repeats in exons 1 and 8 of the myelin oligodendrocyte glycoprotein (MOG) gene. Hum. Genet., 96, 737738.[Web of Science][Medline]
41 Udalova, I.A., Nedospasov, S.A., Webb, G.C., Chaplin, D.D. and Turetskaya, R.L. (1993) Highly informative typing of the human TNF locus using six adjacent polymorphic markers. Genomics, 16, 180186.[Web of Science][Medline]
42 Hsieh, S.L., March, R.E., Khanna, A., Cross, S.J. and Campbell, R.D. (1997) Mapping of 10 novel microsatellites in the MHC class III region: application to the study of autoimmune disease. J. Rheumatol., 24, 220222.[Web of Science][Medline]
43 Horn, G.T., Bugawan, T.L., Long, C.M. and Erlich, H.A. (1988) Allelic sequence variation of the HLADQ loci: relationship to serology and to insulin-dependent diabetes susceptibility. Proc. Natl Acad. Sci. USA, 85, 60126016.
44 Wordsworth, B.P., Allsopp, C.E., Young, R.P. and Bell, J.I. (1990) HLA-DR typing using DNA amplification by the polymerase chain reaction and sequential hybridization to sequence-specific oligonucleotide probes. Immunogenetics, 32, 413418.[Web of Science][Medline]
45 Bugawan, T.L., Begovich, A.B. and Erlich, H.A. (1991) Rapid HLADPB typing using enzymatically amplified DNA and nonradioactive sequence-specific oligonucleotide probes. Immunogenetics, 34, 413.[Web of Science][Medline]
46 Djilali-Saiah, I., Benini, V., Schmitz, J., Timsit, J., Assan, R., Boitard, C., Bach, J.F. and Caillat-Zucman, S. (1996) Absence of primary association between DM gene polymorphism and insulin-dependent diabetes mellitus or celiac disease. Hum. Immunol., 49, 2227.[Medline]
47 Powis, S.H., Tonks, S., Mockridge, I., Kelly, A.P., Bodmer, J.G. and Trowsdale, J. (1993) Alleles and haplotypes of the MHC-encoded ABC transporters TAP1 and TAP2. Immunogenetics, 37, 373380.[Web of Science][Medline]
48 Vinasco, J., Beraun, Y., Nieto, A., Fraile, A., Pareja, E., Mataran, L. and Martin, J. (1997) Heat shock protein 70 gene polymorphisms in rheumatoid arthritis. Tissue Antigens, 50, 7173.[Web of Science][Medline]
49 Sham, P.C. and Curtis, D. (1995) An extended transmission/disequilibrium test (TDT) for multi-allele marker loci. Ann. Hum. Genet., 59, 323336.[Web of Science][Medline]
50 Koeleman, B.P., Dudbridge, F., Cordell, H.J. and Todd, J.A. (2000) Adaptation of the extended transmission/disequilibrium test to distinguish disease associations of multiple loci: the Conditional Extended Transmission/Disequilibrium Test. Ann. Hum. Genet., 64, 207213.[Web of Science][Medline]
51 Thomson, G. (1995) Mapping disease genes: family-based association studies. Am. J. Hum. Genet., 57, 487498.[Web of Science][Medline]
52 Fernandez-Vina, M.A., Gao, X.J., Moraes, M.E., Moraes, J.R., Salatiel, I., Miller, S., Tsai, J., Sun, Y.P., An, J.B., Layrisse, Z. et al. (1991) Alleles at four HLA class II loci determined by oligonucleotide hybridization and their associations in five ethnic groups. Immunogenetics, 34, 299312.[Web of Science][Medline]
53 Lewontin, R.C. (1988) On measures of gametic disequilibrium. Genetics, 120, 849852.
54 Guo, S.W. and Thompson, E.A. (1992) Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics, 48, 361372.[Web of Science][Medline]
55 Hedrick, P.W. (1987) Gametic disequilibrium measures: proceed with caution. Genetics, 117, 331341.
56 Lautenberger, J.A., Stephens, J.C., OBrien, S.J. and Smith, M.W. (2000) Significant admixture linkage disequilibrium across 30 cM around the FY locus in African Americans. Am. J. Hum. Genet., 66, 969978.[Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
E Cocco, C Sardu, R Murru, J Frau, L Lorefice, E Mamusa, P Contu, and M. Marrosu Multiple sclerosis risk: interaction between human leukocyte antigen and the environment in Sardinian population Multiple Sclerosis, September 1, 2009; 15(9): 1030 - 1036. [Abstract] [PDF] |
||||
![]() |
M. R. Lincoln, S. V. Ramagopalan, M. J. Chao, B. M. Herrera, G. C. DeLuca, S.-M. Orton, D. A. Dyment, A D. Sadovnick, and G. C. Ebers Epistasis among HLA-DRB1, HLA-DQA1, and HLA-DQB1 loci determines multiple sclerosis susceptibility PNAS, May 5, 2009; 106(18): 7542 - 7547. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Sawcer The complex genetics of multiple sclerosis: pitfalls and prospects Brain, December 1, 2008; 131(12): 3118 - 3131. [Abstract] [Full Text] [PDF] |
||||
![]() |
D A Dyment, M Z Cader, B M Herrera, S V Ramagopalan, S M Orton, M Chao, C J Willer, A D Sadovnick, N Risch, and G C Ebers A genome scan in a single pedigree with a high prevalence of multiple sclerosis J. Neurol. Neurosurg. Psychiatry, February 1, 2008; 79(2): 158 - 162. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Chao, M. C.N.M. Barnardo, G.-z. Lui, M. R. Lincoln, S. V. Ramagopalan, B. M. Herrera, D. A. Dyment, A. D. Sadovnick, and G. C. Ebers Transmission of class I/II multi-locus MHC haplotypes and multiple sclerosis susceptibility: accounting for linkage disequilibrium Hum. Mol. Genet., August 15, 2007; 16(16): 1951 - 1958. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Schmidt, D. Williamson, and A. Ashley-Koch HLA-DR15 Haplotype and Multiple Sclerosis: A HuGE Review Am. J. Epidemiol., May 15, 2007; 165(10): 1097 - 1109. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. F. Barcellos, S. Sawcer, P. P. Ramsay, S. E. Baranzini, G. Thomson, F. Briggs, B. C.A. Cree, A. B. Begovich, P. Villoslada, X. Montalban, et al. Heterogeneity at the HLA-DRB1 locus and risk for multiple sclerosis Hum. Mol. Genet., September 15, 2006; 15(18): 2813 - 2824. [Abstract] [Full Text] [PDF] |
||||
![]() |
T Fukazawa, S Kikuchi, R Miyagishi, Y Miyazaki, I Yabe, T Hamada, and H Sasaki HLA-DPB1*0501 is not uniquely associated with opticospinal multiple sclerosis in Japanese patients. Important role of DPB1*0301 Multiple Sclerosis, February 1, 2006; 12(1): 19 - 23. [Abstract] [PDF] |
||||
![]() |
D. A. Dyment, B. M. Herrera, M. Z. Cader, C. J. Willer, M. R. Lincoln, A. D. Sadovnick, N. Risch, and G. C. Ebers Complex interactions among MHC haplotypes in multiple sclerosis: susceptibility and resistance Hum. Mol. Genet., July 15, 2005; 14(14): 2019 - 2026. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Cocco, A. Sotgiu, G. Costa, M. R. Murru, C. Mancosu, R. Murru, M. Lai, P. Contu, and M. G. Marrosu HLA-DR,DQ and APOE genotypes and gender influence in Sardinian primary progressive MS Neurology, February 8, 2005; 64(3): 564 - 566. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Brassat, G. Salemi, L. F. Barcellos, G. McNeill, P. Proia, S. L. Hauser, J. R. Oksenberg, and G. Savettieri The HLA locus and multiple sclerosis in Sicily Neurology, January 25, 2005; 64(2): 361 - 363. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. G. Marrosu, C. Motzo, R. Murru, R. Lampis, G. Costa, P. Zavattari, D. Contu, E. Fadda, E. Cocco, and F. Cucca The co-inheritance of type 1 diabetes and multiple sclerosis in Sardinia cannot be explained by genotype variation in the HLA region alone Hum. Mol. Genet., December 1, 2004; 13(23): 2919 - 2924. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Teuscher, J. Y. Bunn, P. D. Fillmore, R. J. Butterfield, J. F. Zachary, and E. P. Blankenhorn Gender, Age, and Season at Immunization Uniquely Influence the Genetic Control of Susceptibility to Histopathological Lesions and Clinical Signs of Experimental Allergic Encephalomyelitis: Implications for the Genetics of Multiple Sclerosis Am. J. Pathol., November 1, 2004; 165(5): 1593 - 1602. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. G. Marrosu, C. Sardu, E. Cocco, G. Costa, M. R. Murru, C. Mancosu, R. Murru, M. Lai, and P. Contu Bias in parental transmission of the HLA-DR3 allele in Sardinian multiple sclerosis Neurology, September 28, 2004; 63(6): 1084 - 1086. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Dyment, A. D. Sadovnick, C. J. Willer, H. Armstrong, Z. M. Cader, S. Wiltshire, B. Kalman, N. Risch, and G. C. Ebers An extended genome scan in 442 Canadian multiple sclerosis-affected sibships: a report from the Canadian Collaborative Study Group Hum. Mol. Genet., May 15, 2004; 13(10): 1005 - 1015. [Abstract] [Full Text] [PDF] |
||||
![]() |
F Capon, I K Toal, J C Evans, M H Allen, S Patel, D Tillman, D Burden, J N W N Barker, and R C Trembath Haplotype analysis of distantly related populations implicates corneodesmosin in psoriasis susceptibility J. Med. Genet., June 1, 2003; 40(6): 447 - 452. [Full Text] [PDF] |
||||
![]() |
B. A. Rybicki, M. J. Maliarik, L. M. Poisson, R. Sheffer, K. M. Chen, M. Major, G. A. Chase, and M. C. Iannuzzi The Major Histocompatibility Complex Gene Region and Sarcoidosis Susceptibility in African Americans Am. J. Respir. Crit. Care Med., February 1, 2003; 167(3): 444 - 449. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. G. De la Concha, M. Fernandez-Arquero, L. Gual, P. Vigil, A. Martinez, E. Urcelay, A. Ferreira, M. C. Garcia-Rodriguez, and G. Fontan MHC Susceptibility Genes to IgA Deficiency Are Located in Different Regions on Different HLA Haplotypes J. Immunol., October 15, 2002; 169(8): 4637 - 4643. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||












