Pathogen-driven selection and worldwide HLA class I diversity

Population Genetics and Natural Selection in Rheumatic Disease

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Natural selection is an important influence on human genetic variation. Because immune and inflammatory function genes are enriched for signals of positive selection, the prevalence of rheumatic disease-risk alleles seen in different populations is partially the result of differing selective pressures (eg, due to pathogens). This review summarizes the genetic regions associated with susceptibility to different rheumatic diseases and concomitant evidence for natural selection, including known agents of selection exerting selective pressure in these regions.

Kourami: graph-guided assembly for novel human leukocyte antigen allele discovery

Accurate typing of human leukocyte antigen (HLA) is important because HLA genes play important roles in immune responses and disease genesis. Previously available computational methods are database-matching approaches and their outputs are inherently limited by the completeness of already known types, making them unsuitable for discovery of novel alleles. We have developed a graph-guided assembly technique for classical HLA genes, which can construct allele sequences given high-coverage whole-genome sequencing data. Our method delivers highly accurate HLA typing, comparable to the current state-of-the-art methods. Using various data, we also demonstrate that our method can type novel alleles.

Genetics of Natural Killer Cells in Human Health, Disease, and Survival

Natural killer (NK) cells have vital functions in human immunity and reproduction. In the innate and adaptive immune responses to infection, particularly by viruses, NK cells respond by secreting inflammatory cytokines and killing infected cells. In reproduction, NK cells are critical for genesis of the placenta, the organ that controls the supply of oxygen and nutrients to the growing fetus. Controlling NK cell functions are interactions of HLA class I with inhibitory NK cell receptors. First evolved was the conserved interaction of HLA-E with CD94:NKG2A; later established were diverse interactions of HLA-A, -B, and -C with killer cell immunoglobulin-like receptors. Characterizing the latter interactions is rapid evolution, which distinguishes human populations and all species of higher primate. Driving this evolution are the different and competing selections imposed by pathogens on NK cell-mediated immunity and by the constraints of human reproduction on NK cell-mediated placentation. Promoting rapid evolution is independent segregation of polymorphic receptors and ligands throughout human populations.

AmpliSAS and AmpliHLA: Web Server Tools for MHC Typing of Non-Model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are web server tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de-novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human-specific version, it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Here we describe four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

A genomic perspective on HLA evolution

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

Q&A: What are pathogens, and what have they done to and for us?

Microbes are found on us, within us and around us. They inhabit virtually every environment on the planet and the bacteria carried by an average human, mostly in their gut, outnumber human cells. The vast majority of microbes are harmless to us, and many play essential roles in plant, animal and human health. Others, however, are either obligate or facultative pathogens exerting a spectrum of deleterious effects on their hosts. Infectious diseases have historically represented the most common cause of death in humans until recently, exceeding by far the toll taken by wars or famines. From the dawn of humanity and throughout history, infectious diseases have shaped human evolution, demography, migrations and history.

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.

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite-mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection ("trans-species polymorphism"), or alternatively favouring the independent emergence of functionally similar alleles post-speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut-parasite communities in four sympatric lemurs. We tested for parasite-MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co-ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC-parasite co-evolution should be envisaged at the community level. We further show that balancing selection maintaining trans-species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.

Detecting signatures of past pathogen selection on human HLA loci: are there needles in the haystack?

Human leucocyte antigens (HLAs) are responsible for the display of peptide fragments for recognition by T-cell receptors. The gene family encoding them is thus integral to human adaptive immunity, and likely to be under strong pathogen selection. Despite this, it has proved difficult to demonstrate specific examples of pathogen-HLA coevolution. Selection from multiple pathogens simultaneously could explain why the evolutionary signatures of particular pathogens on HLAs have proved elusive. Here, we present an individual-based model of HLA evolution in the presence of two mortality-causing pathogens. We demonstrate that it is likely that individual pathogen species causing high mortality have left recognizable signatures on the HLA genomic region, despite more than one pathogen being present. Such signatures are likely to exist at the whole-population level, and involve haplotypic combinations of HLA genes rather than single loci.

A quantitative and qualitative comparison of illumina MiSeq and 454 amplicon sequencing for genotyping the highly polymorphic major histocompatibility complex (MHC) in a non-model species

Background

High-throughput sequencing enables high-resolution genotyping of extremely duplicated genes. 454 amplicon sequencing (454) has become the standard technique for genotyping the major histocompatibility complex (MHC) genes in non-model organisms. However, illumina MiSeq amplicon sequencing (MiSeq), which offers a much higher read depth, is now superseding 454. The aim of this study was to quantitatively and qualitatively evaluate the performance of MiSeq in relation to 454 for genotyping MHC class I alleles using a house sparrow (Passer domesticus) dataset with pedigree information. House sparrows provide a good study system for this comparison as their MHC class I genes have been studied previously and, consequently, we had prior expectations concerning the number of alleles per individual.

Results

We found that 454 and MiSeq performed equally well in genotyping amplicons with low diversity, i.e. amplicons from individuals that had fewer than 6 alleles. Although there was a higher rate of failure in the 454 dataset in resolving amplicons with higher diversity (6–9 alleles), the same genotypes were identified by both 454 and MiSeq in 98% of cases.

Conclusions

We conclude that low diversity amplicons are equally well genotyped using either 454 or MiSeq, but the higher coverage afforded by MiSeq can lead to this approach outperforming 454 in amplicons with higher diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-017-2654-1) contains supplementary material, which is available to authorized users.

Fast and accurate HLA typing from short-read next-generation sequence data with xHLA

The HLA gene complex on human chromosome 6 is one of the most polymorphic regions in the human genome and contributes in large part to the diversity of the immune system. Accurate typing of HLA genes with short-read sequencing data has historically been difficult due to the sequence similarity between the polymorphic alleles. Here, we introduce an algorithm, xHLA, that iteratively refines the mapping results at the amino acid level to achieve 99-100% four-digit typing accuracy for both class I and II HLA genes, taking only [Formula: see text]3 min to process a 30× whole-genome BAM file on a desktop computer.

Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions

Background

Recent work suggests that gene duplications may play an important role in the evolution of immunity genes. Passerine birds, and in particular Sylvioidea warblers, have highly duplicated major histocompatibility complex (MHC) genes, which are key in immunity, compared to other vertebrates. However, reasons for this high MHC gene copy number are yet unclear. High-throughput sequencing (HTS) allows MHC genotyping even in individuals with extremely duplicated genes. This HTS data can reveal evidence of selection, which may help to unravel the putative functions of different gene copies, i.e. neofunctionalization. We performed exhaustive genotyping of MHC class I in a Sylvioidea warbler, the sedge warbler, Acrocephalus schoenobaenus, using the Illumina MiSeq technique on individuals from a wild study population.

Results

The MHC diversity in 863 genotyped individuals by far exceeds that of any other bird species described to date. A single individual could carry up to 65 different alleles, a large proportion of which are expressed (transcribed). The MHC alleles were of three different lengths differing in evidence of selection, diversity and divergence within our study population. Alleles without any deletions and alleles containing a 6 bp deletion showed characteristics of classical MHC genes, with evidence of multiple sites subject to positive selection and high sequence divergence. In contrast, alleles containing a 3 bp deletion had no sites subject to positive selection and had low divergence.

Conclusions

Our results suggest that sedge warbler MHC alleles that either have no deletion, or contain a 6 bp deletion, encode classical antigen presenting MHC molecules. In contrast, MHC alleles containing a 3 bp deletion may encode molecules with a different function. This study demonstrates that highly duplicated MHC genes can be characterised with HTS and that selection patterns can be useful for revealing neofunctionalization. Importantly, our results highlight the need to consider the putative function of different MHC genes in future studies of MHC in relation to disease resistance and fitness.

Electronic supplementary material

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

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

Coloniality and migration are related to selection on MHC genes in birds

The major histocompatibility complex (MHC) plays a key role in pathogen recognition as a part of the vertebrate adaptive immune system. The great diversity of MHC genes in natural populations is maintained by different forms of balancing selection and its strength should correlate with the diversity of pathogens to which a population is exposed and the rate of exposure. Despite this prediction, little is known about how life-history characteristics affect selection at the MHC. Here, we examined whether the strength of balancing selection on MHC class II genes in birds (as measured with nonsynonymous nucleotide substitutions, dN) was related to their social or migratory behavior, two life-history characteristics correlated with pathogen exposure. Our comparative analysis indicated that the rate of nonsynonymous substitutions was higher in colonial and migratory species than solitary and resident species, suggesting that the strength of balancing selection increases with coloniality and migratory status. These patterns could be attributed to: (1) elevated transmission rates of pathogens in species that breed in dense aggregations, or (2) exposure to a more diverse fauna of pathogens and parasites in migratory species. Our study suggests that differences in social structure and basic ecological traits influence MHC diversity in natural vertebrate populations.

Primary immunodeficiencies suggest redundancy within the human immune system

Pathogen-driven evolution has shaped the complexity of the human immune system. Our genome contains at least 1854 gene products involved in immune responses. However, the redundancy and robustness of the immune system need further characterization. One way to examine this redundancy is through the study of monogenic primary immunodeficiencies (PIDs) associated with infections. Causal mutations affecting innate immunity genes are, in relative terms, close to seven times less frequent than those affecting adaptive immunity genes in PIDs. Loss-of-function mutations of innate immunity genes encoding pattern-recognition receptors (PRRs) and associated pathways rarely cause susceptibility to infections, which suggests that PRR pathways are partially redundant in the immune responses to infection. This dispensability has also been observed for constitutive products of the immune system, such as secretory immunoglobulin A, and for innate immune cells, such as natural killer and innate lymphoid cell subsets, which are not essential for viability. This Review discusses these findings in the context of their implications for the identification of previously unknown classes of PIDs and assessment of the susceptibility to infection associated with various targeted immunotherapies.

Understanding rare and common diseases in the context of human evolution

The wealth of available genetic information is allowing the reconstruction of human demographic and adaptive history. Demography and purifying selection affect the purge of rare, deleterious mutations from the human population, whereas positive and balancing selection can increase the frequency of advantageous variants, improving survival and reproduction in specific environmental conditions. In this review, I discuss how theoretical and empirical population genetics studies, using both modern and ancient DNA data, are a powerful tool for obtaining new insight into the genetic basis of severe disorders and complex disease phenotypes, rare and common, focusing particularly on infectious disease risk.

Genetic Contribution of MHC Class II Genes in Susceptibility to West Nile Virus Infection

WNV is a zoonotic neurotropic flavivirus that has recently emerged globally as a significant cause of viral encephalitis. The last five years, 624 incidents of WNV infection have been reported in Greece. The risk for severe WNV disease increases among immunosuppressed individuals implying thus the contribution of the MHC locus to the control of WNV infection. In order to investigate a possible association of MHC class II genes, especially HLA-DPA1, HLA-DQA1, HLA-DRB1, we examined 105 WNV patients, including 68 cases with neuroinvasive disease and 37 cases with mild clinical phenotype, collected during the period from 2010 to2013, and 100 control individuals selected form the Greek population. Typing was performed for exon 2 for all three genes. DQA1*01:01 was considered to be "protective" against WNV infection (25.4% vs 40.1%, P = 0.004) while DQA1*01:02 was associated with increased susceptibility (48.0% vs 32.1%, P = 0.003). Protection against neuroinvasion was associated with the presence of DRB1*11:02 (4.99% vs 0.0%, P = 0.018). DRB1*16:02 was also absent from the control cohort (P = 0.016). Three additional population control groups were used in order to validate our results. No statistically significant association with the disease was found for HLA-DPA alleles. The results of the present study provide some evidence that MHC class II is involved in the response to WNV infection, outlining infection "susceptibility" and "CNS-high-risk" candidates. Furthermore, three new alleles were identified while the frequency of all alleles in the study was compared with worldwide data. The characterization of the MHC locus could help to estimate the risk for severe WNV cases in a country.

Spatial pattern of genetic diversity and selection in the MHC class II DRB of three Neotropical bat species

Background

Although bats are natural reservoirs of many pathogens, few studies have been conducted on the genetic variation and detection of selection in major histocompatibility complex (MHC) genes. These genes are critical for resistance and susceptibility to diseases, and host–pathogen interactions are major determinants of their extensive polymorphism. Here we examined spatial patterns of diversity of the expressed MHC class II DRB gene of three sympatric Neotropical bats, Carollia perspicillata and Desmodus rotundus (Phyllostomidae), and Molossus molossus (Molossidae), all of which use the same environments (e.g., forests, edge habitats, urban areas). Comparison with neutral marker (mtDNA D-loop) diversity was performed at the same time.

Results

Twenty-three DRB alleles were identified in 19 C. perspicillata, 30 alleles in 35 D. rotundus and 20 alleles in 28 M. molossus. The occurrence of multiple DRB loci was found for the two Phyllostomidae species. The DRB polymorphism was high in all sampling sites and different signatures of positive selection were detected depending on the environment. The patterns of DRB diversity were similar to those of neutral markers for C. perspicillata and M. molossus. In contrast, these patterns were different for D. rotundus for which a geographical structure was highlighted. A heterozygote advantage was also identified for this species. No recombination or gene conversion event was found and phylogenetic relationships showed a trans-species mode of evolution in the Phyllostomids.

Conclusions

This study of MHC diversity demonstrated the strength of the environment and contrasting pathogen pressures in shaping DRB diversity. Differences between positively selected sites identified in bat species highlighted the potential role of gut microbiota in shaping immune responses. Furthermore, multiple geographic origins and/or population admixtures observed in C. perspicillata and M. molossus populations acted as an additional force in shaping DRB diversity. In contrast, DRB diversity of D. rotundus was shaped by environment rather than demographic history.

Electronic supplementary material

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

MYD88 and functionally related genes are associated with multiple infections in a model population of Kenyan village dogs

The purpose of this study was to seek associations between immunity-related molecular markers and endemic infections in a model population of African village dogs from Northern Kenya with no veterinary care and no selective breeding. A population of village dogs from Northern Kenya composed of three sub-populations from three different areas (84, 50 and 55 dogs) was studied. Canine distemper virus (CDV), Hepatozoon canis, Microfilariae (Acantocheilonema dracunculoides, Acantocheilonema reconditum) and Neospora caninum were the pathogens studied. The presence of antibodies (CDV, Neospora), light microscopy (Hepatozoon) and diagnostic PCR (Microfilariae) were the methods used for diagnosing infection. Genes involved in innate immune mechanisms, NOS3, IL6, TLR1, TLR2, TLR4, TLR7, TLR9, LY96, MYD88, and three major histocompatibility genes class II genes were selected as candidates. Single nucleotide polymorphism (SNP) markers were detected by Sanger sequencing, next generation sequencing and PCR-RFLP. The Fisher´s exact test for additive and non-additive models was used for association analyses. Three SNPs within the MYD88 gene and one TLR4 SNP marker were associated with more than one infection. Combined genotypes and further markers identified by next generation sequencing confirmed associations observed for individual genes. The genes associated with infection and their combinations in specific genotypes match well our knowledge on their biological role and on the role of the relevant biological pathways, respectively. Associations with multiple infections observed between the MYD88 and TLR4 genes suggest their involvement in the mechanisms of anti-infectious defenses in dogs.

Modulatory effects of perforin gene dosage on pathogen-associated blood-brain barrier (BBB) disruption

Background

CD8 T cell-mediated blood-brain barrier (BBB) disruption is dependent on the effector molecule perforin. Human perforin has extensive single nucleotide variants (SNVs), the significance of which is not fully understood. These SNVs can result in reduced, but not ablated, perforin activity or expression. However, complete loss of perforin expression or activity results in the lethal disease familial hemophagocytic lymphohistiocytosis type 2 (FHL 2). In this study, we address the hypothesis that a single perforin allele can alter the severity of BBB disruption in vivo using a well-established model of CNS vascular permeability in C57Bl/6 mice. The results of this study provide insight into the significance of perforin SNVs in the human population.

Methods

We isolated the effect a single perforin allele has on CNS vascular permeability through the use of perforin-heterozygous (perforin+/−) C57BL/6 mice in the peptide-induced fatal syndrome (PIFS) model of immune-mediated BBB disruption. Seven days following Theiler’s murine encephalomyelitis virus (TMEV) CNS infection, neuroinflammation and TMEV viral control were assessed through flow cytometric analysis and quantitative real-time PCR of the viral genome, respectively. Following immune-mediated BBB disruption, gadolinium-enhanced T1-weighted MRI, with 3D volumetric analysis, and confocal microscopy were used to define CNS vascular permeability. Finally, the open field behavior test was used to assess locomotor activity of mice following immune-mediated BBB disruption.

Results

Perforin-null mice had negligible CNS vascular permeability. Perforin-WT mice have extensive CNS vascular permeability. Interestingly, perforin-heterozygous mice had an intermediate level of CNS vascular permeability as measured by both gadolinium-enhanced T1-weighted MRI and fibrinogen leakage in the brain parenchyma. Differences in BBB disruption were not a result of increased CNS immune infiltrate. Additionally, TMEV was controlled in a perforin dose-dependent manner. Furthermore, a single perforin allele is sufficient to induce locomotor deficit during immune-mediated BBB disruption.

Conclusions

Perforin modulates BBB disruption in a dose-dependent manner. This study demonstrates a potentially advantageous role for decreased perforin expression in reducing BBB disruption. This study also provides insight into the effect SNVs in a single perforin allele could have on functional deficit in neurological disease.

Alice in microbes' land: adaptations and counter-adaptations of vector-borne parasitic protozoa and their hosts

In the present review, we aim to provide a general introduction to different facets of the arms race between pathogens and their hosts/environment, emphasizing its evolutionary aspects. We focus on vector-borne parasitic protozoa, which have to adapt to both invertebrate and vertebrate hosts. Using Leishmania, Trypanosoma and Plasmodium as main models, we review successively (i) the adaptations and counter-adaptations of parasites and their invertebrate host, (ii) the adaptations and counter-adaptations of parasites and their vertebrate host and (iii) the impact of human interventions (chemotherapy, vaccination, vector control and environmental changes) on these adaptations. We conclude by discussing the practical impact this knowledge can have on translational research and public health.

Excess of Deleterious Mutations around HLA Genes Reveals Evolutionary Cost of Balancing Selection

Deleterious mutations are expected to evolve under negative selection and are usually purged from the population. However, deleterious alleles segregate in the human population and some disease-associated variants are maintained at considerable frequencies. Here, we test the hypothesis that balancing selection may counteract purifying selection in neighboring regions and thus maintain deleterious variants at higher frequency than expected from their detrimental fitness effect. We first show in realistic simulations that balancing selection reduces the density of polymorphic sites surrounding a locus under balancing selection, but at the same time markedly increases the population frequency of the remaining variants, including even substantially deleterious alleles. To test the predictions of our simulations empirically, we then use whole-exome sequencing data from 6,500 human individuals and focus on the most established example for balancing selection in the human genome, the major histocompatibility complex (MHC). Our analysis shows an elevated frequency of putatively deleterious coding variants in nonhuman leukocyte antigen (non-HLA) genes localized in the MHC region. The mean frequency of these variants declined with physical distance from the classical HLA genes, indicating dependency on genetic linkage. These results reveal an indirect cost of the genetic diversity maintained by balancing selection, which has hitherto been perceived as mostly advantageous, and have implications both for the evolution of recombination and also for the epidemiology of various MHC-associated diseases.

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.

Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease

The human major histocompatibility complex (MHC) region has been shown to be associated with numerous diseases. However, it remains a challenge to pinpoint the causal variants for these associations because of the extreme complexity of the region. We thus sequenced the entire 5-Mb MHC region in 20,635 individuals of Han Chinese ancestry (10,689 controls and 9,946 patients with psoriasis) and constructed a Han-MHC database that includes both variants and HLA gene typing results of high accuracy. We further identified multiple independent new susceptibility loci in HLA-C, HLA-B, HLA-DPB1 and BTNL2 and an intergenic variant, rs118179173, associated with psoriasis and confirmed the well-established risk allele HLA-C*06:02. We anticipate that our Han-MHC reference panel built by deep sequencing of a large number of samples will serve as a useful tool for investigating the role of the MHC region in a variety of diseases and thus advance understanding of the pathogenesis of these disorders.

Co-evolution of MHC class I and variable NK cell receptors in placental mammals

Shaping natural killer (NK) cell functions in human immunity and reproduction are diverse killer cell immunoglobulin-like receptors (KIRs) that recognize polymorphic MHC class I determinants. A survey of placental mammals suggests that KIRs serve as variable NK cell receptors only in certain primates and artiodactyls. Divergence of the functional and variable KIRs in primates and artiodactyls predates placental reproduction. Among artiodactyls, cattle but not pigs have diverse KIRs. Catarrhine (humans, apes, and Old World monkeys) and platyrrhine (New World monkeys) primates, but not prosimians, have diverse KIRs. Platyrrhine and catarrhine systems of KIR and MHC class I are highly diverged, but within the catarrhines, a stepwise co-evolution of MHC class I and KIR is discerned. In Old World monkeys, diversification focuses on MHC-A and MHC-B and their cognate lineage II KIR. With evolution of C1-bearing MHC-C from MHC-B, as informed by orangutan, the focus changes to MHC-C and its cognate lineage III KIR. Evolution of C2 from C1 and fixation of MHC-C drove further elaboration of MHC-C-specific KIR, as exemplified by chimpanzee. In humans, the evolutionary trajectory changes again. Emerging from reorganization of the KIR locus and selective attenuation of KIR avidity for MHC class I are the functionally distinctive KIR A and KIR B haplotypes.

Neanderthal genomics suggests a pleistocene time frame for the first epidemiologic transition

High quality Altai Neanderthal and Denisovan genomes are revealing which regions of archaic hominin DNA have persisted in the modern human genome. A number of these regions are associated with response to infection and immunity, with a suggestion that derived Neanderthal alleles found in modern Europeans and East Asians may be associated with autoimmunity. As such Neanderthal genomes are an independent line of evidence of which infectious diseases Neanderthals were genetically adapted to. Sympathetically, human genome adaptive introgression is an independent line of evidence of which infectious diseases were important for AMH coming in to Eurasia and interacting with Neanderthals. The Neanderthals and Denisovans present interesting cases of hominin hunter-gatherers adapted to a Eurasian rather than African infectious disease package. Independent sources of DNA-based evidence allow a re-evaluation of the first epidemiologic transition and how infectious disease affected Pleistocene hominins. By combining skeletal, archaeological and genetic evidence from modern humans and extinct Eurasian hominins, we question whether the first epidemiologic transition in Eurasia featured a new package of infectious diseases or a change in the impact of existing pathogens. Coupled with pathogen genomics, this approach supports the view that many infectious diseases are pre-Neolithic, and the list continues to expand. The transfer of pathogens between hominin populations, including the expansion of pathogens from Africa, may also have played a role in the extinction of the Neanderthals and offers an important mechanism to understand hominin-hominin interactions well back beyond the current limits for aDNA extraction from fossils alone. Am J Phys Anthropol 160:379-388, 2016. © 2016 Wiley Periodicals, Inc.

The Uromodulin Gene Locus Shows Evidence of Pathogen Adaptation through Human Evolution

Common variants in the UMOD gene encoding uromodulin, associated with risk of hypertension and CKD in the general population, increase UMOD expression and urinary excretion of uromodulin, causing salt-sensitive hypertension and renal lesions. To determine the effect of selective pressure on variant frequency, we investigated the allelic frequency of the lead UMOD variant rs4293393 in 156 human populations, in eight ancient human genomes, and in primate genomes. The T allele of rs4293393, associated with CKD risk, has high frequency in most modern populations and was the one detected in primate genomes. In contrast, we identified only the derived, C allele in Denisovan and Neanderthal genomes. The distribution of the UMOD ancestral allele did not follow the ancestral susceptibility model observed for variants associated with salt-sensitive hypertension. Instead, the global frequencies of the UMOD alleles significantly correlated with pathogen diversity (bacteria, helminths) and prevalence of antibiotic-resistant urinary tract infections (UTIs). The inverse correlation found between urinary levels of uromodulin and markers of UTIs in the general population substantiates the link between UMOD variants and protection against UTIs. These data strongly suggest that the UMOD ancestral allele, driving higher urinary excretion of uromodulin, has been kept at a high frequency because of its protective effect against UTIs.

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.

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 gene[rate] tools, highlighting the specific aspects and the adaptations required to analyze ambiguous nominal data in population genetics.

A computational method for identification of vaccine targets from protein regions of conserved human leukocyte antigen binding

Background

Computational methods for T cell-based vaccine target discovery focus on selection of highly conserved peptides identified across pathogen variants, followed by prediction of their binding of human leukocyte antigen molecules. However, experimental studies have shown that T cells often target diverse regions in highly variable viral pathogens and this diversity may need to be addressed through redefinition of suitable peptide targets.

Methods

We have developed a method for antigen assessment and target selection for polyvalent vaccines, with which we identified immune epitopes from variable regions, where all variants bind HLA. These regions, although variable, can thus be considered stable in terms of HLA binding and represent valuable vaccine targets.

Results

We applied this method to predict CD8⁺ T-cell targets in influenza A H7N9 hemagglutinin and significantly increased the number of potential vaccine targets compared to the number of targets discovered using the traditional approach where low-frequency peptides are excluded.

Conclusions

We developed a webserver with an intuitive visualization scheme for summarizing the T cell-based antigenic potential of any given protein or proteome using human leukocyte antigen binding predictions and made a web-accessible software implementation freely available at http://met-hilab.cbs.dtu.dk/blockcons/.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology

OBJECTIVES

The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations.

MATERIALS AND METHODS

We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database.

RESULTS

Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure.

DISCUSSION

Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

Heterogeneity of dN/dS Ratios at the Classical HLA Class I Genes over Divergence Time and Across the Allelic Phylogeny

The classical class I HLA loci of humans show an excess of nonsynonymous with respect to synonymous substitutions at codons of the antigen recognition site (ARS), a hallmark of adaptive evolution. Additionally, high polymporphism, linkage disequilibrium, and disease associations suggest that one or more balancing selection regimes have acted upon these genes. However, several questions about these selective regimes remain open. First, it is unclear if stronger evidence for selection on deep timescales is due to changes in the intensity of selection over time or to a lack of power of most methods to detect selection on recent timescales. Another question concerns the functional entities which define the selected phenotype. While most analyses focus on selection acting on individual alleles, it is also plausible that phylogenetically defined groups of alleles ("lineages") are targets of selection. To address these questions, we analyzed how dN/dS (ω) varies with respect to divergence times between alleles and phylogenetic placement (position of branches). We find that ω for ARS codons of class I HLA genes increases with divergence time and is higher for inter-lineage branches. Throughout our analyses, we used non-selected codons to control for possible effects of inflation of ω associated to intra-specific analysis, and showed that our results are not artifactual. Our findings indicate the importance of considering the timescale effect when analysing ω over a wide spectrum of divergences. Finally, our results support the divergent allele advantage model, whereby heterozygotes with more divergent alleles have higher fitness than those carrying similar alleles.

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.

The role of host genetic factors in respiratory tract infectious diseases: systematic review, meta-analyses and field synopsis

Host genetic factors have frequently been implicated in respiratory infectious diseases, often with inconsistent results in replication studies. We identified 386 studies from the total of 24,823 studies identified in a systematic search of four bibliographic databases. We performed meta-analyses of studies on tuberculosis, influenza, respiratory syncytial virus, SARS-Coronavirus and pneumonia. One single-nucleotide polymorphism from IL4 gene was significant for pooled respiratory infections (rs2070874; 1.66 [1.29–2.14]). We also detected an association of TLR2 gene with tuberculosis (rs5743708; 3.19 [2.03–5.02]). Subset analyses identified CCL2 as an additional risk factor for tuberculosis (rs1024611; OR = 0.79 [0.72–0.88]). The IL4-TLR2-CCL2 axis could be a highly interesting target for translation towards clinical use. However, this conclusion is based on low credibility of evidence - almost 95% of all identified studies had strong risk of bias or confounding. Future studies must build upon larger-scale collaborations, but also strictly adhere to the highest evidence-based principles in study design, in order to reduce research waste and provide clinically translatable evidence.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

Evolving models of the immunopathogenesis of T cell-mediated drug allergy: The role of host, pathogens, and drug response

Immune-mediated (IM) adverse drug reactions (ADRs) are an underrecognized source of preventable morbidity, mortality, and cost. Increasingly, genetic variation in the HLA loci is associated with risk of severe reactions, highlighting the importance of T-cell immune responses in the mechanisms of both B cell-mediated and primary T cell-mediated IM-ADRs. In this review we summarize the role of host genetics, microbes, and drugs in IM-ADR development; expand on the existing models of IM-ADR pathogenesis to address multiple unexplained observations; discuss the implications of this work in clinical practice today; and describe future applications for preclinical drug toxicity screening, drug design, and development.

Spatial pattern of adaptive and neutral genetic diversity across different biomes in the lesser anteater (Tamandua tetradactyla)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and activation of the adaptive immune response and are thought to be regulated by natural selection, especially due to pathogen‐driven selective pressure. In this study, we investigated the spatial distribution of MHC class IIDRB exon 2 gene diversity of the lesser anteater (Tamandua tetradactyla) across five Brazilian biomes using next‐generation sequencing and compared the MHC pattern with that of neutral markers (microsatellites). We found a noticeable high level of diversity in DRB (60 amino acid alleles in 65 individuals) and clear signatures of historical positive selection acting on this gene. Higher allelic richness and proportion of private alleles were found in rain forest biomes, especially Amazon forest, a megadiverse biome, possibly harboring greater pathogen richness as well. Neutral markers, however, showed a similar pattern to DRB, demonstrating the strength of demography as an additional force to pathogen‐driven selection in shaping MHC diversity and structure. This is the first characterization and description of diversity of a MHC gene for any member of the magna‐order Xenarthra, one of the basal lineages of placental mammals.

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.

Among-lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic- or benthic-like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild-caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild-caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.

Distribution of HLA haplotypes across Japanese Archipelago: similarity, difference and admixture

The human leukocyte antigen (HLA) region is the most polymorphic region in the human genome. The polymorphic nature of the HLA region is thought to have been shaped from balancing selection. The complex migration events during the Out-of-Africa expansion have influenced geographic patterns of HLA allele frequencies and diversities across present-day human populations. Differences in the HLA allele frequency may contribute geographic differences in the susceptibility to many diseases, such as infectious, autoimmune and metabolic diseases. Here we briefly reviewed characteristics of frequency distribution of HLA alleles and haplotypes in Japanese population. A large part of HLA alleles and haplotypes that are common in Japanese are shared with neighboring Asian populations. The differentiations in HLA alleles and haplotypes across Japanese regional populations may provide clues to model for peopling of Japanese Archipelago and for design of genetic association studies. Finally, we introduce recent topics that new HLA alleles derived from ancient admixtures with Neanderthals and Denisovans are thought to have played an important role in the adaptation of modern humans to local pathogens during Out-of-Africa expansion.