*Department of Biology, University of Michigan;
Molecular Parasitology, Humboldt-University at Berlin, Berlin, Germany; and
DOE Joint Genome Institute and Lawrence Livermore National Laboratory, Walnut Creek, California
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although metazoan mitochondrial sequences are known to evolve rapidly, their gene arrangements often remain unchanged over long periods of evolutionary time. Comparison of the relative arrangements of these genes has several merits for reconstructing ancient evolutionary relationships, including their generally slow rate of change and the large number of potential gene orders (Boore and Brown 1998
). With some notable exceptions, gene arrangements are relatively stable within major groups but variable between them, and arrangement comparisons have resolved several evolutionary relationships, including those among arthropods (Boore et al. 1995
; Boore, Lavrov, and Brown 1998
; Dowton 1999
), annelids (Boore and Brown 2000
), gastropods (Kurabayashi and Ueshima 2000), and echinoderms (Smith et al. 1993
), and placed brachiopods with the eutrochozoans (Stechmann and Schlegel 1999
).
The superphyla Deuterostomia (e.g., Chordata, Echinodermata, Hemichordata) and Protostomia (e.g., Arthropoda, Annelida, Mollusca) each appear to have evolved from an ancestor with a coelomic body cavity. Based largely on the view that the coelom evolved once, in their common ancestor, they are traditionally united in the clade Coelomata to the exclusion of the (presumably) earlier branching Platyhelminthes (Hyman 1951
). Platyhelminths share a number of developmental features with protostomes, such as spiral cleavage, the mesoderm originating from a mesentoblast, and the mouth arising from the blastopore. In the traditional phylogenetic scheme, these shared features are interpreted as primitive for all of the Coelomata, with the Deuterostome variations being derived later for that superphylum. However, others view these shared features as indicating that the Platyhelminthes evolved within the Protostomia, requiring the interpretation either that the coeloms of deuterostomes and protostomes are of separate origins or that Platyhelminthes have lost their coelom secondarily. Most who favor this latter phylogeny place Platyhelminthes as an early offshoot within Protostomia (see Brusca and Brusca 1990
), but some recent studies (Aguinaldo et al. 1997
; Balavoine 1997, 1998
; Carranza, Baguñà, and Riutort 1997
) have concluded that Platyhelminthes (or just those other than the Acoela [Ruiz-Trillo et al. 1999
]) are, instead, allied specifically with the group of coelomate protostomes dubbed the Eutrochozoa, e.g., Mollusca and Annelida, but not Arthropoda (Ghiselin 1988
).
Until recently, there were no complete mitochondrial sequences representing the phylum Platyhelminthes; however, these sequences are now available for the cestodes Taenia crassiceps (Le et al. 2000
) and Echinococcus multilocularis (GenBank record AB018440) and for the trematodes Fasciola hepatica (Garey and Wolstenholme 1989
; Le et al. 2000
; GenBank record AF216697), Schistosoma mansoni (Blair et al. 1999
; Le et al. 2000
), Schistosoma japonicum, Schistosoma mekongi, and Paragonimus westermani (Le et al. 2000
). However, these studies have emphasized specific comparative aspects without yet making a detailed genome description. We analyze and describe here the complete mitochondrial genome sequence of a platyhelminth, the tapeworm (i.e., class Cestoda) Hymenolepis diminuta, and compare aspects of this genome with those already published. Comparison of the mitochondrial gene arrangement of H. diminuta with all others available supports the view that platyhelminths are eutrochozoans.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Initially, a short fragment (450 nt) of H. diminuta cob was amplified by PCR using AmpliTaq DNA polymerase (Fisher) and 40 cycles of 94°C for 15 s, 53°C for 30 s, and 72°C for 2 min, followed by an incubation at 72° for 10 min. Amplifying primers (Cytb-424F [GGW TAY GTW YTW CCW TGR GGW CAR AT] and Cytb-876R [GCR TAW GCR AAW ARR AAR TAY CAY TCW GG]) were designed based on sequences well conserved in many distantly related taxa. The amplification product was visualized on a 1% agarose gel by ethidium bromide staining and UV irradiation and then purified by three serial passages through Ultrafree spin columns (Millipore; 30,000 NMWL) and used as a template for dRhodamine dye-terminator cycle sequencing reactions (Perkin-Elmer). Sequence was resolved using an ABI 310 automated DNA sequencer.
Six oligonucleotide primers were then designed for "long-PCR" (Barnes 1994
). The first two (Hym-Cytb-F [GCA CAA GAG TGG TAG TAA AAA TCC C] and Hym-Cytb-R [CGA ACA AGG GTG AAA CCC GTA ACA G]) matched the ends of the cob sequence as determined above. Two others (Hym-12S-F [TTT AGG ACT TGA TAG TAG GGT AGA C] and Hym-12S-R [ATC GTC CTT TAT AAC ACA CCT TCC C]) matched each end of a previously published fragment of H. diminuta rrnS (von Nickisch-Rosenegk, Lucius, and Loos-Frank 1999
). The remaining pair (Hym-16S-F [TTA TAA ATG GCC GCA GTA TAT TGA C] and Hym-16S-R [AGG CAA TTA ATT ATG CTA CCT TYG C]) was designed from an alignment of three cestode DNA fragments (450 bp of Echinococcus multilocularis, Catenotaenia lobata, and Schistocephalus solidus; unpublished data). In each case, these primers face "out" from the fragments. These primers were used in all 12 possible combinations with rTth-XL DNA polymerase (Perkin-Elmer). Reaction conditions were optimized for Mg++ concentration and annealing temperature as necessary. Visualization, purification, and sequencing of these products were as above, except that an ABI 377 automated sequencer was used for most reactions, with additional oligonucleotides being used for "primer walking." All sequences were determined on both strands.
Analysis
Sequences were produced and assembled using the ABI suite of programs (e.g., Sequencing Analysis and Sequence Navigator). Subsequent manipulations used MacVector, version 6.5 (Oxford Molecular Group), and GCG (Wisconsin Package, version 10.0, Genetics Computer Group, Madison, Wis.). Twelve protein-encoding genes and two ribosomal RNA genes were easily identified by comparisons with similar published sequences. Twenty-two tRNA genes were identified by eye based on their potential for forming tRNA-like secondary structures; specific identities were assigned by anticodon sequences.
Mitochondrial gene arrangements were compared by eye for gene adjacencies in all pairwise combinations for H. diminuta, E. multilocularis (GenBank record AB018440), F. hepatica (Garey and Wolstenholme 1989
; Le et al. 2000
; GenBank record AF216697), S. mansoni (Blair et al. 1999
; Le et al. 2000
), S. japonicum, S. mekongi, T. crassiceps (Le et al. 2000
), Lumbricus terrestris (Boore and Brown 1995
), Platynereis dumerilii (Boore and Brown 2000
; GenBank AF178678), Helobdella robusta (partial), Galathealinum brachiosum (partial) (Boore and Brown 2000
), Limulus polyphemus (Lavrov, Boore, and Brown 2000
), Tetilla sp. (partial) (Watkins and Beckenbach 1999
), and the gene arrangement previously inferred to be primitive for Vertebrata (Boore 1999
; Boore, Daehler, and Brown 1999
). Each of the three contending phylogenetic positions of the Platyhelminthes (see above and below) was evaluated by eye for whether these sets of shared gene boundaries were consistent or homoplastic.
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Gene Content and Organization
The 2 rRNA genes, 22 tRNA genes, and 12 of the 13 protein-encoding genes typical of metazoan mtDNAs were identified for H. diminuta by comparisons of their sequence similarity and/or potential secondary structures with those of other metazoans. Their arrangement is shown in figure 1
. One protein gene, atp8, typically found in animal mtDNAs, appears to be missing from this mtDNA (see below), as has been inferred for each of the other completely sequenced flatworm mtDNAs reported (Le et al. 2000
; GenBank record AB018440). Gene arrangements are substantially similar among the studied flatworm mtDNAs; details are presented below.
|
Absence of atp8
After careful analysis of the H. diminuta mtDNA sequence, the only open reading frame (ORF) that is a candidate for atp8 appears to be (mostly) within the 183-nt region otherwise described as a noncoding region (below). Designated ORF52, it is equivalent to 52 codons (plus a stop codon) in length, similar to the typical size of atp8, beginning 34 nt after trnY and overlapping trnS2(uga) by 10 nt (fig. 2
). This ORF starts with TTG, but the next codon is ATA, a more commonly used start codon for mitochondrial protein genes. It ends with TAA, although potential abbreviated stop codons (see below) are present that could truncate this ORF by one, two, or nine residues, thus requiring less or no overlap with the downstream trnS2(uga). However, there are several reasons to doubt that ORF52 is the atp8 gene. Blast searches and direct comparisons with published atp8 sequences and with the sequence of all published flatworm mtDNAs using both nucleotide and inferred amino acids reveal no similarity to this or to any other published mitochondrial gene. This lack of similarity also extends to the hydrophilicity profile, which is a well conserved feature of this protein (Watkins and Beckenbach 1999
). Characteristic amino acid signatures of the Atp8 protein (e.g., WXW near the carboxyl terminus) are not present in ORF52. Comparisons of the entire sequence of H. diminuta with published sequences of atp8 genes, as well as with the corresponding protein sequences, reveal no similarity. The gene for atp8 is therefore assumed to be absent, and it seems unlikely that ORF52 encodes any functional peptide, considering its small size, A+T-richness, and internal repeat (see below).
|
Initiation and Termination Codons
For H. diminuta mtDNA, 10 of the 12 protein genes can be inferred to initiate with ATG without overlap of the upstream gene. The exceptions are nad4, which appears to use ATT, and cox1, which appears to use GTT. In each case, there are no other eligible start codons that are reasonable alternatives, and the inferred amino terminal residues are similar to those of their homologs in other metazoans. The inference for nad4 requires overlap with the upstream nad4L, as is common for these two genes among metazoan mtDNAs. Although GTT is not usually designated as an eligible start codon, cox1 is most extreme among mitochondrial proteins with regard to initiating with unusual start codons (Wolstenholme 1992
).
Abbreviated stop codons (T and TA) are also common among metazoan mitochondrial protein genes, where posttranscriptional polyadenylation forms a complete UAA stop codon (Ojala, Montoya, and Attardi 1981
). For H. diminuta mtDNA, 11 of the 12 ORFs end with a complete stop codon (nine are TAG and two are TAA). Only cox1 is inferred to end with an abbreviated stop codon, although there is an in-frame TAG stop codon that, if used, would require overlap with the downstream tRNA gene by 10 nt.
Genetic Code
The mitochondrial genomes of most studied invertebrates appear to deviate from the "universal" genetic code with regard to the identities of ATA, TGA, and AGR codons (see Wolstenholme 1992
). Further differences have been inferred for Platyhelminthes such that AAA encodes N rather than K, ATA encodes I (as in the universal code) rather than M, and TAA encodes Y rather than termination (Bessho, Ohama, and Osawa 1992
).
The identities of AAA and of ATA in H. diminuta mtDNA appear to be consistent with these flatworm deviations. Pairwise comparisons using the two best conserved genes, cob and cox1, identified the codons corresponding to AAA of H. diminuta for L. terrestris (AAY = 6; AAR = 2), Katharina tunicata (AAY = 4; AAR = 0), and Drosophila yakuba (AAY = 4; AAR = 0) homologs. Similar comparisons for H. diminuta ATA codons determined their matches in L. terrestris (ATY = 9; ATR = 0), K. tunicata (ATY = 6; ATR = 0), and D. yakuba (ATY = 6; ATR = 0). There are weaknesses of such inferences: the pairwise comparisons are nonindependent due to the evolutionary relationships of the taxa, and it is not clear whether mutation bias or selection for amino acid composition can be stronger influences on substitution than the tendency to synonymous substitutions. However, acknowledging these limitations, we tentatively assign AAA and ATA as N and I, respectively, for H. diminuta mtDNA.
The third case that has been suggested as a deviation for flatworm mtDNAs, that in which TAA encodes Y rather than termination, does not seem possible for H. diminuta. No TAA codon is found in-frame for any protein-encoding gene. Two genes, cox2 and nad6, are predicted to end on a TAA stop codon, and one gene, cox1, is predicted to end on an abbreviated stop codon that is completed to UAA posttranscriptionally (fig. 2 ). Each of these genes would greatly overlap their adjacent downstream gene if they were to extend to the first in-frame TAG codon.
Nucleotide and Amino Acid Composition
Hymenolepis diminuta mtDNA is very A+T-rich (71%), even for an animal mitochondrial genome. The distribution of the purines versus the pyrimidines between the strands for each TA and GC pair is highly biased, with the reported strand being very rich in T and G. The values calculated for TA-skew ([T - A]/[T + A]) and GC-skew ([G - C]/[G + C]) (Perna and Kocher 1995
) are 0.29 and 0.33, respectively. (Skew values range from +1 to -1; the value is 0 if the strands have no skew.) The nucleotide composition of the reported strand is 25.4% A, 45.6% T, 19.3% G, and 9.6% C. Each of these values is even more extreme for the subset of nucleotides in the third positions of codons (table 1 ), which are presumably more free to change.
|
|
Noncoding Sequences
Despite its small size, H. diminuta mtDNA contains two noncoding regions of significant size. The larger, of 444 nt, is between nad5 and trnG and consists of 12 identical repeats of a 31-nt sequence (372 nt) plus one repeat with a single-nucleotide difference. Part of this 13th repeat can form a stem-loop with the remaining 41 nt of this noncoding region. This structure has 25 canonical base pairs and a loop of 11 nt (fig. 4
) and terminates at the beginning of trnG with an overlap of 2 nt. The single-nucleotide difference in this 13th unit is a T-to-C transition at nucleotide 16 of the repeat, creating a perfect match at this position of the stem. This is similar to reported repeat sequences of the nematode Meloidogyne javanica, which contains 11 identical repeats of a 63-nt unit, with the last copy modified by a single substitution, and which can also potentially form a stem-loop structure (Okimoto et al. 1991
). Large potential secondary structures can be found also between nad5 and trnG (as for H. diminuta) and between trnY and trnL1 for each of the cestodes E. multilocularis and T. crassiceps (fig. 4
). In some of these cases, there is potential for competing, mutually exclusive hairpins such as those identified in L. terrestris mtDNA (Boore and Brown 1995
). No similar structures could be identified for the noncoding regions of any of the trematode mtDNAs. The function of this region, if any, is unknown; however, similar stem-loop structures have been shown to initiate replication in mammals (e.g., Monnerot, Solignac, and Wolstenholme 1990
), and it is possible that the putative structure in this region serves a similar function in other metazoans.
|
Phylogenetic Reconstruction Using Gene Arrangements
The mitochondrial gene arrangements of two cestodes, T. crassiceps and E. multilocularis (Le et al. 2000
), are identical and differ from that of H. diminuta only by the reversal in order of trnS2 and trnL1. The trematodes S. japonicum, S. mekongi, and F. hepatica (Le et al. 2000
) have these two tRNA genes in the same order as in T. crassiceps and E. multilocularis and otherwise differ from H. diminuta only in the positions of trnE, trnV, and trnW for the schistosomes and the position of trnE for F. hepatica. The African schistosome S. mansoni has several additional gene rearrangements (Le et al. 2000
); these must be derived independently for this lineage based on the similarity among the cestodes and the other schistosomes, whose shared arrangements are parsimoniously inferred as primitive. Furthermore, the arrangement of the genes within the portion spanning from trnP to trnW of H. diminuta is also identical in three other cestodes whose mtDNAs have been only partially sequenced (Diphyllobothrium nihonkaiens [GenBank accession AB006205], Spirometra erinacei [GenBank accession AB006204], and S. proliferum [GenBank accession AB006203]) except for the lack of trnS1(gcu).
Hymenolepis diminuta has two pairs of adjacent genes that have been commonly found together among animal mtDNAs (see Boore 1999
), nad4L-nad4 and trnL1(uag)-trnL2(uaa). Although the significance is doubtful, it also shares single-gene boundaries with several nematodes: trnQ-trnF, as in Caenorhabditis elegans and Ascaris suum (Okimoto et al. 1992
), and cox1-trnT and trnV-trnA, as in M. javanica (Okimoto et al. 1991
). The most noteworthy gene arrangement similarities, however, are to several animals within the proposed higher-level group Eutrochozoa (e.g., mollusks and annelids; Ghiselin 1988
).
The arrangement of the genes in H. diminuta from nad1 through trnS1(gcu) is nearly identical to that found among several annelids, including one previously considered to be of the phylum Pogonophora (Galathealinum brachiosum) (Boore and Brown 1995, 2000
) (fig. 5 ) and, except for trnS1(gcu), in an echiuran (Urechis caupo; unpublished data). Although shared single-gene boundaries are less strong as a phylogenetic character (Boore and Brown 1998
), it is worth noting that all of these taxa also share the gene arrangement trnG-cox3.
|
However, none of the informative gene boundaries are found in S. mansoni mtDNA except for trnG-cox3. This is the basis for the cautionary message of Le et al. (2000). It is important to note that if the rearranged genome of S. mansoni had been the only one sampled from platyhelminths, there would have been a lack of phylogenetic resolution but still no homoplasy in this analysis. Although the generally slow rate of gene rearrangement in animal mtDNAs (Boore 1999
) makes it more likely that the signal of ancient relationships will be preserved, occasionally rapid rearrangements are not expected to be misleading in phylogenetic analysis. Even if there are many changes in some lineages, the great number of potential arrangements and the lack of apparent "hot spots" (Boore and Brown 1998
) make convergence unlikely, especially for multiple rearrangements (notable exceptions to this may be in apparent movements of noncoding sequence [Mindell, Sorenson, and Dimcheff 1998] or in exchange of nearest-neighbor tRNA genes [Boore and Brown 1998
; Boore 1999
]). Cladistic analysis as used here prevents artifactual clustering of taxa based on shared primitive characters (i.e., the "short" branches) but, rather, unites taxa only on the basis of shared, derived gene rearrangements parsimoniously inferred to have occurred only once in their common ancestor.
Arrangements of this set of genes in taxa related to those in figure 5
are consistent with this inference of evolutionary relationships. Echinoderms are thought to be the sister taxon to Chordates on the tree in figure 5C
; all echinoderm mtDNAs studied so far share the gene arrangements cox2-trnK-atp8 and trnL2(uaa)-nad1-trnI with chordates. There is slight similarity in the gene arrangements of H. diminuta with a few mollusks which would also be expected to be included in the clade of Eutrochozoa. In H. diminuta, the gene arrangement is nad1-trnN-trnP, whereas the arrangement is nad1-trnP in Dentalium eboreum, Nautilus sp. (unpublished data), Loligo bleekeri (GenBank AB009838), and Katharina tunicata (but with trnP in opposite orientation; Boore and Brown 1994
). Finally, it may be noteworthy that three of the four tRNA genes in the cluster preceding nad3 in H. diminuta are the same, although in different arrangements, as those preceding nad3 in K. tunicata (Boore and Brown 1994
) and Nautilus sp. (unpublished data).
The position of the Platyhelminthes in the system of Metazoa has been debated for decades (Brusca and Brusca 1990
). Traditionally, the coelom has been viewed as a feature uniting the protostomes and deuterostomes to the exclusion of Platyhelminthes (Hyman 1951
), although this notion is questioned by those who view the mechanisms of schizocoely (as for protostomes) and enterocoely (as for deuterostomes) as fundamentally different and therefore of separate origins. Some investigators view the simplicity of the Platyhelminthes as a secondary reduction from a coelomate ancestor (Ax 1963
; Smith and Tyler 1985
). Since platyhelminths share some of the basic characters of protostomes, such as spiral cleavage, the mesoderm originating (at least in part; see Boyer, Henry, and Martindale 1996
) from the 4d mesentoblast, and the mouth arising from the blastopore, some include them in Protostomia but generally as basal to the other, presumably derived, schizocoelomate phyla (see Brusca and Brusca 1990
).
However, an alternative interpretation is that Platyhelminthes are derived within the Eutrochozoa, a group of protostomes including annelids and mollusks, but not arthropods (Ghiselin 1988
). This view is supported by ultrastructural analysis (Rieger 1986
), by some sequence comparisons (Aguinaldo et al. 1997
; Balavoine 1997
; Carranza, Baguñà, and Riutort 1997
; Ruiz-Trillo et al. 1999
), and by our gene arrangement comparisons here. Further work on those flatworms comprising the Acoela may determine whether this group also occupies such a phylogenetic position or instead is basal among animal life, as has been inferred by DNA sequence comparisons (Ruiz-Trillo et al. 1999
).
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Keywords: Hymenolepis
tapeworm
phylogeny
mitochondria
evolution
genome
2 Address for correspondence and reprints: Markus von Nickisch-Rosenegk, Molecular Parasitology, Humboldt-University at Berlin, Philippstrasse 13, 10115 Berlin, Germany. laforge{at}rz.hu-berlin.de
![]() |
literature cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Aguinaldo, A. M., J. M. Turbeville, L. S. Linford, M. C. Rivera, J. R. Garey, R. A. Raff, and J. A. Lake. 1997. Evidence for a clade of nematodes, arthropods and other moulting animals. Nature 387:489493.
Anderson, S., A. T. Bankier, B. G. Barrell et al. (14 co-authors). 1981. Sequence and organization of the human mitochondrial genome. Nature 290:457465.
Ax, P. 1963. Relationships and phylogeny of the Turbellaria. Pp. 191224 in E. C. Dougherty, ed. The Lower Metazoa. University of California Press, Berkeley, Calif.
Balavoine, G. 1997. The early emergence of Platyhelminthes is contradicted by the agreement between 18S rRNA and Hox genes data. C. R. Acad. Sci. 320:8394.[ISI]
. 1998. Are Platyhelminthes coelomates without a coelom? An argument based on the evolution of Hox genes. Am. Zool. 38:843858.[ISI]
Barnes, W. M. 1994. PCR amplification of up to 35-kb DNA with high fidelity and high yield from bacteriophage templates. Proc. Natl. Acad. Sci. USA 91:22162220.
Bessho, Y., T. Ohama, and S. Osawa. 1992. Planarian mitochondria. II. The unique genetic code as deduced from cytochrome c oxidase subunit I gene sequences. J. Mol. Evol. 34:331335.[ISI][Medline]
Blair, D. T., H. Le, L. Despres, and D. P. McManus. 1999. Mitochondrial genes of Schistosoma mansoni. Parasitology 119:303313.
Boore, J. L. 1999. Animal mitochondrial genomes. Nucleic Acids Res 27:17671780.
Boore, J. L., and W. M. Brown. 1994. Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata. Genetics 138:423443.
. 1995. Complete DNA sequence of the mitochondrial genome of the annelid worm, Lumbricus terrestris. Genetics 141:305319.
. 1998. Big trees from little genomes: Mitochondrial gene order as a phylogenetic tool. Curr. Opin. Genet. Dev. 8:668674.[ISI][Medline]
. 2000. Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa. Mol. Biol. Evol. 17:87106.
Boore, J. L., T. M. Collins, D. Stanton, L. L. Daehler, and W. M. Brown. 1995. Deducing arthropod phylogeny from mitochondrial DNA rearrangements. Nature 376:163165.
Boore, J. L., L. L. Daehler, and W. M. Brown. 1999. Complete sequence, gene arrangement and genetic code of mitochondrial DNA of the cephalochordate Branchiostoma floridae ("amphioxus"). Mol. Biol. Evol. 16:410418.[Abstract]
Boore, J. L., D. Lavrov, and W. M. Brown. 1998. Gene translocation links insects and crustaceans. Nature 393:667668.
Boyer, B. C., J. Q. Henry, and M. Q. Martindale. 1996. Dual origins of mesoderm in a basal member of the spiralian clade: cell lineage studies in the polyclad turbellarian Hoploplana inquilina. Dev. Biol. 179:329338.
Brusca, R. C., and G. J. Brusca. 1990. Invertebrates. Sinauer, Sunderland, Mass.
Carranza, S., J. Baguñà, and M. Riutort. 1997. Is the Platyhelminthes a monophyletic primitive group? An assessment using 18S rDNA sequences. Mol. Biol. Evol. 14:485497.[Abstract]
Dowton, M. 1999. Relationships among the cyclostome braconid (Hymenoptera: Braconidae) subfamilies inferred from a mitochondrial tRNA gene rearrangement. Mol. Phylogenet. Evol. 11:283287.[ISI][Medline]
Garey, J. R., and D. R. Wolstenholme. 1989. Platyhelminth mitochondrial DNA: evidence for early evolutionary origin of a tRNAseragn that contains a dihydrouridine arm replacement loop and of serine-specifying AGA and AGG codons. J. Mol. Evol. 28:374387.[ISI][Medline]
Ghiselin, M. T. 1988. The origin of molluscs in the light of molecular evidence. Oxf. Surv. Evol. Biol. 5:6695.
Hoffmann, R. J., J. L. Boore, and W. M. Brown. 1992. A novel mitochondrial genome organization for the blue mussel, Mytilus edulis. Genetics 131:397412.
Hyman, L. H. 1951. The invertebrates, Vol. 2. Platyhelminthes and Rhynchocoela: the acoelomate bilateria. McGraw-Hill, New York.
Keddie, E. M., T. Higazi, and T. R. Unnasch. 1998. The mitochondrial genome of Onchocerca volvulus: sequence, structure and phylogenetic analysis. Mol. Biochem. Parasitol. 95:111127.[ISI][Medline]
Kurabayashi, A., and R. Ueshima. 2000. Complete sequence of the mitochondrial DNA of the primitive opisthobranch Pupa strigosa: systematic implications of the genome organization. Mol. Biol. Evol. 17:266277.
Lavrov, D., J. L. Boore, and W. M. Brown. 2000. The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus. Mol. Biol. Evol. 17:813824.
Le, T. H., D. Blair, T. Agatsuma et al. (14 co-authors). 2000. Phylogenies inferred from mitochondrial gene ordersa cautionary tale from the parasitic flatworms. Mol. Biol. Evol. 17:11231125.
Mindell, D., M. D. Sorenson, and D. E. Dimcheff. 1998. Multiple independent origins of mitochondrial gene order in birds. Proc. Natl. Acad. Sci. USA 95:1069310697.
Monnerot, M., M. Solignac, and D. R. Wolstenholme. 1990. Discrepancy in divergence of the mitochondrial and nuclear genomes of Drosophila teissieri and Drosophila yakuba. J. Mol. Evol. 30:500508.[ISI][Medline]
Ojala, D., J. Montoya, and G. Attardi. 1981. tRNA punctuation model of RNA processing in human mitochondria. Nature 290:470474.
Okimoto, R., H. M. Chamberlin, J. L. Macfarlane, and D. R. Wolstenholme. 1991. Repeated sequence sets in mitochondrial DNA molecules of root knot nematodes (Meloidogyne): nucleotide sequences, genome location and potential for host race identification. Nucleic Acids Res. 19:16191626.[Abstract]
Okimoto, R., J. L. Macfarlane, D. O. Clary, and D. R. Wolstenholme. 1992. The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics 130:471498.
Perna, N. T., and T. D. Kocher. 1995. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 41:353358.[ISI][Medline]
Rieger, R. M. 1986. Über den Ursprung der Bilateria: die Bedeutung der Ultrastrukturforschung für ein neues Verstehen der Metazoenevolution. Verh. Dtsch. Zool. Ges. 79:3150.
Ruiz-Trillo, I., M. Riutort, D. T. J. Littlewood, E. A. Herniou, and J. Baguñà. 1999. Acoel flatworms: earliest extant bilaterian metazoans, not members of Platyhelminthes. Science 283:19191923.
Shadel, G. S., and D. A. Clayton. 1997. Mitochondrial DNA maintenance in vertebrates. Annu. Rev. Biochem. 66:409435.[ISI][Medline]
Smith, J. III, and S. Tyler. 1985. The acoel turbellarians: kingpins of metazoan evolution or a specialized offshoot? Pp. 123142 in S. Conway Morris, J. D. George, R. Gibson, and H. M. Platt, eds. The origins and relationships of lower invertebrates. Systematics Association Special Vol. 28, Oxford University Press, Oxford, England.
Smith, M. J., A. Arndt, S. Gorski, and E. Fajber. 1993. The phylogeny of echinoderm classes based on mitochondrial gene arrangements. J. Mol. Evol. 36:545554.[ISI][Medline]
Stechmann, A., and M. Schlegel. 1999. Analysis of the complete mitochondrial DNA sequence of the brachiopod Terebratulina retusa places Brachiopoda within the protostomes. Proc. R. Soc. Lond. B Biol. Sci. 260:20432052.
von Nickisch-Rosenegk, M., R. Lucius, and B. Loos-Frank. 1999. Contributions to the phylogeny of the Cyclophyllidea (Cestoda) inferred from mitochondrial 12S rDNA. J. Mol. Evol. 48:586596.[ISI][Medline]
Watkins, R. F., and A. T. Beckenbach. 1999. Partial sequence of a sponge mitochondrial genome reveals sequence similarity to Cnidaria in cytochrome oxidase subunit II and the large ribosomal RNA subunit. J. Mol. Evol. 48:542554.[ISI][Medline]
Wolstenholme, D. R. 1992. Genetic novelties in mitochondrial genomes of multicellular animals. Curr. Opin. Genet. Dev. 2:918925.[Medline]