Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Evolutionary trade-offs constraining the MHC gene expansion: beyond simple TCR depletion model

The immune system is as much shaped by the pressure of pathogens as it is by evolutionary trade-offs that constrain its structure and function. A perfect example comes from the major histocompatibility complex (MHC), molecules that initiate adaptive immune response by presentation of foreign antigens to T cells. The remarkable, population-level polymorphism of MHC genes is assumed to result mainly from a co-evolutionary arms race between hosts and pathogens, while the limited, within-individual number of functional MHC loci is thought to be the consequence of an evolutionary trade-off between enhanced pathogen recognition and excessive T cell depletion during negative selection in the thymus. Certain mathematical models and infection studies suggest that an intermediate individual MHC diversity would thus be optimal. A recent, more direct test of this hypothesis has shown that the effects of MHC diversity on T-cell receptor (TCR) repertoires may differ between MHC classes, supporting the depletion model only for MHC class I. Here, we used the bank vole (Myodes=Cletronomys glareolus), a rodent species with variable numbers of expressed MHC genes, to test how an individual MHC diversity influences the proportions and TCR repertoires of responding T cell subsets. We found a non-linear relationship between MHC diversity and T cell proportions (with intermediate MHC numbers coinciding with the largest T cell proportions), perhaps reflecting an optimality effect of balanced positive and negative thymic selection. The association was strongest for the relationship between MHC class I and splenic CD8+ T cells. The CD8+ TCR richness alone was unaffected by MHC class I diversity, suggesting that MHC class I expansion may be limited by decreasing T cell counts, rather than by direct depletion of TCR richness. In contrast, CD4+ TCR richness was positively correlated with MHC class II diversity, arguing against a universal TCR depletion. It also suggests that different evolutionary forces or trade-offs may limit the within-individual expansion of the MHC class II loci.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Comparative analysis of host resistance to Sarcoptes scabiei var. cuniculi in two different rabbit breeds

BACKGROUND

Scabies, caused by infestation of the mite Sarcoptes scabiei, is one of the most severe ectoparasitic diseases in rabbits. Scabies seriously affects the commercial rabbit breeding, causing severe economic losses. Host resistance to S. scabiei is an important factor in further development of the rabbit industry. In the present study, we compared the host resistance to S. scabiei var. cuniculi of a new breed of domestic rabbit propagated by the Sichuan Animal Sciences Academy (QiXing rabbit, QX) compared with that of a traditional rabbit breed in the domestic rabbit industry (IRA rabbit, IRA).

METHODS

Both QX and IRA rabbits were experimentally infested with live S. scabiei var. cuniculi mites for 48 h. Then, during the course of four-week experimental infestation period, the body weight of rabbits was recorded every two weeks for calculating body-weight variations in comparison to the non-infested control rabbits. Skin lesions in the foot area were assessed on weekly basis and serum samples were tested weekly for the estimation of changes in the total antibody levels (IgG, IgE and IgM). Moreover, DNA extracted from the blood samples was amplified for analysis of the genetic diversity in the major histocompatibility complex, class II, DQ Alpha (MHC-DQA) gene.

RESULTS

Compared to the IRA rabbits, the QX rabbits showed a significantly higher (P < 0.05) relative body weight gain compared to the non-infested control rabbits and significantly lower (P < 0.05) scores for foot skin lesions and higher levels of IgG, IgE and IgM at weeks 1 to 4, week 2 and week 1 post-infestation, respectively. Furthermore, a polymorphism site at position 103 bp of exon two of MHC-DQA gene and a different gene frequency were found between two rabbit breeds, suggesting the genetic basis for the differential host resistance to the S. scabiei var. cuniculi between two rabbit breeds.

CONCLUSIONS

The QX rabbits showed higher host resistance to S. scabiei var. cuniculi compared to the IRA rabbits at the clinical, immunological and genetic levels. These results provide a reference for the breeding of rabbits with adequately improved and sustained host resistance to scabies in the domestic rabbit industry.

Molecular Cloning and Bioinformatics Analysis of DQA Gene from Mink (Neovison vison)

In the present study, we cloned, sequenced, and explored the structural and functional characteristics of the major histocompatibility complex (MHC)-DQA gene from mink (Neovison vison) for the first time. The full-length sequence of DQA gene was 1147-bp-long, contained a coding region of 768-bp, which was predicted to encoding 255 amino acid residues. The comparison between DQA from mink (Neovison vison) and other MHC-DQA molecules from different animal species showed that nucleotide and encoded amino acid sequences of the mink DQA gene exhibited high similarity with the ferret (Mustela pulourius furo). Phylogenetic analysis revealed that mink (Neovison vison) DQA is grouped with that of ferret (Mustela pulourius furo). The cloned sequence contained a 23-amino acid NH₂-terminal signal sequence with the signal peptide cutting site located in amino acids 23–24, and had three Asn-Xaa-Ser/Thr sequons. Three cysteine residues were also identified (Cys-85, Cys-121, and Cys-138). The 218 to 240 amino acids were predicted to be the transmembrane domains. The prediction of the secondary structure revealed three α-helixes and fourteen β-sheets in Neovison vison DQA protein, while random coil was a major pattern. In this study, the whole CDS sequence of Neovison vison DQA gene was successfully cloned, which was valuable for exploring the function and antiviral molecular mechanisms underlying the molecule. The findings of the present study have laid the foundation for the disease resistance and breeding of mink.

Major histocompatibility complex class I diversity limits the repertoire of T cell receptors

Major histocompatibility complex (MHC) genes encode proteins that initiate adaptive immune responses through the presentation of foreign antigens to T cells. The high polymorphism found at these genes, thought to be promoted and maintained by pathogen-mediated selection, contrasts with the limited number of MHC loci found in most vertebrates. Although expressing many diverse MHC genes should broaden the range of detectable pathogens, it has been hypothesized to also cause deletion of larger fractions of self-reactive T cells, leading to a detrimental reduction of the T cell receptor (TCR) repertoire. However, a key prediction of this TCR depletion hypothesis, that the TCR repertoire should be inversely related to the individual MHC diversity, has never been tested. Here, using high-throughput sequencing and advanced sequencing error correction, we provide evidence of such an association in a rodent species with high interindividual variation in the number of expressed MHC molecules, the bank vole (Myodes glareolus). Higher individual diversity of MHC class I, but not class II, was associated with smaller TCR repertoires. Our results thus provide partial support for the TCR depletion model, while also highlighting the complex, potentially MHC class-specific mechanisms by which autoreactivity may trade off against evolutionary expansion of the MHC gene family.

Profiling of the TCRβ repertoire in non-model species using high-throughput sequencing

In recent years, immune repertoire profiling with high-throughput sequencing (HTS) has advanced our understanding of adaptive immunity. However, fast progress in the field applied mostly to human and mouse research, with only few studies devoted to other model vertebrates. We present the first in-depth characterization of the T-cell receptor (TCR) repertoire in a non-model mammal (bank vole, Myodes glareolus), widely used in ecological and evolutionary research. We used RNA from spleens, 5′RACE and HTS to describe V and J segments of TCRβ, qualitatively characterize preferential V–J segment usage and CDR3 length distribution. Overall orthology to murine genes was preserved, with 11 J and 37 V genes found in voles (although 3 V genes lacked a close orthologue). Further, we implemented unique molecular identifiers for quantitative analysis of CDR3 repertoire with stringent error correction. A conservative, lower bound estimation of the TCRβ repertoire was similar to that found for mice (1.7–2.3 × 10⁵ clonotypes). We hope that by providing an easy-to-follow molecular protocol and on-line bioinformatics tools that do not require reference sequences (AmpliTCR and AmpliCDR3), we will encourage HTS immune repertoire profiling in other non-model vertebrates, thus opening new research avenues in e.g. comparative immunology, ecology and evolutionary biology.

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.

De novo transcriptome assembly facilitates characterisation of fast-evolving gene families, MHC class I in the bank vole (Myodes glareolus)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2-4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.

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.

Evidence for selection maintaining MHC diversity in a rodent species despite strong density fluctuations

Strong spatiotemporal variation in population size often leads to reduced genetic diversity limiting the adaptive potential of individual populations. Key genes of adaptive variation are encoded by the immune genes of the major histocompatibility complex (MHC) playing an essential role in parasite resistance. How MHC variation persists in rodent populations that regularly experience population bottlenecks remains an important topic in evolutionary genetics. We analysed the consequences of strong population fluctuations on MHC class II DRB exon 2 diversity in two distant common vole (Microtus arvalis) populations in three consecutive years using a high-throughput sequencing approach. In 143 individuals, we detected 25 nucleotide alleles translating into 14 unique amino acid MHC alleles belonging to at least three loci. Thus, the overall allelic diversity and amino acid distance among the remaining MHC alleles, used as a surrogate for the range of pathogenic antigens that can be presented to T-cells, are still remarkably high. Both study populations did not show significant population differentiation between years, but significant differences were found between sites. We concluded that selection processes seem to be strong enough to maintain moderate levels of MHC diversity in our study populations outcompeting genetic drift, as the same MHC alleles were conserved between years. Differences in allele frequencies between populations might be the outcome of different local parasite pressures and/or genetic drift. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in cyclic populations.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

Major histocompatibility complex class II polymorphisms in forest musk deer (Moschus berezovskii) and their probable association with purulent disease

Genes of the major histocompatibility complex (MHC) family are crucial in immune responses because they present pathogenic peptides to T cells. In this study, we analysed the genetic variation in forest musk deer (Moschus berezovskii) MHC II genes and its potential association with musk deer purulent disease. In total, 53 purulent disease-susceptible and 46 purulent disease-resistant individuals were selected for MHC II exon 2 fragment analysis. Among them, 16 DQ alleles and four additional DR alleles were identified, with DQ exon 2 fragments displaying a low level of polymorphism. The nonsynonymous substitutions exceeded the synonymous substitutions in the peptide-binding sites of DQA2, DQB1 and DQB2. Then, 28 MHC II alleles were used to analyse the distribution patterns of purulent disease between the susceptible and resistant groups. Among them, three alleles (DQA1*01, DQA1*02 and DQA2*04) were found to be resistant, and five alleles (DRB3*07, DQA1*03, DQA1*04, DQA2*05 and DQA2*06) were found to increase susceptibility. Additionally, three haplotypes were found to be putatively associated with musk deer purulent disease. However, these three haplotypes were only found in the resistant or susceptible group, and their frequencies were low. The results from our study support a contributory role of MHC II polymorphisms in the development of purulent disease in forest musk deer.

Positive selection on MHC class II DRB and DQB genes in the bank vole (Myodes glareolus)

The major histocompatibility complex (MHC) class IIB genes show considerable sequence similarity between loci. The MHC class II DQB and DRB genes are known to exhibit a high level of polymorphism, most likely maintained by parasite-mediated selection. Studies of the MHC in wild rodents have focused on DRB, whilst DQB has been given much less attention. Here, we characterised DQB genes in Swedish bank voles Myodes glareolus, using full-length transcripts. We then designed primers that specifically amplify exon 2 from DRB (202 bp) and DQB (205 bp) and investigated molecular signatures of natural selection on DRB and DQB alleles. The presence of two separate gene clusters was confirmed using BLASTN and phylogenetic analysis, where our seven transcripts clustered according to either DQB or DRB homologues. These gene clusters were again confirmed on exon 2 data from 454-amplicon sequencing. Our DRB primers amplify a similar number of alleles per individual as previously published DRB primers, though our reads are longer. Traditional d N/d S analyses of DRB sequences in the bank vole have not found a conclusive signal of positive selection. Using a more advanced substitution model (the Kumar method) we found positive selection in the peptide binding region (PBR) of both DRB and DQB genes. Maximum likelihood models of codon substitutions detected positively selected sites located in the PBR of both DQB and DRB. Interestingly, these analyses detected at least twice as many positively selected sites in DQB than DRB, suggesting that DQB has been under stronger positive selection than DRB over evolutionary time.

Contrasted evolutionary histories of two Toll-like receptors (Tlr4 and Tlr7) in wild rodents (MURINAE)

BACKGROUND

In vertebrates, it has been repeatedly demonstrated that genes encoding proteins involved in pathogen-recognition by adaptive immunity (e.g. MHC) are subject to intensive diversifying selection. On the other hand, the role and the type of selection processes shaping the evolution of innate-immunity genes are currently far less clear. In this study we analysed the natural variation and the evolutionary processes acting on two genes involved in the innate-immunity recognition of Microbe-Associated Molecular Patterns (MAMPs).

RESULTS

We sequenced genes encoding Toll-like receptor 4 (Tlr4) and 7 (Tlr7), two of the key bacterial- and viral-sensing receptors of innate immunity, across 23 species within the subfamily Murinae. Although we have shown that the phylogeny of both Tlr genes is largely congruent with the phylogeny of rodents based on a comparably sized non-immune sequence dataset, we also identified several potentially important discrepancies. The sequence analyses revealed that major parts of both Tlrs are evolving under strong purifying selection, likely due to functional constraints. Yet, also several signatures of positive selection have been found in both genes, with more intense signal in the bacterial-sensing Tlr4 than in the viral-sensing Tlr7. 92% and 100% of sites evolving under positive selection in Tlr4 and Tlr7, respectively, were located in the extracellular domain. Directly in the Ligand-Binding Region (LBR) of TLR4 we identified two rapidly evolving amino acid residues and one site under positive selection, all three likely involved in species-specific recognition of lipopolysaccharide of gram-negative bacteria. In contrast, all putative sites of LBRTLR7 involved in the detection of viral nucleic acids were highly conserved across rodents. Interspecific differences in the predicted 3D-structure of the LBR of both Tlrs were not related to phylogenetic history, while analyses of protein charges clearly discriminated Rattini and Murini clades.

CONCLUSIONS

In consequence of the constraints given by the receptor protein function purifying selection has been a dominant force in evolution of Tlrs. Nevertheless, our results show that episodic diversifying parasite-mediated selection has shaped the present species-specific variability in rodent Tlrs. The intensity of diversifying selection was higher in Tlr4 than in Tlr7, presumably due to structural properties of their ligands.

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high-amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model-averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.

Global major histocompatibility class II β (mh-IIβ)-polymorphism in Arctic charr Salvelinus alpinus

This study explored the use of the gene encoding the β subunit of the major histocompatibility (MH) receptor as a population marker in Arctic charr Salvelinus alpinus. The use of this polymorphic marker allowed differentiation of the S. alpinus lineages previously defined using mitochondrial DNA (mtDNA) but also allowed differentiation between the populations studied within those lineages. The majority of the variation observed here occurred prior to the last glaciation event. Nevertheless, all S. alpinus populations were differentiated using both MH Class II β (mh-IIβ) sequences and allelic frequencies. The fact that all the populations studied presented high rates of non-synonymous: synonymous substitutions and high levels of interpopulation variation, suggested mh-IIβ as an ideal marker to assess differentiation among S. alpinus populations in ways that may represent divergence both by genetic drift and natural adaptation to the local environment.

Insight into gene evolution within Cervidae and Bovidae through genetic variation in MHC-DQA in the black muntjac (Muntiacus crinifrons)

The critical role that the major histocompatibility complex plays in the immune recognition of parasites and pathogens makes its evolutionary dynamics exceptionally relevant to ecology, population biology, and conservation studies. The black muntjac is a rare deer endemic to a small mountainous region in eastern China. We found that this species has two DQA loci through cDNA expression and sequence variation analysis. The level of variation at both DQA loci was found to be extremely low (three alleles for DQA1 and four alleles for DQA2), possibly because of past bottlenecks and the species' relatively solitary behavior pattern. The ratio of d(N)/d(S) in the putative peptide binding region of the DQA2 locus (13.36, P = 0.012) was significantly larger than one but not that of DQA1 (0.94, P = 0.95), suggesting strong positive selection at the DQA2 but not at the DQA1 locus. This difference might reflect different sets of evolutionary selection pressures acting on the two loci. The phylogenetic tree showed that DQA1 alleles from two species of Cervidae and two of Bovidae grouped together, as did the DQA2 alleles. However, different genes from the four species were located in separate branches. These results lead us to suggest that these DQA alleles are derived from primordial DQA genes from a common ancestor and are maintained in Cervidae and Bovidae since their divergence around 25.5-27.8 million years ago.

Signatures of balancing selection are maintained at disease resistance loci following mating system evolution and a population bottleneck in the genus Capsella

Background

Population bottlenecks can lead to a loss of variation at disease resistance loci, which could have important consequences for the ability of populations to adapt to pathogen pressure. Alternatively, current or past balancing selection could maintain high diversity, creating a strong heterogeneity in the retention of polymorphism across the genome of bottlenecked populations. We sequenced part of the LRR region of 9 NBS-LRR disease resistance genes in the outcrossing Capsella grandiflora and the recently derived, bottlenecked selfing species Capsella rubella, and compared levels and patterns of nucleotide diversity and divergence with genome-wide reference loci.

Results

In strong contrast with reference loci, average diversity at resistance loci was comparable between C. rubella and C. grandiflora, primarily due to two loci with highly elevated diversity indicative of past or present balancing selection. Average between-species differentiation was also reduced at the set of R-genes compared with reference loci, which is consistent with the maintenance of ancestral polymorphism.

Conclusions

Historical or ongoing balancing selection on plant disease resistance genes is a likely contributor to the retention of ancestral polymorphism in some regions of the bottlenecked Capella rubella genome.

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

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

The functional gene diversity in natural populations over postglacial areas: the shaping mechanisms behind genetic composition of longnose dace (Rhinichthys cataractae) in northeastern North America

The diversity of functional genes and the related processes are important issues for conservation biology. This is especially relevant for populations that have suffered from demographic reduction as a consequence of the processes of postglacial colonization. In this perspective, the aims of the present study are (1) to quantify the genetic diversity of functional genes and (2) to disentangle the long- and short-term effects of natural selection that shapes genetic diversity from those of drift, mutation, and allopatric fragmentation. This research was conducted using an extensive genetic polymorphism analysis of populations of longnose dace (Rhinichthys cataractae) living over an area once covered by Pleistocene glaciations. The sequence and diversity of one exon of three genes (MHC IIβ, growth hormone, and trypsin) were jointly analyzed with non-coding nuclear loci from 27 populations; these populations were sampled over four major basins of northeastern North America. The survey revealed a surprisingly low allelic richness, especially for the MHC gene, considering the number of individuals and populations sampled. The results suggest that there is a complex mixture of different evolutionary processes shaping the level of polymorphism among longnose dace. While our study underlines the importance of the short-term effects of neutral processes and the major impact of post-glacial colonization on gene diversity, locally dependent balancing selection was detected on MHC. From this perspective, our results support an understanding of the importance of drift on functional gene diversity but also highlight the transient effects of natural selection on allelic composition, even in populations that show drastic reduction of genetic diversity.

Evolutionary genetics of MHC class II beta genes in the brown hare, Lepus europaeus

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na = 11, D_est = 0.071 and Na = 15, D_est = 0.409, respectively), but not for the DPB gene (Na = 4, D_est = 0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω > 1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids.

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.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

Genetic structure and contrasting selection pattern at two major histocompatibility complex genes in wild house mouse populations

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

Associations between MHC genes and Puumala virus infection in Myodes glareolus are detected in wild populations, but not from experimental infection data

We analysed the influence of MHC class II Dqa and Drb genes on Puumala virus (PUUV) infection in bank voles (Myodes glareolus). We considered voles sampled in five European localities or derived from a previous experiment that showed variable infection success of PUUV. The genetic variation observed in the Dqa and Drb genes was assessed by using single-strand conformation polymorphism and pyrosequencing methods, respectively. Patterns were compared with those obtained from 13 microsatellites. We revealed significant genetic differentiation between PUUV-seronegative and -seropositive bank voles sampled in wild populations, at the Drb gene only. The absence of genetic differentiation observed at neutral microsatellites confirmed the important role of selective pressures in shaping these Drb patterns. Also, we found no significant associations between infection success and MHC alleles among laboratory-colonized bank voles, which is explained by a loss of genetic variability that occurred during the captivity of these voles.

A 454 multiplex sequencing method for rapid and reliable genotyping of highly polymorphic genes in large-scale studies

Background

High-throughput sequencing technologies offer new perspectives for biomedical, agronomical and evolutionary research. Promising progresses now concern the application of these technologies to large-scale studies of genetic variation. Such studies require the genotyping of high numbers of samples. This is theoretically possible using 454 pyrosequencing, which generates billions of base pairs of sequence data. However several challenges arise: first in the attribution of each read produced to its original sample, and second, in bioinformatic analyses to distinguish true from artifactual sequence variation. This pilot study proposes a new application for the 454 GS FLX platform, allowing the individual genotyping of thousands of samples in one run. A probabilistic model has been developed to demonstrate the reliability of this method.

Results

DNA amplicons from 1,710 rodent samples were individually barcoded using a combination of tags located in forward and reverse primers. Amplicons consisted in 222 bp fragments corresponding to DRB exon 2, a highly polymorphic gene in mammals. A total of 221,789 reads were obtained, of which 153,349 were finally assigned to original samples. Rules based on a probabilistic model and a four-step procedure, were developed to validate sequences and provide a confidence level for each genotype. The method gave promising results, with the genotyping of DRB exon 2 sequences for 1,407 samples from 24 different rodent species and the sequencing of 392 variants in one half of a 454 run. Using replicates, we estimated that the reproducibility of genotyping reached 95%.

Conclusions

This new approach is a promising alternative to classical methods involving electrophoresis-based techniques for variant separation and cloning-sequencing for sequence determination. The 454 system is less costly and time consuming and may enhance the reliability of genotypes obtained when high numbers of samples are studied. It opens up new perspectives for the study of evolutionary and functional genetics of highly polymorphic genes like major histocompatibility complex genes in vertebrates or loci regulating self-compatibility in plants. Important applications in biomedical research will include the detection of individual variation in disease susceptibility. Similarly, agronomy will benefit from this approach, through the study of genes implicated in productivity or disease susceptibility traits.

Negative relationships between cellular immune response, Mhc class II heterozygosity and secondary sexual trait in the montane water vole

Heterogeneities in immune responsiveness may affect key epidemiological parameters and the dynamics of pathogens. The roles of immunogenetics in these variations remain poorly explored. We analysed the influence of Major histocompatibility complex (Mhc) genes and epigamic traits on the response to phytohaemagglutinin in males from cyclic populations of the montane water vole (Arvicola scherman). Besides, we tested the relevance of lateral scent glands as honest signals of male quality. Our results did not corroborate neither the hypotheses of genome-wide heterozygosity-fitness correlation nor the Mhc heterozygote advantage. We found a negative relationship between Mhc hetetozygosity and response to phytohaemagglutinin, mediated by a specific Mhc homozygous genotype. Our results therefore support the hypothesis of the Arte-Dqa-05 homozygous genotype being a ‘good’ Mhc variant in terms of immunogenetic quality. The development of the scent glands seems to be an honest signal for mate choice as it is negatively correlated with helminth load. The ‘good gene’ hypothesis was not validated as Arte-Dqa-05 homozygous males did not exhibit larger glands. Besides, the negative relationship observed between the size of these glands and the response to phytohaemagglutinin, mainly for Mhc homozygotes, corroborates the immunocompetence handicap hypothesis. The Mhc variants associated with larger glands remain yet to be determined.

High MHC DQB variation and asymmetric allelic distribution in the endangered Yangtze finless porpoise, Neophocaena phocaenoides asiaeorientalis

The endangered Yangtze finless porpoise is found in the middle and lower reaches of the Yangtze River and its adjoining big lakes. To explore the major histocompatibility complex (MHC) genetic diversity and allelic distribution patterns across its range, we investigated variation at DQB exon 2. From 76 porpoises, we identified 18 DQB sequences. The freshwater Yangtze populations had much higher allelic diversity than marine populations. Among these freshwater populations, the middle-reach population had higher allelic diversity than the lower-reach population. The high DQB diversity level, relative to that of a neutral mtDNA locus, suggests that balancing selection is acting at the DQB gene and that rapid evolution and local positive selection play critical roles in generating and retaining high MHC diversity in the freshwater population. As the balancing selection might be driven by environmental pathogens, we suggest that maintaining MHC variation should be a high priority in the conservation and management of this endangered population, especially as an ex situ conservation strategy.

Polymorphism and signatures of selection in the multimammate rat DQB gene

In order to test if DQB is a good candidate marker to investigate the relationship between major histocompatibility complex genes and pathogens in natural populations of Mastomys natalensis, we assessed the polymorphism and evolutionary history of this gene. Twenty-four individuals were genotyped for exon 2 of DQB using capillary electrophoresis single-strand conformation polymorphism, cloning, and sequencing. We found 21 different alleles. Four individuals show three alleles implying a duplication event in the history of this gene. Each distinct sequence translates to give a distinct amino acid sequence and there are strong signals of positive selection on peptide binding sites. Signals of recombination were found in the sequences suggesting that recombination has played a role in generating allelic diversity. Although trans-taxon polymorphism is present at the interspecific level in DQB exon 2 sequences of Mus species, we did not find any evidence of allele sharing among Muridae genera. This indicates a temporal limit of DQB allele sharing in Muridae of less than 8 Mya. In conclusion, although DQB seems to be a good marker to investigate pathogen-driven selection, the polymorphism of gene copy number may restrict its utility in natural populations.

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.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Positive selection on mammalian MHC-DQ genes revisited from a multispecies perspective

Major histocompatibility complex class II DQA and DQB genes have been shown to be under positive selection in certain mammalian species but not in others, fuelling a debate about how their polymorphism has evolved. In this study, we have analysed whether polymorphism in the peptide-binding region (PBR) of DQA (190 sequences, 11 species) and DQB (209 sequences, 7 species) molecules is positively selected by using both approximate (Nei-Gojobori, Li-Wu-Luo and Pamilo-Bianchi-Li) and maximum-likelihood methods. The results obtained with approximate methods were rather inconsistent for DQA, probably due to the high inaccuracy with which d(S) (PBR) is estimated, whereas evidence of positive selection was observed for most of the DQB PBR sequences. A parallel analysis with CodeML allowed us to demonstrate, in a very consistent way, the occurrence of positive selection in the PBR-encoding region of both DQA and DQB genes. Moreover, we have identified several DQA (alpha47, alpha55, alpha56, alpha68, alpha69, alpha76 and alpha79) and DQB (beta9, beta26 and beta57) codons that appear to be under positive selection in different, and often unrelated, mammalian species. Non-synonymous polymorphism at these sites has been evolutionarily conserved meaning that it might have functional consequences on peptide binding.

Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics

We investigated the factors mediating selection acting on two MHC class II genes (DQA and DRB) in water vole (Arvicola scherman) natural populations in the French Jura Mountains. Population genetics showed significant homogeneity in allelic frequencies at the DQA1 locus as opposed to neutral markers (nine microsatellites), indicating balancing selection acting on this gene. Moreover, almost exhaustive screening for parasites, including gastrointestinal helminths, brain coccidia and antibodies against viruses responsible for zoonoses, was carried out. We applied a co-inertia approach to the genetic and parasitological data sets to avoid statistical problems related to multiple testing. Two alleles, Arte-DRB-11 and Arte-DRB-15, displayed antagonistic associations with the nematode Trichuris arvicolae, revealing the potential parasite-mediated selection acting on DRB locus. Selection mechanisms acting on the two MHC class II genes thus appeared different. Moreover, overdominance as balancing selection mechanism was showed highly unlikely in this system.

Diversity and evolutionary history of the MHC DQA gene in leporids

The European rabbit (Oryctolagus cuniculus) is used as a model for many human diseases, yet comparatively little is known of its genetics, particularly at important loci such as the major histocompatibility complex (MHC). This study investigated genetic diversity and evolutionary history of the DQA gene in a range of leporid species by analysing coding sequence diversity of exon 2 and intron 2 in 53 individuals of 16 different species. Fifty leporid DQA alleles were detected, including 13 novel European rabbit alleles. In the rabbit, the highest levels of diversity were observed in wild rabbits from Portugal, with wild rabbits from England and domestic rabbits showing less diversity. Within the sample, several recombination events were detected and trans-specific evolution of alleles was evidenced, both being general characteristics of mammalian MHC genes. Positive selection is implicated as operating on six codons within exon 2, which are also subject to positive selection in other mammals. Some of these positions are putative antigen recognition sites and underline the importance of pathogen-driven selection on these MHC genes.

Selection, diversity and evolutionary patterns of the MHC class II DAB in free-ranging Neotropical marsupials

BACKGROUND

Research on the genetic architecture and diversity of the MHC has focused mainly on eutherian mammals, birds and fish. So far, studies on model marsupials used in laboratory investigations indicated very little or even no variation in MHC class II genes. However, natural levels of diversity and selection are unknown in marsupials as studies on wild populations are virtually absent. We used two endemic South American mouse opossums, Gracilinanus microtarsus and Marmosops incanus, to investigate characteristic features of MHC selection. This study is the first investigation of MHC selection in free-ranging Neotropical marsupials. In addition, the evolutionary history of MHC lineages within the group of marsupials was examined.

RESULTS

G. microtarsus showed extensive levels of MHC diversity within and among individuals as 47 MHC-DAB alleles and high levels of sequence divergence were detected at a minimum of four loci. Positively selected codon sites were identified, of which most were congruent with human antigen binding sites. The diversity in M. incanus was rather low with only eight observed alleles at presumably two loci. However, these alleles also revealed high sequence divergence. Again, positive selection was identified on specific codon sites, all congruent with human ABS and with positively selected sites observed in G. microtarsus. In a phylogenetic comparison alleles of M. incanus interspersed widely within alleles of G. microtarsus with four alleles being present in both species.

CONCLUSION

Our investigations revealed extensive MHC class II polymorphism in a natural marsupial population, contrary to previous assumptions. Furthermore, our study confirms for the first time in marsupials the presence of three characteristic features common at MHC loci of eutherian mammals, birds and fish: large allelic sequence divergence, positive selection on specific sites and trans-specific polymorphism.

MH class IIalpha polymorphism in local and global adaptation of Arctic charr (Salvelinus alpinus L.)

Arctic charr, a highly plastic salmonid that inhabits the circumpolar region, colonized its current environment after the last glaciation. Recent colonization limits the capacity of many techniques to define and characterize constituent populations. As a novel approach, we used the major histocompatibility (MH) class IIalpha gene polymorphism as a marker that would characterize the genetic divergence of global Arctic charr populations caused by drift and by local adaptation to pathogens. We were able to detect significant isolation of all the lineages previously defined by mitochondrial DNA sequencing and also isolation of some populations within those groups. We found that most of the polymorphism of the class IIalpha gene was distributed globally, which indicates ancestral selection; however, in most cases, distinctive allele frequencies and specific haplotypes distinguished each population suggesting that recent selection has also occurred. Although all studied populations showed similar MH class IIalpha polymorphisms, we also found variation in which particular amino acid positions were polymorphic and which were constant in the different populations studied. This variation provides a greater adaptive capacity for the MH class IIalpha receptors in Arctic charr and is yet another illustration of the extraordinary plasticity of the species.

Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents

Parasite-mediated selection on major histocompatibility complex (MHC) genes has mainly been explored at the intraspecific level, although many molecular studies have revealed trans-species polymorphism. Interspecific patterns of MHC diversity might reveal factors responsible for the long-term evolution of MHC polymorphism. We hypothesize that host taxa harbouring high parasite diversity should exhibit high levels of MHC genetic diversity. We test this assumption using data on rodent species and their helminth parasites compiled from the literature. Controlling for similarity due to common descent, we present evidence indicating that high helminth species richness in rodent species is associated with increased MHC class II polymorphism. Our results are consistent with the idea that parasites sharing a long-term coevolutionary history with their hosts are the agents of selection explaining MHC polymorphism.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.