The evolution of gymnosperms redrawn by phytochrome genes: the Gnetatae appear at the base of the gymnosperms

Gymnosperms possess two to four phytochrome types which apparently are the result of successive gene duplications in the genomes of their common ancestors. Phytochromes are nuclear-encoded proteins whose genes, contrary to chloroplast, mitochondrion, and rRNA genes, have hitherto rarely been used to examine gymnosperm phylogenies. Since the individual phytochrome gene types implied phylogenies that were not completely congruent to one another, conflicting branching orders were sorted by the number of gene lineages present in a taxon. The Gnetatae (two gene types) branched at the base of all gymnosperms, a position supported by bootstrap sampling (distance and character state trees, maximum likelihood). The Gnetatae were followed by Ginkgo, Cycadatae, and Pinaceae (three gene types) and the remaining conifers (four gene types). Therefore, in phytochrome trees, the most ancient branch of the conifers (Pinatae) seems to be the Pinaceae. The next split appears to have separated Araucariaceae plus Podocarpaceae from the Taxaceae/Taxodiaceae/Cupressaceae group. Structural arrangements in the plastid genomes (Raubeson and Jansen 1992) corroborate the finding that there is no close connection between Pinaceae and Gnetatae as suggested by some publications. The analyses are based on 60 phytochrome genes (579 positions in an alignment of PCR fragments) from 28 species. According to rough divergence time estimates, the last common ancestor of gymnosperms and angiosperms is likely to have existed in the Carboniferous.     

Evolution of MADS-box gene induction by FLO/LFY genes

Some MADS-box genes function as floral homeotic genes. The Arabidopsis LFY gene is a positive regulator of floral homeotic genes, and homologs of the FLO/LFY gene family in other angiosperms and gymnosperms are likely to have a similar function. To investigate the origin of the floral homeotic gene regulatory cascade involving the FLO/LFY gene, FLO/LFY homologs were cloned from a leptosporangiate fern (Ceratopteris richardii), two eusporangiate ferns (Angiopteris lygodiifolia and Botrychium multifidum var. robustum), three fern allies (Psilotum nudum, Equisetum arvense, and Isoetes asiatica), and a moss (Physcomitrella patens). The FLO/LFY gene phylogenetic tree indicates that both duplication and loss of FLO/LFY homologs occurred during the course of vascular plant evolution. The expression patterns of the Ceratopteris LFY genes (CrLFY1 and 2) were assessed. CrLFY1 expression was prominent in tissues including shoot tips and circinate reproductive leaves, but very weak in other tissues examined. Expression of CrLFY2 was also prominent in tissues, including shoot tips and circinate reproductive leaves. These patterns of expression are dissimilar to that of any Ceratopteris MADS-box gene previously reported, suggesting that the induction of MADS-box genes by FLO/LFY is not established at the stage of ferns. Received: 4 January 2001 / Accepted: 28 February 2001     

The MADS-box floral homeotic gene lineages predate the origin of seed plants: Phylogenetic and molecular clock estimates

Flower development in angiosperms is controlled in part by floral homeotic genes, many of which are members of the plant MADS-box regulatory gene family. The evolutionary history of these developmental genes was reconstructed using 74 loci from 15 dicot, three monocot, and one conifer species. Molecular clock estimates suggest that the different floral homeotic gene lineages began to diverge from one another about 450–500 mya, around the time of the origin of land plants themselves. Received: 31 January 1997 / Accepted: 9 April 1997     

Evolution of the nad3-rps12 Gene Cluster in Angiosperm Mitochondria: Comparison of Edited and Unedited Sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms. Received: 24 January 1996 / Accepted: 5 June 1996     

Characterization of Mitochondrial Small-Subunit Ribosomal RNAs from Holoparasitic Plants

Mitochondrial small-subunit (19S) rDNA sequences were obtained from 10 angiosperms to further characterize sequence divergence levels and structural variation in this molecule. These sequences were derived from seven holoparasitic (nonphotosynthetic) angiosperms as well as three photosynthetic plants. 19S rRNA is composed of a conservative core region (ca. 1450 nucleotides) as well as two variable regions (V1 and V7). In pairwise comparisons of photosynthetic angiosperms to Glycine, the core 19S rDNA sequences differed by less than 1.4%, thus supporting the observation that variation in mitochondrial rDNA is 3–4 times lower than seen in protein coding and rDNA genes of other subcellular organelles. Sequences representing four distinct lineages of nonasterid holoparasites showed significantly increased numbers of substitutions in their core 19S rDNA sequences (2.3–7.6%), thus paralleling previous findings that showed accelerated rates in nuclear (18S) and plastid (16S) rDNA from the same plants. Relative rate tests confirmed the accelerated nucleotide substitution rates in the holoparasites whereas rates in nonparasitic plants were not significantly increased. Among comparisons of both parasitic and nonparasitic plants, transversions outnumbered transitions, in many cases more than two to one. The core 19S rRNA is conserved in sequence and structure among all nonparasitic angiosperms whereas 19S rRNA from members of holoparasitic Balanophoraceae have unique extensions to the V5 and V6 variable domains. Substitution and insertion/deletion mutations characterized the V1 and V7 regions of the nonasterid holoparasites. The V7 sequence of one holoparasite (Scybalium) contained repeat motifs. The cause of substitution rate increases in the holoparasites does not appear to be a result of RNA editing, hence the underlying molecular mechanism remains to be fully documented. Received: 18 May 1997 / Accepted: 11 July 1997     

Plant NADP-dependent isocitrate dehydrogenases are predominantly localized in the cytosol

The isoenzyme patterns of NADP-isocitrate dehydrogenase (NADP-IDH; EC 1.1.1.42) have been investigated in 15 species of higher plants using dietylaminoethyl ion-exchange chromatography and immunological techniques. The obtained results unambiguously demonstrate that the cytosolic enzyme is the predominant form in leaf extracts of all the surveyed plant species. The chloroplastic isoenzyme, previously reported in pea (Pisum sativum L.) leaves (R.D. Chen et al., 1989, Planta 178; 157–163), is a minor form in ferns and dicotyledonous angiosperms and is undetectable in gymnosperms and monocotyledonous angiosperms. Comparison of immunological relatedness suggests that the proteins of cytosolic isoenzymes have been highly conserved in the course of plant evolution. The data support the previously proposed idea that the cytosol is the major site for α-ketoglutarate production to be used for nitrogen assimilation. Received: 18 June 1998 / Accepted: 6 August 1998     

Complex Patterns of Plastid 16S rRNA Gene Evolution in Nonphotosynthetic Green Algae

This study provides a phylogenetic/comparative approach to deciphering the processes underlying the evolution of plastid rRNA genes in genomes under relaxed functional constraints. Nonphotosynthetic green algal taxa that belong to two distinct classes, Chlorophyceae (Polytoma) and Trebouxiophyceae (Prototheca), were investigated. Similar to the situation described previously for plastid 16S rRNA genes in nonphotosynthetic land plants, nucleotide substitution levels, extent of structural variations, and percentage AT values are increased in nonphotosynthetic green algae compared to their closest photosynthetic relatives. However, the mutational processes appear to be different in many respects. First, with the increase in AT content, more transversions are noted in Polytoma and holoparasite angiosperms, while more transitions characterize the evolution of the 16S rDNA sequences in Prototheca. Second, although structural variations do accumulate in both Polytoma and Prototheca (as well as holoparasitic plastid 16S rRNAs), insertions as large as 1.6 kb characterize the plastid 16S rRNA genes in the former, whereas significantly smaller indels (not exceeding 24 bp) seem to be more prevalent in the latter group. The differences in evolutionary rates and patterns within and between lineages might be due to mutations in replication/repair-related genes; slipped-strand mispairing is likely the mechanism responsible for the expansion of insertions in Polytoma plastid 16S rRNA genes. Received: 29 December 2000 / Accepted: 18 May 2001     

Variation in the psbC gene region of gymnosperms and angiosperms as detected by a single restriction site polymorphism

 This paper reports on a PCR-RFLP analysis in a chloroplast DNA region consisting of coding and intergenic spacer sequences of trnS and the adjacent psbC gene. This region was PCR-amplified in 62 woody plant species, predominantly tree species, that represent a broad systematic range in both gymnosperms and dicotyledonous angiosperms. The amplification products were digested by the restriction endonuclease HaeIII (GG↓CC). Fourteen different restriction patterns occurred, 5 of which characterised representatives of the gymnosperms, and 9 angiosperm representatives. A single restriction site polymorphism revealed most of the species to share restriction patterns. Groups formed which showed relationships to plant systematic units. This phenomenon is discussed with regard to the psbC gene and the GGCC motif for tracing species’ relationships on a high taxonomic level of gymnosperms and angiosperms. Received: 5 June 1996/Accepted: 20 September 1996     

Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996     

Ancient and Recent Horizontal Invasions of Drosophilids by P Elements

P elements of two different subfamilies designated as M- and O-type are thought to have invaded host species in the Drosophila obscura group via horizontal transmission from external sources. Sequence comparisons with P elements isolated from other species suggested that the horizontal invasion by the O-type must have been a rather recent event, whereas the M-type invasion should have occurred in the more distant past. To trace the phylogenetic history of O-type elements, additional taxa were screened for the presence of O- and M-type elements using type-specific PCR primers. The phylogeny deduced from the sequence data of a 927-bp section (14 taxa) indicate that O-type elements have undergone longer periods of regular vertical transmission in the lineages of the saltans and willistoni groups of Drosophila. However, starting from a species of the D. willistoni group they were transmitted horizontally into other lineages. First the lineage of the D. affinis subgroup was infected, and finally, in a more recent wave of horizontal spread, species of three different genera were invaded by O-type elements from the D. affinis lineage: Scaptomyza, Lordiphosa, and the sibling species D. bifasciata/D. imaii of the Drosophila obscura subgroup. The O-type elements isolated from these taxa are almost identical (sequence divergence <1%). In contrast, no such striking similarities are observed among M-type elements. Nevertheless, the sequence phylogeny of M-type elements is also not in accordance with the phylogeny of their host species, suggesting earlier horizontal transfer events. The results imply that P elements cross species barriers more frequently than previously thought but require a particular genomic environment and thus seem to be confined to a rather narrow spectrum of host species. Consequently, different P element types acquired by successive horizontal transmission events often coexist within the same genome. Received: 15 May 2000 / Accepted: 19 July 2000     

The Old World sparrows (genus Passer) phylogeography and their relative abundance of nuclear mtDNA pseudogenes

The phylogenetic relationships of genus Passer (Old World sparrows) have been studied with species covering their complete world living range. Mitochondrial (mt) cyt b genes and pseudogenes have been analyzed, the latter being strikingly abundant in genus Passer compared with other studied songbirds. The significance of these Passer pseudogenes is presently unclear. The mechanisms by which mt cyt b genes become pseudogenes after nuclear translocation are discussed together with their mode of evolution, i.e., transition/transversion mitochondrial ratio is decreased in the nucleus, as is the constraint for variability at the three codon positions. However, the skewed base composition according to codon position (in 1^st position the percentage is very similar for the four bases, in 2^nd position there are fewer percentage of A and G and more percentage of T, and in 3^rd codon position fewer percentage of G and T and is very rich in A and C) is maintained in the translocated nuclear pseudogenes. Different nuclear internal mechanisms and/or selective pressures must exist for explaining this nuclear/mitochondrial differential DNA base evolutive variability. Also, the phylogenetic usefulness of pseudogenes for defining relationships between closely related lineages is stressed. The analyses suggest that the primitive genus Passer species comes from Africa, the Cape sparrow being the oldest: P. hispaniolensis italiae is more likely conspecific to P. domesticus than to P. hispaniolensis. Also, Passer species are not included within weavers or Estrildinae or Emberizinae, as previously suggested. European and American Emberizinae sparrows are closely related to each other and seem to be the earliest species that radiated among the studied songbirds (all in the Miocene Epoch). Received: 29 November 2000 / Accepted: 22 March 2001     

Phylogeny, Genome Evolution, and Host Specificity of Single-Stranded RNA Bacteriophage (Family Leviviridae)

Bacteriophage of the family Leviviridae have played an important role in molecular biology where representative species, such as Qβ and MS2, have been studied as model systems for replication, translation, and the role of secondary structure in gene regulation. Using nucleotide sequences from the coat and replicase genes we present the first statistical estimate of phylogeny for the family Leviviridae using maximum-likelihood and Bayesian estimation. Our analyses reveal that the coliphage species are a monophyletic group consisting of two clades representing the genera Levivirus and Allolevivirus. The Pseudomonas species PP7 diverged from its common ancestor with the coliphage prior to the ancient split between these genera and their subsequent diversification. Differences in genome size, gene composition, and gene expression are shown with a high probability to have changed along the lineage leading to the Allolevivirus through gene expansion. The change in genome size of the Allolevivirus ancestor may have catalyzed subsequent changes that led to their current genome organization and gene expression. Received: 3 March 2000 / Accepted: 17 October 2000     

Divergent Evolution of Plant NBS-LRR Resistance Gene Homologues in Dicot and Cereal Genomes

The majority of plant disease resistance genes are members of very large multigene families. They encode structurally related proteins containing nucleotide binding site domains (NBS) and C-terminal leucine rich repeats (LRR). The N-terminal region of some resistance genes contain a short sequence called TIR with homology to the animal innate immunity factors, Toll and interleukin receptor-like genes. Only a few plant resistance genes have been functionally analyzed and the origin and evolution of plant resistance genes remain obscure. We have reconstructed gene phylogeny by exhaustive analysis of available genome and amplified NBS domain sequences. Our study shows that NBS domains faithfully predict whole gene structure and can be divided into two major groups. Group I NBS domains contain group-specific motifs that are always linked with the TIR sequence in the N terminus. Significantly, Group I NBS domains and their associated TIR domains are widely distributed in dicot species but were not detected in cereal databases. Furthermore, Group I specific NBS sequences were readily amplified from dicot genomic DNA but could not be amplified from cereal genomic DNA. In contrast, Group II NBS domains are always associated with putative coiled-coil domains in their N terminus and appear to be present throughout the angiosperms. These results suggest that the two main groups of resistance genes underwent divergent evolution in cereal and dicot genomes and imply that their cognate signaling pathways have diverged as well. Received: 17 May 1999 / Accepted: 25 September 1999     

Evolutionary History of 4.5SI RNA and Indication That It Is Functional

To date, the small nuclear 4.5S_I RNA has only been studied in the rat (Rattus norvegicus). Combining PCR and hybridization analyses, we have revealed 4.5S_I RNA homologues sequences in the genomes of four myomorph rodent families (Muridae, Cricetidae, Spalicidae, and Rhizomyidae), and not in other myomorph families (Dipodidae, Zapodidae, Geomyidae, and Heteromyidae) or sciuromorph and caviomorph rodents. By Northern-hybridization, 4.5S_I RNA has been detected in the common rat (R. norvegicus, Muridae), golden hamster (Mesocricetus auratus, Cricetidae), and Russian mole rat (Spalax microphthalmus, Spalacidae), but not in the related great jerboa (Allactaga jaculus, Dipodidae) or in four non-myomorph rodent species tested. cDNA derived from 4.5S_I RNA of M. auratus and S. microphthalmus has been cloned and sequenced. The hamster RNA is found to differ from rat 4.5S_I RNA by only one nucleotide substitution. For the mole rat, two variants of 4.5S_I RNA are detected: short (S) and long (L) with length 101 and 108 nt, respectively. The L variant differs from the S variant as well as from murid and cricetid 4.5S_I RNAs by both a 7 nt insertion and a varying number of nucleotide substitutions. The sequence similarity between the spalacid S-variant and murid/crecitid variants of 4.5S_I RNA is 90%. Judging from species distribution, 4.5S_I RNA genes emerged during the same period of time as the related short interspersed element B2 arose. This occurred after the divergence of Dipodidae lineage but before the branching of Spalicidae/Rhizomyidae lineage from a common myomorph rodent stem. S variant genes seemed to emerge in a common ancestor of spalacids and rhizomyds whereas L variant genes formed in spalacids following the divergence of these two families. The low rate of evolutionary changes of 4.5S_I RNA, at least, in murids and cricetids (6 × 10⁻⁴ substitutions per site per million years), suggests that this RNA is under selection constraint and have a function. This is a remarkable fact if the recent origin and narrow species distribution range of 4.5S_I RNA genes is taken into account. Genes with narrow species distribution are proposed to be referred to as stenogenes. Received: 11 December 2000 / Accepted: 27 August 2001     

Sequence Peculiarity of Gnetalean Legumin-Like Seed Storage Proteins

The development of seeds as a specialized organ for the nutrition, protection, and dispersal of the next generation was an important step in the evolution of land plants. Seed maturation is accompanied by massive synthesis of storage compounds such as proteins, starch, and lipids. To study the processes of seed storage protein evolution we have partially sequenced storage proteins from maturing seeds of representatives from the gymnosperm genera Gnetum, Ephedra, and Welwitschia—morphologically diverse and unusual taxa that are grouped in most formal systems into the common order Gnetales. Based on partial N-terminal amino acid sequences, oligonucleotide primers were derived and used for PCR amplification and cloning of the corresponding cDNAs. We also describe the structure of the nuclear gene for legumin of Welwitschia mirabilis. This first gnetalean nuclear gene structure contains introns in only two of the four conserved positions previously characterized in other spermatophyte legumin genes. The distinct phylogenetic status of the gnetalean taxa is also reflected in a sequence peculiarity of their legumin genes. A comparative analysis of exon/intron sequences leads to the hypothesis that legumin genes from Gnetales belong to a monophyletic evolutionary branch clearly distinct from that of legumin genes of extant Ginkgoales and Coniferales as well as from all angiosperms. Received: 5 June 1997 / Accepted: 31 March 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 Position of the Subgenus Lordiphosa of the Genus Drosophila (Diptera: Drosophilidae) Inferred from Alcohol Dehydrogenase (Adh) Gene Sequences

We analyzed the phylogenetic relationship between the species of Lordiphosa and other Drosophilidae using alcohol dehydrogenase (Adh) gene sequences. The phylogenetic trees consistently show that the four species Drosophila kurokawai, D. collinella, D. stackelbergi, and D. clarofinis, which include three species groups of Lordiphosa, form a monophyletic clade. This clade is placed as a sister group to the willistoni and saltans groups of Sophophora. On the other hand, three species of Lordiphosa, D. tenuicauda, D. pseudotenuicauda, and D. acutissima, all of which belong to the tenuicauda group, are not shown to be related to the major Lordiphosa lineage. In the phylogenetic trees, these species are included into the clade comprised of Drosophila and Hirtodrosophila, although it remains uncertain whether the tenuicauda group is a monophyletic group or not. These results indicate that Lordiphosa is polyphyletic and that most of the members of the subgenus have a close relationship to the neotropical groups of Sophophora. The above conclusion is compatible with the hypothesis of Okada (Mushi [1963] 37:79–100) and Lastovka and Máca (Acta Ent Bohemoslov [1978] 75:404–420) that Lordiphosa is most closely related to Sophophora; in contrast, our results contradict the hypothesis of Grimaldi (Bull Am Mus Nat Hist [1990] 197:1–139) that Lordiphosa is a sister group to the genus Scaptomyza. Received: 12 May 1999 / Accepted: 14 April 2000     

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     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996