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.     

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.     

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

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

Characterization,Polymorphism, and 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.     

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

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

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. MHC Class I genes in the Tuatara (Sphenodon spp.): evolution of the MHC in an ancient reptilian order

The major histocompatibility complex (MHC) is an extremely dynamic region of the genome, characterized by high polymorphism and frequent gene duplications and rearrangements. This has resulted in considerable differences in MHC organization and evolution among vertebrate lineages, particularly between birds and mammals. As nonavian reptiles are ancestral to both mammals and birds, they occupy an important phylogenetic position for understanding these differences. However, little is known about reptile MHC genes. To address this, we have characterized MHC class I sequences from the tuatara (Sphenodon spp.), the last survivor of an ancient order of reptiles, Sphenodontia. We isolated two different class I cDNA sequences, which share 93% sequence similarity with each other but are highly divergent from other vertebrate MHC genes. Southern blotting and polymerase chain reaction amplification of class I sequences from seven adult tuatara plus a family group indicate that these sequences represent at least two to three loci. Preliminary analysis of variation among individuals from an island population of tuatara indicates that these loci are highly polymorphic. Maximum likelihood analysis of reptile MHC class I sequences indicates that gene duplication has occurred within reptilian orders. However, the evolutionary relationships among sequences from different reptilian orders cannot be resolved, reflecting the antiquity of the major reptile lineages.     

Retained orthologous relationships of the MHC Class I genes during euteleost evolution

Major histocompatibility complex (MHC) class I molecules play a pivotal role in immune defense system, presenting the antigen peptides to cytotoxic CD8+ T lymphocytes. Most vertebrates possess multiple MHC class I loci, but the analysis of their evolutionary relationships between distantly related species has difficulties because genetic events such as gene duplication, deletion, recombination, and/or conversion have occurred frequently in these genes. Human MHC class I genes have been conserved only within the primates for up to 46-66 My. Here, we performed comprehensive analysis of the MHC class I genes of the medaka fish, Oryzias latipes, and found that they could be classified into four groups of ancient origin. In phylogenetic analysis using these genes and the classical and nonclassical class I genes of other teleost fishes, three extracellular domains of the class I genes showed quite different evolutionary histories. The α1 domains generated four deeply diverged lineages corresponding to four medaka class I groups with high bootstrap values. These lineages were shared with salmonid and/or other acanthopterygian class I genes, unveiling the orthologous relationships between the classical MHC class I genes of medaka and salmonids, which diverged approximately 260 Ma. This suggested that the lineages must have diverged in the early days of the euteleost evolution and have been maintained for a long time in their genome. In contrast, the α3 domains clustered by species or fish groups, regardless of classical or nonclassical gene types, suggesting that this domain was homogenized in each species during prolonged evolution, possibly retaining the potential for CD8 binding even in the nonclassical genes. On the other hand, the α2 domains formed no apparent clusters with the α1 lineages or with species, suggesting that they were diversified partly by interlocus gene conversion, and that the α1 and α2 domains evolved separately. Such evolutionary mode is characteristic to the teleost MHC class I genes and might have contributed to the long-term conservation of the α1 domain.     

MHC class I genes of birds of prey: isolation, polymorphism and diversifying selection

The threat of emerging infectious diseases encourages the investigation of functional loci related to host resilience, such as those belonging to the major histocompatibility complex (MHC). Through careful primer design targeting to conserved regions of MHC class I sequences in birds, we successfully amplified a genomic fragment spanning exons 2–4 in three birds of prey. The identification of a highly conserved region within intron 2 allowed cross-amplifying complete exon 3 sequences in diurnal raptors, owls and New World vultures. We found evidence through PCR and cloning for 1–2 polymorphic class I loci, although this is almost certainly an underestimate. Inferences of diversifying selection in the kestrel MHC revealed that the two major regions of exon 3 exhibiting positive selection mostly agree with those described for the human HLA-A2 molecule. In contrast to passerines, where a high incidence of gene duplications and pseudogenes has been commonly documented, birds of prey emerge as nice model species for the investigation of the evolutionary significance and conservation implications of MHC diversity in vertebrates.     

Spectrum of MHC class II variability in Darwin's finches and their close relatives

The study describes >400 major histocompatibility complex (MHC) class II B exon 2 and 114 intron 2 sequences of 36 passerine bird species, 13 of which belong to the group of Darwin's finches (DFs) and the remaining 23 to close or more distant relatives of DFs in Central and South America. The data set is analyzed by a combination of judiciously selected statistical methods. The analysis reveals that reliable information concerning MHC organization, including the assignment of sequences to loci, and evolution, as well as the process of species divergence, can be obtained in the absence of genomic sequence data, if the analysis is taken several steps beyond the standard phylogenetic tree construction approach. The main findings of the present study are these: The MHC class II B region of the passerine birds is as elaborate in its organization, divergence, and genetic diversity as the MHC of the eutherian mammals, specifically the primates. Hence, the reported simplicity of the fowl MHC is an oddity. With the help of appropriate markers, the divergence of the MHC genes can be traced deep in the phylogeny of the bird taxa. Transspecies polymorphism is rampant at many of the bird MHC loci. In this respect, the DFs behave as if they were a single, genetically undifferentiated population. There is thus far no indication of alleles that could be considered species, genus, or even DF group specific. The implication of these findings is that DFs are in the midst of adaptive radiations, in which morphological differentiation into species is running ahead of genetic differentiation in genetic systems such as the MHC or the mitochondrial DNA. The radiations are so young that there has not been enough time to sort out polymorphisms at most of the loci among the morphologically differentiating species. These findings parallel those on Lake Victoria haplochromine fishes. Several of the DF MHC allelic lineages can be traced back to the MHC genes of the species Tiaris obscura, which we identified previously as the closest extant relative of DFs in continental America.     

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non‐model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non‐passerine birds. Non‐passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino‐acid residues showed marked discrepancies with peptide‐binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.     

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

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

Marsupial MHC class II beta genes are not orthologous to the eutherian beta gene families

The major histocompatibility complex (MHC) class II DRB, DQB, DPB, and DOB gene clusters are shared by different eutherian orders. Such an orthologous relationship is not seen between the beta genes of birds and eutherians. A high degree of uncertainty surrounds the evolutionary relationship of marsupial class II beta sequences with eutherian beta gene families. In particular, it has been suggested that marsupials utilize the DRB gene cluster. A cDNA encoding an MHC class II beta molecule was isolated from a brushtail possum mesenteric lymph node cDNA library. This clone is most similar to Macropus rufogriseus DBB. Our analysis suggests that all known marsupial beta-chain genes, excluding DMB, fall into two separate clades, which are distinct from the eutherian DRB, DQB, DPB, or DOB gene clusters. We recommend that the DAB and DBB nomenclature be reinstated. DAB and DBB orthologs are not present in eutherians. It appears that the marsupial and eutherian lineages have retained different gene clusters following gene duplication events early in mammalian evolution.     

Characterisation of class II B MHC genes from a ratite bird,the little spotted kiwi (<Emphasis Type="Italic">Apteryx owenii</Emphasis>)

Major histocompatibility complex (MHC) genes are important for vertebrate immune response and typically display high levels of diversity due to balancing selection from exposure to diverse pathogens. An understanding of the structure of the MHC region and diversity among functional MHC genes is critical to understanding the evolution of the MHC and species resilience to disease exposure. In this study, we characterise the structure and diversity of class II MHC genes in little spotted kiwi Apteryx owenii, a ratite bird representing the basal avian lineage (paleognaths). Results indicate that little spotted kiwi have a more complex MHC structure than that of other non-passerine birds, with at least five class II MHC genes, three of which are expressed and likely to be functional. Levels of MHC variation among little spotted kiwi are extremely low, with 13 birds assayed having nearly identical MHC genotypes (only two genotypes containing four alleles, three of which are fixed). These results suggest that recent genetic drift due to a species-wide bottleneck of at most seven birds has overwhelmed past selection for high MHC diversity in little spotted kiwi, potentially leaving the species highly susceptible to disease.     

Differential modes of MHC class IIB gene evolution in cichlid fishes

Cichlid fishes are emblematic models for the study of adaptive radiation, driven by natural and sexual selection. Parasite mediated selection is an important component in these processes, and the evolution of their immune system therefore merits special attention. In this study, light is shed on the phylogeny of the b family of cichlid major histocompatibility complex (MHC) class IIB genes. Full-length coding sequences were used to reconstruct phylogenies using criteria of maximum parsimony, maximum likelihood and Bayesian inference. All analyses suggest monophyly of the b family of cichlid MHC class IIB genes, although sequences of the cichlid sister taxa are currently not available. Two evolutionary lineages of these genes, respectively encompassing the recently defined genomic regions DBB-DEB-DFB and DCB-DDB, show highly contrasting levels of differentiation. To explore putative causes for these differences, exon 2 sequences were screened for variation in recombination rate and strength of selection. The more diversified lineage of cichlid MHC class IIB b genes was found to have higher levels of both recombination and selection. This is consistent with the observation in other taxa that recombination facilitates the horizontal spread of positively selected sites across MHC loci and hence contributes to fast sequence evolution. In contrast, the lineage that showed low diversification might either be under stabilizing selection or is evolutionary constrained by its low recombination rate. We speculate whether this lineage might include MHC genes with non-classical functions.     

MHC class I loci of the Bar-Headed goose (Anser indicus)

MHC class I proteins mediate functions in anti-pathogen defense. MHC diversity has already been investigated by many studies in model avian species, but here we chose the bar-headed goose, a worldwide migrant bird, as a non-model avian species. Sequences from exons encoding the peptide-binding region (PBR) of MHC class I molecules were isolated from liver genomic DNA, to investigate variation in these genes. These are the first MHC class I partial sequences of the bar-headed goose to be reported. A preliminary analysis suggests the presence of at least four MHC class I genes, which share great similarity with those of the goose and duck. A phylogenetic analysis of bar-headed goose, goose and duck MHC class I sequences using the NJ method supports the idea that they all cluster within the anseriforms clade.     

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.     

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.