Glenys Thomson

Mapping Multiple Sclerosis Susceptibility to the HLA-DR Locus in African Americans

An underlying complex genetic susceptibility exists in multiple sclerosis (MS), and an association with the HLA-DRB1*1501-DQB1*0602 haplotype has been repeatedly demonstrated in high-risk (northern European) populations. It is unknown whether the effect is explained by the HLA-DRB1 or the HLA-DQB1 gene within the susceptibility haplotype, which are in strong linkage disequilibrium (LD). African populations are characterized by greater haplotypic diversity and distinct patterns of LD compared with northern Europeans. To better localize the HLA gene responsible for MS susceptibility, case-control and family-based association studies were performed for DRB1 and DQB1 loci in a large and well-characterized African American data set. A selective association with HLA-DRB1*15 was revealed, indicating a primary role for the DRB1 locus in MS independent of DQB1*0602. This finding is unlikely to be solely explained by admixture, since a substantial proportion of the susceptibility chromosomes from African American patients with MS displayed haplotypes consistent with an African origin.

How selection shapes variation of the human major histocompatibility complex: a review

The nature of polymorphism and molecular sequence variation in the genes of the human major histocompatibility complex (MHC) provides strong support for the idea that these genes are under selection. With the understanding that selection shapes MHC variation new questions have become the focus of study. What is the mode of selection that accounts for MHC polymorphism? Is variation maintained by pathogen pressure or by reproductive mechanisms? Discerning between these requires drawing on information from studies on association between HLA genes and infectious diseases, reproductive success and mating preferences relative to HLA genotypes, and theoretical studies that compare the outcomes of different selection regimes. The pattern that has emerged suggests that several types of selection are plausible for the maintenance of HLA polymorphism.

PTPN22 Genetic Variation: Evidence for Multiple Variants Associated with Rheumatoid Arthritis

The minor allele of the R620W missense single-nucleotide polymorphism (SNP) (rs2476601) in the hematopoietic-specific protein tyrosine phosphatase gene, PTPN22, has been associated with multiple autoimmune diseases, including rheumatoid arthritis (RA). These genetic data, combined with biochemical evidence that this SNP affects PTPN22 function, suggest that this phosphatase is a key regulator of autoimmunity. To determine whether other genetic variants in PTPN22 contribute to the development of RA, we sequenced the coding regions of this gene in 48 white North American patients with RA and identified 15 previously unreported SNPs, including 2 coding SNPs in the catalytic domain. We then genotyped 37 SNPs in or near PTPN22 in 475 patients with RA and 475 individually matched controls (sample set 1) and selected a subset of markers for replication in an additional 661 patients with RA and 1,322 individually matched controls (sample set 2). Analyses of these results predict 10 common (frequency >1%) PTPN22 haplotypes in white North Americans. The sole haplotype found to carry the previously identified W620 risk allele was strongly associated with disease in both sample sets, whereas another haplotype, identical at all other SNPs but carrying the R620 allele, showed no association. R620W, however, does not fully explain the association between PTPN22 and RA, since significant differences between cases and controls persisted in both sample sets after the haplotype data were stratified by R620W. Additional analyses identified two SNPs on a single common haplotype that are associated with RA independent of R620W, suggesting that R620W and at least one additional variant in the PTPN22 gene region influence RA susceptibility.

Association mapping of disease loci, by use of a pooled DNA genomic screen.

Genomic screening to map disease loci by association requires automation, pooling of DNA samples, and 3,000-6,000 highly polymorphic, evenly spaced microsatellite markers. Case-control samples can be used in an initial screen, followed by family-based data to confirm marker associations. Association mapping is relevant to genetic studies of complex diseases in which linkage analysis may be less effective and to cases in which multigenerational data are difficult to obtain, including rare or late-onset conditions and infectious diseases. The method can also be used effectively to follow up and confirm regions identified in linkage studies or to investigate candidate disease loci. Study designs can incorporate disease heterogeneity and interaction effects by appropriate subdivision of samples before screening. Here we report use of pooled DNA amplifications-the accurate determination of marker-disease associations for both case-control and nuclear family-based data-including application of correction methods for stutter artifact and preferential amplification. These issues, combined with a discussion of both statistical power and experimental design to define the necessary requirements for detecting of disease loci while virtually eliminating false positives, suggest the feasibility and efficiency of association mapping using pooled DNA screening.

Balancing selection and heterogeneity across the classical human leukocyte antigen loci: a meta-analytic review of 497 population studies

This paper presents a meta-analysis of high-resolution human leukocyte antigen (HLA) allele frequency data describing 497 population samples. Most of the datasets were compiled from studies published in eight journals from 1990 to 2007; additional datasets came from the International Histocompatibility Workshops and from the AlleleFrequencies.net database. In all, these data represent approximately 66,800 individuals from throughout the world, providing an opportunity to observe trends that may not have been evident at the time the data were originally analyzed, especially with regard to the relative importance of balancing selection among the HLA loci. Population genetic measures of allele frequency distributions were summarized across populations by locus and geographic region. A role for balancing selection maintaining much of HLA variation was confirmed. Further, the breadth of this meta-analysis allowed the ranking of the HLA loci, with DQA1 and HLA-C showing the strongest balancing selection and DPB1 being compatible with neutrality. Comparisons of the allelic spectra reported by studies since 1990 indicate that most of the HLA alleles identified since 2000 are very-low-frequency alleles. The literature-based allele-count data, as well as maps summarizing the geographic distributions for each allele, are available online.

An integrated haplotype map of the human major histocompatibility complex.

Numerous studies have clearly indicated a role for the major histocompatibility complex (MHC) in susceptibility to autoimmune diseases. Such studies have focused on the genetic variation of a small number of classical human-leukocyte-antigen (HLA) genes in the region. Although these genes represent good candidates, given their immunological roles, linkage disequilibrium (LD) surrounding these genes has made it difficult to rule out neighboring genes, many with immune function, as influencing disease susceptibility. It is likely that a comprehensive analysis of the patterns of LD and variation, by using a high-density map of single-nucleotide polymorphisms (SNPs), would enable a greater understanding of the nature of the observed associations, as well as lead to the identification of causal variation. We present herein an initial analysis of this region, using 201 SNPs, nine classical HLA loci, two TAP genes, and 18 microsatellites. This analysis suggests that LD and variation in the MHC, aside from the classical HLA loci, are essentially no different from those in the rest of the genome. Furthermore, these data show that multi-SNP haplotypes will likely be a valuable means for refining association signals in this region.

PyPop: a software framework for population genomics: analyzing large-scale multi-locus genotype data.

Software to analyze multi-locus genotype data for entire populations is useful for estimating haplotype frequencies, deviation from Hardy-Weinberg equilibrium and patterns of linkage disequilibrium. These statistical results are important to both those interested in human genome variation and disease predisposition as well as evolutionary genetics. As part of the 13th International Histocompatibility and Immunogenetics Working Group (IHWG), we have developed a software framework (PyPop). The primary novelty of this package is that it allows integration of statistics across large numbers of data-sets by heavily utilizing the XML file format and the R statistical package to view graphical output, while retaining the ability to inter-operate with existing software. Largely developed to address human population data, it can, however, be used for population based data for any organism. We tested our software on the data from the 13th IHWG which involved data sets from at least 50 laboratories each of up to 1000 individuals with 9 MHC loci (both class I and class II) and found that it scales to large numbers of data sets well.

A review of theoretical aspects of HLA and disease associations

Abstract A number of antigens (alleles) of the human histocompatibility (HLA) system have been shown to be associated with a wide range of chronic diseases. Recently a great deal of interest has been shown in many of the theoretical aspects of HLA disease associations. These range from statistical problems of detecting significant associations to the development of models to determine the mode of inheritance of the disease. These models, as well as other features of HLA disease associations, will be reviewed in this paper.

The HLA Class I A Locus Affects Susceptibility to Type 1 Diabetes

Human leukocyte antigen A (HLA-A) genotypes were determined for samples from 283 multiplex, Caucasian, type 1 diabetes families from the Human Biological Data Interchange (HBDI) using an immobilized probe assay. Distribution of HLA-A alleles transmitted to patients was significantly different from that in affected family-based controls (AFBAC) (p = 0.004). Transmission disequilibrium test (TDT) analysis revealed differential transmission of several HLA-A alleles from parents to affected offspring. HLA class II DRB1 and DQB1 loci were also typed, allowing assignment of HLA-A alleles to haplotypes and calculation of linkage disequilibrium values. Some of the apparent effects of HLA-A alleles on type 1 diabetes susceptibility were attributable to linkage disequilibrium with DR and DQ alleles, although others were not. The differences in frequencies between patients and controls of alleles A*0101, A*2402, and A*3002 could not be explained by linkage disequilibrium alone. Our results suggest an important role for class I antigens in modulating susceptibility to type 1 diabetes.

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10−16). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10−12), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10−13. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.