- Research article
- Open Access
The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features
BMC Genomics volume 7, Article number: 241 (2006)
Sequence data and other characters from mitochondrial genomes (gene translocations, secondary structure of RNA molecules) are useful in phylogenetic studies among metazoan animals from population to phylum level. Moreover, the comparison of complete mitochondrial sequences gives valuable information about the evolution of small genomes, e.g. about different mechanisms of gene translocation, gene duplication and gene loss, or concerning nucleotide frequency biases.
The Peracarida (gammarids, isopods, etc.) comprise about 21,000 species of crustaceans, living in many environments from deep sea floor to arid terrestrial habitats. Ligia oceanica is a terrestrial isopod living at rocky seashores of the european North Sea and Atlantic coastlines.
The study reveals the first complete mitochondrial DNA sequence from a peracarid crustacean. The mitochondrial genome of Ligia oceanica is a circular double-stranded DNA molecule, with a size of 15,289 bp. It shows several changes in mitochondrial gene order compared to other crustacean species. An overview about mitochondrial gene order of all crustacean taxa yet sequenced is also presented. The largest non-coding part (the putative mitochondrial control region) of the mitochondrial genome of Ligia oceanica is unexpectedly not AT-rich compared to the remainder of the genome. It bears two repeat regions (4× 10 bp and 3× 64 bp), and a GC-rich hairpin-like secondary structure. Some of the transfer RNAs show secondary structures which derive from the usual cloverleaf pattern. While some tRNA genes are putative targets for RNA editing, trnR could not be localized at all.
Gene order is not conserved among Peracarida, not even among isopods. The two isopod species Ligia oceanica and Idotea baltica show a similarly derived gene order, compared to the arthropod ground pattern and to the amphipod Parhyale hawaiiensis, suggesting that most of the translocation events were already present the last common ancestor of these isopods. Beyond that, the positions of three tRNA genes differ in the two isopod species. Strand bias in nucleotide frequency is reversed in both isopod species compared to other Malacostraca. This is probably due to a reversal of the replication origin, which is further supported by the fact that the hairpin structure typically found in the control region shows a reversed orientation in the isopod species, compared to other crustaceans.
The metazoan mitochondrial genome is a circular double stranded DNA molecule of about 12–20 kb length. Due to the important role of mitochondria for cell metabolism its gene content is highly conserved and typically contains the same 37 genes: 13 protein-coding genes, two ribosomal genes and 22 transfer RNA genes . In addition one A+T-rich non-coding part is present which contains essential regulatory elements for transcription and replication. It is therefore referred to as the mitochondrial control region . The organization of the mtDNA is compact with very little non-coding sequences between genes, even gene overlaps by a few nucleotides are commonly found (especially at the boundaries between nad4/nad4L and atp6/atp8). As a result the gene order of mitochondrial genomes is relatively stable because rearrangements are likely to disrupt genes. Thus changes in gene order are relatively rare, whereas tRNA genes more frequently change their position than larger protein-coding and rRNA genes . Mollusca [3–6], Brachiopoda [7–10] and Nematoda [11–14] represent phyla where a lot of rearrangements of mitochondrial genomes were reported, whereas in Chordata only few changes in gene order were found . Among arthropods a lot of species have retained the arthropod ground pattern (or a slight modification in Hexapods and Crustacea), while some taxa show frequent genome rearrangements, e.g. Myriapoda [15–17], Hymenoptera [18, 19], Acari [20, 21], and Araneae [22, 23]. Peracarid crustaceans seem to represent another example, as two partially sequenced mitochondrial genomes [24, 25] exhibit strong differences between each other and from the arthropod ground pattern.
Mitochondrial genomes offer a broad range of characters to study phylogenetic relationships of animal taxa. Besides nucleotide and amino acid sequences, tRNA secondary structures , deviations from the universal genetic code [27, 28], as well as changes in the mitochondrial gene order [29, 30] are successfully used as characters in phylogenetic inference. Especially the changes in gene order prove as extremely reliable phylogenetic characters because the probability that homoplastic translocations occur in closely related taxa is very low. Dowton et al.  calculated a chance of 1/2664 for a single gene translocation event occuring independently in two mitochondrial genomes (starting from the same gene order in both). However, this probability could be underestimated according to yet unidentified constraints on modes of gene rearrangements and should be handled with care. With about 21,000 known species peracarids comprise approximately one third of all crustacean species so far described. Within Crustacea the isopods form the largest subtaxon (10,000 species). Isopods show an amazing ecological diversity and morphological flexibility. They are common around the globe, their habitats range from deep sea plains over freshwater wells to terrestrial, even arid deserts. Although to date 30 complete mitochondrial genomes from crustaceans are available – thereof 16 from malacostracan species – a complete mitochondrial sequence of a peracarid species is still missing. Recent sequencing efforts with the amphipod Parhyale hawaiiensis  and the isopod Idotea baltica  produced almost complete genome records, lacking only the control region and some of the tRNA genes. Here we present the first complete sequence of a peracarid mitochondrial genome. Ligia oceanica (Isopoda: Oniscidea) is a terrestrial species living on rocky seashore habitats. It is found from Norway to Iceland in the north, around the british islands and the north sea coasts south to northern Spain and Portugal. We discuss changes in gene order compared to other crustacean taxa and give an overview about genome rearrangements in Crustacea. In addition we compare nucleotide composition of isopod mitochondrial genes and tRNA secondary structure, and describe in detail uncommon features of the mitochondrial control region from Ligia oceanica.
Results and discussion
The complete mitochondrial genome sequence of L. oceanica has an overall length of 15,289 bp [GenBank:DQ442914]. Successfully accomplished PCRs have proven a circular organization of the molecule (Fig. 1). Although this is the general state of metazoan mitochondrial DNA, it is mentioned here, because there is evidence for a linear organization in a related species, Armadillidum vulgare (Isopoda: Oniscidea) . All 13 protein coding subunits which are usually found in metazoan mitochondrial genomes are present, as well the two rRNA subunits (Table 1). In contrast only 21 tRNA genes instead of the typical number of 22 were identified (see below). In addition one major non-coding sequence was detected, which presumably contains the origin of replication and regulatory elements for transcription (mitochondrial control region). There are small gene overlaps at 14 gene borders. The largest has a length of 15 nucleotides (between nad2 and trnC). Some small non-coding sequences exist which occur quite often in arthropod mitochondrial genomes. The largest extends up to 52 nucleotides and is located between trnT and nad5.
The A+T content of the protein coding genes of the L. oceanica mitochondrial genome is with 60.1% (A = 28.6; C = 16.7%; G = 23.2%; T = 31.5%) at the lower end observed for malacostracan species. The values range from a 60.0% minimum given by Cherax destructor  to a 69.3% maximum by Penaeus monodon . Whereas the majority of the 13 protein coding genes show usual start codons for mtDNA, two genes begin with exceptional codons (Table 1). The gene atp8 probably starts with codon GTG. Although there is an ATA codon nine bp downstream from this start codon, alignments with atp8 amino acid sequences from other arthropods suggest the presence of more amino acids in the starting region. GTG is probably also in use as start codon in mitochondrial genes from Idotea baltica (nad1 and cox2) . The second gene with an apparently exceptional codon is cox1, which starts with ACG. Although this seems to be unusual in metazoan mitochondrial genomes, almost all other malacostracan crustaceans studied so far have this start codon for cox1. The only known exception concerns the crab Portunus trituberculatus . Two of the protein coding genes show truncated stop codons. The gene for nad1 terminates with TA whereas cox3 bears a single thymine at its end. This is a well known phenomenon in the mitochondrial genome and is frequently reported for several species. The stop codons are very likely completed by post-transcriptional polyadenylation, so that each transcript finally obtains a functional UAA terminal codon .
In most arthropods there is a strand specific bias in nucleotide frequencies [37, 38]. In detail the (+)strand contains more cytosine and adenine, while the (-)strand consequently is more rich in guanine and thymine. Some taxa show a reversal in that strand bias, among them the isopod Idotea baltica . Strand bias is best reflected in GC skew [37, 39] of mitochondrial genes (Table 2). In Ligia oceanica, as well as in Idotea baltica, GC skew is positive in (+)strand encoded genes, while it is negative in (-)strand encoded genes. This is in contrast to most other malacostracan crustaceans and is probably due to an inversion of the mitochondrial control region, or at least the replication origin . Further evidence comes from sequence analysis of the control region (see below).
The effective number of codons (ENC) is a statistic describing how far codon usage in protein-coding genes departs from the equal usage of all synonymous codons . Its range lies between 20 (when only one codon is used for each amino acid) and 62 (when all synonymous codons are equally in use). The latter departs from the usual value of 61 for nuclear genes as in invertebrate mitochondrial genomes 62 codons are in use (instead being a stop codon, UGA codes for tryptophane in the invertebrate mitochondrial code). The ENC of all published crustacean mitochondrial genomes was determined for all genes (Fig. 2), except nad4L and atp8, because these genes are too short (less than 100 codons) to get proper results. A positive correlation with G+C content in third codon positions was revealed (r2 = 0.3381; p < 0.01). There is no obvious difference seen between malacostracan and other crustaceans. Genes from the two isopod species (Ligia oceanica and Idotea baltica) are of higher G+C content and therefore show a higher than average number of effective codons. For numbers of effective codons for individual species and genes, as well as GenBank accession numbers see additional file 1.
We identified 21 out of normally found 22 transfer-RNA genes in the mitochondrial genome of Ligia oceanica. Despite extensive efforts to find secondary structures in non-coding regions the gene trnR was not found in the mitochondrial genomes sequence. By all means tRNA-Arg is essential for maintaining translation of mitochondrial gene products, so it has to be either imported into the mitochondrion, or its gene exists in the mitochondrial genome, but is subject to extensive RNA editing and therefore not identifiable by now.
Transfer-RNA genes are spread over the entire genome and are located on both strands (Fig. 1, Table 1). 14 of them were identified using tRNAscan-SE 1.21 . The other seven tRNA genes (trnD, trnC, trnE, trnI, trnF, trnS1, trnW) were found by eye inspection of otherwise non-coding regions. Some of the putative secondary structures derive from the usual cloverleaf pattern (Fig. 3): tRNA-Cys and tRNA-Ser(AGY) lack the DHU-arm. The loss of this arm in tRNA-Ser(AGY) was also observed in many other arthropod species, among malacostracan crustaceans Pseudocarcinas gigas and Macrobrachium rosenbergi , Euphausia superba , Cherax destructor , Penaeus monodon , and Portunus trituberculatus . In contrast to that, the derived structure of tRNA-Cys seems to be unique among malacostracan species studied so far. Transfer-RNA-Val and tRNA-Ile miss the TΨC-arm. Again these features are not seen in other malacostracan crustaceans.
A misplaced adenine was recorded in the anticodon loop of tRNA-Val. Its existence has been proven by repeated sequencing of different PCR-products. To assure the functionality of this gene a correctional RNA editing must be presumed in which a single nucleotide is removed. Similar post-transcriptional events with insertion and deletion of single nucleotides are known from the mitochondrial mRNAs of trypanosomes , tRNA editing was demonstrated in the centipede Lithobius forficatus . In addition several mismatches are found in tRNA stems, most of them in the acceptor stem (Fig. 3: tRNA-Gln, tRNA-Ile, tRNA-Leu1, tRNA-Leu2, tRNA-Pro, tRNA-Val, and in the anticodon stem (Fig. 3: tRNA-Ala, tRNA-Asp, tRNA-Thr, tRNA-Tyr). Such mismatches were also reported from other animal mitochondrial tRNAs and are probably further subjects to RNA editing [15, 22, 45].
Control region and repetitive sequences
There is one major non-coding region of 737 bp length located between trnI and trnE. It is assumed to be the mitochondrial control region. At its boundary to trnI it contains two sections with repetitive sequences (Fig. 4). The first consists of a series of four completely matching sequences of 10 bp each and extends into trnI. The second section is formed by a consecutive triplicate 64 bp segment. No similarities of these sequences to any other mitochondrial gene could be identified. Only a few other mitochondrial genomes were shown to contain any repeat region: rabbit mitochondrial genomes show repeated 153 bp motifs in their mitochondrial control region, varying in copy number between different individuals or tissues ; the highly aberrant mitochondrial genome of the brachiopod Lingula anatina posesses ten different unassigned repeated elements ranging in size between 28 bp and 1092 bp and in copy number between 2 and 11 . Also some insects show tandem repetitions in mitochondrial DNA .
Contrary to expectations the A+T content in the control region (55.8%) is lower than in other parts of the genome (protein coding genes: 60.1%). In contrast most other arthropods have an A+T-rich control region. While the repeat region is A+T-rich (70.3%), a 65 bp region near the 3'-end of the control region has an A+T content of only 14.1% (Fig. 4). That region is putatively folded into a hairpin-like structure with a stem consisting of 19 paired nucleotides (two mismatches) and a loop consisting of 11 nucleotides (Fig. 5). This hairpin-like structure highly resembles stem-loop structures known from insect mitochondrial control regions which have stems ranging between 15–30 bp and loops of about 9–15 nucleotides . Similar stem-loop structures were found in other crustacean species, like the mantis shrimp Squilla mantis and the spiny lobster Panulirus japonicus (Fig 5). The flanking sequences around the stem region show conserved motifs: 5'-flanking sequences show a TATA element, while 3'-flanking sequence contains a GACT in Ligia and Squilla, while the GAAAT motif typical for insects is found in Panulirus. It is assumed that these structures are of functional importance in conjunction with the origin of replication . In Ligia the flanking motifs are found in opposite direction and strand compared to that of Squilla and Panulirus. This fact gives direct evidence for an inversion of the control region in isopods (in addition to the reversed strand bias of nucleotide frequency mentioned above).
An overview about gene translocations in Crustacea
Crustacean systematics is far from being settled. While Malacostraca seems to be a well defined clade, the interrelationships between crustacean subtaxa is under debate and even monophyly of Crustacea is doubtful, with respect to the position of Hexapoda, which are probably the next relatives to a crustacean subtaxon (Malacostraca, Branchiopoda or Copepoda) [24, 34, 49–51]. Mitochondrial genome rearrangements may serve as phylogenetic markers which support sistergroup relationships among Crustacea. From 35 species of Crustcea complete or almost complete mitochondrial genome sequences are recorded in GenBank. Gene order is not conserved among these taxa: only 13 species show no changes compared to the pancrustacean ground pattern (Fig. 6 and 8: Pancrustacea ground pattern). Transfer-RNA genes are more often translocated than other genes, probably because of their small size.
Crustacea and Hexapoda (united as Pancrustacea [30, 52, 53] or Tetraconata ) share the same ground pattern in mitochondrial gene order . It differs from the euarthropod ground pattern  by the position of one tRNA gene: trnL2 is located between cox1 and cox2, whereas in Chelicerata, Myriapoda and Onychophora trnL2 is located between nad1 and rrnl, adjacent to trnL1 [29, 30]. Among other data, mitochondrial gene translocations have shown that the enigmatic Remipedia and Pentastomida definitely belong to Pancrustacea, as they show the above mentioned translocation of trnL2 . Only three crustacean species do not show this character: the cephalocarid Hutchinsoniella macrantha , where trnL2 probably is secondarily translocated to another position (Fig. 6, No. 7 from Hutchinsoniella), and the two copepod species Tigriopus japonicus  and Lepeophtheirus salmonis , which underwent a complete shuffling of the mitochondrial genome.
Three species (belonging to Cephalocarida, Branchiura and Pentatomida) share a translocation of trnK to a position between trnR and trnN . Among these, the tongue worm Armillifer armillatus and the fish louse Argulus americanus share one further translocation (trnQ), together with mtDNA sequence analysis supporting a close relationship between Pentastomida and Branchiura [51, 56]. That was already discussed according to sperm morphology [59, 60] and 18S molecular sequence data . trnK is also translocated in all other taxa refered to as members of "Maxillopoda", a systematic unit only weakly based on morphological characters. However, Ostracoda, Copepoda and Cirripedia each show different positions for trnK compared to the above mentioned taxa, so there is no good reason to take this as a homology. Because of contrary results from morphological and sequence based analyses [51, 62] it is also questionable to unite the Cephalocarida with Branchiura and Pentastomida to one clade, solely based on the common translocation of trnK (No. 1 for Hutchinsoniella /Argulus/Armillifer in Fig. 6 and Fig. 7).
The three species of Cirripedia , [GenBank:NC_006293; GenBank:NC_008974] share several translocations of tRNA genes (trnA, trnE, trnP). Another series of events is difficult to reconstruct: trnC and trnY are translocated in all three species to a position between trnS2 and nad1, but in different order and on different strands (referred to as No. 6 from Cirripedia in Fig. 6 and 7). In addition one species shows a triplication of trnC (No. 7 from Pollicipes polymerus in Fig. 6 and 7). Two further differences are reported in Megabalanus volcano: an inversion of a block of five genes (No. 7) and the probable translocation of trnK and trnQ to a position between trnY and trnC (No. 8 in Fig. 6 and 7). An alternative explanation is that trnK/trnQ were primarily translocated to that position seen in Megabalanus and secondarily translocated to the position seen in Pollicipes and Tetraclita. With data from these three species alone, it is not possible to reconstruct a ground pattern of mitochondrial gene order of Cirripedia with respect to the position of trnK, trnQ, trnY and trnC.
A lot of further translocation events are recorded only in single species, making them useless in phylogenetic analysis of the actual data set. Large genome rearrangements involving also protein-coding genes are seen in the branchiuran Argulus americanus , and in the ostracod Vargula hilgendorfi , both accompanied by a duplication of the control region. But the highest degree of genome rearrangement was found in the two copepod species [57, 58], where a complete reshuffling of the mitochondrial genomes has led to a gene order with almost no similarities between the two species and to other crustaceans. Partial genomes from two other copepod species revealed even more rearrangements . In contrast three species of Branchiopoda have retained the pancrustacean ground pattern. The fourth species, Artemia franciscana shows two tRNA gene translocations (trnI, trnQ) .
Among Malacostraca 10 from 19 species have retained the pancrustacean ground pattern: six mantis shrimps (Stomatopoda) and four members of Decapoda (Fig. 8 and 9). Among Decapoda independent translocation events changed gene order in Pagurus longicarpus , Cherax destructor , and Brachyura. The four species of Brachyura share a translocation of trnH, the freshwater crabs Geothelphusa dehaani  and Eriocheir japonica  show further translocations. In the Euphausiacea Euphausia superba a swap between trnL1 and trnL2 seems to have happened, probably preceded by a gene duplication .
Besides the two isopod species (Ligia oceanica, Idotea baltica), only one other peracarid mitochondrial sequence, from the amphipod Parhyale hawaiiensis, was published before . From six gene translocations that must be assumed to get the gene order of Parhyale hawaiiensis, none is shared with the mitochondrial genomes of isopods (Fig. 8 and 9). Only trnI is translocated in both taxa, but as its new position is different in Parhyale hawaiiensis and Ligia oceanica (and not known in Idotea baltica), there is no reason to presume that a translocation had already happened in their common ancestor. Instead we assume an independent translocation of this gene in isopods and amphipods. This implies, that the ground pattern of gene order in Peracarida must be identical to that of Pancrustacea.
Gene translocations in isopods
A comparison of the complete mitochondrial genome of L. oceanica with the ancestral state of pancrustacea  demonstrates several changes in gene order (Fig. 8). All in all 11 genes (cob, nad1, nad5, rrnS, trnI, trnL1, trnF, trnS1, trnT, trnW, trnV) and the control region are found in other relative positions than reported in other malacostracan crustaceans. By reason of parsimony we assume that these positional changes were due to nine gene translocations (Fig. 8, 9; No. 1–9) and a translocation of the major non-coding region (NCR). The genes nad5 and trnF, as well as nad1 and trnL1 retained their adjacent positions, so that they probably were translocated as a block. The other genes were most likely repositioned by single translocation events. Five translocations led also to inversion of genes to the complementary strand: cob, trnT and trnW changed from the (+)strand to the (-)strand, whereas rrnS and the block of nad 5 and trnF were inverted from (-)strand to (+)strand.
From the 11 genes being translocated in Ligia oceanica, seven are also found in the same new positions in the mitochondrial genome of Idotea baltica (Fig. 8) . In addition the mitochondrial control region, trnW, trnI and trnS1 are not found in their original position in Idotea baltica, but will probably be found in the region not sequenced yet – between rrnS and cob, similar to Ligia oceanica. Of all genes translocated in Ligia oceanica only trnV is found in its original position in Idotea baltica. So of nine gene translocation events supposed for Ligia oceanica, eight must already have happended in the common ancestor of both species (Fig. 8, 9; No. 1–8). The derived gene order is probably the ground pattern for an isopod subtaxon comprised of Oniscidea and Valvifera (and probably more taxa). Translocation of trnV probably happened after the split of the oniscidean and valviferan lineages (Fig. 9, No. 9). In contrast, translocation of trnN probably happened in the lineage leading to Idotea baltica (Fig. 9, No.11). It is located in original position in Ligia oceanica, but missing in that position in Idotea baltica. We do not know about the fate of trnR in Ligia oceanica, but in Idotea baltica it was also subject to a translocation (Fig. 9, No. 10).
It is noticeable that in Ligia oceanica all translocated genes found their new position in a segment comprising about one third of the complete genome (between trnA and trnH). This area bears a cluster of tRNA genes in the ancestral gene order of arthropods. It seems to be a "hot spot" of genome rearrangements in arthropods [19, 20].
The first complete mitochondrial genome sequence of a peracarid arthropod, the isopod Ligia oceanica, shows the usual compact and circular organization known from other Metazoa. Gene order is not conserved among peracarids and even not among isopods. In Ligia oceanica 11 genes plus the control region have changed their relative positions in comparison to the pancrustacean ground pattern, implying to be the result of nine gene translocation events. No gene translocation is shared with the amphipod Parhyale hawaiiensis, whereas eight gene translocations were probably already present in the common ancestor of Ligia oceanica and another isopod, Idotea baltica. Both isopod mitochondrial genomes differ by the position of three tRNA genes (trnR, trnV, trnN) and share an inverted strand bias of nucleotide frequencies compared to other malacostracan crustaceans. Reason for this is probably an inversion of the replication origin. This is confirmed by the fact that the typical hairpin-like secondary structure commonly found in mitochondrial control regions is found in opposite orientation compared to other crustacean species.
A broad survey of mitochondrial gene rearrangements in Crustacea reveals a great variation of gene order. Characters derived from gene order (= gene translocations, inversions or duplications) do not solve overall phylogenetic relationships between major crustacean subtaxa. However, they will probably be helpful in analyses of internal phylogeny of some of these subtaxa (Cirripedia, Brachyura, Peracarida), when more data will be provided.
Sample and DNA extraction
A specimen of Ligia oceanica originally collected at the coast of the North Sea island Helgoland (Germany) and preserved in 99% ethanol was utilized for the DNA extraction process. 2–3 pleopods were applied to the DNeasy Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer's protocol to receive the total genomic DNA.
The first partial mitogenomic sequences were obtained by using two insect based primer pairs N4 + 16S2 and CytB + N4(87) . Additional intragenetic parts of cox1, cox3, nd4, nad5, rrnL and rrnS were determined by using six crustacean primer pairs [25, 70]. The remaining gaps were closed by making use of specific primer pairs designed according to the sequences of the aforementioned genes (for primer sequences see additional file 2). Larger PCR products were sequenced by primer walking strategy. To abbreviate this longsome process two longer PCR products were sequenced with the primer pairs S1–S19 from a set of primers which was successfully applied to decapod crustaceans . All primers were purchased from Metabion (München, Germany).
PCR and purification of PCR products
The PCRs were performed with an Eppendorf Mastercycler or Eppendorf Mastercycler Gradient. The cycling was set up with an initial denaturation step at 94°C for 2 minutes, followed by 40 cycles comprising denaturation at 94°C for 30 seconds, annealing at 45–52°C (primer specific) for 1 minute and elongation at 72° for 2–5 minutes depending on the expected length of the PCR product. The process was completed with a final elongation at 68°C for 2 minutes. The reaction volume amounted 50 μl containing 1 μl dNTP mix (Eppendorf), 0.25 μl HotMasterTaq DNA polymerase (5 U/μl; Eppendorf), 5 μl HotMasterTaq buffer (Eppendorf), 1 μl primer mix (10 μM each), 1 μl DNA template and 41.75 μl sterilized distilled water (Eppendorf). The PCR products were separated with a 1% TBE agarose gel, stained with ethidium bromide and inspected subsequently under UV transillumination.
For the purification of the PCR Products the QIAquick PCR Purification Kit (Qiagen) as well as the Blue Matrix DNA Purification Kit (Eurx) were used. Abiding to the manufacturers protocols both kits produced equivalent yields. All PCR products were stored at -20°C until sequencing was performed.
Cloning and transformation
In one case a single PCR fragment, due to low sequence signal quality, ranging from rrnS to cob, which contained the major non-coding region, was sequenced after cloning with the pGEM-T Easy Vector System (Promega). We followed the manufacturers protocol with the exception that half volumes (5 μl) were used for the cloning reaction. For the transformation Escherichia coli XL-gold (Stratagene) were used. Colonies that contained recombinant plasmids were identified with selection plates (LB/ampicillin/IPTG/X-Gal). To verify insertion of the PCR product, few cells were applied to a Colony PCR using the vector primers M13F and M13R. The reaction volume amounted 20 μl with aforementioned PCR ingredients added proportionally. The colony-cycling consisted of an initial denaturation step (10 minutes, 95°C) followed by 25 cycles of denaturation (30 seconds, 94°C), annealing (30 seconds, 46°C) and elongation (4 minutes, 68°C). The colony PCR was closed by a final elongation (3 minutes, 68°C). The result was inspected with an agarose gel under UV transillumination. Positive tested colonies were proliferated in a LB/ampicillin Medium. Subsequently the plasmids were extracted with the Quantum Prep Kit (Bio Rad) and finally stored at -20°C.
Sequencing and sequence analysis
Cycle sequencing reactions were performed with the CEQ DTCS Quick Start Kit (Beckman Coulter) following the manufacturers protocols. The same primers and thermocyclers were used as in PCRs. The temperature profile included 30 cycles comprising denaturation at 94°C for 20 seconds, annealing at 45–52°C (primer specific) for 20 seconds and elongation at 60°C for 4 minutes. Plasmids were preheated additionally before the sequencing reaction (96°C for 1 minute). The separation was executed by a CEQ 8000 cappillary sequencer (Beckman Coulter) and analyzed with the appendant CEQ software (software version: 5.0.360, instrument version: 6.0.2).
Gene annotation and sequence analysis
The alignment of the fragments to complete the whole mitochondrial DNA sequence was done in BioEdit 184.108.40.206 . Each partial sequence was ascertained twice at least to prevent sequencing faults. Ambiguous base pairs were validated manually referring to chromatograms. Gene identification was determined by BLAST search on GenBank databases  and by comparison to the mitochondrial Genome of Drosophila yakuba (NC001322). Boundaries of the protein coding genes were determined with a multiple alignment of other crustacean amino acid sequences. It was assumed that they were specified by the first start and stop codons in frame. Transfer RNA genes were determined with tRNAscan-SE 1.21  or by eye inspection for anti-codon sequences and secondary structures in regions between identified genes. Hairpin structures in non-coding regions were also identified by eye inspection. The control region and RNA genes were assumed to extend to adjacent genes, due to the lack of resources for a better determination of their boundaries. Nucleotide frequencies of protein coding and RNA genes were calculated with the DAMBE software package , effective number of codons was determined according to  with INCA 1.20 . The complete genome sequence is submitted to NCBI GenBank [GenBank:DQ442914].
- atp6 and 8:
genes encoding ATPase subunit 6 and 8
- cox1-3 :
genes encoding cytochrome oxidase subnunits I-III
- cob :
gene encoding cytochrome b
effective number of codons
- nad1-6 and nad4L:
genes encoding NADH dehydrogenase subunits 1–6 and 4L
polymerase chain reaction
- rrnl :
large rRNA subunit (gene)
- rrnS :
small rRNA subunit (gene)
- trnX (where X is replaced by one letter amino acid code):
Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992, 141: 173-216.
Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.
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.
Boore JL, Medina M, Rosenberg LA: Complete sequences of the highly rearranged molluscan mitochondrial genomes of the Scaphopod Graptacme eborea and the bivalve Mytilus edulis. Mol Biol Evol. 2004, 21: 1492-1503. 10.1093/molbev/msh090.
Dreyer H, Steiner G: The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca). Mol Phylogenet Evol. 2004, 31: 605-617. 10.1016/j.ympev.2003.08.007.
Knudsen B, Kohn AB, Nahir B, McFadden CS, Moroz LL: Complete DNA sequence of the mitochondrial genome of the sea-slug, Aplysia californica: conservation of the gene order in Euthyneura. Mol Phylogenet Evol. 2006, 38: 459-469. 10.1016/j.ympev.2005.08.017.
Endo K, Noguchi Y, Ueshima R, Jacobs HT: Novel repetitive structures, deviant protein-encoding sequences and unidentified ORFs in the mitochondrial genome of the brachiopod Lingula anatina. J Mol Evol. 2005, 61: 36-53. 10.1007/s00239-004-0214-5.
Helfenbein KG, Brown WM, Boore JL: The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa. Mol Biol Evol. 2001, 18: 1734-1744.
Noguchi Y, Endo K, Tajima F, Ueshima R: The mitochondrial genome of the brachiopod Laqueus rubellus. Genetics. 2000, 155: 245-259.
Stechmann A, Schlegel M: Analysis of the complete mitochondrial DNA sequence of the brachiopod Terebratulina retusa places Brachiopoda within the protostomes. P Roy Soc B Biol Sci. 1999, 266: 2043-2052. 10.1098/rspb.1999.0885.
Lavrov DV, Brown WM: 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. Genetics. 2001, 157: 621-637.
He Y, Jones J, Armstrong M, Lamberti F, Moens M: The mitochondrial genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): considerable economization in the length and structural features of encoded genes. J Mol Evol. 2005, 61: 819-833. 10.1007/s00239-005-0102-7.
Hu M, Chilton NB, Gasser RB: The mitochondrial genome of Strongyloides stercoralis (Nematoda) – idiosyncratic gene order and evolutionary implications. Int J Parasitol. 2003, 33: 1393-1408. 10.1016/S0020-7519(03)00130-9.
Keddie EM, Higazi T, Unnasch TR: The mitochondrial genome of Onchocerca volvulus: sequence, structure and phylogenetic analysis. Mol Biochem Parasitol. 1998, 95: 111-127. 10.1016/S0166-6851(98)00102-9.
Lavrov DV, 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: 13738-13742. 10.1073/pnas.250402997.
Lavrov DV, Boore JL, Brown WM: Complete mtDNA sequences of two millipedes suggest a new model for mitochondrial gene rearrangements: duplication and nonrandom loss. Mol Biol Evol. 2002, 19: 163-169.
Negrisolo E, Minelli A, Valle G: Extensive gene order rearrangement in the mitochondrial genome of the centipede Scutigera coleoptrata. J Mol Evol. 2004, 58: 413-423. 10.1007/s00239-003-2563-x.
Dowton M, Austin AD: Evolutionary dynamics of a mitochondrial rearrangement "hot spot" in the Hymenoptera. Mol Biol Evol. 1999, 16: 298-309.
Dowton M, Castro LR, Campbell SL, Bargon SD, Austin AD: Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J Mol Evol. 2003, 56: 517-526. 10.1007/s00239-002-2420-3.
Black WC, Roehrdanz RL: Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes. Mol Biol Evol. 1998, 15: 1772-1785.
Shao R, Mitani H, Barker SC, Takahashi M, Fukunaga M: Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum. J Mol Evol. 2005, 60: 764-773. 10.1007/s00239-004-0226-1.
Masta SE, Boore JL: The complete mitochondrial genome sequence of the spider Habronattus oregonensis reveals rearranged and extremely truncated tRNAs. Mol Biol Evol. 2004, 21: 893-902. 10.1093/molbev/msh096.
Qiu Y, Song D, Zhou K, Sun H: The mitochondrial sequences of Heptathela hangzhouensis and Ornithoctonus huwena reveal unique gene arrangements and atypical tRNAs. J Mol Evol. 2005, 60: 57-71. 10.1007/s00239-004-0010-2.
Cook CE, Yue Q, Akam M: Mitochondrial genomes suggest that hexapods and crustaceans are mutually paraphyletic. P Roy Soc B Biol Sci. 2005, 272: 1295-1304. 10.1098/rspb.2004.3042.
Podsiadlowski L, Bartolomaeus T: Major rearrangements characterize the mitochondrial genome of the isopod Idotea baltica (Crustacea: Peracarida). Mol Phylogenet Evol. 2006, 40: 893-899. 10.1016/j.ympev.2006.04.008.
Macey JR, Schulte JA, Larson A: Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards. Syst Biol. 2000, 49: 257-277. 10.1080/10635159950173843.
Castresana J, Feldmaier-Fuchs G, Paabo S: Codon reassignment and amino acid composition in hemichordate mitochondria. Proc Natl Acad Sci USA. 1998, 95: 3703-3707. 10.1073/pnas.95.7.3703.
Telford MJ, Herniou EA, Russell RB, Littlewood DT: Changes in mitochondrial genetic codes as phylogenetic characters: two examples from the flatworms. Proc Natl Acad Sci USA. 2000, 97: 11359-11364. 10.1073/pnas.97.21.11359.
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.
Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392: 667-668. 10.1038/33577.
Dowton M, Castro LR, Austin AD: Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: The examination of genome 'morphology'. Invertebr Syst. 2002, 16: 345-356. 10.1071/IS02003.
Raimond R, Marcade I, Bouchon D, Rigaud T, Bossy JP, Souty-Grosset C: Organization of the large mitochondrial genome in the isopod Armadillidium vulgare. Genetics. 1999, 151: 203-210.
Miller AD, Nguyen TTT, Burridge CP, Austin CM: Complete mitochondrial DNA sequence of the Australian freshwater crayfish, Cherax destructor (Crustacea: Decapoda: Parastacidae): a novel gene order revealed. Gene. 2004, 331: 65-72. 10.1016/j.gene.2004.01.022.
Wilson K, Cahill V, Ballment E, Benzie J: The complete sequence of the mitochondrial genome of the crustacean Penaeus monodon: Are malacostracan crustaceans more closely related to insects than to branchiopods?. Mol Biol Evol. 2000, 17: 863-874.
Yamauchi MM, Miya MU, Nishida M: Complete mitochondrial DNA sequence of the swimming crab, Portunus trituberculatus (Crustacea: Decapoda: Brachyura). Gene. 2003, 311: 129-135. 10.1016/S0378-1119(03)00582-1.
Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290: 470-474. 10.1038/290470a0.
Hassanin A, Leger 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.
Hassanin A: Phylogeny of Arthropoda inferred from mitochondrial sequences: strategies for limiting the misleading effects of multiple changes in pattern and rates of substitution. Mol Phylogenet Evol. 2006, 38: 100-116. 10.1016/j.ympev.2005.09.012.
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.
Wright F: The Effective Number of Codons Used in A Gene. Gene. 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.955.
Miller AD, Murphy NP, Burridge CP, Austin CM: Complete mitochondrial DNA sequences of the decapod crustaceans Pseudocarcinus gigas (Menippidae) and Macrobrachium rosenbergii (Palaemonidae). Mar Biotechnol. 2005, 7: 339-349. 10.1007/s10126-004-4077-8.
Machida RJ, Miya MU, Yamauchi MM, Nishida M, Nishida S: Organization of the mitochondrial genome of Antarctic krill Euphausia superba (Crustacea: Malacostraca). Mar Biotechnol. 2004, 6: 238-250. 10.1007/s10126-003-0016-6.
Benne R: RNA editing in trypanosomes. Eur J Biochem. 1994, 221: 9-23. 10.1111/j.1432-1033.1994.tb18710.x.
Yokobori S, Paabo S: Transfer RNA editing in land snail mitochondria. Proc Natl Acad Sci USA. 1995, 92: 10432-10435. 10.1073/pnas.92.22.10432.
Casane D, Dennebouy N, de RH, Mounolou JC, Monnerot M: Genetic analysis of systematic mitochondrial heteroplasmy in rabbits. Genetics. 1994, 138: 471-480.
Zhang DX, Hewitt GM: Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997, 25: 99-120. 10.1016/S0305-1978(96)00042-7.
Zhang DX, Szymura JM, Hewitt GM: Evolution and structural conservation of the control region of insect mitochondrial DNA. J Mol Evol. 1995, 40: 382-391. 10.1007/BF00164024.
Mallatt J, Giribet G: Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch. Mol Phylogenet Evol. 2006, 40: 772-794. 10.1016/j.ympev.2006.04.021.
Nardi F, Spinsanti G, Boore JL, Carapelli A, Dallai R, Frati F: Hexapod origins: monophyletic or polyphyletic?. Science. 2003, 299: 1887-1889. 10.1126/science.1078607.
Regier JC, Shultz JW, Kambic RE: Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc R Soc Lond B Biol Sci. 2005, 272: 395-401. 10.1098/rspb.2004.2917.
Friedrich M, Tautz D: Ribosomal DNA phylogeny of the major extant arthropod classes and the evolution of myriapods. Nature. 1995, 376: 165-167. 10.1038/376165a0.
Shultz JW, Regier JC: Phylogenetic analysis of arthropods using two nuclear protein-encoding genes supports a crustacean + hexapod clade. P Roy Soc B Biol Sci. 2000, 267: 1011-1019. 10.1098/rspb.2000.1104.
Richter S: The Tetraconata concept: Hexapod-crustacean relationships and the phylogeny of Crustacea. Org Divers Evol. 2002, 2: 217-237. 10.1078/1439-6092-00048.
Staton JL, Daehler LL, Brown WM: Mitochondrial gene arrangement of the horseshoe crab Limulus polyphemus L.: conservation of major features among arthropod classes. Mol Biol Evol. 1997, 14: 867-874.
Lavrov DV, Brown WM, Boore JL: Phylogenetic position of the Pentastomida and (pan)crustacean relationships. Proc R Soc Lond B Biol Sci. 2004, 271: 537-544. 10.1098/rspb.2003.2631.
Machida RJ, Miya MU, Nishida M, Nishida S: Complete mitochondrial DNA sequence of Tigriopus japonicus (Crustacea: Copepoda). Mar Biotechnol. 2002, 4: 406-417. 10.1007/s10126-002-0033-x.
Tjensvoll K, Hodneland K, Nilsen F, Nylund A: Genetic characterization of the mitochondrial DNA from Lepeophtheirus salmonis (Crustacea: Copepoda). A new gene organization revealed. Gene. 2005, 353: 218-230. 10.1016/j.gene.2005.04.033.
Storch V, Jamieson BGM: Further spermatological evidence for including the Pentastomida (Tongue Worms) in the Crustacea. Int J Parasitol. 1992, 22: 95-108. 10.1016/0020-7519(92)90085-Y.
Wingstrand KG: Comparative spermatology of a pentastomid, Raillietella hemidactyli, and a branchiuran crustacean, Argulus foliaceus, with a discussion of pentastomid relationships. Kong Danske Vid Selsk Biol Skri. 1972, 19: 1-72.
Abele LG, Kim W, Felgenhauer BE: Molecular evidence for inclusion of the phylum Pentastomida in the Crustacea. Mol Biol Evol. 1989, 6: 685-691.
Giribet G, Richter S, Edgecombe GD, Wheeler WC: The position of crustaceans within Arthropoda – evidence from nine molecular loci and morphology. Crustacean Issues. 2005, 16: 307-330.
Ogoh K, Ohmiya Y: Complete mitochondrial DNA sequence of the sea-firefly, Vargula hilgendorfii (Crustacea, Ostracoda) with duplicate control regions. Gene. 2004, 327: 131-139. 10.1016/j.gene.2003.11.011.
Machida RJ, Miya MU, Nishida M, Nishida S: Large-scale gene rearrangements in the mitochondrial genomes of two calanoid copepods Eucalanus bungii and Neocalanus cristatus (Crustacea), with notes on new versatile primers for the srRNA and COI genes. Gene. 2004, 332: 71-78. 10.1016/j.gene.2004.01.019.
Valverde JR, Batuecas B, Moratilla C, Marco R, Garesse R: The complete mitochondrial DNA sequence of the crustacean Artemia franciscana. J Mol Evol. 1994, 39: 400-408. 10.1007/BF00160272.
Hickerson MJ, Cunningham CW: Dramatic mitochondrial gene rearrangements in the hermit crab Pagurus longicarpus (Crustacea, Anomura). Mol Biol Evol. 2000, 17: 639-644.
Segawa RD, Aotsuka T: The mitochondrial genome of the Japanese freshwater crab, Geothelphusa dehaani (Crustacea: Brachyura): Evidence for its evolution via gene duplication. Gene. 2005, 355: 28-39. 10.1016/j.gene.2005.05.020.
Sun H, Zhou K, Song D: Mitochondrial genome of the Chinese mitten crab Eriocheir japonica sinenesis (Brachyura: Thoracotremata: Grapsoidea) reveals a novel gene order and two target regions of gene rearrangements. Gene. 2005, 349: 207-217. 10.1016/j.gene.2004.12.036.
Roehrdanz RL, Degrugillier ME, Black WCI: Novel rearrangements of arthropod mitochondrial DNA detected with long-PCR: Applications to arthropod phylogeny and evolution. Mol Biol Evol. 2002, 19: 841-849.
Podsiadlowski L, Bartolomaeus T: Organization of the mitochondrial genome of mantis shrimp Pseudosquilla ciliata (Crustacea: Stomatopoda). Mar Biotechnol. 2005, 7: 618-624. 10.1007/s10126-005-0017-8.
Yamauchi MM, Miya MU, Machida RJ, Nishida M: PCR-based approach for sequencing mitochondrial genomes of decapod crustaceans, with a practical example from kuruma prawn (Marsupenaeus japonicus). Mar Biotechnol. 2004, 6: 419-429. 10.1007/s10126-003-0036-2.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov]
Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.
Supek F, Vlahovicek K: INCA: synonymous codon usage analysis and clustering by means of self-organizing map. Bioinformatics. 2004, 20: 2329-2330. 10.1093/bioinformatics/bth238.
Monforte A, Barrio E, Latorre A: Characterization of the length polymorphism in the A + T-rich region of the Drosophila obscura group species. J Mol Evol. 1993, 36: 214-223. 10.1007/BF00160476.
Richter S, Scholtz G: Phylogenetic analysis of the Malacostraca (Crustacea). J Zool Syst Evol Res. 2001, 39: 113-136. 10.1046/j.1439-0469.2001.00164.x.
The authors thank Prof. Thomas Bartolomaeus for his kind support during our work in his group. We also thank Christine Bergmann (FU Berlin), Dipl.biol. Achim Meyer (Universität Mainz), and Dr. Christoph Bleidorn (Universität Potsdam) for technical support during this study, and Prof. Heinz-Dieter Franke (Biologische Anstalt Helgoland) for providing us with isopods. We also thank two anonymous reviewers for valuable suggestions to improve the manuscript. This study was in part supported by grants from the German science foundation (DFG): Ba 1520/10-1 (to LP) and Po 765/2-1 (to LP).
FK did the majority of the laboratory work and the primary sequence analysis, LP was the initiator and supervisor of this work. Final analyses of data, discussion of results and drawing of the manuscript was done by both authors in equal shares.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kilpert, F., Podsiadlowski, L. The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features. BMC Genomics 7, 241 (2006). https://doi.org/10.1186/1471-2164-7-241
- Mitochondrial Genome
- Mitochondrial Control Region
- Strand Bias
- Ground Pattern
- Gene Translocation