Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus

MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs

West African Mastomys rodents are the primary reservoir of the zoonotic Lassa virus (LASV). The virus causes haemorrhagic Lassa fever and considerable mortality in humans. To date, the role of Mastomys immunogenetics in resistance to, and persistence of, LASV infections is largely unknown. Here, we investigated the role of Major Histocompatibility Complex class I (MHC-I) on LASV infection status (i.e., active vs. cleared infection, determined via PCR and an immunofluorescence assay on IgG antibodies, respectively) in Mastomys natalensis and M. erythroleucus sampled within southwestern Nigeria. We identified more than 190 and 90 MHC-I alleles by Illumina high throughput-sequencing in M. natalensis and M. erythroleucus, respectively, with different MHC allele compositions and frequencies between LASV endemic and non-endemic sites. In M. natalensis, the MHC allele ManaMHC-I*006 was negatively associated with active infections (PCR-positive) and positively associated with cleared infections (IgG-positive) simultaneously, suggesting efficient immune responses that facilitate LASV clearance in animals carrying this allele. Contrarily, alleles ManaMHC-I*008 and ManaMHC-I*021 in M. natalensis, and MaerMHC-I*008 in M. erythroleucus, were positively associated with active infection, implying susceptibility. Alleles associated with susceptibility shared a glutamic acid at the positively selected codon 57, while ManaMHC-I*006 featured an arginine. There was no link between number of MHC alleles per Mastomys individual and LASV prevalence. Thus, specific alleles, but not MHC diversity per se, seem to mediate antibody responses to viremia. We conclude that co-evolution with LASV likely shaped the MHC-I diversity of the main LASV reservoirs in southwestern Nigeria, and that information on reservoir immunogenetics may hold insights into transmission dynamics and zoonotic spillover risks.

Microsatellite markers of the major histocompatibility complex genomic region of domestic camels

We identified and characterized 11 polymorphic microsatellite markers suitable for routine testing (three in the MHC class I sub-region, four in MHC class II and four in the MHC class III sub-region) of dromedaries and Bactrian camels. In total, 38 dromedaries and 33 Bactrian camels were genotyped, and interspecific differences were observed in the numbers of alleles and in allelic frequencies, as well as in the observed heterozygosity. These loci may be used as markers to study the adaptive genetic diversity of the MHC region in Old World camels.

High MHC diversity confers no advantage for phenotypic quality and reproductive performance in a wild bird

1. Genes of the major histocompatibility complex (MHC) encode antigen binding molecules and are an integral part of the acquired immune response of vertebrates. In general, high individual MHC diversity is expected to increase fitness by broadening the spectrum of pathogens recognized by the immune system, in accordance with the heterozygote advantage mechanism. On the other hand, the optimality hypothesis assumes that individuals with optimal (intermediate), rather than maximum diversity of the MHC will achieve the highest fitness because of inherent costs associated with expressing diverse MHC alleles. 2. Here, we tested for associations between individual diversity of the MHC class I and class II genes (binding antigens of intra- and extra-cellular pathogens, respectively) and a range of fitness-related traits (condition, ornament expression and reproduction) in an urban population of the Eurasian coot Fulica atra. 3. Contrary to our expectation, we found that high within-individual allelic diversity of MHC genes (both class I and II) was associated with poorer condition (lower blood haemoglobin concentrations), weaker expression of the putative ornament (smaller frontal shield), later onset of breeding and smaller clutches. An analysis of functional MHC allele clusters (supertypes) provided further support for negative associations of MHC diversity with phenotypic quality and reproductive performance, but most of these relationships could not be explained by the presence of specific maladaptive supertypes. Finally, we found little empirical support for the optimality hypothesis in the Eurasian coot. 4. Our results suggest that the costs of high MHC diversity outweighed any benefits associated with broad MHC repertoire, which could be driven by depauperate pathogen diversity in an urban landscape. To the best of our knowledge, this is one of the first studies providing consistent evidence for negative associations of MHC diversity with a range of fitness-related traits in a natural avian population.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele–pathogen interactions. Specific MHC allele–pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

Pathogen-associated selection on innate immunity genes (TLR4, TLR7) in a neotropical rodent in landscapes differing in anthropogenic disturbance

Toll-like receptors (TLRs) form part of the innate immune system and can recognize structurally conserved pathogen-associated molecular pattern (PAMP) molecules. Their functional importance in the resistance to pathogens has been documented in laboratory experimental settings and in humans. TLR diversity, however, has been rarely investigated in wildlife species. How the genetic diversity of TLRs is associated with various pathogens and how it is shaped by habitat disturbance are understudied. Therefore, we investigated the role of genetic diversity in the functionally important parts of TLR4 and TLR7 genes in resistance towards gastrointestinal nematodes and Hepacivirus infection. We chose a generalist study species, the rodent Proechimys semispinosus, because it is highly abundant in three Panamanian landscapes that differ in their degree of anthropogenic modification. We detected only two TLR7 haplotypes that differed by one synonymous single-nucleotide polymorphism (SNP) position. The TLR4 variability was higher, and we detected four TLR4 haplotypes that differed at one synonymous SNP and at three amino acid positions within the leucine-rich repeat region. Only TLR4 haplotypes had different frequencies in each landscape. Using generalized linear models, we found evidence that nematode loads and virus prevalence were influenced by both specific TLR4 haplotypes and landscape. Here, the variable “landscape” served as a surrogate for the important influential ecological factors distinguishing landscapes in our study, i.e. species diversity and host population density. Individuals carrying the common TLR4_Ht1 haplotype were less intensely infected by the most abundant strongyle nematode. Individuals carrying the rare TLR4_Ht3 haplotype were all Hepacivirus-positive, where those carrying the rare haplotype TLR4_Ht4 were less often infected by Hepacivirus than individuals with other haplotypes. Our study highlights the role of TLR diversity in pathogen resistance and the importance of considering immune genetic as well as ecological factors in order to understand the effects of anthropogenic changes on wildlife health.

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.

MHC class IIα polymorphism and its association with resistance/susceptibility to Vibrio harveyi in golden pompano (Trachinotus ovatus)

The major histocompatibility complex (MHC) plays an important role in the vertebrate immune response to antigenic peptides, and it is essential for recognizing foreign pathogens in organisms. In this study, MHC class IIα (Trov-MHC IIα) from the golden pompano (Trachinotus ovatus) was first cloned and identified. The gene structure of Trov-MHC IIα was contained four exons and three introns. High levels of polymorphism were found in the exon 2 of Trov-MHC IIα. A total of 29 different MHC class IIα alleles with high polymorphism were identified from 80 individuals. The ratio of non-synonymous substitutions (dN) to synonymous substitutions (dS) was 3.157 (>1) in the peptide binding regions (PBRs) of Trov-MHC IIα, suggesting positive balancing selection. Six alleles were selected to analyze the association between alleles and resistance/susceptibility to Vibrio harveyi in golden pompano. The results showed that Trov-DAA*6401 and Trov-DAA*6702 alleles were associated with the resistance to V. harveyi in golden pompano, while alleles Trov-DAA*6304 and Trov-DAA*7301 were associated with the susceptibility to V. harveyi in golden pompano. This study confirmed the association between alleles of MHC class IIα and disease resistance, and also detected some alleles which might be correlated with high V. harveyi-resistance. These disease resistance-related MHC alleles could be used as potential genetic markers for molecular marker-assisted selective breeding in the golden pompano.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Identification, polymorphism and expression of MHC class Iα in golden pompano, Trachinotus ovatus

The classical major histocompatibility complex class I (MHC I) plays a vital role in the immune system. In this study, we cloned and identified golden pompano (Trachinotus ovatus) MHC Iα (Trov-MHC Iα), which encodes 351 amino acid residues including a leader peptide, α1, α2, α3 domain, a transmembrane region and a cytoplasmic domain. Twenty six different sequences, which encoded various numbers of amino acid residues ranging from 348 to 354, were obtained from 12 individuals. Highly genetic polymorphism was found in the Trov-MHC Iα, especially in the α1 and α2 domains. Meanwhile, in the α1 and α2 domains, 21 positive selected positions were revealed by site models, indicating the diversity of Trov-MHC Iα may be mainly generated by positive selection. Moreover, quantitative real-time reverse transcription PCR and western blotting analyses demonstrated that Trov-MHC Iα was ubiquitously expressed in the nine tested tissues and more highly expressed in intestine, head kidney, gill, and spleen. In the head kidney and spleen, Trov-MHC Iα was significantly upregulated under LPS or poly I:C stimulation. The results of this study provide valuable insight into molecular polymorphism, evolutionary mechanism, expression and function of MHC Iα in the immune system of golden pompano.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.

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.

Identification of 48 full-length MHC-DAB functional alleles in miiuy croaker and evidence for positive selection

Major histocompatibility complex (MHC) molecules play a vital role in the immune response and are a highly polymorphic gene superfamily in vertebrates. As the molecular marker associated with polymorphism and disease susceptibility/resistance, the polymorphism of MHC genes has been investigated in many tetrapods and teleosts. Most studies were focused on the polymorphism of the second exon, which encodes the peptide-binding region (PBR) in the α1- or β1-domain, but few studies have examined the full-length coding region. To comprehensive investigate the polymorphism of MHC gene, we identified 48 full-length miiuy croaker (Miichthys miiuy) MHC class IIB (Mimi-DAB) functional alleles from 26 miiuy croaker individuals. All of the alleles encode 34 amino acid sequences, and a high level of polymorphism was detected in Mimi-DAB alleles. The rate of non-synonymous substitutions (dN) occurred at a significantly higher frequency than that of synonymous substitutions (dS) in the PBR, and this result suggests that balancing selection maintains polymorphisms at the Mimi-DAB locus. Phylogenetic analysis based on the full-length and exon 2 sequences of Mimi-DAB alleles both showed that the Mimi-DAB alleles were clustered into two major groups. A total of 19 positive selected sites were identified on the Mimi-DAB alleles after testing for positive selection, and 14 sites were predicted to be associated with antigen-binding sites, which suggests that most of selected sites are significant for disease resistance. The polymorphism of Mimi-DAB alleles provides an important resource for analyzing the association between the polymorphism of MHC gene and disease susceptibility/resistance, and for researching the molecular selective breeding of miiuy croaker with enhanced disease resistance.

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 Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.

Co-option and evolution of non-olfactory proteinaceous pheromones in a terrestrial lungless salamander

Gene co-option is a major force in the evolution of novel biological functions. In plethodontid salamanders, males deliver proteinaceous courtship pheromones to the female olfactory system or transdermally to the bloodstream. Molecular studies identified three families of highly duplicated, rapidly evolving pheromones (PRF, PMF, and SPF). Analyses for Plethodon salamanders revealed pheromone mixtures of primarily PRF and PMF. The current study demonstrates that in Desmognathus ocoee--a plesiomorphic species with transdermal delivery--SPF is the major pheromone component representing >30% of total protein. Chromatographic profiles of D. ocoee pheromones were consistent from May through October. LC/MS-MS analysis suggested uniform SPF isoform expression between individual male D. ocoee. A gene ancestry for SPF with the Three-Finger Protein superfamily was supported by intron-exon boundaries, but not by the disulfide bonding pattern. Further analysis of the pheromone mixture revealed paralogs to peptide hormones that contained mutations in receptor binding regions, such that these novel molecules may alter female physiology by acting as hormone agonists/antagonists. Cumulatively, gene co-option, duplication, and neofunctionalization have permitted recruitment of additional gene families for pheromone activity. Such independent co-option events may be playing a key role in salamander speciation by altering male traits that influence reproductive success.

Maintaining microendemic primate species along an environmental gradient - parasites as drivers for species differentiation

Understanding the drivers of species adaptations to changing environments on the one hand and the limits for hybridization on the other hand is among the hottest questions in evolutionary biology. Parasites represent one of the major selective forces driving host evolution and at least those with free-living stages are at the same time dependent on the ecological conditions of their host's habitat. Local immunological adaptations of host species to varying parasite pressure are therefore expected and might represent the genetic basis for ecological speciation and the maintenance of recently diverged species. Madagascar provides one of the rare examples where two partially sympatric primate species (Microcebus griseorufus, M. murinus) and their hybrids, as well as an allopatric species (M. cf rufus) live in close proximity along a very steep environmental gradient ranging from southern dry spiny bush to gallery forest to evergreen eastern humid rain forest, thus mimicking the situation encountered during extensions and retreats of vegetation formations under changing climatic conditions. This system was used to study parasite infection and immune gene (MHC) adaptations to varying parasite pressure that might provide selective advantages to pure species over hybrids. Parasite burdens increased with increasing humidity. M. griseorufus, M. murinus, and their hybrids but not M. rufus shared the same MHC alleles, indicating either retention of ancestral polymorphism or recent gene flow. The hybrids had much higher prevalence of intestinal parasites than either of the parent species living under identical environmental conditions. The different representation of parasites can indicate a handicap for hybrids that maintains species identities.

Major Histocompatibility Complex, demographic, and environmental predictors of antibody presence in a free-ranging mammal

Major Histocompatibility Complex (MHC) variability plays a key role in pathogen resistance, but its relative importance compared to environmental and demographic factors that also influence resistance is unknown. We analyzed the MHC II DRB exon 2 for 165 raccoons (Procyon lotor) in Missouri (USA). For each animal we also determined the presence of immunoglobulin G (IgG) and immunoglobulin M (IgM) antibodies to two highly virulent pathogens, canine distemper virus (CDV) and parvovirus. We investigated the role of MHC polymorphism and other demographic and environmental factors previously associated with predicting seroconversion. In addition, using an experimental approach, we studied the relative importance of resource availability and contact rates. We found important associations between IgG antibody presence and several MHC alleles and supertypes but not between IgM antibody presence and MHC. No effect of individual MHC diversity was found. For CDV, supertype S8, one allele within S8 (Prlo-DRB(∗)222), and a second allele (Prlo-DRB(∗)204) were positively associated with being IgG+, while supertype S4 and one allele within the supertype (Prlo-DRB(∗)210) were negatively associated with being IgG+. Age, year, and increased food availability were also positively associated with being IgG+, but allele Prlo-DRB(∗)222 was a stronger predictor. For parvovirus, only one MHC allele was negatively associated with being IgG+ and age and site were stronger predictors of seroconversion. Our results show that negative-frequency dependent selection is likely acting on the raccoon MHC and that while the role of MHC in relation to other factors depends on the pathogen of interest, it may be one of the most important factors predicting successful immune response.

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.

MHC-I affects infection intensity but not infection status with a frequent avian malaria parasite in blue tits

Host resistance against parasites depends on three aspects: the ability to prevent, control and clear infections. In vertebrates the immune system consists of innate and adaptive immunity. Innate immunity is particularly important for preventing infection and eradicating established infections at an early stage while adaptive immunity is slow, but powerful, and essential for controlling infection intensities and eventually clearing infections. Major Histocompatibility Complex (MHC) molecules are central in adaptive immunity, and studies on parasite resistance and MHC in wild animals have found effects on both infection intensity (parasite load) and infection status (infected or not). It seems MHC can affect both the ability to control infection intensities and the ability to clear infections. However, these two aspects have rarely been considered simultaneously, and their relative importance in natural populations is therefore unclear. Here we investigate if MHC class I genotype affects infection intensity and infection status with a frequent avian malaria infection Haemoproteus majoris in a natural population of blue tits Cyanistes caeruleus. We found a significant negative association between a single MHC allele and infection intensity but no association with infection status. Blue tits that carry a specific MHC allele seem able to suppress H. majoris infection intensity, while we have no evidence that this allele also has an effect on clearance of the H. majoris infection, a result that is in contrast with some previous studies of MHC and avian malaria. A likely explanation could be that the clearance rate of avian malaria parasites differs between avian malaria lineages and/or between avian hosts.

Parasite-mediated selection of major histocompatibility complex variability in wild brandt's voles (Lasiopodomys brandtii) from Inner Mongolia, China

Background

Genes of the major histocompatibility complex (MHC) exhibit high levels of variability, which is believed to have arisen through pathogen-mediated selection. We investigated the relationship between parasite load and genetic diversity at selectively neutral, non-coding markers (microsatellites) and adaptive genetic variation at a functionally important part of the MHC in six independent natural populations of Brandt’s voles (Lasiopodomys brandtii) from two regions of the Xilingol Grassland area of Inner Mongolia.

Results

Two-hundred and fifty-two voles were screened for gastrointestinal parasites, and were assessed for genetic variation. Parasite screening was done through non-invasive fecal egg counts, while allelic diversity was determined via single-stranded conformation polymorphism and DNA sequencing. We detected eight distinct helminth egg morphotypes. A total of 10 microsatellite loci were genotyped and 19 unique MHC class II B alleles were isolated. The rate of nonsynonymous substitutions (d_N) exceeded the rate of synonymous substitutions (d_S) at putative antigen binding sites of DRB. Neutral and adaptive genetic diversity differed between the six vole populations. To test the main pathogen-driven selection hypotheses for the maintenance of host MHC diversity and parasite species-specific co-evolutionary effects, multivariate approaches (generalized linear mixed models) were used to test for associations between the MHC class II DRB genotype and infections with nematodes. We found no evidence for heterozygote advantage, and overall heterozygosity was lower than expected in the MHC alleles. We identified an association between the parasite load and specific MHC alleles in the voles, and this pattern varied between geographic regions.

Conclusions

The results suggest that MHC variability in Brandt’s voles is maintained by rare allele advantage and fluctuating selection, but the data failed to show any heterozygote advantage effect. Our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations. In addition, our study contributes to the understanding of what maintains MHC diversity, of host-pathogen coevolution and of how genetic diversity is maintained in voles.

Identification and characterization of the major histocompatibility complex class II DQB (MhcMath-DQB1) alleles in Tibetan macaques (Macaca thibetana)

Tibetan macaque (Macaca thibetana), an endangered primate species endemic to China, have been used as experimental animal model for various human diseases. Major histocompatibility complex (MHC) genes play a crucial role in the susceptibility and/or resistance to many human diseases, but little is known about Tibetan macaques. To gain an insight into the MHC background and to facilitate the experimental use of Tibetan macaques, the second exon of Mhc-DQB1 gene was sequenced in a cohort of wild Tibetan macaques living in the Sichuan province of China. A total of 23 MhcMath-DQB1 alleles were identified for the first time, illustrating a marked allelic polymorphism at the DQB1 locus for these macaques. Most of the sequences (74%) observed in this study belong to DQB1*06 (9 alleles) and DQB1*18 (8 alleles) lineages, and the rest (26%) belong to DQB1*15 (3 alleles) and DQB1*17 (3 alleles) lineages. The most frequent alleles detected among these macaques were MhcMath-DQB1*15:02:02 (17.9%), followed by Math-DQB1*06:06, 17:03 and 18:01, which were detected in 9 (16.1%) of the monkeys, respectively. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the peptide-binding region, suggesting balancing selection for maintaining polymorphisms at the MHC class II DQB1 locus. Phylogenetic analyses confirms the trans-species model of evolution of the Mhc-DQB1 genes in non-human primates, and in particular, the extensive allele sharing is observed between Tibetan and other macaque species.

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.

Domestication does not narrow MHC diversity in Sus scrofa

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

The functional importance of sequence versus expression variability of MHC alleles in parasite resistance

Understanding selection processes driving the pronounced allelic polymorphism of the major histocompatibility complex (MHC) genes and its functional associations to parasite load have been the focus of many recent wildlife studies. Two main selection scenarios are currently debated which explain the susceptibility or resistance to parasite infections either by the effects of (1) specific MHC alleles which are selected frequency-dependent in space and time or (2) a heterozygote or divergent allele advantage. So far, most studies have focused only on structural variance in co-evolutionary processes although this might not be the only trait subject to natural selection. In the present study, we analysed structural variance stretching from exon1 through exon3 of MHC class II DRB genes as well as genotypic expression variance in relation to the gastrointestinal helminth prevalence and infection intensity in wild yellow-necked mice (Apodemus flavicollis). We found support for the functional importance of specific alleles both on the sequence and expression level. By resampling a previously investigated study population we identified specific MHC alleles affected by temporal shifts in parasite pressure and recorded associated changes in allele frequencies. The allele Apfl-DRB*23 was associated with resistance to infections by the oxyurid nematode Syphacia stroma and at the same time with susceptibility to cestode infection intensity. In line with our expectation, MHC mRNA transcript levels tended to be higher in cestode-infected animals carrying the allele Apfl-DRB*23. However, no support for a heterozygote or divergent allele advantage on the sequence or expression level was detected. The individual amino acid distance of genotypes did not explain individual differences in parasite loads and the genetic distance had no effect on MHC genotype expression. For ongoing studies on the functional importance of expression variance in parasite resistance, allele-specific expression data would be preferable.

Comparison of 454 pyrosequencing methods for characterizing the major histocompatibility complex of nonmodel species and the advantages of ultra deep coverage

Characterization and population genetic analysis of multilocus genes, such as those found in the major histocompatibility complex (MHC) is challenging in nonmodel vertebrates. The traditional method of extensive cloning and Sanger sequencing is costly and time-intensive and indirect methods of assessment often underestimate total variation. Here, we explored the suitability of 454 pyrosequencing for characterizing multilocus genes for use in population genetic studies. We compared two sample tagging protocols and two bioinformatic procedures for 454 sequencing through characterization of a 185-bp fragment of MHC DRB exon 2 in wolverines (Gulo gulo) and further compared the results with those from cloning and Sanger sequencing. We found 10 putative DRB alleles in the 88 individuals screened with between two and four alleles per individual, suggesting amplification of a duplicated DRB gene. In addition to the putative alleles, all individuals possessed an easily identifiable pseudogene. In our system, sequence variants with a frequency below 6% in an individual sample were usually artefacts. However, we found that sample preparation and data processing procedures can greatly affect variant frequencies in addition to the complexity of the multilocus system. Therefore, we recommend determining a per-amplicon-variant frequency threshold for each unique system. The extremely deep coverage obtained in our study (approximately 5000×) coupled with the semi-quantitative nature of pyrosequencing enabled us to assign all putative alleles to the two DRB loci, which is generally not possible using traditional methods. Our method of obtaining locus-specific MHC genotypes will enhance population genetic analyses and studies on disease susceptibility in nonmodel wildlife species.

Diversity and selection of MHC class IIb gene exon3 in Chinese alligator

Our study used MHC class IIb gene exon3 complete sequence as markers to investigate genetic variability, selection and population differentiation in Chinese alligator. In this study, 282 bp MHC IIb exon3 complete sequence was got, none of the sequences contained insertions/deletions or stop codons, suggesting that all sequences might come from functional molecules in the genome. The neighbor-joining (NJ) tree revealed that Xuangzhou and Changxing populations were genetically close related, while Wild population showed the most diverse from the other. Gene flow (Nm) was very higher than one, suggesting that inter-group gene flow may have occurred. Furthermore, the results showed that MHC IIb gene might be a good molecular marker, we think that this technology could be used for Chinese alligator breeding and releasing in future.

Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Populations of Arctic foxes (Vulpes lagopus) have been isolated on two of the Commander Islands (Bering and Mednyi) from the circumpolar distributed mainland population since the Pleistocene. In 1970-1980, an epizootic outbreak of mange caused a severe population decline on Mednyi Island. Genes of the major histocompatibility complex (MHC) play a primary role in infectious disease resistance. The main objectives of our study were to compare contemporary variation of MHC class II in mainland and island Arctic foxes, and to document the effects of the isolation and the recent bottleneck on MHC polymorphism by analyzing samples from historical and contemporary Arctic foxes. In 184 individuals, we found 25 unique MHC class II DRB and DQB alleles, and identified evidence of balancing selection maintaining allelic lineages over time at both loci. Twenty different MHC alleles were observed in mainland foxes and eight in Bering Island foxes. The historical Mednyi population contained five alleles and all contemporary individuals were monomorphic at both DRB and DQB. Our data indicate that despite positive and diversifying selection leading to elevated rates of amino acid replacement in functionally important antigen-binding sites, below a certain population size, balancing selection may not be strong enough to maintain genetic diversity in functionally important genes. This may have important fitness consequences and might explain the high pathogen susceptibility in some island populations. This is the first study that compares MHC diversity before and after a bottleneck in a wild canid population using DNA from museum samples.

Evaluation of two approaches to genotyping major histocompatibility complex class I in a passerine-CE-SSCP and 454 pyrosequencing

Genes of the highly dynamic major histocompatibility complex (MHC) are directly linked to individual fitness and are of high interest in evolutionary ecology and conservation genetics. Gene duplication and positive selection usually lead to high levels of polymorphism in the MHC region, making genotyping of MHC a challenging task. Here, we compare the performance of two methods for MHC class I genotyping in a passerine with highly duplicated MHC class I genes: capillary electrophoresis-single-strand conformation polymorphism (CE-SSCP) analysis and 454 GS FLX Titanium pyrosequencing. According to our findings, the number of MHC variants (called alleles for simplicity) detected by CE-SSCP is significantly lower than detected by 454. To resolve discrepancies between the two methods, we cloned and Sanger sequenced a MHC class I amplicon for an individual with high number of alleles. We found a perfect congruence between cloning/Sanger sequencing results and 454. Thus, in case of multi-locus amplification, CE-SSCP considerably underestimates individual MHC diversity. However, numbers of alleles detected by both methods are significantly correlated, although the correlation is weak (r = 0.32). Thus, in systems with highly duplicated MHC, 454 provides more reliable information on individual diversity than CE-SSCP.

Heligmosomoides polygyrus infection is associated with lower MHC class II gene expression in Apodemus flavicollis: indication for immune suppression?

Due to their key role in recognizing foreign antigens and triggering the subsequent immune response the genes of the major histocompatibility complex (MHC) provide a potential target for parasites to attack in order to evade detection and expulsion from the host. A diminished MHC gene expression results in less activated T cells and might serve as a gateway for pathogens and parasites. Some parasites are suspected to be immune suppressors and promote co-infections of other parasites even in other parts of the body. In our study we found indications that the gut dwelling nematode Heligmosomoides polygyrus might exert a systemic immunosuppressive effect in yellow-necked mice (Apodemus flavicollis). The amount of hepatic MHC class II DRB gene RNA transcripts in infected mice was negatively associated with infection intensity with H. polygyrus. The hepatic expression of immunosuppressive cytokines, such as transforming growth factor β and interleukin 10 was not associated with H. polygyrus infection. We did not find direct positive associations of H. polygyrus with other helminth species. But the prevalence and infection intensity of the nematodes Syphacia stroma and Trichuris muris were higher in multiple infected individuals. Furthermore, our data indicated antagonistic effects in the helminth community of A. flavicollis as cestode infection correlated negatively with H. polygyrus and helminth species richness. Our study shows that expression analyses of immune relevant genes can also be performed in wildlife, opening new aspects and possibilities for future ecological and evolutionary research.

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.

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.

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.

Major histocompatibility complex variation and age-specific endoparasite load in subadult European rabbits

Genes of the major histocompatibility complex (MHC) play a fundamental role in the vertebrate immune response and are amongst the most polymorphic genes in vertebrate genomes. It is generally agreed that the highly polymorphic nature of the MHC is maintained through host-parasite co-evolution. Two nonexclusive mechanisms of selection are supposed to act on MHC genes: superiority of MHC heterozygous individuals (overdominance) and an advantage for rare MHC alleles. However, the precise mechanisms and their relative importance are still unknown. Here, we examined MHC dependent parasite load in European rabbits (Oryctolagus cuniculus) from a distinct population with low MHC diversity (three alleles, six genotypes). Using a multivariate approach, we tested for associations of individual MHC class II DRB constitution and the rabbits' intestinal burden with nematodes and coccidia. Rabbits having a particular allele showed lower infestations with hepatic coccidia (E. stiedai). However, a comparison of all six genotypes in the population revealed that carriers of this allele only benefit when they are heterozygous, and furthermore, MHC heterozygosity in general did not affect individual parasite load. In conclusion, this study suggests an immunogenetic basis of European rabbit resistance to hepatic coccidiosis, which can strongly limit survival to maturity in this species. Our study gives a complex picture of MHC-parasite correlations, unveiling the limits of the classical hypotheses of how MHC polymorphism is maintained in natural systems.

Characterisation of MHC class II DRB genes in the northern tree shrew (Tupaia belangeri)

Genes of the major histocompatibility complex (MHC) mainly code for proteins of the immune system of jawed vertebrates. In particular, MHC class I and II cell surface proteins are crucial for the self/non-self discrimination of the adaptive immune system and are the most polymorphic genes in vertebrates. Positive selection, gene duplications and pseudogenes shape the face of the MHC and reflect a highly dynamic evolution. Here, we present for the first time data of the highly polymorphic MHC class II DRB exon 2 of a representative of the mammalian order scandentia, the northern tree shrew Tupaia belangeri. We found up to eight different alleles per individual and determined haplotype constitution by intensively studying their inheritance. The alleles were assigned to four putative loci, all of which were polymorphic. Only the most polymorphic locus was subject to positive selection within the antigen binding sites and only alleles of this locus were transcribed.

Natural selection coupled with intragenic recombination shapes diversity patterns in the major histocompatibility complex class II genes of the giant panda

Ample variations of the major histocompatibility complex (MHC) genes are essential for vertebrates to adapt to various environmental conditions. In this study, we investigated the genetic variations and evolutionary patterns of seven functional MHC class II genes (one DRA, two DRB, two DQA, and two DQB) of the giant panda. The results showed the presence of two monomorphic loci (DRA and DQB2) and five polymorphic loci with different numbers of alleles (seven at DRB1, six at DRB3, seven at DQA1, four at DQA2, six at DQB1). The presence of balancing selection in the giant panda was supported by the following pieces of evidence: (1) The observed heterozygosity was higher than expected. (2) Amino acid heterozygosity was significantly higher at antigen-binding sites (ABS) compared with non-ABS sequences. (3) The selection parameter omega (d(N)/d(S)) was significantly higher at ABS compared with non-ABS sequences. (4) Approximately 95.45% of the positively selected codons (P>0.95) were located at or adjacent to an ABS. Furthermore, this study showed that (1) The Qinling subspecies exhibited high omega values across each locus (all >1), supporting its extensive positive selection. (2) The Sichuan subspecies displayed small omega at DRB1 (omega<0.72) and DQA2 (omega<0.48), suggesting that these sites underwent strong purifying selection. (3) Intragenic recombination was detected in DRB1, DQA1, and DQB1. The molecular diversity in classic Aime-MHC class II genes implies that the giant panda had evolved relatively abundant variations in its adaptive immunity along the history of host-pathogen co-evolution. Collectively, these findings indicate that natural selection accompanied by recombination drives the contrasting diversity patterns of the MHC class II genes between the two studied subspecies of giant panda.

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.