The strength of selection is consistent across both domains of the MHC class I peptide-binding groove in birds

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.

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.

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

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.

Evolution of innate and adaptive immune genes in a non-model waterbird, the common tern

Toll-like receptors (TLRs) and the Major Histocompatibility Complex (MHC) are the key pathogen-recognition genes of vertebrate immune system and they have a crucial role in the initiation of innate and adaptive immune response, respectively. Recent advancements in sequencing technology sparked research on highly duplicated MHC genes in non-model species, but TLR variation in natural vertebrate populations has remained little studied and comparisons of polymorphism across both TLRs and MHC are scarce. Here, we aimed to compare variation across innate (four TLR loci) and adaptive (MHC class I and class II) immune genes in a non-model avian species, the common tern Sterna hirundo. We detected relatively high allelic richness at TLR genes (9-48 alleles per locus), which was similar to or even higher than the estimated per locus allelic richness at the MHC (24-30 alleles at class I and 13-16 alleles at class II under uniform sample sizes). Despite this, the total number of MHC alleles across all duplicated loci (four class I and three class II) was much higher and MHC alleles showed greater sequence divergence than TLRs. Positive selection targeted relatively more sites at the MHC than TLRs, but the strength of selection (dN/dS ratios) at TLRs was higher when compared to MHC class I. There were also differences in the signature of positive selection and recombination (gene conversion) between MHC class I and II (stronger signature at class II), suggesting that mechanisms maintaining variation at the MHC may vary between both classes. Our study indicates that allelic richness of both innate and adaptive immune receptors may be maintained at relatively high levels in viable avian populations and we recommend a transition from the traditional gene-specific to multi-gene approach in studying molecular evolution of vertebrate immune system.

Polymorphism and varying selection within the MHC class I of four Anas species

Ducks (Anatidae) are often vectors for the spread of pathogens because of their long-distance migrations. These migrations also expose ducks to a wide variety of pathogens in their wintering and breeding grounds, and, as a consequence, we might expect strong selection on their immune genes. Here, we studied exons 2 and 3 of the MHC class I in four species of Anas ducks (A. platyrhynchos, A. poecilorhyncha, A. formosa, and A. querquedula) using Illumina-sequencing. Both exons 2 and 3 code for the peptide-binding region of class I molecules; however, most previous studies of birds have only focused on exon 3. Here, we found stronger positive selection on exon 2 than exon 3, as indicated by more species with d_N/d_S > 1 and higher Wu-Kabat values. There was little evidence that divergence time influenced polymorphism, the numbers of identical alleles (partial α1 or α2 regions) among four Anas, or selection, suggesting that these widespread species might share similar levels of selection from pathogens. The high similarity of allele numbers, positively selected sites (PSS), conserved motifs, and variable protein sites (VPS) supported the persistence of trans-species polymorphism in Anas for at least 10 million years. Our study revealed exon 2 as a relatively unexplored source of variation in avian MHC class I, which should be considered in future studies.