Host-Parasite Interactions and Purifying Selection in a Microsporidian Parasite of Honey Bees

To clarify the mechanisms of Nosema ceranae parasitism, we deep-sequenced both honey bee host and parasite mRNAs throughout a complete 6-day infection cycle. By time-series analysis, 1122 parasite genes were significantly differently expressed during the reproduction cycle, clustering into 4 expression patterns. We found reactive mitochondrial oxygen species modulator 1 of the host to be significantly down regulated during the entire infection period. Our data support the hypothesis that apoptosis of honey bee cells was suppressed during infection. We further analyzed genome-wide genetic diversity of this parasite by comparing samples collected from the same site in 2007 and 2013. The number of SNP positions per gene and the proportion of non-synonymous substitutions per gene were significantly reduced over this time period, suggesting purifying selection on the parasite genome and supporting the hypothesis that a subset of N. ceranae strains might be dominating infection.     

Sequence diversity of five <Emphasis Type="Italic">Medicago sativa</Emphasis> genes involved in agronomic traits to set up allele mining in breeding

Lucerne (Medicago sativa) is an autotetraploid forage legume, whose breeding mainly relies on phenotypic recurrent selection to create synthetic populations. Allele mining could be a strategy to take advantage of the large genetic diversity of the species. This paper reports the polymorphism of five genes impacting the phenotype for selected traits: CAD1 and CCoaOMT (digestibility), CONSTANS-like (forage yield), NHX1 (salt tolerance), and WXP1 (drought tolerance). Complete genes were sequenced for 387 genotypes of 43 cultivated accessions and 20 genotypes of wild accessions. Wild versus cultivated polymorphism were compared, adaptive evolution was evaluated by comparing M. sativa and Medicago truncatula sequences and variants of the cultivated pool were characterized. We showed that the wild pool was more variable than the cultivated pool, with 36.6% and 8.4% of the variants that were specific to the wild pool and the cultivated pool, respectively. This result confirmed a bottleneck effect during domestication and selection. We also found that CAD1, CCoaOMT, and NHX1 were under a strong purifying selection and contained few non-synonymous variants (8, 5, and 8, respectively), while CONSTANS-like and WXP1 were under a less pronounced purifying selection and were more polymorphic (45 and 91 non-synonymous variants, respectively). This result suggests that adaptive evolution could be an indicator of expected polymorphism. The application of allele mining strategy for plant breeding in autotetraploid species is discussed.     

The Role of Selection in Shaping Diversity of Natural M. tuberculosis Populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to −9.5×10⁻⁴ (95% CI −1.1×10⁻³ to −6.8×10⁻⁴); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620–1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

The axon degeneration gene SARM1 is evolutionarily distinct from other TIR domain-containing proteins

Many forms of neurodegenerative disease are characterized by Wallerian degeneration, an active program of axonal destruction. Recently, the important player which enacts Wallerian degeneration was discovered, the multidomain protein SARM1. Since the SARM1 protein has classically been thought of as an innate immune molecule, its role in Wallerian degeneration has raised questions on the evolutionary forces acting on it. Here, we synthesize a picture of SARM1’s evolution through various organisms by examining the molecular and genetic changes of SARM1 and the genes around it. Using proteins that possess domains homologous to SARM1, we established distances and Ka/Ks values through 5671 pairwise species–species comparisons. We demonstrate that SARM1 diverged across species in a pattern similar to other SAM domain-containing proteins. This is surprising, because it was expected that SARM1 would behave more like its TIR domain relatives. Going along with this divorce from TIR, we also noted that SARM1’s TIR is under stronger purifying selection than the rest of the TIR domain-containing proteins (remaining highly conserved). In addition, SARM1’s synteny analysis reveals that the surrounding gene cluster is highly conserved, functioning as a potential nexus of gene functionality across species. Taken together, SARM1 demonstrates a unique evolutionary pattern, separate from the TIR domain protein family.     

Evolutionary history of the SARS-CoV-2 Gamma variant of concern (P.1): a perfect storm

Our goal was to describe in more detail the evolutionary history of Gamma and two derived lineages (P.1.1 and P.1.2), which are part of the arms race that SARS-CoV-2 wages with its host. A total of 4,977 sequences of the Gamma strain of SARS-CoV-2 from Brazil were analyzed. We detected 194 sites under positive selection in 12 genes/ORFs: Spike, N, M, E, ORF1a, ORF1b, ORF3, ORF6, ORF7a, ORF7b, ORF8, and ORF10. Some diagnostic sites for Gamma lacked a signature of positive selection in our study, but these were not fixed, apparently escaping the action of purifying selection. Our network analyses revealed branches leading to expanding haplotypes with sites under selection only detected when P.1.1 and P.1.2 were considered. The P.1.2 exclusive haplotype H_5 originated from a non-synonymous mutational step (H3509Y) in H_1 of ORF1a. The selected allele, 3509Y, represents an adaptive novelty involving ORF1a of P.1. Finally, we discuss how phenomena such as epistasis and antagonistic pleiotropy could limit the emergence of new alleles (and combinations thereof) in SARS-COV-2 lineages, maintaining infectivity in humans, while providing rapid response capabilities to face the arms race triggered by host immuneresponses.     

Population genetic structure of the Plasmodium vivax circumsporozoite protein (Pvcsp) in Sri Lanka

Molecular methods elucidate evolutionary and ecological processes in parasites, where interaction between hosts and parasites enlighten the evolution of parasite lifestyles and host defenses. Population genetics of Plasmodium vivax parasites accurately describe transmission dynamics of the parasites and evaluation of malaria control measures. As a first generation vaccine candidate against malaria, the Circumsporozoite Protein (CSP) has demonstrated significant potential in P. falciparum. Extensive polymorphism hinders the development of a potent malaria vaccine. Hence, the genetic diversity of Pvcsp was investigated for the first time in 60 Sri Lankan clinical isolates by obtaining the nucleotide sequence of the central repeat (CR) domain and examining the polymorphism of the peptide repeat motifs (PRMs), the genetic diversity indices and phylogenetic relationships. PCR amplicons determined size polymorphism of 610, 700 and 710 bp in Pvcsp of Sri Lanka where all amino acid sequences obtained were of the VK210 variant, consisting variable repeats of 4 different PRMs. The two most abundant PRMs of the CR domain, GDRADGQPA and GDRAAGQPA consisted ~ 2-4 repeats, while GNRAAGQPA was unique to the island. Though, different nucleotide sequences termed repeat allotypes (RATs) were observed for each PRM, these were synonymous contributing to a less polymorphic CR domain. The genetic diversity of Pvcsp in Sri Lanka was due to the number of repetitive peptide repeat motifs, point mutations, and intragenic recombination. The 19 amino acid haplotypes defined were exclusive to Sri Lanka, whereas the 194 Pvcsp sequences of global isolates generated 57 more distinct a.a. haplotypes of the VK210 variant. Strikingly, the CR domain of both VK210 and VK247 variants was under purifying selection interpreting the scarcity of CSP non-synonymous polymorphisms. Insights to the distribution of RATs in the CR region with geographic clustering of the P. vivax VK210 variant were revealed. The cladogram reiterated this unique geographic clustering of local (VK210) and global isolates (VK210 and VK247), which was further validated by the elevated fixation index values of the VK210 variant.     

European wild boars and domestic pigs display different polymorphic patterns in the Toll-like receptor (TLR) 1, TLR2, and TLR6 genes

During the last decade, the Toll-like receptors (TLRs) have been extensively studied, and their immense importance in innate immunity is now being unveiled. Here, we report pronounced differences—probably reflecting the domestication process and differences in selective pressure—between wild boars and domestic pigs regarding single nucleotide polymorphisms (SNPs) in TLR genes. The open reading frames of TLR1, TLR2, and TLR6 were sequenced in 25 wild boars, representing three populations, and in 15 unrelated domestic pigs of Hampshire, Landrace, and Large White origin. In total, 20, 27, and 26 SNPs were detected in TLR1, TLR2, and TLR6, respectively. In TLR1 and TLR2, the numbers of SNPs detected were significantly lower (P?≤?0.05, P?≤?0.01) in the wild boars than in the domestic pigs. In the wild boars, one major high frequency haplotype was found in all three genes, while the same pattern was exhibited only by TLR2 in the domestic pigs. The relative frequency of non-synonymous (dN) and synonymous (dS) SNPs was lower for the wild boars than for the domestic pigs in all three genes. In addition, differences in diversity between the genes were revealed: the mean heterozygosity at the polymorphic positions was markedly lower in TLR2 than in TLR1 and TLR6. Because of its localization—in proximity of the bound ligand—one of the non-synonymous SNPs detected in TLR6 may represent species-specific function on the protein level. Furthermore, the codon usage pattern in the genes studied deviated from the general codon usage pattern in Sus scrofa.     

Pervasive purifying selection characterizes the evolution of I2 homologs

We sampled 384 sequences related to the Solanum pimpinellifolium (=Lycopersicon pimpinellifolium) disease resistance (R) gene 12 from six species, potato, S. demissum, tomato, eggplant, pepper, and tobacco. These species represent increasing phylogenetic distance from potato to tobacco, within the family Solanaceae. Using sequence data from the nucleotide binding site (NBS) region of this gene, we tested models of gene family evolution and inferred patterns of selection acting on the NBS gene region and I2 gene family. We find that the I2 family has diversified within the family Solanaceae for at least 14 million years and evolves through a slow birth-and-death process requiring approximately 12 million years to homogenize gene copies within a species. Analyses of selection resolved a general pattern of strong purifying selection acting on individual codon positions within the NBS and on NBS lineages through time. Surprisingly, we find nine codon positions strongly affected by positive selection and six pairs of codon positions demonstrating correlated amino acid substitutions. Evolutionary analyses serve as bioinformatic tools with which to sort through the vast R gene diversity in plants and find candidates for new resistance specificities or to identify specific amino acid positions important for biochemical function. The slow birth-and-death evolution of I2 genes suggests that some NBS-leucine rich repeat-mediated resistances may not be overcome rapidly by virulence evolution and that the natural diversity of R genes is a potentially valuable source for durable resistance.     

Diversity and evolutionary relationship of nucleotide binding site-encoding disease-resistance gene analogues in sweet potato (<Emphasis Type="Italic">Ipomoea batatas</Emphasis> Lam.)

Most plant disease-resistance genes (R-genes) isolated so far encode proteins with a nucleotide binding site (NBS) domain and belong to a superfamily. NBS domains related to R-genes show a highly conserved backbone of an amino acid motif, which makes it possible to isolate resistance gene analogues (RGAs) by degenerate primers. Degenerate primers based on the conserved motif (P-loop and GLPL) of the NBS domain from R -genes were used to isolate RGAs from the genomic DNA of sweet potato cultivar Qingnong no.2. Five distinct clusters of RGAs (22 sequences) with the characteristic NBS representing a highly diverse sample were identified in sweet potato genomic DNA. Sequence identity among the 22 RGA nucleotide sequences ranged from 41.2% to 99.4%, while the deduced amino acid sequence identity from the 22 RGAs ranged from 20.6%to 100%. The analysis of sweet potato RGA sequences suggested mutation as the primary source of diversity. The phylogenetic analyses for RGA nucleotide sequences and deduced amino acids showed that RGAs from sweet potato were classified into two distinct groups--toll and interleukin receptor-1 (TIR)-NBS-LRR and non-TIR-NBS-LRR. The high degree of similarity between sweet potato RGAs and NBS sequences derived from R-genes cloned from tomato, tobacco, flax and potato suggest an ancestral relationship. Further studies showed that the ratio of non-synonymous to synonymous substitution within families was low. These data obtained from sweet potato suggest that the evolution of NBS-encoding sequences in sweet potato occur by the gradual accumulation of mutations leading to purifying selection and slow rates of divergence within distinct R-gene families.     

Molecular Variation in AVP and AVPR1a in New World Monkeys (Primates,Platyrrhini): Evolution and Implications for Social Monogamy

The neurohypophysial hormone arginine vasopressin (AVP) plays important roles in fluid regulation and vascular resistance. Differences in AVP receptor expression, particularly mediated through variation in the noncoding promoter region of the primary receptor for AVP (AVPR1a), may play a role in social phenotypes, particularly social monogamy, in rodents and humans. Among primates, social monogamy is rare, but is common among New World monkeys (NWM). AVP is a nonapeptide and generally conserved among eutherian mammals, although a recent paper demonstrated that some NWM species possess a novel form of the related neuropeptide hormone, oxytocin. We therefore characterized variation in the AVP and AVPR1a genes in 22 species representing every genus in the three major platyrrhine families (Cebidae, Atelidae and Pitheciidae). For AVP, a total of 16 synonymous substitutions were detected in 15 NWM species. No non-synonymous substitutions were noted, hence, AVP is conserved in NWM. By contrast, relative to the human AVPR1a, 66 predicted amino acids (AA) substitutions were identified in NWM. The AVPR1a N-terminus (ligand binding domain), third intracellular (G-protein binding domain), and C-terminus were variable among species. Complex evolution of AVPR1a is also apparent in NWM. A molecular phylogenetic tree inferred from AVPR1a coding sequences revealed some consensus taxonomic separation by families, but also a mixed group composed of genera from all three families. The overall dN/dS ratio of AVPR1a was 0.11, but signals of positive selection in distinct AVPR1a regions were observed, including the N-terminus, in which we identified six potential positive selection sites. AA substitutions at positions 241, 319, 399 and 409 occurred uniquely in marmosets and tamarins. Our results enhance the appreciation of genetic diversity in the mammalian AVP/AVPR1a system, and set the stage for molecular modeling of the neurohypophyseal hormones and social behavior in primates.     

Adaptive Protein Evolution in Animals and the Effective Population Size Hypothesis

The rate at which genomes adapt to environmental changes and the prevalence of adaptive processes in molecular evolution are two controversial issues in current evolutionary genetics. Previous attempts to quantify the genome-wide rate of adaptation through amino-acid substitution have revealed a surprising diversity of patterns, with some species (e.g. Drosophila) experiencing a very high adaptive rate, while other (e.g. humans) are dominated by nearly-neutral processes. It has been suggested that this discrepancy reflects between-species differences in effective population size. Published studies, however, were mainly focused on model organisms, and relied on disparate data sets and methodologies, so that an overview of the prevalence of adaptive protein evolution in nature is currently lacking. Here we extend existing estimators of the amino-acid adaptive rate by explicitly modelling the effect of favourable mutations on non-synonymous polymorphism patterns, and we apply these methods to a newly-built, homogeneous data set of 44 non-model animal species pairs. Data analysis uncovers a major contribution of adaptive evolution to the amino-acid substitution process across all major metazoan phyla—with the notable exception of humans and primates. The proportion of adaptive amino-acid substitution is found to be positively correlated to species effective population size. This relationship, however, appears to be primarily driven by a decreased rate of nearly-neutral amino-acid substitution because of more efficient purifying selection in large populations. Our results reveal that adaptive processes dominate the evolution of proteins in most animal species, but do not corroborate the hypothesis that adaptive substitutions accumulate at a faster rate in large populations. Implications regarding the factors influencing the rate of adaptive evolution and positive selection detection in humans vs. other organisms are discussed.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

Systematic Analysis and Comparison of Nucleotide-Binding Site Disease Resistance Genes in a Diploid Cotton Gossypium raimondii

Plant disease resistance genes are a key component of defending plants from a range of pathogens. The majority of these resistance genes belong to the super-family that harbors a Nucleotide-binding site (NBS). A number of studies have focused on NBS-encoding genes in disease resistant breeding programs for diverse plants. However, little information has been reported with an emphasis on systematic analysis and comparison of NBS-encoding genes in cotton. To fill this gap of knowledge, in this study, we identified and investigated the NBS-encoding resistance genes in cotton using the whole genome sequence information of Gossypium raimondii. Totally, 355 NBS-encoding resistance genes were identified. Analyses of the conserved motifs and structural diversity showed that the most two distinct features for these genes are the high proportion of non-regular NBS genes and the high diversity of N-termini domains. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication of disease resistance genes could play an important role in cotton by leading to an increase in the functional diversity of the cotton NBS-encoding genes. Analyses of phylogenetic comparisons indicated that, in cotton, the NBS-encoding genes with TIR domain not only have their own evolution pattern different from those of genes without TIR domain, but also have their own species-specific pattern that differs from those of TIR genes in other plants. Analyses of the correlation between disease resistance QTL and NBS-encoding resistance genes showed that there could be more than half of the disease resistance QTL associated to the NBS-encoding genes in cotton, which agrees with previous studies establishing that more than half of plant resistance genes are NBS-encoding genes.     

Relationship Among Intron Length,Gene Expression,and Nucleotide Diversity in the Pacific Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

Crassostrea gigas is a model mollusk, but its genetic features have not been studied comprehensively. In this study, we used whole-genome resequencing data to identify and characterize nucleotide diversity and population recombination rate in a diverse collection of 21 C. gigas samples. Our analyses revealed that C. gigas harbors both extremely high genetic diversity and recombination rates across the whole genome as compared with those of the other taxa. The noncoding regions, introns, intergenic spacers, and untranslated regions (UTRs) showed a lower level diversity than the synonymous sites. The larger introns tended to have lower diversity. Moreover, we found a negative association of the non-synonymous diversity with gene expression, which suggested that purifying selection played an important role in shaping genetic diversity. The nucleotide diversity at the 100- and 50-kb levels was positively correlated with population recombination rates, which was expected if the diversity was shaped by purifying selection or hitchhiking of advantageous mutants. Our work gives a general picture of the oyster’s polymorphism pattern and its association with recombination rates.     

Isolation and diversity analysis of resistance gene analogues (RGAs) from cultivated and wild strawberries

Degenerate oligonucleotide primers, designed based on conserved regions of Nucleotide Binding Site (NBS) domains from previously cloned plant resistance genes, were used to isolate Resistance Gene Analogues (RGAs) from wild and cultivated strawberries. Seven distinct families of RGAs of the NBS-LRR type were identified from two related wild species, Fragaria vesca and F. chiloensis, and six different Fragaria × ananassa cultivars. With one exception (GAV-3), the deduced amino acid sequences of strawberry RGAs showed strong similarity to TIR (Toll Interleukin I Receptor)-type R genes from Arabidopsis, tobacco and flax, suggesting the existence of common ancestors. GAV-3 seemed to be more closely related to the non-TIR type. Further studies showed that the recombination level and the ratio of non-synonymous to synonymous substitutions within families were low. These data suggest that NBS-encoding sequences of RGAs in strawberry are subject to a gradual accumulation of mutations leading to purifying selection, rather than to a diversifying process. The present paper is the first report on RGAs in strawberry.Communicated by M.-A. Grandbastien     

Complete mitochondrial genome of the <Emphasis Type="Italic">Pleuronichthys lighti</Emphasis> (Pleuronectiformes,Pleuronectidae) with phylogenetic consideration

Mitochondrial genome has been used to shed light on many fields of both basic and applied research, including the study of molecular evolution. The complete mitochondrial genome sequence of 17368 bp nucleotides from the Pleuronichthys lighti was determined. It was a circular double-stranded DNA molecule with identical set of 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes as well as a non-coding control region. Stand asymmetry in the nucleotide composition was reflected in the codon usage of genes oriented in opposite directions. In the control region, we identified the extended termination associated sequence domain, the central conserved sequence block domain and the conserved sequence block domain, and two complete repeat region. They were “TTACAATA” and “TGTTGTAA”, respectively. All known 12 mitochondrial genomes of Pleuronectinae fishes were downloaded and analyzed; there were 5570 variable sites in the consensus sequences of 15241 base pairs, calculation of total sites were 35.5%. The highest sequence divergence was 50% (ATP8) and the Kimura-2-parameter genetic distance was 0.235 (ND6), whereas the COIII had the lowest sequence divergence (28.8%) and genetic distance (0.128); the protein coding genes were mainly acted by purifying selection which was detected by selection tests. Analysis of confidence and the information content for per nucleotide revealed ND5, ATP6, COI and ND4 genes were suitable molecular markers for phylogenetic study of Pleuronectinae fishes. Phylogenetic analysis using Bayesian computational algorithms based on COI genes provided support for the taxonomic status of P. lighti, which was consistent with the traditional taxonomy.     

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD⁴⁺ T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d _S) and non-synonymous (d _N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d _N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.