The complete sequence of the mitochondrial genome of Nautilus macromphalus (Mollusca: Cephalopoda)
© Boore; licensee BioMed Central Ltd. 2006
Received: 11 January 2006
Accepted: 19 July 2006
Published: 19 July 2006
Mitochondria contain small genomes that are physically separate from those of nuclei. Their comparison serves as a model system for understanding the processes of genome evolution. Although complete mitochondrial genome sequences have been reported for more than 600 animals, the taxonomic sampling is highly biased toward vertebrates and arthropods, leaving much of the diversity yet uncharacterized.
The mitochondrial genome of the bellybutton nautilus, Nautilus macromphalus, a cephalopod mollusk, is 16,258 nts in length and 59.5% A+T, both values that are typical of animal mitochondrial genomes. It contains the 37 genes that are almost universally found in animal mtDNAs, with 15 on one DNA strand and 22 on the other. The arrangement of these genes can be derived from that of the distantly related Katharina tunicata (Mollusca: Polyplacophora) by a switch in position of two large blocks of genes and transpositions of four tRNA genes. There is strong skew in the distribution of nucleotides between the two strands, and analysis of this yields insight into modes of transcription and replication. There is an unusual number of non-coding regions and their function, if any, is not known; however, several of these demark abrupt shifts in nucleotide skew, and there are several identical sequence elements at these junctions, suggesting that they may play roles in transcription and/or replication. One of the non-coding regions contains multiple repeats of a tRNA-like sequence. Some of the tRNA genes appear to overlap on the same strand, but this could be resolved if the polycistron were cleaved at the beginning of the downstream gene, followed by polyadenylation of the product of the upstream gene to form a fully paired structure.
Nautilus macromphalus mtDNA contains an expected gene content that has experienced few rearrangements since the evolutionary split between cephalopods and polyplacophorans. It contains an unusual number of non-coding regions, especially considering that these otherwise often are generated by the same processes that produce gene rearrangements. The skew in nucleotide composition between the two strands is strong and associated with the direction of transcription in various parts of the genomes, but a comparison with K. tunicata implies that mutational bias during replication also plays a role. This appears to be yet another case where polyadenylation of mitochondrial tRNAs restores what would otherwise be an incomplete structure.
Animal mitochondrial DNA (mtDNA) is nearly always a closed circular molecule and, with a few exceptions [e.g. [1–4]], contains the same 37 genes, specifying 13 proteins, two ribosomal RNAs, and 22 tRNAs . Sequences of these diminutive genomes have been broadly used to address phylogenetic questions ranging from the population [6, 7] to the interphylum [8–11] levels and to model many processes of genome evolution [12, 13]. Although there are exceptions, most mtDNAs contain no introns and are between 14 and 17 kb. Typically there are few intergenic nucleotides except for a single large non-coding region generally thought to contain elements that regulate replication and transcription . Occasionally non-coding regions have been found that contain repeated elements  or contain pseudogenes [12, 16] or that may be remnants of duplicated regions, perhaps those that mediate gene rearrangements [12, 16, 17]. Gene rearrangements tend to be uncommon and to occur in a saltatory manner [see ]. The "universal" genetic code has been modified in many animal lineages, to include the use of alternative start codons and abbreviated stop codons [18, 19]. In some mtDNAs there is pronounced skew in nucleotide composition, often with one strand being rich in G and T and the other in A and C . Post-transcriptional modification of nucleotides has been observed for tRNAs [21, 22].
Little study has been done to date on mollusk mtDNAs compared to those of vertebrates or arthropods , but it is already apparent that mollusks exhibit much variation in the features of their mitochondrial genomes, including losses and gains of genes , atypically large amounts of duplicated or non-coding nucleotides [15, 24], highly rearranged genomes [2, 25], and an unusual pattern of passage termed doubly uniparental inheritance [26, 27]. This is furthered here by reporting and comparing the features of the mitochondrial genome of the first nautiloid to be so studied, N. macromphalus (Mollusca: Cephalopoda).
Nautiloids were once abundant and diverse in the Paleozoic seas, but only a handful of species remain. They are part of the molluscan class Cephalopoda, which otherwise contains octopi, squid, and cuttlefish. They are the earliest diverging lineage of this group and are often considered to be "living fossils" since living forms seem to have changed little from their ancient ancestors. They live in spiral-shaped shells which are filled with gas to control buoyancy and they move about by squirting jets of water. They are carnivorous, using their many grooved tentacles to grasp prey and pass it to their mouth, where a beak-like jaw tears it and passes it to the shredding radula. They live throughout the Southwest Pacific Ocean, at depths as great as 610 meters, and traverse a great range, as shallow as 90 meters, apparently in search of prey.
Complete mtDNA sequences have been determined for 23 mollusks (see Additional file 1), including a representative (Katharina tunicata)  of a basal group (Polyplacophora). This sampling includes seven other cephalopods: Octopus vulgaris , Loligo bleekeri , Todarodes pacificus [29, 31], O. ocellatus, Sepioteuthis lessoniana, Watasenia scintillans, and Sepia officinalis . Comparisons of the features of the N. macromphalus mtDNA with those of some other mollusks are presented here.
Results and discussion
Gene content and organization
In total, there are now available complete mtDNA sequences from eight cephalopod species to compare. In addition to N. macromphalus and O. vulgaris, these are the squids Loligo bleekeri , Todarodes pacificus [29, 31], Watasenia scintillans , and Sepioteuthis lessoniana , the octopus O. ocellatus , and the cuttlefish Sepia officinalis . O. ocellatus shares an identical gene arrangement with O. vulgaris. Two of the squids, L. bleekeri and S. lessoniana, share a nearly identical gene arrangement (differing only by a transposition of one block of genes: trnI, -rrnL, -trnV, -rrnS, -trnW [minus symbol indicates opposite transcriptional orientation]). This gene arrangement, plus another separately rearranged in S. officinalis, are highly derived and each shares only a few blocks of colinearity with the more conserved gene order of N. macromphalus mtDNA. All of these cephalopod mtDNAs have the same gene content except for W. scintillans and T. pacificus, the two representatives of the group Oegopsida. These two mtDNAs have a nearly identical gene arrangement, differing only in the position of trnM, that is highly rearranged from those of other mollusks, and contain duplicated copies of cox1, cox2, cox3, atp6, atp8, and trnD, such that they contain genes for a total of 18 proteins, 2 rRNAs, and 23 tRNAs. In all of these studied cephalopod mtDNAs, all genes retain the same transcriptional orientation, that is, all rearrangements are transpositions and none are inversions. Akasaki et al.  provide a comprehensive and well reasoned review of this pattern of arrangements, including proposals for mechanism of rearrangement, the role of the many large, non-coding regions, and evidence for concerted evolution of duplicated genes.
Gene initiation and termination
Mitochondrial genomes often use a variety of non-standard initiation codons , but N. macromphalus mtDNA has only one type of deviation; three genes (nad3, nad4, and nad5) start with GTG and all others use the standard ATG (Additional file 4). Seven genes have unambiguous termination codons, either TAG (atp6,cox1, nad5) or TAA (atp8, cox3, nad1, nad2). In four cases (cox2, cob, nad3, nad4) genes are probably abbreviated to a single T or to TA such that the excision of the adjacent, downstream tRNA from the polycistronic message leaves an mRNA that is polyadenylated to complete a TAA stop codon. However, in each of these cases, a complete stop codon is available if there is, alternatively, overlap of only one or two nucleotides with the downstream tRNA. Perhaps these act as a "backup" for cases where translation precedes message cleavage. The other two cases are more ambiguous. nad4L could have an abbreviated stop codon, but is inferred to overlap nad4 by seven nucleotides to the first legitimate stop codon, since overlap of this pair has been commonly observed for other mtDNAs, where they are thought to be translated as a bicistron. nad6 is inferred to overlap cob by eight nucleotides, perhaps suggesting that these are processed also as a bicistron, but could instead end on an abbreviated stop codon if there were some signal for message cleavage (i.e., other than a tRNA) that we do not recognize. Inferred in this way, all protein-encoding genes have lengths nearly identical to those of K. tunicata mtDNA (Additional file 2).
There are three cases where tRNA genes appear to overlap, and these potential structures suggest how this is resolved. trnL1(nag) appears to overlap trnL2(yaa) by only the former's discriminator nucleotide (A). trnQ appears to overlap trnW by two nucleotides. trnK appears to overlap trnA by four nucleotides, GGCT. These are well-paired in the potential structure of tRNA(A), but these four correspond to two G-T pairs, one mismatch, and the discriminator nucleotide of tRNA(K). It appears for each case that cleavage to form a complete downstream tRNA followed by (poly)adenylation of the upstream tRNA (as has been demonstrated for some mitochondrial tRNAs ) would yield fully formed, well-paired structures for all. This is illustrated in Figure 3 by lower case, parenthetical letter "a" appended to the genome-encoded nucleotide to indicate likely nucleotides in the actual transcript.
Usually T is in the first anticodon position for tRNAs that recognize either four-fold degenerate codon families or to specifically recognize NNR codons; G is usually in this position only to specifically recognize NNY codons. (Due to the convention of always drawing RNAs from 5' to 3' in orientation, the first nucleotide listed for an anticodon pairs with the last nucleotide of a codon.) All but two of the N. macromphalus mitochondrial tRNAs follow this pattern. One exception is tRNA(M), which has the anticodon CAT (to recognize both ATG and ATA), as is almost universally the case for all animal mitochondrial systems. In some cases the C is known to be post-transcriptionally modified to 5-formylcytidine to enable the necessary pairing with the ATA codon . However, it is less common that the tRNA(S) expected to recognize codon AGN has a GCT anticodon, since this requires the G to pair with all four nucleotides in the wobble position of AGN codons. It is clear the AGA and AGG codons are being used and are not stop codons (as is the case in vertebrate mtDNAs), since they appear in the reading frames of protein encoding genes 117 times. GCT is used as the tRNA(S) anticodon for all of the cephalopods with complete mtDNA sequences (above), and it is likely that this anticodon is modified post-transcriptionally for all, as is known to occur for the Loligo bleekeri tRNA(S), for which the G is modified to 7-methylguanosine .
Number of nucleotides at gene boundaries. Negative numbers refer to overlapping nucleotides.
cox1 – cox2
trnH – nad4
cox2 – trnD
nad4 – nad4L
trnD – atp8
nad4L – trnS2(nga)
atp8 – trnF
trnS2(nga) – cob
trnF – trnL2(yaa)
cob – nad6
trnL2(yaa) – trnL1(nag)
nad6 – trnP
trnL1(nag) – rrnL
trnP – nad1
rrnL – trnV
nad1 – trnE
trnV – rrnS
trnE – cox3
rrnS – trnM
cox3 – trnK
trnM – trnC
trnK – trnA
trnC – trnY
trnA – trnR
trnY – trnW
trnR – trnN
trnW – trnQ
trnN – trnI
trnQ – trnT
trnI – nad3
trnT – trnG
nad3 – trnS1(nct)
trnG – atp6
trnS1(nct) – nad2
atp6 – nad5
nad2 – cox1
nad5 – trnH
In the largest non-coding region, between trnQ and trnT, and beginning adjacent to a (CA)13 run (see below), there are six repeats of a 62 nucleotide element followed by a partial repeat of 39 nucleotides. Within this are five overlapping regions that have potential for forming tRNA-like structures (Fig. 3). The anticodon portion of these structures is AGT, which would pair with codon ACT (or perhaps ACN) to specify threonine. However, having A in this anticodon position would be very unusual and there is little sequence similarity to trnT (or any other tRNA).
Tandem repeats of CA are common, with (CA)3 in each of the intergenic regions of trnA-trnR and trnG-atp6 and an especially noteworthy (CA)13 in the region between trnQ and trnT. Homopolymer runs of T10, nine C9, and A20 are in the regions trnQ-trnT, trnG-atp6, and trnE-cox3, respectively. Non-coding, non-functional portions of mtDNA are generally eliminated rapidly , presumably due to selection for small size at the point of entry into the primordial germ plasm during embryogenesis , but whether these or any particular motif plays any role in regulating replication of transcription awaits experimentation.
Base composition and codon usage
Base composition of a chiton and cephalopod mtDNAs. This refers in all cases to the strand that is in the sense orientation for cox1
Codon usage for the 13 mitochondrial proteins of Nautilus macromphalus. The total number of codons is 3711. Stop codons were not included in this count. Here the plus strand refers arbitrarily to the one that contains cox1.
All protein encoding genes
Plus strand only
Minus strand only
All protein encoding genes
Plus strand only
Minus strand only
The mitochondrial genome of the chiton, K. tunicata, contrasts with this. Although the gene arrangement is quite similar, here the pattern of bias is opposite in two different respects. First, it is the plus-strand genes that have strong skew in nucleotide composition, with the minus-strand genes being nearly neutral for this bias. Secondly, the bias for these is strongly toward G+T for the reported strand. Here again, the sharp transitions in base composition are flanked by non-coding regions at least 20 nucleotides in length, which could potentially serve as signaling elements for transcription or replication.
Such skews with one strand being rich in A+C and the other rich in G+T are common for mitochondrial genomes . (See  for a review of the proposed causes and an analysis specific to mtDNAs.) This is thought to be due predominantly to the commonality of deamination of adenine and cytosine nucleotides in single-stranded DNA [37–39] which appears transiently during replication and transcription. The relative contribution of these two processes remains unclear , as each accounts for one strand being displaced by the nascent DNA or RNA, respectively. (Although this is not without controversy in the case of mitochondrial replication [41–44]). Deaminated adenine forms hypoxanthine, which pairs with cytosine (rather than thymine) and deaminated cytosine forms uracil, which pairs with adenine (rather than guanine). Therefore, the displaced strand, existing in single-stranded form for sometimes protracted periods, tends to become rich in G+T (the analogs of hypoxanthine and uracil) and its complementary strand, therefore, becomes rich in A+C.
Since N. macromphalus and K. tunicata mtDNAs each have sharp boundaries in base compositional bias that correspond so precisely to shifts in transcriptional orientation, it appears that lesions in the displaced strand during transcription are an important contribution. On the other hand, the contrast in the bias being strong for the minus-strand genes of N. macromphalus and for the plus-strand genes of K. tunicata shows that some other factor must be at work.
According to the more long-standing and broadly accepted model of mtDNA replication  (but see [42–44]), and demonstrated for the few cases where it has been studied, replication of mtDNA is very slow and very asymmetrical, with one strand in single stranded form for a protracted period, so this may be an important factor in strand compositional bias. The nucleotide skew between the two mitochondrial strands is expected to be a combination of various factors, and one could imagine a model whereby a reversal between N. macromphalus and K. tunicata in which strand is leading during replication could account for their differing skew patterns. If replication in K. tunicata mtDNA were to proceed first in the rightward direction according to Figure 4, then the bias introduced during replication would make the reported strand rich in G+T. This would be reinforced by biases introduced during transcription in the regions of the plus-strand genes, causing especially high bias, and countered by the biases introduced during transcription in the regions of the minus-strand genes, causing them to approximately cancel out. If N. macromphalus mtDNA replication were to proceed in the opposite direction, right-to-left as in Figure 4, then the effect would be the opposite, with skew generated by mutational bias during replication reinforcing that from transcription of minus-strand genes and opposing that from transcription of plus-strand genes, and accounting for the patterns shown in Figure 4.
It is not clear whether the isolated trnG is transcribed individually or is part of the transcription unit that otherwise ends at nad5. Separating trnG and nad5 is a single plus-strand gene, atp6, and it is possible that this is transcribed in reverse as part of the larger transcription unit, with this antisense message excised and degraded. When considering only the composition of the third positions of four-fold degenerate codon families, G and T comprise 0.04 and 0.24 of atp6, values nearly identical to 0.08 and 0.26 for the other plus-strand genes collectively. However, A and C are 0.24 and 0.48 for atp6 vs. 0.41 and 0.24, respectively, for the other plus-strand genes. Perhaps this indicates a modifying force for mutational bias, perhaps the regular reverse transcription of the gene. On the other hand, trnG is flanked by large blocks of non-coding sequence which could potentially be signals for initiating and terminating transcription for this individual gene.
Table 2 compares the mtDNA size, base composition, A+T-richness, and strand skews for K. tunicata and the eight cephalopod species with complete mtDNA sequences. The other cephalopods all have strand skew measures that are in the same direction, but of lesser magnitude, than N. macromphalus, and all of these cephalopods have strand skews in opposite direction from that of the outgroup K. tunicata. There also appears to be a trend for larger mtDNAs in the cephalopods and, for the octopus and squid lineages, for greater A+T-richness after the split of that leading to the nautiloids.
Potential signaling elements
An attempt was made to find potential regulatory sequence elements by comparing all pairs of non-coding regions that are greater than 13 nucleotides in length (Additional file 3) for any blocks 10 nucleotides or longer with identity at least 80% while considering both strands. In addition to the homopolymer runs and dinucleotide repeats discussed above and underlined in Additional file 3, three elements were identified, all associated with reversals in transcriptional orientation. In the largest non-coding region between the oppositely oriented trnQ and trnT is the sequence TTAAAACAA, also found in the region between atp6 and nad5. Although both are at a point where transcriptional orientation reverses, the first case is of genes arranged head-to-head but the second of genes tail-to-tail. Also in the trnQ-trnT region is the sequence CCNATTTTA which is also found in the region between trnT and trnG; again in the first case the genes are head-to-head and in the second tail-to-tail. The sequence ATAACAAAACTA occurs in the region between trnE and cox3 and also between trnG and atp6, in each case pairs of genes arranged head-to-head. (There are three cases total where genes are arranged head-to-head, these two plus trnQ-trnT. There are three cases total where genes are arranged tail-to-tail, trnT-trnG, atp6-nad5, and atp8-trnF.) None of these sequences are present in K. tunicata mtDNA and none are present in the non-coding regions of any of the other studied cephalopods.
A comparison was also made between each non-coding region of N. macromphalus and each non-coding region of all of the other cephalopod mtDNAs greater than 19 nucleotides in length for any blocks of length 20 or greater matching at least 70%. Although some matches were found, none were consistent across all (or even most) species. Lastly, a search was made for all available cephalopod mtDNAs for long stretches of alternating CA or TA, suggested to play a role in regulation of replication and/or transcription. Of note is that N. macromphalus has several regions of alternating CA (Additional file 3), the longest of which is (CA)13. Only two of the other cephalopods, L. bleekeri and T. pacificus have any as long as (CA)4. In contrast, while N. macromphalus has no regions of alternating TA longer than (TA)3 (which occurs in eight places), each of the other cephalopods has many such regions at least of length (TA)9 (T. pacificus and W. scintillans), and some as long as (TA)22 (S. lessoniana). (The longest alternating TA for L. bleekeri is (TA)14, for O. ocellatus is (TA)16, and for both S. officinalis and O. vulgaris is (TA)12.) Of course, it is possible that actual regulatory elements may be more complex and difficult to identify.
To date, complete mtDNA sequences had been determined for 23 mollusks, a very small sampling compared to those available for vertebrates or arthropods . Even these few studies have revealed that mollusks' mtDNAs have much variation in their features, including losses and gains of genes , unusually large amounts of duplicated or non-coding nucleotides [15, 24], numerous gene rearrangements [2, 25], and doubly uniparental inheritance [26, 27]. By contrast, the mtDNA of the cephalopod Nautilus macromphalus is fairly typical in many respects, with a size, gene content, and A+T-richness similar to those most common for animal mtDNAs. There have been only a few gene arrangements in this lineage even since its divergence from the basal mollusk group Polyplacophora, and these rearrangements can be confidently polarized among the two lineages by comparing them to mtDNAs of less related animals.
There is strong skew in the distribution of nucleotides between the two strands and it appears that biases in mutational spectrum during both transcription and replication are responsible for this. Compared with most animal mtDNAs, there are a large number of non-coding regions. Although their functions, if any, are not known, the fact that several are at positions of abrupt shift in nucleotide skew and that some contain identical sequence elements suggests that they may contain regulatory signals for transcription and/or replication. This appears to be another example where polyadenylation of tRNAs creates part of the amino-acyl acceptor stem. These, and other features can be interpreted in detail for the systems of these diminutive genomes, and further sampling of complete mtDNA sequences across the tree of life promises to provide insights into general aspects of genome evolution.
Testis tissue, stored for longer than a decade at -80°C, but without any record of which species of Nautilus had been sampled, was the gift of Wesley Brown. Fortunately, GenBank contains a short fragment of mitochondrial rRNA for each of the six species of the genus, and this was used for specific identification. The 401 nucleotides in common for this sample and the GenBank records were compared to determine that only two positions differ with the record of N. macromphalus (this 0.5% difference is presumably due to intraspecies polymorphism), whereas all others differ by from 16 to 24 positions; therefore, it appears that this sample was from N. macromphalus.
Mitochondrial DNA was isolated from approximately 1 g of this tissue by first grinding in liquid nitrogen using a mortar and pestle. This powder was dissolved in 14 ml of homogenization buffer (210 mM mannitol, 70 mM sucrose, 50 mM Tris HCl-pH 75, 3 mM CaCl2) and processed using a Tissuemizer T-25 (Tekmar) with three strokes of five seconds each. Membranes were lysed by adding 1/10 volume of 20% SDS and incubating for 20 min at RT. A 1/6 volume of saturated CsCl in water was added and this mixture incubated on ice for 15 min. Debris was pelleted at 17,000 × G for 10 minutes at 4°C. Propidium iodide was added to the collected supernatant to a final concentration of 500 μg/ml and the CsCl concentration was adjusted to a density of 1.57 g/ml. Nuclear and mitochondrial DNA were separated by density gradient centrifugation in a VTi65 rotor at 55,000 × G for 15 hours at 21°C. Although no mitochondrial band was visible in the gradient, the region from about 2–10 mm below the nuclear band was collected using a needle. This was then extracted multiple times with water-saturated butanol to remove the propidium iodide and dialyzed against TE for 24 hours with three buffer changes to remove the CsCl, leaving the sample in a 100 μl volume.
This product was used in PCR as in  to amplify first several short fragments of cox1, rrnL, and cob using primers found in [45–47]. The fragment of cox1 was cloned into pBluescript (Stratagene) that had been digested with Eco RV, T-tailed using Taq polymerase, and gel purified using GeneClean (QBiogene). A successful recombinant clone was selected and DNA prepared using standard techniques. The other fragments were purified by three serial passages through an Ultrafree (NMWL 30,000) spin column (Millipore) and sequenced directly. The sequences of these fragments were determined using an ABI377 automated DNA sequencer with BigDye chemistry (Applied Biosystems) according to supplier's instructions.
Primers were designed to known sequences for use in long PCR  with rTth-XL polymerase (Applied Biosystems) according to supplier's instructions, sometimes combined with primers to conserved mtDNA regions. Generously overlapping fragments were amplified from cox1-nad1 (using conserved nad1 primer CCTGATACTAATTCAGATTCTCCTTC), nad1-cob, cob-rrnL, and rrnL – cox1 (using conserved primer 16SARL ), jointly comprising the entire mtDNA. Because there was no information available for the gene arrangement, many combinations of primers were tried, but only these reactions gave bright, singular bands during electrophoretic analysis. Sequence was determined for each as above, then by primer walking through each fragment. To ensure accuracy, all sequence was determined from both strands. Sequencing reads were assembled manually and quality verified by eye using Sequence Navigator (Applied Biosystems).
Gene annotation and analysis
Genes encoding rRNAs and proteins were identified by matching nucleotide or inferred amino acid sequences to those of K. tunicata mtDNA  through the use of MacVector (Accelrys). Since it is not possible to precisely determine the ends of rRNA genes by sequence data alone, they were assumed to extend to the boundaries of flanking genes. Each protein gene was inferred to begin at an eligible initiation codon nearest to the beginning of its alignment with homologous genes that does not cause overlap with the preceding gene. In five cases, an abbreviated stop codon was inferred where cleavage of a downstream tRNA from the transcript would leave a partial codon of T or TA, such that subsequent mRNA polyadenylation could generate a TAA stop codon; however, in each of these cases, if the reading frame extended through the first legitimate stop codon there would be only a short overlap with the downstream gene. Genes for tRNAs were identified by eye, generically by their ability to fold into a cloverleaf structure and specifically by anticodon sequence. Subsequent analyses, such as counting anticodon usage, calculating nucleotide frequencies and strand skew values, and identifying repeated elements, were performed using MacVector (Accelrys).
- cox1 :
cox2,cox3, cytochrome oxidase subunit I, II, and III protein genes
- cob :
cytochrome b gene
- atp6 :
atp8, ATP synthase subunit 6 and 8 genes
- nad1 :
nad2,nad3,nad4,nad4L,nad5,nad6, NADH dehydrogenase subunit 1–6, 4L genes
- trnA :
trnC,trnD,trnE,trnF,trnG,trnH,trnI,trnK,trnL1(nag),trnL2(yaa),trnM,trnN,trnP,trnQ,trnR,trnS1(nct),trnS2(nga),trnT,trnV,trnW,trnY, transfer RNA genes designated by the one-letter code for the specified amino acid
- in cases where there is more than one tRNA for a particular amino acid:
they are differentiated by numeral and with anticodon (which is maximally ambiguous, e.g. "nag" rather than "tag", to allow recognizing homology with those of other organisms).
I am grateful to Wesley Brown for Nautilus tissue and for many years of guidance and encouragement. This work was supported by funding from the National Science Foundation (DEB-9807100, EAR-0342392, DEB-0089624) and was performed partly under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231.
- Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR: The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992, 130: 471-498.PubMedPubMed CentralGoogle Scholar
- Hoffmann RJ, Boore JL, Brown WM: A novel mitochondrial genome organization for the blue mussel, Mytilus edulis. Genetics. 1992, 131: 397-412.PubMedPubMed CentralGoogle Scholar
- Beagley CT, Okimoto R, Wolstenholme DR: The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): Introns, a paucity of tRNA genes, and a near-standard genetic code. Genetics. 1998, 148: 1091-1108.PubMedPubMed CentralGoogle Scholar
- Helfenbein KG, Fourcade HM, Vanjani RG, Boore JL: The mitochondrial genome of Paraspadella gotoi is highly reduced and reveals that chaetognaths are a sister-group to protostomes. Proc Natl Acad Sci USA. 2004, 101 (29): 10639-10643. 10.1073/pnas.0400941101.PubMedPubMed CentralView ArticleGoogle Scholar
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.PubMedPubMed CentralView ArticleGoogle Scholar
- Nyakaana S, Arctander P, Siegismund H: Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci. Heredity. 2002, 89 (2): 90-98. 10.1038/sj.hdy.6800110.PubMedView ArticleGoogle Scholar
- Ingman M, Kaessmann H, Pääbo S, Gyllensten U: Mitochondrial genome variation and the origin of modern humans. Nature. 2001, 408: 708-713.Google Scholar
- Smith MJ, Arndt A, Gorski S, Fajber E: The phylogeny of echinoderm classes based on mitochondrial gene arrangements. J Mol Evol. 1993, 36: 545-554. 10.1007/BF00556359.PubMedView ArticleGoogle Scholar
- Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM: Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements. Nature. 1995, 376: 163-165. 10.1038/376163a0.PubMedView ArticleGoogle Scholar
- Boore JL, Brown WM: Big trees from little genomes: Mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998, 8 (6): 668-674. 10.1016/S0959-437X(98)80035-X.PubMedView ArticleGoogle Scholar
- Helfenbein KG, Boore JL: The mitochondrial genome of Phoronis architecta – Comparisons demonstrate that phoronids are lophotrochozoan protostomes. Mol Biol Evol. 2004, 21 (1): 153-157. 10.1093/molbev/msh011.PubMedView ArticleGoogle Scholar
- Mueller RL, Boore JL: Molecular mechanisms of extensive mitochondrial gene rearrangement in plethodontid salamanders. Mol Biol Evol. 2005, 22: 2104-2112. 10.1093/molbev/msi204.PubMedView ArticleGoogle Scholar
- Helfenbein KG, Brown WM, Boore JL: The complete mitochondrial genome of a lophophorate, the brachiopod Terebratalia transversa. Mol Biol Evol. 2001, 18 (9): 1734-1744.PubMedView ArticleGoogle Scholar
- Shadel GS, Clayton DA: Mitochondrial DNA maintenance in vertebrates. Annu Rev Biochem. 1997, 66: 409-435. 10.1146/annurev.biochem.66.1.409.PubMedView ArticleGoogle Scholar
- Rigaa A, Monnerot M, Sellos D: Molecular cloning and complete nucleotide sequence of the repeated unit and flanking gene of the scallop Pecten maximus mitochondrial DNA: Putative replication origin features. J Mol Evol. 1995, 41: 189-195. 10.1007/BF00170672.PubMedView ArticleGoogle Scholar
- Arndt A, Smith MJ: Mitochondrial gene rearrangement in the sea cucumber genus Cucumaria. Mol Biol Evol. 1998, 15 (8): 1009-1016.PubMedView ArticleGoogle Scholar
- Boore JL: The duplication/random loss model for gene rearrangement exemplified by mitochondrial genomes of deuterostome animals. Comparative Genomics. Edited by: Sankoff D, Nadeau J. 2000, Computational Biology Series, Kluwer Academic Publishers, Dordrecht, Netherlands, 1: 133-147.View ArticleGoogle Scholar
- Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290: 470-474. 10.1038/290470a0.PubMedView ArticleGoogle Scholar
- Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Intl Rev Cytology. 1992, 141: 173-216.View ArticleGoogle Scholar
- Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995, 41: 353-358. 10.1007/BF01215182.PubMedView ArticleGoogle Scholar
- Lavrov D, Brown WM, Boore JL: A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci USA. 2000, 97 (25): 13738-13742. 10.1073/pnas.250402997.PubMedPubMed CentralView ArticleGoogle Scholar
- Yokobori S, Pääbo S: Transfer RNA editing in land snail mitochondria. Proc Natl Acad Sci. 1995, 92 (22): 10432-10435. 10.1073/pnas.92.22.10432.PubMedPubMed CentralView ArticleGoogle Scholar
- Vallès Y, Boore JL: Lophotrochozoan mitochondrial genomes. Integrative Comp Biol. 2006,Google Scholar
- Fuller KM, Zouros E: Dispersed length polymorphism of mitochondrial DNA in the scallop Placopecten magellanicus (Gmelin). Curr Genet. 1993, 23: 365-369. 10.1007/BF00310901.PubMedView ArticleGoogle Scholar
- Boore JL, Medina M, Rosenberg LA: Complete sequences of two highly rearranged molluscan mitochondrial genomes, those of the scaphopod Graptacme eborea and of the bivalve Mytilus edulis. Mol Biol Evol. 2004, 21 (8): 1492-1503. 10.1093/molbev/msh090.PubMedView ArticleGoogle Scholar
- Stewart DT, Saavedra C, Stanwood RR, Ball AO, Zouros E: Male and female mitochondrial DNA lineages in the Blue Mussel (Mytilus edulis) species group. Mol Biol Evol. 1995, 12: 735-747.PubMedGoogle Scholar
- Passamonti M, Boore JL, Scali V: Molecular evolution and recombination in gender-associated mitochondrial DNAs of the Manila clam Tapes philippinarum. Genetics. 2003, 164: 603-611.PubMedPubMed CentralGoogle Scholar
- Boore JL, Brown WM: Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata. Genetics. 1994, 138: 423-443.PubMedPubMed CentralGoogle Scholar
- Yokobori S, Fukuda N, Nakamura M, Aoyama T, Oshima T: Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes. Mol Biol Evol. 2004, 21 (11): 2034-2046. 10.1093/molbev/msh227.PubMedView ArticleGoogle Scholar
- Tomita K, Yokobori S, Oshima T, Ueda T, Watanabe K: The cephalopod Loligo bleekeri mitochondrial genome: Multiplied noncoding regions and transposition of tRNA genes. J Mol Evol. 2002, 54 (4): 486-500. 10.1007/s00239-001-0039-4.PubMedView ArticleGoogle Scholar
- Akasaki T, Nikaido M, Tsuchiya K, Segawa S, Hasegawa M, Okada N: Extensive mitochondrial gene arrangements in coleoid Cephalopoda and their phylogenetic implications. Mol Phylogenet Evol. 2006, 38: 648-658. 10.1016/j.ympev.2005.10.018.PubMedView ArticleGoogle Scholar
- Tomita K, Ueda T, Ishiwa S, Crain PF, McCloskey JA, Watanabe K: Codon reading patterns in Drosophila melanogaster mitochondria based on their tRNA sequences: a unique wobble rule in animal mitochondria. Nucleic Acids Res. 1999, 27: 4291-4297. 10.1093/nar/27.21.4291.PubMedPubMed CentralView ArticleGoogle Scholar
- Tomita K, Ueda T, Watanabe K: 7-Methylguanosine at the anticodon wobble position of squid mitochondrial tRNA(Ser)GCU: molecular basis for assignment of AGA/AGG codons as serine in invertebrate mitochondria. Biochim Biophys Acta. 1998, 1399 (1): 78-82.PubMedView ArticleGoogle Scholar
- Ashley MV, Laipis PJ, Hauswirth WW: Rapid segregation of heteroplasmic bovine mitochondria. Nucleic Acids Res. 1989, 17 (18): 7325-7331.PubMedPubMed CentralView ArticleGoogle Scholar
- Mignotte F, Tourte M, Mounolou J-C: Segregation of mitochondria in the cytoplasm of Xenopus vitellogenic oocytes. Biol of the Cell. 1987, 60: 97-102.View ArticleGoogle Scholar
- Hassanin A, Léger N, Deutsch J: Evidence for multiple reversals of asymmetric mutational constraints during the evolution of the mitochondrial genome of Metazoa, and consequences for phylogenetic inferences. Syst Biol. 2005, 54: 277-298. 10.1080/10635150590947843.PubMedView ArticleGoogle Scholar
- Beletskii A, Bhagwat AS: Transcription-induced mutations: Increase in C to T mutations in the nontranscribed strand during transcription in Escherichia coli. Proc Natl Acad Sci. 1996, 93: 13919-13924. 10.1073/pnas.93.24.13919.PubMedPubMed CentralView ArticleGoogle Scholar
- Frederico LA, Kunkel TA, Shaw BA: A sensitive genetic assay for the detection of cytosine deamination: Determination of rate constant and the activation energy. Biochemistry. 1990, 29: 2532-2537. 10.1021/bi00462a015.PubMedView ArticleGoogle Scholar
- Sancar A, Sancar GB: DNA repair enzymes. Annu Rev Biochem. 1998, 57: 29-67. 10.1146/annurev.bi.57.070188.000333.View ArticleGoogle Scholar
- Francino MP, Ochman H: Strand asymmetries in DNA evolution. Trends Genet. 1997, 13 (6): 240-245. 10.1016/S0168-9525(97)01118-9.PubMedView ArticleGoogle Scholar
- Clayton DA: Replication and transcription of vertebrate mitochondrial DNA. Annu Rev Cell Biol. 1991, 7: 453-478. 10.1146/annurev.cb.07.110191.002321.PubMedView ArticleGoogle Scholar
- Holt IJ, Lorimer HE, Jacobs HT: Coupled leading- and lagging-strand synthesis of mammalian mitochondrial DNA. Cell. 2000, 100: 515-524. 10.1016/S0092-8674(00)80688-1.PubMedView ArticleGoogle Scholar
- Bogenhagen DF, Clayton DA: The mitochondrial DNA replication bubble has not burst. Trends Biochem Sci. 2003, 28: 357-360. 10.1016/S0968-0004(03)00132-4.PubMedView ArticleGoogle Scholar
- Yasukawa T, Yang MY, Jacobs HT, Holt IJ: A bidirectional origin of replication maps to the major noncoding region of human mitochondrial DNA. Mol Cell. 2005, 18: 651-662. 10.1016/j.molcel.2005.05.002.PubMedView ArticleGoogle Scholar
- Boore JL, Brown WM: 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. 2000, 17 (1): 87-106.PubMedView ArticleGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech. 1994, 3: 294-299.Google Scholar
- Palumbi SR: Nucleic acids II: The polymerase chain reaction. Molecular Systematics. Edited by: Hillis DM, Moritz C, Mable BK. 1996, Sinauer Associates, Sunderland, Massachusetts, USA, 205-247.Google Scholar
- Barnes WM: PCR amplification of up to 35-kb DNA with high fidelity and high yield from bacteriophage templates. Proc Natl Acad Sci USA. 1994, 91: 2216-2220. 10.1073/pnas.91.6.2216.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.