Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O₂ transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O₂ affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.Subject terms: Evolutionary genetics, Evolutionary biology     

Polymorphism and divergence at the 5' flanking region of the sex- determining locus, Sry, in mice

We have investigated patterns of evolution in the nonrecombining portion of the Y chromosome in mice by comparing levels of polymorphism within Mus domesticus with levels of divergence between M. domesticus and M. spretus. A 1,277-bp fragment of noncoding sequence flanking the sex determining locus (Sry) was PCR amplified, and 1,063 bases were sequenced and compared among 20 M. domesticus and 1 M. spretus. Two polymorphic base substitutions and two polymorphic insertion/deletion sites were identified within M. domesticus; nucleotide diversity was estimated to be 0.1%. Divergence between M. domesticus and M. spretus for this region (1.9%) was slightly lower than the average divergence of single-copy nuclear DNA for these species. Comparison of levels of polymorphism and divergence at Sry with levels of polymorphism and divergence in the mitochondrial DNA control region provided no evidence of a departure from the expectations of neutral molecular evolution. These findings are consistent with the presumed lack of function for much of the Y chromosome.      

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

The genomic rate of adaptive amino acid substitution in Drosophila

The proportion of amino acid substitutions driven by adaptive evolution can potentially be estimated from polymorphism and divergence data by an extension of the McDonald-Kreitman test. We have developed a maximum-likelihood method to do this and have applied our method to several data sets from three Drosophila species: D. melanogaster, D. simulans, and D. yakuba. The estimated number of adaptive substitutions per codon is not uniformly distributed among genes, but follows a leptokurtic distribution. However, the proportion of amino acid substitutions fixed by adaptive evolution seems to be remarkably constant across the genome (i.e., the proportion of amino acid substitutions that are adaptive appears to be the same in fast-evolving and slow-evolving genes; fast-evolving genes have higher numbers of both adaptive and neutral substitutions). Our estimates do not seem to be significantly biased by selection on synonymous codon use or by the assumption of independence among sites. Nevertheless, an accurate estimate is hampered by the existence of slightly deleterious mutations and variations in effective population size. The analysis of several Drosophila data sets suggests that approximately 25% +/- 20% of amino acid substitutions were driven by positive selection in the divergence between D. simulans and D. yakuba.     

Positive and negative selection on noncoding DNA close to protein-coding genes in wild house mice

During the past two decades, evidence has accumulated of adaptive evolution within protein-coding genes in a variety of species. However, with the exception of Drosophila and humans, little is known about the extent of adaptive evolution in noncoding DNA. Here, we study regions upstream and downstream of protein-coding genes in the house mouse Mus musculus castaneus, a species that has a much larger effective population size (N(e)) than humans. We analyze polymorphism data for 78 genes from 15 wild-caught M. m. castaneus individuals and divergence to a closely related species, Mus famulus. We find high levels of nucleotide diversity and moderate levels of selective constraint in upstream and downstream regions compared with nonsynonymous sites of protein-coding genes. From the polymorphism data, we estimate the distribution of fitness effects (DFE) of new mutations and infer that most new mutations in upstream and downstream regions behave as effectively neutral and that only a small fraction is strongly negatively selected. We also estimate the fraction of substitutions that have been driven to fixation by positive selection (α) and the ratio of adaptive to neutral divergence (ω(α)). We find that α for upstream and downstream regions (～ 10%) is much lower than α for nonsynonymous sites (～ 50%). However, ω(α) estimates are very similar for nonsynonymous sites (～ 10%) and upstream and downstream regions (～ 5%). We conclude that negative selection operating in upstream and downstream regions of M. m. castaneus is weak and that the low values of α for upstream and downstream regions relative to nonsynonymous sites are most likely due to the presence of a higher proportion of neutrally evolving sites and not due to lower absolute rates of adaptive substitution.     

Purification and Characterization of the Coniferyl Aldehyde Dehydrogenase from Pseudomonas sp. Strain HR199 and Molecular Characterization of the Gene

The coniferyl aldehyde dehydrogenase (CALDH) of Pseudomonas sp. strain HR199 (DSM7063), which catalyzes the NAD⁺-dependent oxidation of coniferyl aldehyde to ferulic acid and which is induced during growth with eugenol as the carbon source, was purified and characterized. The native protein exhibited an apparent molecular mass of 86,000 ± 5,000 Da, and the subunit mass was 49.5 ± 2.5 kDa, indicating an α₂ structure of the native enzyme. The optimal oxidation of coniferyl aldehyde to ferulic acid was obtained at a pH of 8.8 and a temperature of 26°C. The K_m values for coniferyl aldehyde and NAD⁺ were about 7 to 12 μM and 334 μM, respectively. The enzyme also accepted other aromatic aldehydes as substrates, whereas aliphatic aldehydes were not accepted. The NH₂-terminal amino acid sequence of CALDH was determined in order to clone the encoding gene (calB). The corresponding nucleotide sequence was localized on a 9.4-kbp EcoRI fragment (E94), which was subcloned from a Pseudomonas sp. strain HR199 genomic library in the cosmid pVK100. The partial sequencing of this fragment revealed an open reading frame of 1,446 bp encoding a protein with a relative molecular weight of 51,822. The deduced amino acid sequence, which is reported for the first time for a structural gene of a CALDH, exhibited up to 38.5% amino acid identity (60% similarity) to NAD⁺-dependent aldehyde dehydrogenases from different sources.     

Balance between promiscuity and specificity in phage λ host range

As hosts acquire resistance to viruses, viruses must overcome that resistance to re-establish infectivity, or go extinct. Despite the significant hurdles associated with adapting to a resistant host, viruses are evolutionarily successful and maintain stable coevolutionary relationships with their hosts. To investigate the factors underlying how pathogens adapt to their hosts, we performed a deep mutational scan of the region of the λ tail fiber tip protein that mediates contact with the receptor on λ’s host, Escherichia coli. Phages harboring amino acid substitutions were subjected to selection for infectivity on wild type E. coli, revealing a highly restrictive fitness landscape, in which most substitutions completely abrogate function. A subset of positions that are tolerant of mutation in this assay, but diverse over evolutionary time, are associated with host range expansion. Imposing selection for phage infectivity on three λ-resistant hosts, each harboring a different missense mutation in the λ receptor, reveals hundreds of adaptive variants in λ. We distinguish λ variants that confer promiscuity, a general ability to overcome host resistance, from those that drive host-specific infectivity. Both processes may be important in driving adaptation to a novel host.Subject terms: Bacteriophages, Molecular evolution, Viral genetics     

Cold-Adapted β-Galactosidase from the Antarctic Psychrophile Pseudoalteromonas haloplanktis

The β-galactosidase from the Antarctic gram-negative bacterium Pseudoalteromonas haloplanktis TAE 79 was purified to homogeneity. The nucleotide sequence and the NH₂-terminal amino acid sequence of the purified enzyme indicate that the β-galactosidase subunit is composed of 1,038 amino acids with a calculated M_r of 118,068. This β-galactosidase shares structural properties with Escherichia coli β-galactosidase (comparable subunit mass, 51% amino sequence identity, conservation of amino acid residues involved in catalysis, similar optimal pH value, and requirement for divalent metal ions) but is characterized by a higher catalytic efficiency on synthetic and natural substrates and by a shift of apparent optimum activity toward low temperatures and lower thermal stability. The enzyme also differs by a higher pI (7.8) and by specific thermodynamic activation parameters. P. haloplanktis β-galactosidase was expressed in E. coli, and the recombinant enzyme displays properties identical to those of the wild-type enzyme. Heat-induced unfolding monitored by intrinsic fluorescence spectroscopy showed lower melting point values for both P. haloplanktis wild-type and recombinant β-galactosidase compared to the mesophilic enzyme. Assays of lactose hydrolysis in milk demonstrate that P. haloplanktis β-galactosidase can outperform the current commercial β-galactosidase from Kluyveromyces marxianus var. lactis, suggesting that the cold-adapted β-galactosidase could be used to hydrolyze lactose in dairy products processed in refrigerated plants.     

Genetic analysis of genome-scale recombination rate evolution in house mice

The rate of meiotic recombination varies markedly between species and among individuals. Classical genetic experiments demonstrated a heritable component to population variation in recombination rate, and specific sequence variants that contribute to recombination rate differences between individuals have recently been identified. Despite these advances, the genetic basis of species divergence in recombination rate remains unexplored. Using a cytological assay that allows direct in situ imaging of recombination events in spermatocytes, we report a large (∼30%) difference in global recombination rate between males of two closely related house mouse subspecies (Mus musculus musculus and M. m. castaneus). To characterize the genetic basis of this recombination rate divergence, we generated an F2 panel of inter-subspecific hybrid males (n = 276) from an intercross between wild-derived inbred strains CAST/EiJ (M. m. castaneus) and PWD/PhJ (M. m. musculus). We uncover considerable heritable variation for recombination rate among males from this mapping population. Much of the F2 variance for recombination rate and a substantial portion of the difference in recombination rate between the parental strains is explained by eight moderate- to large-effect quantitative trait loci, including two transgressive loci on the X chromosome. In contrast to the rapid evolution observed in males, female CAST/EiJ and PWD/PhJ animals show minimal divergence in recombination rate (∼5%). The existence of loci on the X chromosome suggests a genetic mechanism to explain this male-biased evolution. Our results provide an initial map of the genetic changes underlying subspecies differences in genome-scale recombination rate and underscore the power of the house mouse system for understanding the evolution of this trait.     

Contributions of Protein-Coding and Regulatory Change to Adaptive Molecular Evolution in Murid Rodents

The contribution of regulatory versus protein change to adaptive evolution has long been controversial. In principle, the rate and strength of adaptation within functional genetic elements can be quantified on the basis of an excess of nucleotide substitutions between species compared to the neutral expectation or from effects of recent substitutions on nucleotide diversity at linked sites. Here, we infer the nature of selective forces acting in proteins, their UTRs and conserved noncoding elements (CNEs) using genome-wide patterns of diversity in wild house mice and divergence to related species. By applying an extension of the McDonald-Kreitman test, we infer that adaptive substitutions are widespread in protein-coding genes, UTRs and CNEs, and we estimate that there are at least four times as many adaptive substitutions in CNEs and UTRs as in proteins. We observe pronounced reductions in mean diversity around nonsynonymous sites (whether or not they have experienced a recent substitution). This can be explained by selection on multiple, linked CNEs and exons. We also observe substantial dips in mean diversity (after controlling for divergence) around protein-coding exons and CNEs, which can also be explained by the combined effects of many linked exons and CNEs. A model of background selection (BGS) can adequately explain the reduction in mean diversity observed around CNEs. However, BGS fails to explain the wide reductions in mean diversity surrounding exons (encompassing ∼100 Kb, on average), implying that there is a substantial role for adaptation within exons or closely linked sites. The wide dips in diversity around exons, which are hard to explain by BGS, suggest that the fitness effects of adaptive amino acid substitutions could be substantially larger than substitutions in CNEs. We conclude that although there appear to be many more adaptive noncoding changes, substitutions in proteins may dominate phenotypic evolution.     

Wheat Domestication Accelerated Evolution and Triggered Positive Selection in the β-Xylosidase Enzyme of Mycosphaerella graminicola

Plant cell wall degrading enzymes (PCWDEs) of plant pathogens are receiving increasing interest for their potential to trigger plant defense reactions. In an antagonistic co-evolutionary arms race between host and pathogen, PCWDEs could be under strong selection. Here, we tested the hypothesis that PCWDEs in the fungal wheat pathogen Mycosphaerella graminicola have been positively selected by analyzing ratios of non-synonymous and synonymous nucleotide changes in the genes encoding these enzymes. Analyses of five PCWDEs demonstrated that one (β-xylosidase) has been under strong positive selection and experienced an accelerated rate of evolution. In contrast, PCWDEs in the closest relatives of M. graminicola collected from wild grasses did not show evidence for selection or deviation from a molecular clock. Since the genealogical divergence of M. graminicola from these latter species coincided with the onset of agriculture, we hypothesize that the recent domestication of the host plant and/or agricultural practices triggered positive selection in β-xylosidase and that this enzyme played a key role in the emergence of a host-specialized pathogen.     

Nucleotide Sequence Evolution at the κ-Casein Locus: Evidence for Positive Selection within the Family Bovidae

κ-Casein is a mammalian milk protein involved in a number of important physiological processes. In the gut, the ingested protein is split into an insoluble peptide (para κ-casein) and a soluble hydrophilic glycopeptide (caseinomacropeptide). Caseinomacropeptide is responsible for increased efficiency of digestion, prevention of neonate hypersensitivity to ingested proteins, and inhibition of gastric pathogens. Variation within this peptide has significant effects associated with important traits such as milk production. The nucleotide sequences for regions of κ-casein exon and intron four were determined for representatives of the artiodactyl family Bovidae. The pattern of nucleotide substitution in κ-casein sequences for distantly related bovid taxa demonstrates that positive selection has accelerated their divergence at the amino acid sequence level. This selection has differentially influenced the molecular evolution of the two κ-casein split peptides and is focused within a 34-codon region of caseinomacropeptide.     

Nucleotide Variation, Linkage Disequilibrium and Founder-Facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (K_ST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.     

Evolutionary Relationships among Five Subspecies of MUS MUSCULUS Based on Restriction Enzyme Cleavage Patterns of Mitochondrial DNA

The intra- and intersubspecific genetic distances between five subspecies of Mus musculus were estimated from restriction enzyme cleavage patterns or maps of mitochondrial DNA (mtDNA). The European subspecies, M. m. domesticus and Asian subspecies, M. m. bactrianus, M. m. castaneus, M. m. molossinus and M. m. urbanus were examined. For each subspecies, except M. m. urbanus, at least two local races from widely separated localities were examined. Intrasubspecific heterogeneity was found in the mtDNA cleavage patterns of M. m. bactrianus and M. m. castaneus. M. m. molossinus and M. m. domesticus, however, revealed no intrasubspecific heterogeneity. Four of the subspecies had distinct cleavage patterns. The fifth, M. m. urbanus, had cleavage patterns identical to those of M. m. castaneus with several enzymes. Estimates of genetic distances between the various races and subspecies were obtained by comparing cleavage maps of the mtDNAs with various restriction enzymes. Nucleotide sequence divergences of mtDNA between local races were estimated to be less than 0.4% in M. m. bactrianus and less than 0.3% in M. m. castaneus. The times of divergence of both subspecies were calculated to be 0.1–0.2 x 10⁶ years. These values suggest that the intrasubspecific divergence began some 0.1–0.2 x 10⁶ years ago. On the other hand, nucleotide sequence divergences between European subspecies M. m. domesticus and Asian subspecies M. m. bactrianus and M. m. castaneus were 7.1% and 5.8%, respectively. The times of divergence were calculated to be 2.1–2.6 x 10⁶ years. Further, the nucleotide sequence divergence and time of divergence between M. m. molossinus and the other two Asian subspecies were comparable to those between M. m. molossinus and M. m. domesticus (about 3% and 1 x 10⁶ years, respectively). These results suggest that M. m. molossinus is situated in a unique evolutionary position among Asian subspecies.     

The Evolution of Vp1 Gene in Enterovirus C Species Sub-Group That Contains Types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99

Genus Enterovirus (Family Picornaviridae,) consists of twelve species divided into genetically diverse types by their capsid protein VP1 coding sequences. Each enterovirus type can further be divided into intra-typic sub-clusters (genotypes). The aim of this study was to elucidate what leads to the emergence of novel enterovirus clades (types and genotypes). An evolutionary analysis was conducted for a sub-group of Enterovirus C species that contains types Coxsackievirus A21 (CVA-21), CVA-24, Enterovirus C95 (EV-C95), EV-C96 and EV-C99. VP1 gene datasets were collected and analysed to infer the phylogeny, rate of evolution, nucleotide and amino acid substitution patterns and signs of selection. In VP1 coding gene, high intra-typic sequence diversities and robust grouping into distinct genotypes within each type were detected. Within each type the majority of nucleotide substitutions were synonymous and the non-synonymous substitutions tended to cluster in distinct highly polymorphic sites. Signs of positive selection were detected in some of these highly polymorphic sites, while strong negative selection was indicated in most of the codons. Despite robust clustering to intra-typic genotypes, only few genotype-specific ‘signature’ amino acids were detected. In contrast, when different enterovirus types were compared, there was a clear tendency towards fixation of type-specific ‘signature’ amino acids. The results suggest that permanent fixation of type-specific amino acids is a hallmark associated with evolution of different enterovirus types, whereas neutral evolution and/or (frequency-dependent) positive selection in few highly polymorphic amino acid sites are the dominant forms of evolution when strains within an enterovirus type are compared.     

Uptake of the β-Lactam Precursor α-Aminoadipic Acid in Penicillium chrysogenum Is Mediated by the Acidic and the General Amino Acid Permease

External addition of the β-lactam precursor α-aminoadipic acid to the filamentous fungus Penicillium chrysogenum leads to an increased intracellular α-aminoadipic acid concentration and an increase in penicillin production. The exact route for α-aminoadipic acid uptake is not known, although the general amino acid and acidic amino acid permeases have been implicated in this process. Their corresponding genes, PcGAP1 and PcDIP5, of P. chrysogenum were cloned and functionally expressed in a mutant of Saccharomyces cerevisiae (M4276) in which the acidic amino acid and general amino acid permease genes (DIP5 and GAP1, respectively) are disrupted. Transport assays show that both PcGap1 and PcDip5 mediated the uptake of α-aminoadipic acid, although PcGap1 showed a higher affinity for α-aminoadipic acid than PcDip5 (K_m values, 230 and 800 μM, respectively). Leucine strongly inhibits α-aminoadipic acid transport via PcGap1 but not via PcDip5. This difference was exploited to estimate the relative contribution of each transport system to the α-aminoadipic acid flux in β-lactam-producing P. chrysogenum. The transport measurements demonstrate that both PcGap1 and PcDip5 contribute to the α-aminoadipic acid flux.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

The Role of Gene Duplication and Unconstrained Selective Pressures in the Melanopsin Gene Family Evolution and Vertebrate Circadian Rhythm Regulation

Melanopsin is a photosensitive cell protein involved in regulating circadian rhythms and other non-visual responses to light. The melanopsin gene family is represented by two paralogs, OPN4x and OPN4m, which originated through gene duplication early in the emergence of vertebrates. Here we studied the melanopsin gene family using an integrated gene/protein evolutionary approach, which revealed that the rhabdomeric urbilaterian ancestor had the same amino acid patterns (DRY motif and the Y and E conterions) as extant vertebrate species, suggesting that the mechanism for light detection and regulation is similar to rhabdomeric rhodopsins. Both OPN4m and OPN4x paralogs are found in vertebrate genomic paralogons, suggesting that they diverged following this duplication event about 600 million years ago, when the complex eye emerged in the vertebrate ancestor. Melanopsins generally evolved under negative selection (ω = 0.171) with some minor episodes of positive selection (proportion of sites = 25%) and functional divergence (θ_I = 0.349 and θ_II = 0.126). The OPN4m and OPN4x melanopsin paralogs show evidence of spectral divergence at sites likely involved in melanopsin light absorbance (200F, 273S and 276A). Also, following the teleost lineage-specific whole genome duplication (3R) that prompted the teleost fish radiation, type I divergence (θ_I = 0.181) and positive selection (affecting 11% of sites) contributed to amino acid variability that we related with the photo-activation stability of melanopsin. The melanopsin intracellular regions had unexpectedly high variability in their coupling specificity of G-proteins and we propose that Gq/11 and Gi/o are the two G-proteins most-likely to mediate the melanopsin phototransduction pathway. The selection signatures were mainly observed on retinal-related sites and the third and second intracellular loops, demonstrating the physiological plasticity of the melanopsin protein group. Our results provide new insights on the phototransduction process and additional tools for disentangling and understanding the links between melanopsin gene evolution and the specializations observed in vertebrates, especially in teleost fish.     

Reduced nucleotide variability at an androgen-binding protein locus (Abpa) in house mice: evidence for positive natural selection.

Previous work has shown that the gene for the alpha subunit of androgen-binding protein, Abpa, may be involved in premating isolation between different subspecies of the house mouse, Mus musculus. We investigated patterns of DNA sequence variation at Abpa within and between species of mice to test several predictions of a model of neutral molecular evolution. Intraspecific variation among 10 Mus musculus domesticus alleles was compared with divergence between M. m. domesticus and M. caroli for Abpa and two X-linked genes, Glra2 and Amg. No variation was observed at Abpa within M. m. domesticus. The ratio of polymorphism to divergence was significantly lower at Abpa than at Glra2 and Amg, despite the fact that all three genes experience similar rates of recombination. Interspecific comparisons among M. m. domesticus, Mus musculus musculus, Mus musculus castaneus, Mus spretus, Mus spicilegus, and Mus caroli revealed that the ratio of nonsynonymous substitutions to synonymous substitutions on a per-site basis (Ka/Ks) was generally greater than one. The combined observations of no variation at Abpa within M. m: domesticus and uniformly high Ka/Ks values between species suggest that positive directional selection has acted recently at this locus.