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

Phylogenetic analyses of endoparasitic Acanthocephala based on mitochondrial genomes suggest secondary loss of sensory organs

The metazoan taxon Syndermata (Monogononta, Bdelloidea, Seisonidea, Acanthocephala) comprises species with vastly different lifestyles. The focus of this study is on the phylogeny within the syndermatan subtaxon Acanthocephala (thorny-headed worms, obligate endoparasites). In order to investigate the controversially discussed phylogenetic relationships of acanthocephalan subtaxa we have sequenced the mitochondrial (mt) genomes of Echinorhynchus truttae (Palaeacanthocephala), Paratenuisentis ambiguus (Eoacanthocephala), Macracanthorhynchus hirudinaceus (Archiacanthocephala), and Philodina citrina (Bdelloidea). In doing so, we present the largest molecular phylogenetic dataset so far for this question comprising all major subgroups of Acanthocephala. Alongside with publicly available mt genome data of four additional syndermatans as well as 18 other lophotrochozoan (spiralian) taxa and one outgroup representative, the derived protein-coding sequences were used for Maximum Likelihood as well as Bayesian phylogenetic analyses. We achieved entirely congruent results, whereupon monophyletic Archiacanthocephala represent the sister taxon of a clade comprising Eoacanthocephala and monophyletic Palaeacanthocephala (Echinorhynchida). This topology suggests the secondary loss of lateral sensory organs (sensory pores) within Palaeacanthocephala and is further in line with the emergence of apical sensory organs in the stem lineage of Archiacanthocephala.

Comparative mitogenomic analysis of damsel bugs representing three tribes in the family Nabidae (Insecta: Hemiptera)

BACKGROUND

Nabidae, a family of predatory heteropterans, includes two subfamilies and five tribes. We previously reported the complete mitogenome of Alloeorhynchus bakeri, a representative of the tribe Prostemmatini in the subfamily Prostemmatinae. To gain a better understanding of architecture and evolution of mitogenome in Nabidae, mitogenomes of five species representing two tribes (Gorpini and Nabini) in the subfamily Nabinae were sequenced, and a comparative mitogenomic analysis of three nabid tribes in two subfamilies was carried out.

METHODOLOGY/PRINCIPAL FINDINGS

Nabid mitogenomes share a similar nucleotide composition and base bias, except for the control region, where differences are observed at the subfamily level. In addition, the pattern of codon usage is influenced by the GC content and consistent with the standard invertebrate mitochondrial genetic code and the preference for A+T-rich codons. The comparison among orthologous protein-coding genes shows that different genes have been subject to different rates of molecular evolution correlated with the GC content. The stems and anticodon loops of tRNAs are extremely conserved, and the nucleotide substitutions are largely restricted to TψC and DHU loops and extra arms, with insertion-deletion polymorphisms. Comparative analysis shows similar rates of substitution between the two rRNAs. Long non-coding regions are observed in most Gorpini and Nabini mtDNAs in-between trnI-trnQ and/or trnS2-nad1. The lone exception, Nabis apicalis, however, has lost three tRNAs. Overall, phylogenetic analysis using mitogenomic data is consistent with phylogenies constructed mainly form morphological traits.

CONCLUSIONS/SIGNIFICANCE

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the family Nabidae. Nucleotide diversity and mitogenomic traits are phylogenetically informative at subfamily level. Furthermore, inclusion of a broader range of samples representing various taxonomic levels is critical for the understanding of mitogenomic evolution in damsel bugs.

The complete mitochondrial DNA sequence of Crotalus horridus (timber rattlesnake)

The complete mitogenome of the timber rattlesnake (Crotalus horridus) was completed using Sanger sequencing. It is 17,260 bp with 13 protein-coding genes, 21 tRNAs, two rRNAs and two control regions. Gene synteny is consistent with other snakes with the exception of a missing redundant tRNA (Ser) . This mitogenome should prove to be a useful addition of a well-known member of the Viperidae snake family.

The mitochondrial genome of Paraminabea aldersladei (Cnidaria: Anthozoa: Octocorallia) supports intramolecular recombination as the primary mechanism of gene rearrangement in octocoral mitochondrial genomes

Sequencing of the complete mitochondrial genome of the soft coral Paraminabea aldersladei (Alcyoniidae) revealed a unique gene order, the fifth mt gene arrangement now known within the cnidarian subclass Octocorallia. At 19,886 bp, the mt genome of P. aldersladei is the second largest known for octocorals; its gene content and nucleotide composition are, however, identical to most other octocorals, and the additional length is due to the presence of two large, noncoding intergenic regions. Relative to the presumed ancestral octocoral gene order, in P. aldersladei a block of three protein-coding genes (nad6–nad3–nad4l) has been translocated and inverted. Mapping the distribution of mt gene arrangements onto a taxonomically comprehensive phylogeny of Octocorallia suggests that all of the known octocoral gene orders have evolved by successive inversions of one or more evolutionarily conserved blocks of protein-coding genes. This mode of genome evolution is unique among Metazoa, and contrasts strongly with that observed in Hexacorallia, in which extreme gene shuffling has occurred among taxonomic orders. Two of the four conserved gene blocks found in Octocorallia are, however, also conserved in the linear mt genomes of Medusozoa and in one group of Demospongiae. We speculate that the rate and mechanism of gene rearrangement in octocorals may be influenced by the presence in their mt genomes of mtMutS, a putatively active DNA mismatch repair protein that may also play a role in mediating intramolecular recombination.

Replication-Associated Mutational Pressure (RMP) Governs Strand-Biased Compositional Asymmetry (SCA) and Gene Organization in Animal Mitochondrial Genomes

The nucleotide composition of the light (L-) and heavy (H-) strands of animal mitochondrial genomes is known to exhibit strand-biased compositional asymmetry (SCA). One of the possibilities is the existence of a replication-associated mutational pressure (RMP) that may introduce characteristic nucleotide changes among mitochondrial genomes of different animal lineages. Here, we discuss the influence of RMP on nucleotide and amino acid compositions as well as gene organization. Among animal mitochondrial genomes, RMP may represent the major force that compels the evolution of mitochondrial protein-coding genes, coupled with other process-based selective pressures, such as on components of translation machinery— tRNAs and their anticodons. Through comparative analyses of sequenced mitochondrial genomes among diverse animal lineages and literature reviews, we suggest a strong RMP effect, observed among invertebrate mitochondrial genes as compared to those of vertebrates, that is either a result of positive selection on the invertebrate or a relaxed selective pressure on the vertebrate mitochondrial genes.

A unique tRNA gene family and a novel, highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS(-AGY) and trnQ genes, a trnM isomerism, but it lacks trnS(-CUN). We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP→trnS(-AGY)) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, "Oxidored_q1_N", which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.

New features of Asian Crassostrea oyster mitochondrial genomes: a novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.

The complete mitochondrial genome of the giant pill millipede, Sphaerotheriidae sp. (Myriapoda: Diplopoda: Sphaerotheriida)

The complete mitochondrial genome (mitogenome) of the millipede Sphaerotheriidae sp. has been studied. The genome is 14,970 bp long and contains the typical complement of 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. Gene order in Sphaerotheriidae sp. mitogenome is assumed to represent the myriapod ground pattern, which is shared by myriapod-chelicerate clade.

The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes)

An exponential growth of mitochondrial genome information has brought significant progress in understanding the organismal phylogeny and mitochondrial genome evolution for many metazoans including platyhelminth groups. In this study, we determined the complete mitochondrial genome sequence for Benedenia hoshinai, an ectoparasitic monogenean species, and compared it with its congener Benedenia seriolae. The complete mitochondrial genome is 13,554 bp in length and contains 12 protein-coding genes (lacking the atp8 gene), 2 rRNA genes, and 22 tRNA genes, all encoded in the same direction as found in all other platyhelminth species sequenced to date. The gene arrangement of B. hoshinai mtDNA is almost identical to B. seriolae, differing only by the translocation of trnT between cox1 and rrnL. It is unclear whether the shared position of trnT between B. hoshinai and Gyrodactylus represents evidence for their phylogenetic affinity; testing this hypothesis requires further mitogenomic evidence.

Characterization of the complete mitochondrial genome sequence of Spirometra erinaceieuropaei (Cestoda: Diphyllobothriidae) from China

Sparganosis, caused by the plerocercoid larvae of members of the genus Spirometra, can cause significant public health problem and considerable economic losses. In the present study, the complete mitochondrial DNA (mtDNA) sequence of Spirometra erinaceieuropaei from China was determined, characterized and compared with that of S. erinaceieuropaei from Japan. The gene arrangement in the mt genome sequences of S. erinaceieuropaei from China and Japan is identical. The identity of the mt genomes was 99.1% between S. erinaceieuropaei from China and Japan, and the complete mtDNA sequence of S. erinaceieuropaei from China is slightly shorter (2 bp) than that from Japan. Phylogenetic analysis of S. erinaceieuropaei with other representative cestodes using two different computational algorithms [Bayesian inference (BI) and maximum likelihood (ML)] based on concatenated amino acid sequences of 12 protein-coding genes, revealed that S. erinaceieuropaei is closely related to Diphyllobothrium spp., supporting classification based on morphological features. The present study determined the complete mtDNA sequences of S. erinaceieuropaei from China that provides novel genetic markers for studying the population genetics and molecular epidemiology of S. erinaceieuropaei in humans and animals.

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.

Pseudoscorpion mitochondria show rearranged genes and genome-wide reductions of RNA gene sizes and inferred structures, yet typical nucleotide composition bias

BACKGROUND

Pseudoscorpions are chelicerates and have historically been viewed as being most closely related to solifuges, harvestmen, and scorpions. No mitochondrial genomes of pseudoscorpions have been published, but the mitochondrial genomes of some lineages of Chelicerata possess unusual features, including short rRNA genes and tRNA genes that lack sequence to encode arms of the canonical cloverleaf-shaped tRNA. Additionally, some chelicerates possess an atypical guanine-thymine nucleotide bias on the major coding strand of their mitochondrial genomes.

RESULTS

We sequenced the mitochondrial genomes of two divergent taxa from the chelicerate order Pseudoscorpiones. We find that these genomes possess unusually short tRNA genes that do not encode cloverleaf-shaped tRNA structures. Indeed, in one genome, all 22 tRNA genes lack sequence to encode canonical cloverleaf structures. We also find that the large ribosomal RNA genes are substantially shorter than those of most arthropods. We inferred secondary structures of the LSU rRNAs from both pseudoscorpions, and find that they have lost multiple helices. Based on comparisons with the crystal structure of the bacterial ribosome, two of these helices were likely contact points with tRNA T-arms or D-arms as they pass through the ribosome during protein synthesis.The mitochondrial gene arrangements of both pseudoscorpions differ from the ancestral chelicerate gene arrangement. One genome is rearranged with respect to the location of protein-coding genes, the small rRNA gene, and at least 8 tRNA genes. The other genome contains 6 tRNA genes in novel locations. Most chelicerates with rearranged mitochondrial genes show a genome-wide reversal of the CA nucleotide bias typical for arthropods on their major coding strand, and instead possess a GT bias. Yet despite their extensive rearrangement, these pseudoscorpion mitochondrial genomes possess a CA bias on the major coding strand. Phylogenetic analyses of all 13 mitochondrial protein-coding gene sequences consistently yield trees that place pseudoscorpions as sister to acariform mites.

CONCLUSION

The well-supported phylogenetic placement of pseudoscorpions as sister to Acariformes differs from some previous analyses based on morphology. However, these two lineages share multiple molecular evolutionary traits, including substantial mitochondrial genome rearrangements, extensive nucleotide substitution, and loss of helices in their inferred tRNA and rRNA structures.

Novel gene rearrangements in the mitochondrial genome of a webspinner, Aposthonia japonica (Insecta: Embioptera)

Webspinners (order Embioptera) are polyneopteran insects characterized by enlarged foretarsi with silk glands, whose silk is used to produce galleries in which the insects live gregariously. The phylogenetic position of webspinners has been debated. In the present study, an almost complete mitochondrial DNA (mtDNA) sequence of Embioptera is reported for the first time. The mtDNA of a webspinner, Aposthonia japonica , has the 13 protein-coding genes (PCGs) generally found in metazoan mtDNA sequences. There is a translocation of a large region including atp6, atp8, cox3, nad3, and nad5 as well as a duplication of the 12S rRNA gene. The rearrangement does not seem to affect nucleotide composition, although amino acid composition in some parts of the mtDNA is biased compared with other Polyneoptera species. Based on phylogenetic analyses using nucleotide sequences of all PCGs concatenated with two rRNA genes and the amino acid sequences of all PCGs, A. japonica is sister to Verophasmatodea, a suborder of typical stick and leaf insects.

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.

Nucleotide composition of CO1 sequences in Chelicerata (Arthropoda): detecting new mitogenomic rearrangements

Here we study the evolution of nucleotide composition in third codon-positions of CO1 sequences of Chelicerata, using a phylogenetic framework, based on 180 taxa and three markers (CO1, 18S, and 28S rRNA; 5,218 nt). The analyses of nucleotide composition were also extended to all CO1 sequences of Chelicerata found in GenBank (1,701 taxa). The results show that most species of Chelicerata have a positive strand bias in CO1, i.e., in favor of C nucleotides, including all Amblypygi, Palpigradi, Ricinulei, Solifugae, Uropygi, and Xiphosura. However, several taxa show a negative strand bias, i.e., in favor of G nucleotides: all Scorpiones, Opisthothelae spiders and several taxa within Acari, Opiliones, Pseudoscorpiones, and Pycnogonida. Several reversals of strand-specific bias can be attributed to either a rearrangement of the control region or an inversion of a fragment containing the CO1 gene. Key taxa for which sequencing of complete mitochondrial genomes will be necessary to determine the origin and nature of mtDNA rearrangements involved in the reversals are identified. Acari, Opiliones, Pseudoscorpiones, and Pycnogonida were found to show a strong variability in nucleotide composition. In addition, both mitochondrial and nuclear genomes have been affected by higher substitution rates in Acari and Pseudoscorpiones. The results therefore indicate that these two orders are more liable to fix mutations of all types, including base substitutions, indels, and genomic rearrangements.

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

BACKGROUND

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

PRINCIPAL FINDINGS

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

CONCLUSIONS/SIGNIFICANCE

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.

Complete mitochondrial genomes of two cockroaches, Blattella germanica and Periplaneta americana, and the phylogenetic position of termites

The mitochondrial genomes are one of the most information-rich markers in phylogenetics. The relationships within superorder Dictyoptera have been debated in the literature. However, the closely related termites (Isoptera) are retained as unranked taxon within the order Blattaria (cockroaches). In this work, we sequenced the complete mitogenomes of two cockroaches, reconstructed the molecular phylogeny and attempted to infer the phylogenetic position of termites in Blattaria more reliably. The complete mtDNA nucleotide sequences of the peridomestic American cockroach (Periplaneta americana L.) and the domestic German cockroach (Blattella germanica L.) are 15,025 and 15,584 bp in size, respectively. The genome shares the gene order and orientation with previously known Blattaria mitogenomes. Most tRNAs could be folded into the typical cloverleaf secondary structure, but the tRNA-Ser (AGN) of P. americana appears to be missing the dihydrouridine arm. Using nucleotide and amino acid sequences as phylogenetic markers, we proposed that termites should be treated as a superfamily (Termitoidea) of cockroaches. We suggested that Polyphagoidea was the sister group of Termitoidea in Blattaria and supported that the suborder Caelifera is more closely related to the Phasmatodea than to the suborder Ensifera of Orthoptera.

The complete mitochondrial genome of the rodent intra-arterial nematodes Angiostrongylus cantonensis and Angiostrongylus costaricensis

The two rodent intra-arterial nematodes, Angiostrongylus cantonensis and Angiostrongylus costaricensis, can cause human ill-health. The present study aimed to characterize and compare the mitochondrial (mt) genomes of these two species, and clarify their phylogenetic relationship and the position in the phylum Nematoda. The complete mt genomes of A. cantonensis and A. costaricensis are 13,497 and 13,585 bp in length, respectively. Hence, they are the smallest in the class of Chromadorea characterized thus far. Like many nematode species in the class of Chromadorea, they encode 12 proteins, 22 transfer RNAs, and two ribosomal RNAs. All genes are located on the same strand. Nucleotide identity of the two mt genomes is 81.6%, ranging from 77.7% to 87.1% in individual gene pairs. Our mt genome-wide analysis identified three major gene arrangement patterns (II-1, II-2, and II-3) from 48 nematode mt genomes. Both patterns II-1 and II-2 are distinct from pattern II-3, which covers the Spirurida, supporting a closer relationship between Ascaridida and Strongylida rather than Spirurida. Thymine (T) was highly concentrated on coding strands in Chromadorea, but balanced between the two strands in Enoplea, probably due to the gene arrangement pattern. Interestingly, the gene arrangement pattern of mt genomes and phylogenetic analysis based on concatenated amino acids indicated a closer relationship between the order Ascaridida and Rhabditida rather than Spirurida as indicated in previous studies. These discrepancies call for new research, reassessing the position of the order of Ascaridida in the phylogenetic tree. Once consolidated, the findings are important for population genetic studies and target identification.

The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae)

Many of true bugs are important insect pests to cultivated crops and some are important vectors of human diseases, but few cladistic analyses have addressed relationships among the seven infraorders of Heteroptera. The Enicocephalomorpha and Nepomorpha are consider the basal groups of Heteroptera, but the basal-most lineage remains unresolved. Here we report the mitochondrial genome of the unique-headed bug Stenopirates sp., the first mitochondrial genome sequenced from Enicocephalomorpha. The Stenopirates sp. mitochondrial genome is a typical circular DNA molecule of 15, 384 bp in length, and contains 37 genes and a large non-coding fragment. The gene order differs substantially from other known insect mitochondrial genomes, with rearrangements of both tRNA genes and protein-coding genes. The overall AT content (82.5%) of Stenopirates sp. is the highest among all the known heteropteran mitochondrial genomes. The strand bias is consistent with other true bugs with negative GC-skew and positive AT-skew for the J-strand. The heteropteran mitochondrial atp8 exhibits the highest evolutionary rate, whereas cox1 appears to have the lowest rate. Furthermore, a negative correlation was observed between the variation of nucleotide substitutions and the GC content of each protein-coding gene. A microsatellite was identified in the putative control region. Finally, phylogenetic reconstruction suggests that Enicocephalomorpha is the sister group to all the remaining Heteroptera.

The complete mitochondrial genome sequence of Oncicola luehei (Acanthocephala: Archiacanthocephala) and its phylogenetic position within Syndermata

In the present study, we determined the complete mitochondrial genome sequence of Oncicola luehei (14,281bp), the first archiacanthocephalan representative and the second complete sequence from the phylum Acanthocephala. The complete genome contains 36 genes including 12 protein coding genes, 22 transfer RNA (tRNA) genes and 2 ribosomal RNA genes (rrnL and rrnS) as reported for other syndermatan species. All genes are encoded on the same strand. The overall nucleotide composition of O. luehei mtDNA is 37.7% T, 29.6% G, 22.5% A, and 10.2% C. The overall A+T content (60.2%) is much lower, compared to other syndermatan species reported so far, due to the high frequency (18.3%) of valine encoded by GTN in its protein-coding genes. Results from phylogenetic analyses of amino acid sequences for 10 protein-coding genes from 41 representatives of major metazoan groups including O. luehei supported monophyly of the phylum Acanthocephala and of the clade Syndermata (Acanthocephala+Rotifera), and the paraphyly of the clade Eurotatoria (classes Bdelloidea+Monogononta from phylum Rotifera). Considering the position of the acanthocephalan species within Syndermata, it is inferred that obligatory parasitism characteristic of acanthocephalans was acquired after the common ancestor of acanthocephalans diverged from its sister group, Bdelloidea. Additional comparison of complete mtDNA sequences from unsampled acanthocephalan lineages, especially classes Polyacanthocephala and Eoacanthocephala, is required to test if mtDNA provides reliable information for the evolutionary relationships and pattern of life history diversification found in the syndermatan groups.

MitoZoa 2.0: a database resource and search tools for comparative and evolutionary analyses of mitochondrial genomes in Metazoa

The MITOchondrial genome database of metaZOAns (MitoZoa) is a public resource for comparative analyses of metazoan mitochondrial genomes (mtDNA) at both the sequence and genomic organizational levels. The main characteristics of the MitoZoa database are the careful revision of mtDNA entry annotations and the possibility of retrieving gene order and non-coding region (NCR) data in appropriate formats. The MitoZoa retrieval system enables basic and complex queries at various taxonomic levels using different search menus. MitoZoa 2.0 has been enhanced in several aspects, including: a re-annotation pipeline to check the correctness of protein-coding gene predictions; a standardized annotation of introns and of precursor ORFs whose functionality is post-transcriptionally recovered by RNA editing or programmed translational frameshifting; updates of taxon-related fields and a BLAST sequence similarity search tool. Database novelties and the definition of standard mtDNA annotation rules, together with the user-friendly retrieval system and the BLAST service, make MitoZoa a valuable resource for comparative and evolutionary analyses as well as a reference database to assist in the annotation of novel mtDNA sequences. MitoZoa is freely accessible at http://www.caspur.it/mitozoa.

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.

Probable presence of an ubiquitous cryptic mitochondrial gene on the antisense strand of the cytochrome oxidase I gene

BACKGROUND

Mitochondria mediate most of the energy production that occurs in the majority of eukaryotic organisms. These subcellular organelles contain a genome that differs from the nuclear genome and is referred to as mitochondrial DNA (mtDNA). Despite a disparity in gene content, all mtDNAs encode at least two components of the mitochondrial electron transport chain, including cytochrome c oxidase I (Cox1).

PRESENTATION OF THE HYPOTHESIS

A positionally conserved ORF has been found on the complementary strand of the cox1 genes of both eukaryotic mitochondria (protist, plant, fungal and animal) and alpha-proteobacteria. This putative gene has been named gau for gene antisense ubiquitous in mtDNAs. The length of the deduced protein is approximately 100 amino acids. In vertebrates, several stop codons have been found in the mt gau region, and potentially functional gau regions have been found in nuclear genomes. However, a recent bioinformatics study showed that several hypothetical overlapping mt genes could be predicted, including gau; this involves the possible import of the cytosolic AGR tRNA into the mitochondria and/or the expression of mt antisense tRNAs with anticodons recognizing AGR codons according to an alternative genetic code that is induced by the presence of suppressor tRNAs. Despite an evolutionary distance of at least 1.5 to 2.0 billion years, the deduced Gau proteins share some conserved amino acid signatures and structure, which suggests a possible conserved function. Moreover, BLAST analysis identified rare, sense-oriented ESTs with poly(A) tails that include the entire gau region. Immunohistochemical analyses using an anti-Gau monoclonal antibody revealed strict co-localization of Gau proteins and a mitochondrial marker.

TESTING THE HYPOTHESIS

This hypothesis could be tested by purifying the gau gene product and determining its sequence. Cell biological experiments are needed to determine the physiological role of this protein.

IMPLICATIONS OF THE HYPOTHESIS

Studies of the gau ORF will shed light on the origin of novel genes and their functions in organelles and could also have medical implications for human diseases that are caused by mitochondrial dysfunction. Moreover, this strengthens evidence for mitochondrial genes coded according to an overlapping genetic code.

Complete mitochondrial DNA sequence of the European flat oyster Ostrea edulis confirms Ostreidae classification

BACKGROUND

Because of its typical architecture, inheritance and small size, mitochondrial (mt) DNA is widely used for phylogenetic studies. Gene order is generally conserved in most taxa although some groups show considerable variation. This is particularly true in the phylum Mollusca, especially in the Bivalvia. During the last few years, there have been significant increases in the number of complete mitochondrial sequences available. For bivalves, 35 complete mitochondrial genomes are now available in GenBank, a number that has more than doubled in the last three years, representing 6 families and 23 genera. In the current study, we determined the complete mtDNA sequence of O. edulis, the European flat oyster. We present an analysis of features of its gene content and genome organization in comparison with other Ostrea, Saccostrea and Crassostrea species.

RESULTS

The Ostrea edulis mt genome is 16 320 bp in length and codes for 37 genes (12 protein-coding genes, 2 rRNAs and 23 tRNAs) on the same strand. As in other Ostreidae, O. edulis mt genome contains a split of the rrnL gene and a duplication of trnM. The tRNA gene set of O. edulis, Ostrea denselamellosa and Crassostrea virginica are identical in having 23 tRNA genes, in contrast to Asian oysters, which have 25 tRNA genes (except for C. ariakensis with 24). O. edulis and O. denselamellosa share the same gene order, but differ from other Ostreidae and are closer to Crassostrea than to Saccostrea. Phylogenetic analyses reinforce the taxonomic classification of the 3 families Ostreidae, Mytilidae and Pectinidae. Within the Ostreidae family the results also reveal a closer relationship between Ostrea and Saccostrea than between Ostrea and Crassostrea.

CONCLUSIONS

Ostrea edulis mitogenomic analyses show a high level of conservation within the genus Ostrea, whereas they show a high level of variation within the Ostreidae family. These features provide useful information for further evolutionary analysis of oyster mitogenomes.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

The complete mitochondrial genome of Calicogorgia granulosa (Anthozoa: Octocorallia): potential gene novelty in unidentified ORFs formed by repeat expansion and segmental duplication

Mitochondrial genomes of many nonbilaterian animals show high diversity of genome size and gene content, revealing many intergenic regions (IGRs), diverse repeats and additional genes. Here we present a new complete mitogenome of the cnidarian sea fan species, Calicogorgia granulosa (Anthozoa: Octocorallia) and its novel genomic features. The 20,246 bp of the complete mitogenome, which is the largest among the nine octocorals sequenced to date, contains 13 protein coding genes, 2 rRNAs and a tRNA within its circular form of mitochondrial DNA. We found an identical segmental duplication (S1 and S2, 913 bp) composed of an ORF (672 bp) coding for a hypothetical protein within which Direct Variant Repeat (DVR) expansions reside in-frame to the coding sequence. Additionally, the duplicated segmental DNA showed no variation in nucleotide sequences both between S1 and S2 and across multiple individual samples. Upon these observations, we discuss plausible causes for the intramitochondrial segmental duplication and the absence of sequence variation, and a need for further investigation of the novel ORF as well. In conclusion the present mitogenome of C. granulosa adds more information to our understanding of the diversity and evolution of mitogenomes of nonbilaterian animals.

The architecture and complete sequence of mitochondrial genome of an assassin bug Agriosphodrus dohrni (Hemiptera: Reduviidae)

The 16, 470 bp nucleotide sequence of the mitochondrial genome (mitogenome) of an assassin bug from the reduviid subfamily Harpactorinae, Agriosphodrus dohrni, has been revealed. The entire genome encodes for two ribosomal RNA genes (rrnL and rrnS), 22 transfer RNA (tRNA) genes, 13 protein-coding genes, and a control region. The nucleotide composition is biased toward adenine and thymine (A+T = 72.2%). Comparative analysis with two other reduviid species Triatoma dimidiata and Valentia hoffmanni, exhibited highly conserved genome architectures including genome contents, gene order, nucleotide composition, codon usage, amino acid composition, as well as genome asymmetry. All protein-coding genes use standard mitochondrial initiation codons (methionine and isoleucine), except that nad1 starts with GTG. All tRNAs have the classic clover-leaf structure, except that the dihydrouridine (DHU) arm of tRNA(Ser(AGN)) forms a simple loop. Secondary structure comparisons of the two mitochondrial ribosomal subunits among sequenced assassin bugs show that the sequence and structure of rrnL is more conservative than that of rrnS. The presence of structural elements in the control region is also discussed, with emphasis on their implications in the regulation of replication and/or transcription of the reduviid mitogenome. The phylogenetic analyses indicated that within Reduviidae, Harpactorinae is a sister group to the Salyavatinae + Triatominae clade.

Transcriptome analysis in Concholepas concholepas (Gastropoda, Muricidae): mining and characterization of new genomic and molecular markers

The marine gastropod Concholepas concholepas, locally known as the "loco", is the main target species of the benthonic Chilean fisheries. Genetic and genomic tools are necessary to study the genome of this species in order to understand the molecular basis of its development, growth, and other key traits to improve the management strategies and to identify local adaptation to prevent loss of biodiversity. Here, we use pyrosequencing technologies to generate the first transcriptomic database from adult specimens of the loco. After trimming, a total of 140,756 Expressed Sequence Tag sequences were achieved. Clustering and assembly analysis identified 19,219 contigs and 105,435 singleton sequences. BlastN analysis showed a significant identity with Expressed Sequence Tags of different gastropod species available in public databases. Similarly, BlastX results showed that only 895 out of the total 124,654 had significant hits and may represent novel genes for marine gastropods. From this database, simple sequence repeat motifs were also identified and a total of 38 primer pairs were designed and tested to assess their potential as informative markers and to investigate their cross-species amplification in different related gastropod species. This dataset represents the first publicly available 454 data for a marine gastropod endemic to the southeastern Pacific coast, providing a valuable transcriptomic resource for future efforts of gene discovery and development of functional markers in other marine gastropods.

Transcriptome map of plant mitochondria reveals islands of unexpected transcribed regions

Background

Plant mitochondria contain a relatively large amount of genetic information, suggesting that their functional regulation may not be as straightforward as that of metazoans. We used a genomic tiling array to draw a transcriptomic atlas of Oryza sativa japonica (rice) mitochondria, which was predicted to be approximately 490-kb long.

Results

Whereas statistical analysis verified the transcription of all previously known functional genes such as the ones related to oxidative phosphorylation, a similar extent of RNA expression was frequently observed in the inter-genic regions where none of the previously annotated genes are located. The newly identified open reading frames (ORFs) predicted in these transcribed inter-genic regions were generally not conserved among flowering plant species, suggesting that these ORFs did not play a role in mitochondrial principal functions. We also identified two partial fragments of retrotransposon sequences as being transcribed in rice mitochondria.

Conclusion

The present study indicated the previously unexpected complexity of plant mitochondrial RNA metabolism. Our transcriptomic data (Oryza sativa Mitochondrial rna Expression Server: OsMES) is publicly accessible at [http://bioinf.mind.meiji.ac.jp/cgi-bin/gbrowse/OsMes/#search].

Evidence for a fourteenth mtDNA-encoded protein in the female-transmitted mtDNA of marine Mussels (Bivalvia: Mytilidae)

Background

A novel feature for animal mitochondrial genomes has been recently established: i.e., the presence of additional, lineage-specific, mtDNA-encoded proteins with functional significance. This feature has been observed in freshwater mussels with doubly uniparental inheritance of mtDNA (DUI). The latter unique system of mtDNA transmission, which also exists in some marine mussels and marine clams, is characterized by one mt genome inherited from the female parent (F mtDNA) and one mt genome inherited from the male parent (M mtDNA). In freshwater mussels, the novel mtDNA-encoded proteins have been shown to be mt genome-specific (i.e., one novel protein for F genomes and one novel protein for M genomes). It has been hypothesized that these novel, F- and M-specific, mtDNA-encoded proteins (and/or other F- and/or M-specific mtDNA sequences) could be responsible for the different modes of mtDNA transmission in bivalves but this remains to be demonstrated.

Methodology/Principal Findings

We investigated all complete (or nearly complete) female- and male-transmitted marine mussel mtDNAs previously sequenced for the presence of ORFs that could have functional importance in these bivalves. Our results confirm the presence of a novel F genome-specific mt ORF, of significant length (>100aa) and located in the control region, that most likely has functional significance in marine mussels. The identification of this ORF in five Mytilus species suggests that it has been maintained in the mytilid lineage (subfamily Mytilinae) for ∼13 million years. Furthermore, this ORF likely has a homologue in the F mt genome of Musculista senhousia, a DUI-containing mytilid species in the subfamily Crenellinae. We present evidence supporting the functionality of this F-specific ORF at the transcriptional, amino acid and nucleotide levels.

Conclusions/Significance

Our results offer support for the hypothesis that “novel F genome-specific mitochondrial genes” are involved in key biological functions in bivalve species with DUI.

Organellar RNA editing

RNA editing is a term used for a number of mechanistically different processes that alter the nucleotide sequence of RNA molecules to differ from the gene sequence. RNA editing occurs in a wide variety of organisms and is particularly frequent in organelle transcripts of eukaryotes. The discontiguous phylogenetic distribution of mRNA editing, the mechanistic differences observed in different organisms, and the nonhomologous editing machinery described in different taxonomic groups all suggest that RNA editing has appeared independently several times. This raises questions about the selection pressures acting to maintain editing that are yet to be completely resolved. Editing tends to be frequent in organisms with atypical organelle genomes and acts to correct the effect of DNA mutations that would otherwise compromise the synthesis of functional proteins. Additional functions of editing in generating protein diversity or regulating gene expression have been proposed but so far lack widespread experimental evidence, at least in organelles.

Mutation and selection on the wobble nucleotide in tRNA anticodons in marine bivalve mitochondrial genomes

Background

Animal mitochondrial genomes typically encode one tRNA for each synonymous codon family, so that each tRNA anticodon essentially has to wobble to recognize two or four synonymous codons. Several factors have been hypothesized to determine the nucleotide at the wobble site of a tRNA anticodon in mitochondrial genomes, such as the codon-anticodon adaptation hypothesis, the wobble versatility hypothesis, the translation initiation and elongation conflict hypothesis, and the wobble cost hypothesis.

Principal Findings

In this study, we analyzed codon usage and tRNA anticodon wobble sites of 29 marine bivalve mitochondrial genomes to evaluate features of the wobble nucleotides in tRNA anticodons. The strand-specific mutation bias favors G and T on the H strand in all the 29 marine bivalve mitochondrial genomes. A bias favoring G and T is also visible in the third codon positions of protein-coding genes and the wobble sites of anticodons, rejecting that codon usage bias drives the wobble sites of tRNA anticodons or tRNA anticodon bias drives the evolution of codon usage. Almost all codon families (98.9%) from marine bivalve mitogenomes support the wobble versatility hypothesis. There are a few interesting exceptions involving tRNA^Trp with an anticodon CCA fixed in Pectinoida species, tRNA^Ser with a GCU anticodon fixed in Mytiloida mitogenomes, and the uniform anticodon CAU of tRNA^Met translating the AUR codon family.

Conclusions/Significance

These results demonstrate that most of the nucleotides at the wobble sites of tRNA anticodons in marine bivalve mitogenomes are determined by wobble versatility. Other factors such as the translation initiation and elongation conflict, and the cost of wobble translation may contribute to the determination of the wobble nucleotide in tRNA anticodons. The finding presented here provides valuable insights into the previous hypotheses of the wobble nucleotide in tRNA anticodons by adding some new evidence.

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.

The mitochondrial genome of Colossendeis megalonyx supports a basal position of Colossendeidae within the Pycnogonida

We present the almost complete (16,007 bp) mitochondrial genome of a Colossendeis megalonyx specimen from the Southern Ocean and discuss gene order and tRNA structure in a comparative phylogenetic context. Our data suggest a basal position of the colossendeid lineage corroborating earlier phylogenetic studies but disagreeing with results of a recently published study that supported a highly derived sister-group relationship of Colossendeidae and Nymphonidae. Our results, together with BLAST searches and phylogenetic comparisons, indicate that the specimen presented as Colossendeis sp. in a series of recent studies had been misidentified. It has now been identified as a nymphonid species.

Novel protein genes in animal mtDNA: a new sex determination system in freshwater mussels (Bivalvia: Unionoida)?

Mitochondrial (mt) function depends critically on optimal interactions between components encoded by mt and nuclear DNAs. mitochondrial DNA (mtDNA) inheritance (SMI) is thought to have evolved in animal species to maintain mito-nuclear complementarity by preventing the spread of selfish mt elements thus typically rendering mtDNA heteroplasmy evolutionarily ephemeral. Here, we show that mtDNA intraorganismal heteroplasmy can have deterministic underpinnings and persist for hundreds of millions of years. We demonstrate that the only exception to SMI in the animal kingdom, that is, the doubly uniparental mtDNA inheritance system in bivalves, with its three-way interactions among egg mt-, sperm mt- and nucleus-encoded gene products, is tightly associated with the maintenance of separate male and female sexes (dioecy) in freshwater mussels. Specifically, this mother-through-daughter and father-through-son mtDNA inheritance system, containing highly differentiated mt genomes, is found in all dioecious freshwater mussel species. Conversely, all hermaphroditic species lack the paternally transmitted mtDNA (=possess SMI) and have heterogeneous macromutations in the recently discovered, novel protein-coding gene (F-orf) in their maternally transmitted mt genomes. Using immunoelectron microscopy, we have localized the F-open reading frame (ORF) protein, likely involved in specifying separate sexes, in mitochondria and in the nucleus. Our results support the hypothesis that proteins coded by the highly divergent maternally and paternally transmitted mt genomes could be directly involved in sex determination in freshwater mussels. Concomitantly, our study demonstrates novel features for animal mt genomes: the existence of additional, lineage-specific, mtDNA-encoded proteins with functional significance and the involvement of mtDNA-encoded proteins in extra-mt functions. Our results open new avenues for the identification, characterization, and functional analyses of ORFs in the intergenic regions, previously defined as "noncoding," found in a large proportion of animal mt genomes.

Rampant gene rearrangement and haplotype hypervariation among nematode mitochondrial genomes

Rare syntenic conservation, sequence duplication, and the use of both DNA strands to encode genes are signature architectural features defining mitochondrial genomes of enoplean nematodes. These characteristics stand in contrast to the more conserved mitochondrial genome sizes and transcriptional organizations of mitochondrial DNAs (mtDNAs) derived from chromadorean nematodes. To address the frequency of gene rearrangement within nematode mitochondrial DNA (mtDNA), mitochondrial genome variation has been characterized within a more confined enoplean taxonomic unit, the family Mermithidae. The complete nucleotide sequences of the mosquito parasitic nematodes Romanomermis culicivorax, R. nielseni, and R. iyengari mtDNA have been determined. Duplicated expanses encompassing different regions of the mitochondrial genomes were found in each of these congeners. These mtDNA shared few rRNA and protein gene junctions, indicating extensive gene rearrangement within the Romanomermis lineage. Rapid structural changes are also observed at the conspecific level where no two individual nematodes carry the same haplotype. Rolling circle amplification was used to isolate complete mitochondrial genomes from individuals in local populations of Thaumamermis cosgrovei, a parasite of terrestrial isopods. Mitochondrial DNA length variants ranging from 19 to 34 kb are observed, but haplotypes are not shared between any two individuals. The complete nucleotide sequences of three haplotypes have been determined, revealing a constant region encoding most mitochondrial genes and a hypervariable segment that contains intact and pseudogene copies of several mitochondrial genes, duplicated to different copy numbers, resulting in mtDNA size variation. Constant rearrangement generates new T. cosgrovei mtDNA forms.