HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Immunity and mental illness: findings from a Danish population-based immunogenetic study of seven psychiatric and neurodevelopmental disorders

Human leukocyte antigen (HLA) genes encode proteins with important roles in the regulation of the immune system. Many studies have also implicated HLA genes in psychiatric and neurodevelopmental disorders. However, these studies usually focus on one disorder and/or on one HLA candidate gene, often with small samples. Here, we access a large dataset of 65,534 genotyped individuals consisting of controls (N = 19,645) and cases having one or more of autism spectrum disorder (N = 12,331), attention deficit hyperactivity disorder (N = 14,397), schizophrenia (N = 2401), bipolar disorder (N = 1391), depression (N = 18,511), anorexia (N = 2551) or intellectual disability (N = 3175). We imputed participants' HLA alleles to investigate the involvement of HLA genes in these disorders using regression models. We found a pronounced protective effect of DPB1*1501 on susceptibility to autism (p = 0.0094, OR = 0.72) and intellectual disability (p = 0.00099, OR = 0.41), with an increased protective effect on a comorbid diagnosis of both disorders (p = 0.003, OR = 0.29). We also identified a risk allele for intellectual disability, B*5701 (p = 0.00016, OR = 1.33). Associations with both alleles survived FDR correction and a permutation procedure. We did not find significant evidence for replication of previously-reported associations for autism or schizophrenia. Our results support an implication of HLA genes in autism and intellectual disability, which requires replication by other studies. Our study also highlights the importance of large sample sizes in HLA association studies.

Multiplicative fitness, rapid haplotype discovery, and fitness decay explain evolution of human MHC

The major histocompatibility complex (MHC) is a central component of the vertebrate immune system and hence evolves in the regime of a host-pathogen evolutionary race. The MHC is associated with quantitative traits which directly affect fitness and are subject to selection pressure. The evolution of haplotypes at the MHC HLA (HLA) locus is generally thought to be governed by selection for increased diversity that is manifested in overdominance and/or negative frequency-dependent selection (FDS). However, recently, a model combining purifying selection on haplotypes and balancing selection on alleles has been proposed. We compare the predictions of several population dynamics models of haplotype frequency evolution to the distributions derived from 6.59-million-donor HLA typings from the National Marrow Donor Program registry. We show that models that combine a multiplicative fitness function, extremely high haplotype discovery rates, and exponential fitness decay over time produce the best fit to the data for most of the analyzed populations. In contrast, overdominance is not supported, and population substructure does not explain the observed haplotype frequencies. Furthermore, there is no evidence of negative FDS. Thus, multiplicative fitness, rapid haplotype discovery, and rapid fitness decay appear to be the major factors shaping the HLA haplotype frequency distribution in the human population.

Towards uniformity in the definition of acceptable mismatches for highly sensitized patients

The Eurotransplant (The Eurotransplant International Foundation) acceptable mismatch programme has been shown to be a successful tool to enhance transplantation of highly sensitized patients(HSPs). However, patients with rare HLA phenotypes in relation to the Eurotransplant donor population remain on the waiting list. EUROSTAM is an European Union funded project to explore the feasibility of a Europe‐wide acceptable mismatch programme enabling transplantation of HSPs with rare HLA phenotypes within their own organ exchange organization. The present study, which forms part of the EUROSTAM project, assesses the differences in the practices of the laboratories in different countries with respect to their HLA antibody profiling and risk adverseness. In the serum exchange exercises of 18 samples, a high level of variability has been shown in both assays and interpretation of results. In the data exchange exercise when all participants were given the same Luminex raw data for analysis, a high degree of consensus was reached where the median fluorescent intensity values of beads were <500 and >2000 for standard single antigen bead assays, or <500 and >5000 for assignment of acceptable mismatches. The risk adverseness analysis has showed distinct patterns of attitudes towards the perceived risks based on HLA antibody assay results, most probably influenced by the local protocols of the clinical transplant programme of each laboratory. In order to ensure fairness and maintain consistencies of organ exchange among partner transplant centres, a centralized facility will be instrumental for a uniform definition of acceptable mismatches.

Deciphering the fine nucleotide diversity of full HLA class I and class II genes in a well-documented population from sub-Saharan Africa

With the aim to understand how next‐generation sequencing (NGS) improves both our assessment of genetic variation within populations and our knowledge on HLA molecular evolution, we sequenced and analysed 8 HLA loci in a well‐documented population from sub‐Saharan Africa (Mandenka). The results of full‐gene NGS‐MiSeq sequencing compared with those obtained by traditional typing techniques or limited sequencing strategies showed that segregating sites located outside exon 2 are crucial to describe not only class I but also class II population diversity. A comprehensive analysis of exons 2, 3, 4 and 5 nucleotide diversity at the 8 HLA loci revealed remarkable differences among these gene regions, notably a greater variation concentrated in the antigen recognition sites of class I exons 3 and some class II exons 2, likely associated with their peptide‐presentation function, a lower diversity of HLA‐C exon 3, possibly related to its role as a KIR ligand, and a peculiar molecular diversity of HLA‐A exon 2, revealing demographic signals. Based on full‐length HLA sequences, we also propose that the most frequent DRB1 allele in the studied population, DRB1*13:04, emerged from an allelic conversion involving 3 potential alleles as donors and DRB1*11:02:01 as recipient. Finally, our analysis revealed a high occurrence of the DRB1*13:04‐DQA1*05:05:01‐DQB1*03:19 haplotype, possibly resulting from a selective sweep due to protection to Onchorcerca volvulus, a prevalent pathogen in West Africa. This study unveils highly relevant information on the molecular evolution of HLA genes in relation to their immune function, calling for similar analyses in other populations living in contrasting environments.

Genetic diversity of the Thao people of Taiwan using Y-chromosome, mitochondrial DNA and HLA gene systems

BACKGROUND

Despite attempts in retracing the history of the Thao people in Taiwan using folktales, linguistics, physical anthropology, and ethnic studies, their history remains incomplete. The heritage of Thao has been associated with the Pazeh Western plains peoples and several other mountain peoples of Taiwan. In the last 400 years, their culture and genetic profile have been reshaped by East Asian migrants. They were displaced by the Japanese and the construction of a dam and almost faced extinction. In this paper, genetic information from mitochondrial DNA (mtDNA), Histoleucocyte antigens (HLA), and the non-recombining Y chromosome of 30 Thao individuals are compared to 836 other Taiwan Mountain and Plains Aborigines (TwrIP & TwPp), 384 Non-Aboriginal Taiwanese (non-TwA) and 149 Continental East Asians.

RESULTS

The phylogeographic analyses of mtDNA haplogroups F4b and B4b1a2 indicated gene flow between Thao, Bunun, and Tsou, and suggested a common ancestry from 10,000 to 3000 years ago. A claim of close contact with the heavily Sinicized Pazeh of the plains was not rejected and suggests that the plains and mountain peoples most likely shared the same Austronesian agriculturist gene pool in the Neolithic.

CONCLUSIONS

Having been moving repeatedly since their arrival in Taiwan between 6000 and 4500 years ago, the Thao finally settled in the central mountain range. They represent the last plains people whose strong bonds with their original culture allowed them to preserve their genetic heritage, despite significant gene flow from the mainland of Asia. Representing a considerable contribution to the genealogical history of the Thao people, the findings of this study bear on ongoing anthropological and linguistic debates on their origin.

Avian Toll-like receptor allelic diversity far exceeds human polymorphism: an insight from domestic chicken breeds

Immune genes show remarkable levels of adaptive variation shaped by pathogen-mediated selection. Compared to humans, however, population polymorphism in animals has been understudied. To provide an insight into immunogenetic diversity in birds, we sequenced complete protein-coding regions of all Toll-like receptor (TLR) genes with direct orthology between mammals and birds (TLR3, TLR4, TLR5 and TLR7) in 110 domestic chickens from 25 breeds and compared their variability with a corresponding human dataset. Chicken TLRs (chTLRs) exhibit on average nine-times higher nucleotide diversity than human TLRs (hTLRs). Increased potentially functional non-synonymous variability is found in chTLR ligand-binding ectodomains. While we identified seven sites in chTLRs under positive selection and found evidence for convergence between alleles, no selection or convergence was detected in hTLRs. Up to six-times more alleles were identified in fowl (70 chTLR4 alleles vs. 11 hTLR4 alleles). In chTLRs, high numbers of alleles are shared between the breeds and the allelic frequencies are more equal than in hTLRs. These differences may have an important impact on infectious disease resistance and host-parasite co-evolution. Though adaptation through high genetic variation is typical for acquired immunity (e.g. MHC), our results show striking levels of intraspecific polymorphism also in poultry innate immune receptors.

The genetic architecture of type 1 diabetes mellitus

Type 1 diabetes mellitus (T1D) is a complex autoimmune disorder characterised by loss of the insulin-producing pancreatic beta cells in genetically predisposed individuals, ultimately resulting in insulin deficiency and hyperglycaemia. T1D is most common among children and young adults, and the incidence is on the rise across the world. The aetiology of T1D is hypothesized to involve genetic and environmental factors that result in the T-cell mediated destruction of pancreatic beta cells. There is a strong genetic risk to T1D; with genome-wide association studies (GWAS) identifying over 60 susceptibility regions within the human genome which are marked by single nucleotide polymorphisms (SNPs). Here, we review what is currently known about the genetics of T1D. We argue that advancing our understanding of the aetiology and pathogenesis of T1D will require the integration of genome biology (omics-data) with GWAS data, thereby making it possible to elucidate the putative gene regulatory networks modulated by disease-associated SNPs. This approach has a potential to revolutionize clinical management of T1D in an era of precision medicine.

Fine mapping of MHC region in lung cancer highlights independent susceptibility loci by ethnicity

The basis for associations between lung cancer and major histocompatibility complex genes is not completely understood. Here the authors further consider genetic variation within the MHC region in lung cancer patients and identify independent associations within HLA genes that explain MHC lung cancer associations in Europeans and Asian populations.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA–tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

Human Leukocyte Antigen-A, B, C, DRB1, and DQB1 Allele and Haplotype Frequencies in a Subset of 237 Donors in the South African Bone Marrow Registry

Human leukocyte antigen- (HLA-) A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 allele and haplotype frequencies were studied in a subset of 237 volunteer bone marrow donors registered at the South African Bone Marrow Registry (SABMR). Hapl-o-Mat software was used to compute allele and haplotype frequencies from individuals typed at various resolutions, with some alleles in multiple allele code (MAC) format. Four hundred and thirty-eight HLA-A, 235 HLA-B, 234 HLA-DRB1, 41 HLA-DQB1, and 29 HLA-C alleles are reported. The most frequent alleles were A∗02:02g (0.096), B∗07:02g (0.082), C∗07:02g (0.180), DQB1∗06:02 (0.157), and DRB1∗15:01 (0.072). The most common haplotype was A∗03:01g~B∗07:02g~C∗07:02g~DQB1∗06:02~DRB1∗15:01 (0.067), which has also been reported in other populations. Deviations from Hardy-Weinberg equilibrium were observed in A, B, and DRB1 loci, with C~DQB1 being the only locus pair in linkage disequilibrium. This study describes allele and haplotype frequencies from a subset of donors registered at SABMR, the only active bone marrow donor registry in Africa. Although the sample size was small, our results form a key resource for future population studies, disease association studies, and donor recruitment strategies.

HLA class I (-A, -B, -C) and class II (-DR, -DQ) polymorphism in the Mauritanian population

BACKGROUND

HLA antigens have been widely studied for their role in transplantation biology, human diseases and population diversity. The aim of this study was to provide the first profile of HLA class I and class II alleles in the Mauritanian population.

METHODS

HLA typing was carried in 93 healthy Mauritanian blood donors, using single specific primer amplification (PCR-SSP).

RESULTS

Occurrences of the main HLA class I (-A, -B, -C) and class II (-DR, -DQ) antigens in the general population showed that out of the 17 HLA-A allele groups detected, five main HLA-A allele groups: A*02 (18.42%), A*01 (14.04%), A*23 (14.04%), A*30 (13.16%) and A*29 (12.28%) were the most common identified along other 12 relatively minor allele groups. Twenty three allele groups were observed in the locus B of which B*07 (13.46%) was the most prevalent followed by B*15, B*35, B*08 and B*27 all, with a frequency between 7 to 8%. Three prevalent HLA-C allele groups (C*02: 35.09%, C*07: 20.19% and C*06: 13.6%) were detected. The main HLA class II observed allele groups were: DRB1*13 (27.42%), DRB1*03 (24.73%), DRB1*11 (13.98%), DQB1*03 (36.03%), DQB1*02 (22.06%) and DQB1*05 (18.8%). Except for few haplotype in class I (A*02-B*07: 4.45%, A*02-C02: 10%, A*23-C*02: 8.8%, B*07-C*02: 8.8%, B*15-C*02: 8.8%) and in class II (DRB1*13-DQB1*06: 11.94%, DRB1*03-DQB1*02:11.19% and DRB1*03-DQB1*03: 10.45%), the majority of locus combination were in the range of 2-3%. A single predominant haplotype C*02-DRB1*03 (16.67%) was found.

CONCLUSIONS

These results, in agreement with previous data using different tissues markers, underlined the ethnic heterogeneity of the Mauritanian population.

Single molecule real-time DNA sequencing of HLA genes at ultra-high resolution from 126 International HLA and Immunogenetics Workshop cell lines

The hyperpolymorphic HLA genes play important roles in disease and transplantation and act as genetic markers of migration and evolution. A panel of 107 B-lymphoblastoid cell lines (B-LCLs) was established in 1987 at the 10th International Histocompatibility Workshop as a resource for the immunogenetics community. These B-LCLs are well characterised and represent diverse ethnicities and HLA haplotypes. Here we have applied Pacific Biosciences' Single Molecule Real-Time (SMRT) DNA sequencing to HLA type 126 B-LCL, including the 107 International HLA and Immunogenetics Workshop (IHIW) cells, to ultra-high resolution. Amplicon sequencing of full-length HLA class I genes (HLA-A, -B and -C) and partial length HLA class II genes (HLA-DRB1, -DQB1 and -DPB1) was performed. We typed a total of 931 HLA alleles, 895 (96%) of which were consistent with the typing in the IPD-IMGT/HLA Database (Release 3.27.0, January 20, 2017), with 595 (64%) typed at a higher resolution. Discrepant types, including novel alleles (n = 10) and changes in zygosity (n = 13), as well as previously unreported types (n = 34) were observed. In addition, patterns of linkage disequilibrium were distinguished by four-field resolution typing of HLA-B and HLA-C. By improving and standardising the HLA typing of these B-LCLs, we have ensured their continued usefulness as a resource for the immunogenetics community in the age of next generation DNA sequencing.

HLA-A, -B, -C and -DRB1 allele and haplotype frequencies in the Macedonian population based on a family study

AIM

The aim of this study was to determine HLA allele and 2-, 3- and 4-loci haplotype frequencies in a sample from Macedonian population with defined haplotypes based on family history.

MATERIAL AND METHODS

We analysed 286 unrelated individuals with Macedonian origin, parents of patients who needed stem cell transplantation, in the period of 01.01.2003 till 31.12.2016. Allele and haplotype frequencies, as well as Hardy-Weinberg equilibrium were calculated using the Arlequin3.5 software. Population comparison was calculated using the PHYLIP software.

RESULTS

We identified 18 HLA-A, 26 HLA-B, 13 HLA-C and 13 HLA-DRB1 allele group families. The most frequent allele groups in our population were HLA-A*02 (29.0%), HLA-A*24 (13.8%), HLA-B*35 (16.1%), HLA-B*51 (14.7%), HLA-B*18 (14.7%), HLA-C*07 (27.9%), HLA-DRB1*11 (25.5%) and HLA-DRB1*16 (14.8%). The most frequent four loci haplotype was HLA-A*01-B*08-C*07-DRB1*03 (2.7%). Our comparison showed that the Macedonian population is closely related to the neighbouring countries in the Balkan Peninsula.

CONCLUSION

This study provides data about the HLA diversity in the Macedonian population, which can be very important in the process of unrelated donor search, and in addition yields control group for future disease association studies in our population.

The HLA-B landscape of Africa: Signatures of pathogen-driven selection and molecular identification of candidate alleles to malaria protection

Human leukocyte antigen (HLA) genes play a key role in the immune response to infectious diseases, some of which are highly prevalent in specific environments, like malaria in sub-Saharan Africa. Former case-control studies showed that one particular HLA-B allele, B*53, was associated with malaria protection in Gambia, but this hypothesis was not tested so far within a population genetics framework. In this study, our objective was to assess whether pathogen-driven selection associated with malaria contributed to shape the HLA-B genetic landscape of Africa. To that aim, we first typed the HLA-A and -B loci in 484 individuals from 11 populations living in different environments across the Sahel, and we analysed these data together with those available for 29 other populations using several approaches including linear modelling on various genetic, geographic and environmental parameters. In addition to relevant signatures of populations' demography and migrations history in the genetic differentiation patterns of both HLA-A and -B loci, we found that the frequencies of three HLA alleles, B*53, B*78 and A*74, were significantly associated with Plasmodium falciparum malaria prevalence, suggesting their increase through pathogen-driven selection in malaria-endemic environments. The two HLA-B alleles were further identified, by high-throughput sequencing, as B*53:01:01 (in putative linkage disequilibrium with one HLA-C allele, C*04:01:01:01) and B*78:01 in all but one individuals tested, making them appropriate candidates to malaria protection. These results highlight the role of environmental factors in the evolution of the HLA polymorphism and open key perspectives for functional studies focusing on HLA peptide-binding properties.

High-Resolution HLA-A Typing in Normal Iranian Population

Background:

Human leukocyte antigen (HLA) gene is a highly polymorphic region. HLA typing is required to match patients and donors for transplantation; therefore, development of HLA registries is necessary for finding HLA match donors. HLA system is highly informative, and numerous studies have been conducted on HLA allele distribution in different populations.

Methods:

In this study, 100 unrelated Iranian individuals were typed for HLA-A locus using sequence-based typing method. Samples were subjected to the PCR, followed by Sanger sequencing and software analysis.

Results:

A*02:01 (13%) and A*24:02 (12%) were the two most frequent alleles, while A*01:14, A*02:05, A*02:11, A*02:34, A*02:50, A*11:04, A*23:02, A*24:34, A*25:01, A*26:09, A*26:43, A*29:67, A*30:54, A*31:02, A*31:66, A*32:03, A*32:04, A*33:03, and A*66:15 alleles had the least frequencies (1%).

Conclusion:

This is the first report of HLA-A allele level typing in a randomized population of Iran and can be useful for development of national registries of HLA-typed volunteer marrow donors and local cord blood banks.

Genotype-Specific Evolution of Hepatitis E Virus

Hepatitis E virus (HEV) is the most common cause of acute viral hepatitis globally. HEV comprises four genotypes with different geographic distributions and host ranges. We utilize this natural case-control study for investigating the evolution of zoonotic viruses compared to single-host viruses, using 244 near-full-length HEV genomes. Genome-wide estimates of the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS ratio) located a region of overlapping reading frames, which is subject to positive selection in genotypes 3 and 4. The open reading frames (ORFs) involved have functions related to host-pathogen interaction, so genotype-specific evolution of these regions may reflect their fitness. Bayesian inference of evolutionary rates shows that genotypes 3 and 4 have significantly higher rates than genotype 1 across all ORFs. Reconstruction of the phylogenies of zoonotic genotypes demonstrates significant intermingling of isolates between hosts. We speculate that the genotype-specific differences may result from cyclical adaptation to different hosts in genotypes 3 and 4.IMPORTANCE Hepatitis E virus (HEV) is increasingly recognized as a pathogen that affects both the developing and the developed world. While most often clinically mild, HEV can be severe or fatal in certain demographics, such as expectant mothers. Like many other viral pathogens, HEV has been classified into several distinct genotypes. We show that most of the HEV genome is evolutionarily constrained. One locus of positive selection is unusual in that it encodes two distinct protein products. We are the first to detect positive selection in this overlap region. Genotype 1, which infects humans only, appears to be evolving differently from genotypes 3 and 4, which infect multiple species, possibly because genotypes 3 and 4 are unable to achieve the same fitness due to repeated host jumps.

The HLA-A, -B and -DRB1 polymorphism in a large dataset of South Brazil bone marrow donors from Rio Grande do Sul

Human leukocyte antigen (HLA) genes are very informative in population genetics studies and their variability has been widely used to reconstruct the history of geographic and/or demographic expansions of human populations. The characterization of HLA diversity at the population level is also fundamental in clinical studies, particularly for bone marrow transplantation programs. In this study, we investigated the HLA molecular variation in Rio Grande do Sul, South Brazil, in order to identify possible regional differences across this state. More than 97,000 bone marrow donors were typed at the HLA- A, -B and -DRB1 loci and analyzed by considering two kinds of subdivisions based on both self-identified ethnicity and place of residence: (a) the official geographic subdivision defined by the Brazilian Institute of Geography and Statistics and (b) known information about the colonization history of the state. HLA allele and haplotype frequencies were estimated and compared among the defined subgroups. The results indicate a lack of correlation between genetic variation and geography and thus no clear HLA genetic structure based on geographic criteria. On the other hand, major differences were observed regarding ethnicity. In addition, local populations from Rio Grande do Sul were found to be genetically similar to their corresponding parental European populations from Germany, Italy and Portugal, as documented by historical data. Overall, this study provides a thorough characterization of the HLA genetic variation in Rio Grande do Sul and a better understanding of its demographic history, being most useful for the development of more efficient strategies in bone marrow donors' recruitment.

Mapping the HLA diversity of the Iberian Peninsula

The polymorphism of HLA genes can be used to reconstruct human peopling history. However, this huge diversity impairs successful matching in stem cell transplantation, a situation which has led to the recruitment of millions of donors worldwide. In parallel to the increase of recruitment, registries are progressively relying on information from population genetics to optimize their donor pools in terms of HLA variability. In this study, the HLA data of 65,000 Spanish bone marrow donors were analyzed together with 60,000 Portuguese individuals to provide a comprehensive HLA genetic map of the Iberian Peninsula. The frequencies of many alleles were shown to vary continuously across the Peninsula, either increasing or decreasing from the Mediterranean coast to the Atlantic domain or from the Strait of Gibraltar to the Pyrenees and Bay of Biscay. Similar patterns were observed for several haplotypes. In addition, within some regions neighboring provinces share a close genetic similarity. These results outline the genetic landscape of the Iberian Peninsula, and confirm that the analysis of the HLA polymorphism may reveal relevant signatures of past demographic events even when data from donor registries are used. This conclusion stimulates future developments of the Spanish registry, presented here for the first time.

Evidence of extensive non-allelic gene conversion among LTR elements in the human genome

Long Terminal Repeats (LTRs) are nearly identical DNA sequences found at either end of Human Endogenous Retroviruses (HERVs). The high sequence similarity that exists among different LTRs suggests they could be substrate of ectopic gene conversion events. To understand the extent to which gene conversion occurs and to gain new insights into the evolutionary history of these elements in humans, we performed an intra-species phylogenetic study of 52 LTRs on different unrelated Y chromosomes. From this analysis, we obtained direct evidence that demonstrates the occurrence of ectopic gene conversion in several LTRs, with donor sequences located on both sex chromosomes and autosomes. We also found that some of these elements are characterized by an extremely high density of polymorphisms, showing one of the highest nucleotide diversities in the human genome, as well as a complex patchwork of sequences derived from different LTRs. Finally, we highlighted the limits of current short-read NGS studies in the analysis of genetic diversity of the LTRs in the human genome. In conclusion, our comparative re-sequencing analysis revealed that ectopic gene conversion is a common event in the evolution of LTR elements, suggesting complex genetic links among LTRs from different chromosomes.

HLA class I molecular variation and peptide-binding properties suggest a model of joint divergent asymmetric selection

The main function of HLA class I molecules is to present pathogen-derived peptides to cytotoxic T lymphocytes. This function is assumed to drive the maintenance of an extraordinary amount of polymorphism at each HLA locus, providing an immune advantage to heterozygote individuals capable to present larger repertories of peptides than homozygotes. This seems contradictory, however, with a reduced diversity at individual HLA loci exhibited by some isolated populations. This study shows that the level of functional diversity predicted for the two HLA-A and HLA-B genes considered simultaneously is similar (almost invariant) between 46 human populations, even when a reduced diversity exists at each locus. We thus propose that HLA-A and HLA-B evolved through a model of joint divergent asymmetric selection conferring all populations an equivalent immune potential. The distinct pattern observed for HLA-C is explained by its functional evolution towards killer cell immunoglobulin-like receptor (KIR) activity regulation rather than peptide presentation.

Forward-in-Time, Spatially Explicit Modeling Software to Simulate Genetic Lineages Under Selection

SELECTOR is a software package for studying the evolution of multiallelic genes under balancing or positive selection while simulating complex evolutionary scenarios that integrate demographic growth and migration in a spatially explicit population framework. Parameters can be varied both in space and time to account for geographical, environmental, and cultural heterogeneity. SELECTOR can be used within an approximate Bayesian computation estimation framework. We first describe the principles of SELECTOR and validate the algorithms by comparing its outputs for simple models with theoretical expectations. Then, we show how it can be used to investigate genetic differentiation of loci under balancing selection in interconnected demes with spatially heterogeneous gene flow. We identify situations in which balancing selection reduces genetic differentiation between population groups compared with neutrality and explain conflicting outcomes observed for human leukocyte antigen loci. These results and three previously published applications demonstrate that SELECTOR is efficient and robust for building insight into human settlement history and evolution.

Impact of host genetic polymorphisms on vaccine induced antibody response

Many host- and vaccine-specific factors modulate an antibody response. Host genetic polymorphisms, in particular, modulate the immune response in multiple ways on different scales. This review article describes how information on host genetic polymorphisms and corresponding immune cascades may be used to generate personalized vaccine strategies to optimize the antibody response.

Variation in NAT2 acetylation phenotypes is associated with differences in food-producing subsistence modes and ecoregions in Africa

Background

Dietary changes associated to shifts in subsistence strategies during human evolution may have induced new selective pressures on phenotypes, as currently held for lactase persistence. Similar hypotheses exist for arylamine N-acetyltransferase 2 (NAT2) mediated acetylation capacity, a well-known pharmacogenetic trait with wide inter-individual variation explained by polymorphisms in the NAT2 gene. The environmental causative factor (if any) driving its evolution is as yet unknown, but significant differences in prevalence of acetylation phenotypes are found between hunter-gatherer and food-producing populations, both in sub-Saharan Africa and worldwide, and between agriculturalists and pastoralists in Central Asia. These two subsistence strategies also prevail among sympatric populations of the African Sahel, but knowledge on NAT2 variation among African pastoral nomads was up to now very scarce. Here we addressed the hypothesis of different selective pressures associated to the agriculturalist or pastoralist lifestyles having acted on the evolution of NAT2 by sequencing the gene in 287 individuals from five pastoralist and one agriculturalist Sahelian populations.

Results

We show that the significant NAT2 genetic structure of African populations is mainly due to frequency differences of three major haplotypes, two of which are categorized as decreased function alleles (NAT2*5B and NAT2*6A), particularly common in populations living in arid environments, and one fast allele (NAT2*12A), more frequently detected in populations living in tropical humid environments. This genetic structure does associate more strongly with a classification of populations according to ecoregions than to subsistence strategies, mainly because most Sahelian and East African populations display little to no genetic differentiation between them, although both regions hold nomadic or semi-nomadic pastoralist and sedentary agriculturalist communities. Furthermore, we found significantly higher predicted proportions of slow acetylators in pastoralists than in agriculturalists, but also among food-producing populations living in the Sahelian and dry savanna zones than in those living in humid environments, irrespective of their mode of subsistence.

Conclusion

Our results suggest a possible independent influence of both the dietary habits associated with subsistence modes and the chemical environment associated with climatic zones and biomes on the evolution of NAT2 diversity in sub-Saharan African populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0543-6) contains supplementary material, which is available to authorized users.

Intrahaplotypic Variants Differentiate Complex Linkage Disequilibrium within Human MHC Haplotypes

Distinct regions of long-range genetic fixation in the human MHC region, known as conserved extended haplotypes (CEHs), possess unique genomic characteristics and are strongly associated with numerous diseases. While CEHs appear to be homogeneous by SNP analysis, the nature of fine variations within their genomic structure is unknown. Using multiple, MHC-homozygous cell lines, we demonstrate extensive sequence conservation in two common Asian MHC haplotypes: A33-B58-DR3 and A2-B46-DR9. However, characterization of phase-resolved MHC haplotypes revealed unique intra-CEH patterns of variation and uncovered 127 single nucleotide variants (SNVs) which are missing from public databases. We further show that the strong linkage disequilibrium structure within the human MHC that typically confounds precise identification of genetic features can be resolved using intra-CEH variants, as evidenced by rs3129063 and rs448489, which affect expression of ZFP57, a gene important in methylation and epigenetic regulation. This study demonstrates an improved strategy that can be used towards genetic dissection of diseases.

Computer simulation of human leukocyte antigen genes supports two main routes of colonization by human populations in East Asia

Background

Recent genetic studies have suggested that the colonization of East Asia by modern humans was more complex than a single origin from the South, and that a genetic contribution via a Northern route was probably quite substantial.

Results

Here we use a spatially-explicit computer simulation approach to investigate the human migration hypotheses of this region based on one-route or two-route models. We test the likelihood of each scenario by using Human Leukocyte Antigen (HLA) − A, −B, and − DRB1 genetic data of East Asian populations, with both selective and demographic parameters considered. The posterior distribution of each parameter is estimated by an Approximate Bayesian Computation (ABC) approach.

Conclusions

Our results strongly support a model with two main routes of colonization of East Asia on both sides of the Himalayas, with distinct demographic histories in Northern and Southern populations, characterized by more isolation in the South. In East Asia, gene flow between populations originating from the two routes probably existed until a remote prehistoric period, explaining the continuous pattern of genetic variation currently observed along the latitude. A significant although dissimilar level of balancing selection acting on the three HLA loci is detected, but its effect on the local genetic patterns appears to be minor compared to those of past demographic events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0512-0) contains supplementary material, which is available to authorized users.

HLA supertype variation across populations: new insights into the role of natural selection in the evolution of HLA-A and HLA-B polymorphisms

Supertypes are groups of human leukocyte antigen (HLA) alleles which bind overlapping sets of peptides due to sharing specific residues at the anchor positions-the B and F pockets-of the peptide-binding region (PBR). HLA alleles within the same supertype are expected to be functionally similar, while those from different supertypes are expected to be functionally distinct, presenting different sets of peptides. In this study, we applied the supertype classification to the HLA-A and HLA-B data of 55 worldwide populations in order to investigate the effect of natural selection on supertype rather than allelic variation at these loci. We compared the nucleotide diversity of the B and F pockets with that of the other PBR regions through a resampling procedure and compared the patterns of within-population heterozygosity (He) and between-population differentiation (G ST) observed when using the supertype definition to those estimated when using randomized groups of alleles. At HLA-A, low levels of variation are observed at B and F pockets and randomized He and G ST do not differ from the observed data. By contrast, HLA-B concentrates most of the differences between supertypes, the B pocket showing a particularly high level of variation. Moreover, at HLA-B, the reassignment of alleles into random groups does not reproduce the patterns of population differentiation observed with supertypes. We thus conclude that differently from HLA-A, for which supertype and allelic variation show similar patterns of nucleotide diversity within and between populations, HLA-B has likely evolved through specific adaptations of its B pocket to local pathogens.

A limit to the divergent allele advantage model supported by variable pathogen recognition across HLA-DRB1 allele lineages

Genetic diversity in human leukocyte antigen (HLA) molecules is thought to have arisen from the co‐evolution between host and pathogen and maintained by balancing selection. Heterozygote advantage is a common proposed scenario for maintaining high levels of diversity in HLA genes, and extending from this, the divergent allele advantage (DAA) model suggests that individuals with more divergent HLA alleles bind and recognize a wider array of antigens. While the DAA model seems biologically suitable for driving HLA diversity, there is likely an upper threshold to the amount of sequence divergence. We used peptide‐binding and pathogen‐recognition capacity of DRB1 alleles as a model to further explore the DAA model; within the DRB1 locus, we examined binding predictions based on two distinct phylogenetic groups (denoted group A and B) previously identified based on non‐peptide‐binding region (PBR) nucleotide sequences. Predictions in this study support that group A allele and group B allele lineages have contrasting binding/recognition capacity, with only the latter supporting the DAA model. Furthermore, computer simulations revealed an inconsistency in the DAA model alone with observed extent of polymorphisms, supporting that the DAA model could only work effectively in combination with other mechanisms. Overall, we support that the mechanisms driving HLA diversity are non‐exclusive. By investigating the relationships among HLA alleles, and pathogens recognized, we can provide further insights into the mechanisms on how humans have adapted to infectious diseases over time.

Human Leukocyte Antigen Diversity: A Southern African Perspective

Despite the increasingly well-documented evidence of high genetic, ethnic, and linguistic diversity amongst African populations, there is limited data on human leukocyte antigen (HLA) diversity in these populations. HLA is part of the host defense mechanism mediated through antigen presentation to effector cells of the immune system. With the high disease burden in southern Africa, HLA diversity data is increasingly important in the design of population-specific vaccines and the improvement of transplantation therapeutic interventions. This review highlights the paucity of HLA diversity data amongst southern African populations and defines a need for information of this kind. This information will support disease association studies, provide guidance in vaccine design, and improve transplantation outcomes.

A genome-wide association study of myasthenia gravis

IMPORTANCE

Myasthenia gravis is a chronic, autoimmune, neuromuscular disease characterized by fluctuating weakness of voluntary muscle groups. Although genetic factors are known to play a role in this neuroimmunological condition, the genetic etiology underlying myasthenia gravis is not well understood.

OBJECTIVE

To identify genetic variants that alter susceptibility to myasthenia gravis, we performed a genome-wide association study.

DESIGN, SETTING, AND PARTICIPANTS

DNA was obtained from 1032 white individuals from North America diagnosed as having acetylcholine receptor antibody-positive myasthenia gravis and 1998 race/ethnicity-matched control individuals from January 2010 to January 2011. These samples were genotyped on Illumina OmniExpress single-nucleotide polymorphism arrays. An independent cohort of 423 Italian cases and 467 Italian control individuals were used for replication.

MAIN OUTCOMES AND MEASURES

We calculated P values for association between 8,114,394 genotyped and imputed variants across the genome and risk for developing myasthenia gravis using logistic regression modeling. A threshold P value of 5.0×10(-8) was set for genome-wide significance after Bonferroni correction for multiple testing.

RESULTS

In the overall case-control cohort, we identified association signals at CTLA4 (rs231770; P=3.98×10(-8); odds ratio, 1.37; 95% CI, 1.25-1.49), HLA-DQA1 (rs9271871; P=1.08×10(-8); odds ratio, 2.31; 95% CI, 2.02-2.60), and TNFRSF11A (rs4263037; P=1.60×10(-9); odds ratio, 1.41; 95% CI, 1.29-1.53). These findings replicated for CTLA4 and HLA-DQA1 in an independent cohort of Italian cases and control individuals. Further analysis revealed distinct, but overlapping, disease-associated loci for early- and late-onset forms of myasthenia gravis. In the late-onset cases, we identified 2 association peaks: one was located in TNFRSF11A (rs4263037; P=1.32×10(-12); odds ratio, 1.56; 95% CI, 1.44-1.68) and the other was detected in the major histocompatibility complex on chromosome 6p21 (HLA-DQA1; rs9271871; P=7.02×10(-18); odds ratio, 4.27; 95% CI, 3.92-4.62). Association within the major histocompatibility complex region was also observed in early-onset cases (HLA-DQA1; rs601006; P=2.52×10(-11); odds ratio, 4.0; 95% CI, 3.57-4.43), although the set of single-nucleotide polymorphisms was different from that implicated among late-onset cases.

CONCLUSIONS AND RELEVANCE

Our genetic data provide insights into aberrant cellular mechanisms responsible for this prototypical autoimmune disorder. They also suggest that clinical trials of immunomodulatory drugs related to CTLA4 and that are already Food and Drug Administration approved as therapies for other autoimmune diseases could be considered for patients with refractory disease.

HLA genetic diversity in Hungarians and Hungarian Gypsies: complementary differentiation patterns and demographic signals revealed by HLA-A, -B and -DRB1 in Central Europe

Systematic analyses of human leukocyte antigen (HLA) profiles in different populations may increase the efficiency of bone marrow donor selection and help reconstructing human peopling history. We typed HLA-A, -B, and -DRB1 allele groups in two bone marrow donor cohorts of 2402 Hungarians and 186 Hungarian Gypsies and compared them with several Central-European, Spanish Gypsy, and Indian populations. Our results indicate that different European Gypsy populations share a common origin but diverged genetically as a consequence of founder effect and rapid genetic drift, whereas other European populations are related genetically in relation to geography. This study also suggests that while HLA-A accurately depicts the effects of genetic drift, HLA-B, and -DRB1 conserve more signatures of ancient population relationships, as a result of balancing selection.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

HLA-A, B and DRB1 genetic heterogeneity in Quebec

Recent studies have shown that under specific conditions such as high sample sizes and Hardy-Weinberg equilibrium, bone marrow donor registry data can be used to describe HLA molecular variation across a specific geographic area, thus providing excellent data sets to infer human migrations history. The province of Quebec is known to have experienced a complex history of settlement, characterized by multiple migrations and demographic changes. We thus analysed the data of more than 13 000 unrelated individuals acting as volunteer bone marrow donors who were molecularly typed for HLA-A, B and DRB1 polymorphisms in the Héma-Quebec registry. HLA allelic and haplotypic frequencies were estimated and compared among regions. The results indicate that, despite an overall low genetic diversity in Quebec, genetic variation is correlated with geography, compatible with isolation-by-distance across the province. However, some localities also harbour contrasting genetic profiles, that is a highly diversified genetic pool in the two main urban centres (Montréal and Laval) and a more pronounced genetic divergence of two specific regions characterized by a peculiar peopling history (Saguenay-Lac-St-Jean and Gaspésie-Îles-De-La-Madeleine). In agreement with other independent molecular markers, the observations based on HLA data thus account for the main demographic mechanisms that shaped the genetic structure of the present day Quebecer population. In addition, the detailed analysis of the Héma-Quebec registry provides key genetic information on which an efficient bone marrow transplantation recruitment strategy can be settled.

Seven novel HLA alleles reflect different mechanisms involved in the evolution of HLA diversity: description of the new alleles and review of the literature

The human leukocyte antigen (HLA) loci are among the most polymorphic genes in the human genome. The diversity of these genes is thought to be generated by different mechanisms including point mutation, gene conversion and crossing-over. During routine HLA typing, we discovered seven novel HLA alleles which were probably generated by different evolutionary mechanisms. HLA-B*41:21, HLA-DQB1*02:10 and HLA-DQA1*01:12 likely emerged from the common alleles of their groups by point mutations, all of which caused non-synonymous amino acid substitutions. In contrast, a deletion of one nucleotide leading to a frame shift with subsequent generation of a stop codon is responsible for the appearance of a null allele, HLA-A*01:123N. Whereas HLA-B*35:231 and HLA-B*53:31 were probably products of intralocus gene conversion between HLA-B alleles, HLA-C*07:294 presumably evolved by interlocus gene conversion between an HLA-C and an HLA-B allele. Our analysis of these novel alleles illustrates the different mechanisms which may have contributed to the evolution of HLA polymorphism.

The relevance of HLA sequencing in population genetics studies

Next generation sequencing (NGS) is currently being adapted by different biotechnological platforms to the standard typing method for HLA polymorphism, the huge diversity of which makes this initiative particularly challenging. Boosting the molecular characterization of the HLA genes through efficient, rapid, and low-cost technologies is expected to amplify the success of tissue transplantation by enabling us to find donor-recipient matching for rare phenotypes. But the application of NGS technologies to the molecular mapping of the MHC region also anticipates essential changes in population genetic studies. Huge amounts of HLA sequence data will be available in the next years for different populations, with the potential to change our understanding of HLA variation in humans. In this review, we first explain how HLA sequencing allows a better assessment of the HLA diversity in human populations, taking also into account the methodological difficulties it introduces at the statistical level; secondly, we show how analyzing HLA sequence variation may improve our comprehension of population genetic relationships by facilitating the identification of demographic events that marked human evolution; finally, we discuss the interest of both HLA and genome-wide sequencing and genotyping in detecting functionally significant SNPs in the MHC region, the latter having also contributed to the makeup of the HLA molecular diversity observed today.

Performance of HLA allele prediction methods in African Americans for class II genes HLA-DRB1, -DQB1, and -DPB1

Background

The expense of human leukocyte antigen (HLA) allele genotyping has motivated the development of imputation methods that use dense single nucleotide polymorphism (SNP) genotype data and the region’s haplotype structure, but the performance of these methods in admixed populations (such as African Americans) has not been adequately evaluated. We compared genotype-based—derived from both genome-wide genotyping and targeted sequencing—imputation results to existing allele data for HLA–DRB1, −DQB1, and –DPB1.

Results

In European Americans, the newly-developed HLA Genotype Imputation with Attribute Bagging (HIBAG) method outperformed HLA*IMP:02. In African Americans, HLA*IMP:02 performed marginally better than HIBAG pre-built models, but HIBAG models constructed using a portion of our African American sample with both SNP genotyping and four-digit HLA class II allele typing had consistently higher accuracy than HLA*IMP:02. However, HIBAG was significantly less accurate in individuals heterozygous for local ancestry (p ≤0.04). Accuracy improved in models with equal numbers of African and European chromosomes. Variants added by targeted sequencing and SNP imputation further improved both imputation accuracy and the proportion of high quality calls.

Conclusion

Combining the HIBAG approach with local ancestry and dense variant data can produce highly-accurate HLA class II allele imputation in African Americans.

HLA variation reveals genetic continuity rather than population group structure in East Asia

Genetic differences between Northeast Asian (NEA) and Southeast Asian (SEA) populations have been observed in numerous studies. At the among-population level, despite a clear north-south differentiation observed for many genetic markers, debates were led between abrupt differences and a continuous pattern. At the within-population level, whether NEA or SEA populations have higher genetic diversity is also highly controversial. In this study, we analyzed a large set of HLA data from East Asia in order to map the genetic variation among and within populations in this continent and to clarify the distribution pattern of HLA lineages and alleles. We observed a genetic differentiation between NEA and SEA populations following a continuous pattern from north to south, and we show a significant and continuous decrease of HLA diversity by the same direction. This continuity is shaped by clinal distributions of many HLA lineages and alleles with increasing or decreasing frequencies along the latitude. These results bring new evidence in favor of the "overlapping model" proposed previously for East Asian peopling history, whereby modern humans migrated eastward from western Eurasia via two independent routes along each side of the Himalayas and, later, overlapped in East Asia across open land areas. Our study strongly suggests that intensive gene flow between NEA and SEA populations occurred and shaped the latitude-related continuous pattern of genetic variation and the peculiar HLA lineage and allele distributions observed in this continent. Probably for a very long period, the exact duration of these events remains to be estimated.

Population-specific differences in gene conversion patterns between human SUZ12 and SUZ12P are indicative of the dynamic nature of interparalog gene conversion

Nonallelic homologous gene conversion (NAHGC) resulting from interparalog recombination without crossover represents an important influence on the evolution of duplicated sequences in the human genome. In 17q11.2, different paralogous sequences mediate large NF1 deletions by nonallelic homologous recombination with crossover (NAHR). Among these paralogs are SUZ12 and its pseudogene SUZ12P which harbour the breakpoints of type-2 (1.2-Mb) NF1 deletions. Such deletions are caused predominantly by mitotic NAHR since somatic mosaicism with normal cells is evident in most patients. Investigating whether SUZ12 and SUZ12P have also been involved in NAHGC, we observed gene conversion tracts between these paralogs in both Africans (AFR) and Europeans (EUR). Since germline type-2 NF1 deletions resulting from meiotic NAHR are very rare, the vast majority of the gene conversion tracts in SUZ12 and SUZ12P are likely to have resulted from mitotic recombination during premeiotic cell divisions of germ cells. A higher number of gene conversion tracts were noted within SUZ12 and SUZ12P in AFR as compared to EUR. Further, the distinctive signature of NAHGC (a high number of SNPs per paralog and a high number of shared SNPs between paralogs), a characteristic of many actively recombining paralogs, was observed in both SUZ12 and SUZ12P but only in AFR and not in EUR. A novel polymorphic 2.3-kb deletion in SUZ12P was identified which exhibited a high allele frequency in EUR. We postulate that this interparalog structural difference, together with low allelic recombination rates, could have caused a reduction in NAHGC between SUZ12 and SUZ12P during human evolution.

Assessing the HLA diversity of cord blood units collected from a birth clinic caring for pregnant women in an ethnically diverse metropolitan area

BACKGROUND

New strategies are emerging in cord blood banking where focusing on birth clinics caring for a high number of mothers belonging to ethnic minorities could offer new possibilities for allotransplantation both for patients of European origin and for patients from ethnic minorities or mixed ancestries.

STUDY DESIGN AND METHODS

Marseilles Cord Blood Bank works with one university birth clinic caring for a culturally and sociologically diverse population. Stringent French legal restrictions apply to recording the geographic origin of parents. To circumvent this limitation and evaluate the contribution of newly banked cord blood units (CBUs) to increasing HLA diversity, we applied an algorithm that allows for the determination of parents' putative haplotypes and thus grossly deduce information on their ancestry. Generic resolution HLA-A, HLA-B, and allelic resolution HLA-DRB1 genotyping for 328 CBUs and 2691 unrelated donors (UDs) between January 2009 and May 2012 were performed. Enrichment from international CBU registry with nonreferenced generic HLA-A, HLA-B, and allelic HLA-DRB1 phenotypes was compared between CBUs identified with one or two non-European haplotypes and CBUs identified with two European haplotypes.

RESULTS

Marseilles CBUs display an increased proportion of HLA antigens frequently expressed in African populations compared to UDs. Whereas 93% of 199 CBUs identified with one or two non-European haplotypes enrich international CBU registry, this result is reduced to 42% for the 129 CBUs identified with two European haplotypes.

CONCLUSION

This study supports a new method to assess HLA diversity. However, such an increased of HLA diversity raises questions about frequencies of CBUs released and clinical relevance from their uses.

Human leukocyte antigen DQ2/8 prevalence in non-celiac patients with gastrointestinal diseases

AIM: To investigate the prevalence of human leukocyte antigen (HLA) DQ2/8 alleles in Southern Italians with liver and gastrointestinal (GI) diseases outside of celiac disease.

METHODS: HLA DQ2/8 status was assessed in 443 patients from three ambulatory gastroenterology clinics in Southern Italy (University of Federico II, Naples, Loreto Crispi Hospital, Ruggi D’Aragona Hospital, Salerno). Patients were grouped based on disease status [pre-post transplant liver disease, esophageal/gastric organic and functional diseases, irritable bowel syndrome (IBS) and inflammatory bowel disease (IBD)] and DQ2/8 alleles, which correspond to a celiac disease genetic risk gradient. Subject allele frequencies were compared to healthy Italian controls.

RESULTS: One hundred and ninety-six out of four hundred and forty-three (44.2%) subjects, median age 56 years and 42.6% female, were DQ2/8 positive. When stratifying by disease we found that 86/188 (45.7%) patients with liver disease were HLA DQ2/8 positive, 39/73 (53.4%) with functional upper GI diseases and 19/41 (46.3%) with organic upper GI diseases were positive. Furthermore, 38/105 (36.2%) patients with IBS and 14/36 (38.9%) with IBD were HLA DQ2/8 positive (P = 0.21). Compared to healthy controls those with functional upper GI diseases disorders had a 1.8 times higher odds of DQ2/8 positivity. Those with liver disease had 1.3 times the odds, albeit not statistically significant, of DQ2/8 positivity. Both those with IBS and IBD had a lower odds of DQ2/8 positivity compared to healthy controls.

CONCLUSION: The proportion of individuals HLA DQ2/8 positive is higher in those with liver/upper functional GI disease and lower in IBS/IBD as compared to general population estimates.

TNFA promoter alleles--frequencies and linkage with classical HLA genes in a South German Caucasian population

The tumor necrosis factor alpha (TNFA) promoter region exhibits several polymorphisms, which have been hypothesized to influence gene expression, thereby associating positively or negatively with inflammatory conditions. Many studies have focused on single nucleotide polymorphisms (SNPs) taking not into account additive or inverse effects between different SNPs. We typed 1,021 healthy Caucasian volunteer stem cell donors for their TNFA promoter as well as their HLA-A,-C,-B,-DRB1 loci. Using statistical methods, we reconstructed TNFA promoter alleles and analyzed their frequency and linkage with HLA genes. We show that the number of TNFA promoter alleles frequent enough to be analyzed in clinical studies is limited and that a strong linkage with classical HLA genes is present, especially for the extended HLA-haplotype HLA-A*01:01/HLA-C*07:01/HLA-B*08:01/TNFA promoter allele 3/HLA*DRB1*03:01. Taking into account SNP frequency information, it is possible to quite accurately deduce TNFA promoter alleles by generic Sanger sequencing, obviating the need for elaborating allele-specific sequencing. This information may enable investigators to consider the complete TNFA regulatory region in a phase-separated manner in contrast to previous approaches examining only one or few isolated SNPs.

16(th) IHIW: analysis of HLA population data, with updated results for 1996 to 2012 workshop data (AHPD project report)

We present here the results of the Analysis of HLA Population Data (AHPD) project of the 16th International HLA and Immunogenetics Workshop (16IHIW) held in Liverpool in May-June 2012. Thanks to the collaboration of 25 laboratories from 18 different countries, HLA genotypic data for 59 new population samples (either well-defined populations or donor registry samples) were gathered and 55 were analysed statistically following HLA-NET recommendations. The new data included, among others, large sets of well-defined populations from north-east Europe and West Asia, as well as many donor registry data from European countries. The Gene[rate] computer tools were combined to create a Gene[rate] computer pipeline to automatically (i) estimate allele frequencies by an expectation-maximization algorithm accommodating ambiguities, (ii) estimate heterozygosity, (iii) test for Hardy-Weinberg equilibrium (HWE), (iv) test for selective neutrality, (v) generate frequency graphs and summary statistics for each sample at each locus and (vi) plot multidimensional scaling (MDS) analyses comparing the new samples with previous IHIW data. Intrapopulation analyses show that HWE is rarely rejected, while neutrality tests often indicate a significant excess of heterozygotes compared with neutral expectations. The comparison of the 16IHIW AHPD data with data collected during previous workshops (12th-15th) shows that geography is an excellent predictor of HLA genetic differentiations for HLA-A, -B and -DRB1 loci but not for HLA-DQ, whose patterns are probably more influenced by natural selection. In Europe, HLA genetic variation clearly follows a north to south-east axis despite a low level of differentiation between European, North African and West Asian populations. Pacific populations are genetically close to Austronesian-speaking South-East Asian and Taiwanese populations, in agreement with current theories on the peopling of Oceania. Thanks to this project, HLA genetic variation is more clearly defined worldwide and better interpreted in relation to human peopling history and HLA molecular evolution.

A sequence-based approach demonstrates that balancing selection in classical human leukocyte antigen (HLA) loci is asymmetric

Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1, HLA-B and HLA-DPA1, with significant deviation (P ≤ 1.1 × 10(-4)) in about half of the populations. There were significant results at all loci except HLA-DQB1/DQA1. We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

The heterogeneous HLA genetic makeup of the Swiss population

This study aims at investigating the HLA molecular variation across Switzerland in order to determine possible regional differences, which would be highly relevant to several purposes: optimizing donor recruitment strategies in hematopoietic stem cell transplantation (HSCT), providing reliable reference data in HLA and disease association studies, and understanding the population genetic background(s) of this culturally heterogeneous country. HLA molecular data of more than 20,000 HSCT donors from 9–13 recruitment centers of the whole country were analyzed. Allele and haplotype frequencies were estimated by using new computer tools adapted to the heterogeneity and ambiguity of the data. Non-parametric and resampling statistical tests were performed to assess Hardy-Weinberg equilibrium, selective neutrality and linkage disequilibrium among different loci, both in each recruitment center and in the whole national registry. Genetic variation was explored through genetic distance and hierarchical analysis of variance taking into account both geographic and linguistic subdivisions in Switzerland. The results indicate a heterogeneous genetic makeup of the Swiss population: first, allele frequencies estimated on the whole national registry strongly deviate from Hardy-Weinberg equilibrium, by contrast with the results obtained for individual centers; second, a pronounced differentiation is observed for Ticino, Graubünden, and, to a lesser extent, Wallis, suggesting that the Alps represent(ed) a barrier to gene flow; finally, although cultural (linguistic) boundaries do not represent a main genetic differentiation factor in Switzerland, the genetic relatedness between population from south-eastern Switzerland and Italy agrees with historical and linguistic data. Overall, this study justifies the maintenance of a decentralized donor recruitment structure in Switzerland allowing increasing the genetic diversity of the national—and hence global—donor registry. It also indicates that HLA data of local donor recruitment centers can be used as reference data in both epidemiological and population genetic studies focusing on the genetic history of present European populations.

Influence of HLA-DRB1 and HLA-DQB1 alleles on IgG antibody response to the P. vivax MSP-1, MSP-3α and MSP-9 in individuals from Brazilian endemic area

BACKGROUND

The antibody response generated during malaria infections is of particular interest, since the production of specific IgG antibodies is required for acquisition of clinical immunity. However, variations in antibody responses could result from genetic polymorphism of the HLA class II genes. Given the increasing focus on the development of subunit vaccines, studies of the influence of class II alleles on the immune response in ethnically diverse populations is important, prior to the implementation of vaccine trials.

METHODS AND FINDINGS

In this study, we evaluated the influence of HLA-DRB1* and -DQB1* allelic groups on the naturally acquired humoral response from Brazilian Amazon individuals (n = 276) against P. vivax Merozoite Surface Protein-1 (MSP-1), MSP-3α and MSP-9 recombinant proteins. Our results provide information concerning these three P. vivax antigens, relevant for their role as immunogenic surface proteins and vaccine candidates. Firstly, the studied population was heterogeneous presenting 13 HLA-DRB1* and 5 DQB1* allelic groups with a higher frequency of HLA-DRB1*04 and HLA-DQB1*03. The proteins studied were broadly immunogenic in a naturally exposed population with high frequency of IgG antibodies against PvMSP1-19 (86.7%), PvMSP-3 (77%) and PvMSP-9 (76%). Moreover, HLA-DRB1*04 and HLA-DQB1*03 alleles were associated with a higher frequency of IgG immune responses against five out of nine antigens tested, while HLA-DRB1*01 was associated with a high frequency of non-responders to repetitive regions of PvMSP-9, and the DRB1*16 allelic group with the low frequency of responders to PvMSP3 full length recombinant protein.

CONCLUSIONS

HLA-DRB1*04 alleles were associated with high frequency of antibody responses to five out of nine recombinant proteins tested in Rondonia State, Brazil. These features could increase the success rate of future clinical trials based on these vaccine candidates.

Evaluation of the imputation performance of the program IMPUTE in an admixed sample from Mexico City using several model designs

BACKGROUND

We explored the imputation performance of the program IMPUTE in an admixed sample from Mexico City. The following issues were evaluated: (a) the impact of different reference panels (HapMap vs. 1000 Genomes) on imputation; (b) potential differences in imputation performance between single-step vs. two-step (phasing and imputation) approaches; (c) the effect of different INFO score thresholds on imputation performance and (d) imputation performance in common vs. rare markers.

METHODS

The sample from Mexico City comprised 1,310 individuals genotyped with the Affymetrix 5.0 array. We randomly masked 5% of the markers directly genotyped on chromosome 12 (n=1,046) and compared the imputed genotypes with the microarray genotype calls. Imputation was carried out with the program IMPUTE. The concordance rates between the imputed and observed genotypes were used as a measure of imputation accuracy and the proportion of non-missing genotypes as a measure of imputation efficacy.

RESULTS

The single-step imputation approach produced slightly higher concordance rates than the two-step strategy (99.1% vs. 98.4% when using the HapMap phase II combined panel), but at the expense of a lower proportion of non-missing genotypes (85.5% vs. 90.1%). The 1,000 Genomes reference sample produced similar concordance rates to the HapMap phase II panel (98.4% for both datasets, using the two-step strategy). However, the 1000 Genomes reference sample increased substantially the proportion of non-missing genotypes (94.7% vs. 90.1%). Rare variants (<1%) had lower imputation accuracy and efficacy than common markers.

CONCLUSIONS

The program IMPUTE had an excellent imputation performance for common alleles in an admixed sample from Mexico City, which has primarily Native American (62%) and European (33%) contributions. Genotype concordances were higher than 98.4% using all the imputation strategies, in spite of the fact that no Native American samples are present in the HapMap and 1000 Genomes reference panels. The best balance of imputation accuracy and efficiency was obtained with the 1,000 Genomes panel. Rare variants were not captured effectively by any of the available panels, emphasizing the need to be cautious in the interpretation of association results for imputed rare variants.

Strategies to work with HLA data in human populations for histocompatibility, clinical transplantation, epidemiology and population genetics: HLA-NET methodological recommendations

HLA-NET (a European COST Action) aims at networking researchers working in bone marrow transplantation, epidemiology and population genetics to improve the molecular characterization of the HLA genetic diversity of human populations, with an expected strong impact on both public health and fundamental research. Such improvements involve finding consensual strategies to characterize human populations and samples and report HLA molecular typings and ambiguities; proposing user-friendly access to databases and computer tools and defining minimal requirements related to ethical aspects. The overall outcome is the provision of population genetic characterizations and comparisons in a standard way by all interested laboratories. This article reports the recommendations of four working groups (WG1-4) of the HLA-NET network at the mid-term of its activities. WG1 (Population definitions and sampling strategies for population genetics’ analyses) recommends avoiding outdated racial classifications and population names (e.g. ‘Caucasian’) and using instead geographic and/or cultural (e.g. linguistic) criteria to describe human populations (e.g. ‘pan-European’). A standard ‘HLA-NET POPULATION DATA QUESTIONNAIRE’ has been finalized and is available for the whole HLA community. WG2 (HLA typing standards for population genetics analyses) recommends retaining maximal information when reporting HLA typing results. Rather than using the National Marrow Donor Program coding system, all ambiguities should be provided by listing all allele pairs required to explain each genotype, according to the formats proposed in ‘HLA-NET GUIDELINES FOR REPORTING HLA TYPINGS’. The group also suggests taking into account a preliminary list of alleles defined by polymorphisms outside the peptide-binding sites that may affect population genetic statistics because of significant frequencies. WG3 (Bioinformatic strategies for HLA population data storage and analysis) recommends the use of programs capable of dealing with ambiguous data, such as the ‘gene[rate]’ computer tools to estimate frequencies, test for Hardy–Weinberg equilibrium and selective neutrality on data containing any number and kind of ambiguities. WG4 (Ethical issues) proposes to adopt thorough general principles for any HLA population study to ensure that it conforms to (inter)national legislation or recommendations/guidelines. All HLA-NET guidelines and tools are available through its website http://hla-net.eu.