Polymorphism, haplotype composition, and selection in the Mhc-DRB of wild baboons

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Relationship between genome-wide and MHC class I and II genetic diversity and complementarity in a nonhuman primate

Although mate choice is expected to favor partners with advantageous genetic properties, the relative importance of genome‐wide characteristics, such as overall heterozygosity or kinship, versus specific loci, is unknown. To disentangle genome‐wide and locus‐specific targets of mate choice, we must first understand congruence in global and local variation within the same individual. This study compares genetic diversity, both absolute and relative to other individuals (i.e., complementarity), assessed across the genome to that found at the major histocompatibility complex (MHC), a hyper‐variable gene family integral to immune system function and implicated in mate choice across species. Using DNA from 22 captive olive baboons (Papio anubis), we conducted double digest restriction site‐associated DNA sequencing to estimate genome‐wide heterozygosity and kinship, and sequenced two class I and two class II MHC loci. We found that genome‐wide diversity was not associated with MHC diversity, and that diversity at class I MHC loci was not correlated with diversity at class II loci. Additionally, kinship was a significant predictor of the number of MHC alleles shared between dyads at class II loci. Our results provide further evidence of the strong selective pressures maintaining genetic diversity at the MHC in comparison to other randomly selected sites throughout the genome. Furthermore, our results indicate that class II MHC disassortative mate choice may mediate inbreeding avoidance in this population. Our study suggests that mate choice favoring genome‐wide genetic diversity is not always synonymous with mate choice favoring MHC diversity, and highlights the importance of controlling for kinship when investigating MHC‐associated mate choice.

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.

Patterns of MHC-dependent sexual selection in a free-living population of sheep

The MHC is one of the most polymorphic gene clusters in vertebrates and play an essential role in adaptive immunity. Apart from pathogen-mediated selection, sexual selection can also contribute to the maintenance of MHC diversity. MHC-dependent sexual selection could occur via several mechanisms but at present there is no consensus as to which of these mechanisms are involved and their importance. Previous studies have often suffered from limited genetic and behavioural data and small sample size, and were rarely able to examine all the mechanisms together, determine whether signatures of MHC-based non-random mating are independent of genomic effects or differentiate whether MHC-dependent sexual selection takes place at the pre- or post-copulatory stage. In this study, we use Monte Carlo simulation to investigate evidence for non-random MHC-dependent mating patterns by all three mechanisms in a free-living population of Soay sheep. Using 1710 sheep diplotyped at the MHC class IIa region and genome-wide SNPs, together with field observations of consorts, we found sexual selection against a particular haplotype in males at the pre-copulatory stage and sexual selection against female MHC heterozygosity during the rut. We also found MHC-dependent disassortative mating at the post-copulatory stage, along with strong evidence of inbreeding avoidance at both stages. However, results from generalized linear mixed models suggest that the pattern of MHC-dependent disassortative mating could be a by-product of inbreeding avoidance. Our results therefore suggest that while multiple apparent mechanisms of non-random mating with respect to the MHC may occur, some of them have alternative explanations.

Genetic variation at MHC class II loci influences both olfactory signals and scent discrimination in ring-tailed lemurs

BACKGROUND

Diversity at the Major Histocompatibility Complex (MHC) is critical to health and fitness, such that MHC genotype may predict an individual's quality or compatibility as a competitor, ally, or mate. Moreover, because MHC products can influence the components of bodily secretions, an individual's body odors may signal its MHC composition and influence partner identification or mate choice. Here, we investigated MHC-based signaling and recipient sensitivity by testing for odor-gene covariance and behavioral discrimination of MHC diversity and pairwise dissimilarity in a strepsirrhine primate, the ring-tailed lemur (Lemur catta).

METHODS

First, we coupled genotyping of the MHC class II gene, DRB, with gas chromatography-mass spectrometry of genital gland secretions to investigate if functional genetic diversity is signaled by the chemical diversity of lemur scent secretions. We also assessed if the chemical similarity between individuals correlated with their MHC-DRB similarity. Next, we assessed if lemurs discriminated this chemically encoded, genetic information in opposite-sex conspecifics.

RESULTS

We found that both sexes signaled overall MHC-DRB diversity and pairwise MHC-DRB similarity via genital secretions, but in a sex- and season-dependent manner. Additionally, the sexes discriminated absolute and relative MHC-DRB diversity in the genital odors of opposite-sex conspecifics, suggesting that lemur genital odors function to advertise genetic quality.

CONCLUSIONS

In summary, genital odors of ring-tailed lemurs provide honest information about an individual's absolute and relative MHC quality. Complementing evidence in humans and Old World monkeys, we suggest that reliance on scent signals to communicate MHC quality may be important across the primate lineage.

Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca)

Few major histocompatibility complex (MHC)-based mate choice studies include all MHC genes at the inter-individual, sperm-egg, and mother-fetus recognition levels. We tested three hypotheses of female mate choice in a 17-year study of the giant panda (Ailuropoda melanoleuca) while using ten functional MHC loci (four MHC class I loci: Aime-C, Aime-F, Aime-I, and Aime-L; six MHC class II loci: Aime-DRA, Aime-DRB3, Aime-DQA1, Aime-DQA2, Aime-DQB1, and Aime-DQB2); five super haplotypes (SuHa, SuHaI, SuHaII, DQ, and DR); and, seven microsatellites. We found female choice for heterozygosity at Aime-C, Aime-I, and DQ and for disassortative mate choice at Aime-C, DQ, and DR at the inter-individual recognition level. High mating success occurred in MHC-dissimilar mating pairs. No significant results were found based on any microsatellite parameters, suggesting that MHCs were the mate choice target and there were no signs of inbreeding avoidance. Our results indicate Aime-DQA1- and Aime-DQA2-associated disassortative selection at the sperm-egg recognition level and a possible Aime-C- and Aime-I-associated assortative maternal immune tolerance mechanism. The MHC genes were of differential importance at the different recognition levels, so all of the functional MHC genes should be included when studying MHC-dependent reproductive mechanisms.

Lack of evidence for selection favouring MHC haplotypes that combine high functional diversity

High rates of gene duplication and the highest levels of functional allelic diversity in vertebrate genomes are the main hallmarks of the major histocompatibility complex (MHC), a multigene family with a primordial role in pathogen recognition. The usual tight linkage among MHC gene duplicates may provide an opportunity for the evolution of haplotypes that associate functionally divergent alleles and thus grant the transmission of optimal levels of diversity to coming generations. Even though such associations may be a crucial component of disease resistance, this hypothesis has been given little attention in wild populations. Here, we leveraged pedigree data from a barn owl (Tyto alba) population to characterize MHC haplotype structure across two MHC class I (MHC-I) and two MHC class IIB (MHC-IIB) duplicates, in order to test the hypothesis that haplotypes' genetic diversity is higher than expected from randomly associated alleles. After showing that MHC loci are tightly linked within classes, we found limited evidence for shifts towards MHC haplotypes combining high diversity. Neither amino acid nor functional within-haplotype diversity were significantly higher than in random sets of haplotypes, regardless of MHC class. Our results therefore provide no evidence for selection towards high-diversity MHC haplotypes in barn owls. Rather, high rates of concerted evolution may constrain the evolution of high-diversity haplotypes at MHC-I, while, in contrast, for MHC-IIB, fixed differences among loci may provide barn owls with already optimized functional diversity. This suggests that at the MHC-I and MHC-IIB respectively, different evolutionary dynamics may govern the evolution of within-haplotype diversity.

Genetic wealth, population health: Major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring-tailed lemur (Lemur catta) as a model, we compared two populations under long-term study: a relatively "open," wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003-2013) and a "closed," captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980-2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC-DRB diversity and, across populations, we compared the number of unique MHC-DRB alleles and their distributions. Wild individuals possessed more MHC-DRB alleles than did captive individuals, and overall, the wild population had more unique MHC-DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive-breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.

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

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

A quick and robust MHC typing method for free-ranging and captive primate species

Gene products of the major histocompatibility complex (MHC) of human and non-human primates play a crucial role in adaptive immunity, and most of the relevant genes not only show a high degree of variability (polymorphism) but also copy number variation (CNV) is observed. Due to this diversity, MHC proteins influence the capability of individuals to cope with various pathogens. MHC and/or MHC-linked gene products such as odorant receptor genes are thought to influence mate choice and reproductive success. Therefore, MHC typing of wild and captive primate populations is considered to be useful in conservation biology, which is, however, often hampered by the need of invasive and time-consuming methods. All intact Mhc-DRB genes in primates appear to possess a complex and highly divergent microsatellite, DRB-STR. A panel of 154 pedigreed olive baboons (Papio anubis) was examined for their DRB content by DRB-STR analysis of genomic DNA. Using the same methodology on DNA of feces samples, DRB variability of a silvery gibbon population (Hylobates moloch) (N = 24), an endangered species, could successfully be studied. In both species, length determination of the DRB-STR resulted in the definition of unique genotyping patterns that appeared to be specific for a certain chromosome. Moreover, the different STR lengths were shown to segregate with the allelic variation of the respective gene. The results obtained expand data gained previously on DRB-STR typing in macaques, great apes, and humans and strengthen the conclusion that this protocol is applicable in molecular ecology, conservation biology, and colony management, especially of endangered primate species.

No evidence of prenatal diversifying selection at locus or supertype levels in the dog MHC class II loci

BACKGROUND

Despite decades of studying, the mechanisms maintaining high diversity in the genes of the Major Histocompatibility Complex (MHC) are still puzzling scientists. In addition to pathogen recognition and other functions, MHC molecules may act prenatally in mate choice and in maternal-foetal interactions. These interactions are potential selective mechanisms that increase genetic diversity in the MHC. During pregnancy, immune response has a dual role: the foetus represents foreign tissue compared to mother, but histo-incompatibility is required for successful pregnancy. We have studied the prenatal selection in MHC class II loci (DLA-DQA1, DLA-DQB1 and DLA-DRB1) in domestic dogs by comparing the observed and expected offspring genotype proportions in 110 dog families. Several potential selection targets were addressed, including the peptide-binding site, the MHC locus, three-locus haplotype and supertype levels. For the supertype analysis, the first canine supertype classification was created based on in silico analysis of peptide-binding amino-acid polymorphism.

RESULTS

In most loci and levels, no deviation from the expected genotype frequencies was observed. However, one peptide-binding site in DLA-DRB1 had an excess of heterozygotes among the offspring. In addition, if the father shared a DLA-DRB1 allele with the mother, that allele was inherited by the offspring more frequently than expected, suggesting the selective advantage of a histo-compatible foetus, in contrast to our expectations.

CONCLUSIONS

We conclude that there is some evidence of post-copulatory selection at nucleotide site level in the MHC loci of pet dogs. But due to no indication of selection at locus, three-locus, or supertype levels, we estimated that the prenatal selection coefficient is less than 0.3 in domestic dogs and very likely other factors are more important in maintaining the genetic diversity in MHC loci.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

Major histocompatibility complex and mate choice in the polygynous primate: the Sichuan snub-nosed monkey (Rhinopithecus roxellana)

The highly polymorphic genes within the major histocompatibility complex (MHC) not only play a major role in immunity resistance, but also seem to provide hints for mate choice in some animal populations. In the present study we investigated MHC-related mate choice in a small natural population (group size 40-55 individuals) of a polygynous primate, the Sichuan snub-nosed monkey (Rhinopithecus roxellana). We found that there was no evidence either for MHC-disassortative mating, or for females to mate with males based on MHC heterozygosity or specific alleles. Nevertheless, of the 11 alleles identified, we found that the frequencies of 2 alleles, Rhro-DRB2 (P < 0.01) and Rhro-DRB5 (P < 0.05) were higher in offspring than in their parents. These findings suggest that MHC-DRB in this population of R. roxellana is unlikely to be associated with mating preferences. Limited female opportunities for mate choice are likely due, in part, to the harem breeding structure present in R. roxellana, and the relatively small number of resident adult males in our study band (N = 4-6). In addition, we suggest that differences in the frequency of particular alleles across generations may be linked to parasite resistance in a fluctuating environment; however, confirmation of this finding requires further study.

Signature Patterns of MHC Diversity in Three Gombe Communities of Wild Chimpanzees Reflect Fitness in Reproduction and Immune Defense against SIVcpz

Major histocompatibility complex (MHC) class I molecules determine immune responses to viral infections. These polymorphic cell-surface glycoproteins bind peptide antigens, forming ligands for cytotoxic T and natural killer cell receptors. Under pressure from rapidly evolving viruses, hominoid MHC class I molecules also evolve rapidly, becoming diverse and species-specific. Little is known of the impact of infectious disease epidemics on MHC class I variant distributions in human populations, a context in which the chimpanzee is the superior animal model. Population dynamics of the chimpanzees inhabiting Gombe National Park, Tanzania have been studied for over 50 years. This population is infected with SIVcpz, the precursor of human HIV-1. Because HLA-B is the most polymorphic human MHC class I molecule and correlates strongly with HIV-1 progression, we determined sequences for its ortholog, Patr-B, in 125 Gombe chimpanzees. Eleven Patr-B variants were defined, as were their frequencies in Gombe's three communities, changes in frequency with time, and effect of SIVcpz infection. The growing populations of the northern and central communities, where SIVcpz is less prevalent, have stable distributions comprising a majority of low-frequency Patr-B variants and a few high-frequency variants. Driving the latter to high frequency has been the fecundity of immigrants to the northern community, whereas in the central community, it has been the fecundity of socially dominant individuals. In the declining population of the southern community, where greater SIVcpz prevalence is associated with mortality and emigration, Patr-B variant distributions have been changing. Enriched in this community are Patr-B variants that engage with natural killer cell receptors. Elevated among SIVcpz-infected chimpanzees, the Patr-B*06:03 variant has striking structural and functional similarities to HLA-B*57, the human allotype most strongly associated with delayed HIV-1 progression. Like HLA-B*57, Patr-B*06:03 correlates with reduced viral load, as assessed by detection of SIVcpz RNA in feces.

Selection on MHC class II supertypes in the New Zealand endemic Hochstetter's frog

BACKGROUND

The New Zealand native frogs, family Leiopelmatidae, are among the most archaic in the world. Leiopelma hochstetteri (Hochstetter's frog) is a small, semi-aquatic frog with numerous, fragmented populations scattered across New Zealand's North Island. We characterized a major histocompatibility complex (MHC) class II B gene (DAB) in L. hochstetteri from a spleen transcriptome, and then compared its diversity to neutral microsatellite markers to assess the adaptive genetic diversity of five populations ("evolutionarily significant units", ESUs).

RESULTS

L. hochstetteri possessed very high MHC diversity, with 74 DAB alleles characterized. Extremely high differentiation was observed at the DAB locus, with only two alleles shared between populations, a pattern that was not reflected in the microsatellites. Clustering analysis on putative peptide binding residues of the DAB alleles indicated four functional supertypes, all of which were represented in 4 of 5 populations, albeit at different frequencies. Otawa was an exception to these observations, with only two DAB alleles present.

CONCLUSIONS

This study of MHC diversity highlights extreme population differentiation at this functional locus. Supertype differentiation was high among populations, suggesting spatial and/or temporal variation in selection pressures. Low DAB diversity in Otawa may limit this population's adaptive potential to future pathogenic challenges.

454 screening of individual MHC variation in an endemic island passerine

Genes of the major histocompatibility complex (MHC) code for receptors that are central to the adaptive immune response of vertebrates. These genes are therefore important genetic markers with which to study adaptive genetic variation in the wild. Next-generation sequencing (NGS) has increasingly been used in the last decade to genotype the MHC. However, NGS methods are highly prone to sequencing errors, and although several methodologies have been proposed to deal with this, until recently there have been no standard guidelines for the validation of putative MHC alleles. In this study, we used the 454 NGS platform to screen MHC class I exon 3 variation in a population of the island endemic Berthelot's pipit (Anthus berthelotii). We were able to characterise MHC genotypes across 309 individuals with high levels of repeatability. We were also able to determine alleles that had low amplification efficiencies, whose identification within individuals may thus be less reliable. At the population level we found lower levels of MHC diversity in Berthelot's pipit than in its widespread continental sister species the tawny pipit (Anthus campestris), and observed trans-species polymorphism. Using the sequence data, we identified signatures of gene conversion and evidence of maintenance of functionally divergent alleles in Berthelot's pipit. We also detected positive selection at 10 codons. The present study therefore shows that we have an efficient method for screening individual MHC variation across large datasets in Berthelot's pipit, and provides data that can be used in future studies investigating spatio-temporal patterns and scales of selection on the MHC.

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.

Maintaining functional major histocompatibility complex diversity under inbreeding: the case of a selfing vertebrate

Major histocompatibility complex (MHC) genes encode proteins that present pathogen-derived antigens to T-cells, initiating the adaptive immune response in vertebrates. Although populations with low MHC diversity tend to be more susceptible to pathogens, some bottlenecked populations persist and even increase in numbers despite low MHC diversity. Thus, the relative importance of MHC diversity versus genome-wide variability for the long-term viability of populations after bottlenecks and/or under high inbreeding is controversial. We tested the hypothesis that genome-wide inbreeding (estimated using microsatellites) should be more critical than MHC diversity alone in determining pathogen resistance in the self-fertilizing fish Kryptolebias marmoratus by analysing MHC diversity and parasite loads in natural and laboratory populations with different degrees of inbreeding. Both MHC and neutral diversities were lost after several generations of selfing, but we also found evidence of parasite selection acting on MHC diversity and of non-random loss of alleles, suggesting a possible selective advantage of those individuals with functionally divergent MHC, in accordance with the hypothesis of divergent allele advantage. Moreover, we found that parasite loads were better explained by including MHC diversity in the model than by genome-wide (microsatellites) heterozygosity alone. Our results suggest that immune-related overdominance could be the key in maintaining variables rates of selfing and outcrossing in K. marmoratus and other mixed-mating species.

Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing

The critical role of major histocompatibility complex (MHC) genes in disease resistance, along with their putative function in sexual selection, reproduction and chemical ecology, make them an important genetic system in evolutionary ecology. Studying selective pressures acting on MHC genes in the wild nevertheless requires population-wide genotyping, which has long been challenging because of their extensive polymorphism. Here, we report on large-scale genotyping of the MHC class II loci of the grey mouse lemur (Microcebus murinus) from a wild population in western Madagascar. The second exons from MHC-DRB and -DQB of 772 and 672 individuals were sequenced, respectively, using a 454 sequencing platform, generating more than 800,000 reads. Sequence analysis, through a stepwise variant validation procedure, allowed reliable typing of more than 600 individuals. The quality of our genotyping was evaluated through three independent methods, namely genotyping the same individuals by both cloning and 454 sequencing, running duplicates, and comparing parent-offspring dyads; each displaying very high accuracy. A total of 61 (including 20 new) and 60 (including 53 new) alleles were detected at DRB and DQB genes, respectively. Both loci were non-duplicated, in tight linkage disequilibrium and in Hardy-Weinberg equilibrium, despite the fact that sequence analysis revealed clear evidence of historical selection. Our results highlight the potential of 454 sequencing technology in attempts to investigate patterns of selection shaping MHC variation in contemporary populations. The power of this approach will nevertheless be conditional upon strict quality control of the genotyping data.

Characterization and 454 pyrosequencing of major histocompatibility complex class I genes in the great tit reveal complexity in a passerine system

BACKGROUND

The critical role of Major Histocompatibility Complex (Mhc) genes in disease resistance and their highly polymorphic nature make them exceptional candidates for studies investigating genetic effects on survival, mate choice and conservation. Species that harbor many Mhc loci and high allelic diversity are particularly intriguing as they are potentially under strong selection and studies of such species provide valuable information as to the mechanisms maintaining Mhc diversity. However comprehensive genotyping of complex multilocus systems has been a major challenge to date with the result that little is known about the consequences of this complexity in terms of fitness effects and disease resistance.

RESULTS

In this study, we genotyped the Mhc class I exon 3 of the great tit (Parus major) from two nest-box breeding populations near Oxford, UK that have been monitored for decades. Characterization of Mhc class I exon 3 was adopted and bidirectional sequencing was carried using the 454 sequencing platform. Full analysis of sequences through a stepwise variant validation procedure allowed reliable typing of more than 800 great tits based on 214,357 reads; from duplicates we estimated the repeatability of typing as 0.94. A total of 862 alleles were detected, and the presence of at least 16 functional loci was shown - the highest number characterized in a wild bird species. Finally, the functional alleles were grouped into 17 supertypes based on their antigen binding affinities.

CONCLUSIONS

We found extreme complexity at the Mhc class I of the great tit both in terms of allelic diversity and gene number. The presence of many functional loci was shown, together with a pseudogene family and putatively non-functional alleles; there was clear evidence that functional alleles were under strong balancing selection. This study is the first step towards an in-depth analysis of this gene complex in this species, which will help understanding how parasite-mediated and sexual selection shape and maintain host genetic variation in nature. We believe that study systems like ours can make important contributions to the field of evolutionary biology and emphasize the necessity of integrating long-term field-based studies with detailed genetic analysis to unravel complex evolutionary processes.

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.

MHC, mate choice and heterozygote advantage in a wild social primate

Preferences for mates carrying dissimilar genes at the major histocompatibility complex (MHC) may help animals increase offspring pathogen resistance or avoid inbreeding. Such preferences have been reported across a range of vertebrates, but have rarely been investigated in social species other than humans. We investigated mate choice and MHC dynamics in wild baboons (Papio ursinus). MHC Class II DRB genes and 16 microsatellite loci were genotyped across six groups (199 individuals). Based on the survey of a key segment of the gene-rich MHC, we found no evidence of mate choice for MHC dissimilarity, diversity or rare MHC genotypes. First, MHC dissimilarity did not differ from random expectation either between parents of the same offspring or between immigrant males and females from the same troop. Second, female reproductive success was not influenced by MHC diversity or genotype frequency. Third, population genetic structure analysis revealed equally high genotypic differentiation among troops, and comparable excess heterozygosity within troops for juveniles, at both Mhc-DRB and neutral loci. Nevertheless, the age structure of Mhc-DRB heterozygosity suggested higher longevity for heterozygotes, which should favour preferences for MHC dissimilarity. We propose that high levels of within-group outbreeding, resulting from group-living and sex-biased dispersal, might weaken selection for MHC-disassortative mate choice.

A female signal reflects MHC genotype in a social primate

BACKGROUND

Males from many species are believed to advertise their genetic quality through striking ornaments that attract mates. Yet the connections between signal expression, body condition and the genes associated with individual quality are rarely elucidated. This is particularly problematic for the signals of females in species with conventional sex roles, whose evolutionary significance has received little attention and is poorly understood. Here we explore these questions in the sexual swellings of female primates, which are among the most conspicuous of mammalian sexual signals and highly variable in size, shape and colour. We investigated the relationships between two components of sexual swellings (size and shape), body condition, and genes of the Major Histocompatibility Complex (MHC) in a wild baboon population (Papio ursinus) where males prefer large swellings.

RESULTS

Although there was no effect of MHC diversity on the sexual swelling components, one specific MHC supertype (S1) was associated with poor body condition together with swellings of small size and a particular shape. The variation in swelling characteristics linked with the possession of supertype S1 appeared to be partially mediated by body condition and remained detectable when taking into account the possession of other supertypes.

CONCLUSIONS

These findings suggest a pathway from immunity genes to sexual signals via physical condition for the first time in females. They further indicate that mechanisms of sexual selection traditionally assigned to males can also operate in females.

A prime-boost immunisation regimen using recombinant BCG and Pr55(gag) virus-like particle vaccines based on HIV type 1 subtype C successfully elicits Gag-specific responses in baboons

Mycobacterium bovis BCG is considered an attractive live bacterial vaccine vector. In this study, we investigated the immune response of baboons to a primary vaccination with recombinant BCG (rBCG) constructs expressing the gag gene from a South African HIV-1 subtype C isolate, and a boost with HIV-1 subtype C Pr55(gag) virus-like particles (Gag VLPs). Using an interferon enzyme-linked immunospot assay, we show that although these rBCG induced only a weak or an undetectable HIV-1 Gag-specific response on their own, they efficiently primed for a Gag VLP boost, which strengthened and broadened the immune responses. These responses were predominantly CD8+ T cell-mediated and recognised similar epitopes as those targeted by humans with early HIV-1 subtype C infection. In addition, a Gag-specific humoral response was elicited. These data support the development of HIV-1 vaccines based on rBCG and Pr55(gag) VLPs.