Skip to main content

Evolutionary history of black grouse major histocompatibility complex class IIB genes revealed through single locus sequence-based genotyping



Gene duplications are frequently observed in the Major Histocompatibility Complex (MHC) of many species, and as a consequence loci belonging to the same MHC class are often too similar to tell apart. In birds, single locus genotyping of MHC genes has proven difficult due to concerted evolution homogenizing sequences at different loci. But studies on evolutionary history, mode of selection and heterozygosity correlations on the MHC cannot be performed before it is possible to analyse duplicated genes separately. In this study we investigate the architecture and evolution of the MHC class IIB genes in black grouse. We developed a sequence-based genotyping method for separate amplification of the two black grouse MHC class IIB genes BLB1 and BLB2. Based on this approach we are able to study differences in structure and selection between the two genes in black grouse and relate these results to the chicken MHC structure and organization.


Sequences were obtained from 12 individuals and separated into alleles using the software PHASE. We compared nucleotide diversity measures and employed selection tests for BLB1 and BLB2 to explore their modes of selection. Both BLB1 and BLB2 are transcribed and display classic characteristics of balancing selection as predicted for expressed MHC class IIB genes. We found evidence for both intra- and interlocus recombination or gene conversion, as well as indication for positive but differential selection at both loci. Moreover, the two loci appear to be linked. Phylogenetic analyses revealed orthology of the black grouse MHC class IIB genes to the respective BLB loci in chicken.


The results indicate that the duplication of the BLB gene occurred before the species divergence into black grouse, chicken and pheasant. Further, we conclude that BLB1 and BLB2 in black grouse are subjected to homogenizing concerted evolution due to interlocus genetic exchange after species divergence. The loci are in linkage disequilibrium, which is in line with the theory of tightly coevolving genes within the MHC under the minimal essential MHC hypothesis. Our results support the conclusion that MHC form and function in birds derived from studies on the domesticated chicken are not artefacts of the domestication process.


The duplication of genes is a common process in the evolution and adaptation of most organisms [1]. Understanding the molecular mechanisms that shaped duplicated genes is essential to study the evolutionary history and selection processes acting on complex gene families in a larger context. A very prominent example of a multigene family with frequent genetic and genomic duplications is the Major Histocompatibility Complex (MHC), which holds a key role in the vertebrate adaptive immune response in pathogen recognition and defense [2]. Thereby, MHC genes are linked to individual fitness and the outcome of infections (reviewed in e.g. [3]). The genes of the MHC are known to be the most polymorphic and among the best studied functional genes in vertebrates [4]. In humans, for instance, close to a thousand different alleles are known from the MHC HLA-DRB1 locus [5]. The extreme polymorphism commonly found at MHC genes is believed to be maintained by means of balancing selection, driven by selection through diverse pathogens and leading to maintenance of allelic variation within populations [68]. Balancing selection acts on variation at the MHC by means of heterozygote advantage (dominant as well as overdominant selection), rare-allele advantage (negative frequency-dependent selection) and/or selection that fluctuates over space and time [3]. Gene conversion and recombination [9], trans-species evolution [10], MHC-dependent sexual selection (reviewed in [4]) and selection against hitchhiking recessive deleterious mutations [11] have also been shown to play a role in shaping MHC diversity.

Duplications of MHC loci are frequently observed, and as a consequence loci belonging to the same MHC class are, in many species, too similar to distinguish by their allelic sequences [1216]. However, studies on the evolutionary history, the mode of selection and heterozygosity correlations on the MHC cannot be performed before it is possible to analyse duplicated loci separately [4]. Single locus amplification of polymorphic MHC genes in birds has been a goal for evolutionary biologists for over a decade [12, 17], but has not been achieved until recently. Very few bird species studied so far have only a single MHC class IIB locus (i.e. the green-rumped parrotlet [18], penguins [19] and kestrels [20]). Among the numerous birds with duplicated MHC class IIB loci, single locus amplification has so far only been reported in red jungle fowl/domestic chicken (Gallus gallus) [21, 22], captive turkey (Meleagris gallapavo) [23] and Barn owl (Tyto alba) [24].

The MHC shows pronounced differences in genomic organization and in number of MHC loci between vertebrate lineages [25], especially when mammalian and non-mammalian species are compared [26, 27]. Different hypotheses have emerged for how the complex MHC family has evolved. The prevailing consensus is that MHC evolution is characterized by repeated gene duplication (birth) and gene loss (death), whereby the loci evolve under a birth-and-death model [2, 28, 29]. Phylogenetic relationships can reveal whether gene duplication arose post speciation or pre speciation. Post speciation duplication leads to a pattern where paralogous MHC loci within a species are more closely related than orthologous MHC loci between species [30, 31]. Pre speciation duplication, in contrast, took place in an ancestral species and results in orthologous MHC loci from different species being more closely related than paralogous MHC loci within species. However, these patterns can likewise be generated by different processes. A situation identical to what to be expected under a post speciation duplication event can be caused by concerted evolution acting on early duplicated genes [26, 32]. Under concerted evolution, gene fragments are frequently exchanged between paralogous loci, and thereby homogenized (interlocus genetic exchange by recombination or gene conversion [9, 33]). However, if convergent selection acts on either different loci within a species or between different species, those loci will functionally converge and mask the actual evolutionary relationship [32, 34]. These different processes are not mutually exclusive, which adds a further level of difficulty to the interpretation of evolutionary patterns at the MHC.

The domestic chicken (Gallus gallus domesticus) is the major avian MHC model species with a history of more than 60 years of MHC studies [35]. The core of MHC in the chicken is called the classical MHC or the BF/BL region which, among others, consists of MHC class I and class II genes, TAP genes, Tapasin and one MHC class IV (BG) gene. This gives a total of 19 genes in 92 kb [36] situated on a single chromosome compared to the human MHC region which comprises 128 expressed genes as well as pseudo genes on a stretch of 3600 kb on several chromosomes [37]. The chicken MHC also appears more compact than what has been observed in the passerine bird which is best known, the zebra finch (Taeniopygia guttata), in which the MHC region is more complex and found on at least four chromosomes [38]. Because of the compact and small structure of the chicken MHC compared to the mammalian MHC, and the presence of few expressed genes, it has been named “a minimal essential MHC”. In addition, the chicken MHC is arranged differently than the mammalian MHC [36]. MHC class I and class II genes are tightly linked [31, 39], and TAP genes has been demonstrated to co-evolve with MHC class I (BF) genes [31]. It has also been suggested that other MHC genes may co-evolve [30, 40]. Based on experiments, recombination within BF/BL in the chicken is rarely observed [41] but sequencing showed evidence of gene conversion and recombination in shaping the chicken MHC [39]. Furthermore, examinations of BF/BL in pedigreed families in the closely related galliforms turkey and Japanese quail (Coturnix japonica) have revealed direct evidence of gene conversion and recombination [42, 43].

In the chicken, both MHC class I (BF) and MHC class IIB (BLB) comprise two loci but only one locus in each class is highly expressed. While the BLB2 locus is considered to be dominantly expressed [44, 45], the BLB1 locus has been suggested not be involved in peptide binding and to be neutral to selection [30, 46]. Nevertheless, interlocus genetic exchange has been observed between BLB1 and BLB2 in a farmed population of pheasant (Phasianus colchicus) and in the domestic chicken, i.e. a case of concerted evolution [47]. In this study we test whether BLB1 and BLB2 evolved independently or in a concerted way in a wild close relative to both the pheasant and the chicken: the black grouse (Tetrao tetrix). Disentangling the evolutionary history of BLB1 and BLB2 in black grouse will add valuable information about how these genes evolved in the chicken and related species. To date, it is largely unknown if the uniquely compact MHC in the chicken [36] is a result of the long domestication process or was present before domestication.

Domestic chicken and black grouse belong to the same avian family (Phasianidae; Galliformes) and both possess two MHC class II B (BLB) loci that surround the Tapasin gene (Figure 1) [44, 48]. Although the MHC class II B region is well studied in the domestic chicken, no primers could be developed for independent amplification of either MHC class II B minor (BLB1) or MHC class II B major (BLB2) [49]. The chicken BLB1 and BLB2 are nearly identical in sequence and can only be amplified separately with PCR primers anchored outside the BLB genes along with a nested non-locus-specific PCR [22, 44].

Figure 1

Schematic figure of the genomic focus region of the MHC. Genes are illustrated by white arrows indicating their orientation, with the names above the drawing. The exon-intron structure is given for BLB1 and BLB2 to show the target region and primer positions. The positions and amplification directions of all primers used in this study are indicated by small one-headed arrows. Below the upper part of the drawing, the primers for the long-range PCR amplifications of BLB1 and BLB2 (C275, NBG262 and preBLB2F, BLBex3R) are given. Further down, the target region Exon 2 is enlarged, with the positions of the nested PCR primers (BLB1: RNAR1a, RNAF1a, NBG262; BLB2: Vorinex, RNAF1a, RNAR1a, Postex) given below. Note that BLB1 and BLB2 are orientated in opposite directions.

In this study, we present an approach to genotype the two MHC class II B loci separately in a wild galliform, the black grouse. The procedure comprises a long-range anchored PCR method followed by nested PCR reactions and sequence-based genotyping. Based on this approach we first aim to characterize and contrast MHC class II B diversity and mode of selection at the BLB1 and BLB2 loci in the black grouse. Second, we explore whether the two loci underlie concerted evolution or evolve independently. And third, we relate and compare the results to the MHC class II B structure and organization in the chicken.



We studied locus-specific MHC class IIB variation in twelve black grouse individuals, eleven of which were genotyped non locus-specific in our previous studies [15, 50]. We included one more individual (JHGO 213) and genomic DNA was extracted from tissue with a salt extraction protocol [51]. In addition, we used sequence data obtained in a fosmid sequencing study on the same individual JHGO 213 [48] [GenBank JQ028669]. Sequences from other species used in the data analysis derived from GenBank, accession numbers are given in the respective figure legends.

Locus-specific amplification of BLB1 and BLB2

In order to amplify the two black grouse MHC class IIB genes BLB1 and BLB2 separately, we applied a long-range anchored PCR method combined with nested non locus-specific PCRs. For each BLB locus, we anchored one of the PCR primers outside of the respective BLB gene for a locus-specific long-range PCR (Figure 1). In a second step, these long-range PCR fragments were used as template in nested PCRs with non locus-specific primers to amplify parts of the BLB exon 2.

We selectively amplified a long fragment of the BLB1 gene with a modification of primers designed for the chicken BLB minor[44]. We used the forward primer C275 located in the Blec1 gene adjacent to BLB1 gene, and the reverse primer NBG262, slightly modified from C262 situated within BLB1 exon 2 (Figure 1) (see Additional file 1 for primer sequences). These primers yielded a 1932 bp PCR product. All attempts to design primers amplifying the complete BLB1 exon 2 based on chicken genome sequences or on black grouse MHC fosmid data failed. In the 20 μl BLB1 long-range reaction, approximately 100 ng DNA, 0.3 μM of each primer C275 and NBG262, 0.8 μM dNTP, 0.6 μl DMSO, 1x Phusion GC buffer and 1U Phusion High-Fidelity DNA polymerase (Finnzymes, Espoo, Finland) were used. The PCR programme was initiated with 40 s at 98°C, followed by 35 cycles of 10 s at 98°C, 20 s at 62°C and 68 s at 72°C and ended with 6 min at 72°C.

Based on the black grouse MHC fosmid data [GenBank JQ028669] [48] we designed primers for a selective amplification of the BLB2 gene. We imported the BLB2 sequence and flanking regions to the NCBI web based primer design program Primer BLAST [52] and marked areas that were conserved between the black grouse BLB2, the chicken BLB major and the turkey BLB2 (Meleagris gallopavo), but were different from BLB1. The resulting primer pair forward preBLB2F and reverse BLBex3R yielded a 1359 bp PCR product ranging from the intragenic region between the Tapasin gene and BLB2 to the BLB2 exon 3 (Figure 1) (see Additional file 1 for primer sequences). For amplification of the BLB2 we used approximately 100 ng DNA, 0.3 μM of each primer preBLB2F and BLBex3R, 0.3 μM KAPA dNTP, 1x KAPA GC buffer and 0.5 U KAPA HiFi DNA polymerase (KAPA BIOSYSTEMS, Boston, United States) in a 25 μl reaction. The PCR programme was initiated with 5 min at 98°C, followed by 35 cycles of 20 s at 98°C, 15 s at 62°C and 90 s at 72°C and ended with a final extension for 5 min at 72°C.

Long-range PCR products for BLB1 and BLB2 revealed each one strong single band on a 1.5% agarose gel. A sterile toothpick was inserted into each band, taking up a small amount of the PCR product, and transferred into 50 μl ddH20. This solution was stored overnight at -20°C before using it in the nested PCRs as template as described in the following.

Nested PCR to amplify exon 2 and sequencing

The primers RNAF1a and RNAR1a [15] (Additional file 1) were used in a nested PCR to amplify 125 bp on the BLB1 and BLB2 long-range PCR products as template. 1 μl of the template solution, 0.48 μM of each primer, 0.6 mM of dNTP, 3 mM MgCl2, 1 x buffer and 0.75 U BioTaq DNA polymerase (DNA Technology, Aarhus, Denmark) were used in a 25 μl PCR reaction. The PCR programme was initiated at 94°C for 5 min, followed by 30 cycles at 94°C for 1 min and 30 s at both 64.9°C and 72°C, and ended for 10 min at 72°C. We also performed another nested PCR amplifying a longer part of the exon 2 on the BLB2 long-range PCR product only. A second primer pair, Vorinex2 and Postex2, developed for the closely related willow grouse (Lagopus lagopus), was used to amplify 251 bp of exon 2 (see Additional file 1). 1 μl of the template solution, 0.4 μM of each primer, 0.48 mM of dNTP, 1 x KCl buffer and 1 U Taq DNA polymerase (Fermentas, St. Leon-Rot, Germany) were used in a 25 μl PCR reaction. The programme was initiated at 94°C for 5 min, followed by 30 cycles at 94°C for 1 min and 30 s at both 66°C and 72°C before a final extension for 10 min at 72°C. Using the primer pair Vorinex2/Postex2 a longer portion of BLB2 exon 2 sequences were amplified that were not known from before, so every new 251 bp sequence was subjected to verifying PCR runs (both BLB2 long-range PCR and nested PCR) and sequencing. For BLB1, the 251 bp nested fragment could not be amplified due to the position of the nested primer Vorinex2 outside of the BLB1 long-range PCR product.

The nested PCR products were cleaned with an Exo-Sap reaction (Exonuclease I-Shrimp Alkaline Phospate, Fermentas, St. Leon-Rot, Germany). Cycle sequencing reactions were performed in both directions and the sequencing products were subjected to a post-reaction clean-up (GE Healthcare, Uppsala, Sweden), following the protocols of the manufacturer. All direct sequencing was performed both forward and reverse on a MegaBACE 1000 DNA analysing system (GE Healthcare, Uppsala, Sweden).

Confirmation of MHC alleles

PCR artefacts are a large issue in MHC studies and cautions to reduce the formation of artificial alleles are important to address [53]. To decrease the probability of PCR artefacts we only regarded alleles present in two independent PCR reactions as confirmed [54]. For the long-range PCR reactions for BLB1 and BLB2 we used different enzymes (Phusion High-Fidelity DNA polymerase and KAPA HiFi DNA), which could produce different levels of PCR artefacts. We were able to compare the sequences derived by the long-range PCR presented here, with previously cloned sequences using the non locus-specific primers RNAF1a and RNAR1a directly on genomic DNA [15, 50]. Between 16 and 26 clones were sequenced per individual in the previous studies. We found that enzymes in the long-range PCRs of BLB1 and BLB2 worked equivalently well in amplifying alleles corresponding to those previously found (Table 1).

Table 1 Allele designations for the twelve black grouse individuals included in this study

Data analysis

Identification of sequences and analysis of sequence variation

CodonCode Aligner version 3.7.1 was used for sequence editing and performing ambiguity codes for heterozygous positions. Allelic phase for heterozygote sequences was determined by computational inference with the PHASE haplotype reconstruction algorithm in DnaSP version 5.10.05 [55]. The PHASE algorithm has been proven reliable for determining the allelic phase in a number of previous MHC studies (e.g. [56, 57]). The parameters used were 1000 iterations, thinning 10, 100 burn in, and a recombination model. For the haplotype reconstruction of the 125 bp fragment we included sequences that were previously verified by cloning. Sequence variation statistics were calculated in DnaSP version 5.10.05 [55] including the number of segregating sites, nucleotide diversity (π) and average number of nucleotide differences (theta k).

Recombination, gene conversion and linkage

The power of any method for the detection of recombination or gene conversion strongly depends on the number of sequences included, the recombination rate and the number of sites differentiating the recombinant sequences [58, 59]. Therefore it is suggested to apply an array of different methods to analyse recombination. We used a set of seven methods to detect recombination signals and putative recombinant sequences. The minimum number of recombination events (Rm) according to the four-gamete test by Hudson & Kaplan [60] was calculated in DnaSP version 5.10.05 [55]. This method tests for overall evidence of recombination in the alignment. We used the RDP3 package [61] to apply the methods RDP [62], Maxchi [63], Chimaera [64] and GENECONV [65]. These methods are designed to detect recombination breakpoint locations; the first three use a dynamic sliding window, whereas GENECONV searches for unusual long regions of identity between sequences. Moreover, we used the method GARD [66] implemented in Datamonkey [67], which likewise identifies recombination breakpoints, searching all possible partitions in a probabilistic way. Recombination rate ρ and mutation rate θ were calculated using LDhat recombination rate scan [68].

We tested for linkage between the BLB1 and the BLB2 locus applying a likelihood ratio test of linkage disequilibrium implemented in Arlequin 3.5 [69]. In this test, the likelihood of the sample evaluated under the hypothesis of no association between loci (linkage equilibrium) is compared to the likelihood of the sample when association is allowed (linkage disequilibrium).

Analysis of selection

Positive selection is the spread to fixation of an allele that increases the fitness of individuals [70]. Initially, a locus undergoing balancing selection seems to be subjected to positive selection, but then experiences negative selection on alleles that became too frequent (under frequency dependent selection). Testing for long-term balancing selection in genetic regions therefore often includes testing for positive selection, as inferred from the ratio of non-synonymous/synonymous (dN/dS or ω) substitution rates, for instance [71]. Therefore, the outcome of methods can reflect a positive selection process, although in our context, testing for positive selection within a species is translated into balancing selection.

Tajima’s D, Fu and Li’s D and F tests were analysed in DnaSP version 5.10.05 [55] and Fisher exact tests were employed to test for differences from neutrality. In addition to the averaged Tajima’s D, we performed sliding window Tajima’s D with a window size of 11 bp and a step size of 1 bp for the 125 bp fragments of the black grouse BLB1 and BLB2. A negative Tajima’s D (as Fu and Li’s D, Fu and Li’s F) is a sign of negative selection or population expansion and a positive value is a sign of positive/balancing selection or a population bottleneck. Averaged synonymous (dS) and non-synonymous (dN) substitutions per synonymous and non-synonymous site were calculated in MEGA 5.05 [72]. The Nei-Gojobori method with Jukes Cantor corrections and 5000 bootstrap replicates was used to calculate the overall average of dN/dS (ω) for all sites, the peptide binding region (PBR) and non-PBR according to [73]. Deviations from neutrality (dN = dS) were assessed by a Z-test. Differences in dN and dS between PBR and non-PBR were tested for significance using Mann–Whitney U-tests.

A different approach to detect molecular evidence of positive selection is to calculate ω per codon using the maximum likelihood method CODEML implemented in PAML version 4.6 [74]. The programme estimates heterogeneous ω among sites applying different models of codon evolution. We compared M2a (allowing ω to vary between ω <1 (conserved), ω = 0 (neutral) and ω > 1 (positive)) with M1a (ω allowed to vary between ω <1 and ω = 0) and M8 (ω can vary in beta distribution of 0 and 1 including ω > 1) with M7 (ω allowed to vary in a beta distribution of 0 and 1) using a log-likelihood test. Significant positive selection was inferred if twice the difference in log-likelihood values between the two models was greater than the χ2 critical value for the given degrees of freedom. The Bayes empirical Bayes (BEB) approach was used to identify significantly positively selected codon sites. Tree files used in PAML analyses were generated using a maximum likelihood approach in PhyML3.0 [75], under the F81 model of nucleotide substitution and estimated gamma shape parameter. Models of nucleotide substitution and the distribution of rate variation across nucleotide sites (gamma) were estimated in jModelTest 0.1.1 [76]. The best molecular evolution models were selected by the Aikaike information criterion (AIC) (see Additional file 2).

Reconstruction of phylogenetic relationships

We constructed phylogenetic networks based on the BLB exon 2 sequences and on the 3rd codon positions of exon 2 only to illustrate the phylogenetic relationship among the black grouse BLB1 and BLB2 alleles and possible orthology to the chicken BLB1 and BLB2. The networks were built with the software SplitsTree4, using the neighbour net method [77]. Phylogenetic trees based on exon 3 sequences and on 3′ untranslated regions (UTRs) were constructed using the Neighbour-joining method with bootstrapping (1000 replicates) implemented in MEGA 5.0 [72]. For BLB orthology between the black grouse and the chicken, we particularly looked at the 3′UTRs, since orthology between the pheasant and the chicken BLB1 and BLB2 has previously been demonstrated with the 3′UTRs [47].


Locus-specific genotyping of black grouse BLB

We successfully amplified locus-specific long-range PCR products for both BLB1 and BLB2 in all twelve black grouse individuals. In subsequent nested PCR reactions, we amplified and sequenced a 125 bp PBR-rich fragment of the exon 2 for both loci, and almost the whole exon 2 for the BLB2 locus (251 out of 270 bp). Positions and sequences of the different primers are given in Figure 1 and Additional file 1. All heterozygous sequences could be assigned to haplotypes using the software PHASE implemented in DnaSP [55]. At the BLB1 locus we detected 6 unique 125 bp alleles (BLB1125, Table 1, Figure 2) and at the BLB2 locus 7 unique 125 bp alleles (BLB2125) and 9 unique 251 bp alleles (BLB2251). None of these sequences contained an indel or a stop codon, we thus assume that they are functional.

Figure 2

Alignment of amino acid sequences. PBR positions from Tong et al. [73] are marked with a +. The shaded amino acids are the PAML derived positively selected codon positions, blue for BLB2 and pink for BLB1.

At the BLB1125 locus seven out of twelve individuals were heterozygous, whereas at the BLB2 all but one individual (BLB2125), respective all individuals (BLB2251) were heterozygous (Table 1). For allele BLB2125-01 the 251 bp fragment revealed additional variation towards the 5′-end of exon 2. This resulted in a subdivision into BLB2251-01A, BLB2251-01B and BLB2251-01C, each differing by one amino acid exchange (Figure 2). We found one 125 bp-allele that occurred at both BLB loci, namely BLB125-04. There is indication, however, that BLB125-04 differs at the 5′-end of exon 2 between the two loci [see Additional file 3, BLB2251-04B and BLB202-04-cDNA, the latter from our previous study [15].

All individuals displayed unique MHC class II B genotypes. Combining the data for BLB1 and BLB2, 12 individuals had a minimum of 14 unique BLB alleles and each individual carried between three and four alleles. In comparison, 11 out of the 12 individuals previously analysed by non locus-specific BLB cloning [15, 50] revealed 11 unique BLB alleles and individuals carried between two and four alleles (Table 1). This was a minor underestimation of both the number of unique alleles and the number of alleles per individual due to the applied method compared to the here presented more accurate single locus amplifications.

There is evidence that both BLB loci are expressed in the black grouse. In an earlier study we found four BLB sequences in one individual by non locus-specific cloning based on cDNA (D870, Table 1, [15]). Our current study reveals that two of these sequences stem from BLB1 and the other two from BLB2.

Sequence analyses: gene diversity, recombination and linkage

Although the locus BLB1125 showed a lower number of unique alleles than BLB2125, it revealed more segregating sites, a higher number of average pairwise differences theta k and a higher nucleotide diversity π (Table 2). As expected, the longer fragment BLB2251 revealed more alleles, more segregating sites, a higher theta k but a lower π than the shorter BLB2125. For BLB1125, 7 out of 8 positions involved in antigen binding [73] were polymorphic, whereas 6 out of 8 and 12 out of 15 positions were variable for BLB2125 and BLB2251, respectively (Figure 2).

Table 2 Gene diversity measures

We used a set of different methods to detect recombination signals and putative recombinant sequences, and found evidence for recombination in each of the fragments (BLB1125, BLB2125, BLB2251, BLB1&2125), at least with the four-gamete test and the method Maxchi (Table 3). These recombination signals lead to a number of recombinant sequences including both intra- and interlocus recombination or gene conversion (Table 4). Considering both loci together (BLB1&2125) the estimated population recombination rate ρ was 11.42, and the mutation rate θ was 9.27, leading to a recombination-mutation ratio of 1.23. A ρ/θ ratio of > 1 is an indication for recombination being more prevalent in the dataset than point mutations [78]. For the loci separately, however, the estimated mutation rates were much higher than ρ, leading to ρ/θ ratios < 1. The results of the different tests combined strongly suggest that recombination or gene conversion has occurred in the black grouse MHC class II.

Table 3 Number of recombination events calculated by different methods
Table 4 Recombinant sequences and their potential parental sequences at the black grouse BLB1 and BLB2 (calculated with Maxchi)

We observed significant linkage disequilibrium between the BLB1 and BLB2 (likelihood-ratio-test, LnLikelihood LD: -46.087, LnLikelihood LE: -62.120, exact p = 0.0239), indicating that the two loci are linked.

Analyses of selection

Positive values for Tajima’s D, Fu & Li’s D and Fu & Li’s F were detected for BLB2125 and BLB2251 (Table 2), which are an indication for positive or balancing selection acting on the BLB2 locus. On the other hand, we found these values negative for BLB1125 and the combination of BLB1&2125, which is generally a sign of negative or purifying selection. None of these values were significant, however, so that they have to be interpreted as indication rather than evidence. A sliding window analysis of Tajima’s D along the BLB1125 and BLB2125 further disentangled this pattern of differences in selection at specific sites between the two loci (Figure 3).

Figure 3

Sliding window Tajimas’ D, for the 125 bp exon 2 fragments of the black grouse BLB1 and BLB2 (window size 11 bp, step size 1 bp). The threshold for P < 0.05 is shown by the dotted line.

We calculated the relative rates of non-synonymous (d N ) and synonymous (d S ) substitutions for all sites, the PBR and the non-PBR according to Tong et al. (2006) (Table 5). For the PBR, the values for the ratio d N /d S were considerably larger than 1 in all fragments, which was significant in the case of BLB2251. This is considered to be evidence for positive selection acting on the PBR. Moreover, both d N and d S were significantly higher at the PBR sites compared to the non-PBR sites for all fragments (Mann–Whitney U-test, all p < 0.01), a pattern which has recently been explained to be created by gene conversion combined with positive selection [78]. Synonymous mutations hitchhike alongside beneficial non-synonymous mutations and thereby occur more frequent than under neutral expectations.

Table 5 Relative rates of non-synonymous ( d N ) and synonymous ( d S ) substitutions with standard errors calculated for the two MHC class IIB loci, averaged over all sites, the peptide binding region (PBR) and non-PBR according to Tong et al.[73]

Positive selection on specific codon sites was detected using the maximum likelihood method CODEML implemented in PAML4.6 [74]. Two pairs of models were applied: M2a versus M1a, and M8 versus M7. The models M2a and M8, which allow for positive selection, fitted the data significantly better than the neutral models for both BLB loci (Table 6 and Additional file 4). Several codons were identified as significantly positively selected (ω > 1, Table 6 and Figure 2). Four selected codons were identified with both M2a and M8 for BLB1125 and four and five codons were identified with M2a and M8, respectively, for BLB2125. Codons 71, 74 and 86 were identified as positively selected in both BLB1125 and BLB2125. In addition, position 67 was identified in BLB1125 and position 61 and 85 in BLB2125. The longer fragment BLB2251 confirmed the same positively selected codons as calculated for the fragment BLB2125, and revealed one additional position at the 5′-end of the molecule (position 9). All identified positively selected codon positions were identical or directly situated next to a peptide binding site as identified by Tong et al. [73].

Table 6 Summary table for CODEML

Reconstruction of the phylogenetic relationships

In order to take a closer look at the orthology of the black grouse BLB1 and BLB2 with the corresponding BLB loci in the chicken and pheasant, we constructed a phylogenetic tree of the BLB1 and BLB2 3′UT regions (Figure 4), which are untranslated DNA sequences containing regulatory sequences. It has been shown that the BLB loci in both chicken and pheasant can be clearly distinguished by this part of the gene [47], in terms of nucleotide divergence and length. We found that the 3′UTR’s of the three species cluster according to gene and not according to species, forming one distinct BLB1 cluster and one BLB2 cluster, and similar differences in length. This suggests orthology of the BLB1 gene between the black grouse, chicken and pheasant, and orthology of the BLB2 gene between the three species.

Figure 4

Neighbour joining tree for the 3′UTR. BLB1 and BLB2 sequences derived from black grouse (fosmid individual JHGO 213 [48], chicken [GenBank AB268588] and pheasant [GenBank AJ224349]. The scale bar represents substitutions per site.

A contrasting phylogenetic relationship is shown in the neighbour network of the exon 2 sequences of BLB1 and BLB2 (Figure 5a). Here, sequences cluster according to species and not according to gene. All black grouse sequences cluster together, with BLB1 and BLB2 sequences intermingling. The black grouse cluster is distinct from the sequences of the domestic chicken and red jungle fowl, which form a mixed cluster, again BLB1 and BLB2 intermingling. The same picture is apparent in a neighbour network based on the third codon positions of exon 2, at which variation is mostly synonymous and therefore likely neutral. Further, exon 3 sequences of BLB1 and BLB2 repeat the same clustering according to species (see Additional file 5).

Figure 5

Neighbour network for BLB1 and BLB2. Sequences of the black grouse (pink and blue), chicken (dark and light green) [GenBank AB426144, AB426150-51, AJ248577, M29763] and red jungle fowl (brown and orange) [GenBank AM489767 - AM489776] for the a) 125 nucleotide sequence of exon 2 and b) only 3rd codon positions in the 125 nucleotide sequence of exon 2. Long-eared Owl Asot-DAB1*02 [GenBank EF641225] was used as an outgroup.


The chicken MHC stands out with an unparalleled simple and compact architecture compared to mammals or other avian groups like the passerines [26, 36]. The tight linkage of MHC genes and reduced recombination rates were suggested to have resulted in a close co-evolution of genes within MHC haplotypes ([40] but see [39]), leading to strong disease-associations of specific MHC haplotypes in consequence. In the evolutionary history of the avian MHC it is still not resolved whether this unique MHC structure in the chicken is a domestication artefact or was present before domestication as a galliform feature.

We here present a study on diversity, selection patterns and the co-evolutionary history of two transcribed MHC class IIB genes in a wild galliform bird species, the black grouse. We found evidence for both intra- and interlocus recombination or gene conversion, as well as indication for positive selection on the PBR at both loci. However, we also detected differences in the selection between the two loci, as positive selection was indicated at the BLB2 locus and purifying selection at the BLB1 locus. We were able to relate our findings to the structure of the chicken MHC class II by resolving orthology to the chicken BLB loci.

Similar to the work of [22] on red jungle fowl, we amplified the two black grouse BLB loci separately in twelve individuals. Eleven of the individuals included in this study were cloned in our previous studies with non locus-specific primers that gave a 125 bp BLB exon 2 product [15, 50]. The resulting sequences showed that individuals carried between two and four BLB alleles, without locus assignment. This is a slight underestimation in number of alleles per individual compared to the here presented study, which revealed between three and four alleles per individual for BLB1 and BLB2 combined. Part of this additional variation was detected through the longer BLB2 fragment, which revealed further variable positions towards the 5′-end of the exon 2. To another part, the underestimation by the cloning/sequencing approach may be a result of too few analysed clones (16–26 per individual). On the other hand, for individual D248 the cloning procedure revealed one allele (BLB*125-04) which could not be confirmed by the single-locus amplification. This additional allele is likely to be a contamination during the cloning process. Considering that cloning is more time-consuming than direct sequencing, can add artefacts derived from mismatch repair or the formation of chimeric sequences [53] and is more contamination-prone, the here presented single locus sequence-based typing method clearly outmatches the cloning/sequencing approach in both accuracy and effort.

Another finding in this study concerns the expression of the two BLB loci in black grouse. We earlier amplified BLB cDNA sequences for individual D870 [15], which we could show now derived from both BLB1 and BLB2. Thus, both BLB1 and BLB2 are transcribed loci in black grouse. In domestic chicken BLB2 is considered to be dominantly expressed [44, 45], while BLB1 is less expressed and has even been suggested to be neutral to selection and not involved in peptide binding [30, 46]. This ultimately leads to the question whether the BLB1 and BLB2 loci in black grouse show a similar situation or differ in the degree and mode of selection compared to the chicken. We therefore compared diversity and selection patterns between BLB1 and BLB2 to explore their respective modes of selection.

We observed significant linkage between BLB1 and BLB2, in line with the hypothesis of Kaufman [30, 40] that MHC genes in chicken form a tightly linked cluster of co-evolving genes. Similarly, Agudo et al. [79] showed in a study on Egyptian vultures (Neophron percnopterus) the existence of linkage groups, containing pairs of MHC alleles in strong linkage disequilibrium. They interpret the presence of linkage groups with similar MHC alleles as indication of concerted evolution acting on the MHC gene duplicates. However, concerted evolution does not necessarily homogenize a gene cluster. Gene conversion can both homogenize and diversify among paralogues, as was shown for the gene family hsp70 in Drosophila[80].

Intuitively one might think that gene conversion and recombination should tend to break up linkage disequilibrium. But this does not seem to be this simple. The human olfactory receptor (OR) gene cluster is another gene family suggested to evolve under concerted evolution with gene conversion events [81]. Linkage disequilibrium is highly significant in the centromeric part of the gene cluster, whereas no linkage disequilibrium could be observed in the telomeric part of the cluster. These considerable differences suggest the presence of several recombination hotspots within the OR cluster.

In terms of number of different alleles, MHC diversity was similar at the two BLB loci, however, BLB2 revealed higher heterozygosity and was less variable in terms of nucleotide diversity than BLB1. Signals of positive selection were found on both the BLB1 and BLB2 locus. Elevated dN/dS ratios were obtained at the PBR of both loci, though only significant for the longer BLB2 fragment. Maximum likelihood analyses confirmed that models allowing for positive selection fitted our data significantly better than neutral models. With these models, four identical codons were identified to be under significant positive selection in BLB1 and BLB2, and two additional positions were identified for the BLB2 locus. All of these positively selected codon positions were either congruent or directly adjacent to a peptide binding site as identified by Tong et al. [73] for human MHC class II molecules. This congruence emphasises the structural similarity of codon positions involved in peptide binding throughout evolutionary lineages. However, it has been suggested the methods used here might be prone to overestimate positive selection, particularly when recombination rate is high [82].

A different pattern of selection was found, however, regarding neutrality tests as Tajima’s D. We found a negative Tajima’s D value for BLB1, indicative of purifying selection, and a positive Tajima’s D for BLB2 indicating positive or balancing selection, although not significant in both cases. In a sliding window analysis of Tajima’s D comparing BLB1 and BLB2 it became more apparent that these different selection footprints are due to differences at some specific nucleotide sites only, while at other positions positive selection between the two loci coincides. These opposite but vague signals of differential selection at BLB1 and BLB2 made us explore this further. Worley et al. [22] found likewise a negative Tajima’s D for BLB1 and positive value for BLB2 (both n.s.) in a captive population of red jungle fowl but the authors points out that a negative Tajima’s D could also be a result of a population bottleneck. As a further comparison, we obtained domestic chicken BLB1 and BLB2 sequences from GenBank and calculated Tajima’s D the same way as we did for the black grouse sequences. This analysis repeated the observation of positive Tajima’s D for BLB2 and negative for BLB1 (data not shown). We also observed that theta k and π were higher in BLB1 than BLB2 in the domestic chicken, which was likewise observed for black grouse. In conclusion, we interpret this as a repeated pattern of differential selection on BLB1 and BLB2 in galliform birds.

Detecting recombination or gene conversion events using statistical methods can be highly problematic, in particular when small fragments are transferred or the gene conversion rate is too high. Hence, we have to keep in mind that a lack of evidence is no evidence for absence. To minimize the risk of missing the footprint of gene conversion, it is recommended to use multiple statistical methods for detection [58]. Following this recommendation we applied a set of seven statistical methods, which have been evaluated to be powerful and accurate for different conditions and scenarios [59]. While two of the methods (RDP and GENECONV) failed to detect most recombination signals, the other methods detected recombination or gene conversion both within and 'between the two BLB loci. As an additional hint for the occurrence of gene conversion, we detected significantly higher d N and d S values at the PBR sites compared to the non-PBR sites at both loci. This situation is likely to be created by a combination of positive selection and gene conversion [78, 83]. Gene conversion events transferring advantageous non-synonymous substitutions at the PBR are positively selected, and in doing so synonymous substitutions are carried along. This way, synonymous substitutions appear more often at the PBR than expected by point mutations under neutrality. We conclude that both intra- and interlocus genetic exchange play an important role in shaping the black grouse MHC class II.

The 3′UTRs between BLB1 versus BLB2 differ in domestic chicken as well as between the corresponding DAB1 and DAB2 loci in pheasants in both length and nucleotide composition, so that the sequences cluster together as orthologous genes rather than according to species [47]. We have observed similar differences in the 3′UTRs between black grouse BLB1 and BLB2 and could prove orthology of the BLB1 gene between the black grouse, chicken and pheasant, and orthology of the BLB2 gene between the three species. The origin of the pheasant like birds Phasianoidea is estimated to approximately 40 million years ago [84] and chicken and black grouse/turkey diverged approximately 30 million years ago [85]. It seems that the duplication of BLB1 and BLB2 is a case of pre-speciation duplication that has arisen in the ancestral species before the split into chicken, pheasant and black grouse. The observed differences in the length of the 3′UTRs between the two loci may also reflect a difference in function [86].

In all other phylogenetic reconstructions based on the exon 2, the third codon positions of exon 2 and the exon 3, sequences clustered species-specific and not locus-specific. This is indication for frequent interlocus genetic exchange homogenizing sequences between the paralogous loci, as suggested by Wittzell et al. [47]. The fact that the whole exon 2 sequences and third codon positions of the exon 2 networks present similar phylogenetic relations contains more indication for an early duplication in the evolutionary history of BLB1 and BLB2. Under convergent selection on the two loci, the two phylogenies would reveal deviating patterns [34]. A phylogenetic tree based on the codon sequences should reflect functional similarities, whereas the third codon positions are expected to mirror the neutral gene history. Under the early duplication hypothesis, in contrast, a phylogeny based on the third codon positions will match the phylogenetic relation based on the whole codon sequences, as in our case. Note that in this analysis all third codon positions are considered and not only the synonymous positions at the PBR, which are possibly under a selective sweep of the positions under selection, as discussed earlier.

In summary, we infer that the BLB gene duplicated before the species divergence into chicken, black grouse and pheasant and thus is a case of pre speciation duplication. Further, we conclude that BLB1 and BLB2 in black grouse are subjected to homogenizing concerted evolution due to inter-genetic exchange between loci after species divergence. Different selection patterns indicated for BLB1 and BLB2 may be a sign of different immunogenetic functions. Both BLB1 and BLB2 have been under balancing selection during their evolutionary history. Nevertheless, it is likewise possible that, at present, balancing selection may be operating directly only on BLB2. Because of the tight linkage between the loci, BLB1 would be hitchhiking with BLB2.


We have presented a powerful single locus genotyping method for amplifying MHC class IIB loci in black grouse. This method will allow exploring correlations between MHC heterozygosity and reproductive success, mate choice and disease resistance in this and related species.

It is evident that both BLB1 and BLB2 are expressed and under balancing selection in black grouse. Even though we found some differences between the loci in selection patterns, the similarities between loci were apparent. In fact, our data show inter-locus genetic exchange between BLB1 and BLB2. The loci are in linkage disequilibrium, which is in line with what have been stated regarding tightly coevolving genes within MHC under the minimal essential MHC hypothesis [31, 36]. Our results support the conclusion that MHC form and function in birds derived from studies on the domesticated chicken are not artefacts of the domestication process. However, the data from black grouse do suggest more recombination than previously observed in chicken (but see [39]). These data are important for understanding how the MHC of birds and other non-mammalian vertebrates have evolved.

The duplicated genes BLB1 and BLB2 are a case of early duplication and have co-evolved in a concerted way by interlocus gene exchange not only in chicken but also in black grouse. This supports the so far non-tested hypothesis of concerted evolution of BLB1 and BLB2 in galliformes [47].


  1. 1.

    Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290 (5494): 1151-1155. 10.1126/science.290.5494.1151.

    Article  CAS  PubMed  Google Scholar 

  2. 2.

    Klein J: Natural history of the major histocompatibility complex. 1986, New York: John Wiley and sons

    Google Scholar 

  3. 3.

    Spurgin LG, Richardson DS: How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings. Proc R Soc B-Biolo Scie. 2010, 277 (1684): 979-988. 10.1098/rspb.2009.2084.

    Article  CAS  Google Scholar 

  4. 4.

    Piertney S, Oliver M: The evolutionary ecology of the major histocompatibility complex. Heredity. 2006, 96 (1): 7-21.

    CAS  PubMed  Google Scholar 

  5. 5.

    Robinson J, Mistry K, McWilliam H, Lopez R, Parham P, Marsh SGE: The IMGT/HLA database. Nucleic Acids Res. 2011, 39 (suppl 1): D1171-D1176.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  6. 6.

    Apanius V, Penn D, Slev PR, Ruff LR, Potts WK: The nature of selection on the major histocompatibility complex. Crit Rev Immunol. 1997, 17 (2): 179-224. 10.1615/CritRevImmunol.v17.i2.40.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Bernatchez L, Landry C: MHC studies in nonmodel vertebrates: what have we learned about natural selection in 15 years?. J Evolution Biol. 2003, 16 (3): 363-377. 10.1046/j.1420-9101.2003.00531.x.

    Article  CAS  Google Scholar 

  8. 8.

    Hughes AL, Yeager M: Natural selection at major histocompatibility complex loci of vertebrates. Annu Rev Genet. 1998, 32: 415-435. 10.1146/annurev.genet.32.1.415.

    Article  CAS  PubMed  Google Scholar 

  9. 9.

    Martinsohn JT, Sousa AB, Guethlein LA, Howard JC: The gene conversion hypothesis of MHC evolution: a review. Immunogenetics. 1999, 50 (3): 168-200. 10.1007/s002510050593.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Klein J: Origin of major histocompatibility complex polymorphism: the transspecies hypothesis. Hum Immunol. 1987, 19: 155-162. 10.1016/0198-8859(87)90066-8.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    van Oosterhout C: A new theory of MHC evolution: beyond selection on the immune genes. Proc R Soc B-Biol Scie. 2009, 276 (1657): 657-665. 10.1098/rspb.2008.1299.

    Article  CAS  Google Scholar 

  12. 12.

    Edwards SV, Gasper J, March M: Genomics and polymorphism of Agph-DAB1, an MHC class II B gene in red-winged blackbirds (Agelaius phoeniceus). Mol Biol Evol. 1998, 15 (3): 236-250. 10.1093/oxfordjournals.molbev.a025921.

    Article  CAS  PubMed  Google Scholar 

  13. 13.

    Ekblom R, Grahn M, Höglund J: Patterns of polymorphism in the MHC class II of a non-passerine bird, the great snipe (Gallinago media). Immunogenetics. 2003, 54 (10): 734-741.

    CAS  PubMed  Google Scholar 

  14. 14.

    Miller HC, Allendorf F, Daugherty CH: Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.). Mol Ecol. 2010, 19 (18): 3894-3908. 10.1111/j.1365-294X.2010.04771.x.

    Article  PubMed  Google Scholar 

  15. 15.

    Strand T, Westerdahl H, Höglund J, Alatalo RV, 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 (9): 725-734. 10.1007/s00251-007-0234-6.

    Article  CAS  PubMed  Google Scholar 

  16. 16.

    Westerdahl H, Wittzell H, von Schantz T: Polymorphism and transcription of Mhc class I genes in a passerine bird, the great reed warbler. Immunogenetics. 1999, 49 (3): 158-170. 10.1007/s002510050477.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Miller HC, Lambert DM: Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae). Mol Ecol. 2004, 13 (12): 3709-3721. 10.1111/j.1365-294X.2004.02368.x.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Hughes CR, Miles S, Walbroehl JM: Support for the minimal essential MHC hypothesis: a parrot with a single, highly polymorphic MHC class IIB gene. Immunogenetics. 2008, 60 (5): 219-231. 10.1007/s00251-008-0287-1.

    Article  CAS  PubMed  Google Scholar 

  19. 19.

    Tsuda TT, Tsuda M, Naruse T, Kawata H, Ando A, Shiina T, Fukuda M, Kurita M, LeMaho I, Kulski JK: Phylogenetic analysis of penguin (Spheniscidae) species based on sequence variation in MHC class II genes. Immunogenetics. 2001, 53 (8): 712-716. 10.1007/s002510100369.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Alcaide M, Edwards S, Negro J: Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey. J Mol Evol. 2007, 65: 541-554. 10.1007/s00239-007-9033-9.

    Article  CAS  PubMed  Google Scholar 

  21. 21.

    O’Neill A, Livant E, Ewald S: The chicken < i > BF1 (classical MHC class I) gene shows evidence of selection for diversity in expression and in promoter and signal peptide regions. Immunogenetics. 2009, 61 (4): 289-302. 10.1007/s00251-008-0354-7.

    Article  PubMed  Google Scholar 

  22. 22.

    Worley K, Gillingham M, Jensen P, Kennedy LJ, Pizzari T, Kaufman J, Richardson DS: Single locus typing of MHC class I and class IIB loci in a population of red jungle fowl. Immunogenetics. 2008, 60 (5): 233-247. 10.1007/s00251-008-0288-0.

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    Reed K, Bauer M, Monson M, Benoit B, Chaves L, O’Hare T, Delany M: Defining the Turkey MHC: identification of expressed class I- and class IIB-like genes independent of the MHC-B. Immunogenetics. 2011, 63: 753-771. 10.1007/s00251-011-0549-1.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Burri R, Niculita-Hirzel H, Roulin A, Fumagalli L: Isolation and characterization of major histocompatibility complex (MHC) class II B genes in the Barn owl (Aves: Tyto alba). Immunogenetics. 2008, 60 (9): 543-550. 10.1007/s00251-008-0308-0.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Kelley J, Walter L, Trowsdale J: Comparative genomics of major histocompatibility complexes. Immunogenetics. 2005, 56 (10): 683-695. 10.1007/s00251-004-0717-7.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Hess CM, Edwards SV: The Evolution of the Major Histocompatibility Complex in Birds. Bioscience. 2002, 52 (5): 423-431. 10.1641/0006-3568(2002)052[0423:TEOTMH]2.0.CO;2.

    Article  Google Scholar 

  27. 27.

    Kulski J, Shiina T, Anzai T, Kohara S, Inoko H: Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man. Immunol Rev. 2002, 190 (1): 95-122. 10.1034/j.1600-065X.2002.19008.x.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Colloquium on Genetics and the Origin of Species: Jan 30-Feb 01 1997; Irvine, Ca: Natl Acad Sciences. 1997, 94: 7799-7806.

    CAS  Google Scholar 

  29. 29.

    Nei M, Rooney A: Concerted and birth-and-death evolution of multigene families. Ann Rev Gen. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.

    Article  CAS  Google Scholar 

  30. 30.

    Kaufman J: The avian MHC. Avian Immunology. Edited by: Davison F, Kaspers B, Schat KA. 2008, Oxford: Academic

    Google Scholar 

  31. 31.

    Walker BA, Hunt LG, Sowa AK, Skjødt K, Göbel TW, Lehner PJ, Kaufman J: The dominantly expressed class I molecule of the chicken MHC is explained by coevolution with the polymorphic peptide transporter (TAP) genes. Proc Natl Acad Sci. 2011, 108 (20): 8396-8401. 10.1073/pnas.1019496108.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  32. 32.

    Burri R, Hirzel H, Salamin N, Roulin A, Fumagalli L: Evolutionary patterns of MHC class IIB in owls and their implications for the understanding of avian MHC evolution. Mol Biol Evol. 2008, 25 (6): 1180-1191. 10.1093/molbev/msn065.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Fawcett J, Innan H: Neutral and Non-Neutral Evolution of Duplicated Genes with Gene Conversion. Genes. 2011, 2 (1): 191-209.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  34. 34.

    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 (10): 2360-2374. 10.1093/molbev/msq120.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Briles WE, McGibbon WH, Irwin MR: On multiple alleles effecting cellular antigens in the chicken. Genetics. 1950, 35 (6): 633-652.

    PubMed Central  CAS  PubMed  Google Scholar 

  36. 36.

    Kaufman J, Milne S, Gobel TWF, Walker BA, Jacob JP, Auffray C, Zoorob R, Beck S: The chicken B locus is a minimal essential major histocompatibility complex. Nature. 1999, 401 (6756): 923-925. 10.1038/44856.

    Article  CAS  PubMed  Google Scholar 

  37. 37.

    The MHC sequencing consortium: Complete sequence and gene map of a human major histocompatibility complex. Nature. 1999, 401 (6756): 921-10.1038/44853.

    Article  Google Scholar 

  38. 38.

    Balakrishnan C, Ekblom R, Volker M, Westerdahl H, Godinez R, Kotkiewicz H, Burt D, Graves T, Griffin D, Warren W: Gene duplication and fragmentation in the zebra finch major histocompatibility complex. BMC Biol. 2010, 8 (1): 29-10.1186/1741-7007-8-29.

    PubMed Central  Article  PubMed  Google Scholar 

  39. 39.

    Hosomichi K, Miller MM, Goto RM, Wang YJ, Suzuki S, Kulski JK, Nishibori M, Inoko H, Hanzawa K, Shiina T: Contribution of mutation, recombination, and gene conversion to chicken MHC-B haplotype diversity. J Immunol. 2008, 181 (5): 3393-3399.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  40. 40.

    Kaufman J: Co-evolving genes in MHC haplotypes: the “rule” for nonmammalian vertebrates?. Immunogenetics. 1999, 50 (3–4): 228-236.

    Article  CAS  PubMed  Google Scholar 

  41. 41.

    Skjodt K, Koch C, Crone M, Simonsen M: Analysis of chickens for recombination within the MHC (B-complex). Tissue Antigens. 1985, 25 (5): 278-282.

    Article  CAS  PubMed  Google Scholar 

  42. 42.

    Chaves LD, Faile GM, Krueth SB, Hendrickson JA, Reed KM: Haplotype variation, recombination, and gene conversion within the turkey MHC-B locus. Immunogenetics. 2010, 62 (7): 465-477. 10.1007/s00251-010-0451-2.

    Article  CAS  PubMed  Google Scholar 

  43. 43.

    Shiina T, Shimizu S, Hosomichi K, Kohara S, Watanabe S, Hanzawa K, Beck S, Kulski JK, Inoko H: Comparative genomic analysis of two avian (quail and chicken) MHC regions. J Immunol. 2004, 172 (11): 6751-6763.

    Article  CAS  PubMed  Google Scholar 

  44. 44.

    Jacob JP, Milne S, Beck S, Kaufman J: The major and a minor class II beta-chain (B-LB) gene flank the Tapasin gene in the B-F/B-L region of the chicken major histocompatibility complex. Immunogenetics. 2000, 51 (2): 138-147. 10.1007/s002510050022.

    Article  CAS  PubMed  Google Scholar 

  45. 45.

    Pharr GT, Dodgson JB, Hunt HD, Bacon LD: Class II MHC cDNAs in 15I(5) B-congenic chickens. Immunogenetics. 1998, 47 (5): 350-354. 10.1007/s002510050369.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    Shaw I, Powell TJ, Marston DA, Baker K, van Hateren A, Riegert P, Wiles MV, Milne S, Beck S, Kaufman J: 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 (9): 5744-5752.

    Article  CAS  PubMed  Google Scholar 

  47. 47.

    Wittzell H, Madsen T, Westerdahl H, Shine R, von Schantz T: MHC variation in birds and reptiles. Genetica. 1999, 104 (3): 301-309.

    Article  CAS  Google Scholar 

  48. 48.

    Wang B, Ekblom R, Strand T, Portela-Bens S, Hoglund J: Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC. BMC Genomics. 2012, 13 (1): 553-10.1186/1471-2164-13-553.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  49. 49.

    Li XX, Han LX, Han JL: No Specific Primer Can Independently Amplify the Complete Exon 2 of Chicken BLB1 or BLB2 Genes. Int J Poultry Scie. 2010, 9 (2): 192-197.

    Article  CAS  Google Scholar 

  50. 50.

    Strand T, Höglund J: Genotyping of black grouse MHC class II B using reference Strand-Mediated Conformational Analysis (RSCA). BMC Res Notes. 2011, 4 (1): 183-10.1186/1756-0500-4-183.

    PubMed Central  Article  PubMed  Google Scholar 

  51. 51.

    Paxton RJ, Thoren PA, Tengo J, Estoup A, Pamilo P: Mating structure and nestmate relatedness in a communal bee, Andrena jacobi (Hymenoptera, Andrenidae), using microsatellites. Mol Ecol. 1996, 5 (4): 511-519. 10.1111/j.1365-294X.1996.tb00343.x.

    Article  CAS  PubMed  Google Scholar 

  52. 52.

    Rozen S, Skaletsky H: Primer3 on the WWW for General Users and for Biologist Programmers Bioinformatics Methods and Protocols. Edited by: Misener S, Krawetz SA. 2000, Humana Press, 132: 365-386.

    Google Scholar 

  53. 53.

    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 (1–2): 117-123.

    Article  CAS  PubMed  Google Scholar 

  54. 54.

    Westerdahl H, Wittzell H, Schantz T, Bensch S: MHC class I typing in a songbird with numerous loci and high polymorphism using motif-specific PCR and DGGE. Heredity. 2004, 92 (6): 534-542. 10.1038/sj.hdy.6800450.

    Article  CAS  PubMed  Google Scholar 

  55. 55.

    Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.

    Article  CAS  PubMed  Google Scholar 

  56. 56.

    Alcaide M, Rodriguez A, Negro J: Sampling strategies for accurate computational inferences of gametic phase across highly polymorphic major histocompatibility complex loci. BMC Res Notes. 2011, 4 (1): 151-10.1186/1756-0500-4-151.

    PubMed Central  Article  PubMed  Google Scholar 

  57. 57.

    Bos D, Turner S, Dewoody J: Haplotype inference from diploid sequence data: evaluating performance using non-neutral MHC sequences. Hereditas. 2007, 144: 228-234. 10.1111/j.2007.0018-0661.01994.x.

    Article  PubMed  Google Scholar 

  58. 58.

    Mansai SP, Innan H: The power of the methods for detecting interlocus gene conversion. Genetics. 2010, 184 (2): 517-527. 10.1534/genetics.109.111161.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  59. 59.

    Martin DP, Lemey P, Posada D: Analysing recombination in nucleotide sequences. Mol Ecol Resour. 2011, 11 (6): 943-955. 10.1111/j.1755-0998.2011.03026.x.

    Article  PubMed  Google Scholar 

  60. 60.

    Hudson RR, Kaplan NL: Statistical properties of the number of recombination events in the history of a sample of DNA-sequences. Genetics. 1985, 111 (1): 147-164.

    PubMed Central  CAS  PubMed  Google Scholar 

  61. 61.

    Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P: RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010, 26 (19): 2462-2463. 10.1093/bioinformatics/btq467.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  62. 62.

    Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000, 16 (6): 562-563. 10.1093/bioinformatics/16.6.562.

    Article  CAS  PubMed  Google Scholar 

  63. 63.

    Maynard Smith J: Analyzing the mosaic structure of genes. J Mol Evol. 1992, 34 (2): 126-129.

    Google Scholar 

  64. 64.

    Posada D, Crandall KA: Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. Proc Natl Acad Sci USA. 2001, 98 (24): 13757-13762. 10.1073/pnas.241370698.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  65. 65.

    Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology. 1999, 265 (2): 218-225. 10.1006/viro.1999.0056.

    Article  CAS  PubMed  Google Scholar 

  66. 66.

    Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SDW: Automated phylogenetic detection of recombination using a genetic algorithm. Mol Biol Evol. 2006, 23 (10): 1891-1901. 10.1093/molbev/msl051.

    Article  PubMed  Google Scholar 

  67. 67.

    Delport W, Poon AFY, Frost SDW, Pond SLK: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 2010, 26 (19): 2455-2457. 10.1093/bioinformatics/btq429.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  68. 68.

    McVean GAT, Myers SR, Hunt S, Deloukas P, Bentley DR, Donnelly P: The fine-scale structure of recombination rate variation in the human genome. Science. 2004, 304 (5670): 581-584. 10.1126/science.1092500.

    Article  CAS  PubMed  Google Scholar 

  69. 69.

    Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10 (3): 564-567. 10.1111/j.1755-0998.2010.02847.x.

    Article  PubMed  Google Scholar 

  70. 70.

    Hurst LD: Genetics and the understanding of selection. Nat Rev Genet. 2009, 10 (2): 83-93.

    Article  CAS  PubMed  Google Scholar 

  71. 71.

    Garrigan D, Hedrick PW, Mitton J: Perspective: Detecting adaptive molecular polymorphism: Lessons from the MHC. Evolution. 2003, 57 (8): 1707-1722.

    Article  CAS  PubMed  Google Scholar 

  72. 72.

    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-2739. 10.1093/molbev/msr121.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  73. 73.

    Tong J, Bramson J, Kanduc D, Chow S, Sinha A, Ranganathan S: Modeling the bound conformation of Pemphigus Vulgaris-associated peptides to MHC Class II DR and DQ Alleles. Immunome Res. 2006, 2 (1): 1-10.1186/1745-7580-2-1.

    PubMed Central  Article  PubMed  Google Scholar 

  74. 74.

    Yang ZH: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.

    Article  CAS  PubMed  Google Scholar 

  75. 75.

    Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.

    Article  CAS  PubMed  Google Scholar 

  76. 76.

    Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.

    Article  CAS  PubMed  Google Scholar 

  77. 77.

    Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23 (2): 254-267.

    Article  CAS  PubMed  Google Scholar 

  78. 78.

    Spurgin LG, van Oosterhout C, Illera JC, Bridgett S, Gharbi K, Emerson BC, Richardson DS: Gene conversion rapidly generates major histocompatibility complex diversity in recently founded bird populations. Mol Ecol. 2011, 20 (24): 5213-5225. 10.1111/j.1365-294X.2011.05367.x.

    Article  CAS  PubMed  Google Scholar 

  79. 79.

    Agudo R, Alcaide M, Rico C, Lemus JA, Blanco G, Hiraldo F, DonÁZar JA: Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection. Mol Ecol. 2011, 20 (11): 2329-2340. 10.1111/j.1365-294X.2011.05107.x.

    Article  PubMed  Google Scholar 

  80. 80.

    Bettencourt B, Feder M: Rapid Concerted Evolution via Gene Conversion at the Drosophila hso70 Genes. J Mol Evol. 2002, 54: 569-586. 10.1007/s00239-001-0044-7.

    Article  CAS  PubMed  Google Scholar 

  81. 81.

    Sharon D, Gilad Y, Glusman G, Khen M, Lancet D, Kalush F: Identification and characterization of coding single-nucleotide polymorphisms within a human olfactory receptor gene cluster. Gene. 2000, 260 (1–2): 87-94.

    Article  CAS  PubMed  Google Scholar 

  82. 82.

    Anisimova M, Nielsen R, Yang Z: Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics. 2003, 164 (3): 1229-1236.

    PubMed Central  CAS  PubMed  Google Scholar 

  83. 83.

    Ohta T: Gene conversion vs point mutation in generating variability at the antigen recognition site of Major Histocompatibility Complex loci. J Mol Evol. 1995, 41 (2): 115-119.

    Article  CAS  PubMed  Google Scholar 

  84. 84.

    Helm-Bychowski KM, Wilson AC: Rates of nuclear DNA evolution in pheasant-like birds: evidence from restriction maps. Proc Natl Acad Sci. 1986, 83 (3): 688-692. 10.1073/pnas.83.3.688.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  85. 85.

    Pereira SL, Baker AJ: A molecular timescale for galliform birds accounting for uncertainty in time estimates and heterogeneity of rates of DNA substitutions across lineages and sites. Mol Phylogenet Evol. 2006, 38 (2): 499-509. 10.1016/j.ympev.2005.07.007.

    Article  CAS  PubMed  Google Scholar 

  86. 86.

    Stark A, Brennecke J, Bushati N, Russell RB, Cohen SM: Animal MicroRNAs Confer Robustness to Gene Expression and Have a Significant Impact on 3′UTR Evolution. Cell. 2005, 123 (6): 1133-1146. 10.1016/j.cell.2005.11.023.

    Article  CAS  PubMed  Google Scholar 

Download references


Thanks to Heli Siitari, Rauno Alatalo, and Sten Svartaas for providing black grouse samples. We are grateful to Gunilla Engström for assistance in the lab, to Robert Ekblom and three anonymous reviewers for valuable comments on the manuscript and Reto Burri for helpful discussions. Thanks also to Eleanor Jones and Radoslav Kozma for language improvements. This project was funded through the Swedish Research Council (VR) to JH.

Author information



Corresponding author

Correspondence to Jacob Höglund.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TS and JH conceived and designed the study, TS and BW carried out the molecular genetic studies. TS, YML and BW performed the statistical analyses and interpretation of the data. The manuscript was written by TS and YML, with input from JH and BW. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: MHC primers and corresponding amplification lengths. The fragment column highlights the locus and region amplified with each primer pair. The PCR product length is given without primers. (DOCX 16 KB)

Additional file 2: Best evolutionary models estimated by Aikaike’s information criterion in jModelTest. The models were used to construct phylogenetic trees. (DOCX 14 KB)

Additional file 3: Alignment of black grouse MHC class II B exon 2 nucleotide sequences. Locus-designated sequences derived in the present study are indicated with 125 bp and 251 bp. Sequences indicated as BLB* are cloned in previous studies and not designated to locus, for example 202 bp cDNA sequences (from individual D870, see Table 1). Sequences likely to be from the same allele are grouped by grey shading. (DOCX 40 KB)

Supplement for the test for positive selection.

Additional file 4: Likelihood values and parameter estimates for the different models calculated with CODEML implemented in PAML 4.6 [74]. M = model, lnL = Log-likelihood value. In the neutral model M0, ω is equivalent to averaged dN/dS. Dark grey shade highlights the significantly best models. (DOCX 18 KB)

Neighbour Joining tree for exon 3.

Additional file 5: BLB1 and BLB2 sequences derived from black grouse (fosmid individual JHGO 213 [48] [GenBank JQ028669] and chicken [GenBank AB268588]. Goose [GenBank EU999169] was used as an outgroup. (DOCX 48 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Strand, T., Wang, B., Meyer-Lucht, Y. et al. Evolutionary history of black grouse major histocompatibility complex class IIB genes revealed through single locus sequence-based genotyping. BMC Genet 14, 29 (2013).

Download citation


  • Balancing selection
  • Concerted evolution
  • Early duplication
  • BLB1
  • BLB2
  • Galliformes
  • 3′UTR
  • Tetrao tetrix