MHC polymorphism under host-pathogen coevolution

Heterozygote advantage can explain the extraordinary diversity of immune genes

The majority of highly polymorphic genes are related to immune functions and with over 100 alleles within a population, genes of the major histocompatibility complex (MHC) are the most polymorphic loci in vertebrates. How such extraordinary polymorphism arose and is maintained is controversial. One possibility is heterozygote advantage (HA), which can in principle maintain any number of alleles, but biologically explicit models based on this mechanism have so far failed to reliably predict the coexistence of significantly more than 10 alleles. We here present an eco-evolutionary model showing that evolution can result in the emergence and maintenance of more than 100 alleles under HA if the following two assumptions are fulfilled: first, pathogens are lethal in the absence of an appropriate immune defence; second, the effect of pathogens depends on host condition, with hosts in poorer condition being affected more strongly. Thus, our results show that HA can be a more potent force in explaining the extraordinary polymorphism found at MHC loci than currently recognised.

Europeans and Americans of European origin show differences between their biological pathways related to the major histocompatibility complex

In this study, we analysed biological pathway diversity among Europeans and Northern Americans of European origin, the groups of people that share a common genetic ancestry but live in different geographic regions. We used a novel complex approach for analysing genomic data: we studied the total effects of multiple weak selection signals, accumulated from independent SNPs within a pathway. We found significant differences between immunity-related biological pathways from the two groups. All identified pathways included genes belonging to the major histocompatibility complex (MHC) system, which plays an important role in adaptive immune responses. We suggest that the ways of evolution were different for the MHC-I and MHC-II gene groups at least in Europeans and Americans of European origin. We hypothesise that the observed variability between the two populations was triggered by selection pressures due to the different pathogen landscapes and pathogen loads on the two continents. Our findings can be important for epidemic prevention and control, as well as for analysing processes related to allergies, organ transplantation, and autoimmune diseases.

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

Organisation and evolution of the major histocompatibility complex class I genes in cetaceans

A quarter of marine mammals are at risk of extinction, with disease and poor habitat quality contributing to population decline. Investigation of the Major Histocompatibility Complex (MHC) provides insight into species' capacity to respond to immune and environmental challenges. The eighteen available cetacean chromosome level genomes were used to annotate MHC Class I loci, and to reconstruct the phylogenetic relationship of the described loci. The highest number of loci was observed in the striped dolphin (Stenella coeruleoalba), while the least was observed in the pygmy sperm whale (Kogia breviceps) and rough toothed dolphin (Steno bredanensis). Of the species studied, Mysticetes had the most pseudogenes. Evolutionarily, MHC Class I diverged before the speciation of cetaceans. Yet, locus one was genomically and phylogenetically similar in many species, persisting over evolutionary time. This characterisation of MHC Class I in cetaceans lays the groundwork for future population genetics and MHC expression studies.

Diverse cytomegalovirus US11 antagonism and MHC-A evasion strategies reveal a tit-for-tat coevolutionary arms race in hominids

Recurrent, ancient arms races between viruses and hosts have shaped both host immunological defense strategies as well as viral countermeasures. One such battle is waged by the glycoprotein US11 encoded by the persisting human cytomegalovirus. US11 mediates degradation of major histocompatibility class I (MHC-I) molecules to prevent CD8+ T-cell activation. Here, we studied the consequences of the arms race between US11 and primate MHC-A proteins, leading us to uncover a tit-for-tat coevolution and its impact on MHC-A diversification. We found that US11 spurred MHC-A adaptation to evade viral antagonism: In an ancestor of great apes, the MHC-A A2 lineage acquired a Pro184Ala mutation, which confers resistance against the ancestral US11 targeting strategy. In response, US11 deployed a unique low-complexity region (LCR), which exploits the MHC-I peptide loading complex to target the MHC-A2 peptide-binding groove. In addition, the global spread of the human HLA-A*02 allelic family prompted US11 to employ a superior LCR strategy with an optimally fitting peptide mimetic that specifically antagonizes HLA-A*02. Thus, despite cytomegaloviruses low pathogenic potential, the increasing commitment of US11 to MHC-A has significantly promoted diversification of MHC-A in hominids.

Secrets and lies of host-microbial interactions: MHC restriction and trans-regulation of T cell trafficking conceal the role of microbial agents on the edge between health and multifactorial/complex diseases

Here we critically discuss data supporting the view that microbial agents (pathogens, pathobionts or commensals alike) play a relevant role in the pathogenesis of multifactorial diseases, but their role is concealed by the rules presiding over T cell antigen recognition and trafficking. These rules make it difficult to associate univocally infectious agents to diseases' pathogenesis using the paradigm developed for canonical infectious diseases. (Cross-)recognition of a variable repertoire of epitopes leads to the possibility that distinct infectious agents can determine the same disease(s). There can be the need for sequential infection/colonization by two or more microorganisms to develop a given disease. Altered spreading of infectious agents can determine an unwanted activation of T cells towards a pro-inflammatory and trafficking phenotype, due to differences in the local microenvironment. Finally, trans-regulation of T cell trafficking allows infectious agents unrelated to the specificity of T cell to modify their homing to target organs, thereby driving flares of disease. The relevant role of microbial agents in largely prevalent diseases provides a conceptual basis for the evaluation of more specific therapeutic approaches, targeted to prevent (vaccine) or cure (antibiotics and/or Biologic Response Modifiers) multifactorial diseases.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

Integration of the immune memory into the pathogen-driven MHC polymorphism hypothesis

Major histocompatibility complex (MHC) genes (referred to as human leukocyte antigen or HLA in humans) are a key component of vertebrate immune systems, coding for proteins which present antigens to T-cells. These genes are outstanding in their degree of polymorphism, with important consequences for human and animal health. The polymorphism is thought to arise from selection pressures imposed by pathogens on MHC allomorphs, which differ in their antigen-binding capacity. However, the existing theory has not considered MHC selection in relation to the formation of immune memory. In this paper, we argue that this omission limits our understanding of the evolution of MHC polymorphism and its role in disease. We review recent evidence that has emerged from the massive research effort related to the SARS-CoV-2 pandemics, and which provides new evidence for the role of MHC in shaping immune memory. We then discuss why the inclusion of immune memory within the existing theory may have non-trivial consequence for our understanding of the evolution of MHC polymorphism. Finally, we will argue that neglecting immune memory hinders our interpretation of empirical findings, and postulate that future studies focusing on pathogen-driven MHC selection would benefit from stratifying the available data according to the history of infection (and vaccination, if relevant).

The maintenance of genetic diversity under host-parasite coevolution in finite, structured populations

As a corollary to the Red Queen hypothesis, host-parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host-parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Cytokine gene polymorphism and parasite susceptibility in free-living rodents: Importance of non-coding variants

Associations between genetic variants and susceptibility to infections have long been studied in free-living hosts so as to infer the contemporary evolutionary forces that shape the genetic polymorphisms of immunity genes. Despite extensive studies of proteins interacting with pathogen-derived ligands, such as MHC (major histocompatilbility complex) or TLR (Toll-like receptors), little is known about the efferent arm of the immune system. Cytokines are signalling molecules that trigger and modulate the immune response, acting as a crucial link between innate and adaptive immunity. In the present study we investigated how genetic variation in cytokines in bank voles Myodes glareolus affects their susceptibility to infection by parasites (nematodes: Aspiculuris tianjensis, Heligmosomum mixtum, Heligmosomoides glareoli) and microparasites (Cryptosporidium sp, Babesia microti, Bartonella sp.). We focused on three cytokines: tumour necrosis factor (TNF), lymphotoxin alpha (LTα), and interferon beta (IFNβ1). Overall, we identified four single nucleotide polymorphisms (SNPs) associated with susceptibility to nematodes: two located in LTα and two in IFNβ1. One of those variants was synonymous, another located in an intron. Each SNP associated with parasite load was located in or next to a codon under selection, three codons displayed signatures of positive selection, and one of purifying selection. Our results indicate that cytokines are prone to parasite-driven selection and that non-coding variants, although commonly disregarded in studies of the genetic background of host-parasite co-evolution, may play a role in susceptibility to infections in wild systems.

Human leukocyte antigen class II gene diversity tunes antibody repertoires to common pathogens

Allelic diversity of human leukocyte antigen (HLA) class II genes may help maintain humoral immunity against infectious diseases. In this study, we investigated germline genetic variation in classical HLA class II genes and employed a systematic, unbiased approach to explore the relative contribution of this genetic variation in the antibody repertoire to various common pathogens. We leveraged a well-defined cohort of 800 adults representing the general Arab population in which genetic material is shared because of the high frequency of consanguineous unions. By applying a high-throughput method for large-scale antibody profiling to this well-defined cohort, we were able to dissect the overall effect of zygosity for classical HLA class II genes, as well as the effects associated with specific HLA class II alleles, haplotypes and genotypes, on the antimicrobial antibody repertoire breadth and antibody specificity with unprecedented resolution. Our population genetic studies revealed that zygosity of the classical HLA class II genes is a strong predictor of antibody responses to common human pathogens, suggesting that classical HLA class II gene heterozygosity confers a selective advantage. Moreover, we demonstrated that multiple HLA class II alleles can have additive effects on the antibody repertoire to common pathogens. We also identified associations of HLA-DRB1 genotypes with specific antigens. Our findings suggest that HLA class II gene polymorphisms confer specific humoral immunity against common pathogens, which may have contributed to the genetic diversity of HLA class II loci during hominine evolution.

Comparative assessment of immunological tolerance in fish with natural immunodeficiency

Natural occurrences of immunodeficiency by definition should lead to compromised immune function. The major histocompatibility complexes (MHC) are key components of the vertebrate adaptive immune system, charged with mediating allorecognition and antigen presentation functions. To this end, the genomic loss of the MHC II pathway in Syngnathus pipefishes raises questions regarding their immunological vigilance and allorecognition capabilities. Utilising allograft and autograft fin-transplants, we compared the allorecognition immune responses of two pipefish species, with (Nerophis ophidion) and without (Syngnathus typhle) a functional MHC II. Transcriptome-wide assessments explored the immunological tolerance and potential compensatory measures occupying the role of the absent MHC II. Visual observations suggested a more acute rejection response in N. ophidion allografts compared with S. typhle allografts. Differentially expressed genes involved in innate immunity, angiogenesis and tissue recovery were identified among transplantees. The intriguing upregulation of the cytotoxic T-cell implicated gzma in S. typhle allografts, suggests a prominent MHC I related response, which may compensate for the MHC II and CD4 loss. MHC I related downregulation in N. ophidion autografts hints at an immunological tolerance related reaction. These findings may indicate alternative measures evolved to cope with the MHC II genomic loss enabling the maintenance of appropriate tolerance levels. This study provides intriguing insights into the immune and tissue recovery mechanisms associated with syngnathid transplantation, and can be a useful reference for future studies focusing on transplantation transcriptomics in non-model systems.

How HLA diversity is apportioned: influence of selection and relevance to transplantation

In his 1972 paper ‘The apportionment of human diversity’, Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low F_ST values. However, when we analyse haplotypes defined by these SNPs (which define ‘HLA alleles’), we find marked differences in frequencies between geographic regions. These differences are not reflected in the F_ST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting F_ST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

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.

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer

BACKGROUND

Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range.

RESULTS

We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables.

CONCLUSIONS

Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

Evolutionary pressures rendered by animal husbandry practices for avian influenza viruses to adapt to humans

Commercial poultry operations produce and crowd billions of birds every year, which is a source of inexpensive animal protein. Commercial poultry is intensely bred for desirable production traits, and currently presents very low variability at the major histocompatibility complex. This situation dampens the advantages conferred by the MHC’s high genetic variability, and crowding generates immunosuppressive stress. We address the proteins of influenza A viruses directly and indirectly involved in host specificities. We discuss how mutants with increased virulence and/or altered host specificity may arise if few class I alleles are the sole selective pressure on avian viruses circulating in immunocompromised poultry. This hypothesis is testable with peptidomics of MHC ligands. Breeding strategies for commercial poultry can easily and inexpensively include high variability of MHC as a trait of interest, to help save billions of dollars as a disease burden caused by influenza and decrease the risk of selecting highly virulent strains.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

Urbanization processes drive divergence at the major histocompatibility complex in a common waterbird

Urban sprawl is one of the most common landscape alterations occurring worldwide, and there is a growing list of species that are recognised to have adapted to urban life. To be successful, processes of urban colonization by wildlife require a broad spectrum of phenotypic (e.g., behavioural or physiological) adjustments, but evidence for genetic adaptations is much scarcer. One hypothesis proposes that different pathogen-driven selective pressures between urban and non-urban landscapes leads to adaptations in host immune genes. Here, we examined urbanization-related differentiation at the key pathogen-recognition genes of vertebrate adaptive immunity-the major histocompatibility complex (MHC)-in a common waterbird, the Eurasian coot (Fulica atra). Samples were collected from an old urban population (established before the 1950s), a new urban population (established in the 2000s), and two rural populations from central Poland. We found strong significant divergence (as measured with Jost's D) at the MHC class II between the old urban population and the remaining (new urban and rural) populations. Also, there was a moderate, but significant divergence at the MHC between the new urban population and two rural populations, while no divergence was found between the two rural populations. The total number of MHC alleles and the number of private (population-specific) MHC alleles was lower in old urban populations, as compared to the rural ones. These patterns of differentiation at the MHC were not consistent with patterns found for neutral genetic markers (microsatellites), which showed few differences between the populations. Our results indicate that MHC allele composition depended on the level of anthropogenic disturbance and the time which passed since urban colonization, possibly due to the processes of genotype sorting and local adaptation. As such, our study contributes to the understanding of genetic mechanisms associated with urbanization processes in wildlife.

Comparative Analysis of SLA-1, SLA-2, and DQB1 Genetic Diversity in Locally-Adapted Kenyan Pigs and Their Wild Relatives, Warthogs

Swine leukocyte antigen (SLA) plays a central role in controlling the immune response by discriminating self and foreign antigens and initiating an immune response. Studies on SLA polymorphism have demonstrated associations between SLA allelic variants, immune response, and disease resistance. The SLA polymorphism is due to host-pathogen co-evolution resulting in improved adaptation to diverse environments making SLA a crucial genomic region for comparative diversity studies. Although locally-adapted African pigs have small body sizes, they possess increased resilience under harsh environmental conditions and robust immune systems with reported tolerance to some diseases, including African swine fever. However, data on the SLA diversity in these pigs are not available. We characterized the SLA of unrelated locally-adapted domestic pigs from Homa Bay, Kenya, alongside exotic pigs and warthogs. We undertook SLA comparative diversity of the functionally expressed SLA class I (SLA-1, SLA-2) and II (DQB1) repertoires in these three suids using the reverse transcription polymerase chain reaction (RT-PCR) sequence-based typing (SBT) method. Our data revealed higher genetic diversity in the locally-adapted pigs and warthogs compared to the exotic pigs. The nucleotide substitution rates were higher in the peptide-binding regions of the SLA-1, SLA-2, and DQB1 loci, indicative of adaptive evolution. We obtained high allele frequencies in the three SLA loci, including some breed-specific private alleles, which could guide breeders to increase their frequency through selection if confirmed to be associated with enhanced resilience. Our study contributes to the growing body of knowledge on genetic diversity in free-ranging animal populations in their natural environment, availing the first DQB1 gene data from locally-adapted Kenyan pigs.

Diversity of MHC IIB genes and parasitism in hybrids of evolutionarily divergent cyprinoid species indicate heterosis advantage

The genes of the major histocompatibility complex (MHC) are an essential component of the vertebrate immune system and MHC genotypes may determine individual susceptibility to parasite infection. In the wild, selection that favors MHC variability can create situations in which interspecies hybrids experience a survival advantage. In a wild system of two naturally hybridizing leuciscid fish, we assessed MHC IIB genetic variability and its potential relationships to hosts' ectoparasite communities. High proportions of MHC alleles and parasites were species-specific. Strong positive selection at specific MHC codons was detected in both species and hybrids. MHC allele expression in hybrids was slightly biased towards the maternal species. Controlling for a strong seasonal effect on parasite communities, we found no clear associations between host-specific parasites and MHC alleles or MHC supertypes. Hybrids shared more MHC alleles with the more MHC-diverse parental species, but expressed intermediate numbers of MHC alleles and positively selected sites. Hybrids carried significantly fewer ectoparasites than either parent species, suggesting a hybrid advantage via potential heterosis.

Hybridization with mountain hares increases the functional allelic repertoire in brown hares

Brown hares (Lepus europaeus Pallas) are able to hybridize with mountain hares (L. timidus Linnaeus) and produce fertile offspring, which results in cross-species gene flow. However, not much is known about the functional significance of this genetic introgression. Using targeted sequencing of candidate loci combined with mtDNA genotyping, we found the ancestral genetic diversity in the Finnish brown hare to be small, likely due to founder effect and range expansion, while gene flow from mountain hares constitutes an important source of functional genetic variability. Some of this variability, such as the alleles of the mountain hare thermogenin (uncoupling protein 1, UCP1), might have adaptive advantage for brown hares, whereas immunity-related MHC alleles are reciprocally exchanged and maintained via balancing selection. Our study offers a rare example where an expanding species can increase its allelic variability through hybridization with a congeneric native species, offering a route to shortcut evolutionary adaptation to the local environmental conditions.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Identifying HLA DRB1-DQB1 alleles associated with Chlamydia trachomatis infection and in silico prediction of potentially-related peptides

HLA class II (HLA-II) genes' polymorphism influences the immune response to Chlamydia trachomatis (Ct), it is considered a sexually transmitted infection. However, associations between HLA-II alleles and Ct-infection have been little explored in humans; this study was thus aimed at determining HLA-DRB1-DQB1 alleles/haplotypes' effect on Ct-infection outcome in a cohort of Colombian women. Cervical sample DNA was used as template for detecting Ct by PCR and typing HLA-DRB1-DQB1 alleles/haplotypes by Illumina MiSeq sequencing. Survival models were adjusted for identifying the alleles/haplotypes' effect on Ct-outcome; bioinformatics tools were used for predicting secreted bacterial protein T- and B-cell epitopes. Sixteen HLA-DRB1 alleles having a significant effect on Ct-outcome were identified in the 262 women analysed. DRB1*08:02:01G and DRB1*12:01:01G were related to infection-promoting events. Only the DQB1*05:03:01G allele related to clearance/persistence events was found for HLA-DQB1. HLA-DRB1 allele homozygous women were associated with events having a lower probability of clearance and/or early occurrence of persistence. Twenty-seven peptides predicted in silico were associated with protective immunity against Ct; outer membrane and polymorphic membrane protein-derived peptides had regions having dual potential for being T- or B-cell epitopes. This article describes HLA-DRB1-DQB1 alleles/haplotypes related to Ct-infection resolution and the peptides predicted in silico which might probably be involved in host immune response. The data provides base information for developing future studies leading to the development of effective prevention measures against Ct-infection.

Systematic screening of viral and human genetic variation identifies antiretroviral resistance and immune escape link

Background

Considering the remaining threat of drug-resistantmutations (DRMs) to antiretroviral treatment (ART) efficacy, we investigated how the selective pressure of human leukocyte antigen (HLA)-restricted cytotoxic T lymphocytes drives certain DRMs' emergence and retention.

Methods

We systematically screened DRM:HLA class I allele combinations in 3997 ART-naïve Swiss HIV Cohort Study (SHCS) patients. For each pair, a logistic regression model preliminarily tested for an association with the DRM as the outcome. The three HLA:DRM pairs remaining after multiple testing adjustment were analyzed in three ways: cross-sectional logistic regression models to determine any HLA/infection time interaction, survival analyses to examine if HLA type correlated with developing specific DRMs, and via NetMHCpan to find epitope binding evidence of immune escape.

Results

Only one pair, RT-E138:HLA-B18, exhibited a significant interaction between infection duration and HLA. The survival analyses predicted two pairs with an increased hazard of developing DRMs: RT-E138:HLA-B18 and RT-V179:HLA-B35. RT-E138:HLA-B18 exhibited the greatest significance in both analyses (interaction term odds ratio [OR] 1.169 [95% confidence interval (CI) 1.075-1.273]; p-value<0.001; survival hazard ratio 12.211 [95% CI 3.523-42.318]; p-value<0.001). The same two pairs were also predicted by netMHCpan to have epitopic binding.

Conclusions

We identified DRM:HLA pairs where HLA presence is associated with the presence or emergence of the DRM, indicating that the selective pressure for these mutations alternates direction depending on the presence of these HLA alleles.

Funding

Funded by the Swiss National Science Foundation within the framework of the SHCS, and the University of Zurich, University Research Priority Program: Evolution in Action: From Genomes Ecosystems, in Switzerland.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback?

Major Histocompatibility Complex (MHC) genes code for proteins that recognize foreign protein antigens to initiate T-cell-mediated adaptive immune responses. They are often the most polymorphic genes in vertebrate genomes. How evolution maintains this diversity remains of debate. Three main hypotheses seek to explain the maintenance of MHC diversity by invoking pathogen-mediated selection: heterozygote advantage, frequency-dependent selection, and fluctuating selection across landscapes or through time. Here, we use a large-scale field parasite survey in a stickleback metapopulation to test predictions derived from each of these hypotheses. We identify over 1000 MHC IIβ variants (alleles spanning paralogous genes) and find that many of them covary positively or negatively with parasite load, suggesting that these genes contribute to resistance or susceptibility. However, despite our large sample-size, we find no evidence for the widely cited stabilizing selection on MHC heterozygosity, in which individuals with an intermediate number of MHC variants have the lowest parasite burden. Nor do we observe a rare-variant advantage, or widespread fluctuating selection across populations. In contrast, we find that MHC diversity is best predicted by neutral genome-wide heterozygosity and between-population genomic divergence, suggesting neutral processes are important in shaping the pattern of metapopulation MHC diversity. Thus, although MHC IIβ is highly diverse and relevant to the type and intensity of macroparasite infection in these populations of stickleback, the main models of MHC evolution still provide little explanatory power in this system.

Coevolution fails to maintain genetic variation in a host-parasite model with constant finite population size

Coevolutionary negative frequency-dependent selection has been hypothesized to maintain genetic variation in host and parasites. Despite the extensive literature pertaining to host-parasite coevolution, the temporal dynamics of genetic variation have not been examined in a matching-alleles model (MAM) with a finite population size relative to the expectation under neutral genetic drift alone. The dynamics of the MA coevolution in an infinite population, in fact, suggests that genetic variation in these coevolving populations behaves neutrally. By comparing host heterozygosity to the expectation in a single-species model of neutral genetic drift we find that while this is also largely true in finite populations two additional phenomena arise. First, reciprocal natural selection acting on stochastic perturbations in host and pathogen allele frequencies results in a slight increase or decrease in genetic variation depending on the parameter conditions. Second, following the fixation of an allele in the parasite, selection in the MAM becomes directional, which then rapidly erodes genetic variation in the host. Hence, rather than maintain it, we find that, on average, matching-alleles coevolution depletes genetic variation.

Genetic origins and diversity of bushpigs from Madagascar (Potamochoerus larvatus, family Suidae)

The island of Madagascar, situated off the southeast coast of Africa, shows the first evidence of human presence ~ 10,000 years ago; however, other archaeological data indicates a settlement of the modern peoples of the island distinctly more recent, perhaps > 1500 years ago. Bushpigs of the genus Potamochoerus (family Suidae), are today widely distributed in Madagascar and presumed to have been introduced from Africa at some stage by human immigrants to the island. However, disparities about their origins in Madagascar have been presented in the literature, including the possibility of endemic subspecies, and few empirical data are available. Furthermore, the separation of bushpigs in Madagascar from their mainland relatives may have favoured the evolution of a different repertoire of immune genes first due to a founder effect and then as a response to distinct pathogens compared to their ancestors. Molecular analysis confirmed the species status of the bushpig in Madagascar as P. larvatus, likely introduced from the central region of southern Africa, with no genetic evidence for the recognition of eastern and western subspecies as suggested from previous cranial morphology examination. Investigation of the immunologically important SLA-DQB1 peptide-binding region showed a different immune repertoire of bushpigs in Madagascar compared to those on the African mainland, with seventeen exon-2 haplotypes unique to bushpigs in Madagascar (2/28 haplotypes shared). This suggests that the MHC diversity of the Madagascar populations may have enabled Malagasy bushpigs to adapt to new environments.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Balancing selection versus allele and supertype turnover in MHC class II genes in guppies

Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but F_ST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

HIV Transmission Chains Exhibit Greater HLA-B Homogeneity Than Randomly Expected

BACKGROUND

HIV's capacity to escape immune recognition by human leukocyte antigen (HLA) is a core component of HIV pathogenesis. A better understanding of the distribution of HLA class I in HIV-infected patients would improve our knowledge of pathogenesis in relation to the host HLA type and could better improve therapeutic strategies against HIV.

MATERIALS AND METHODS

Three hundred one to 325 transmission pairs and 469-496 clusters were identified for analysis among Swiss HIV Cohort Study (SHCS) participants using HIV pol sequences from the drug resistance database. HLA class I data were compiled at 3 specificity levels: 4-digit, 2-digit alleles, and HLA-B supertype. The analysis tabulated HLA-I homogeneity as 2 measures: the proportion of transmission pairs, which are HLA concordant, and the average percentage of allele matches within all clusters. These measures were compared with the mean value across randomizations with randomly assorted individuals.

RESULTS

We repeated the analysis for different HLA classification levels and separately for HLA-A, -B, and -C. Subanalyses by the risk group were performed for HLA-B. HLA-B showed significantly greater homogeneity in the transmission chains (2-digit clusters: 0.291 vs. 0.251, P value = 0.009; supertype clusters: 0.659 vs. 0.611, P value = 0.002; supertype pairs: 0.655 vs. 0.608, P value = 0.014). Risk group restriction caused the effect to disappear for men-who-have-sex-with-men but not for other risk groups. We also examined if protective HLA alleles B27 and B57 were under- or overrepresented in the transmission chains, although this yielded no significant pattern.

CONCLUSIONS

The HLA-B alleles of patients within HIV-1 transmission chains segregate in homogenous clusters/pairs, potentially indicating preferential transmission among HLA-B concordant individuals.

In silico study of the association of the HLA-A*31:01 allele (human leucocyte antigen allele 31:01) with neuroantigenic epitopes of PLP (proteolipid protein), MBP (myelin basic protein) and MOG proteins (myelin oligodendrocyte glycoprotein) for studying the multiple sclerosis disease pathogenesis

The main pathologic hallmark of multiple sclerosis is a demyelinating plaque that contains a prominent immunologic response dominated by T cells of the immune system. PLP (proteolipid protein), MPB (myelin basic protein), and Myelin oligodendrocyte glycoprotein (MOG) proteins are important autoantigens for the demyelinating of CNS in multiple sclerosis. There is good evidence indicating that T CD8⁺ cells and MHC class I molecules play an important role in this disease. The HLA-A*31:01 allele of MHC class I is a member of HLA-A3 superfamily and there is no clear report concerning the relationship of this allele with MS. Feeling this gap, we studied the possible association of the HLA-A*31:01 with MS by prediction of neuroantigenic epitopes of human MBP, PLP, and MOG proteins of myelin sheath using in silico methods. PLP did not show any neuroantigenic epitope, but the two epitopes of MBP and seven epitopes of MOG for HLA-A*31:01 were determined via bioinformatics servers. In silico study of the nine epitope showed that MOG195-204 (LIICYNWLHR) peptide of the membrane-associated/cytoplasmic part of human MOG has suitable binding affinity to the HLA-A*31:01 allele as a potential neuroantigenic epitope. Further investigations of this peptide revealed that the binding of C-terminal residue of this peptide has a more significant effect on binding to this allele than the N-terminal part of the peptide. Altogether, this combination of "LIICYNWLHR/A*31:01 allele "may play an important role in MS pathogenesis and this complex is suggested for further studies such as T cell receptor.Communicated by Ramaswamy H. Sarma.

Mating preferences can drive expansion or contraction of major histocompatibility complex gene family

Major histocompatibility complex (MHC)-based mating rules can evolve as a way to avoid inbreeding or to increase offspring immune competence. While the role of mating preference in shaping the MHC diversity in vertebrates has been acknowledged, its impact on individual MHC diversity has not been considered. Here, we use computer simulations to investigate how simple mating rules favouring MHC-dissimilar partners affect the evolution of the number of MHC variants in individual genomes, accompanying selection for resistance to parasites. We showed that the effect of such preferences could sometimes be dramatic. If preferences are aimed at avoiding identical alleles, the equilibrium number of MHC alleles is much smaller than under random mating. However, if the mating rule minimizes the ratio of shared to different alleles in partners, MHC number is higher than under random mating. Additionally, our simulations revealed that a negative correlation between the numbers of MHC variants in mated individuals can arise from simple rules of MHC-disassortative mating. Our results reveal unexpected potential of MHC-based mating preferences to drive MHC gene family expansions or contractions and highlight the need to study the mechanistic basis of such preferences, which is currently poorly understood.

Diversity of HLA Class I and Class II blocks and conserved extended haplotypes in Lacandon Mayans

Here we studied HLA blocks and haplotypes in a group of 218 Lacandon Maya Native American using a high-resolution next generation sequencing (NGS) method. We assessed the genetic diversity of HLA class I and class II in this population, and determined the most probable ancestry of Lacandon Maya HLA class I and class II haplotypes. Importantly, this Native American group showed a high degree of both HLA homozygosity and linkage disequilibrium across the HLA region and also lower class II HLA allelic diversity than most previously reported populations (including other Native American groups). Distinctive alleles present in the Lacandon population include HLA-A*24:14 and HLA-B*40:08. Furthermore, in Lacandons we observed a high frequency of haplotypes containing the allele HLA-DRB1*04:11, a relatively frequent allele in comparison with other neighboring indigenous groups. The specific demographic history of the Lacandon population including inbreeding, as well as pathogen selection, may have elevated the frequencies of a small number of HLA class II alleles and DNA blocks. To assess the possible role of different selective pressures in determining Native American HLA diversity, we evaluated the relationship between genetic diversity at HLA-A, HLA-B and HLA-DRB1 and pathogen richness for a global dataset and for Native American populations alone. In keeping with previous studies of such relationships we included distance from Africa as a covariate. After correction for multiple comparisons we did not find any significant relationship between pathogen diversity and HLA genetic diversity (as measured by polymorphism information content) in either our global dataset or the Native American subset of the dataset. We found the expected negative relationship between genetic diversity and distance from Africa in the global dataset, but no relationship between HLA genetic diversity and distance from Africa when Native American populations were considered alone.

Episodic positive diversifying selection on key immune system genes in major avian lineages

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

A Glimpse of the Peptide Profile Presentation by Xenopus laevis MHC Class I: Crystal Structure of pXela-UAA Reveals a Distinct Peptide-Binding Groove

The African clawed frog, Xenopus laevis, is a model species for amphibians. Before metamorphosis, tadpoles do not efficiently express the single classical MHC class I (MHC-I) molecule Xela-UAA, but after metamorphosis, adults express this molecule in abundance. To elucidate the Ag-presenting mechanism of Xela-UAA, in this study, the Xela-UAA structure complex (pXela-UAAg) bound with a peptide from a synthetic random peptide library was determined. The amino acid homology between the Xela-UAA and MHC-I sequences of different species is <45%, and these differences are fully reflected in the three-dimensional structure of pXela-UAAg. Because of polymorphisms and interspecific differences in amino acid sequences, pXela-UAAg forms a distinct peptide-binding groove and presents a unique peptide profile. The most important feature of pXela-UAAg is the two-amino acid insertion in the α2-helical region, which forms a protrusion of ∼3.8 Å that is involved in TCR docking. Comparison of peptide-MHC-I complex (pMHC-I) structures showed that only four amino acids in β2-microglobulin that were bound to MHC-I are conserved in almost all jawed vertebrates, and the most unique feature in nonmammalian pMHC-I molecules is that the AB loop bound β2-microglobulin. Additionally, the binding distance between pMHC-I and CD8 molecules in nonmammals is different from that in mammals. These unique features of pXela-UAAg provide enhanced knowledge of T cell immunity and bridge the knowledge gap regarding the coevolutionary progression of the MHC-I complex from aquatic to terrestrial species.

The Graft-Versus-Leukemia Effect in AML

Allogeneic hematopoietic stem cell transplantation (allo-SCT) is the most established and commonly used cellular immunotherapy in cancer care. It is the most potent anti-leukemic therapy in patients with acute myeloid leukemia (AML) and is routinely used with curative intent in patients with intermediate and poor risk disease. Donor T cells, and possibly other immune cells, eliminate residual leukemia cells after prior (radio)chemotherapy. This immune-mediated response is known as graft-versus-leukemia (GvL). Donor alloimmune responses can also be directed against healthy tissues, which is known as graft-versus-host disease (GvHD). GvHD and GvL often co-occur and, therefore, a major barrier to exploiting the full immunotherapeutic benefit of donor immune cells against patient leukemia is the immunosuppression required to treat GvHD. However, curative responses to allo-SCT and GvHD do not always occur together, suggesting that these two immune responses could be de-coupled in some patients. To make further progress in successfully promoting GvL without GvHD, we must transform our limited understanding of the cellular and molecular basis of GvL and GvHD. Specifically, in most patients we do not understand the antigenic basis of immune responses in GvL and GvHD. Identification of antigens important for GvL but not GvHD, and vice versa, could impact on donor selection, allow us to track GvL immune responses and begin to specifically harness and strengthen anti-leukemic immune responses against patient AML cells, whilst minimizing the toxicity of GvHD.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia-one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty-three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population-for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human-induced changes to the coastal ecosystem it inhabits.

Immunogenetic response of the bananaquit in the face of malarial parasites

BACKGROUND

In the arms race between hosts and parasites, genes involved in the immune response are targets for natural selection. Toll-Like Receptor (TLR) genes play a role in parasite detection as part of the innate immune system whereas Major Histocompatibility Complex (MHC) genes encode proteins that display antigens as part of the vertebrate adaptive immune system. Thus, both gene families are under selection pressure from pathogens. The bananaquit (Coereba flaveola) is a passerine bird that is a common host of avian malarial parasites (Plasmodium sp. and Haemoproteus sp.). We assessed molecular variation of TLR and MHC genes in a wild population of bananaquits and identified allelic associations with resistance/susceptibility to parasitic infection to address hypotheses of avian immune response to haemosporidian parasites.

RESULTS

We found that allele frequencies are associated with infection status at the immune loci studied. A consistent general trend showed the infected groups possessed more alleles at lower frequencies, and exhibited unique alleles, compared to the uninfected group.

CONCLUSIONS

Our results support the theory of natural selection favoring particular alleles for resistance while maintaining overall genetic diversity in the population, a mechanism which has been demonstrated in some systems in MHC previously but understudied in TLRs.

Evolution of major histocompatibility complex gene copy number

MHC genes, which code for proteins responsible for presenting pathogen-derived antigens to the host immune system, show remarkable copy-number variation both between and within species. However, the evolutionary forces driving this variation are poorly understood. Here, we use computer simulations to investigate whether evolution of the number of MHC variants in the genome can be shaped by the number of pathogen species the host population encounters (pathogen richness). Our model assumed that while increasing a range of pathogens recognised, expressing additional MHC variants also incurs costs such as an increased risk of autoimmunity. We found that pathogen richness selected for high MHC copy number only when the costs were low. Furthermore, the shape of the association was modified by the rate of pathogen evolution, with faster pathogen mutation rates selecting for increased host MHC copy number, but only when pathogen richness was low to moderate. Thus, taking into account factors other than pathogen richness may help explain wide variation between vertebrate species in the number of MHC genes. Within population, variation in the number of unique MHC variants carried by individuals (INV) was observed under most parameter combinations, except at low pathogen richness. This variance gave rise to positive correlations between INV and host immunocompetence (proportion of pathogens recognised). However, within-population variation in host immunocompetence declined with pathogen richness. Thus, counterintuitively, pathogens can contribute more to genetic variance for host fitness in species exposed to fewer pathogen species, with consequences to predictions from “Hamilton-Zuk” theory of sexual selection.

Highly polymorphic genes of the Major Histocompatibility Complex (MHC) code for proteins responsible for presenting antigens to lymphocytes, thus initiating adaptive immune response. The polymorphism is driven by coevolution with parasites which are selected to evade recognition by MHC proteins. Expressing many MHC molecules could ensure that an individual could present antigens of most pathogen species encountered, but this comes at a cost, such as enhanced negative selection on lymphocytes leading to holes in T-cell receptor repertoire. Our simulations showed that evolution of the number of MHC genes in the genome is driven by a complex interaction between three factors we explored: pathogen richness, the intrinsic cost of expressing additional MHC variants, and pathogen mutation rate. In contrast to verbal arguments, our results indicate that pathogen richness does not always selects for MHC gene family expansion. Taking into account factors other than pathogen richness, in particular costs of expressing additional MHC variants which are still poorly understood, may help explain striking interspecific variation in the number of MHC genes. Counterintuitively, our results also demonstrated that opportunity for selection on immunocompetence should decrease with MHC gene family expansion.

Mechanisms of Plastic Rescue in Novel Environments

Adaptive phenotypic plasticity provides a mechanism of developmental rescue in novel and rapidly changing environments. Understanding the underlying mechanism of plasticity is important for predicting both the likelihood that a developmental response is adaptive and associated life-history trade-offs that could influence patterns of subsequent evolutionary rescue. Although evolved developmental switches may move organisms toward a new adaptive peak in a novel environment, such mechanisms often result in maladaptive responses. The induction of generalized physiological mechanisms in new environments is relatively more likely to result in adaptive responses to factors such as novel toxins, heat stress, or pathogens. Developmental selection forms of plasticity, which rely on within-individual selective processes, such as shaping of tissue architecture, trial-and-error learning, or acquired immunity, are particularly likely to result in adaptive plasticity in a novel environment. However, both the induction of plastic responses and the ability to be plastic through developmental selection come with significant costs, resulting in delays in reproduction, increased individual investment, and reduced fecundity. Thus, we might expect complex interactions between plastic responses that allow survival in novel environments and subsequent evolutionary responses at the population level.

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.