Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

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.

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.

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

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

Association between Y-Maze Acquisition Learning and Major Histocompatibility Complex Class II Polymorphisms in Mice

Objective

To explore the association between the acquisition process in the Y-maze and H-2 class II polymorphisms in mice.

Methods

Mice were trained for 5 consecutive days in the Y-maze. The value of the slope of the latent period was considered an indication for the acquisition process. A slope < 0 indicated learning during the training and a slope > 0 indicated no learning. The H-2 polymorphism was determined with PCR amplification, and the correlation between the alleles and the acquisition process was analyzed.

Results

The overall percentage of mice that learned was 46.1%. The percentage of mice that had learned with MudoEb5 (37.9%) was significantly lower than that of mice without MudoEb5 (61.1%; P < 0.05). The percentage of mice that had learned with MudoEb7 (26.1%) was significantly lower than that of mice without MudoEb7 (51.9%; P < 0.05).

Conclusions

The major histocompatibility complex (MHC) and other alleles may be involved in the acquisition process. There may be a biological basis for learning in mice.

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

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

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.

Balancing selection and genetic drift create unusual patterns of MHCIIβ variation in Galápagos mockingbirds

The extracellular subunit of the major histocompatibility complex MHCIIβ plays an important role in the recognition of pathogens and the initiation of the adaptive immune response of vertebrates. It is widely accepted that pathogen-mediated selection in combination with neutral micro-evolutionary forces (e.g. genetic drift) shape the diversity of MHCIIβ, but it has proved difficult to determine the relative effects of these forces. We evaluated the effect of genetic drift and balancing selection on MHCIIβ diversity in 12 small populations of Galápagos mockingbirds belonging to four different species, and one larger population of the Northern mockingbird from the continental USA. After genotyping MHCIIβ loci by high-throughput sequencing, we applied a correlational approach to explore the relationships between MHCIIβ diversity and population size by proxy of island size. As expected when drift predominates, we found a positive effect of population size on the number of MHCIIβ alleles present in a population. However, the number of MHCIIβ alleles per individual and number of supertypes were not correlated with population size. This discrepancy points to an interesting feature of MHCIIβ diversity dynamics: some levels of diversity might be shaped by genetic drift while others are independent and possibly maintained by balancing selection.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Spatial and temporal variation at major histocompatibility complex class IIB genes in the endangered Blakiston's fish owl

INTRODUCTION

Quantifying intraspecific genetic variation in functionally important genes, such as those of the major histocompatibility complex (MHC), is important in the establishment of conservation plans for endangered species. The MHC genes play a crucial role in the vertebrate immune system and generally show high levels of diversity, which is likely due to pathogen-driven balancing selection. The endangered Blakiston's fish owl (Bubo blakistoni) has suffered marked population declines on Hokkaido Island, Japan, during the past several decades due to human-induced habitat loss and fragmentation. We investigated the spatial and temporal patterns of genetic diversity in MHC class IIβ genes in Blakiston's fish owl, using massively parallel pyrosequencing.

RESULTS

We found that the Blakiston's fish owl genome contains at least eight MHC class IIβ loci, indicating recent gene duplications. An analysis of sequence polymorphism provided evidence that balancing selection acted in the past. The level of MHC variation, however, was low in the current fish owl populations in Hokkaido: only 19 alleles were identified from 174 individuals. We detected considerable spatial differences in MHC diversity among the geographically isolated populations. We also detected a decline of MHC diversity in some local populations during the past decades.

CONCLUSIONS

Our study demonstrated that the current spatial patterns of MHC variation in Blakiston's fish owl populations have been shaped by loss of variation due to the decline and fragmentation of populations, and that the short-term effects of genetic drift have counteracted the long-term effects of balancing selection.

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.

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

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

Selection and trans-species polymorphism of major histocompatibility complex class II genes in the order Crocodylia

Major Histocompatibility Complex (MHC) class II genes encode for molecules that aid in the presentation of antigens to helper T cells. MHC characterisation within and between major vertebrate taxa has shed light on the evolutionary mechanisms shaping the diversity within this genomic region, though little characterisation has been performed within the Order Crocodylia. Here we investigate the extent and effect of selective pressures and trans-species polymorphism on MHC class II α and β evolution among 20 extant species of Crocodylia. Selection detection analyses showed that diversifying selection influenced MHC class II β diversity, whilst diversity within MHC class II α is the result of strong purifying selection. Comparison of translated sequences between species revealed the presence of twelve trans-species polymorphisms, some of which appear to be specific to the genera Crocodylus and Caiman. Phylogenetic reconstruction clustered MHC class II α sequences into two major clades representing the families Crocodilidae and Alligatoridae. However, no further subdivision within these clades was evident and, based on the observation that most MHC class II α sequences shared the same trans-species polymorphisms, it is possible that they correspond to the same gene lineage across species. In contrast, phylogenetic analyses of MHC class II β sequences showed a mixture of subclades containing sequences from Crocodilidae and/or Alligatoridae, illustrating orthologous relationships among those genes. Interestingly, two of the subclades containing sequences from both Crocodilidae and Alligatoridae shared specific trans-species polymorphisms, suggesting that they may belong to ancient lineages pre-dating the divergence of these two families from the common ancestor 85-90 million years ago. The results presented herein provide an immunogenetic resource that may be used to further assess MHC diversity and functionality in Crocodylia.

Evolution of MHC class I in the order Crocodylia

The major histocompatibility complex (MHC) is a dynamic genomic region with an essential role in the adaptive immunity of jawed vertebrates. The evolution of the MHC has been dominated by gene duplication and gene loss, commonly known as the birth-and-death process. Evolutionary studies of the MHC have mostly focused on model species. However, the investigation of this region in non-avian reptiles is still in its infancy. To provide insights into the evolutionary mechanisms that have shaped the diversity of this region in the Order Crocodylia, we investigated MHC class I exon 3, intron 3, and exon 4 across 20 species of the families Alligatoridae and Crocodilidae. We generated 124 DNA sequences and identified 31 putative functional variants as well as 14 null variants. Phylogenetic analyses revealed three gene groups, all of which were present in Crocodilidae but only one in Alligatoridae. Within these groups, variants generally appear to cluster at the genus or family level rather than in species-specific groups. In addition, we found variation in gene copy number and some indication of interlocus recombination. These results suggest that MHC class I in Crocodylia underwent independent events of gene duplication, particularly in Crocodilidae. These findings enhance our understanding of MHC class I evolution and provide a preliminary framework for comparative studies of other non-avian reptiles as well as diversity assessment within Crocodylia.

MHC class IIB exon 2 polymorphism in the Grey partridge (Perdix perdix) is shaped by selection, recombination and gene conversion

Among bird species, the most studied major histocompatibility complex (MHC) is the chicken MHC. Although the number of studies on MHC in free-ranging species is increasing, the knowledge on MHC variation in species closely related to chicken is required to understand the peculiarities of bird MHC evolution. Here we describe the variation of MHC class IIB (MHCIIB) exon 2 in a population of the Grey partridge (Perdix perdix), a species of high conservation concern throughout Europe and an emerging galliform model in studies of sexual selection. We found 12 alleles in 108 individuals, but in comparison to other birds surprisingly many sites show signatures of historical positive selection. Individuals displayed between two to four alleles both on genomic and complementary DNA, suggesting the presence of two functional MHCIIB loci. Recombination and gene conversion appear to be involved in generating MHCIIB diversity in the Grey partridge; two recombination breakpoints and several gene conversion events were detected. In phylogenetic analysis of galliform MHCIIB, the Grey partridge alleles do not cluster together, but are scattered through the tree instead. Thus, our results indicate that the Grey partridge MHCIIB is comparable to most other galliforms in terms of copy number and population polymorphism.

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.

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.

MHC class I of saltwater crocodiles (Crocodylus porosus): polymorphism and balancing selection

Saltwater crocodiles are in high demand for the production of luxury fashion items. However, their susceptibility to disease incurs substantial losses and it is hoped to be able to genetically select these animals for disease resistance. So far, this has only been enabled by phenotypic selection. Investigating the major histocompatibility complex (MHC) could provide insight into the ability of an individual to respond to pathogens acting as a selective pressure on the host. Here, we assessed genetic diversity and a role of selection in shaping the diversity of MHC class I exon 3 among 42 saltwater crocodiles from nine river basins in the Northern Territory, Australia. We generated 640 sequences using cloning and sequencing methods and identified 43 MHC variants among them. Phylogenetic analyses clustered these variants into two major clades, which may suggest two gene lineages. We found the number of variants within an individual varying between one and seven, indicating that there are at least four gene loci in this species. Selection detection analyses revealed an elevated ratio of nonsynonymous to synonymous substitutions (mean = 1.152 per codon), suggesting balancing selection. Population differentiation analyses revealed that the MHC did not show structuring among the river basins, and there were some shared variants among them. This may be a result of possible gene flow and/or similar selection pressures among populations. These findings provide background knowledge to identify potential MHC markers, which could be used for selecting genetically variable individuals for future disease associations. All MHC class I exon 3 sequences reported in this paper were submitted to the GenBank database with following accession numbers: HQ008785-HQ008789, HQ008791-HQ008798, HQ008808-HQ008815, HQ008824, HQ008826-HQ008830, HQ008835, HQ008839, HQ008842-HQ008850, and JX023536-JX023540.

Evolutionary genetics of MHC class II beta genes in the brown hare, Lepus europaeus

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na = 11, D_est = 0.071 and Na = 15, D_est = 0.409, respectively), but not for the DPB gene (Na = 4, D_est = 0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω > 1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids.

Genotyping of black grouse MHC class II B using reference Strand-Mediated Conformational Analysis (RSCA)

BACKGROUND

The Major Histocompatibility Complex (MHC) is a cluster of genes involved in the vertebrate immune system and includes loci with an extraordinary number of alleles. Due to the complex evolution of MHC genes, alleles from different loci within the same MHC class can be very similar and therefore difficult to assign to separate loci. Consequently, single locus amplification of MHC genes is hard to carry out in species with recently duplicated genes in the same MHC class, and multiple MHC loci have to be genotyped simultaneously. Since amplified alleles have the same length, accurate genotyping is difficult. Reference Strand-Mediated Conformational Analysis (RSCA), which is increasingly used in studies of natural populations with multiple MHC genes, is a genotyping method capable to provide high resolution and accuracy in such cases.

FINDINGS

We adapted the RSCA method to genotype multiple MHC class II B (BLB) genes in black grouse (Tetrao tetrix), a non-model galliform bird species, using a 96-Capillary Array Electrophoresis, the MegaBACE™ 1000 DNA Analysing System (GE Healthcare). In this study we used fluorescently labelled reference strands from both black grouse and hazel grouse and observed good agreement between RSCA and cloning/sequencing since 71 alleles were observed by cloning/sequencing and 76 alleles by RSCA among the 24 individuals included in the comparison. At the individual level however, there was a trend towards more alleles scored with RSCA (1-6 per individual) than cloning/sequencing (1-4 per individual). In 63% of the pair-wise comparison, the identical allele was scored in RSCA as in cloning/sequencing. Nine out of 24 individuals had the same number of alleles in RSCA as in cloning/sequencing. Our RSCA protocol allows a faster RSCA genotyping than presented in many other RSCA studies.

CONCLUSIONS

In this study, we have developed the RSCA typing method further to work on a 96-Capillary Array Electrophoresis (MegaBACE™ 1000). Our RSCA protocol can be applied to fast and reliable screening of MHC class II B diversity of black grouse populations. This will facilitate future large-scale population studies of black grouse and other galliformes species with multiple inseparable MHC loci.

454 sequencing reveals extreme complexity of the class II Major Histocompatibility Complex in the collared flycatcher

Background

Because of their functional significance, the Major Histocompatibility Complex (MHC) class I and II genes have been the subject of continuous interest in the fields of ecology, evolution and conservation. In some vertebrate groups MHC consists of multiple loci with similar alleles; therefore, the multiple loci must be genotyped simultaneously. In such complex systems, understanding of the evolutionary patterns and their causes has been limited due to challenges posed by genotyping.

Results

Here we used 454 amplicon sequencing to characterize MHC class IIB exon 2 variation in the collared flycatcher, an important organism in evolutionary and immuno-ecological studies. On the basis of over 152,000 sequencing reads we identified 194 putative alleles in 237 individuals. We found an extreme complexity of the MHC class IIB in the collared flycatchers, with our estimates pointing to the presence of at least nine expressed loci and a large, though difficult to estimate precisely, number of pseudogene loci. Many similar alleles occurred in the pseudogenes indicating either a series of recent duplications or extensive concerted evolution. The expressed alleles showed unambiguous signals of historical selection and the occurrence of apparent interlocus exchange of alleles. Placing the collared flycatcher's MHC sequences in the context of passerine diversity revealed transspecific MHC class II evolution within the Muscicapidae family.

Conclusions

454 amplicon sequencing is an effective tool for advancing our understanding of the MHC class II structure and evolutionary patterns in Passeriformes. We found a highly dynamic pattern of evolution of MHC class IIB genes with strong signals of selection and pronounced sequence divergence in expressed genes, in contrast to the apparent sequence homogenization in pseudogenes. We show that next generation sequencing offers a universal, affordable method for the characterization and, in perspective, genotyping of MHC systems of virtually any complexity.

A locus-wide approach to assessing variation in the avian MHC: the B-locus of the wild turkey

Studies of major histocompatibility complex (MHC) diversity in non-model vertebrates typically focus on structure and sequence variation in the antigen-presenting loci: the highly variable and polymorphic class I and class IIB genes. Although these studies provide estimates of the number of genes and alleles/locus, they often overlook variation in functionally related and co-inherited genes important in the immune response. This study utilizes the sequence of the MHC B-locus derived from a commercial turkey to investigate MHC variation in wild birds. Sequences were obtained for nine interspersed MHC amplicons (non-class I/II) from each of 40 birds representing 3 subspecies of wild turkey (Meleagris gallopavo). Analysis of aligned sequences identified 238 single-nucleotide variants approximately one-third of which had minor allele frequencies >0.2 in the sampled birds. PHASE analysis identified 70 prospective MHC haplotypes in the wild turkeys, whereas a combined analysis with commercial birds identified almost 100 haplotypes in the species. Denaturing gradient gel electrophoresis (DGGE) of the class IIB loci was used to test the efficacy of single-nucleotide polymorphism (SNP) haplotyping to capture locus-wide variation. Diversity in SNP haplotypes and haplotype sharing among individuals was directly reflected in the DGGE patterns. Utilization of a reference haplotype to sequence interspersed regions of the MHC has significant advantages over other methods of surveying diversity while identifying high-frequency SNPs for genotyping. SNP haplotyping provides a means to identify both divergent haplotypes and homozygous individuals for assessment of immunological variation in wild and domestic populations.

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca

Background

Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts.

Findings

Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (β chain). Taking the pied flycatcher Ficedula hypoleuca as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species.

Conclusions

Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.