Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana)
© Luo et al.; licensee BioMed Central Ltd. 2012
Received: 28 February 2012
Accepted: 5 October 2012
Published: 19 October 2012
Small, isolated populations often experience loss of genetic variation due to random genetic drift. Unlike neutral or nearly neutral markers (such as mitochondrial genes or microsatellites), major histocompatibility complex (MHC) genes in these populations may retain high levels of polymorphism due to balancing selection. The relative roles of balancing selection and genetic drift in either small isolated or bottlenecked populations remain controversial. In this study, we examined the mechanisms maintaining polymorphisms of MHC genes in small isolated populations of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) by comparing genetic variation found in MHC and microsatellite loci. There are few studies of this kind conducted on highly endangered primate species.
Two MHC genes were sequenced and sixteen microsatellite loci were genotyped from samples representing three isolated populations. We isolated nine DQA1 alleles and sixteen DQB1 alleles and validated expression of the alleles. Lowest genetic variation for both MHC and microsatellites was found in the Shennongjia (SNJ) population. Historical balancing selection was revealed at both the DQA1 and DQB1 loci, as revealed by excess non-synonymous substitutions at antigen binding sites (ABS) and maximum-likelihood-based random-site models. Patterns of microsatellite variation revealed population structure. FST outlier analysis showed that population differentiation at the two MHC loci was similar to the microsatellite loci.
MHC genes and microsatellite loci showed the same allelic richness pattern with the lowest genetic variation occurring in SNJ, suggesting that genetic drift played a prominent role in these isolated populations. As MHC genes are subject to selective pressures, the maintenance of genetic variation is of particular interest in small, long-isolated populations. The results of this study may contribute to captive breeding and translocation programs for endangered species.
KeywordsBalancing selection Conservation genetics Gene drift MHC Rhinopithecus roxellana
Understanding how levels of genetic variation influence the survival of threatened species is of fundamental interest to evolutionary and conservation biologists because many natural populations are threatened by intense reduction and fragmentation of habitat, leading to isolation, declining populations, and decreasing genetic diversity [1, 2]. Loss of genetic diversity may increase the risk of extinction due to decreased reproductive fitness, decreased adaptive flexibility, and increased disease susceptibility . Clarifying the mechanism that determines genetic variation in small, isolated populations is therefore essential for their conservation . An important assumption in conservation genetics is that small, isolated populations are more sensitive to genetic drift and inbreeding [5, 6]. 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. 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 .
However, some functionally important genes that are maintained by balancing selection, such as major histocompatibility complex (MHC) genes, may have a different evolutionary pattern compared with neutral markers. The multi-gene MHC family is found in vertebrates, codes for cell surface glycoproteins, and is important in animal conservation due to its role in resisting pathogens . Compared to nearly neutral markers such as microsatellite loci or mitochondrial DNA, which are informative for phylogenetic and phylogeographic reconstructions , MHC variability is believed to determine the capability of individuals to resist continuously evolving pathogens and parasites. Consequently, MHC variability is a reflection of the processes that are related to adaptive evolution within and between populations . Thus, most variation at MHC loci reflects the effects of balancing selection , which is the main mechanism for retaining high MHC genetic diversity. Balancing selection includes frequency-dependent selection, overdominance and diversifying selection and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages in mammals . Besides balancing selection, intragenic recombination has also been suggested as one evolutionary mechanism for generation of MHC sequence diversity [11, 12]. Still, several significant questions remain. Current results conflict or are unclear with regard to the relative roles of balancing selection and genetic drift in maintaining MHC polymorphism in small, isolated, or severely bottlenecked populations [13, 14]. Studies of guppies, Mexican wolves and Namibian leopards found that even when genetic diversity at neutral markers was poor, polymorphisms in MHC were still maintained by balancing selection [13–17]. In contrast, other research has found low levels of detectable polymorphisms for MHC genes in populations with lower diversity in neutral markers, including studies on fallow deer (Cervus dama) , northern elephant seals (Mirounga angustirostris) , great crested newt (Triturus cristatus) , and the black-footed rock-wallaby (Petrogale lateralis lateralis) . These results suggest that compared to genetic drift balancing selection is relatively weak in small populations, leading to reduced variation at some MHC loci . To elucidate the role of balancing selection and genetic drift in populations, FST outlier analysis  is widely used because demographic processes affect neutral loci and lead to population differentiation (measured by FST) . The FST values are computed for all genes to distinguish the genes under selection from those non-selective genes .
Genetic variation at microsatellites
Population genetic parameters for populations estimated from microsatellite data
No. of Loci
MHC variation, expression analyses, and recombination analysis
For each locus, we analysed 894 clones from 64 individuals that represented three monkey populations. We obtained 13 different DQA1 fragments with lengths that ranged from 436 bp to 444 bp including exon 2 and partial intron 2 and intron 3 [GenBank: JQ217094-JQ217106]. After performing BLAST with the HLA-DQA1 and Macaca mulatta DQA1 exon 2, we removed the intron and obtained 9 DQA1 exon 2 sequences with an equal length of 249 bp. These sequences were labelled as Rhro-DQA1*01–09 [GenBank: JQ217107-JQ217115] according to the nomenclature of Klein et al. and were used for analysis. At the DQB1 locus, we identified 16 unique DQB1 exon 2 sequences that we called Rhro-DQB1*01–16 [GenBank: JQ217116- JQ217131].
Fragments of MHC loci were successfully amplified from the cDNA of 2 blood samples that were also used in DNA amplifications. At the DQA1 locus, 189 bp fragments of DQA1 exon 2 were obtained from 2 blood samples, and these fragments were part of the alleles obtained from DNA in the same individual [EMBL: HE616682- HE616683]. At the DQB1 locus longer fragments were obtained, including the entire exon 2, partial exon 1, and exon 3; however, these fragments were only obtained from 1 blood sample [GenBank: JQ217132-JQ217133]. Shorter fragments of DQB1 exon 2 were obtained from 2 blood samples [EMBL: HE616684- HE616686]. No sequences contained stop codons. The sequences from RNA were identical with that from DNA, which supports the hypothesis that these MHC PCR products were transcribed and expressed.
Summary of MHC variation in R. roxellana populations
Parsimony- informative sites
Over mean distance
Population structure and phylogenetic analysis
Average non-synonymous substitutions per non-synonymous site ( d N ) and synonymous substitutions per synonymous sites ( d S )
Results of maximum-likelihood models for exon 2 of the DQA1 gene
Positively selected sites
p0= 0.723 (p1= 0.277)
p0= 0.911, p1= 0.033(p2= 0.056) ω2= 5.148
p0= 0.736, p1= 0.244 (p2= 0.018) ω1= 1.298, ω2= 8.174
p= 0.008, q= 0.028
M8(beta and omega)
p0= 0.974 (p1= 0.025) p= 0.008, q= 0.023, ω= 7.292
Results of maximum-likelihood models for exon 2 of the DQB1 gene
Positively selected sites
p0= 0.869 (p1= 0.131)
p0=0.818, p1=0.176(p2= 0.006) ω2= 8.276
p1= 0.059 (p2= 0.018) ω1= 0.076, ω2= 3.390
p= 0.006, q= 0.036
M8(beta and omega)
p0= 0.994 (p1= 0.006) p= 0.005, q= 0.020, ω= 8.643
Habitat isolation generates barriers to gene flow among populations that often result in loss of genetic diversity through genetic drift and inbreeding . For vertebrates, genetic variation is of special importance in MHC genes due to the significant role they play in immune functions . MHC diversity is presumed to improve parasite resistance, reproductive success, and population viability , and has been studied in species such as western gorilla (Gorilla gorilla) , brown bear (Ursus arctos) , Ethiopian wolf (Canis simensis) , European bison (Bison bonasus) , Bengal tiger (Panthera tigris tigris) , and Namibian leopard (Panthera pardus pardus) . In this study, a relatively high level of MHC variation was found in the golden snub-nosed monkey, with 9 DQA1 and 16 DQB1 alleles found in 64 individuals. However, we found lower MHC variation in the SNJ population, which also showed lower genetic diversity in microsatellites and mitochondrial genes [27, 34]. Small sample size cannot account for the reduced variability, because sample collected across the entire distribution of the SNJ population showed similar patterns of genetic diversity .
For each MHC locus, regardless of how many clones were sequenced from an individual, no more than 2 alleles were observed in an individual, a strong indicator that we amplified single loci in all cases. We assumed that all alleles were from a single functional gene. The assumption that our sequences were from functional genes was supported by three findings: 1) sites inferred to have been exposed to significant selection, most of which were ABS sites, indicating historical selection at functional genes; 2) no reading frame or stop codon disruptions found in any alleles; and 3) the sequences identified from cDNA were parts of or included sequences obtained from DNA.
Historical balancing selection
Golden snub-nosed monkey MHC genes reflect historical balancing selection in that an excess of non-synonymous substitutions was mainly concentrated in the ABS (Table 3, Additional file 4: Table S3, Additional file 5: Table S4). According to neutrality theory , the synonymous nucleotide substitution rate is larger than the non-synonymous substitution rate because a change in amino acid sequence has a greater possibility of being deleterious. The elevated rate of non-synonymous substitutions at the ABS provided clear evidence of positive selection [43, 44] shaping genetic variation . The p-value was not significant at DQA1 ABS may due to a weaker recombination within DQA1 as a lower recombination rate has been shown before. Without higher recombination and stronger selection, some DQA1 alleles might be lost when population size decreased. The sharing of MHC alleles among populations also indicates that MHC alleles may have been conserved by selection . Second, random site models analysis proved the existence of historical selection based on the maximum likelihood method, which revealed that, for MHC genes, the models including selection (M2a, M3, and M8) match MHC alleles better than models without selection (Table 4,and 5). Under the M2a and M8 models, some sites of the 2 MHC loci were under significant selection pressure. Furthermore, trans-species evolution of the MHC alleles revealed historical balancing selection. Under balancing selection, some MHC alleles or allelic lineages are reported in other species, which indicates that they are ancestral alleles .
Genes vary in terms of the level of selection, and both DQA1 and DQB1 revealed different patterns of selection. Each population had unique DQB1 alleles, while not every population had unique DQA1 alleles (Table 2). Population divergence, measured as pairwise F ST , was larger in DQB1 than in DQA1 except for that between SG and QL (Figure 3). Similar results were reported in water vole, where balancing selection pressure was different at MHC genes in continuous populations . MHC genes are assumed to be closely linked . In our study, however, linkage disequilibrium between DQA1 and DQB1 was only observed in the SG population. This weak linkage disequilibrium and the different selection pressures on these loci could be a result of recombination, which is common at the MHC genes [49, 50]. The higher recombination rate that was found in DQB1 genes may explain their larger allelic richness compared to DQA1 genes. Recombinants maintained by selection may counteract the linkage of closely linked genes  and play an important adaptive role in DQB1 evolution. In the present study, historical selection was found, but this does not conclusively indicate that balancing selection is acting on current populations. First, an excess of non-synonymous mutations requires a long time to accumulate . Once present, this pattern would take a long time to vanish after the disappearance of selection . Hence, we investigated whether selection continues to play a major role at present.
Patterns of selection and drift
Although selection historically maintained MHC diversity, recent population isolation and fragmentation has increased the role of genetic drift in shaping patterns of MHC variation in snub-nosed monkeys. First, compared with neutral forces, balancing selection is supposed to diminish population differentiation as measured by conventional pairwise FST[22, 52]. Thus, the population structure of genes under balancing selection should not be pronounced . However, in the present study, half of the pairwise FST values were greater than 0.05, and two FST values at DQB were even greater than those at microsatellites (Figure 3). Second, FST outlier analysis showed that the structure in the MHC loci was within the neutrality level for all populations and for each population. Considering all populations, one microsatellite (D14S306) showed a FST value that was lower than the neutral level, indicating its linkage with other genes under selection . Lastly, a positive correlation was found between allelic richness in MHC and microsatellites. The QL population had the highest allelic richness in microsatellites and MHC, while SNJ had the lowest. This positive correlation indicates that genetic drift plays a significant role in maintaining MHC diversity for snub-nosed monkeys . Maintenance of MHC variation through balancing selection may be hampered in small, isolated populations because of their lower effective recombination rate . In all, our results indicate that even though selection acts on MHC, it is overwhelmed by genetic drift in small, isolated populations.
The positive correlation in the allelic richness of MHC and microsatellites, together with other evidence, indicates that genetic drift has a great influence on the maintenance of MHC variations in small, isolated populations of snub-nosed monkeys [7, 25, 55]. No evidence showed that MHC polymorphism had increased in populations that contained low neutral variation [27, 42]. Under neutral evolution theory , alleles are expected to be neutral when s < 1/2 Ne (s = selection coefficient, Ne = effective population size). Therefore, the smaller Ne becomes, the greater the likelihood of genetic drift . The SNJ population is subject to more genetic drift than the other two populations as found in a previous study . Other animals whose patterns of MHC polymorphism have been contributed to drift over selection include the great crested newt (Triturus cristatus) , black-footed rock-wallaby (Petrogale lateralis lateralis) , tuatara (Sphenodon spp.) , and the Egyptian vulture (Neophron percnopterus) . These results suggest that selection on MHC is not strong enough to counteract drift that results from population fragmentation, isolation and bottleneck.
Variation at MHC loci is widely accepted as being maintained by balancing selection (reviewed in 2), even with a low level of neutral variability in some species [29, 59]. This shows the importance of balancing selection for maintaining variation in the field and exposes the problem of using neutral genes as substitutes for variation in fitness-related genes . However, in small, isolated populations or bottlenecked populations, balancing selection is overwhelmed by drift . In this study, we found the same genetic variation pattern both at neutral and MHC markers, suggesting that genetic drift was stronger than selection, thus leading to a reduction in MHC diversity in the most isolated populations. Such findings may contribute to the conservation of endangered species such as snub-nosed monkeys both in captive breeding and translocation programs. Though the relationships between MHC with mate choice and pregnancy outcomes still remain controversial , it could be helpful to examine MHC variation in captive individuals that are involved in mating programs. In the field, translocation of individuals from demographically and genetically healthy populations to populations that suffer from reduced genetic diversity can improve the chances of genetic recovery [61, 62]. Recently, researchers have found that inbred populations could thrive with the import of migrants as part of a conservation management program [63–65]. Further, the introduction of individuals from western populations of golden snub-nosed monkeys could restore genetic diversity to the relatively homogenous SNJ groups, reducing the likelihood of inbreeding depression.
Sample collection and DNA extraction
We collected 64 R. roxellana samples (muscle, skin, and blood) from the current population (SG population = 25, SNJ population = 17, and QL population = 22) (Figure 1 and Table 1). All collections complied with the relevant animal welfare institutions and laws of China. Also, the Institute of Zoology (Chinese Academy of Sciences) provided ethical approval for this study. Muscle and skin samples were collected from carcasses that were provided by local museums and nature reserves. Skin samples were stored dry, muscle samples were stored in 95% ethanol, and blood samples were collected while trapping individuals for physical examination and were stored in a refrigerator at −80°C. During DNA extraction, benches and plasticware were washed with 10% bleach and sterile water and were exposed to UV light for 30 min prior to treatment to prevent contamination during DNA extraction. During the subsequent polymerase chain reaction (PCR), 8 extraction controls were used without any positive amplification.
Samples were genotyped at 16 microsatellite loci . Thirty-five cycles of PCR amplification were executed at the same time for up to 3 loci with reliable genotypes, with combinations selected by fluorescent dye (HEX, ISMRA, or FAM), Tm, and fragment size using the QIAGEN Mutliplex PCR kit following the manufacturer’s protocols at optimized annealing temperatures (55°C). An ABI 377 prism automated sequencer was used to resolve products, which were then analysed using GeneScan v3.1.2 and Genotyper 2.5 (Applied Biosystems).
MHC amplification, cloning, sequencing, and expression analyses
The PCR was carried out in a 50 μL solution including 10 mM Tris–HCl (pH 8.4), 50 mM KCl, 2.5 mM MgCl2, 0.4 μM each primer, 0.2 mM each dNTP, 1.0 unit Hotstart-Taq DNA polymerase (Takara), and 10–100 ng DNA template. The amplification profile consisted of 5 min at 94°C, followed by 35 cycles of 30 s at 94°C, 30 s at 56°C (DQA1) or 58°C (DQB1), and 30 s at 72°C, with a final extension of 10 min at 72°C. The following primers were used for amplifying exon 2: 5'DQA1-AAGCCCA TAATATT TGAAAGTCAGT and 3'DQA1-TATGTGATTTTAGAGATGGGAGATG, or 5'DQB1-TGTAAAC GACGGCCAGTTCCCCGCAGAGGATTTCGTG and 3'DQB1-TGCTCTAGAGGGCGACGACGACGCCTCACCTC . The DQA1 primers were designed based on BLAST sequences from relative species. According to the manufacturer’s protocols, a Wizard PCR Preps DNA Purification Kit (Promega) was used to purify the PCR products. Purified PCR products were cloned using the pMD-18T vector (Takara) following manufacturer’s instructions. Ten to twenty clones containing inserts from each individual were sequenced on an ABI 377 or ABI-PRISMTM 3100 Genetic Analyzer (Applied Biosystems Inc.) with the Prism BigDyeTM Terminator Ready Reaction Kit (Applied Biosystems Inc.).
Expression analyses were conducted to validate the expression of the obtained MHC II-DQA and DQB alleles. The RNA was isolated from 2.5 mL whole blood using Trizol (Invitrogen) following the manufacturer’s instructions. To ensure that the genomic DNA was removed from the isolated RNA, a second DNA digestion was performed using the DNase I RNase-free Set (Promega). The cDNA synthesis was obtained using 200 U of M-MuLV Reverse Transcriptase (Promega) in a 25 μL reaction tube containing 1μL Oligo(dT)12–18 primer (0.5 μg/μL; Invitrogen), 2 μg total RNA as a template, and 25 U ribonuclease inhibitor, RNase-free water, 5 μL dNTP mix (10 mM), and 5 μL 5× reaction buffer. The reactions were incubated at 42°C for 60 min. The cDNA was acquired from 2 samples and was amplified by PCR using the primer sets (Additional file 6: Table S5). The primer sets were exon-spanning to detect the amplification of genomic DNA contaminants based on the BLAST results of human and Macaca mulatta exons. All amplified cDNA products were analysed by cloning and sequencing as described above.
Microsatellite loci diversity analyses
Genetic diversity was measured as observed (HO) and expected heterozygosities (HE) . The above analyses were computed using Arlequin v3 . Inbreeding coefficients (FIS) and allelic richness (AR) analyses were performed in Fstat version 2.9.3 . The AR calculations were set to the smallest population size correction. Pairwise FST among populations and Hardy–Weinberg equilibrium were analysed across all loci for each population in an exact probability test in Genepop 4.0 . Using the Bonferroni correction, significance values were adjusted for multiple comparisons. Linkage disequilibrium (LD) was calculated across all loci with the Weir correlation coefficient  in Genepop 4.0 . To determine which LD values were significant, a permutation test was used. The presence of stuttering, null alleles, and allele dropout were examined using Micro-Checker .
MHC diversity analyses
Sequences were aligned and translated into the amino acid sequences in Mega 4 . The gene identity was verified through homology with publicized MHC alleles of other species using BLASTN (http://blast.ncbi.nlm.nih.gov/Blast.cgi) from NCBI. Because some obtained sequences are false alleles potentially corresponded to PCR amplification artifacts , a new sequence was assumed to be an allele when it was identified from 3 separate amplications from the same sample or from at least 3 different samples . We used Mega 4 to detect the variable and parsimony-informative sites (sites with at least 2 different nucleotides or amino acids), to compute the mean number of nucleotide differences, and to derive the overall mean genetic distances of nucleotide sequences. Standard estimate errors were detected using 1,000 bootstrap replicates. Hardy–Weinberg equilibrium tests, pairwise FST values, and LD were calculated in Genepop. The AR and FIS of DQA1 and DQB1 were compared among populations using the comparison among groups of samples option in Fstat version 2.9.3 . The AR estimates were adjusted to the smallest sample size.
Phylogenetic analysis and recombination analysis
To detect the genetic structures of microsatellite loci, Bayesian clustering was used in Structure 2.1 . Because of the subtle population structure, correlated allele frequencies between populations and the admixture model were chosen . The clusters (K) tested range was set from 1 to 10, and 10 independent runs were executed for each analysis. The Markov Chain Monte Carlo (MCMC) iterations lengths and burn-ins were set at 5,000,000 and 200,000, respectively. For most situations, ΔK had a mode at true K; thus, true K was selected using the ΔK statistic, which was calculated according to the rate of the log probability change of the data between successive K values .
The relationship of the MHC allele phylogenesis was constructed using a maximum likelihood method in PhyML 3.0 . Before phylogenetic analysis, the most appropriate evolution models of sequence were estimated based on Akaike Information Criterion (AIC) in Modeltest 3.7 . The K80 + G model was suggested to be the optimum model for MHC sequences with the gamma shape parameter α = 0.6532. The reliability of the obtained tree topology structure was carried out with 1,000 bootstrap replications. In addition, intraspecific phylogenetic structures were inferred using median-joining networks in Network 4.600 . To estimate the rate of population recombination ρ (ρ = 4Ner), the composite-likelihood method  in LDhat  was used. The ρ is calculated using the crossing over for an effective population size (Ne) and rate per generation (r), and is estimated without prior information . Even for the sequences that evolved under balancing selection in the existence of recombination events, the LDhat still works efficiently .
Detecting historical selection
Two methods were used to detect historical selection. First, we calculated the rates of non-synonymous and synonymous substitutions at all amino acid sites, ABS sites, and non-ABS in Mega 4 by the Nei-Gojobori method with the Jukes-Cantor correction  and 1,000 bootstrap replicates to obtain standard errors. The putative ABS and non-ABS locations were derived according to the structure of human MHC II . Compare with the method that used HLA-DR1  structure to derive ABS sites, Reche and Reinherz’s research provided a better way to avoid the emergence of uncertain gaps while blasting DQ1 genes with DR1 genes, and was also adopted by other studies [87, 88]. Historical selection evidence was obtained with Codeml in the Paml 4 package . This procedure examined heterogeneity in ω (ω = dN/dS)  among codons based on the maximum likelihood method, with positive selection indicated by ω = dN/dS > 1. Six models (M0, M2a, M3, M7, M8) allowing for different selection intensity among sites were tested [91, 92]. The sites revealed to be under selection were compared with ABS sites.
Detecting recent selection
The FST outlier analysis  was used to examine whether MHC genes had significantly different population differentiation measures compared with microsatellite loci. Loci outside the confidence intervals of neutral FST values were likely under selection. The FST was estimated for DQA1, DQB1, and 16 microsatellite loci using Lositan [22, 93]. Simulations involving 100 demes were sampled, and HE and FST were printed for each of the 100,000 simulations. The infinite allele model was used in simulations that were performed using Fdist2. The aforementioned FST outlier analysis was also carried out on each population respectively. For each population, simulation parameters were similar using the same parameter set.
This project was funded by National Key Technology R&D Program (Key technology research on germplasm resources conservation of Shennongjia golden snub-nosed monkey), Natural Science Foundation of China (31130061, 30970427), and Project of Public Benefit (201104073). Thanks to B.P. Ren, Z.F. Chang, B.H. Yang, B.S. Wang, L. Zhang, and L.F. Zhu for advice and lab assistance. Thanks to Edanz Group China and Alicia Krzton for English editing. Special thanks to the editor and two innominate reviewers for their kind suggestions for improving this paper.
- Peacock MM, Smith AT: The effects of habitat fragmentation on dispersal patterns, mating behavior and genetic variation in a pica (Ochotona princeps) metapopulation. Oecol. 1997, 112: 524-533. 10.1007/s004420050341.View ArticleGoogle Scholar
- Sommer S: The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005, 2: 16-10.1186/1742-9994-2-16.PubMedPubMed CentralView ArticleGoogle Scholar
- Allendorf FW, Luikart G: Units of conservation. In: Conservation and the genetics of populations. 2007, Oxford: Blackwell Publishing, 380-420.Google Scholar
- Frankham R: Genetics and extinction. Biol Conserv. 2005, 126: 131-140. 10.1016/j.biocon.2005.05.002.View ArticleGoogle Scholar
- Frankham R, Ballou JD, Briscoe DA: Introduction to Conservation Genetics (2nd edn). 2010, Cambridge: Cambridge University Press, 1-618.View ArticleGoogle Scholar
- Ouborg NJ, Pertoldi C, Loeschcke V, Bijlsma R, Hedrick PW: Conservation genetics in transition to conservation genomics. Trends Genet. 2010, 26: 177-187. 10.1016/j.tig.2010.01.001.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Morin PA, Luikart G, Wayne RK, the SNP workshop group: SNPs in ecology, evolution and conservation. Trends Ecol Evo. 2004, 19: 208-216. 10.1016/j.tree.2004.01.009.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Garrigan D, Hedrick PW: Perspective: detecting adaptive molecular polymorphism: lessons from the MHC. Evolution. 2003, 57: 1707-1722.PubMedView ArticleGoogle Scholar
- Andersson L, Mikko S: Generation of MHC class II diversity by intra- and intergenic recombination. Immunol Rev. 1995, 143: 5-12. 10.1111/j.1600-065X.1995.tb00667.x.PubMedView ArticleGoogle Scholar
- Bergstrom TF, Josefsson A, Erlich HA, Gyllensten U: Recent origin of HLA-DRB1 alleles and implications for human evolution. Nat Genet. 1998, 18: 237-242. 10.1038/ng0398-237.PubMedView ArticleGoogle Scholar
- Van Oosterhout C, Joyce DA, Cummings SM, Blais J, Barson NJ, Ramnarine IW, Mohammed RS, Persad N, Cable J: Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata). Evolution. 2006, 60: 2562-2574.PubMedView ArticleGoogle Scholar
- Sutton JT, Nakagawa S, Robertson BC, Jamieson IG: Disentangling the roles of natural selection and genetic drift in shaping variation at MHC immunity genes. Mol Ecol. 2011, 20: 4408-4420. 10.1111/j.1365-294X.2011.05292.x.PubMedView ArticleGoogle Scholar
- Miller KM, Kaukinen KH, Beacham TD, Withler RE: Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon. Genetica. 2001, 111: 237-257. 10.1023/A:1013716020351.PubMedView ArticleGoogle Scholar
- Hedrick PW, Lee RN, Buchanan C: Canine parvovirus enteritis, canine distemper, and major histocompatibility complex genetic variation in Mexican wolves. J Wildlife Dis. 2003, 39: 909-913.View ArticleGoogle Scholar
- Castro-Prieto A, Wachter B, Melzheimer J, Thalwitzer S, Sommer S: Diversity and evolutionary patterns of immune genes in free-ranging Namibian leopards (Panthera pardus pardus). Heredity. 2011, 102: 653-665. 10.1093/jhered/esr097.View ArticleGoogle Scholar
- Mikko S, Roed K, Schmutz S, Andersson L: Monomorphism and polymorphism at MHC DRB loci in domestic and wild ruminants. Immunol Rev. 1999, 167: 169-178. 10.1111/j.1600-065X.1999.tb01390.x.PubMedView ArticleGoogle Scholar
- Weber DS, Stewart BS, Schienman J, Lehman N: Major histocompatibility complex variation at three class II loci in the northern elephant seal. Mol Ecol. 2004, 13: 711-718. 10.1111/j.1365-294X.2004.02095.x.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Mason RAB, Browning TL, Eldridge DB: Reduced MHC class II diversity in island compared to mainland populations of the black-footed rock-wallaby (Petrogale lateralis lateralis). Conserv Genet. 2011, 12: 91-103. 10.1007/s10592-009-9993-y.View ArticleGoogle Scholar
- Biedrzycka A, Radwan J: Population fragmentation and major histocompatibility complex variation in the spotted suslik, Spermophilus suslicus. Mol Ecol. 2008, 17: 4801-4811. 10.1111/j.1365-294X.2008.03955.x.PubMedView ArticleGoogle Scholar
- Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure. Proc R Soc B. 1996, 263: 1619-1626. 10.1098/rspb.1996.0237.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Aguilar A, Garza JC: A comparison of variability and population structure for major histocompatibiltiy complex and microsatellite loci in California coastal steelhead (Oncorhynchus mykiss Walbaum). Mol Ecol. 2006, 15: 923-937. 10.1111/j.1365-294X.2006.02843.x.PubMedView ArticleGoogle Scholar
- Zhang RZ: The primates of China: biogeography and conservation status—past, present and future. 2002, Beijing: China Forestry Publishing HouseGoogle Scholar
- Li M, Liu ZJ, Gou JX, Ren BP, Pan RL, Su YJ, Funk SM, Wei FW: Phylogeography and population structure of the golden monkeys (Rhinopithecus roxellana): inferred from mitochondrial DNA sequences. Am J Primatol. 2007, 69: 1195-1209. 10.1002/ajp.20425.PubMedView ArticleGoogle Scholar
- Pan D, Hu HX, Meng SJ, Men ZM, Fu YX, Zhang YP: A high polymorphism level in Rhinopithecus roxellana. Int J Primatol. 2009, 30: 337-351. 10.1007/s10764-009-9347-0.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Boessenkool S, Taylor SS, Tepolt CK, Komdeur J, Jamieson IG: Large mainland populations of South Island robins retain greater genetic diversity than offshore island refuges. Conserv Genet. 2007, 8: 705-714. 10.1007/s10592-006-9219-5.View ArticleGoogle Scholar
- Hedrick PW, Parker KM, Miller EL, Miller PS: Major histocompatiblity complex variation in the endangered Przewalski's horse. Genetics. 2001, 152: 1701-1710.Google Scholar
- 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.PubMedGoogle Scholar
- Chang ZF, Liu ZJ, Yang JY, Li M, Vigilant L: Noninvasive genetic assessment of the population trend and sex ratio of the Shennongjia population of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Chinese Sci Bull. 2011, 57: 1135-1141.View ArticleGoogle Scholar
- Luo MF, Liu ZJ, Pan HJ, Zhao L, Li M: Historical geographic dispersal of the golden snub-nosed monkey (Rhinopithecus roxellana) and the influence of climatic oscillations. Am J Primatol. 2012, 74: 91-101. 10.1002/ajp.21006.PubMedView ArticleGoogle Scholar
- Bijlsma R, Bundgaard J, Boerema AC: Does inbreeding affect the extinction risk of small populations? Predictions from Drosophila. J Evolution Biol. 2000, 13: 502-514. 10.1046/j.1420-9101.2000.00177.x.View ArticleGoogle Scholar
- Thoß M, Ilmonen P, Musolf K, Penn DJ: Major histocompatibility complex heterozygosity enhances reproductive success. Mol Ecol. 2011, 20: 1546-1557. 10.1111/j.1365-294X.2011.05009.x.PubMedView ArticleGoogle Scholar
- Lukas D, Bradley BJ, Nsubuga AM, Doran-Sheehy D, Robbins MM, Vigilant L: Major histocompatibility complex and microsatellite variation in two populations of wild gorillas. Mol Ecol. 2004, 13: 3389-3402. 10.1111/j.1365-294X.2004.02353.x.PubMedView ArticleGoogle Scholar
- Goda1 N, Mano T, Kosintsev P, Vorobie A, Masuda R: Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae. Tissue Antigens. 2010, 76: 404-410. 10.1111/j.1399-0039.2010.01528.x.View ArticleGoogle Scholar
- Kennedy LJ, Randall DA, Knobel D, Brown JJ, Fooks AR, Argaw K, Shiferaw F, Ollier WER, Sillero-Zubiri C, Macdonald DW, Laurenson MK: Major histocompatibility complex diversity in the endangered Ethiopian wolf (Canis simensis). Tissue Antigens. 2010, 77: 118-125.View ArticleGoogle Scholar
- Radwan J, Demiaszkiewicz AW, Kowalczyk R, Lachowicz J, Kawalko A, Wojcik JM, Pyziel AM, Babik W: An evaluation of two potential risk factors, MHC diversity and host density, for infection by an invasive nematode Ashworthius sidemi in endangered European bison (Bison bonasus). Biol Conserv. 2010, 143: 2049-2053. 10.1016/j.biocon.2010.05.012.View ArticleGoogle Scholar
- Pokorny I, Sharma R, Goyal SP, Mishra S, Tiedemann R: MHC class I and MHC class II DRB gene variability in wild and captive Bengal tigers (Panthera tigris tigris). Immunogenetics. 2010, 62: 667-679. 10.1007/s00251-010-0475-7.PubMedView ArticleGoogle Scholar
- Kimura M: Preponderance of synonymous changes as evidence for the neutrality theory of molecular evolution. Nature. 1977, 267: 275-276. 10.1038/267275a0.PubMedView ArticleGoogle Scholar
- Hughes AL, Nei M: Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988, 335: 167-170. 10.1038/335167a0.PubMedView ArticleGoogle Scholar
- Hughes AL, Nei M: Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection. Proc Natl Acad Sci USA. 1989, 86: 958-962. 10.1073/pnas.86.3.958.PubMedPubMed CentralView ArticleGoogle Scholar
- Hughes AL: Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level. Heredity. 2007, 99: 364-373. 10.1038/sj.hdy.6801031.PubMedView ArticleGoogle Scholar
- Klein J, Sato A, Nikolaidis N: MHC, TSP, and the origin of species: From immunogenetics to evolutionary genetics. Annu Rev Genet. 2007, 41: 281-304. 10.1146/annurev.genet.41.110306.130137.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Li HS, Ligons DL, Rose NR: Genetic complexity of autoimmune myocarditis. Autoimm reviews. 2008, 7: 168-173. 10.1016/j.autrev.2007.11.010.View ArticleGoogle Scholar
- Schaschl H, Suchentruk F, Hammer S, Goodman SJ: Recombination and the origin of sequence diversity in the DRB MHC class II locus in chamois (Rupicapra spp.). Immunogenetics. 2005, 57: 108-115. 10.1007/s00251-005-0784-4.PubMedView ArticleGoogle Scholar
- Xu SX, Ren WH, Zhou XM, Zhou KY, Yang G: Sequence polymorphism and geographical variation at a positively selected MHC-DRB gene in the finless porpoise (Neophocaena phocaenoides): Implication for recent differentiation of the Yangtze finless porpoise?. J Mol Ecol. 2010, 71: 6-22. 10.1007/s00239-010-9357-8.View ArticleGoogle Scholar
- Kauppi L, Jeffreys AJ, Keeney S: Where the crossovers are: recombination distributions in mammals. Nat Rev Genet. 2004, 5: 413-424. 10.1038/nrg1346.PubMedView ArticleGoogle Scholar
- Bonneaud C, Chastel O, Federici P, Westerdahl , Sorci G: Complex Mhc-based mate choice in a wild passerine. Proc R Soc B. 2006, 273: 1111-1116. 10.1098/rspb.2005.3325.PubMedPubMed CentralView ArticleGoogle Scholar
- Muirhead CA: Consequences of population structure on genes under balancing selection. Evolution. 2001, 55: 1532-1541.PubMedView ArticleGoogle Scholar
- Richman AD, Herrera LG, Nash D: Evolution of MHC class II E beta diversity within the genus Peromysus. Genetics. 2003, 164: 289-297.PubMedPubMed CentralView ArticleGoogle Scholar
- Miller HC, Lambert DM: Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae). Mol Ecol. 2004, 13: 3709-3721. 10.1111/j.1365-294X.2004.02368.x.PubMedView ArticleGoogle Scholar
- Kimura M: The Neutral Theory of Molecular Evolution. 1983, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Maureen BP, Thomas FT: Genetic variation of the major histocompatibility complex (MHC class II b gene) in the threatened Gila trout, Oncorhynchus gilae gilae. Conserv Genet. 2008, 9: 257-270. 10.1007/s10592-007-9336-9.View ArticleGoogle Scholar
- Agudo R, Alcaide M, Rico C, Lemus JA, Blanco G, Hiraldo F, Donázar JA: Major histocompatibility complex varian in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection. Mol Ecol. 2011, 20: 2329-2340. 10.1111/j.1365-294X.2011.05107.x.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Ujvari B, Belov K: Major histocompatibility complex (MHC) markers in conservation biology. Int J Mol Sci. 2011, 12: 5168-5186. 10.3390/ijms12085168.PubMedPubMed CentralView ArticleGoogle Scholar
- Tallmon DA, Luikart G, Waples RS: The alluring simplicity and complex reality of genetic rescue. Trends Ecol Evol. 2004, 19: 489-496. 10.1016/j.tree.2004.07.003.PubMedView ArticleGoogle Scholar
- Edmands S: Between a rock and a hard place: Evaluating the relative risks of inbreeding and outbreeding for conservation and management. Mol Ecol. 2007, 16: 463-475.PubMedView ArticleGoogle Scholar
- Vilà C, Sundqvist AK, Flagstad O, Seddon J, Björnerfeldt S, Kojola I, Casulli A, Sand H, Wabakken P, Ellegren H: Rescue of a severely bottlenecked wolf (Canis lupus) population by a single immigrant. Proc R Soc Lond B. 2003, 270: 91-97. 10.1098/rspb.2002.2184.View ArticleGoogle Scholar
- Madsen T, Ujvari B, Olsson M: Novel genes continue to enhance population growth in adders (Vipera berus). Biol Conserv. 2004, 120: 145-147. 10.1016/j.biocon.2004.01.022.View ArticleGoogle Scholar
- Bouzat JL, Johnson JE, Toepfer JE, Simpson SA, Esker TL, Westemeier RL: Beyond the beneficial effects of translocations in an effective tool for the genetic restoration of isolated populations. Conserv Genet. 2009, 10: 191-201. 10.1007/s10592-008-9547-8.View ArticleGoogle Scholar
- Doxiadis GGM, Otting N, de Groot NG, Groot N, Rouweler AJM, Noort R, Verschoor EJ, Bontjer I, Bontrop RE: Evolutionary stability of MHC class II haplotypes in diverse rhesus macaque populations. Immunogenetics. 2003, 55: 540-551. 10.1007/s00251-003-0590-9.PubMedView ArticleGoogle Scholar
- Nei M: Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978, 89: 583-590.PubMedPubMed CentralGoogle Scholar
- Excoffier L, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform. 2005, 1: 47-50.Google Scholar
- Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). 2001, Available from http://www.unil.ch/izea/softwares/fstat.html. Updated from Goudet (1995).Google Scholar
- Rousset F: Genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Res. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.View ArticleGoogle Scholar
- Weir BS: Inferences about linkage disequilibrium. Biometrics. 1979, 35: 235-254. 10.2307/2529947.PubMedView ArticleGoogle Scholar
- Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA 4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Kennedy LJ, Ryvar R, Gaskell RM, Addie D, Willoughby K, Carter S, Thomson W, Ollier W, Radford A: Sequence analysis of MHC DRB alleles in domestic cats from the United Kingdom. Immunogenetics. 2002, 54: 348-352. 10.1007/s00251-002-0465-5.PubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Falush D, Stephens M, Pritchard J: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Guindon S, Gascuel O: A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA: Modeltest: Testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48. 10.1093/oxfordjournals.molbev.a026036.PubMedView ArticleGoogle Scholar
- Hudson RR: Two-locus sampling distributions and their application. Genetics. 2001, 159: 1805-1817.PubMedPubMed CentralGoogle Scholar
- McVean G, Awadalla P, Fearnhead P: A coalescennt-based method for detecting and estimating recombination from gene sequences. Genetics. 2002, 160: 1231-1241.PubMedPubMed CentralGoogle Scholar
- Richman AD, Herrera LG, Nash D, Schierup MH: Relative roles of mutation and recombination in generating allelic polymorphism at an MHC class II locus in Peromyscus maniculatus. Genet Res. 2003, 82: 89-99. 10.1017/S0016672303006347.PubMedView ArticleGoogle Scholar
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and non-synonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMedGoogle Scholar
- Reche PA, Reinherz EL: Sequence variability analysis of human class I and class II MHC molecules: functional and structural correlates of amino acid polymorphisms. J Mol Biol. 2003, 331: 623-641. 10.1016/S0022-2836(03)00750-2.PubMedView ArticleGoogle Scholar
- Brown JH, Jardetzky TS, Gorga JC, Stern LJ, Urban RG, Strominger JL, Wiley DC: 3-Dimensional structure of the human class-II histocompatibility antigen HLA-DR1. Nature. 1993, 364: 33-39. 10.1038/364033a0.PubMedView ArticleGoogle Scholar
- Kuduk K, Johanet A, Allaine D, Cohas A, Radwan J: Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). J Evolutiona Biol. 2012, 10.1111/j.1420-9101.2012.02537.x.Google Scholar
- Kloch A, Baran K, Buczek M, Konarzewski M, Radwan J: MHC influences infection with parasites and wintersurvival in the root vole Microtus oeconomus. Evol Ecol. 2012, 10.1007/s10682-012-9611-1.Google Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-736.PubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AMK: Codon substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.PubMedPubMed CentralGoogle Scholar
- Yang Z, Wong WSW, Nielsen R: Bayes Empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1111. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
- Antao T, Lopes A, Lopes RJ, Beja-Pereira A, Luikart G: LOSITAN: a workbench to detect molecular adaptation based on a Fst-outlier method. BMC Bioinformatics. 2008, 9: 323-10.1186/1471-2105-9-323.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.