Stéphane Buhler

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

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajimas tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies. Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.

Immunogenetics as a tool in anthropological studies

The genes coding for the main molecules involved in the human immune system – immunoglobulins, human leucocyte antigen (HLA) molecules and killer‐cell immunoglobulin‐like receptors (KIR) – exhibit a very high level of polymorphism that reveals remarkable frequency variation in human populations. ‘Genetic marker’ (GM) allotypes located in the constant domains of IgG antibodies have been studied for over 40 years through serological typing, leading to the identification of a variety of GM haplotypes whose frequencies vary sharply from one geographic region to another. An impressive diversity of HLA alleles, which results in amino acid substitutions located in the antigen‐binding region of HLA molecules, also varies greatly among populations. The KIR differ between individuals according to both gene content and allelic variation, and also display considerable population diversity. Whereas the molecular evolution of these polymorphisms has most likely been subject to natural selection, principally driven by host–pathogen interactions, their patterns of genetic variation worldwide show significant signals of human geographic expansion, demographic history and cultural diversification. As current developments in population genetic analysis and computer simulation improve our ability to discriminate among different – either stochastic or deterministic – forces acting on the genetic evolution of human populations, the study of these systems shows great promise for investigating both the peopling history of modern humans in the time since their common origin and human adaptation to past environmental (e.g. pathogenic) changes. Therefore, in addition to mitochondrial DNA, Y‐chromosome, microsatellites, single nucleotide polymorphisms and other markers, immunogenetic polymorphisms represent essential and complementary tools for anthropological studies.

Analysis of the HLA population data (AHPD) submitted to the 15th International Histocompatibility/Immunogenetics Workshop by using the Gene[rate] computer tools accommodating ambiguous data (AHPD project report)

During the 15th International Histocompatibility and Immunogenetics Workshop (IHIWS), 14 human leukocyte antigen (HLA) laboratories participated in the Analysis of HLA Population Data (AHPD) project where 18 new population samples were analyzed statistically and compared with data available from previous workshops. To that aim, an original methodology was developed and used (i) to estimate frequencies by taking into account ambiguous genotypic data, (ii) to test for Hardy-Weinberg equilibrium (HWE) by using a nested likelihood ratio test involving a parameter accounting for HWE deviations, (iii) to test for selective neutrality by using a resampling algorithm, and (iv) to provide explicit graphical representations including allele frequencies and basic statistics for each series of data. A total of 66 data series (1-7 loci per population) were analyzed with this standard approach. Frequency estimates were compliant with HWE in all but one population of mixed stem cell donors. Neutrality testing confirmed the observation of heterozygote excess at all HLA loci, although a significant deviation was established in only a few cases. Population comparisons showed that HLA genetic patterns were mostly shaped by geographic and/or linguistic differentiations in Africa and Europe, but not in America where both genetic drift in isolated populations and gene flow in admixed populations led to a more complex genetic structure. Overall, a fruitful collaboration between HLA typing laboratories and population geneticists allowed finding useful solutions to the problem of estimating gene frequencies and testing basic population diversity statistics on highly complex HLA data (high numbers of alleles and ambiguities), with promising applications in either anthropological, epidemiological, or transplantation studies.

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.

HLA class II genetic diversity in southern Tunisia and the Mediterranean area

North Africa is populated by many Arab‐ and Berber‐speaking populations whose genetic history is still poorly understood. In this study, we analyse the HLA‐DRB1 and DQB1 molecular diversity in three populations from the south of Tunisia — Berbers from Jerba, Berbers from Matmata and Arabs from Gabes — and we compare them to a large set of populations from the whole Mediterranean region. Among the three populations studied, the Berbers from Jerba are the most peculiar, as they diverge significantly from other North Africans while being genetically highly diversified and close to populations from the Near East. Thus, Jerba may have been a crossing point, in historical times, where colonization from the eastern Mediterranean area left significant genetic traces. By contrast, the populations from Matmata and Gabes are genetically similar to other Arab and Berber‐speaking populations from different areas of the Maghrib, despite some peculiar allele and haplotype frequencies. At a larger scale, northwest Africa and southwest Europe are closely related according to these polymorphisms, while the populations from the eastern Mediterranean area are much more differentiated. The close genetic relatedness found for HLA among populations of the western Mediterranean region challenges previous results based on Y chromosome analyses, where the Gibraltar Strait appeared to constitute a main genetic barrier.

Positive association of DRB1*04 and DRB1*15 alleles with Fya immunization in a Southern European population

BACKGROUND: Anti‐Fya has been implicated in hemolytic transfusion reactions. However, not all Fy(a−) patients develop anti‐Fya after transfusion with 1 unit of blood [Fy(a+)]. This study was designed to identify HLA‐DRB1 alleles associated with a predisposition to Fya immunization after blood transfusion.

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.

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.

HLA-G UTR Haplotype Conservation in the Malian Population: Association with Soluble HLA-G

The HLA-G molecule plays an important role in immunomodulation. In a previous study carried out on a southern French population our team showed that HLA-G haplotypes, defined by SNPs in the coding region and specific SNPs located in 5′URR and 3′UTR regulatory regions, are associated with differential soluble HLA-G expression (sHLA-G). Furthermore, the structure of these HLA-G haplotypes appears to be conserved in geographically distant populations. The aim of our study is to confirm these expectations in a sub-Saharan African population and to explore additional factors, such as HLA-A alleles, that might influence sHLA-G expression. DNA and plasma samples were collected from 229 Malians; HLA-G and HLA-A genotyping were respectively performed by the Snap Shot® method and by Luminex™ technology. sHLA-G dosage was performed using an ELISA kit. HLA-G and HLA-A allelic and haplotypic frequencies were estimated using an EM algorithm from the Gene[Rate] program. Associations between genetic and non genetic parameters with sHLA-G were performed using a non-parametric test with GRAPH PAD Prism 5. Our results reveal a good conservation of the HLA-G UTR haplotype structure in populations with different origins and demographic histories. These UTR haplotypes appear to be involved in different sHLA-G expression patterns. Specifically, the UTR-2 haplotype was associated with low sHLA-G levels, displaying a dominant negative effect. Furthermore, an allelic effect of both HLA-G and HLA-A, as well as non genetic parameters, such as age and gender possibly linked to osteogenesis and sexual hormones, also seem to be involved in the modulation of sHLA-G. These data suggest that further investigation in larger cohorts and in populations from various ethnical backgrounds is necessary not only to detect new functional polymorphism in HLA-G regulatory regions, but also to reveal the extent of biological phenomena that influence sHLA-G secretion and this might therefore have an impact on transplantation practice.

A New HLA Map of Europe: Regional Genetic Variation and Its Implication for Peopling History, Disease-Association Studies and Tissue Transplantation

Objectives: HLA genes are highly polymorphic in human populations as a result of diversifying selection related to their immune function. However, HLA geographic variation worldwide suggests that demographic factors also shaped their evolution. We here analyzed in detail HLA genetic variation in Europe in order to identify signatures of migration history and/or natural selection. Methods: Relationships between HLA diversity and geography were analyzed at 7 loci through several approaches including linear regression on gene diversity and haplotype frequencies. Regional variation was also assessed on HLA multi-locus phenotypes through structure analysis. Deviation from neutrality was tested by resampling. Results: Geographic distance was a strong predictor of HLA variation at 5 loci (A, B, C, DRB1 and DPB1) in Europe, and latitude significantly shaped HLA gene diversity and haplotype frequencies. Whereas the main level of genetic diversity was found within populations, both HLA gene frequencies and phenotypic profiles revealed regional variation, Southeast Europe, Great Britain and Finland being the most distinctive. Effects of natural selection were suggested at the DQ loci. Conclusions: HLA regional variation was observed in Europe and can be related to population history, locus HLA-A providing by far the strongest signals. This new HLA map of Europe represents an invaluable reference for disease-association studies and tissue transplantation.