The Duplication/Random Loss Model for Gene Rearrangement Exemplified by Mitochondrial Genomes of Deuterostome Animals

Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

Selected Lark Mitochondrial Genomes Provide Insights into the Evolution of Second Control Region with Tandem Repeats in Alaudidae (Aves, Passeriformes)

The control region (CR) regulates the replication and transcription of the mitochondrial genome (mitogenome). Some avian mitogenomes possess two CRs, and the second control region (CR2) may enhance replication and transcription; however, the CR2 in lark mitogenome appears to be undergoing loss and is accompanied by tandem repeats. Here, we characterized six lark mitogenomes from Alaudala cheleensis, Eremophila alpestris, Alauda razae, and Calandrella cinerea and reconstructed the phylogeny of Passerida. Through further comparative analysis among larks, we traced the evolutionary process of CR2. The mitochondrial gene orders were conserved in all published lark mitogenomes, with Cytb-trnT-CR1-trnP-ND6-trnE-remnant CR2 with tandem repeat-trnF-rrnS. Phylogenetic analysis revealed Alaudidae and Panuridae are sister groups at the base of Sylvioidea, and sporadic losses of CR2 may occur in their common ancestor. CR sequence and phylogeny analysis indicated CR2 tandem repeats were generated within CR2, originating in the ancestor of all larks, rather than inherited from CR1. The secondary structure comparison of tandem repeat units within and between species suggested slipped-strand mispairing and DNA turnover as suitable models for explaining the origin and evolution of these repeats. This study reveals the evolutionary process of the CR2 containing tandem repeat in Alaudidae, providing reference for understanding the evolutionary characteristics and dynamics of tandem repeats.

A Mysterious Asian Firefly Genus, Oculogryphus Jeng, Engel & Yang (Coleoptera, Lampyridae): The First Complete Mitochondrial Genome and Its Phylogenetic Implications

The firefly genus Oculogryphus Jeng, Engel & Yang, 2007 is a rare-species group endemic to Asia. Since its establishment, its position has been controversial but never rigorously tested. To address this perplexing issue, we are the first to present the complete mitochondrial sequence of Oculogryphus, using the material of O. chenghoiyanae Yiu & Jeng, 2018 determined through a comprehensive morphological identification. Our analyses demonstrate that its mitogenome exhibits similar characteristics to that of Stenocladius, including a rearranged gene order between trnC and trnW, and a long intergenic spacer (702 bp) between the two rearranged genes, within which six remnants (29 bp) of trnW were identified. Further, we incorporated this sequence into phylogenetic analyses of Lampyridae based on different molecular markers and datasets using ML and BI analyses. The results consistently place Oculogryphus within the same clade as Stenocladius in all topologies, and the gene rearrangement is a synapomorphy for this clade. It suggests that Oculogryphus should be classified together with Stenocladius in the subfamily Ototretinae at the moment. This study provides molecular evidence confirming the close relationship between Oculogryphus and Stenocladius and discovers a new phylogenetic marker helpful in clarifying the monophyly of Ototretinae, which also sheds a new light on firefly evolution.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

Multiple independent origins of duplicated mitochondrial control regions indicate an apomorphy in the Thysanoptera (Insecta)

The mitochondrial genome (mitogenome) of thrips is characterized by the presence of control region (CR) duplication. However, the evolution pattern of duplicated CRs in thrips is still unclear. In this study, the multiple independent origins of duplicated CR indicated that the CR duplication was not an ancestral state for Thysanoptera. The macroevolutionary pattern suggested that the earliest CR duplication event occurred in the middle Cretaceous (94.85 Ma) coincided with rearrangement events forming the ancestors of Aeolothripidae, but much later than that forming the ancestors of the suborder Terebrantia. The mitogenome with duplicated CRs showed a higher rate of gene rearrangement. The sequence similarity of the CR copies and divergence time were negatively correlated, indicating age-related deterioration of mitochondrial function. No significant differences were found in the mitochondrial DNA, the P123 and P4FD between the single and multiple-CR charactered mitogenomes, which suggested that the duplicated CRs may not affect the replication process in thrip mitogenome. The mitogenomes with duplicated CRs (mean: 0.0088 subs/s/my) show a significantly increased evolutionary rate than that with a single one (mean: 0.0058 subs/s/my). However, it seems that this higher evolutionary rate did not have adaptive mechanisms in Terebrantia. We speculated that the duplicated CRs may cause a more intense production of energy by mitochondria, and an accelerated mutation and substitution rate is expected in such mitogenomes. Our study provided new insights into the presence of CR duplications and their evolution in the mitogenomes of thrips.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Detecting gene breakpoints in noisy genome sequences using position-annotated colored de-Bruijn graphs

BACKGROUND

Identifying the locations of gene breakpoints between species of different taxonomic groups can provide useful insights into the underlying evolutionary processes. Given the exact locations of their genes, the breakpoints can be computed without much effort. However, often, existing gene annotations are erroneous, or only nucleotide sequences are available. Especially in mitochondrial genomes, high variations in gene orders are usually accompanied by a high degree of sequence inconsistencies. This makes accurately locating breakpoints in mitogenomic nucleotide sequences a challenging task.

RESULTS

This contribution presents a novel method for detecting gene breakpoints in the nucleotide sequences of complete mitochondrial genomes, taking into account possible high substitution rates. The method is implemented in the software package DeBBI. DeBBI allows to analyze transposition- and inversion-based breakpoints independently and uses a parallel program design, allowing to make use of modern multi-processor systems. Extensive tests on synthetic data sets, covering a broad range of sequence dissimilarities and different numbers of introduced breakpoints, demonstrate DeBBI 's ability to produce accurate results. Case studies using species of various taxonomic groups further show DeBBI 's applicability to real-life data. While (some) multiple sequence alignment tools can also be used for the task at hand, we demonstrate that especially gene breaks between short, poorly conserved tRNA genes can be detected more frequently with the proposed approach.

CONCLUSION

The proposed method constructs a position-annotated de-Bruijn graph of the input sequences. Using a heuristic algorithm, this graph is searched for particular structures, called bulges, which may be associated with the breakpoint locations. Despite the large size of these structures, the algorithm only requires a small number of graph traversal steps.

The mitochondrial genome of Binodoxys acalephae (Hymenoptera: Braconidae) with unique gene rearrangement and phylogenetic implications

BACKGROUND

Species in the subfamily Aphidiinae from the Braconidae of Hymenoptera are endoparasitic wasps that exclusively utilize aphids as hosts. Some Aphidiinae species are widely used as biological agents. However, there were only one species with determined complete mitochondrial genome from this subfamily.

METHODS AND RESULTS

In this study, we sequenced and annotated the mitochondrial genome (mitogenome) of Binodoxys acalephae, which was 15,116 bp in size and contained 37 genes. The start codon of 13 protein-coding genes was ATN, and the complete stop codon TAA and TAG was widely assigned to 11 protein-coding genes. The lrRNA contains 43 stem-loop structures, and srRNA contains 25 stem-loop structures. Translocation and inversion of tRNA genes was found to be dominant in B. acalephae. In contrast to Aphidius gifuensis from the same subfamily Aphidiinae, inverted tRNALeu1 was translocated to the gene cluster between tRNALeu2 and COX2, and the control region between tRNAIle and tRNAMet was deleted in the mitogenome of B. acalephae. Within Braconidae, gene clusters tRNATrp-tRNACys-tRNATyr and CR-tRNAIle-tRNAGln-tRNAMet were hotspots for gene rearrangement. Phylogenetic analysis showed that both Bayesian and maximum-likelihood methods recovered the monophyly of Aphidiinae and suggested that Aphidiinae formed sister clades with the remaining subfamilies. The phylogenetic analyses of nine subfamilies supported the monophyly of Cyclostomes and Noncyclostomes in Braconidae.

CONCLUSION

The arrangement of mitochondrial genes and the phylogenetic relationships among nine Braconidae subfamilies were constructed better to understand the diversity and evolution of Aphidiinae mitogenomes.

Comparative Mitogenomic Analysis of Two Snake Eels Reveals Irregular Gene Rearrangement and Phylogenetic Implications of Ophichthidae

The family Ophichthidae has the largest number and the most various species (about 359 valid species) in the order Anguilliformes worldwide. Both morphological and molecular characteristics have been used to assess their taxonomic status. However, due to the ambiguous morphological features, molecular data such as mitochondrial DNA sequences have been implemented for the correct identification and classification of these fishes. In this study, the gene arrangement and structure characteristics of two Ophichthidae mitochondrial genomes were investigated for the first time. The total mitogenome lengths of O. evermanni and O. erabo were 17,759 bp and 17,856 bp, respectively. Comparing with the ancestral mitochondrial gene order, the irregular gene rearrangement happened between ND6 and tRNA-Pro (P) genes with another similar control region emerging between tRNA-Thr (T) and ND6 genes, which could be explained by the tandem duplication and random loss (TDRL) model appropriately. ML phylogenetic tree demonstrated that the family Ophichthidae was monophyletic origin, but genus Ophichthus might be polyphyletic because of the confused cluster relationships among different species.

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.

Full-Length Transcriptome Profiling of Coridius chinensis Mitochondrial Genome Reveals the Transcription of Genes with Ancestral Arrangement in Insects

Coridius chinensis (Hemiptera: Dinidoridae) is a medicinal insect. Its mitochondrial gene arrangement is consistent with that of Drosophila melanogaster and Erthesina fullo, the two insects with well-studied mitochondrial transcription. To investigate whether the structural consistency of mitochondrial genes leads to similarities in transcription and post-transcriptional processing, we improved the gene annotation and constructed a quantitative transcription map for the C. chinensis mitochondrial genome (mitogenome) using full-length transcriptome sequencing. The size of this mitogenome was 16,214 bp and the proposed model of mitochondrial transcription was similar to that of Drosophila. Both strands were nearly entirely transcribed except for the antisense genes downstream of trnS2 on N strand. The expression of cytochrome c subunit genes is higher than that of NADH-dehydrogenase subunit genes. The post-transcriptional cleavage process followed the "tRNA punctuation" model, and both the "reverse cleavage" model in Drosophila and "forward cleavage" model in E. fullo were found in C. chinensis. In addition, we found that long non-coding RNAs from the control region contained tandem repeats. Polyadenylation was performed after CCA triplet at the 3' end of tRNA. The isoform diversity of lrRNA was identified. Our study sheds light on the transcriptional regulation and RNA processing of insect mitogenomes with the putative ancestral gene arrangement.

A comparison of gene organisations and phylogenetic relationships of all 22 squamate species listed in South Korea using complete mitochondrial DNA

Studies using complete mitochondrial genome data have the potential to increase our understanding on gene organisations and evolutionary species relationships. In this study, we compared complete mitochondrial genomes between all 22 squamate species listed in South Korea. In addition, we constructed Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI) phylogenetic trees using 13 mitochondrial protein-coding genes. The mitochondrial genes for all six species in the suborder Sauria followed the same organisation as the sequenced Testudines (turtle) outgroup. In contrast, 16 snake species in the suborder Serpentes contained some gene organisational variations. For example, all snake species contained a second control region (CR2), while three species in the family Viperidae had a translocated tRNA-Pro gene region. In addition, the snake species, Elapheschrenckii, carried a tRNA-Pro pseudogene. We were also able to identify a translocation of a tRNA-Asn gene within the five tRNA (WANCY gene region) gene clusters for two true sea snake species in the subfamily Hydrophiinae. Our BI phylogenetic tree was also well fitted against currently known Korean squamate phylogenetic trees, where each family and genus unit forms monophyletic clades and the suborder Sauria is paraphyletic to the suborder Serpentes. Our results may form the basis for future northeast Asian squamate phylogenetic studies.

Mitochondrial genomes of 10 Mantidae species and their phylogenetic implications

This article aims to present a phylogenetic evaluation of Mantidae based on a mitochondrial genome (mitogenome) data set. The mitogenome of 10 Mantidae species were sequenced using next-generation sequencing. The length of nine the complete mitogenomes ranged from 15,371 bp in Tenodera aridifolia to 16,063 bp in Hierodula longa. Mantidae mitogenomes have 37 genes and control region with two exceptions: five trnR copies in Statilia maculata, and H. zhangi was incomplete missing trnI, trnQ, trnM and a portion of the control region. There was a large noncoding region (LNC) between trnM and nad2 in H. chinensis, H. longa, H. maculata and Titanodula sp. Most of protein-coding genes (PCGs) used the typical start ATN codon and TAA/TAG stop codons. All tRNAs fold into the typical clover-leaf secondary structure except trnS1 which lacks a dihydrouracil (DHU) arm. Nucleotide diversity and Ka/Ks analysis of 13 PCGs showed that atp8 had the highest variability and fastest evolutionary rate. Phylogenetic relationships among 42 Mantidae species were reconstructed using the 13 PCGs and two rRNA genes using Bayesian Inference (BI) and Maximum Likelihood (ML) methods. Of the seven mantid subfamilies included in this analysis, only four had multiple exemplars, and of those only Mantinae and Vatinae formed monophyletic groups in BI and ML trees. Consistent with previous studies, the monophyly of the Hierudulinae and Tenoderinae were not been supported. The present results imply that it is necessary to combine nuclear molecular markers and external characteristic to understand the phylogenetic relationships within Mantidae.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

Background

Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.

Results

P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.

Conclusions

Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08607-4.

Novel duplication remnant in the first complete mitogenome of Hemitriakis japanica and the unique phylogenetic position of family Triakidae

In this study, we firstly determined the complete mitogenome of the Japanese topeshark (Hemitriakis japonica), which belong to the family Triakidae and was assessed as Endangered A2d on the IUCN Red List in 2021. The mitogenome is 17,301 bp long, has a high AT content (60.0%), and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, a control region and specially a 594 bp-long non-coding region between Cytb gene and tRNA-Thr gene. The novel non-coding region share high sequence similarity with segments of the former and latter genes, so it was recognized as a duplication remnant. In addition, the Cytb gene and tRNA-Thr gene tandemly duplicated twice while accompanied by being deleted once at least. This is the first report of mitogenomic gene-arrangement in Triakidae. The phylogenetic trees were constructed using Bayesian inference (BI) and maximum likelihood (ML) methods based on the mitogenomic data of 51 shark species and two outgroups. In summary, basing on a novel type of gene rearrangements in houndshark mitogenome, the possibly rearranged process was analyzed and contributed further insight of shark mitogenomes evolution and phylogeny.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

Mitogenomes of Accipitriformes and Cathartiformes Were Subjected to Ancestral and Recent Duplications Followed by Gradual Degeneration

The rearrangement of 37 genes with one control region, firstly identified in Gallus gallus mitogenome, is believed to be ancestral for all Aves. However, mitogenomic sequences obtained in recent years revealed that many avian mitogenomes contain duplicated regions that were omitted in previous genomic versions. Their evolution and mechanism of duplication are still poorly understood. The order of Accipitriformes is especially interesting in this context because its representatives contain a duplicated control region in various stages of degeneration. Therefore, we applied an appropriate PCR strategy to look for duplications within the mitogenomes of the early diverged species Sagittarius serpentarius and Cathartiformes, which is a sister order to Accipitriformes. The analyses revealed the same duplicated gene order in all examined taxa and the common ancestor of these groups. The duplicated regions were subjected to gradual degeneration and homogenization during concerted evolution. The latter process occurred recently in the species of Cathartiformes as well as in the early diverged lineages of Accipitriformes, that is, Sagittarius serpentarius and Pandion haliaetus. However, in other lineages, that is, Pernis ptilorhynchus, as well as representatives of Aegypiinae, Aquilinae, and five related subfamilies of Accipitriformes (Accipitrinae, Circinae, Buteoninae, Haliaeetinae, and Milvinae), the duplications were evolving independently for at least 14-47 Myr. Different portions of control regions in Cathartiformes showed conflicting phylogenetic signals indicating that some sections of these regions were homogenized at a frequency higher than the rate of speciation, whereas others have still evolved separately.

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

Mitochondrial genomes of two Polydora (Spionidae) species provide further evidence that mitochondrial architecture in the Sedentaria (Annelida) is not conserved

Contrary to the early evidence, which indicated that the mitochondrial architecture in one of the two major annelida clades, Sedentaria, is relatively conserved, a handful of relatively recent studies found evidence that some species exhibit elevated rates of mitochondrial architecture evolution. We sequenced complete mitogenomes belonging to two congeneric shell-boring Spionidae species that cause considerable economic losses in the commercial marine mollusk aquaculture: Polydora brevipalpa and Polydora websteri. The two mitogenomes exhibited very similar architecture. In comparison to other sedentarians, they exhibited some standard features, including all genes encoded on the same strand, uncommon but not unique duplicated trnM gene, as well as a number of unique features. Their comparatively large size (17,673 bp) can be attributed to four non-coding regions larger than 500 bp. We identified an unusually large (putative) overlap of 14 bases between nad2 and cox1 genes in both species. Importantly, the two species exhibited completely rearranged gene orders in comparison to all other available mitogenomes. Along with Serpulidae and Sabellidae, Polydora is the third identified sedentarian lineage that exhibits disproportionally elevated rates of mitogenomic architecture rearrangements. Selection analyses indicate that these three lineages also exhibited relaxed purifying selection pressures.

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'.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

The First Complete Mitochondrial Genome of Lachninae Species and Comparative Genomics Provide New Insights into the Evolution of Gene Rearrangement and the Repeat Region

Complete mitochondrial genomes are valuable resources for different research fields such as genomics, molecular evolution and phylogenetics. The subfamily Lachninae represents one of the most ancient evolutionary lineages of aphids. To date, however, no complete Lachninae mitogenome is available in public databases. Here we report the Stomaphis sinisalicis mitogenome, representing the first complete mitogenome of Lachninae. The S. sinisalicis mitogenome is consist of 13 protein-coding genes (PCGs), two rRNA genes (rRNAs), 22 tRNA genes (tRNAs), a control region and a large tandem repeat region. Strikingly, the mitogenome exhibits a novel, highly rearranged gene order between trnE and nad1 compared with that of other aphids. The presence of repeat region in the basal Lachninae may further indicate it is probably an ancestral feature of aphid mitogenomes. Collectively, this study provides new insights on mitogenome evolution and valuable data for future comparative studies across different insect lineages.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Laelia suffusa (Lepidoptera: Erebidae, Lymantriinae)

In this study, the complete mitochondrial genome of a white tussock moth, Laelia suffusa (Walker, 1855) (Lepidoptera: Erebidae, Lymantriinae), was sequenced and annotated. The genome sequence was 15,502 bp in length and comprised 13 PCGs, 2 rRNAs, 22 tRNAs, and a single noncoding control region (CR). The nucleotide composition of the genome was highly A + T biased, accounting for 79.04% of the whole genome and with a slightly positive AT skewness (0.015). Comparing the gene order with the basal species of Lepidoptera, a typical trnM rearrangement was detected in the mitogenome of L. suffusa. Besides, the trnM rearrangement was found at the head of trnI and trnQ, rather than at the back. The 13 PCGs used ATN as their start codons, except for the cox1 which used CGA. Out of the 22 tRNAs, only 1 tRNA (trnS1) failed to fold in a typical cloverleaf secondary structure. The conserved motif 'ATAGA + poly-T' was detected at the start of the control region which was similar to other Lepidoptera species. In total, 10 overlapping regions and 19 intergenic spacers were identified, ranging from 1 to 41 and 2 to 73 bp, respectively. Phylogenetic analysis showed that Lymantriinae was a monophyletic group with a high support value and L. suffusa was closely related to tribe Orgyiini (Erebidae, Lymantriinae). Moreover, the phylogenetic relationship of Noctuoidea (Lepidoptera) species was reconstructed using two datasets (13 PCGs and 37 genes) and these supported the topology of (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))).

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Rearrangement and evolution of mitochondrial genomes in Thysanoptera (Insecta)

Prior to this study, complete mitochondrial genomes from Order Thysanoptera were restricted to a single family, the Thripidae, resulting in a biased view of their evolution. Here we present the sequences for the mitochondrial genomes of four additional thrips species, adding three extra families and an additional subfamily, thus greatly improving taxonomic coverage. Thrips mitochondrial genomes are marked by high rates of gene rearrangement, duplications of the control region and tRNA mutations. Derived features of mitochondrial tRNAs in thrips include gene duplications, anticodon mutations, loss of secondary structures and high gene translocation rates. Duplicated control regions are found in the Aeolothripidae and the 'core' Thripinae clade but do not appear to promote gene rearrangement as previously proposed. Phylogenetic analysis of thrips mitochondrial sequence data supports the monophyly of two suborders, a sister-group relationship between Stenurothripidae and Thripidae, and suggests a novel set of relationships between thripid genera. Ancestral state reconstructions indicate that genome rearrangements are common, with just eight gene blocks conserved between any thrips species and the ancestral insect mitochondrial genome. Conversely, 71 derived rearrangements are shared between at least two species, and 24 of these are unambiguous synapomorphies for clades identified by phylogenetic analysis. While the reconstructed sequence of genome rearrangements among the protein-coding and ribosomal RNA genes could be inferred across the phylogeny, direct inference of phylogeny from rearrangement data in MLGO resulted in a highly discordant set of relationships inconsistent with both sequence-based phylogenies and previous morphological analysis. Given the demonstrated rates of genomic evolution within thrips, extensive sampling is needed to fully understand these phenomena across the order.

Mesozoic mitogenome rearrangements and freshwater mussel (Bivalvia: Unionoidea) macroevolution

Using a new fossil-calibrated mitogenome-based approach, we identified macroevolutionary shifts in mitochondrial gene order among the freshwater mussels (Unionoidea). We show that the early Mesozoic divergence of the two Unionoidea clades, Margaritiferidae and Unionidae, was accompanied by a synchronous split in the gene arrangement in the female mitogenome (i.e., gene orders MF1 and UF1). Our results suggest that this macroevolutionary jump was completed within a relatively short time interval (95% HPD 201-226 Ma) that coincided with the Triassic-Jurassic mass extinction. Both gene orders have persisted within these clades for ~200 Ma. The monophyly of the so-called "problematic" Gonideinae taxa was supported by all the inferred phylogenies in this study using, for the first time, the M- and F-type mitogenomes either singly or combined. Within Gonideinae, two additional splits in the gene order (UF1 to UF2, UF2 to UF3) occurred in the Mesozoic and have persisted for ~150 and ~100 Ma, respectively. Finally, the mitogenomic results suggest ancient connections between freshwater basins of East Asia and Europe near the Cretaceous-Paleogene boundary, probably via a continuous paleo-river system or along the Tethys coastal line, which are well supported by at least three independent but almost synchronous divergence events.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

The complete mitochondrial genome of Pyxicephalus adspersus: high gene rearrangement and phylogenetics of one of the world's largest frogs

The family Pyxicephalidae including two subfamilies (Cacosterninae and Pyxicephalinae) is an ecologically important group of frogs distributed in sub-Saharan Africa. However, its phylogenetic position among the Anura has remained uncertain. The present study determined the complete mitochondrial genome sequence of Pyxicephalus adspersus, the first representative mitochondrial genome from the Pyxicephalinae, and reconstructed the phylogenetic relationships within Ranoidae using 10 mitochondrial protein-coding genes of 59 frog species. The P. adspersus mitochondrial genome showed major gene rearrangement and an exceptionally long length that is not shared with other Ranoidae species. The genome is 24,317 bp in length, and contains 15 protein-coding genes (including extra COX3 and Cyt b genes), four rRNA genes (including extra 12S rRNA and 16S rRNA genes), 29 tRNA genes (including extra tRNALeu (UAG), tRNALeu (UUR), tRNAThr , tRNAPro , tRNAPhe , tRNAVal , tRNAGln genes) and two control regions (CRs). The Dimer-Mitogenome and Tandem duplication and random loss models were used to explain these gene arrangements. Finally, both Bayesian inference and maximum likelihood analyses supported the conclusion that Pyxicephalidae was monophyletic and that Pyxicephalidae was the sister clade of (Petropedetidae + Ptychadenidae).

Mitochondrial genome reorganization provides insights into the relationship between oribatid mites and astigmatid mites (Acari: Sarcoptiformes: Oribatida)

Oribatida s.l. represents one of the most species-rich mite lineages, including two recognized groups: oribatid mites (Oribatida s.s., non-astigmatan oribatids) and astigmatid mites (Astigmata). However, the relationship between these two groups has been debated. Here, we sequenced the complete mitochondrial (mt) genome of one oribatid mite and one astigmatid mite, retrieved complete mt genomes of three oribatid mites, and compared them with two other oribatid mites and 12 astigmatid mites sequenced previously. We find that gene orders in the mt genomes of both oribatid mites and astigmatid mites are rearranged relative to the hypothetical ancestral arrangement of the arthropods. Based on the shared derived gene clusters in each mt genome group, rearranged mt genomes are roughly divided into two groups corresponding to each mite group (oribatid mites or astigmatid mites). Phylogenetic results show that Astigmata nested in Oribatida. The monophyly of Astigmata is recovered, while paraphyly of Oribatida s.s. is observed. Our results show that rearranged gene orders in the mt genomes characterize various lineages of oribatid mites and astigmatid mites, and have potential phylogenetic information for resolving the high-level (cohort or supercohort) phylogeny of Oribatida.

Characterization of the complete mitochondrial genome of Brentisentisyangtzensis Yu & Wu, 1989 (Acanthocephala, Illiosentidae)

The mitogenome of Brentisentisyangtzensis is 13,864 bp in length and has the circular structure typical of metazoans. It contains 36 genes: 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and 12 protein-encoding genes (PCGs). All genes are transcribed from the same strand. Thirteen overlapping regions were found in the mitochondrial genome. The overall A+T content of B.yangtzensis is 68.3% versus 31.7% of G+C content (A = 27.8%, T = 40.5%, C = 9.0%, G = 22.7%). B.yangtzenensis (Illiosentidae) and Leptorhynchoidesthecatus (Rhadinorhynchidae) form a sister clade, showing the relatively close relationship between the Illiosentidae and the Rhadinorhynchidae. The mitochondrial gene arrangements of acanthocephalan species are relatively conserved, with only a few translocations of tRNAs (trnS1, trnS2, trnV, and trnK) detected. An identical gene order was found both in a sister clade (Centrorhynchusaluconis and Plagiorhynchustransversus) and across different classes (B.yangtzensis (Palaeacanthocephala), Acanthosentischeni (Eoacanthocephala) and Macracanthorhynchushirudinaceus (Archiacanthocephala), Oncicolaluehei and L.thecatus (Palaeacanthocephala)). More studies and more sequences of acanthocephalan species are needed to gain a clear understanding of the phylogenetic relationships.

High and Variable Rates of Repeat-Mediated Mitochondrial Genome Rearrangement in a Genus of Plants

For 30 years, it has been clear that angiosperm mitochondrial genomes evolve rapidly in sequence arrangement (i.e., synteny), yet absolute rates of rearrangement have not been measured in any plant group, nor is it known how much these rates vary. To investigate these issues, we sequenced and reconstructed the rearrangement history of seven mitochondrial genomes in Monsonia (Geraniaceae). We show that rearrangements (occurring mostly as inversions) not only take place at generally high rates in these genomes but also uncover significant variation in rearrangement rates. For example, the hyperactive mitochondrial genome of Monsonia ciliata has accumulated at least 30 rearrangements over the last million years, whereas the branch leading to M. ciliata and its sister species has sustained rearrangement at a rate that is at least ten times lower. Furthermore, our analysis of published data shows that rates of mitochondrial genome rearrangement in seed plants vary by at least 600-fold. We find that sites of rearrangement are highly preferentially located in very close proximity to repeated sequences in Monsonia. This provides strong support for the hypothesis that rearrangement in angiosperm mitochondrial genomes occurs largely through repeat-mediated recombination. Because there is little variation in the amount of repeat sequence among Monsonia genomes, the variable rates of rearrangement in Monsonia probably reflect variable rates of mitochondrial recombination itself. Finally, we show that mitochondrial synonymous substitutions occur in a clock-like manner in Monsonia; rates of mitochondrial substitutions and rearrangements are therefore highly uncoupled in this group.

Gene annotation errors are common in the mammalian mitochondrial genomes database

Background

Although animal mitochondrial DNA sequences are known to evolve rapidly, their gene arrangements often remain unchanged over long periods of evolutionary time. Therefore, comparisons of mitochondrial genomes may result in significant insights into the evolution both of organisms and of genomes. Mammalian mitochondrial genomes recently published in the GenBank database of NCBI show numerous rearrangements in various regions of the genome, from which it may be inferred that the mammalian mitochondrial genome is more dynamic than expected. However, it is alternatively possible that these are errors of annotation and, if so, are misleading our interpretations. In order to verify these possible errors of annotation, we performed a comparative genomic analysis of mammalian mitochondrial genomes available in the NCBI database.

Results

Using a combination of bioinformatics methods to carefully examine the mitochondrial gene arrangements in 304 mammalian species, we determined that there are only two sets of gene arrangements, one that is shared by all of the marsupials and another that is shared by all of the monotremes and eutherians, with these two arrangements differing only by the positions of tRNA genes in the region commonly designated as “WANCY” for the genes it comprises. All of the 68 other cases of reported gene rearrangements are errors. We note that there are also numerous errors of impossibly short, incorrect gene annotations, cases where genomes that are reported as complete are actually missing portions of the sequence, and genes that are clearly present but were not annotated in these records.

Conclusions

We judge that the application of simple bioinformatic tools in the verification of gene annotation, particularly for organelle genomes, would be a very useful enhancement for the curation of genome sequences submitted to GenBank.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5447-1) contains supplementary material, which is available to authorized users.

Mitochondrial genomes of two diplectanids (Platyhelminthes: Monogenea) expose paraphyly of the order Dactylogyridea and extensive tRNA gene rearrangements

BACKGROUND

Recent mitochondrial phylogenomics studies have reported a sister-group relationship of the orders Capsalidea and Dactylogyridea, which is inconsistent with previous morphology- and molecular-based phylogenies. As Dactylogyridea mitochondrial genomes (mitogenomes) are currently represented by only one family, to improve the phylogenetic resolution, we sequenced and characterized two dactylogyridean parasites, Lamellodiscus spari and Lepidotrema longipenis, belonging to a non-represented family Diplectanidae.

RESULTS

The L. longipenis mitogenome (15,433 bp) contains the standard 36 flatworm mitochondrial genes (atp8 is absent), whereas we failed to detect trnS1, trnC and trnG in L. spari (14,614 bp). Both mitogenomes exhibit unique gene orders (among the Monogenea), with a number of tRNA rearrangements. Both long non-coding regions contain a number of different (partially overlapping) repeat sequences. Intriguingly, these include putative tRNA pseudogenes in a tandem array (17 trnV pseudogenes in L. longipenis, 13 trnY pseudogenes in L. spari). Combined nucleotide diversity, non-synonymous/synonymous substitutions ratio and average sequence identity analyses consistently showed that nad2, nad5 and nad4 were the most variable PCGs, whereas cox1, cox2 and cytb were the most conserved. Phylogenomic analysis showed that the newly sequenced species of the family Diplectanidae formed a sister-group with the Dactylogyridae + Capsalidae clade. Thus Dactylogyridea (represented by the Diplectanidae and Dactylogyridae) was rendered paraphyletic (with high statistical support) by the nested Capsalidea (represented by the Capsalidae) clade.

CONCLUSIONS

Our results show that nad2, nad5 and nad4 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in the Diplectanidae. The unique gene order pattern further suggests discontinuous evolution of mitogenomic gene order arrangement in the Class Monogenea. This first report of paraphyly of the Dactylogyridea highlights the need to generate more molecular data for monogenean parasites, in order to be able to clarify their relationships using large datasets, as single-gene markers appear to provide a phylogenetic resolution which is too low for the task.

Genome Rearrangement with ILP

The weighted Genome Sorting Problem (wGSP) is to find a minimum-weight sequence of rearrangement operations that transforms a given gene order into another given gene order using rearrangement operations that are associated with a predefined weight. This paper presents a polynomial sized Integer Linear Program -called GeRe-ILP- for solving the wGSP for the following three types of rearrangement operations: inversion , transposition, and inverse transposition. GeRe-ILP uses variables and constraints for gene orders of length . It is studied experimentally on simulated data how different weighting schemes influence the reconstructed scenarios. The influences of the length of the gene orders and of the size of the reconstructed scenarios on the runtime of GeRe-ILP are studied as well.

The mitochondrial genomes of the mesozoans Intoshia linei, Dicyema sp. and Dicyema japonicum

The Dicyemida and Orthonectida are two groups of tiny, simple, vermiform parasites that have historically been united in a group named the Mesozoa. Both Dicyemida and Orthonectida have just two cell layers and appear to lack any defined tissues. They were initially thought to be evolutionary intermediates between protozoans and metazoans but more recent analyses indicate that they are protostomian metazoans that have undergone secondary simplification from a complex ancestor. Here we describe the first almost complete mitochondrial genome sequence from an orthonectid, Intoshia linei, and describe nine and eight mitochondrial protein-coding genes from Dicyema sp. and Dicyema japonicum, respectively. The 14,247 base pair long I. linei sequence has typical metazoan gene content, but is exceptionally AT-rich, and has a unique gene order. The data we have analysed from the Dicyemida provide very limited support for the suggestion that dicyemid mitochondrial genes are found on discrete mini-circles, as opposed to the large circular mitochondrial genomes that are typical of the Metazoa. The cox1 gene from dicyemid species has a series of conserved, in-frame deletions that is unique to this lineage. Using cox1 genes from across the genus Dicyema, we report the first internal phylogeny of this group.

The complete mitochondrial genome of Vanessa indica and phylogenetic analyses of the family Nymphalidae

Vanessa indica is a small butterfly lacking historical molecular and biological research. Vanessa indica belongs to the family Nymphalidae (Lepidoptera: Papilionoidea), which is the largest group of butterflies and are nearly ubiquitous. However, after more than a century of taxonomic and molecular studies, there is no consensus for family classification, and the phylogenetic relationships within Nymphalidae are controversial. The first objective was to sequence and characterize the complete mitochondrial genome of V. indica. The most important objective was to completely reconstruct the phylogenetic relationships for family members within Nymphalidae. The mitochondrial genomic DNA (mtDNA) of V. indica was extracted and amplified by polymerase chain reaction. The complete mitochondrial sequence was annotated and characterized by analyzing sequences with SeqMan program. The phylogenetic analyses were conducted on thirteen protein coding genes (PCGs) in 95 mtDNA of Nymphalidae downloaded from GenBank for reference using the maximum likelihood method and Bayesian inference to ensure the validity of the results. The complete mitogenome was a circular molecule with 15,191 bp consisting of 13 protein coding genes, two ribosomal RNA genes (16S rRNA and 12S rRNA), 22 transfer RNA (tRNA) genes, and an A + T-rich region (D-loop). The nucleotide composition of the genome was highly biased for A + T content, which accounts for 80.0% of the nucleotides. All the tRNAs have putative secondary structures that are characteristic of mitochondrial tRNAs, except tRNA^Ser(AGN). All the PCGs started with ATN codons, except cytochrome c oxidase subunit 1 (COX1), which was found to start with an unusual CGA codon. Four genes were observed to have unusual codons: COX1 terminated with atypical TT and the other three genes terminated with a single T. The A + T rich region of 327 bp consisted of repetitive sequences, including a ATAGA motif, a 19-bp poly-T stretch, and two microsatellite-like regions (TA)₈. The phylogenetic analyses consistently placed Biblidinae as a sister cluster to Heliconiinae and Calinaginae as a sister clade to Satyrinae. Moreover, the phylogenetic tree identified Libytheinae as a monophyletic group within Nymphalidae. The complete mitogenome of V. indica was 15,191 bp with mitochondrial characterizations common for lepidopteran species, which enriched the mitochondria data of Nymphalid species. And the phylogenetic analysis revealed different classifications and relationships than those previously described. Our results are significant because they would be useful in further understanding of the evolutionary biology of Nymphalidae.

Multiple independent structural dynamic events in the evolution of snake mitochondrial genomes

BACKGROUND

Mitochondrial DNA sequences have long been used in phylogenetic studies. However, little attention has been paid to the changes in gene arrangement patterns in the snake's mitogenome. Here, we analyzed the complete mitogenome sequences and structures of 65 snake species from 14 families and examined their structural patterns, organization and evolution. Our purpose was to further investigate the evolutionary implications and possible rearrangement mechanisms of the mitogenome within snakes.

RESULTS

In total, eleven types of mitochondrial gene arrangement patterns were detected (Type I, II, III, III-A, III-B, III-B1, III-C, III-D, III-E, III-F, III-G), with mitochondrial genome rearrangements being a major trend in snakes, especially in Alethinophidia. In snake mitogenomes, the rearrangements mainly involved three processes, gene loss, translocation and duplication. Within Scolecophidia, the O_L was lost several times in Typhlopidae and Leptotyphlopidae, but persisted as a plesiomorphy in the Alethinophidia. Duplication of the control region and translocation of the tRNA^Leu gene are two visible features in Alethinophidian mitochondrial genomes. Independently and stochastically, the duplication of pseudo-Pro (P*) emerged in seven different lineages of unequal size in three families, indicating that the presence of P* was a polytopic event in the mitogenome.

CONCLUSIONS

The WANCY tRNA gene cluster and the control regions and their adjacent segments were hotspots for mitogenome rearrangement. Maintenance of duplicate control regions may be the source for snake mitogenome structural diversity.

Probing recalcitrant problems in polyclad evolution and systematics with novel mitochondrial genome resources

For their apparent morphological simplicity, the Platyhelminthes or "flatworms" are a diverse clade found in a broad range of habitats. Their body plans have however made them difficult to robustly classify. Molecular evidence is only beginning to uncover the true evolutionary history of this clade. Here we present nine novel mitochondrial genomes from the still undersampled orders Polycladida and Rhabdocoela, assembled from short Illumina reads. In particular we present for the first time in the literature the mitochondrial sequence of a Rhabdocoel, Bothromesostoma personatum (Typhloplanidae, Mesostominae). The novel mitochondrial genomes examined generally contained the 36 genes expected in the Platyhelminthes, with all possessing 12 of the 13 protein-coding genes normally found in metazoan mitochondrial genomes (ATP8 being absent from all Platyhelminth mtDNA sequenced to date), along with two ribosomal RNA genes. The majority presented possess 22 transfer RNA genes, and a single tRNA gene was absent from two of the nine assembled genomes. By comparison of mitochondrial gene order and phylogenetic analysis of the protein coding and ribosomal RNA genes contained within these sequences with those of previously sequenced species we are able to gain a firm molecular phylogeny for the inter-relationships within this clade. Our phylogenetic reconstructions, using both nucleotide and amino acid sequences under several models and both Bayesian and Maximum Likelihood methods, strongly support the monophyly of Polycladida, and the monophyly of Acotylea and Cotylea within that clade. They also allow us to speculate on the early emergence of Macrostomida, the monophyly of a "Turbellarian-like" clade, the placement of Rhabditophora, and that of Platyhelminthes relative to the Lophotrochozoa (=Spiralia). The data presented here therefore represent a significant advance in our understanding of platyhelminth phylogeny, and will form the basis of a range of future research in the still-disputed classifications within this taxon.

Combinatorics of Tandem Duplication Random Loss Mutations on Circular Genomes

The tandem duplication random loss operation (TDRL) is an important genome rearrangement operation in metazoan mitochondrial genomes. A TDRL consists of a duplication of a contiguous set of genes in tandem followed by a random loss of one copy of each duplicated gene. This paper presents an analysis of the combinatorics of TDRLs on circular genomes, e.g., the mitochondrial genome. In particular, results on TDRLs for circular genomes and their linear representatives are established. Moreover, the distance between gene orders with respect to linear TDRLs and circular TDRLs is studied. An analysis of the available animal mitochondrial gene orders shows the practical relevance of the theoretical results.

Genome Rearrangement Analysis: Cut and Join Genome Rearrangements and Gene Cluster Preserving Approaches

Genome rearrangements are mutations that change the gene content of a genome or the arrangement of the genes on a genome. Several years of research on genome rearrangements have established different algorithmic approaches for solving some fundamental problems in comparative genomics based on gene order information. This review summarizes the literature on genome rearrangement analysis along two lines of research. The first line considers rearrangement models that are particularly well suited for a theoretical analysis. These models use rearrangement operations that cut chromosomes into fragments and then join the fragments into new chromosomes. The second line works with rearrangement models that reflect several biologically motivated constraints, e.g., the constraint that gene clusters have to be preserved. In this chapter, the border between algorithmically "easy" and "hard" rearrangement problems is sketched and a brief review is given on the available software tools for genome rearrangement analysis.

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders

Background

Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication.

Results

We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years.

Conclusions

The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.

Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera

BACKGROUND

Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae.

RESULTS

We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny.

CONCLUSIONS

This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.