Extensive MHC Class II B Gene Duplication in a Passerine, the Common Yellowthroat (Geothlypis trichas)

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.

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.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

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

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

Trans‐specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters ( Etheostoma )

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans‐specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep‐sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

Contemporary evolution of the innate immune receptor gene TLR3 in an isolated vertebrate population

Understanding where genetic variation exists, and how it influences fitness within populations is important from an evolutionary and conservation perspective. Signatures of past selection suggest that pathogen-mediated balancing selection is a key driver of immunogenetic variation, but studies tracking contemporary evolution are needed to help resolve the evolutionary forces and mechanism at play. Previous work in a bottlenecked population of Seychelles warblers (Acrocephalus sechellensis) show that functional variation has been maintained at the viral-sensing Toll-like receptor 3 (TLR3) gene, including one nonsynonymous SNP, resulting in two alleles. Here, we characterise evolution at this TLR3 locus over a 25-year period within the original remnant population of the Seychelles warbler, and in four other derived, populations. Results show a significant and consistent temporal decline in the frequency of the TLR3^C allele in the original population, and that similar declines in the TLR3^C allele frequency occurred in all the derived populations. Individuals (of both sexes) with the TLR3^CC genotype had lower survival, and males - but not females - that carry the TLR3^C allele had significantly lower lifetime reproductive success than those with only the TLR3^A allele. These results indicate that positive selection on the TLR3^A allele, caused by an as yet unknown agent, is driving TLR3 evolution in the Seychelles warbler. No evidence of heterozygote advantage was detected. However, whether the positive selection observed is part of a longer-term pattern of balancing selection (through fluctuating selection or rare-allele advantage) cannot be resolved without tracking the TLR3^C allele over an extended time period.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

Association between louse abundance and MHC II supertypes in Galápagos mockingbirds

Major histocompatibility complex class II (MHC II) is an essential molecule triggering the adaptive immune response by the presentation of pathogens to helper T cells. The association between individual MHC II variants and various parasites has become a frequent finding in studies of vertebrate populations. However, although bird ectoparasites have a significant effect on their host's fitness, and the host's immune system can regulate ectoparasitic infections, no study has yet investigated the association between MHC II polymorphism and ectoparasite infection in the populations of free-living birds. Here, we test whether an association exists between the abundance of a chewing louse (Myrsidea nesomimi) and MHC II polymorphism of its hosts, the Galápagos mockingbirds (Mimus). We have found that the presence of two MHC II supertypes (functionally differentiated clusters) was significantly associated with louse abundance. This pattern supports the theory that a co-evolutionary interaction stands behind the maintenance of MHC polymorphism. Moreover, we have found a positive correlation between louse abundance and heterophil/lymphocyte ratio (an indicator of immunological stress) that serves as an additional piece of evidence that ectoparasite burden is affected by immunological state of Galápagos mockingbirds.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

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.

"Balancing" balancing selection? Assortative mating at the major histocompatibility complex despite molecular signatures of balancing selection

In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the set of pathogens to which an individual's adaptive immune system can respond. MHC genes are extraordinarily polymorphic, often showing elevated nonsynonymous relative to synonymous sequence variation and sharing presumably ancient polymorphisms between lineages. These patterns likely reflect pathogen-mediated balancing selection, for example, rare-allele or heterozygote advantage. Such selection is often reinforced by disassortative mating at MHC. We characterized exon 2 of MHC class II, corresponding to the hypervariable peptide-binding region, in song sparrows (Melospiza melodia). We compared nonsynonymous to synonymous sequence variation in order to identify positively selected sites; assessed evidence for trans-species polymorphisms indicating ancient balancing selection; and compared MHC similarity of socially mated pairs to expectations under random mating. Six codons showed elevated ratios of nonsynonymous to synonymous variation, consistent with balancing selection, and we characterized several alleles similar to those occurring in at least four other avian families. Despite this evidence for historical balancing selection, mated pairs were significantly more similar at MHC than were randomly generated pairings. Nonrandom mating at MHC thus appears to partially counteract, not reinforce, pathogen-mediated balancing selection in this system. We suggest that in systems where individual fitness does not increase monotonically with MHC diversity, assortative mating may help to avoid excessive offspring heterozygosity that could otherwise arise from long-standing balancing selection.

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

Background

Genes of the Major Histocompatibility Complex (MHC) are essential for adaptive immune response in vertebrates, as they encode receptors that recognize peptides derived from the processing of intracellular (MHC class I) and extracellular (MHC class II) pathogens. High MHC diversity in natural populations is primarily generated and maintained by pathogen-mediated diversifying and balancing selection. It is, however, debated whether selection at the MHC can counterbalance the effects of drift in bottlenecked populations. The aim of this study was to assess allelic diversity of MHC genes in a recently bottlenecked bird of prey, the White-tailed Eagle Haliaeetus albicilla, as well as to compare mechanisms that shaped the evolution of MHC class I and class II in this species.

Results

We showed that significant levels of MHC diversity were retained in the core Central European (Polish) population of White-tailed Eagles. Ten MHC class I and 17 MHC class II alleles were recovered in total and individual birds showed high average MHC diversity (3.80 and 6.48 MHC class I and class II alleles per individual, respectively). Distribution of alleles within individuals provided evidence for the presence of at least three class I and five class II loci the White-tailed Eagle, which suggests recent duplication events. MHC class II showed greater sequence polymorphism than MHC class I and there was much stronger signature of diversifying selection acting on MHC class II than class I. Phylogenetic analysis provided evidence for trans-species similarity of class II, but not class I, sequences, which is likely consistent with stronger balancing selection at MHC class II.

Conclusions

Relatively high MHC diversity retained in the White-tailed Eagles from northern Poland reinforces high conservation value of local eagle populations. At the same time, our study is the first to demonstrate contrasting patterns of allelic diversity and selection at MHC class I and class II in an accipitrid species, supporting the hypothesis that different mechanisms can shape evolutionary trajectories of MHC class I and class II genes.

Electronic supplementary material

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

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

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

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.

Chemical composition of preen wax reflects major histocompatibility complex similarity in songbirds

In jawed vertebrates, genes of the major histocompatibility complex (MHC) play a key role in immunity by encoding cell-surface proteins that recognize and bind non-self antigens. High variability at MHC suggests that these loci may also function in social signalling such as mate choice and kin recognition. This requires that MHC genotype covaries with some perceptible phenotypic trait. In mammals and fish, MHC is signalled chemically through volatile and non-volatile peptide odour cues, facilitating MHC-dependent mate choice and other behaviours. In birds, despite evidence for MHC-dependent mating, candidate mechanisms for MHC signalling remain largely unexplored. However, feather preen wax has recently been implicated as a potential source of odour cues. We examined whether the chemical composition of preen wax correlates with MHC class IIβ genotypes of wild song sparrows (Melospiza melodia). Pairwise chemical distance reflected amino acid distance at MHC for male-female dyads, although not for same-sex dyads. Chemical diversity did not reflect MHC diversity. We used gas chromatography-mass spectrometry (GC-MS) to characterize preen wax compounds, and identified four wax esters that best reflect MHC similarity. Provided songbirds can detect variation in preen wax composition, this cue may allow individuals to assess MHC compatibility of potential mates.

Population Differences at MHC Do Not Explain Enhanced Resistance of Song Sparrows to Local Parasites

Infectious disease represents an emerging threat to natural populations, particularly when hosts are more susceptible to novel parasites (allopatric) than to parasites from the local area (sympatric). This pattern could arise through evolutionary processes (host populations become adapted to their local parasites and genetically differentiated from other populations at immune-related loci) and/or through ecological interactions (host individuals develop resistance to local parasites through previous exposure). The relative importance of these candidate mechanisms remains unclear. In jawed vertebrates, genes of the major histocompatibility complex (MHC) play a fundamental role in immunity and are compelling candidates for spatially varying selection. We recently showed that song sparrows (Melospiza melodia) are more susceptible to allopatric than to sympatric strains of malaria (Plasmodium). In the current study, to determine whether population differences at MHC explain this pattern, we characterized the peptide-binding regions of MHC (classes I and II) of birds that did or did not become infected in the previous experiment. We recovered up to 4 alleles per individual at class I, implying at least 2 loci, and up to 26 alleles per individual at class II, implying at least 13 loci. Individuals with more class I alleles were less likely to become infected by Plasmodium, consistent with parasite-mediated balancing selection. However, we found no evidence for population genetic differentiation at either class of MHC, based on 36 individuals sequenced. Resistance to sympatric parasites previously described for this system likely stems from individuals' prior immune experience, not from population differentiation and locally protective alleles at MHC.

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

BACKGROUND

Gene duplication has led to a most remarkable adaptation involved in vertebrates' host-pathogen arms-race, the major histocompatibility complex (MHC). However, MHC duplication history is as yet poorly understood in non-mammalian vertebrates, including birds.

RESULTS

Here, we provide evidence for the evolution of two ancient avian MHC class IIB (MHCIIB) lineages by a duplication event prior to the radiation of all extant birds >100 million years ago, and document the role of concerted evolution in eroding the footprints of the avian MHCIIB duplication history.

CONCLUSIONS

Our results suggest that eroded footprints of gene duplication histories may mimic birth-death evolution and that in the avian MHC the presence of the two lineages may have been masked by elevated rates of concerted evolution in several taxa. Through the presence of a range of intermediate evolutionary stages along the homogenizing process of concerted evolution, the avian MHCIIB provides a remarkable illustration of the erosion of multigene family duplication history.

Polymorphism in the major histocompatibility complex (MHC class II B) genes of the Rufous-backed Bunting (Emberiza jankowskii)

Genetic diversity is one of the pillars of conservation biology research. High genetic diversity and abundant genetic variation in an organism may be suggestive of capacity to adapt to various environmental changes. The major histocompatibility complex (MHC) is known to be highly polymorphic and plays an important role in immune function. It is also considered an ideal model system to investigate genetic diversity in wildlife populations. The Rufous-backed Bunting (Emberiza jankowskii) is an endangered species that has experienced a sharp decline in both population and habitat size. Many historically significant populations are no longer present in previously populated regions, with only three breeding populations present in Inner Mongolia (i.e., the Aolunhua, Gahaitu and Lubei557 populations). Efforts focused on facilitating the conservation of the Rufous-backed Bunting (Emberiza jankowskii) are becoming increasingly important. However, the genetic diversity of E. jankowskii has not been investigated. In the present study, polymorphism in exon 2 of the MHCIIB of E. jankowskii was investigated. This polymorphism was subsequently compared with a related species, the Meadow Bunting (Emberiza cioides). A total of 1.59 alleles/individual were detected in E. jankowskii and 1.73 alleles/individual were identified in E. cioides. The maximum number of alleles per individual from the three E. jankowskii populations suggest the existence of at least three functional loci, while the maximum number of alleles per individual from the three E. cioides populations suggest the presence of at least four functional loci. Two of the alleles were shared between the E. jankowskii and E. cioides. Among the 12 unique alleles identified in E. jankowskii, 10.17 segregating sites per allele were detected, and the nucleotide diversity was 0.1865. Among the 17 unique alleles identified in E. cioides, eight segregating sites per allele were detected, and the nucleotide diversity was 0.1667. Overall, compared to other passerine birds, a relatively low level of MHC polymorphism was revealed in E. jankowskii, which was similar to that in E. cioides. Positive selection was detected by PAML/SLAC/FEL analyses in the region encoding the peptide-binding region in both species, and no recombination was detected. Phylogenetic analysis showed that the alleles from E. jankowskii and E. cioides belong to the same clade and the two species shared similar alleles, suggesting the occurrence of a trans-species polymorphism between the two Emberiza species.

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.

MHC class II β exon 2 variation in pardalotes (Pardalotidae) is shaped by selection, recombination and gene conversion

The high levels of polymorphism and allelic diversity which characterise genes in the major histocompatibility complex (MHC) are thought to be generated and maintained through the combined effects of different evolutionary processes. Here, we characterised exon 2 of the MHC class II β genes in two congeneric passerine species, the spotted (Pardalotus punctatus) and striated pardalote (Pardalotus striatus). We estimated the levels of allelic diversity and tested for signatures of recombination, gene conversion and balancing selection to determine if these processes have influenced MHC variation in the two species. Both species showed high levels of polymorphism and allelic diversity, as well as evidence of multiple gene loci and putative pseudogenes based on the presence of stop codons. We found higher levels of MHC diversity in the striated pardalote than the spotted pardalote, based on the levels of individual heterozygosity, sequence divergence and number of polymorphic sites. The observed differences may reflect variable selection pressure on the species, resulting from differences in patterns of movement among populations. We identified strong signatures of historical balancing selection, recombination and gene conversion at the sequence level, indicating that MHC variation in the two species has been shaped by a combination of processes.

Early Duplication of a Single MHC IIB Locus Prior to the Passerine Radiations

A key characteristic of MHC genes is the persistence of allelic lineages over macroevolutionary periods, often through multiple speciation events. This phenomenon, known as trans-species polymorphism (TSP), is well documented in several major taxonomic groups, but has less frequently been observed in birds. The order Passeriformes is arguably the most successful terrestrial vertebrate order in terms of diversity of species and ecological range, but the reasons for this success remain unclear. Passerines exhibit the most highly duplicated MHC genes of any major vertebrate taxonomic group, which may generate increased immune response relative to other avian orders with fewer MHC loci. Here, we describe phylogenetic patterns of the MHC IIB in the passerine family Corvidae. Our results indicate wide-spread TSP within this family, with at least four supported MHC IIB allelic lineages that predate speciation by many millions of years. Markov chain Monte Carlo simulations indicate that divergence of these lineages occurred near the time of the divergence of the Passeriformes and other avian orders. We suggest that the current MHC diversity observed in passerines is due in part to the multiple duplication of a single MHC locus, DAB1, early in passerine evolution and that subsequent duplications of these paralogues have contributed to the enormous success of this order by increasing their ability to recognize and mount immune responses to novel pathogens.

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC-disassortative mate choice. However, many species lack this expected pattern of MHC-disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus-specific primers for high-throughput sequencing of two expressed MHC Class II B genes in Leach's storm-petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene-specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC-dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection.

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.

Characterization of class II β chain major histocompatibility complex genes in a family of Hawaiian honeycreepers: 'amakihi (Hemignathus virens)

Hawaiian honeycreepers (Drepanidinae) have evolved in the absence of mosquitoes for over five million years. Through human activity, mosquitoes were introduced to the Hawaiian archipelago less than 200 years ago. Mosquito-vectored diseases such as avian malaria caused by Plasmodium relictum and Avipoxviruses have greatly impacted these vulnerable species. Susceptibility to these diseases is variable among and within species. Due to their function in adaptive immunity, the role of major histocompatibility complex genes (Mhc) in disease susceptibility is under investigation. In this study, we evaluate gene organization and levels of diversity of Mhc class II β chain genes (exon 2) in a captive-reared family of Hawaii 'amakihi (Hemignathus virens). A total of 233 sequences (173 bp) were obtained by PCR+1 amplification and cloning, and 5720 sequences were generated by Roche 454 pyrosequencing. We report a total of 17 alleles originating from a minimum of 14 distinct loci. We detected three linkage groups that appear to represent three distinct haplotypes. Phylogenetic analysis revealed one variable cluster resembling classical Mhc sequences (DAB) and one highly conserved, low variability cluster resembling non-classical Mhc sequences (DBB). High net evolutionary divergence values between DAB and DBB resemble that seen between chicken BLB system and YLB system genes. High amino acid identity among non-classical alleles from 12 species of passerines (DBB) and four species of Galliformes (YLB) was found, suggesting that these non-classical passerine sequences may be related to the Galliforme YLB sequences.

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.

Molecular characterization of classical and nonclassical MHC class I genes from the golden pheasant (Chrysolophus pictus)

Classical major histocompatibility complex (MHC) class I allelic polymorphism is essential for competent antigen presentation. To improve the genotyping efforts in the golden pheasant, it is necessary to differentiate more accurately between classical and nonclassical class I molecules. In our study, all MHC class I genes were isolated from one golden pheasant based on two overlapping PCR amplifications. In total, six full-length class I nucleotide sequences (A-F) were identified, and four were novel. Two (A and C) belonged to the IA1 gene, two (B and D) were alleles derived from the IA2 gene through transgene amplification, and two (E and F) comprised a third novel locus, IA3 that was excluded from the core region of the golden pheasant MHC-B. IA1 and IA2 exhibited the broad expression profiles characteristic of classical loci, while IA3 showed no expression in multiple tissues and was therefore defined as a nonclassical gene. Phylogenetic analysis indicated that the three IA genes in the golden pheasant share a much closer evolutionary relationship than the corresponding sequences in other galliform species. This observation was consistent with high sequence similarity among them, which likely arises from the homogenizing effect of recombination. Our careful distinction between the classical and nonclassical MHC class I genes in the golden pheasant lays the foundation for developing locus-specific genotyping and establishing a good molecular marker system of classical MHC I loci.

Characterization of a non-classical MHC class II gene in the vulnerable Chinese egret (Egretta eulophotes)

Genes of the major histocompatibility complex (MHC) are valuable makers of adaptive genetic variation in evolutionary ecology research, yet the non-classical MHC genes remain largely unstudied in wild vertebrates. In this study, we have characterized the non-classical MHC class II gene, Egeu-DAB4, in the vulnerable Chinese egret (Ciconiiformes, Ardeidae, Egretta eulophotes). Gene expression analyses showed that Egeu-DAB4 gene had a restricted tissue expression pattern, being expressed in seven examined tissues including the liver, heart, kidney, esophagus, stomach, gallbladder, and intestine, but not in muscle. With respect to polymorphism, only one allele of exon 2 was obtained from Egeu-DAB4 using asymmetric PCR, indicating that Egeu-DAB4 is genetically monomorphic in exon 2. Comparative analyses showed that Egeu-DAB4 had an unusual sequence, with amino acid differences suggesting that its function may differ from those of classical MHC genes. Egeu-DAB4 gene was only found in 30.56-36.56 % of examined Chinese egret individuals. Phylogenetic analysis showed a closer relationship between Egeu-DAB4 and the DAB2 genes in nine other ardeid species. These new findings provide a foundation for further studies to clarify the immunogenetics of non-classical MHC class II gene in the vulnerable Chinese egret and other ciconiiform birds.

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.

Different ornaments signal male health and MHC variation in two populations of a warbler

Male traits that signal health and vigour are used by females to choose better quality mates, but in some cases the male trait selected by females differs among populations. Multiple male traits can be maintained through female mate choice if both traits are equally honest indicators of male quality, but tests of this prediction are rare. By choosing males based on such traits, females could gain direct benefits from males (assistance with parental care), but when females choose extra-pair mates based on these traits, females gain only male sperm, and potentially indirect genetic benefits for their offspring. In common yellowthroats (Geothylpis trichas), female choice of extra-pair mates targets two different plumage ornaments: the black mask in a Wisconsin population and the yellow bib in a New York population. Previously, we found that the black mask in Wisconsin is related to greater major histocompatibility complex (MHC) class II variation, which in turn signals better survival and disease resistance. In this study, we examined the signalling function of the yellow bib in New York to test whether it signals the same aspects of male quality as the black mask in Wisconsin. As predicted, we found that the yellow bib in New York is most closely associated with MHC variation, which also signals survival and resistance to blood parasites. Thus, the ornament preferred by females differs between the two populations, but the different ornaments signal similar aspects of male health and genetic quality, specifically information regarding MHC variation and potential indirect genetic benefits to females.

Evidence for multiple MHC class II β loci in New Zealand's critically endangered kakapo, Strigops habroptilus

Immunologically important genes of the major histocompatibility complex (MHC) have been characterized in a number of avian species with the general finding of considerable variation in size and structural organization among organisms. A range of nonpasserines which represent early-diverging Neoave lineages have been described as having only one MHC class II β locus potentially leading to the conclusion that this is the ancestral condition. Here, we examine the monotypic, early-diverging, critically endangered kakapo, Strigops habroptilus, for allelic variation at MHC class II β exon 2, as part of species' recovery efforts. We found two to four confirmed sequence variants per individual indicating the presence of more than one MHC class II β locus. Given the kakapo's basal evolutionary status, evidence for multiple MHC class II β loci seems to counter the proposed mono-locus history of modern birds. However, MHC gene duplication, maintenance, and loss among and within bird species may confound avian relationships making it difficult to elucidate the ancestral state. This study adds essential data for disentangling the course of MHC structural evolution in birds.

Characterization, polymorphism and selection of major histocompatibility complex (MHC) DAB genes in vulnerable Chinese egret (Egretta eulophotes)

The major histocompatibility complex (MHC) is an excellent molecular marker for the studies of evolutionary ecology and conservation genetics because it is a family of highly polymorphic genes that play a key role in vertebrate immune response. In this study, the functional genes of MHC Class II B (DAB) were isolated for the first time in a vulnerable species, the Chinese egret (Egrettaeulophotes). Using a full length DNA and cDNA produced by PCR and RACE methods, four potential MHC DAB loci were characterized in the genome of this egret and all four were expressed in liver and blood. At least four copies of the MHC gene complex were similar to two copies of the minimal essential MHC complex of chicken, but are less complex than the multiple copies expressed in passerine species. In MHC polymorphism, 19 alleles of exon 2 were isolated from 48 individuals using PCR. No stop codons or frameshift mutations were found in any of the coding regions. The signatures of positive selection detected in potential peptide-binding regions by Bayesian analysis, suggesting that all of these genes were functional. These data will provide the fundamental basis for further studies to elucidate the mechanisms and significance of MHC molecular adaptation in vulnerable Chinese egret and other ardeids.

Characteristics of MHC class I genes in house sparrows Passer domesticus as revealed by long cDNA transcripts and amplicon sequencing

In birds the major histocompatibility complex (MHC) organization differs both among and within orders; chickens Gallus gallus of the order Galliformes have a simple arrangement, while many songbirds of the order Passeriformes have a more complex arrangement with larger numbers of MHC class I and II genes. Chicken MHC genes are found at two independent loci, classical MHC-B and non-classical MHC-Y, whereas non-classical MHC genes are yet to be verified in passerines. Here we characterize MHC class I transcripts (α1 to α3 domain) and perform amplicon sequencing using a next-generation sequencing technique on exon 3 from house sparrow Passer domesticus (a passerine) families. Then we use phylogenetic, selection, and segregation analyses to gain a better understanding of the MHC class I organization. Trees based on the α1 and α2 domain revealed a distinct cluster with short terminal branches for transcripts with a 6-bp deletion. Interestingly, this cluster was not seen in the tree based on the α3 domain. 21 exon 3 sequences were verified in a single individual and the average numbers within an individual were nine and five for sequences with and without a 6-bp deletion, respectively. All individuals had exon 3 sequences with and without a 6-bp deletion. The sequences with a 6-bp deletion have many characteristics in common with non-classical MHC, e.g., highly conserved amino acid positions were substituted compared with the other alleles, low nucleotide diversity and just a single site was subject to positive selection. However, these alleles also have characteristics that suggest they could be classical, e.g., complete linkage and absence of a distinct cluster in a tree based on the α3 domain. Thus, we cannot determine for certain whether or not the alleles with a 6-bp deletion are non-classical based on our present data. Further analyses on segregation patterns of these alleles in combination with dating the 6-bp deletion through MHC characterization across the genus Passer may solve this matter in the future.

MHC class II B diversity in blue tits: a preliminary study

In this study, we partly characterize major histocompatibility complex (MHC) class II B in the blue tit (Cyanistes caeruleus). A total of 22 individuals from three different European locations: Spain, The Netherlands, and Sweden were screened for MHC allelic diversity. The MHC genes were investigated using both PCR-based methods and unamplified genomic DNA with restriction fragment length polymorphism (RFLP) and southern blots. A total of 13 different exon 2 sequences were obtained independently from DNA and/or RNA, thus confirming gene transcription and likely functionality of the genes. Nine out of 13 alleles were found in more than one country, and two alleles appeared in all countries. Positive selection was detected in the region coding for the peptide binding region (PBR). A maximum of three alleles per individual was detected by sequencing and the RFLP pattern consisted of 4-7 fragments, indicating a minimum number of 2-4 loci per individual. A phylogenetic analysis, demonstrated that the blue tit sequences are divergent compared to sequences from other passerines resembling a different MHC lineage than those possessed by most passerines studied to date.

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

Genes of the Major Histocompatibility Complex (MHC) have become an important marker for the investigation of adaptive genetic variation in vertebrates because of their critical role in pathogen resistance. However, despite significant advances in the last few years the characterization of MHC variation in non-model species still remains a challenging task due to the redundancy and high variation of this gene complex. Here we report the utility of a single pair of primers for the cross-amplification of the third exon of MHC class I genes, which encodes the more polymorphic half of the peptide-binding region (PBR), in oscine passerines (songbirds; Aves: Passeriformes), a group especially challenging for MHC characterization due to the presence of large and complex MHC multigene families. In our survey, although the primers failed to amplify exon 3 from two suboscine passerine birds, they amplified exon 3 of multiple MHC class I genes in all 16 species of oscine songbirds tested, yielding a total of 120 sequences. The 16 songbird species belong to 14 different families, primarily within the Passerida, but also in the Corvida. Using a conservative approach based on the analysis of cloned amplicons (n = 16) from each species, we found between 3 and 10 MHC sequences per individual. Each allele repertoire was highly divergent, with the overall number of polymorphic sites per species ranging from 33 to 108 (out of 264 sites) and the average number of nucleotide differences between alleles ranging from 14.67 to 43.67. Our survey in songbirds allowed us to compare macroevolutionary dynamics of exon 3 between songbirds and non-passerine birds. We found compelling evidence of positive selection acting specifically upon peptide-binding codons across birds, and we estimate the strength of diversifying selection in songbirds to be about twice that in non-passerines. Analysis using comparative methods suggest weaker evidence for a higher GC content in the 3rd codon position of exon 3 in non-passerine birds, a pattern that contrasts with among-clade GC patterns found in other avian studies and may suggests different mutational mechanisms. Our primers represent a useful tool for the characterization of functional and evolutionarily relevant MHC variation across the hyperdiverse songbirds.

In silico peptide-binding predictions of passerine MHC class I reveal similarities across distantly related species, suggesting convergence on the level of protein function

The major histocompatibility complex (MHC) genes are the most polymorphic genes found in the vertebrate genome, and they encode proteins that play an essential role in the adaptive immune response. Many songbirds (passerines) have been shown to have a large number of transcribed MHC class I genes compared to most mammals. To elucidate the reason for this large number of genes, we compared 14 MHC class I alleles (α1-α3 domains), from great reed warbler, house sparrow and tree sparrow, via phylogenetic analysis, homology modelling and in silico peptide-binding predictions to investigate their functional and genetic relationships. We found more pronounced clustering of the MHC class I allomorphs (allele specific proteins) in regards to their function (peptide-binding specificities) compared to their genetic relationships (amino acid sequences), indicating that the high number of alleles is of functional significance. The MHC class I allomorphs from house sparrow and tree sparrow, species that diverged 10 million years ago (MYA), had overlapping peptide-binding specificities, and these similarities across species were also confirmed in phylogenetic analyses based on amino acid sequences. Notably, there were also overlapping peptide-binding specificities in the allomorphs from house sparrow and great reed warbler, although these species diverged 30 MYA. This overlap was not found in a tree based on amino acid sequences. Our interpretation is that convergent evolution on the level of the protein function, possibly driven by selection from shared pathogens, has resulted in allomorphs with similar peptide-binding repertoires, although trans-species evolution in combination with gene conversion cannot be ruled out.

MHC class I variation in a natural blue tit population (Cyanistes caeruleus)

The major histocompatibility complex (MHC) is central to the vertebrate immune system and its highly polymorphic genes are considered to influence several life-history traits of individuals. To characterize the MHC in a natural population of blue tits (Cyanistes caeruleus) we investigated the class I exon 3 diversity of more than 900 individuals. We designed two pairs of motif-specific primers that reliably amplify independent subsets of MHC alleles. Applying denaturing gradient gel electrophoresis (DGGE) we obtained 48 independently inherited units of unique band patterns (DGGE-haplogroups), which were validated in a segregation analysis within 105 families. In a second approach, we extensively sequenced 6 unrelated individuals to confirm that DGGE-haplogroup composition reflects individual allelic variation. The highest number of different DGGE-haplogroups in a single individual corresponded in 19 MHC exon 3 sequences, suggesting a minimum of 10 amplified MHC class I loci in the blue tit. In total, we identified 50 unique functional and 3 non-functional sequences. Functional sequences showed high levels of recombination and strong positive selection in the antigen binding region, whereas nucleotide diversity was comparatively low in the range of all passerine species. Finally, in a phylogenetic comparison of passerine MHC class I exon 3 sequences we discuss conflicting evolutionary signals possibly due to recent gene duplication, recombination events and concerted evolution. Our results indicate that the described method is suitable to effectively explore the MHC diversity and its ecological impacts in blue tits in future studies.

Social and extra-pair mating in relation to major histocompatibility complex variation in common yellowthroats

Females are thought to gain better-quality genes for their offspring by mating with particular males. Genes of the major histocompatibility complex (MHC) play a critical role in adaptive immunity, and several studies have examined female mate choice in relation to MHC variation. In common yellowthroats, females prefer males that have larger black facial masks, an ornament associated with MHC variation, immune function and condition. Here we also tested whether mating patterns are directly correlated with MHC diversity or similarity. Using pyrosequencing, we found that the presence of extra-pair young in the brood was not related to male MHC diversity or similarity between the female and her within-pair mate. Furthermore, extra-pair sires did not differ in overall diversity from males they cuckolded, or in their similarity to the female. MHC diversity is extremely high in this species, and it may limit the ability of females to assess MHC variation in males. Thus, mating may be based on ornaments, such as mask size, which are better indicators of overall male health and genetic quality.

Variation at innate immunity Toll-like receptor genes in a bottlenecked population of a New Zealand robin

Toll-like receptors (TLRs) are an ancient family of genes encoding transmembrane proteins that bind pathogen-specific molecules and initiate both innate and adaptive aspects of the immune response. Our goal was to determine whether these genes show sufficient genetic diversity in a bottlenecked population to be a useful addition or alternative to the more commonly employed major histocompatibility complex (MHC) genotyping in a conservation genetics context. We amplified all known avian TLR genes in a severely bottlenecked population of New Zealand's Stewart Island robin (Petroica australis rakiura), for which reduced microsatellite diversity was previously observed. We genotyped 17-24 birds from a reintroduced island population (including the 12 founders) for nine genes, seven of which were polymorphic. We observed a total of 24 single-nucleotide polymorphisms overall, 15 of which were non-synonymous, representing up to five amino-acid variants at a locus. One locus (TLR1LB) showed evidence of past directional selection. Results also confirmed a passerine duplication of TLR7. The levels of TLR diversity that we observe are sufficient to justify their further use in addressing conservation genetic questions, even in bottlenecked populations.

Ancestral polymorphism at the major histocompatibility complex (MHCIIß) in the Nesospiza bunting species complex and its sister species (Rowettia goughensis)

BACKGROUND

The major histocompatibility complex (MHC) is an important component of the vertebrate immune system and is frequently used to characterise adaptive variation in wild populations due to its co-evolution with pathogens. Passerine birds have an exceptionally diverse MHC with multiple gene copies and large numbers of alleles compared to other avian taxa. The Nesospiza bunting species complex (two species on Nightingale Island; one species with three sub-species on Inaccessible Island) represents a rapid adaptive radiation at a small, isolated archipelago, and is thus an excellent model for the study of adaptation and speciation. In this first study of MHC in Nesospiza buntings, we aim to characterize MHCIIß variation, determine the strength of selection acting at this gene region and assess the level of shared polymorphism between the Nesospiza species complex and its putative sister taxon, Rowettia goughensis, from Gough Island.

RESULTS

In total, 23 unique alleles were found in 14 Nesospiza and 2 R. goughensis individuals encoding at least four presumably functional loci and two pseudogenes. There was no evidence of ongoing selection on the peptide binding region (PBR). Of the 23 alleles, 15 were found on both the islands inhabited by Nesospiza species, and seven in both Nesospiza and Rowettia; indications of shared, ancestral polymorphism. A gene tree of Nesospiza MHCIIß alleles with several other passerine birds shows three highly supported Nesospiza-specific groups. All R. goughensis alleles were shared with Nesospiza, and these alleles were found in all three Nesospiza sequence groups in the gene tree, suggesting that most of the observed variation predates their phylogenetic split.

CONCLUSIONS

Lack of evidence of selection on the PBR, together with shared polymorphism across the gene tree, suggests that population variation of MHCIIß among Nesospiza and Rowettia is due to ancestral polymorphism rather than local selective forces. Weak or no selection pressure could be attributed to low parasite load at these isolated Atlantic islands. The deep divergence between the highly supported Nesospiza-specific sequence Groups 2 and 3, and the clustering of Group 3 close to the distantly related passerines, provide strong support for preserved ancestral polymorphism, and present evidence of one of the rare cases of extensive ancestral polymorphism in birds.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

Maintenance of MHC Class IIB diversity in a recently established songbird population

We examined variation at MHC Class IIB genes in a recently established population of dark-eyed juncos (Junco hyemalis) in a coastal urban environment in southern California, USA relative to an ancestral-range population from a nearby species-typical montane environment. The founding population is estimated to have been quite small, but we predicted that variation at the major histocompatibility complex (MHC) among the founders would nevertheless be preserved owing to the high functional significance of MHC. Previous studies of MHC in songbirds have had varying degrees of success in isolating loci, as passerines show extensive MHC gene duplication. In order to compare diversity in the two populations, we employed two published approaches to sequencing MHC Class II exon 2: direct sequencing with exon-based primers, and traditional cloning and sequencing with intron-based primers. Results from both methods show that the colonist population has maintained high levels of variation. Our results also indicate varying numbers of alleles across individuals, corroborating evidence for gene duplication in songbird MHC. While future studies in songbirds may need to take a genomic approach to fully understand the structure of MHC in this lineage, our results show that it is possible to use traditional methods to reveal functional variation across populations.