How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

Temporal challenges in detecting balancing selection from population genomic data

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

Correlated evolution of oxidative physiology and MHC-based immunosurveillance in birds

Maintenance and activation of the immune system incur costs, not only in terms of substrates and energy but also via collateral oxidative damage to host cells or tissues during immune response. So far, associations between immune function and oxidative damage have been primarily investigated at intra-specific scales. Here, we hypothesized that pathogen-driven selection should favour the evolution of effective immunosurveillance mechanisms (e.g. major histocompatibility complex, MHC) and antioxidant defences to mitigate oxidative damage resulting from immune function. Using phylogenetically informed comparative approaches, we provided evidence for the correlated evolution of host oxidative physiology and MHC-based immunosurveillance in birds. Species selected for more robust MHC-based immunosurveillance (higher gene copy numbers and allele diversity) showed stronger antioxidant defences, although selection for MHC diversity still showed a positive evolutionary association with oxidative damage to lipids. Our results indicate that historical pathogen-driven selection for highly duplicated and diverse MHC could have promoted the evolution of efficient antioxidant mechanisms, but these evolutionary solutions may be insufficient to keep oxidative stress at bounds. Although the precise nature of mechanistic links between the MHC and oxidative stress remains unclear, our study suggests that a general evolutionary investment in immune function may require co-adaptations at the level of host oxidative metabolism.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

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.

The Narrow Footprint of Ancient Balancing Selection Revealed by Heterokaryon Incompatibility Genes in Aspergillus fumigatus

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2,142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Only rare classical MHC-I alleles are highly expressed in the European house sparrow

The exceptional polymorphism observed within genes of the major histocompatibility complex (MHC), a core component of the vertebrate immune system, has long fascinated biologists. The highly polymorphic classical MHC class-I (MHC-I) genes are maintained by pathogen-mediated balancing selection (PMBS), as shown by many sites subject to positive selection, while the more monomorphic non-classical MHC-I genes show signatures of purifying selection. In line with PMBS, at any point in time, rare classical MHC alleles are more likely than common classical MHC alleles to confer a selective advantage in host-pathogen interactions. Combining genomic and expression data from the blood of wild house sparrows Passer domesticus, we found that only rare classical MHC-I alleles were highly expressed, while common classical MHC-I alleles were lowly expressed or not expressed. Moreover, highly expressed rare classical MHC-I alleles had more positively selected sites, indicating exposure to stronger PMBS, compared with lowly expressed classical alleles. As predicted, the level of expression was unrelated to allele frequency in the monomorphic non-classical MHC-I alleles. Going beyond previous studies, we offer a fine-scale view of selection on classical MHC-I genes in a wild population by revealing differences in the strength of PMBS according to allele frequency and expression level.

Immunogenetic-pathogen networks shrink in Tome's spiny rat, a generalist rodent inhabiting disturbed landscapes

Anthropogenic disturbance may increase the emergence of zoonoses. Especially generalists that cope with disturbance and live in close contact with humans and livestock may become reservoirs of zoonotic pathogens. Yet, whether anthropogenic disturbance modifies host-pathogen co-evolutionary relationships in generalists is unknown. We assessed pathogen diversity, neutral genome-wide diversity (SNPs) and adaptive MHC class II diversity in a rodent generalist inhabiting three lowland rainforest landscapes with varying anthropogenic disturbance, and determined which MHC alleles co-occurred more frequently with 13 gastrointestinal nematodes, blood trypanosomes, and four viruses. Pathogen-specific selection pressures varied between landscapes. Genome-wide diversity declined with the degree of disturbance, while MHC diversity was only reduced in the most disturbed landscape. Furthermore, pristine forest landscapes had more functional important MHC-pathogen associations when compared to disturbed forests. We show co-evolutionary links between host and pathogens impoverished in human-disturbed landscapes. This underscores that parasite-mediated selection might change even in generalist species following human disturbance which in turn may facilitate host switching and the emergence of zoonoses.

Limited associations between MHC diversity and reproductive success in a bird species with biparental care

The selective pressure from pathogens on individuals can have direct consequences on reproduction. Genes from the major histocompatibility complex (MHC) are central to the vertebrate adaptive immune system and pathogen resistance. In species with biparental care, each sex has distinct reproductive roles and levels of investment, and due to a trade-off with immunity, one can expect different selective regimes acting upon the MHC of each parent. Here, we addressed whether couples combine each other's variation at MHC loci to increase their breeding success. Specifically, we used a 23-year dataset from a barn owl population (Tyto alba) to understand how MHC class Iα and IIβ functional divergence and supertypes of each parent were associated with clutch size and fledging success. We did not detect associations between MHC diversity and supertypes with the clutch size or with the fledging success. In addition, to understand the relative contribution from the MHC of the genetic parents and the social parents, we analyzed the fledging success using only a cross-fostered dataset. We found several associations of weak-to-moderate effect sizes between the father's MHC and fledging success: (i) lower MHC-Iα divergence in the genetic father increases fledging success, which might improve paternal care during incubation, and (ii) one and two MHC-IIβ DAB2 supertypes in the social father decrease and increase, respectively, fledging success, which may affect the paternal care after hatching. Furthermore, fledging success increased when both parents did not carry MHC-IIβ DAB1 supertype 2, which could suggest conditional effects of this supertype. Although our study relied on a substantial dataset, we showed that the associations between MHC diversity and reproductive success remain scarce and of complex interpretation in the barn owl. Moreover, our results highlighted the need to incorporate more than one proxy of reproductive success and several MHC classes to capture more complex associations.

Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response

BACKGROUND

The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer.

METHODS

RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles.

RESULTS

The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls.

CONCLUSIONS

The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.

Evolutionary trade-offs constraining the MHC gene expansion: beyond simple TCR depletion model

The immune system is as much shaped by the pressure of pathogens as it is by evolutionary trade-offs that constrain its structure and function. A perfect example comes from the major histocompatibility complex (MHC), molecules that initiate adaptive immune response by presentation of foreign antigens to T cells. The remarkable, population-level polymorphism of MHC genes is assumed to result mainly from a co-evolutionary arms race between hosts and pathogens, while the limited, within-individual number of functional MHC loci is thought to be the consequence of an evolutionary trade-off between enhanced pathogen recognition and excessive T cell depletion during negative selection in the thymus. Certain mathematical models and infection studies suggest that an intermediate individual MHC diversity would thus be optimal. A recent, more direct test of this hypothesis has shown that the effects of MHC diversity on T-cell receptor (TCR) repertoires may differ between MHC classes, supporting the depletion model only for MHC class I. Here, we used the bank vole (Myodes=Cletronomys glareolus), a rodent species with variable numbers of expressed MHC genes, to test how an individual MHC diversity influences the proportions and TCR repertoires of responding T cell subsets. We found a non-linear relationship between MHC diversity and T cell proportions (with intermediate MHC numbers coinciding with the largest T cell proportions), perhaps reflecting an optimality effect of balanced positive and negative thymic selection. The association was strongest for the relationship between MHC class I and splenic CD8+ T cells. The CD8+ TCR richness alone was unaffected by MHC class I diversity, suggesting that MHC class I expansion may be limited by decreasing T cell counts, rather than by direct depletion of TCR richness. In contrast, CD4+ TCR richness was positively correlated with MHC class II diversity, arguing against a universal TCR depletion. It also suggests that different evolutionary forces or trade-offs may limit the within-individual expansion of the MHC class II loci.

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.

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

AmpliSAS and AmpliHLA are tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human specific version; it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Both tools are available in AmpliSAT web-server as well as scripts for local/server installation. Here we describe the installation and deployment of AmpliSAS and AmpliHLA Perl scripts and dependencies on a local or a server computer. We will show how to run them in the command line using as examples four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

Assessing micro-macroparasite selective pressures and anthropogenic disturbance as drivers of immune gene diversity in a Neotropical wild cat

Anthropogenic environmental change is reducing available habitat for wild species, providing novel selection pressures such as infectious diseases and causing species to interact in new ways. The potential for emerging infectious diseases and zoonoses at the interface between humans, domestic animals, and wild species is a key global concern. In vertebrates, diversity at the major histocompatibility complex MHC is critical to disease resilience, and its study in wild populations provides insights into eco-evolutionary dynamics that human activities alter. In natural populations, variation at MHC loci is partly maintained by balancing selection, driven by pathogenic selective pressures. We hypothesize that MHC genetic diversity differs between guigna populations inhabiting human-dominated landscapes (higher pathogen pressures) versus more natural habitats (lower pathogen pressures). We predict that MHC diversity in guignas would be highest in human-dominated landscapes compared with continuous forest habitats. We also expected to find higher MHC diversity in guignas infected with micro and macro parasites (higher parasite load) versus non infected guignas. We characterized for the first time the genetic diversity at three MHC class I and II exons in 128 wild guignas (Leopardus guigna) across their distribution range in Chile (32-46° S) and Argentina, representing landscapes with varying levels of human disturbance. We integrated MHC sequence diversity with multiple measures of anthropogenic disturbance and both micro and macro parasite infection data. We also assessed signatures of positive selection acting on MHC genes. We found significantly higher MHC class I diversity in guignas inhabiting landscapes where houses were present, and with lower percentage of vegetation cover, and also in animals with more severe cardiorespiratory helminth infection (richness and intensity) and micro-macroparasite co-infection. This comprehensive, landscape-level assessment further enhances our knowledge on the evolutionary dynamics and adaptive potential of vertebrates in the face of emerging infectious disease threats and increasing anthropogenic impacts.

The effects of urban life on animal immunity: Adaptations and constraints

Land transformation, including urbanization, is a dominant form of anthropogenic change to the global environment at the dawn of the Anthropocene epoch. More and more species are brought into direct contact with humans, being either required to develop broad-scale adaptations to urban environment or filtered out from urbanized areas. While behavioural or physiological adaptations are at the forefront of urban biology research, there is accumulating evidence for divergent pathogen pressure across urbanization gradients, requiring adjustments in host immune function. At the same time, host immunity may be constrained by unfavourable components of an urban environment, such as poor-quality food resources, disturbance, or pollution. Here, I reviewed existing evidence for adaptations and constrains in the immune system of urban animals, focusing on the recent implementation of metabarcoding, genomic, transcriptomic, and epigenomic approaches in urban biology research. I show that spatial variation in pathogen pressure across urban and non-urban landscapes is highly complex and may be context-dependent, but there is solid evidence for pathogen-driven immunostimulation in urban-dwelling animals. I also show that genes coding for molecules directly involved in interactions with pathogens are the prime candidates for immunogenetic adaptations to urban life. Evidence emerging from landscape genomics and transcriptomics show that immune adaptations to urban life may have a polygenic nature, but immune traits may not be among the key biological functions experiencing broad-scale microevolutionary changes in response to urbanization. Finally, I provided recommendations for future research, including i) a better integration of different 'omic' approaches to obtain a more complete picture of immune adaptations to urban life in non-model animal taxa, ii) quantification of fitness landscapes for immune phenotypes and genotypes across urbanization gradient, and iii) much broader taxonomic coverage (including invertebrates) necessary to draw more robust conclusions on how general (or taxa-specific) are immune responses of animals to urbanization.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

The quest for good genes: Epigamic traits, fitness, MHC and multilocus heterozygosity in the guppy

The 'good genes' hypothesis for the evolution of male secondary sexual traits poses that female preferences for such traits are driven by indirect genetic benefits. However, support for the hypothesis remains ambiguous, and, in particular, the genetic basis for the benefits has rarely been investigated. Here, we use seminatural populations of Trinidadian guppies to investigate whether sexually selected traits (orange, black and iridescent colouration, gonopodium length and body size) predict fitness measured as the number of grandoffspring, a metric that integrates across fitness components and sexes. Furthermore, we tested whether two potential sources of genetic benefits-major histocompatibility complex (MHC) genotypes and multilocus heterozygosity (MLH)-are significant predictors of fitness and of the size of sexually selected traits. We found a significant, nonlinear effect of the area of black pigmentation and male body size on the number of grandoffspring, suggesting stabilizing selection on black area, and nonlinear selection favouring small body size. MLH was heritable (h2  = 0.14) and significantly predicted the number of grandoffspring, indicating the potential for genetic benefits based on heterozygosity. We also found support for local heterozygosity effects, which may reflect a noneven distribution of genetic load across the genome. MHC genotype was not significantly associated with any tested fitness component, or with the load of Gyrodactylus parasites. Neither MHC nor MLH was significant predictor of sexually selected traits. Overall, our results highlight the role of heterozygosity in determining fitness, but do not provide support for male sexually selected traits being indicators of genetic quality.

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Adaptive genetic variations are key for understanding evolutionary processes influencing invasions. However, we have limited knowledge on how adaptive genetic diversity in invasive species responds to new pathogenic environments. Here, we compared variations in immune major histocompatibility complex (MHC) class-II β gene and neutral loci in relation to pathogenic chytrid fungus (Batrachochytrium dendrobatidis, Bd) infection across invasive and native populations of American bullfrog between China and United States (US). Chinese invasive populations show a 60% reduction in neutral cytb variations relative to US native populations, and there were similar MHC variation and functional diversity between them. One MHC allele private to China was under recent positive selection and associated with decreased Bd infection, partly explaining the lower Bd prevalence for Chinese populations than for native US populations. These results suggest that pathogen-mediated selection favors adaptive MHC variations and functional diversity maintenance against serious bottlenecks during the early invasions (within 15 generations) of bullfrogs.

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.

The characterization of toll-like receptor repertoire in Pinna nobilis after mass mortality events suggests adaptive introgression

The fan mussel Pinna nobilis is currently on the brink of extinction due to a multifactorial disease mainly caused to the highly pathogenic parasite Haplosporidium pinnae, meaning that the selection pressure outweighs the adaptive potential of the species. Hopefully, rare individuals have been observed somehow resistant to the parasite, stretching the need to identify the traits underlying this better fitness. Among the candidate to explore at first intention are fast-evolving immune genes, of which toll-like receptor (TLR). In this study, we examined the genetic diversity at 14 TLR loci across P. nobilis, Pinna rudis and P. nobilis × P. rudis hybrid genomes, collected at four physically distant regions, that were found to be either resistant or sensitive to the parasite H. pinnae. We report a high genetic diversity, mainly observed at cell surface TLRs compared with that of endosomal TLRs. However, the endosomal TLR-7 exhibited unexpected level of diversity and haplotype phylogeny. The lack of population structure, associated with a high genetic diversity and elevated dN/dS ratio, was interpreted as balancing selection, though both directional and purifying selection were detected. Interestingly, roughly 40% of the P. nobilis identified as resistant to H. pinnae were introgressed with P. rudis TLR. Specifically, they all carried a TLR-7 of P. rudis origin, whereas sensitive P. nobilis were not introgressed, at least at TLR loci. Small contributions of TLR-6 and TLR-4 single-nucleotide polymorphisms to the clustering of resistant and susceptible individuals could be detected, but their specific role in resistance remains highly speculative. This study provides new information on the diversity of TLR genes within the P. nobilis species after MME and additional insights into adaptation to H. pinnae that should contribute to the conservation of this Mediterranean endemic species.

Neutral Forces and Balancing Selection Interplay to Shape the Major Histocompatibility Complex Spatial Patterns in the Striped Hamster in Inner Mongolia: Suggestive of Broad-Scale Local Adaptation

BACKGROUND

The major histocompatibility complex (MHC) plays a key role in the adaptive immune response to pathogens due to its extraordinary polymorphism. However, the spatial patterns of MHC variation in the striped hamster remain unclear, particularly regarding the relative contribution of the balancing selection in shaping MHC spatial variation and diversity compared to neutral forces.

METHODS

In this study, we investigated the immunogenic variation of the striped hamster in four wild populations in Inner Mongolia which experience a heterogeneous parasitic burden. Our goal was to identify local adaptation by comparing the genetic structure at the MHC with that at seven microsatellite loci, taking into account neutral processes.

RESULTS

We observed significant variation in parasite pressure among sites, with parasite burden showing a correlation with temperature and precipitation. Molecular analysis revealed a similar co-structure between MHC and microsatellite loci. We observed lower genetic differentiation at MHC loci compared to microsatellite loci, and no correlation was found between the two.

CONCLUSIONS

Overall, these results suggest a complex interplay between neutral evolutionary forces and balancing selection in shaping the spatial patterns of MHC variation. Local adaptation was not detected on a small scale but may be applicable on a larger scale.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Immune genotypes, immune responses, and survival in a wild bird population

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

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.

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

Habitat fragmentation matters more than habitat loss: The case of host-parasite interactions

While ecologists agree that habitat loss has a substantial negative effect on biodiversity it is still very much a matter of debate whether habitat fragmentation has a lesser effect and whether this effect is positive or negative for biodiversity. Here, we assess the relative influence of tropical forest loss and fragmentation on the prevalence of vector-borne blood parasites of the genera Plasmodium and Haemoproteus in six forest bird species. We also determine whether habitat loss and fragmentation are associated with a rise or fall in prevalence. We sample more than 4000 individual birds from 58 forest sites in Guadeloupe and Martinique. Considering 34 host-parasite combinations independently and a fine characterization of the amount and spatial configuration of habitat, we use partial least square regressions to disentangle the relative effects of forest loss, forest fragmentation, landscape heterogeneity, and local weather conditions on spatial variability of parasite prevalence. Then we test for the magnitude and the sign of the effect of each environmental descriptor. Strikingly, we show that forest fragmentation explains twice as much of the variance in prevalence as habitat loss or landscape heterogeneity. In addition, habitat fragmentation leads to an overall rise in prevalence in Guadeloupe, but its effect is variable in Martinique. Both habitat loss and landscape heterogeneity exhibit taxon-specific effects. Our results suggest that habitat loss and fragmentation may have contrasting effects between tropical and temperate regions and that inter-specific interactions may not respond in the same way as more commonly used biodiversity metrics such as abundance and diversity.

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

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

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

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.

Partitioning variance in immune traits in a zooplankton host-Fungal parasite system

Host immune traits arise from both genetic and environmental sources of variation. When immune traits have a strong genetic basis, the presence and severity of disease in a population may influence the distribution of those traits. Our study addressed how two immune-related traits (gut penetrability and the hemocyte response) are shaped by genetic and environmental sources of variation, and how the presence of a virulent disease altered the relative frequency of these traits in natural populations. Daphnia dentifera hosts were sampled from five Indiana lakes between June and December 2017 before and during epidemics of their fungal pathogen, Metschnikowia bicuspidata. Collected Daphnia were experimentally exposed to Metschnikowia and assayed for their gut penetrability, hemocyte response, and multi-locus genotype. Mixed-effects models were constructed to partition variance in immune traits between genetic and environmental sources. We then isolated the genetic sources to produce genotype-specific estimates of immune traits for each multi-locus genotype. Finally, we assessed the relative frequency and dynamics of genotypes during epidemics and asked whether genotypes with more robust immune responses increased in frequency during epidemics. Although genotype was an important source of variation for both gut penetrability and the hemocyte response, environmental factors (e.g., resource availability, Metschnikowia prevalence, and co-infection) still explained a large portion of observed variation, suggesting a high degree of flexibility in Daphnia immune traits. Additionally, no significant associations were detected between a genotype's immune traits and its frequency in a population. Our study highlights the power of variance partitioning in understanding the factors driving variation in Daphnia traits and motivates further research on immunological flexibility and the ecological drivers of immune variation.

Complex patterns shape immune genes diversity during invasion of common raccoon in Europe - Selection in action despite genetic drift

Rapid adaptation is common in invasive populations and is crucial to their long-term success. The primary target of selection in the invasive species' new range is standing genetic variation. Therefore, genetic drift and natural selection acting on existing variation are key evolutionary processes through which invaders will evolve over a short timescale. In this study, we used the case of the raccoon Procyon lotor invasion in Europe to identify the forces shaping the diversity of immune genes during invasion. The genes involved in the defence against infection should be under intense selection pressure in the invasive range where novel pathogens are expected to occur. To disentangle the selective and demographic processes shaping the adaptive immune diversity of its invasive and expanding populations, we have developed species-specific single-nucleotide polymorphism markers located in the coding regions of targeted immune-related genes. We characterised the genetic diversity of 110 functionally important immune genes in two invasive and one native raccoon genetic clusters, each presenting a different demographic history. Despite the strong effect of demographic processes in the invasive clusters, we detected a subset of genes exhibiting the diversity pattern suggestive of selection. The most likely process shaping the variation in those genes was balancing selection. The selected genes belong to toll-like receptors and cytokine-related genes. Our results suggest that the prevalence of selection depends on the level of diversity, that is - less genetically diverse invasive population from the Czech Republic displayed fewer signs of selection. Our results highlight the role of standing genetic variation in adapting to new environment. Understanding the evolutionary mechanisms behind invasion success would enable predicting how populations may respond to environmental change.

Heterogeneity in disease resistance and the impact of antibiotics in the US

We hypothesize that the impact of antibiotics is moderated by a population's inherent (genetic) resistance to infectious disease. Using the introduction of sulfa drugs in 1937, we show that US states that are more genetically susceptible to infectious disease saw larger declines in their bacterial mortality rates following the introduction of sulfa drugs in 1937. This suggests area-level genetic endowments of disease resistance and the discovery of medical technologies have acted as substitutes in determining levels of health across the US. We also document immediate effects of sulfa drug exposure to the age of the workforce and cumulative effects on educational attainment for cohorts exposed to sulfa drugs in early life.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

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.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower β diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Genetic diversity of the endangered Mongolian saiga antelope Saiga tatarica mongolica (Artiodactyla: Bovidae) provides insights into conservation

The Saiga antelope (Saiga tatarica) is one of the few megafauna species from the mammoth steppe still living today. Currently, saiga are classified as critically endangered, persisting only in small areas of Central Asian steppe and desert ecosystems. The species is divided into two subspecies: Saiga tatarica mongolica and Saiga tatarica tatarica. In this study, we have for the first time characterized the genetic diversity of the Mongolian saiga (S. t. mongolica) using both mitochondrial DNA and microsatellite markers. We also analysed S. t. tatarica specimens in order to genetically compare both subspecies. The mitochondrial control region was sequenced for a total of 89 individuals: 20 skin, 53 umbilical cord, three placentae and a muscle sample from S. t. mongolica, and a total of 12 hair samples from S. t. tatarica. Additionally, 19 microsatellites developed for saiga antelope were also screened. Our results revealed that the Mongolian saiga presents very low genetic diversity at the mitochondrial level, with no shared mitochondrial haplotype between the two subspecies. Low genetic diversity is also present at the autosomal level, with most loci having low heterozygosity (Ho/He) and a low number of alleles per locus. Despite the low genetic diversity, we found no separation between the subpopulations in Mongolia, indicating that conservation corridors are actually promoting contact between different herds. Our results validate current conservation efforts and inform the implementation of new measures to increase the viability of the S. t. mongolica subspecies.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Selection balancing at innate immune genes: adaptive polymorphism maintenance in Toll-like receptors

Balancing selection is a classic mechanism for maintaining variability in immune genes involved in host–pathogen interactions. However, it remains unclear how widespread the mechanism is across immune genes other than the major histocompatibility complex (MHC). Although occasional reports suggest that balancing selection (heterozygote advantage, negative frequency-dependent selection, and fluctuating selection) may act on other immune genes, the current understanding of the phenomenon in non-MHC immune genes is far from solid. In this review, we focus on Toll-like receptors (TLRs), innate immune genes directly involved in pathogen recognition and immune response activation, as there is a growing body of research testing the assumptions of balancing selection in these genes. After reviewing infection- and fitness-based evidence, along with evidence based on population allelic frequencies and heterozygosity levels, we conclude that balancing selection maintains variation in TLRs, though it tends to occur under specific conditions in certain evolutionary lineages rather than being universal and ubiquitous. Our review also identifies key gaps in current knowledge and proposes promising areas for future research. Improving our understanding of host–pathogen interactions and balancing selection in innate immune genes are increasingly important, particularly regarding threats from emerging zoonotic diseases.

High MHC diversity confers no advantage for phenotypic quality and reproductive performance in a wild bird

1. Genes of the major histocompatibility complex (MHC) encode antigen binding molecules and are an integral part of the acquired immune response of vertebrates. In general, high individual MHC diversity is expected to increase fitness by broadening the spectrum of pathogens recognized by the immune system, in accordance with the heterozygote advantage mechanism. On the other hand, the optimality hypothesis assumes that individuals with optimal (intermediate), rather than maximum diversity of the MHC will achieve the highest fitness because of inherent costs associated with expressing diverse MHC alleles. 2. Here, we tested for associations between individual diversity of the MHC class I and class II genes (binding antigens of intra- and extra-cellular pathogens, respectively) and a range of fitness-related traits (condition, ornament expression and reproduction) in an urban population of the Eurasian coot Fulica atra. 3. Contrary to our expectation, we found that high within-individual allelic diversity of MHC genes (both class I and II) was associated with poorer condition (lower blood haemoglobin concentrations), weaker expression of the putative ornament (smaller frontal shield), later onset of breeding and smaller clutches. An analysis of functional MHC allele clusters (supertypes) provided further support for negative associations of MHC diversity with phenotypic quality and reproductive performance, but most of these relationships could not be explained by the presence of specific maladaptive supertypes. Finally, we found little empirical support for the optimality hypothesis in the Eurasian coot. 4. Our results suggest that the costs of high MHC diversity outweighed any benefits associated with broad MHC repertoire, which could be driven by depauperate pathogen diversity in an urban landscape. To the best of our knowledge, this is one of the first studies providing consistent evidence for negative associations of MHC diversity with a range of fitness-related traits in a natural avian population.