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

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.

Individual copy number variation and extensive diversity between major MHC-DAB1 allelic lineages in the European bitterling

Polymorphism of the major histocompatibility complex (MHC), DAB1 gene was characterized for the first time in the European bitterling (Rhodeus amarus), a freshwater fish employed in studies of host-parasite coevolution and mate choice, taking advantage of newly designed primers coupled with high-throughput amplicon sequencing. Across 221 genotyped individuals, we detected 1–4 variants per fish, with 28% individuals possessing 3–4 variants. We identified 36 DAB1 variants, and they showed high sequence diversity mostly located within predicted antigen-binding sites, and both global and codon-specific excess of non-synonymous mutations. Despite deep divergence between two major allelic lineages, functional diversity was surprisingly low (3 supertypes). Overall, these findings suggest the role of positive and balancing selection in promotion and long-time maintenance of DAB1 polymorphism. Further investigations will clarify the role of pathogen-mediated selection to drive the evolution of DAB1 variation.

Diversity of MHC IIB genes and parasitism in hybrids of evolutionarily divergent cyprinoid species indicate heterosis advantage

The genes of the major histocompatibility complex (MHC) are an essential component of the vertebrate immune system and MHC genotypes may determine individual susceptibility to parasite infection. In the wild, selection that favors MHC variability can create situations in which interspecies hybrids experience a survival advantage. In a wild system of two naturally hybridizing leuciscid fish, we assessed MHC IIB genetic variability and its potential relationships to hosts' ectoparasite communities. High proportions of MHC alleles and parasites were species-specific. Strong positive selection at specific MHC codons was detected in both species and hybrids. MHC allele expression in hybrids was slightly biased towards the maternal species. Controlling for a strong seasonal effect on parasite communities, we found no clear associations between host-specific parasites and MHC alleles or MHC supertypes. Hybrids shared more MHC alleles with the more MHC-diverse parental species, but expressed intermediate numbers of MHC alleles and positively selected sites. Hybrids carried significantly fewer ectoparasites than either parent species, suggesting a hybrid advantage via potential heterosis.

Spatiotemporal adaptive evolution of an MHC immune gene in a frog-fungus disease system

Genetic diversity of major histocompatibility complex (MHC) genes is linked to reduced pathogen susceptibility in amphibians, but few studies also examine broad spatial and temporal patterns of MHC and neutral genetic diversity. Here, we characterized range-wide MHC diversity in the Northern leopard frog, Rana pipiens, a species found throughout North America that is experiencing disease-related declines. We used previously sequenced neutral markers (mitochondrial DNA and microsatellites), sequenced an expressed MHC class IIß gene fragment, and measured infection prevalence and intensity of the global fungal pathogen Batrachochytrium dendrobatidis (Bd) across 14 populations. Four populations were sampled across two decades, enabling temporal comparisons of selection and demography. We recovered 37 unique MHC alleles, including 17 that were shared across populations. Phylogenetic and population genetic patterns between MHC and neutral markers were incongruent, and five MHC codon positions associated with peptide binding were under positive selection. MHC heterozygosity, but not neutral marker heterozygosity, was a significant factor explaining spatial patterns of Bd prevalence, whereas only environmental variables predicted Bd intensity. MHC allelic richness (AR) decreased significantly over time but microsatellite-based AR did not, highlighting a loss of functional immunogenetic diversity that may be associated with Bd selective pressures. MHC supertype 4 was significantly associated with an elevated risk of Bd infection, whereas one supertype 2 allele was associated with a nearly significant reduced risk of Bd. Taken together, these results provide evidence that positive selection contributes to MHC class IIß evolution in R. pipiens and suggest that functional MHC differences across populations may contribute to disease adaptation.

Polymorphism of MHC class IIB in an acheilognathid species, Rhodeus sinensis shaped by historical selection and recombination

BACKGROUND

Rhodeus sinensis is a bitterling species occurring throughout the numerous freshwater systems on the East Asia. Here, we analyzed the diversity of the MHC class IIB (DAB) genes from this species, which may offer meaningful insights into evolutionary processes in this species as well as other bitterlings.

RESULTS

Using cDNA and gDNA samples from 50 individuals, we discovered classical 140 allelic sequences that could be allocated into either DAB1 (Rhsi-DAB1) or DAB3 (Rhsi-DAB3). DAB sequences completely lacking the intron, but identical or similar to Rhsi-DAB1, were also discovered from our gDNA samples, and this intron loss likely originated from the retrotransposition events of processed mDNA. The β1 domain was the most polymorphic in both Rhsi-DAB1 and -DAB3. Putative peptide biding residues (PBRs) in Rhsi-DAB1, but not in Rhsi-DAB3, exhibited a significant dN/dS, presumably indicating that different selection pressures have acted on those two DABs. Recombination between different alleles seemed to have contributed to the increase of diversity in Rhsi-DABs. Upon phylogenetic analysis, Rhsi-DAB1 and -DAB3 formed independent clusters. Several alleles from other species of Cypriniformes were embedded in the clade of Rhsi-DAB1, whereas Rhsi-DAB3 clustered with alleles from the wider range of taxa (Cyprinodontiformes), indicating that these two Rhsi-DABs have taken different historical paths.

CONCLUSIONS

A great deal of MHC class IIB allelic diversity was found in R. sinensis, and gene duplication, selection and recombination may have contributed to this diversity. Based on our data, it is presumed that such historical processes have commonly or differently acted on the polymorphism of Rhsi-DAB1 and -DAB3.

MHC class IIA polymorphisms and their association with resistance-susceptibility to Streptococcus agalactiae in Nile tilapia, Oreochromis niloticus

The association between major histocompatibility complex (MHC) class IIA polymorphisms and the severity of infection by Streptococcus agalactiae was investigated using 40 susceptible and 40 resistant individuals of Nile tilapia Oreochromis niloticus. Twenty-five alleles were identified from 80 individuals, which belong to 22 major allele types. High polymorphism of mhcIIa gene and at least two loci were discovered in O. niloticus. In peptide-binding region (PBR) and non-PBR, the ratio of nonsynonymous substitution (dN) to synonymous substitution (dS) was 1.294 (>1) and 1.240 (>1), suggesting that the loci are evolving under positive balancing selection. Association analysis showed that the allele, orni-daa*0501, was significantly associated with resistance to S. agalactiae, while the alleles, orni-daa*1101, orni-daa*1301, orni-daa*1401 and orni-daa*1201, were associated with susceptibility to S. agalactiae. To confirm these correlations, another independent challenge experiment was performed in the Huizhou population of the O. niloticus. The frequency distribution showed that the orni-daa*1101 allele was significantly more frequent in the Huizhou-Susceptible group (HZ-SG) than in the Huizhou-Resistant group (HZ-RG) (P < 0.05), which was consistent with the first challenge. However, orni-daa*0501 did not present in HZ-SG and HZ-RG and the distribution frequencies of the orni-daa*1201, orni-daa*1301 and orni-daa*1401 alleles were not significantly more frequent in HZ-SG than in HZ-RG. These results indicate that the orni-daa*1101 allele confers susceptibility to S. agalactia infection. These results suggest that the diversity of exon 2 of mcaIIa alleles could be used to explore the association between disease susceptibility or resistance and the multiformity of mcaIIa and to achieve the molecular-assisted selection of O. niloticus with enhanced disease resistance.

High polymorphism in MHC-DRB genes in golden snub-nosed monkeys reveals balancing selection in small, isolated populations

Background

Maintaining variation in immune genes, such as those of the major histocompatibility complex (MHC), is important for individuals in small, isolated populations to resist pathogens and parasites. The golden snub-nosed monkey (Rhinopithecus roxellana), an endangered primate endemic to China, has experienced a rapid reduction in numbers and severe population fragmentation over recent years. For this study, we measured the DRB diversity among 122 monkeys from three populations in the Qinling Mountains, and estimated the relative importance of different agents of selection in maintaining variation of DRB genes.

Results

We identified a total of 19 DRB sequences, in which five alleles were novel. We found high DRB variation in R. roxellana and three branches of evidence suggesting that balancing selection has contributed to maintaining MHC polymorphism over the long term in this species: i) different patterns of both genetic diversity and population differentiation were detected at MHC and neutral markers; ii) an excess of non-synonymous substitutions compared to synonymous substitutions at antigen binding sites, and maximum-likelihood-based random-site models, showed significant positive selection; and iii) phylogenetic analyses revealed a pattern of trans-species evolution for DRB genes.

Conclusions

High levels of DRB diversity in these R. roxellana populations may reflect strong selection pressure in this species. Patterns of genetic diversity and population differentiation, positive selection, as well as trans-species evolution, suggest that pathogen-mediated balancing selection has contributed to maintaining MHC polymorphism in R. roxellana over the long term. This study furthers our understanding of the role pathogen-mediated balancing selection has in maintaining variation in MHC genes in small and fragmented populations of free-ranging vertebrates.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1148-7) contains supplementary material, which is available to authorized users.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Extensive shared polymorphism at non-MHC immune genes in recently diverged North American prairie grouse

Gene polymorphisms shared between recently diverged species are thought to be widespread and most commonly reflect introgression from hybridization or retention of ancestral polymorphism through incomplete lineage sorting. Shared genetic diversity resulting from incomplete lineage sorting is usually maintained for a relatively short period of time, but under strong balancing selection it may persist for millions of years beyond species divergence (balanced trans-species polymorphism), as in the case of the major histocompatibility complex (MHC) genes. However, balancing selection is much less likely to act on non-MHC immune genes. The aim of this study was to investigate the patterns of shared polymorphism and selection at non-MHC immune genes in five grouse species from Centrocercus and Tympanuchus genera. For this purpose, we genotyped five non-MHC immune genes that do not interact directly with pathogens, but are involved in signaling and regulate immune cell growth. In contrast to previous studies with MHC, we found no evidence for balancing selection or balanced trans-species polymorphism among the non-MHC immune genes. No haplotypes were shared between genera and in most cases more similar allelic variants sorted by genus. Between species within genera, however, we found extensive shared polymorphism, which was most likely attributable to introgression or incomplete lineage sorting following recent divergence and large ancestral effective population size (i.e., weak genetic drift). Our study suggests that North American prairie grouse may have attained relatively low degree of reciprocal monophyly at nuclear loci and reinforces the rarity of balancing selection in non-MHC immune genes.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Molecular characterization of MHC class IIB genes of sympatric Neotropical cichlids

BACKGROUND

The Major Histocompatibility Complex (MHC) is a key component of the adaptive immune system of all vertebrates and consists of the most polymorphic genes known to date. Due to this complexity, however, MHC remains to be characterized in many species including any Neotropical cichlid fish. Neotropical crater lake cichlids are ideal models to study evolutionary processes as they display one of the most convincing examples of sympatric and repeated parallel radiation events within and among isolated crater lakes.

RESULTS

Here, we characterized the genes of MHC class IIB chain of the Midas cichlid species complex (Amphilophus cf. citrinellus) including fish from five lakes in Nicaragua. We designed 19 new specific primers anchored in a stepwise fashion in order to detect all alleles present. We obtained 866 genomic DNA (gDNA) sequences from thirteen individuals and 756 additional sequences from complementary DNA (cDNA) of seven of those individuals. We identified 69 distinct alleles with up to 25 alleles per individual. We also found considerable intron length variation and mismatches of alleles detected in cDNA and gDNA suggesting that some loci have undergone pseudogenization. Lastly, we created a model of protein structure homology for each allele and identified their key structural components.

CONCLUSIONS

Overall, the Midas cichlid has one of the most diverse repertoires of MHC class IIB genes known, which could serve as a powerful tool to elucidate the process of divergent radiations, colonization and speciation in sympatry.

Early Duplication of a Single MHC IIB Locus Prior to the Passerine Radiations

A key characteristic of MHC genes is the persistence of allelic lineages over macroevolutionary periods, often through multiple speciation events. This phenomenon, known as trans-species polymorphism (TSP), is well documented in several major taxonomic groups, but has less frequently been observed in birds. The order Passeriformes is arguably the most successful terrestrial vertebrate order in terms of diversity of species and ecological range, but the reasons for this success remain unclear. Passerines exhibit the most highly duplicated MHC genes of any major vertebrate taxonomic group, which may generate increased immune response relative to other avian orders with fewer MHC loci. Here, we describe phylogenetic patterns of the MHC IIB in the passerine family Corvidae. Our results indicate wide-spread TSP within this family, with at least four supported MHC IIB allelic lineages that predate speciation by many millions of years. Markov chain Monte Carlo simulations indicate that divergence of these lineages occurred near the time of the divergence of the Passeriformes and other avian orders. We suggest that the current MHC diversity observed in passerines is due in part to the multiple duplication of a single MHC locus, DAB1, early in passerine evolution and that subsequent duplications of these paralogues have contributed to the enormous success of this order by increasing their ability to recognize and mount immune responses to novel pathogens.

Adaptive major histocompatibility complex (MHC) and neutral genetic variation in two native Baltic Sea fishes (perch Perca fluviatilis and zander Sander lucioperca) with comparisons to an introduced and disease susceptible population in Australia (P. fluviatilis): assessing the risk of disease epidemics

This study assessed the major histocompatibility complex (MHC) and neutral genetic variation and structure in two percid species, perch Perca fluviatilis and zander Sander lucioperca, in a unique brackish ecosystem, the Baltic Sea. In addition, to assess the importance of MHC diversity to disease susceptibility in these populations, comparisons were made to an introduced, disease susceptible, P. fluviatilis population in Australia. Eighty-three MHC class II B exon 2 variants were amplified: 71 variants from 92 P. fluviatilis samples, and 12 variants from 82 S. lucioperca samples. Microsatellite and MHC data revealed strong spatial genetic structure in S. lucioperca, but not P. fluviatilis, across the Baltic Sea. Both microsatellite and MHC data showed higher levels of genetic diversity in P. fluviatilis from the Baltic Sea compared to Australia, which may have facilitated the spread of an endemic virus, EHNV in the Australian population. The relatively high levels of genetic variation in the Baltic Sea populations, together with spatial genetic structure, however, suggest that there currently seems to be little risk of disease epidemics in this system. To ensure this remains the case in the face of ongoing environmental changes, fisheries and habitat disturbance, the conservation of local-scale genetic variation is recommended.

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.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Characterization of MHC class IIB for four endangered Australian freshwater fishes obtained from ecologically divergent populations

Genetic diversity is an essential aspect of species viability, and assessments of neutral genetic diversity are regularly implemented in captive breeding and conservation programs. Despite their importance, information from adaptive markers is rarely included in such programs. A promising marker of significance in fitness and adaptive potential is the major histocompatibility complex (MHC), a key component of the adaptive immune system. Populations of Australian freshwater fishes are generally declining in numbers due to human impacts and the introduction of exotic species, a scenario of particular concern for members of the family Percichthyidae, several of which are listed as nationally vulnerable or endangered, and hence subject to management plans, captive breeding, and restoration plans. We used a next-generation sequencing approach to characterize the MHC IIB locus and provide a conservative description of its levels of diversity in four endangered percichthyids: Gadopsis marmoratus, Macquaria australasica, Nannoperca australis, and Nannoperca obscura. Evidence is presented for a duplicated MHC IIB locus, positively selected sites and recombination of MHC alleles. Relatively moderate levels of diversity were detected in the four species, as well as in different ecotypes within each species. Phylogenetic analyses revealed genus specific clustering of alleles and no allele sharing among species. There were also no shared alleles observed between two ecotypes within G. marmoratus and within M. australasica, which might be indicative of ecologically-driven divergence and/or long divergence times. This represents the first characterization and assessment of MHC diversity for Percichthyidae, and also for Australian freshwater fishes in general, providing key genetic resources for a vertebrate group of increasing conservation concern.

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.

Loss of genetic diversity at an MHC locus in the endangered Tokyo bitterling Tanakia tanago (Teleostei: Cyprinidae)

Genetic diversity at a major histocompatibility complex (MHC) class II B gene was examined for two wild and three captive populations of the endangered Tokyo bitterling Tanakia tanago. A specific primer set was first developed to amplify the MHC II B exon 2 locus. Using single strand conformation polymorphism (SSCP) and sequencing analysis, 16 DAB3 alleles were detected with 56 nucleotide substitutions in the 276-bp region. In the putative antigen-binding sites of exon 2, the rate of nonsynonymous substitutions was significantly higher than that of synonymous substitutions (dN/dS = 2.79), indicating positive selection on the retention of polymorphism. The population from the Handa Natural Habitat Conservation Area and that from the Tone River system exhibited low variation (one and three alleles, respectively), whereas the captive population that originated from a mix of three distinct populations had the highest amounts of variation (14 alleles). The levels of heterozygosity at the MHC varied considerably among populations and showed significant correlations with those at putative neutral microsatellite markers, suggesting that genetic drift following a bottleneck has affected MHC variability in some populations. Comparisons between endangered and non-endangered fish species in previous reports and the present results indicate that the number of MHC alleles per population is on average 70% lower in endangered species than non-endangered species. Considering the functional consequence of this locus, attention should be paid to captive and wild endangered fish populations in terms of further loss of MHC alleles.

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

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

MH II-DAB gene expression in grass carp Ctenopharyngodon idella (Valenciennes) after infection with the ciliate parasite, Ichthyophthirius multifiliis

The grass carp, Ctenopharyngodon idella (Valenciennes), is one of the most extensively aquacultured freshwater fish in China. However, because of the lack of effective control measures and the high-density culture environment, considerable economic losses are caused by infection of C. idella with the parasitic ciliate, Ichthyophthirius multifiliis. The major histocompatibility (MH) DAB gene belongs to antigen-presented genes in the class II genomic region, which is associated with parasite resistance. To understand the relationship of the DAB gene with I. multifiliis infection in grass carp, the expression profiles of MH II-DAB were studied in tissues using real-time quantitative polymerase chain reaction. The results showed that expression of the MH II-DAB gene was up-regulated in head kidney after I. multifiliis infection, and the expression peak appeared earlier in the study (case) group than in the control group. The obvious up-regulation peak of MH II-DAB gene was found at days 2 and 4 in skin; at 12 h to day 4 in spleen; at 12 h and days 1 and 6 in gill; and at day 10 in blood, whereas the MH II-DAB gene was down-regulated in liver and intestines after I. multifiliis infection. These results have implications for better understanding C. idella resistance to I. multifiliis infection.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1'665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context. Using next-generation sequencing, the present manuscript identifies the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of a cyprinid fish: the European minnow (Phoxinus phoxinus). We highlight that the relative importance of neutral versus adaptive processes in shaping immune competence may differ between duplicates as a consequence of alternative selective regimes or different genomic contexts.

Genomic Porosity between Invasive Chondrostoma nasus and Endangered Endemic Parachondrostoma toxostoma (Cyprinidae): The Evolution of MHC IIB Genes

Two cyprinid species, Parachondrostoma toxostoma, an endemic threatened species, and Chondrostoma nasus, an invasive species, live in sympatry in southern France and form two sympatric zones where the presence of intergeneric hybrids is reported. To estimate the potential threat to endemic species linked to the introduction of invasive species, we focused on the DAB genes (functional MHC IIB genes) because of their adaptive significance and role in parasite resistance. More specifically, we investigated (1) the variability of MHC IIB genes, (2) the selection pattern shaping MHC polymorphism, and (3) the extent to which trans-species evolution and intergeneric hybridization affect MHC polymorphism. In sympatric areas, the native species has more diversified MHC IIB genes when compared to the invasive species, probably resulting from the different origins and dispersal of both species. A similar level of MHC polymorphism was found at population level in both species, suggesting similar mechanisms generating MHC diversity. In contrast, a higher number of DAB-like alleles per specimen were found in invasive species. Invasive species tended to express the alleles of two DAB lineages, whilst native species tended to express the alleles of only the DAB3 lineage. Hybrids have a pattern of MHC expression intermediate between both species. Whilst positive selection acting on peptide binding sites (PBS) was demonstrated in both species, a slightly higher number of positively selected sites were identified in C. nasus, which could result from parasite-mediated selection. Bayesian clustering analysis revealed a similar pattern of structuring for the genetic variation when using microsatellites or the MHC approach. We confirmed the importance of trans-species evolution for MHC polymorphism. In addition, we demonstrated bidirectional gene flow for MHC IIB genes in sympatric areas. The positive significant correlation between MHC and microsatellites suggests that demographic factors may contribute to MHC variation on a short time scale.

Molecular cloning and expression analysis of major histocompatibility complex class IIB gene of the Whitespotted bambooshark (Chiloscyllium plagiosum)

Major histocompatibility complex (MHC) plays an important role in the immune response to antigenic peptides in vertebrates. In this study, the full length of MHC IIB cDNA was isolated from the Whitespotted bambooshark (Chiloscyllium plagiosum) by homology cloning, and the rapid amplification of cDNA ends polymerase chain reaction. As a result, the MHC IIB cDNA is 1,407 bp, which contains an open reading frame (ORF) of 831 bp encoding a protein of 276 amino acids. Furthermore, seven alleles of the complete MHC IIB ORF were detected and the variable sites were mainly located in the immunoglobulin-like (β2) region. Tissue distribution analysis showed that MHC IIB can be detected in all the ten tissues examined, with the highest expression in the spleen and gill. Challenge of C. plagiosum with the pathogenic bacteria, Vibrio harveyi, resulted in significant changes in the expression of MHC IIB mRNA in the three immune-related tissues (gill, liver and spleen). These results show that the MHC IIB plays an important role in response to bacterial infection in elasmobranches.

Domestication does not narrow MHC diversity in Sus scrofa

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.

Female rose bitterling prefer MHC-dissimilar males: experimental evidence

The role of genetic benefits in female mate choice remains a controversial aspect of sexual selection theory. In contrast to "good allele" models of sexual selection, "compatible allele" models of mate choice predict that females prefer mates with alleles complementary to their own rather than conferring additive effects. While correlative results suggest complementary genetic effects to be plausible, direct experimental evidence is scarce. A previous study on the Chinese rose bitterling (Rhodeus ocellatus) demonstrated a positive correlation between female mate choice, offspring growth and survival, and the functional dissimilarity between the Major Histocompatibility Complex (MHC) alleles of males and females. Here we directly tested whether females used cues associated with MHC genes to select genetically compatible males in an experimental framework. By sequentially pairing females with MHC similar and dissimilar males, based on a priori known MHC profiles, we showed that females discriminated between similar and dissimilar males and deposited significantly more eggs with MHC dissimilar males. Notably, the degree of dissimilarity was an important factor for female decision to mate, possibly indicating a potential threshold value of dissimilarity for decision making, or of an indirect effect of the MHC.

MHC polymorphism and disease resistance to Vibrio anguillarum in 8 families of half-smooth tongue sole (Cynoglossus semilaevis)

Background

Genes in the major histocompatibility complex (MHC) have a critical role in both the innate and adaptive immune responses because of their involvement in presenting foreign peptides to T cells. However, the nature has remained largely unknown.

Results

We examined the genetic variation in MHC class IIB in half-smooth tongue sole (Cynoglossus semilaevis) after challenge with vibrio anguillarum. Two thousand and four hundred fry from 12 half-smooth tongue sole families were challenged with Vibrio anguillarum. To determine any association between alleles and resistance or susceptibility to V. anguillarum, 160 individuals from four high-resistance (HR, < 40.55% mortality) families and four low-resistance (LR, > 73.27% mortality) families were selected for MHC IIB exon2 gene sequence analysis. The MHC IIB exon2 genes of tongue sole displayed a high level of polymorphism and were discovered at least four loci. Meanwhile, the d_N/d_S [the ratio of non-synonymous (d_N) substitutions to synonymous (d_S) substitutions] in the peptide-binding region (PBR) was higher than that in the non-peptide-binding region (non-PBR). Eighty-eight alleles were discovered among 160 individuals, and 13 out of 88 alleles were used to analyze the distribution pattern between the resistant and susceptible families. Certain alleles presented in HR and LR with a different frequency, while other alleles were discovered in only the HR or LR families, not both. Five alleles, Cyse-DBB*6501, Cyse-DBB*4002, Cyse-DBB*6102, Cyse-DBB*5601 and Cyse-DBB*2801, were found to be associated with susceptibility to V. anguillarum with a frequency of 1.25%, 1.25%, 1.25%, 1.25% and 2.5% in the HR families, and 35%, 33.75%, 27.5%, 16.25%, 15% in the LR families (p < 0.01, 0.01, 0.01, 0.01, 0.01), respectively. Four alleles, Cyse-DBB*3301, Cyse-DBB*4701, Cyse-DBB*6801 and Cyse-DBB*5901, were found to be associated with resistance to V. anguillarum, with a frequency of 13.75%, 11.25%, 11.25%, 8.75% in the HR families and 1.25%, 1.25%, 1.25%, 1.25% and 1.25% in the LR families (p < 0.01, 0.05, 0.05 and p = 0.064), respectively.

Conclusions

Elucidation of the role of MHC II B genes in half-smooth tongue sole should prove to be helpful to the in-depth development of marker-assisted selective breeding in half-smooth tongue sole.

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

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

Adaptive molecular evolution of the Major Histocompatibility Complex genes, DRA and DQA, in the genus Equus

BACKGROUND

Major Histocompatibility Complex (MHC) genes are central to vertebrate immune response and are believed to be under balancing selection by pathogens. This hypothesis has been supported by observations of extremely high polymorphism, elevated nonsynonymous to synonymous base pair substitution rates and trans-species polymorphisms at these loci. In equids, the organization and variability of this gene family has been described, however the full extent of diversity and selection is unknown. As selection is not expected to act uniformly on a functional gene, maximum likelihood codon-based models of selection that allow heterogeneity in selection across codon positions can be valuable for examining MHC gene evolution and the molecular basis for species adaptations.

RESULTS

We investigated the evolution of two class II MHC genes of the Equine Lymphocyte Antigen (ELA), DRA and DQA, in the genus Equus with the addition of novel alleles identified in plains zebra (E. quagga, formerly E. burchelli). We found that both genes exhibited a high degree of polymorphism and inter-specific sharing of allele lineages. To our knowledge, DRA allelic diversity was discovered to be higher than has ever been observed in vertebrates. Evidence was also found to support a duplication of the DQA locus. Selection analyses, evaluated in terms of relative rates of nonsynonymous to synonymous mutations (dN/dS) averaged over the gene region, indicated that the majority of codon sites were conserved and under purifying selection (dN <dS). However, the most likely evolutionary codon models allowed for variable rates of selection across codon sites at both loci and, at the DQA, supported the hypothesis of positive selection acting on specific sites.

CONCLUSIONS

Observations of elevated genetic diversity and trans-species polymorphisms supported the conclusion that balancing selection may be acting on these loci. Furthermore, at the DQA, positive selection was occurring at antigen binding sites, suggesting that a few selected residues may play a significant role in equid immune function. Future studies in natural equid populations will be valuable for understanding the functional significance of the uniquely diverse DRA locus and for elucidating the mechanism maintaining diversity at these MHC loci.

MHC class IIalpha gene polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus)

Association between polymorphism of the major histocompatibility complex (MHC) genes and disease resistance has been documented for few teleosts. In this study, we first investigate the genetic variation at the MHC IIalpha gene in Japanese flounder (Paralichthys olivaceus) for survival after challenge with bacterial infection. To explore the specific allele associated with disease resistance, about 6000 individuals from 60 families challenged with Vibrio anguillarum, which causes significantly different mortality in flounder families. 15-20 individuals from each of six high-resistance (HR) and six low-resistance (LR) families were screened for their MHC class IIalpha genotypes using sequence analysis. High polymorphism of MHC IIalpha gene and at least two loci were discovered in Japanese flounder and the rate of d(N) occurred at a significantly higher frequency than that of d(S) in PBR and non-PBR, the balancing selection hypothesis could explain the high polymorphism of MHC IIalpha gene in flounder. From the 212 individuals, a total of 55 MHC class IIalpha alleles were identified, and 9 alleles were used to study association between alleles and resistance to disease. Three alleles, Paol-DAA*1301, Paol-DAA*1401 and Paol-DAA*2201 were significantly associated with resistance against V. anguillarum, and Paol-DAA*1001 and Paol-DAA*1501 alleles were significantly associated with increased susceptibility to V. anguillarum. This study confirmed the association between alleles of MHC class IIalpha gene and disease resistance or susceptibility to bacterial infection in flounder, and the disease resistance-related MHC markers could be used for molecular marker-assisted selective breeding in the flounder.

Evolution of an MHC class Ia gene fragment in four North American Morone species

A nucleotide sequence analysis of a fragment of a Morone MHC class Ia gene detected high levels of polymorphism in striped bass Morone saxatilis, white perch Morone americana and yellow bass Morone mississippiensis. Extremely low levels of MHC diversity, however, were detected in white bass Morone chrysops, suggesting the possibility of a severe population bottleneck for this species.

Major histocompatibility complex variation at class II DQA locus in the brown hare (Lepus europaeus)

The major histocompatability complex (MHC) is a multigene family of receptors that bind and present antigenic peptides to T-cells. Genes of the MHC are characterized by an outstanding genetic polymorphism, which is considered to be maintained by positive selection. Sites involved in peptide binding form binding pockets (P) that are collectively termed the peptide-binding region (PBR). In this study, we examined the level of MHC genetic diversity within and among natural populations of brown hare (Lepus europaeus) from Europe and Anatolia choosing for analysis of the second exon of the DQA locus, one of the most polymorphic class II loci. We aimed at an integrated population genetic analysis of L. europeaus by (i) correlating MHC polymorphism to genetic variability and phylogenetic status estimated previously from maternally (mtDNA) and biparentally (allozymes, microsatellites) inherited loci; and (ii) comparing full-length exon amino acid polymorphism with functional polymorphism in the PBR and the binding pockets P1, P6 and P9. A substantial level of DQA exon 2 polymorphism was detected with two completely different set of alleles between the Anatolian and European populations. However, the phylogeny of full-length exon 2 Leeu-DQA alleles did not show a strong phylogeographic signal. The presence of balancing selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous in the PBR and a trans-species pattern of evolution detected after phylogenetic reconstruction. The differentiating patterns detected between genetic and functional polymorphism, i.e. the number and the distribution of pocket variants within and among populations, indicated a hierarchical action of selection pressures.

Polymorphism and selection in the major histocompatibility complex DRA and DQA genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.

Resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 is influenced by major histocompatibility (MH) class II B gene polymorphism

The role of MH class II B (Cyca-DAB1-like) genes in resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 (CyHV-3), also known as koi herpesvirus (KHV) was analysed. The material consisted of 934 fish from six carp crosses. Fish were challenged with CyHV-3 at an age of 7 and 10 months. During challenge experiments the peak of mortality caused by CyHV-3 was observed at days 8-12 p.i. and the overall cumulative mortality reached 79.9%. Among six Cyca-DAB1-like genotypes, revealed by PCR-RF-SSCP analysis, one genotype (E) was found associated with higher resistance to CyHV-3. Three other genotypes (B, H and J) could be linked to higher susceptibility to CyHV-3. Analysis of the alleles that compose the Cyca-DAB1-like genotypes linked one particular allele (Cyca-DAB1*05) to significantly increased, and two alleles (Cyca-DAB1*02 and Cyca-DAB1*06) to significantly decreased resistance to CyHV-3. Our data indicate that MH class II B genes could be used as potential genetic markers in breeding of common carp for resistance to this virus.

Parental relatedness and major histocompatibility effects on early embryo survivorship in Atlantic salmon

Salmon have provided key insights into the relative influence of natural and sexual selection on major histocompatibility complex (MHC) variation. Natural selection on salmon MHC genes has been demonstrated in pathogen studies, and there is evidence of MHC-based mate choice (sexual selection). We tested whether parental MHC genes affect survivorship of juvenile Atlantic salmon (Salmo salar) by quantifying the influence of parental genome-wide relatedness and MHC genotype on survivorship to the swim-up stage. Thirteen microsatellite loci were used to estimate the influence of genome-wide relatedness between parents on offspring survivorship and MHC genotypes were determined by sequencing part of the class IIbeta gene. Our results revealed no significant relationship between early offspring survivorship and genome-wide relatedness, predicted MHC heterozygosity, or MHC allelic similarity. Overall, our data are consistent with the contention that excess MHC heterozygosity in Atlantic salmon juveniles is due to sexual selection as well as differential survival of offspring due to MHC genotype.