The human Major Histocompatibility Complex as a paradigm in genomics research

Host genetics and viral load in primary HIV-1 infection: clear evidence for gene by sex interactions

Research in the past two decades has generated unequivocal evidence that host genetic variations substantially account for the heterogeneous outcomes following human immunodeficiency virus type 1 (HIV-1) infection. In particular, genes encoding human leukocyte antigens (HLA) have various alleles, haplotypes, or specific motifs that can dictate the set-point (a relatively steady state) of plasma viral load (VL), although rapid viral evolution driven by innate and acquired immune responses can obscure the long-term relationships between HLA genotypes and HIV-1-related outcomes. In our analyses of VL data from 521 recent HIV-1 seroconverters enrolled from eastern and southern Africa, HLA-A*03:01 was strongly and persistently associated with low VL in women (frequency = 11.3 %, P < 0.0001) but not in men (frequency = 7.7 %, P = 0.66). This novel sex by HLA interaction (P = 0.003, q = 0.090) did not extend to other frequent HLA class I alleles (n = 34), although HLA-C*18:01 also showed a weak association with low VL in women only (frequency = 9.3 %, P = 0.042, q > 0.50). In a reduced multivariable model, age, sex, geography (clinical sites), previously identified HLA factors (HLA-B*18, B*45, B*53, and B*57), and the interaction term for female sex and HLA-A*03:01 collectively explained 17.0 % of the overall variance in geometric mean VL over a 3-year follow-up period (P < 0.0001). Multiple sensitivity analyses of longitudinal and cross-sectional VL data yielded consistent results. These findings can serve as a proof of principle that the gap of "missing heritability" in quantitative genetics can be partially bridged by a systematic evaluation of sex-specific associations.

Associations of HLA alleles with specific language impairment

Background

Human leukocyte antigen (HLA) loci have been implicated in several neurodevelopmental disorders in which language is affected. However, to date, no studies have investigated the possible involvement of HLA loci in specific language impairment (SLI), a disorder that is defined primarily upon unexpected language impairment. We report association analyses of single-nucleotide polymorphisms (SNPs) and HLA types in a cohort of individuals affected by language impairment.

Methods

We perform quantitative association analyses of three linguistic measures and case-control association analyses using both SNP data and imputed HLA types.

Results

Quantitative association analyses of imputed HLA types suggested a role for the HLA-A locus in susceptibility to SLI. HLA-A A1 was associated with a measure of short-term memory (P = 0.004) and A3 with expressive language ability (P = 0.006). Parent-of-origin effects were found between HLA-B B8 and HLA-DQA1*0501 and receptive language. These alleles have a negative correlation with receptive language ability when inherited from the mother (P = 0.021, P = 0.034, respectively) but are positively correlated with the same trait when paternally inherited (P = 0.013, P = 0.029, respectively). Finally, case control analyses using imputed HLA types indicated that the DR10 allele of HLA-DRB1 was more frequent in individuals with SLI than population controls (P = 0.004, relative risk = 2.575), as has been reported for individuals with attention deficit hyperactivity disorder (ADHD).

Conclusion

These preliminary data provide an intriguing link to those described by previous studies of other neurodevelopmental disorders and suggest a possible role for HLA loci in language disorders.

The association of HLA-DRB1 alleles and drug use with HIV infection in a Chinese Han Cohort

Objective

To investigate the relationship between the polymorphism of human leukocyte antigen (HLA)-DRB1 and the susceptibility and repellency of drug use combined with HIV infection in Chinese.

Methods

A total of 213 unrelated healthy people, 41 HIV-infected drug users, 24 HIV-uninfected drug users, and 64 HIV-infected non-drug users were recruited. Their HLA-DRB1 allele frequencies were analyzed by PCR-SSP and allele distribution was analyzed.

Results

Compared with healthy controls, in drug users, the frequencies of HLA-DRB1 *0401-041, *1001 were significantly higher; in HIV-infected patients, the frequencies of HLA-DRB1 *0101-0103, *0401-0411, *1001 were significantly higher, while the frequencies of DRB1 *1501-1502, *1101-1105, *1301-1302, DRB4, DRB5 were significantly lower; in HIV-infected drug users, the frequencies of HLA-DRB1 *0101-0103, *0401-0411, *0801-0806, *1001, *1401/1404/1405 were significantly higher, while the frequencies of DRB1 *1301/1302, 1501-1502, DRB5 were significantly lower.

Conclusion

There is close relationship between the polymorphism of HLA-DRB1 alleles and drug use with HIV infection, which plays an important role in elucidating the pathogenesis and providing the basis for therapeutics and prophylaxis of patients with drug use and HIV infection.

Effects of highly conserved major histocompatibility complex (MHC) extended haplotypes on iron and low CD8+ T lymphocyte phenotypes in HFE C282Y homozygous hemochromatosis patients from three geographically distant areas

Hereditary Hemochromatosis (HH) is a recessively inherited disorder of iron overload occurring commonly in subjects homozygous for the C282Y mutation in HFE gene localized on chromosome 6p21.3 in linkage disequilibrium with the human leukocyte antigen (HLA)-A locus. Although its genetic homogeneity, the phenotypic expression is variable suggesting the presence of modifying factors. One such genetic factor, a SNP microhaplotype named A-A-T, was recently found to be associated with a more severe phenotype and also with low CD8⁺T-lymphocyte numbers. The present study aimed to test whether the predictive value of the A-A-T microhaplotype remained in other population settings. In this study of 304 HH patients from 3 geographically distant populations (Porto, Portugal 65; Alabama, USA 57; Nord-Trøndelag, Norway 182), the extended haplotypes involving A-A-T were studied in 608 chromosomes and the CD8⁺ T-lymphocyte numbers were determined in all subjects. Patients from Porto had a more severe phenotype than those from other settings. Patients with A-A-T seemed on average to have greater iron stores (p = 0.021), but significant differences were not confirmed in the 3 separate populations. Low CD8⁺ T-lymphocytes were associated with HLA-A*03-A-A-T in Porto and Alabama patients but not in the greater series from Nord-Trøndelag. Although A-A-T may signal a more severe iron phenotype, this study was unable to prove such an association in all population settings, precluding its use as a universal predictive marker of iron overload in HH. Interestingly, the association between A-A-T and CD8⁺ T-lymphocytes, which was confirmed in Porto and Alabama patients, was not observed in Nord-Trøndelag patients, showing that common HLA haplotypes like A*01–B*08 or A*03–B*07 segregating with HFE/C282Y in the three populations may carry different messages. These findings further strengthen the relevance of HH as a good disease model to search for novel candidate loci associated with the genetic transmission of CD8⁺ T-lymphocyte numbers.

Major histocompatibility complex genomics and human disease

Over several decades, various forms of genomic analysis of the human major histocompatibility complex (MHC) have been extremely successful in picking up many disease associations. This is to be expected, as the MHC region is one of the most gene-dense and polymorphic stretches of human DNA. It also encodes proteins critical to immunity, including several controlling antigen processing and presentation. Single-nucleotide polymorphism genotyping and human leukocyte antigen (HLA) imputation now permit the screening of large sample sets, a technique further facilitated by high-throughput sequencing. These methods promise to yield more precise contributions of MHC variants to disease. However, interpretation of MHC-disease associations in terms of the functions of variants has been problematic. Most studies confirm the paramount importance of class I and class II molecules, which are key to resistance to infection. Infection is likely driving the extreme variation of these genes across the human population, but this has been difficult to demonstrate. In contrast, many associations with autoimmune conditions have been shown to be specific to certain class I and class II alleles. Interestingly, conditions other than infections and autoimmunity are also associated with the MHC, including some cancers and neuropathies. These associations could be indirect, owing, for example, to the infectious history of a particular individual and selective pressures operating at the population level.

Molecular characterization, balancing selection, and genomic organization of the tree shrew (Tupaia belangeri) MHC class I gene

The major histocompatibility complex (MHC) class I genes play a pivotal role in the adaptive immune response among vertebrates. Accordingly, in numerous mammals the genomic structure and molecular characterization of MHC class I genes have been thoroughly investigated. To date, however, little is known about these genes in tree shrews, despite the increasingly popularity of its usage as an animal model. To address this deficiency, we analyzed the structure and characteristic of the tree shrew MHC class I genes (Tube-MHC I) and performed a comparative gene analysis of the tree shrew and other mammal species. We found that the full-length cDNA sequence of the tree shrew MHC class I is 1074bp in length. The deduced peptide is composed of 357 amino acids containing a leader peptide, an α1 and α2 domain, an α3 domain, a transmembrane domain and a cytoplasmic domain. Among these peptides, the cysteines, CD8(+) interaction and N-glycosylation sites are all well conserved. Furthermore, the genomic sequence of the tree shrew MHC class I gene was identified to be 3180bp in length, containing 8 exons and 7 introns. In 21 MHC class I sequences, we conducted an extensive study of nucleotide substitutions. The results indicated that in the peptide binding region (PBR) the rate of non-synonymous substitutions (dN) to synonymous substitutions (dS) was greater than 1, suggesting balancing selection at the PBR. These findings provide valuable contributions in furthering our understanding of the structure, molecular polymorphism, and function of the MHC class I genes in tree shrews, further improving their utility as an animal model in biomedical research.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Genomic modulators of the immune response

Our understanding of immunity has historically been informed by studying heritable mutations in both the adaptive and innate immune responses, including primary immunodeficiency and autoimmune diseases. Recent advances achieved through the application of genomic and epigenomic approaches are reshaping the study of immune dysfunction and opening up new avenues for therapeutic interventions. Moreover, applying genomic techniques to resolve functionally important genetic variation between individuals is providing new insights into immune function in health. This review describes progress in the study of rare variants and primary immunodeficiency diseases arising from whole-exome sequencing (WES), and discusses the application, success, and challenges of applying genome-wide association studies (GWAS) to disorders of immune function and how they may inform more rational use of therapeutics. In addition, the application of expression quantitative-trait mapping to immune phenotypes, progress in understanding MHC disease associations, and insights into epigenetic mechanisms at the interface of immunity and the environment are reviewed.

Personal receptor repertoires: olfaction as a model

Background

Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires.

Results

By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively “resurrected” from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs.

Conclusions

These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.

Meiotic recombinations within major histocompatibility complex of human embryos

We studied the rate and pattern of recombinations within the extended major histocompatibility complex (MHC) locus of the human embryos obtained during preimplantation genetic diagnosis (PGD) for HLA compatibility. Recombinant allele frequency was on average 5.33 %, and recombination rate was 0.44 cM/Mb in the 12.2 Mb of the extended MHC locus. Recombination rate varied up to 14-fold (0.19-2.73 cM/Mb) between cases, and maternal recombination rate was on average 3.8 times higher than paternal alleles. More than 69 % of the recombination hot spots were clustered within the extended class II region where the recombination rate was 5.4 times more than that in extended class I region. These findings indicate the potential of PGD to study the mechanisms of linkage disequilibrium within MHC locus of human embryos, demonstrate the recombination characteristics within extended MHC loci of human embryos in comparison to sperm and family studies, and point to the significance of design and interpretation of PGD for HLA compatibility to avoid misdiagnosis because of meiotic recombinations.

Expression QTL analysis of top loci from GWAS meta-analysis highlights additional schizophrenia candidate genes

There is genetic evidence that schizophrenia is a polygenic disorder with a large number of loci of small effect on disease susceptibility. Genome-wide association studies (GWASs) of schizophrenia have had limited success, with the best finding at the MHC locus at chromosome 6p. A recent effort of the Psychiatric GWAS consortium (PGC) yielded five novel loci for schizophrenia. In this study, we aim to highlight additional schizophrenia susceptibility loci from the PGC study by combining the top association findings from the discovery stage (9394 schizophrenia cases and 12 462 controls) with expression QTLs (eQTLs) and differential gene expression in whole blood of schizophrenia patients and controls. We examined the 6192 single-nucleotide polymorphisms (SNPs) with significance threshold at P<0.001. eQTLs were calculated for these SNPs in a sample of healthy controls (n=437). The transcripts significantly regulated by the top SNPs from the GWAS meta-analysis were subsequently tested for differential expression in an independent set of schizophrenia cases and controls (n=202). After correction for multiple testing, the eQTL analysis yielded 40 significant cis-acting effects of the SNPs. Seven of these transcripts show differential expression between cases and controls. Of these, the effect of three genes (RNF5, TRIM26 and HLA-DRB3) coincided with the direction expected from meta-analysis findings and were all located within the MHC region. Our results identify new genes of interest and highlight again the involvement of the MHC region in schizophrenia susceptibility.

Profiles of epigenetic histone post-translational modifications at type 1 diabetes susceptible genes

Both genetic and environmental factors are implicated in type 1 diabetes (T1D). Because environmental factors can trigger epigenetic changes, we hypothesized that variations in histone post-translational modifications (PTMs) at the promoter/enhancer regions of T1D susceptible genes may be associated with T1D. We therefore evaluated histone PTM variations at known T1D susceptible genes in blood cells from T1D patients versus healthy nondiabetic controls, and explored their connections to T1D. We used the chromatin immunoprecipitation-linked to microarray approach to profile key histone PTMs, including H3-lysine 4 trimethylation (H3K4me3), H3K27me3, H3K9me3, H3K9 acetylation (H3K9Ac), and H4K16Ac at genes within the T1D susceptible loci in lymphocytes, and H3K4me3, H3K9me2, H3K9Ac, and H4K16Ac at the insulin-dependent diabetes mellitus 1 region in monocytes of T1D patients and healthy controls separately. We screened for potential variations in histone PTMs using computational methods to compare datasets from T1D and controls. Interestingly, we observed marked variations in H3K9Ac levels at the upstream regions of HLA-DRB1 and HLA-DQB1 within the insulin-dependent diabetes mellitus 1 locus in T1D monocytes relative to controls. Additional experiments with THP-1 monocytes demonstrated increased expression of HLA-DRB1 and HLA-DQB1 in response to interferon-γ and TNF-α treatment that were accompanied by changes in H3K9Ac at the same promoter regions as that seen in the patient monocytes. These results suggest that the H3K9Ac status of HLA-DRB1 and HLA-DQB1, two genes highly associated with T1D, may be relevant to their regulation and transcriptional response toward external stimuli. Thus, the promoter/enhancer architecture and chromatin status of key susceptible loci could be important determinants in their functional association to T1D susceptibility.

Living with the enemy or uninvited guests: functional genomics approaches to investigating host resistance or tolerance traits to a protozoan parasite, Theileria annulata, in cattle

Many breeds of cattle with long histories of living in areas of endemic disease have evolved mechanisms that enable them to co-exist with specific pathogens. Understanding the genes that control tolerance and resistance could provide new strategies to improve the health and welfare of livestock. Around one sixth of the world cattle population is estimated to be at risk from one of the most debilitating tick-borne diseases of cattle, caused by the protozoan parasite, Theileria annulata. The parasite mainly infects cells of the myeloid lineage which are also the main producers of inflammatory cytokines. If an infectious or inflammatory insult is sufficiently great, inflammatory cytokines produced by macrophages enter the circulation and induce an acute phase proteins (APP) response. The Bos taurus Holstein breed produces higher and more prolonged levels of inflammatory cytokine induced APP than the Bos indicus Sahiwal breed in response to experimental infection with T. annulata. The Sahiwal exhibits significantly less pathology and survives infection, unlike the Holstein breed. Therefore, we hypothesised that the causal genes were likely to be expressed in macrophages and control the production of inflammatory cytokines. A functional genomics approach revealed that the transcriptome profile of the B. taurus macrophages was more associated with an inflammatory programme than the B. indicus macrophages. In particular the most differentially expressed gene was a member of the signal regulatory protein (SIRP) family. These are mainly expressed on myeloid cell surfaces and control inflammatory responses. Other differentially expressed genes included bovine major histocompatibility complex (MHC) (BoLA) class II genes, particularly BoLA DQ, and transforming growth factor (TGF)B2. We are now exploring whether sequence and functional differences in the bovine SIRP family may underlie the resistance or tolerance to T. annulata between the breeds. Potentially, our research may also have more general implications for the control of inflammatory processes against other pathogens. Genes controlling the balance between pathology and protection may determine how livestock can survive in the face of infectious onslaught. Next generation sequencing and RNAi methodologies for livestock species will bring new opportunities to link diversity at the genome level to functional differences in health traits in livestock species.

Analytical methods for disease association studies with immunogenetic data

Disease association studies involving highly polymorphic immunogenetic data may involve analyses at one or many units of analysis, including amino acid, allele, genotype and haplotype levels, as well as consideration of gene-gene or gene-environment interactions. The selection of the appropriate statistical tests is critical and will be dependent on the nature of the dataset (e.g., case-control vs. family data) as well as the specific research hypotheses being tested. This paper describes the various study and analysis categories used for such analyses, including the advantages and limitations of such techniques.

A community standard for immunogenomic data reporting and analysis: proposal for a STrengthening the REporting of Immunogenomic Studies statement

Modern high-throughput HLA and KIR typing technologies are generating a wealth of immunogenomic data with the potential to revolutionize the fields of histocompatibility and immune-related disease association and population genetic research, much as SNP-based approaches have revolutionized association research. The STrengthening the REporting of Genetic Association studies (STREGA) statement provides community-based data reporting and analysis standards for genomic disease-association studies, identifying specific areas in which adoption of reporting guidelines can improve the consistent interpretation of genetic studies. While aspects of STREGA can be applied to immunogenomic studies, HLA and KIR research requires additional consideration, as the high levels of polymorphism associated with immunogenomic data pose unique methodological and computational challenges to the synthesis of information across datasets. Here, we outline the principle challenges to consistency in immunogenomic studies, and propose that an immunogenomic-specific analog to the STREGA statement, a STrengthening the REporting of Immunogenomic Studies (STREIS) statement, be developed as part of the 16th International HLA and Immunogenetics Workshop. We propose that STREIS extends at least four of the 22 elements of the STREGA statement to specifically address issues pertinent to immunogenomic data: HLA and KIR nomenclature, data-validation, ambiguity resolution, and the analysis of highly polymorphic genetic systems. As with the STREGA guidelines, the intent behind STREIS is not to dictate the design of immunogenomic studies, but to ensure consistent and transparent reporting of research, facilitating the synthesis of HLA and KIR data across studies.

Reconsidering the association between the major histocompatibility complex and bipolar disorder

Bipolar disorder (BD) is a cyclical and chronic affective disorder, globally recognized as an important public health problem and characterized by mood changes with recurring phases such as mania and depression. It is considered a complex disease, depending on the interaction of genetic and environmental triggers (stressors factors), but with a poorly known pathogenesis. Recent studies have implicated immune factors in the pathogenesis of BD and more particularly associated with different human major histocompatibility complex (MHC) regions. A major consortium study have recently linked BD to hundreds of variations with stronger associations in the MHC region, such as the rs3130297 SNP, located in the NOTCH4 gene, with an additional overlapping association with schizophrenia. This short review focuses on studies that investigated the association between bipolar disorder and the MHC, and the involvement of the immune system in the pathogenesis of the disease, in order to provide further information for additional diagnostic and therapeutic strategies. Fully understanding the etiology and pathophysiology of BD is extremely important to define new approaches for intervention and prevention, maybe through the modulation of the immune system.

Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex

The human major histocompatibility complex (MHC) on chromosome 6p21 is a paradigm for genomics, showing remarkable polymorphism and striking association with immune and non-immune diseases. The complex genomic landscape of the MHC, notably strong linkage disequilibrium, has made resolving causal variants very challenging. A promising approach is to investigate gene expression levels considered as tractable intermediate phenotypes in mapping complex diseases. However, how transcription varies across the MHC, notably relative to specific haplotypes, remains unknown. Here, using an original hybrid tiling and splice junction microarray that includes alternate allele probes, we draw the first high-resolution strand-specific transcription map for three common MHC haplotypes (HLA-A1-B8-Cw7-DR3, HLA-A3-B7-Cw7-DR15, and HLA-A26-B18-Cw5-DR3-DQ2) strongly associated with autoimmune diseases including type 1 diabetes, systemic lupus erythematosus, and multiple sclerosis. We find that haplotype-specific differences in gene expression are common across the MHC, affecting 96 genes (46.4%), most significantly the zing finger protein gene ZFP57. Differentially expressed probes are correlated with polymorphisms between haplotypes, consistent with cis effects that we directly demonstrate for ZFP57 in a cohort of healthy volunteers (P = 1.2 × 10(-14)). We establish that alternative splicing is significantly more frequent in the MHC than genome-wide (72.5% vs. 62.1% of genes, P ≤ 1 × 10(-4)) and shows marked haplotypic differences. We also unmask novel and abundant intergenic transcription involving 31% of transcribed blocks identified. Our study reveals that the renowned MHC polymorphism also manifests as transcript diversity, and our novel haplotype-based approach marks a new step toward identification of regulatory variants involved in the control of MHC-associated phenotypes and diseases.

Accurate HLA type inference using a weighted similarity graph

Background

The human leukocyte antigen system (HLA) contains many highly variable genes. HLA genes play an important role in the human immune system, and HLA gene matching is crucial for the success of human organ transplantations. Numerous studies have demonstrated that variation in HLA genes is associated with many autoimmune, inflammatory and infectious diseases. However, typing HLA genes by serology or PCR is time consuming and expensive, which limits large-scale studies involving HLA genes. Since it is much easier and cheaper to obtain single nucleotide polymorphism (SNP) genotype data, accurate computational algorithms to infer HLA gene types from SNP genotype data are in need. To infer HLA types from SNP genotypes, the first step is to infer SNP haplotypes from genotypes. However, for the same SNP genotype data set, the haplotype configurations inferred by different methods are usually inconsistent, and it is often difficult to decide which one is true.

Results

In this paper, we design an accurate HLA gene type inference algorithm by utilizing SNP genotype data from pedigrees, known HLA gene types of some individuals and the relationship between inferred SNP haplotypes and HLA gene types. Given a set of haplotypes inferred from the genotypes of a population consisting of many pedigrees, the algorithm first constructs a weighted similarity graph based on a new haplotype similarity measure and derives constraint edges from known HLA gene types. Based on the principle that different HLA gene alleles should have different background haplotypes, the algorithm searches for an optimal labeling of all the haplotypes with unknown HLA gene types such that the total weight among the same HLA gene types is maximized. To deal with ambiguous haplotype solutions, we use a genetic algorithm to select haplotype configurations that tend to maximize the same optimization criterion. Our experiments on a previously typed subset of the HapMap data show that the algorithm is highly accurate, achieving an accuracy of 96% for gene HLA-A, 95% for HLA-B, 97% for HLA-C, 84% for HLA-DRB1, 98% for HLA-DQA1 and 97% for HLA-DQB1 in a leave-one-out test.

Conclusions

Our algorithm can infer HLA gene types from neighboring SNP genotype data accurately. Compared with a recent approach on the same input data, our algorithm achieved a higher accuracy. The code of our algorithm is available to the public for free upon request to the corresponding authors.

Regulation of major histocompatibility complex class II gene expression, genetic variation and disease

Major histocompatibility complex (MHC) class II molecules are central to adaptive immune responses and maintenance of self-tolerance. Since the early 1970s, the MHC class II region at chromosome 6p21 has been shown to be associated with a remarkable number of autoimmune, inflammatory and infectious diseases. Given that a full explanation for most MHC class II disease associations has not been reached through analysis of structural variation alone, in this review we examine the role of genetic variation in modulating gene expression. We describe the intricate architecture of the MHC class II regulatory system, indicating how its unique characteristics may relate to observed associations with disease. There is evidence that haplotype-specific variation involving proximal promoter sequences can alter the level of gene expression, potentially modifying the emergence and expression of key phenotypic traits. Although much emphasis has been placed on cis-regulatory elements, we also examine the role of more distant enhancer elements together with the evidence of dynamic inter- and intra-chromosomal interactions and epigenetic processes. The role of genetic variation in such mechanisms may hold profound implications for susceptibility to common disease.