Species-specific evolution of class I MHC genes in iguanas (Order: Squamata; Subfamily: Iguaninae)

Over the last few decades, the major histocompatibility complex (MHC) has emerged as a model for understanding the influence of natural selection on genetic diversity in populations as well as for investigating the genetic basis of host resistance to pathogens. However, many vertebrate taxa remain underrepresented in the field of MHC research, preventing its application to studies of disease, evolution, and conservation genetics in these groups. This is particularly true for squamates, which are by far the most diversified order of non-avian reptiles but have not been the subject of any recent MHC studies. In this paper, we present MHC class I complementary DNA data from three squamate species in the subfamily Iguaninae (iguanas): the Galápagos marine iguana (Amblyrhynchus cristatus), the Galápagos land iguana (Conolophus subcristatus), and the green iguana (Iguana iguana). All sequences obtained are related to the few published class I genes from other squamates. There is evidence for multiple loci in each species, and the conserved alpha-3 domain appears to be evolving in a species-specific manner. Conversely, there is some indication of shared polymorphism between species in the peptide-binding alpha-1 and alpha-2 domains, suggesting that these two regions have different phylogenetic histories. The great similarity between alpha-3 sequences in marine iguanas in particular suggests that concerted evolution is acting to homogenize class I loci within species. However, while less likely, the data are also compatible with a birth and death model of evolution.     

MHC class I variation in a natural blue tit population (Cyanistes caeruleus)

The major histocompatibility complex (MHC) is central to the vertebrate immune system and its highly polymorphic genes are considered to influence several life-history traits of individuals. To characterize the MHC in a natural population of blue tits (Cyanistes caeruleus) we investigated the class I exon 3 diversity of more than 900 individuals. We designed two pairs of motif-specific primers that reliably amplify independent subsets of MHC alleles. Applying denaturing gradient gel electrophoresis (DGGE) we obtained 48 independently inherited units of unique band patterns (DGGE-haplogroups), which were validated in a segregation analysis within 105 families. In a second approach, we extensively sequenced 6 unrelated individuals to confirm that DGGE-haplogroup composition reflects individual allelic variation. The highest number of different DGGE-haplogroups in a single individual corresponded in 19 MHC exon 3 sequences, suggesting a minimum of 10 amplified MHC class I loci in the blue tit. In total, we identified 50 unique functional and 3 non-functional sequences. Functional sequences showed high levels of recombination and strong positive selection in the antigen binding region, whereas nucleotide diversity was comparatively low in the range of all passerine species. Finally, in a phylogenetic comparison of passerine MHC class I exon 3 sequences we discuss conflicting evolutionary signals possibly due to recent gene duplication, recombination events and concerted evolution. Our results indicate that the described method is suitable to effectively explore the MHC diversity and its ecological impacts in blue tits in future studies.     

Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae)

In contrast to mammals, the evolution of MHC genes in birds appears to be characterized by high rates of gene duplication and concerted evolution. To further our understanding of the evolution of passerine MHC genes, we have isolated class II B sequences from two species of New Zealand robins, the South Island robin (Petroica australis australis), and the endangered Chatham Island black robin (Petroica traversi). Using an RT-PCR based approach we isolated four transcribed class II B MHC sequences from the black robin, and eight sequences from the South Island robin. RFLP analysis indicated that all class II B loci were contained within a single linkage group. Analysis of 3-untranslated region sequences enabled putative orthologous loci to be identified in the two species, and indicated that multiple rounds of gene duplication have occurred within the MHC of New Zealand robins. The orthologous relationships are not retained within the coding region of the gene, instead the sequences group within species. A number of putative gene conversion events were identified across the length of our sequences that may account for this. Exon 2 sequences are highly diverse and appear to have diverged under balancing selection. It is also possible that gene conversion involving short stretches of sequence within exon 2 adds to this diversity. Our study is the first report of putative orthologous MHC loci in passerines, and provides further evidence for the importance of gene duplication and gene conversion in the evolution of the passerine MHC.Nucleotide sequence data reported in this paper are available in the GenBank database under the accession numbers AY258333–AY258335, AY428561–AY428570, and AY530534–AY530535     

Evidence of gene orthology and trans‐species polymorphism,but not of parallel evolution,despite high levels of concerted evolution in the major histocompatibility complex of flamingo species

The major histocompatibility complex (MHC) is a cornerstone in the study of adaptive genetic diversity. Intriguingly, highly polymorphic MHC sequences are often not more similar within species than between closely related species. Divergent selection of gene duplicates, balancing selection maintaining trans‐species polymorphism (TSP) that predate speciation and parallel evolution of species sharing similar selection pressures can all lead to higher sequence similarity between species. In contrast, high rates of concerted evolution increase sequence similarity of duplicated loci within species. Assessing these evolutionary models remains difficult as relatedness and ecological similarities are often confounded. As sympatric species of flamingos are more distantly related than allopatric species, flamingos represent an ideal model to disentangle these evolutionary models. We characterized MHC Class I exon 3, Class IIB exon 2 and exon 3 of the six extant flamingo species. We found up to six MHC Class I loci and two MHC Class IIB loci. As all six species shared the same number of MHC Class IIB loci, duplication appears to predate flamingo speciation. However, the high rate of concerted evolution has prevented the divergence of duplicated loci. We found high sequence similarity between all species regardless of codon position. The latter is consistent with balancing selection maintaining TSP, as under this mechanism amino acid sites under pathogen‐mediated selection should be characterized by fewer synonymous codons (due to their common ancestry) than under parallel evolution. Overall, balancing selection maintaining TSP appears to result in high MHC similarity between species regardless of species relatedness and geographical distribution.     

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.     

ERVK9, transposons and the evolution of MHC class I duplicons within the alpha-block of the human and chimpanzee

The genomic sequences within the alpha-block (approximately 288-310 kb) of the human and chimpanzee MHC class I region contains ten MHC class I genes and three MIC gene fragments grouped together within alternating duplicated genomic segments or duplicons. In this study, the chimpanzee and human genomic sequences were analyzed in order to determine whether the remnants of the ERVK9 and other retrotransposon sequences are useful genomic markers for reconstructing the evolutionary history of the duplicated MHC gene families within the alpha-block. A variety of genes, pseudogenes, autologous DNA transposons and retrotransposons such as Alu and ERVK9 were used to categorize the ten duplicons into four distinct structural groups. The phylogenetic relationship of the ten duplicons was examined by using the neighbour joining method to analyze transposon sequence topologies of selected Alu members, LTR16B and Charlie9. On the basis of these structural groups and the phylogeny of the duplicated transposon sequences, a duplication model was reconstructed involving four multipartite tandem duplication steps to explain the organization and evolution of the ten duplicons within the alpha-block of the chimpanzee and human. The phylogenetic analysis and inferred duplication history suggests that the Patr/HLA-F was the first MHC class I gene to have been fixed and not required as a precursor for further duplication within the alpha-block of the ancestral species.     

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.     

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.     

Extremely high MHC class I variation in a population of a long-distance migrant,the Scarlet Rosefinch (<Emphasis Type="Italic">Carpodacus erythrinus</Emphasis>)

Although the number of studies focusing on the major histocompatibility complex (MHC) in non-model vertebrates is increasing, results are often contradictory, and the structure of MHC is still poorly understood in wild species. Here, we describe the structure and diversity of exon 3 of MHC class I in a passerine bird, the Scarlet Rosefinch (Carpodacus erythrinus). Using capillary electrophoresis single-strand conformation polymorphism, we identified 82 different MHC class I variants in one Rosefinch population nesting at one site in the Czech Republic. Thus far, this is the highest intra-populational MHC class I variation observed in birds. We have not found support for ‘minimal essential’ MHC in this species since individuals exhibited between three and nine different exon 3 sequences, indicating that there may be at least five amplified MHC class I genes. By cloning, we obtained and analysed 29 exon sequences and found that all of them could be translated into potentially functional proteins. We also show that strong positive selection appears to be acting mainly, but not only, on previously described antigen-binding sites in MHC class I genes. Furthermore, our results indicate that recombination has played an important role in generating genetic diversity of these genes in the Scarlet Rosefinch; we discuss the significance of this extremely high genetic diversity in light of the life history traits of this species, such as long-distance migration. The sequence data described in this paper are accessible in GenBank data library under accession numbers FJ392762–FJ392790.     

MHC variation in birds and reptiles

The major histocompatibility complex (MHC) has been studied in a multitude of mammals by now, but much less is known about its organisation and variation in other vertebrate species. The mammalian MHC is organised as a single gene cluster, but recent studies on birds suggest that this paradigm of MHC organisation has to be supplemented. The domestic chicken thus possesses two separate gene clusters which both contain MHC class I and class II B genes, and we have shown that the ring‐necked pheasant Phasianus colchicus also has two unlinked clusters of class II B genes. We are studying the effect of the MHC on mate choice, survival and reproductive success in natural populations of birds and reptiles. For this reason, we are developing DNA techniques to determine the animals' MHC genotype. The amplification of the hypervariable exon 3 of the class I gene from songbirds and reptiles has provided us with species specific probes that can be used in Southern blot analysis. The first results indicate very extensive variation in all studied species, that is starlings Sturnus vulgaris, great reed warblers Acrocephalus arundinaceus and water pythons Liasis fuscus. The restriction fragment length polymorphism (RFLP) analysis also suggests that the number of MHC genes is significantly larger in these species than in pheasants and domestic chickens. This revised version was published online in July 2006 with corrections to the Cover Date.     

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat‐specific three‐ or five‐amino‐acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide‐binding groove, as previously assayed by a bat MHCI gene with a three‐amino‐acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon‐specific manner.     

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.     

Proteasome, transporter associated with antigen processing, and class I genes in the nurse shark Ginglymostoma cirratum: evidence for a stable class I region and MHC haplotype lineages.

Cartilaginous fish (e.g., sharks) are derived from the oldest vertebrate ancestor having an adaptive immune system, and thus are key models for examining MHC evolution. Previously, family studies in two shark species showed that classical class I (UAA) and class II genes are genetically linked. In this study, we show that proteasome genes LMP2 and LMP7, shark-specific LMP7-like, and the TAP1/2 genes are linked to class I/II. Functional LMP7 and LMP7-like genes, as well as multiple LMP2 genes or gene fragments, are found only in some sharks, suggesting that different sets of peptides might be generated depending upon inherited MHC haplotypes. Cosmid clones bearing the MHC-linked classical class I genes were isolated and shown to contain proteasome gene fragments. A non-MHC-linked LMP7 gene also was identified on another cosmid, but only two exons of this gene were detected, closely linked to a class I pseudogene (UAA-NC2); this region probably resulted from a recent duplication and translocation from the functional MHC. Tight linkage of proteasome and class I genes, in comparison with gene organizations of other vertebrates, suggests a primordial MHC organization. Another nonclassical class I gene (UAA-NC1) was detected that is linked neither to MHC nor to UAA-NC2; its high level of sequence similarity to UAA suggests that UAA-NC1 also was recently derived from UAA and translocated from MHC. These data further support the principle of a primordial class I region with few class I genes. Finally, multiple paternities in one family were demonstrated, with potential segregation distortions.     

Locus specificity of polymorphic alleles and evolution by a birth-and-death process in mammalian MHC genes

We have conducted an extensive phylogenetic analysis of polymorphic alleles from human and mouse major histocompatibility complex (MHC) class I and class II genes. The phylogenetic tree obtained for 212 complete human class I allele sequences (HLA-A, -B, and -C) has shown that all alleles from the same locus form a single cluster, which is highly supported by bootstrap values, except for one HLA-B allele (HLA-B*7301). Mouse MHC class I loci did not show locus-specific clusters of polymorphic alleles. This was considered to be because of either interlocus genetic exchange or the confusing designation of loci in different haplotypes at the present time. The locus specificity of polymorphic alleles was also observed in human and mouse MHC class II loci. It was therefore concluded that interlocus recombination or gene conversion is not very important for generating MHC diversity, with a possible exception of mouse class I loci. According to the phylogenetic trees of complete coding sequences, we classified human MHC class I (HLA-A, -B, and -C) and class II (DRB1) alleles into three to five major allelic lineages (groups), which were monophyletic with high bootstrap values. Most of these allelic groups remained unchanged even in phylogenetic trees based on individual exons, though this does not exclude the possibility of intralocus recombination involving short DNA segments. These results, together with the previous observation that MHC loci are subject to frequent duplication and deletion, as well as to balancing selection, indicate that MHC evolution in mammals is in agreement with the birth-and-death model of evolution, rather than with the model of concerted evolution.     

Generation and maintenance of diversity in the cattle MHC class I region

Major histocompatibility complex (MHC) class I genes play a crucial role in the immune defence against intracellular pathogens. An important evolutionary strategy is to generate and maintain a high level of diversity in these genes. Humans express three highly polymorphic classical MHC class I genes (HLA-A, HLA-B and HLA-C). In contrast, some species, for example rat and rhesus macaque, maintain diversity by generation of haplotypes that vary considerably with regard to the number and combination of transcribed genes. Cattle appear to use both strategies. We show that various combinations of six apparently classical genes, three of which are highly polymorphic, are transcribed on different haplotypes. Although additional sequences were identified in both cDNA and gDNA, it was not possible to assign them to any of these defined genes. Most were highly divergent or were non-classical class I genes. Thus, we found little evidence for frequent duplication and deletion of classical class I genes as reported in some other species. However, the maintenance of class I diversity in cattle may involve limited gene shuffling and deletion, possibly as a result of unequal crossing-over within the class I region.The first two authors made an equal contribution to this work.     

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next‐generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans‐species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long‐term survival of the species.     

Molecular evolution of Dmrt1 accompanies change of sex-determining mechanisms in reptilia

In reptiles, sex-determining mechanisms have evolved repeatedly and reversibly between genotypic and temperature-dependent sex determination. The gene Dmrt1 directs male determination in chicken (and presumably other birds), and regulates sex differentiation in animals as distantly related as fruit flies, nematodes and humans. Here, we show a consistent molecular difference in Dmrt1 between reptiles with genotypic and temperature-dependent sex determination. Among 34 non-avian reptiles, a convergently evolved pair of amino acids encoded by sequence within exon 2 near the DM-binding domain of Dmrt1 distinguishes species with either type of sex determination. We suggest that this amino acid shift accompanied the evolution of genotypic sex determination from an ancestral condition of temperature-dependent sex determination at least three times among reptiles, as evident in turtles, birds and squamates. This novel hypothesis describes the evolution of sex-determining mechanisms as turnover events accompanied by one or two small mutations.