Ancient haplotypes of the HLA Class II region

Qualitative, rather than quantitative, differences between HLA-DQ alleles affect HLA-DQ immunogenicity in organ transplantation

Prolonging the lifespan of transplanted organs is critical to combat the shortage of this life-saving resource. Chronic rejection, with irreversible demise of the allograft, is often caused by the development of donor-specific HLA antibodies. Currently, enumerating molecular (amino acid) mismatches between recipient and donor is promoted to identify patients at higher risk of developing HLA antibodies, for use in organ allocation, and immunosuppression-minimization strategies. We have counseled against the incorporation of such approaches into clinical use and hypothesized that not all molecular mismatches equally contribute to generation of donor-specific immune responses. Herein, we document statistical shortcomings in previous study design: for example, use of individuals who lack the ability to generate donor-specific-antibodies (HLA identical) as part of the negative cohort. We provide experimental evidence, using CRISPR-Cas9-edited cells, to rebut the claim that the HLAMatchmaker eplets represent "functional epitopes." We further used unique sub-cohorts of patients, those receiving an allograft with two HLA-DQ mismatches yet developing antibodies only to one mismatch (2MM1DSA), to interrogate differential immunogenicity. Our results demonstrate that mismatches of DQα05-heterodimers exhibit the highest immunogenicity. Additionally, we demonstrate that the DQα chain critically contributes to the overall qualities of DQ molecules. Lastly, our data proposes that an augmented risk to develop donor-specific HLA-DQ antibodies is dependent on qualitative (evolutionary and functional) divergence between recipient and donor, rather than the mere number of molecular mismatches. Overall, we propose an immunological mechanistic rationale to explain differential HLA-DQ immunogenicity, with potential ramifications for other pathological processes such as autoimmunity and infections.

Genome-wide association with footrot in hair and wool sheep

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

Narcolepsy type I-associated DNA methylation and gene expression changes in the human leukocyte antigen region

Narcolepsy type 1 (NT1) is caused by a loss of hypothalamic orexin-producing cells, and autoreactive CD4⁺ and CD8⁺ T cells have been suggested to play a role in the autoimmune mechanism. Although NT1 showed a strong association with human leukocyte antigen (HLA)-DQB1*06:02, the responsible antigens remain unidentified. We analyzed array-based DNA methylation and gene expression data for the HLA region in CD4⁺ and CD8⁺ T cells that were separated from the peripheral blood mononuclear cells of Japanese subjects (NT1, N = 42; control, N = 42). As the large number of SNPs in the HLA region might interfere with the affinity of the array probes, we conducted a comprehensive assessment of the reliability of each probe. The criteria were based on a previous study reporting that the presence of frequent SNPs, especially on the 3' side of the probe, makes the probe unreliable. We confirmed that 90.3% of the probes after general filtering in the HLA region do not include frequent SNPs, and are thus suitable for analysis, particularly in Japanese subjects. We then performed an association analysis, and found that several CpG sites in the HLA class II region of the patients were significantly hypomethylated in CD4⁺ and CD8⁺ T cells. This association was not detected when the effect of HLA-DQB1*06:02 was considered, suggesting that the hypomethylation was possibly derived from HLA-DQB1*06:02. Further RNA sequencing revealed reduced expression levels of HLA-DQB1 alleles other than HLA-DQB1*06:02 in the patients with NT1. Our results suggest the involvement of epigenetic and expressional changes in HLA-DQB1 in the pathogenesis of NT1.

HLA-DQ heterodimers in hematopoietic cell transplantation

HLA-DQ heterodimers increase the susceptibility to autoimmune diseases, but their role in hematopoietic cell transplantation is unknown. We tested the hypothesis that outcome after HLA-matched and HLA-DQ-mismatched hematopoietic cell transplantation is influenced by HLA-DQ heterodimers. Heterodimers were defined in 5164 HLA-matched and 520 HLA-DQ-mismatched patients and their transplant donors according to well-established crystallographic criteria. Group 1 (G1) heterodimers are any DQA1*02/03/04/05/06α paired with any DQB1*02/03/04β. Group 2 (G2) heterodimers are DQA1*01α paired with any DQB1*05/06β. Multivariable models identified significantly higher relapse risk in G1G2 and G2G2 compared with G1G1 HLA-matched patients with malignant disease; risk increased with an increasing number of G2 molecules. In HLA-DQ-mismatched transplantation for malignant diseases, matching or mismatching for G2 increased relapse risk. G2 lowered disease-free survival after both HLA-matched and HLA-DQ-mismatched transplantation. A paradigm based on HLA-DQ heterodimers provides a functional definition of the hematopoietic cell transplantation barrier and a means to lower risks for future patients.

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action

BACKGROUND

HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR).

METHODS

For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used.

RESULTS

All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found.

CONCLUSIONS

Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Evaluation of the major histocompatibility complex (MHC) class II as a candidate for sudden acquired retinal degeneration syndrome (SARDS) in Dachshunds

Objective

Sudden acquired retinal degeneration syndrome (SARDS) is one of the leading causes of acute blindness in dogs, with an unknown etiology and no effective treatment. Certain breeds such as Dachshunds are overrepresented among SARDS patients, and therefore, the syndrome is suspected to have a genetic component. The objective of this study was to determine if a genetic locus associated with SARDS in Dachshunds could be identified using a genome‐wide association study (GWAS).

Procedures

Genome‐wide association mapping was performed in 15 SARDS‐affected and 16 unaffected Dachshunds. Genotyping of three classical DLA class II genes (DLA‐DRB1, DLA‐DQA1, and DLA‐DQB1) was performed in 34 SARDS‐affected and 66 unaffected Dachshunds to evaluate for an association in this region.

Results

Although no single nucleotide polymorphisms (SNPs) were of genome‐wide statistical significance (P_Bonferroni < 0.05), 5 of the top 9 SNPs were in the major histocompatibility complex (MHC). Using DLA typing, the allele DLA‐DRB1*09401 was identified as a risk factor for the development of SARDS (P = 0.0032, OR = 4.0). The alleles DLA‐DQB1*00101 (P = 0.0050, OR = 0.31), DLA‐DQA1*00901 (P = 0.0087, OR = 0.33), and a previously identified DLA‐DRB1allele described as “DRB1‐T” (P = 0.0284, OR = 0.37) were identified as protective factors.

Conclusions

Although far from definitive, association of SARDS with alleles of immunologic importance further supports the hypothesis that autoimmunity may play a role in the pathogenesis of SARDS.

Peptide presentation by HLA-DQ molecules is associated with the development of immune tolerance

HLA class II proteins are important elements of human adaptive immune recognition and are associated with numerous infectious and immune-mediated diseases. These highly variable molecules can be classified into DP, DQ and DR groups. It has been proposed that in contrast with DP and DR, epitope binding by DQ variants rather results in immune tolerance. However, the pieces of evidence are limited and controversial. We found that DQ molecules bind more human epitopes than DR. Pathogen-associated epitopes bound by DQ molecules are more similar to human proteins than the ones bound by DR. Accordingly, DQ molecules bind epitopes of significantly different pathogen species. Moreover, the binding of autoimmunity-associated epitopes by DQ confers protection from autoimmune diseases. Our results suggest a special role of HLA-DQ in immune homeostasis and help to better understand the association of HLA molecules with infectious and autoimmune diseases.

Characterizing alleles with large deletions using region specific extraction

Two novel HLA class II alleles, DRB4*03:01N and DQB1*03:276N, containing large deletions were identified during routine typing. Extraction of DNA encompassing the deletions was carried out with a panel of capture oligonucleotides followed by whole genome amplification. Next generation DNA sequencing was then used to characterize the sequences. DRB4*03:01N has a 16 kilobase pair deletion stretching upstream from intron 2 toward centromeric DRB8. DQB1*03:276N has two deletions separated by 844 nucleotides. The first deletion (3.7 kilobase pairs) is upstream of intron 1 and the second deletion removes 3.3 kilobase pairs further upstream towards centromeric DQA2.

HLA frequency in candidates to transplant without compatible cord blood at the National Center of Blood Transfusion (Mexico)

INTRODUCTION

Umbilical Cord Blood Units (UCBU) for transplantation, are a therapeutic possibility for patients with a wide range of oncohaematological diseases and other immunologic disorders. The search of compatible donors for bone marrow transplantation is increasingly difficult for patients of mixed ethnicity. The aim of this work was determine the HLA frequency of candidates for transplantation without compatible UCBU at the National Center of from Blood Transfusion (NCBT) - Mexico.

MATERIAL AND METHODS

A retrospective analysis of candidates to transplant without compatible UCBU was performed in the archives from 2003 to 2016 at the NCBT. HLA class I and II genotyping of candidates was performed by medium resolution methods: Sequence Specific Primer and/or Specific Sequence Oligonucleotide. HLA frequencies were obtained by including individuals without any particular bias to a phenotype and strict statistical and genetic analysis of populations were done. The database in www.allelefrequencies.net was used in order to identify the ethnic origin of the most frequent alleles.

RESULTS

Three hundred and sixty-four candidates without compatible UCBU for transplantation were identified. The most frequent haplotype HLA I and II were: HLA-A*02/02, 02/24, 24/24, 02/68, 01/24 and 24/68; HLA-B*39/39, 35/51, 44/44, 44/40, 35/40 and 35/35; HLA-DRB1*04/04, 13/07, 04/13, 13/13 and 03/11. The ethnic origins of the analyzed data were represented in most cases by Amerindians, Caucasics, Orientals, Asians, Arabs and Africans.

CONCLUSION

This work shows the existence of a broad genetic diversity of candidates for transplantation with UCBU, making it difficult to find compatible units considering donors only from the capital.

Transancestral mapping and genetic load in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P<5 × 10⁻⁸), refined association signals in established regions, extended associations to additional ancestries, and a disentangled complex HLA multigenic effect. The risk allele count (genetic load) exhibits an accelerating pattern of SLE risk, leading us to posit a cumulative hit hypothesis for autoimmune disease. Comparing results across the three ancestries identifies both ancestry-dependent and ancestry-independent contributions to SLE risk. Our results are consistent with the unique and complex histories of the populations sampled, and collectively help clarify the genetic architecture and ethnic disparities in SLE.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong ethnic and gender bias. In a transancestral genetic association study, Langefeld et al. identify 24 novel regions associated with risk to lupus and propose a cumulative hits hypothesis for loci conferring risk to SLE.

An analytical workflow for accurate variant discovery in highly divergent regions

Background

Current variant discovery methods often start with the mapping of short reads to a reference genome; yet, their performance deteriorates in genomic regions where the reads are highly divergent from the reference sequence. This is particularly problematic for the human leukocyte antigen (HLA) region on chromosome 6p21.3. This region is associated with over 100 diseases, but variant calling is hindered by the extreme divergence across different haplotypes.

Results

We simulated reads from chromosome 6 exonic regions over a wide range of sequence divergence and coverage depth. We systematically assessed combinations between five mappers and five callers for their performance on simulated data and exome-seq data from NA12878, a well-studied individual in which multiple public call sets have been generated. Among those combinations, the number of known SNPs differed by about 5 % in the non-HLA regions of chromosome 6 but over 20 % in the HLA region. Notably, GSNAP mapping combined with GATK UnifiedGenotyper calling identified about 20 % more known SNPs than most existing methods without a noticeable loss of specificity, with 100 % sensitivity in three highly polymorphic HLA genes examined. Much larger differences were observed among these combinations in INDEL calling from both non-HLA and HLA regions. We obtained similar results with our internal exome-seq data from a cohort of chronic lymphocytic leukemia patients.

Conclusions

We have established a workflow enabling variant detection, with high sensitivity and specificity, over the full spectrum of divergence seen in the human genome. Comparing to public call sets from NA12878 has highlighted the overall superiority of GATK UnifiedGenotyper, followed by GATK HaplotypeCaller and SAMtools, in SNP calling, and of GATK HaplotypeCaller and Platypus in INDEL calling, particularly in regions of high sequence divergence such as the HLA region. GSNAP and Novoalign are the ideal mappers in combination with the above callers. We expect that the proposed workflow should be applicable to variant discovery in other highly divergent regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3045-z) contains supplementary material, which is available to authorized users.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

Genetic variation and the de novo assembly of human genomes

The discovery of genetic variation and the assembly of genome sequences are both inextricably linked to advances in DNA-sequencing technology. Short-read massively parallel sequencing has revolutionized our ability to discover genetic variation but is insufficient to generate high-quality genome assemblies or resolve most structural variation. Full resolution of variation is only guaranteed by complete de novo assembly of a genome. Here, we review approaches to genome assembly, the nature of gaps or missing sequences, and biases in the assembly process. We describe the challenges of generating a complete de novo genome assembly using current technologies and the impact that being able to perfectly sequence the genome would have on understanding human disease and evolution. Finally, we summarize recent technological advances that improve both contiguity and accuracy and emphasize the importance of complete de novo assembly as opposed to read mapping as the primary means to understanding the full range of human genetic variation.

Cell-surface MHC density profiling reveals instability of autoimmunity-associated HLA

Polymorphisms within HLA gene loci are strongly associated with susceptibility to autoimmune disorders; however, it is not clear how genetic variations in these loci confer a disease risk. Here, we devised a cell-surface MHC expression assay to detect allelic differences in the intrinsic stability of HLA-DQ proteins. We found extreme variation in cell-surface MHC density among HLA-DQ alleles, indicating a dynamic allelic hierarchy in the intrinsic stability of HLA-DQ proteins. Using the case-control data for type 1 diabetes (T1D) for the Swedish and Japanese populations, we determined that T1D risk-associated HLA-DQ haplotypes, which also increase risk for autoimmune endocrinopathies and other autoimmune disorders, encode unstable proteins, whereas the T1D-protective haplotypes encode the most stable HLA-DQ proteins. Among the amino acid variants of HLA-DQ, alterations in 47α, the residue that is located on the outside of the peptide-binding groove and acts as a key stability regulator, showed strong association with T1D. Evolutionary analysis suggested that 47α variants have been the target of positive diversifying selection. Our study demonstrates a steep allelic hierarchy in the intrinsic stability of HLA-DQ that is associated with T1D risk and protection, suggesting that HLA instability mediates the development of autoimmune disorders.

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.

Diplotyper: diplotype-based association analysis

Background

It was previously reported that an association analysis based on haplotype clusters increased power over single-locus tests, and that another association test based on diplotype trend regression analysis outperformed other, more common association approaches. We suggest a novel algorithm to combine haplotype cluster- and diplotype-based analyses.

Methods

Diplotyper combines a novel algorithm designed to cluster haplotypes of interest from a given set of haplotypes with two existing tools: Haploview, for analyses of linkage disequilibrium blocks and haplotypes, and PLINK, to generate all possible diplotypes from given genotypes of samples and calculate linear or logistic regression. In addition, procedures for generating all possible diplotypes from the haplotype clusters and transforming these diplotypes into PLINK formats were implemented.

Results

Diplotyper is a fully automated tool for performing association analysis based on diplotypes in a population. Diplotyper was tested through association analysis of hepatic lipase (LIPC) gene polymorphisms or diplotypes and levels of high-density lipoprotein (HDL) cholesterol.

Conclusions

Diplotyper is useful for identifying more precise and distinct signals over single-locus tests.

Microsatellite variation in the equine MHC

Genes within the major histocompatibility complex (MHC) encode proteins involved in innate and adaptive immune responses. Genetic variation in this region can influence the immune response of an individual animal to challenges from a variety of pathogens; however, a complete documentation of genetic variation in the MHC is lacking for most domestic animals, including horses. To provide additional genetic markers for study of the horse MHC, or ELA (equine lymphocyte antigen), we identified 37 polymorphic microsatellite repeats in ELA and used these variations separately and together with published SNPs to investigate linkage disequilibrium (LD) and haplotype structure in a sample of Thoroughbred horses. ELA SNPs alone detected little LD, but microsatellites, either separately or combined with SNPs, revealed substantially more LD. A subset of markers in very high LD across the breadth of ELA may be predictive of structural polymorphisms or linked epistases that are important drivers of haplotype structure in Thoroughbreds.

Human genetic individuality

The central preoccupation of human genetics is an effort to understand the genotypic basis of human phenotypic diversity. Although recent progress in identifying the genes that, when mutated, underlie major genetic diseases has been rapid, knowledge of the genetic influences on the vast range of variable, and at least partially heritable, traits that constitute the "normal" range of human phenotypic variation lags. Spectacular advances in our knowledge of human genetic variation have laid the groundwork for a synthesis of insights from medical genetics, population genetics, molecular evolution, and the study of human origins that places basic constraints on models of human genetic individuality. Balancing selection, local adaptation, mutation-selection balance, and founder effects have all extensively shaped contemporary genetic variation. Long-term-balancing selection appears largely to reflect the consequences of host-pathogen arms races. Local adaptation has been widespread-and involved responses to a plethora of selective pressures, some identifiable but most unknown. However, it appears to be a combination of mutation-selection balance and founder effects that largely accounts for genetic individuality. If true, this inference has major implications for future research programs in human genetics.

Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird

The extent to which pathogens maintain the extraordinary polymorphism at vertebrate Major Histocompatibility Complex (MHC) genes via balancing selection has intrigued evolutionary biologists for over half a century, but direct tests remain challenging. Here we examine whether a well-characterized epidemic of Mycoplasmal conjunctivitis resulted in balancing selection on class II MHC in a wild songbird host, the house finch (Carpodacus mexicanus). First, we confirmed the potential for pathogen-mediated balancing selection by experimentally demonstrating that house finches with intermediate to high multi-locus MHC diversity are more resistant to challenge with Mycoplasma gallisepticum. Second, we documented sequence and diversity-based signatures of pathogen-mediated balancing selection at class II MHC in exposed host populations that were absent in unexposed, control populations across an equivalent time period. Multi-locus MHC diversity significantly increased in exposed host populations following the epidemic despite initial compromised diversity levels from a recent introduction bottleneck in the exposed host range. We did not observe equivalent changes in allelic diversity or heterozygosity across eight neutral microsatellite loci, suggesting that the observations reflect selection rather than neutral demographic processes. Our results indicate that a virulent pathogen can exert sufficient balancing selection on class II MHC to rescue compromised levels of genetic variation for host resistance in a recently bottlenecked population. These results provide evidence for Haldane's long-standing hypothesis that pathogens directly contribute to the maintenance of the tremendous levels of genetic variation detected in natural populations of vertebrates.

Expanded dog leukocyte antigen (DLA) single nucleotide polymorphism (SNP) genotyping reveals spurious class II associations

The dog leukocyte antigen (DLA) system contains many of the functional genes of the immune system, thereby making it a candidate region for involvement in immune-mediated disorders. A number of studies have identified associations between specific DLA class II haplotypes and canine immune hemolytic anemia, thyroiditis, immune polyarthritis, type I diabetes mellitus, hypoadrenocorticism, systemic lupus erythematosus-related disease complex, necrotizing meningoencephalitis (NME) and anal furunculosis. These studies have relied on sequencing approximately 300 bases of exon 2 of each of the DLA class II genes: DLA-DRB1, DLA-DQA1 and DLA-DQB1. In the present study, an association (odds ratio=4.29) was identified by this method between Weimaraner dogs with hypertrophic osteodystrophy (HOD) and DLA-DRB1∗01501. To fine map the association with HOD, a genotyping assay of 126 coding single nucleotide polymorphisms (SNPs) from across the entire DLA, spanning a region of 2.5 Mb (3,320,000-5,830,000) on CFA12, was developed and tested on Weimaraners with HOD, as well as two additional breeds with diseases associated with DLA class II: Nova Scotia duck tolling retrievers with hypoadrenocorticism and Pug dogs with NME. No significant associations were found between Weimaraners with HOD or Nova Scotia duck tolling retrievers with hypoadrenocorticism and SNPs spanning the DLA region. In contrast, significant associations were found with NME in Pug dogs, although the associated region extended beyond the class II genes. By including a larger number of genes from a larger genomic region, a SNP genotyping assay was generated that provides coverage of the extended DLA region and may be useful in identifying and fine mapping DLA associations in dogs.

Haplotype-based genomic sequencing of a chromosomal polymorphism in the white-throated sparrow (Zonotrichia albicollis)

Inversion polymorphisms have been linked to a variety of fundamental biological and evolutionary processes. Yet few studies have used large-scale genomic sequencing to directly compare the haplotypes associated with the standard and inverted chromosome arrangements. Here we describe the targeted genomic sequencing and comparison of haplotypes representing alternative arrangements of a common inversion polymorphism linked to a suite of phenotypes in the white-throated sparrow (Zonotrichia albicollis). More than 7.4 Mb of genomic sequence was generated and assembled from both the standard (ZAL2) and inverted (ZAL2(m)) arrangements. Sequencing of a pair of inversion breakpoints led to the identification of a ZAL2-specific segmental duplication, as well as evidence of breakpoint reusage. Comparison of the haplotype-based sequence assemblies revealed low genetic differentiation outside versus inside the inversion indicative of historical patterns of gene flow and suppressed recombination between ZAL2 and ZAL2(m). Finally, despite ZAL2(m) being maintained in a near constant state of heterozygosity, no signatures of genetic degeneration were detected on this chromosome. Overall, these results provide important insights into the genomic attributes of an inversion polymorphism linked to mate choice and variation in social behavior.

Trans-species polymorphism, HLA-disease associations and the evolution of the MHC

Currently, the paradigm is that major histocompatibility complex (MHC) polymorphism is maintained by balancing selection on the immune genes. However, other evolutionary forces besides selection also play a role in the population genetics of this multigene family. van Oosterhout proposed a new theory of MHC evolution called associative balancing complex (ABC) evolution.1 This theory incorporates the effects of the evolutionary forces in the entire MHC region (peri-MHC), and it proposes that recessive deleterious mutations can accumulate in the peri-MHC in a process similar to Muller's ratchet.2 These mutations are not easily purged because epistasis and high gene diversity in the MHC reduce the efficacy of natural selection. Because natural selection is less efficient, it could also make the MHC prone to the onslaught of genomic parasites such as retroviruses and transposable elements (TEs). The accumulated genetic load has important consequences for the evolution of this immune gene family, and it can reinforce linkage disequilibria and help to maintain the MHC polymorphism. ABC evolution offers new insights into some of the most puzzling aspects of the MHC, including the occurrence of identical MHC sequences in diverged species (i.e., trans-species polymorphism). It may also explain why the large numbers of disease-associated mutations are not removed by natural selection, and why the genes that protect vertebrates against infectious diseases are associated to such a wide variety of genetic disorders.

Darwin and genetics

Darwin's theory of natural selection lacked an adequate account of inheritance, making it logically incomplete. We review the interaction between evolution and genetics, showing how, unlike Mendel, Darwin's lack of a model of the mechanism of inheritance left him unable to interpret his own data that showed Mendelian ratios, even though he shared with Mendel a more mathematical and probabilistic outlook than most biologists of his time. Darwin's own "pangenesis" model provided a mechanism for generating ample variability on which selection could act. It involved, however, the inheritance of characters acquired during an organism's life, which Darwin himself knew could not explain some evolutionary situations. Once the particulate basis of genetics was understood, it was seen to allow variation to be passed intact to new generations, and evolution could then be understood as a process of changes in the frequencies of stable variants. Evolutionary genetics subsequently developed as a central part of biology. Darwinian principles now play a greater role in biology than ever before, which we illustrate with some examples of studies of natural selection that use DNA sequence data and with some recent advances in answering questions first asked by Darwin.

Polymorphic major histocompatibility complex class II Alu insertions at five loci and their association with HLA-DRB1 and -DQB1 in Japanese and Caucasians

We investigated polymorphic Alu insertion (POALIN) frequencies at five loci in the major histocompatibility complex (MHC) class II genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-DRB1 and -DQB1 genes for 100 Japanese, 174 Australian Caucasians and 67 HLA reference cell lines obtained from different ethnic groups. The POALINs varied in frequency between 11% and 57% with significant differences between the Japanese and Caucasians at three loci. One POALIN locus deviated significantly from Hardy-Weinberg equilibrium (HWE) and four POALIN loci were in significant linkage disequilibrium and had a high percentage association with a variety of HLA-DRB1 or -DQB1 two-digit alleles. Inferred haplotype analysis among two-locus, five-locus and seven-locus haplotype structures showed maximum differences between the Japanese and Caucasians with the seven-locus haplotypes. The most common multilocus haplotype in Caucasians was DRB1*1501/DQB1*0602/AluDQ1/AluDRB1/AluORF10/AluDPB2 (6.7%), whereas the second most common allele HLA-DRB1*15 (17.5%) in Japanese was associated with three or four Alu insertions. The HLA class II POALINs also differentiated within and between HLA-DRB1 super-haplotypes DR1, DR8, DR51, DR52 and DR53. This is the first comparative population study of multilocus POALINs in the HLA class II region, which shows that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic markers for the identification of allele and haplotype lineages and variations within the same and/or different populations.

A new theory of MHC evolution: beyond selection on the immune genes

The major histocompatibility complex (MHC) is a dense region of immune genes with high levels of polymorphism, which are arranged in haplotype blocks. Traditional models of balancing selection (i.e. overdominance and negative frequency dependence) were developed to study the population genetics of single genes. However, the MHC is a multigene family surrounded by linked (non-neutral) polymorphisms, and not all of its features are well explained by these models. For example, (i) the high levels of polymorphism in small populations, (ii) the unexpectedly large genetic differentiation between populations, (iii) the shape of the allelic genealogy associated with trans-species evolution, and (iv) the close associations between particular MHC (human leucocyte antigen, HLA) haplotypes and the approximately 100 pathologies in humans. Here, I propose a new model of MHC evolution named Associative Balancing Complex evolution that can explain these phenomena. The model proposes that recessive deleterious mutations accumulate as a 'sheltered load' nearby MHC genes. These mutations can accumulate because (i) they are rarely expressed as homozygotes given the high MHC gene diversity and (ii) purifying selection is inefficient with low recombination rates (cf. Muller's ratchet). Once fixed, these mutations add to balancing selection and further reinforce linkage through epistatic selection against recombinants.

A comprehensive evaluation of SNP genotype imputation

Genome-wide association studies have contributed significantly to the genetic dissection of complex diseases. In order to increase the power of existing marker sets even further, methods have been proposed to predict individual genotypes at un-typed loci from other marker sets by imputation, usually employing HapMap data as a reference. Although various imputation algorithms have been used in practice already, a comprehensive evaluation and comparison of these approaches, using genome-wide SNP data from one and the same population is still lacking. We therefore investigated four publicly available programs for genotype imputation (BEAGLE, IMPUTE, MACH, and PLINK) using data from 449 German individuals genotyped in our laboratory for three genome-wide SNP sets [Affymetrix 5.0 (500 k), Affymetrix 6.0 (1,000 k), and Illumina 550 k]. We observed that HapMap-based imputation in a northern European population is powerful and reliable, even in highly variable genomic regions such as the extended MHC on chromosome 6p21. However, while genotype predictions were found to be highly accurate with all four programs, the number of SNPs for which imputation was actually carried out ('imputation efficacy') varied substantially. BEAGLE, IMPUTE, and MACH yielded nearly identical trade-offs between imputation accuracy and efficacy whereas PLINK performed consistently poorer. We nevertheless recommend either MACH or BEAGLE for practical use because these two programs are more user-friendly and generally require less memory than IMPUTE.

Genomic imprinting in the development and evolution of psychotic spectrum conditions

I review and evaluate genetic and genomic evidence salient to the hypothesis that the development and evolution of psychotic spectrum conditions have been mediated in part by alterations of imprinted genes expressed in the brain. Evidence from the genetics and genomics of schizophrenia, bipolar disorder, major depression, Prader-Willi syndrome, Klinefelter syndrome, and other neurogenetic conditions support the hypothesis that the etiologies of psychotic spectrum conditions commonly involve genetic and epigenetic imbalances in the effects of imprinted genes, with a bias towards increased relative effects from imprinted genes with maternal expression or other genes favouring maternal interests. By contrast, autistic spectrum conditions, including Kanner autism, Asperger syndrome, Rett syndrome, Turner syndrome, Angelman syndrome, and Beckwith-Wiedemann syndrome, commonly engender increased relative effects from paternally expressed imprinted genes, or reduced effects from genes favouring maternal interests. Imprinted-gene effects on the etiologies of autistic and psychotic spectrum conditions parallel the diametric effects of imprinted genes in placental and foetal development, in that psychotic spectrum conditions tend to be associated with undergrowth and relatively-slow brain development, whereas some autistic spectrum conditions involve brain and body overgrowth, especially in foetal development and early childhood. An important role for imprinted genes in the etiologies of psychotic and autistic spectrum conditions is consistent with neurodevelopmental models of these disorders, and with predictions from the conflict theory of genomic imprinting.

Human MHC architecture and evolution: implications for disease association studies

Major histocompatibility complex (MHC) variation is a key determinant of susceptibility and resistance to a large number of infectious, autoimmune and other diseases. Identification of the MHC variants conferring susceptibility to disease is problematic, due to high levels of variation and linkage disequilibrium. Recent cataloguing and analysis of variation over the complete MHC has facilitated localization of susceptibility loci for autoimmune diseases, and provided insight into the MHC's evolution. This review considers how the unusual genetic characteristics of the MHC impact on strategies to identify variants causing, or contributing to, disease phenotypes. It also considers the MHC in relation to novel mechanisms influencing gene function and regulation, such as epistasis, epigenetics and microRNAs. These developments, along with recent technological advances, shed light on genetic association in complex disease.

A highly divergent microsatellite facilitating fast and accurate DRB haplotyping in humans and rhesus macaques

The DRB region of the MHC in primate species is known to display abundant region configuration polymorphism with regard to the number and content of genes present per haplotype. Furthermore, depending on the species studied, the different DRB genes themselves may display varying degrees of allelic polymorphism. Because of this combination of diversity (differential gene number) and polymorphism (allelic variation), molecular typing methods for the primate DRB region are cumbersome. All intact DRB genes present in humans and rhesus macaques appear to possess, however, a complex and highly divergent microsatellite. Microsatellite analysis of a sizeable panel of outbred rhesus macaques, covering most of the known Mamu-DRB haplotypes, resulted in the definition of unique genotyping patterns that appear to be specific for a given haplotype. Subsequent examination of a representative panel of human cells illustrated that this approach also facilitates high-resolution HLA-DRB typing in an easy, quick, and reproducible fashion. The genetic composition of this complex microsatellite is shown to be in concordance with the phylogenetic relationships of various HLA-DRB and Mamu-DRB exon 2 gene/lineage sequences. Moreover, its length variability segregates with allelic variation of the respective gene. This simple protocol may find application in a variety of research avenues such as transplantation biology, disease association studies, molecular ecology, paternity testing, and forensic medicine.

Full-length sequence analysis of the HLA-DRB1 locus suggests a recent origin of alleles

The HLA region harbors some of the most polymorphic loci in the human genome. Among them is the class II locus HLA-DRB1, with more than 400 known alleles. The age of the polymorphism and the rate at which new alleles are generated at HLA loci has caused much controversy over the years. Previous studies have mostly been restricted to the 270 base pairs that constitute the second exon and represent the most variable part of the gene. Here, we investigate the evolutionary history of the HLA-DRB1 locus on the basis of an analysis of 15 genomic full-length alleles (10-15 kb). In addition, the variation in 49 complete coding sequences and 322 exon 2 sequences were analyzed. When excluding exon 2 from the analysis, the diversity at the synonymous sites was found to be similar to the intron diversity. The overall diversity in noncoding region was also similar to the genome average. The DRB1*03 lineage has been found in human, chimpanzee, bonobo, gorilla, and orangutan. An ancestral "proto HLA-DRB1*03 lineage" appeared to have diverged in the last 5 million years into the human-specific lineages *08, *11, *13, and *14. With exception to exon 2, both the coding- and the noncoding diversity suggests a recent origin (<1 million years ago) for most of the alleles at the HLA-DRB1 locus. Sites encoding for amino acids involved in antigen binding [antigen recognizing sites (ARS)] appear to have a more ancient origin. Taken together, the recent origin of most alleles, the high diversity between allelic lineages, and the ancient origin of sequence motifs in exon 2, is consistent with a relatively rapid generation of novel alleles by gene conversion like events.

Genetic control of alternative splicing in the TAP2 gene: possible implication in the genetics of type 1 diabetes

The transporter 2, ATP-binding cassette, subfamily B (TAP2) is involved in the transport of antigenic peptides to HLA molecules. Coding TAP2 polymorphisms shows a strong association with type 1 diabetes, but it is not clear whether this association may be entirely due to linkage disequilibrium with HLA DR and DQ. Functionally, rat Tap2 nonsynonymous single-nucleotide polymorphisms (nsSNPs) confer differential selectivity for antigenic peptides, but this was not shown to be the case for human TAP2 nsSNPs. In the human, differential peptide selectivity is rather conferred by two splicing isoforms with alternative carboxy terminals. Here, we tested the hypothesis that alleles at the coding SNPs favor different splicing isoforms, thus determining peptide selectivity indirectly. This may be the basis for independent contribution to the type 1 diabetes association. In RNA from heterozygous lymphoblastoid lines, we measured the relative abundance of each SNP haplotype in each isoform. In isoform NM_000544, the G (Ala) allele at 665 Thr>Ala (rs241447) is more than twice as abundant as A (Thr) (GA = 2.2 +/- 0.4, P = 1.5 x 10(-4)), while isoform NM_018833 is derived almost exclusively from chromosomes carrying A (AG = 18.1 +/- 5.6, P = 2.04 x 10(-7)). In 889 Canadian children with type 1 diabetes, differential transmission of parental TAP2 alleles persisted (P = 0.011) when analysis was confined to chromosomes carrying only DQ*02 alleles, which mark a conserved DR-DQ haplotype, thus eliminating most of the variation at DR-DQ. Thus, we present evidence of TAP2 association with type 1 diabetes that is independent of HLA DR-DQ and describe a plausible functional mechanism based on allele dependence of splicing into isoforms known to have differential peptide selectivities.

Ancient polymorphism and functional variation in the primate MHC-DQA1 5' cis-regulatory region

Precise regulation of MHC gene expression is critical to vertebrate immune surveillance and response. Polymorphisms in the 5' proximal promoter region of the human class II gene HLA-DQA1 have been shown to influence its transcriptional regulation and may contribute to the pathogenesis of autoimmune diseases. We investigated the evolutionary history of this cis-regulatory region by sequencing the DQA1 5' proximal promoter region in eight nonhuman primate species. We observed unexpectedly high levels of sequence variation and multiple strong signatures of balancing selection in this region. Specifically, the considerable DQA1 promoter region diversity was characterized by abundant shared (or trans-species) polymorphism and a pronounced lack of fixed differences between species. The majority of transcription factor binding sites in the DQA1 promoter region were polymorphic within species, and these binding site polymorphisms were commonly shared among multiple species despite evidence for negative selection eliminating a significant fraction of binding site mutations. We assessed the functional consequences of intraspecific promoter region diversity using a cell line-based reporter assay and detected significant differences among baboon DQA1 promoter haplotypes in their ability to drive transcription in vitro. The functional differentiation of baboon promoter haplotypes, together with the significant deviations from neutral sequence evolution, suggests a role for balancing selection in the evolution of DQA1 transcriptional regulation in primates.

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.