Department of Biology, University of Michigan; and DOE Joint Genome Institute and Lawrence Berkeley National Laboratory, Walnut Creek, California
Correspondence: E-mail: jlboore{at}lbl.gov.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: Phoronis architecta lophophorates mitochondrial genome metazoan phylogeny gene arrangements evolution
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The number of complete mitochondrial DNA (mtDNA) sequences is growing rapidly. Among the triploblast animals, this circular molecule ranges in size from about 13 to 18 kb and almost always contains the same 37 genes: 13 for proteins, 22 for tRNAs, and 2 for rRNAs (Boore 1999). For some groups of animals, the arrangement of these genes is very stable over long periods of time; for example, all 37 genes of human and shark mtDNA are identically arranged (Boore 1999). For other lineages, rearrangements are much more rapid; for example, other than the short block trnL1-trnL2-nad1, there are no gene boundaries in common between the mitochondrial genomes of two mollusks, a chiton (class Polyplacophora) Katharina tunicata (Boore and Brown 1994) and a mussel (class Bivalvia) Mytilus edulis (Hoffman, Boore, and Brown 1992).
In recent years, the comparison of complete mitochondrial genomes has become a very powerful tool for reconstructing phylogeny. The comparisons of complete mitochondrial genome sequences have been shown to have much greater resolving power than for individual gene sequences (Cao et al. 2000; Ingman et al. 2001; Miya, Kawaguchi, and Nishida 2001). In addition, the relative arrangement of these 37 genes can be an especially powerful tool for phylogenetic analysis (Boore et al. 1995; Boore and Brown 1998, 2000). The great number of potential gene rearrangements makes it very unlikely that different animal lineages would independently adopt the same gene order or that any gene would move back to a previous location, so gene arrangements that are shared in a derived condition (see below) are a very strong indicator of relatedness (Boore and Brown 1998).
Among lophophorate taxa, the only complete mitochondrial genome sequences available are for three brachiopods, Terebratulina retusa (Stechmann and Schlegel 1999), Laqueus rubellus (Noguchi et al., 2000), and Terebratalia transversa (Helfenbein, Brown, and Boore 2001). Although similarity of gene arrangement to non-lophophorate taxa was noted in these studies, no formal phylogenetic analysis was performed. Phoronida is a second lophophorate phylum, comprising coelomate, vermiform animals living in chitinous tubes. Here we report the sequence of nearly all (i.e., 14,018 nts) of the mitochondrial genome of the phoronid Phoronis architecta, from which we can determine the arrangements of all 13 of the expected protein-encoding genes, both of the rRNA genes, and 16 of the expected 22 tRNA genes.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Sequence Analysis
Sequences were produced and assembled using the ABI suite of programs (e.g., Sequencing Analysis, Sequence Navigator, Autoassembler). Subsequent manipulations used MacVector 6.5 (Oxford Molecular Group). Protein encoding genes were identified by BLAST matching to other animal mtDNAs, with start codons inferred as being an eligible in-frame start codon corresponding at least to the extent of alignment that does not overlap the upstream gene. Abbreviated stop codons, presumably completed by polyadenylation of the mRNA after cleavage, were inferred in cases where extensive gene overlap beyond the extent of homologous gene-sequence matching would be required to reach the first complete stop codon. Ribosomal RNAs were identified by BLAST matching. Transfer RNA genes were identified generically by their potential secondary structures and specifically by anticodon sequence, in some cases using the tRNAscan-SE Search Server (www.genetics.wustl.edu/eddy/tRNAscan-SE).
Sequences and Alignments
We chose a broad representation of taxa from available complete mtDNA sequences for this study. All of the mtDNAs of the triploblast taxa contain unambiguous copies of the 37 genes mentioned above. The following sequences were retrieved from GenBank: Terebratulina retusa (accession number NC_000941), Laqueus rubellus (NC_002322), Terebratalia transversa (NC_003086), Platynereis dumerilii (NC_000931), Lumbricus terrestris (NC_001673), Loligo bleekeri (NC_002507), Katharina tunicata (NC_001636), Drosophila yakuba (NC_001322), Daphnia pulex (NC_000844), Locusta migratoria (NC_001712), Limulus polyphemus (NC_003057), Paracentrotus lividus (NC_001572), Asterina pectinifera (NC_001627), Florometra serratissima (NC_001878), Balanoglossus carnosus (NC_001887), Homo sapiens (NC_001807), Mustelus manazo (NC_000890), Metridium senile (NC_000933), Podospora anserina (NC_001329), and Pichia canadensis (NC_001762).
Inferred amino acid sequences of nine mitochondrial protein genes (atp6, cob, cox1, cox2, cox3, nad1, nad3, nad4, and nad5) from 18 ingroup and three outgroup taxa were aligned using ClustalW. Pairwise alignments were done in slow mode with an open gap penalty of 10 and an extend gap penalty of 1 using the BLOSUM similarity matrix. The multiple alignment used the same parameters as the pairwise plus a "delay divergent sequences" setting of 40%. Ambiguously aligned positions at the termini of each gene alignment were trimmed. Gaps were included in the analysis and scored as "missing characters." (See Supplementary Material online.)
Phylogenetic Reconstruction
Parsimony analysis of aligned inferred amino acid sequences employed a heuristic search with 10,000 random additions in PAUP* version 4.01b10 (Swofford 1998). Bootstrap values were calculated using 500 replicates with an addition sequence of 10 random replicates. Bremer branch support values (i.e., decay indices; Bremer 1994) were calculated using Tree Rot (Sorenson 1996). Maximum-likelihood analysis of the same alignment employed 10,000 quartet puzzling steps, an mtREV24 model of substitution, and 8 Gamma rate categories in TreePuzzle version 5.0 (Strimmer and von Haeseler 1996).
For phylogenetic analysis of gene arrangements, a gene adjacency matrix was constructed for 18 animal mtDNAs as in Boore et al. (1995) in their analysis of arthropod gene arrangements. Briefly, this scores 74 characters as being "upstream of" and "downstream of" (according to transcriptional orientation) each of the 37 genes. For each character there are 72 possible character states, each being the 5' or 3' end of one of the other 36 genes. This matrix is then analyzed using parsimony criteria with PAUP* version 4.01b10 (Swofford 1998). This is essentially identical to the technique termed "maximum parsimony of multiple encodings" (MPME) by Wang et al. (2002). Boundaries not identified for the Phoronis mtDNA were scored as missing data.
The bootstrap consensus tree of the gene adjacency matrix was calculated using 1,000 replicates with an addition sequence of 10 random replicates. No outgroups were used in this analysis because no sequenced mtDNA from an appropriate taxon has the same gene content as the triploblast animals studied here. All of the gene arrangements used in this study can be found listed in the Mitochondrial Gene Arrangement Source Guide (follow the Evolutionary Genomics link from http://www.jgi.doe.gov/).
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Remarkably, the arrangement of the 31 genes in the P. architecta mtDNA is nearly identical to that of the mollusk Katharina tunicata (Boore and Brown 1994). Only three of these 31 genes are arranged differently between these two mitochondrial genomes (fig. 1). The six unaccounted for tRNA genes of P. architecta are clustered in the K. tunicata mtDNA and are between rrnS and cox3, inviting the speculation that they may also have similar or identical arrangement in the phoronid mtDNA.
|
Our cladistic analysis of the gene adjacency matrix for 18 animals resulted in 108 trees of equal length. The strict consensus of the 108 trees includes the lophophorates among the protostomes in a lophotrochozoan clade with 98% bootstrap support (see fig. 2). This data set includes 73 parsimony informative characters, i.e., gene junctions that are shared by members of one clade, but for which other taxa share an identically rearranged state. For example, preceding trnK of Katharina and Phoronis is cox3, whereas preceding trnK of Limulus, many vertebrates, acorn worm, and echinoderms is cox2. For another example, Katharina, Phoronis, and Limulus all share the gene arrangement trnR-trnN, whereas vertebrates, acorn worm, and echinoderms share the arrangement trnR-nad4L.
|
|
Although there is strong evidence in the mtDNA data presented here for a protostome clade containing phoronids, brachiopods, annelids, and mollusks, there is not enough phylogenetic information in the data set to resolve the relationships within the group. Of particular interest is the question of lophophorate monophyly. Halanych and co-workers found evidence of a clade of phoronids and articulate brachiopods (Halanych et al. 1995), but this result was severely questioned (Conway Morris et al. 1996). Some evidence for monophyly of brachiopods + phoronids has been presented (e.g., Cohen 2000), but more analyses of nuclear genes are needed to resolve this question and other questions of relationships within the Lophotrochozoa.
From the results of these analyses, it appears that the morphological and embryological features thought to indicate a close lophophorate-deuterostome relationship such as the mouth deriving from a secondary body opening rather than from the blastopore, the presence of ciliary bands formed by monociliate cells, the position of the nerve concentration/brain, the presence of a tripartite coelom, and other features (Nielsen, Scharff, and Eibye-Jacobsen 1996; Brusca and Brusca 2003) have either arisen independently or have been lost multiple times within the triploblastic animals. Indeed, in so strongly corroborating the phylogenetic position of phoronids and brachiopods as protostomes with this study, we would agree that morphological traits (e.g., egg cleavage patterns, modes of coelom formation, mouth formation, and trimery) are either unreliable characters for deep metazoan phylogeny, too open to interpretation as to their character states, or not polarizable, as has been suggested previously (De Rosza 2001). That the morphological characters commonly used for metazoan systematics have not evolved as conservatively as had been believed is also suggested by a recent study of the Myzostomida (Eeckhaut et al. 2000).
Resolution of the relationships within the Lophotrochozoa is an important and seemingly difficult problem in the study of animal evolution. It may be possible to resolve these questions using multiple nuclear-encoded genes. To this end, and to gather data for many additional analyses, we recommend that large-scale genome sequencing be performed on members of this clade that until now has been overlooked by the genome-sequencing community.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Hervé Philippe, Associate Editor
![]() |
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:489-493.[CrossRef][ISI][Medline]
Boore, J. L. 1999. Animal mitochondrial genomes. Nucleic Acids Res. 27:1767-1780.
Boore, J. L., and W. M. Brown. 1994. Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata. Genetics 138:423-443.
Boore, J. L., and W. M. Brown. 1998. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr. Opin. Genet. Dev. 8:668-674.[CrossRef][ISI][Medline]
Boore, J. L., and W. M. Brown. 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:87-106.
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:163-165.[CrossRef][ISI][Medline]
Bremer, K. 1994. Branch support and tree stability. Cladistics 10:295-304.[CrossRef][ISI]
Brusca, R. C., and G. J. Brusca. 2003. Invertebrates, 2nd edition. Sinauer Associates, Sunderland, Mass.
Cao, Y., M. Fujiwara, M. Nikaido, N. Okada, and M. Hasegawa. 2000. Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data. Gene 259:149-158.[CrossRef][ISI][Medline]
Cohen, B. L. 2000. Monophyly of brachiopods and phoronids: reconciliation of molecular evidence with Linnaean classification (the subphylum Phoroniformea nov.). Proc. R. Soc. Lond. Ser. B Biol. Sci. 267:225-231.[CrossRef][ISI][Medline]
Cohen, B. L., A. Gawthrop, and T. Cavalier-Smith. 1998. Molecular phylogeny of brachiopods and phoronids based on nuclear-encoded small subunit ribosomal RNA gene sequences. Phil. Trans. R. Soc. Lond. Ser. B. 353:2039-2061.[CrossRef][ISI]
Cohen, B. L., S. Stark, A. B. Gawthrop, M. E. Burke, and C. W. Thayer. 1998. Comparison of articulate brachiopod nuclear and mitochondrial gene trees leads to a clade-based redefinition of protostomes Protostomozoa and deuterostomes Deuterostomozoa. Proc. R. Soc. Lond. Ser. B Biol. Sci. 265:475-482.[CrossRef][ISI][Medline]
Conway Morris, C., B. L. Cohen, A. B. Gawthrop, T. Cavalier-Smith, and B. Winnepenninckx. 1996. Lophophorate phylogeny. Science 272:282.[ISI][Medline]
De Rosza, R. 2001. Molecular data indicate the protostome affinity of brachiopods. Syst. Biol. 50:848-859.[CrossRef][ISI][Medline]
Eeckhaut, I., D. McHugh, P. Mardulyn, R. Tiedemann, D. Monteyne, M. Jangoux, and M. C. Milinkovitch. 2000. Myzostomida: a link between trochozoans and flatworms? Proc. R. Soc. Lond. Ser. B Biol. Sci. 267:1383-1392.[CrossRef][ISI][Medline]
Endo, K. 2001. The phylogenetic position of brachiopods inferred from mitochondrial gene orders. Pp. 142149 in C. H. C. Brunton, L. R. M. Cocks, and S. L. Long, eds. Brachiopods past and present. Taylor & Francis, New York, N.Y.
Field, K. G., G. J. Olsen, D. J. Lane, S. J. Giovannoni, M. T. Ghiselin, E. C. Raff, M. R. Pace, and R. A. Raff. 1988. Molecular phylogeny of the animal kingdom. Science 239:748-753.[ISI][Medline]
Halanych, K. M., J. D. Bacheller, A. M. Aguinaldo, S. M. Liva, D. M. Hillis, and J. A. Lake. 1995. Evidence from 18S ribosomal DNA that the lophophorates are protostome animals. Science 267:1641-1643.[ISI][Medline]
Helfenbein, K. G., W. M. Brown, and J. L. Boore. 2001. The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa. Mol. Biol. Evol. 18:1734-1744.
Hoffmann, R. J., J. L. Boore, and W. M. Brown. 1992. A novel mitochondrial genome organization for the blue mussel, Mytilus edulis. Genetics 131:397-412.
Hyman, L. 1959. The Invertebrates, Vol. 5,. Smaller Coelomate Groups, McGraw-Hill, New York.
Ingman, M., H. Kaessmann, S. Pääbo, and U. Gyllensten. 2001. Mitochondrial genome variation and the origin of modern humans. Nature 408:708-713.[CrossRef][ISI]
Mackey, L. Y., B. Winnepenninckx, R. DeWachter, T. Backeljau, P. Emschermann, and J. R. Garey. 1996. 18S rRNA suggests that Entoprocta are protostomes, unrelated to Ectoprocta. J. Mol. Evol. 42:552-559.[ISI][Medline]
Maddison, W. P., and D. R. Maddison. 2000. MacClade4: Analysis of Phylogeny and Character Evolution. Sinauer Associates, Sunderland, Mass.
Miya, M., A. Kawaguchi, and M. Nishida. 2001. Mitogenomic exploration of higher teleostean phylogenies: a case study for moderate-scale evolutionary genomics with 38 newly determined complete mitochondrial DNA sequences. Mol. Biol. Evol. 18:1993-2009.
Nielsen, C., N. Scharff, and D. Eibye-Jacobsen. 1996. Cladistic analyses of the animal kingdom. Biol. J. Linn. Soc. Lond. 57:385-410.[CrossRef][ISI]
Noguchi, Y., K. Endo, F. Tajima, and R. Ueshima. 2000. The mitochondrial genome of the brachiopod Laqueus rubellus. Genetics 155:245-259.
Sorenson, M. D. 1996. TreeRot. University of Michigan, Ann Arbor.
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. Ser. B Biol. Sci. 260:2043-2052.
Strimmer, K., and A. von Haeseler. 1996. Quartet Puzzling: a quartet maximum-likelihood method for reconstructing tree topologies. Mol. Biol. Evol. 13:964-969.
Swofford, D. L. 1998. PAUP*, Phylogenetic Analysis Using Parsimony (*and other methods). Version 4. Sinauer Associates, Sunderland, Mass.
Wang, L-S., R. K. Jansen, B. M. E. Moret, L. A. Raubeson, and T. Warnow. 2002. Phylogenetic methods for the analysis of genome rearrangement data: an empirical study. Pp. 524535 in Proceedings of Fifth Pacific Symp. of Biocomputing PSB '02, Hawaii, USA.