Organization of the large mitochondrial genome in the isopod Armadillidium vulgare

Rapid evolution, rearrangements and whole mitogenome duplication in the Australian stingless bees Tetragonula (Hymenoptera: Apidae): A steppingstone towards understanding mitochondrial function and evolution

The extreme conservation of mitochondrial genomes in metazoans poses a significant challenge to understanding mitogenome evolution. However, the presence of variation in gene order or genome structure, found in a small number of taxa, can provide unique insights into this evolution. Previous work on two stingless bees in the genus Tetragonula (T. carbonaria and T. hockingsi) revealed highly divergent CO1 regions between them and when compared to the bees from the same tribe (Meliponini), indicating rapid evolution. Using mtDNA isolation and Illumina sequencing, we elucidated the mitogenomes of both species. In both species, there has been a duplication of the whole mitogenome to give a total genome size of 30,666 bp in T. carbonaria; and 30,662 bp in T. hockingsi. These duplicated genomes present a circular structure with two identical and mirrored copies of all 13 protein coding genes and 22 tRNAs, with the exception of a few tRNAs that are present as single copies. In addition, the mitogenomes are characterized by rearrangements of two block of genes. We believe that rapid evolution is present in the whole Indo-Malay/Australasian group of Meliponini but is extraordinarily elevated in T. carbonaria and T. hockingsi, probably due to founder effect, low effective population size and the mitogenome duplication. All these features - rapid evolution, rearrangements, and duplication - deviate significantly from the vast majority of the mitogenomes described so far, making the mitogenomes of Tetragonula unique opportunities to address fundamental questions of mitogenome function and evolution.

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.

The first complete mitochondrial genome of a parasitic isopod supports Epicaridea Latreille, 1825 as a suborder and reveals the less conservative genome of isopods

The complete mitochondrial genome sequence of the holoparasitic isopod Gyge ovalis (Shiino, 1939) has been determined. The mitogenome is 14,268 bp in length and contains 34 genes: 13 protein-coding genes, two ribosomal RNA, 19 tRNA and a control region. Three tRNA genes (trnE, trnI and trnS1) are missing. Most of the tRNA genes show secondary structures which derive from the usual cloverleaf pattern except for trnC which is characterised by the loss of the DHU-arm. Compared to the isopod ground pattern and Eurydice pulchra Leach, 1815 (suborder Cymothoida Wägele, 1989), the genome of G. ovalis shows few differences, with changes only around the control region. However, the genome of G. ovalis is very different from that of non-cymothoidan isopods and reveals that the gene order evolution in isopods is less conservative compared to other crustaceans. Phylogenic trees were constructed using maxiumum likelihood and Bayesian inference analyses based on 13 protein-coding genes. The results do not support the placement of G. ovalis with E. pulchra and Bathynomus sp. in the same suborder; rather, G. ovalis appears to have a closer relationship to Ligia oceanica (Linnaeus, 1767), but this result suggests a need for more data and further analysis. Nevertheless, these results cast doubt that Epicaridea Latreille, 1825 can be placed as an infraorder within the suborder Cymothoida, and Epicaridea appears to also deserve subordinal rank. Further development of robust phylogenetic relationships across Isopoda Latreille, 1817 will require more genetic data from a greater diversity of taxa belonging to all isopod suborders.

Untangling Heteroplasmy, Structure, and Evolution of an Atypical Mitochondrial Genome by PacBio Sequencing

The highly compact mitochondrial (mt) genome of terrestrial isopods (Oniscidae) presents two unusual features. First, several loci can individually encode two tRNAs, thanks to single nucleotide polymorphisms at anticodon sites. Within-individual variation (heteroplasmy) at these loci is thought to have been maintained for millions of years because individuals that do not carry all tRNA genes die, resulting in strong balancing selection. Second, the oniscid mtDNA genome comes in two conformations: a ∼14 kb linear monomer and a ∼28 kb circular dimer comprising two monomer units fused in palindrome. We hypothesized that heteroplasmy actually results from two genome units of the same dimeric molecule carrying different tRNA genes at mirrored loci. This hypothesis, however, contradicts the earlier proposition that dimeric molecules result from the replication of linear monomers-a process that should yield totally identical genome units within a dimer. To solve this contradiction, we used the SMRT (PacBio) technology to sequence mirrored tRNA loci in single dimeric molecules. We show that dimers do present different tRNA genes at mirrored loci; thus covalent linkage, rather than balancing selection, maintains vital variation at anticodons. We also leveraged unique features of the SMRT technology to detect linear monomers closed by hairpins and carrying noncomplementary bases at anticodons. These molecules contain the necessary information to encode two tRNAs at the same locus, and suggest new mechanisms of transition between linear and circular mtDNA. Overall, our analyses clarify the evolution of an atypical mt genome where dimerization counterintuitively enabled further mtDNA compaction.

The highly rearranged mitochondrial genomes of the crabs Maja crispata and Maja squinado (Majidae) and gene order evolution in Brachyura

We sequenced the mitochondrial genomes of the spider crabs Maja crispata and Maja squinado (Majidae, Brachyura). Both genomes contain the whole set of 37 genes characteristic of Bilaterian genomes, encoded on both α- and β-strands. Both species exhibit the same gene order, which is unique among known animal genomes. In particular, all the genes located on the β-strand form a single block. This gene order was analysed together with the other nine gene orders known for the Brachyura. Our study confirms that the most widespread gene order (BraGO) represents the plesiomorphic condition for Brachyura and was established at the onset of this clade. All other gene orders are the result of transformational pathways originating from BraGO. The different gene orders exhibit variable levels of genes rearrangements, which involve only tRNAs or all types of genes. Local homoplastic arrangements were identified, while complete gene orders remain unique and represent signatures that can have a diagnostic value. Brachyura appear to be a hot-spot of gene order diversity within the phylum Arthropoda. Our analysis, allowed to track, for the first time, the fully evolutionary pathways producing the Brachyuran gene orders. This goal was achieved by coupling sophisticated bioinformatic tools with phylogenetic analysis.

Multiple Conserved Heteroplasmic Sites in tRNA Genes in the Mitochondrial Genomes of Terrestrial Isopods (Oniscidea)

Mitochondrial genome structure and organization are relatively conserved among metazoans. However, in many isopods, especially the terrestrial isopods (Oniscidea), the mitochondrial genome consists of both ∼14-kb linear monomers and ∼28-kb circular dimers. This unusual organization is associated with an ancient and conserved constitutive heteroplasmic site. This heteroplasmy affects the anticodon of a tRNA gene, allowing this single locus to function as a “dual” tRNA gene for two different amino acids. Here, we further explore the evolution of these unusual mitochondrial genomes by assembling complete mitochondrial sequences for two additional Oniscidean species, Trachelipus rathkei and Cylisticus convexus. Strikingly, we find evidence of two additional heteroplasmic sites that also alter tRNA anticodons, creating additional dual tRNA genes, and that are conserved across both species. These results suggest that the unique linear/circular organization of isopods’ mitochondrial genomes may facilitate the evolution of stable mitochondrial heteroplasmies, and, conversely, once such heteroplasmies have evolved, they constrain the multimeric structure of the mitochondrial genome in these species. Finally, we outline some possible future research directions to identify the factors influencing mitochondrial genome evolution in this group.

If the cap fits, wear it: an overview of telomeric structures over evolution

Genome organization into linear chromosomes likely represents an important evolutionary innovation that has permitted the development of the sexual life cycle; this process has consequently advanced nuclear expansion and increased complexity of eukaryotic genomes. Chromosome linearity, however, poses a major challenge to the internal cellular machinery. The need to efficiently recognize and repair DNA double-strand breaks that occur as a consequence of DNA damage presents a constant threat to native chromosome ends known as telomeres. In this review, we present a comparative survey of various solutions to the end protection problem, maintaining an emphasis on DNA structure. This begins with telomeric structures derived from a subset of prokaryotes, mitochondria, and viruses, and will progress into the typical telomere structure exhibited by higher organisms containing TTAGG-like tandem sequences. We next examine non-canonical telomeres from Drosophila melanogaster, which comprise arrays of retrotransposons. Finally, we discuss telomeric structures in evolution and possible switches between canonical and non-canonical solutions to chromosome end protection.

Inverted repeats and genome architecture conversions of terrestrial isopods mitochondrial DNA

Mitochondrial DNA (mtDNA) is usually depicted as a circular molecule, however, there is increasing evidence that linearization of mtDNA evolved independently many times in organisms such as fungi, unicellular eukaryotes, and animals. Recent observations in various models with linear mtDNA revealed the presence of conserved inverted repeats (IR) at both ends that, when they become single-stranded, may be able to fold on themselves to create telomeric-hairpins involved in genome architecture conversions. The atypical mtDNA of terrestrial isopods (Crustacea: Oniscidea) composed of linear monomers and circular dimers is an interesting model to study genome architecture conversions. Here, we present the mtDNA control region sequences of two species of the genus Armadillidium: A. vulgare and A. pelagicum. All features of arthropods mtDNA control regions are present (origin of replication, poly-T stretch, GA and TA-rich blocks and one variable domain), plus a conserved IR. This IR can potentially fold into a hairpin structure and is present in two different orientations among the A. vulgare populations: either in one sense or in its reverse complement. This polymorphism, also observed in a single individual (heteroplasmy), might be a signature of genome architecture conversions from linear to circular monomeric mtDNA via successive opening and closing of the molecules.

Complete mitogenome sequences of four flatfishes (Pleuronectiformes) reveal a novel gene arrangement of L-strand coding genes

Background

Few mitochondrial gene rearrangements are found in vertebrates and large-scale changes in these genomes occur even less frequently. It is difficult, therefore, to propose a mechanism to account for observed changes in mitogenome structure. Mitochondrial gene rearrangements are usually explained by the recombination model or tandem duplication and random loss model.

Results

In this study, the complete mitochondrial genomes of four flatfishes, Crossorhombus azureus (blue flounder), Grammatobothus krempfi, Pleuronichthys cornutus, and Platichthys stellatus were determined. A striking finding is that eight genes in the C. azureus mitogenome are located in a novel position, differing from that of available vertebrate mitogenomes. Specifically, the ND6 and seven tRNA genes (the Q, A, C, Y, S₁, E, P genes) encoded by the L-strand have been translocated to a position between tRNA-T and tRNA-F though the original order of the genes is maintained.

Conclusions

These special features are used to suggest a mechanism for C. azureus mitogenome rearrangement. First, a dimeric molecule was formed by two monomers linked head-to-tail, then one of the two sets of promoters lost function and the genes controlled by the disabled promoters became pseudogenes, non-coding sequences, and even were lost from the genome. This study provides a new gene-rearrangement model that accounts for the events of gene-rearrangement in a vertebrate mitogenome.

Mitochondrial DNA of Clathrina clathrus (Calcarea, Calcinea): six linear chromosomes, fragmented rRNAs, tRNA editing, and a novel genetic code

Sponges (phylum Porifera) are a large and ancient group of morphologically simple but ecologically important aquatic animals. Although their body plan and lifestyle are relatively uniform, sponges show extensive molecular and genetic diversity. In particular, mitochondrial genomes from three of the four previously studied classes of Porifera (Demospongiae, Hexactinellida, and Homoscleromorpha) have distinct gene contents, genome organizations, and evolutionary rates. Here, we report the mitochondrial genome of Clathrina clathrus (Calcinea, Clathrinidae), a representative of the fourth poriferan class, the Calcarea, which proves to be the most unusual. Clathrina clathrus mitochondrial DNA (mtDNA) consists of six linear chromosomes 7.6-9.4 kb in size and encodes at least 37 genes: 13 protein codings, 2 ribosomal RNAs (rRNAs), and 24 transfer RNAs (tRNAs). Protein genes include atp9, which has now been found in all major sponge lineages, but no atp8. Our analyses further reveal the presence of a novel genetic code that involves unique reassignments of the UAG codons from termination to tyrosine and of the CGN codons from arginine to glycine. Clathrina clathrus mitochondrial rRNAs are encoded in three (srRNA) and ≥6 (lrRNA) fragments distributed out of order and on several chromosomes. The encoded tRNAs contain multiple mismatches in the aminoacyl acceptor stems that are repaired posttranscriptionally by 3'-end RNA editing. Although our analysis does not resolve the phylogenetic position of calcareous sponges, likely due to their high rates of mitochondrial sequence evolution, it confirms mtDNA as a promising marker for population studies in this group. The combination of unusual mitochondrial features in C. clathrus redefines the extremes of mtDNA evolution in animals and further argues against the idea of a "typical animal mtDNA."

Genetic aspects of mitochondrial genome evolution

Many years of extensive studies of metazoan mitochondrial genomes have established differences in gene arrangements and genetic codes as valuable phylogenetic markers. Understanding the underlying mechanisms of replication, transcription and the role of the control regions which cause e.g. different gene orders is important to assess the phylogenetic signal of such events. This review summarises and discusses, for the Metazoa, the general aspects of mitochondrial transcription and replication with respect to control regions as well as several proposed models of gene rearrangements. As whole genome sequencing projects accumulate, more and more observations about mitochondrial gene transfer to the nucleus are reported. Thus occurrence and phylogenetic aspects concerning nuclear mitochondrial-like sequences (NUMTS) is another aspect of this review.

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.

Widespread atypical mitochondrial DNA structure in isopods (Crustacea, Peracarida) related to a constitutive heteroplasmy in terrestrial species

Metazoan mitochondrial DNA (mtDNA) is generally composed of circular monomeric molecules. However, a few exceptions do exist and among them two terrestrial isopods Armadillidium vulgare and Porcellionides pruinosus have an atypical mtDNA composed of linear monomers associated with circular "head-to-head" dimers: a very unusual structure for animal mtDNA genome. To assess the distribution of this atypical mtDNA among isopods, we performed RFLP and Southern blot analyses on mtDNA of 16 terrestrial (Oniscidea family) and two aquatic isopod species: the marine Sphaeroma serratum (suborder Flabellifera, sister group of Oniscidea) and the freshwater Asellus aquaticus (Asellota, early derived taxon of isopod). The atypical mtDNA structure was observed in 15 terrestrial isopod species and A. aquaticus, suggesting a wide distribution of atypical mtDNA among isopods. However, a typical metazoan mtDNA structure was detected in the marine isopod S. serratum and the Oniscidea Ligia oceanica . Our results suggest two possible scenarios: an early origin of the atypical mtDNA in isopods followed by reversion to the typical ancestral mtDNA structure for several species, or a convergent appearance of the atypical mtDNA structure in two isopod suborders. We compare this distribution of the atypical mtDNA structure with the presence of a heteroplasmy also observed in the mtDNA of several terrestrial isopod species. We discuss if this transmitted heteroplasmy is vectored by the atypical mtDNA and its impact on the maintenance of the atypical mtDNA in isopods.

Evolution of linear mitochondrial genomes in medusozoan cnidarians

In nearly all animals, mitochondrial DNA (mtDNA) consists of a single circular molecule that encodes several subunits of the protein complexes involved in oxidative phosphorylation as well as part of the machinery for their expression. By contrast, mtDNA in species belonging to Medusozoa (one of the two major lineages in the phylum Cnidaria) comprises one to several linear molecules. Many questions remain on the ubiquity of linear mtDNA in medusozoans and the mechanisms responsible for its evolution, replication, and transcription. To address some of these questions, we determined the sequences of nearly complete linear mtDNA from 24 species representing all four medusozoan classes: Cubozoa, Hydrozoa, Scyphozoa, and Staurozoa. All newly determined medusozoan mitochondrial genomes harbor the 17 genes typical for cnidarians and map as linear molecules with a high degree of gene order conservation relative to the anthozoans. In addition, two open reading frames (ORFs), polB and ORF314, are identified in cubozoan, schyphozoan, staurozoan, and trachyline hydrozoan mtDNA. polB belongs to the B-type DNA polymerase gene family, while the product of ORF314 may act as a terminal protein that binds telomeres. We posit that these two ORFs are remnants of a linear plasmid that invaded the mitochondrial genomes of the last common ancestor of Medusozoa and are responsible for its linearity. Hydroidolinan hydrozoans have lost the two ORFs and instead have duplicated cox1 at each end of their mitochondrial chromosome(s). Fragmentation of mtDNA occurred independently in Cubozoa and Hydridae (Hydrozoa, Hydroidolina). Our broad sampling allows us to reconstruct the evolutionary history of linear mtDNA in medusozoans.

The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods

The complete mitochondrial (mt) genome sequence of the Australian fresh water isopod Eophreatoicus sp.-14 has been determined. The new species is a member of the taxon Phreatoicidea, a clade of particular interest, as it is often regarded as the sister group to all other Isopoda. Although the overall genome organization of Eophreatoicus sp.-14 conforms to the typical state of Metazoa--it is a circular ring of DNA hosting the usual 37 genes and one major non-coding region--it bears a number of derived characters that fall within the scope of "genome morphology". Earlier studies have indicated that the isopod mitochondrial gene order is not as conserved as that of other crustaceans. Indeed, the mt genome of Eophreatoicus sp.-14 shows an inversion of seven genes (including cox1), which is as far as we know unique. Even more interesting is the derived arrangement of nad1, trnL(CUN), rrnS, control region, cob, trnT, nad5 and trnF that is shared by nearly all available isopod mt genomes. A striking feature is the close proximity of the rearranged genes to the mt control region. Inferable gene translocation events are, however, more suitable to trace the evolution of mt genomes. Genes like nad1/trnL(CUN) and nad5/trnF, which retained their adjacent position after being rearranged, were most likely translocated together. A very good example for the need to understand the mechanisms of translocations is the remolding of trnL(UUR) to trnL(CUN). Both tRNA genes are adjacent and have a high sequence similarity, probably the result of a gene duplication and subsequent anticodon mutation. Modified secondary structures were found in three tRNAs of Eophreatoicus sp.-14, which are all characterized by the loss of the DHU-arm. This is common to crustaceans for tRNA Serine(AGY), while the arm-loss in tRNA Cysteine within Malacostraca is only shared by other isopods. Modification of the third tRNA, Isoleucine, is not known from any other related species. Nucleotide frequencies of genes have been found to be indirectly correlated to the orientation of the mitochondrial replication process. In Eophreatoicus sp.-14 and in other Isopoda the associated nucleotide bias is inversed to the state of other Malacostraca. This is a strong indication for an inversion of the control region that most likely evolved in the isopod ancestor.

A thirty million year-old inherited heteroplasmy

Due to essentially maternal inheritance and a bottleneck effect during early oogenesis, newly arising mitochondrial DNA (mtDNA) mutations segregate rapidly in metazoan female germlines. Consequently, heteroplasmy (i.e. the mixture of mtDNA genotypes within an organism) is generally resolved to homoplasmy within a few generations. Here, we report an exceptional transpecific heteroplasmy (predicting an alanine/valine alloacceptor tRNA change) that has been stably inherited in oniscid crustaceans for at least thirty million years. Our results suggest that this heteroplasmy is stably transmitted across generations because it occurs within mitochondria and therefore escapes the mtDNA bottleneck that usually erases heteroplasmy. Consistently, at least two oniscid species possess an atypical trimeric mitochondrial genome, which provides an adequate substrate for the emergence of a constitutive intra-mitochondrial heteroplasmy. Persistence of a mitochondrial polymorphism on such a deep evolutionary timescale suggests that balancing selection may be shaping mitochondrial sequence evolution in oniscid crustaceans.

Structure and Evolution of the Atypical Mitochondrial Genome of Armadillidium vulgare (Isopoda, Crustacea)

The crustacean isopod Armadillidium vulgare is characterized by an unusual approximately 42-kb-long mitochondrial genome consisting of two molecules co-occurring in mitochondria: a circular approximately 28-kb dimer formed by two approximately 14-kb monomers fused in opposite polarities and a linear approximately 14-kb monomer. Here we determined the nucleotide sequence of the fundamental monomeric unit of A. vulgare mitochondrial genome, to gain new insight into its structure and evolution. Our results suggest that the junction zone between monomers of the dimer structure is located in or near the control region. Direct sequencing indicated that the nucleotide sequences of the different monomer units are virtually identical. This suggests that gene conversion and/or replication processes play an important role in shaping nucleotide sequence variation in this mitochondrial genome. The only heteroplasmic site we identified predicts an alloacceptor tRNA change from tRNA(Ala) to tRNA(Val). Therefore, in A. vulgare, tRNA(Ala) and tRNA(Val) are found at the same locus in different monomers, ensuring that both tRNAs are present in mitochondria. The presence of this heteroplasmic site in all sequenced individuals suggests that the polymorphism is selectively maintained, probably because of the necessity of both tRNAs for maintaining proper mitochondrial functions. Thus, our results provide empirical evidence for the tRNA gene recruitment model of tRNA evolution. Moreover, interspecific comparisons showed that the A. vulgare mitochondrial gene order is highly derived compared to the putative ancestral arthropod type. By contrast, an overall high conservation of mitochondrial gene order is observed within crustacean isopods.

Mitochondrial genome haplotype hypervariation within the isopod parasitic nematode Thaumamermis cosgrovei

Characterization of mitochondrial genomes from individual Thaumamermis cosgrovei nematodes, obligate parasites of the isopod Armadillidium vulgare, revealed that numerous mtDNA haplotypes, ranging in size from 19 to 34 kb, are maintained in several spatially separated isopod populations. The magnitude and frequency of conspecific mtDNA size variation is unprecedented among all studied size-polymorphic metazoan mitochondrial genomes. To understand the molecular basis of this hypervariation, complete nucleotide sequences of two T. cosgrovei mtDNA haplotypes were determined. A hypervariable segment, residing between the atp6 and rrnL genes, contributes exclusively to T. cosgrovei mtDNA size variation. Within this region, mtDNA coding genes and putative nonfunctional sequences have accumulated substitutions and are duplicated and rearranged to varying extents. Hypervariation at this level has enabled a first insight into the life history of T. cosgrovei. In five A. vulgare hosts infected with multiple nematodes, four carried nematodes with identical mtDNA haplotypes, suggesting that hosts may become infected by ingesting a recently hatched egg clutch or become parasitized by individuals from the same brood prior to dispersal of siblings within the soil.

Comparative genomics of large mitochondria in placozoans

The first sequenced mitochondrial genome of a placozoan, Trichoplax adhaerens, challenged the conventional wisdom that a compact mitochondrial genome is a common feature among all animals. Three additional placozoan mitochondrial genomes representing highly divergent clades have been sequenced to determine whether the large Trichoplax mtDNA is a shared feature among members of the phylum Placozoa or a uniquely derived condition. All three mitochondrial genomes were found to be very large, 32- to 37-kb, circular molecules, having the typical 12 respiratory chain genes, 24 tRNAs, rnS, and rnL. They share with the Trichoplax mitochondrial genome the absence of atp8, atp9, and all ribosomal protein genes, the presence of several cox1 introns, and a large open reading frame containing an intron group I LAGLIDADG endonuclease domain. The differences in mtDNA size within Placozoa are due to variation in intergenic spacer regions and the presence or absence of long open reading frames of unknown function. Phylogenetic analyses of the 12 respiratory chain genes support the monophyly of Placozoa. The similarities in composition and structure between the three mitochondrial genomes reported here and that of Trichoplax's mtDNA suggest that their uncompacted state is a shared ancestral feature to other nonmetazoans while their gene content is a derived feature shared only among the Metazoa.

The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features

Background

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.

Results

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.

Conclusion

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.

Novel Repetitive Structures, Deviant Protein-Encoding Sequences andUnidentified ORFs in the Mitochondrial Genome of the BrachiopodLingula anatina

Complete sequence determination of the brachiopod Lingula anatina mtDNA (28,818 bp) revealed an organization that is remarkably atypical for an animal mt-genome. In addition to the usual set of 37 animal mitochondrial genes, which make up only 57% (16,555 bp) of the entire sequence, the genome contains lengthy unassigned sequences. All the genes are encoded in the same DNA strand, generally in a compact way, whereas the overall gene order is highly divergent in comparison with known animal mtDNA. Individual genes are generally longer and deviate considerably in sequence from their homologues in other animals. The genome contains two major repeat regions, in which 11 units of unassigned sequences and six genes (atp8, trnM, trnQ, trnV, and part of cox2 and nad2) are found in repetition, in the form of nested direct repeats of unparalleled complexity. One of the repeat regions contains unassigned repeat units dispersed among several unique sequences, novel repetitive structure for animal mtDNAs. Each of those unique sequences contains an open reading frame for a polypeptide between 80 and 357 amino acids long, potentially encoding a functional molecule, but none of them has been identified with known proteins. In both repeat regions, tRNA genes or tRNA gene-like sequences flank major repeated units, supporting the view that those structures play a role in the mitochondrial gene rearrangements. Although the intricate repeated organization of this genome can be explained by recurrent tandem duplications and subsequent deletions mediated by replication errors, other mechanisms, such as nonhomologous recombinations, appear to explain certain structures more easily.

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

We determined the complete 16,354 bp mitochondrial DNA sequence of the Chinese mitten crab Eriocheir japonica sinesnesis. It consists of 13 protein coding genes (PCGs), 2 rRNAs, and 22 tRNAs, typical of metazoan mitochondrial genomes. With respect to the ancestral crustacean mt gene order, a PCG, the rRNAs, and 12 of the tRNAs appear to be rearranged. This rearrangement is novel in the arthropods and suggests an accelerated rate of mt genome rearrangement in this brachyuran lineage based on the relative rate of gene rearrangement. It is typical in arthropods that all of the rearranged genes or gene blocks take place at both nad3-nad5 and nad5-nad4 gene junctions. Such occurrence additionally revealed two target regions of frequent rearrangement in mitochondrial genomes of decapods, even in that of the non-hexapod arthropods according to our comparative studies among 32 taxa. Additionally, selective constraint on sharing the single introducing location is apparent for most of the rearrangements that occurred at the nad3-nad5 gene junction of these taxa. The gene arrangement features at both gene junctions allow the reconstruction of relationships among the advanced decapods. These features are therefore characteristic molecular markers in phylogenetic inference. The genomic organization differences at both gene junctions provide new evidence of extremely divergence between Heterotremata and Thoracotremata in brachyuran crabs. A duplication/loss (random and nonrandom)-recombination model was proposed to account for the generation of the gene order in E. japonica sinesnesis under the guide of intergenic spacers.

Sequencing and comparing whole mitochondrial genomes of animals

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, which can be especially powerful. We describe here the protocols commonly used for physically isolating mitochondrial DNA (mtDNA), for amplifying these by polymerase chain reaction (PCR) or rolling circle amplification (RCA), for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences with determining and comparing complete mitochondrial DNA sequences.

Mitochondria of Protists

Over the past several decades, our knowledge of the origin and evolution of mitochondria has been greatly advanced by determination of complete mitochondrial genome sequences. Among the most informative mitochondrial genomes have been those of protists (primarily unicellular eukaryotes), some of which harbor the most gene-rich and most eubacteria-like mitochondrial DNAs (mtDNAs) known. Comparison of mtDNA sequence data has provided insights into the radically diverse trends in mitochondrial genome evolution exhibited by different phylogenetically coherent groupings of eukaryotes, and has allowed us to pinpoint specific protist relatives of the multicellular eukaryotic lineages (animals, plants, and fungi). This comparative genomics approach has also revealed unique and fascinating aspects of mitochondrial gene expression, highlighting the mitochondrion as an evolutionary playground par excellence.

Novel rearrangements of arthropod mitochondrial DNA detected with long-PCR: applications to arthropod phylogeny and evolution

Rearrangements of mitochondrial DNA gene order have been suggested as a tool for defining the pattern of evolutionary divergence in arthropod taxa. We have employed a combination of highly conserved insect-based polymerase chain reaction (PCR) primers with long-PCR to survey 14 noninsect arthropods for mitochondrial gene rearrangements. The size of the amplified fragments was used to order the primer containing genes. Five chelicerates exhibit amplicons that are consistent with the presumptive ancestral arthropod mtDNA gene order. These five species comprise two soft ticks, two prostriate hard ticks, and an opilionid. Six other chelicerates, all metastriate hard ticks, have a different arrangement that was originally discovered by this procedure and has been previously detailed in a complete mtDNA sequence. Three new arthropod mtDNA gene arrangements are described here. They were discovered in a terrestrial crustacean (Isopoda) and two myriapods (Chilopoda, centipede; Diplopoda, millipede). These rearrangements include major realignments of some of the large coding regions and two possible new positions for the tRNA(Met) (M) gene in arthropods. The long-PCR approach affords an opportunity to quickly screen divergent taxa for major rearrangements. Taxa exhibiting rearrangements can be targeted for DNA sequencing of gene boundaries to establish the details of the mtDNA organization.

Complete mtDNA sequences of two millipedes suggest a new model for mitochondrial gene rearrangements: duplication and nonrandom loss

We determined the complete mitochondrial DNA (mtDNA) sequences of the millipedes Narceus annularus and Thyropygus sp. (Arthropoda: Diplopoda) and identified, in both genomes, all 37 genes typical for metazoan mtDNA. The arrangement of these genes is identical in the two millipedes, but differs from others found in arthropod mtDNAs in the location of at least four genes or gene blocks. This novel gene arrangement is unusual for animal mtDNA in that genes with identical transcriptional polarities are clustered in the genome, and the two clusters are separated by two noncoding regions. The only exception to this pattern is the gene for cysteine tRNA, which is located in the part of the genome that otherwise contains all genes with the opposite transcriptional polarity. We suggest that a mechanism involving complete mtDNA duplication followed by the loss of genes, predetermined by their transcriptional polarity and location in the genome, could generate this gene arrangement from the one ancestral for arthropods. The proposed mechanism has important implications for phylogenetic inferences that are drawn on the basis of gene arrangement comparisons.

One-step PCR amplification of complete arthropod mitochondrial genomes

A new PCR primer set which enables one-step amplification of complete arthropod mitochondrial genomes was designed from two conserved 16S rDNA regions for the long PCR technique. For this purpose, partial 16S rDNAs amplified with universal primers 16SA and 16SB were newly sequenced from six representative arthropods: Armadillidium vulgare and Macrobrachium nipponense (Crustacea), Anopheles sinensis (Insecta), Lithobius forficatus and Megaphyllum sp. (Myriapoda), and Limulus polyphemus (Chelicerata). The genomic locations of two new primers, HPK16Saa and HPK16Sbb, correspond to positions 13314-13345 and 12951-12984, respectively, in the Drosophila yakuba mitochondrial genome. The usefulness of the primer set was experimentally examined and confirmed with five of the representative arthropods, except for A. vulgare, which has a linearized mitochondrial genome. With this set, therefore, we could easily and rapidly amplify complete mitochondrial genomes with small amounts of arthropod DNA. Although the primers suggested here were examined only with arthropod groups, a possibility of successful application to other invertebrates is very high, since the high degree of sequence conservation is shown on the primer sites in other invertebrates. Thus, this primer set can serve various research fields, such as molecular evolution, population genetics, and molecular phylogenetics based on DNA sequences, RFLP, and gene rearrangement of mitochondrial genomes in arthropods and other invertebrates.

Mitochondrial DNA polymorphism, sex ratio distorters and population genetics in the isopod Armadillidium vulgare

Two maternally inherited sex ratio distorters (SRD) impose female-biased sex ratios on the wood louse Armadillidium vulgare by feminizing putative males. These SRD are (i) an intracytoplasmic bacterium of the genus Wolbachia, and (ii) another non-Mendelian element of unknown nature: the f element. Mitochondrial DNA variation was investigated in A. vulgare field populations to trace the evolution of host-SRD relationships and to investigate the effect of SRD on host cytoplasmic polymorphism. The Wolbachia endosymbionts showed no polymorphism in their ITS2 sequence and were associated with two closely related mitochondrial types. This situation probably reflects a single infection event followed by a slight differentiation of mitochondria. There was no association between the f element and a given mitochondrial type, which may confirm the fact that this element can be partially paternally transmitted. The spreading of a maternally inherited SRD in a population should reduce the mitochondrial diversity by a hitchhiking process. In A. vulgare, however, a within-population mtDNA polymorphism was often found, because of the deficient spread of Wolbachia and the partial paternal inheritance of the f element. The analysis of molecular variance indicated that A. vulgare populations are genetically structured, but without isolation by distance.

Mitochondrial DNA variability and wolbachia infection in two sibling woodlice species

Several morphological races and subspecies have been described and later included within the terrestrial isopod species Porcellionides pruinosus. During our study of this species, we have worked on specimens from France, Greece, Tunisia and Reunion island. Laboratory crosses have revealed two separate groups of populations: French populations (four localities) in one group, and those from Tunisia, Reunion island and Greece in the other. French individuals were reproductively isolated from those of the other populations. We have undertaken a survey of mitochondrial DNA (mtDNA) polymorphism in these seven populations. We observed two groups of mitotypes corresponding to the two groups of populations. Interfertility experiments between populations and the mitochondrial genetic distances between mitotypes both suggest the presence of two different species, one in France and one in Greece, Tunisia and Reunion island. The two species harbour, respectively, two different Wolbachia lines. Another feature of the molecular genetic analysis was the apparent mitochondrial monomorphism in the French populations and the low variability in the other three populations. The result can be related to the possibility of Wolbachia-induced genetic hitchhiking in these populations.