Patterns of nucleotide polymorphism and divergence in the odorant-binding protein genes OS-E and OS-F: analysis in the melanogaster species subgroup of Drosophila

The Olfactory Specific-E and -F genes (OS-E and OS-F) belong to the odorant-binding protein gene family, which includes the general odorant-binding proteins and the pheromone-binding proteins. In Drosophila melanogaster, these genes are arranged in tandem in a genomic region near the centromere of chromosome arm 3R. We examined the pattern of DNA sequence variation in an approximately 7-kb genomic region encompassing the two OS genes in four species of the melanogaster subgroup of Drosophila and in a population sample of D. melanogaster. We found that both the OS-E and the OS-F gene are present in all surveyed species. Nucleotide divergence estimates would support that the two genes are functional, although they diverge in their functional constraint. The pattern of nucleotide variation in D. melanogaster also differed between genes. Variation in the OS-E gene region exhibited an unusual and distinctive pattern: (i) a relatively high number of fixed amino acid replacements in the encoded protein and (ii) a peak of nucleotide polymorphism around the OS-E gene. These results are unlikely under the neutral model and suggest the action of natural selection in the evolution of the two odorant-binding protein genes.     

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.     

Historical Selection, Amino Acid Polymorphism and Lineage-Specific Divergence at the G6pd Locus in Drosophila Melanogaster and D. Simulans

The nucleotide diversity across 1705 bp of the G6pd gene is studied in 50 Drosophila melanogaster and 12 D. simulans lines. Our earlier report contrasted intraspecific polymorphism and interspecific differences at silent and replacement sites in these species. This report expands the number of European and African lines and examines the pattern of polymorphism with respect to the common A/B allozymes. In D. melanogaster the silent nucleotide diversity varies 2.8-fold across localities. The B allele sequences are two- to fourfold more variable than the derived A allele, and differences between allozymes are twice as among B alleles. There is strong linkage disequilibrium across the G6pd region. In both species the level of silent polymorphism increases from the 5' to 3' ends, while there is no comparable pattern in level of silent site divergence or fixation. The neutral model is not rejected in either species. Using D. yakuba as an outgroup, the D. melanogaster lineage shows a twofold greater rate of silent fixation, but less than half the rate of amino acid replacement. Lineage-specific differences in mutation fixation are inconsistent with neutral expectations and suggest the interaction of species-specific population size differences with both weakly advantageous and deleterious selection.     

Recurrent deletion and gene presence/absence polymorphism: telomere dynamics dominate evolution at the tip of 3L in Drosophila melanogaster and D. simulans

Although Drosophila melanogaster has been the subject of intensive analysis of polymorphism and divergence, little is known about the distribution of variation at the most distal regions of chromosomes arms. Here we report a survey of genetic variation on the tip of 3L in D. melanogaster and D. simulans. Levels of single nucleotide polymorphism in the most distal euchromatic sequence are approximately one order of magnitude less than that typically observed in genomic regions of normal crossing over, consistent with what might be expected under models of linked selection in regions of low crossing over. However, despite this reduced level of nucleotide variation, we found abundant deletion polymorphism. These deletions create at least three gene presence/absence polymorphisms within D. melanogaster: the putative G-protein coupled receptor mthl-8 (which is the most distal known or predicted gene on 3L) and the unannotated mRNAs AY060886 and BT006009. Strikingly, D. simulans is also segregating deletions that cause mthl8 presence/absence polymorphism. Breakpoint sequencing and tests of correlations with segregating SNPs in D. melanogaster suggest that each deletion is unique. Cloned breakpoint sequences revealed the presence of Het-A elements just distal to unique, canonical euchromatic sequences. This pattern suggests a model in which repeated telomeric deficiencies cause deletions of euchromatic sequence followed by subsequent "healing" by retrotranposition of Het-A elements. These data reveal the dominance of telomeric dynamics on the evolution of closely linked sequences in Drosophila.     

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.     

The Rosy Region of Drosophila Melanogaster and Drosophila Simulans. I. Contrasting Levels of Naturally Occurring DNA Restriction Map Variation and Divergence

A 40-kb region around the rosy and snake loci was analyzed for restriction map variation among 60 lines of Drosophila melanogaster and 30 lines of Drosophila simulans collected together at a single locality in Raleigh, North Carolina. DNA sequence variation in D. simulans was estimated to be 6.3 times greater than in D. melanogaster (heterozygosities per nucleotide of 1.9% vs. 0.3%). This result stands in marked contrast to results of studies of phenotypic variation including proteins (allozymes), morphology and chromosome arrangements which are generally less variable and less geographically differentiated in D. simulans. Intraspecific polymorphism is not distributed uniformly over the 40-kb region. The level of heterozygosity per nucleotide varies more than 12-fold across the region in D. simulans, being highest over the hsc2 gene. Similar, though less extreme, variation in heterozygosity is also observed in D. melanogaster. Average interspecific divergence (corrected for intraspecific polymorphism) averaged 3.8%. The pattern of interspecific divergence over the 40-kb region shows some disparities with the spatial distribution of intraspecific variation, but is generally consistent with selective neutrality predictions: the most polymorphic regions within species are generally the most divergent between species. Sequence-length polymorphism is observed for D. melanogaster to be at levels comparable to other gene regions in this species. In contrast, no sequence length variation was observed among D. simulans chromosomes (limit of resolution approximately 100 bp). These data indicate that transposable elements play at best a minor role in the generation of naturally occurring genetic variation in D. simulans compared to D. melanogaster. We hypothesize that differences in species effective population size are the major determinant of the contrasting levels and patterns of DNA sequence and insertion/deletion variation that we report here and the patterns of allozyme and morphological variation and differentiation reported by other workers for these two species.     

Polymorphism and Divergence at a Drosophila Pseudogene Locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcpψ) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcpψ fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcpψ revealed lower levels of polymorphism ( & = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcpψ within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcpψ lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Lack of nucleotide polymorphism in the Y-linked sperm flagellar dynein gene Dhc-Yh3 of Drosophila melanogaster and D. simulans

We studied levels of intra- and interspecific nucleotide variation associated with a Y-linked gene in five members of the Drosophila melanogaster subgroup. Using published sequence for 348 bp of the Dhc-Yh3 gene, and degenerate PCR primers designed from comparisons of the sea urchin and Chlamydomonas flagellar dynein genes, we recovered a 1738-bp region in D. melanogaster. Analyses of sequence variation in a worldwide collection of 11 lines of D. melanogaster and 10 lines of D. simulans found only a single silent polymorphism in the latter species. The synonymous site divergence per site for Dhc-Yh3 is comparable to values for X and autosomal genes. Assuming a Wright-Fisher population model, the lack of variation is statistically less than expected using appropriately reduced estimates of theta from the X and autosomes. Because the Y chromosome encodes only six known genes, genetic hitchhiking associated with background selection is unlikely to explain this low variation. Conversely, adaptive hitchhiking, as associated with sex-ratio chromosomes, or a large variance in male fertility may reduce the polymorphism on the Y chromosome. Codon bias is very low, as seen for other genes in regions of low recombination.     

Molecular evolutionary analysis of a histone gene repeating unit from Drosophila simulans

A repeating unit of the histone gene cluster from Drosophila simulans containing the H1, H2A, H2B and H4 genes (the H3 gene region has already been analyzed) was cloned and analyzed. A nucleotide sequence of about 4.6 kbp was determined to study the nucleotide divergence and molecular evolution of the histone gene cluster. Comparison of the structure and nucleotide sequence with those of Drosophila melanogaster showed that the four histone genes were located at identical positions and in the same directions. The proportion of different nucleotide sites was 6.3% in total. The amino acid sequence of H1 was divergent, with a 5.1% difference. However, no amino acid change has been observed for the other three histone proteins. Analysis of the GC contents and the base substitution patterns in the two lineages, D. melanogaster and D. simulans, with a common ancestor showed the following. 1) A strong negative correlation was found between the GC content and the nucleotide divergence in the whole repeating unit. 2) The mode of molecular evolution previously found for the H3 gene was also observed for the whole repeating unit of histone genes; the nucleotide substitutions were stationary in the 3' and spacer regions, and there was a directional change of the codon usage to the AT-rich codons. 3) No distinct difference in the mode or pattern of molecular evolution was detected for the histone gene repeating unit in the D. melanogaster and D. simulans lineages. These results suggest that selectional pressure for the coding regions of histones, which eliminate A and T, is less effective in the D. melanogaster and D. simulans lineages than in the other GC-rich species.     

Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

Excess Polymorphism at the Adh Locus in DROSOPHILA MELANOGASTER

The evolutionary history of a region of DNA encompassing the Adh locus is studied by comparing patterns of variation in Drosophila melanogaster and its sibling species, D. simulans. An unexpectedly high level of silent polymorphism in the Adh coding region relative to the 5' and 3' flanking regions in D. melanogaster is revealed by a populational survey of restriction polymorphism using a four-cutter filter hybridization technique as well as by direct sequence comparisons. In both of these studies, a region of the Adh gene encompassing the three coding exons exhibits a frequency of polymorphism equal to that of a 4-kb 5' flanking region. In contrast, an interspecific sequence comparison shows a two-fold higher level of divergence in the 5' flanking sequence compared to the structural locus. Analysis of the patterns of variation suggest an excess of polymorphism within the D. melanogaster Adh locus, rather than lack of polymorphism in the 5' flanking region. An approach is outlined for testing neutral theory predictions about patterns of variation within and between species. This approach indicates that the observed patterns of variation are incompatible with an infinite site neutral model.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.     

Molecular Evolution of Drosophila Metallothionein Genes

The metallothionein genes of Drosophila melanogaster, Mtn and Mto, may play an important role in heavy metal detoxification. Several different tandem duplications of Mtn have been shown to increase cadmium and copper tolerance, as well as Mtn expression. In order to investigate the possibility of increased selection for duplications of these genes in natural populations exposed to high levels of heavy metals, we compared the frequencies of such duplications among flies collected from metal-contaminated and non-contaminated orchards in Pennsylvania, Tennessee and Georgia. Restriction enzyme analysis was used to screen 1666 wild third chromosomes for Mtn duplications and a subset (327) of these lines for Mto duplications. The frequency of pooled Mtn duplications found ranged from 0% to 20%, and was not significantly higher at the contaminated sites. No Mto duplications were identified. Estimates of sequence diversity at the Mtn locus among a subsample (92) of the duplication survey were obtained using four-cutter analysis. This analysis revealed a low level of polymorphism, consistent with both selection at the Mtn locus, and a fairly recent origin for the duplications. To further examine this hypothesis, we sequenced an Mtn allele of Drosophila simulans and measured the amount of nucleotide sequence divergence between D. simulans and its sibling species D. melanogaster. The levels of silent nucleotide polymorphism and divergence in the Mtn region were compared with those in the Adh region, using the neutrality test of R.R. Hudson, M. Kreitman and M. Aguadé.     

Nucleotide variation of the Est-6 gene region in natural populations of Drosophila melanogaster

We have investigated nucleotide polymorphism in the Est-6 gene region in four samples of Drosophila melanogaster derived from natural populations of East Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). There are two divergent sequence types in the North and South American samples, which are not perfectly (North America) or not at all (South America) associated with the Est-6 allozyme variation. Less pronounced or no sequence dimorphism occurs in the European and African samples, respectively. The level of nucleotide diversity is highest in the African sample, lower (and similar to each other) in the samples from Europe and North America, and lowest in the sample from South America. The extent of linkage disequilibrium is low in Africa (1.23% significant associations), but much higher in non-African populations (22.59, 21.45, and 37.68% in Europe, North America, and South America, respectively). Tests of neutrality with recombination are significant in non-African samples but not significant in the African sample. We propose that demographic history (bottleneck and admixture of genetically different populations) is the major factor shaping the nucleotide patterns in the Est-6 gene region. However, positive selection modifies the pattern: balanced selection creates elevated levels of nucleotide variation around functionally important (target) polymorphic sites (RsaI-/RsaI+ in the promoter region and F/S in the coding region) in both African and non-African samples; and directional selection, acting during the geographic expansion phase of D. melanogaster, creates an excess of very similar sequences (RsaI- and S allelic lineages, in the promoter and coding regions, respectively) in the non-African samples.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

A high-frequency null mutant of an odorant-binding protein gene, Obp57e, in Drosophila melanogaster

We have found a null mutant of an odorant-binding protein, Obp57e, in Drosophila melanogaster. This frameshift mutation, which is a 10-bp deletion in the coding region, is at a high frequency in the Kyoto population and is also present in Taiwan and Africa. We have sequenced a 1.5-kb region including the tandemly duplicated gene, Obp57d, from 16 inbred lines sampled in Kyoto, Japan. The analyses showed a peak of nucleotide diversity and strong linkage disequilibrium around this mutation. This pattern suggests an elevated mutation rate or an influence of balancing selection in this region. The level of nucleotide divergence between D. melanogaster and D. simulans does not support the former possibility. Thus, this presence/absence polymorphism may be due to balancing selection, which takes advantage of the relatively weak functional constraint in members of a large gene family. In addition, the Obp57d gene region showed an excess of high-frequency-derived mutants that is consistent with a pattern predicted under positive natural selection.