Jose Manuel Nunes

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.

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.

Haplotype diversity generated by ancient recombination-like events in the MHC of Indian rhesus macaques

The Mamu-A, Mamu-B, and Mamu-DRB genes of the rhesus macaque show several levels of complexity such as allelic heterogeneity (polymorphism), copy number variation, differential segregation of genes/alleles present on a haplotype (diversity) and transcription level differences. A combination of techniques was implemented to screen a large panel of pedigreed Indian rhesus macaques (1,384 individuals representing the offspring of 137 founding animals) for haplotype diversity in an efficient and inexpensive manner. This approach allowed the definition of 140 haplotypes that display a relatively low degree of region variation as reflected by the presence of only 17 A, 18 B and 22 DRB types, respectively, exhibiting a global linkage disequilibrium comparable to that in humans. This finding contrasts with the situation observed in rhesus macaques from other geographic origins and in cynomolgus monkeys from Indonesia. In these latter populations, nearly every haplotype appears to be characterised by a unique A, B and DRB region. In the Indian population, however, a reshuffling of existing segments generated “new” haplotypes. Since the recombination frequency within the core MHC of the Indian rhesus macaques is relatively low, the various haplotypes were most probably produced by recombination events that accumulated over a long evolutionary time span. This idea is in accord with the notion that Indian rhesus macaques experienced a severe reduction in population during the Pleistocene due to a bottleneck caused by geographic changes. Thus, recombination-like processes appear to be a way to expand a diminished genetic repertoire in an isolated and relatively small founder population.

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 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.

Allele Frequency Estimation from Ambiguous Data: Using Resampling Schema in Validating Frequency Estimates and in Selective Neutrality Testing

Abstract The development of molecular typing techniques applied to the study of population genetic diversity originates data with increasing precision but at the cost of some ambiguities. As distinct techniques may produce distinct kinds of ambiguities, a crucial issue is to assess the differences between frequency distributions estimated from data produced by alternative techniques for the same sample. To that aim, we developed a resampling scheme that allows evaluating, by statistical means, the significance of the difference between two frequency distributions. The same approach is then shown to be applicable to test selective neutrality when only sample frequencies are known. The use of these original methods is presented here through an application to the genetic study of a Munda human population sample, where three different HLA loci were typed using two different molecular methods (reverse PCR-SSO typing on microbeads arrays based on Luminex technology and PCR-SSP typing), as described in details in the companion article by Riccio et al. [The Austroasiatic Munda population from India and its enigmatic origin: An HLA diversity study. Hum. Biol. 38:405–435 (2011)]. The differences between the frequency estimates of the two typing techniques were found to be smaller than those resulting from sampling. Overall, we show that using a resampling scheme in validating frequency estimates is effective when alternative frequency estimates are available. Moreover, resampling appears to be the unique way to test selective neutrality when only frequency data are available to describe the genetic structure of populations.

Common and well-documented HLA alleles over all of Europe and within European sub-regions: A catalogue from the European Federation for Immunogenetics

A catalogue of common and well‐documented (CWD) human leukocyte antigen (HLA), previously established by the American Society for Histocompatibility and Immunogenetics (ASHI), is widely used as indicator for typing ambiguities to be resolved in tissue transplantation or for checking the universality of any HLA allele in the world. However, European population samples, which are characterized by a substantial level of genetic variation, are underrepresented in the ASHI catalogue. Therefore, the Population Genetics Working Group of the European Federation for Immunogenetics (EFI) has facilitated data collection for an European CWD catalogue.

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 (GST) 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 GST 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.

Using UNIFORMAT and GENE[RATE] to Analyze Data with Ambiguities in Population Genetics:

Some genetic systems frequently present ambiguous data that cannot be straightforwardly analyzed with common methods of population genetics. Two possibilities arise to analyze such data: one is the arbitrary simplification of the data and the other is the development of methods adapted to such ambiguous data. In this article, we present an attempt at such a development, the UNIFORMAT grammar and the GENEE[RATE] tools, highlighting the specific aspects and the adaptations required to analyze ambiguous nominal data in population genetics.