Computational prediction of MHC II-antigen binding supports divergent allele advantage and explains trans-species polymorphism

HLA gene polymorphism is a modifier of age-related breast cancer penetrance in carriers of BRCA1 pathogenic alleles

PURPOSE

Female carriers of germline BRCA1 mutations almost invariably develop breast cancer (BC); however, the age at onset is a subject of variation. We hypothesized that the age-related penetrance of BRCA1 mutations may depend on inherited variability in the host immune system.

METHODS

Next-generation sequencing was utilized for genotyping of HLA class I/II genes (HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1/3/4/5) in patients with BRCA1-associated BC with early (< / = 38 years, n = 215) and late (> / = 58 years, n = 108) age at onset.

RESULTS

HLA-DQB1*06:03P prevalence was higher in the late-onset group due to the excess of allele carriers [25/108 (23.1%) vs. 22/215 (10.2%); OR 2.96, p < 0.001]. For all HLA-I loci, there was a trend toward an increase in the number of homozygotes in the early-onset group. This trend reached statistical significance for the HLA-A [14.4% vs. 6.5%, p = 0.037; OR 2.4, p = 0.042]. The frequencies of HLA-DPB1, HLA-DQB1, and HLA-DRB1/3/4/5 homozygous genotypes did not differ between young-onset and late-onset patients. The maximum degree of homozygosity detected in this study was 6 out of 7 HLA class I/II loci; all six carriers of these genotypes were diagnosed with BC at the age < / = 38 years [OR 6.97, p = 0.187].

CONCLUSION

HLA polymorphism may play a role in modifying the penetrance of BRCA1 pathogenic variants. Certain HLA alleles or HLA homozygosity may modify the risk of BC in BRCA1 carriers.

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.

Climate heterogeneity, season variation, and sexual dimorphism modulate the association between MHC II diversity and parasite variation in striped hamster

Parasite-mediated selection is widely believed to play a crucial role in maintaining the diversity of the major histocompatibility complex (MHC) genes, which is thought to be maintained through heterozygote advantage, rare-allele advantage, and fluctuating selection. However, the relationship between parasite pressure and MHC diversity has yielded inconsistent findings. These inconsistencies may arise from the influence of environmental factors and individual variations in traits on host-parasite interactions. To address these issues, our study extensively investigated populations of striped hamsters inhabiting regions characterized by environmental heterogeneity. The primary objective was to examine the universality of parasite-mediated selection mechanisms. Our observations revealed the presence of multiple parasite infections, accompanied by spatial and temporal variations in parasite communities and infection patterns among individual hamsters. Specifically, the temperature was found to influence all four parasite indices, while the presence of gamasid mites and parasite richness decreased with increasing precipitation. We also noted significant seasonal variation in parasite dynamics. Moreover, a significant sexual dimorphism was observed with males exhibiting a considerably higher parasite burden compared to their female counterparts. Lastly, we identified the maintenance of MHC polymorphism in striped hamsters as being driven by the heterozygote advantage and fluctuating selection mechanisms. This study underscores the significance of ecological processes in comprehending host-parasite systems and highlights the necessity of considering environmental factors and individual traits when elucidating the mechanisms underlying MHC diversity mediated by parasites.

Allogeneic CD4 T Cells Sustain Effective BK Polyomavirus-Specific CD8 T Cell Response in Kidney Transplant Recipients

Introduction

BK polyomavirus-associated nephropathy (BKPyVAN) is a significant complication in kidney transplant recipients (KTRs), associated with a higher level of plasmatic BK polyomavirus (BKPyV) replication and leading to poor graft survival.

Methods

We prospectively followed-up with 100 KTRs with various degrees of BKPyV reactivation (no BKPyV reactivation, BKPyV-DNAuria, BKPyV-DNAemia, and biopsy-proven BKPyVAN [bp-BKPyVAN], 25 patients per group) and evaluated BKPyV-specific T cell functionality and phenotype.

Results

We demonstrate that bp-BKPyVAN is associated with a loss of BKPyV-specific T cell proliferation, cytokine secretion, and cytotoxic capacities. This severe functional impairment is associated with an overexpression of lymphocyte inhibitory receptors (programmed cell death 1 [PD1], cytotoxic T lymphocyte-associated protein 4, T cell immunoreceptor with Ig and ITIM domains, and T cell immunoglobulin and mucin domain-containing-3), highlighting an exhausted-like phenotype of BKPyV-specific CD4 and CD8 T cells in bp-BKPyVAN. This T cell dysfunction is associated with low class II donor-recipient human leukocyte antigen (HLA) divergence. In contrast, in the context of higher class II donor-recipient HLA (D/R-HLA) divergence, allogeneic CD4 T cells can provide help that sustains BKPyV-specific CD8 T cell responses. In vitro, allogeneic HLA-mismatched CD4 T cells rescue BKPyV-specific CD8 T cell responses.

Conclusion

Our findings suggest that in KTRs, allogeneic CD4 T cells can help to maintain an effective BKPyV-specific CD8 T cell response that better controls BKPyV replication in the kidney allograft and may protect against BKPyVAN.

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.

The structure of songbird MHC class I reveals antigen binding that is flexible at the N-terminus and static at the C-terminus

Long-distance migratory animals such as birds and bats have evolved to withstand selection imposed by pathogens across the globe, and pathogen richness is known to be particularly high in tropical regions. Immune genes, so-called Major Histocompatibility Complex (MHC) genes, are highly duplicated in songbirds compared to other vertebrates, and this high MHC diversity has been hypothesised to result in a unique adaptive immunity. To understand the rationale behind the evolution of the high MHC genetic diversity in songbirds, we determined the structural properties of an MHC class I protein, Acar3, from a long-distance migratory songbird, the great reed warbler Acrocephalus arundinaceus (in short: Acar). The structure of Acar3 was studied in complex with pathogen-derived antigens and shows an overall antigen presentation similar to human MHC class I. However, the peptides bound to Acar3 display an unusual conformation: Whereas the N-terminal ends of the peptides display enhanced flexibility, the conformation of their C-terminal halves is rather static. This uncommon peptide-binding mode in Acar3 is facilitated by a central Arg residue within the peptide-binding groove that fixes the backbone of the peptide at its central position, and potentially permits successful interactions between MHC class I and innate immune receptors. Our study highlights the importance of investigating the immune system of wild animals, such as birds and bats, to uncover unique immune mechanisms which may neither exist in humans nor in model organisms.

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).

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

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.

Assessing the causes and consequences of gut mycobiome variation in a wild population of the Seychelles warbler

BACKGROUND

Considerable research has focussed on the importance of bacterial communities within the vertebrate gut microbiome (GM). However, studies investigating the significance of other microbial kingdoms, such as fungi, are notably lacking, despite their potential to influence host processes. Here, we characterise the fungal GM of individuals living in a natural population of Seychelles warblers (Acrocephalus sechellensis). We evaluate the extent to which fungal GM structure is shaped by environment and host factors, including genome-wide heterozygosity and variation at key immune genes (major histocompatibility complex (MHC) and Toll-like receptor (TLR)). Importantly, we also explore the relationship between fungal GM differences and subsequent host survival. To our knowledge, this is the first time that the genetic drivers and fitness consequences of fungal GM variation have been characterised for a wild vertebrate population.

RESULTS

Environmental factors, including season and territory quality, explain the largest proportion of variance in the fungal GM. In contrast, neither host age, sex, genome-wide heterozygosity, nor TLR3 genotype was associated with fungal GM differences in Seychelles warblers. However, the presence of four MHC-I alleles and one MHC-II allele was associated with changes in fungal GM alpha diversity. Changes in fungal richness ranged from between 1 and 10 sequencing variants lost or gained; in some cases, this accounted for 20% of the fungal variants carried by an individual. In addition to this, overall MHC-I allelic diversity was associated with small, but potentially important, changes in fungal GM composition. This is evidenced by the fact that fungal GM composition differed between individuals that survived or died within 7 months of being sampled.

CONCLUSIONS

Our results suggest that environmental factors play a primary role in shaping the fungal GM, but that components of the host immune system-specifically the MHC-may also contribute to the variation in fungal communities across individuals within wild populations. Furthermore, variation in the fungal GM can be associated with differential survival in the wild. Further work is needed to establish the causality of such relationships and, thus, the extent to which components of the GM may impact host evolution. Video Abstract.

Origin and persistence of polymorphism in loci targeted by disassortative preference: a general model

The emergence and persistence of polymorphism within populations generally requires specific regimes of natural or sexual selection. Here, we develop a unified theoretical framework to explore how polymorphism at targeted loci can be generated and maintained by either disassortative mating choice or balancing selection due to, for example, heterozygote advantage. To this aim, we model the dynamics of alleles at a single locus A in a population of haploid individuals, where reproductive success depends on the combination of alleles carried by the parents at locus A. Our theoretical study of the model confirms that the conditions for the persistence of a given level of allelic polymorphism depend on the relative reproductive advantages among pairs of individuals. Interestingly, equilibria with unbalanced allelic frequencies were shown to emerge from successive introduction of mutants. We then investigate the role of the function linking allelic divergence to reproductive advantage on the evolutionary fate of alleles within the population. Our results highlight the significance of the shape of this function for both the number of alleles maintained and their level of genetic divergence. Large number of alleles are maintained with substantial replacement of alleles, when disassortative advantage slowly increases with allelic differentiation . In contrast, few highly differentiated alleles are predicted to be maintained when genetic differentiation has a strong effect on disassortative advantage. These opposite effects predicted by our model explain how disassortative mate choice may lead to various levels of allelic differentiation and polymorphism, and shed light on the effect of mate preferences on the persistence of balanced and unbalanced polymorphism in natural population.

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.

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.

Class I/Class II HLA Evolutionary Divergence Ratio Is an Independent Marker Associated With Disease-Free and Overall Survival After Allogeneic Hematopoietic Stem Cell Transplantation for Acute Myeloid Leukemia

Class I Human Leukocyte Antigen (HLA) evolutionary divergence (HED) is a metric which reflects immunopeptidome diversity and has been associated with immune checkpoint inhibitor responses in solid tumors. Its impact and interest in allogeneic hematopoietic stem cell transplantation (HCT) have not yet been thoroughly studied. This study analyzed the clinical and immune impact of class I and II HED in 492 acute myeloid leukemia (AML) recipients undergoing HCT. The overall cohort was divided into a training (n=338) and a testing (n=132) set. Univariate cox screening found a positive impact of a high class I HED and a negative impact of a high class II HED on both disease-free (DFS) and overall survival (OS). These results were combined in a unique marker, class I/class II HED ratio, and assessed in the testing cohort. The final multivariate cox model confirmed the positive impact of a high versus low class I/class II HED ratio on both DFS (Hazard Ratio (HR) 0.41 [95% CI 0.2-0.83]; p=0.01) and OS (HR 0.34 [0.19-0.59]; p<0.001), independently of HLA matching and other HCT parameters. No significant association was found between the ratio and graft-versus-host disease (GvHD) nor with neutrophil and platelet recovery. A high class I HED was associated with a tendency for an increase in NK, CD8 T-cell, and B cell recovery at 12 months. These results introduce HED as an original and independent prognosis marker reflecting immunopeptidome diversity and alloreactivity after HCT.

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.

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.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

Genetic variation and selection in the major histocompatibility complex Class II gene in the Guizhou pony

The Guizhou pony (GZP) is an indigenous species of equid found in the mountains of the Guizhou province in southwest China. We selected four regions of the equine leukocyte antigen (ELA), including DQA, DRA, DQB, and DRB, and used them to assess the diversity of the major histocompatibility complex (MHC) class II gene using direct sequencing technology. DRA had the lowest d_N/d_S ratio (0.560) compared with the other three loci, indicating that DRA was conserved and could be conserved after undergoing selective processes. Nine DQA, five DQB, nine DRA, and seven DRB codons were under significant positive selection at the antigen binding sites (ABS), suggesting that the selected residues in ABS may play a significant role in the innate immune system of the GZP. Two GZP alleles were shared with Przewalski’s horse, and six older GZP haplotypes had a better relationship with other horse species by one or two mutational steps, indicating that the GZP may be a natural ancient variety of equid. The specific diversity of ABS and the numbers of unique haplotypes in the evolutionary process affords this species a better genetic fitness and ability to adapt to the native environment.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

HLA Evolutionary Divergence as a Prognostic Marker for AML Patients Undergoing Allogeneic Stem Cell Transplantation

The diversity of human leukocyte antigens (HLAs) enables the presentation of immense repertoires of peptides, including tumor-associated antigens (TAAs). As a surrogate for immunopeptidome diversity, the HLA evolutionary divergence (HED) between individual HLA alleles might directly define the ability to present TAAs, a prerequisite for graft-versus-leukemia effects. We therefore analyzed the impact of HED on survival within a cohort of 171 acute myeloid leukemia (AML) patients after matched donor allogeneic hematopoietic stem cell transplantation (HSCT). Low HED (<25th percentile) of HLA class I (HED^{class I}) or HLA-DR antigens (HED^DR) was a strong determinant for adverse overall survival after allogeneic HSCT (OS), with a hazard ratio for death of 1.9 (95% CI 1.2–3.2) and 2.1 (95% CI 1.3–3.4), respectively. Defining a cutoff value for the combined HED^total (HED^{class I} and HED^DR), the respective 5 year OS was 29.7% and 64.9% in patients with low and high HED^total (p < 0.001), respectively. Furthermore, the risk of relapse was significantly higher in patients with low HED^total (hazard ratio (HR) 2.2, 95% CI 1.3–3.6) and event-free survival (EFS) was significantly reduced (5 year EFS 25.7% versus 54.4%, p < 0.001). We here introduce HED, a fundamental metric of immunopeptidome diversity, as a novel prognostic factor for AML patients undergoing allogeneic HSCT.

Host-parasite interactions in non-native invasive species are dependent on the levels of standing genetic variation at the immune locus

Background

Parasites may mediate the success of biological invasions through their effect on host fitness and thus, on host population growth and stability. However, a release from the pressure of parasites is strongly related to the genetic differentiation of the host. In invasive host populations, the number of available genetic variants, allowing them to ‘fight’ the infection, are likely to be influenced by founder events and genetic drift. The level standing genetic variation of invasive populations may be crucial in successfully adapting to new environments and resisting diseases. We studied invasive populations of raccoon that experienced a random reduction in genetic diversity during the establishment and evaluated the relationship between host immune genetic diversity and intestinal parasites infection.

Results

We distinguished two different genetic clusters that are characterized by different sets of functionally relevant MHC-DRB alleles. Both clusters were characterized by considerably different allele-parasite associations and different levels of parasite infection. The specific resistance MHC-DRB alleles explained the lower prevalence of Digenea parasites. An increased infection intensity was related to the presence of two MHC-DRB alleles. One of these alleles significantly decreased in frequency over time, causing a decrease of Digenea abundance in raccoons in consecutive years.

Conclusions

Our findings suggest that intestinal parasites can exert selective pressure on an invasive host with lowered levels of immune genetic diversity and contribute to promoting local adaptation over time. The random genetic drift that created the two different genetic clusters in the invasive raccoon range imposed completely different MHC-parasite associations, strongly associated with the infection status of populations. Our findings underline the role of standing genetic variation in shaping host-parasite relationships and provide empirical support that functional genetic variation may be, at least partly, responsible for differences in the success of invasive populations.

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.

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.

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.

Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds

Pathogen communities can vary substantially between geographical regions due to different environmental conditions. However, little is known about how host immune systems respond to environmental variation across macro-ecological and evolutionary scales. Here, we select 37 species of songbird that inhabit diverse environments, including African and Palaearctic residents and Afro-Palaearctic migrants, to address how climate and habitat have influenced the evolution of key immune genes, the major histocompatibility complex class I (MHC-I). Resident species living in wetter regions, especially in Africa, had higher MHC-I diversity than species living in drier regions, irrespective of the habitats they occupy. By contrast, no relationship was found between MHC-I diversity and precipitation in migrants. Our results suggest that the immune system of birds has evolved greater pathogen recognition in wetter tropical regions. Furthermore, evolving transcontinental migration appears to have enabled species to escape wet, pathogen-rich areas at key periods of the year, relaxing selection for diversity in immune genes and potentially reducing immune system costs.

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.

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.

Plumage microbiota covaries with the major histocompatibility complex in blue petrels

To increase fitness, a wide range of vertebrates preferentially mate with partners that are dissimilar at the major histocompatibility complex (MHC) or that have high MHC diversity. Although MHC often can be assessed through olfactory cues, the mechanism by which MHC genes influence odour remains largely unclear. MHC class IIB molecules, which enable recognition and elimination of extracellular bacteria, have been suggested to influence odour indirectly by shaping odour-producing microbiota, i.e. bacterial communities. However, there is little evidence of the predicted covariation between an animal's MHC genotype and its bacterial communities in scent-producing body surfaces. Here, using high-throughput sequencing, we tested the covariation between MHC class IIB genotypes and feather microbiota in the blue petrel (Halobaena caerulea), a seabird with highly developed olfaction that has been suggested to rely on oduor cues during an MHC-based mate choice. First, we show that individuals with similar MHC class IIB profiles also have similar bacterial assemblages in their feathers. Then, we show that individuals with high MHC diversity have less diverse feather microbiota and also a reduced abundance of a bacterium of the genus Arsenophonus, a genus in which some species are symbionts of avian ectoparasites. Our results, showing that feather microbiota covary with MHC, are consistent with the hypothesis that individual MHC genotype may shape the semiochemical-producing microbiota in birds.

Cryptic haplotype-specific gamete selection yields offspring with optimal MHC immune genes

Females choose specific mates in order to produce fitter offspring. However, several factors interfere with females' control over fertilization of their eggs, including sneaker males and phenotypically unpredictable allele segregation during meiosis. Mate choice at the individual level thus provides only a poor approximation for obtaining the best genetic match. Consequently, postcopulatory sperm selection by female oocytes has been proposed as a mechanism to achieve complementary combinations of parental haplotypes. Here, using controlled in vitro fertilization of three‐spined stickleback eggs, we find haplotype‐specific fertilization bias toward gametes with complementary major histocompatibility complex (MHC) immunogenes. The resulting zygote (and thus offspring) genotypes exhibit an intermediate level of individual MHC diversity that was previously shown to confer highest pathogen resistance. Our finding of haplotype‐specific gamete selection thus represents an intriguing mechanism for fine‐tuned optimization of the offspring's immune gene composition and an evolutionary advantage in the Red Queen dynamics of host‐parasite coevolution.

Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans

Leprosy, a chronic infectious disease caused by Mycobacterium leprae (M. leprae), was very common in Europe till the 16th century. Here, we perform an ancient DNA study on medieval skeletons from Denmark that show lesions specific for lepromatous leprosy (LL). First, we test the remains for M. leprae DNA to confirm the infection status of the individuals and to assess the bacterial diversity. We assemble 10 complete M. leprae genomes that all differ from each other. Second, we evaluate whether the human leukocyte antigen allele DRB1*15:01, a strong LL susceptibility factor in modern populations, also predisposed medieval Europeans to the disease. The comparison of genotype data from 69 M. leprae DNA-positive LL cases with those from contemporary and medieval controls reveals a statistically significant association in both instances. In addition, we observe that DRB1*15:01 co-occurs with DQB1*06:02 on a haplotype that is a strong risk factor for inflammatory diseases today.

Leprosy, caused by infection with Mycobacterium leprae, was common in Europe in the Middle Ages. Here, Krause-Kyora et al. analyze ancient DNA from a medieval Danish leprosarium to assemble 10 complete bacterial genomes and perform association analysis of the DRB1*15:01 allele with risk of leprosy infection.

Lack of evidence for selection favouring MHC haplotypes that combine high functional diversity

High rates of gene duplication and the highest levels of functional allelic diversity in vertebrate genomes are the main hallmarks of the major histocompatibility complex (MHC), a multigene family with a primordial role in pathogen recognition. The usual tight linkage among MHC gene duplicates may provide an opportunity for the evolution of haplotypes that associate functionally divergent alleles and thus grant the transmission of optimal levels of diversity to coming generations. Even though such associations may be a crucial component of disease resistance, this hypothesis has been given little attention in wild populations. Here, we leveraged pedigree data from a barn owl (Tyto alba) population to characterize MHC haplotype structure across two MHC class I (MHC-I) and two MHC class IIB (MHC-IIB) duplicates, in order to test the hypothesis that haplotypes' genetic diversity is higher than expected from randomly associated alleles. After showing that MHC loci are tightly linked within classes, we found limited evidence for shifts towards MHC haplotypes combining high diversity. Neither amino acid nor functional within-haplotype diversity were significantly higher than in random sets of haplotypes, regardless of MHC class. Our results therefore provide no evidence for selection towards high-diversity MHC haplotypes in barn owls. Rather, high rates of concerted evolution may constrain the evolution of high-diversity haplotypes at MHC-I, while, in contrast, for MHC-IIB, fixed differences among loci may provide barn owls with already optimized functional diversity. This suggests that at the MHC-I and MHC-IIB respectively, different evolutionary dynamics may govern the evolution of within-haplotype diversity.

Evolutionary genetics of immunological supertypes reveals two faces of the Red Queen

Red Queen host-parasite co-evolution can drive adaptations of immune genes by positive selection that erodes genetic variation (Red Queen arms race) or results in a balanced polymorphism (Red Queen dynamics) and long-term preservation of genetic variation (trans-species polymorphism). These two Red Queen processes are opposite extremes of the co-evolutionary spectrum. Here we show that both Red Queen processes can operate simultaneously by analysing the major histocompatibility complex (MHC) in guppies (Poecilia reticulata and P. obscura) and swamp guppies (Micropoecilia picta). Sub-functionalisation of MHC alleles into 'supertypes' explains how polymorphisms persist during rapid host-parasite co-evolution. Simulations show the maintenance of supertypes as balanced polymorphisms, consistent with Red Queen dynamics, whereas alleles within supertypes are subject to positive selection in a Red Queen arms race. Building on the divergent allele advantage hypothesis, we show that functional aspects of allelic diversity help to elucidate the evolution of polymorphic genes involved in Red Queen co-evolution.

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

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

Odour-based discrimination of similarity at the major histocompatibility complex in birds

Many animals are known to preferentially mate with partners that are dissimilar at the major histocompatibility complex (MHC) in order to maximize the antigen binding repertoire (or disease resistance) in their offspring. Although several mammals, fish or lizards use odour cues to assess MHC similarity with potential partners, the ability of birds to assess MHC similarity using olfactory cues has not yet been explored. Here we used a behavioural binary choice test and high-throughput-sequencing of MHC class IIB to determine whether blue petrels can discriminate MHC similarity based on odour cues alone. Blue petrels are seabirds with particularly good sense of smell, they have a reciprocal mate choice and are known to preferentially mate with MHC-dissimilar partners. Incubating males preferentially approached the odour of the more MHC-dissimilar female, whereas incubating females showed opposite preferences. Given their mating pattern, females were, however, expected to show preference for the odour of the more MHC-dissimilar male. Further studies are needed to determine whether, as in women and female mice, the preference varies with the reproductive cycle in blue petrel females. Our results provide the first evidence that birds can use odour cues only to assess MHC dissimilarity.

Spatio-temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi-arid regions. Fish-parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species' adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio-temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio-temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.

Parallel or convergent evolution in human population genomic data revealed by genotype networks

BACKGROUND

Genotype networks are representations of genetic variation data that are complementary to phylogenetic trees. A genotype network is a graph whose nodes are genotypes (DNA sequences) with the same broadly defined phenotype. Two nodes are connected if they differ in some minimal way, e.g., in a single nucleotide.

RESULTS

We analyze human genome variation data from the 1,000 genomes project, and construct haploid genotype (haplotype) networks for 12,235 protein coding genes. The structure of these networks varies widely among genes, indicating different patterns of variation despite a shared evolutionary history. We focus on those genes whose genotype networks show many cycles, which can indicate homoplasy, i.e., parallel or convergent evolution, on the sequence level.

CONCLUSION

For 42 genes, the observed number of cycles is so large that it cannot be explained by either chance homoplasy or recombination. When analyzing possible explanations, we discovered evidence for positive selection in 21 of these genes and, in addition, a potential role for constrained variation and purifying selection. Balancing selection plays at most a small role. The 42 genes with excess cycles are enriched in functions related to immunity and response to pathogens. Genotype networks are representations of genetic variation data that can help understand unusual patterns of genomic variation.