- Research article
- Open Access
New data from basal Australian songbird lineages show that complex structure of MHC class II β genes has early evolutionary origins within passerines
© Balasubramaniam et al. 2016
Received: 10 October 2015
Accepted: 10 May 2016
Published: 21 May 2016
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.
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.
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.
The major histocompatibility complex (MHC) is a complex multigene family that regulates the function of the adaptive branch of the vertebrate immune system. The MHC is comprised of many different types of genes, such as classical and non-classical MHC genes, as well as non-MHC genes such as natural killer cells and tapasin . Classical and non-classical MHC genes may be distinguished from each other in that the former are polymorphic and highly expressed, whereas the latter exhibit lower levels of expression and are often monomorphic . MHC genes may also be grouped as class I, II or III genes based on the types of molecules they encode. Class I genes encode receptors that are presented on the surfaces of most nucleated cells and primarily facilitate immune responses to intracellular pathogens, whereas class II genes are only found on a subset of cells and are associated with immune responses to extracellular pathogens . Class III genes encode molecules involved in the complement component of the immune response rather than the adaptive response .
The high levels of polymorphism characteristic of classical MHC class I and II genes may be largely attributed to the peptide-binding region (PBR), which displays great diversity of alleles and extensive sequence variation among different alleles [3, 4]. The genes of this region encode proteins that form the molecular groove where peptides are bound and then presented to T-cells, resulting in the appropriate immune response being triggered . Each allele typically responds to a category of potential antigens, so individuals and populations with greater MHC diversity may be better able to cope with a range of infections [5, 6].
A number of models of multigene family evolution, such as the birth-and-death model  and the accordion model , have been proposed to explain the generally high levels of polymorphism and the large and variable number of classical MHC loci among individuals and species. According to the birth-and-death model, new genes are created via repeated duplication (birth) and then either retained in the genome or eventually lost (death) after becoming non-functional through deleterious mutations . The accordion model, on the other hand, posits that the number of MHC genes expands and contracts in response to fluctuating selection pressure . These models are not mutually exclusive, and both models predict the redundancy that occurs in the genome of most species, whereby multiple duplicated gene loci exist in both class I and II gene regions . One explanation for this redundancy is that once alleles are generated, they may remain in the genome even in the absence of selection pressure . Alleles may be maintained over long evolutionary time scales through speciation events, which may then resolve as trans-species polymorphisms (TSP) in phylogenetic reconstructions, where similar alleles are present in groups that diverged over millions of years ago such as rodents and primates [10, 11]. On a smaller scale, genetic variation in the MHC may also be generated among alleles and between loci via the processes of point mutation, recombination and gene conversion, and maintained through the influence of balancing selection.
Among birds, the most studied and well-described MHC is that of the chicken (Gallus gallus), which only contains one dominantly expressed molecule of class I and one of class II [6, 7]. The chicken MHC is small and densely arranged; a compact organisation of genes with short introns, small physical size, lack of redundancy, and low overall number of classical class I and II genes with few pseudogenes, has led to the chicken MHC being described as a “minimal essential MHC” [7, 8]. The small size and simplicity of the chicken MHC allows for the co-evolution of genes as haplotypes over considerable periods of time, with no recombination between class I and II genes having been detected over thousands of experiments [7, 9]. The result of this is that the most striking associations between MHC haplotype and resistance or susceptibility to disease have been characterised in chickens [10, 11].
Similarly simple genetic organisation and small numbers of expressed classical MHC loci have been noted in a few galliform species, such as the pheasant (Phasianus colchicus) , black grouse (Tetrao tetrix) , and grey partridge (Perdix perdix) , although the quail (Coturnix japonica) has been found to have higher levels of gene duplication . Variable levels of MHC complexity have been characterised among other avian groups, such as penguins, seabirds, raptors, rails, cranes, and waders [15–20], and passerine species typically show multiple loci with extensive gene duplication, evidence of recombination, high levels of polymorphism, and the presence of pseudogenes [21–25]. This suggests the minimal essential model is not applicable to all avian species, as the passerine MHC appears to be larger and more complex than that described in chickens. The passerine MHC, along with a number of non-passerine lineages, has been characterised as having high rates of concerted evolution and/or recent duplications [2, 12], leading to few reports of orthologous relationships among birds [12, 26–29]. Thus, although the overall function of the MHC is conserved across birds, there appear to be considerable differences in the genomic organisation of the MHC between lineages .
Sample information and MHC allelic diversity
Brown treecreeper Climacteris picumnus
Yellow-tufted honeyeater Lichenostomus melanops
White-plumed honeyeater Lichenostomus penicillatus
Fuscous honeyeater Lichenostomus fuscus
Superb fairy-wren Malurus cyaneus
Spotted pardalote Pardalotus punctatus
Striated pardalote Pardalotus striatus
Results and discussion
In order to assess levels of polymorphism and allelic diversity in the MHC class II genes of basal passerines, we chose the following species: Brown Treecreeper Climacteris picumnus, Superb Fairy-wren Malurus cyaneus¸ Spotted Pardalote Pardalotus punctatus, Striated Pardalote P. striatus, White-plumed honeyeater Lichenostomus penicillatus, Fuscous honeyeater L. fuscus, and Yellow-tufted Honeyeater L. melanops. These seven species are from four families (Climacteridae, Maluridae, Pardalotidae, Meliphagidae) that comprise a prominent component of the Corvida. We amplified a 159 bp region that includes the highly diverse PBR of the MHC class II β exon 2 using cloning and 454 sequencing methods. The high level of coverage provided by 454 sequencing makes it a convenient method for assessing genetic diversity in highly polymorphic, multilocus systems, such as the MHC [25, 35, 36]. We used 454 sequencing and cloning with Sanger sequencing complementarily to assess levels of genetic diversity in seven basal passerine species.
The two methods resulted in different numbers of alleles being identified in each species. We isolated 98 MHC class II β exon 2 sequences from 21 individuals using cloning with Sanger sequencing (see Additional file 1), with the total number of alleles detected in each species ranging from 7 in the Spotted Pardalote to 19 in the Striated Pardalote. Of the 98 sequences isolated, 17 were found in more than one species, resulting in 81 unique alleles across the seven species. Sample sizes were the same in every species and all samples were from the same study region, so the variation in allele numbers across species is unlikely to be a product of sampling design. Within species, the number of alleles identified in each individual ranged from 5 to 14.
454 sequencing was conducted on 13 samples in total (see Additional file 1), across all seven study species. Due to small volumes of available blood samples, different individuals were sequenced using 454 compared to Sanger sequencing, and from a different location in the study area (see Methods). Using 454 sequencing, we obtained 16,257 reads of the 159-bp region showing a complete match to the forward and reverse primers. After a stepwise filtering procedure, we were left with 66 % of the original reads (10,783), with coverage depth ranging from 114 to 570 reads per allele (see Additional file 2). This remaining 66 % of reads comprised a total of 146 alleles, of which 12 were found in more than one species. Of the remaining 134 unique alleles recovered through 454 sequencing, 38 were also identified through cloning. Similar to the sequences isolated through cloning, numbers of alleles varied widely across species and individuals, ranging from 11 in the White-plumed Honeyeater to 26 in the Brown Treecreeper. The isolation of more than two alleles from each individual indicates the presence of multiple MHC class II β loci in each of these species, from a minimum of five loci in the Fuscous Honeyeater to nine loci in the Brown Treecreeper. The conservative filtering process we used would have potentially excluded a larger number of true alleles from individuals with higher coverage compared to individuals with lower coverage, so levels of MHC allelic diversity may be underestimated in some species.
Despite fewer individuals having been screened through 454 sequencing compared to cloning, larger numbers of alleles were identified in total and in each individual (see Additional file 1). Although 454 sequencing has the advantage of avoiding some artefacts associated with bacterial cloning, it is still vulnerable to artefacts arising during PCR and has been shown to produce a higher percentage of sequencing errors than Sanger sequencing [37–40]. In contrast, Sanger sequencing probably represents a smaller proportion of the allelic diversity of a species, compared with 454 sequencing, with a prevalence of locally common alleles while less common (or rare) alleles are missed in cloning. Our use of both these methods suggests 454 sequencing is a much more efficient, although possibly less accurate, way to assess levels of allelic diversity in species and populations.
Polymorphism and allelic diversity of MHC loci
Mean nucleotide and amino acid distances
Phylogenetic relationships among MHC alleles
The basal position of the Superb Fairy-wren sequences may result from recombination generating a pattern of mixed ancestry, which causes recombinant sequences to fall outside the parental clusters . The intermingling of other Superb Fairy-wren sequences within clades in the phylogenetic tree containing other genera may indicate a trans-species mode of evolution, a pattern which has been widely documented in MHC class II β sequences of both passerines and non-passerines [28, 46–50]. Trans-species polymorphism (TSP) is a mechanism by which identical alleles occur in related species by being passed on from ancestral to descendant species . This can be difficult to differentiate from convergent evolution, which can produce similar alleles in different species through independent evolutionary pathways. Being sympatric, the seven species in this study may display signatures of convergent evolution if they are subject to similar parasite or pathogen pressure. To differentiate between signatures of convergent evolution and TSP, we compared phylogenetic reconstructions based on non-synonymous and synonymous substitutions only. If the species are subject to similar selective pressures, the phylogeny based on non-synonymous substitutions should show mixed clustering of sequences across families whereas the phylogeny based on synonymous substitutions should show alleles grouping according to family or species. If, however, the similarity of alleles across species is a result of TSP, the two phylogenies should show similar patterns of allelic clustering across families [48, 52]. A comparison of phylogenies based on non-synonymous and synonymous substitutions showed the latter pattern of clustering, with the majority of clusters in both phylogenies showing alleles from more than one family (see Additional file 3). Mixed clusters predominantly comprised alleles from honeyeaters, pardalotes and fairy-wrens, whereas alleles from treecreepers clustered intraspecifically for the most part. Based on our analysis of exon 2, the phylogenetic signatures support a TSP explanation of the clustering pattern rather than one of convergent evolution. Although this result suggests the presence of TSP, it must be noted that distinguishing between phylogenetic patterns generated by TSP and other non-TSP patterns and processes is problematic when the loci are unknown. The evolution of exon 2 is complex and it is possible that, for example, rates of allele sharing based only on this region are an overestimation, thereby leading to incorrect inferences of TSP [52, 53]. Comparison of phylogenies based on intron and exon regions, as well as other related species with allopatric distributions, and identification of alleles associated with individual loci, would clarify the roles of convergent evolution, TSP, and other processes in generating these patterns of clustering.
Fifteen alleles in the dataset occurred in at least two different species, and there was variation among genera in levels of TSP. In 11 instances, alleles occurred across genera between fairy-wrens and another species. In all other instances TSP occurred within genera, either between the two species of pardalotes or three species of honeyeaters. Treecreepers did not show any instances of TSP with other study species, which may be related to treecreepers being in a separate evolutionary lineage to the other study species (Climacteridae; Fig. 1). Retention of alleles through TSP between treecreepers and the other basal lineages in this study would therefore need to occur over longer evolutionary timeframes of at least 60 million years .
Inspection of amino acid alignments revealed seven sites with residues unique to each of the Climacteridae, Pardalotidae and Meliphagidae, distinguishing sequences in the three families from each other (see Additional file 4). Fairy-wren sequences, being intermingled within these groups, showed the residues particular to whichever group they fell in with. This pattern of clustering by species is not uncommon in birds and could be explained by either a) recent duplication, where genes diverge following speciation and then duplicate to form similar copies of the gene, or b) concerted evolution, where duplication occurs prior to speciation but the duplicated genes become homogenised through gene conversion . Both these processes result in intraspecific loci being more similar to each other than to orthologous loci in other species. Among birds, patterns of orthology may be obscured because of rapid gene duplication and homogenisation of MHC class II genes .
In this study we detected evidence of multiple gene loci, high levels of polymorphism and allelic diversity, and the presence of pseudogenes among classical MHC class II β exon 2 of seven passerine species. Our analyses strongly support a monophyletic grouping of passerines from the Corvida and Passerida, signifying an early evolutionary origin of complex MHC class II β structure in the Passeriformes. Phylogenetic analyses based on non-synonymous and synonymous substitutions supported a trans-species explanation for the mixed clustering of alleles across the different families in this study. However, we cannot rule out the possibility that convergent evolution may have also played a role in generating these phylogenetic signatures. Phylogenetic analyses based on sequences from other MHC class II exon and intron regions may clarify the relative roles of TSP and convergent evolution in generating these signatures. A comparison of 454 sequencing and cloning methods suggested the former is a much more efficient way to assess levels of allelic diversity in species and populations. Continued characterisation of this gene region, as well as other non-coding regions of the MHC, in species from different phylogenetic levels may improve our understanding of the rates of gene conversion, recombination, diversification and homogenisation occurring in the avian MHC. Such research will also provide insights into the evolutionary significance of the disparity in MHC complexity between passerine and non-passerine species.
Sample collection and study design
We collected blood samples from one to three individuals of seven Australian passerine species from Victoria, Australia, as described in Amos et al. (; see Additional file 1). These seven species (Brown Treecreeper Climacteris picumnus, Superb Fairy-wren Malurus cyaneus, Spotted Pardalote Pardalotus punctatus, Striated Pardalote P. striatus, White-plumed honeyeater Lichenostomus penicillatus, Fuscous honeyeater L. fuscus, Yellow-tufted Honeyeater L. melanops) are from four families (Climacteridae, Maluridae, Pardalotidae, Meliphagidae) that comprise a prominent component of basal passerine lineages.
As part of a pilot study we tested previously published MHC class II β primers on the study species (Additional file 5), from which we identified a single set of primers that consistently amplified across all the study species. These primers amplified a 159 bp region of the MHC class II β exon 2 and were used for all our genetic work. To assess levels of polymorphism and allelic diversity in the seven study species, we amplified the 159 bp region of the MHC class II β exon 2 using cloning and 454 sequencing methods. This 159 bp fragment includes the highly diverse PBR as well as fit within the length restrictions for 454 sequencing. The PBR and non-PBR were inferred assuming functional congruence to the human HLA-DR1 molecule .
DNA extraction, PCR, cloning, and sequencing
Blood samples were digested overnight using Proteinase K and DNA was extracted with standard phenol-chloroform protocols for Sanger sequencing , while DNA was extracted from blood samples using a salting-out protocol for 454 sequencing . The resulting DNA was suspended in Tris-EDTA. For Sanger sequencing the degenerate primers 326 and 325  were used to amplify a 159 bp region of the MHC class II β gene, spanning the majority of exon 2. For 454 sequencing we amplified the same 159 bp region of the MHC class II β gene using HPLC-purified fusion primers. The forward fusion primer (5′-CCATCTCATCCCTGCGTGTCTCCGACTCAGNNNNNNNNNNGAGTGYCAYTAYYTNAAYGGYAC-3′) comprised the 454 GS FLX Titanium Primer A, a 10 bp Multiplex Identifier (MID) tag (indicated with Ns) to differentiate among individuals, and the 326 primer sequence (in bold). The reverse fusion primer (5′-CCTATCCCCTGTGTGCCTTGGCAGTCTCAGNNNNNNNNNNGTAGTTGTGNCKGCAGTANSTGTCCAC-3′) similarly comprised the 454 GS FLX Titanium Primer B, a 10 bp MID tag (indicated with Ns) and the 325 primer sequence (in bold).
For both Sanger and 454 sequencing standard 25 μl PCRs were conducted with approximately 25 ng of genomic DNA. Sanger sequencing PCRs: 1 x Taq buffer (Promega), 2.5 mM MgCl2, 0.2 mM of each dNTP, 0.5 units platinum Taq polymerase (Invitrogen), 0.02 mg/ml BSA (Sigma-Aldrich) and 0.4 mM of each primer. 454 sequencing PCRs: 1 x GoTaq Colourless Master Mix (Promega), 2 mM MgCl2, and 0.4 mM of each primer. A negative control was included in each PCR and a touch-down protocol was used for all amplifications. The thermal profile consisted of an initial 5 min denaturation at 95 °C, followed by 10 cycles denaturation at 95 °C for 30 s, annealing at 60–44 °C for 30 s, with the annealing temperature decreasing by 4 °C every 2 cycles, and extension at 72 °C for 90 s. This was followed by 35 cycles of denaturation at 95 °C for 45 s, annealing at 45 °C for 45 s, and extension at 72 °C for 90 s, followed by a final extension at 72 °C for 10 min. PCR products were visualised on 1.2 % agarose gels to confirm amplification.
Cloning of PCR products for Sanger sequencing was undertaken by ligating PCR product into the pCR®II-TOPO bacterial plasmid (TOPO TA Cloning® Kit Dual Promoter, Invitrogen) and transformed into TOP10F’ chemically competent cells following the manufacturer’s protocol for the TOPO TA Cloning® Kit Dual Promoter (Invitrogen). Recombinant clones were detected by blue/white screening and selected clones were suspended in 30 μl of ddH2O for a minimum of 1 h. Prior to sequencing, clones were screened for the expected insert size using 2 μl of bacterial water in a PCR containing M13 forward and reverse primers. For each individual, at least 20 amplified clones of the expected insert size were purified with 1 μl ExoSAP-IT (USB), and sequenced commercially (Macrogen, Korea) using M13 primers in Sanger sequencing. To validate each allelic sequence, DNA from each individual was amplified and cloned twice. Cloned sequences were retained if they occurred in at least two independent PCRs. The retained sequences were edited, assembled and aligned using Geneious v. 6 .
PCR products for 454 sequencing were purified using the Agencourt AMPure XP purification kit (Beckman Coulter) according to the manufacturer’s instructions. Purified products were pooled in equimolar concentrations and sequenced commercially (Macrogen, Korea) on an eighth plate on a 454 GS-FLX run.
Bioinformatics following 454 sequencing
Following 454 sequencing, we used the jMHC software  to extract sequences and assign reads to individuals. Only reads showing a complete match to the forward and reverse primers were retained, ensuring that all assigned reads covered the whole amplicon. At the level of the whole dataset, we retained only sequences of exactly 159 bp which occurred in at least two independent PCRs, each represented by at least three reads . Remaining sequences were assigned to individuals based on MID tags and aligned in Geneious v. 6 . At the individual level, we discarded any alleles with coverage lower than 10 % of the allele with the highest coverage within each individual in order to remove less reliable sequence variants . Despite the potential for excluding true alleles in individuals with very high coverage, we believe this conservative approach provides a sufficiently comprehensive assessment of the levels of diversity across species to assess patterns of MHC evolution in passerines.
Allelic diversity and phylogenetic analyses
We calculated intraspecific nucleotide and amino acid diversity using the Kimura 2-parameter and Poisson models respectively in MEGA v.5 and estimated standard errors through 1000 bootstrap replicates . We reconstructed the evolutionary relationships among MHC class II β exon 2 sequences from a) the basal passerine species in this study, and b) a wider range of avian taxa, including species from Accipitridae, Apterygidae, Ardeidae, Galliformes, Procellariformes, Spheniscidae, Strigidae, and non-Australian Passeriformes (see Additional file 6). We estimated the best-fit model of evolutionary change using jModelTest2 ([GTR + Γ; [64, 65]) based on Akaike’s Information Criterion and constructed a Bayesian phylogeny in MRBAYES v. 3.2  with a Nile crocodile sequence (Crocodylus niloticus FJ886734) as an outgroup. Trees were sampled every 1,000 generations over 50 million generations of two simultaneous runs, with one cold and three heated Markov Coupled Monte-Carlo chains. The first 25 % of trees were discarded as burn-in and the remaining used to construct a consensus tree which was visualised in FigTree v. 1.4 . We evaluated the monophyly of the passerine sequences in our dataset by comparing estimates of marginal likelihood in natural log units for positively and negatively constrained topologies  in MRBAYES v. 3.2. We performed topology testing, using a Bayesian stepping stone approach for 5 million generations and assessed support for the constrained topology over the negative constraint using Bayes Factors . To assess phylogenetic relationships based on putatively neutral and adaptive genetic variation, we constructed Neighbour-Joining trees in MEGA v. 5  using the Nei-Gojobori method with Jukes-Cantor correction based on a) non-synonymous (dN) and b) synonymous substitutions (dS). Bootstrap tests of trees were conducted using 5,000 replicates. As exon 2 of the MHC class II β genes has been demonstrated to undergo high rates of recombination and gene conversion [2, 70, 71], we complementarily utilised phylogenetic networks in addition to phylogenetic trees to visualise relationships among MHC sequences. We used the neighbour-net algorithm based on Jukes-Cantor distances in Splitstree v. 4  to examine the relationships among MHC sequences from the seven basal passerine species.
Samples were collected under permits from the Victorian Department of Environment and Primary Industries (numbers 10004294 under the Wildlife Act 1975 and the National Parks Act 1975, and NWF10455 under section 52 of the Forest Act 1958), the Australian Bird and Bat Banding Schemes and under approval and monitoring of Monash University ethics processes (BSCI/2007/07).
Availability of supporting data
The sequence data supporting the results of this article are available in the Figshare digital repository and can be accessed at doi:10.4225/49/571353BCF0125.
We thank all the members of the Birds Linkage team, especially A. Lill and N. Takeuchi, and numerous volunteers for assistance with fieldwork.
This work was supported by the Australian Research Council (LP0776322), Department of Environment and Primary Industries, Museum Victoria, Parks Victoria, the North Central and Goulburn Broken Catchment Management Authorities, University of Melbourne (SB) and the Australian Academy of Sciences (JM).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Klein J. Natural history of the major histocompatibility complex. New York: Wiley; 1986.Google Scholar
- Hess CM, Edwards SV. The Evolution of the Major Histocompatibility Complex in Birds. BioScience. 2002;52:423.View ArticleGoogle Scholar
- Hughes AL, Yeager M. Natural selection at major histocompatibility complex loci of vertebrates. Annu Rev Genet. 1998;32:415–35.View ArticlePubMedGoogle Scholar
- Garrigan D, Hedrick PW. Perspective: detecting adaptive molecular polymorphism: lessons from the MHC. Evol Int J Org Evol. 2003;57:1707–22.View ArticleGoogle Scholar
- Klein J, Ono H, Klein D, O’hUigin C. The Accordion Model of Mhc Evolution. In: Gergely J, Benczúr M, Erdei A, Falus A, Füst G, Medgyesi G, editors. Prog. Immunol, vol. VIII. Berlin Heidelberg: Springer; 1993. p. 137–43.Google Scholar
- Kaufman J, Völk H, Wallny HJ. A “minimal essential Mhc” and an “unrecognized Mhc”: two extremes in selection for polymorphism. Immunol Rev. 1995;143:63–88.View ArticlePubMedGoogle Scholar
- Kaufman J, Milne S, Göbel TW, Walker BA, Jacob JP, Auffray C, et al. The chicken B locus is a minimal essential major histocompatibility complex. Nature. 1999;401:923–5.View ArticlePubMedGoogle Scholar
- Shaw I, Powell TJ, Marston DA, Baker K, van Hateren A, Riegert P, et al. Different evolutionary histories of the two classical class I genes BF1 and BF2 illustrate drift and selection within the stable MHC haplotypes of chickens. J Immunol. 2007;178:5744–52.View ArticlePubMedGoogle Scholar
- Skjødt K, Koch C, Crone M, Simonsen M. Analysis of chickens for recombination within the MHC (B-complex). Tissue Antigens. 1985;25:278–82.View ArticlePubMedGoogle Scholar
- Plachy J, Pink JR, Hála K. Biology of the chicken MHC (B complex). Crit Rev Immunol. 1992;12:47–79.PubMedGoogle Scholar
- Schat K. Immunity in Marek’s disease and other tumors. In: Toivanen A, Toivanen P, editors. Avian Immunol. Basis Pract. Boca Raton: CRC Press; 1987.Google Scholar
- Wittzell H, Bernot A, Auffray C, Zoorob R. Concerted evolution of two Mhc class II B loci in pheasants and domestic chickens. Mol Biol Evol. 1999;16:479–90.View ArticlePubMedGoogle Scholar
- Promerová M, Králová T, Bryjová A, Albrecht T, Bryja J. MHC class IIB exon 2 polymorphism in the grey partridge (Perdix perdix) is shaped by selection, recombination and gene conversion. PLoS ONE. 2013;8:e69135.View ArticlePubMedPubMed CentralGoogle Scholar
- Hosomichi K, Shiina T, Suzuki S, Tanaka M, Shimizu S, Iwamoto S, et al. The major histocompatibility complex (Mhc) class IIB region has greater genomic structural flexibility and diversity in the quail than the chicken. BMC Genomics. 2006;7:322.View ArticlePubMedPubMed CentralGoogle Scholar
- Kikkawa EF, Tsuda TT, Naruse TK, Sumiyama D, Fukuda M, Kurita M, et al. Analysis of the sequence variations in the Mhc DRB1-like gene of the endangered Humboldt penguin (Spheniscus humboldti). Immunogenetics. 2005;57:99–107.View ArticlePubMedGoogle Scholar
- Alcaide M, Edwards SV, Negro JJ. Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey. J Mol Evol. 2007;65:541–54.View ArticlePubMedGoogle Scholar
- Silva MC, Edwards SV. Structure and Evolution of a New Avian MHC Class II B Gene in a Sub-Antarctic Seabird, the Thin-Billed Prion (Procellariiformes: Pachyptila belcheri). J Mol Evol. 2009;68:279–91.View ArticlePubMedGoogle Scholar
- Ekblom R, Grahn M, Hoglund J. Patterns of polymorphism in the MHC class II of a non-passerine bird, the great snipe (Gallinago media). Immunogenetics. 2003;54:734–41.PubMedGoogle Scholar
- Alcaide M, Muñoz J, Martínez-de la Puente J, Soriguer R, Figuerola J. Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae). Ecol. Evol. 2014;4:688–98.Google Scholar
- Kohyama TI, Akiyama T, Nishida C, Takami K, Onuma M, Momose K, et al. Isolation and characterization of major histocompatibility complex class II B genes in cranes. Immunogenetics. 2015;67:705–10.View ArticlePubMedGoogle Scholar
- Balakrishnan CN, Ekblom R, Völker M, Westerdahl H, Godinez R, Kotkiewicz H, et al. Gene duplication and fragmentation in the zebra finch major histocompatibility complex. BMC Biol. 2010;8:29.View ArticlePubMedPubMed CentralGoogle Scholar
- Bollmer JL, Dunn PO, Whittingham LA, Wimpee C. Extensive MHC class II B gene duplication in a passerine, the common Yellowthroat (Geothlypis trichas). J Hered. 2010;101:448–60.View ArticlePubMedGoogle Scholar
- Sato A, Tichy H, Grant PR, Grant BR, Sato T, O’hUigin C. Spectrum of MHC class II variability in Darwin’s finches and their close relatives. Mol Biol Evol. 2011;28:1943–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Promerová M, Albrecht T, Bryja J. Extremely high MHC class I variation in a population of a long-distance migrant, the Scarlet Rosefinch (Carpodacus erythrinus). Immunogenetics. 2009;61:451–61.View ArticlePubMedGoogle Scholar
- Sepil I, Moghadam HK, Huchard E, Sheldon BC. Characterization and 454 pyrosequencing of Major Histocompatibility Complex class I genes in the great tit reveal complexity in a passerine system. BMC Evol Biol. 2012;12:68.View ArticlePubMedPubMed CentralGoogle Scholar
- Strand T, Westerdahl H, Höglund J, V Alatalo R, Siitari H. The Mhc class II of the Black grouse (Tetrao tetrix) consists of low numbers of B and Y genes with variable diversity and expression. Immunogenetics. 2007;59:725–34.View ArticlePubMedGoogle Scholar
- Miller HC, Lambert DM. Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae). Immunogenetics. 2004;56:178–91.PubMedGoogle Scholar
- Burri R, Hirzel HN, Salamin N, Roulin A, Fumagalli L. Evolutionary patterns of MHC class II B in owls and their implications for the understanding of avian MHC evolution. Mol Biol Evol. 2008;25:1180–91.View ArticlePubMedGoogle Scholar
- Burri R, Salamin N, Studer RA, Roulin A, Fumagalli L. Adaptive divergence of ancient gene duplicates in the avian MHC class II β. Mol Biol Evol. 2010;27:2360–74.View ArticlePubMedGoogle Scholar
- Kelley J, Walter L, Trowsdale J. Comparative genomics of major histocompatibility complexes. Immunogenetics. 2005;56:683–95.View ArticlePubMedGoogle Scholar
- Alcaide M, Liu M, Edwards SV. Major histocompatibility complex class I evolution in songbirds: universal primers, rapid evolution and base compositional shifts in exon 3. PeerJ. 2013;1:e86.View ArticlePubMedPubMed CentralGoogle Scholar
- Sutton JT, Robertson BC, Grueber CE, Stanton J-AL, Jamieson IG. Characterization of MHC class II B polymorphism in bottlenecked New Zealand saddlebacks reveals low levels of genetic diversity. Immunogenetics. 2013;65:619–33.View ArticlePubMedGoogle Scholar
- Edwards SV, Wakeland EK, Potts WK. Contrasting histories of avian and mammalian Mhc genes revealed by class II B sequences from songbirds. Proc Natl Acad Sci U S A. 1995;92:12200–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Sibley C, Ahlquist J. Phylogeny and classification of the birds of the world. New Haven: Yale University Press; 1990.Google Scholar
- Babik W, Taberlet P, Ejsmond MJA, Radwan J. New generation sequencers as a tool for genotyping of highly polymorphic multilocus MHC system. Mol Ecol Resour. 2009;9:713–9.View ArticlePubMedGoogle Scholar
- Galan M, Guivier E, Caraux G, Charbonnel N, Cosson J-F. A 454 multiplex sequencing method for rapid and reliable genotyping of highly polymorphic genes in large-scale studies. BMC Genomics. 2010;11:296.View ArticlePubMedPubMed CentralGoogle Scholar
- Moore MJ, Dhingra A, Soltis PS, Shaw R, Farmerie WG, Folta KM, et al. Rapid and accurate pyrosequencing of angiosperm plastid genomes. BMC Plant Biol. 2006;6:17.View ArticlePubMedPubMed CentralGoogle Scholar
- Longeri M, Zanotti M, Damiani G. Recombinant DRB sequences produced by mismatch repair of heteroduplexes during cloning in Escherichia coli. Eur J Immunogenet. 2002;29:517–23.View ArticlePubMedGoogle Scholar
- Lenz TL, Becker S. Simple approach to reduce PCR artefact formation leads to reliable genotyping of MHC and other highly polymorphic loci–implications for evolutionary analysis. Gene. 2008;427:117–23.View ArticlePubMedGoogle Scholar
- Burri R, Promerová M, Goebel J, Fumagalli L. PCR-based isolation of multigene families: lessons from the avian MHC class IIB. Mol Ecol Resour. 2014;14:778–88.View ArticlePubMedGoogle Scholar
- Babik W, Durka W, Radwan J. Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber). Mol Ecol. 2005;14:4249–57.View ArticlePubMedGoogle Scholar
- Oppelt C, Wutzler R, von Holst D. Characterisation of MHC class II DRB genes in the northern tree shrew (Tupaia belangeri). Immunogenetics. 2010;62:613–22.View ArticlePubMedGoogle Scholar
- Bonneaud C, Sorci G, Morin V, Westerdahl H, Zoorob R, Wittzell H. Diversity of Mhc class I and IIB genes in house sparrows (Passer domesticus). Immunogenetics. 2004;55:855–65.View ArticlePubMedGoogle Scholar
- Anmarkrud JA, Johnsen A, Bachmann L, Lifjeld JT. Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes: MHC class II diversity in bluethroats. J Evol Biol. 2010;23:1206–17.View ArticlePubMedGoogle Scholar
- Huson DH, Bryant D. Application of Phylogenetic Networks in Evolutionary Studies. Mol Biol Evol. 2006;23:254–67.View ArticlePubMedGoogle Scholar
- Seddon JM, Ellegren H. A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution. Proc R Soc Lond B Biol Sci. 2004;271:2283–91.View ArticleGoogle Scholar
- Radwan J, Kawałko A, Wójcik JM, Babik W. MHC-DRB3 variation in a free-living population of the European bison, Bison bonasus. Mol Ecol. 2007;16:531–40.View ArticlePubMedGoogle Scholar
- Eimes JA, Townsend AK, Sepil I, Nishiumi I, Satta Y. Patterns of evolution of MHC class II genes of crows (Corvus) suggest trans-species polymorphism. PeerJ [Internet]. 2015;3. Available from: https://peerj.com/articles/853. Accessed 17 May 2016.
- Strandh M, Lannefors M, Bonadonna F, Westerdahl H. Characterization of MHC class I and II genes in a subantarctic seabird, the blue petrel, Halobaena caerulea (Procellariiformes). Immunogenetics. 2011;63:653–66.View ArticlePubMedGoogle Scholar
- Miller HC, Bowker-Wright G, Kharkrang M, Ramstad K. Characterisation of class II B MHC genes from a ratite bird, the little spotted kiwi (Apteryx owenii). Immunogenetics. 2011;63:223–33.View ArticlePubMedGoogle Scholar
- Klein J, Sato A, Nagl S, O’hUigín C. Molecular Trans-Species Polymorphism. Annu Rev Ecol Syst. 1998;29:1-C1.View ArticleGoogle Scholar
- Kriener K, O’hUigin C, Tichy H, Klein J. Convergent evolution of major histocompatibility complex molecules in humans and New World monkeys. Immunogenetics. 2000;51:169–78.View ArticlePubMedGoogle Scholar
- Takahata N, Satta Y. Selection, convergence, and intragenic recombination in HLA diversity. Genetica. 1998;102–103:157–69.View ArticlePubMedGoogle Scholar
- Barker FK, Cibois A, Schikler P, Feinstein J, Cracraft J. Phylogeny and diversification of the largest avian radiation. Proc Natl Acad Sci U S A. 2004;101:11040–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Amos JN, Harrisson KA, Radford JQ, White M, Newell G, Nally RM, et al. Species-and sex-specific connectivity effects of habitat fragmentation in a suite of woodland birds. Ecology. 2014;95:1556–68.View ArticlePubMedGoogle Scholar
- Brown J, Jardetzky T, Gorga J, Stern L, Urban R, Strominger J, et al. Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1. Nature. 1993;364:33–9.View ArticlePubMedGoogle Scholar
- Sambrook J, Fritsch E, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, New York: Harbor Laboratory Press; 1989.Google Scholar
- Harrisson K, Pavlova A, Amos JN, Takeuchi N, Lill A, Radford JQ, et al. Fine-scale effects of habitat loss and fragmentation despite large-scale gene flow for some regionally declining woodland bird species. Landsc Ecol. 2012;27:813–27.View ArticleGoogle Scholar
- Edwards S, Grahn M, Potts W. Dynamics of Mhc evolution in birds and crocodilians: amplification of class II genes with degenerate primers. Mol Ecol. 1995;4:719–30.View ArticlePubMedGoogle Scholar
- Drummond A, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, et al. Geneious [Internet]. 2011. Available from: http://www.geneious.com/. Accessed 17 May 2016.Google Scholar
- Stuglik MT, Radwan J, Babik W. jMHC: software assistant for multilocus genotyping of gene families using next-generation amplicon sequencing. Mol Ecol Resour. 2011;11:739–42.View ArticlePubMedGoogle Scholar
- Stiebens VA, Merino SE, Chain FJJ, Eizaguirre C. Evolution of MHC class I genes in the endangered loggerhead sea turtle (Caretta caretta) revealed by 454 amplicon sequencing. BMC Evol Biol. 2013;13:95.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Darriba D, Taboada G, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.View ArticlePubMedPubMed CentralGoogle Scholar
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Rambaut A. FigTree [Internet]. 2009. Available from: http://tree.bio.ed.ac.uk/software/figtree. Accessed 17 May 2016.Google Scholar
- Xie W, Lewis PO, Fan Y, Kuo L, Chen M-H. Improving Marginal Likelihood Estimation for Bayesian Phylogenetic Model Selection. Syst Biol. 2011;60:150–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Kass RE, Raftery AE. Bayes Factors. J Am Stat Assoc. 1995;90:773.View ArticleGoogle Scholar
- Schaschl H, Wandeler P, Suchentrunk F, Obexer-Ruff G, Goodman SJ. Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates. Heredity. 2006;97:427–37.View ArticlePubMedGoogle Scholar
- Nei M, Gu X, Sitnikova T. Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci U S A. 1997;94:7799–806.View ArticlePubMedPubMed CentralGoogle Scholar
- Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491:444–8.View ArticlePubMedGoogle Scholar