Evolution by recombination and transspecies polymorphism in the MHC class I gene of Xenopus laevis

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

An amphibian emerging infectious disease (EID), chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), originated in Asia but primarily led to declines and extinctions in amphibian populations outside of Asia. Host major histocompatibility complex (MHC) molecules exhibit high polymorphism, and the evolution of MHC can be influenced by recombination and pathogens. Previous studies have indicated that host MHC class II is associated with Bd resistance. In this study, I conducted recombination and selection tests on functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) and an Australasian Bd-susceptible species (Litoria caerulea). Recombination at the same site was identified in both species, supporting the hypothesis that recombination contributes to MHC IIß1 diversity in amphibians. Positive selection was observed in MHC IIß1 alleles in both species. In L. caerulea, at least four amino acid sites were identified under significant positive selection in the MHC IIß1, whereas these sites were either negatively selected or conserved in B. gargarizans. This suggests these sites might be selected for Bd resistance. Hydrophobicity was detected in certain amino acid sites relating to Bd resistance, suggesting this physicochemical property may be a factor selected to counteract Bd infection. These findings of this study provide an evolutionary basis for understanding how amphibian MHC IIß1 may undergo selection in response to chytrid infection.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

Genetic diversity of major histocompatibility complex class I genes in Zootoca vivipara

The Major Histocompatibility Complex (MHC), as a family of highly polymorphic genes associated with immunity in the genome of the vertebrate, has become an important indicator for assessing the evolutionary potential of wildlife. In order to better protect Zootoca vivipara in the Greater Khingan Range and Lesser Khingan Range, to understand the genetic structure of Z. vivipara, and to explore the mechanism and phylogenetic relationship of the gene polymorphisms, the MHC molecular marker method was used to analyze Z. vivipara population. Forty-seven alleles were obtained from four populations. The four populations were highly polymorphic, rich in genetic information, and had significant genetic diversity. There were certain inbreeding phenomena. There was a high degree of genetic differentiation among populations, which was caused by genetic drift and natural selection. The sequence undergoes genetic duplication and recombination. The existence of trans-species polymorphism was found in the constructed phylogenetic tree. The present study provides a theoretical basis for species protection of Z. vivipara.

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora, Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1-3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52-57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.

The MHC Class Ia Genes in Chenfu's Treefrog (Zhangixalus chenfui) Evolved via Gene Duplication, Recombination, and Selection

Simple Summary

Amphibians, the first terrestrial vertebrates, provide materials for adaptive evolutionary studies, such as the evolution of the major histocompatibility complex (MHC). To date, various MHC evolutionary mechanisms have been identified in frogs, but more research is needed to determine the evolutionary mechanisms of the frog MHC. The main purpose of this study was to evaluate polymorphisms in the MHC class Ia genes of the Chenfu’s Treefrog. The MHC class Ia genes of the Chenfu’s Treefrog have high polymorphism. The mechanisms responsible for the formation of the polymorphisms include gene duplication, recombination, and selection.

Abstract

The molecular mechanisms underlying the evolution of adaptive immunity-related proteins can be deduced by a thorough examination of the major histocompatibility complex (MHC). Currently, in vertebrates, there is a relatively large amount of research on MHCs in mammals and birds. However, research related to amphibian MHC genes and knowledge about the evolutionary patterns is limited. This study aimed to isolate the MHC class I genes from Chenfu’s Treefrog (Zhangixalus chenfui) and reveal the underlying evolutionary processes. A total of 23 alleles spanning the coding region of MHC class Ia genes were identified in 13 individual samples. Multiple approaches were used to test and identify recombination from the 23 alleles. Amphibian MHC class Ia alleles, from NCBI, were used to construct the phylogenetic relationships in MEGA. Additionally, the partition strategy was adopted to construct phylogenetic relationships using MrBayes and MEGA. The sites of positive selection were identified by FEL, PAML, and MEME. In Chenfu’s Treefrog, we found that: (1) recombination usually takes place between whole exons of MHC class Ia genes; (2) there are at least 3 loci for MHC class Ia, and (3) the diversity of genes in MHC class Ia can be attributed to recombination, gene duplication, and positive selection. We characterized the evolutionary mechanisms underlying MHC class Ia genes in Chenfu’s Treefrog, and in so doing, broadened the knowledge of amphibian MHC systems.

Influence of immunogenetics, sex and body condition on the cutaneous microbial communities of two giant salamanders

The complex association between hosts and microbial symbionts requires the implementation of multiple approaches to evaluate variation in host physiology. Within amphibians, heterogeneity in immunogenetic traits and cutaneous microbiota is associated with variation in disease resistance. Ozark (Cryptobranchus alleganiensis bishopi) and eastern hellbenders (C. a. alleganiensis) provide a model system to assess variation in host traits and microbial communities. Ozark hellbenders have experienced declines throughout their range, are federally endangered and experience wound retardation that is absent in the eastern subspecies. Previous microbial investigations indicate differentiation in the composition of the skin microbiota of both hellbender subspecies, but it is not clear whether these patterns are concurrent with diversity in the major histocompatibility complex (MHC) genes. We characterized the MHC IIB and the skin microbiota of hellbenders in Missouri, where both subspecies co-occur though not sympatric. We compared the microbiota composition and MHC diversity between both subspecies and investigated whether individual-level MHC diversity, sex and body condition were associated with microbiota composition. Overall, MHC IIB diversity was lower in Ozark hellbenders compared to the eastern subspecies. Multivariate statistical comparisons identified microbiota differentiation between Ozark and eastern hellbenders. MHC IIB allele presence/absence, allele divergence, body composition and sex defined grouping of hellbender microbiotas within populations. Differentiation of the cutaneous microbiotas and MHC IIB genes between eastern and Ozark hellbenders suggests that differences exist in immunity between the two subspecies. This study demonstrates how simultaneous assessments of host genetic traits and microbiotas can inform patterns of microbial community structure in natural systems.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

Amphibians, as the first terrestrial vertebrates, offer a window into early major histocompatibility complex (MHC) evolution. We characterized the MHC class I of two Korean amphibians, the Asiatic toad (Bufo gargarizans) and the Japanese tree frog (Hyla japonica). We found at least four transcribed MHC class I (MHC I) loci, the highest number confirmed in any anuran to date. Furthermore, we identified MHC I transcripts in terrestrial adults, and possibly in aquatic larvae, of both species. We conducted a phylogenetic analysis based on MHC I sequence data and found that B. gargarizans and H. japonica cluster together in the superfamily Nobleobatrachia. We further identified three supertypes shared by the two species. Our results reveal substantial variation in the number of MHC I loci in anurans and suggest that certain supertypes have particular physiochemical properties that may confer pathogen resistance.

Historical gene flow and profound spatial genetic structure among golden pheasant populations suggested by multi-locus analysis

Major histocompatibility complex (MHC) is a good marker system for geographical genetics since they are functional genes in the immune system that are likely to affect the fitness of the individual, and the survival and evolutionary potential of a population in a changing environment. Golden pheasant (Chrysolophus pictus) is a wild Phasianidae distributed in central and north China. In this study, we used a locus-specific genotyping technique for MHC IIB genes of golden pheasant. Combining with microsatellites (simple sequence repeat, SSR) and mitochondrial DNA (mtDNA) D-loop region, we investigated the demographic history and illuminate genetic structure of this bird in detail. SYR (south of Yangtze river) - NYR (north of Yangtze river) lineages, separated by Yangtze River, were defined in genetic structure of MHC IIB. NYR was supposed as refuge during glacial period, suggested by diversity parameters and more ancient alleles in this region. Based on this hypothesis, there was gene flow from NYR to SYR, which was proved by three pieces of evidence: (1) distinct demographic histories of SYR (kept stable) and NYR (experienced expansion); (2) specific affiliation of LC in genetic structure of SSR and MHC genes; (3) significant gene flow from NYR to SYR. Moreover, we also found balancing selection by combination of three Grouping A2's regions (SC, QL and North) into one in Grouping B4 (NYR) and no pattern of isolation by distance (IBD) found in MHC IIB, whereas for SSR we found a relatively strong and significant IBD. Several mechanisms in the evolution of MHC IIB genes, including recombination, historically positive selection, trans-species evolution and concerted evolution, were shown by molecular and phylogenetic analysis. Overall these results suggest the Yangtze River was inferred to be a geological barrier for this avian and NYR might experience population expansion, which invaded into a neighboring region. This study contributes to the understanding of the effects of geographic features on contemporary patterns of genetic variation in the golden pheasant in China, and helps us to define the adaptive unite (AU) for this avian.

MHC class II β exon 2 variation in pardalotes (Pardalotidae) is shaped by selection, recombination and gene conversion

The high levels of polymorphism and allelic diversity which characterise genes in the major histocompatibility complex (MHC) are thought to be generated and maintained through the combined effects of different evolutionary processes. Here, we characterised exon 2 of the MHC class II β genes in two congeneric passerine species, the spotted (Pardalotus punctatus) and striated pardalote (Pardalotus striatus). We estimated the levels of allelic diversity and tested for signatures of recombination, gene conversion and balancing selection to determine if these processes have influenced MHC variation in the two species. Both species showed high levels of polymorphism and allelic diversity, as well as evidence of multiple gene loci and putative pseudogenes based on the presence of stop codons. We found higher levels of MHC diversity in the striated pardalote than the spotted pardalote, based on the levels of individual heterozygosity, sequence divergence and number of polymorphic sites. The observed differences may reflect variable selection pressure on the species, resulting from differences in patterns of movement among populations. We identified strong signatures of historical balancing selection, recombination and gene conversion at the sequence level, indicating that MHC variation in the two species has been shaped by a combination of processes.

Characterisation of major histocompatibility complex class I genes in Japanese Ranidae frogs

The major histocompatibility complex (MHC) is a key component of adaptive immunity in all jawed vertebrates, and understanding the evolutionary mechanisms that have shaped these genes in amphibians, one of the earliest terrestrial tetrapods, is important. We characterised MHC class I variation in three common Japanese Rana species (Rana japonica, Rana ornativentris and Rana tagoi tagoi) and identified a total of 60 variants from 21 individuals. We also found evolutionary signatures of gene duplication, recombination and balancing selection (including trans-species polymorphism), all of which drive increased MHC diversity. A unique feature of MHC class I from these three Ranidae species includes low synonymous differences per site (d_S) within species, which we attribute to a more recent diversification of these sequences or recent gene duplication. The resulting higher d_N/d_S ratio relative to other anurans studied could be related to stronger selection pressure at peptide binding sites. This is one of the first studies to investigate MHC in Japanese amphibians and permits further exploration of the polygenetic factors associated with resistance to infectious diseases.

Contrasting evolutionary histories of MHC class I and class II loci in grouse--effects of selection and gene conversion

Genes of the major histocompatibility complex (MHC) encode receptor molecules that are responsible for recognition of intracellular and extracellular pathogens (class I and class II genes, respectively) in vertebrates. Given the different roles of class I and II MHC genes, one might expect the strength of selection to differ between these two classes. Different selective pressures may also promote different rates of gene conversion at each class. Despite these predictions, surprisingly few studies have looked at differences between class I and II genes in terms of both selection and gene conversion. Here, we investigated the molecular evolution of MHC class I and II genes in five closely related species of prairie grouse (Centrocercus and Tympanuchus) that possess one class I and two class II loci. We found striking differences in the strength of balancing selection acting on MHC class I versus class II genes. More than half of the putative antigen-binding sites (ABS) of class II were under positive or episodic diversifying selection, compared with only 10% at class I. We also found that gene conversion had a stronger role in shaping the evolution of MHC class II than class I. Overall, the combination of strong positive (balancing) selection and frequent gene conversion has maintained higher diversity of MHC class II than class I in prairie grouse. This is one of the first studies clearly demonstrating that macroevolutionary mechanisms can act differently on genes involved in the immune response against intracellular and extracellular pathogens.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

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.

Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera

Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner (csd) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

Patterns of genetic differentiation at MHC class I genes and microsatellites identify conservation units in the giant panda

Background

Evaluating patterns of genetic variation is important to identify conservation units (i.e., evolutionarily significant units [ESUs], management units [MUs], and adaptive units [AUs]) in endangered species. While neutral markers could be used to infer population history, their application in the estimation of adaptive variation is limited. The capacity to adapt to various environments is vital for the long-term survival of endangered species. Hence, analysis of adaptive loci, such as the major histocompatibility complex (MHC) genes, is critical for conservation genetics studies. Here, we investigated 4 classical MHC class I genes (Aime-C, Aime-F, Aime-I, and Aime-L) and 8 microsatellites to infer patterns of genetic variation in the giant panda (Ailuropoda melanoleuca) and to further define conservation units.

Results

Overall, we identified 24 haplotypes (9 for Aime-C, 1 for Aime-F, 7 for Aime-I, and 7 for Aime-L) from 218 individuals obtained from 6 populations of giant panda. We found that the Xiaoxiangling population had the highest genetic variation at microsatellites among the 6 giant panda populations and higher genetic variation at Aime-MHC class I genes than other larger populations (Qinling, Qionglai, and Minshan populations). Differentiation index (FST)-based phylogenetic and Bayesian clustering analyses for Aime-MHC-I and microsatellite loci both supported that most populations were highly differentiated. The Qinling population was the most genetically differentiated.

Conclusions

The giant panda showed a relatively higher level of genetic diversity at MHC class I genes compared with endangered felids. Using all of the loci, we found that the 6 giant panda populations fell into 2 ESUs: Qinling and non-Qinling populations. We defined 3 MUs based on microsatellites: Qinling, Minshan-Qionglai, and Daxiangling-Xiaoxiangling-Liangshan. We also recommended 3 possible AUs based on MHC loci: Qinling, Minshan-Qionglai, and Daxiangling-Xiaoxiangling-Liangshan. Furthermore, we recommend that a captive breeding program be considered for the Qinling panda population.

Characteristics of MHC class I genes in house sparrows Passer domesticus as revealed by long cDNA transcripts and amplicon sequencing

In birds the major histocompatibility complex (MHC) organization differs both among and within orders; chickens Gallus gallus of the order Galliformes have a simple arrangement, while many songbirds of the order Passeriformes have a more complex arrangement with larger numbers of MHC class I and II genes. Chicken MHC genes are found at two independent loci, classical MHC-B and non-classical MHC-Y, whereas non-classical MHC genes are yet to be verified in passerines. Here we characterize MHC class I transcripts (α1 to α3 domain) and perform amplicon sequencing using a next-generation sequencing technique on exon 3 from house sparrow Passer domesticus (a passerine) families. Then we use phylogenetic, selection, and segregation analyses to gain a better understanding of the MHC class I organization. Trees based on the α1 and α2 domain revealed a distinct cluster with short terminal branches for transcripts with a 6-bp deletion. Interestingly, this cluster was not seen in the tree based on the α3 domain. 21 exon 3 sequences were verified in a single individual and the average numbers within an individual were nine and five for sequences with and without a 6-bp deletion, respectively. All individuals had exon 3 sequences with and without a 6-bp deletion. The sequences with a 6-bp deletion have many characteristics in common with non-classical MHC, e.g., highly conserved amino acid positions were substituted compared with the other alleles, low nucleotide diversity and just a single site was subject to positive selection. However, these alleles also have characteristics that suggest they could be classical, e.g., complete linkage and absence of a distinct cluster in a tree based on the α3 domain. Thus, we cannot determine for certain whether or not the alleles with a 6-bp deletion are non-classical based on our present data. Further analyses on segregation patterns of these alleles in combination with dating the 6-bp deletion through MHC characterization across the genus Passer may solve this matter in the future.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Selection and phylogenetics of salmonid MHC class I: wild brown trout (Salmo trutta) differ from a non-native introduced strain

We tested how variation at a gene of adaptive importance, MHC class I (UBA), in a wild, endemic Salmo trutta population compared to that in both a previously studied non-native S. trutta population and a co-habiting Salmo salar population (a sister species). High allelic diversity is observed and allelic divergence is much higher than that noted previously for co-habiting S. salar. Recombination was found to be important to population-level divergence. The α1 and α2 domains of UBA demonstrate ancient lineages but novel lineages are also identified at both domains in this work. We also find examples of recombination between UBA and the non-classical locus, ULA. Evidence for strong diversifying selection was found at a discrete suite of S. trutta UBA amino acid sites. The pattern was found to contrast with that found in re-analysed UBA data from an artificially stocked S. trutta population.

In silico peptide-binding predictions of passerine MHC class I reveal similarities across distantly related species, suggesting convergence on the level of protein function

The major histocompatibility complex (MHC) genes are the most polymorphic genes found in the vertebrate genome, and they encode proteins that play an essential role in the adaptive immune response. Many songbirds (passerines) have been shown to have a large number of transcribed MHC class I genes compared to most mammals. To elucidate the reason for this large number of genes, we compared 14 MHC class I alleles (α1-α3 domains), from great reed warbler, house sparrow and tree sparrow, via phylogenetic analysis, homology modelling and in silico peptide-binding predictions to investigate their functional and genetic relationships. We found more pronounced clustering of the MHC class I allomorphs (allele specific proteins) in regards to their function (peptide-binding specificities) compared to their genetic relationships (amino acid sequences), indicating that the high number of alleles is of functional significance. The MHC class I allomorphs from house sparrow and tree sparrow, species that diverged 10 million years ago (MYA), had overlapping peptide-binding specificities, and these similarities across species were also confirmed in phylogenetic analyses based on amino acid sequences. Notably, there were also overlapping peptide-binding specificities in the allomorphs from house sparrow and great reed warbler, although these species diverged 30 MYA. This overlap was not found in a tree based on amino acid sequences. Our interpretation is that convergent evolution on the level of the protein function, possibly driven by selection from shared pathogens, has resulted in allomorphs with similar peptide-binding repertoires, although trans-species evolution in combination with gene conversion cannot be ruled out.

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.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline MHC-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.

A new theory of MHC evolution: beyond selection on the immune genes

The major histocompatibility complex (MHC) is a dense region of immune genes with high levels of polymorphism, which are arranged in haplotype blocks. Traditional models of balancing selection (i.e. overdominance and negative frequency dependence) were developed to study the population genetics of single genes. However, the MHC is a multigene family surrounded by linked (non-neutral) polymorphisms, and not all of its features are well explained by these models. For example, (i) the high levels of polymorphism in small populations, (ii) the unexpectedly large genetic differentiation between populations, (iii) the shape of the allelic genealogy associated with trans-species evolution, and (iv) the close associations between particular MHC (human leucocyte antigen, HLA) haplotypes and the approximately 100 pathologies in humans. Here, I propose a new model of MHC evolution named Associative Balancing Complex evolution that can explain these phenomena. The model proposes that recessive deleterious mutations accumulate as a 'sheltered load' nearby MHC genes. These mutations can accumulate because (i) they are rarely expressed as homozygotes given the high MHC gene diversity and (ii) purifying selection is inefficient with low recombination rates (cf. Muller's ratchet). Once fixed, these mutations add to balancing selection and further reinforce linkage through epistatic selection against recombinants.

Comparative and developmental study of the immune system in Xenopus

Xenopus laevis is the model of choice for evolutionary, comparative, and developmental studies of immunity, and invaluable research tools including MHC-defined clones, inbred strains, cell lines, and monoclonal antibodies are available for these studies. Recent efforts to use Silurana (Xenopus) tropicalis for genetic analyses have led to the sequencing of the whole genome. Ongoing genome mapping and mutagenesis studies will provide a new dimension to the study of immunity. Here we review what is known about the immune system of X. laevis integrated with available genomic information from S. tropicalis. This review provides compelling evidence for the high degree of similarity and evolutionary conservation between Xenopus and mammalian immune systems. We propose to build a powerful and innovative comparative biomedical model based on modern genetic technologies that takes take advantage of X. laevis and S. tropicalis, as well as the whole Xenopus genus. Developmental Dynamics 238:1249-1270, 2009. (c) 2009 Wiley-Liss, Inc.

Correlated evolution of nucleotide substitution rates and allelic variation in Mhc-DRB lineages of primates

BACKGROUND

The major histocompatibility complex (MHC) is a key model of genetic polymorphism. Selection pressure by pathogens or other microevolutionary forces may result in a high rate of non-synonymous substitutions at the codons specifying the contact residues of the antigen binding sites (ABS), and the maintenance of extreme MHC allelic variation at the population/species level. Therefore, selection forces favouring MHC variability for any reason should cause a correlated evolution between substitution rates and allelic polymorphism. To investigate this prediction, we characterised nucleotide substitution rates and allelic polymorphism (i.e. the number of alleles detected in relation to the number of animals screened) of several Mhc class II DRB lineages in 46 primate species, and tested for a correlation between them.

RESULTS

First, we demonstrate that species-specific and lineage-specific evolutionary constraints favour species- and lineage-dependent substitution rate at the codons specifying the ABS contact residues (i.e. certain species and lineages can be characterised by high substitution rate, while others have low rate). Second, we show that although the degree of the non-synonymous substitution rate at the ABS contact residues was systematically higher than the degree of the synonymous substitution rate, these estimates were strongly correlated when we controlled for species-specific and lineage-specific effects, and also for the fact that different studies relied on different sample size. Such relationships between substitution rates of different types could even be extended to the non-contact residues of the molecule. Third, we provide statistical evidence that increased substitution rate along a MHC gene may lead to allelic variation, as a high substitution rate can be observed in those lineages in which many alleles are maintained. Fourth, we show that the detected patterns were independent of phylogenetic constraints. When we used phylogenetic models that control for similarity between species, due to common descent, and focused on variations within a single lineage (DRB1*03), the positive relationship between different substitution rates and allelic polymorphisms was still robust. Finally, we found the same effects to emerge in the analyses that eliminated within-species variation in MHC traits by using strictly single population-level studies. However, in a set of contrasting analyses, in which we focused on the non-functional DRB6 locus, the correlation between substitution rates and allelic variation was not prevalent.

CONCLUSION

Our results indicate that positive selection for the generation of allelic polymorphism acting on the functional part of the protein has consequences for the nucleotide substitution rate along the whole exon 2 sequence of the Mhc-DRB gene. Additionally, we proved that an increased substitution rate can promote allelic variation within lineages. Consequently, the evolution of different characteristics of genetic polymorphism is not independent.

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.

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.

Contrasting mode of evolution between the MHC class I genomic region and class II region in the three-spined stickleback (Gasterosteus aculeatus L.; Gasterosteidae: Teleostei)

Major histocompatibility complex (MHC) class I molecules display peptides on cell surfaces for subsequent T-cell recognition and are involved in the immune response against intracellular pathogens. In this study, a BAC library was created from a single three-spined stickleback and screened for clones containing MHC class I genes. In a 163.2-kb genomic sequence segment of a single clone, we identified three MHC class I genes in the same transcriptional orientation. Two class I genes are potentially expressed and functional. In one class I gene, the transmembrane region is missing and could therefore present a pseudogene. Alternatively, it presents a functional gene that encodes a soluble MHC class Ib molecule. Despite genomic similarities to the MHC class II region, which is characterized by interlocus recombination, we did not find any evidence for this kind of recombination in the class I genes. It thus seems that interlocus recombination may play a rather minor role in generating class I diversity in stickleback and that the class I region displays a higher genomic stability (i.e., lower local recombination rate). In addition, two non-MHC genes (Oct-2 beta and Na(+),K(+)-ATPasealpha3) have been identified in the analyzed class I region. The Oct-2 beta gene is a transcription factor that is expressed primarily in B lymphocytes, in activated T-cells, and in neuronal cells. The Na(+),K(+)-ATPasealpha3 gene is primarily expressed in the brain and heart and mediates catalytic activities. Both genes are located on the same linkage group together with the MHC class I genes in the zebra fish. In humans, however, homologues of Oct-2 beta and ATPasealpha3 lie outside the MHC region, which indicates that the concentration of immune genes found in mammalian genomes is a derived state.

Polymorphism, natural selection, and structural modeling of class Ia MHC in the African clawed frog (Xenopus laevis)

In the African clawed frog (Xenopus laevis), two deeply divergent allelic lineages of multiple genes of the class I MHC region have been discovered. For the MHC class I UAA locus, functional differences and the molecular basis for lineages maintenance are unknown. Alleles of linked class I region genes also exhibit strong disequilibrium with specific MHC alleles, but the underlying cause is not clear. We use MHC class Ia sequence data to estimate substitution rates and investigate structural differences between allelic lineages from protein models. Results indicate the operation of natural selection, and differences in the steric properties in the F pocket of the peptide-binding region among lineages. Variability in this pocket likely enables allelic lineages to bind very different sets of peptides and to interact differently with MHC chaperones in the endoplasmic reticulum. These results constitute evidence of the molecular evolutionary basis for 1) the maintenance of allelic lineages, 2) functional differences among lineages, and 3) strong linkage disequilibrium of allelic variants of class I region genes in X. laevis.