Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.