Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

Phylomitogenomic analyses strongly support the sister relationship of the Chaetognatha and Protostomia

The phylum Chaetognatha (arrow worms) comprises a group of small marine predators that constitute a critical component of the zooplankton community throughout the world's oceans. Various phylogenetic affiliations have been proposed for the Chaetognatha, for which there are at least nine possible phylogenetic positions. Resolving the phylogenetic position of the chaetognaths is a key in understanding the fundamental developmental features of bilaterians. In comparison with the typical gene content of metazoan mitogenomes, two protein‐coding genes (atp6 and atp8) are absent from all chaetognaths. The two mitogenomes sequenced from Sagitta crassa and Zonosagitta nagae in this study nevertheless contain two and four tRNA genes, respectively, in contrast to those of the other five chaetognaths reported where only one tRNA gene (trnMet) is present, thus invalidating the view that all chaetognath mitogenomes have a single tRNA gene. A conserved major gene order shared by all chaetognaths could be partially identified in many protostome mitogenomes, but not in any ancestral mitogenome gene arrangement of the four deuterostome groups. Phylogenetic analysis of the deduced amino acid sequences of protein‐coding genes from 85 mitogenomes of 19 groups suggests the Chaetognatha to be a sister group to the protostomes, a result consistent with evidences from the developmental pattern and other molecular analyses.     

Complete mitochondrial genome of Membranipora grandicella (Bryozoa: Cheilostomatida) determined with next-generation sequencing: The first representative of the suborder Malacostegina

Next-generation sequencing (NGS) has proven a valuable platform for fast and easy obtaining of large numbers of sequences at relatively low cost. In this study we use shot-gun sequencing method on Illumina HiSeq 2000, to obtain enough sequences for the assembly of the bryozoan Membranipora grandicella (Bryozoa: Cheilostomatida) mitochondrial genome, which is the first representative of the suborder Malacostegina. The complete mitochondrial genome is 15,861 bp in length, which is relatively larger than other studied bryozoans. The mitochondrial genome contains 13 protein-coding genes, 2 ribosomal RNAs and 20 transfer RNAs. To investigate the phylogenetic position and the inner relationships of the phylum Bryozoa, phylogenetic trees were constructed with amino acid sequences of 11 PCGs from 30 metazoans. Two superclades of protostomes, namely Lophotrochozoa and Ecdysozoa, are recovered as monophyletic with strong support in both ML and Bayesian analyses. Somewhat to surprise, Bryozoa appears as the sister group of Chaetognatha with moderate or high support. The relationship among five bryozoans is Tubulipora flabellaris + (M. grandicella + (Flustrellidra hispida + (Bugula neritina + Watersipora subtorquata))), which supports for the view that Cheilostomatida is not a natural, monophyletic clade. NGS proved to be a quick and easy method for sequencing a complete mitochondrial genome.     

Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans

Using "long-PCR," we amplified in overlapping fragments the complete mitochondrial genome of the tapeworm Hymenolepis diminuta (Platyhelminthes: Cestoda) and determined its 13,900-nt sequence. The gene content is the same as that typically found for animal mitochondrial DNA (mtDNA) except that atp8 appears to be lacking, a condition found previously for several other animals. Despite the small size of this mtDNA, there are two large noncoding regions, one of which contains 13 repeats of a 31-nt sequence and a potential stem-loop structure of 25 bp with an 11-member loop. Large potential secondary structures were identified also for the noncoding regions of two other cestode mtDNAS: Comparison of the mitochondrial gene arrangement of H. diminuta with those previously published supports a phylogenetic position of flatworms as members of the Eutrochozoa, rather than placing them basal to either a clade of protostomes or a clade of coelomates.     

Complete mtDNA sequence of the North American freshwater mussel,Lampsilis ornata (Unionidae): an examination of the evolution and phylogenetic utility of mitochondrial genome organization in Bivalvia (Mollusca)

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. Two factors inhibiting our understanding the evolution of gene rearrangement in bivalves are inadequate taxonomic sampling and failure to examine gene order in a phylogenetic framework. Here, we report the first complete nucleotide sequence (16,060 bp) of the mitochondrial (mt) genome of a North American freshwater bivalve, Lampsilis ornata (Mollusca: Paleoheterodonta: Unionidae). Gene order and mt genome content is examined in a comparative phylogenetic framework for Lampsilis and five other bivalves, representing five families. Mitochondrial genome content is shown to vary by gene duplication and loss among taxa and between male and female mitotypes within a species. Although mt gene arrangement is highly variable among bivalves, when optimized on an independently derived phylogenetic hypothesis, it allows for the reconstruction of ancestral gene order states and indicates the potential phylogenetic utility of the data. However, the interpretation of reconstructed ancestral gene order states must take in to account both the accuracy of the phylogenetic estimation and the probability of character state change across the topology, such as the presence/absence of atp8 in bivalve lineages. We discuss what role, if any, doubly uniparental inheritance (DUI) and recombination between sexual mitotypes may play in influencing gene rearrangement of the mt genome in some bivalve lineages.     

Genome structure and gene content in protist mitochondrial DNAs.

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.     

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.     

Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails

Using DNA sequence data from multiple genes (often from more than one genome compartment) to reconstruct phylogenetic relationships has become routine. Augmenting this approach with genomic structural characters (e.g., intron gain and loss, changes in gene order) as these data become available from comparative studies already has provided critical insight into some long-standing questions about the evolution of land plants. Here we report on the presence of a group II intron located in the mitochondrial atp1 gene of leptosporangiate and marattioid ferns. Primary sequence data for the atp1 gene are newly reported for 27 taxa, and results are presented from maximum likelihood-based phylogenetic analyses using Bayesian inference for 34 land plants in three data sets: (1) single-gene mitochondrial atp1 (exon+intron sequences); (2) five combined genes (mitochondrial atp1 [exon only]; plastid rbcL, atpB, rps4; nuclear SSU rDNA); and (3) same five combined genes plus morphology. All our phylogenetic analyses corroborate results from previous fern studies that used plastid and nuclear sequence data: the monophyly of euphyllophytes, as well as of monilophytes; whisk ferns (Psilotidae) sister to ophioglossoid ferns (Ophioglossidae); horsetails (Equisetopsida) sister to marattioid ferns (Marattiidae), which together are sister to the monophyletic leptosporangiate ferns. In contrast to the results from the primary sequence data, the genomic structural data (atp1 intron distribution pattern) would seem to suggest that leptosporangiate and marattioid ferns are monophyletic, and together they are the sister group to horsetails--a topology that is rarely reconstructed using primary sequence data.     

Extreme mitochondrial evolution in the ctenophore Mnemiopsis leidyi: Insight from mtDNA and the nuclear genome

Recent advances in sequencing technology have led to a rapid accumulation of mitochondrial DNA (mtDNA) sequences, which now represent the wide spectrum of animal diversity. However, one animal phylum--Ctenophora--has, to date, remained completely unsampled. Ctenophores, a small group of marine animals, are of interest due to their unusual biology, controversial phylogenetic position, and devastating impact as invasive species. Using data from the Mnemiopsis leidyi genome sequencing project, we Polymerase Chain Reaction (PCR) amplified and analyzed its complete mitochondrial (mt-) genome. At just over 10 kb, the mt-genome of M. leidyi is the smallest animal mtDNA ever reported and is among the most derived. It has lost at least 25 genes, including atp6 and all tRNA genes. We show that atp6 has been relocated to the nuclear genome and has acquired introns and a mitochondrial targeting presequence, while tRNA genes have been genuinely lost, along with nuclear-encoded mt-aminoacyl tRNA synthetases. The mt-genome of M. leidyi also displays extremely high rates of sequence evolution, which likely led to the degeneration of both protein and rRNA genes. In particular, encoded rRNA molecules possess little similarity with their homologs in other organisms and have highly reduced secondary structures. At the same time, nuclear encoded mt-ribosomal proteins have undergone expansions, likely to compensate for the reductions in mt-rRNA. The unusual features identified in M. leidyi mtDNA make this organism an interesting system for the study of various aspects of mitochondrial biology, particularly protein and tRNA import and mt-ribosome structures, and add to its value as an emerging model species. Furthermore, the fast-evolving M. leidyi mtDNA should be a convenient molecular marker for species- and population-level studies.     

Mitochondrial Genome of <Emphasis Type="Italic">Savalia savaglia</Emphasis> (Cnidaria,Hexacorallia) and Early Metazoan Phylogeny

Mitochondrial genomes have recently become widely used in animal phylogeny, mainly to infer the relationships between vertebrates and other bilaterians. However, only 11 of 723 complete mitochondrial genomes available in the public databases are of early metazoans, including cnidarians (Anthozoa, mainly Scleractinia) and sponges. Although some cnidarians (Medusozoa) are known to possess atypical linear mitochondrial DNA, the anthozoan mitochondrial genome is circular and its organization is similar to that of other metazoans. Because the phylogenetic relationships among Anthozoa as well as their relation to other early metazoans still need to be clarified, we tested whether sequencing the complete mitochondrial genome of Savalia savaglia, an anthozoan belonging to the order Zoantharia (=Zoanthidea), could be useful to infer such relationships. Compared to other anthozoans, S. savaglia’s genome is unusually long (20,766 bp) due to the presence of several noncoding intergenic regions (3691 bp). The genome contains all 13 protein coding genes commonly found in metazoans, but like other Anthozoa it lacks most of the tRNAs. Phylogenetic analyses of S. savaglia mitochondrial sequences show Zoantharia branching closely to other Hexacorallia, either as a sister group to Actiniaria or as a sister group to Actiniaria and Scleractinia. The close relationships suggested between Zoantharia and Actiniaria are reinforced by strong similarities in their gene order and the presence of similar introns in the COI and ND5 genes. Our study suggests that mitochondrial genomes can be a source of potentially valuable information on the phylogeny of Hexacorallia and may provide new insights into the evolution of early metazoans. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Axel Meyer]     

Use of atp6 in fungal phylogenetics: an example from the boletales

Complete nucleotide sequences have been determined for atp6 from Suillus luteus and cox3 from Suillus sinuspaulianus (Boletales, Hymenomycetes, Basidiomycota), which code for ATPase subunit 6 and cytochrome oxidase subunit 3, respectively. These sequences were used to design PCR primers for the amplification of partial atp6 and cox3 sequences from other members of the Boletales and outgroup taxa. In atp6 and cox3 from Russula rosacea, one of the outgroup taxa, we observed a number of in-frame TGA(trp) codons, which imply a Neurospora crassa-type mitochondrial code in R. rosacea and possibly in basidiomycetes in general. Interestingly, however, most basidiomycetes other than R. rosacea appear to strongly prefer the TGG(trp) codon, which is unusual, given the strong A + T bias in fungal mitochondrial genomes. Pairwise comparisons were performed between atp6 sequences from increasingly divergent fungal lineages, and results show that all three codon positions become saturated in substitutions after an estimated divergence time of approx 300 Ma. This means that atp6 is likely to provide phylogenetic resolution within fungal classes but not at higher taxonomic levels. Also, because of the strong A + T bias in fungal mitochondrial genomes, A/T transversions were found to be more common than any other type of substitution, resulting in transversions being about two to three times more common in most pairwise sequence comparisons. Finally, atp6 sequences were used to infer phylogenetic relationships between 27 taxa from the Boletales and 4 outgroup taxa. Analyses were performed (i) on nucleotide sequence data using parsimony (successive approximation) as well as maximum likelihood methods and (ii) on deduced amino acid sequences using distance methods based on empirical substitution probabilities. Results from the various analyses are largely concordant with each other as well as with prior analyses of partial mitochondrial large-subunit rDNA (mtLSU rDNA). Analysis of the combined atp6 and mtLSU rDNA sequences results in increased bootstrap support for several key branches. Relationships that have been resolved for the first time in the current analysis are discussed.     

The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

The partial mitochondrial genome of the Cephalothrix rufifrons (Nemertea, Palaeonemertea): characterization and implications for the phylogenetic position of Nemertea

A continuous 10.1kb fragment of the Cephalothrix rufifrons (Nemertea, Palaeonemertea) mitochondrial genome was sequenced and characterized to further assess organization of protostome mitochondrial genomes and evaluate the phylogenetic potential of gene arrangement and amino acid characters. The genome is A-T rich (72%), and this biased base composition is partly reflected in codon usage. Inferred tRNA secondary structures are typical of those reported for other metazoan mitochondrial DNAs. The arrangement of the 26 genes contained in the fragment exhibits marked similarity to those of many protostome taxa, most notably molluscs with highly conserved arrangements and a phoronid. Separate and simultaneous phylogenetic analyses of inferred amino acid sequences and gene adjacencies place the nemertean within the protostomes among coelomate lophotrochozoan taxa, but do not find a well-supported sister taxon link.     

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

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

Complete Mitochondrial Genome of Haplorchis taichui and Comparative Analysis with Other Trematodes

Mitochondrial genomes have been extensively studied for phylogenetic purposes and to investigate intra- and interspecific genetic variations. In recent years, numerous groups have undertaken sequencing of platyhelminth mitochondrial genomes. Haplorchis taichui (family Heterophyidae) is a trematode that infects humans and animals mainly in Asia, including the Mekong River basin. We sequenced and determined the organization of the complete mitochondrial genome of H. taichui. The mitochondrial genome is 15,130 bp long, containing 12 protein-coding genes, 2 ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). Like other trematodes, it does not encode the atp8 gene. All genes are transcribed from the same strand. The ATG initiation codon is used for 9 protein-coding genes, and GTG for the remaining 3 (nad1, nad4, and nad5). The mitochondrial genome of H. taichui has a single long non-coding region between trnE and trnG. H. taichui has evolved as being more closely related to Opisthorchiidae than other trematode groups with maximal support in the phylogenetic analysis. Our results could provide a resource for the comparative mitochondrial genome analysis of trematodes, and may yield genetic markers for molecular epidemiological investigations into intestinal flukes.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Mitochondrial diversity of early-branching metazoa is revealed by the complete mt genome of a haplosclerid demosponge

The first mitochondrial (mt) genomes of demosponges have recently been sequenced and appear to be markedly different from published eumetazoan mt genomes. Here we show that the mt genome of the haplosclerid demosponge Amphimedon queenslandica has features that it shares with both demosponges and eumetazoans. Although the A. queenslandica mt genome has typical demosponge features, including size, long noncoding regions, and bacterialike rRNA genes, it lacks atp9, which is found in the other demosponges sequenced to date. We found strong evidence of a recent transposon-mediated transfer of atp9 to the nuclear genome. In addition, A. queenslandica bears an incomplete tRNA set, unusual amino acid deletion patterns, and a putative control region. Furthermore, the arrangement of mt rRNA genes differs from that of other demosponges. These genes evolve at significantly higher rates than observed in other demosponges, similar to previously observed nuclear rRNA gene rates in other haplosclerid demosponges.