Multiple Recurrent Evolution of Trophic Types in Northeastern Atlantic and Mediterranean Seabreams (Sparidae, Percoidei)

Seabreams are among the most valuable fish, not only for small-scale and semiindustrial fisheries but also for aquaculture throughout the Mediterranean. Nevertheless, their phylogenetic relationships are not at all clear. The current taxonomy is based solely on trophic morphology and rests on the assumption that each trophic type evolved only once from a less specialized ancestral condition. We analyzed a 486-bp segment of the mitochondrial 16S rDNA of all 24 seabream species described for the northeastern Atlantic and the Mediterranean to elucidate their generic and subfamily-level relationships. Three major mitochondrial lineages, each comprising species of different feeding strategy and dentition, were found that do not agree with the present taxonomic assignments. Most of the investigated genera were resolved paraphyletically, indicating that the structure and arrangement of oral teeth must have repeatedly evolved from a less specialized ancestral condition. Further, the genus Sparus was resolved as distantly related to the genus Pagrus, in that it was assigned to a different major mitochondrial lineage. Oblada melanura was consistently placed within the Diplodus radiation as sister group to Diplodus puntazzo. Our phylogenetic hypothesis thus suggests multiple independent origins of similar trophic specializations within the Sparidae and indicates that the currently recognized three or four subfamilies need to be redefined. Received: 5 October 1999 / Accepted: 9 November 1999     

Partial Sequence of a Sponge Mitochondrial Genome Reveals Sequence Similarity to Cnidaria in Cytochrome Oxidase Subunit II and the Large Ribosomal RNA Subunit

A 2550-bp portion of the mitochondrial genome of a Demosponge, genus Tetilla, was amplified from whole genomic DNA extract and sequenced. The sequence was found to code for the 3′ end of the 16S rRNA gene, cytochrome c oxidase subunit II, a lysine tRNA, ATPase subunit 8, and a 5′ portion of ATPase subunit 6. The Porifera cluster distinctly within the eumetazoan radiation, as a sister group to the Cnidaria. Also, the mitochondrial genetic code of this sponge is likely identical to that found in the Cnidaria. Both the full COII DNA and protein sequences and a portion of the 16S rRNA gene were found to possess a striking similarity to published Cnidarian mtDNA sequences, allying the Porifera more closely to the Cnidaria than to any other metazoan phylum. The gene arrangement, COII—tRNA^Lys—ATP8—ATP6, is observed in many Eumetazoan phyla and is apparently ancestral in the metazoa. Received: 24 November 1997 / Accepted: 14 September 1998     

Extreme Endemic Radiation of the Malagasy Vangas (Aves: Passeriformes)

Phylogenetic relationships of the family Vangidae and representatives of several other passeriform families were inferred from 882 base positions of mitochondrial DNA sequences of 12S and 16S rRNA genes. Results indicated the monophyly of the Vangidae, which includes the genus Tylas, hitherto often placed in the family Pycnonotidae. Our results also revealed the Malagasy endemic Newtonia, a genus never previously assigned to the Vangidae, to be a member of this family. These results suggest the occurrence of an extensive in situ radiation of this family within Madagascar, and that the extant high diversity of this family is not the result of multiple colonizations from outside. The extremely high morphological and ecological diversification of the family seems to have been enhanced through the use and ultimate occupancy of vacant niches in this island. Received: 8 September 2000 / Accepted: 13 February 2001     

Multiple Substitutions Affect the Phylogenetic Utility of Cytochrome b and 12S rDNA Data: Examining a Rapid Radiation in Leporid (Lagomorpha) Evolution

Partial sequences of two mitochondrial genes, the 12S ribosomal gene (739 bp) and the cytochrome b gene (672 bp), were analyzed in hopes of reconstructing the evolutionary relationships of 11 leporid species, representative of seven genera. However, partial cytochrome b sequences were of little phylogenetic value in this study. A suite of pairwise comparisons between taxa revealed that at the intergeneric level, the cytochrome b gene is saturated at synonymous coding positions due to multiple substitution events. Furthermore, variation at the nonsynonymous positions is limited, rendering the cytochrome b gene of little phylogenetic value for assessing the relationships between leporid genera. If the cytochrome b data are analyzed without accounting for these two classes of nucleotides (i.e., synonymous and nonsynonymous sites), one may incorrectly conclude that signal exists in the cytochrome b data. The mitochondrial 12S rRNA gene, on the other hand, has not experienced excessive saturation at either stem or loop positions. Phylogenies reconstructed from the 12S rDNA data support hypotheses based on fossil evidence that African rock rabbits (Pronolagus) are outside of the main leporid stock and that leporids experienced a rapid radiation. However, the molecular data suggest that this radiation event occurred in the mid-Miocene several millions of years earlier than the Pleistocene dates suggested by paleontological evidence. Received: 23 April 1998 / Accepted: 14 May 1998     

Mitochondrial DNA Phylogeny of the Family Cichlidae: Monophyly and Fast Molecular Evolution of the Neotropical Assemblage

A mitochondrial DNA (mtDNA) phylogeny of cichlid fish is presented for the most taxonomically inclusive data set compiled to date (64 taxa). 16S rDNA data establish with confidence relationships among major lineages of cichlids, with a general pattern congruent with previous morphological studies and less inclusive molecular phylogenies based on nuclear genes. Cichlids from Madagascar and India are the most basal groups of the family Cichlidae and sister to African–Neotropical cichlids. The cichlid phylogeny suggests drift-vicariance events, consistent with the fragmentation of Gondwana, to explain current biogeographic distributions. Important phylogenetic findings include the placement of the controversial genus Heterochromis basal among African cichlids, the South American genus Retroculus as the most basal taxon of the Neotropical cichlid assemblage, and the close relationship of the Neotropical genera Cichla with Astronotus rather than with the crenicichlines. Based on a large number of South American genera, the Neotropical cichlids are defined as a monophyletic assemblage and shown to harbor significantly higher levels of genetic variation than their African counterparts. Relative rate tests suggest that Neotropical cichlids have experienced accelerated rates of molecular evolution. But these high evolutionary rates were significantly higher among geophagine cichlids. Received: 18 September 1998 / Accepted: 16 December 1998     

Sequence of PRAT Satellite DNA ``Frozen' in Some Coleopteran Species

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera (Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50–60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.     

Sociality and the Rate of rDNA Sequence Evolution in Wasps (Vespidae) and Honeybees (Apis)

Sequence data of mitochondrial 16S ribosomal DNA (mt-rDNA) and nuclear 28S ribosomal DNA (nuc-rDNA) were compared in two honeybee species (Apis mellifera and Apis dorsata) and a selection of 22 wasp species (Vespidae) with different levels of sociality. The averge substitution rates in mt-rDNA and nuc-rDNA were almost-equal in solitary species. In species with larger nests, however, the difference between the nuclear and the mitochondrial substitution rate significantly increased. The average substitution ratio, ψ (nucleotide substitutions in mt-rDNA/nucleotide substitutions in nuc-rDNA) was 1.48 ± 0.12 (SE) among the solitary Eumeninae, 3.70 ± 0.15 among five primitive social Stenogastrinae species, 3.24 ± 0.20 among five Polistinae species, 5.76 ± 0.33 among nine highly eusocial Vespinae, and 12.7 in the two Apis species. The high egg-laying rate and the effective population size skew between the sexes may contribute to the rise of the substitution ratio in the highly eusocial species. Drift and bottleneck effects in the mitochondrial DNA pool during speciation events as well as polyandry may further enhance this phenomenon. Received: 12 January 1998 / Accepted: 28 April 1998     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998     

The Peperomia Mitochondrial coxI Group I Intron: Timing of Horizontal Transfer and Subsequent Evolution of the Intron

The Peperomia polybotrya coxI gene intron is the only currently reported group I intron in a vascular plant mitochondrial genome and it likely originated by horizontal transfer from a fungal donor. We provide a clearer picture of the horizontal transfer and a portrayal of the evolution of the group I intron since it was gained by the Peperomia mitochondrial genome. The intron was transferred recently in terms of plant evolution, being restricted to the single genus Peperomia among the order Piperales. Additional support is presented for the suggestion that a recombination/repair mechanism was used by the intron for integration into the Peperomia mitochondrial genome, as a perfect 1:1 correspondence exists between the intron's presence in a species and the presence of divergent nucleotide markers flanking the intron insertion site. Sequencing of coxI introns from additional Peperomia species revealed that several mutations have occurred in the intron since the horizontal transfer, but sequence alterations have not caused frameshifts or created stop codons in the intronic open reading frame. In addition, two coxI pseudogenes in Peperomia cubensis were discovered that lack a large region of coxI exon 2 and contain a truncated version of the group I intron that likely cannot be spliced out. Received: 29 May 1997 / Accepted: 1 November 1997     

Gene Order Breakpoint Evidence in Animal Mitochondrial Phylogeny

Multiple genome rearrangement methodology facilitates the inference of animal phylogeny from gene orders on the mitochondrial genome. The breakpoint distance is preferable to other, highly correlated but computationally more difficult, genomic distances when applied to these data. A number of theories of metazoan evolution are compared to phylogenies reconstructed by ancestral genome optimization, using a minimal total breakpoints criterion. The notion of unambiguously reconstructed segments is introduced as a way of extracting the invariant aspects of multiple solutions for a given ancestral genome; this enables a detailed reconstruction of the evolution of non-tRNA mitochondrial gene order. Received: 15 July 1998 / Accepted: 5 March 1999     

Phylogeny of Japanese Papilionid Butterflies Inferred from Nucleotide Sequences of the Mitochondrial ND5 Gene

Phylogenetic relationships among the Japanese papilionid butterflies were analyzed by comparing 783 nucleotide sequences of the mitochondrial gene encoding NADH dehydrogenase subunit 5 (ND5). Phylogenetic trees of the representative species from each family in the superfamily Papilionoidea revealed that the species of the family Papilionidae and those of all other families formed distinct clusters, with a few species of the family Hesperiidae (Hesperioidea) as an outgroup. In the phylogenetic trees of most Japanese species of the family Papilionidae with Nymphalis xanthomelas (Nymphalidae) as an outgroup, the tribe Parnassiini (Parnassiinae) formed a cluster, and the rest formed the other cluster in which the tribe Zerynthiini (Parnassiinae) and the subfamily Papilioninae formed different subclusters. In the Papilioninae cluster, the tribes Troidini and Graphiini formed a subcluster, and the tribe Papilionini formed the other subcluster. These results generally agree with the traditional classification of the papilionid butterflies based on their morphological characteristics and support the proposed evolutionary genealogy of the butterflies based on their morphology, behavior, and larval host plants, except that the tribes Parnasiini and Zerynthiini (both Parnassiinae) are not in the same cluster. Received: 16 March 1998 / Accepted: 28 April 1998     

Molecular Data from the 16S rRNA Gene for the Phylogeny of Pectinidae (Mollusca: Bivalvia)

The phylogenetic relationships among the species belonging to the family Pectinidae are still an issue of debate. The mitochondrial DNA sequences from the large ribosomal RNA gene may be of great value for systematic and phylogenetic studies within families. Partial sequences of the 16S rRNA gene were obtained for the scallop species Adamussium colbecki, Aequipecten opercularis, Chlamys glabra, C. islandica, C. varia, and Pecten jacobeus and compared with the published sequence of Pecten maximus. The present molecular data show that Chlamys are polyphyletic and do not support the assignment of these species to the two subfamilies Chlamydinae and Pectininae. Moreover, the minimal genetic distance between P. maximus and P. jacobeus suggests that they could belong to the same species. Received: 24 May 1999 / Accepted: 1 September 1999     

Phylogeny and Biogeography of Wood-Feeding Cockroaches, Genus Salganea Stål (Blaberidae: Panesthiinae), in Southeast Asia Based on Mitochondrial DNA Sequences

Molecular phylogenetic relationships among 25 species of the wood-feeding cockroach belonging to the genus Salganea St?l (Panesthiinae; Blaberidae) in Southeast Asia were analyzed based on the DNA sequence of the complete mitochondrial cytochrome oxidase II (COII) gene. Most basal relationships among species of Salganea are poorly resolved by both neighbor-joining and nonweighted parsimony analyses, suggesting the possibility of a hard polytomy due to a rapid and potentially simultaneous radiation early in the history of the genus. For more apical relationships, however, some interesting phylogenetic relationships were recognized. The monophyly of the two species groups, morio and foveolata, the former of which is distributed mainly in the Sunda lands (containing the Malay Peninsula, Sumatra, Java, and Borneo), whereas the latter is Sulawesi endemic, was strongly supported. Based on the inferred phylogenetic patterns and recent palaeogeographic scenario for Southeast Asia, it is suggested that a radiation of Salganea species occurred in Southeast Asia presumably in the early Tertiary, and several barriers against dispersal and gene flow, such as the formation of straits or high mountains, have arisen from the middle Tertiary. Received: 4 April 2001 / Accepted: 20 April 2001     

Neutral and Nonneutral Mitochondrial Genetic Variation in Deep-Sea Clams from the Family Vesicomyidae

Nucleotide sequences at two mitochondrial genes from 57 individuals representing eight species of deep-sea clams (Vesicomyidae) were examined for variation consistent with the neutral model of molecular evolution. One gene, cytochrome oxidase subunit I (COI), deviated from the expectations of neutrality by containing an excess of intraspecific nonsynonymous polymorphism. Additionally, one species, Calyptogena kilmeri, showed a significant excess of rare polymorphism specifically at the COI locus. In contrast, a second mitochondrial gene, the large-subunit 16S ribosomal RNA gene (16S), showed little deviation from neutrality either between or within species. Together, COI and 16S show no deviation from neutral expectations by the HKA test, produce congruent phylogenetic relationships between species, and show correlated numbers of fixed differences between species and polymorphism within species. These patterns of both neutral and nonneutral evolution within the mitochondrial genome are most consistent with a model where intraspecific nonsynonymous polymorphism at COI is near neutrality. In addition to examining the forces of molecular evolution, we extend hypotheses about interspecific relationships within this family for geographical locations previously unexamined by molecular methods including habitats near the Middle Atlantic, the Aleutian Trench, and Costa Rica. Received: 10 March 1999 / Accepted: 13 September 1999     

Evolution of trophic types in emperor fishes ( Lethrinus,Lethrinidae, Percoidei) based on cytochrome B gene sequence variation

Three trophic categories exist within emperor fishes, genus Lethrinus, relating to body form and dentition type. One group contains low-bodied, high speed, stalking predators with conical teeth. Another group comprises high-bodied, slow speed carnivores with molariform teeth capable of crushing hard-shelled benthic prey. A third group is also high-bodied but with conical teeth feeding mostly on small or soft-shelled benthic prey. Inferring the evolution of these trophic types within Lethrinus using morphology is problematic since these characters are typically correlated with feeding mode and are potentially homoplasious. We use mitochondrial DNA sequences, to independently determine a phylogenetic hypothesis for Lethrinus, which are not dependent on morphological characters relating to trophic categories. We analyzed complete cytochrome b gene sequences (1140 bp) for 20 species of Lethrinus, representing the three trophic types, and for 13 outgroup species, including four other representatives of the Lethrinidae. A monophyletic Lethrinidae did not resolve, but the monophyly of Lethrinus is well supported. In addition, two major clades within Lethrinus are well supported. One of these clades exclusively contains low-bodied species with conical teeth while the other clade only comprises the high-bodied species with molariform teeth. A high-bodied species with conical teeth, Lethrinus miniatus, appears most ancestral and sister to all other Lethrinus species. We hypothesize that this generalist trophic type was the evolutionary precursor to both of the other primary trophic types.     

Identification of Multiple Gypsy LTR-Retrotransposon Lineages in Vertebrate Genomes

Gypsy LTR-retrotransposons have been identified in the genomes of many organisms, but only a small number of vertebrate examples have been reported to date. Here we show that members of this family are likely to be widespread in many vertebrate classes with the possible exceptions of mammals and birds. Phylogenetic analyses demonstrate that although there are several distinct lineages of vertebrate gypsy LTR-retrotransposons, the majority clusters into one monophyletic clade. Groups of fungal, plant, and insect elements were also observed, suggesting horizontal transfer between phyla may be infrequent. However, in contrast to this, there was little evidence to support sister relationships between elements derived from vertebrate and insect hosts. In fact, the majority of the vertebrate elements appeared to be most closely related to a group of gypsy LTR-retrotransposons present within fungi. This implies either that at least one horizontal transmission between these two phyla has occurred previously or that a gypsy LTR-retrotransposon lineage has been lost from insect taxa. Received: 22 December 1998 / Accepted: 6 April 1999     

Characterization of a Novel Class of Interspersed LTR Elements in Primate Genomes: Structure, Genomic Distribution, and Evolution

Retrovirus-like sequences and their solitary (solo) long terminal repeats (LTRs) are common repetitive elements in eukaryotic genomes. We reported previously that the tandemly arrayed genes encoding U2 snRNA (the RNU2 locus) in humans and apes contain a solo LTR (U2-LTR) which was presumably generated by homologous recombination between the two LTRs of an ancestral provirus that is retained in the orthologous baboon RNU2 locus. We have now sequenced the orthologous U2-LTRs in human, chimpanzee, gorilla, orangutan, and baboon and examined numerous homologs of the U2-LTR that are dispersed throughout the human genome. Although these U2-LTR homologs have been collectively referred to as LTR13 in the literature, they do not display sequence similarity to any known retroviral LTRs; however, the structure of LTR13 closely resembles that of other retroviral LTRs with a putative promoter, polyadenylation signal, and a tandemly repeated 53-bp enhancer-like element. Genomic blotting indicates that LTR13 is primate-specific; based on sequence analysis, we estimate there are about 2,500 LTR13 elements in the human genome. Comparison of the primate U2-LTR sequences suggests that the homologous recombination event that gave rise to the solo U2-LTR occurred soon after insertion of the ancestral provirus into the ancestral U2 tandem array. Phylogenetic analysis of the LTR13 family confirms that it is diverse, but the orthologous U2-LTRs form a coherent group in which chimpanzee is closest to the humans; orangutan is a clear outgroup of human, chimpanzee, and gorilla; and baboon is a distant relative of human, chimpanzee, gorilla, and orangutan. We compare the LTR13 family with other known LTRs and consider whether these LTRs might play a role in concerted evolution of the primate RNU2 locus. Received: 29 September 1997 / Accepted: 16 January 1998     

Evolution of Orthologous Intronless and Intron-Bearing Globin Genes in Two Insect Species

While globin genes ctt-2β and ctt-9.1 in Chironomus thummi thummi each have a single intron, all of the other insect globin genes reported so far are intronless. We analyzed four globin genes linked to the two intron-bearing genes in C. th. thummi. Three have a single intron at the same position as ctt-2β and ctt-9.1; the fourth is intronless and lies between intron bearing genes. Finally, in addition to its intron, one gene (ctt-13RT) was recently interrupted by retrotransposition. Phylogenetic analyses show that the six genes in C. th. thummi share common ancestry with five globin genes in the distantly related species C. tentans, and that a 5-gene ancestral cluster predates the divergence of the two species. One gene in the ancestral cluster gave rise to ctn-ORFB in C. tentans, and duplicated in C. th. thummi to create ctt-11 and ctt-12. From parsimonious calculations of evolutionary distances since speciation, ctt-11, ctt-12, and ctn-ORFB evolved rapidly, while ctn-ORFE in C. tentans evolved slowly compared to other globin genes in the clusters. While these four globins are under selective pressure, we suggest that most chironomid globin genes were not selected for their unique function. Instead, we propose that high gene copy number itself was selected because conditions favored organisms that could synthesize more hemoglobin. High gene copy number selection to produce more of a useful product may be the basis of forming multigene families, all of whose members initially accumulate neutral substitutions while retaining essential function. Maintenance of a large family of globin genes not only ensured high levels of hemoglobin production, but may have facilitated the extensive divergence of chironomids into as many as 5000 species. Received: 31 December 1996 / Accepted: 16 May 1997     

Phylogenetic Relationships Within the Genus Equus and the Evolution of α and θ Globin Genes

Sequences of the α1, α2 and θ globin genes from six equid species have been determined to investigate relationships within the genus Equus. Analyses using standard phylogenetic methods, or an approach designed to account for the effects of gene conversion between the α genes, gave broadly similar results and show that the horses diverged from the zebra/ass ancestor ∼2.4 million years ago and that the zebra and ass species arose in a rapid radiation ∼0.9 million years ago. These results from the α genes are corroborated by θ gene data and are in contrast to mitochondrial DNA studies of the phylogeny of this genus, which suggest a more gradual set of speciation events. Received: 22 April 1997 / Accepted: 20 July 1998     

Recent Selection on Synonymous Codon Usage in Drosophila

Evidence from a variety of sources indicates that selection has influenced synonymous codon usage in Drosophila. It has generally been difficult, however, to distinguish selection that acted in the distant past from ongoing selection. However, under a neutral model, polymorphisms usually reflect more recent mutations than fixed differences between species and may, therefore, be useful for inferring recent selection. If the ancestral state is preferred, selection should shift the frequency distribution of derived states/site toward lower values; if the ancestral is unpreferred, selection should increase the number of derived states/site. Polymorphisms were classified as ancestrally preferred or unpreferred for several genes of D. simulans and D. melanogaster. A computer simulation of coalescence was employed to derive the expected frequency distributions of derived states/site under various modifications of the Wright–Fisher neutral model, and distributions of test statistics (t and Mann–Whitney U) were derived by appropriate sampling. One-tailed tests were applied to transformed frequency data to assess whether the two frequency distributions deviated from neutral expectations in the direction predicted by selection on codon usage. Several genes from D. simulans appear to be subject to recent selection on synonymous codons, including one gene with low codon bias, esterase-6. Selection may also be acting in D. melanogaster. Received: 15 April 1998 / Accepted: 13 May 1999