Analysis of the complete mitochondrial DNA sequence of the brachiopod terebratulina retusa places Brachiopoda within the protostomes

Brachiopod phylogeny is still a controversial subject. Analyses using nuclear 18SrRNA and mitochondrial 12SrDNA sequences place them within the protostomes but some recent interpretations of morphological data support a relationship with deuterostomes. In order to investigate brachiopod affinities within the metazoa further, we compared the gene arrangement on the brachiopod mitochondrial genome with several metazoan taxa. The complete (15 451 bp) mitochondrial DNA (mtDNA) sequence of the articulate brachiopod Terebratulina retusa was determined from two overlapping long polymerase chain reaction products. All the genes are encoded on the same strand and gene order comparisons showed that.only one major rearrangement is required to interconvert the T. retusa and Katharina tunicata (Mollusca: Polvplacophora) mitochondrial genomes. The partial mtDNA sequence of the prosobranch mollusc Littorina saxatilis shows complete congruence with the T. rehtusa gene arrangement with regard to the ribosomal and protein coding genes. This high similarity in gene arrangement is the first to be reported within the protostomes. Sequence analyses of mitochondrial protein coding genes also support a close relationship of the brachiopod with molluscs and annelids, thus supporting the clade Lophotrochozoa. Though being highly informative, sequence analyses of the mitochondrial protein coding genes failed to resolve the branching order within the lophotrochozoa.     

The mitochondrial genome of the brachiopod Laqueus rubellus

The complete nucleotide sequence of the 14,017-bp mitochondrial (mt) genome of the articulate brachiopod Laqueus rubellus is presented. Being one of the smallest of known mt genomes, it has an extremely compact gene organization. While the same 13 polypeptides, two rRNAs, and 22 tRNAs are encoded as in most other animal mtDNAs, lengthy noncoding regions are absent, with the longest apparent intergenic sequence being 54 bp in length. Gene-end sequence overlaps are prevalent, and several stop codons are abbreviated. The genes are generally shorter, and three of the protein-coding genes are the shortest among known homologues. All of the tRNA genes indicate size reduction in either or both of the putative TPsiC and DHU arms compared with standard tRNAs. Possession of a TV (TPsiC arm-variable loop) replacement loop is inferred for tRNA(R) and tRNA(L-tag). The DHU arm appears to be unpaired not only in tRNA(S-tct) and tRNA(S-tga), but also in tRNA(C), tRNA(I), and tRNA(T), a novel condition. All the genes are encoded in the same DNA strand, which has a base composition rich in thymine and guanine. The genome has an overall gene arrangement drastically different from that of any other organisms so far reported, but contains several short segments, composed of 2-3 genes, which are found in other mt genomes. Combined cooccurrence of such gene assortments indicates that the Laqueus mt genome is similar to the annelid Lumbricus, the mollusc Katharina, and the octocoral Sarcophyton mt genomes, each with statistical significance. Widely accepted schemes of metazoan phylogeny suggest that the similarity with the octocoral could have arisen through a process of convergent evolution, while it appears likely that the similarities with the annelid and the mollusc reflect phylogenetic relationships.     

Complete sequences of the highly rearranged molluscan mitochondrial genomes of the Scaphopod Graptacme eborea and the bivalve Mytilus edulis

We have determined the complete sequence of the mitochondrial genome of the scaphopod mollusk Graptacme eborea (14,492 nts) and completed the sequence of the mitochondrial genome of the bivalve mollusk Mytilus edulis (16,740 nts). (The name Graptacme eborea is a revision of the species formerly known as Dentalium eboreum.) G. eborea mtDNA contains the 37 genes that are typically found and has the genes divided about evenly between the two strands, but M. edulis contains an extra trnM and is missing atp8, and it has all genes on the same strand. Each has a highly rearranged gene order relative to each other and to all other studied mtDNAs. G. eborea mtDNA has almost no strand skew, but the coding strand of M. edulis mtDNA is very rich in G and T. This is reflected in differential codon usage patterns and even in amino acid compositions. G. eborea mtDNA has fewer noncoding nucleotides than any other mtDNA studied to date, with the largest noncoding region only 24 nt long. Phylogenetic analysis using 2,420 aligned amino acid positions of concatenated proteins weakly supports an association of the scaphopod with gastropods to the exclusion of Bivalvia, Cephalopoda, and Polyplacophora, but it is generally unable to convincingly resolve the relationships among major groups of the Lophotrochozoa, in contrast to the good resolution seen for several other major metazoan groups.     

The Mitochondrial Genomes of Two Nematodes, Caenorhabditis Elegans and Ascaris Suum

The nucleotide sequences of the mitochondrial DNA (mtDNA) molecules of two nematodes, Caenorhabditis elegans [13,794 nucleotide pairs (ntp)], and Ascaris suum (14,284 ntp) are presented and compared. Each molecule contains the genes for two ribosomal RNAs (s-rRNA and l-rRNA), 22 transfer RNAs (tRNAs) and 12 proteins, all of which are transcribed in the same direction. The protein genes are the same as 12 of the 13 protein genes found in other metazoan mtDNAs: Cyt b, cytochrome b; COI-III, cytochrome c oxidase subunits I-III; ATPase6, Fo ATPase subunit 6; ND1-6 and 4L, NADH dehydrogenase subunits 1-6 and 4L: a gene for ATPase subunit 8, common to other metazoan mtDNAs, has not been identified in nematode mtDNAs. The C. elegans and A. suum mtDNA molecules both include an apparently noncoding sequence that contains runs of AT dinucleotides, and direct and inverted repeats (the AT region: 466 and 886 ntp, respectively). A second, apparently noncoding sequence in the C. elegans and A. suum mtDNA molecules (109 and 117 ntp, respectively) includes a single, hairpin-forming structure. There are only 38 and 89 other intergenic nucleotides in the C. elegans and A. suum mtDNAs, and no introns. Gene arrangements are identical in the C. elegans and A. suum mtDNA molecules except that the AT regions have different relative locations. However, the arrangement of genes in the two nematode mtDNAs differs extensively from gene arrangements in all other sequenced metazoan mtDNAs. Unusual features regarding nematode mitochondrial tRNA genes and mitochondrial protein gene initiation codons, previously described by us, are reviewed. In the C. elegans and A. suum mt-genetic codes, AGA and AGG specify serine, TGA specifies tryptophan and ATA specifies methionine. From considerations of amino acid and nucleotide sequence similarities it appears likely that the C. elegans and A. suum ancestral lines diverged close to the time of divergence of the cow and human ancestral lines, about 80 million years ago.     

Complete DNA sequence of the mitochondrial genome of the ascidian Halocynthia roretzi (Chordata, Urochordata)

The complete nucleotide sequence of the 14,771-bp-long mitochondrial (mt) DNA of a urochordate (Chordata)-the ascidian Halocynthia roretzi-was determined. All the Halocynthia mt-genes were found to be located on a single strand, which is rich in T and G rather than in A and C. Like nematode and Mytilus edulis mtDNAs, that of Halocynthia encodes no ATP synthetase subunit 8 gene. However, it does encode an additional tRNA gene for glycine (anticodon TCT) that enables Halocynthia mitochondria to use AGA and AGG codons for glycine. The mtDNA carries an unusual tRNA(Met) gene with a TAT anticodon instead of the usual tRNA(Met)(CAT) gene. As in other metazoan mtDNAs, there is not any long noncoding region. The gene order of Halocynthia mtDNA is completely different from that of vertebrate mtDNAs except for tRNA(His)-tRNA(Ser)(GCU), suggesting that evolutionary change in the mt-gene structure is much accelerated in the urochordate line compared with that in vertebrates. The amino acid sequences of Halocynthia mt-proteins deduced from their gene sequences are quite different from those in other metazoans, indicating that the substitution rate in Halocynthia mt-protein genes is also accelerated.     

Complete Sequence of the Mitochondrial DNA of the Annelid Worm Lumbricus Terrestris

We have determined the complete nucleotide (nt) sequence of the mitochondrial genome of an oligochaete annelid, the earthworm Lumbricus terrestris. This genome contains the 37 genes typical of metazoan mitochondrial DNA (mtDNA), including ATPase8, which is missing from some invertebrate mtDNAs. ATPase8 is not immediately upstream of ATPase6, a condition found previously only in the mtDNA of snails. All genes are transcribed from the same DNA strand. The largest noncoding region is 384 nt and is characterized by several homopolymer runs, a tract of alternating TA pairs, and potential secondary structures. All protein-encoding genes either overlap the adjacent downstream gene or end at an abbreviated stop codon. In Lumbricus mitochondria, the variation of the genetic code that is typical of most invertebrate mitochondrial genomes is used. Only the codon ATG is used for translation initiation. Lumbricus mtDNA is A + T rich, which appears to affect the codon usage pattern. The DHU arm appears to be unpaired not only in tRNA(ser(AGN)), as is typical for metazoans, but perhaps also in tRNA(ser(UCN)), a condition found previously only in a chiton and among nematodes. Relating the Lumbricus gene organization to those of other major protostome groups requires numerous rearrangements.     

A Novel Mitochondrial Genome Organization for the Blue Mussel, Mytilus Edulis

The sequence of 13.9 kilobases (kb) of the 17.1-kb mitochondrial genome of Mytilus edulis has been determined, and the arrangement of all genes has been deduced. Mytilus mitochondrial DNA (mtDNA) contains 37 genes, all of which are transcribed from the same DNA strand. The gene content of Mytilus is typically metazoan in that it includes genes for large and small ribosomal RNAs, for a complete set of transfer RNAs and for 12 proteins. The protein genes encode the cytochrome b apoenzyme, cytochrome c oxidase (CO) subunits I-III, NADH dehydrogenase (ND) subunits 1-6 and 4L, and ATP synthetase (ATPase) subunit 6. No gene for ATPase subunit 8 could be found. The reading frames for the ND1, COI, and COIII genes contain long extensions relative to those genes in other metazoan mtDNAs. There are 23 tRNA genes, one more than previously found in any metazoan mtDNA. The additional tRNA appears to specify methionine, making Mytilus mtDNA unique in having two tRNA(Met) genes. Five lengthy unassigned intergenic sequences are present, four of which vary in length from 79 to 119 nucleotides and the largest of which is 1.2 kb. The base compositions of these are unremarkable and do not differ significantly from that of the remainder of the mtDNA. The arrangement of genes in Mytilus mtDNA is remarkably unlike that found in any other known metazoan mtDNA.     

Nearly complete mitochondrial genome of Polyascus gregaria and the phylogenetic relationships among maxillopodans

We determined for the first time the nearly complete mitochondrial genome sequence of the entozoic Polyascus gregaria, a representative of Rhizocephala, Cirripedia. The nearly complete mitogenome was 15, 465 bp in length, consisting of 11 protein-coding genes, two rRNA genes, 22 tRNA genes and one major incomplete noncoding region. In total there are 73 overlapping nucleotides and 17 spacers between genes. All genes sequenced in P. gregaria mtDNA (including RNAs) were encoded on the same strand of the DNA, and the gene arrangement differed from that of other metazoan animals. The mitochondrial genome rearrangements included translocation of at least 8 genes and even inversion of the coding polarity of at least 2 genes. Comparative analysis of the gene orders with other maxillopodan mtDNAs showed that the unique characteristics of the thoracican cirripeds lineage were not observed in this representative of rhizocephalan. Phylogenetic analyses supported a close affinity of Rhizocephala to Thoracica. By adding the mitochondrial genomes from 4 copepods, the reciprocally monophyletic cirripeds and copepods clustered as sister groups, refusing the close relationship between Cirripedia and Remipedia. However, the monophyly of Maxillopoda was not supported in this study.     

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

Summary The nucleotide sequence of a segment of the mitochondrial DNA (mtDNA) molecule of the liver flukeFasciola hepatica (phylum Platyhelminthes, class Trematoda) has been determined, within which have been identified the genes for tRNA^ala, tRNA^asp, respiratory chain NADH dehydrogenase subunit I (ND1), tRNA^asn, tRNA^pro, tRNA^ile, tRNA^lys, ND3, tRNA^serAGN, tRNA^trp, and cytochromec oxidase subunit I (COI). The 11 genes are arranged in the order given and are all transcribed from the same strand of the molecule. The overall order of theF. hepatica mitochondrial genes differs from what is found in other metazoan mtDNAs. All of the sequenced tRNA genes except the one for tRNA^serAGN can be folded into a secondary structure with four arms resembling most other metazoan mitochondrial tRNAs, rather than the tRNAs that contain a TψC arm replacement loop, found in nematode mtDNAs. TheF. hepatica mitochondrial tRNA^serAGN gene contains a dihydrouridine arm replacement loop, as is the case in all other metazoan mtDNAs examined to date. AGA and AGG are found in theF. hepatica mitochondrial protein genes and both codons appear to specify serine. These findings concerningF. hepatica mtDNA indicate that both a dihydrouridine arm replacement loop-containing tRNA^serAGN gene and the use of AGA and AGG codons to specify serine must first have occurred very early in, or before, the evolution of metazoa.     

Complete DNA Sequence of the Mitochondrial Genome of the Black Chiton, Katharina Tunicata

The DNA sequence of the 15,532-base pair (bp) mitochondrial DNA (mtDNA) of the chiton Katharina tunicata has been determined. The 37 genes typical of metazoan mtDNA are present: 13 for protein subunits involved in oxidative phosphorylation, 2 for rRNAs and 22 for tRNAs. The gene arrangement resembles those of arthropods much more than that of another mollusc, the bivalve Mytilus edulis. Most genes abut directly or overlap, and abbreviated stop codons are inferred for four genes. Four junctions between adjacent pairs of protein genes lack intervening tRNA genes; however, at each of these junctions there is a sequence immediately adjacent to the start codon of the downstream gene that is capable of forming a stem-and-loop structure. Analysis of the tRNA gene sequences suggests that the D arm is unpaired in tRNA(ser(AGN)), which is typical of metazoan mtDNAs, and also in tRNA(ser(UCN)), a condition found previously only in nematode mtDNAs. There are two additional sequences in Katharina mtDNA that can be folded into structures resembling tRNAs; whether these are functional genes is unknown. All possible codons except the stop codons TAA and TAG are used in the protein-encoding genes, and Katharina mtDNA appears to use the same variation of the mitochondrial genetic code that is used in Drosophila and Mytilus. Translation initiates at the codons ATG, ATA and GTG. A + T richness appears to have affected codon usage patterns and, perhaps, the amino acid composition of the encoded proteins. A 142-bp non-coding region between tRNA(glu) and CO3 contains a 72-bp tract of alternating A and T.     

The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus

We determined the complete 14,985-nt sequence of the mitochondrial DNA of the horseshoe crab Limulus polyphemus (Arthropoda: Xiphosura). This mtDNA encodes the 13 protein, 2 rRNA, and 22 tRNA genes typical for metazoans. The arrangement of these genes and about half of the sequence was reported previously; however, the sequence contained a large number of errors, which are corrected here. The two strands of Limulus mtDNA have significantly different nucleotide compositions. The strand encoding most mitochondrial proteins has 1. 25 times as many A's as T's and 2.33 times as many C's as G's. This nucleotide bias correlates with the biases in amino acid content and synonymous codon usage in proteins encoded by different strands and with the number of non-Watson-Crick base pairs in the stem regions of encoded tRNAs. The sizes of most mitochondrial protein genes in Limulus are either identical to or slightly smaller than those of their Drosophila counterparts. The usage of the initiation and termination codons in these genes seems to follow patterns that are conserved among most arthropod and some other metazoan mitochondrial genomes. The noncoding region of Limulus mtDNA contains a potential stem-loop structure, and we found a similar structure in the noncoding region of the published mtDNA of the prostriate tick Ixodes hexagonus. A simulation study was designed to evaluate the significance of these secondary structures; it revealed that they are statistically significant. No significant, comparable structure can be identified for the metastriate ticks Rhipicephalus sanguineus and Boophilus microplus. The latter two animals also share a mitochondrial gene rearrangement and an unusual structure of mt-tRNA(C) that is exactly the same association of changes as previously reported for a group of lizards. This suggests that the changes observed are not independent and that the stem-loop structure found in the noncoding regions of Limulus and Ixodes mtDNA may play the same role as that between trnN and trnC in vertebrates, i.e., the role of lagging strand origin of replication.     

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features.

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.     

Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAS and has a gene arrangement relatable to those of coelomate metazoans

The complete mitochondrial DNA (mtDNA) of the nematode Trichinella spiralis has been amplified in four overlapping fragments and 16,656 bp of its sequence has been determined. This sequence contains the 37 genes typical of metazoan mtDNAs, including a putative atp8, which is absent from all other nematode mtDNAs examined. The genes are transcribed from both mtDNA strands and have an arrangement relatable to those of coelomate metazoans, but not to those of secernentean nematodes. All protein genes appear to initiate with ATN codons, typical for metazoans. Neither TTG nor GTT start codons, inferred for several genes of other nematodes, were found. The 22 T. spiralis tRNA genes fall into three categories: (i) those with the potential to form conventional "cloverleaf" secondary structures, (ii) those with TPsiC arm + variable arm replacement loops, and (iii) those with DHU-arm replacement loops. Mt-tRNA(R) has a 5'-UCG-3' anticodon, as in most other metazoans, instead of the very unusual 5'-ACG-3' present in the secernentean nematodes. The sequence also contains a large repeat region that is polymorphic in size at the population and/or individual level.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

The mitochondrial genome of the gymnosperm Cycas taitungensis contains a novel family of short interspersed elements, Bpu sequences, and abundant RNA editing sites

The mtDNA of Cycas taitungensis is a circular molecule of 414,903 bp, making it 2- to 6-fold larger than the known mtDNAs of charophytes and bryophytes, but similar to the average of 7 elucidated angiosperm mtDNAs. It is characterized by abundant RNA editing sites (1,084), more than twice the number found in the angiosperm mtDNAs. The A + T content of Cycas mtDNA is 53.1%, the lowest among known land plants. About 5% of the Cycas mtDNA is composed of a novel family of mobile elements, which we designated as "Bpu sequences." They share a consensus sequence of 36 bp with 2 terminal direct repeats (AAGG) and a recognition site for the Bpu 10I restriction endonuclease (CCTGAAGC). Comparison of the Cycas mtDNA with other plant mtDNAs revealed many new insights into the biology and evolution of land plant mtDNAs. For example, the noncoding sequences in mtDNAs have drastically expanded as land plants have evolved, with abrupt increases appearing in the bryophytes, and then in the seed plants. As a result, the genomic organizations of seed plant mtDNAs are much less compact than in other plants. Also, the Cycas mtDNA appears to have been exempted from the frequent gene loss observed in angiosperm mtDNAs. Similar to the angiosperms, the 3 Cycas genes nad1, nad2, and nad5 are disrupted by 5 group II intron squences, which have brought the genes into trans-splicing arrangements. The evolutionary origin and invasion/duplication mechanism of the Bpu sequences in Cycas mtDNA are hypothesized and discussed.     

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.     

Mitochondrial genomes of two demosponges provide insights into an early stage of animal evolution

Mitochondrial DNA (mtDNA) of multicellular animals (Metazoa) is typically a small ( approximately 16 kbp), circular-mapping molecule that encodes 37 tightly packed genes. The structures of mtDNA-encoded transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs) are usually highly unorthodox, and proteins are translated with multiple deviations from the standard genetic code. In contrast, mtDNA of the choanoflagellate Monosiga brevicollis, the closest unicellular relative of animals, is four times larger, contains 1.5 times as many genes, and lacks mentioned peculiarities of animal mtDNA. To investigate the evolutionary transition that led to the specific organization of metazoan mtDNA, we determined complete mitochondrial sequences from the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. We found that poriferan mtDNAs resemble those of other animals in their compact organization, lack of introns, and a well-conserved animal-like gene order. Yet, they contain several extra genes, encode bacterial-like rRNAs and tRNAs, and use a minimally derived genetic code. Our findings suggest that the evolution of the typical metazoan mtDNA has been a multistep process in which the compact genome organization and the reduced gene content were established prior to the reduction of tRNA and rRNA structures and the introduction of multiple changes of the translation code.