High levels of diversity characterize mandrill (Mandrillus sphinx) Mhc-DRB sequences

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.

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.

Testing for post-copulatory selection for major histocompatibility complex genotype in a semi-free-ranging primate population

A large body of evidence suggests that major histocompatibility complex (MHC) genotype influences mate choice. However, few studies have investigated MHC-mediated post-copulatory mate choice under natural, or even semi-natural, conditions. We set out to explore this question in a large semi-free-ranging population of mandrills (Mandrillus sphinx) using MHC-DRB genotypes for 127 parent-offspring triads. First, we showed that offspring MHC heterozygosity correlates positively with parental MHC dissimilarity suggesting that mating among MHC dissimilar mates is efficient in increasing offspring MHC diversity. Second, we compared the haplotypes of the parental dyad with those of the offspring to test whether post-copulatory sexual selection favored offspring with two different MHC haplotypes, more diverse gamete combinations, or greater within-haplotype diversity. Limited statistical power meant that we could only detect medium or large effect sizes. Nevertheless, we found no evidence for selection for heterozygous offspring when parents share a haplotype (large effect size), genetic dissimilarity between parental haplotypes (we could detect an odds ratio of ≥1.86), or within-haplotype diversity (medium-large effect). These findings suggest that comparing parental and offspring haplotypes may be a useful approach to test for post-copulatory selection when matings cannot be observed, as is the case in many study systems. However, it will be extremely difficult to determine conclusively whether post-copulatory selection mechanisms for MHC genotype exist, particularly if the effect sizes are small, due to the difficulty in obtaining a sufficiently large sample.

Odour signals major histocompatibility complex genotype in an Old World monkey

The major histocompatibility complex (MHC) is an extraordinarily diverse cluster of genes that play a key role in the immune system. MHC gene products are also found in various body secretions, leading to the suggestion that MHC genotypes are linked to unique individual odourtypes that animals use to assess the suitability of other individuals as potential mates or social partners. We investigated the relationship between chemical odour profiles and genotype in a large, naturally reproducing population of mandrills, using gas chromatography-mass spectrometry and MHC genotyping. Odour profiles were not linked to the possession of particular MHC supertypes. Sex influenced some measures of odour diversity and dominance rank influenced some measures of odour diversity in males, but not in females. Odour similarity was strongly related to similarity at the MHC, and, in some cases, to pedigree relatedness. Our results suggest that odour provides both a cue of individual genetic quality and information against which the receiver can compare its own genotype to assess genetic similarity. These findings provide a potential mechanism underlying mate choice for genetic diversity and MHC similarity as well as kin selection.

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.

Opposites attract: MHC-associated mate choice in a polygynous primate

We investigated reproduction in a semi-free-ranging population of a polygynous primate, the mandrill, in relation to genetic relatedness and male genetic characteristics, using neutral microsatellite and major histocompatibility complex (MHC) genotyping. We compared genetic dissimilarity to the mother and genetic characteristics of the sire with all other potential sires present at the conception of each offspring (193 offspring for microsatellite genetics, 180 for MHC). The probability that a given male sired increased as pedigree relatedness with the mother decreased, and overall genetic dissimilarity and MHC dissimilarity with the mother increased. Reproductive success also increased with male microsatellite heterozygosity and MHC diversity. These effects were apparent despite the strong influence of dominance rank on male reproductive success. The closed nature of our study population is comparable to human populations for which MHC-associated mate choice has been reported, suggesting that such mate choice may be especially important in relatively isolated populations with little migration to introduce genetic variation.

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

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

Correlated evolution of nucleotide substitution rates and allelic variation in Mhc-DRB lineages of primates

BACKGROUND

The major histocompatibility complex (MHC) is a key model of genetic polymorphism. Selection pressure by pathogens or other microevolutionary forces may result in a high rate of non-synonymous substitutions at the codons specifying the contact residues of the antigen binding sites (ABS), and the maintenance of extreme MHC allelic variation at the population/species level. Therefore, selection forces favouring MHC variability for any reason should cause a correlated evolution between substitution rates and allelic polymorphism. To investigate this prediction, we characterised nucleotide substitution rates and allelic polymorphism (i.e. the number of alleles detected in relation to the number of animals screened) of several Mhc class II DRB lineages in 46 primate species, and tested for a correlation between them.

RESULTS

First, we demonstrate that species-specific and lineage-specific evolutionary constraints favour species- and lineage-dependent substitution rate at the codons specifying the ABS contact residues (i.e. certain species and lineages can be characterised by high substitution rate, while others have low rate). Second, we show that although the degree of the non-synonymous substitution rate at the ABS contact residues was systematically higher than the degree of the synonymous substitution rate, these estimates were strongly correlated when we controlled for species-specific and lineage-specific effects, and also for the fact that different studies relied on different sample size. Such relationships between substitution rates of different types could even be extended to the non-contact residues of the molecule. Third, we provide statistical evidence that increased substitution rate along a MHC gene may lead to allelic variation, as a high substitution rate can be observed in those lineages in which many alleles are maintained. Fourth, we show that the detected patterns were independent of phylogenetic constraints. When we used phylogenetic models that control for similarity between species, due to common descent, and focused on variations within a single lineage (DRB1*03), the positive relationship between different substitution rates and allelic polymorphisms was still robust. Finally, we found the same effects to emerge in the analyses that eliminated within-species variation in MHC traits by using strictly single population-level studies. However, in a set of contrasting analyses, in which we focused on the non-functional DRB6 locus, the correlation between substitution rates and allelic variation was not prevalent.

CONCLUSION

Our results indicate that positive selection for the generation of allelic polymorphism acting on the functional part of the protein has consequences for the nucleotide substitution rate along the whole exon 2 sequence of the Mhc-DRB gene. Additionally, we proved that an increased substitution rate can promote allelic variation within lineages. Consequently, the evolution of different characteristics of genetic polymorphism is not independent.

Comparative genetics of a highly divergent DRB microsatellite in different macaque species

The DRB region of the major histocompatibility complex (MHC) of cynomolgus and rhesus macaques is highly plastic, and extensive copy number variation together with allelic polymorphism makes it a challenging enterprise to design a typing protocol. All intact DRB genes in cynomolgus monkeys (Mafa) appear to possess a compound microsatellite, DRB-STR, in intron 2, which displays extensive length polymorphism. Therefore, this STR was studied in a large panel of animals, comprising pedigreed families as well. Sequencing analysis resulted in the detection of 60 Mafa-DRB exon 2 sequences that were unambiguously linked to the corresponding microsatellite. Its length is often allele specific and follows Mendelian segregation. In cynomolgus and rhesus macaques, the nucleotide composition of the DRB-STR is in concordance with the phylogeny of exon 2 sequences. As in humans and rhesus monkeys, this protocol detects specific combinations of different DRB-STR lengths that are unique for each haplotype. In the present panel, 22 Mafa-DRB region configurations could be defined, which exceeds the number detected in a comparable cohort of Indian rhesus macaques. The results suggest that, in cynomolgus monkeys, even more frequently than in rhesus macaques, new haplotypes are generated by recombination-like events. Although both macaque species are known to share several identical DRB exon 2 sequences, the lengths of the corresponding microsatellites often differ. Thus, this method allows not only fast and accurate DRB haplotyping but may also permit discrimination between highly related macaque species.

Preservation of a pseudogene by gene conversion and diversifying selection

Interlocus gene conversion is considered a crucial mechanism for generating novel combinations of polymorphisms in duplicated genes. The importance of gene conversion between duplicated genes has been recognized in the major histocompatibility complex and self-incompatibility genes, which are likely subject to diversifying selection. To theoretically understand the potential role of gene conversion in such situations, forward simulations are performed in various two-locus models. The results show that gene conversion could significantly increase the number of haplotypes when diversifying selection works on both loci. We find that the tract length of gene conversion is an important factor to determine the efficacy of gene conversion: shorter tract lengths can more effectively generate novel haplotypes given the gene conversion rate per site is the same. Similar results are also obtained when one of the duplicated genes is assumed to be a pseudogene. It is suggested that a duplicated gene, even after being silenced, will contribute to increasing the variability in the other locus through gene conversion. Consequently, the fixation probability and longevity of duplicated genes increase under the presence of gene conversion. On the basis of these findings, we propose a new scenario for the preservation of a duplicated gene: when the original donor gene is under diversifying selection, a duplicated copy can be preserved by gene conversion even after it is pseudogenized.

Sequence variability at three MHC loci of finless porpoises (Neophocaena phocaenoides)

Major histocompatibility complex (MHC) class II DQB and DRA genes and class I gene of finless porpoises (Neophocaena phocaenoides) were investigated by single-strand conformation polymorphism and sequence analysis. The DRA, DQB, and MHC-I loci each contained 5, 14, and 34 unique sequences, respectively, and considerable sequence variation was found at the MHC-I and DQB loci. Gene duplication was manifested as three to five distinct sequences at each of the DQB and MHC-I loci from some individuals, and these sequences at each of the two loci separately clustered into four groups (cluster A, B, C, and D) based on the phylogenetic trees. Phylogenetic reconstruction revealed a trans-species pattern of evolution. Relatively high rates of non-synonymous (dN) vs synonymous (dS) substitution in the peptide-binding region (PBR) suggested balancing selection for maintaining polymorphisms at the MHC-I and DQB loci. In contrast, one single locus with little sequence variation was detected in the DRA gene, and no non-synonymous substitutions in the PBR indicated no balancing selection on this gene.

MHC class I A region diversity and polymorphism in macaque species

The HLA-A locus represents a single copy gene that displays abundant allelic polymorphism in the human population, whereas, in contrast, a nonhuman primate species such as the rhesus macaque (Macaca mulatta) possesses multiple HLA-A-like (Mamu-A) genes, which parade varying degrees of polymorphism. The number and combination of transcribed Mamu-A genes present per chromosome display diversity in a population of Indian animals. At present, it is not clearly understood whether these different A region configurations are evolutionarily stable entities. To shed light on this issue, rhesus macaques from a Chinese population and a panel of cynomolgus monkeys (Macaca fascicularis) were screened for various A region-linked variations. Comparisons demonstrated that most A region configurations are old entities predating macaque speciation, whereas most allelic variation (>95%) is of more recent origin. The latter situation contrasts the observations of the major histocompatibility complex class II genes in rhesus and cynomolgus macaques, which share a high number of identical alleles (>30%) as defined by exon 2 sequencing.