Allelic diversity and selection at the MHC class I and class II in a bottlenecked bird of prey, the White-tailed Eagle

Partial molecular characterization, expression pattern and polymorphism analysis of MHC I genes in Chinese domestic goose (Anser cygnoides)

Major histocompatibility complex (MHC) allelic polymorphism is critically important for mediating antigen presentation in vertebrates. Presently, there are insufficient studies of MHC genetic diversity in domestic Anseriform birds. In this study, we analyzed the expression profile of MHC I genes and screened for MHC I exon 2 polymorphism in one domestic goose population from China using Illumina MiSeq sequencing. The results showed that four MHC I alleles (Ancy-IE2*09/*11/*13/*21) in one goose were identified based on cDNA cloning and sequencing using four primer combinations, and the varying number of cDNA clones implied that these four classical sequences showed differential expression patterns. Through next-generation sequencing, 27 alleles were obtained from 68 geese with 3-10 putative alleles per individual, indicating at least the existence of 5 MHC I loci in the goose. The marked excess of the non-synonymous over the synonymous substitution in the peptide-binding region (PBR) along 27 alleles and five positively selected sites (PSSs) detected around the PBR indicated that balancing selection might be the major force in shaping high MHC variation in the goose. Additionally, IA alleles displaying lower polymorphism were subject to less positive selection pressure than non-IA alleles with a higher level of polymorphism.

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.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The F_ST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

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

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

Genetic Diversity and Differentiation of MHC Class I Genes in Red-Crowned Crane Populations

The red-crowned crane (Grus japonensis) has been demoted to “vulnerable species” because its populations have apparently stabilized in Japan and Korea. Low variation and genetic drift may cause damage to the nascent recovery of the G. japonensis population. The major histocompatibility complex (MHC) is one of the most polymorphic gene families in the vertebrate genome and can reflect information on the adaptive evolution of endangered species. In this study, variations in MHC I exon 3 of captive G. japonensis in China were assessed and compared with those in cranes from Japan. Forty MHC alleles of 274 base pairs were isolated from 32 individuals from two captive populations in China. There was high variability in the nucleotide and amino acid composition, showing the proportion of polymorphic sites of 18.98 and 32.97%, respectively. Comparative analyses of the Chinese and Japanese populations based on 222 base pair sequences revealed more alleles and higher variation in the Chinese population. The lack of significant geographical differentiation of G. japonensis was supported by the genetic differentiation coefficient (0.04506) between the Chinese and Japanese populations. Positive selection of antigen-binding sites was observed, which contributed to maintaining the diversity of MHC class I genes. Phylogenetic analysis suggested the persistence of trans-species polymorphisms among MHC class I genes in Gruidae species. Our results may contribute to optimizing the management of G. japonensis populations and population recovery of this threatened species.

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.

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.

MHC Reflects Fine-Scale Habitat Structure in White-Tailed Eagles, Haliaeetus albicilla

The major histocompatibility complex (MHC) genes code for key immune receptors responsible for recognition of intra- and extracellular pathogens (MHC class I and class II, respectively). It was hypothesized that MHC polymorphism can be maintained via fluctuating selection resulting from between-habitat variation in pathogen regimes. We examined associations between MHC class I and class II genes and habitat structure in an apex avian predator, the white-tailed eagle, Haliaeetus albicilla. We genotyped MHC class I and class II genes in ca. 150 white-tailed eagle chicks from nearly 100 nesting territories distributed across 3 distinct populations in Poland. Habitat structure was quantified at the level of foraging territories and directly at the nest sites. We found strong support for associations of habitat traits with diversity and allelic composition at the MHC class II. Forest area within territory and forest productivity were identified as the major habitat predictors of MHC class II polymorphism, whereas other habitat traits (distance to nearest open water, grassland, and water area within territory or understory presence) showed fewer associations with class II alleles. In contrast, there was little support for associations between MHC class I genes and habitat structure. All significant associations were apparent at the within-population level rather than between populations. Our results suggest that extracellular (rather than intracellular) pathogens may exert much stronger selective pressure on the white-tailed eagle. Associations of habitat structure with MHC class II may reflect fluctuating (balancing) selection, which maintains MHC diversity within populations.

Allelic diversity and patterns of selection at the major histocompatibility complex class I and II loci in a threatened shorebird, the Snowy Plover (Charadrius nivosus)

Background

Understanding the structure and variability of adaptive loci such as the major histocompatibility complex (MHC) genes is a primary research goal for evolutionary and conservation genetics. Typically, classical MHC genes show high polymorphism and are under strong balancing selection, as their products trigger the adaptive immune response in vertebrates. Here, we assess the allelic diversity and patterns of selection for MHC class I and class II loci in a threatened shorebird with highly flexible mating and parental care behaviour, the Snowy Plover (Charadrius nivosus) across its broad geographic range.

Results

We determined the allelic and nucleotide diversity for MHC class I and class II genes using samples of 250 individuals from eight breeding population of Snowy Plovers. We found 40 alleles at MHC class I and six alleles at MHC class II, with individuals carrying two to seven different alleles (mean 3.70) at MHC class I and up to two alleles (mean 1.45) at MHC class II. Diversity was higher in the peptide-binding region, which suggests balancing selection. The MHC class I locus showed stronger signatures of both positive and negative selection than the MHC class II locus. Most alleles were present in more than one population. If present, private alleles generally occurred at very low frequencies in each population, except for the private alleles of MHC class I in one island population (Puerto Rico, lineage tenuirostris).

Conclusion

Snowy Plovers exhibited an intermediate level of diversity at the MHC, similar to that reported in other Charadriiformes. The differences found in the patterns of selection between the class I and II loci are consistent with the hypothesis that different mechanisms shape the sequence evolution of MHC class I and class II genes. The rarity of private alleles across populations is consistent with high natal and breeding dispersal and the low genetic structure previously observed at neutral genetic markers in this species.

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

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

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

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

Evolutionary relationship and population structure of domestic Bovidae animals based on MHC-linked and neutral autosomal microsatellite markers

Major histocompatibility complex (MHC) genes are critical for disease resistance or susceptibility responsible for host-pathogen interactions determined mainly by extensive polymorphisms in the MHC genes. Here, we examined the diversity and phylogenetic pattern of MHC haplotypes reconstructed using three MHC-linked microsatellite markers in 55 populations of five Bovidae species and compared them with those based on neutral autosomal microsatellite markers (NAMs). Three-hundred-and-forty MHC haplotypes were identified in 1453 Bovidae individuals, suggesting significantly higher polymorphism and heterozygosity compared with those based on NAMs. The ambitious boundaries in population differentiation (phylogenetic network, pairwise F_ST and STRUCTURE analyses) within and between species assessed using the MHC haplotypes were different from those revealed by NAMs associated closely with speciation, geographical distribution, domestication and management histories. In addition, the mean F_ST was significantly correlated negatively with the number of observed alleles (N_A), observed (H_O) and expected (H_E) heterozygosity and polymorphism information content (PIC) (P < 0.05) in the MHC haplotype dataset while there was no correction of the mean F_ST estimates (P> 0.05) between the MHC haplotype and NAMs datasets. Analysis of molecular variance (AMOVA) revealed a lower percentage of total variance (PTV) between species/groups based on the MHC-linked microsatellites than NAMs. Therefore, it was inferred that individuals within populations accumulated as many MHC variants as possible to increase their heterozygosity and thus the survival rate of their affiliated populations and species, which eventually reduced population differentiation and thereby complicated their classification and phylogenetic relationship inference. In summary, host-pathogen coevolution and heterozygote advantage, rather than demographic history, act as key driving forces shaping the MHC diversity within the populations and determining the interspecific MHC diversity.

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.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.