Major histocompatibility complex class I evolution in songbirds: universal primers, rapid evolution and base compositional shifts in exon 3

Associations among MHC genes, latitude, and avian malaria infections in the rufous-collared sparrow (Zonotrichia capensis)

The major histocompatibility complex (MHC) is a genetic region in jawed vertebrates that contains key genes involved in the immune response. Associations between the MHC and avian malaria infections in wild birds have been observed and mainly explored in the Northern Hemisphere, while a general lack of information remains in the Southern Hemisphere. Here, we investigated the associations between the MHC genes and infections with Plasmodium and Haemoproteus blood parasites along a latitudinal gradient in South America. We sampled 93 rufous-collared sparrows (Zonotrichia capensis) individuals from four countries, Colombia, Ecuador, Peru, and Chile, and estimated MHC-I and MHC-II allele diversity. We detected between 1-4 (MHC-I) and 1-6 (MHC-II) amino acidic alleles per individual, with signs of positive selection. We obtained generalized additive mixed models to explore the associations between MHC-I and MHC-II diversity and latitude. We also explored the relationship between infection status and latitude/biome. We found a non-linear association between the MHC-II amino acidic allele diversity and latitude. Individuals from north Chile presented a lower MHC genetic diversity than those from other locations. We also found an association between deserts and xeric shrublands and a lower prevalence of Haemoproteus parasites. Our results support a lower MHC genetic in arid or semi-arid habitats in the region with the lower prevalence of Haemoproteus parasites.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

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.

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.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

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.

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.

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

Background

The Major Histocompatibility Complex (MHC) codes for the key vertebrate immune receptors responsible for pathogen recognition. Foreign antigens are recognized via their compatibility to hyper-variable region of the peptide-binding groove (PBR), which consists of two separate protein domains. Specifically, the PBR of the MHC class I receptors, which recognize intra-cellular pathogens, has two α domains encoded by exon 2 (α₁) and exon 3 (α₂) of the same gene. Most research on avian MHC class I polymorphism has traditionally focused exclusively on exon 3 and comparisons of selection between the two domains have been hampered by the scarcity of molecular data for exon 2. Thus, it is not clear whether the two domains vary in their specificity towards different antigens and whether they are subject to different selective pressure.

Results

Here, we took advantage of rapidly accumulating genomic resources to test for the differences in selection patterns between both MHC class I domains of the peptide-binding groove in birds. For this purpose, we compiled a dataset of MHC class I exon 2 and 3 sequences for 120 avian species from 46 families. Our phylogenetically-robust approach provided strong evidence for highly consistent levels of selection on the α₁ and α₂ domains. There were strong correlations in all selection measures (number of positively/negatively selected residues and dN/dS ratios) between both PBR exons. Similar positive associations were found for the level of amino acid polymorphism across the two domains.

Conclusions

We conclude that the strength of selection and the level of polymorphism are highly consistent between both peptide-binding domains (α₁ and α₂) of the avian MHC class I.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01812-x.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

Episodic positive diversifying selection on key immune system genes in major avian lineages

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

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.

Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba's short-tailed bat (Carollia perspicillata)

The major histocompatibility complex (MHC) is one of the most diverse genetic regions under pathogen-driven selection because of its central role in antigen binding and immunity. The highest MHC variability, both in terms of the number of individual alleles and gene copies, has so far been found in passerine birds; this is probably attributable to passerine adaptation to both a wide geographic range and a diverse array of habitats. If extraordinary high MHC variation and duplication rates are adaptive features under selection during the evolution of ecologically and taxonomically diverse species, then similarly diverse MHC architectures should be found in bats. Bats are an extremely species-rich mammalian group that is globally widely distributed. Many bat species roost in multitudinous groups and have high contact rates with pathogens, conspecifics, and allospecifics. We have characterized the MHC class I diversity in 116 Panamanian Seba's short-tailed bats (Carollia perspicillata), a widely distributed, generalist, neotropical species. We have detected a remarkable individual and population-level diversity of MHC class I genes, with between seven and 22 alleles and a unique genotype in each individual. This diversity is comparable with that reported in passerine birds and, in both taxonomic groups, further variability has evolved through length polymorphisms. Our findings support the hypothesis that, for species with a geographically broader range, high MHC class I variability is particularly adaptive. Investigation of the details of the underlying adaptive processes and the role of the high MHC diversity in pathogen resistance are important next steps for a better understanding of the role of bats in viral evolution and as carriers of several deadly zoonotic viruses.

Immune genes are hotspots of shared positive selection across birds and mammals

Consistent patterns of positive selection in functionally similar genes can suggest a common selective pressure across a group of species. We use alignments of orthologous protein-coding genes from 39 species of birds to estimate parameters related to positive selection for 11,000 genes conserved across birds. We show that functional pathways related to the immune system, recombination, lipid metabolism, and phototransduction are enriched for positively selected genes. By comparing our results with mammalian data, we find a significant enrichment for positively selected genes shared between taxa, and that these shared selected genes are enriched for viral immune pathways. Using pathogen-challenge transcriptome data, we show that genes up-regulated in response to pathogens are also enriched for positively selected genes. Together, our results suggest that pathogens, particularly viruses, consistently target the same genes across divergent clades, and that these genes are hotspots of host-pathogen conflict over deep evolutionary time.

Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine

Pathogens are one of the main forces driving the evolution and maintenance of the highly polymorphic genes of the vertebrate major histocompatibility complex (MHC). Although MHC proteins are crucial in pathogen recognition, it is still poorly understood how pathogen-mediated selection promotes and maintains MHC diversity, and especially so in host species with highly duplicated MHC genes. Sedge warblers (Acrocephalus schoenobaenus) have highly duplicated MHC genes, and using data from high-throughput MHC genotyping, we were able to investigate to what extent avian malaria parasites explain temporal MHC class I supertype fluctuations in a long-term study population. We investigated infection status and infection intensities of two different strains of Haemoproteus, that is avian malaria parasites that are known to have significant fitness consequences in sedge warblers. We found that prevalence of avian malaria in carriers of specific MHC class I supertypes was a significant predictor of their frequency changes between years. This finding suggests that avian malaria infections partly drive the temporal fluctuations of the MHC class I supertypes. Furthermore, we found that individuals with a large number of different supertypes had higher resistance to avian malaria, but there was no evidence for an optimal MHC class I diversity. Thus, the two studied malaria parasite strains appear to select for a high MHC class I supertype diversity. Such selection may explain the maintenance of the extremely high number of MHC class I gene copies in sedge warblers and possibly also in other passerines where avian malaria is a common disease.

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non-model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non-passerine birds. Non-passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino-acid residues showed marked discrepancies with peptide-binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.

Genotyping strategy matters when analyzing hypervariable major histocompatibility complex-Experience from a passerine bird

Genotyping of classical major histocompatibility complex (MHC) genes is challenging when they are hypervariable and occur in multiple copies. In this study, we used several different approaches to genotype the moderately variable MHC class I exon 3 (MHCIe3) and the highly polymorphic MHC class II exon 2 (MHCIIβe2) in the bluethroat (Luscinia svecica). Two family groups (eight individuals) were sequenced in replicates at both markers using Ion Torrent technology with both a single- and a dual-indexed primer structure. Additionally, MHCIIβe2 was sequenced on Illumina MiSeq. Allele calling was conducted by modifications of the pipeline developed by Sommer et al. (BMC Genomics, 14, 2013, 542) and the software AmpliSAS. While the different genotyping strategies gave largely consistent results for MHCIe3, with a maximum of eight alleles per individual, MHCIIβe2 was remarkably complex with a maximum of 56 MHCIIβe2 alleles called for one individual. Each genotyping strategy detected on average 50%-82% of all MHCIIβe2 alleles per individual, but dropouts were largely allele-specific and consistent within families for each strategy. The discrepancies among approaches indicate PCR biases caused by the platform-specific primer tails. Further, AmpliSAS called fewer alleles than the modified Sommer pipeline. Our results demonstrate that allelic dropout is a significant problem when genotyping the hypervariable MHCIIβe2. As these genotyping errors are largely nonrandom and method-specific, we caution against comparing genotypes across different genotyping strategies. Nevertheless, we conclude that high-throughput approaches provide a major advance in the challenging task of genotyping hypervariable MHC loci, even though they may not reveal the complete allelic repertoire.

A quantitative and qualitative comparison of illumina MiSeq and 454 amplicon sequencing for genotyping the highly polymorphic major histocompatibility complex (MHC) in a non-model species

Background

High-throughput sequencing enables high-resolution genotyping of extremely duplicated genes. 454 amplicon sequencing (454) has become the standard technique for genotyping the major histocompatibility complex (MHC) genes in non-model organisms. However, illumina MiSeq amplicon sequencing (MiSeq), which offers a much higher read depth, is now superseding 454. The aim of this study was to quantitatively and qualitatively evaluate the performance of MiSeq in relation to 454 for genotyping MHC class I alleles using a house sparrow (Passer domesticus) dataset with pedigree information. House sparrows provide a good study system for this comparison as their MHC class I genes have been studied previously and, consequently, we had prior expectations concerning the number of alleles per individual.

Results

We found that 454 and MiSeq performed equally well in genotyping amplicons with low diversity, i.e. amplicons from individuals that had fewer than 6 alleles. Although there was a higher rate of failure in the 454 dataset in resolving amplicons with higher diversity (6–9 alleles), the same genotypes were identified by both 454 and MiSeq in 98% of cases.

Conclusions

We conclude that low diversity amplicons are equally well genotyped using either 454 or MiSeq, but the higher coverage afforded by MiSeq can lead to this approach outperforming 454 in amplicons with higher diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-017-2654-1) contains supplementary material, which is available to authorized users.

Extreme MHC class I diversity in the sedge warbler (Acrocephalus schoenobaenus); selection patterns and allelic divergence suggest that different genes have different functions

Background

Recent work suggests that gene duplications may play an important role in the evolution of immunity genes. Passerine birds, and in particular Sylvioidea warblers, have highly duplicated major histocompatibility complex (MHC) genes, which are key in immunity, compared to other vertebrates. However, reasons for this high MHC gene copy number are yet unclear. High-throughput sequencing (HTS) allows MHC genotyping even in individuals with extremely duplicated genes. This HTS data can reveal evidence of selection, which may help to unravel the putative functions of different gene copies, i.e. neofunctionalization. We performed exhaustive genotyping of MHC class I in a Sylvioidea warbler, the sedge warbler, Acrocephalus schoenobaenus, using the Illumina MiSeq technique on individuals from a wild study population.

Results

The MHC diversity in 863 genotyped individuals by far exceeds that of any other bird species described to date. A single individual could carry up to 65 different alleles, a large proportion of which are expressed (transcribed). The MHC alleles were of three different lengths differing in evidence of selection, diversity and divergence within our study population. Alleles without any deletions and alleles containing a 6 bp deletion showed characteristics of classical MHC genes, with evidence of multiple sites subject to positive selection and high sequence divergence. In contrast, alleles containing a 3 bp deletion had no sites subject to positive selection and had low divergence.

Conclusions

Our results suggest that sedge warbler MHC alleles that either have no deletion, or contain a 6 bp deletion, encode classical antigen presenting MHC molecules. In contrast, MHC alleles containing a 3 bp deletion may encode molecules with a different function. This study demonstrates that highly duplicated MHC genes can be characterised with HTS and that selection patterns can be useful for revealing neofunctionalization. Importantly, our results highlight the need to consider the putative function of different MHC genes in future studies of MHC in relation to disease resistance and fitness.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0997-9) contains supplementary material, which is available to authorized users.

Characterization of MHC class I in a long distance migratory wader, the Icelandic black-tailed godwit

The major histocompatibility complex (MHC) encodes proteins that are central for antigen presentation and pathogen elimination. MHC class I (MHC-I) genes have attracted a great deal of interest among researchers in ecology and evolution and have been partly characterized in a wide range of bird species. So far, the main focus has been on species within the bird orders Galliformes and Passeriformes, while Charadriiformes remain vastly underrepresented with only two species studied to date. These two Charadriiformes species exhibit striking differences in MHC-I characteristics and MHC-I diversity. We therefore set out to study a third species within Charadriiformes, the Icelandic subspecies of black-tailed godwits (Limosa limosa islandica). This subspecies is normally confined to parasite-poor environments, and we hence expected low MHC diversity. MHC-I was partially characterized first using Sanger sequencing and then using high-throughput sequencing (MiSeq) in 84 individuals. We verified 47 nucleotide alleles in open reading frame with classical MHC-I characteristics, and each individual godwit had two to seven putatively classical MHC alleles. However, in contrast to previous MHC-I data within Charadriiformes, we did not find any evidence of alleles with low sequence diversity, believed to represent non-classical MHC genes. The diversity and divergence of the godwits MHC-I genes to a large extent fell between the previous estimates within Charadriiformes. However, the MHC genes of the migratory godwits had few sites subject to positive selection, and one possible explanation could be a low exposure to pathogens.

Very high MHC Class IIB diversity without spatial differentiation in the mediterranean population of greater Flamingos

Background

Selective pressure from pathogens is thought to shape the allelic diversity of major histocompatibility complex (MHC) genes in vertebrates. In particular, both local adaptation to pathogens and gene flow are thought to explain a large part of the intraspecific variation observed in MHC allelic diversity. To date, however, evidence that adaptation to locally prevalent pathogens maintains MHC variation is limited to species with limited dispersal and, hence, reduced gene flow. On the one hand high gene flow can disrupt local adaptation in species with high dispersal rates, on the other hand such species are much more likely to experience spatial variation in pathogen pressure, suggesting that there may be intense pathogen mediated selection pressure operating across breeding sites in panmictic species. Such pathogen mediated selection pressure operating across breeding sites should therefore be sufficient to maintain high MHC diversity in high dispersing species in the absence of local adaptation mechanisms. We used the Greater Flamingo, Phoenicopterus roseus, a long-lived colonial bird showing a homogeneous genetic structure of neutral markers at the scale of the Mediterranean region, to test the prediction that higher MHC allelic diversity with no population structure should occur in large panmictic populations of long-distance dispersing birds than in other resident species.

Results

We assessed the level of allelic diversity at the MHC Class IIB exon 2 from 116 individuals born in four different breeding colonies of Greater Flamingo in the Mediterranean region. We found one of the highest allelic diversity (109 alleles, 2 loci) of any non-passerine avian species investigated so far relative to the number of individuals and loci genotyped. There was no evidence of population structure between the four major Mediterranean breeding colonies.

Conclusion

Our results suggest that local adaptation at MHC Class IIB in Greater Flamingos is constrained by high gene flow and high MHC diversity appears to be maintained by population wide pathogen-mediated selection rather than local pathogen-mediated selection. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in species with large panmictic populations and high dispersal rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0905-3) contains supplementary material, which is available to authorized users.

Coloniality and migration are related to selection on MHC genes in birds

The major histocompatibility complex (MHC) plays a key role in pathogen recognition as a part of the vertebrate adaptive immune system. The great diversity of MHC genes in natural populations is maintained by different forms of balancing selection and its strength should correlate with the diversity of pathogens to which a population is exposed and the rate of exposure. Despite this prediction, little is known about how life-history characteristics affect selection at the MHC. Here, we examined whether the strength of balancing selection on MHC class II genes in birds (as measured with nonsynonymous nucleotide substitutions, dN) was related to their social or migratory behavior, two life-history characteristics correlated with pathogen exposure. Our comparative analysis indicated that the rate of nonsynonymous substitutions was higher in colonial and migratory species than solitary and resident species, suggesting that the strength of balancing selection increases with coloniality and migratory status. These patterns could be attributed to: (1) elevated transmission rates of pathogens in species that breed in dense aggregations, or (2) exposure to a more diverse fauna of pathogens and parasites in migratory species. Our study suggests that differences in social structure and basic ecological traits influence MHC diversity in natural vertebrate populations.

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next-generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans-species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long-term survival of the species.

Egernia stokesii (gidgee skink) MHC I positively selected sites lack concordance with HLA peptide binding regions

Genes of the major histocompatibility complex (MHC) play an important role in vertebrate disease resistance, kin recognition and mate choice. Mammalian MHC is the most widely characterised of all vertebrates, and attention is often given to the peptide binding regions of the MHC because they are presumed to be under stronger selection than non-peptide binding regions. For vertebrates where the MHC is less well understood, researchers commonly use the amino acid positions of the peptide binding regions of the human leukocyte antigen (HLA) to infer the peptide binding regions within the MHC sequences of their taxon of interest. However, positively selected sites within MHC have been reported to lack correspondence with the HLA in fish, frogs, birds and reptiles including squamates. Despite squamate diversity, the MHC has been characterised in few snakes and lizards. The Egernia group of scincid lizards is appropriate for investigating mechanisms generating MHC variation, as their inclusion will add a new lineage (i.e. Scincidae) to studies of selection on the MHC. We aimed to identify positively selected sites within the MHC of Egernia stokesii and then determine if these sites corresponded with the peptide binding regions of the HLA. Six positively selected sites were identified within E. stokesii MHC I, only two were homologous with the HLA. E. stokesii positively selected sites corresponded more closely to non-lizard than other lizard taxa. The characterisation of the MHC of more intermediate taxa within the squamate order is necessary to understand the evolution of the MHC across all vertebrates.

New data from basal Australian songbird lineages show that complex structure of MHC class II β genes has early evolutionary origins within passerines

Background

The major histocompatibility complex (MHC) plays a crucial role in the adaptive immune system and has been extensively studied across vertebrate taxa. Although the function of MHC genes appears to be conserved across taxa, there is great variation in the number and organisation of these genes. Among avian species, for instance, there are notable differences in MHC structure between passerine and non-passerine lineages: passerines typically have a high number of highly polymorphic MHC paralogs whereas non-passerines have fewer loci and lower levels of polymorphism. Although the occurrence of highly polymorphic MHC paralogs in passerines is well documented, their evolutionary origins are relatively unexplored. The majority of studies have focussed on the more derived passerine lineages and there is very little empirical information on the diversity of the MHC in basal passerine lineages. We undertook a study of MHC diversity and evolutionary relationships across seven species from four families (Climacteridae, Maluridae, Pardalotidae, Meliphagidae) that comprise a prominent component of the basal passerine lineages. We aimed to determine if highly polymorphic MHC paralogs have an early evolutionary origin within passerines or are a more derived feature of the infraorder Passerida.

Results

We identified 177 alleles of the MHC class II β exon 2 in seven basal passerine species, with variation in numbers of alleles across individuals and species. Overall, we found evidence of multiple gene loci, pseudoalleles, trans-species polymorphism and high allelic diversity in these basal lineages. Phylogenetic reconstruction of avian lineages based on MHC class II β exon 2 sequences strongly supported the monophyletic grouping of basal and derived passerine species.

Conclusions

Our study provides evidence of a large number of highly polymorphic MHC paralogs in seven basal passerine species, with strong similarities to the MHC described in more derived passerine lineages rather than the simpler MHC in non-passerine lineages. These findings indicate an early evolutionary origin of highly polymorphic MHC paralogs in passerines and shed light on the evolutionary forces shaping the avian MHC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0681-5) contains supplementary material, which is available to authorized users.

Contrasting evolutionary histories of MHC class I and class II loci in grouse--effects of selection and gene conversion

Genes of the major histocompatibility complex (MHC) encode receptor molecules that are responsible for recognition of intracellular and extracellular pathogens (class I and class II genes, respectively) in vertebrates. Given the different roles of class I and II MHC genes, one might expect the strength of selection to differ between these two classes. Different selective pressures may also promote different rates of gene conversion at each class. Despite these predictions, surprisingly few studies have looked at differences between class I and II genes in terms of both selection and gene conversion. Here, we investigated the molecular evolution of MHC class I and II genes in five closely related species of prairie grouse (Centrocercus and Tympanuchus) that possess one class I and two class II loci. We found striking differences in the strength of balancing selection acting on MHC class I versus class II genes. More than half of the putative antigen-binding sites (ABS) of class II were under positive or episodic diversifying selection, compared with only 10% at class I. We also found that gene conversion had a stronger role in shaping the evolution of MHC class II than class I. Overall, the combination of strong positive (balancing) selection and frequent gene conversion has maintained higher diversity of MHC class II than class I in prairie grouse. This is one of the first studies clearly demonstrating that macroevolutionary mechanisms can act differently on genes involved in the immune response against intracellular and extracellular pathogens.

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata)

BACKGROUND

Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system.

RESULTS

Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies.

CONCLUSIONS

Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences.

Characterization of major histocompatibility complex class I loci of the lark sparrow (Chondestes grammacus) and insights into avian MHC evolution

The major histocompatibilty complex (MHC) has become increasingly important in the study of the immunocapabilities of non-model vertebrates due to its direct involvement in the immune response. The characterization of MHC class I loci in the lark sparrow (Chondestes grammacus) revealed multiple MHC class I loci with elevated genetic diversity at exon 3, evidence of differential selection between the peptide binding region (PBR) and non-PBR, and the presence of multiple pseudogenes with limited divergence. The minimum number of functional MHC class I loci was estimated at four. Sequence analysis revealed d N /d S ratios significantly less than one at non-PBR sites, indicative of negative selection, whereas PBR sites associated with antigen recognition showed ratios greater than 1 but non-significant. GenBank surveys and phylogenetic analyses of previously reported avian MHC class I sequences revealed variable signatures of evolutionary processes acting upon this gene family, including gene duplication and potential concerted evolution. An increase in the number of class I loci across species coincided with an increase in pseudogene prevalence, revealing the importance of gene duplication in the expansion of multigene families and the creation of pseudogenes.

Patterns of evolution of MHC class II genes of crows (Corvus) suggest trans-species polymorphism

A distinguishing characteristic of genes that code for the major histocompatibility complex (MHC) is that alleles often share more similarity between, rather than within species. There are two likely mechanisms that can explain this pattern: convergent evolution and trans-species polymorphism (TSP), in which ancient allelic lineages are maintained by balancing selection and retained by descendant species. Distinguishing between these two mechanisms has major implications in how we view adaptation of immune genes. In this study we analyzed exon 2 of the MHC class IIB in three passerine bird species in the genus Corvus: jungle crows (Corvus macrorhynchos japonensis) American crows (C. brachyrhynchos) and carrion crows (C. corone orientalis). Carrion crows and American crows are recently diverged, but allopatric, sister species, whereas carrion crows and jungle crows are more distantly related but sympatric species, and possibly share pathogens linked to MHC IIB polymorphisms. These patterns of evolutionary divergence and current geographic ranges enabled us to test for trans-species polymorphism and convergent evolution of the MHC IIB in crows. Phylogenetic reconstructions of MHC IIB sequences revealed several well supported interspecific clusters containing all three species, and there was no biased clustering of variants among the sympatric carrion crows and jungle crows. The topologies of phylogenetic trees constructed from putatively selected sites were remarkably different than those constructed from putatively neutral sites. In addition, trees constructed using non-synonymous substitutions from a continuous fragment of exon 2 had more, and generally more inclusive, supported interspecific MHC IIB variant clusters than those constructed from the same fragment using synonymous substitutions. These phylogenetic patterns suggest that recombination, especially gene conversion, has partially erased the signal of allelic ancestry in these species. While clustering of positively selected amino acids by supertyping revealed a single supertype shared by only jungle and carrion crows, a pattern consistent with convergence, the overall phylogenetic patterns we observed suggest that TSP, rather than convergence, explains the interspecific allelic similarity of MHC IIB genes in these species of crows.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.

454 screening of individual MHC variation in an endemic island passerine

Genes of the major histocompatibility complex (MHC) code for receptors that are central to the adaptive immune response of vertebrates. These genes are therefore important genetic markers with which to study adaptive genetic variation in the wild. Next-generation sequencing (NGS) has increasingly been used in the last decade to genotype the MHC. However, NGS methods are highly prone to sequencing errors, and although several methodologies have been proposed to deal with this, until recently there have been no standard guidelines for the validation of putative MHC alleles. In this study, we used the 454 NGS platform to screen MHC class I exon 3 variation in a population of the island endemic Berthelot's pipit (Anthus berthelotii). We were able to characterise MHC genotypes across 309 individuals with high levels of repeatability. We were also able to determine alleles that had low amplification efficiencies, whose identification within individuals may thus be less reliable. At the population level we found lower levels of MHC diversity in Berthelot's pipit than in its widespread continental sister species the tawny pipit (Anthus campestris), and observed trans-species polymorphism. Using the sequence data, we identified signatures of gene conversion and evidence of maintenance of functionally divergent alleles in Berthelot's pipit. We also detected positive selection at 10 codons. The present study therefore shows that we have an efficient method for screening individual MHC variation across large datasets in Berthelot's pipit, and provides data that can be used in future studies investigating spatio-temporal patterns and scales of selection on the MHC.

PCR-based isolation of multigene families: lessons from the avian MHC class IIB

The amount of sequence data available today highly facilitates the access to genes from many gene families. Primers amplifying the desired genes over a range of species are readily obtained by aligning conserved gene regions, and laborious gene isolation procedures can often be replaced by quicker PCR-based approaches. However, in the case of multigene families, PCR-based approaches bear the often ignored risk of incomplete isolation of family members. This problem is most prominent in gene families with highly variable and thus unpredictable number of gene copies among species, such as in the major histocompatibility complex (MHC). In this study, we (i) report new primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds and (ii) share our experience with isolating MHCIIB genes from an unprecedented number of avian species from all over the avian phylogeny. We report important and usually underappreciated problems encountered during PCR-based multigene family isolation and provide a collection of measures to help significantly improving the chance of successfully isolating complete multigene families using PCR-based approaches.

Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next-generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide-binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high-frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non-passerine bird.