Sequence variation at the major histocompatibility complex locus DQ beta in beluga whales (Delphinapterus leucas) [published erratum appears in Mol Biol Evol 1995 Nov;12(6):1174]

Genetic variation at the Major Histocompatibility Complex locus DQ beta was analyzed in 233 beluga whales (Delphinapterus leucas) from seven populations: St. Lawrence Estuary, eastern Beaufort Sea, eastern Chukchi Sea, western Hudson Bay, eastern Hudson Bay, southeastern Baffin Island, and High Arctic and in 12 narwhals (Monodon monoceros) sympatric with the High Arctic beluga population. Variation was assessed by amplification of the exon coding for the peptide binding region via the polymerase chain reaction, followed by either cloning and DNA sequencing or single-stranded conformation polymorphism analysis. Five alleles were found across the beluga populations and one in the narwhal. Pairwise comparisons of these alleles showed a 5:1 ratio of nonsynonymous to synonymous substitutions per site leading to eight amino acid differences, five of which were nonconservative substitutions, centered around positions previously shown to be important for peptide binding. Although the amount of allelic variation is low when compared with terrestrial mammals, the nature of the substitutions in the peptide binding sites indicates an important role for the DQ beta locus in the cellular immune response of beluga whales. Comparisons of allele frequencies among populations show the High Arctic population to be different (P < or = .005) from the other beluga populations surveyed. In these other populations an allele, Dele-DQ beta*0101-2, was found in 98% of the animals, while in the High Arctic it was found in only 52% of the animals. Two other alleles were found at high frequencies in the High Arctic population, one being very similar to the single allele found in narwhal.      

Allelic and haplotype variation of major histocompatibility complex class II DRB1 and DQB loci in the St Lawrence beluga (Delphinapterus leucas)

In order to assess levels of major histocompatibility complex (Mhc) variation within the St Lawrence beluga (Delphinapterus leucas) the variation at the beluga Mhc DRB1 class II locus was assessed by single-strand conformation polymorphism (SSCP) analysis of the peptide-binding region for 313 whales collected from 13 sampling locations across North America. In addition, samples from west Greenland and the St Lawrence were also typed at the DQB locus, allowing comparison to a previous study and assessment of linkage disequilibrium of alleles at the two loci. Comparisons of DRB1 and DQB allele frequencies among all sampling locations indicated genetic structure (α < 0.005). Most of this structure resulted from differences between the different wintering groups. Significant genetic structure (α = 0.05) exists among each pair of the following groups at both the DRB1 and DQB loci; St Lawrence, Hudson Strait, Bering Sea, Cunningham Inlet, and Davis Strait (minus Cunningham Inlet), except the St Lawrence and Hudson Strait for the DQB locus. In the St Lawrence population, six of the eight DRB1 alleles are present representing all five known allelic lineages. Evidence of linkage disequilibrium between the DRB1 and DQB is present in two sampling locations, the St Lawrence and Nuussuaq (α = 0.05). Analysis of probable DRB1–DQB haplotypes among groups of beluga suggests a haplotype reduction in the St Lawrence.     

Characterization of Mhc-DRB allelic diversity in white-tailed deer (Odocoileus virginianus) provides insight into Mhc-DRB allelic evolution within Cervidae

 Although white-tailed deer (Odocoileus virginianus) are one of North America's best studied mammals, no information is available concerning allelic diversity at any locus of the major histocompatibility complex in this taxon. Using the polymerase chain reaction, single-stranded conformation polymorphism analysis, and DNA sequencing techniques, 15 DRB exon 2 alleles were identified among 150 white-tailed deer from a single population in southeastern Oklahoma. These alleles represent a single locus and exhibit a high degree of nucleotide and amino acid polymorphism, with most amino acid variation occurring at positions forming the peptide binding sites. Furthermore, twenty-seven amino acid residues unique to white-tailed deer DRB alleles were detected, with 19 of these occurring at residues forming contact points of the peptide binding region. Significantly higher rates of nonsynonymous than synonymous substitutions were detected among these DRB alleles. In contrast to other studies of Artiodactyla DRB sequences, interallelic recombination does not appear to be playing a significant role in the generation of allelic diversity at this locus in white-tailed deer. To examine evolution of white-tailed deer (Odvi-DRB) alleles within Cervidae, we performed a phylogenetic analysis of all published red deer (Ceel-DRB), roe deer (Caca-DRB), and moose (Alal-DRB) DRB alleles. The phylogenetic tree clearly shows a trans-species persistence of DRB lineages among these taxa. Moreover, this phylogenetic tree provides insight into evolution of DRB allelic lineages within Cervidae and may aid in assignment of red deer DRB alleles to specific loci. Received: 25 June 1998 / Revised: 2 September 1998     

Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex, Capra pyrenaica

During the last two centuries, the Spanish ibex (Capra pyrenaica) has shown a significant demographic decline as a result of the progressive destruction of its natural habitat, disease epidemics, and uncontrolled hunting. Partial sequencing of the class II MHC DRB1 gene revealed that the Spanish ibex has remarkably low levels of genetic variation at this locus, with only six different DRB1 alleles and an observed heterozygosity of 0.429-0.579. The rates of nonsynonymous vs synonymous substitutions were significantly different in the peptide-binding region (dN/dS=5.347, P=0.002), a feature that indicates that the DRB1 gene is under positive selection. A phylogenetic analysis of the Spanish ibex and a set of domestic goat DRB1 alleles revealed that the reported sequences represent four major allelic lineages. The limited allelic repertoire of the DRB1 gene in the Spanish ibex is likely the direct result of the recent history of population bottlenecks and marked demographic decline of this species. A genetic survey of 13 microsatellite loci was consistent with this idea. The Spanish ibex subspecies C. p. hispanica and C. p. victoriae consistently showed considerably lower levels of microsatellite heterozygosity (Ho=0.184-0.231) and allelic diversity (mean number of alleles per locus=2-2.4) than those reported in other wild ruminants. This study demonstrates the significance of both natural selection and the demographic history of populations in determining patterns of genetic variation at MHC loci. In addition, our results emphasize the importance of locally adapted populations for the preservation of genetic diversity.     

Sequence variation and gene duplication at MHC DQB loci of baiji (Lipotes vexillifer), a Chinese river dolphin

The major histocompatibility complex (MHC) is a fundamental part of the vertebrate immune system, and the high variability in many MHC genes is thought to play an important role in the recognition of parasites. Baiji (Lipotes vexillifer) is one of the most endangered species in the world. Its wild population has declined to fewer than 100 individuals and has a very high risk of becoming extinct in the near future. In this study we present a first step in the molecular characterization of a DQB-like locus of baiji by nucleotide sequence analysis of the polymorphic exon 2 segments. In the examined 172 bp sequences from a group of 18 incidentally captured or stranded individuals, 48 variable sites were determined and 43 alleles were identified, many of which were represented by only one clone. Three to seven alleles were found in each individual, suggesting gene duplications. No deletion, insertion, or exceptional stop codon was detected, suggesting these alleles function in vivo. Phylogenetic reconstruction using neighbor joining grouped the 43 alleles into two distinct lineages, differing by seven nucleotides and four amino acids. Substitutions of amino acids tend to be clustered around sites postulated to be responsible for selective peptide recognition. In the peptide-binding region (PBR) of the DQB locus, the average number of nonsynonymous substitutions per site is greater than that of synonymous substitutions per site (0.1962 versus 0.0256, respectively). Nucleotide and amino acid sequences both showed a relatively high level of similarity (nucleotides 90.6%; amino acids 80.6%) to those of beluga whale (Delphinapterus leucas) and narwhal (Monodon monoceros). The high level of baiji MHC polymorphism revealed in the present study has not been reported in other cetaceans and could be a consequence of the small baiji population adapting to freshwater with a relatively high level of pathogens.     

Analysis of polymorphisms in the major expressed class I locus (B-FIV) of the chicken

 We analyzed the polymorphic nature of eleven alleles expressed by the major class I locus (B-FIV) in chickens. Similar to mammalian class I loci, the nucleotide substitutions with high variability occur in exons 2 and 3 encoding the α1 and α2 domains. However, the nonsynonymous to synonymous ratio of nucleotide substitutions in exon 3 encoding the α helix and β sheets is reversed compared with HLA. The region of exon 3 encoding the α2 helix demonstrates a much lower nonsynonymous to synonymous ratio, suggesting evolutionary selection of a more conserved α2 helix in B-FIV compared with HLA. Amino acid residues with high Wu-Kabat variability are typically located in positions predicted to impact antigen presentation. B-FIV amino acid residues predicted to interact with the CDR1α region of the T-cell receptor (Tcr) demonstrate less variability than in mouse and human class I alleles. The combination of a reduced nonsynonymous to synonymous ratio in exon 3 encoding the α2 helix and the limited variability in CDR1α contact residues is discussed with regard to concerted evolution between a minimal major histocompatibility complex and compaction of Tcr variable gene segments in the chicken. Received: 18 Juli 1997 /  Revised: 21 November 1997     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

Diversifying Selection Governs Sequence Polymorphism in the Major Adhesin Proteins FimA, PapA, and SfaA of Escherichia coli

Fimbriae or pili are essential adherence factors usually found in pathogenic bacteria to aid colonization of host cells. Three major structural pilin genes, fimA, sfaA, and papA, from Escherichia coli natural isolates were examined and nucleotide sequence data revealed elevated levels of both synonymous and nonsynonymous site variation at these loci. Examination of synonymous site variation shows a fivefold increase in fimA sites, relative to the housekeeping gene mdh; and similarly the sfaA and papA genes have increased synonymous sites variation relative to fimA. Nonsynonymous site variation is also elevated at all three loci but, in particular, at the papA locus (k _N= 0.44). The k _N/k _S ratio for the three genes are among the highest yet reported for E. coli genes. Regional variation in nucleotide polymorphism within each of the genes reveal hypervariable segments where nonsynonymous substitutions exceed synonymous substitutions. We propose that at the fimA, papA, and sfaA genes, diversifying selection has brought about the increase levels of polymorphism. Received: 7 August 1997 / Accepted: 8 March 1998     

Molecular considerations in the evolution of bacterial genes

Summary Synonymous and nonsynonymous substitution rates at the loci encoding glyceraldehyde-3-phosphate dehydrogenase (gap) and outer membrane protein 3A (ompA) were examined in 12 species of enteric bacteria. By examining homologous sequences in species of varying degrees of relatedness and of known phylogenetic relationships, we analyzed the patterns of synonymous and nonsynonymous substitutions within and among these genes. Although both loci accumulate synonymous substitutions at reduced rates due to codon usage bias, portions of thegap andompA reading frames show significant deviation in synonymous substitution rates not attributable to local codon bias. A paucity of synonymous substitutions in portions of theompA gene may reflect selection for a novel mRNA secondary structure. In addition, these studies allow comparisons of homologous protein-coding sequences (gap) in plants, animals, and bacteria, revealing differences in evolutionary constraints on this glycolytic enzyme in these lineages.     

Major histocompatibility complex class II genetic variation in forest musk deer (Moschus berezovskii) in China

The major histocompatibility complex (MHC) plays an important role in the immune system of vertebrates. We used the second exon of four MHC class II genes (DRA, DQA1, DQA2 and DRB3) to assess the overall MHC variation in forest musk deer (Moschus berezovskii). We also compared the MHC variation in captive and wild populations. We observed 22 alleles at four loci (four at DRA, four at DQA1, four at DQA2 and 10 at DRB3), 15 of which were newly identified alleles. Results suggest that forest musk deer maintain relatively high MHC variation, which may result from balancing selection. Moreover, considerable diversity was observed at the DRA locus. We found a high frequency of Mobe‐DRA*02, Mobe‐DQA1*01 and Mobe‐DQA2*05 alleles, which may be important for pathogen resistance. A Ewens–Watterson test showed that the DRB3 locus in the wild population had experienced recent balancing selection. We detected a small divergence at the DRA locus, suggesting the effect of weak positive selection on the DRA gene. Alternatively, this locus may be young and not yet adapted a wide spectrum of alleles for pathogen resistance. The significant heterozygosity deficit observed at the DQA1 and DRB3 loci in the captive population and at all four loci in the wild population may be the result of a population bottleneck. Additionally, MHC genetic diversity was higher in the wild population than in the captive, suggesting that the wild population may have the ability to respond to a wider range of pathogens.     

Low synonymous site variation at thelacY locus inEscherichia coli suggests the action of positive selection

We have determined the nucleotide sequences of sevenlacY alleles isolated from natural isolates ofEscherichia coli. Nucleotide heterozygosity estimates for this locus were compared to those obtained from previous studies of intraspecific variation at chromosomal loci, revealing thatlacY has unusually low synonymous site variation. The average pairwise heterozygosity of synonymous sites (Ks=0.0112+/-0.0100) is the second lowest reported and the lowest for loci that have an equivalent level of nonsynonymous variation. We consider several hypotheses to explain how different forces in evolution could act to create the observed pattern of polymorphism, including selection for translational efficiency and positive selection. Our analysis most strongly supports the hypothesis that positive selection has acted on thelacY locus inE. coli.     

Evolution of MHC class II SLA‐DRB1 locus in the Croatian wild boar (Sus scrofa) implies duplication and weak signals of positive selection

The wild boar is an ancestor of the domestic pig and an important game species with the widest geographical range of all ungulates. Although a large amount of data are available on major histocompatibility complex (MHC) variability in domestic pigs, only a few studies have been performed on wild boars. Due to their crucial role in appropriate immune responses and extreme polymorphism, MHC genes represent some of the best candidates for studying the processes of adaptive evolution. Here, we present the results on the variability and evolution of the entire MHC class II SLA‐DRB1 locus exon 2 in 133 wild boars from Croatia. Using direct sequencing and cloning methods, we identified 20 SLA‐DRB1 alleles, including eight new variants, with notable divergence. In some individuals, we documented functional locus duplication, and SLA‐DRB1*04:10 was identified as the allele involved in the duplication. The expression of a duplicated locus was confirmed by cloning and sequencing cDNA‐derived amplicons. Based on individual genotypes, we were able to assume that alleles SLA‐DRB1*04:10 and SLA‐DRB1*06:07 are linked as an allelic combination that co‐evolves as a two‐locus haplotype. Our investigation of evolutionary processes at the SLA‐DRB1 locus confirmed the role of intralocus recombination in generating allelic variability, whereas tests of positive selection based on the dN/dS (non‐synonymous/synonymous substitution rate ratio) test revealed atypically weak and ambiguous signals.     

Positive Selection in the N-Terminal Extramembrane Domain of Lung Surfactant Protein C (SP-C) in Marine Mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (ω) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Richard Kliman     

Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles

Seventy-five individuals of Bubalus bubalis belonging to four different breeds, three of river buffalo and one of swamp buffalo, were studied for polymorphism in MHC DRB (Bubu-DRB) and DRA (Bubu-DRA) loci. Eight alleles of Bubu-DRB were found, and all alleles in the swamp type were shared with the three river breeds. All alleles sampled from the breed of European origin (Mediterranean) were present in breeds sampled in Brazil, thus variability of this locus may have been preserved to a great extent in the more recently founded Brazilian population. Bubu-DRB alleles contained higher proportions of synonymous vs. non-synonymous substitutions in the non-peptide-binding sites (PBS) region, in contrast to the pattern of variation found in BoLA-DRB3, the orthologous locus in cattle. This indicated that either the first domain exon (exon 2) of Bubu-DRB has not undergone as much recombination and/or gene conversion as in cattle alleles, or Bubu-DRB may be more ancient than BoLA-DRB3 alleles. Phylogenetic analysis of DRB alleles from Bubalus, Syncerus c. caffer, the Cape buffalo, and domestic cattle demonstrated transspecies polymorphism. Water buffalo contained two alleles of DRA that differed from each other in two amino acid positions, including one in the PBS (alpha22) that was also shared with Anoa depressicornis, the anoa. Discovery of variation in DRA was surprising as the first domain of DRA is a highly conserved polypeptide in mammals in general and especially in ruminants, where no other substitution in PBS was seen.     

Allelic diversity at the Mhc-DQA locus in cotton rats (Sigmodon hispidus) and a comparison of DQA sequences within the family Muridae (Mammalia: Rodentia)

 The cotton rat (Sigmodon hispidus) is a common murid rodent of the southern United States, Mexico, and Central America. Using single-stranded conformation polymorphism analysis and DNA sequencing techniques, 11 DQA exon 2 alleles were detected among 180 S. hispidus from Caddo County, Oklahoma, USA. The alleles represent a single locus exhibiting a high level of polymorphism. Nucleotide and amino acid distance values among DQA alleles of S. hispidus were higher than those within Mus musculus and species of Rattus. Although the distribution of polymorphic amino acid residues among alleles of S. hispidus was similiar to that of Mus and Rattus, some residues of the α-helix region were more variable in S. hispidus. Comparisons of nonsynonymous and synonymous substitutions indicated a trend toward higher numbers of nonsynonymous substitutions; however, this difference was not significant statistically among S. hispidus alleles. To examine evolution of DQA alleleswithin Muridae, we performed a phylogenetic analysis that included DQA alleles from S. hispidus, Peromyscus leucopus, M. musculus, R. norvegicus, and six Australian species of Rattus. Results depicted monophyly for each genus, and this concordance between species and gene trees represents a lack of evidence for trans-species persistence of alleles among these genera. Received: 22 October 1998 / Revised: 17 March 1999     

Sequence variability analysis on major histocompatibility complex class II DRB alleles in three felines

The variation of the exon 2 of the major histo-compatibility complex (MHC) class II gene DRB locus in three feline species were examined on clouded leopard (Neofelis nebulosa), leopard (Panthera pardus) and Amur tiger (Panthera tigris altaica). A pair of degenerated primers was used to amplify DRB locus covering almost the whole exon 2. Exon 2 encodes the β1 domain which is the most variable fragments of the MHC class II molecule. Single-strand conformational polymorphism (SSCP) analysis was applied to detect different MHC class II DRB haplotypes. Fifteen recombinant plasmids for each individual were screened out, isolated, purified and sequenced finally. Totally eight distinct haplotypes of exon 2 were obtained in four individuals. Within 237 bp nucleotide sequences from four samples, 30 variable positions were found, and 21 putative peptide-binding positions were disclosed in 79 amino acid residues. The ratio of nonsynonymous substitutions (d _N ) was much higher than that of synonymous substitutions (d _S ), which indicated that balancing selection probably maintain the variation of exon 2. MEGA neighbor joining (NJ) and PAUP maximum parsimony (MP) methods were used to reconstruct phylogenetic trees among species, respectively. Results displayed a more close relationship between leopard and tiger; however, clouded leopard has a comparatively distant relationship form the other two. __________ Translated from Zoological Research, 2006, 27(2): 181-C188 [译自:动物学研究]     

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.     

Major histocompatibility complex variation and evolution at a single, expressed DQA locus in two genera of elephants

Genes of the vertebrate major histocompatibility complex (MHC) are crucial to defense against infectious disease, provide an important measure of functional genetic diversity, and have been implicated in mate choice and kin recognition. As a result, MHC loci have been characterized for a number of vertebrate species, especially mammals; however, elephants are a notable exception. Our study is the first to characterize patterns of genetic diversity and natural selection in the elephant MHC. We did so using DNA sequences from a single, expressed DQA locus in elephants. We characterized six alleles in 30 African elephants (Loxodonta africana) and four alleles in three Asian elephants (Elephas maximus). In addition, for two of the African alleles and three of the Asian alleles, we characterized complete coding sequences (exons 1–5) and nearly complete non-coding sequences (introns 2–4) for the class II DQA loci. Compared to DQA in other wild mammals, we found moderate polymorphism and allelic diversity and similar patterns of selection; patterns of non-synonymous and synonymous substitutions were consistent with balancing selection acting on the peptides involved in antigen binding in the second exon. In addition, balancing selection has led to strong trans-species allelism that has maintained multiple allelic lineages across both genera of extant elephants for at least 6 million years. We discuss our results in the context of MHC diversity in other mammals and patterns of evolution in elephants.     

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1–2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.