Intensity of natural selection at the major histocompatibility complex loci

Graduate Student Literature Review: The DRB3 gene of the bovine major histocompatibility complex-Discovery, diversity, and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine major histocompatibility complex (MHC), also known as the bovine leukocyte antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. The BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and Creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

Early human migration determines the risk of being attacked by wolves: ethnic gene diversity on the development of systemic lupus erythematosus

The prevalence of systemic lupus erythematosus (SLE) varies significantly based on ethnicity rather than geographic distribution; thus, the prevalence is higher in Asian, Hispanic, and Black African populations than in European populations. The risk of developing lupus nephritis (LN) is the highest among Asian populations. Therefore, we hypothesize that human genetic diversity between races has occurred through the early human migration and human genetic adaptation to various environments, with a particular focus on pathogens. Additionally, we compile the currently available evidence on the ethnic gene diversity of SLE and how it relates to disease severity. The human leukocyte antigen (HLA) locus is well established as associated with susceptibility to SLE; specific allele distributions have been observed across diverse populations. Notably, specific amino acid residues within these HLA loci demonstrate significant associations with SLE risk. The non-HLA genetic loci associated with SLE risk also varies across diverse ancestries, implicating distinct immunological pathways, such as the type-I interferon and janus kinase-signal transducers and activators of transcription (JAK-STAT) pathways in Asians, the type-II interferon signaling pathway in Europeans, and B cell activation pathway in Africans. Furthermore, assessing individual genetic susceptibility using genetic risk scores (GRS) for SLE helps to reveal the diverse prevalence, age of onset, and clinical phenotypes across different ethnicities. A higher GRS increases the risk of LN and the severity of SLE. Therefore, understanding ethnic gene diversity is crucial for elucidating disease mechanisms and SLE severity, which could enable the development of novel drugs specific to each race.

Joint estimation of selection intensity and mutation rate under balancing selection with applications to HLA

A composite likelihood method is introduced for jointly estimating the intensity of selection and the rate of mutation, both scaled by the effective population size, when there is balancing selection at a single multi-allelic locus in an isolated population at demographic equilibrium. The performance of the method is tested using simulated data. Average estimated mutation rates and selection intensities are close to the true values but there is considerable variation about the averages. Allowing for both population growth and population subdivision does not result in qualitative differences but the estimated mutation rates and selection intensities do not in general reflect the current effective population size. The method is applied to three class I (HLA-A, HLA-B and HLA-C) and two class II loci (HLA-DRB1 and HLA-DQA1) in the 1000 Genomes populations. Allowing for asymmetric balancing selection has only a slight effect on the results from the symmetric model. Mutations that restore symmetry of the selection model are preferentially retained because of the tendency of natural selection to maximize average fitness. However, slight differences in selective effects result in much longer persistence time of some alleles. Trans-species polymorphism (TSP), which is characteristic of major-histocompatibility loci in vertebrates, is more likely when there are small differences in allelic fitness than when complete symmetry is assumed. Therefore, variation in allelic fitness expands the range of parameter values consistent with observations of TSP.

Challenging ancient DNA results about putative HLA protection or susceptibility to Yersina 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.

Fate of MHCII in salmonids following 4WGD

Major histocompatibility complex (MHC) genes are key players in the adaptive immunity providing a defense against invading pathogens. Although the basic structures are similar when comparing mammalian and teleost MHC class II (MHCII) molecules, there are also clear-cut differences. Based on structural requirements, the teleosts non-classical MHCII molecules do not comply with a function similar to the human HLA-DM and HLA-DO, i.e., assisting in peptide loading and editing of classical MHCII molecules. We have previously studied the evolution of teleost class II genes identifying various lineages and tracing their phylogenetic occurrence back to ancient ray-finned fishes. We found no syntenic MHCII regions shared between cyprinids, salmonids, and neoteleosts, suggesting regional instabilities. Salmonids have experienced a unique whole genome duplication 94 million years ago, providing them with the opportunity to experiment with gene duplicates. Many salmonid genomes have recently become available, and here we set out to investigate how MHCII has evolved in salmonids using Northern pike as a diploid sister phyla, that split from the salmonid lineage prior to the fourth whole genome duplication (4WGD) event. We identified 120 MHCII genes in pike and salmonids, ranging from 11 to 20 genes per species analyzed where DB-group genes had the most expansions. Comparing the MHC of Northern pike with that of Atlantic salmon and other salmonids species provides a tale of gene loss, translocations, and genome rearrangements.

Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar

Background

Myanmar cattle populations predominantly consist of native cattle breeds (Pyer Sein and Shwe), characterized by their geographical location and coat color, and the Holstein-Friesian crossbreed, which is highly adapted to the harsh tropical climates of this region. Here, we analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus that has been linked to the immune response, in Myanmar cattle populations.

Methods

Blood samples (n = 294) were taken from two native breeds (Pyer Sein, n = 163 and Shwe Ni, n = 69) and a cattle crossbreed (Holstein-Friesian, n = 62) distributed across six regions of Myanmar (Bago, n = 38; Sagaing, n = 77; Mandalay, n = 46; Magway, n = 46; Kayin, n = 43; Yangon, n = 44). In addition, a database that included 2428 BoLA-DRB3 genotypes from European (Angus, Hereford, Holstein, Shorthorn, Overo Negro, Overo Colorado, and Jersey), Zebuine (Nellore, Brahman and Gir), Asian Native from Japan and Philippine and Latin-American Creole breeds was also included. Furthermore, the information from the IPD–MHC database was also used in the present analysis. DNA was genotyped using the sequence-based typing method. DNA electropherograms were analyzed using the Assign 400ATF software.

Results

We detected 71 distinct alleles, including three new variants for the BoLA-DRB3 gene. Venn analysis showed that 11 of these alleles were only detected in Myanmar native breeds and 26 were only shared with Asian native and/or Zebu groups. The number of alleles ranged from 33 in Holstein-Friesians to 58 in Pyer Seins, and the observed versus unbiased expected heterozygosity were higher than 0.84 in all the three the populations analyzed. The F_ST analysis showed a low level of genetic differentiation between the two Myanmar native breeds (F_ST = 0.003), and between these native breeds and the Holstein-Friesians (F_ST <  0.021). The average F_ST value for all the Myanmar Holstein-Friesian crossbred and Myanmar native populations was 0.0136 and 0.0121, respectively. Principal component analysis (PCA) and tree analysis showed that Myanmar native populations grouped in a narrow cluster that diverged clearly from the Holstein-Friesian populations. Furthermore, the BoLA-DRB3 allele frequencies suggested that while some Myanmar native populations from Bago, Mandalay and Yangon regions were more closely related to Zebu breeds (Gir and Brahman), populations from Kayin, Magway and Sagaing regions were more related to the Philippines native breeds. On the contrary, PCA showed that the Holstein-Friesian populations demonstrated a high degree of dispersion, which is likely the result of the different degrees of native admixture in these populations.

Conclusion

This study is the first to report the genetic diversity of the BoLA-DRB3 gene in two native breeds and one exotic cattle crossbreed from Myanmar. The results obtained contribute to our understanding of the genetic diversity and distribution of BoLA-DRB3 gene alleles in Myanmar, and increases our knowledge of the worldwide variability of cattle BoLA-DRB3 genes, an important locus for immune response and protection against pathogens.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

Divergent Allele Advantage Provides a Quantitative Model for Maintaining Alleles with a Wide Range of Intrinsic Merits

The Major Histocompatibility Complex (MHC) is the most genetically diverse region of the genome in most vertebrates. Some form of balancing selection is necessary to account for the extreme diversity, but the precise mechanism of balancing selection is unknown. Due to the way MHC molecules determine immune recognition, overdominance (also referred to as heterozygote advantage) has been suggested as the main driving force behind this unrivalled diversity. However, both theoretical results and simulation models have shown that overdominance in its classical form cannot maintain large numbers of alleles unless all alleles confer unrealistically similar levels of fitness. There is increasing evidence that heterozygotes containing genetically divergent alleles allow for broader antigen presentation to immune cells, providing a selective mechanism for MHC polymorphism. By framing competing models of overdominance within a general framework, we show that a model based on Divergent Allele Advantage (DAA) provides a superior mechanism for maintaining alleles with a wide range of intrinsic merits, as intrinsically less-fit MHC alleles that are more divergent can survive under DAA. Specifically, our results demonstrate that a quantitative mechanism built from the DAA hypothesis is able to maintain polymorphism in the MHC. Applying such a model to both livestock breeding and conservation could provide a better way of identifying superior heterozygotes, and quantifying the advantages of genetic diversity at the MHC.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Whereas specialist major histocompatibility complex (MHC) molecules initiate immune response against only relatively few pathogens, generalists provide protection against a broad range. Accordingly, this study shows that the geographical distribution of generalist MHC alleles in human populations reflects exposure to diverse infectious diseases.

Variation in the human genome influences our susceptibility to infectious diseases, but the causal link between disease and underlying mutation often remains enigmatic. Major histocompatibility complex II (MHC class II) molecules shape both our immune response against pathogens and our tolerance of self-peptides. The genomic region that encodes MHC molecules is of particular interest, as it is home to more genetic disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Here, we propose that MHC class II molecules can be categorized into two major types; specialists initiate effective immune response against only relatively few pathogens, while generalists provide protection against a broad range of pathogens. As support, we demonstrate that generalist MHC class II variants are more prevalent in human populations residing in pathogen-rich areas.

Peptide presentation by HLA-DQ molecules is associated with the development of immune tolerance

HLA class II proteins are important elements of human adaptive immune recognition and are associated with numerous infectious and immune-mediated diseases. These highly variable molecules can be classified into DP, DQ and DR groups. It has been proposed that in contrast with DP and DR, epitope binding by DQ variants rather results in immune tolerance. However, the pieces of evidence are limited and controversial. We found that DQ molecules bind more human epitopes than DR. Pathogen-associated epitopes bound by DQ molecules are more similar to human proteins than the ones bound by DR. Accordingly, DQ molecules bind epitopes of significantly different pathogen species. Moreover, the binding of autoimmunity-associated epitopes by DQ confers protection from autoimmune diseases. Our results suggest a special role of HLA-DQ in immune homeostasis and help to better understand the association of HLA molecules with infectious and autoimmune diseases.

Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations

Background

Bovine leukocyte antigens (BoLAs) are used extensively as markers of disease and immunological traits in cattle. However, until now, characterization of BoLA gene polymorphisms in Zebu breeds using high resolution typing methods has been poor. Here, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 421 cattle (116 Bolivian Nellore, 110 Bolivian Gir, and 195 Peruvian Nellore-Brahman). Data from 1416 Taurine and Zebu samples were also included in the analysis.

Results

We identified 46 previously reported alleles and no novel variants. Of note, 1/3 of the alleles were detected only in Zebu cattle. Comparison of the degree of genetic variability at the population and sequence levels with genetic distance in the three above mentioned breeds and nine previously reported breeds revealed that Zebu breeds had a gene diversity score higher than 0.86, a nucleotide diversity score higher than 0.06, and a mean number of pairwise differences greater than 16, being similar to those estimated for other cattle breeds. A neutrality test revealed that only Nellore-Brahman cattle showed the even gene frequency distribution expected under a balanced selection scenario. The F_ST index and the exact G test showed significant differences across all cattle populations (F_ST = 0.057; p <  0.001). Neighbor-joining trees and principal component analysis identified two major clusters: one comprising mainly European Taurine breeds and a second comprising Zebu breeds. This is consistent with the historical and geographical origin of these breeds. Some of these differences may be explained by variation of amino acid motifs at antigen-binding sites.

Conclusions

The results presented herein show that the historical divergence between Taurine and Zebu cattle breeds is a result of origin, selection, and adaptation events, which would explain the observed differences in BoLA-DRB3 gene diversity between the two major bovine types. This allelic information will be important for investigating the relationship between the major histocompatibility complex and disease, and contribute to an ongoing effort to catalog bovine MHC allele frequencies according to breed and location.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0618-7) contains supplementary material, which is available to authorized users.

Selection outweighs drift at a fine scale: Lack of MHC differentiation within a family living lizard across geographically close but disconnected rocky outcrops

The highly polymorphic genes of the major histocompatibility complex (MHC) are involved in disease resistance, mate choice and kin recognition. Therefore, they are widely used markers for investigating adaptive variation. Although selection is the key driver, gene flow and genetic drift also influence adaptive genetic variation, sometimes in opposing ways and with consequences for adaptive potential. To further understand the processes that generate MHC variation, it is helpful to compare variation at the MHC with that at neutral genetic loci. Differences in MHC and neutral genetic variation are useful for inferring the relative influence of selection, gene flow and drift on MHC variation. To date, such investigations have usually been undertaken at a broad spatial scale. Yet, evolutionary and ecological processes can occur at a fine spatial scale, particularly in small or fragmented populations. We investigated spatial patterns of MHC variation among three geographically close, naturally discrete, sampling sites of Egernia stokesii, an Australian lizard. The MHC of E. stokesii has recently been characterized, and there is evidence for historical selection on the MHC. We found E. stokesii MHC weakly differentiated among sites compared to microsatellites, suggesting selection, acting similarly at each site, has outweighed any effects of low gene flow or of genetic drift on E. stokesii MHC variation. Our findings demonstrate the strength of selection in shaping patterns of MHC variation or consistency at a fine spatial scale.

The investigation of the origin of Southern Tunisians using HLA genes

The south of Tunisia is characterized by marked ethnic diversity, highlighted by the coexistence of native Berbers with Blacks, Jews and Arab-speaking populations. Despite this heterogeneity, genetic anthropology studies investigating the origin of current Southern Tunisians were rarely reported. We examined human leukocyte antigen (HLA) class I (A, B) and class II (DRB1, DQB1) gene profiles of 250 unrelated Southern Tunisians, and compared them with those of Arab-speaking communities, along with Mediterranean and sub-Sahara African populations using genetic distances, neighbor-joining dendrograms, correspondence and haplotype analysis. In total, 137 HLA alleles were detected, which comprised 32 HLA-A, 52 HLA-B, 32 DRB1 and 21 DQB1 alleles. The most frequent alleles were HLA-A*02:01(18.02%), HLA-B*50:01 (9.11%), HLA-DRB1*07:01 (22.06%) and HLA-DQB1*02:01 (17.21%). All pairs of HLA loci show significant linkage disequilibrium. The four loci depict negative F_nd (the normalized deviate of the homozygosity) values indicating an overall trend to balancing selection. Southern Tunisians appear to be closely related to others Tunisian populations including Berbers, North Africans and Iberians. On the contrary, Southern Tunisians were distinct from Palestinian, Lebanese and Jordanian Middle Eastern Arab-speaking population, despite the deep Arab incursions and Arabization that affected Southern Tunisia. In addition, Southern Tunisians were distant from many sub-Saharan communities, evidenced by genetic distance analysis. Collectively, this indicates a limited genetic contribution of Arab invasion and Black caravans on the makeup of Southern Tunisian gene pool.

Elucidating the origin of HLA-B*73 allelic lineage: Did modern humans benefit by archaic introgression?

A previous study reported that some of the human leukocyte antigen (HLA) alleles and haplotypes in present-day humans were acquired by admixture with archaic humans; specifically, an exceptionally diverged HLA-B*73 allele was proposed to be transmitted from Denisovans, although the DNA sequence of HLA-B*73 has not been detected in the Denisovan genome. Here, we argue against the hypothesis that HLA-B*73 introgressed from Denisovans into early modern humans. A phylogenetic analysis revealed that HLA-B*73:01 formed a monophyletic group with a chimpanzee MHC-B allele, strongly suggesting that the HLA-B*73 allelic lineage has been maintained in humans as well as in chimpanzees since the divergence of humans and chimpanzees. The global distribution of HLA-B*73 allele showed that the population frequency of HLA-B*73 in west Asia (0.24 %)—a possible site of admixture with Denisovans—is lower than that in Europe (0.72 %) and in south Asia (0.69 %). Furthermore, HLA-B*73 is not observed in Melanesia even though the Melanesian genome contains the highest proportion of Denisovan ancestry in present-day human populations. Single nucleotide polymorphisms in HLA-A*11-HLA-C*12:02 or HLA-A*11-C*15 haplotypes, one of which was assumed to be transmitted together with HLA-B*73 from Denisovans by the study of Abi-Rached and colleagues, were not differentiated from those in other HLA-A-C haplotypes in modern humans. These results do not support the introgression hypothesis. Thus, we conclude that it is highly likely that HLA-B*73 allelic lineage has been maintained in the direct ancestors of modern humans.

TE-domestication and horizontal transfer in a putative Nef-AP1mu mimic of HLA-A cytoplasmic domain re-trafficking

Genes of the major histocompatibility complex (MHC; also called HLA in human) are polymorphic elements in the genomes of sharks to humans. Class-I and class-II MHC loci appear responsible for much of the genetic linkage to myriad disease states via the capacity to bind short (~8-15 a.a.) peptides of a given pathogen's proteome, or in some cases, the altered proteomes of cancerous cells, and even (in autoimmunity) certain nominal 'self' peptides (Janeway, 2004).(1) Unfortunately, little is known about how the canonical structure of the MHC-I/-II peptide-presenting gene evolved, particularly since beyond ~500 Mya (sharks) no paralogs exist.(2,3) We previously reported that HLA-A isotype alleles with the α1-helix, R65 motif, are wide-spread in phylogeny, but that the α 2-helix, H151R motif, has apparently segregated out of most species. Surprisingly, an uncharacterized orf in T. syrichta (Loc-103275158) encoded R151, but within a truncated A-23 like gene containing 5'- and 3'- footprints of the transposon (TE), tigger-1; the extant tarsier A-23 allele is totally missing exon-3 and part-of exon-4; together, suggesting TE-mediated inactivation of an intact/ancestral A-23 allele (Murray, 2015a).(4) The unique Loc-103275158 orf encodes a putative 15-exon transcript with no apparent paralogs throughout phylogeny. However, an HLA-A11 like gene in M. leucophaeus with a shortened C-terminal domain, and an HLA-A like orf in C. atys with two linked α1/α2/α3 domains, both contain a second transmembrane segment, which is conserved in Loc-103275158. Thus, we could model the putative protein with its Nef-like tail domain docked to its MHC-I like α3 domain (i.e., on the same side of a membrane). This modeled tertiary structure is strikingly similar to the solved structure of the Nef:MHC-I CD:AP1mu transporter (Jia, 2012).(5) Nef:AP1mu binds the CD of MHC-I in trafficking MHC-I away from the trans-golgi and into the endocytic pathway in HIV-1 infected cells. The CD loop of the Loc-103275158 provisional protein conserved the nominal MHC-I CD tyrosine phosphorylation site, and it has an N-terminal SH3 domain that we docked in one conformation to its internal Nef-like domain. Here, we suggest that phosphorylation of the protein's CD-loop signals an exchange between the internal Nef-like domain and a lentiviral-Nef for binding the N-terminal SH3 domain - freeing the Nef-like domain to bind MHC-I CD. Since the 5'-tigger sequence encodes part of the pseudo α1/α2 MHC-I domain, and the 3'-tigger part of the Nef-like domain, we speculate that transposition proceeded phylogenetically disparate horizontal transfers, involving adjacent 5'- and 3'- parasitic footprints, which we also found in the Loc-103275158 orf.

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.

Patterns of MHC-DRB1 polymorphism in a post-glacial island canid, the Newfoundland red fox (Vulpes vulpes deletrix), suggest balancing selection at species and population timescales

As the only native insular Newfoundland canid between the extinction of the wolf in the 1930s and the recent arrival of coyotes, the red fox (Vulpes vulpes deletrix Bangs 1898) poses interesting questions about genetic distinctiveness and the post-glacial colonization history of the island’s depauperate mammalian fauna. Here, we characterized genetic variability at the major histocompatibility complex (MHC) class II DR β1 domain (DRB1) locus in 28 red foxes from six sampling localities island-wide and compared it with mitochondrial control region (CR) diversity and DRB1 diversity in other canids. Our goals were to describe novel DRB1 alleles in a new canid population and to make inferences about the role of selection in maintaining their diversity. As in numerous studies of vertebrates, we found an order-of-magnitude higher nucleotide diversity at the DRB1 locus compared with the CR and significantly positive nonsynonymous-to-synonymous substitution ratios, indicative of selection in the distant past. Although the evidence is weaker, the Ewens-Watterson test of neutrality and the geographical distribution of variation compared with the CR suggest a role for selection over the evolutionary timescale of populations. We report the first genetic data from the DRB1 locus in the red fox and establish baseline information regarding immunogenetic variation in this island canid population which should inform continued investigations of population demography, adaptive genetic diversity, and wildlife disease in red foxes and related species.

Balancing selection and recombination as evolutionary forces caused population genetic variations in golden pheasant MHC class I genes

BACKGROUND

The major histocompatibility complex (MHC) genes are vital partners in the acquired immune processes of vertebrates. MHC diversity may be directly associated with population resistance to infectious pathogens. Here, we screened for polymorphisms in exons 2 and 3 of the IA1 and IA2 genes in 12 golden pheasant populations across the Chinese mainland to characterize their genetic variation levels, to understand the effects of historical positive selection and recombination in shaping class I diversity, and to investigate the genetic structure of wild golden pheasant populations.

RESULTS

Among 339 individual pheasants, we identified 14 IA1 alleles in exon 2 (IA1-E2), 11 IA1-E3 alleles, 27 IA2-E2 alleles, and 28 IA2-E3 alleles. The non-synonymous substitution rate was significantly greater than the synonymous substitution rate at sequences in the IA2 gene encoding putative peptide-binding sites but not in the IA1 gene; we also found more positively selected sites in IA2 than in IA1. Frequent recombination events resulted in at least 9 recombinant IA2 alleles, in accordance with the intermingling pattern of the phylogenetic tree. Although some IA alleles are widely shared among studied populations, large variation occurs in the number of IA alleles across these populations. Allele frequency analysis across 2 IA loci showed low levels of genetic differentiation among populations on small geographic scales; however, significant genetic differentiation was observed between pheasants from the northern and southern regions of the Yangtze River. Both STRUCTURE analysis and F-statistic (F ST ) value comparison classified those populations into 2 major groups: the northern region of the Yangtze River (NYR) and the southern region of the Yangtze River (SYR).

CONCLUSIONS

More extensive polymorphisms in IA2 than IA1 indicate that IA2 has undergone much stronger positive-selection pressure during evolution. Moreover, the recombination events detected between the genes and the intermingled phylogenetic pattern indicate that interlocus recombination accounts for much of the allelic variation in IA2. Analysis of the population differentiation implied that homogenous balancing selection plays an important part in maintaining an even distribution of MHC variations. The natural barrier of the Yangtze River and heterogeneous balancing selection might help shape the NYR-SYR genetic structure in golden pheasants.

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.

High genetic diversity in the endangered and narrowly distributed amphibian species Leptobrachium leishanense

Threatened species typically have a small or declining population size, which make them highly susceptible to loss of genetic diversity through genetic drift and inbreeding. Genetic diversity determines the evolutionary potential of a species; therefore, maintaining the genetic diversity of threatened species is essential for their conservation. In this study, we assessed the genetic diversity of the adaptive major histocompatibility complex (MHC) genes in an endangered and narrowly distributed amphibian species, Leptobrachium leishanense in Southwest China. We compared the genetic variation of MHC class I genes with that observed in neutral markers (5 microsatellite loci and cytochrome b gene) to elucidate the relative roles of genetic drift and natural selection in shaping the current MHC polymorphism in this species. We found a high level of genetic diversity in this population at both MHC and neutral markers compared with other threatened amphibian species. Historical positive selection was evident in the MHC class I genes. The higher allelic richness in MHC markers compared with that of microsatellite loci suggests that selection rather than genetic drift plays a prominent role in shaping the MHC variation pattern, as drift can affect all the genome in a similar way but selection directly targets MHC genes. Although demographic analysis revealed no recent bottleneck events in L. leishanense, additional population decline will accelerate the dangerous status for this species. We suggest that the conservation management of L. leishanense should concentrate on maximizing the retention of genetic diversity through preventing their continuous population decline. Protecting their living habitats and forbidding illegal hunting are the most important measures for conservation of L. leishanense.

Epigenetic regulation of differential HLA-A allelic expression levels

MHC class I expression levels influence the strength of immune responses and represent another variable in determining outcome to disease beyond peptide binding alone. Identification of the HLA loci that vary in allelic expression levels and delineating the mechanism responsible for expression variation may provide the opportunity to modify their expression therapeutically. We have examined the expression levels of allelic lineages at the HLA-A locus in a sample of 216 European Americans using a real-time polymerase chain reaction assay, which amplifies all HLA-A lineages specifically with equal efficiency, and observed a gradient of expression that associates with HLA-A allelic lineage (R = 0.6, P = 5 × 10(-25)). DNA methylation of the HLA-A gene appears to contribute to the variation in HLA-A mRNA expression levels, as a significant inverse correlation was observed between HLA-A mRNA expression levels in untreated cells and the degree to which expression is increased after treatment of the cells with a DNA methyltransferase inhibitor (R = 0.6, P = 2.8 × 10(-6)). Further, deep-sequencing and immunoprecipitation assays revealed allelic lineage-specific methylation patterns within the HLA-A promoter region where increased DNA methylation levels correlated significantly with reduced HLA-A expression levels (R = 0.89, P = 3.7 × 10(-9)). These data demonstrate HLA-A allelic lineage-specific variation in expression levels, and DNA methylation as a likely factor in contributing to this variation.

Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele

Bovine leukocyte antigens (BoLAs) are used extensively as markers for bovine disease and immunological traits. In this study, we estimated BoLA-DRB3 allele frequencies using 888 cattle from 10 groups, including seven cattle breeds and three crossbreeds: 99 Red Angus, 100 Black Angus, 81 Chilean Wagyu, 49 Hereford, 95 Hereford × Angus, 71 Hereford × Jersey, 20 Hereford × Overo Colorado, 113 Holstein, 136 Overo Colorado, and 124 Overo Negro cattle. Forty-six BoLA-DRB3 alleles were identified, and each group had between 12 and 29 different BoLA-DRB3 alleles. Overo Negro had the highest number of alleles (29); this breed is considered in Chile to be an 'Old type' European Holstein Friesian descendant. By contrast, we detected 21 alleles in Holstein cattle, which are considered to be a 'Present type' Holstein Friesian cattle. Chilean cattle groups and four Japanese breeds were compared by neighbor-joining trees and a principal component analysis (PCA). The phylogenetic tree showed that Red Angus and Black Angus cattle were in the same clade, crossbreeds were closely related to their parent breeds, and Holstein cattle from Chile were closely related to Holstein cattle in Japan. Overall, the tree provided a thorough description of breed history. It also showed that the Overo Negro breed was closely related to the Holstein breed, consistent with historical data indicating that Overo Negro is an 'Old type' Holstein Friesian cattle. This allelic information will be important for investigating the relationship between major histocompatibility complex (MHC) and disease.

Global distribution of malaria-resistant MHC-HLA alleles: the number and frequencies of alleles and malaria risk

Background

The major histocompatibility complex (MHC) is the most polymorphic genetic region in vertebrates, but the origin of such genetic diversity remains unresolved. Several studies have demonstrated at the within-population level that individuals harbouring particular alleles can be less or more susceptible to malaria, but these do not allow strong generalization.

Methods

Here worldwide data on the frequencies of several hundred MHC alleles of the human leucocyte antigen (HLA) system in relation to malaria risk at the between-population level were analysed in a phylogenetic framework, and results for different alleles were quantitatively summarized in a meta-analysis.

Results

There was an overall positive relationship between malaria pressure and the frequency of several HLA alleles indicating that allele frequencies increase in countries with strong malaria pressure. Nevertheless, considerable heterogeneity was observed across alleles, and some alleles showed a remarkable negative relationship with malaria risk. When heterogeneities were partitioned into different organization groups of the MHC, the strongest positive relationships were detected for alleles of the HLA-A and HLA-B loci, but there were also differences between MHC supertypes that constitute functionally distinct nucleotide sequences. Finally, the number of MHC alleles that are maintained within countries was also related to malaria risk.

Conclusion

Therefore, malaria represents a key selection pressure for the human MHC and has left clear evolutionary footprints on both the frequencies and the number of alleles observed in different countries.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-349) contains supplementary material, which is available to authorized users.

The relevance of HLA sequencing in population genetics studies

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

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.

Nonsynonymous substitution rate heterogeneity in the peptide-binding region among different HLA-DRB1 lineages in humans

An extraordinary diversity of amino acid sequences in the peptide-binding region (PBR) of human leukocyte antigen [HLA; human major histocompatibility complex (MHC)] molecules has been maintained by balancing selection. The process of accumulation of amino acid diversity in the PBR for six HLA genes (HLA-A, B, C, DRB1, DQB1, and DPB1) shows that the number of amino acid substitutions in the PBR among alleles does not linearly correlate with the divergence time of alleles at the six HLA loci. At these loci, some pairs of alleles show significantly less nonsynonymous substitutions at the PBR than expected from the divergence time. The same phenomenon was observed not only in the HLA but also in the rat MHC. To identify the cause for this, DRB1 sequences, a representative case of a typical nonlinear pattern of substitutions, were examined. When the amino acid substitutions in the PBR were placed with maximum parsimony on a maximum likelihood tree based on the non-PBR substitutions, heterogeneous rates of nonsynonymous substitutions in the PBR were observed on several branches. A computer simulation supported the hypothesis that allelic pairs with low PBR substitution rates were responsible for the stagnation of accumulation of PBR nonsynonymous substitutions. From these observations, we conclude that the nonsynonymous substitution rate at the PBR sites is not constant among the allelic lineages. The deceleration of the rate may be caused by the coexistence of certain pathogens for a substantially long time during HLA evolution.

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.

Current perspectives on the intensity of natural selection of MHC loci

Polymorphism of genes in the major histocompatibility complex (MHC) is believed to be maintained by balancing selection. However, direct evidence of selection has proven difficult to demonstrate. In 1994, Satta and colleagues estimated the selection intensity of the human MHC (human leukocyte antigen (HLA)) loci; however, at that time the number of HLA sequences was limited. By comparing five different methods, this study demonstrated the best way to calculate the selection coefficient, through a computer simulation study. Since the study, many HLA nucleotide sequences have been made available. Our new analysis takes advantage of these newly available sequences and compares new estimates with those of the previous study. Generally, our new results are consistent with those of the 1994 study. Our results show that, even after 20 years of exhaustive sequencing of human HLA, the number of dominant HLA alleles, on which our original estimate of selection intensity depended, appears to be conserved. Indeed, according to the frequency distribution for each HLA allele, most sequences in the database were minor or private alleles; therefore, we conclude that the selection intensities of HLA loci are at most 4.4 % even though the HLA is the prominent example on which the natural selection has been operating.

Pathogen-driven selection in the human genome

Infectious diseases and epidemics have always accompanied and characterized human history, representing one of the main causes of death. Even today, despite progress in sanitation and medical research, infections are estimated to account for about 15% of deaths. The hypothesis whereby infectious diseases have been acting as a powerful selective pressure was formulated long ago, but it was not until the availability of large-scale genetic data and the development of novel methods to study molecular evolution that we could assess how pervasively infectious agents have shaped human genetic diversity. Indeed, recent evidences indicated that among the diverse environmental factors that acted as selective pressures during the evolution of our species, pathogen load had the strongest influence. Beside the textbook example of the major histocompatibility complex, selection signatures left by pathogen-exerted pressure can be identified at several human loci, including genes not directly involved in immune response. In the future, high-throughput technologies and the availability of genetic data from different populations are likely to provide novel insights into the evolutionary relationships between the human host and its pathogens. Hopefully, this will help identify the genetic determinants modulating the susceptibility to infectious diseases and will translate into new treatment strategies.

MHC class II DQB diversity in the Japanese black bear, Ursus thibetanus japonicus

Background

The major histocompatibility complex (MHC) genes are one of the most important genetic systems in the vertebrate immune response. The diversity of MHC genes may directly influence the survival of individuals against infectious disease. However, there has been no investigation of MHC diversity in the Asiatic black bear (Ursus thibetanus). Here, we analyzed 270-bp nucleotide sequences of the entire exon 2 region of the MHC DQB gene by using 188 samples from the Japanese black bear (Ursus thibetanus japonicus) from 12 local populations.

Results

Among 185 of 188 samples, we identified 44 MHC variants that encoded 31 different amino acid sequences (allotypes) and one putative pseudogene. The phylogenetic analysis suggests that MHC variants detected from the Japanese black bear are derived from the DQB locus. One of the 31 DQB allotypes, Urth-DQB*01, was found to be common to all local populations. Moreover, this allotype was shared between the black bear on the Asian continent and the Japanese black bear, suggesting that Urth-DQB*01 might have been maintained in the ancestral black bear population for at least 300,000 years. Our findings, from calculating the ratio of non-synonymous to synonymous substitutions, indicate that balancing selection has maintained genetic variation of peptide-binding residues at the DQB locus of the Japanese black bear. From examination of genotype frequencies among local populations, we observed a considerably lower level of observed heterozygosity than expected.

Conclusions

The low level of observed heterozygosity suggests that genetic drift reduced DQB diversity in the Japanese black bear due to a bottleneck event at the population or species level. The decline of DQB diversity might have been accelerated by the loss of rare variants that have been maintained by negative frequency-dependent selection. Nevertheless, DQB diversity of the black bear appears to be relatively high compared with some other endangered mammalian species. This result suggests that the Japanese black bears may also retain more potential resistance against pathogens than other endangered mammalian species. To prevent further decline of potential resistance against pathogens, a conservation policy for the Japanese black bear should be designed to maintain MHC rare variants in each local population.

MHC variability supports dog domestication from a large number of wolves: high diversity in Asia

The process of dog domestication is still somewhat unresolved. Earlier studies indicate that domestic dogs from all over the world have a common origin in Asia. So far, major histocompatibility complex (MHC) diversity has not been studied in detail in Asian dogs, although high levels of genetic diversity are expected at the domestication locality. We sequenced the second exon of the canine MHC gene DLA-DRB1 from 128 Asian dogs and compared our data with a previously published large data set of MHC alleles, mostly from European dogs. Our results show that Asian dogs have a higher MHC diversity than European dogs. We also estimated that there is only a small probability that new alleles have arisen by mutation since domestication. Based on the assumption that all of the currently known 102 DLA-DRB1 alleles come from the founding wolf population, we simulated the number of founding wolf individuals. Our simulations indicate an effective population size of at least 500 founding wolves, suggesting that the founding wolf population was large or that backcrossing has taken place.

Emerging prion disease drives host selection in a wildlife population

Infectious diseases are increasingly recognized as an important force driving population dynamics, conservation biology, and natural selection in wildlife populations. Infectious agents have been implicated in the decline of small or endangered populations and may act to constrain population size, distribution, growth rates, or migration patterns. Further, diseases may provide selective pressures that shape the genetic diversity of populations or species. Thus, understanding disease dynamics and selective pressures from pathogens is crucial to understanding population processes, managing wildlife diseases, and conserving biological diversity. There is ample evidence that variation in the prion protein gene (PRNP) impacts host susceptibility to prion diseases. Still, little is known about how genetic differences might influence natural selection within wildlife populations. Here we link genetic variation with differential susceptibility of white-tailed deer to chronic wasting disease (CWD), with implications for fitness and disease-driven genetic selection. We developed a single nucleotide polymorphism (SNP) assay to efficiently genotype deer at the locus of interest (in the 96th codon of the PRNP gene). Then, using a Bayesian modeling approach, we found that the more susceptible genotype had over four times greater risk of CWD infection; and, once infected, deer with the resistant genotype survived 49% longer (8.25 more months). We used these epidemiological parameters in a multi-stage population matrix model to evaluate relative fitness based on genotype-specific population growth rates. The differences in disease infection and mortality rates allowed genetically resistant deer to achieve higher population growth and obtain a long-term fitness advantage, which translated into a selection coefficient of over 1% favoring the CWD-resistant genotype. This selective pressure suggests that the resistant allele could become dominant in the population within an evolutionarily short time frame. Our work provides a rare example of a quantifiable disease-driven selection process in a wildlife population, demonstrating the potential for infectious diseases to alter host populations. This will have direct bearing on the epidemiology, dynamics, and future trends in CWD transmission and spread. Understanding genotype-specific epidemiology will improve predictive models and inform management strategies for CWD-affected cervid populations.

Distinct evolutionary strategies of human leucocyte antigen loci in pathogen-rich environments

Human leucocyte antigen (HLA) loci have a complex evolution where both stochastic (e.g. genetic drift) and deterministic (natural selection) forces are involved. Owing to their extraordinary level of polymorphism, HLA genes are useful markers for reconstructing human settlement history. However, HLA variation often deviates significantly from neutral expectations towards an excess of genetic diversity. Because HLA molecules play a crucial role in immunity, this observation is generally explained by pathogen-driven-balancing selection (PDBS). In this study, we investigate the PDBS model by analysing HLA allelic diversity on a large database of 535 populations in relation to pathogen richness. Our results confirm that geographical distances are excellent predictors of HLA genetic differentiation worldwide. We also find a significant positive correlation between genetic diversity and pathogen richness at two HLA class I loci (HLA-A and -B), as predicted by PDBS, and a significant negative correlation at one HLA class II locus (HLA-DQB1). Although these effects are weak, as shown by a loss of significance when populations submitted to rapid genetic drift are removed from the analysis, the inverse relationship between genetic diversity and pathogen richness at different loci indicates that HLA genes have adopted distinct evolutionary strategies to provide immune protection in pathogen-rich environments.

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

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.

Cheetah paradigm revisited: MHC diversity in the world's largest free-ranging population

For more than two decades, the cheetah (Acinonyx jubatus) has been considered a paradigm of disease vulnerability associated with low genetic diversity, particularly at the immune genes of the major histocompatibility complex (MHC). Cheetahs have been used as a classic example in numerous conservation genetics textbooks as well as in many related scientific publications. However, earlier studies used methods with low resolution to quantify MHC diversity and/or small sample sizes. Furthermore, high disease susceptibility was reported only for captive cheetahs, whereas free-ranging cheetahs show no signs of infectious diseases and a good general health status. We examined whether the diversity at MHC class I and class II-DRB loci in 149 Namibian cheetahs was higher than previously reported using single-strand conformation polymorphism analysis, cloning, and sequencing. MHC genes were examined at the genomic and transcriptomic levels. We detected ten MHC class I and four class II-DRB alleles, of which nine MHC class I and all class II-DRB alleles were expressed. Phylogenetic analyses and individual genotypes suggested that the alleles belong to four MHC class I and three class II-DRB putative loci. Evidence of positive selection was detected in both MHC loci. Our study indicated that the low number of MHC class I alleles previously observed in cheetahs was due to a smaller sample size examined. On the other hand, the low number of MHC class II-DRB alleles previously observed in cheetahs was further confirmed. Compared with other mammalian species including felids, cheetahs showed low levels of MHC diversity, but this does not seem to influence the immunocompetence of free-ranging cheetahs in Namibia and contradicts the previous conclusion that the cheetah is a paradigm species of disease vulnerability.

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

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

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

Antagonistic coevolution with parasites maintains host genetic diversity: an experimental test

Genetic variation in natural populations is a prime prerequisite allowing populations to respond to selection, but is under constant threat from forces that tend to reduce it, such as genetic drift and many types of selection. Haldane emphasized the potential importance of parasites as a driving force of genetic diversity. His theory has been taken for granted ever since, but despite numerous studies showing correlations between genetic diversity and parasitism, Haldane's hypothesis has rarely been tested experimentally for unambiguous support. We experimentally staged antagonistic coevolution between the host Tribolium castaneum and its natural microsporidian parasite, Nosema whitei, to test for the relative importance of two separate evolutionary forces (drift and parasite-induced selection) on the maintenance of genetic variation. Our results demonstrate that coevolution with parasites indeed counteracts drift as coevolving populations had significantly higher levels of heterozygosity and allelic diversity. Genetic drift remained a strong force, strongly reducing genetic variation and increasing genetic differentiation in small populations. To our surprise, differentiation between the evolving populations was smaller when they coevolved with parasites, suggesting parallel balancing selection. Hence, our results experimentally vindicate Haldane's original hypothesis 60 years after its conception.

Human genetic differentiation across the Strait of Gibraltar

BACKGROUND

The Strait of Gibraltar is a crucial area in the settlement history of modern humans because it represents a possible connection between Africa and Europe. So far, genetic data were inconclusive about the fact that this strait constitutes a barrier to gene flow, as previous results were highly variable depending on the genetic locus studied. The present study evaluates the impact of the Gibraltar region in reducing gene flow between populations from North-Western Africa and South-Western Europe, by comparing formally various genetic loci. First, we compute several statistics of population differentiation. Then, we use an original simulation approach in order to infer the most probable evolutionary scenario for the settlement of the area, taking into account the effects of both demography and natural selection at some loci.

RESULTS

We show that the genetic patterns observed today in the region of the Strait of Gibraltar may reflect an ancient population genetic structure which has not been completely erased by more recent events such as Neolithic migrations. Moreover, the differences observed among the loci (i.e. a strong genetic boundary revealed by the Y-chromosome polymorphism and, at the other extreme, no genetic differentiation revealed by HLA-DRB1 variation) across the strait suggest specific evolutionary histories like sex-mediated migration and natural selection. By considering a model of balancing selection for HLA-DRB1, we here estimate a coefficient of selection of 2.2% for this locus (although weaker in Europe than in Africa), which is in line with what was estimated from synonymous versus non-synonymous substitution rates. Selection at this marker thus appears strong enough to leave a signature not only at the DNA level, but also at the population level where drift and migration processes were certainly relevant.

CONCLUSIONS

Our multi-loci approach using both descriptive analyses and Bayesian inferences lead to better characterize the role of the Strait of Gibraltar in the evolution of modern humans. We show that gene flow across the Strait of Gibraltar occurred at relatively high rates since pre-Neolithic times and that natural selection and sex-bias migrations distorted the demographic signal at some specific loci of our genome.

Structure and evolution of a new avian MHC class II B gene in a sub-Antarctic seabird, the thin-billed prion (Procellariiformes: Pachyptila belcheri)

The major histocompatibility complex encodes molecules that present foreign peptides to T cells of the immune system. The peptide binding region (PBR) of these molecules is among the most polymorphic regions found in vertebrate taxa. Genomic cloning approaches are improving our understanding of the evolution of this multigene family in nonmodel avian groups. By building a cosmid library, a new MHC class II B gene, Pabe-DAB1, was isolated and characterized at the genomic level in a sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri). Pabe-DAB1 exhibits the hallmark structural features of functional MHC class II loci. Direct sequencing of the PBR encoding exon in a panel of prions revealed significantly higher levels of genetic diversity compared to two noncoding neutral loci, with most alleles differing by at least one replacement substitution in the peptide binding codons. We estimated evolutionary dynamics for Pabe-DAB1 using a variety of Bayesian and other approaches. Evidence for balancing selection comes from a spatially variable ratio of nonsynonymous-to-synonymous substitutions (mean d (N)/d (S) = 2.87) in the PBR, with sites predicted to be functionally relevant exhibiting the highest omega values. We estimate the population recombination rate to be approximately 0.3 per site per generation, indicating an important role for recombination in generating polymorphism at this locus. Pabe-DAB1 is among the few avian class II loci characterized at the genomic level and with a known intron-exon structure, a feature that greatly facilitated the amplification and sequencing of a single MHC locus in this species.

Patterns of nucleotide sequence variation in ICAM1 and TNF genes in twelve ethnic groups of India: roles of demographic history and natural selection

We have studied DNA sequence variation in and around the genes ICAM1 and TNF, which play functional and correlated roles in inflammatory processes and immune cell responses, in 12 diverse ethnic groups of India, with a view to investigating the relative roles of demographic history and natural selection in shaping the observed patterns of variation. The total numbers of single nucleotide polymorphisms (SNPs) detected at the ICAM1 and TNF loci were 29 and 12, respectively. Haplotype and allele frequencies differed significantly across populations. The site frequency spectra at these loci were significantly different from those expected under neutrality, and showed an excess of intermediate-frequency variants consistent with balancing selection. However, as expected under balancing selection, there was no significant reduction of F(ST) values compared to neutral autosomal loci. Mismatch distributions were consistent with population expansion for both loci. On the other hand, the phylogenetic network among haplotypes for the TNF locus was similar to expectations under population expansion, while that for the ICAM1 was as expected under balancing selection. Nucleotide diversity at the ICAM1 locus was an order of magnitude lower in the promoter region, compared to the introns or exons, but no such difference was noted for the TNF gene. Thus, we conclude that the pattern of nucleotide variation in these genes has been modulated by both demographic history and selection. This is not surprising in view of the known allelic associations of several polymorphisms in these genes with various diseases, both infectious and noninfectious.

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.

Detecting the Genetic Signature of Natural Selection in Human Populations: Models, Methods, and Data

Patterns of DNA sequence variation in the genome contain a record of past selective events. The ability to collect increasingly large data sets of polymorphisms has allowed investigators to perform hypothesis-driven studies of candidate genes as well as genome-wide scans for signatures of adaptations. This genetic approach to the study of natural selection has identified many signals consistent with predictions from anthropological studies. Selective pressures related to variation in climate, diet, and pathogen exposure have left strong marks on patterns of human variation. Additional signals of adaptations are observed in genes involved in chemosensory perception and reproduction. Several ongoing projects aim to sequence the complete genome of 1000 individuals from different human populations. These large-scale projects will provide data for more complete genome scans of selection, but more focused studies aimed at testing specific hypotheses will continue to hold an important place in elucidating the history of adaptations in humans.