100 million years of multigene family evolution: origin and evolution of the avian MHC class IIB

Macroevolution of avian T cell receptor C segments using genomic data

All jawed vertebrates have four T cell receptor (TCR) chains expressed by thymus-derived lymphocytes that play a significant role in animal immune defense. However, avian TCR studies have been limited to a few species, although their co-functional major histocompatibility complexes (MHCs) have been studied for decades, showing various copy numbers and polymorphisms. Here, using public genome data, we characterized the copy numbers, the phylogenic relationship and selection of T cell receptor complex (TCR-C) segments, and the genomic organization of TCR loci across birds. Various numbers of C segments were found in the TCRα/TCRδ, TCRβ, and TCRγ loci, and phylogenetic analysis reflected both ancient gene duplication events (two Cβ segments and Cδ segments divergent into CδI and CδII) and contemporary evolution (lineage-specific and species-specific characteristics). Most passerines lack CδII segments and a second TRD locus, except Meliphagidae and Maluridae. A relatively stable structure was verified in four TCR loci of birds, except for the arrangement of V segment groups. In this study, we explored the phylogenetic relationships of TCR-C segments across avians for the first time. We inferred gene duplication and loss events during the evolution process. The finding of diverse TCR germline repertoires provides a better understanding of the immune systems of birds.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Major histocompatibility complex genomic investigation of endangered Chinese alligator provides insights into the evolution of tetrapod major histocompatibility complex and survival of critically bottlenecked species

Background

The major histocompatibility complex (MHC) gene family, a vital immune gene family in vertebrates, helps animals defend against pathogens. The polymorphism of MHC genes is important for a species and is considered to be caused by the numerous alleles of MHC antigen-presenting genes. However, the mechanism of this process is unclear due to the lack of data on the MHC structure. The evolutionary trajectories of the tetrapod MHC are also unclear because of insufficient studies on the reptile MHC architecture. Here, we studied the Chinese alligator (Alligator sinensis), which experienced a population bottleneck, but the population increased rapidly in the past 30 years and is proposed to have a unique MHC system to face pathogenic challenges.

Results

We successfully constructed a 2 Mb MHC region using bacterial artificial chromosome (BAC) library and genome data of the Chinese alligator and checked the antigen-presenting genes using transcriptome data and the rapid amplification of cDNA ends (RACE) technique. The MHC architecture reported here uncovers adjacent Class I and Class II subregions and a unique CD1 subregion. This newly added information suggested that the Class I-II structure pattern was more ancient in tetrapods and helped reconstruct the evolution of the MHC region architecture. We also found multiple groups of MHC class I (MHC-I) (12 duplicated loci, belonging to three groups, two of which were novel) and MHC class II (MHC-II) (11 duplicated loci, belonging to two groups) inside the 2 Mb MHC region, and there were three more duplicated MHC-I loci outside it. These highly duplicated antigen-presenting loci had differences in expression, amino acid length of antigen-presenting exons, and splice signal of exon and intron, which together promoted the polymorphism of duplicated genes. Although the MHC antigen-presenting genes were identified as monomorphic or oligomorphic in our previous population study, the loci with high copy numbers and many differences can make up for this loss, presenting another mechanism for polymorphism in antigen presentation. These MHC-I and MHC-IIB loci with low polymorphism for each locus, but high numbers in all, may also contribute to MHC antigen-presenting binding variability in a population.

Conclusion

To summarize, the fine MHC region architecture of reptiles presented in this study completes the evolutionary trajectories of the MHC structure in tetrapods, and these distinctive MHC gene groups in the Chinese alligator may have helped this species to expand rapidly in the past recent years.

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.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

Reconstructing Macroevolutionary Patterns in Avian MHC Architecture With Genomic Data

The Major Histocompatibility Complex (MHC) is a hyper-polymorphic genomic region, which forms a part of the vertebrate adaptive immune system and is crucial for intra- and extra-cellular pathogen recognition (MHC-I and MHC-IIA/B, respectively). Although recent advancements in high-throughput sequencing methods sparked research on the MHC in non-model species, the evolutionary history of MHC gene structure is still poorly understood in birds. Here, to explore macroevolutionary patterns in the avian MHC architecture, we retrieved contigs with antigen-presenting MHC and MHC-related genes from available genomes based on third-generation sequencing. We identified: 1) an ancestral avian MHC architecture with compact size and tight linkage between MHC-I, MHC-IIA/IIB and MHC-related genes; 2) three major patterns of MHC-IIA/IIB unit organization in different avian lineages; and 3) lineage-specific gene translocation events (e.g., separation of the antigen-processing TAP genes from the MHC-I region in passerines), and 4) the presence of a single MHC-IIA gene copy in most taxa, showing evidence of strong purifying selection (low dN/dS ratio and low number of positively selected sites). Our study reveals long-term macroevolutionary patterns in the avian MHC architecture and provides the first evidence of important transitions in the genomic arrangement of the MHC region over the last 100 million years of bird evolution.

Structure, Organization, and Evolution of Satellite DNAs: Insights from the Drosophila repleta and D. virilis Species Groups

The fact that satellite DNAs (satDNAs) in eukaryotes are abundant genomic components, can perform functional roles, but can also change rapidly across species while being homogenous within a species, makes them an intriguing and fascinating genomic component to study. It is also becoming clear that satDNAs represent an important piece in genome architecture and that changes in their structure, organization, and abundance can affect the evolution of genomes and species in many ways. Since the discovery of satDNAs more than 50 years ago, species from the Drosophila genus have continuously been used as models to study several aspects of satDNA biology. These studies have been largely concentrated in D. melanogaster and closely related species from the Sophophora subgenus, even though the vast majority of all Drosophila species belong to the Drosophila subgenus. This chapter highlights some studies on the satDNA structure, organization, and evolution in two species groups from the Drosophila subgenus: the repleta and virilis groups. We also discuss and review the classification of other abundant tandem repeats found in these species in the light of the current information available.

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

First Evidence of a Relationship Between Female Major Histocompatibility Complex Diversity and Eggshell Bacteria in House Sparrows (Passer domesticus)

Bacteria are known to exert positive and negative influences on animals’ health and fitness. Bacteria, in particular those inhabiting the skin and inner organs of vertebrates, are horizontally or vertically transmitted. Specifically, mothers of bird species can transfer bacterial strains to their offspring when the egg is passing the reproductive tract, as the eggshell rubs against the wall of the uterus. In this context, the female immune system might play an important role in influencing the vertical transmission of bacteria. Here, we investigate the relationship between the major histocompatibility complex (MHC) and cultivable eggshell bacteria originating putatively from the female urogenital tract in a captive population of house sparrows (Passer domesticus). We predict that females with a more variable MHC will transfer fewer bacteria onto the eggshells. Our results show a negative relationship between the number of functional MHC class I alleles and bacteria originating in the urinary tract and growing on a selective medium. This is the first study to find a correlation between female MHC diversity and eggshell bacteria.

Long-Read Genome Assemblies Reveal Extraordinary Variation in the Number and Structure of MHC Loci in Birds

Our knowledge of the Major Histocompatibility Complex (MHC) in birds is limited because it often consists of numerous duplicated genes within individuals that are difficult to assemble with short read sequencing technology. Long-read sequencing provides an opportunity to overcome this limitation because it allows the assembly of long regions with repetitive elements. In this study, we used genomes based on long-read sequencing to predict the number and location of MHC loci in a broad range of bird taxa. From the long-read-based genomes of 34 species, we found that there was extremely large variation in the number of MHC loci between species. Overall, there were greater numbers of both class I and II loci in passerines than nonpasserines. The highest numbers of loci (up to 193 class II loci) were found in manakins (Pipridae), which had previously not been studied at the MHC. Our results provide the first direct evidence from passerine genomes of this high level of duplication. We also found different duplication patterns between species. In some species, both MHC class I and II genes were duplicated together, whereas in most species they were duplicated independently. Our study shows that the analysis of long-read-based genomes can dramatically improve our knowledge of MHC structure, although further improvements in chromosome level assembly are needed to understand the evolutionary mechanisms producing the extraordinary interspecific variation in the architecture of the MHC region.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

A Developmental Switch Generating Phenotypic Plasticity Is Part of a Conserved Multi-gene Locus

Switching between alternative complex phenotypes is often regulated by "supergenes," polymorphic clusters of linked genes such as in butterfly mimicry. In contrast, phenotypic plasticity results in alternative complex phenotypes controlled by environmental influences rather than polymorphisms. Here, we show that the developmental switch gene regulating predatory versus non-predatory mouth-form plasticity in the nematode Pristionchus pacificus is part of a multi-gene locus containing two sulfatases and two α-N-acetylglucosaminidases (nag). We provide functional characterization of all four genes, using CRISPR-Cas9-based reverse genetics, and show that nag genes and the previously identified eud-1/sulfatase have opposing influences. Members of the multi-gene locus show non-overlapping neuronal expression and epistatic relationships. The locus architecture is conserved in the entire genus Pristionchus. Interestingly, divergence between paralogs is counteracted by gene conversion, as inferred from phylogenies and genotypes of CRISPR-Cas9-induced mutants. Thus, we found that physical linkage accompanies regulatory linkage between switch genes controlling plasticity in P. pacificus.

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 Copy Number at the MHC Varies across the Avian Tree of Life

The evolution of the major histocompatibility complex (MHC) is shaped by frequent gene duplications and deletions, which generate extensive variation in the number of loci (gene copies) between different taxa. Here, we collected estimates of copy number at the MHC for over 250 bird species from 68 families. We found contrasting patterns of copy number evolution between MHC class I and class IIB, which encode receptors for intra- and extracellular pathogens, respectively. Across the avian evolutionary tree, there was evidence of accelerated evolution and stabilizing selection acting on copy number at class I, while copy number at class IIB was primarily influenced by fluctuating selection and drift. Reconstruction of MHC copy number variation showed ancestrally low numbers of MHC loci in nonpasserines and evolution toward larger numbers of loci in passerines. Different passerine lineages had the highest duplication rates for MHC class I (Sylvioidea) and class IIB (Muscicapoidea and Passeroidea). We also found support for the correlated evolution of MHC copy number and life-history traits such as lifespan and migratory behavior. These results suggest that MHC copy number evolution in birds has been driven by life histories and differences in exposure to intra- and extracellular pathogens.

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

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