Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

Characterization of antisera raised against stickleback (Gasterosteus aculeatus) MHC class I and class II molecules

The three-spined stickleback (Gasterosteus aculeatus) is an important model organism for investigations on the maintenance of polymorphism of the major histocompatibility complex (MHC) of vertebrates. Analysis of functional aspects of MHC diversity in stickleback would benefit from the availability of MHC specific reagents. Here we characterize antisera raised against recombinant fusion proteins of stickleback MHC class I alpha and class II alpha and beta. Western blot analysis using recombinant proteins confirmed the specificity of the antisera. In brain and muscle preparations, neither of the MHC types was detectable. High levels of each MHC receptor type were observed in gills and spleen and lower levels in head kidneys. In histological sections of gills, epithelial cells of primary and secondary lamellae stained positive with MHC class I antiserum, while single, scattered cells stained positive for MHC class II. In sections of spleen and head kidney, considerable numbers of cells positive for either MHC type were detected. Molecular weight shift in SDS-PAGE after deglycosylation of MHC class I alpha and class II beta confirmed the predicted glyco-protein character of the molecules. The majority of MHC II alpha was not glycosylated; only a small fraction of MHC II alpha was susceptible to deglycosylation. This suggests differential expression of the two stickleback MHC II alpha genes (Gaac-DAA, Gaac-DBA) only one of which (Gaac-DBA) has a site for N-linked glycosylation.     

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.     

Diversity and evolutionary patterns of immune genes in free-ranging Namibian leopards (Panthera pardus pardus)

The genes of the major histocompatibility complex (MHC) are a key component of the mammalian immune system and have become important molecular markers for fitness-related genetic variation in wildlife populations. Currently, no information about the MHC sequence variation and constitution in African leopards exists. In this study, we isolated and characterized genetic variation at the adaptively most important region of MHC class I and MHC class II-DRB genes in 25 free-ranging African leopards from Namibia and investigated the mechanisms that generate and maintain MHC polymorphism in the species. Using single-stranded conformation polymorphism analysis and direct sequencing, we detected 6 MHC class I and 6 MHC class II-DRB sequences, which likely correspond to at least 3 MHC class I and 3 MHC class II-DRB loci. Amino acid sequence variation in both MHC classes was higher or similar in comparison to other reported felids. We found signatures of positive selection shaping the diversity of MHC class I and MHC class II-DRB loci during the evolutionary history of the species. A comparison of MHC class I and MHC class II-DRB sequences of the leopard to those of other felids revealed a trans-species mode of evolution. In addition, the evolutionary relationships of MHC class II-DRB sequences between African and Asian leopard subspecies are discussed.     

Simple approach to reduce PCR artefact formation leads to reliable genotyping of MHC and other highly polymorphic loci--implications for evolutionary analysis

Genetic variation in coding regions is of strong interest for biologists as it represents an important factor that drives evolution. To analyse polymorphic loci, researchers usually rely on commonly used typing techniques such as cloning, SSCP, DGGE or RSCA. However, there are potential pitfalls in screening multi-allelic templates, which are mainly the formation of sequence chimeras during PCR amplification, and mosaic sequences during cloning. One of the most challenging genomic regions to explore is the Major Histocompatibility Complex (MHC), which codes for peptide-binding proteins of the vertebrate's adaptive immune system and is well known for its exceptional polymorphism. We compared the effect of two different PCR amplification approaches in a study of the MHC class IIB genes of the three-spined stickleback (Gasterosteus aculeatus). One approach used standard PCR conditions and the other a combination of several measures to eliminate PCR artefacts. In both approaches, the amplicons obtained were cloned and sequenced. In the first, established approach, 24% of the clones represented artefacts, while in the second approach the number of artefacts were reduced ten-fold. Furthermore, it enabled easy differentiation between real alleles and artificial sequences. We also analysed the potential effects of such artefacts in genetic analysis and evolutionary interpretation, and found a slight reduction in the signature of positive selection and an increase in recombination events. Consequently, we strongly recommend to apply the new PCR approach described in this study when genotyping MHC or other polymorphic genes.     

Characterization of MHC class I and II genes in a subantarctic seabird,the blue petrel, <Emphasis Type="Italic">Halobaena caerulea</Emphasis> (Procellariiformes)

The great polymorphism observed in the major histocompatibility complex (MHC) genes is thought to be maintained by pathogen-mediated selection possibly combined with MHC-disassortative mating, guided by MHC-determined olfactory cues. Here, we partly characterize the MHC class I and II B of the blue petrel, Halobaena caerulea (Procellariiformes), a bird with significant olfactory abilities that lives under presumably low pathogen burdens in Subantarctica. Blue petrels are long-lived, monogamous birds which suggest the necessity of an accurate mate choice process. The species is ancestral to songbirds (Passeriformes; many MHC loci), although not to gamefowls (Galliformes; few MHC loci). Considering the phylogenetic relationships and the low subantarctic pathogen burden, we expected few rather than many MHC loci in the blue petrel. However, when we analysed partial MHC class I and class II B cDNA and gDNA sequences we found evidence for as many as at least eight MHC class I loci and at least two class II B loci. These class I and II B sequences showed classical MHC characteristics, e.g. high nucleotide diversity, especially in putative peptide-binding regions where signatures of positive selection was detected. Trans-species polymorphism was found between MHC class II B sequences of the blue petrel and those of thin-billed prion, Pachyptila belcheri, two species that diverged ∼25 MYA. The observed MHC allele richness in the blue petrel may well serve as a basis for mate choice, especially since olfactory discrimination of MHC types may be possible in this species.     

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 variation in MHC class II expression and interactions with MHC sequence polymorphism in three-spined sticklebacks

Genes of the major histocompatibility complex (MHC) have been studied for several decades because of their pronounced allelic polymorphism. Structural allelic polymorphism is, however, not the only source of variability subjected to natural selection. Genetic variation may also exist in gene expression patterns. Here, we show that in a natural population of three-spined sticklebacks (Gasterosteus aculeatus) the expression of MHC class IIB genes was positively correlated with parasite load, which indicates increased immune activation of the MHC when infections are frequent. To experimentally study MHC expression, we used laboratory-bred sticklebacks that were exposed to three naturally occurring species of parasite. We found strong differences in MHC class IIB expression patterns among fish families, which were consistent over two generations, thus demonstrating a genetic component. The average number of MHC class IIB sequence variants within families was negatively correlated to the MHC expression level suggesting compensatory up-regulation in fish with a low (i.e. suboptimal) MHC sequence variability. The observed differences among families and the negative correlation with individual sequence diversity imply that MHC expression is evolutionary relevant for the onset and control of the immune response in natural populations.     

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

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

Parasite diversity,patterns of MHC II variation and olfactory based mate choice in diverging three-spined stickleback ecotypes

Ecological speciation has been the subject of intense research in evolutionary biology but the genetic basis of the actual mechanism driving reproductive isolation has rarely been identified. The extreme polymorphism of the major histocompatibility complex (MHC), probably maintained by parasite-mediated selection, has been proposed as a potential driver of population divergence. We performed an integrative field and experimental study using three-spined stickleback river and lake ecotypes. We characterized their parasite load and variation at MHC class II loci. Fish from lakes and rivers harbor contrasting parasite communities and populations possess different MHC allele pools that could be the result of a combined action of genetic drift and parasite-mediated selection. We show that individual MHC class II diversity varies among populations and is lower in river ecotypes. Our results suggest the action of homogenizing selection within habitat type and diverging selection between habitat types. Finally, reproductive isolation was suggested by experimental evidence: in a flow channel design females preferred assortatively the odor of their sympatric male. This demonstrates the role of olfactory cues in maintaining reproductive isolation between diverging fish ecotypes.     

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.     

Historical and contemporary selection of teleost MHC genes: did we leave the past behind?

The extreme polymorphism of antigen‐presenting genes of the major histocompatibility complex (MHC) has spurred intense research unparalleled for any other gene family. This applies also to teleosts where sequence information is available for 3559 MHC class I and class II allelic variants from 137 species. This review summarizes current knowledge on the origin and maintenance of diversity at classical MHC loci. Most studies identified positive selection (i.e. elevated rates of non‐synonymous over synonymous substitutions, dN/dS) as a sign of balancing selection. A meta‐analysis on nine species with sufficient numbers of class I and class II sequences revealed that recombination rate and intensity of positive selection were positively correlated, suggesting that recombination and gene conversion played a significant role in shaping the allelic repertoire. Processes that create diversity over long timescales need to be complemented by contemporary balancing selection, either through overdominance or frequency‐dependent selection, in order to explain the high allelic diversity observed today. While some evidence for overdominance exists for a few taxa (mainly salmonids) by correlating parasite infection data or survival to MHC genotypes, field or experimental data on negative frequency‐dependent selection are lacking altogether, even though some fish species are particularly suitable as model systems. Theoretical predictions suggest that negative frequency‐dependent selection is necessary to maintain the existing polymorphism. Hence, future empirical studies should focus on detecting signals that differentiate between mechanisms of contemporary selection rather than repeatedly showing historical selection events.     

Multiple parasites are driving major histocompatibility complex polymorphism in the wild

Parasite mediated selection may result in arms races between host defence and parasite virulence. In particular, simultaneous infections from multiple parasite species should cause diversification (i.e. balancing selection) in resistance genes both at the population and the individual level. Here, we tested these ideas in highly polymorphic major histocompatibility complex (MHC) genes from three-spined sticklebacks (Gasterosteus aculeatus L.). In eight natural populations, parasite diversity (15 different species), and MHC class IIB diversity varied strongly between habitat types (lakes vs. rivers vs. estuaries) with lowest values in rivers. Partial correlation analysis revealed an influence of parasite diversity on MHC class IIB variation whereas general genetic diversity assessed at seven microsatellite loci was not significantly correlated with parasite diversity. Within individual fish, intermediate, rather than maximal allele numbers were associated with minimal parasite load, supporting theoretical models of self-reactive T-cell elimination. The optimal individual diversity matched those values female fish try to achieve in their offspring by mate choice. We thus present correlative evidence supporting the 'allele counting' strategy for optimizing the immunocompetence in stickleback offspring.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

Molecular characterization,balancing selection,and genomic organization of the tree shrew (Tupaia belangeri) MHC class I gene

The major histocompatibility complex (MHC) class I genes play a pivotal role in the adaptive immune response among vertebrates. Accordingly, in numerous mammals the genomic structure and molecular characterization of MHC class I genes have been thoroughly investigated. To date, however, little is known about these genes in tree shrews, despite the increasingly popularity of its usage as an animal model. To address this deficiency, we analyzed the structure and characteristic of the tree shrew MHC class I genes (Tube-MHC I) and performed a comparative gene analysis of the tree shrew and other mammal species. We found that the full-length cDNA sequence of the tree shrew MHC class I is 1074 bp in length. The deduced peptide is composed of 357 amino acids containing a leader peptide, an α1 and α2 domain, an α3 domain, a transmembrane domain and a cytoplasmic domain. Among these peptides, the cysteines, CD8⁺ interaction and N-glycosylation sites are all well conserved. Furthermore, the genomic sequence of the tree shrew MHC class I gene was identified to be 3180 bp in length, containing 8 exons and 7 introns. In 21 MHC class I sequences, we conducted an extensive study of nucleotide substitutions. The results indicated that in the peptide binding region (PBR) the rate of non-synonymous substitutions (dN) to synonymous substitutions (dS) was greater than 1, suggesting balancing selection at the PBR. These findings provide valuable contributions in furthering our understanding of the structure, molecular polymorphism, and function of the MHC class I genes in tree shrews, further improving their utility as an animal model in biomedical research.     

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.