Skip to main content

Genetic diversity and differentiation of the rhesus macaque (Macaca mulatta) population in western Sichuan, China, based on the second exon of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles



Rhesus macaques living in western Sichuan, China, have been separated into several isolated populations due to habitat fragmentation. Previous studies based on the neutral or nearly neutral markers (mitochondrial DNA or microsatellites) showed high levels of genetic diversity and moderate genetic differentiation in the Sichuan rhesus macaques. Variation at the major histocompatibility complex (MHC) loci is widely accepted as being maintained by balancing selection, even with a low level of neutral variability in some species. However, in small and isolated or bottlenecked populations, balancing selection may be overwhelmed by genetic drift. To estimate microevolutionary forces acting on the isolated rhesus macaque populations, we examined genetic variation at Mhc-DQB1 loci in 119 wild rhesus macaques from five geographically isolated populations in western Sichuan, China, and compared the levels of MHC variation and differentiation among populations with that previously observed at neutral microsatellite markers.


23 Mamu-DQB1 alleles were identified in 119 rhesus macaques in western Sichuan, China. These macaques exhibited relatively high levels of genetic diversity at Mamu-DQB1. The Hanyuan population presented the highest genetic variation, whereas the Heishui population was the lowest. Analysis of molecular variance (AMOVA) and pairwise FST values showed moderate genetic differentiation occurring among the five populations at the Mhc-DQB1 locus. Non-synonymous substitutions occurred at a higher frequency than synonymous substitutions in the peptide binding region. Levels of MHC variation within rhesus macaque populations are concordant with microsatellite variation. On the phylogenetic tree for the rhesus and crab-eating macaques, extensive allele or allelic lineage sharing is observed betweenthe two species.


Phylogenetic analyses confirm the apparent trans-species model of evolution of the Mhc-DQB1 genes in these macaques. Balancing selection plays an important role in sharing allelic lineages between species, but genetic drift may share balancing selection dominance to maintain MHC diversity. Great divergence at neutral or adaptive markers showed that moderate genetic differentiation had occurred in rhesus macaque populations in western Sichuan, China, due to the habitat fragmentation caused by long-term geographic barriers and human activity. The Heishui population should be paid more attention for its lowest level of genetic diversity and relatively great divergence from others.


How genetic variation is maintained within populations of endangered species is a central issue in evolutionary and conservation biology [1]. Diversity as a measure of individual variation within a population is widely thought to reflect the number of different genotypes in the population, taking into account their frequencies [2]. Neutral markers are often a measure of genetic diversity, but neutral genetic variations do not always correlate with variations at adaptively important genes [3]. Therefore, the importance of genetic diversity at markers of adaptive significance has been increasingly recognized [4, 5]. The multi-gene major histocompatibility complex (MHC) family is found in vertebrates, coding for cell surface glycoproteins and is important in animal conservation due to its role in resisting pathogens [6]. The classical MHC genes of rhesus macaques can be divided into MHC class I and II genes. The MHC class I genes include mainly -A and -B alleles, and MHC class II genes include mainly -DM, −DO, −DP, −DQ, and -DR alleles. MHC genes are well-known examples of genes of adaptive significance and are particularly relevant to conservation [7]. Consequently, MHC variability is a reflection of the processes that are related to adaptive evolution within and between populations [8], and the generation and maintenance of allelic polymorphisms in the MHC genes is a major topic in evolutionary genetics. The high levels of polymorphisms usually observed at MHC genes are most likely to be maintained by balancing selection, and driven largely by exposure to a diversity of pathogens [911], but some studies on MHC variation in natural populations show that balancing selection does not always maintain high levels of MHC variation relative to neutral markers [12, 13]. In particular in small populations, the strength of balancing selection on MHC genes may be weak relative to other microevolutionary forces such as genetic drift [6, 14]. Genetic drift is the random fluctuation of allele frequencies over time, thus, adaptive alleles may be lost, and deleterious alleles could be fixed in the population. Many natural populations are threatened by intense reduction and fragmentation of habitat, which leads to population isolation, decline in the number and the loss of genetic diversity [15]. The small population size and fixation of deleterious alleles leads to inbreeding depression and reduction of individual fitness, which decreases viability and compromises a population’s evolutionary adaptive potential [16]. The loss of genetic diversity may increase the risk of extinction due to decreased reproductive fitness, decreased adaptive flexibility, and increased disease susceptibility [17]. An important assumption in conservation genetics is that small and isolated populations are more sensitive to genetic drift and inbreeding [18, 19]. Clarifying the mechanism that determines genetic variation in small and isolated populations is therefore essential for species conservation [19].

Western Sichuan, located along the southeastern edge of the Qinghai-Tibetan Plateau (QTP) and the northern section of Hengduan Mountains, is the global biodiversity hotspots. Affected by the Pliocene uplifting of the QTP, western Sichuan formed a complex and diverse topography with high mountains and deep valleys. The region, from southwest to northeast, has four major mountains, namely, Liangshan, Xiangling, Qionglai, and Mingshan and four large rivers, namely, Jinshajiang, Yalongjiang, Daduhe, Mingjiang and Jialingjiang (Figure 1). Additionally, there are many human settlements and roads in the region. Thus, many wild animals and plants – endangered species in particular – have actually formed fragmented distributions in western Sichuan due to the barrier effect of these natural and man-made impediments. The rhesus macaque (Macaca mulatta) is mainly found throughout most of southern Asia crossing eastern Afghanistan, Bangladesh, Bhutan, northern and central India, central and southern China, Lao PDR, Myanmar, Nepal, northern Pakistan, northern Thailand, and Vietnam, and is a endangered species due to the sharp reduction in the number throughout its wider distribution in recent years. This species is listed in CITES Appendix II, Schedule III of the Bangladesh Wildlife (Preservation) Act (1974), Schedule I of Part I of the Indian Wildlife (Protection) Act (amended up to 2002), and Category II of the Chinese Wildlife Protection Act (1989). According to Jiang [20], Chinese rhesus macaques were divided into six subspecies. Rhesus macaques in western Sichuan are classified as a separate subspecies and designated as M. mulatta lasiotis, which was isolated from the other five subspecies by human activities and natural barriers (Yangtze River) and has been fragmentized in western Sichuan mainly due to the special topography in this region (Figure 1). Therefore, western Sichuan is a suitable place to explore how microevolutionary forces such as balancing selection and genetic drift maintain the variation pattern of rhesus macaque populations in western China. Understanding evolutionary processes that maintain genetic diversity in natural populations is an essential goal of population genetics [21, 22]. Neutral or nearly neutral markers such as mitochondrial or microsatellite DNA are informative for phylogenetic and phylogeographic reconstructions [23]. Genetic diversity analysis using mitochondrial DNA markers showed that significant genetic differentiation had occurred among different isolated populations of Sichuan rhesus macaques [24]. Subsequent microsatellite loci analysis also demonstrated that the genetic differentiation has occurred between a local (Heishui) population and other populations [25]. Another study, using a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, revealed two maternal super haplogroups that exist in Chinese rhesus macaques, one in eastern China and another in western China, and further analysis showed that orogenesis likely drove the divergence of western populations in China [26].

Figure 1

Distribution of sampling localitys in western Sichuan. Each filled circle representing a population and filled triangle representing Gongga Mountains (GGM). JL, Jiulong; HY, Hanyuan; DB, Danba; XJ, Xiaojin; HS, Heishui. The map in this figure was produced by author Qiuxia Dai.

The rhesus macaque, cynomolgus macaque, common marmoset and owl monkey are non-human primates that are often used as experimental model animals for biomedical research, such as transplantation studies and vaccine development against measles, Ebola, or other infectious diseases including acquired immune deficiency syndrome (AIDS) [27, 28]. For use in immune-related studies, it is necessary that the MHC background of these animals be characterized. Until now, various MHC genes including class I and II have been identified and characterized in different degrees for these non-human primate animals. For example, 14 Aona-DQB1 alleles were identified that separated the two lineages Aona-DQB1*22 and Aona-DQB1*23 in a group of 19 unrelated owl monkeys (Aotus nancymaae) [29]; Otting et al. [30] sequenced exon 2 of the Mhc-DQB gene in each of a random panel of 60 non-pedigreed cynomolgous macaques (Macaca fascicularis), resulting in the detection of 23 Mafa-DQB1 alleles that had not previously been published and confirmed the trans-species model of evolution of the Mhc-DQB lineages, in which a group of major allele lineages is shared by different species of non-human primates in the phylogeny. Qiu et al. [31] identified 37 MhcMamu-DQB1 alleles in 105 Chinese rhesus macaques by cloning and sequencing, illustrating a marked allelic polymorphism at DQB1 in these monkeys. Notably, in recent years more MHC class II DQB alleles have been identified and characterized in two closely related monkeys, rhesus and cynomolgous macaques, due to their common use as experimental animal models in biomedical research.

The above work provided us with the inspiration and base to explore the genetic variation patterns of rhesus macaques in western China based on MHC genes with adaptive significance because of their role in pathogen resistance. The aims of this study included the following: 1) to investigate polymorphism levels of MHC class II DQB1 locus in five isolated populations of rhesus macaques in western Sichuan, China; 2) to estimate the degree of genetic differentiation among or within populations; 3) to determine the relative roles of microevolutionary forces such as balancing selection or genetic drift in populations; and 4) to confirm the trans-species model of evolution of Mhc-DQB1 lineages between rhesus and cynomolgous macaques.


Frequency distribution of Mamu-DQB1alleles

A total of 119 rhesus macaque individuals in western Sichuan were typed into 49 genotypes through the polymerase chain reaction (PCR) and single-strand conformation polymorphism (SSCP) analysis based on the polymorphism of the second exon of the Mamu-DQB1 gene, and then cloning and sequencing of the second exon of the Mamu-DQB1 gene were conducted for the samples representing the 49 genotypes. Consequently, 23 Mamu-DQB1 alleles were identified in Sichuan rhesus macaques (Table 1). Each Mamu-DQB1 sequence obtained was 274 bp (excluding primers) in length and spanned the complete sequence of exon 2. No more than two alleles were detected in each individual, and no stop codon or insertion/deletion was observed in Mamu-DQB1. By blast searches in the IPD-MHC Database (, the 23 allelic sequences were found to be completely identical to some designated alleles, showing that these alleles also existed in other rhesus macaque individuals in previous study.

Table 1 Distribution of allelic frequencies for Mamu-DQB1 in five rhesus macaque populations

The number and frequency of the Mamu-DQB1 alleles varied substantially in five different geographical populations (Table 1 and Figure 2). The highest frequency was Mamu-DQB*0614 among these rhesus macaques, which was found in 31 (13.03%) of the 119 monkeys; the next most common frequency was Mamu-DQB*1503, which was detected in 30 (12.61%) of these monkeys. The third frequency was Mamu-DQB*1703, which was detected in 29 (24.37%) of these monkeys. Three of the five geographical populations had their unique alleles. Mamu-DQB*1502, Mamu-DQB*1601, Mamu-DQB*170601, Mamu-DQB*1804, and Mamu-DQB*1819 were found only in the Hanyuan population, Mamu-DQB*061301 and Mamu-DQB*1818 were found only in the Danba population, and Mamu-DQB*0607 was found only in the Jiulong population. Mamu-DQB*0614 was found in all populations at different frequency from 1.79% to 21.43%. Mamu-DQB*061301 was detected only in one individual.

Figure 2

Histogram of MhcMamu-DQB1 frequency distribution in populations. The y-axis represents the alliic frequency and x-axis represents 23 MhcMmamu-DQB1 alleles, No. 1–23 presents each allele name in Table 1; JL, Jiulong; HY, Hanyuan; DB, Danba; XJ, Xiaojin; HS, Heishui. Red asterisk (*) represents group-specific alleles in this population.

Genetic diversity of the rhesus macaque populations

In order to estimate genetic diversity of the rhesus macaques in western Sichuan based on the adaptive marker, several genetic diversity parameters such as haplotype diversity (h), nucleotide diversity (π), expected heterozygosity (He) and allelic richness (A R ) were calculated (Table 2). The values of h, He, and A R were 0.718 ~ 0.891, 0.780 ~ 0.875, and 5.733 ~ 11.938, respectively, indicating a relatively high level of genetic diversity at Mamu-DQB1 locus. The analysis of nucleotide diversity revealed that most of the nucleotide variation was located in the antigen-binding site (ABS) region, and there were at least 20 nucleotide variations on average in each population. The Hanyuan population exhibited the highest genetic diversity: 15 allele haplotypes were observed in 28 animal individuals. The Heishui population exhibited the lowest genetic diversity: only 7 allele haplotypes were observed in 30 animal individuals. Surprisingly, the inbreeding coefficient (F is ) in the Hanyuan population was also the highest, whereas it was lowest in Heishui population.

Table 2 Genetic diversity of Mamu-DQB1 within rhesus macaque populations

Phylogenetic analysis

The phylogenetic tree was constructed for 156 Mhc-DQB1 gene sequences, including 69 Mamu-DQB1, 85 Mafa-DQB1 and 2 Aona-DQB1 alleles by using the neighbor-joining method (Figure 3), and aona-DQB1 alleles were used as the out-group. The 156 sequences can be searched in the IPD-MHC Database ( As is evident from the tree, the majority of Mhc-DQB1 alleles tended to cluster together according to lineage group rather than species or sample origins, revealing a trans-species model of evolution. Twenty-three Mamu-DQB1 alleles were classed into five lineages: DQB1*06, *15, *16, *17, and *18, of which lineage DQB1*06 and *18 were most dominant (9/23 = 39% and 7/23 = 30%, respectively). Almost all previously reported major Mamu-DQB1 lineages were detected (DQB1*06, *15, *17, *18), except for DQB1*24 lineages that were not found in the western Sichuan macaques examined in this study. It is worth noting that rhesus macaques shared only 14 identical alleles with crab-eating macaques on the phylogenetic tree, but there were 10 alleles to be detected in this study (see Table 1).

Figure 3

Phylogenetic tree of the MhcMamu- and MhcMafa-DQB1 alleles using neighbour-joing methods. 23 Mamu-DQB1 sequences detected in this study are shown with a solid round spot. 46 Mamu-DQB1 alleles not detected in animals of this study are showed with a circle. 85 Mafa-DQB1 sequences are showed without anything. 2 Aona-DQB1 (Aona-DQB1*220101 and Aona-DQB1*220102) were selected as outgroup. Except for 23 Mamu-DQB1 sequences detected in this study, all other DQB1 sequences were retrieved from the IPD-MHC Database ( Numbers at the node are bootstrap value in statistic analysis, only values above 50% are showed.

Selective pressure analysis

To analyze the selective pressure on Mamu-DQB1 exon 2 sequences in western Sichuan rhesus populations, the average non-synonymous and synonymous substitution rates (d N /d S , ω) for all sites, antigen-binding sites (ABS), and non-antigen binding sites (non-ABS) were calculated (Table 3, Dataset A: 23 sequences). The highest ω ratio value was found in ABS and was greater than one (ω = 1.458), whereas the lowest ω ratio value was found in non-ABS and less than one (ω = 0.976), but these values did not significantly deviate from the neutral expectation in the Z-test (p > 0.05). To better understand the substitution pattern of Mamu-DQB1 exon 2 sequences, 46 Mamu-DQB1 exon 2 sequences from the previous studies were added to our sequence data to expand our analysis (Table 3, Dataset B: 65 sequences). The same analysis with Dataset A was performed. We found that all ω ratio values for dataset B were greater than that of Dataset A and greater than one, especially the ω ratio value in the ABS of Dataset B was found to be significantly greater than one (ω = 2.041, p < 0.05). These results indicated that the ABS of Mamu-DQB1 exon 2 sequences had been subjected to a strong balancing selection, but the selection signal was not significant in the Mamu-DQB1 exon 2 sequences of western Sichuan rhesus macaques.

Table 3 Mean numbers of nucleotide substitutions per non-synonymous site ( d N ) and synonymous site ( d S ) and the d N /d S ratios (ω)

Genetic differentiation among five rhesus-macaque populations

From the analysis of molecular variance (AMOVA, Table 4) for the five populations, 85.29% of the total variation in our dataset was observed within populations and 14.71% was observed among populations, showing that the genetic variation primarily came from within population. Although the genetic differentiation within a population was significantly higher than among populations, there was still a moderate differentiation among populations because of a larger fixation index (FST) (0.14712, p < 0.01). To further assess the genetic differentiation among populations, the pairwise comparison FST and gene flow (Nm) values among populations were calculated (Table 5). The degree of genetic differentiation was smallest between the Danba and Xiaojin populations, whereas the Heishui and Jiulong populations had the largest differentiation. The Heishui population had relatively little gene exchange with others.

Table 4 Analysis of molecular variance (AMOVA) of Mamu-DQB1 sequences in five rhesus macaque populations
Table 5 Genetic differentiation among five rhesus macaque populations

To further explore genetic differentiation, we compared He and A R between microsatellite and Mamu-DQB1 markers (Figure 4). The richness of MHC sequences per population provides a measure of allelic diversity within each population; the A R values of both the microsatellite sites and Mamu-DQB1 are significantly positively correlated (r2 = 0.961, df = 4, p = 0.003). However, there is not a significant correlation between microsatellite sites and Mamu-DQB1 heterozygosity per population (r2 = 0.2293, df = 4, p = 0.4145). No significant correlation to pairwise FST between microsatellite sites and Mamu-DQB1 (r2 = 0.2198, df = 9, p = 0.1716). The FST values of Mamu-DQB1 showed higher divergence among populations when compared with microsatellite loci (except for the Danba-Hanyuan and Danba-Xiaojin populations), and all of the FST values (except for that of the Danba-Xiaojin population) at the Mamu-DQB1 locus were higher than 0.05, indicating that divergence existed among populations, especially between the Heishui and Jiulong populations, which had the largest F ST . Overally, the population difference was more obvious at the Mamu-DQB1 locus than the microsatelites: six F ST values from the Mamu-DQB1 locus were far larger than that of the microsatellites in ten pariwise population comparisons (Table 5).

Figure 4

Correlations between MhcMamu-DQB1 allele and microsatellite locus variation. (A) Mean microsatellite heterozygosity vs. MHC heterozygosity. (B) Mean microsatellite allelic richness vs MHC allelic richness per population.


Genetic diversity of Mamu-MhcDQB1locus

In various studies, the maintenance of a high level of genetic diversity is assumed to be a preliminary strategy for the population conservation [3234]. More specifically, high genetic diversity has often been shown to be positively correlated with indicator values of individual fitness such as the mating success of animals and adaptability to changing environmental conditions [35, 36]. In accordance with a previous study in rhesus monkeys, Chinese-origin rhesus macaques have a greater genetic variation at the Mamu-DQB1 locus and higher nucleotide diversity of mtDNA when compared with Indian-origin rhesus macaques [31, 37]. For example, only 16 Mamu-DQB1 alleles were identified from 107 Indian origin individuals born in captivity [38], whereas 23 Mamu-DQB1 alleles were identified from 119 Chinese rhesus macaques in this study.

The high level of polymorphism usually observed in MHC genes is likely to be maintained by balancing selection, driven largely by exposure to a diversity of pathogens [911]. In the present study, the Hanyuan population had conspicuously higher A R , haplotype diversity (h) and He values than any other populations (Table 2), exhibiting higher genetic diversity. This implied a strong balancing selection force that may have occurred in this population. In addition, the Hanyuan population lives in the Xiangling Mountains which is a transition zone between the north and south of two distinct climates, the wet and cold Qionglai mountains in the north and the dry and hot Liangshan mountains in the south. Thus, the climate transition zone might also be a genetic intersection area of species, leading to higher genetic diversity than other nearby areas. However, the Heishui population exhibited a lower level of genetic diversity than other populations. In particular, Mamu-DQB1*0614 and Mamu-DQB1*1703 are significantly more frequent in this population (Table 1). It is implied that a founder effect or genetic drift may have occurred in Heishui population. In addition, previous analyses of MHC variation in species that span large geographical ranges have emphasized the role of diversifying selection on MHC genes [3942], thus we can not rule out the possibility of strong selection for the two alleles (Mamu-DQB1*0614 and Mamu-DQB1*1703) that are present at considerably higher frequencies in this population. Nevertheless, it is possible that there are additional MHC alleles that were not detected in Heishui population because of its smaller sample size.

Disassortative mate choice may also contribute to the maintenance of MHC diversity [43, 44]. Choosing partners with maximal or optimal genetic dissimilarity may help to avoid inbreeding or excessive outbreeding, respectively. To increase offspring diversity at key functional loci in a natural population, a recent study indicated that inbreeding avoidance was along with disassortative mate choice at DRB, but not at DQB locus [45]. Our results also seem to support this point. For example, the Hanyuan and Jiulong populations had high levels of genetic diversity, but they also had the highest inbreeding coefficients (Fis) at the Mamu-DQB1 locus (Table 2).

Balancing selection and trans-species polymorphism

Balancing selection is expected to preserve high levels of polymorphisms at MHC loci by retaining alleles during species diversification events [46, 47]. The d N /d S test for high rates of non-synonymous (d N ) compared to synonymous (d S ) substitutions (ω = d N /d S ) is the most common measure used to detect balancing selection acting on protein-coding genes. Polymorphisms within the ABS of MHC genes are thought to be maintained by balancing selection imposed by pathogens/parasites. Fast-evolving parasites may adapt to the most common host genotype and escape presentation of their antigens to the adaptive immune system of the host [48]. In this study, the analysis of the nucleotide acid substitution pattern showed that Mamu-DQB1 exon 2 sequences had been subjected to a strong balancing selection, but the selection signal was not significant in the Mamu-DQB1 exon 2 sequences of western Sichuan rhesus macaques. This may well be caused by genetic drift acting differently on the different populations, e.g., the Heishui population may has been subject to a genetic drift, leading to a low level of genetic diversity. Alternatively, this may also be due to the smaller sample size surveyed here.

Under balancing selection, some allelic lineages exhibit unusual longevity that predates speciation events, leading to “trans-species polymorphisms” (TSP) used to describe polymorphisms among species [49]. In an extreme case of TSP, different species share identical alleles. Identical MHC gene sequences can be found in different species that have been separated by millions of years of independent evolution [50]. This is remarkable, given that these genes are involved in a co-evolutionary arms race with parasites [51]. In the present study, the phylogenetic analysis has illustrated that Mhc-DQB1 allele sequences from rhesus and cynomolgus macaques tended to cluster together according to allelic lineages rather than species or sample origins, leading to the sharing of allelic lineages by different species. The result clearly confirmed the trans-species model of evolution of the Mhc-DQB1 lineages. In particular, 14 Mamu-DQB1 alleles matched exactly with the Mafa-DQB1 alleles, suggesting that they may have occurred because of trans-species conservation after the divergence of the rhesus and cynomolgus macaque 1.8–2.0 million years ago (mya) [52, 53, 28]. Alternatively, the identical nucleotide sequences in the rhesus and cynomolgus macaques may have arisen by occasional interbreeding between the two species. In fact, a recent study on genome sequencing and comparison of the cynomolgus and Chinese rhesus macaques showed that the cynomolgus macaque genome has been shaped by introgression after hybridization with the Chinese rhesus macaque [54]. The allele sharing between the rhesus and cynomolgus macaques has also been found in other MHC class II genes such as DQA1, DPB1, and DRB1 [38]. However, it is unusual that the sharing phenomenon has not been by now detected in any MHC class I genes, which is a synthetic group in the same chromosome with class II genes.

Genetic differentiation among populations

AMOVA of the five populations indicated that variations mainly occurred among populations and that the Fst was moderate (0.14712) (Table 4). Fst represents the level of genetic differentiation among populations [55, 56]: an Fst of 0 to 0.05 represents “little differentiation,” 0.05 to 0.25 “moderate differentiation,” and values greater than 0.25 “very great differentiation” [57]. Therefore, our results revealed that moderate genetic differentiation had occurred among the five rhesus macaque populations. To further assess genetic differentiation among the five populations, we also estimated Fst and Nm. Nm (also known as gene migration) is the transfer of genes from one population to another. Generally, if Nm < 1, genetic drift will result in substantial local differentiation; if Nm > 1, gene flow between populations is higher and the extent of genetic differentiation is smaller. When Nm > 4, gene exchange is more frequent and genetic differentiation is much smaller [58, 59]. In this study, all values of Nm were larger than 1 (1.533–6.175). Especially, the values of Nm among three populations (Danba, Xiaojin, and Hanyuan) were even larger than 4, suggesting that gene exchange was more frequent among the three populations.

One hallmark of animal mtDNA is strict maternal inheritance. Previous study about mtDNA indicated poor gene exchange occurred in these five populations with the Nm values less than 1, thus gene exchange may occur mainly in nuclear genes. In other words, genetic exchanges mainly come from the migration of paternal macaques, which is consistent with male-biased migration behavior in various macaques. However, the levels of differentiation among the five populations were moderate. The differentiation were highest between the Jiulong and Heishui populations, and lowest between the Danba and Xiaojin populations. We speculated that genetic differentiation between these populations has arisen due to the geographical barrier of the high mountains, large rivers and long geographic distances between populations. For example, Jiulong and Danba populations are located between the Yalongjiang and Dadu Rivers, but they are separated by the Gongga Mountains. It is difficult for rhesus macaques to cross these geographical barriers because of their vast width, depth, length, and flow. In fact, it has been shown that large rivers formed a natural barrier that obstructs gene exchange between many mammal populations [60, 61].

Differences in parasite abundance and diversity in environments may result in selection for different sets of alleles in different populations. Loci under balancing selection are predicted to own a higher variation within population, but lower differentiation between populations when compared with neutral or nearly neutral markers, assuming that the selective pressures in the two populations are similar [62]. However, genes under directional selection present increased divergence between populations relative to neutral or nearly neutral markers if selective pressures differ between populations, or they may present decreased divergence relative to neutral or nearly neutral markers if selection pressures are similar [63]. Thus, compared with neutral forces, balancing selection is supposed to diminish population differentiation as measured by conventional pairwise F ST . Hence, the population structure under balancing selection should not be distinct. In the present study, nine out of ten of the pairwise F ST values of Mamu-DQB1 were greater than 0.05 and eight pairwise FST values were even greater than that of the microsatellites (Table 2). A significant positive correlation was found between allelic richness in Mamu-DQB1 and microsatellites; a similar positive correlation was also found on pairwise F ST values, implying that genetic drift had played a significant role in maintaining MHC diversity for rhesus macaques. However, as mentioned above, the sharing of allelic lineages between species observed on our phylogenetic tree indicated a “trans-species polymorphism” under balancing selection, but the rhesus and cynomolgus macaques have diverged from each only 1.8–2.0 mya [53]. Therefore, the balancing selection acting on Mamu-DQB1 may be historical.


It is widely accepted that variation at MHC loci is maintained by balancing selection, even with a low level of neutral variability in some species [64, 65], indicating the importance of balancing selection for maintaining genetic variation in the wild and also exposing the limitations of the neutral marker as substitutes for variation in fitness-related genes [66]. However, balancing selection may be overwhelmed by other microevolutionary forces such as genetic drift in small and isolated or bottlenecked populations [67]. In this study, we detected evidence of a weaker balancing selection and found similar genetic variation pattern both at neutral and adaptive markers, suggesting that genetic drift was stronger than selection, thus leading to a reduction of MHC diversity in the local populations. AMOVA results also revealed that a moderate genetic differentiation had occurred in several segregated populations in western Sichuan, China, due to habitat fragmentation caused by long-term geographic barriers and human activities. The Heishui population should be paid more attention for its lowest level of genetic diversity and relatively high divergence from others. This study provides deep insight into genetic variation of rhesus macaque populations in the fragmented habitat and helps to make an effective strategy for species conservation and facilitate the use of rhesus macaque in biomedical research. However, we have to acknowledge that the findings in this study are still very limited because of the smaller sample size as well as a single MHC gene locus. Therefore, to clarify thoroughly the genetic variation pattern of rhesus macaque populations in western China, further research is still needed using other MHC genes, other immune system loci and other markers, as well as more comprehensive sampling.


Samples and Genomic DNA extraction

The whole blood samples were collected from 119 western Sichuan rhesus macaques at the National Experimental Rhesus Research Centre of Sichuan Agricultural University. All of these animals were the first-captured founder population from five different geographical populations in the wild, Jiulong (n = 20), Hanyuan (n = 28), Danba (n = 14), Xiaojin (n = 27), and Heishui (n = 30), located in western Sichuan, China (Figure 1). For each population, the samples are selected randomly to try to ensure that no correlation existed between animal individuals, i.e., they were not from the same family. All samples were collected in strict compliance with the Chinese Wildlife Conservation Act. Also, the Sichuan Agricultural University Ethics Committee provided ethical approval for this study. Genomic DNA was extracted from samples according to the standard phenol-chloroform method [68]. Products were examined with 1.0% agarose gel electrophoresis and were visualized under ultraviolet light.

PCR and SSCP typing

A pair of specific primers (Mamu-DQB1F 5’–GCCTGACTGACTGGACGGTGATTC-3’ and Mamu-DQB1R 5’-GGGGCGACAACGCTCACCTC-3’) were used to amplify an approximately 317-bp DNA fragment containing the complete second exon of the Mamu-DQB1 gene. For polymerase chain reaction (PCR) amplification of Mamu-DQB1, all PCRs were performed in a reaction volume of 10 μl, containing 5 μl of 2X ES Taq MasterMix (CWBIO, Beijing, China), 0.3 μl (10 pm/μl) of each DQB1 primer, roughly 20 ng template DNA, and distilled, deionized H2O to 10 μl. All PCR amplifications were carried out using a MyCyclerTM Thermal Cycler (BIO-RAD, USA), with the first step at 94°C for 5 minutes, 30 cycles of 94°C for 30 seconds, 62°C for 25 seconds, 72°C for 45 seconds, and a final extension at 72°C for 10 minutes. Products were examined with 1% agarose gel electrophoresis and visualized under ultraviolet light.

All positive PCR amplification products were screened for polymorphisms in DQB1 exon 2 by SSCP analysis. Two μl of the PCR sample was mixed with 8 μl of the denaturating buffer (95% formamide, 5% 6X loading buffer) before heating for 10 min at 98°C. The samples after denaturation were immediately cooled in ice water at least 10 min, then loaded in a 2-mm-thick 12% acrylamide gel (no glycerol contained) and were electrophorized in 0.5X TBE at 130 V (15 h) at room temperature. The resulting bands were visualized by silver staining. PCR products that showed the same SSCP pattern were counted and typed.

Cloning and sequencing

Representative samples from each population were amplified with a reaction volume of 50 μl by PCR; the product was fractionated by electrophoresis in 18% agarose gels and purified using a Gel Extraction Kit (Sangon Biotech, Shanghai, China). Purified pieces of DNA were then cloned using the pMD19-T Vector kit (Takara, Dalian, China). Briefly, ~2–3 μl of each positive purified DNA was added into each mixture containing 2.5 μl of Solution I and 0.5 μl of pMD19-T Vector. The mixtures were incubated at 16°C for 40 min. For transformation, we added 25 μl of Escherichia coli DH5α-competent cells to the incubated ligation mixture, and then we put this on ice for 30 minutes, heat-shocked the cells at 42°C for 52 seconds in a water bath, immediately put this on ice for at least 3 minutes, and added 800 μl 37°C LB media. This mixture was then incubated at 37°C with shaking at approximately 150 rpm for 80 min, and recovery cells were mixed with 7 μl IPTG and 40 μl X-Gal, then 160 μl of the mixture was plated onto LB/ampicillin plates and incubated at 37°C overnight (14–16 hours). After blue⁄white selection, a number of positive clones were picked and grown, and 10–12 clones were chosen and sequenced by the sequencing company (Invitrogen, Shanghai, China). To avoid errors arising from PCR-based recombination, apart from the optimization of the PCR program, we have taken several measures to determine a real sequence: 1) a sequence must be detected in at least 2 animal individuals; and/or 2) if a sequence was found in only one animal individual, it must be confirmed by at least 2 independent PCRs.

Data analysis

Mamu-DQB1diversity analyses

The sequences obtained from this study were initially adjusted using DNAStar software (DNASTAR Inc., Madison, WI). Multiple sequence alignments were created using the Clustal W program in DAMBE software [69]. Nucleotide diversity (π) and haplotype diversity (h) were calculated using DnaSP (Version 5) [70]. Allelic frequency distribution, expected (He) and observed (Ho) heterozygosities, pairwise FST and genetic flow (Nm) among populations were estimated with GenAlEx 6.5 [71]. Inbreeding coefficient (F IS ) and allelic richness (A R ) per population were calculated with the FSTAT program, version 2.9.3 [72]. Allelic richness is a measure of the number of alleles independent of sample size. Hence, it can be compared between samples of different sizes.

Selection pressure and phylogenetic analysis

Non-synonymous (d N ) and synonymous (d S ) substitution rates at Mamu-DQB1 amino acid sites were calculated according to the Nei–Gjobori method in Mega 5.05 [73]; the Z-test for selection was used to evaluate whether d N was significantly greater than d S for entire sequences. The ABS and non-ABS were distinguished based on HLA-DQB1 structure [74]. Phylogenetic trees were reconstructed for Mhc-DQB1 allelic sequences using Mega 5.05, based on the neighbor-joining method with Kimura 2-para model. The reliability of the neighbor-joining tree was estimated by bootstrapping, in which 1000 bootstrap pseudosamples were used. To expand our data and compare with the closely related species, an additional 46 Mamu-DQB1, 85 Mafa-DQB1 (Macaca fascicularis), and 2 Aona-DQB1 (Aotus nancymaae) alleles were also retrieved from the Immuno Polymorphism Database. Thus, a total of 156 Mhc-DQB1 gene sequences including 23 sequences identified in this study were used in the phylogenetic analysis, and Aona-DQB1 was used as an out-group. Information for all these sequences is shown in Additional file1: Table S1.

Mamu-DQB1genetic differentiation analyses

The degree of population subdivision was estimated on the basis of hierarchical analysis of molecular variance (AMOVA) using the Arlequin version 3.1 package [75]. This package was used to compare the component of genetic diversity from the variance among the sampled major geographical regions and among subpopulations within each region. In order to measure the role of balancing selection in maintaining MHC variation among the population, the microsatellite data sets of Li et al. [25] that analyzed microsatellite variation at 30 loci in six populations (including all populations in this study) were used, and correlations between Mhc-DQB1 and microsatellite data were assessed using the Pearson product–moment correlation in SPSS 16.0.

Availability of supporting data

The data sets supporting the results of this article are available in the Dryad repository, doi:10.5061/dryad.h8b93 [76].


  1. 1.

    Sommer S: Effects of habitat fragmentation and changes of dispersal behaviour after a recent population decline on the genetic variability of noncoding and coding DNA of a monogamous Malagasy rodent. Mol Ecol. 2003, 12 (10): 2845-2851. 10.1046/j.1365-294X.2003.01906.x.

    PubMed  CAS  Article  Google Scholar 

  2. 2.

    Gregorius HR: The relationship between the concepts of genetic diversity and differentiation. Theor Appl Genet. 1987, 74 (3): 397-401. 10.1007/BF00274724.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Aguilar A, Roemer G, Debenham S, Binns M, Garcelon D, Wayne R: High MHC diversity maintained by balancing selection in an otherwise genetically monomorphic mammal. Proc Natl Acad Sci USA. 2004, 101: 3490-3494. 10.1073/pnas.0306582101.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  4. 4.

    Kohn MH, Murphy WJ, Ostrander EA, Wayne RK: Genomics and conservation genetics. Trends Ecol Evol. 2006, 21: 629-637. 10.1016/j.tree.2006.08.001.

    PubMed  Article  Google Scholar 

  5. 5.

    Bonin A, Nicole F, Pompanon F, Miaud C, Taberlet P: Population adaptive index: a new method to help measure intraspecific genetic diversity and prioritize populations for conservation. Conserv Biol. 2007, 21: 697-708. 10.1111/j.1523-1739.2007.00685.x.

    PubMed  Article  Google Scholar 

  6. 6.

    Miller HC, Allendorf F, Daugherty CH: Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.). Mol Ecol. 2010, 19: 3894-3908. 10.1111/j.1365-294X.2010.04771.x.

    PubMed  Article  Google Scholar 

  7. 7.

    Siddle HV, Kreiss A, Eldridge MDB, Noonan E, Clarke CJ, Pyecroft S, Woods GM, Belov K: Transmission of a fatal clonal tumor by biting occurs due to depleted MHC diversity in a threatened carnivorous marsupial. Proc Natl Acad Sci USA. 2007, 104: 16221-16226. 10.1073/pnas.0704580104.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  8. 8.

    Spurgin LG, Richardson DS: How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings. Proc R Soc B. 2010, 277: 979-988. 10.1098/rspb.2009.2084.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  9. 9.

    Bernatchez L, Landry C: MHC studies in nonmodel vertebrates: what have we learned about natural selection in15 years?. J Evol Biol. 2003, 16: 363-377. 10.1046/j.1420-9101.2003.00531.x.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Piertney SB, Oliver MK: The evolutionary ecology of the major histocompatibility complex. Heredity. 2006, 96: 7-21.

    PubMed  CAS  Google Scholar 

  11. 11.

    Alcaide M, Lemus JA, Blanco G, Tella J, Serrand D, Negro JJ, Rodriguez A, Garcia-Montijano M: MHC diversity and differential exposure to pathogens in kestrels (Aves: Falconidae). Mol Ecol. 2010, 19: 691-705. 10.1111/j.1365-294X.2009.04507.x.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Campos JL, Posada D, Moran P: Genetic variation at MHC, mitochondrial and microsatellite loci in isolated populations of Brown trout (Salmo trutta). Conserv Genet. 2006, 7: 515-530. 10.1007/s10592-005-9063-z.

    CAS  Article  Google Scholar 

  13. 13.

    Babik W, Pabijan M, Arntzen JW, Cogalniceanu D, Durka W, Radwan J: Long-term survival of a urodele amphibian despite depleted major histocompatibility complex variation. Mol Ecol. 2009, 18: 769-781. 10.1111/j.1365-294X.2008.04057.x.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Luo MF, Pan HJ, Liu ZJ, Li M: Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana). BMC Evol Biol. 2012, 12: 207-10.1186/1471-2148-12-207.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Sommer S: The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005, 2: 1-18. 10.1186/1742-9994-2-1.

    Article  Google Scholar 

  16. 16.

    Oubory NJ, Pertoldi C, Loeschcke V, Bijlsma R, Hedrick PW: Conservation genetics in transition to conservation genomics. Trends Genet. 2010, 26 (4): 177-187. 10.1016/j.tig.2010.01.001.

    Article  Google Scholar 

  17. 17.

    Allendorf FW, Luikart G: Conservation and the Genetics of Populations. 2007, Malden, MA: Blackwell Publishing

    Google Scholar 

  18. 18.

    Frankham R: Challenges and opportunities of genetic approaches to biological conservation. Biol Conserv. 2010, 143 (9): 1919-1927. 10.1016/j.biocon.2010.05.011.

    Article  Google Scholar 

  19. 19.

    Frankham R: Genetics and extinction. Biol Conserv. 2005, 126 (2): 131-140. 10.1016/j.biocon.2005.05.002.

    Article  Google Scholar 

  20. 20.

    Jiang XL, Wang YX, Ma SL: Taxonomic revision and distributin of subspecies of rhesus monke (macaca mulatta) in China. Zool Res. 1991, 12: 241-247.

    Google Scholar 

  21. 21.

    Lewontin RC: The Genetic Basis of Evolutionary Change. 1974, New York: Columbia University Press

    Google Scholar 

  22. 22.

    Rockman MV: The QTN program and the alleles that matter for evolution: all that's gold does not glitter. Evolution. 2012, 66 (1): 1-17. 10.1111/j.1558-5646.2011.01486.x.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Morin PA, Luikart G: Wayne RK, the SNP workshop group: SNPs in ecology, evolution and conservation. Trends Ecol Evol. 2004, 19 (4): 208-216. 10.1016/j.tree.2004.01.009.

    Article  Google Scholar 

  24. 24.

    Xu HL, Wang YT, Cheng AC, Yao YF, Ni QY, Zeng W, Bi FJ, Yang ZX, Chen XY: Polymorphism of MHC-DPB1 gene exon 2 in rhesus macaques (Macaca mulatta). Yi Chuan. 2010, 32: 588-598. 10.3724/SP.J.1005.2010.00588.

    PubMed  CAS  Article  Google Scholar 

  25. 25.

    Li DY, Xu HL, Trask JS, Zhu Q, Cheng AC, Smith DG, George D, Zhang L: Genetic diversity and population structure in wild Sichuan rhesus macaques. Mol Biol Rep. 2013, 40: 3033-3041. 10.1007/s11033-012-2377-2.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Wu SJ, Luo J, Li QQ, Wang YQ, Robert WM, Blair C, Wu SF, Yue BS, Zhang YP: Ecological genetics of Chinese rhesus macaque in response to mountain building: all things are not equal. PLoS One. 2013, 8: 2-

    Google Scholar 

  27. 27.

    Leuchte N, Berry N, Köhler B, Almond N, LeGrand R, Thorstensson R, Titti F, Sauermann U: MhcDRB-sequences from cynomolgus macaques (Macaca fascicularis) of different origin. Tissue Antigens. 2004, 63 (6): 529-537. 10.1111/j.0001-2815.2004.0222.x.

    PubMed  CAS  Article  Google Scholar 

  28. 28.

    Sano K, Shiina T, Kohara S, Yanagiya K, Hosomichi K, Shimizu S, Anzai T, Watanabe A, Ogasawara K, Torii R, Kulski JK, Inoko H: Novel cynomolgus macaque MHC-DPB1 polymorphisms in three South-East Asian populations. Tissue Antigens. 2006, 67 (4): 297-306. 10.1111/j.1399-0039.2006.00577.x.

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Diaz D, Naegeli M, Rodriguez R, Nino-Vasquez JJ, Moreno A, Patarroyo ME, Pluschke G, Daubenberger CA: Sequence and diversity of MHC DQA and DQB genes of the owl monkey Aotus nancymaae. Immunogenetics. 2000, 51 (7): 528-537. 10.1007/s002510000189.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Otting N, de Groot NG, Doxiadis GG, Bontrop RE: Extensive Mhc-DQB variation in humans and non-human primate species. Immunogenetics. 2002, 54 (4): 230-239. 10.1007/s00251-002-0461-9.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Qiu CL, Yang GB, Yu K, Li Y, Li LX, Liu Q, Zhao H, Xing H, Shao YM: Characterization of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles in a cohort of Chinese rhesus macaques (Macaca mulatta). Hum Immunol. 2008, 69: 513-521. 10.1016/j.humimm.2008.05.014.

    PubMed  CAS  Article  Google Scholar 

  32. 32.

    Hansson B, Westerberg L: On the correlation between heterozygosity and fitness in natural populations. Mol Ecol. 2002, 11: 2467-2474. 10.1046/j.1365-294X.2002.01644.x.

    PubMed  Article  Google Scholar 

  33. 33.

    Reed DH, Frankham R: Correlation between fitness and genetic diversity conserve. Biogeosciences. 2003, 17: 230-237.

    Google Scholar 

  34. 34.

    Leimu R, Mutikainen P, Koricheva J, Fischer M: How general are positive relationships between plant population size, fitness and genetic variation?. J Ecol. 2006, 94: 942-952. 10.1111/j.1365-2745.2006.01150.x.

    Article  Google Scholar 

  35. 35.

    Joron M, Brakefield PM: Captivity masks inbreeding effects on male mating success in butterflies. Nature. 2003, 424: 191-194. 10.1038/nature01713.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Booy G, Hendriks RJJ, Smulders MJM, Van Groenendael JM, Vosman B: Genetic diversity and the survival of populations Plant. Biogeosciences. 2000, 2: 379-395.

    Google Scholar 

  37. 37.

    Smith DG, McDonough J: Mitochondrial DNA variation in Chinese and Indian rhesus macaques (Macaca mulatta). Am J Primatol. 2005, 65 (1): 1-25. 10.1002/ajp.20094.

    PubMed  CAS  Article  Google Scholar 

  38. 38.

    Doxiadis GG, Otting N, de Groot NG, de Groot N, Rouweler AJ, Noort R: Evolutionary stability of MHC class II haplotypes in diverse rhesus macaque populations. Immunogenetics. 2003, 55: 540-10.1007/s00251-003-0590-9.

    PubMed  CAS  Article  Google Scholar 

  39. 39.

    Miller KM, Kaukinen KH, Beacham TD, Withler RE: Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon. Genetics. 2001, 111: 237-257.

    CAS  Google Scholar 

  40. 40.

    Bryja J, Charbonnel N, Berthier K, Galan M, Cosson JF: Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles. Mol Ecol. 2007, 16: 5084-5097. 10.1111/j.1365-294X.2007.03584.x.

    PubMed  CAS  Article  Google Scholar 

  41. 41.

    Ekblom R, Saether SA, Jacobsson P, Fiske P, Sahlman T, Grahn M, Kalas JA, Hoglund J: Spatial pattern of MHC class II variation in the great snipe (Gallinago media). Mol Ecol. 2007, 16: 1439-1451. 10.1111/j.1365-294X.2007.03281.x.

    PubMed  Article  Google Scholar 

  42. 42.

    Alcaide M, Edwards SV, Negro JJ, Serrano D, Tella JL: Extensive polymorphism and geographical variation at a positively selected MHC class IIB gene of the lesser kestrel (Falco naumanni). Mol Ecol. 2008, 17: 2652-2665. 10.1111/j.1365-294X.2008.03791.x.

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Penn DJ: The scent of genetic compatibility: sexual selection and the major histocompatibility complex. Ethology. 2002, 108: 1-21. 10.1046/j.1439-0310.2002.00768.x.

    Article  Google Scholar 

  44. 44.

    Milinski M: The major histocompatibility complex, sexual selection, and mate choice. Annu Rev Ecol Evol Syst. 2006, 37: 159-186. 10.1146/annurev.ecolsys.37.091305.110242.

    Article  Google Scholar 

  45. 45.

    Huchard E, Baniel A, Schliehe-Diecks S, Kappeler PM: MHC-disassortative mate choice and inbreeding avoidance in a solitary primate. Mol Ecol. 2013, 22 (15): 4071-4086. 10.1111/mec.12349.

    PubMed  Article  Google Scholar 

  46. 46.

    Klein J: Origin of major histocompatibility complex polymorphism: the trans-species hypothesis. Hum Immunol. 1987, 19: 155-162. 10.1016/0198-8859(87)90066-8.

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Kamath PL, Getz WM: Adaptive molecular evolution of the major histocompatibility complex genes, DRA and DQA, in the genus Equus. BMC Evol Biol. 2011, 11: 128-10.1186/1471-2148-11-128.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  48. 48.

    Radwan J, Biedrzycka A, Babik W: Does reduced MHC diversity decrease viability of vertebrate populations?. Biol Conserv. 2010, 143: 537-544. 10.1016/j.biocon.2009.07.026.

    Article  Google Scholar 

  49. 49.

    Wei K, Zhang ZH, Wang XF, Zhang WP, Xu X, Shen FJ, Yue BS: Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline Mhc-DRB peptide-binding region. Immunogenetics. 2010, 62: 307-317. 10.1007/s00251-010-0440-5.

    PubMed  CAS  Article  Google Scholar 

  50. 50.

    Yao YF, Zhao JJ, Dai QX, Li JY, Zhou L, Wang YT, Ni QY, Zhang MW, Xu HL: Identification and characterization of the major histocompatibility complex class II DQB (MhcMath-DQB1) alleles in Tibetan macaques (Macaca thibetana). Tissue Antigens. 2013, 82 (2): 113-121. 10.1111/tan.12145.

    PubMed  CAS  Article  Google Scholar 

  51. 51.

    van Oosterhout C: Trans-species polymorphism, HLA-disease associations and the evolution of the MHC. Commun Integr Biol. 2009, 2: 408-410. 10.4161/cib.2.5.8765.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  52. 52.

    Figureroa F, Gunther E, Klein J: MHC polymorphism pre-dating speciation. Nature. 1988, 335: 265-267. 10.1038/335265a0.

    Article  Google Scholar 

  53. 53.

    Hayasaka K, Fujii K, Horai S: Molecular phylogeny of macaques: implications of nucleotide sequences from an 896 base pair region of mitochondrial DNA. Mol Biol Evol. 1996, 13: 1044-1053. 10.1093/oxfordjournals.molbev.a025655.

    PubMed  CAS  Article  Google Scholar 

  54. 54.

    Yan G, Zhang G, Fang X, Zhang Y, Li C, Ling F, Cooper DN, Li Q, Li Y, van Gool AJ, Du H, Chen J, Chen R, Zhang P, Huang Z, Thompson JR, Meng Y, Bai Y, Wang J, Zhuo M, Wang T, Huang Y, Wei L, Li J, Wang Z, Hu H, Yang P, Le L, Stenson PD, Li B, et al: Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques. Nat Biotechnol. 2011, 29 (11): 1019-1023. 10.1038/nbt.1992.

    PubMed  CAS  Article  Google Scholar 

  55. 55.

    Wright S: The genetical structure of populations. Ann Eugenics. 1949, 15: 323-354. 10.1111/j.1469-1809.1949.tb02451.x.

    Article  Google Scholar 

  56. 56.

    Yasukochi Y, Nishida S, Han SH, Kurosaki T, Yoneda M, Koike H: Genetic structure of the Asiatic black bear in Japan using mitochondrial DNA analysis. J Hered. 2009, 100: 297-308. 10.1093/jhered/esn097.

    PubMed  CAS  Article  Google Scholar 

  57. 57.

    Wright S: The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution. 1965, 19: 395-420. 10.2307/2406450.

    Article  Google Scholar 

  58. 58.

    Millar C, Libby W: Strategies for conserving clinal, ccotypic, and disjunct population diversity in widespread species. Genetics and Conservation of Rare Plants. Edited by: Fald DA, Holsinger KE. 1991, New York: Oxford University Press, 149-170.

    Google Scholar 

  59. 59.

    Slatkin M: Gene flow and the geographic structure of natural populations. Science. 1987, 236 (4803): 787-10.1126/science.3576198.

    PubMed  CAS  Article  Google Scholar 

  60. 60.

    Zhu L, Zhang S, Gu X, Wei F: Significant genetic boundaries and spatial dynamics of giant pandas occupying fragmented habitat across southwest China. Mol Ecol. 2011, 20: 1122-1132. 10.1111/j.1365-294X.2011.04999.x.

    PubMed  Article  Google Scholar 

  61. 61.

    Yao YF, Zhong LJ, Li BF, Li JY, Ni QY, Xu HL: Genetic variation between two Tibetan macaque (Macaca thibetana) populations in the eastern China based on mitochondrial DNA control region sequences. Mitochondrial DNA. 2013, 24 (3): 267-275. 10.3109/19401736.2012.748040.

    PubMed  CAS  Article  Google Scholar 

  62. 62.

    Schierup MH, Hein J: Consequences of recombination on traditional phylogenetic analysis. Genetics. 2000, 156 (2): 879-891.

    PubMed  CAS  PubMed Central  Google Scholar 

  63. 63.

    Luikart G, England PR, Tallmon D, Jordan S, Taberlet P: The power and promise of population genomics: from genotyping to genome typing. Nat Rev Genet. 2003, 4: 981-994.

    PubMed  CAS  Article  Google Scholar 

  64. 64.

    Chang ZF, Luo MF, Liu ZJ, Yang JY, Xiang ZF, Li M, Vigilant L: Human influence on the population decline and loss of genetic diversity in a small and isolated population of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Genetica. 2012, 140: 105-114. 10.1007/s10709-012-9662-9.

    PubMed  Article  Google Scholar 

  65. 65.

    Oliver MK, Lambin X, Cornulier T, Piertney SB: Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations. Mol Ecol. 2009, 18: 80-92.

    PubMed  CAS  Google Scholar 

  66. 66.

    Klein J, Bontrop RE, Dawkins RL, Erlich HA, Gyllensten UB, Heise ER, Jones PP, Parham P, Wakeland EK, Watkins DI: Nomenclature for the major histocompatibility complexes of different species: A proposal. Immunogenetics. 1990, 31: 217-219.

    PubMed  CAS  Google Scholar 

  67. 67.

    Mason RAB, Browning TL, Eldridge DB: Reduced MHC class II diversity in island compared to mainland populations of the black-footed rockwallaby (Petrogale lateralis lateralis). Conserv Genet. 2011, 12: 91-103. 10.1007/s10592-009-9993-y.

    Article  Google Scholar 

  68. 68.

    Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning, A Laboratory Manual. 1989, New York: Cold Spring Harbor Laboratory Press

    Google Scholar 

  69. 69.

    Xia X, Xie Z: DAMBE: data analysis in molecular biology andevolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.

    PubMed  CAS  Article  Google Scholar 

  70. 70.

    Rozas R, Sánchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.

    PubMed  CAS  Article  Google Scholar 

  71. 71.

    Peakall R, Peter E: Smouse: genetic analysis in excel, population genetic software for teaching and research. Bioinformatics. 2012, 28 (19): 2537-2539. 10.1093/bioinformatics/bts460.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  72. 72.

    Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). 2001, Available from

    Google Scholar 

  73. 73.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA 5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  74. 74.

    Amills M, Ramirez O, Tomas A, Obexer-Ruff G, Vidal O: Positive selection on mammalian MHC-DQ genes revisited from a multispecies perspective. Genes Immun. 2008, 9: 651-658. 10.1038/gene.2008.62.

    PubMed  CAS  Article  Google Scholar 

  75. 75.

    Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform. 2005, 1: 47-50.

    CAS  Google Scholar 

  76. 76.

    Yao YF, Dai QX, Li J, Ni QY, Zhang MW, Xu HL: Data from: Genetic diversity and differentiation of the rhesus macaque (Macaca mulatta) population in western Sichuan, China, based on the second exon of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles. Dryad Digital Repository. 2014, doi:10.5061/dryad.h8b93

    Google Scholar 

Download references


We thank An-Chun Cheng, Ming-Shu Wang and Wen Zeng for helping to collect samples; This work was supported by the National Natural Science Foundation of China (No. 31370407, 30970383), Natural Science Foundation of Educational Commission of Sichuan Province of China (No. 08ZA076), the academic leader training projects in Sichuan Province of China (No. 01320206) and the Innovative Research Team in University of Sichuan Bureau of Education.

Author information



Corresponding author

Correspondence to Huai-Liang Xu.

Additional information

Competing interests

The authors have declared that no competing interests.

Authors’ contributions

H-LX conceived the study and helped to draft the manuscript; Y-FY and Q-XD performed the experiments and wrote the manuscript; JL, Q-YN and M-WZ analyzed the data. All authors read and approved the final manuscript.

Yong-Fang Yao, Qiu-Xia Dai contributed equally to this work.

Electronic supplementary material


Additional file 1: Table S1: Mhc-DQB1 alleles of different primate species used in this study A total of 156 Mhc-DQB1 alleles from three primate species (Macaca mulatta, Macaca fascicularis, Aotus nancymaae) were used in this study, including 69 Mamu-DQB1, 85 Mafa-DQB1 and 2 Aona-DQB1 alleles. Alleles detected in this study are shown in bold. IPD Acc No. is accession number in the IPD-MHC Database ( (PDF 13 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yao, Y., Dai, Q., Li, J. et al. Genetic diversity and differentiation of the rhesus macaque (Macaca mulatta) population in western Sichuan, China, based on the second exon of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles. BMC Evol Biol 14, 130 (2014).

Download citation


  • Genetic diversity
  • Genetic differentiation
  • Macaca mulatta
  • Major histocompatibility complex
  • Trans-species evolution