Challenging Ancient DNA Results About Putative HLA Protection or Susceptibility to Yersinia pestis.

In a recent article, Immel et al. (Immel A, Key FM, Szolek A, Barquera R, Robinson MK, Harrison GF, Palmer WH, Spyrou MA, Susat J, Krause-Kyora B, et al. 2021. Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes. Mol Biol Evol. 38:4059-4076) extracted DNA from 36 individuals dead from plague in Ellwangen, Southern Germany, during the 16th century. By comparing their human leukocyte antigen (HLA) genotypes with those of 50 present-day Ellwangen inhabitants, the authors reported a significant decrease of HLA-B*51:01 and HLA-C*06:02 and a significant increase of HLA-DRB1*13:01/13:02 frequencies from ancient to modern populations. After comparing these frequencies with a larger sample of 8,862 modern Germans and performing simulations of natural selection, they concluded that these changes had been driven by natural selection. In an attempt to provide more evidence on such stimulating results, we explored the HLA frequency patterns over all of Europe, we predicted binding affinities of HLA-B/C/DRB1 alleles to 106,515 Yersinia pestis-derived peptides, and we performed forward simulations of HLA genetic profiles under neutrality. Our analyses do not sustain the conclusions of HLA protection or susceptibility to plague based on ancient DNA.

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules.

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution.

BACKGROUND: Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations. RESULTS: Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations. CONCLUSIONS: We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

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

The above article originally published with an incomplete bibliographic information for Bitarello et al. (2016) and presented correctly in this article.

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.

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.

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.

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.

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.

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.

HLA-DRB1 and HLA-DQB1 allele associations in an Albanian patient population with rheumatoid arthritis: correlations with the specific autoantibody markers and inter-population DRB1 allele frequency variability.

The prevalence of rheumatoid arthritis and its specific autoantibodies varies in different populations. This variability depends on the genetic polymorphism of the immune response genes among which the HLA system plays a major role. In this context, we studied the HLA-DRB1 and HLA-DQB1 first-level allele frequencies in 100 Albanian patients with rheumatoid arthritis (RA), and taking into account their rheumatoid factor (RF) and anticitrullinated peptide antibodies (ACPA) serologic subgroups, we compared them with the respective frequencies in a population of 191 Albanian individuals without known pathology. No differences were found between the controls and the RA patient group as a whole, but three statistically significant differences were found: an increase in DRB1*04 among ACPA+, RF+ and ACPA+/RF+ patients, a significant decrease in DRB1*11 among ACPA+/RF+ and also a decrease in DRB1*13 among RF+ patient subgroups. Comparing allele frequencies of putatively associated RA alleles in different European populations revealed a significant negative correlation between the RA predisposing DRB1*04 and protective DRB1*11 allele frequencies. A statistically significant correlation was also found between RA prevalence rates and DRB1*04 as well as DRB1*11 frequencies. The relatively low frequencies of DRB1*04 and high DRB1*11 in the Albanian population might explain the rather low positivity rate of ACPA and RF antibodies among the Albanian RA patients. These specific association patterns suggest that this first study of RA in an Albanian population should be followed up to include second level or higher definition of HLA alleles and to compare RA patterns among European populations.

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.

The most frequent HLA alleles around the world: A fundamental synopsis.

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

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.

Population designations in biomedical research: Limitations and perspectives.

In biomedical research, population differences are of central interest. Variations in the frequency and severity of diseases and in treatment effects among human subpopulation groups are common in many medical conditions. Unfortunately, the practices in terms of subpopulation labeling do not exhibit the level of rigor one would expect in biomedical research, especially when studying multifactorial diseases such as cancer or atherosclerosis. The reporting of population differences in clinical research is characterized by large disparities in practices, and fraught with methodological issues and inconsistencies. The actual designations such as "Black" or "Asian" refer to broad and heterogeneous groups, with a great discrepancy among countries. Moreover, the use of obsolete concepts such as "Caucasian" is unfortunate and imprecise. The use of adequate labeling to reflect the scientific hypothesis needs to be promoted. Furthermore, the use of "race/ethnicity" as a unique cause of human heterogeneity may distract from investigating other factors related to a medical condition, particularly if this label is employed as a proxy for cultural habits, diet, or environmental exposure. In addition, the wide range of opinions among researchers does not facilitate the attempts made for resolving this heterogeneity in labeling. "Race," "ethnicity," "ancestry," "geographical origin," and other similar concepts are saturated with meanings. Even if the feasibility of a global consensus on labeling seems difficult, geneticists, sociologists, anthropologists, and ethicists should help develop policies and practices for the biomedical field.

A registry-based population study of the HLA in Québec, Canada.

As part of the worldwide effort to better characterize HLA diversity in populations, we have studied the population of Québec in Canada. This province has been defined by a complex history with multiple founder effects and migration patterns. We analyzed the typing data of 3806 individuals registered in Héma-Québec's Registry, which covered most administrative regions in Québec. Typing information was resolved at the second field level of resolution by next-generation sequencing (NGS) or by Sanger sequencing. We used the HLA-net.eu GENE[RATE] tools to estimate allele and two-locus haplotype frequencies for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1, as well as Hardy-Weinberg equilibrium (HWE), selective neutrality, and linkage disequilibrium. The chord genetic distance was also calculated between administrative regions and was visualized using non-metric multidimensional scaling (NMDS) analysis. While most individual regions were in HWE, HWE was rejected for the province considered as a whole. Some regions exhibited signatures of selection, mostly toward an excess of heterozygotes. Allele and haplotype frequencies revealed outlier regions that strongly differed from the other regions. NMDS plots also showed differences between regions. The administrative regions of the province of Québec displayed heterogeneity in their HLA profiles. This heterogeneity was attributable to differing allele and haplotype specificities by region. In particular, regions 02-Saguenay-Lac-Saint-Jean and 01-Bas-St-Laurent diverged from the rest of the regions. The urban regions 06-Montréal and 13-Laval were very diversified in their HLA profiles. Together, these results will help optimize donor recruitment strategies in Québec.

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.

Allele frequency estimation from ambiguous data: using resampling schema in validating frequency estimates and in selective neutrality testing.

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.

The Austroasiatic Munda population from India and Its enigmatic origin: a HLA diversity study.

The Austroasiatic linguistic family disputes its origin between two geographically distant regions of Asia, India, and Southeast Asia, respectively. As genetic studies based on classical and gender-specific genetic markers provided contradictory results to this debate thus far, we investigated the HLA diversity (HLA-A, -B, and -DRB1 loci) of an Austroasiatic Munda population from Northeast India and its relationships with other populations from India and Southeast Asia. Because molecular methods currently used to test HLA markers often provide ambiguous results due to the high complexity of this polymorphism, we applied two different techniques (reverse PCR-SSO typing on microbeads arrays based on Luminex technology, and PCR-SSP typing) to type the samples. After validating the resulting frequency distributions through the original statistical method described in our companion article ( Nunes et al. 2011 ), we compared the HLA genetic profile of the sampled Munda to those of other Asiatic populations, among which Dravidian and Indo-European-speakers from India and populations from East and Southeast Asia speaking languages belonging to different linguistic families. We showed that the Munda from Northeast India exhibit a peculiar genetic profile with a reduced level of HLA diversity compared to surrounding Indian populations. They also exhibit less diversity than Southeast Asian populations except at locus DRB1. Several analyses using genetic distances indicate that the Munda are much more closely related to populations from the Indian subcontinent than to Southeast Asian populations speaking languages of the same Austroasiatic linguistic family. On the other hand, they do not share a closer relationship with Dravidians compared with Indo-Europeans, thus arguing against the idea that the Munda share a common and ancient Indian origin with Dravidians. Our results do not favor either a scenario where the Munda would be representative of an ancestral Austroasiatic population giving rise to an eastward Austroasiatic expansion to Southeast Asia. Rather, their peculiar genetic profile is better explained by a decrease in genetic diversity through genetic drift from an ancestral population having a genetic profile similar to present-day Austroasiatic populations from Southeast Asia (thus suggesting a possible southeastern origin), followed by intensive gene flow with neighboring Indian populations. This conclusion is in agreement with archaeological and linguistic information. The history of the Austroasiatic family represents a fascinating example where complex interactions among culturally distinct human populations occurred in the past.

Challenging views on the peopling history of East Asia: the story according to HLA markers.

The peopling of East Asia by the first modern humans is strongly debated from a genetic point of view. A north-south genetic differentiation observed in this geographic area suggests different hypotheses on the origin of Northern East Asian (NEA) and Southern East Asian (SEA) populations. In this study, the highly polymorphic HLA markers were used to investigate East Asian genetic diversity. Our database covers a total of about 127,000 individuals belonging to 84 distinct Asian populations tested for HLA-A, -B, -C, -DPB1, and/or -DRB1 alleles. Many Chinese populations are represented, which have been sampled in the last 30 years but rarely taken into account in international research due to their data published in Chinese. By using different statistical methods, we found a significant correlation between genetics and geography and relevant genetic clines in East Asia. Additionally, HLA alleles appear to be unevenly distributed: some alleles observed in NEA populations are widespread at the global level, while some alleles observed in SEA populations are virtually unique in Asia. The HLA genetic variation in East Asia is also characterized by a decrease of diversity from north to south, although a reverse pattern appears when one only focuses on alleles restricted to Asia. These results reflect a more complex migration history than that illustrated by the "southern-origin" hypothesis, as genetic contribution of ancient human migrations through a northern route has probably been quite substantial. We thus suggest a new overlapping model where northward and southward opposite migrations occurring at different periods overlapped.