Complete mtDNA of Ciona intestinalis reveals extensive gene rearrangement and the presence of an atp8 and an extra trnM gene in ascidians

Evolutionary Insights from the Mitochondrial Genome of Oikopleura dioica: Sequencing Challenges, RNA Editing, Gene Transfers to the Nucleus, and tRNA Loss

Sequencing the mitochondrial genome of the tunicate Oikopleura dioica is a challenging task due to the presence of long poly-A/T homopolymer stretches, which impair sequencing and assembly. Here, we report on the sequencing and annotation of the majority of the mitochondrial genome of O. dioica by means of combining several DNA and amplicon reads obtained by Illumina and MinIon Oxford Nanopore Technologies with public RNA sequences. We document extensive RNA editing, since all homopolymer stretches present in the mitochondrial DNA correspond to 6U-regions in the mitochondrial RNA. Out of the 13 canonical protein-coding genes, we were able to detect eight, plus an unassigned open reading frame that lacked sequence similarity to canonical mitochondrial protein-coding genes. We show that the nad3 gene has been transferred to the nucleus and acquired a mitochondria-targeting signal. In addition to two very short rRNAs, we could only identify a single tRNA (tRNA-Met), suggesting multiple losses of tRNA genes, supported by a corresponding loss of mitochondrial aminoacyl-tRNA synthetases in the nuclear genome. Based on the eight canonical protein-coding genes identified, we reconstructed maximum likelihood and Bayesian phylogenetic trees and inferred an extreme evolutionary rate of this mitochondrial genome. The phylogenetic position of appendicularians among tunicates, however, could not be accurately determined.

Mitochondrial genomic analyses provide new insights into the "missing" atp8 and adaptive evolution of Mytilidae

Background

Mytilidae, also known as marine mussels, are widely distributed in the oceans worldwide. Members of Mytilidae show a tremendous range of ecological adaptions, from the species distributed in freshwater to those that inhabit in deep-sea. Mitochondria play an important role in energy metabolism, which might contribute to the adaptation of Mytilidae to different environments. In addition, some bivalve species are thought to lack the mitochondrial protein-coding gene ATP synthase F0 subunit 8. Increasing studies indicated that the absence of atp8 may be caused by annotation difficulties for atp8 gene is characterized by highly divergent, variable length.

Results

In this study, the complete mitochondrial genomes of three marine mussels (Xenostrobus securis, Bathymodiolus puteoserpentis, Gigantidas vrijenhoeki) were newly assembled, with the lengths of 14,972 bp, 20,482, and 17,786 bp, respectively. We annotated atp8 in the sequences that we assembled and the sequences lacking atp8. The newly annotated atp8 sequences all have one predicted transmembrane domain, a similar hydropathy profile, as well as the C-terminal region with positively charged amino acids. Furthermore, we reconstructed the phylogenetic trees and performed positive selection analysis. The results showed that the deep-sea bathymodiolines experienced more relaxed evolutionary constraints. And signatures of positive selection were detected in nad4 of Limnoperna fortunei, which may contribute to the survival and/or thriving of this species in freshwater.

Conclusions

Our analysis supported that atp8 may not be missing in the Mytilidae. And our results provided evidence that the mitochondrial genes may contribute to the adaptation of Mytilidae to different environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08940-8.

Comparative mitogenomic analysis of the superfamily Tellinoidea (Mollusca: Bivalvia): Insights into the evolution of the gene rearrangements

The superfamily Tellinoidea is widespread and contains approximately 180 living species, which is one of the most diverse and representative groups among the bivalves. In order to extend our knowledge on evolution of tellinoidean species, we newly determined five tellinoidean mitochondrial genomes (mitogenomes). The newly determined mitogenome vary in size from 16,333 to 16,986 bp. The results show that the genome size and genome organization are conserved in tellinoideans. However, gene arrangement and the location of the major non-coding region (NCR) show diversity. The atp8 gene presents in all the five new mitogenomes. Two trnK and trnP genes were detected in Gari togata mitogenome. Phylogenetic analysis supports the monophyly of Tellinoidea, however, it's family Psammobiidae is polyphyletic. CREx analysis suggests that the gene order of Nuttallia olivacea is assumed as the most primitive condition of Tellinoidea. We map the gene order onto the phylogeny and infers the possible gene rearrangement scenarios among tellinoidean mitogenomes. The mitochondrial gene rearrangement is a useful information that help reassessing the phylogeny of Tellinoidea. Phylogenetic relationship and gene arrangement analyses suggest that a careful review for the current taxonomy of the family Psammobiidae is required.

Identification and population genetic comparison of three ascidian species based on mtDNA sequences

Ascidians are sessile marine chordate invertebrates found along seashores worldwide and are typically regarded as invasive organisms. Knowledge concerning their global genetic structure and subsequent invasive potential is limited. Here, we identified three ascidians-Ciona robusta, Ciona savignyi, and Styela clava from the northeast region of China using morphological characteristics and mitochondrial cytochrome c oxidase subunit I (cox1) as genetic marker. We additionally used phylogenetics to aid in the identification of these three species. The results of a population genetic analysis showed that among the three species, the level of haplotype diversity was particularly high within C. savignyi, and nucleotide diversity varied moderately. We divided the three species separately into native and invasive populations using 170 cox1 sequences from global resources to explore population genetic structure and invasive potential. Although in the network analysis Ciona spp. formed haplogroups of native and invasive populations, some haplotypes were still shared. We found that the haplotypes did not cluster within the network of S. clava. Our AMOVA results also showed that Ciona spp. had a weak genetic structure, and less genetic differentiation was present in S. clava. These data suggest that there are extensive incursions of these three ascidians into different geographical regions. Global comparisons of ascidian populations will help in the understanding of their population genetic structure and invasive potential, hence providing important insights regarding conservation as well as management.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

Extensive mitochondrial gene rearrangements in Ctenophora: insights from benthic Platyctenida

BACKGROUND

Complete mitochondrial (mt) genomes have been sequenced for thousands of animals and represent a molecule of choice for many evolutionary studies. Nevertheless, some animal groups have remained under-sampled. Ctenophora (comb jellies) is one such example, with only two complete mt sequences determined hitherto for this phylum, which encompasses ca. 150-200 described species. This lack of data derives from the extremely fast mt evolutionary rate in this lineage, complicating primer design and DNA amplification. Indeed, in the two ctenophore mt genomes sequenced to date, i.e. those of Mnemiopsis leidyi (order Lobata) and Pleurobrachia bachei (order Cydippida), both rRNA and protein coding genes exhibit an extraordinary size reduction and have highly derived sequences. Additionally, all tRNAs, and the atp6 and atp8 genes are absent. In order to determine whether these characteristics are shared by other ctenophores, we obtained the complete mt genomes of three benthic ctenophores belonging to the so far unsampled order of Platyctenida: Coeloplana loyai, Coeloplana yulianicorum and Vallicula multiformis.

RESULTS

The mt genomes of benthic ctenophores reveal the same peculiarities found in Mnemiopsis and Pleurobrachia, demonstrating that the fast evolutionary rate is a general trait of the ctenophore mt genomes. Our results also indicate that this high evolutionary rate not only affects the nucleotide substitution but also gene rearrangements. Indeed, gene order was highly rearranged among representatives of the different taxonomic orders in which it was close to random, but also quite variable within Platyctenida, in which the genera Coeloplana and Vallicula share only four conserved synteny blocks. However, the two congeneric Coeloplana species display exactly the same gene order. Because of the extreme evolutionary rate, our phylogenetic analyses were unable to resolve the phylogenetic position of ctenophores within metazoans or the relationships among the different Ctenophora orders. Comparative sequence-analyses allowed us to correct the annotation of the Pleurobrachia mt genome, confirming the absence of tRNAs, the presence of both rRNA genes, and the existence of a reassignment of codon TGA from tryptophan to serine for this species.

CONCLUSIONS

Since Platyctenida is an early diverging lineage among Ctenophora, our findings suggest that the mt traits described above are ancestral characteristics of this phylum.

Digging deeper: new gene order rearrangements and distinct patterns of codons usage in mitochondrial genomes among shrimps from the Axiidea, Gebiidea and Caridea (Crustacea: Decapoda)

BACKGROUND

Whole mitochondrial DNA is being increasingly utilized for comparative genomic and phylogenetic studies at deep and shallow evolutionary levels for a range of taxonomic groups. Although mitogenome sequences are deposited at an increasing rate into public databases, their taxonomic representation is unequal across major taxonomic groups. In the case of decapod crustaceans, several infraorders, including Axiidea (ghost shrimps, sponge shrimps, and mud lobsters) and Caridea (true shrimps) are still under-represented, limiting comprehensive phylogenetic studies that utilize mitogenomic information.

METHODS

Sequence reads from partial genome scans were generated using the Illumina MiSeq platform and mitogenome sequences were assembled from these low coverage reads. In addition to examining phylogenetic relationships within the three infraorders, Axiidea, Gebiidea, and Caridea, we also investigated the diversity and frequency of codon usage bias and mitogenome gene order rearrangements.

RESULTS

We present new mitogenome sequences for five shrimp species from Australia that includes two ghost shrimps, Callianassa ceramica and Trypaea australiensis, along with three caridean shrimps, Macrobrachium bullatum, Alpheus lobidens, and Caridina cf. nilotica. Strong differences in codon usage were discovered among the three infraorders and significant gene order rearrangements were observed. While the gene order rearrangements are congruent with the inferred phylogenetic relationships and consistent with taxonomic classification, they are unevenly distributed within and among the three infraorders.

DISCUSSION

Our findings suggest potential for mitogenome rearrangements to be useful phylogenetic markers for decapod crustaceans and at the same time raise important questions concerning the drivers of mitogenome evolution in different decapod crustacean lineages.

Mitochondrial Genomes of Kinorhyncha: trnM Duplication and New Gene Orders within Animals

Many features of mitochondrial genomes of animals, such as patterns of gene arrangement, nucleotide content and substitution rate variation are extensively used in evolutionary and phylogenetic studies. Nearly 6,000 mitochondrial genomes of animals have already been sequenced, covering the majority of animal phyla. One of the groups that escaped mitogenome sequencing is phylum Kinorhyncha-an isolated taxon of microscopic worm-like ecdysozoans. The kinorhynchs are thought to be one of the early-branching lineages of Ecdysozoa, and their mitochondrial genomes may be important for resolving evolutionary relations between major animal taxa. Here we present the results of sequencing and analysis of mitochondrial genomes from two members of Kinorhyncha, Echinoderes svetlanae (Cyclorhagida) and Pycnophyes kielensis (Allomalorhagida). Their mitochondrial genomes are circular molecules approximately 15 Kbp in size. The kinorhynch mitochondrial gene sequences are highly divergent, which precludes accurate phylogenetic inference. The mitogenomes of both species encode a typical metazoan complement of 37 genes, which are all positioned on the major strand, but the gene order is distinct and unique among Ecdysozoa or animals as a whole. We predict four types of start codons for protein-coding genes in E. svetlanae and five in P. kielensis with a consensus DTD in single letter code. The mitochondrial genomes of E. svetlanae and P. kielensis encode duplicated methionine tRNA genes that display compensatory nucleotide substitutions. Two distant species of Kinorhyncha demonstrate similar patterns of gene arrangements in their mitogenomes. Both genomes have duplicated methionine tRNA genes; the duplication predates the divergence of two species. The kinorhynchs share a few features pertaining to gene order that align them with Priapulida. Gene order analysis reveals that gene arrangement specific of Priapulida may be ancestral for Scalidophora, Ecdysozoa, and even Protostomia.

Mitochondrial gene order variation in the brachiopod Lingula anatina and its implications for mitochondrial evolution in lophotrochozoans

Vertebrate mitochondrial (mt) genomes display highly conserved gene order and relatively low evolutionary rates. However, these features are variable in marine invertebrates. Here we present the mt genome of the lingulid brachiopod, Lingula anatina, from Amami Island, Japan, as part of the nuclear genome project. We obtain ~2000-fold coverage of the 17.9-kb mt genome using Illumina sequencing, and we identify hypervariable regions within the same individual. Transcriptome analyses show that mt transcripts are polycistronic and expressed differentially. Unexpectedly, we find that the mt gene order of Amami Lingula is completely shuffled compared to that of a specimen from Yanagawa, suggesting that there may be cryptic species. Using breakpoint distance analyses with 101 metazoan mt genomes, we show that the evolutionary history of mt gene order among lophotrochozoans is unique. Analyses of non-synonymous substitution rates reveal that mt protein-coding genes of Lingula have experienced rapid evolution comparable to that expected for interspecific comparisons. Whole genome phylogenetic analyses suggest that mt genomes have limited value for inferring the phylogenetic positions of lophotrochozoans because of their high evolutionary rates in brachiopods and bivalves.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities

BACKGROUND

Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform.

RESULTS

The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations.

CONCLUSIONS

Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

Ascidian mitogenomics: comparison of evolutionary rates in closely related taxa provides evidence of ongoing speciation events

Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (dN) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri, where most pairs show a dN value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1. Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

BACKGROUND

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

RESULTS

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

CONCLUSIONS

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Host control of symbiont natural product chemistry in cryptic populations of the tunicate Lissoclinum patella

Natural products (secondary metabolites) found in marine invertebrates are often thought to be produced by resident symbiotic bacteria, and these products appear to play a major role in the symbiotic interaction of bacteria and their hosts. In these animals, there is extensive variation, both in chemistry and in the symbiotic bacteria that produce them. Here, we sought to answer the question of what factors underlie chemical variation in the ocean. As a model, we investigated the colonial tunicate Lissoclinum patella because of its rich and varied chemistry and its broad geographic range. We sequenced mitochondrial cytochrome c oxidase 1 (COXI) genes, and found that animals classified as L. patella fall into three phylogenetic groups that may encompass several cryptic species. The presence of individual natural products followed the phylogenetic relationship of the host animals, even though the compounds are produced by symbiotic bacteria that do not follow host phylogeny. In sum, we show that cryptic populations of animals underlie the observed chemical diversity, suggesting that the host controls selection for particular secondary metabolite pathways. These results imply novel approaches to obtain chemical diversity from the oceans, and also demonstrate that the diversity of marine natural products may be greatly impacted by cryptic local extinctions.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.

Tunicate mitogenomics and phylogenetics: peculiarities of the Herdmania momus mitochondrial genome and support for the new chordate phylogeny

Background

Tunicates represent a key metazoan group as the sister-group of vertebrates within chordates. The six complete mitochondrial genomes available so far for tunicates have revealed distinctive features. Extensive gene rearrangements and particularly high evolutionary rates have been evidenced with regard to other chordates. This peculiar evolutionary dynamics has hampered the reconstruction of tunicate phylogenetic relationships within chordates based on mitogenomic data.

Results

In order to further understand the atypical evolutionary dynamics of the mitochondrial genome of tunicates, we determined the complete sequence of the solitary ascidian Herdmania momus. This genome from a stolidobranch ascidian presents the typical tunicate gene content with 13 protein-coding genes, 2 rRNAs and 24 tRNAs which are all encoded on the same strand. However, it also presents a novel gene arrangement, highlighting the extreme plasticity of gene order observed in tunicate mitochondrial genomes. Probabilistic phylogenetic inferences were conducted on the concatenation of the 13 mitochondrial protein-coding genes from representatives of major metazoan phyla. We show that whereas standard homogeneous amino acid models support an artefactual sister position of tunicates relative to all other bilaterians, the CAT and CAT+BP site- and time-heterogeneous mixture models place tunicates as the sister-group of vertebrates within monophyletic chordates. Moreover, the reference phylogeny indicates that tunicate mitochondrial genomes have experienced a drastic acceleration in their evolutionary rate that equally affects protein-coding and ribosomal-RNA genes.

Conclusion

This is the first mitogenomic study supporting the new chordate phylogeny revealed by recent phylogenomic analyses. It illustrates the beneficial effects of an increased taxon sampling coupled with the use of more realistic amino acid substitution models for the reconstruction of animal phylogeny.

The mitochondrial genome structure of Xenoturbella bocki (phylum Xenoturbellida) is ancestral within the deuterostomes

Background

Mitochondrial genome comparisons contribute in multiple ways when inferring animal relationships. As well as primary sequence data, rare genomic changes such as gene order, shared gene boundaries and genetic code changes, which are unlikely to have arisen through convergent evolution, are useful tools in resolving deep phylogenies. Xenoturbella bocki is a morphologically simple benthic marine worm recently found to belong among the deuterostomes. Here we present analyses comparing the Xenoturbella bocki mitochondrial gene order, genetic code and control region to those of other metazoan groups.

Results

The complete mitochondrial genome sequence of Xenoturbella bocki was determined. The gene order is most similar to that of the chordates and the hemichordates, indicating that this conserved mitochondrial gene order might be ancestral to the deuterostome clade. Using data from all phyla of deuterostomes, we infer the ancestral mitochondrial gene order for this clade. Using inversion and breakpoint analyses of metazoan mitochondrial genomes, we test conflicting hypotheses for the phylogenetic placement of Xenoturbella and find a closer affinity to the hemichordates than to other metazoan groups. Comparative analyses of the control region reveal similarities in the transcription initiation and termination sites and origin of replication of Xenoturbella with those of the vertebrates. Phylogenetic analyses of the mitochondrial sequence indicate a weakly supported placement as a basal deuterostome, a result that may be the effect of compositional bias.

Conclusion

The mitochondrial genome of Xenoturbella bocki has a very conserved gene arrangement in the deuterostome group, strikingly similar to that of the hemichordates and the chordates, and thus to the ancestral deuterostome gene order. Similarity to the hemichordates in particular is suggested by inversion and breakpoint analysis. Finally, while phylogenetic analyses of the mitochondrial sequences support a basal deuterostome placement, support for this decreases with the use of more sophisticated models of sequence evolution.

Comparative mitogenomic analyses of three scallops (Bivalvia: Pectinidae) reveal high level variation of genomic organization and a diversity of transfer RNA gene sets

Background

It can be seen from the available mollusk mitogenomes that the family Pectinidae exhibits the most variation in genome organization. In this study, comparative mitogenomic analyses were performed for three scallops from the subfamily Chlamydinae (Pectinidae), with the goal of characterizing the degree of variability of mitogenome organization and other characteristics among species from the same subfamily and exploring their possible evolution route.

Findings

The complete or nearly complete mtDNA sequences of scallop Mimachlamys nobilis (17 935 bp), Mizuhopecten yessoensis (20 964 bp) and Chlamys farreri (17 035 bp) were determined using long PCR amplification and primer walking sequencing strategy. Highly variable size difference of the three genomes resulted primarily from length and number variations of non-coding regions, and the major difference in gene content of the three scallop species are due to varying tRNA gene sets. Only 21, 16, and 17 tRNA genes were detected in the mitogenomes of M. nobilis, M. yessoensis and C. farreri, respectively. Remarkably, no trnS gene could be identified in any of the three scallops. A newly-detected trnA-like sequence within the mitogenome of M. yessoensis seems to exemplify the functional loss of a tRNA gene, and the duplication of trnD in M. yessoensis raises a fundamental question of whether the retention of the tRNA gene copy of 2-tRNAs is easier than that of 4-tRNAs. Analysis of putative evolutionary pathways of gene rearrangement indicates that transposition of neighboring gene blocks may play an important role in the evolution of mitogenomes in scallops. Parsimonious analysis of the genomic variations implies that the mitogenomes of M. yessoensis and C. farreri are likely to derive independently from a common ancestor that was closely related to M. nobilis.

Conclusion

Comparative mitogenomic analyses among three species from the subfamily Chlamydinae show that the three genomes exhibit a high level of genomic variation and a diversity of tRNA gene sets, characterized by extensive translocation of genes. These features provide useful clues and information for evolutionary analysis of scallop mitogenomes.

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

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The mitochondrial genome of Phallusia mammillata and Phallusia fumigata (Tunicata, Ascidiacea): high genome plasticity at intra-genus level

Background

Within Chordata, the subphyla Vertebrata and Cephalochordata (lancelets) are characterized by a remarkable stability of the mitochondrial (mt) genome, with constancy of gene content and almost invariant gene order, whereas the limited mitochondrial data on the subphylum Tunicata suggest frequent and extensive gene rearrangements, observed also within ascidians of the same genus.

Results

To confirm this evolutionary trend and to better understand the evolutionary dynamics of the mitochondrial genome in Tunicata Ascidiacea, we have sequenced and characterized the complete mt genome of two congeneric ascidian species, Phallusia mammillata and Phallusia fumigata (Phlebobranchiata, Ascidiidae). The two mtDNAs are surprisingly rearranged, both with respect to one another and relative to those of other tunicates and chordates, with gene rearrangements affecting both protein-coding and tRNA genes. The new data highlight the extraordinary variability of ascidian mt genome in base composition, tRNA secondary structure, tRNA gene content, and non-coding regions (number, size, sequence and location). Indeed, both Phallusia genomes lack the trnD gene, show loss/acquisition of DHU-arm in two tRNAs, and have a G+C content two-fold higher than other ascidians. Moreover, the mt genome of P. fumigata presents two identical copies of trnI, an extra tRNA gene with uncertain amino acid specificity, and four almost identical sequence regions. In addition, a truncated cytochrome b, lacking a C-terminal tail that commonly protrudes into the mt matrix, has been identified as a new mt feature probably shared by all tunicates.

Conclusion

The frequent occurrence of major gene order rearrangements in ascidians both at high taxonomic level and within the same genus makes this taxon an excellent model to study the mechanisms of gene rearrangement, and renders the mt genome an invaluable phylogenetic marker to investigate molecular biodiversity and speciation events in this largely unexplored group of basal chordates.

Mitogenomics reveals two cryptic species in Ciona intestinalis

Individual mitochondrial genes or genomic features are commonly used as phylogenetic markers at many taxonomic levels. We used a mitogenomics approach to demonstrate the existence of two cryptic species in the ascidian Ciona intestinalis, a model chordate whose status as a single species has recently been questioned. Comprehensive comparative analysis of the mitochondrial genome of the two cryptic species revealed significant differences in gene order, size and number of noncoding regions, compositional features and divergence of protein-coding genes.

Initial stage of genetic mapping inCiona intestinalis

The use of classic genetics is emerging in the ascidian Ciona intestinalis; recent advances in genomics and high-quality developmental and evolutionary studies have made this animal an attractive model for research purposes. Genetic mapping in Ciona will likely make a major contribution to ascidian genomics and developmental biology by providing support for genome assembly and annotation and for the isolation of genes with particular mutations, while construction of genetic maps advances classic genetics in this species. Two major issues must be overcome before fine genetic maps can be constructed: the choice of proper genetic backgrounds and the establishment of laboratory strains. A high degree of polymorphism is useful for genetic mapping if we consider particular combinations of genetic backgrounds and techniques, although it is necessary to pay attention to the confused classification of C. intestinalis. Thus, it is preferred to establish laboratory strains instead of using samples with various genetic backgrounds. As these issues are unresolved, only amplified fragment length polymorphism-based maps have been created, while bulk segregant analysis is expected to isolate markers flanking mutant loci. However, rich genomic resources should facilitate the next stage of genetic map construction based on type I markers using coding sequences. The meiotic events that occur in crossing experiments for mapping purposes should shed light on population genetics and speciation issues. The results of such investigations may provide feedback for comparative genomics and developmental genetics in the near future.

Ribosomal RNA genes and deuterostome phylogeny revisited: More cyclostomes, elasmobranchs, reptiles, and a brittle star

This is an expanded study of the relationships among the deuterostome animals based on combined, nearly complete 28S and 18S rRNA genes (>3925 nt.). It adds sequences from 20 more taxa to the approximately 45 sequences used in past studies. Seven of the new taxa were sequenced here (brittle star Ophiomyxa, lizard Anolis, turtle Chrysemys, sixgill shark Hexanchus, electric ray Narcine, Southern Hemisphere lamprey Geotria, and Atlantic hagfish Myxine for 28S), and the other 13 were from GenBank and the literature (from a chicken, dog, rat, human, three lungfishes, and several ray-finned fishes, or Actinopterygii). As before, our alignments were based on secondary structure but did not account for base pairing in the stems of rRNA. The new findings, derived from likelihood-based tree-reconstruction methods and by testing hypotheses with parametric bootstrapping, include: (1) brittle star joins with sea star in the echinoderm clade, Asterozoa; (2) with two hagfishes and two lampreys now available, the cyclostome (jawless) fishes remain monophyletic; (3) Hexanchiform sharks are monophyletic, as Hexanchus groups with the frilled shark, Chlamydoselachus; (4) turtle is the sister taxon of all other amniotes; (5) bird is closer to the lizard than to the mammals; (6) the bichir Polypterus is in a monophyletic Actinopterygii; (7) Zebrafish Danio is the sister taxon of the other two teleosts we examined (trout and perch); (8) the South American and African lungfishes group together to the exclusion of the Australian lungfish. Other findings either upheld those of the previous rRNA-based studies (e.g., echinoderms and hemichordates group as Ambulacraria; orbitostylic sharks; batoids are not derived from any living lineage of sharks) or were obvious (monophyly of mammals, gnathostomes, vertebrates, echinoderms, etc.). Despite all these findings, the rRNA data still fail to resolve the relations among the major groups of deuterostomes (tunicates, Ambulacraria, cephalochordates and vertebrates) and of gnathostomes (chondrichthyans, lungfishes, coelacanth, actinopterygians, amphibians, and amniotes), partly because tunicates and lungfishes are rogue taxa that disrupt the tree. Nonetheless, parametric bootstrapping showed our RNA-gene data are only consistent with these dominant hypotheses: (1) deuterostomes consist of Ambulacraria plus Chordata, with Chordata consisting of tunicates and 'vertebrates plus cephalochordates'; and (2) lungfishes are the closest living relatives of tetrapods.

The mechanisms of codon reassignments in mitochondrial genetic codes

Many cases of nonstandard genetic codes are known in mitochondrial genomes. We carry out analysis of phylogeny and codon usage of organisms for which the complete mitochondrial genome is available, and we determine the most likely mechanism for codon reassignment in each case. Reassignment events can be classified according to the gain-loss framework. The “gain” represents the appearance of a new tRNA for the reassigned codon or the change of an existing tRNA such that it gains the ability to pair with the codon. The “loss” represents the deletion of a tRNA or the change in a tRNA so that it no longer translates the codon. One possible mechanism is codon disappearance (CD), where the codon disappears from the genome prior to the gain and loss events. In the alternative mechanisms the codon does not disappear. In the unassigned codon mechanism, the loss occurs first, whereas in the ambiguous intermediate mechanism, the gain occurs first. Codon usage analysis gives clear evidence of cases where the codon disappeared at the point of the reassignment and also cases where it did not disappear. CD is the probable explanation for stop to sense reassignments and a small number of reassignments of sense codons. However, the majority of sense-to-sense reassignments cannot be explained by CD. In the latter cases, by analysis of the presence or absence of tRNAs in the genome and of the changes in tRNA sequences, it is sometimes possible to distinguish between the unassigned codon and the ambiguous intermediate mechanisms. We emphasize that not all reassignments follow the same scenario and that it is necessary to consider the details of each case carefully.

Conflict between translation initiation and elongation in vertebrate mitochondrial genomes

The strand-biased mutation spectrum in vertebrate mitochondrial genomes results in an AC-rich L-strand and a GT-rich H-strand. Because the L-strand is the sense strand of 12 protein-coding genes out of the 13, the third codon position is overall strongly AC-biased. The wobble site of the anticodon of the 22 mitochondrial tRNAs is either U or G to pair with the most abundant synonymous codon, with only one exception. The wobble site of Met-tRNA is C instead of U, forming the Watson-Crick match with AUG instead of AUA, the latter being much more frequent than the former. This has been attributed to a compromise between translation initiation and elongation; i.e., AUG is not only a methionine codon, but also an initiation codon, and an anticodon matching AUG will increase the initiation rate. However, such an anticodon would impose selection against the use of AUA codons because AUA needs to be wobble-translated. According to this translation conflict hypothesis, AUA should be used relatively less frequently compared to UUA in the UUR codon family. A comprehensive analysis of mitochondrial genomes from a variety of vertebrate species revealed a general deficiency of AUA codons relative to UUA codons. In contrast, urochordate mitochondrial genomes with two tRNA(Met) genes with CAU and UAU anticodons exhibit increased AUA codon usage. Furthermore, six bivalve mitochondrial genomes with both of their tRNA-Met genes with a CAU anticodon have reduced AUA usage relative to three other bivalve mitochondrial genomes with one of their two tRNA-Met genes having a CAU anticodon and the other having a UAU anticodon. We conclude that the translation conflict hypothesis is empirically supported, and our results highlight the fine details of selection in shaping molecular evolution.

Huntingtin gene evolution in Chordata and its peculiar features in the ascidian Ciona genus

Background

To gain insight into the evolutionary features of the huntingtin (htt) gene in Chordata, we have sequenced and characterized the full-length htt mRNA in the ascidian Ciona intestinalis, a basal chordate emerging as new invertebrate model organism. Moreover, taking advantage of the availability of genomic and EST sequences, the htt gene structure of a number of chordate species, including the cogeneric ascidian Ciona savignyi, and the vertebrates Xenopus and Gallus was reconstructed.

Results

The C. intestinalis htt transcript exhibits some peculiar features, such as spliced leader trans-splicing in the 98 nt-long 5' untranslated region (UTR), an alternative splicing in the coding region, eight alternative polyadenylation sites, and no similarities of both 5' and 3'UTRs compared to homologs of the cogeneric C. savignyi. The predicted protein is 2946 amino acids long, shorter than its vertebrate homologs, and lacks the polyQ and the polyP stretches found in the the N-terminal regions of mammalian homologs. The exon-intron organization of the htt gene is almost identical among vertebrates, and significantly conserved between Ciona and vertebrates, allowing us to hypothesize an ancestral chordate gene consisting of at least 40 coding exons.

Conclusion

During chordate diversification, events of gain/loss, sliding, phase changes, and expansion of introns occurred in both vertebrate and ascidian lineages predominantly in the 5'-half of the htt gene, where there is also evidence of lineage-specific evolutionary dynamics in vertebrates. On the contrary, the 3'-half of the gene is highly conserved in all chordates at the level of both gene structure and protein sequence. Between the two Ciona species, a fast evolutionary rate and/or an early divergence time is suggested by the absence of significant similarity between UTRs, protein divergence comparable to that observed between mammals and fishes, and different distribution of repetitive elements.

The complete sequences and gene organisation of the mitochondrial genomes of the heterodont bivalves Acanthocardia tuberculata and Hiatella arctica--and the first record for a putative Atpase subunit 8 gene in marine bivalves

BACKGROUND

Mitochondrial (mt) gene arrangement is highly variable among molluscs and especially among bivalves. Of the 30 complete molluscan mt-genomes published to date, only one is of a heterodont bivalve, although this is the most diverse taxon in terms of species numbers. We determined the complete sequence of the mitochondrial genomes of Acanthocardia tuberculata and Hiatella arctica, (Mollusca, Bivalvia, Heterodonta) and describe their gene contents and genome organisations to assess the variability of these features among the Bivalvia and their value for phylogenetic inference.

RESULTS

The size of the mt-genome in Acanthocardia tuberculata is 16.104 basepairs (bp), and in Hiatella arctica 18.244 bp. The Acanthocardia mt-genome contains 12 of the typical protein coding genes, lacking the Atpase subunit 8 (atp8) gene, as all published marine bivalves. In contrast, a complete atp8 gene is present in Hiatella arctica. In addition, we found a putative truncated atp8 gene when re-annotating the mt-genome of Venerupis philippinarum. Both mt-genomes reported here encode all genes on the same strand and have an additional trnM. In Acanthocardia several large non-coding regions are present. One of these contains 3.5 nearly identical copies of a 167 bp motive. In Hiatella, the 3' end of the NADH dehydrogenase subunit (nad)6 gene is duplicated together with the adjacent non-coding region. The gene arrangement of Hiatella is markedly different from all other known molluscan mt-genomes, that of Acanthocardia shows few identities with the Venerupis philippinarum. Phylogenetic analyses on amino acid and nucleotide levels robustly support the Heterodonta and the sister group relationship of Acanthocardia and Venerupis. Monophyletic Bivalvia are resolved only by a Bayesian inference of the nucleotide data set. In all other analyses the two unionid species, being to only ones with genes located on both strands, do not group with the remaining bivalves.

CONCLUSION

The two mt-genomes reported here add to and underline the high variability of gene order and presence of duplications in bivalve and molluscan taxa. Some genomic traits like the loss of the atp8 gene or the encoding of all genes on the same strand are homoplastic among the Bivalvia. These characters, gene order, and the nucleotide sequence data show considerable potential of resolving phylogenetic patterns at lower taxonomic levels.

Darwin's dilemma: the realities of the Cambrian 'explosion'

The Cambrian 'explosion' is widely regarded as one of the fulcrum points in the history of life, yet its origins and causes remain deeply controversial. New data from the fossil record, especially of Burgess Shale-type Lagerstätten, indicate, however, that the assembly of bodyplans is not only largely a Cambrian phenomenon, but can already be documented in fair detail. This speaks against a much more ancient origin of the metazoans, and current work is doing much to reconcile the apparent discrepancies between the fossil record, including the Ediacaran assemblages of latest Neoproterozoic age and molecular 'clocks'. Hypotheses to explain the Cambrian 'explosion' continue to be generated, but the recurrent confusion of cause and effect suggests that the wrong sort of question is being asked. Here I propose that despite its step-like function this evolutionary event is the inevitable consequence of Earth and biospheric change.

The phylogeny of Mediterranean tortoises and their close relatives based on complete mitochondrial genome sequences from museum specimens

As part of an ongoing project to generate a mitochondrial database for terrestrial tortoises based on museum specimens, the complete mitochondrial genome sequences of 10 species and a approximately 14kb sequence from an eleventh species are reported. The sampling of the present study emphasizes Mediterranean tortoises (genus Testudo and their close relatives). Our new sequences are aligned, along with those of two testudinoid turtles from GenBank, Chrysemys picta and Mauremys reevesii, yielding an alignment of 14,858 positions, of which 3238 are parsimony informative. We develop a phylogenetic taxonomy for Testudo and related species based on well-supported, diagnosable clades. Several well-supported nodes are recovered, including the monophyly of a restricted Testudo, T. kleinmanni+T. marginata (the Chersus clade), and the placement of the enigmatic African pancake tortoise (Malacochersus tornieri) within the predominantly Palearctic greater Testudo group (Testudona tax. nov.). Despite the large amount of sequence reported, there is low statistical support for some nodes within Testudona and so we do not propose names for those groups. A preliminary and conservative estimation of divergence times implies a late Miocene diversification for the testudonan clade (6-10 million years ago), matching their first appearance in the fossil record. The multi-continental distribution of testudonan turtles can be explained by the establishment of permanent connections between Europe, Africa, and Asia at this time. The arrival of testudonan turtles to Africa occurred after one or more initial tortoise invasions gave rise to the diverse (>25 species) 'Geochelone complex.' Two unusual genomic features are reported for the mtDNA of one tortoise, M. tornieri: (1) nad4 has a shift of reading frame that we suggest is resolved by translational frameshifting of the mRNA on the ribosome during protein synthesis and (2) there are two copies of the control region and trnF, with the latter having experienced multiple-nucleotide substitutions in a pattern suggesting that each is being maintained by selection.

Complete mitochondrial DNA sequence of the eastern oyster Crassostrea virginica

The complete mitochondrial genome of the eastern oyster Crassostrea virginica (GenBank accession number AY905542) is 17,243 bp in length and contains 2 ribosomal genes, 12 protein-coding genes, and 23 transfer RNAs. The arrangement of protein-coding genes is identical to that of the congeneric Pacific oyster C. gigas, but tRNA genes show several duplications and extensive rearrangements between the species. Unique features in C. virginica include an additional trnM gene, the absence of an ATPase subunit 8 (atp8) gene, and an inferred translational frameshift within the cytochrome b (cob) gene. In both species the large subunit ribosomal RNA gene is encoded by 2 separate regions of the mitochondrial genome, the first reported case of a split ribosomal RNA gene in a metazoan. Translation of protein-coding genes in both species is initiated with methionine, with the exception of cob, which uses leucine. In C. virginica translation of all protein-coding genes (except possibly cob) terminates with TAA, with polyadenylation completing the primary transcript in cytochrome oxidase subunit III (cox3) and NADH dehydrogenase subunit 4L (nad4L), whereas C. gigas employs stop codons TAA and TAG equally. Interspecific divergence of mitochondrially encoded proteins is considerable, with amino acid identities ranging from 47% to 92%. A single major noncoding region representing the putative control region is found in both species.

Genetic characterization of the mitochondrial DNA from Lepeophtheirus salmonis (Crustacea; Copepoda). A new gene organization revealed

The mitochondrial DNA (mtDNA) from the salmon louse, Lepeophtheirus salmonis, is 15445 bp. It includes the genes coding for cytochrome B (Cyt B), ATPase subunit 6 and 8 (A6 and A8), NADH dehydrogenase subunits 1-6 and 4L (ND1, ND2, ND3, ND4, ND4L, ND5 and ND6), cytochrome c oxidase subunits I-III (COI, COII and COIII), two rRNA genes (12S rRNA and 16S rRNA) and 22 tRNAs. Two copies of tRNA-Lys are present in the mtDNA of L. salmonis, while tRNA-Cys was not identified. Both DNA strands contain coding regions in the salmon louse, in contrast to the other copepod characterized Tigriopus japonicus, but only a few genes overlap. In vertebrates, ND4 and ND4L are transcribed as one bicistronic mRNA, and are therefore localized together. The same organization is also found in crustaceans, with the exceptions of T. japonicus, Neocalanus cristatus and L. salmonis that deviate from this pattern. Another exception of the L. salmonis mtDNA is that A6 and A8 do not overlap, but are separated by several genes. The protein-coding genes have a bias towards AT-rich codons. The mitochondrial gene order in L. salmonis differs significantly from the copepods T. japonicus, Eucalanus bungii, N. cristatus and the other 13 crustaceans previously characterized. Furthermore, the mitochondrial rRNA genes are encoded on opposite strands in L. salmonis. This has not been found in any other arthropods, but has been reported in two starfish species. In a phylogenetic analysis, using an alignment of mitochondrial protein sequences, L. salmonis groups together with T. japonicus, being distant relatives to the other crustaceans.

Forces maintaining organellar genomes: is any as strong as genetic code disparity or hydrophobicity?

It remains controversial why mitochondria and chloroplasts retain the genes encoding a small subset of their constituent proteins, despite the transfer of so many other genes to the nucleus. Two candidate obstacles to gene transfer, suggested long ago, are that the genetic code of some mitochondrial genomes differs from the standard nuclear code, such that a transferred gene would encode an incorrect amino acid sequence, and that the proteins most frequently encoded in mitochondria are generally very hydrophobic, which may impede their import after synthesis in the cytosol. More recently it has been suggested that both these interpretations suffer from serious "false positives" and "false negatives": genes that they predict should be readily transferred but which have never (or seldom) been, and genes whose transfer has occurred often or early, even though this is predicted to be very difficult. Here I consider the full known range of ostensibly problematic such genes, with particular reference to the sequences of events that could have led to their present location. I show that this detailed analysis of these cases reveals that they are in fact wholly consistent with the hypothesis that code disparity and hydrophobicity are much more powerful barriers to functional gene transfer than any other. The popularity of the contrary view has led to the search for other barriers that might retain genes in organelles even more powerfully than code disparity or hydrophobicity; one proposal, concerning the role of proteins in redox processes, has received widespread support. I conclude that this abandonment of the original explanations for the retention of organellar genomes has been premature. Several other, relatively minor, obstacles to gene transfer certainly exist, contributing to the retention of relatively many organellar genes in most lineages compared to animal mtDNA, but there is no evidence for obstacles as severe as code disparity or hydrophobicity. One corollary of this conclusion is that there is currently no reason to suppose that engineering nuclear versions of the remaining mammalian mitochondrial genes, a feat that may have widespread biomedical relevance, should require anything other than sequence alterations obviating code disparity and causing modest reductions in hydrophobicity without loss of enzymatic function.

Complete nucleotide sequence of the mitochondrial genome of Doliolum nationalis with implications for evolution of urochordates

The evolutionary history of the diverse lifestyles adopted by urochordates has attracted intense interest because it may effect the evolutionary history of vertebrates. Here, we report the complete mitochondrial (mt) DNA sequence of the pelagic thaliacean doliolid Doliolum nationalis. The doliolid mt genome shares the unusual tRNAs of trnM(uau) and trnG(ucu) with other ascidians, such as Halocynthia and Ciona. On the other hand, the gene order of the doliolid mt genome is significantly different from that of any ascidian species or vertebrate reported to date. Phylogenetic analyses of the amino acid sequences of 12 protein-coding genes strongly support the sister-grouping of doliolids and the Phlebobranch ascidian Ciona, with the Stolidobranch ascidian alocynthia as the outgroup, thereby providing strong support for the paraphyly of ascidians, as has been suggested by 18S rDNA studies. Given the paraphyletic nature of ascidians, it seems likely that the common ancestor of ascidians and thaliaceans was sessile, as are the present-day ascidians, and that the thaliaceans subsequently evolved a pelagic lifestyle.