Factors shaping genetic variation in the MHC of natural non-human primate populations

The Cynomolgus Macaque MHC Polymorphism in Experimental Medicine

Among the non-human primates used in experimental medicine, cynomolgus macaques (Macaca fascicularis hereafter referred to as Mafa) are increasingly selected for the ease with which they are maintained and bred in captivity. Macaques belong to Old World monkeys and are phylogenetically much closer to humans than rodents, which are still the most frequently used animal model. Our understanding of the Mafa genome has progressed rapidly in recent years and has greatly benefited from the latest technical advances in molecular genetics. Cynomolgus macaques are widespread in Southeast Asia and numerous studies have shown a distinct genetic differentiation of continental and island populations. The major histocompatibility complex of cynomolgus macaque (Mafa MHC) is organized in the same way as that of human, but it differs from the latter by its high degree of classical class I gene duplication. Human polymorphic MHC regions play a pivotal role in allograft transplantation and have been associated with more than 100 diseases and/or phenotypes. The Mafa MHC polymorphism similarly plays a crucial role in experimental allografts of organs and stem cells. Experimental results show that the Mafa MHC class I and II regions influence the ability to mount an immune response against infectious pathogens and vaccines. MHC also affects cynomolgus macaque reproduction and impacts on numerous biological parameters. This review describes the Mafa MHC polymorphism and the methods currently used to characterize it. We discuss some of the major areas of experimental medicine where an effect induced by MHC polymorphism has been demonstrated.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Differential recombination dynamics within the MHC of macaque species

A panel of 15 carefully selected microsatellites (short tandem repeats, STRs) has allowed us to study segregation and haplotype stability in various macaque species. The STRs span the major histocompatibility complex (MHC) region and map in more detail from the centromeric part of the Mhc-A to the DR region. Two large panels of Indian rhesus and Indonesian/Indochinese cynomolgus macaques have been subjected to pedigree analysis, allowing the definition of 161 and 36 different haplotypes and the physical mapping of 10 and 5 recombination sites, respectively. Although most recombination sites within the studied section of the Indian rhesus monkey MHC are situated between the Mhc-A and Mhc-B regions, the resulting recombination rate for this genomic segment is low and similar to that in humans. In contrast, in Indonesian/Indochinese macaques, two recombination sites, which appear to be absent in rhesus macaques, map between the class III and II regions. As a result, the mean recombination frequency of the core MHC, Mhc-A to class II, is higher in Indonesian/Indochinese cynomolgus than in Indian rhesus macaques, but as such is comparable to that in humans. The present communication demonstrates that the dynamics of recombination 'hot/cold spots' in the MHC, as well as their frequencies, may differ substantially between highly related macaque species.

Study of MHC class II region polymorphism in the Filipino cynomolgus macaque population

The cynomolgus macaque (Macaca fascicularis) is currently used as an animal model in various fields of immunology especially in the development of innovative vaccines for the prevention and treatment of infectious diseases. The polymorphism of the major histocompatibility complex (MHC) influences the development of adaptive immune responses and it is crucial to characterize the polymorphism of cynomolgus MHC genes. We present here a systematic study of the MHC class II haplotypes in the Filipino macaque population. By the study of a large sample of Filipino animals (N = 353), we have characterized 18 MHC class II haplotypes by means of genotyping seven microsatellites. The animals were DRB genotyped by means of PCR-SSO or DGGE-sequencing on genomic amplified fragments. We cloned and sequenced the complementary DNA (cDNA) of DQA, DQB, DPA, and DPB genes of 117 animals. Combining the microsatellite genotyping and cDNA characterized in the 117 animals, we defined genetic association between the cDNA and the microsatellites and characterized 18 MHC class II haplotypes. For 104 animals out of the 353 studied, the presence of a recombinant haplotype was highly probable. Thirty-four percent of recombination was located in 256 kb segment between D6S2876 and D6S2747 microsatellites, a region encompassing several hot spots of recombination in the human MHC.

In vivo imaging in NHP models of malaria: challenges, progress and outlooks

Animal models of malaria, mainly mice, have made a large contribution to our knowledge of host-pathogen interactions and immune responses, and to drug and vaccine design. Non-human primate (NHP) models for malaria are admittedly under-used, although they are probably closer models than mice for human malaria; in particular, NHP models allow the use of human pathogens (Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae and Plasmodium knowlesi). NHPs, whether natural hosts or experimentally challenged with a simian Plasmodium, can also serve as robust pre-clinical models. Some simian parasites are closely related to a human counterpart, with which they may share a common ancestor, and display similar major features with the human infection and pathology. NHP models allow longitudinal studies, from the early events following sporozoite inoculation to the later events, including analysis of organs and tissues, particularly liver, spleen, brain and bone marrow. NHP models have one other significant advantage over mouse models: NHPs are our closest relatives and thus their biology is very similar to ours. Recently developed in vivo imaging tools have provided insight into malaria parasite infection and disease in mouse models. One advantage of these tools is that they limit the need for invasive procedures, such as tissue biopsies. Many such technologies are now available for NHP studies and provide new opportunities for elucidating host/parasite interactions. The aim of this review is to bring the malaria community up to date on what is currently possible and what soon will be, in terms of in vivo imaging in NHP models of malaria, to consider the pros and the cons of the various techniques, and to identify challenges.

Differential response of the cynomolgus macaque gut microbiota to Shigella infection

Little is known about the role of gut microbiota in response to live oral vaccines against enteric pathogens. We examined the effect of immunization with an oral live-attenuated Shigella dysenteriae 1 vaccine and challenge with wild-type S. dysenteriae 1 on the fecal microbiota of cynomolgus macaques using 16 S rRNA analysis of fecal samples. Multi-dimensional cluster analysis identified different bacterial community types within macaques from geographically distinct locations. The fecal microbiota of Mauritian macaques, observed to be genetically distinct, harbored a high-diversity community and responded differently to Shigella immunization, as well as challenge compared to the microbiota in non-Mauritian macaques. While both macaque populations exhibited anti-Shigella antibody responses, clinical shigellosis was observed only among non-Mauritian macaques. These studies highlight the importance of further investigation into the possible protective role of the microbiota against enteric pathogens and consideration of host genetic backgrounds in conducting vaccine studies.

Positive selection in the major histocompatibility complex class III region of cynomolgus macaques (Macaca fascicularis) of the Philippines origin

The cynomolgus macaque (Macaca fascicularis) is a model of choice among primates for the study of local adaptation processes because of its mixed and wide insular and continental distribution. In a previous study, by using 12 markers [5 microsatellites located in the major histocompatibility complex (MHC) region and 7 outside MHC], we have detected a signal of positive selection on the microsatellite DRACA located inside the gene DRA. In order to refine the location of this signal of positive selection in the MHC region, we studied the genetic diversity of 36 markers (18 microsatellites spread across the MHC region and 18 autosomal microsatellites outside MHC) in a sample of 254 individuals from four populations (Vietnam, Java, the Philippines, and Mauritius). We estimated for each locus the deviation of F(st) from a neutral model by using two methods based on contrasted demographic models. The two approaches showed a signal of positive selection in the MHC class III region that is much more significant than the one previously reported for the marker DRACA which could have been influenced by a hitchhiking effect due to its proximity with the class III region.

Mitochondrial DNA diversity among three subpopulations of cynomolgus macaques (Macaca fascicularis) originating from the Indochinese region

The cynomolgus macaque (Macaca fascicularis) has emerged as an important experimental animal model for biomedical research in various domains, necessitating the more extensive characterization of the genetic backgrounds influencing the macaque's response to drugs and sensitivity to experimental disease. The diversity of the variable mitochondrial DNA (mtDNA) D-loop region has been analyzed phylogenetically among geographically isolated populations or within subdivisions of the same regional population. However, the genetic differences among several substructures originating from a common population have not yet been investigated. By sequencing fragments of the mtDNA D-loop region from two subpopulations from the Indochinese region (Cambodian-Chinese and Vietnamese) along with two native Indonesian and Filipino populations, we identified 87 mtDNA D-loop haplotypes, of which 67 are new. The phylogenetic relationship suggests that the Indochinese haplotypes are intermingled in comparison to the distinct divergence of the Indonesian and Filipino lineages. The subpopulations were shown by estimation of evolutionary divergence and Wright's F-statistic (Fst) to have little genetic differentiation. Altogether, the subpopulations may be used in biomedical research, even though a slight difference is observed in haplotype frequencies among them. Therefore, genetic diversity analyses will be necessary for the elucidation of genetic differences among the populations, as well as to obtain a better understanding of genetic diversity for biomedical research. This will involve the selection of macaques and the monitoring of genetic heterogeneity among and within breeding facilities.

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.

Impact of MHC class II polymorphism on blood counts of CD4+ T lymphocytes in macaque

While the number of peripheral blood T lymphocytes and of their two main subsets (CD4+CD8- and CD4-CD8+) varies little in a given healthy individual, substantial variation is observed between individuals. It was proposed that these counts could be influenced by MHC polymorphisms because of the well-established role of MHC molecules in thymic T lymphocyte maturation and presentation of antigenic peptides to peripheral T lymphocytes. To test this hypothesis, we have chosen the crab-eating macaque (Macaca fascicularis), an animal model phylogenetically close to man. We selected the Philippine macaque population because of a restriction of the MHC polymorphism in this islander population. Peripheral blood lymphocytes were counted with an automated analyzer and T lymphocyte subsets were assessed by immunolabeling and flow cytometry. The MHC polymorphism was investigated in 200 unrelated subjects using 14 microsatellites markers distributed across the MHC and the DRB locus that was genotyped by denaturing gradient gel electrophoresis and sequencing. All markers were in Hardy-Weinberg equilibrium. Allelic associations were tested with the UNPHASED software. We revealed a significant influence of the MHC class II region on CD4+ T lymphocyte blood count with the largest effect associated with a two-locus haplotypes combining the DRACA allele 274 and the DRB haplotype #8a (p < 8 × 10(-7)). Our data should stimulate a similar association study of the CD4+ T cell counts in humans.

Erratum to: rhesus macaque MHC class I molecules show differential subcellular localizations

The MHC class I gene family of rhesus macaques is characterised by considerable gene duplications. While a HLA-C-orthologous gene is absent, the Mamu-A and in particular the Mamu-B genes have expanded, giving rise to plastic haplotypes with differential gene content. Although some of the rhesus macaque MHC class I genes are known to be associated with susceptibility/resistance to infectious diseases, the functional significance of duplicated Mamu-A and Mamu-B genes and the expression pattern of their encoded proteins are largely unknown. Here, we present data of the subcellular localization of AcGFP-tagged Mamu-A and Mamu-B molecules. We found strong cell surface and low intracellular expression for Mamu-A1, Mamu-A2 and Mamu-A3-encoded molecules as well as for Mamu-B*01704, Mamu-B*02101, Mamu-B*04801, Mamu-B*06002 and Mamu-B*13401. In contrast, weak cell surface and strong intracellular expression was seen for Mamu-A4*1403, Mamu-B*01202, Mamu-B*02804, Mamu-B*03002, Mamu-B*05704, Mamu-I*010201 and Mamu-I*0121. The different expression patterns were assigned to the antigen-binding alpha1 and alpha2 domains, suggesting failure of peptide binding is responsible for retaining 'intracellular' Mamu class I molecules in the endoplasmic reticulum. These findings indicate a diverse functional role of the duplicated rhesus macaque MHC class I genes.

Rhesus macaque MHC class I molecules show differential subcellular localizations

The MHC class I gene family of rhesus macaques is characterised by considerable gene duplications. While a HLA-C-orthologous gene is absent, the Mamu-A and in particular the Mamu-B genes have expanded, giving rise to plastic haplotypes with differential gene content. Although some of the rhesus macaque MHC class I genes are known to be associated with susceptibility/resistance to infectious diseases, the functional significance of duplicated Mamu-A and Mamu-B genes and the expression pattern of their encoded proteins are largely unknown. Here, we present data of the subcellular localization of AcGFP-tagged Mamu-A and Mamu-B molecules. We found strong cell surface and low intracellular expression for Mamu-A1, Mamu-A2 and Mamu-A3-encoded molecules as well as for Mamu-B*01704, Mamu-B*02101, Mamu-B*04801, Mamu-B*06002 and Mamu-B*13401. In contrast, weak cell surface and strong intracellular expression was seen for Mamu-A4*1403, Mamu-B*01202, Mamu-B*02804, Mamu-B*03002, Mamu-B*05704, Mamu-I*010201 and Mamu-I*0121. The different expression patterns were assigned to the antigen-binding α1 and α2 domains, suggesting failure of peptide binding is responsible for retaining ‘intracellular’ Mamu class I molecules in the endoplasmic reticulum. These findings indicate a diverse functional role of the duplicated rhesus macaque MHC class I genes.

Study of cynomolgus monkey (Macaca fascicularis) DRA polymorphism in four populations

To describe the polymorphism of the DRA gene in Macaca fascicularis, we have studied 141 animals either at cDNA level (78 animals from Mauritius, the Philippines, and Vietnam) or genomic level (63 animals from the Philippines, Indonesia, and Vietnam). In total, we characterized 22 cDNA DRA alleles, 13 of which had not been described until now. In the Mauritius population, we confirmed the presence of three DRA alleles. In the Philippine and Vietnam populations, we observed 11 and 14 DRA alleles, respectively. Only two alleles were present in all three populations. All DRA alleles but one differ from the consensus sequence by one to three mutations, most being synonymous; so, only seven DR alpha proteins were deduced from the 22 cDNA alleles. One DRA cDNA allele, Mafa-DRA*02010101, differs from all other alleles by 11 to 14 mutations of which only four are non-synonymous. The two amino acid changes inside the peptide groove of Mafa-DRA*02010101 are highly conservative. The very low proportion of non-synonymous/synonymous mutations is compatible with a purifying selection which is comparable to all previous observations concerning the evolution of the DRA gene in mammals. Homologues of the allele Mafa-DRA*02010101 are also found in two other Asian macaques (Macaca mulatta and Macaca nemestrina). The forces able to maintain this highly divergent allele in three different macaque species remain hypothetical.

MHC class I A loci polymorphism and diversity in three Southeast Asian populations of cynomolgus macaque

Cynomolgus macaques (Macaca fascicularis, Mafa) have emerged as important animal models for biomedical research, necessitating a more extensive characterization of their major histocompatibility complex polymorphic regions. The current information on the polymorphism or diversity of the polygenetic Mafa class I A loci is limited in comparison to the more commonly studied rhesus macaque Mafa class I A loci. Therefore, in this paper, to better elucidate the degree and types of polymorphisms and genetic differences of Mafa-A1 among three native Southeast Asian populations (Indonesian, Vietnamese, and Filipino) and to investigate how the allele differences between macaques and humans might have evolved to affect their respective immune responses, we identified 83 Mafa-A loci-derived alleles by DNA sequencing of which 66 are newly described. Most alleles are unique to each population, but seven of the most frequent alleles were identical in sequence to some alleles in other macaque species. We also revealed (1) the large and dynamic genetic and structural differences and similarities in allelic variation by analyzing the population allele frequencies, Hardy-Weinberg's equilibrium, heterozygosity, nucleotide diversity profiles, and phylogeny, (2) the difference in genetic structure of populations by Wright's FST statistic and hierarchical analysis of molecular variance, and (3) the different demographic and selection pressures on the three populations by performing Tajima's D test of neutrality. The large level of diversity and polymorphism at the Mafa-A1 was less evident in the Filipino than in the Vietnam or the Indonesian populations, which may have important implications in animal capture, selection, and breeding for medical research.

Characterization of 47 MHC class I sequences in Filipino cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) provide increasingly common models for infectious disease research. Several geographically distinct populations of these macaques from Southeast Asia and the Indian Ocean island of Mauritius are available for pathogenesis studies. Though host genetics may profoundly impact results of such studies, similarities and differences between populations are often overlooked. In this study we identified 47 full-length MHC class I nucleotide sequences in 16 cynomolgus macaques of Filipino origin. The majority of MHC class I sequences characterized (39 of 47) were unique to this regional population. However, we discovered eight sequences with perfect identity and six sequences with close similarity to previously defined MHC class I sequences from other macaque populations. We identified two ancestral MHC haplotypes that appear to be shared between Filipino and Mauritian cynomolgus macaques, notably a Mafa-B haplotype that has previously been shown to protect Mauritian cynomolgus macaques against challenge with a simian/human immunodeficiency virus, SHIV(89.6P). We also identified a Filipino cynomolgus macaque MHC class I sequence for which the predicted protein sequence differs from Mamu-B*17 by a single amino acid. This is important because Mamu-B*17 is strongly associated with protection against simian immunodeficiency virus (SIV) challenge in Indian rhesus macaques. These findings have implications for the evolutionary history of Filipino cynomolgus macaques as well as for the use of this model in SIV/SHIV research protocols.

Balancing and directional selection at exon-2 of the MHC DQB1 locus among populations of odontocete cetaceans

The diversity of exon-2 (peptide-binding region) of the DQB1 locus (Class II, major histocompatibility complex, MHC) was investigated on an extended sample of populations of three focal cetacean species (two sibling delphinid species and another in the same family). We tested the hypothesis that dolphin populations with a worldwide distribution across different habitats and geographic regions will be under differential selective pressure by comparing DQB1 variation with variation at neutral markers and by investigating putative functional residues within the exon-2 sequence at the population level. Variation at the DQB1 locus was not correlated to neutral differentiation (assessed by comparison with microsatellite DNA markers), and overall F(ST) values were significantly lower for the MHC locus, consistent with expectations for balancing selection. Measures of heterozygosity and d(n)/d(s) ratios were also consistent with balancing selection. However, outliers in the F(ST) comparisons and the analysis of putative functional residues suggested incidences of directional selection in local populations.

Assessing natural introgression in 2 biomedical model species, the rhesus macaque (Macaca mulatta) and the long-tailed macaque (Macaca fascicularis)

Rhesus macaque (Macaca mulatta) and long-tailed macaque (Macaca fascicularis) are the 2 most commonly used primate model species in biomedical sciences. Although morphological studies have revealed a weak hybridization at the interspecific contact zone, in the north of Indochina, a molecular study has suggested an ancient introgression from rhesus to long-tailed macaque into the Indo-Chinese peninsula. However, the gene flow between these 2 taxa has never been quantified using genetic data and theoretical models. In this study, we have examined genetic variation within and between the parapatric Chinese rhesus macaque and Indo-Chinese long-tailed macaque populations, using 13 autosomal, 5 sex-linked microsatellite loci and mitochondrial DNA sequence data. From these data, we assessed genetic structure and estimated gene flow using a Bayesian clustering approach and the "Isolation with Migration" model. Our results reveal a weak interspecific genetic differentiation at both autosomal and sex-linked loci, suggesting large population sizes and/or gene flow between populations. According to the Bayesian clustering, Chinese rhesus macaque is a highly homogeneous gene pool that contributes strongly to the current Indo-Chinese long-tailed macaque genetic makeup, whether or not current admixture is assumed. Coalescent simulations, which integrated the characteristics of the loci, pointed out 1) a higher effective population size in rhesus macaque, 2) no mitochondrial gene flow, and 3) unilateral and male-mediated nuclear gene flow of approximately 10 migrants per generation from rhesus to long-tailed macaque. These patterns of genetic structure and gene flow suggest extensive ancient introgression from Chinese rhesus macaque into the Indo-Chinese long-tailed macaque population.

MHC class I characterization of Indonesian cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) are quickly becoming a useful model for infectious disease and transplantation research. Even though cynomolgus macaques from different geographic regions are used for these studies, there has been limited characterization of full-length major histocompatibility complex (MHC) class I immunogenetics of distinct geographic populations. Here, we identified 48 MHC class I cDNA nucleotide sequences in eleven Indonesian cynomolgus macaques, including 41 novel Mafa-A and Mafa-B sequences. We found seven MHC class I sequences in Indonesian macaques that were identical to MHC class I sequences identified in Malaysian or Mauritian macaques. Sharing of nucleotide sequences between these geographically distinct populations is also consistent with the hypothesis that Indonesia was a source of the Mauritian macaque population. In addition, we found that the Indonesian cDNA sequence Mafa-B7601 is identical throughout its peptide binding domain to Mamu-B03, an allele that has been associated with control of Simian immunodeficiency virus (SIV) viremia in Indian rhesus macaques. Overall, a better understanding of the MHC class I alleles present in Indonesian cynomolgus macaques improves their value as a model for disease research, and it better defines the biogeography of cynomolgus macaques throughout Southeast Asia.