Ancestral polymorphisms of MHC class II genes: divergent allele advantage

Comparison of human leukocyte antigen in patients with paroxysmal nocturnal hemoglobinuria of different clone sizes

Glycosylphosphatidylinositol-anchored protein-deficient hematopoietic stem and progenitor cell development caused by PIGA mutations cannot fully explain the pathogenesis of paroxysmal nocturnal hemoglobinuria (PNH). Herein, patients newly diagnosed with PNH at our hospital between April 2019 and April 2021 were recruited. The human leukocyte antigen (HLA) class I and II loci were analyzed, and patients were stratified by PNH clone sizes: small (< 50%) and large (≥ 50%). In 40 patients (29 males; 72.5%), the median PNH clone size was 72%. Thirteen (32.5%) and twenty-seven (67.5%) patients harbored small and large PNH clones, respectively. DRB1*15:01 and DQB1*06:02 had higher frequencies in patients with PNH than in healthy controls (adjusted P-value = 4.10 × 10-4 and 4.10 × 10-4, respectively). Whole HLA class I and II allele contributions differed (P = 0.046 and 0.065, not significant difference) when comparing patients with small and large PNH clones. B*13:01 and C*04:01 allelic frequencies were significantly higher in patients with small clones (P = 0.032 and P = 0.032, respectively). Patients with small clones had higher class II HLA evolutionary divergence (HED) (P = 0.041) and global class I and II HED (P = 0.019). In the entire cohort, 17 HLA aberrations were found in 11 (27.5%) patients. No significant differences in HLA aberrations were found between patients with small or large clones. In conclusion, patients with small clones tended to have a higher frequency of immune attack-associated alleles. A higher HED in patients with small clones may reflect a propensity for T cell-mediated autoimmunity. HLA aberrations were similar between patients with small and large clones.

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

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

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Impact of HLA class I functional divergence on HIV control

Heterozygosity of Human leukocyte antigen (HLA) class I genes is linked to beneficial outcomes after HIV infection, presumably through greater breadth of HIV epitope presentation and cytotoxic T cell response. Distinct allotype pairs, however, differ in the extent to which they bind shared sets of peptides. We developed a functional divergence metric that measures pairwise complementarity of allotype-associated peptide binding profiles. Greater functional divergence for pairs of HLA-A and/or HLA-B allotypes was associated with slower AIDS progression and independently with enhanced viral load control. The metric predicts immune breadth at the peptide level rather than gene level and redefines HLA heterozygosity as a continuum differentially affecting disease outcome. Functional divergence may affect response to additional infections, vaccination, immunotherapy, and other diseases where HLA heterozygote advantage occurs.

HLA Genes: A Hallmark of Functional Genetic Variation and Complex Evolution

The major histocompatibility complex (MHC) with its highly polymorphic HLA genes represents one of the most intensely studied genomic regions in the genome. MHC proteins play a key role in antigen-specific immunity and are associated with a wide range of complex diseases. Despite decades of research and many advances in the field, the characterization and interpretation of its genetic and genomic variability remain challenging. Here an overview is provided of the MHC, the nature of its exceptional variability, and the complex evolutionary processes assumed to drive this variability. Highlighted are also recent advances in the field that promise to improve our understanding of the variability in the MHC and in antigen-specific immunity more generally.

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.

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

Correlative CD4 and CD8 T-cell immunodominance in humans and mice: Implications for preclinical testing

Antigen-specific T-cell recognition is restricted by Major Histocompatibility Complex (MHC) molecules, and differences between CD4 and CD8 immunogenicity in humans and animal species used in preclinical vaccine testing are yet to be fully understood. In this study, we addressed this matter by analyzing experimentally identified epitopes based on published data curated in the Immune Epitopes DataBase (IEDB) database. We first analyzed SARS-CoV-2 spike (S) and nucleoprotein (N), which are two common targets of the immune response and well studied in both human and mouse systems. We observed a weak but statistically significant correlation between human and H-2b mouse T-cell responses (CD8 S specific (r = 0.206, p = 1.37 × 10-13); CD4 S specific (r = 0.118, p = 2.63 × 10-5) and N specific (r = 0.179, p = 2.55 × 10-4)). Due to intrinsic differences in MHC molecules across species, we also investigated the association between the immunodominance of common Human Leukocyte Antigen (HLA) alleles for which HLA transgenic mice are available, namely, A*02:01, B*07:02, DRB1*01:01, and DRB1*04:01, and found higher significant correlations for both CD8 and CD4 (maximum r = 0.702, p = 1.36 × 10-31 and r = 0.594, p = 3.04-122, respectively). Our results further indicated that some regions are commonly immunogenic between humans and mice (either H-2b or HLA transgenic) but that others are human specific. Finally, we noted a significant correlation between CD8 and CD4 S- (r = 0.258, p = 7.33 × 1021) and N-specific (r = 0.369, p = 2.43 × 1014) responses, suggesting that discrete protein subregions can be simultaneously recognized by T cells. These findings were confirmed in other viral systems, providing general guidance for the use of murine models to test T-cell immunogenicity of viral antigens destined for human use.

Tutorial: a statistical genetics guide to identifying HLA alleles driving complex disease

The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.

MHC-B Diversity and Signs of Respiratory Illness in Wild, East African Chimpanzees ( Pan troglodytes schweinfurthii )

Many traits, intrinsic and extrinsic to an organism, contribute to interindividual variation in immunity in wild habitats. The vertebrate Major Histocompatibility Complex (MHC) includes genes encoding antigen-presenting molecules that are highly variable, and that variation often predicts susceptibility/resistance to and recovery from pathogen infection. I compare MHC-B variation at two long-term chimpanzee research sites, Kibale National Park in Uganda and Gombe National Park in Tanzania. Using decades of respiratory health data available for these chimpanzees, I test hypotheses associated with maintenance of diversity at MHC loci, including heterozygote, divergent allele, and rare allele advantage hypotheses, and predictions for unique function of MHC-B in great apes. I found, despite confirmation of recent shared ancestry between Kibale and Gombe chimpanzees, including an overlapping MHC-B allele repertoire and similar MHC-B phenotype compositions, chimpanzees from the two research sites experienced differences in the occurrence of respiratory signs and had different associations of MHC-B diversity with signs of respiratory illness. Kibale chimpanzees with heterozygous genotypes and different peptide-binding supertypes were observed less often with respiratory signs than those homozygous or possessing the same supertypes, but this same association was not observed among Gombe chimpanzees. Gombe chimpanzees with specific MHC-B phenotypes that enable engagement of Natural Killer (NK) cells were observed more often with respiratory signs than chimpanzees with other phenotypes, but this was not observed at Kanyawara. This study emphasizes local adaptation in shaping genetic and phenotypic traits in different infectious disease contexts, even among close genetic relatives of the same subspecies, and highlights utility for continued and simultaneous tracking of host immune genes and specific pathogens in wild species.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

Association of HLA diversity with the risk of 25 cancers in the UK Biobank

BACKGROUND

The human leukocyte antigen (HLA) is a highly polymorphic region, and HLA diversity may play a role in presenting tumour-associated peptides and inducing immune responses. However, the effect of HLA diversity on cancers has not been fully assessed. We aimed to explore the role of HLA diversity on cancer development.

METHODS

A pan-cancer analysis was performed to evaluate the effect of HLA diversity, measured by HLA heterozygosity and HLA evolutionary divergence (HED), on the susceptibility of 25 cancers in the UK Biobank.

FINDINGS

We observed that the diversity of HLA class II locus was associated with a lower risk of lung cancer (OR_hetero = 0.94, 95% CI = 0.90-0.97, P = 1.29 × 10^-4) and head and neck cancer (OR_hetero = 0.91, 95% CI = 0.86-0.96, P = 1.56 × 10^-3). Besides, a lower risk of non-Hodgkin lymphoma was associated with an increased diversity of HLA class I (OR_hetero = 0.92, 95% CI = 0.87-0.98, P = 8.38 × 10^-3) and class II locus (OR_hetero = 0.89, 95% CI = 0.86-0.92, P = 1.65 × 10^-10). A lower risk of Hodgkin lymphoma was associated with the HLA class I diversity (OR_hetero = 0.85, 95% CI = 0.75-0.96, P = 0.011). The protective effect of HLA diversity was mainly observed in pathological subtypes with higher tumour mutation burden, such as lung squamous cell carcinoma (P = 9.39 × 10^-3) and diffuse large B cell lymphoma (P_{class I} = 4.12 × 10^-4; P_{class Ⅱ} = 4.71 × 10^-5), as well as the smoking subgroups of lung cancer (P = 7.45 × 10^-5) and head and neck cancer (P = 4.55 × 10^-3).

INTERPRETATION

We provided a systematic insight into the effect of HLA diversity on cancers, which might help to understand the etiological role of HLA on cancer development.

FUNDING

This study was supported by grants from the National Natural Science Foundation of China (82273705, 82003520); the Basic and Applied Basic Research Foundation of Guangdong Province, China (2021B1515420007); the Science and Technology Planning Project of Guangzhou, China (201804020094); Sino-Sweden Joint Research Programme (81861138006); the National Natural Science Foundation of China (81973131, 81903395, 81803319, 81802708).

Large scale phenotypic characterisation of Hierophis viridiflavus (Squamata: Serpentes): climatic and environmental drivers suggest the role of evolutionary processes in a polymorphic species

Colour variability is largely widespread in the animal world as it is tightly associated with fitness and survivorship. Therefore, the drivers and implications of such variability have been of great interest for zoologists in the past decades. Reptiles are excellent models to investigate colour variations and expression under different conditions. Here, we focused on melanism occurrence in the two main lineages of Hierophis viridiflavus at the scale of the species distribution, by extracting available data from iNaturalist, a citizen science network, with the aim of detecting any pure effect of climate or local habitat on colour expression. Our analyses highlighted that habitat does not explain differences in phenotypes, whereas marked effects of geographic and climatic variables were detected. However, the observed climatic effects could be a proxy of the geographical distribution of the two groups, and thus the high occurrence of bright colourations in western populations of the eastern lineage could be addressed to an ongoing event of asymmetric gene flow in contact zones. The current distribution of phenotypes could be the outcome of the evolutionary history of the species combined with the geological history of the Mediterranean region. This investigation, though, is only preliminary and molecular analyses on highly variable regions of the genome are mandatory to address this issue.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

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.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

Novel insights on aquatic mammal MHC evolution: Evidence from manatee DQB diversity

The low diversity in marine mammal major histocompatibility complex (MHC) appears to support the hypothesis of reduced pathogen selective pressure in aquatic systems compared to terrestrial environments. However, the lack of characterization of the aquatic and evolutionarily distant Sirenia precludes drawing more generalized conclusions. Therefore, we aimed to characterize the MHC DQB diversity of two manatee species and compare it with those reported for marine mammals. Our results identified 12 and 6 alleles in T. inunguis and T. manatus, respectively. Alleles show high rates of nonsynonymous substitutions, suggesting loci are evolving under positive selection. Among aquatic mammals, Pinnipeda DQB had smaller numbers of alleles, higher synonymous substitution rate, and a dN/dS ratio closer to 1, suggesting it may be evolving under more relaxed selection compared to fully aquatic mammals. This contradicts one of the predictions of the hypothesis that aquatic environments impose reduced pathogen pressure to mammalian immune system. These results suggest that the unique evolutionary trajectories of mammalian MHC may impose challenges in drawing ecoevolutionary conclusions from comparisons across distant vertebrate lineages.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat

Astroviruses (AstVs) infect numerous mammalian species including reservoirs such as bats. Peptides encoded by the genes of the highly polymorphic Major Histocompatibility Complex (MHC) form the first line of host defence against pathogens. Aside from direct involvement in mounting adaptive immune responses, MHC class II genes are hypothesized to regulate gut commensal diversity and shape the production of immune-modulatory substances by microbes, indirectly affecting host susceptibility. Despite initial empirical evidence for the link between host MHC and the microbiota, associations among these factors remain largely unknown. To fill this gap, we examined MHC allelic diversity and constitution, the gut bacterial community and abundance pattern of a wild population of the neotropical bat (Artibeus jamaicensis) challenged by AstV infections. First, we show an age-dependent relationship between the host MHC class II diversity and constitution and the gut microbiota in AstV uninfected bats. Crucially, these associations changed in AstV infected bats. Additionally, we identify changes in abundance of specific bacterial taxa linked to the presence of certain MHC supertypes and AstV infection. We suggest changes in the microbiota to be either a result of AstV infection or the MHC-mediated modulation of microbial communities. The latter could subsequently affect microbe-mediated immunity and resistance against AstV infection. Our results underscore that the reciprocal nature of host immune genetics, gut microbial diversity and pathogen infection requires attention, which is particularly important given its repercussions for disease susceptibility and severity in wild animal populations with a history of zoonotic spillover and frequent human contact.

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.

Immunogenetic variation shapes the gut microbiome in a natural vertebrate population

BACKGROUND

The gut microbiome (GM) can influence many biological processes in the host, impacting its health and survival, but the GM can also be influenced by the host's traits. In vertebrates, Major Histocompatibility Complex (MHC) genes play a pivotal role in combatting pathogens and are thought to shape the host's GM. Despite this-and the documented importance of both GM and MHC variation to individual fitness-few studies have investigated the association between the GM and MHC in the wild.

RESULTS

We characterised MHC class I (MHC-I), MHC class II (MHC-II) and GM variation in individuals within a natural population of the Seychelles warbler (Acrocephalus sechellensis). We determined how the diversity and composition of the GM varied with MHC characteristics, in addition to environmental factors and other host traits. Our results show that the presence of specific MHC alleles, but not MHC diversity, influences both the diversity and composition of the GM in this population. MHC-I alleles, rather than MHC-II alleles, had the greatest impact on the GM. GM diversity was negatively associated with the presence of three MHC-I alleles (Ase-ua3, Ase-ua4, Ase-ua5), and one MHC-II allele (Ase-dab4), while changes in GM composition were associated with the presence of four different MHC-I alleles (Ase-ua1, Ase-ua7, Ase-ua10, Ase-ua11). There were no associations between GM diversity and TLR3 genotype, but GM diversity was positively correlated with genome-wide heterozygosity and varied with host age and field period.

CONCLUSIONS

These results suggest that components of the host's immune system play a role in shaping the GM of wild animals. Host genotype-specifically MHC-I and to a lesser degree MHC-II variation-can modulate the GM, although whether this occurs directly, or indirectly through effects on host health, is unclear. Importantly, if immune genes can regulate host health through modulation of the microbiome, then it is plausible that the microbiome could also influence selection on immune genes. As such, host-microbiome coevolution may play a role in maintaining functional immunogenetic variation within natural vertebrate populations. Video abstract.

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.

The similarity of class II HLA genotypes defines patterns of autoreactivity in idiopathic bone marrow failure disorders

Idiopathic aplastic anemia (IAA) is a rare autoimmune bone marrow failure (BMF) disorder initiated by a human leukocyte antigen (HLA)-restricted T-cell response to unknown antigens. As in other autoimmune disorders, the predilection for certain HLA profiles seems to represent an etiologic factor; however, the structure-function patterns involved in the self-presentation in this disease remain unclear. Herein, we analyzed the molecular landscape of HLA complexes of a cohort of 300 IAA patients and almost 3000 healthy and disease controls by deeply dissecting their genotypic configurations, functional divergence, self-antigen binding capabilities, and T-cell receptor (TCR) repertoire specificities. Specifically, analysis of the evolutionary divergence of HLA genotypes (HED) showed that IAA patients carried class II HLA molecules whose antigen-binding sites were characterized by a high level of structural homology, only partially explained by specific risk allele profiles. This pattern implies reduced HLA binding capabilities, confirmed by binding analysis of hematopoietic stem cell (HSC)-derived self-peptides. IAA phenotype was associated with the enrichment in a few amino acids at specific positions within the peptide-binding groove of DRB1 molecules, affecting the interface HLA-antigen-TCR β and potentially constituting the basis of T-cell dysfunction and autoreactivity. When analyzing associations with clinical outcomes, low HED was associated with risk of malignant progression and worse survival, underlying reduced tumor surveillance in clearing potential neoantigens derived from mechanisms of clonal hematopoiesis. Our data shed light on the immunogenetic risk associated with IAA etiology and clonal evolution and on general pathophysiological mechanisms potentially involved in other autoimmune disorders.

Evidence of MHC class I and II influencing viral and helminth infection via the microbiome in a non-human primate

Until recently, the study of major histocompability complex (MHC) mediated immunity has focused on the direct link between MHC diversity and susceptibility to parasite infection. However, MHC genes can also influence host health indirectly through the sculpting of the bacterial community that in turn shape immune responses. We investigated the links between MHC class I and II gene diversity gut microbiome diversity and micro- (adenovirus, AdV) and macro- (helminth) parasite infection probabilities in a wild population of non-human primates, mouse lemurs of Madagascar. This setup encompasses a plethora of underlying interactions between parasites, microbes and adaptive immunity in natural populations. Both MHC classes explained shifts in microbiome composition and the effect was driven by a few select microbial taxa. Among them were three taxa (Odoribacter, Campylobacter and Prevotellaceae-UCG-001) which were in turn linked to AdV and helminth infection status, correlative evidence of the indirect effect of the MHC via the microbiome. Our study provides support for the coupled role of MHC diversity and microbial flora as contributing factors of parasite infection.

The selective pressure of the major histocompatibility complex (MHC) on microbial communities, and the potential role of this interaction in driving parasite resistance has been largely neglected. Using a natural population of the primate Microcebus griseorufus, we provide correlative evidence of two outstanding findings: that MHCI and MHCII diversity shapes the composition of the gut microbiota; and that select taxa associated with MHC diversity predicted adenovirus and helminth infection status. Our study highlights the importance of incorporating the microbiome when investigating parasite-mediated MHC selection.

MHC-II distance between parents predicts sex allocation decisions in a genetically monogamous bird

Theory predicts that parental heritable characteristics should shape sex allocation decisions when their effects on reproduction or survival are offspring sex-dependent. Numerous studies have questioned to what extent characteristics displayed by one of the parents matched theoretical expectations. This contrasts with the handful of studies that investigated whether compatibility between parents could also trigger selective pressures for sex allocation adjustments. We studied the genetically monogamous black-legged kittiwake (Rissa tridactyla), where previous data revealed that female chicks suffered higher fitness costs from low diversity at genes of the major histocompatibility complex (MHC) than male chicks. We predicted, and found in our dataset, that MHC-similar parents, producing low MHC-diverse offspring, should avoid the production of females. The relation between MHC-distance between parents (i.e. the functional distinctness of their MHC alleles) and offspring sex was not linear, such that MHC-dissimilar parents also overproduced sons. Overall, our results suggest that the genetically monogamous black-legged kittiwake parents flexibly adapt their reproduction and circumvent the costs of suboptimal pairing by manipulating offspring sex.

Associations between MHC class II variation and phenotypic traits in a free-living sheep population

Pathogen-mediated selection (PMS) is thought to maintain the high level of allelic diversity observed in the major histocompatibility complex (MHC) class II genes. A comprehensive way to demonstrate contemporary selection is to examine associations between MHC variation and individual fitness. As individual fitness is hard to measure, many studies examine associations between MHC variation and phenotypic traits, including direct or indirect measures of adaptive immunity thought to contribute to fitness. Here, we tested associations between MHC class II variation and five phenotypic traits measured in free-living sheep captured in August: weight, strongyle faecal egg count, and plasma IgA, IgE and IgG immunoglobulin titres against the gastrointestinal nematode parasite Teladorsagia circumcincta. We found no association between MHC class II variation and weight or strongyle faecal egg count. We did, however, find associations between MHC class II variation and immunoglobulin levels which varied with isotype, age and sex. Our results suggest associations between MHC and phenotypic traits are more likely to be found for traits more closely associated with pathogen defence than integrative traits such as bodyweight and highlight the association between MHC variation and antibodies in wild populations.

Impact of HLA type, age and chronic viral infection on peripheral T-cell receptor sharing between unrelated individuals

The human adaptive immune system must generate extraordinary diversity to be able to respond to all possible pathogens. The T-cell repertoire derives this high diversity through somatic recombination of the T-cell receptor (TCR) locus, a random process that results in repertoires that are largely private to each individual. However, factors such as thymic selection and T-cell proliferation upon antigen exposure can affect TCR sharing among individuals. By immunosequencing the TCRβ variable region of 426 healthy individuals, we find that, on average, fewer than 1% of TCRβ clones are shared between individuals, consistent with largely private TCRβ repertoires. However, we detect a significant correlation between increased HLA allele sharing and increased number of shared TCRβ clones, with each additional shared HLA allele contributing to an increase in ~0.01% of the total shared TCRβ clones, supporting a key role for HLA type in shaping the immune repertoire. Surprisingly, we find that shared antigen exposure to CMV leads to fewer shared TCRβ clones, even after controlling for HLA, indicative of a largely private response to major viral antigenic exposure. Consistent with this hypothesis, we find that increased age is correlated with decreased overall TCRβ clone sharing, indicating that the pattern of private TCRβ clonal expansion is a general feature of the T-cell response to other infectious antigens as well. However, increased age also correlates with increased sharing among the lowest frequency clones, consistent with decreased repertoire diversity in older individuals. Together, all of these factors contribute to shaping the TCRβ repertoire, and understanding their interplay has important implications for the use of T cells for therapeutics and diagnostics.

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.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele–pathogen interactions. Specific MHC allele–pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Unique pathogen peptidomes facilitate pathogen-specific selection and specialization of MHC alleles

A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes—the pool of short peptides derived from the pathogen’s proteome—so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.

Patterns of MHC-dependent sexual selection in a free-living population of sheep

The MHC is one of the most polymorphic gene clusters in vertebrates and play an essential role in adaptive immunity. Apart from pathogen-mediated selection, sexual selection can also contribute to the maintenance of MHC diversity. MHC-dependent sexual selection could occur via several mechanisms but at present there is no consensus as to which of these mechanisms are involved and their importance. Previous studies have often suffered from limited genetic and behavioural data and small sample size, and were rarely able to examine all the mechanisms together, determine whether signatures of MHC-based non-random mating are independent of genomic effects or differentiate whether MHC-dependent sexual selection takes place at the pre- or post-copulatory stage. In this study, we use Monte Carlo simulation to investigate evidence for non-random MHC-dependent mating patterns by all three mechanisms in a free-living population of Soay sheep. Using 1710 sheep diplotyped at the MHC class IIa region and genome-wide SNPs, together with field observations of consorts, we found sexual selection against a particular haplotype in males at the pre-copulatory stage and sexual selection against female MHC heterozygosity during the rut. We also found MHC-dependent disassortative mating at the post-copulatory stage, along with strong evidence of inbreeding avoidance at both stages. However, results from generalized linear mixed models suggest that the pattern of MHC-dependent disassortative mating could be a by-product of inbreeding avoidance. Our results therefore suggest that while multiple apparent mechanisms of non-random mating with respect to the MHC may occur, some of them have alternative explanations.

The HLA diversity of the Anthony Nolan register

While the success of allogeneic stem cell transplantation depends on a high degree of HLA compatibility between donor and patient, finding a suitable donor remains challenging due to the hyperpolymorphic nature of HLA genes. We calculated high-resolution allele, haplotype and phenotype frequencies for HLA-A, -C, -B, -DRB1 and -DQB1 for 10 subpopulations of the Anthony Nolan (AN) register using an in-house expectation-maximisation (EM) algorithm run on mixed resolution HLA data, covering 676 155 individuals. Sample sizes range from 599 410 for British/Irish North West European (BINWE) individuals, the largest subpopulation in the United Kingdom to 1105 for the British Bangladeshi population. Calculation of genetic distance between the subpopulations based on haplotype frequencies shows three broad clusters, each following a major continental group: European, African and Asian. We further analysed the HLA haplotype and phenotype diversity of each subpopulation, and found that 35.52% of BINWE individuals ranging to 98.34% of Middle Eastern individuals on the register had a unique phenotype within their subpopulation. These analyses and the allele, haplotype and phenotype frequency data of the subpopulation on the AN register are a valuable resource in understanding the HLA diversity in the United Kingdom and can be used to improve the accuracy of match likelihoods and to inform future donor recruitment strategies.

Within-litter covariance of allele-specific MHC heterozygosity, coccidian endoparasite load and growth is modulated by sibling differences in starting mass

Although littermates in altricial mammals usually experience highly similar environmental conditions during early life, considerable differences in growth and health can emerge among them. In a study on subadults of a European rabbit (Oryctolagus cuniculus) population with low MHC polymorphism, we tested whether litter-sibling differences in endoparasitic coccidia load and body mass at the end of the vegetation period were associated with within-litter differences in starting body mass (measured around 2 weeks prior to weaning) and in immune-genetic (MHC class II DRB) constitution. We hypothesized that siblings with a lighter starting mass might be more susceptible to endoparasite infections and thus, negative effects of a more unfavourable MHC constitution might be particularly pronounced in such individuals. Within-litter comparisons revealed that animals with a lighter starting mass reached a relatively lower body mass in autumn. Furthermore, there were indications for an allele-specific heterozygote advantage, as animals with heterozygous combinations of the allele Orcu-DRB*4 had relatively lower hepatic coccidia loads than their littermates with certain homozygous allele combinations. Consistent with our hypothesis, significantly higher hepatic coccidia loads and tendentially lower autumn body masses in homozygous compared to heterozygous individuals for the allele Orcu-DRB*4 were evident in initially lighter but not in heavier siblings, suggesting synergistic effects between an unfavourable MHC constitution and a light starting mass. Taken together, these effects might lead to notable differences in fitness among litter siblings, as a low body mass and a high endoparasite burden are key factors limiting young rabbits' survival during winter.

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.

Sex and hatching order modulate the association between MHC-II diversity and fitness in early-life stages of a wild seabird

Genes of the major histocompatibility complex (MHC) play a pivotal role in parasite resistance, and their allelic diversity has been associated with fitness variations in several taxa. However, studies report inconsistencies in the direction of this association, with either positive, quadratic or no association being described. These discrepancies may arise because the fitness costs and benefits of MHC diversity differ among individuals depending on their exposure and immune responses to parasites. Here, we investigated in black-legged kittiwake (Rissa tridactyla) chicks whether associations between MHC class-II diversity and fitness vary with sex and hatching order. MHC-II diversity was positively associated with growth and tick clearance in female chicks, but not in male chicks. Our data also revealed a positive association between MHC-II diversity and survival in second-hatched female chicks (two eggs being the typical clutch size). These findings may result from condition-dependent parasite infections differentially impacting sexes in relation to hatching order. We thus suggest that it may be important to account for individual heterogeneities in traits that potentially exert selective pressures on MHC diversity in order to properly predict MHC-fitness associations.

Adult survival in migratory caribou is negatively associated with MHC functional diversity

Genes of the major histocompatibility complex (MHC) are involved in acquired immunity in vertebrates. Only a few studies have investigated the fitness consequences of MHC gene diversity in wild populations. Here, we looked at the association between annual survival and body mass and MHC-DRB exon 2 (MHC-DRB) genetic diversity, obtained from high-throughput sequencing, in two declining migratory caribou (Rangifer tarandus) herds. To disentangle the potential direct and general effects of MHC-DRB genetic diversity, we compared different indices of diversity that were either based on DNA-sequence variation or on physicochemical divergence of the translated peptides, thereby covering a gradient of allelic-to-functional diversity. We found that (1) body mass was not related to MHC-DRB diversity or genotype, and (2) adult survival probability was negatively associated with point accepted mutation distance, a corrected distance that considers the likelihood of each amino acid substitution to be accepted by natural selection. In addition, we found no evidence of fluctuating selection over time on MHC-DRB diversity. We concluded that direct effects were involved in the negative relationship between MHC functional diversity and survival, although the mechanism underlying this result remains unclear. A possible explanation could be that individuals with higher MHC diversity suffer higher costs of immunity (immunopathology). Our results suggest that genetic diversity is not always beneficial even in genes that are likely to be strongly shaped by balancing selection.

Comparing raccoon major histocompatibility complex diversity in native and introduced ranges: Evidence for the importance of functional immune diversity for adaptation and survival in novel environments

The adaptive potential of invasive species is related to the genetic diversity of the invader, which is influenced by genetic drift and natural selection. Typically, the genetic diversity of invaders is studied with neutral genetic markers; however, the expectation of reduced diversity has not been consistently supported by empirical studies. Here, we describe and interpret genetic diversity at both neutral microsatellite loci and the immune-related MHC-DRB locus of native and invasive populations of raccoon to better understand of how drift and selection impact patterns of genetic diversity during the invasion process. We found that despite the loss of many MHC (major histocompatibility complex) alleles in comparison with native populations, functional MHC supertypes are preserved in the invasive region. In the native raccoon population, the number of supertypes within individuals was higher than expected under a neutral model. The high level of individual functional divergence may facilitate the adaptation to local conditions in the invasive range. In the invasive populations, we also detected increased population structure at microsatellites compared to the MHC locus, further suggesting that balancing selection is acting on adaptively important regions of the raccoon genome. Finally, we found that alleles known to exhibit resistance to rabies in the native range, Prlo-DRB*4, Prlo-DRB*16 and Prlo-DRB*102, were the most common alleles in the European populations, suggesting directional selection is acting on this locus. Our research shows empirical support for the importance of functional immune diversity for adaptation and survival in novel environments.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.

Divergent Allele Advantage at Human MHC Genes: Signatures of Past and Ongoing Selection

The highly polymorphic genes of the major histocompatibility complex (MHC) play a key role in adaptive immunity. Divergent allele advantage, a mechanism of balancing selection, is proposed to contribute to their exceptional polymorphism. It assumes that MHC genotypes with more divergent alleles allow for broader antigen-presentation to immune effector cells, by that increasing immunocompetence. However, the direct correlation between pairwise sequence divergence and the corresponding repertoire of bound peptides has not been studied systematically across different MHC genes. Here, we investigated this relationship for five key classical human MHC genes (human leukocyte antigen; HLA-A, -B, -C, -DRB1, and -DQB1), using allele-specific computational binding prediction to 118,097 peptides derived from a broad range of human pathogens. For all five human MHC genes, the genetic distance between two alleles of a heterozygous genotype was positively correlated with the total number of peptides bound by these two alleles. In accordance with the major antigen-presentation pathway of MHC class I molecules, HLA-B and HLA-C alleles showed particularly strong correlations for peptides derived from intracellular pathogens. Intriguingly, this bias coincides with distinct protein compositions between intra- and extracellular pathogens, possibly suggesting adaptation of MHC I molecules to present specifically intracellular peptides. Eventually, we observed significant positive correlations between an allele’s average divergence and its population frequency. Overall, our results support the divergent allele advantage as a meaningful quantitative mechanism through which pathogen-mediated selection leads to the evolution of MHC diversity.

Not all birds have a single dominantly expressed MHC-I gene: Transcription suggests that siskins have many highly expressed MHC-I genes

Passerine birds belong to the most species rich bird order and are found in a wide range of habitats. The extremely polymorphic adaptive immune system of passerines, identified through their major histocompatibility complex class I genes (MHC-I), may explain some of this extreme radiation. Recent work has shown that passerines have higher numbers of MHC-I gene copies than other birds, but little is currently known about expression and function of these gene copies. Non-passerine birds have a single highly expressed MHC-I gene copy, a pattern that seems unlikely in passerines. We used high-throughput sequencing to study MHC-I alleles in siskins (Spinus spinus) and determined gene expression, phylogenetic relationships and sequence divergence. We verified between six and 16 MHC-I alleles per individual and 97% of these were expressed. Strikingly, up to five alleles per individual had high expression. Out of 88 alleles 18 were putatively non-classical with low sequence divergence and expression, and found in a single phylogenetic cluster. The remaining 70 alleles were classical, with high sequence divergence and variable degrees of expression. Our results contradict the suggestion that birds only have a single dominantly expressed MHC-I gene by demonstrating several highly expressed MHC-I gene copies in a passerine.

Evolutionary divergence of HLA class I genotype impacts efficacy of cancer immunotherapy

Functional diversity of the highly polymorphic human leukocyte antigen class I (HLA-I) genes underlies successful immunologic control of both infectious disease and cancer. The divergent allele advantage hypothesis dictates that an HLA-I genotype with two alleles with sequences that are more divergent enables presentation of more diverse immunopeptidomes^1-3. However, the effect of sequence divergence between HLA-I alleles-a quantifiable measure of HLA-I evolution-on the efficacy of immune checkpoint inhibitor (ICI) treatment for cancer remains unknown. In the present study the germline HLA-I evolutionary divergence (HED) of patients with cancer treated with ICIs was determined by quantifying the physiochemical sequence divergence between HLA-I alleles of each patient's genotype. HED was a strong determinant of survival after treatment with ICIs. Even among patients fully heterozygous at HLA-I, patients with an HED in the upper quartile respond better to ICIs than patients with a low HED. Furthermore, HED strongly impacts the diversity of tumor, viral and self-immunopeptidomes and intratumoral T cell receptor clonality. Similar to tumor mutation burden, HED is a fundamental metric of diversity at the major histocompatibility complex-peptide complex, which dictates ICI efficacy. The data link divergent HLA allele advantage to immunotherapy efficacy and unveil how ICI response relies on the evolved efficiency of HLA-mediated immunity.

MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach

Proteins encoded by extraordinarily polymorphic major histocompatibility complex (MHC) genes are involved in the adaptive immune response. Balancing selection is believed to maintain MHC polymorphism in the long term, although neutral processes also play a role in shaping MHC diversity. However, the relative contribution of these processes is poorly understood. Here we characterized MHC class II variation of a low-dispersal, pond-breeding newt (Triturus carnifex) over a restricted, geographically structured area. We aimed to (1) evaluate the contribution of selection and neutral processes to shaping MHC diversity at two geographic scales, and (2) test for signatures of divergent allele advantage (DAA), which is a potentially important mechanism of balancing selection. The dominant role of selection in shaping MHC variation was suggested by the lack of correlation between MHC and neutral (microsatellite) variation. Although most variation occurred within populations for both types of markers, they differed in the extent of structuring at the two spatial scales. MHC structuring was more pronounced at local scales, suggesting the role of local selection, while structuring was not detectable at a larger scale, possibly due to the effect of balancing selection. Microsatellites showed the opposite pattern. As expected under DAA, the observed genotypes combined more sequence diversity than expected under a random association of alleles. Thus, DAA may contribute to maintaining MHC polymorphism, which is ancient, as supported by signatures of historical positive selection and trans-species polymorphism. Our results point to the importance of a multi-scale approach in studying MHC variation, especially in low-dispersal taxa, which are genetically structured at fine spatial scales.

Massive introgression of major histocompatibility complex (MHC) genes in newt hybrid zones

Variation in the vertebrate major histocompatibility complex (MHC) genes is crucial for fighting pathogen assault. Because new alleles confer a selective advantage, MHC should readily introgress between species, even under limited hybridization. Using replicated transects through two hybrid zones between strongly reproductively isolated European newts, Lissotriton montandoni and L. vulgaris, we demonstrated recent and ongoing MHC class I and II introgression in the Carpathian region. The extent of introgression correlated with the age of contact. In the older zone, MHC similarity between species within transects exceeded similarity between transects within species, implying pervasive introgression - a massive exchange of MHC genes, not limited to specific variants. In simulations, the observed pattern emerged under the combined action of balancing selection and hybridization, but not when these processes acted separately. Thus, massive introgression at advanced stages of divergence can introduce novel and restore previously lost MHC variation, boosting the adaptive potential of hybridizing taxa. In consequence, MHC genes may be the last to stop introgressing between incipient species.

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.

Characterisation of major histocompatibility complex class IIa haplotypes in an island sheep population

The ovine MHC class IIa is known to consist of six to eight loci located in close proximity on chromosome 20, forming haplotypes that are typically inherited without recombination. Here, we characterise the class IIa haplotypes within the Soay sheep (Ovis aries) on St. Kilda to assess the diversity present within this unmanaged island population. We used a stepwise sequence-based genotyping strategy to identify alleles at seven polymorphic MHC class IIa loci in a sample of 118 Soay sheep from four cohorts spanning 15 years of the long-term study on St. Kilda. DRB1, the most polymorphic MHC class II locus, was characterised first in all 118 sheep and identified six alleles. Using DRB1 homozygous animals, the DQA (DQA1, DQA2 and DQA2-like) and DQB (DQB1, DQB2 and DQB2-like) loci were sequenced, revealing eight haplotypes. Both DQ1/DQ2 and DQ2/DQ2-like haplotype configurations were identified and a single haplotype carrying three DQB alleles. A test sample of 94 further individuals typed at the DRB1 and DQA loci found no exceptions to the eight identified haplotypes and a haplotype homozygosity of 21.3%. We found evidence of historic positive selection at DRB1, DQA and DQB. The limited variation at MHC class IIa loci in Soay sheep enabled haplotype characterisation but showed that no single locus could capture the full extent of the expressed variation in the region.

Selection, drift, and introgression shape MHC polymorphism in lizards

The major histocompatibility complex (MHC) has long served as a model for the evolution of adaptive genetic diversity in wild populations. Pathogen-mediated selection is thought to be a main driver of MHC diversity, but it remains elusive to what degree selection shapes MHC diversity in complex biogeographical scenarios where other evolutionary processes (e.g. genetic drift and introgression) may also be acting. Here we focus on two closely related green lizard species, Lacerta trilineata and L. viridis, to address the evolutionary forces acting on MHC diversity in populations with different biogeographic structure. We characterized MHC class I exon 2 and exon 3, and neutral diversity (microsatellites), to study the relative importance of selection, drift, and introgression in shaping MHC diversity. As expected, positive selection was a significant force shaping the high diversity of MHC genes in both species. Moreover, introgression significantly increased MHC diversity in mainland populations, with a primary direction of gene flow from L. viridis to L. trilineata. Finally, we found significantly fewer MHC alleles in island populations, but maintained MHC sequence and functional diversity, suggesting that positive selection counteracted the effect of drift. Overall, our data support that different evolutionary processes govern MHC diversity in different biogeographical scenarios: positive selection occurs broadly while introgression acts in sympatry and drift when the population sizes decrease.

Deciphering the fine nucleotide diversity of full HLA class I and class II genes in a well-documented population from sub-Saharan Africa

With the aim to understand how next‐generation sequencing (NGS) improves both our assessment of genetic variation within populations and our knowledge on HLA molecular evolution, we sequenced and analysed 8 HLA loci in a well‐documented population from sub‐Saharan Africa (Mandenka). The results of full‐gene NGS‐MiSeq sequencing compared with those obtained by traditional typing techniques or limited sequencing strategies showed that segregating sites located outside exon 2 are crucial to describe not only class I but also class II population diversity. A comprehensive analysis of exons 2, 3, 4 and 5 nucleotide diversity at the 8 HLA loci revealed remarkable differences among these gene regions, notably a greater variation concentrated in the antigen recognition sites of class I exons 3 and some class II exons 2, likely associated with their peptide‐presentation function, a lower diversity of HLA‐C exon 3, possibly related to its role as a KIR ligand, and a peculiar molecular diversity of HLA‐A exon 2, revealing demographic signals. Based on full‐length HLA sequences, we also propose that the most frequent DRB1 allele in the studied population, DRB1*13:04, emerged from an allelic conversion involving 3 potential alleles as donors and DRB1*11:02:01 as recipient. Finally, our analysis revealed a high occurrence of the DRB1*13:04‐DQA1*05:05:01‐DQB1*03:19 haplotype, possibly resulting from a selective sweep due to protection to Onchorcerca volvulus, a prevalent pathogen in West Africa. This study unveils highly relevant information on the molecular evolution of HLA genes in relation to their immune function, calling for similar analyses in other populations living in contrasting environments.