Changing identities: tRNA duplication and remolding within animal mitochondrial genomes

Incorporation of microgastropoda species provides novel insights into phylogeny of Trochoidea (Gastropoda: Vetigastropoda)

Trochoidea is the richest and most diverse group within Vetigastropoda, serving as one of the main focuses on studies of marine ecology and systematics. Both morphological and molecular studies have sought to resolve the phylogenetic framework of Trochoidea; however, the phylogenetic relationships among some lineages remain controversial. In order to explore the phylogenetic relationships within Trochoidea, we sequenced the mitochondrial genomes of 9 trochoids and analyzed them with data from 38 previously published mitochondrial genomes and 27 transcriptomic data representing 11 families within this group. The mitochondrial genomes of all Trochidae exhibited a consistent gene arrangement and showed conserved genome size and nucleotide composition; however, Colloniidae and Phasianellidae showed higher levels of gene order rearrangement. Furthermore, three-nucleotide insertions were observed in the cox1 gene of Colloniidae and nad4L gene of Phasianellidae, while three-nucleotide deletions were detected in nad4 gene of Trochidae. Both maximum likelihood and Bayesian inference analyses supported the monophyly of all families within the Trochoidea, except Tegulidae, and placed Liotiidae as sister to the rest of Trochoidea with poor to moderate support. Areneidae was recovered as the sister group to a clade including Phasianellidae and Colloniidae. These findings challenge the traditional classification of this family based on both morphological and molecular data. Our study provides new insights into the phylogeny of Trochoidea, especially with the incorporation of micromolluscs taxa, and highlight the significance of incorporating microgastropoda taxa into molecular phylogenetic reconstructions of gastropod subgroups.

Comparative Analysis of Metopograpsus quadridentatus (Crustacea: Decapoda: Grapsidae) Mitochondrial Genome Reveals Gene Rearrangement and Phylogeny

The taxonomy of the genus Metopograpsus is still a matter of debate because its species have limited morphological differences. Mitochondrial genomes, which are highly informative and easily accessible genetic markers, have been widely used to study molecular evolution and clarify relationships among species. This study aims to investigate the mitochondrial genome of Metopograpsus quadridentatus, a species with unique ecological and evolutionary significance. By analyzing the mitochondrial genome, we seek to address taxonomic uncertainties and provide insights into the evolutionary history of this species. In this study, we sequenced and analyzed the mitochondrial genome of M. quadridentatus to investigate its gene rearrangement patterns and its place within Brachyura. We compared the gene order of M. quadridentatus with that of 40 other Brachyuran species and created phylogenetic trees based on the nucleotide and amino acid sequences of 13 protein-coding genes (PCGs). We found that the mitochondrial gene arrangement of M. quadridentatus is mostly unchanged, similar to the original crustacean pattern, except for the movement of the trnH gene. Notably, the gene orders of several families, such as Eriphiidae, Grapsidae, Camptandriidae, Dotillidae, Plagusiidae, Ocypodidae, and Gecarcinidae, are the same. Phylogenetic analyses consistently placed M. quadridentatus within the genus Metopograpsus and the family Grapsidae, confirming its current taxonomic classification. These results offer important insights into evolutionary relationships and gene order conservation within Brachyura. Our study highlights the significance of mitochondrial genomes in resolving taxonomic uncertainties within the genus Metopograpsus.

Mitochondrial Genome of Grapsus albolineatus and Insights into the Phylogeny of Brachyura

Brachyura is among the most diverse groups of crustaceans, with over 7000 described species. Crab mitogenomes are important for understanding molecular evolution and phylogenetic relationships. Grapsus albolineatus exhibits specific rearrangements compared with the Pancrustacean ground pattern and other Brachyura species. The gene arrangement of G. albolineatus is similar to that of ancestral crustaceans, barring that of the translocated trnH gene. In phylogenetic analyses, the Bayesian inference estimation was observed to be superior to the maximum likelihood estimation when the nodal support values were compared. Considering the results of the gene rearrangement pattern and phylogenetic analysis, we speculate that G. albolineatus belongs to Grapsidae. Our comparative study indicated that mitogenomes are a useful phylogenetic tool at the subfamily level within Brachyura. The findings indicate that mitogenomes could be a useful tool for systematics in other Brachyuran species.

Aberrant Mitochondrial tRNA Genes Appear Frequently in Animal Evolution

Mitochondrial tRNAs have acquired a diverse portfolio of aberrant structures throughout metazoan evolution. With the availability of more than 12,500 mitogenome sequences, it is essential to compile a comprehensive overview of the pattern changes with regard to mitochondrial tRNA repertoire and structural variations. This, of course, requires reanalysis of the sequence data of more than 250,000 mitochondrial tRNAs with a uniform workflow. Here, we report our results on the complete reannotation of all mitogenomes available in the RefSeq database by September 2022 using mitos2. Based on the individual cases of mitochondrial tRNA variants reported throughout the literature, our data pinpoint the respective hotspots of change, i.e. Acanthocephala (Lophotrochozoa), Nematoda, Acariformes, and Araneae (Arthropoda). Less dramatic deviations of mitochondrial tRNAs from the norm are observed throughout many other clades. Loss of arms in animal mitochondrial tRNA clearly is a phenomenon that occurred independently many times, not limited to a small number of specific clades. The summary data here provide a starting point for systematic investigations into the detailed evolutionary processes of structural reduction and loss of mitochondrial tRNAs as well as a resource for further improvements of annotation workflows for mitochondrial tRNA annotation.

Complete mitochondrial genome and phylogenetic analysis of Mancinella alouina

BACKGROUND

The Muricidae family in the Class Gastropoda comprises numerous species with a vast range of morphological features and a worldwide presence. The phylogeny of the Muricidae has been analyzed in previous studies; however, the evolutionary relationships among the main branches of the Muricidae remain unknown.

METHODS AND RESULTS

In the present study, the mitochondrial genome of Mancinella alouina was sequenced. The mitochondrial genome was found to be 16,671 bp in length and made up of 37 genes (13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes). The genome has an A-T-rich region (66.5% A + T content) and all of the PCGs use the ATN start codon and the TAG or TAA stop codons. The mitochondrial gene arrangement of Mancinella alouina is similar to that of other Muricidae, except for Ocinebrellus inornatus and Ceratostoma burnetti. On the basis of a flexible molecular clock model, time-calibrated phylogenetic results indicate that the genus Mancinella diverged roughly 18.09 Mya, and that the family Muricidae emerged in the Late Cretaceous.

CONCLUSIONS

This study reveals the structural and sequence information features of the mitochondrial genome of Mancinella alouina. This study provides evidence for the relationships within the family Muricidae at the molecular level, and infer the divergence time. The results of phylogenetic analyses strongly support the current classification.

Characterization of the Complete Mitochondrial Genome of Agelas nakamurai from the South China Sea

The Agelas genus sponges are widely distributed and provide shelter for organisms that inhabit reefs. However, there is a lack of research on the genetic diversity of the Agelas sponges. Additionally, only one Agelas mitochondrial genome has been documented, leaving the characteristics of the Agelas genus's mitogenome in need of further clarification. To address this research gap, we utilized Illumina HiSeq4000 sequencing and de novo assembly to ascertain the complete mitochondrial genome of Agelas sp. specimens, sourced from the South China Sea. Our analysis of the cox1 barcoding similarity and phylogenetic relationship reveals that taxonomically, the Agelas sp. corresponds to Agelas nakamurai. The mitogenome of Agelas nakamurai is 20,885 bp in length, encoding 14 protein-coding genes, 24 transfer RNA genes, and 2 ribosomal RNA genes. Through a comparison of the mitochondrial genes, we discovered that both Agelas nakamurai and Agelas schmidti have an identical gene arrangement. Furthermore, we observed a deletion in the trnD gene and duplication and remodeling of the trnL gene in the Agelas nakamurai's mitogenome. Our evolutionary analysis also identified lineage-specific positive selection sites in the nad3 and nad5 genes of the Agelas sponges' mitogenome. These findings shed light on the gene rearrangement events and positive selection sites in the mitogenome of Agelas nakamurai, providing valuable molecular insights into the evolutionary processes of this genus.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

The compact mitogenome of Ceratocystiopsis pallidobrunnea

Ceratocystiopsis is a fungal genus that has been assigned to the Ophiostomatales, fungi known for their association with various bark beetles and other arthropods. The mitochondrial genome of Ceratocystiopsis pallidobrunnea has been characterized and compared with other members of the genus Ceratocystiopsis and Ophiostomatales. At 29 022 bp, the mitogenome of C. pallidobrunnea is the smallest reported so far for this genus. Gene arrangement was observed to be conserved for this group of fungi and mitogenome variation appears to be mostly due to the absence and presence of introns. The long-term goal is to apply mitogenomes to resolve taxonomic issues within the Ophiostomatales and within the various genera that comprise the Ophiostomataceae.

A novel gene order and remolded tRNAs revealed in the mitogenome of Asian gecarcinucid freshwater crabs (Brachyura, Gecarcinucidae)

Here we report the first mitochondrial genomes (mitogenomes) of four species of gecarcinucid freshwater crabs (FWCs) in two genera, two from China (Somanniathelphusa hainanensis and S. yangshanensis), one from Laos (Esanthelphusa dugasti), and one from Myanmar (Esanthelphusa keyini). A novel gecarcinucid mitochondrial gene order (GMGO2) that was only found in E. dugasti that contains a total of 42 genes, including one pseudogene, two remolded tRNAs and two duplicated tRNAs. The GMGO2 of E. dugasti was compared with the brachyuran ground-pattern mitochondrial gene order (BMGO), revealing the rearrangements of the positions of 10 tRNAs, two PCGs, and one mNCR. The three other gecarcinucids in this study were all found to possess a previously reported gecarcinucid mitochondrial gene order (GMGO1). The phylogenetic tree reconstructed using the secondary structures of 22 tRNAs of the mitogenomes of 41 species of FWCs provides insights into the evolution of the mitogenome of E. dugasti (GMGO2) which includes remolded and duplicated tRNAs.

Molluscan mitochondrial genomes break the rules

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

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Hidden cases of tRNA gene duplication and remolding in mitochondrial genomes of amphipods

The evolution of tRNA genes in mitochondrial (mt) genomes is a complex process that includes duplications, degenerations, and transpositions, as well as a specific process of identity change through mutations in the anticodon (tRNA gene remolding or tRNA gene recruitment). Using amphipod-specific tRNA models for annotation, we show that tRNA duplications are more common in the mt genomes of amphipods than what was revealed by previous annotations. Seventeen cases of tRNA gene duplications were detected in the mt genomes of amphipods, and ten of them were tRNA genes that underwent remolding. The additional tRNA gene findings were verified using phylogenetic analysis and genetic distance analysis. The majority of remolded tRNA genes (seven out of ten cases) were found in the mt genomes of endemic amphipod species from Lake Baikal. All additional mt tRNA genes arose independently in the Baikalian amphipods, indicating the unusual plasticity of tRNA gene evolution in these species assemblages. The possible reasons for the unusual abundance of additional tRNA genes in the mt genomes of Baikalian amphipods are discussed. The amphipod-specific tRNA models developed for MiTFi refine existing predictions of tRNA genes in amphipods and reveal additional cases of duplicated tRNA genes overlooked by using less specific Metazoa-wide models. The application of these models for mt tRNA gene prediction will be useful for the correct annotation of mt genomes of amphipods and probably other crustaceans.

The first divergence time estimation of the subfamily Stenogastrinae (Hymenoptera: Vespidae) based on mitochondrial phylogenomics

In this study, the mitochondrial genomes of three Stenogastrinae species, Eustenogaster scitula, Liostenogaster nitidipennis and Parishnogaster mellyi were sequenced and annotated, and a total of 16 vespid mtgenomes are comparatively analyzed. Our results indicate that codon usage bias is mainly dominated by mutational pressure, and affected only slightly by natural selection. Selective pressure analysis of protein-coding genes (PCGs) shows that the highest evolutionary rate is present in NADH complex I, and the lowest in cox1. Compared with the reported mtgenomes of other Vespidae, in Stenogastrinae, trnH is shifted to a new position. Phylogenetic analyses are performed using Bayesian method and Maximum Parsimony. Phylogenetic analysis further confirms that the Stenogastrinae is the sister group of all remaining Vespidae. Divergence time of Stenogastrinae from other Vespidae is estimated at ~ 166 Mya. Our results also support that eusociality evolved twice in the family Vespidae.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.

Complete mitochondrial genome of the first deep-sea spongicolid shrimp Spongiocaris panglao (Decapoda: Stenopodidea): Novel gene arrangement and the phylogenetic position and origin of Stenopodidea

Stenopodidea Claus, 1872 (Crustacea: Decapoda) is one of the major groups of decapods crustaceans. Hitherto, only one complete mitochondrial genome (mitogenome) from the family Stenopodidae is available for the infraorder Stenopodidea. Here, we determined the complete mitogenome of Spongiocaris panglao de Grave and Saito, 2016 using Illumina sequencing, representing the first species from the family Spongicolidae. The 15,909 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. Although the overall genome organization is typical for metazoans, the mitogenome of S. panglao shows some derived characters. A + T content of 77.42% in S. pamglao mitogenome is second-highest among the dacapods described to date. The trnR gene exhibit modified secondary structure with the TψC loop completely missing, which might be a putative autapomorphy of S. pamglao mitogenome. Compared with the shallow-water stenopodidean species S. hispidus, the control region of S. pamglao exhibits three characteristics: larger size, higher A + T content, and more tandem repeat sequences. The gene order exhibited difference from the ancestral mitogenome pattern of the Pancrustacea, with 5 tRNA genes rearrangement. The result from BI was agreed with most morphological characters and molecular evidences, revealing that Stenopodidea and Reptantia had the closest relationship, as the sister group of Caridea. Still, the alternative hypothesis supported from ML topology cannot be completely rejected based on the current data. Estimated times revealed that the two stenopodideans families Stenopodidae and Spongicolidae diverged from each other around 122 Mya. The divergence time of spongicolid shrimp is in good agreement with the origin of their hexactinellid hosts (78-144 Mya).

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

Mitochondrial Genomes of Two Bombycoidea Insects and Implications for Their Phylogeny

The mitochondrial genome (mt genome) provides important information for understanding molecular evolution and phylogenetics. As such, the two complete mt genomes of Ampelophaga rubiginosa and Rondotia menciana were sequenced and annotated. The two circular genomes of A. rubiginosa and R. menciana are 15,282 and 15,636 bp long, respectively, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The nucleotide composition of the A. rubiginosa mt genome is A + T rich (81.5%) but is lower than that of R. menciana (82.2%). The AT skew is slightly positive and the GC skew is negative in these two mt genomes. Except for cox1, which started with CGA, all other 12PCGs started with ATN codons. The A + T-rich regions of A. rubiginosa and R. menciana were 399 bp and 604 bp long and consist of several features common to Bombycoidea insects. The order and orientation of A. rubiginosa and R. menciana mitogenomes with the order trnM-trnI-trnQ-nad2 is different from the ancestral insects in which trnM is located between trnQ and nad2 (trnI-trnQ-trnM-nad2). Phylogenetic analyses indicate that A. rubiginosa belongs in the Sphingidae family, and R. menciana belongs in the Bombycidae family.

Complete mitochondrial genome of Clistocoeloma sinensis (Brachyura: Grapsoidea): Gene rearrangements and higher-level phylogeny of the Brachyura

Deciphering the animal mitochondrial genome (mitogenome) is very important to understand their molecular evolution and phylogenetic relationships. In this study, the complete mitogenome of Clistocoeloma sinensis was determined. The mitogenome of C. sinensis was 15,706 bp long, and its A+T content was 75.7%. The A+T skew of the mitogenome of C. sinensis was slightly negative (−0.020). All the transfer RNA genes had the typical cloverleaf structure, except for the trnS1 gene, which lacked a dihydroxyuridine arm. The two ribosomal RNA genes had 80.2% A+T content. The A+T-rich region spanned 684 bp. The gene order within the complete mitogenome of C. sinensis was identical to the pancrustacean ground pattern except for the translocation of trnH. Additionally, the gene order of trnI-trnQ-trnM in the pancrustacean ground pattern becomes trnQ-trnI-trnM in C. sinensis. Our phylogenetic analysis showed that C. sinensis and Sesarmops sinensis cluster together with high nodal support values, indicating that C. sinensis and S. sinensis have a sister group relationship. The results support that C. sinensis belongs to Grapsoidea, Sesarmidae. Our findings also indicate that Varunidae and Sesarmidae species share close relationships. Thus, mitogenomes are likely to be valuable tools for systematics in other groups of Crustacea.

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

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

The complete mitochondrial genome of Sesarmops sinensis reveals gene rearrangements and phylogenetic relationships in Brachyura

Mitochondrial genome (mitogenome) is very important to understand molecular evolution and phylogenetics. Herein, in this study, the complete mitogenome of Sesarmops sinensis was reported. The mitogenome was 15,905 bp in size, and contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The AT skew and the GC skew are both negative in the mitogenomes of S. sinensis. The nucleotide composition of the S. sinensis mitogenome was also biased toward A + T nucleotides (75.7%). All tRNA genes displayed a typical mitochondrial tRNA cloverleaf structure, except for the trnS1 gene, which lacked a dihydroxyuridine arm. S. sinensis exhibits a novel rearrangement compared with the Pancrustacean ground pattern and other Brachyura species. Based on the 13 PCGs, the phylogenetic analysis showed that S. sinensis and Sesarma neglectum were clustered on one branch with high nodal support values, indicating that S. sinensis and S. neglectum have a sister group relationship. The group (S. sinensis + S. neglectum) was sister to (Parasesarmops tripectinis + Metopaulias depressus), suggesting that S. sinensis belongs to Grapsoidea, Sesarmidae. Phylogenetic trees based on amino acid sequences and nucleotide sequences of mitochondrial 13 PCGs using BI and ML respectively indicate that section Eubrachyura consists of four groups clearly. The resulting phylogeny supports the establishment of a separate subsection Potamoida. These four groups correspond to four subsections of Raninoida, Heterotremata, Potamoida, and Thoracotremata.

Mitochondrial Genomes Provide Insights into the Phylogeny of Lauxanioidea (Diptera: Cyclorrhapha)

The superfamily Lauxanioidea is a significant dipteran clade including over 2500 known species in three families: Lauxaniidae, Celyphidae and Chamaemyiidae. We sequenced the first five (three complete and two partial) lauxanioid mitochondrial (mt) genomes, and used them to reconstruct the phylogeny of this group. The lauxanioid mt genomes are typical of the Diptera, containing all 37 genes usually present in bilaterian animals. A total of three conserved intergenic sequences have been reported across the Cyclorrhapha. The inferred secondary structure of 22 tRNAs suggested five substitution patterns among the Cyclorrhapha. The control region in the Lauxanioidea has apparently evolved very fast, but four conserved structural elements were detected in all three complete mt genome sequences. Phylogenetic relationships based on the mt genome data were inferred by Maximum Likelihood and Bayesian methods. The traditional relationships between families within the Lauxanioidea, (Chamaemyiidae + (Lauxaniidae + Celyphidae)), were corroborated; however, the higher-level relationships between cyclorrhaphan superfamilies are mostly poorly supported.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

Updated mitochondrial phylogeny of Pteriomorph and Heterodont Bivalvia, including deep-sea chemosymbiotic Bathymodiolus mussels, vesicomyid clams and the thyasirid clam Conchocele cf. bisecta

The mitochondrial genomes of bivalves have often been used for comparative genomics and for resolving phylogenetic relationships. More than 100 bivalve complete mitochondrial genomes have been sequenced to date. However, few mitochondrial genomes have been reported for deep-sea chemosymbiotic bivalves, which belong to the subclasses Pteriomorphia and Heterodonta. In the present study, we sequenced the mitochondrial genomes of eight deep-sea chemosymbiotic bivalve species: three species of Bathymodiolus mussels (B. japonicus, B. platifrons, and B. septemdierum), four species of vesicomyid clams (Abyssogena mariana, A. phaseoliformis, Isorropodon fossajaponicum, and Phreagena okutanii, all of which were formerly classified in the genus Calyptogena), and one species of thyasirid clam (Conchocele cf. bisecta). With a few exceptions, these mitochondrial genomes contained genes that are typical of metazoans: 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. The major non-coding region with a high A+T content of each genome, which contained tandem repeats and hairpins, was hypothesized to function as a control region. The phylogenetic trees of Pteriomorphia and Heterodonta were reconstructed based on the concatenated sequences of 14 shared genes. Bathymodiolus formed a monophyletic clade with asymbiotic Mytilidae mussels, the vesicomyid clams formed a monophyly that was sister to the Veneridae, and C. cf. bisecta branched basally in the Heterodonta. It is known that the gene orders of mitochondrial genomes vary among bivalves. To examine whether gene order variation exhibits phylogenetic signals, tree topologies based on the minimum number of gene rearrangements were reconstructed for two clades (superfamily Tellinoidea, which includes the Psammobiidae, Semelidae, Solecurtidae, and Tellinidae; and the clade comprising the Myidae, Mactridae, Arcticidae, Vesicomyidae, and Veneridae) with high statistical support in sequence-based phylogenies. The resulting tree topologies were almost identical to those of the sequence-based trees. Our present findings suggest that the evolution of bivalves could be precisely traced back through the analysis of mitochondrial genomes, and that such an analysis could contribute to understanding bivalve evolution and diversity.

Orthologs, turn-over, and remolding of tRNAs in primates and fruit flies

Background

Transfer RNAs (tRNAs) are ubiquitous in all living organism. They implement the genetic code so that most genomes contain distinct tRNAs for almost all 61 codons. They behave similar to mobile elements and proliferate in genomes spawning both local and non-local copies. Most tRNA families are therefore typically present as multicopy genes. The members of the individual tRNA families evolve under concerted or rapid birth-death evolution, so that paralogous copies maintain almost identical sequences over long evolutionary time-scales. To a good approximation these are functionally equivalent. Individual tRNA copies thus are evolutionary unstable and easily turn into pseudogenes and disappear. This leads to a rapid turnover of tRNAs and often large differences in the tRNA complements of closely related species. Since tRNA paralogs are not distinguished by sequence, common methods cannot not be used to establish orthology between tRNA genes.

Results

In this contribution we introduce a general framework to distinguish orthologs and paralogs in gene families that are subject to concerted evolution. It is based on the use of uniquely aligned adjacent sequence elements as anchors to establish syntenic conservation of sequence intervals. In practice, anchors and intervals can be extracted from genome-wide multiple sequence alignments. Syntenic clusters of concertedly evolving genes of different families can then be subdivided by list alignments, leading to usually small clusters of candidate co-orthologs. On the basis of recent advances in phylogenetic combinatorics, these candidate clusters can be further processed by cograph editing to recover their duplication histories. We developed a workflow that can be conceptualized as stepwise refinement of a graph of homologous genes. We apply this analysis strategy with different types of synteny anchors to investigate the evolution of tRNAs in primates and fruit flies. We identified a large number of tRNA remolding events concentrated at the tips of the phylogeny. With one notable exception all phylogenetically old tRNA remoldings do not change the isoacceptor class.

Conclusions

Gene families evolving under concerted evolution are not amenable to classical phylogenetic analyses since paralogs maintain identical, species-specific sequences, precluding the estimation of correct gene trees from sequence differences. This leaves conservation of syntenic arrangements with respect to “anchor elements” that are not subject to concerted evolution as the only viable source of phylogenetic information. We have demonstrated here that a purely synteny-based analysis of tRNA gene histories is indeed feasible. Although the choice of synteny anchors influences the resolution in particular when tight gene clusters are present, and the quality of sequence alignments, genome assemblies, and genome rearrangements limits the scope of the analysis, largely coherent results can be obtained for tRNAs. In particular, we conclude that a large fraction of the tRNAs are recent copies. This proliferation is compensated by rapid pseudogenization as exemplified by many very recent alloacceptor remoldings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2927-4) contains supplementary material, which is available to authorized users.

Duplication and Remolding of tRNA Genes in the Mitochondrial Genome of Reduvius tenebrosus (Hemiptera: Reduviidae)

Most assassin bugs are predators that act as important natural enemies of insect pests. Mitochondrial (mt) genomes of these insects are double-strand circular DNAs that encode 37 genes. In the present study, we explore the duplication and rearrangement of tRNA genes in the mt genome of Reduvius tenebrosus, the first mt genome from the subfamily Reduviinae. The gene order rearranges from CR (control region)-trnI-trnQ-trnM-ND2 to CR-trnQ-trnI2-trnI1-trnM-ND2. We identified 23 tRNA genes, including 22 tRNAs commonly found in insects and an additional trnI (trnI2), which has high sequence similarity to trnM. We found several pseudo genes, such as pseudo-trnI, pseudo-CR, and pseudo-ND2, in the hotspot region of gene rearrangement (between the control region and ND2). These features provided evidence that this novel gene order could be explained by the tandem duplication/random loss (TDRL) model. The tRNA duplication/anticodon mutation mechanism further explains the presence of trnI2, which is remolded from a duplicated trnM in the TDRL process (through an anticodon mutation of CAT to GAT). Our study also raises new questions as to whether the two events proceed simultaneously and if the remolded tRNA gene is fully functional. Significantly, the duplicated tRNA gene in the mitochondrial genome has evolved independently at least two times within assassin bugs.

Rearrangement of mitochondrial tRNA genes in flat bugs (Hemiptera: Aradidae)

The typical insect mitochondrial (mt) genome organization, which contains a single chromosome with 37 genes, was found in the infraorder Pentatomomorpha (suborder Heteroptera). The arrangement of mt genes in these true bugs is usually the same as the ancestral mt gene arrangement of insects. Rearrangement of transfer RNA (tRNA) genes, however, has been found in two subfamilies of flat bugs (Mezirinae and Calisiinae, family Aradidae). In this study, we sequenced the complete mt genomes of four species from three other subfamilies (Aradinae, Carventinae and Aneurinae). We found tRNA gene rearrangement in all of these four species. All of the rearranged tRNA genes are located between the mitochondrial control region and cox1, indicating this region as a hotspot for gene rearrangement in flat bugs; the rearrangement is likely caused by events of tandem duplication and random deletion of genes. Furthermore, our phylogenetic and dating analyses indicated that the swap of positions between trnQ and trnI occurred ~162 million years ago (MYA) in the most recent common ancestor of the five subfamilies of flat bugs investigated to date, whereas the swap of positions between trnC and trnW occurred later in the lineage leading to Calisiinae, and the translocation of trnC and trnY occurred later than 134 MYA in the lineage leading to Aradinae.

Mitochondrial genomes of praying mantises (Dictyoptera, Mantodea): rearrangement, duplication, and reassignment of tRNA genes

Insect mitochondrial genomes (mitogenomes) contain a conserved set of 37 genes for an extensive diversity of lineages. Previously reported dictyopteran mitogenomes share this conserved mitochondrial gene arrangement, although surprisingly little is known about the mitogenome of Mantodea. We sequenced eight mantodean mitogenomes including the first representatives of two families: Hymenopodidae and Liturgusidae. Only two of these genomes retain the typical insect gene arrangement. In three Liturgusidae species, the trnM genes have translocated. Four species of mantis (Creobroter gemmata, Mantis religiosa, Statilia sp., and Theopompa sp.-HN) have multiple identical tandem duplication of trnR, and Statilia sp. additionally includes five extra duplicate trnW. These extra trnR and trnW in Statilia sp. are erratically arranged and form another novel gene order. Interestingly, the extra trnW is converted from trnR by the process of point mutation at anticodon, which is the first case of tRNA reassignment for an insect. Furthermore, no significant differences were observed amongst mantodean mitogenomes with variable copies of tRNA according to comparative analysis of codon usage. Combined with phylogenetic analysis, the characteristics of tRNA only possess limited phylogenetic information in this research. Nevertheless, these features of gene rearrangement, duplication, and reassignment provide valuable information toward understanding mitogenome evolution in insects.

Towards a comprehensive picture of alloacceptor tRNA remolding in metazoan mitochondrial genomes

Remolding of tRNAs is a well-documented process in mitochondrial genomes that changes the identity of a tRNA. It involves a duplication of a tRNA gene, a mutation that changes the anticodon and the loss of the ancestral tRNA gene. The net effect is a functional tRNA that is more closely related to tRNAs of a different alloacceptor family than to tRNAs with the same anticodon in related species. Beyond being of interest for understanding mitochondrial tRNA function and evolution, tRNA remolding events can lead to artifacts in the annotation of mitogenomes and thus in studies of mitogenomic evolution. Therefore, it is important to identify and catalog these events. Here we describe novel methods to detect tRNA remolding in large-scale data sets and apply them to survey tRNA remolding throughout animal evolution. We identify several novel remolding events in addition to the ones previously mentioned in the literature. A detailed analysis of these remoldings showed that many of them are derived from ancestral events.

Comparative mitogenomics of plant bugs (Hemiptera: Miridae): identifying the AGG codon reassignments between serine and lysine

Insect mitochondrial genomes are very important to understand the molecular evolution as well as for phylogenetic and phylogeographic studies of the insects. The Miridae are the largest family of Heteroptera encompassing more than 11,000 described species and of great economic importance. For better understanding the diversity and the evolution of plant bugs, we sequence five new mitochondrial genomes and present the first comparative analysis of nine mitochondrial genomes of mirids available to date. Our result showed that gene content, gene arrangement, base composition and sequences of mitochondrial transcription termination factor were conserved in plant bugs. Intra-genus species shared more conserved genomic characteristics, such as nucleotide and amino acid composition of protein-coding genes, secondary structure and anticodon mutations of tRNAs, and non-coding sequences. Control region possessed several distinct characteristics, including: variable size, abundant tandem repetitions, and intra-genus conservation; and was useful in evolutionary and population genetic studies. The AGG codon reassignments were investigated between serine and lysine in the genera Adelphocoris and other cimicomorphans. Our analysis revealed correlated evolution between reassignments of the AGG codon and specific point mutations at the antidocons of tRNALys and tRNASer(AGN). Phylogenetic analysis indicated that mitochondrial genome sequences were useful in resolving family level relationship of Cimicomorpha. Comparative evolutionary analysis of plant bug mitochondrial genomes allowed the identification of previously neglected coding genes or non-coding regions as potential molecular markers. The finding of the AGG codon reassignments between serine and lysine indicated the parallel evolution of the genetic code in Hemiptera mitochondrial genomes.

tRNA anticodon shifts in eukaryotic genomes

Embedded in the sequence of each transfer RNA are elements that promote specific interactions with its cognate aminoacyl tRNA-synthetase. Although many such "identity elements" are known, their detection is difficult since they rely on unique structural signatures and the combinatorial action of multiple elements spread throughout the tRNA molecule. Since the anticodon is often a major identity determinant itself, it is possible to switch between certain tRNA functional types by means of anticodon substitutions. This has been shown to have occurred during the evolution of some genomes; however, the scale and relevance of "anticodon shifts" to the evolution of the tRNA multigene family is unclear. Using a synteny-conservation-based method, we detected tRNA anticodon shifts in groups of closely related species: five primates, 12 Drosophila, six nematodes, 11 Saccharomycetes, and 61 Enterobacteriaceae. We found a total of 75 anticodon shifts: 31 involving switches of identity (alloacceptor shifts) and 44 between isoacceptors that code for the same amino acid (isoacceptor shifts). The relative numbers of shifts in each taxa suggest that tRNA gene redundancy is likely the driving factor, with greater constraint on changes of identity. Sites that frequently covary with alloacceptor shifts are located at the extreme ends of the molecule, in common with most known identity determinants. Isoacceptor shifts are associated with changes in the midsections of the tRNA sequence. However, the mutation patterns of anticodon shifts involving the same identities are often dissimilar, suggesting that alternate sets of mutation may achieve the same functional compensation.

tRNA genes rapidly change in evolution to meet novel translational demands

Changes in expression patterns may occur when organisms are presented with new environmental challenges, for example following migration or genetic changes. To elucidate the mechanisms by which the translational machinery adapts to such changes, we perturbed the tRNA pool of Saccharomyces cerevisiae by tRNA gene deletion. We then evolved the deletion strain and observed that the genetic adaptation was recurrently based on a strategic mutation that changed the anticodon of other tRNA genes to match that of the deleted one. Strikingly, a systematic search in hundreds of genomes revealed that anticodon mutations occur throughout the tree of life. We further show that the evolution of the tRNA pool also depends on the need to properly couple translation to protein folding. Together, our observations shed light on the evolution of the tRNA pool, demonstrating that mutation in the anticodons of tRNA genes is a common adaptive mechanism when meeting new translational demands.

DOI:http://dx.doi.org/10.7554/eLife.01339.001

Genes contain the blueprints for the proteins that are essential for countless biological functions and processes, and the path that leads from a particular gene to the corresponding protein is long and complex. The genetic information stored in the DNA must first be transcribed to produce a messenger RNA molecule, which then has to be translated to produce a string of amino acids that fold to form a protein. The translation step is performed by a molecular machine called the ribosome, with transfer RNA molecules bringing the amino acids that are needed to make the protein.

The information in messenger RNA is stored as a series of letters, with groups of three letters called codons representing the different amino acids. Since there are four letters—A, C, G and U—it is possible to form 64 different codons. And since there are only 20 amino acids, two or more different codons can specify the same amino acid (for example, AGU and AGC both specify serine), and two or more different transfer RNA molecules can take this amino acid to the ribosome. Moreover, some codons are found more often than others in the messenger RNA molecules, so the genes that encode the related transfer RNA molecules are more common than the genes for other transfer RNA molecules.

Environmental pressures mean that organisms must adapt to survive, with some genes and proteins increasing in importance, and others becoming less important. Clearly the relative numbers of the different transfer RNA molecules will also need to change to reflect these evolutionary changes, but the details of how this happens were not understood.

Now Yona et al. have explored this issue by studying yeast cells that lack a gene for one of the less common transfer RNA molecules (corresponding to the codon AGG, which specifies the amino acid arginine). At first this mutation resulted in slower growth of the yeast cells, but after being allowed to evolve over 200 generations, the rate of growth matched that of a normal strain with all transfer RNA genes. Yona et al. found that the gene for a more common transfer RNA molecule, corresponding to the codon AGA, which also specifies arginine, had mutated to AGG. As a result, the mutated yeast was eventually able to produce proteins as quickly as wild type yeast. Moreover, further experiments showed that the levels of some transfer RNAs are kept deliberately low in order to slow down the production of proteins so as to ensure that the proteins assume their correct structure.

But does the way these cells evolved in the lab resemble what happened in nature? To address this question Yona et al. examined a database of transfer RNA sequences from more than 500 species, and found evidence for the same codon-based switching mechanism in many species across the tree of life.

DOI:http://dx.doi.org/10.7554/eLife.01339.002

Ancestral genome organization: an alignment approach

We present a comparative genomics approach for inferring ancestral genome organization and evolutionary scenarios, based on present-day genomes represented as ordered gene sequences with duplicates. We develop our methodology for a model of evolution restricted to duplication and loss, and then show how to extend it to other content-modifying operations, and to inversions. From a combinatorial point of view, the main consequence of ignoring rearrangements is the possibility of formulating the problem as an alignment problem. On the other hand, duplications and losses are asymmetric operations that are applicable to one of the two aligned sequences. Consequently, an ancestral genome can directly be inferred from a duplication-loss scenario attached to a given alignment. Although alignments are a priori simpler to handle than rearrangements, we show that a direct approach based on dynamic programming leads, at best, to an efficient heuristic. We present an exact pseudo-boolean linear programming algorithm to search for the optimal alignment along with an optimal scenario of duplications and losses. Although exponential in the worst case, we show low running times on real datasets as well as synthetic data. We apply our algorithm (*) in a phylogenetic context to the evolution of stable RNA (tRNA and rRNA) gene content and organization in Bacillus genomes. Our results lead to various biological insights, such as rates of ribosomal RNA proliferation among lineages, their role in altering tRNA gene content, and evidence of tRNA class conversion.

Evolution of extensively fragmented mitochondrial genomes in the lice of humans

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. The human body louse, Pediculus humanus, however, has its mt genes on 20 minichromosomes. We sequenced the mt genomes of two other human lice: the head louse, P. capitis, and the pubic louse, Pthirus pubis. Comparison among the three human lice revealed the presence of fragmented mt genomes in their most recent common ancestor, which lived ∼7 Ma. The head louse has exactly the same set of mt minichromosomes as the body louse, indicating that the number of minichromosomes, and the gene content and gene arrangement in each minichromosome have remained unchanged since the body louse evolved from the head louse ∼107,000 years ago. The pubic louse has the same pattern of one protein-coding or rRNA gene per minichromosome (except one minichromosome with two protein-coding genes, atp6 and atp8) as the head louse and the body louse. This pattern is apparently ancestral to all human lice and has been stable for at least 7 Myr. Most tRNA genes of the pubic louse, however, are on different minichromosomes when compared with their counterparts in the head louse and the body louse. It is evident that rearrangement of four tRNA genes (for leucine, arginine and glycine) was due to gene-identity switch by point mutation at the third anticodon position or by homologous recombination, whereas rearrangement of other tRNA genes was by gene translocation between minichromosomes, likely caused by minichromosome split via gene degeneration and deletion.

DNA barcoding methods for invertebrates

Invertebrates comprise approximately 34 phyla, while vertebrates represent one subphylum and insects a (very large) class. Thus, the clades excepting vertebrates and insects encompass almost all of animal diversity. Consequently, the barcoding challenge in invertebrates is that of barcoding animals in general. While standard extraction, cleaning, PCR methods, and universal primers work for many taxa, taxon-specific challenges arise because of the shear genetic and biochemical diversity present across the kingdom, and because problems arising as a result of this diversity, and solutions to them, are still poorly characterized for many metazoan clades. The objective of this chapter is to emphasize general approaches, and give practical advice for overcoming the diverse challenges that may be encountered across animal taxa, but we stop short of providing an exhaustive inventory. Rather, we encourage researchers, especially those working on poorly studied taxa, to carefully consider methodological issues presented below, when standard approaches perform poorly.

Mitogenomics does not resolve deep molluscan relationships (yet?)

The origin of molluscs among lophotrochozoan metazoans is unresolved and interclass relationships are contradictory between morphology-based, multi-locus, and recent phylogenomic analyses. Within the "Deep Metazoan Phylogeny" framework, all available molluscan mitochondrial genomes were compiled, covering 6 of 8 classes. Genomes were reannotated, and 13 protein coding genes (PCGs) were analyzed in various taxon settings, under multiple masking and coding regimes. Maximum Likelihood based methods were used for phylogenetic reconstructions. In all cases, molluscs result mixed up with lophotrochozoan outgroups, and most molluscan classes with more than single representatives available are non-monophyletic. We discuss systematic errors such as long branch attraction to cause aberrant, basal positions of fast evolving ingroups such as scaphopods, patellogastropods and, in particular, the gastropod subgroup Heterobranchia. Mitochondrial sequences analyzed either as amino acids or nucleotides may perform well in some (Cephalopoda) but not in other palaeozoic molluscan groups; they are not suitable to reconstruct deep (Cambrian) molluscan evolution. Supposedly "rare" mitochondrial genome level features have long been promoted as phylogenetically informative. In our newly annotated data set, features such as genome size, transcription on one or both strands, and certain coupled pairs of PCGs show a homoplastic, but obviously non-random distribution. Apparently congruent (but not unambiguous) signal for non-trivial subclades, e.g. for a clade composed of pteriomorph and heterodont bivalves, needs confirmation from a more comprehensive bivalve sampling. We found that larger clusters not only of PCGs but also of rRNAs and even tRNAs can bear local phylogenetic signal; adding trnG-trnE to the end of the ancestral cluster trnM-trnC-trnY-trnW-trnQ might be synapomorphic for Mollusca. Mitochondrial gene arrangement and other genome level features explored and reviewed herein thus failed as golden bullets, but are promising as additional characters or evidence supporting deep molluscan clades revealed by other data sets. A representative and dense sampling of molluscan subgroups may contribute to resolve contentious interclass relationships in the future, and is vital for exploring the evolution of especially diverse mitochondrial genomes in molluscs.

Genetic aspects of mitochondrial genome evolution

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

Next generation sequencing and comparative analyses of Xenopus mitogenomes

Background

Mitochondrial genomes comprise a small but critical component of the total DNA in eukaryotic organisms. They encode several key proteins for the cell’s major energy producing apparatus, the mitochondrial respiratory chain. Additonally, their nucleotide and amino acid sequences are of great utility as markers for systematics, molecular ecology and forensics. Their characterization through nucleotide sequencing is a fundamental starting point in mitogenomics. Methods to amplify complete mitochondrial genomes rapidly and efficiently from microgram quantities of tissue of single individuals are, however, not always available. Here we validate two approaches, which combine long-PCR with Roche 454 pyrosequencing technology, to obtain two complete mitochondrial genomes from individual amphibian species.

Results

We obtained two new xenopus frogs (Xenopus borealis and X. victorianus) complete mitochondrial genome sequences by means of long-PCR followed by 454 of individual genomes (approach 1) or of multiple pooled genomes (approach 2), the mean depth of coverage per nucleotide was 9823 and 186, respectively. We also characterised and compared the new mitogenomes against their sister taxa; X. laevis and Silurana tropicalis, two of the most intensely studied amphibians. Our results demonstrate how our approaches can be used to obtain complete amphibian mitogenomes with depths of coverage that far surpass traditional primer-walking strategies, at either the same cost or less. Our results also demonstrate: that the size, gene content and order are the same among xenopus mitogenomes and that S. tropicalis form a separate clade to the other xenopus, among which X. laevis and X. victorianus were most closely related. Nucleotide and amino acid diversity was found to vary across the xenopus mitogenomes, with the greatest diversity observed in the Complex 1 gene nad4l and the least diversity observed in Complex 4 genes (cox1-3). All protein-coding genes were shown to be under strong negative (purifying selection), with genes under the strongest pressure (Complex 4) also being the most highly expressed, highlighting their potentially crucial functions in the mitochondrial respiratory chain.

Conclusions

Next generation sequencing of long-PCR amplicons using single taxon or multi-taxon approaches enabled two new species of Xenopus mtDNA to be fully characterized. We anticipate our complete mitochondrial genome amplification methods to be applicable to other amphibians, helpful for identifying the most appropriate markers for differentiating species, populations and resolving phylogenies, a pressing need since amphibians are undergoing drastic global decline. Our mtDNAs also provide templates for conserved primer design and the assembly of RNA and DNA reads following high throughput “omic” techniques such as RNA- and ChIP-seq. These could help us better understand how processes such mitochondrial replication and gene expression influence xenopus growth and development, as well as how they evolved and are regulated.

MITOS: improved de novo metazoan mitochondrial genome annotation

About 2000 completely sequenced mitochondrial genomes are available from the NCBI RefSeq data base together with manually curated annotations of their protein-coding genes, rRNAs, and tRNAs. This annotation information, which has accumulated over two decades, has been obtained with a diverse set of computational tools and annotation strategies. Despite all efforts of manual curation it is still plagued by misassignments of reading directions, erroneous gene names, and missing as well as false positive annotations in particular for the RNA genes. Taken together, this causes substantial problems for fully automatic pipelines that aim to use these data comprehensively for studies of animal phylogenetics and the molecular evolution of mitogenomes. The MITOS pipeline is designed to compute a consistent de novo annotation of the mitogenomic sequences. We show that the results of MITOS match RefSeq and MitoZoa in terms of annotation coverage and quality. At the same time we avoid biases, inconsistencies of nomenclature, and typos originating from manual curation strategies. The MITOS pipeline is accessible online at http://mitos.bioinf.uni-leipzig.de.

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.

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.

Adding resolution to ordinal level relationships of tapeworms (Platyhelminthes: Cestoda) with large fragments of mtDNA

The construction of a stable phylogeny for the Cestoda, indicating the interrelationships of recognised orders and other major lineages, has proceeded iteratively since the group first received attention from phylogenetic systematists. Molecular analyses using nuclear ribosomal RNA gene fragments from the small (ssrDNA) and large (lsrDNA) subunits have been used to test competing evolutionary scenarios based on morphological data but could not arbitrate between some key conflicting hypotheses. To the ribosomal data, we have added a contiguous fragment of mitochondrial (mt) genome data (mtDNA) of partial nad1-trnN-trnP-trnI-trnK-nad3-trnS-trnW-cox1-trnT-rrnL-trnC-partial rrnS, spanning 4034-4447 bp, where new data for this region were generated for 18 species. Bayesian analysis of mtDNA and rDNA as nucleotides, and where appropriate as amino acids, demonstrated that these two classes of genes provide complementary signal across the phylogeny. In all analyses, except when using mt amino acids only, the Gyrocotylidea is sister group to all other Cestoda (Nephroposticophora), and Amphilinidea forms the sister group to the Eucestoda. However, an earliest-diverging position of Amphilinidea is strongly supported in the mt amino acid analysis. Amphilinidea exhibit a unique tRNA arrangement (nad1-trnI-trnL2-trnP-trnK-trnV-trnA-trnN-nad3), whereas Gyrocotylidea shares that of the derived lineages, providing additional evidence of the uniqueness of amphilinid genes and genomes. The addition of mtDNA to the rDNA genes supported the Caryophyllidea as the sister group to (Spathebothriidea+remaining Eucestoda), a hypothesis consistently supported by morphology. This relationship suggests a history of step-wise evolutionary transitions from simple monozoic, unsegmented tapeworms to the more familiar polyzoic, externally segmented (strobilate) forms. All our data partitions recovered Haplobothriidea as the sister group to Diphyllobothriidae. The sister-group relationship between Diphyllidea and Trypanorhyncha, as previously established using rDNA, is not supported by the mt data, although it is supported by the combined mt and rDNA analysis. With regards to the more derived taxa, in all except the mt amino acid analysis, the following topology is supported: (Bothriocephalidea (Litobothriidea (Lecanicephalidea (Rhinebothriidea (Tetraphyllidea, (Acanthobothrium, Proteocephalidea), (Nippotaeniidea, Mesocestoididae, Tetrabothriidea, Cyclophyllidea)))))), where the Tetraphyllidea are paraphyletic. Evidence from the mt data provides strong (nucleotides) to moderate (amino acids) support for Tetraphyllidea forming a group to the inclusion of Proteocephalidea, with the latter consistently forming the sister group to Acanthobothrium. The interrelationships among Nippotaeniidea, Mesocestoididae, Tetrabothriidea and Cyclophyllidea remain ambiguous and require further systematic attention. Mitochondrial and nuclear rDNA data provide conflicting signal for certain parts of the cestode tree. In some cases mt data offer results in line with morphological evidence, such as the interrelationships of the early divergent lineages. Also, Tetraphyllidea, although remaining paraphyletic with the inclusion of the Proteocephalidea, does not include the most derived cestodes; a result which has consistently been obtained with rDNA.

Complete mitochondrial DNA sequence of Stenopus hispidus (Crustacea: Decapoda: Stenopodidea) and a novel tRNA gene cluster

As a phylogenetically valuable decapoda, a complete mitochondrial genome from Stenopodidea has not been reported to date. Here, we determined the complete mitochondrial DNA sequence of Stenopus hispidus (Olivier, 1811). The 15,528 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs) and two ribosomal RNA (rRNA) genes plus the putative control region (CR). This finding is similar to other metazoan animals but with the exception of 23 transfer RNA (tRNA) genes, which contain an additional tRNA-Gln compared with other crustaceans. With respect to the pancrustacean ground pattern mitochondria gene order, 5 tRNAs appear to be rearranged (tRNAs-Leu (CUN), Arg, Glu, Gln, and Met), one of which has also undergone inversion (tRNA-Leu (CUN)). Phylogenetic analyses reveal Stenopodidea and Reptantia form a clade sister to Caridea, which agrees with Abele and Felgenhauer's (1986) hypothesis. This topology contrasts with previous results based on morphological and some molecular data.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

The mitochondrial genome of Sinentomon erythranum (Arthropoda: Hexapoda: Protura): an example of highly divergent evolution

BACKGROUND

The phylogenetic position of the Protura, traditionally considered the most basal hexapod group, is disputed because it has many unique morphological characters compared with other hexapods. Although mitochondrial genome information has been used extensively in phylogenetic studies, such information is not available for the Protura. This has impeded phylogenetic studies on this taxon, as well as the evolution of the arthropod mitochondrial genome.

RESULTS

In this study, the mitochondrial genome of Sinentomon erythranum was sequenced, as the first proturan species to be reported. The genome contains a number of special features that differ from those of other hexapods and arthropods. As a very small arthropod mitochondrial genome, its 14,491 nucleotides encode 37 typical mitochondrial genes. Compared with other metazoan mtDNA, it has the most biased nucleotide composition with T = 52.4%, an extreme and reversed AT-skew of -0.351 and a GC-skew of 0.350. Two tandemly repeated regions occur in the A+T-rich region, and both could form stable stem-loop structures. Eighteen of the 22 tRNAs are greatly reduced in size with truncated secondary structures. The gene order is novel among available arthropod mitochondrial genomes. Rearrangements have involved in not only small tRNA genes, but also PCGs (protein-coding genes) and ribosome RNA genes. A large block of genes has experienced inversion and another nearby block has been reshuffled, which can be explained by the tandem duplication and random loss model. The most remarkable finding is that trnL2(UUR) is not located between cox1 and cox2 as observed in most hexapod and crustacean groups, but is between rrnL and nad1 as in the ancestral arthropod ground pattern. The "cox1-cox2" pattern was further confirmed in three more representative proturan species. The phylogenetic analyses based on the amino acid sequences of 13 mitochondrial PCGs suggest S. erythranum failed to group with other hexapod groups.

CONCLUSIONS

The mitochondrial genome of S. erythranum shows many different features from other hexapod and arthropod mitochondrial genomes. It underwent highly divergent evolution. The "cox1-cox2" pattern probably represents the ancestral state for all proturan mitogenomes, and suggests a long evolutionary history for the Protura.

The mitogenome of Paphia euglypta (Bivalvia: Veneridae) and comparative mitogenomic analyses of three venerids

Extraordinary variation has been found in mitochondrial (mt) genome inheritance, gene content and arrangement among bivalves. However, only few bivalve mt genomes have been comparatively analyzed to infer their evolutionary scenarios. In this study, the complete mt genome of the venerid Paphia euglypta (Bivalvia: Veneridae) was firstly studied and, secondly, it was comparatively analyzed with other venerids (e.g., Venerupis philippinarum and Meretrix petechialis) to better understand the mt genome evolution within a family. Though several common features such as the AT content, codon usage of protein-coding genes, and AT/GC skew are shared by the three venerids, a high level of variability is observed in genome size, gene content, gene order, arrangements and primary sequence of nucleotides or amino acids. Most of the gene rearrangement can be explained by the "tandem duplication and random loss" model. From the observed rearrangement patterns, we speculate that block interchange between adjacent genes may be common in the evolution of mt genomes in venerids. Furthermore, this study presents several new findings in mt genome annotation of V. philippinarum and M. petechialis, and hence we have reannotated the genome of these two species as: (1) the ORF of the formerly annotated cox2 gene in V. philippinarum is deduced by using a truncated "T" codon and a second cox2 gene is identified; (2) the trnS-AGN gene is identified and marked in the mt genome of both venerids. Thus, this study demonstrated a high variability of mt genomes in the Veneridae, and showed the importance of comparative mt genome analysis to interpret the evolution of the bivalve mt genome.

Gene recruitment--a common mechanism in the evolution of transfer RNA gene families

The evolution of alloacceptor transfer RNAs (tRNAs) has been traditionally thought to occur vertically and reflect the evolution of the genetic code. Yet there have been several indications that a tRNA gene could evolve horizontally, from a copy of an alloacceptor tRNA gene in the same genome. Earlier, we provided the first unambiguous evidence for the occurrence of such "tRNA gene recruitment" in nature--in the mitochondrial (mt) genome of the demosponge Axinella corrugata. Yet the extent and the pattern of this process in the evolution of tRNA gene families remained unclear. Here we analyzed tRNA genes from 21 mt genomes of demosponges as well as nuclear genomes of rhesus macaque, chimpanzee and human. We found four new cases of alloacceptor tRNA gene recruitment in mt genomes and eleven cases in the nuclear genomes. In most of these cases we observed a single nucleotide substitution at the middle position of the anticodon, which resulted in the change of not only the tRNA's amino-acid identity but also the class of the amino-acyl tRNA synthetases (aaRSs) involved in amino-acylation. We hypothesize that the switch to a different class of aaRSs may have prevented the conflict between anticodon and amino-acid identities of recruited tRNAs. Overall our results suggest that gene recruitment is a common phenomenon in tRNA multigene family evolution and should be taken into consideration when tRNA evolutionary history is reconstructed.