Evolution of Copy Number at the MHC Varies across the Avian Tree of Life

Many loci make light work: High individual diversity despite low population diversity and random mating at class I MHC in a Critically Endangered island songbird

Multi-cellular organisms are under constant attack from parasites, making immune defence a critical aspect of fitness. In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the breadth of pathogens to which individuals can respond. Having many MHC alleles can confer better protection against infectious disease, and balancing selection at MHC is widespread. Indeed, MHC loci are famously variable, with some populations harbouring thousands of alleles (Biedrzycka et al., 2018; Robinson, Soormally, Hayhurst, & Marsh, 2016). MHC has also long fascinated behavioural ecologists because mate choice-for example, preferring MHC-dissimilar partners-may amplify the effects of natural selection (Penn & Potts, 1999). But despite keen interest in the evolutionary ecology of MHC, extensive duplication (Minias, Pikus, Whittingham, & Dunn, 2019) has made these genes challenging to study. In a From the Cover article in this issue of Molecular Ecology, Stervander, Dierickx, Thorley, Brooke, and Westerdahl (2020) characterizes class I MHC in a Critically Endangered songbird, relating genotype to mate choice and survivorship. By inferring copy number and patterns of allelic co-segregation, the authors pave the way to elucidating the genomic architecture of MHC in this bottlenecked population. These insights help reconcile apparently counterintuitive findings: no effect of MHC genotype on mate choice or survival, and high MHC diversity within individuals despite low diversity at the population level. The latter finding is cause for optimism regarding conservation prospects. Moreover, these results suggest that ancient duplication events can have longstanding effects on the adaptive landscapes of natural and sexual selection.

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

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

Avian major histocompatibility complex copy number variation is associated with helminth richness

Genes of the major histocompatibility complex (MHC) play a key role in the adaptive immunity of vertebrates, as they encode receptors responsible for antigen recognition. Evolutionary history of the MHC proceeded through numerous gene duplications, which increase the spectrum of pathogens recognized by individuals. Although pathogen-mediated selection is believed to be a primary driver of MHC expansion over evolutionary times, empirical evidence for this association is virtually lacking. Here, we used an extensive dataset on MHC class II copy number variation in non-passerine birds to test for an evolutionary correlation with helminth parasite richness. As expected, our phylogenetically-informed modelling revealed a positive association between MHC copy number and total helminth richness, even after controlling for a broad spectrum of ecological and life-history traits. Thus, total helminth richness appears to be the most important correlate of MHC copy number, supporting a leading role of pathogen-mediated selection in the evolution of MHC in birds. Our results provide some of the first, although correlative, evidence linking parasitism to interspecific variation in MHC copy number among birds.

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC-I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC-I genotype/diversity on mate choice or survival. However, we demonstrate that MHC-I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2-8 linked MHC-I loci. Within-locus homozygosity is high, contributing to low population-wide diversity. Conversely, each individual had comparably many alleles, 6-16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within-individual MHC-I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species' adaptive potential.

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.

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.

Molecular evolution of MHC class II DRB exon 2 in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Raccoon dogs, Nyctereutes procyonoides, are native to East Asia, but have been introduced into western Russia and eastern Europe. To determine allelic diversity and elucidate the evolution of major histocompatibility complex (MHC) class II genes in the raccoon dog, we analysed a 237-bp region of DRB exon 2 from 36 individuals of native and introduced populations from Japan and Russia. We detected 23 DRB alleles (Nypr-DRBs), 22 of which were novel. Some alleles were found across the species’ range, while others were geographically restricted. For both native and introduced populations, the ratio of non-synonymous to synonymous substitution rates for codons at predicted antigen-binding sites was significantly greater than 2, indicating that Nypr-DRBs have evolved under positive selection. Mixed effect model evolution analysis and an algorithm to detect recombination showed five positively selected codons and one recombination breakpoint, respectively. Overall, our results suggest that the diversity of MHC class II DRB in N. procyonoides was influenced and maintained by recombination, pathogen-driven positive selection, geographical barriers and the founder effect. A Bayesian phylogenetic tree revealed no evidence of trans-species polymorphism (TSP), but instead showed monophyly for the Nypr-DRB alleles within a larger clade of canid sequences. The lack of TSP may have been due to long-term divergence of raccoon dogs from other canids, or to their having encountered different sets of pathogens due to occupying a different ecological niche.

Endless forms of sexual selection

In recent years, the field of sexual selection has exploded, with advances in theoretical and empirical research complementing each other in exciting ways. This perspective piece is the product of a "stock-taking" workshop on sexual selection and sexual conflict. Our aim is to identify and deliberate on outstanding questions and to stimulate discussion rather than provide a comprehensive overview of the entire field. These questions are organized into four thematic sections we deem essential to the field. First we focus on the evolution of mate choice and mating systems. Variation in mate quality can generate both competition and choice in the opposite sex, with implications for the evolution of mating systems. Limitations on mate choice may dictate the importance of direct vs. indirect benefits in mating decisions and consequently, mating systems, especially with regard to polyandry. Second, we focus on how sender and receiver mechanisms shape signal design. Mediation of honest signal content likely depends on integration of temporally variable social and physiological costs that are challenging to measure. We view the neuroethology of sensory and cognitive receiver biases as the main key to signal form and the 'aesthetic sense' proposed by Darwin. Since a receiver bias is sufficient to both initiate and drive ornament or armament exaggeration, without a genetically correlated or even coevolving receiver, this may be the appropriate 'null model' of sexual selection. Thirdly, we focus on the genetic architecture of sexually selected traits. Despite advances in modern molecular techniques, the number and identity of genes underlying performance, display and secondary sexual traits remains largely unknown. In-depth investigations into the genetic basis of sexual dimorphism in the context of long-term field studies will reveal constraints and trajectories of sexually selected trait evolution. Finally, we focus on sexual selection and conflict as drivers of speciation. Population divergence and speciation are often influenced by an interplay between sexual and natural selection. The extent to which sexual selection promotes or counteracts population divergence may vary depending on the genetic architecture of traits as well as the covariance between mating competition and local adaptation. Additionally, post-copulatory processes, such as selection against heterospecific sperm, may influence the importance of sexual selection in speciation. We propose that efforts to resolve these four themes can catalyze conceptual progress in the field of sexual selection, and we offer potential avenues of research to advance this progress.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

Divergence between genes but limited allelic polymorphism in two MHC class II A genes in Leach's storm-petrels Oceanodroma leucorhoa

The major histocompatibility complex (MHC) is critical to host-pathogen interactions. Class II MHC is a heterodimer, with α and β subunits encoded by different genes. The peptide-binding groove is formed by the first domain of both subunits (α₁ and β₁), but studies of class II variation or natural selection focus primarily on the β subunit and II B genes. We explored MHC II A in Leach's storm-petrel, a seabird with two expressed, polymorphic II B genes. We found two II A genes, Ocle-DAA and Ocle-DBA, in contrast to the single II A gene in chicken and duck. In exon 2 which encodes the α₁ domain, the storm-petrel II A genes differed strongly from each other but showed little within-gene polymorphism in 30 individuals: just one Ocle-DAA allele, and three Ocle-DBA alleles differing from each other by single non-synonymous substitutions. In a comparable sample, the two II B genes had nine markedly diverged alleles each. Differences between the α₁ domains of Ocle-DAA and Ocle-DBA showed signatures of positive selection, but mainly at non-peptide-binding site (PBS) positions. In contrast, positive selection within and between the II B genes corresponded to putative PBS codons. Phylogenetic analysis of the conserved α₂ domain did not reveal deep or well-supported lineages of II A genes in birds, in contrast to the pronounced differentiation of DQA, DPA, and DRA isotypes in mammals. This uncertain homology complicates efforts to compare levels of functional variation and modes of evolution of II A genes across taxa.

Evolution of major histocompatibility complex gene copy number

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

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

Genetic Diversity and Differentiation at Structurally Varying MHC Haplotypes and Microsatellites in Bottlenecked Populations of Endangered Crested Ibis

Investigating adaptive potential and understanding the relative roles of selection and genetic drift in populations of endangered species are essential in conservation. Major histocompatibility complex (MHC) genes characterized by spectacular polymorphism and fitness association have become valuable adaptive markers. Herein we investigate the variation of all MHC class I and II genes across seven populations of an endangered bird, the crested ibis, of which all current individuals are offspring of only two pairs. We inferred seven multilocus haplotypes from linked alleles in the Core Region and revealed structural variation of the class II region that probably evolved through unequal crossing over. Based on the low polymorphism, structural variation, strong linkage, and extensive shared alleles, we applied the MHC haplotypes in population analysis. The genetic variation and population structure at MHC haplotypes are generally concordant with those expected from microsatellites, underlining the predominant role of genetic drift in shaping MHC variation in the bottlenecked populations. Nonetheless, some populations showed elevated differentiation at MHC, probably due to limited gene flow. The seven populations were significantly differentiated into three groups and some groups exhibited genetic monomorphism, which can be attributed to founder effects. We therefore propose various strategies for future conservation and management.