Complete Mitochondrial DNA Sequence of Conger myriaster (Teleostei: Anguilliformes): Novel Gene Order for Vertebrate Mitochondrial Genomes and the Phylogenetic Implications for Anguilliform Families

The complete nucleotide sequence of the mitochondrial genome was determined for a conger eel, Conger myriaster (Elopomorpha: Anguilliformes), using a PCR-based approach that employs a long PCR technique and many fish-versatile primers. Although the genome [18,705 base pairs (bp)] contained the same set of 37 mitochondrial genes [two ribosomal RNA (rRNA), 22 transfer RNA (tRNA), and 13 protein-coding genes] as found in other vertebrates, the gene order differed from that recorded for any other vertebrates. In typical vertebrates, the ND6, tRNA^Glu, and tRNA^Pro genes are located between the ND5 gene and the control region, whereas the former three genes, in C. myriaster, have been translocated to a position between the control region and the tRNA^Phe gene that are contiguously located at the 5′ end of the 12S rRNA gene in typical vertebrates. This gene order is similar to the recently reported gene order in four lineages of birds in that the latter lack the ND6, tRNA^Glu, and tRNA^Pro genes between the ND5 gene and the control region; however, the relative position of the tRNA^Pro to the ND6–tRNA^Glu genes in C. myriaster was different from that in the four birds, which presumably resulted from different patterns of tandem duplication of gene regions followed by gene deletions in two distantly related groups of organisms. Sequencing of the ND5–cyt b region in 11 other anguilliform species, representing 11 families, plus one outgroup species, revealed that the same gene order as C. myriaster was shared by another 4 families, belonging to the suborder Congroidei. Although the novel gene orders of four lineages of birds were indicated to have multiple independent origins, phylogenetic analyses using nucleotide sequences from the mitochondrial 12S rRNA and cyt b genes suggested that the novel gene orders of the five anguilliform families had originated in a single ancestral species. Received: 13 July 2000 / Accepted: 30 November 2000     

Mitochondrial DNA of the coral sarcophyton glaucum contains a gene for a homologue of bacterial muts: A possible case of gene transfer from the nucleus to the mitochondrion

The nucleotide sequences of two segments of 6,737 ntp and 258 ntp of the 18.4-kb circular mitochondrial (mt) DNA molecule of the soft coral Sarcophyton glaucum (phylum Cnidaria, class Anthozoa, subclass Octocorallia, order Alcyonacea) have been determined. The larger segment contains the 3′ 191 ntp of the gene for subunit 1 of the respiratory chain NADH dehydrogenase (ND1), complete genes for cytochrome b (Cyt b), ND6, ND3, ND4L, and a bacterial MutS homologue (MSH), and the 5′ terminal 1,124 ntp of the gene for the large subunit rRNA (l-rRNA). These genes are arranged in the order given and all are transcribed from the same strand of the molecule. The smaller segment contains the 3′ terminal 134 ntp of the ND4 gene and a complete tRNA^f-Met gene, and these genes are transcribed in opposite directions. As in the hexacorallian anthozoan, Metridium senile, the mt-genetic code of S. glaucum is near standard: that is, in contrast to the situation in mt-genetic codes of other invertebrate phyla, AGA and AGG specify arginine, and ATA specifies isoleucine. However, as appears to be universal for metazoan mt-genetic codes, TGA specifies tryptophan rather than termination. Also, as in M. senile the mt-tRNA^f-Met gene has primary and secondary structural features resembling those of Escherichia coli initiator tRNA, including standard dihydrouridine and TψC loop sequences, and a mismatched nucleotide pair at the top of the amino-acyl stem. The presence of a mutS gene homologue, which has not been reported to occur in any other known mtDNA, suggests that there is mismatch repair activity in S. glaucum mitochondria. In support of this, phylogenetic analysis of MutS family protein sequences indicates that the S. glaucum mtMSH protein is more closely related to the nuclear DNA-encoded mitochondrial mismatch repair protein (MSH1) of the yeast Saccharomyces cerevisiae than to eukaryotic homologues involved in nuclear function, or to bacterial homologues. Regarding the possible origin of the S. glaucum mtMSH gene, the phylogenetic analysis results, together with comparative base composition considerations, and the absence of an MSH gene in any other known mtDNA best support the hypothesis that S. glaucum mtDNA acquired the mtMSH gene from nuclear DNA early in the evolution of octocorals. The presence of mismatch repair activity in S. glaucum mitochondria might be expected to influence the rate of evolution of this organism's mtDNA. Received: 13 January 1997 / Accepted: 23 September 1997     

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     

The Cephalopod Loligo bleekeri Mitochondrial Genome: Multiplied Noncoding Regions and Transposition of tRNA Genes

We previously reported the sequence of a 9260-bp fragment of mitochondrial (mt) DNA of the cephalopod Loligo bleekeri [J. Sasuga et al. (1999) J. Mol. Evol. 48:692–702]. To clarify further the characteristics of Loligo mtDNA, we have sequenced an 8148-bp fragment to reveal the complete mt genome sequence. Loligo mtDNA is 17,211 bp long and possesses a standard set of metazoan mt genes. Its gene arrangement is not identical to any other metazoan mt gene arrangement reported so far. Three of the 19 noncoding regions longer than 10 bp are 515, 507, and 509 bp long, and their sequences are nearly identical, suggesting that multiplication of these noncoding regions occurred in an ancestral Loligo mt genome. Comparison of the gene arrangements of Loligo, Katharina tunicata, and Littorina saxatilis mt genomes revealed that 17 tRNA genes of the Loligo mt genome are adjacent to noncoding regions. A majority (15 tRNA genes) of their counterparts is found in two tRNA gene clusters of the Katharina mt genome. Therefore, the Loligo mt genome (17 tRNA genes) may have spread over the genome, and this may have been coupled with the multiplication of the noncoding regions. Maximum likelihood analysis of mt protein genes supports the clade Mollusca + Annelida + Brachiopoda but fails to infer the relationships among Katharina, Loligo, and three gastropod species. Received: 9 May 2001 / Accepted: 3 October 2001     

Current Intraspecific Dynamics of Sequence Evolution Differs from Long-Term Trends and Can Account for the AT-Richness of Honeybee Mitochondrial DNA

The very high AT content of hymenopteran mtDNA has warranted speculation about nucleotide substitution processes in this group. Here we investigate the pattern of honeybee, Apis mellifera, mtDNA nucleotide polymorphisms inferred from phylogeny in terms of differences between the ATPase6, COI, COII, COIII, cytochrome b, and ND2 genes and strand asymmetry in mutation rates. The observed transition/transversion ratios and the distribution of nonsynonymous substitutions between regions differed significantly. The pattern of differences between genes leading to these heterogeneities (the ATPase6 and COIII genes group apart from the rest) differed markedly from that predicted on the basis of long-term evolutionary change and may indicate differences between current and long-term dynamics of sequence evolution. Also, there is strong strand asymmetry in substitutions, which probably results in a mutability of G and C sufficiently high to account for the AT-richness of honeybee mtDNA. Received: 21 October 1998 / Accepted: 27 January 1999     

Sequence Analysis of the Mitochondrial Genome of Sarcophyton glaucum: Conserved Gene Order Among Octocorals

The nucleotide sequence for an 11,715-bp segment of the mitochondrial genome of the octocoral Sarcophyton glaucum is presented, completing the analysis of the entire genome for this anthozoan member of the phylum Cnidaria. The genome contained the same 13 protein-coding and 2 ribosomal RNA genes as in other animals. However, it also included an unusual mismatch repair gene homologue reported previously and codes for only a single tRNA gene. Intermediate in length compared to two other cnidarians (17,443 and 18,911 bp), this organellar genome contained the smallest amount of noncoding DNA (428, compared to 1283 and 781 nt, respectively), making it the most compact one found for the phylum to date. The mitochondrial genes of S. glaucum exhibited an identical arrangement to that found in another octocoral, Renilla kolikeri, with five protein-coding genes in the same order as has been found in insect and vertebrate mitochondrial genomes. Although gene order appears to be highly conserved among octocorals, compared to the hexacoral, Metridium senile, few similarities were found. Like other metazoan mitochondrial genomes, the A + T composition was elevated and a general bias against codons ending in G or C was observed. However, an exception to this was the infrequent use of TGA compared to TGG to code for tryptophan. This divergent codon bias is unusual but appears to be a conserved feature among two rather distantly related anthozoans. Received: 27 January 1998 / Accepted: 25 May 1998     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.     

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     

The Mitochondrial Genome of Chlorogonium elongatum Inferred from the Complete Sequence

The 22,704-bp circular mitochondrial DNA (mtDNA) of the chlamydomonad alga Chlorogonium elongatum was completely cloned and sequenced. The genome encodes seven proteins of the respiratory electron transport chain, subunit 1 of the cytochrome oxidase complex (cox1), apocytochrome b (cob), five subunits of the NADH dehydrogenase complex (nad1, nad2, nad4, nad5, and nad6), a set of three tRNAs (Q, W, M), and the large (LSU)- and small (SSU)-subunit ribosomal RNAs. Six group-I introns were found, two each in the cox1, cob, and nad5 genes. In each intron an open reading frame (ORF) related to maturases or endonucleases was identified. Both the LSU and the SSU rRNA genes are split into fragments intermingled with each other and with other genes. Although the average A + T content is 62.2%, GC-rich clusters were detected in intergenic regions, in variable domains of the rRNA genes, and in introns and intron-encoded ORFs. A comparison of the genome maps reveals that C. elongatum and Chlamydomonas eugametos mtDNAs are more closely related to one another than either is to Chlamydomonas reinhardtii mtDNA. Received: 3 November 1997 / Accepted: 12 January 1998     

A Translocated Mitochondrial Cytochrome b Pseudogene in Voles (Rodentia: Microtus)

A full-length cytochrome b pseudogene was found in rodents; it has apparently been translocated from a mitochondrion to the nuclear genome in the subfamily Arvicolinae. The pseudogene (ψcytb) differed from its mitochondrial counterpart at 201 of 1143 sites (17.6%) and by four indels. Cumulative evidence suggests that the pseudogene has been translocated to the nucleus. Phylogenetic reconstruction indicates that the pseudogene arose before the diversification of M. arvalis/M. rossiaemeridionalis from M. oeconomus, but after the divergence of the peromyscine/sigmodontine/arvicoline clades some ∼10 MYA. Published rates of divergence between mitochondrial genes and their nuclear pseudogenes suggest that the translocation of this mitochondrial gene to the nuclear genome occurred some 6 MYA, in agreement with the phylogenetic evidence. Received: 16 January 1998 / Accepted: 18 July 1998     

Evolutionary Shifts in Three Major Structural Features of the Mitochondrial Genome Among Iguanian Lizards

A phylogenetic tree for major lineages of iguanian lizards is estimated from 1,488 aligned base positions (858 informative) of newly reported mitochondrial DNA sequences representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Two well-supported groups are defined, the Acrodonta and the Iguanidae (sensu lato). This phylogenetic hypothesis is used to investigate evolutionary shifts in mitochondrial gene order, origin for light-strand replication, and secondary structure of tRNA^Cys. These three characters shift together on the branch leading to acrodont lizards. Plate tectonics and the fossil record indicate that these characters changed in the Jurassic. We propose that changes to the secondary structure of tRNA^Cys may destroy function of the origin for light-strand replication which, in turn, may facilitate shifts in gene order. Received: 28 May 1996 / Accepted: 27 December 1996     

Mitochondrial DNA of Hydra attenuata (Cnidaria): A Sequence That Includes an End of One Linear Molecule and the Genes for l-rRNA, tRNAf-Met, tRNATrp, COII, and ATPase8

The 3231-nucleotide-pair (ntp) sequence of one end of one of the two linear mitochondrial (mt) DNA molecules of Hydra attenuata (phylum Cnidaria, class Hydrozoa, order Anthomedusae) has been determined. This segment contains complete genes for tRNA^f-Met, l-rRNA, tRNA^Trp, subunit 2 of cytochrome c oxidase (COII), subunit 8 of ATP synthetase (ATPase8), and the 5′ 136 ntp of ATPase6. These genes are arranged in the order given and are transcribed from the same strand of the molecule. As in two other cnidarians, the hexacorallian anthozoan Metridium senile and the octocorallian anthozoan Sarcophyton glaucum, the mt-genetic code of H. attenuata is near standard. The only modification appears to be that TGA specifies tryptophan rather than termination. Also as in M. senile and S. glaucum, the encoded H. attenuata mt-tRNA^f-Met has primary and secondary structural features resembling those of Escherichia coli initiator tRNA^t-Met. As the encoded mt-tRNA^Trp cannot be folded into a totally orthodox secondary structure, two alternative forms are suggested. The encoded H. attenuata mt-l-rRNA is 1738 nt, which is 451 nt shorter than the M. senile mt-l-rRNA. Comparisons of secondary structure models of these two mt-l-rRNAs indicate that most of the size difference results from loss of nucleotides in the H. attenuata molecule at a minimum of 46 locations, which includes elimination of six distinct helical elements. Received: 9 March 2000 / Accepted: 24 July 2000     

A Reexamination of the Phylogenetic Position of Callimico (Primates) Incorporating New Mitochondrial DNA Sequence Data

The New World monkeys are divided into two main groups, Callitrichidae and Cebidae. Callimico goeldii shares traits with both the Cebidae and the Callitrichidae. Recent morphological phyletic studies generally place Callimico as the most basal member of the Callitrichidae. In contrast, genetic studies (immunological, restriction fragment, and sequence data) have consistently placed Callimico somewhere within the Callitrichidae, not basal to this clade. A DNA sequence data set from the terminal 236 codons of the mitochondrial ND4 gene and the tRNA^His, tRNA^Ser, and tRNA^Leu genes was generated to clarify the position of Callimico. The sequences of 887 base pairs were analyzed by maximum-parsimony, neighbor-joining, and maximum-likelihood methods. The results of these various methods are generally congruent and place Callimico within the Callitrichidae between the marmosets (Callithrix and Cebuella) and the tamarins (Saguinus and Leontopithecus). Combined analyses of all suitable nuclear and mitochondrial gene sequences confirm the position of Callimico between the marmosets and the tamarins. As available molecular evidence indicates that Callimico is more closely related to the marmosets than to the tamarins, a reconsideration of the morphological evidence in light of the consensus tree from DNA sequence analyses is warranted. The marmosets and tamarins share four morphological characters (loss of the third molar, loss of the hypocone, reduced body size, reproductive twinning). Dwarfism may have evolved repeatedly among the Callitrichidae. It is well-known that the loss of a character can occur many times independently. The reproduction of marmosets and tamarins is extremely specialized and it is difficult to imagine that this complex and unique twinning system evolved separately in marmosets and tamarins. However, it is possible that a secondary reversal to single offspring took place in Callimico. Received: 20 March 1997 / Accepted: 17 December 1997     

Plant Mitochondrial RNA Editing

RNA editing affects messenger RNAs and transfer RNAs in plant mitochondria by site-specific exchange of cytidine and uridine bases in both seed and nonseed plants. Distribution of the phenomenon among bryophytes has been unclear since RNA editing has been detected in some but not all liverworts and mosses. A more detailed understanding of RNA editing in plants required extended data sets for taxa and sequences investigated. Toward this aim an internal region of the mitochondrial nad5 gene (1104 nt) was analyzed in a large collection of bryophytes and green algae (Charales). The genomic nad5 sequences predict editing in 30 mosses, 2 hornworts, and 7 simple thalloid and leafy liverworts (Jungermanniidae). No editing is, however, required in seven species of the complex thalloid liverworts (Marchantiidae) and the algae. RNA editing among the Jungermanniidae, on the other hand, reaches frequencies of up to 6% of codons being modified. Predictability of RNA editing from the genomic sequences was confirmed by cDNA analysis in the mosses Schistostega pennata and Rhodobryum roseum, the hornworts Anthoceros husnotii and A. punctatus, and the liverworts Metzgeria conjugata and Moerckia flotoviana. All C-to-U nucleotide exchanges predicted to reestablish conserved codons were confirmed. Editing in the hornworts includes the removal of genomic stop codons by frequent reverse U-to-C edits. Expectedly, no RNA editing events were identified by cDNA analysis in the marchantiid liverworts Ricciocarpos natans, Corsinia coriandra, and Lunularia cruciata. The findings are discussed in relation to models on the phylogeny of land plants. Received: 2 April 1998 / Accepted: 4 August 1998     

Mitochondrial Genes Collectively Suggest the Paraphyly of Crustacea with Respect to Insecta

Complete sequences of seven protein coding genes from Penaeus notialis mitochondrial DNA were compared in base composition and codon usage with homologous genes from Artemia franciscana and four insects. The crustacean genes are significantly less A + T-rich than their counterpart in insects and the pattern of codon usage (ratio of G + C-rich versus A + T-rich codon) is less biased. A phylogenetic analysis using amino acid sequences of the seven corresponding polypeptides supports a sister-taxon status for mollusks–annelid and arthropods. Furthermore, a distance matrix-based tree and two most-parsimonious trees both suggest that crustaceans are paraphyletic with respect to insects. This is also supported by the inclusion of Panulirus argus COII (complete) and COI and COIII (partial) sequence data. From analysis of single and combined genes to infer phylogenies, it is observed that obtained from single genes are not well supported in most topologies cases and notably differ from that of the tree based on all seven genes. Received: 25 August 1998 / Accepted: 8 March 1999     

Primary Structure and Phylogenetic Relationships of a Malate Dehydrogenase Gene from Giardia lamblia

The lactate and malate dehydrogenases comprise a complex protein superfamily with multiple enzyme homologues found in eubacteria, archaebacteria, and eukaryotes. In this study we describe the sequence and phylogenetic relationships of a malate dehydrogenase (MDH) gene from the amitochondriate diplomonad protist, Giardia lamblia. Parsimony, distance, and maximum-likelihood analyses of the MDH protein family solidly position G. lamblia MDH within a eukaryote cytosolic MDH clade, to the exclusion of chloroplast, mitochondrial, and peroxisomal homologues. Furthermore, G. lamblia MDH is specifically related to a homologue from Trichomonas vaginalis. This MDH topology, together with published phylogenetic analyses of β-tubulin, chaperonin 60, valyl-tRNA synthetase, and EF-1α, suggests a sister-group relationship between diplomonads and parabasalids. Since these amitochondriate lineages contain genes encoding proteins which are characteristic of mitochondria and α-proteobacteria, their shared ancestry suggests that mitochondrial properties were lost in the common ancestor of both groups. Received: 14 September 1998 / Accepted: 29 December 1998     

The complete mitochondrial genome of the marbled rockfish <Emphasis Type="Italic">Sebastiscus marmoratus</Emphasis> (Scorpaeniformes,Scorpaenidae): Genome characterization and phylogenetic considerations

The complete mitochondrial genome sequence of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae) was determined and phylogenetic analysis was conducted to elucidate the evolutionary relationship of the marbled rockfish with other Sebastinae species. This mitochondrial genome, consisting of 17301 bp, is highly similar to that of most other vertebrates, containing the same gene order and an identical number of genes or regions, including 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and one putative control region. Most of the genes are encoded on the H-strand, while the ND6 and seven tRNA genes (for Gln, Ala, Asn, Tyr, Ser (UCA), Glu, and Pro) are encoded on the L-strand. The reading frame of two pairs of genes overlapped on the same strand (the ATPase 8 and 6 genes overlapped by ten nucleotides; ND4L and ND4 genes overlapped by seven nucleotides). The possibly nonfunctional light-strand replication origin folded into a typical stem-loop secondary structure and a conserved motif (5′-GCCGG-3′) was found at the base of the stem within the tRNA^Cys gene. An extent termination-associated sequence (ETAS) and conserved sequence blocks (CSB) were identified in the control region, except for CSB-1; unusual long tandem repeats were found at the 3′ end of the control region. Phylogenetic analyses supported the view that Sebastinae comprises four genera (Sebastes, Hozukius, Helicolenus, and Sebastiscus).     

The Dentin Matrix Protein 1 Gene of Prototherian and Metatherian Mammals

Mineralization of tooth dentin (the deposition of hydroxyapatite crystals in and around collagen type I fibers of the extracellular matrix) requires the involvement of several genes, among them the gene coding for the dentin matrix protein 1, DMP1. We determined the exon–intron organization of the cattle DMP1 gene and used this information to amplify by the polymerase chain reaction homologous gene fragments from the genomic DNA of two species of metatherian (marsupial) mammals and one prototherian (monotreme) species. The translated proto- and metatherian protein sequences are highly divergent from the eutherian sequences but retain the general characteristics of the DMP1 (high acidity, serine-richness, multiple glycosylation sites, and the presence of the RGD cell attachment tripeptide). They therefore appear to be functional even though, evolutionarily, teeth are in a regression phase in prototherians. It is possible, therefore, that DMP1 is also involved in other functions besides dentinogenesis. The DMP1 gene appears to evolve rapidly and apparently tolerates non-frame-shifting insertions/deletions throughout the coding sequence. Received: 22 February 1998 / Accepted: 11 July 1998     

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