Dynamics of Mhc evolution in birds and crocodilians: amplification of class II genes with degenerate primers

Molecular Tools for Lynx spp. qPCR Identification and STR-Based Individual Identification of Eurasian Lynx (Lynx lynx) in Forensic Casework

The Eurasian lynx (Lynx lynx) is listed in CITES Appendix II and is protected under the Bern Convention and the EU Habitats Directive, yet it remains a frequent target of wildlife crime, highlighting the urgent need for reliable identification methods. This study focuses on determination and DNA quantification of the Lynx spp. using quantitative real-time PCR (qPCR). The Llynx Qplex quantification multiplex system effectively distinguishes Lynx spp. from other Feliformia species by targeting mitochondrial and nuclear markers. Additionally, we present the results of the developmental validation of the Llyn STRplex system for individual identification and databasing using six STR loci. This study followed ISFG recommendations for non-human DNA testing and developmental validation guidelines. Both systems demonstrate high sensitivity (5 pg genomic DNA for Llynx Qplex and 30 pg of mtDNA for Llyn STRplex) and high specificity to Lynx spp., confirmed by testing against 16 related Feliformia species. Robustness was evaluated, showing sensitivity to temperature variation, and both repeatability and reproducibility were successfully tested across replicates and conditions. Given that forensic casework often involves degraded and limited biological material, molecular tools must be both sensitive and specific to ensure accurate results. Developing precise and efficient tools is essential for supporting investigations of wildlife crime involving the Eurasian lynx, as well as efforts aimed at conserving the species.

Expansion of MHC-IIB Has Constrained the Evolution of MHC-IIA in Passerines

The major histocompatibility complex (MHC) is central in adaptive immunity, with the highly polymorphic MHC genes encoding antigen-presenting molecules. Two MHC class II (MHC-II) loci, DA1 and DA2, predate the radiation of extant birds and persist throughout much of the avian phylogeny. Within each locus, the MHC-II molecules are encoded by A-genes (DAA) and B-genes (DAB), which are arranged in A-B dyads. However, in passerines (order Passeriformes), the DA2 locus has been lost, and the ancestral A-B dyad at the DA1 locus has been replaced by a putatively single A-gene (DAA1) and an array of highly polymorphic B-genes (DAB1). In this study, we genotyped the DAA1 gene of 15 passerine species and confirmed that passerines possess just one copy of DAA1. We then compared selection patterns in DAA1 between passerines and nonpasserines and found that exon 2, which encodes the antigen-presenting domain, has been subject to weaker positive selection and stronger negative selection in passerines compared with nonpasserines. Additional comparisons showed that the patterns of selection in the passerine DAA1 gene are unlikely to be related to the loss of the DA2 locus. Instead, our findings suggest that the expansion of DAB1 (MHC-IIB) has imposed an evolutionary constraint on the passerine DAA1 (MHC-IIA) gene. We speculate that this constraint may be the result of each DAA1 chain forming heterodimers with many different DAB1 chains.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

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

Background

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

Results

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

Conclusion

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

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

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

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

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

New data from basal Australian songbird lineages show that complex structure of MHC class II β genes has early evolutionary origins within passerines

Background

The major histocompatibility complex (MHC) plays a crucial role in the adaptive immune system and has been extensively studied across vertebrate taxa. Although the function of MHC genes appears to be conserved across taxa, there is great variation in the number and organisation of these genes. Among avian species, for instance, there are notable differences in MHC structure between passerine and non-passerine lineages: passerines typically have a high number of highly polymorphic MHC paralogs whereas non-passerines have fewer loci and lower levels of polymorphism. Although the occurrence of highly polymorphic MHC paralogs in passerines is well documented, their evolutionary origins are relatively unexplored. The majority of studies have focussed on the more derived passerine lineages and there is very little empirical information on the diversity of the MHC in basal passerine lineages. We undertook a study of MHC diversity and evolutionary relationships across seven species from four families (Climacteridae, Maluridae, Pardalotidae, Meliphagidae) that comprise a prominent component of the basal passerine lineages. We aimed to determine if highly polymorphic MHC paralogs have an early evolutionary origin within passerines or are a more derived feature of the infraorder Passerida.

Results

We identified 177 alleles of the MHC class II β exon 2 in seven basal passerine species, with variation in numbers of alleles across individuals and species. Overall, we found evidence of multiple gene loci, pseudoalleles, trans-species polymorphism and high allelic diversity in these basal lineages. Phylogenetic reconstruction of avian lineages based on MHC class II β exon 2 sequences strongly supported the monophyletic grouping of basal and derived passerine species.

Conclusions

Our study provides evidence of a large number of highly polymorphic MHC paralogs in seven basal passerine species, with strong similarities to the MHC described in more derived passerine lineages rather than the simpler MHC in non-passerine lineages. These findings indicate an early evolutionary origin of highly polymorphic MHC paralogs in passerines and shed light on the evolutionary forces shaping the avian MHC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0681-5) contains supplementary material, which is available to authorized users.

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

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

Patterns of evolution of MHC class II genes of crows (Corvus) suggest trans-species polymorphism

A distinguishing characteristic of genes that code for the major histocompatibility complex (MHC) is that alleles often share more similarity between, rather than within species. There are two likely mechanisms that can explain this pattern: convergent evolution and trans-species polymorphism (TSP), in which ancient allelic lineages are maintained by balancing selection and retained by descendant species. Distinguishing between these two mechanisms has major implications in how we view adaptation of immune genes. In this study we analyzed exon 2 of the MHC class IIB in three passerine bird species in the genus Corvus: jungle crows (Corvus macrorhynchos japonensis) American crows (C. brachyrhynchos) and carrion crows (C. corone orientalis). Carrion crows and American crows are recently diverged, but allopatric, sister species, whereas carrion crows and jungle crows are more distantly related but sympatric species, and possibly share pathogens linked to MHC IIB polymorphisms. These patterns of evolutionary divergence and current geographic ranges enabled us to test for trans-species polymorphism and convergent evolution of the MHC IIB in crows. Phylogenetic reconstructions of MHC IIB sequences revealed several well supported interspecific clusters containing all three species, and there was no biased clustering of variants among the sympatric carrion crows and jungle crows. The topologies of phylogenetic trees constructed from putatively selected sites were remarkably different than those constructed from putatively neutral sites. In addition, trees constructed using non-synonymous substitutions from a continuous fragment of exon 2 had more, and generally more inclusive, supported interspecific MHC IIB variant clusters than those constructed from the same fragment using synonymous substitutions. These phylogenetic patterns suggest that recombination, especially gene conversion, has partially erased the signal of allelic ancestry in these species. While clustering of positively selected amino acids by supertyping revealed a single supertype shared by only jungle and carrion crows, a pattern consistent with convergence, the overall phylogenetic patterns we observed suggest that TSP, rather than convergence, explains the interspecific allelic similarity of MHC IIB genes in these species of crows.

Preen secretions encode information on MHC similarity in certain sex-dyads in a monogamous seabird

Animals are known to select mates to maximize the genetic diversity of their offspring in order to achieve immunity against a broader range of pathogens. Although several bird species preferentially mate with partners that are dissimilar at the major histocompatibility complex (MHC), it remains unknown whether they can use olfactory cues to assess MHC similarity with potential partners. Here we combined gas chromatography data with genetic similarity indices based on MHC to test whether similarity in preen secretion chemicals correlated with MHC relatedness in the black-legged kittiwake (Rissa tridactyla), a species that preferentially mates with genetically dissimilar partners. We found that similarity in preen secretion chemicals was positively correlated with MHC relatedness in male-male and male-female dyads. This study provides the first evidence that preen secretion chemicals can encode information on MHC relatedness and suggests that odor-based mechanisms of MHC-related mate choice may occur in birds.

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

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

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

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

Evolution of MHC class I in the order Crocodylia

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

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

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

MHC class II DQB diversity in the Japanese black bear, Ursus thibetanus japonicus

Background

The major histocompatibility complex (MHC) genes are one of the most important genetic systems in the vertebrate immune response. The diversity of MHC genes may directly influence the survival of individuals against infectious disease. However, there has been no investigation of MHC diversity in the Asiatic black bear (Ursus thibetanus). Here, we analyzed 270-bp nucleotide sequences of the entire exon 2 region of the MHC DQB gene by using 188 samples from the Japanese black bear (Ursus thibetanus japonicus) from 12 local populations.

Results

Among 185 of 188 samples, we identified 44 MHC variants that encoded 31 different amino acid sequences (allotypes) and one putative pseudogene. The phylogenetic analysis suggests that MHC variants detected from the Japanese black bear are derived from the DQB locus. One of the 31 DQB allotypes, Urth-DQB*01, was found to be common to all local populations. Moreover, this allotype was shared between the black bear on the Asian continent and the Japanese black bear, suggesting that Urth-DQB*01 might have been maintained in the ancestral black bear population for at least 300,000 years. Our findings, from calculating the ratio of non-synonymous to synonymous substitutions, indicate that balancing selection has maintained genetic variation of peptide-binding residues at the DQB locus of the Japanese black bear. From examination of genotype frequencies among local populations, we observed a considerably lower level of observed heterozygosity than expected.

Conclusions

The low level of observed heterozygosity suggests that genetic drift reduced DQB diversity in the Japanese black bear due to a bottleneck event at the population or species level. The decline of DQB diversity might have been accelerated by the loss of rare variants that have been maintained by negative frequency-dependent selection. Nevertheless, DQB diversity of the black bear appears to be relatively high compared with some other endangered mammalian species. This result suggests that the Japanese black bears may also retain more potential resistance against pathogens than other endangered mammalian species. To prevent further decline of potential resistance against pathogens, a conservation policy for the Japanese black bear should be designed to maintain MHC rare variants in each local population.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird

The extent to which pathogens maintain the extraordinary polymorphism at vertebrate Major Histocompatibility Complex (MHC) genes via balancing selection has intrigued evolutionary biologists for over half a century, but direct tests remain challenging. Here we examine whether a well-characterized epidemic of Mycoplasmal conjunctivitis resulted in balancing selection on class II MHC in a wild songbird host, the house finch (Carpodacus mexicanus). First, we confirmed the potential for pathogen-mediated balancing selection by experimentally demonstrating that house finches with intermediate to high multi-locus MHC diversity are more resistant to challenge with Mycoplasma gallisepticum. Second, we documented sequence and diversity-based signatures of pathogen-mediated balancing selection at class II MHC in exposed host populations that were absent in unexposed, control populations across an equivalent time period. Multi-locus MHC diversity significantly increased in exposed host populations following the epidemic despite initial compromised diversity levels from a recent introduction bottleneck in the exposed host range. We did not observe equivalent changes in allelic diversity or heterozygosity across eight neutral microsatellite loci, suggesting that the observations reflect selection rather than neutral demographic processes. Our results indicate that a virulent pathogen can exert sufficient balancing selection on class II MHC to rescue compromised levels of genetic variation for host resistance in a recently bottlenecked population. These results provide evidence for Haldane's long-standing hypothesis that pathogens directly contribute to the maintenance of the tremendous levels of genetic variation detected in natural populations of vertebrates.

Molecular cloning, polymorphism and tissue distribution of the MHC class IIB gene in the Chinese goose (Anser cygnoides)

1. The goose major histocompatibility complex (MHC) class IIB cDNA (Ancy-MHCII) was cloned by homology cloning and rapid amplification of cDNA ends by polymerase chain reaction (RACE-PCR), and the genomic structure and tissue expression were investigated. 2. Three different 5'-RACE sequences (Ancy-MHC II5'-1, Ancy-MHC II5'-2, Ancy-MHC II5'-3), one 3'-RACE sequence (Ancy-MHC II-3') and two different full length Ancy-MHC IIB cDNA sequences (Ancy-CD01, Ancy-CD02), which came from different alleles at one locus or different loci, were determined. 3. The genomic organisation is composed of 6 exons and 5 introns, with a longer intron region than that of the chicken. The alleles encode 259 and 260 amino acids in the mature protein. 4. The number of non-synonymous substitutions (dN) in the peptide-binding region of exon 2 from 8 alleles was higher than that of the synonymous substitutions (dS). 5. Tissue-specific expression of Ancy-MHC II mRNA was detected in an adult goose using RT-PCR. These results showed that Ancy-MHC II mRNA was expressed in the lung, spleen, liver, intestine, heart, kidney, pancreas, brain, skin and muscle. This is consistent with the expression of MHC class IIB in various tissues from the chicken. 6. Sequences from goose, snipe and duck clustered together when compared with known MHC class IIB sequences from the other species, significantly differing from mammals and aquatic species, indicating a pattern consistent with accepted evolutionary pathways.

Isolation and characterization of major histocompatibility class IIβ genes in an endangered North American cyprinid fish, the Rio Grande silvery minnow (Hybognathus amarus)

The major histocompatibility complex (MHC) is a critical component of the adaptive immune response in vertebrates. Due to the role that MHC plays in immunity, absence of variation within these genes may cause species to be vulnerable to emerging diseases. The freshwater fish family Cyprinidae comprises the most diverse and species-rich group of freshwater fish in the world, but some are imperiled. Despite considerable species richness and the long evolutionary history of the family, there are very few reports of MHC sequences (apart from a few model species), and no sequences are reported from endemic North American cyprinids (subfamily Leuciscinae). Here we isolate and characterize the MH Class II beta genes from complementary DNA and genomic DNA of the non-model, endangered Rio Grande silvery minnow (Hybognathus amarus), a North American cyprinid. Phylogenetic reconstruction revealed two groups of divergent MH alleles that are paralogous to previously described loci found in deeply divergent cyprinid taxa including common carp, zebrafish, African large barb and bream. Both groups of alleles were under the influence of diversifying selection yet not all individuals had alleles belonging to both allelic groups. We concluded that the general organization and pattern of variation of MH class II genes in Rio Grande silvery minnow is similar to that identified in other cyprinid fishes studied to date, despite distant evolutionary relationships and evidence of a severe genetic bottleneck.

Sequence-based evidence for major histocompatibility complex-disassortative mating in a colonial seabird

The major histocompatibility complex (MHC) is a polymorphic gene family associated with immune defence, and it can play a role in mate choice. Under the genetic compatibility hypothesis, females choose mates that differ genetically from their own MHC genotypes, avoiding inbreeding and/or enhancing the immunocompetence of their offspring. We tested this hypothesis of disassortative mating based on MHC genotypes in a population of great frigatebirds (Fregata minor) by sequencing the second exon of MHC class II B. Extensive haploid cloning yielded two to four alleles per individual, suggesting the amplification of two genes. MHC similarity between mates was not significantly different between pairs that did (n = 4) or did not (n = 42) exhibit extra-pair paternity. Comparing all 46 mated pairs to a distribution based on randomized re-pairings, we observed the following (i): no evidence for mate choice based on maximal or intermediate levels of MHC allele sharing (ii), significantly disassortative mating based on similarity of MHC amino acid sequences, and (iii) no evidence for mate choice based on microsatellite alleles, as measured by either allele sharing or similarity in allele size. This suggests that females choose mates that differ genetically from themselves at MHC loci, but not as an inbreeding-avoidance mechanism.

Characterization and evolution of MHC class II B genes in Ardeid birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.

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

Background

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

Findings

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

Conclusions

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

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

BACKGROUND

Due to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.

RESULTS

The zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.

CONCLUSION

The zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.

Trans-species polymorphism of the Mhc class II DRB-like gene in banded penguins (genus Spheniscus)

The Major Histocompatibility Complex (Mhc) class II DRB locus of vertebrates is highly polymorphic and some alleles may be shared between closely related species as a result of balancing selection in association with resistance to parasites. In this study, we developed a new set of PCR primers to amplify, clone, and sequence overlapping portions of the Mhc class II DRB-like gene from the 5'UTR end to intron 3, including exons 1, 2, and 3 and introns 1 and 2 in four species (20 Humboldt, six African, five Magellanic, and three Galapagos penguins) of penguin from the genus Spheniscus (Sphe). Analysis of gene sequence variation by the neighbor-joining method of 21 Sphe sequences and 20 previously published sequences from four other penguin species revealed overlapping clades within the Sphe species, but species-specific clades for the other penguin species. The overlap of the DRB-like gene sequence variants between the four Sphe species suggests that, despite their allopatric distribution, the Sphe species are closely related and that some shared DRB1 alleles may have undergone a trans-species inheritance because of balancing selection and/or recent rapid speciation. The new primers and PCR assays that we have developed for the identification of the DRB1 DNA and protein sequence variations appear to be useful for the characterization of the molecular evolution of the gene in closely related Penguin species and might be helpful for the assessment of the genetic health and the management of the conservation and captivity of these endangered species.

RSCA genotyping of MHC for high-throughput evolutionary studies in the model organism three-spined stickleback Gasterosteus aculeatus

Background

In all jawed vertebrates, highly polymorphic genes of the major histocompatibility complex (MHC) encode antigen presenting molecules that play a key role in the adaptive immune response. Their polymorphism is composed of multiple copies of recently duplicated genes, each possessing many alleles within populations, as well as high nucleotide divergence between alleles of the same species. Experimental evidence is accumulating that MHC polymorphism is a result of balancing selection by parasites and pathogens. In order to describe MHC diversity and analyse the underlying mechanisms that maintain it, a reliable genotyping technique is required that is suitable for such highly variable genes.

Results

We present a genotyping protocol that uses Reference Strand-mediated Conformation Analysis (RSCA), optimised for recently duplicated MHC class IIB genes that are typical for many fish and bird species, including the three-spined stickleback, Gasterosteus aculeatus. In addition we use a comprehensive plasmid library of MHC class IIB alleles to determine the nucleotide sequence of alleles represented by RSCA allele peaks. Verification of the RSCA typing by cloning and sequencing demonstrates high congruency between both methods and provides new insight into the polymorphism of classical stickleback MHC genes. Analysis of the plasmid library additionally reveals the high resolution and reproducibility of the RSCA technique.

Conclusion

This new RSCA genotyping protocol offers a fast, but sensitive and reliable way to determine the MHC allele repertoire of three-spined sticklebacks. It therefore provides a valuable tool to employ this highly polymorphic and adaptive marker in future high-throughput studies of host-parasite co-evolution and ecological speciation in this emerging model organism.

Structure and evolution of a new avian MHC class II B gene in a sub-Antarctic seabird, the thin-billed prion (Procellariiformes: Pachyptila belcheri)

The major histocompatibility complex encodes molecules that present foreign peptides to T cells of the immune system. The peptide binding region (PBR) of these molecules is among the most polymorphic regions found in vertebrate taxa. Genomic cloning approaches are improving our understanding of the evolution of this multigene family in nonmodel avian groups. By building a cosmid library, a new MHC class II B gene, Pabe-DAB1, was isolated and characterized at the genomic level in a sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri). Pabe-DAB1 exhibits the hallmark structural features of functional MHC class II loci. Direct sequencing of the PBR encoding exon in a panel of prions revealed significantly higher levels of genetic diversity compared to two noncoding neutral loci, with most alleles differing by at least one replacement substitution in the peptide binding codons. We estimated evolutionary dynamics for Pabe-DAB1 using a variety of Bayesian and other approaches. Evidence for balancing selection comes from a spatially variable ratio of nonsynonymous-to-synonymous substitutions (mean d (N)/d (S) = 2.87) in the PBR, with sites predicted to be functionally relevant exhibiting the highest omega values. We estimate the population recombination rate to be approximately 0.3 per site per generation, indicating an important role for recombination in generating polymorphism at this locus. Pabe-DAB1 is among the few avian class II loci characterized at the genomic level and with a known intron-exon structure, a feature that greatly facilitated the amplification and sequencing of a single MHC locus in this species.

MHC haplotype involvement in avian resistance to an ectoparasite

Research on immune function in evolutionary ecology has frequently focused on avian ectoparasites (e.g., mites and lice). However, host immunogenetics involved with bird resistance to ectoparasites has not been determined. The critical role of the major histocompatibility complex (MHC) in adaptive immunity and high genetic variation found within the MHC make this gene complex useful for exploring the immunogenetic basis for bird resistance to ectoparasites. The objective of this study was to determine if the avian MHC influenced resistance to a blood-feeding ectoparasite. Four congenic lines of chickens, differing only at the MHC, were comparatively infested with a cosmopolitan ectoparasite of birds-northern fowl mite (NFM)-which is also a serious pest species of poultry. Mite infestations were monitored over time and mite densities (weekly and maximum) were compared among lines. Chickens with the MHC haplotype B21 were relatively resistant to NFM, compared with birds in the B15 congenic line (P < 0.02). To test for similar effects in an outbred genetic background, a separate experiment was performed with 107 commercial chickens (white leghorn, W-36 strain) infested with NFM. Hens were genotyped using a MHC microsatellite marker (LEI0258) and associations between MHC haplotype and NFM density were tested. The highest peak NFM populations occurred more often on hens with the B15 haplotype versus the B21 haplotype (P = 0.012), which supported the results of the congenic study. These data indicate the avian MHC influences ectoparasite resistance, which is relevant to disease ecology and avian-ectoparasite interaction.

Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.

Molecular characterization of MHC class II in a nonmodel anuran species, the fire-bellied toad Bombina bombina

While the anuran Xenopus comprises one of the best characterized nonmammalian taxa regarding the major histocompatibility complex (MHC), the organization of this gene complex has never been studied in other anurans, and information on amphibian MHC (other than Xenopus) is generally very scarce. Here, we describe the characterization of the first MHC class II B cDNA sequences from a nonmodel anuran species, the European fire-bellied toad (Bombina bombina). We isolated two transcript sequences differing substantially in amino acid composition and length within the beta2 domain. To investigate the variability of the peptide binding region in this species, we sequenced a 158-bp large fragment from wild B. bombina (n = 20) and identified eight distinct alleles. All substitutions but one were nonsynonymous, and many of the highly polymorphic sites corresponded with amino acid positions known to be involved in antigen binding. The level of variation we found in B. bombina was similar compared to that previously found in a comparable sample of a wild urodelan species, Ambystoma tigrinum, and to that found in Xenopus laevis. Based on the cDNA data and the individual's allelic diversity, we conclude that Bombina possesses at least two class II B loci. With our new beta1 primers, we were able to generate sequences in other species of anurans. We provide here a first phylogenetic analysis of this gene in amphibians.

Evolutionary genetics of Carpodacus mexicanus, a recently colonized host of a bacterial pathogen, Mycoplasma gallisepticum

We present molecular data documenting how introduction to the eastern United States and an epizootic involving a bacterial pathogen has affected the genetic diversity of house finches, a cardueline songbird. Population bottlenecks during introduction can cause loss of genetic variation and may negatively affect a population's ability to adapt to novel stressors such as disease. Although a genome-wide survey using Amplified Fragment Length Polymorphism (AFLP) markers suggests little loss of genetic diversity in introduced populations, an epizootic of bacterial Mycoplasma has nonetheless caused dramatic declines in the eastern US population. Sequence analysis of a candidate gene for pathogen resistance in the Major Histocompatibity Complex (MHC) in pre- and post-epizootic population samples reveals allele frequency shifts since introduction of the pathogen, but similar shifts are also observed in control populations not exposed to the bacteria, and in a neutral non-coding locus. Expression studies using a novel subtractive hybridization approach indicate decreased expression of the class II MHC locus upon exposure to Mycoplasma, a pattern also seen in MHC class I loci in mice infected with cytomegalovirus and consistent with manipulation of the finch immune system by Mycoplasma. These results will be further expanded using experimental studies as well as examination of evolution of the pathogen genome itself.

Major histocompatibility complex variation in the endangered crested ibis Nipponia nippon and implications for reintroduction

The major histocompatibility complex (MHC), with its extraordinary levels of genetic variation, is thought to be an essential aspect of the ability of an organism to recognize different parasites and pathogens. It has also been proposed to regulate reproductive processes in many aspects. Here we examine the genetic variation of the second exon of the MHC class II B genes of the crested ibis, an endangered species known to descend from just two breeding pairs rediscovered in 1981. Only five alleles are identified by single-strand conformation polymorphism (SSCP) analysis of 36 samples taken from both wild and captive populations, and a comparatively low level of divergence between MHC alleles is observed. We suggest that representative sampling of individuals with most of the different MHC allele genotypes to constitute a founder population, together with the monitoring of the pathogen status of candidate sites before release, is of great importance for raising the success rate of reintroduction for the crested ibis.

Patterns of variation in MHC class II beta loci of the little greenbul (Andropadus virens) with comments on MHC evolution in birds

We have isolated major histocompatibility complex (MHC) class II beta loci from the little greenbul (Andropadus virens), an African songbird. We utilized preexisting information about conserved regions of the avian MHC to design primers to amplify a pool of sequences representing multiple loci. From this pool, a unique locus spanning 1109 bp that we designate as Anvi-DAB1 was cloned and sequenced. We designed locus-specific primers based on this sequence information and amplified six alleles from seven individuals. Compared to other A. virens MHC sequences obtained from genomic DNA or cDNA, the variability of sequences from Anvi-DAB1 was low and the ratio of nonsynonymous to synonymous substitution was much less than one, suggesting that Anvi-DAB1 may either be a pseudogene or a nonclassical MHC locus. Phylogenetic analysis revealed that the Anvi-DAB1 locus was highly divergent when compared with other passerine or A. virens genomic or transcribed MHC sequences. The use of conserved MHC primers followed by analysis of cloned sequences allows rapid isolation of MHC loci from exotic species and avoids laborious large-scale cloning and sequencing.

Characterization of MHC class II genes from an ancient reptile lineage, Sphenodon (tuatara)

The organization and evolution of major histocompatibility complex (MHC) genes vary considerably among vertebrate lineages. MHC genes have been well characterized in mammals, birds, amphibians and fish, but little is known about their organization in reptiles, despite the fact that reptiles occupy an important phylogenetic position for understanding the evolutionary history of both mammalian and avian MHC genes. Here we describe the characterization of the first MHC class II B cDNA sequences from a non-avian reptile, the tuatara (Sphenodon spp.). Three class II B sequences were isolated from a tuatara cDNA library, and four additional partial sequences were isolated by reverse transcriptase-polymerase chain reaction. Six of these sequences appear to belong to the same gene family, which we have named SppuDAB. The remaining sequence (named SppuDBB) shares only 43.9% amino acid similarity with SppuDAB and thus appears to represent a separate gene family. SppuDBB may be a non-classical locus as it does not contain all the conserved residues expected of a classical MHC class II gene. Southern blot analysis indicates that only a single copy of SppuDBB exists in tuatara, but that multiple loci related to SppuDAB are present. The SppuDAB sequences have the highest amino acid similarity (57.2-62.4%) with class II B sequences from the spectacled caiman, but only 26.4-48.7% similarity with sequences from other vertebrates. The tuatara sequences do not strongly group with other reptile sequences on a phylogenetic tree, reflecting the antiquity of the Sphenodon lineage and the lack of closely related sequences for comparison.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae)

The Chatham Island black robin, Petroica traversi, is a highly inbred, endangered passerine with extremely low levels of variation at hypervariable neutral DNA markers. In this study we investigated variation in major histocompatibility complex (MHC) class II genes in both the black robin and its nonendangered relative, the South Island robin Petroica australis australis. Previous studies have shown that Petroica have at least four expressed class II B MHC genes. In this study, the sequences of introns flanking exon 2 of these loci were characterized to design primers for peptide-binding region (PBR) sequence analysis. Intron sequences were comprised of varying numbers of repeated units, with highly conserved regions immediately flanking exon 2. Polymerase chain reaction primers designed to this region amplified three or four sequences per black robin individual, and eight to 14 sequences per South Island robin individual. MHC genes are fitness-related genes thought to be under balancing selection, so they may be more likely to retain variation in bottlenecked populations. To test this, we compared MHC variation in the black robin with artificially bottlenecked populations of South Island robin, and with their respective source populations, using restriction fragment length polymorphism analyses and DNA sequencing of the PBR. Our results indicate that the black robin is monomorphic at class II B MHC loci, while both source and bottlenecked populations of South Island robin have retained moderate levels of variation. Comparison of MHC variation with minisatellite DNA variation indicates that genetic drift outweighs balancing selection in determining MHC diversity in the bottlenecked populations. However, balancing selection appears to influence MHC diversity over evolutionary timescales, and the effects of gene conversion are evident.

Major histocompatibility complex variation and mate choice in a lekking bird, the great snipe (Gallinago media)

Genes of the major histocompatibility complex (MHC) play a major part in the activation of the vertebrate immune system. In addition, they also appear to function as cues for mate choice. In mammals especially, several kinds of MHC-dependent mate choice have been hypothesized and observed. These include choice of mates that share no or few alleles with the choosing individual, choice of mates with alleles that differ as much as possible from the choosing individual, choice of heterozygous mates, choice of certain genotypes and choice of rare alleles. We investigated these different aspects of mate choice in relation to MHC in a lekking bird species, the great snipe (Gallinago media). We found no evidence for MHC disassortative mating, no preference for males with many MHC alleles and no preference for rare alleles. However, we did find that some allelic lineages were more often found in males with mating success than in males without mating success. Females do not seem to use themselves as references for the MHC-dependent mate choice, rather they seem to prefer males with certain allele types. We speculate that these alleles may be linked to resistance to common parasites.

Comparative genomics of major histocompatibility complexes

The major histocompatibility complex (MHC) is a gene dense region found in all jawed vertebrates examined to date. The MHC contains a high percentage of immune genes, in particular genes involved in antigen presentation, which are generally highly polymorphic. The region plays an important role in disease resistance. The clustering of MHC genes could be advantageous for co-evolution or regulation, and its study in many species is desirable. Even though some linkage of MHC genes is apparent in all gnathostomes, the genomic organization can differ greatly by species, suggesting rapid evolution of MHC genes after divergence from a common ancestor. Previous reviews of comparative MHC organization have been written when relatively fragmentary sequence and mapping data were available on many species. This review compares maps of MHC gene orders in commonly studied species, where extensive sequencing has been performed.

Natural selection of the major histocompatibility complex (Mhc) in Hawaiian honeycreepers (Drepanidinae)

The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.