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

The complete mitochondrial genome of Calyptogena marissinica (Heterodonta: Veneroida: Vesicomyidae): Insight into the deep-sea adaptive evolution of vesicomyids

The deep-sea chemosynthetic environment is one of the most extreme environments on the Earth, with low oxygen, high hydrostatic pressure and high levels of toxic substances. Species of the family Vesicomyidae are among the dominant chemosymbiotic bivalves found in this harsh habitat. Mitochondria play a vital role in oxygen usage and energy metabolism; thus, they may be under selection during the adaptive evolution of deep-sea vesicomyids. In this study, the mitochondrial genome (mitogenome) of the vesicomyid bivalve Calyptogena marissinica was sequenced with Illumina sequencing. The mitogenome of C. marissinica is 17,374 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes. All of these genes are encoded on the heavy strand. Some special elements, such as tandem repeat sequences, “G(A)_nT” motifs and AT-rich sequences, were observed in the control region of the C. marissinica mitogenome, which is involved in the regulation of replication and transcription of the mitogenome and may be helpful in adjusting the mitochondrial energy metabolism of organisms to adapt to the deep-sea chemosynthetic environment. The gene arrangement of protein-coding genes was identical to that of other sequenced vesicomyids. Phylogenetic analyses clustered C. marissinica with previously reported vesicomyid bivalves with high support values. Positive selection analysis revealed evidence of adaptive change in the mitogenome of Vesicomyidae. Ten potentially important adaptive residues were identified, which were located in cox1, cox3, cob, nad2, nad4 and nad5. Overall, this study sheds light on the mitogenomic adaptation of vesicomyid bivalves that inhabit the deep-sea chemosynthetic environment.

Complete mitochondrial genome of two Thitarodes species (Lepidoptera, Hepialidae), the host moths of Ophiocordyceps sinensis and phylogenetic implications

Thitarodes (Lepidoptera, Hepialidae) is the only genus that hosts to the Ophiocordyceps sinensis, a traditional Chinese medicine considered as a powerful medicinal supplement. In this study, the complete mitochondrial genomes (mitogenomes) of two species, T. damxungensis and T. pui, have been sequenced, which are 15,928 bp and 15,362 bp in size respectively, and both contain 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and an AT-rich region. Like other hepialoids, the gene arrangement of the mitogenomes of T. damxungensis and T. pui is identical to the ancestral arrangement but differs from those of other lepidopteran species on account of the different arrangements of trnM, trnI, and trnQ. The size of AT-rich region is 545 bp in T. damxungensis and 1030 bp in T. pui. Tandem repetition in the AT-rich region is responsible for the length difference of the A + T-rich region in both species. In Hepialidae, the phylogenetic study based on the dataset of the sequences that combined the protein-coding genes and RNA genes suggested that the species T. yunnanensis should still belong to the genus Thitarodes rather than Ahamns, which is different from the results based on the traditional phylogeny.

Complete mitochondrial genomes of Bittacus strigosus and Panorpa debilis and genomic comparisons of Mecoptera

Mitochondrial genomes play a significant role in reconstructing phylogenetic relationships and revealing molecular evolution of insects. However, only four mitochondrial genomes were reported in Mecoptera to date. Here, we obtained two new complete mitochondrial genomes of the hangingfly Bittacus strigosus Hagen, 1861 and the scorpionfly Panorpa debilis Westwood, 1846. The results show that the complete mitogenome sequences of B. strigosus and P. debilis are 15,825 and 17,018 bp, respectively, both containing 37 genes and one control region. The mecopteran mitogenomes are highly similar in A + T bias, AT-skew, and GC-skew. Tandem repeats of the control region were discovered in Mecoptera for the first time. The sliding window, genetic distance, and Ka/Ks ratio analyses indicate the purifying selection of 13 protein-coding genes, the lowest evolutionary rate of cox1, and the highest sequence variability of atp8. Considering the sufficiently large size, fast evolution, and high ratio of Ka/Ks, nad4L and nad6 are regarded as potential markers for future phylogenetic analyses, population genetics, and species delimitations in Mecoptera. The phylogenetic relationships were reconstructed for four families of Mecoptera based on all six available mitogenomes using Bayesian inference and maximum likelihood methods. The phylogeny is presented as Boreidae + (Nannochoristidae + (Bittacidae + Panorpidae)).

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

The tick Haemaphysalis longicornis is widely distributed in eastern Asia, New Zealand and Australia, and is well-known as a vector of multiple zoonotic pathogens. This species exhibits two reproductive strategies, bisexual and obligate parthenogenetic reproduction. Hence, in the current study, the complete mitochondrial genomes of the bisexual and parthenogenetic populations were assembled and analyzed, and the expression of the mitochondrial protein-coding genes was evaluated and compared between the two reproductive populations. The results indicated that the length of the mitochondrial genomes of the two reproductive populations is 14,694 and 14,693 bp in the bisexual and parthenogenetic populations, respectively. The AT content in the mitochondrial genome of the bisexual and obligate parthenogenetic population reached 77.22 and 77.34%, respectively. The phylogenetic tree was constructed combining 13 protein-coding genes, which showed that the genetic distance between the bisexual and parthenogenetic populations was less than that between the subspecies. The expression of the mitochondrial protein-coding genes was quantitatively analyzed at different feeding status for the bisexual and parthenogenetic populations, and the results showed significant differences in the expression patterns of these genes, suggesting that they might trigger specific energy utilization mechanisms due to their different reproductive strategies and environmental pressures.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Mitochondrial phylogeny and comparative mitogenomics of closely related pine moth pests (Lepidoptera: Dendrolimus)

BACKGROUND

Pine moths, Dendrolimus spp. (Lasiocampidae), are serious economic pests of conifer forests. Six closely related species (Dendrolimus punctatus, D. tabulaeformis, D. spectabilis, D. superans, D. houi, and D. kikuchii) occur in China and cause serious damage to coniferophyte. The complete mito genomes of Dendrolimus genus are significant to resolve the phylogenetic relationship and provide theoretical support in pest control.

METHODS

The complete mitogenomes of three species (D. superans, D. houi, and D. kikuchii) were sequenced based on PCR-amplified with universal primers, which were used to amplify initial fragments. Phylogenetic analyses were carried out with 78 complete mitogenomes of lepidopteran species from 10 superfamilies.

RESULTS

The complete mitochondrial genomes of these three species were 15,417, 15,381, and 15,377 bp in length, separately. The phylogenetic analyses produced consistent results for six Dendrolimus species based on complete mitogenomes, two major clades were formed, one containing D. spectabilis clustered with D. punctatus + D. tabulaeformis, and D. superans as the sister group to this three-taxon clade, the other containing D. kikuchii and D. houi. Comparative analyses of the congeneric mitochondrial genomes were performed, which showed that non-coding regions were more variable than the A+T rich region. The mitochondrial nucleotide diversity was more variable when compared within than among genus, and the concatenated tRNA region was the most conserved and the nd6 genes was the most variable.

Comparative mitogenomic analysis of species in the subfamily Amphinemurinae (Plecoptera: Nemouridae) reveal conserved mitochondrial genome organization

The subfamily Amphinemurinae has five genera in China, with each genus of similar morphology. To gain a better understanding of architecture and evolution of mitogenome in Amphinemurinae, mitogenomes of eight species representing four genera (Amphinemura, Indonemoura, Protonemura and Sphaeronemoura) in the subfamily Amphinemurinae were sequenced, and a comparative mitogenomic analysis of five genera (including a published stonefly genus, Mesonemoura) was carried out. By comparative analysis, we found highly conserved genome organization of ten Amphinemurinae species including genome contents, gene order, nucleotide composition, codon usage, amino acid composition, as well as genome asymmetry. GC content was the most significant factor in determining codon bias among organisms. The Ka/Ks values for all PCGs were far lower than 1, indicating that these genes were evolving under purifying selection. We also found some important conserved stem and loop in the cloverleaf structure of tRNAs, and found conserved helices and loops in each domain of the secondary structure of rRNAs. The presence of structural elements in the control region is also discussed. The phylogenetic analyses indicated that within Amphinemurinae, Sphaeronemoura was assigned the sister group of Mesonemoura. Our analyses inferred a relationship within Euholognatha: ((Nemouridae + Notonemouridae) + (Taeniopterygidae + Capniidae) + Scopuridae) + Leuctridae.

Sequencing and analysis of the complete mitogenome of Doriprismatica atromarginata (Cuvier, 1804)

Doriprismatica atromarginata is a common nudibranch in the tropical and sub-tropical Indo-Pacific region. In this study, the complete mitogenome of D. atromarginata from Vietnam was recorded for the first time. The circular mitogenome had a size of 14,421 bp and consisted of 37 genes (13 protein-coding genes, two ribosomal RNA genes, and 22 tRNA genes). Phylogenetic tree based on amino acid sequences of coding genes demonstrated that D. atromarginata has sister group relationship with the genus Chromodoris.

OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes

Organellar (plastid and mitochondrial) genomes play an important role in resolving phylogenetic relationships, and next-generation sequencing technologies have led to a burst in their availability. The ongoing massive sequencing efforts require software tools for routine assembly and annotation of organellar genomes as well as their display as physical maps. OrganellarGenomeDRAW (OGDRAW) has become the standard tool to draw graphical maps of plastid and mitochondrial genomes. Here, we present a new version of OGDRAW equipped with a new front end. Besides several new features, OGDRAW now has access to a local copy of the organelle genome database of the NCBI RefSeq project. Together with batch processing of (multi-)GenBank files, this enables the user to easily visualize large sets of organellar genomes spanning entire taxonomic clades. The new OGDRAW server can be accessed at https://chlorobox.mpimp-golm.mpg.de/OGDraw.html.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

Long-read sequencing reveals a 4.4 kb tandem repeat region in the mitogenome of Echinococcus granulosus (sensu stricto) genotype G1

Background

Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method.

Methods

We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow.

Results

We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441–445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat.

Conclusions

We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such ‘new’ repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

Comparison of mitochondrial DNA enrichment and sequencing methods from fish tissue

Sparid fish species have different commercial values related to their organoleptic features. Mitochondrial (mt) DNA provides a potential tool to distinguish species, but the enrichment of high-quality mtDNA from total genomic DNA is critical to obtain entire mtDNA sequences. Conventional mtDNA isolation is relatively low-cost and proficient. However, high numbers of PCR cycles can lead to artefacts (10^-6 mutations/bp). We describe a rapid protocol for mtDNA extraction and enrichment from fish tissues, based on conventional miniprep columns and paramagnetic bead-based purification, without the need to employ PCR amplification. This newly described method generates a substrate for next-generation sequencing (NGS) analysis and is likely to have wider applications for mitochondrial studies in other fish families to help ensure traceability and differentiation of fish with high commercial values.

Mitochondrial genomes of the stoneflies Mesonemourametafiligera and Mesonemouratritaenia (Plecoptera, Nemouridae), with a phylogenetic analysis of Nemouroidea

In this study, two new mitochondrial genomes (mitogenomes) of Mesonemourametafiligera and Mesonemouratritaenia from the family Nemouridae (Insecta: Plecoptera) were sequenced. The Mesonemourametafiligera mitogenome was a 15,739 bp circular DNA molecule, which was smaller than that of M.tritaenia (15,778 bp) due to differences in the size of the A+T-rich region. Results show that gene content, gene arrangement, base composition, and codon usage were highly conserved in two species. Ka/Ks ratios analyses of protein-coding genes revealed that the highest and lowest rates were found in ND6 and COI and that all these genes were evolving under purifying selection. All tRNA genes in nemourid mitogenomes had a typical cloverleaf secondary structure, except for tRNA^Ser(AGN) which appeared to lack the dihydrouridine arm. The multiple alignments of nemourid lrRNA and srRNA genes showed that sequences of three species were highly conserved. All the A+T-rich region included tandem repeats regions and stem-loop structures. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results. Amphinemurinae and Nemourinae were sister-groups and the family Nemouridae was placed as sister to Capniidae and Taeniopterygidae.

Mitochondrial Genome Fragmentation Unites the Parasitic Lice of Eutherian Mammals

Organelle genome fragmentation has been found in a wide range of eukaryotic lineages; however, its use in phylogenetic reconstruction has not been demonstrated. We explored the use of mitochondrial (mt) genome fragmentation in resolving the controversial suborder-level phylogeny of parasitic lice (order Phthiraptera). There are approximately 5000 species of parasitic lice in four suborders (Amblycera, Ischnocera, Rhynchophthirina, and Anoplura), which infest mammals and birds. The phylogenetic relationships among these suborders are unresolved despite decades of studies. We sequenced the mt genomes of eight species of parasitic lice and compared them with 17 other species of parasitic lice sequenced previously. We found that the typical single-chromosome mt genome is retained in the lice of birds but fragmented into many minichromosomes in the lice of eutherian mammals. The shared derived feature of mt genome fragmentation unites the eutherian mammal lice of Ischnocera (family Trichodectidae) with Anoplura and Rhynchophthirina to the exclusion of the bird lice of Ischnocera (family Philopteridae). The novel clade, namely Mitodivisia, is also supported by phylogenetic analysis of mt genome and cox1 gene sequences. Our results demonstrate, for the first time, that organelle genome fragmentation is informative for resolving controversial high-level phylogenies.

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

Comparative analysis of the mitochondrial genomes of Orthonectida: insights into the evolution of an invertebrate parasite species

Among invertebrates, only a few groups still have uncertain phylogenetic position, Orthonectida, a small group of rare multi-cellular parasites of marine invertebrates, being one of them. Recent molecular and morphological findings suggest that orthonectids belong to Lophotrochozoa and are close to Annelida. Nevertheless, phylogenetic relationships between orthonectids and annelids are unclear, and the phylogeny within the group itself has never been studied. Sequencing of mitochondrial genomes is used here to clarify this issue. Complete mt genomes of the orthonectids Intoshia variabili and Rhopalura litoralis were characterized and compared with Intoshia linei mt genome. Our results show that Orthonectida mt genomes have undergone reduction and gene loss, and that they have complicated organization revealed in strand asymmetry in nucleotide composition, in some features of intergenic non-coding regions, tRNA duplication and folding. Moreover, all species of Orthonectida have a unique gene order with complicated rearrangement landscape. Significant differences in mitochondrial genomes in the three orthonectid species could be explained by the fact that their host species belong to different taxa (flat worms, nemertines and gastropods). Among the analyzed mt genomes of Orthonectida, I. linei possesses the closest gene order to the ancestral genome. All Orthonectida species are monophyletic, and in the phylogenetic tree are close to Pleistoannelida, and specifically, to Clitellata.

Incorporation of deep-sea and small-sized species provides new insights into gastropods phylogeny

The use of phylogeny with uneven or limited taxon sampling may bias our interpretation of organismal evolution, for instance, the origin(s) of the deep-sea animals. The Mollusca is the second most speciose phylum, in which the Gastropoda forms the largest group. However, the currently proposed hypotheses of gastropod phylogeny are mainly based on part of their taxonomic diversity, notably on the large-sized and shallow-water species. In this study, we aimed at correcting this bias by reconstructing the phylogeny with new mitogenomes of deep-sea gastropods including Anatoma sp., Bathysciadiidae sp., Bayerotrochus teramachii, Calliotropis micraulax, Coccocrater sp., Cocculina subcompressa, Lepetodrilus guaymasensis, Peltospira smaragdina, Perotrochus caledonicus, Pseudococculinidae sp., and Shinkailepas briandi. This dataset provided the first reports of the mitogenomes for the Cocculiniformia, three vetigastropod superfamilies: Pleurotomarioidea, Lepetelloidea, and Scissurelloidea, and the neritimorph family Phenacolepadidae. The addition of deep-sea representatives also allowed us to evaluate the evolution of habitat use in gastropods. Our results showed a strongly supported sister-group relationship between the deep-sea lineages Cocculiniformia and Neomphalina. Within the Vetigastropoda, the Pleurotomarioidea was revealed as the sister-group of the remaining vetigastropods. Although this clade was presently restricted to the deep sea, fossil records showed that it has only recently invaded this habitat, thus suggesting that shallow waters was the ancestral habitat for the Vetigastropoda. The deep-sea Lepetelloidea and Lepetodriloidea formed a well-supported clade, with the Scissurelloidea sister to it, suggesting an early transition from shallow water to deep sea in this lineage. In addition, the switch between different chemosynthetic habitats was also observed in deep-sea gastropod lineages, notably in Neomphalina and Lepetelloidea. In both cases, the biogenic substrates appeared as the putative ancestral habitat, confirming the previously proposed hypothesis of a wooden-step to deep-sea vents scenario of evolution of habitat use for these taxa.

Unusual AT-skew of Sinorhodeus microlepis mitogenome provides new insights into mitogenome features and phylogenetic implications of bitterling fishes

Sinorhodeus microlepis (S. microlepis) is recently described as a new species and represents a new genus Sinorhodeu of the subfamily Acheilognathinae. In this study, we first sequenced the complete mitogenome of S. microlepis and compared with the other 29 bitterling mitogenomes. The S. microlepis mitogenome is 16,591 bp in length and contains 37 genes. Gene distribution pattern is identical among 30 bitterling mitogenomes. A significant linear correlation between A+T% and AT-skew were found among 29 bitterling mitogenomes, except S. microlepis shows unusual AT-skew with slightly negative in tRNAs and PCGs. Bitterling mitogenomes exhibit highly conserved usage bias of start codon, relative synonymous codons and amino acids, overlaps and non-coding intergenic spacers. Phylogenetic trees constructed by 13 PCGs strongly support the polyphyly of the genus Acheilognathus and the paraphyly of Rhodeus and Tanakia. Together with the unusual characters of S. microlepis mitogenomes and phylogenetic trees, S. microlepis should be a sister species to the genus Rhodeu that might diverge about 13.69 Ma (95% HPD: 12.96-14.48 Ma).

Cause or casualty: The role of mitochondrial DNA in aging and age-associated disease

The mitochondrial genome (mtDNA) represents a tiny fraction of the whole genome, comprising just 16.6 kilobases encoding 37 genes involved in oxidative phosphorylation and the mitochondrial translation machinery. Despite its small size, much interest has developed in recent years regarding the role of mtDNA as a determinant of both aging and age-associated diseases. A number of studies have presented compelling evidence for key roles of mtDNA in age-related pathology, although many are correlative rather than demonstrating cause. In this review we will evaluate the evidence supporting and opposing a role for mtDNA in age-associated functional declines and diseases. We provide an overview of mtDNA biology, damage and repair as well as the influence of mitochondrial haplogroups, epigenetics and maternal inheritance in aging and longevity.

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.

Characterization of the complete mitochondrial genome of Simulium (Byssodon) maculatum (Diptera: Simuliidae) and its phylogenetic implications

The mitochondrial (mt) genome of the black fly Simulium (Byssodon) maculatum, a pest of great importance for both humans and livestock, is sequenced and annotated for the first time. The genome structure, gene order and codon usage are typical among Diptera mt genomes. The mt genome is circular and 15,799 bp in length with 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR), and with weakly positive AT-skew (0.02) and negative GC-skew (-0.12). Phylogenetic relationships of 16 species representing five families of Culicomorpha and two outgroups, based on mt genome data, were analyzed using both Maximum Likelihood and Bayesian methods. The monophyly of Culicomorpha is well supported, while Chironomoidea is indicated as a paraphyletic group. The well supported monophyletic Simuliidae is the sister group to Culicidae.

Comparative analysis of the mitochondrial genomes of oriental spittlebug trible Cosmoscartini: insights into the relationships among closely related taxa

Background

Cosmoscartini (Hemiptera: Cercopoidea: Cercopidae) is a large and brightly colored Old World tropical tribe, currently containing over 310 phytophagous species (including some economically important pests of eucalyptus in China) in approximately 17 genera. However, very limited information of Cosmoscartini is available except for some scattered taxonomic studies. Even less is known about its phylogenetic relationship, especially among closely related genera or species. In this study, the detailed comparative genomic and phylogenetic analyses were performed on nine newly sequenced mitochondrial genomes (mitogenomes) of Cosmoscartini, with the purpose of exploring the taxonomic status of the previously defined genus Okiscarta and some closely related species within the genus Cosmoscarta.

Results

Mitogenomes of Cosmoscartini display similar genomic characters in terms of gene arrangement, nucleotide composition, codon usage and overlapping regions. However, there are also many differences in intergenic spacers, mismatches of tRNAs, and the control region. Additionally, the secondary structures of rRNAs within Cercopidae are inferred for the first time.

Based on comparative genomic (especially for the substitution pattern of tRNA secondary structure) and phylogenetic analyses, the representative species of Okiscarta uchidae possesses similar structures with other Cosmoscarta species and is placed consistently in Cosmoscarta. Although Cosmoscarta bimacula is difficult to be distinguished from Cosmoscarta bispecularis by traditional morphological methods, evidence from mitogenomes highly support the relationships of (C. bimacula + Cosmoscarta rubroscutellata) + (C. bispecularis + Cosmoscarta sp.).

Conclusions

This study presents mitogenomes of nine Cosmoscartini species and represents the first detailed comparative genomic and phylogenetic analyses within Cercopidae. It is indicated that knowledge of mitogenomes can be effectively used to resolve phylogenetic relationships at low taxonomic levels. Sequencing more mitogenomes at various taxonomic levels will also improve our understanding of mitogenomic evolution and phylogeny in Cercopidae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5365-7) contains supplementary material, which is available to authorized users.

Relation between mitochondrial DNA hyperdiversity, mutation rate and mitochondrial genome evolution in Melarhaphe neritoides (Gastropoda: Littorinidae) and other Caenogastropoda

Mitochondrial DNA hyperdiversity is primarily caused by high mutation rates (µ) and has potential implications for mitogenome architecture and evolution. In the hyperdiverse mtDNA of Melarhaphe neritoides (Gastropoda: Littorinidae), high mutational pressure generates unusually large amounts of synonymous variation, which is expected to (1) promote changes in synonymous codon usage, (2) reflect selection at synonymous sites, (3) increase mtDNA recombination and gene rearrangement, and (4) be correlated with high mtDNA substitution rates. The mitogenome of M. neritoides was sequenced, compared to closely related littorinids and put in the phylogenetic context of Caenogastropoda, to assess the influence of mtDNA hyperdiversity and high µ on gene content and gene order. Most mitogenome features are in line with the trend in Mollusca, except for the atypical secondary structure of the methionine transfer RNA lacking the TΨC-loop. Therefore, mtDNA hyperdiversity and high µ in M. neritoides do not seem to affect its mitogenome architecture. Synonymous sites are under positive selection, which adds to the growing evidence of non-neutral evolution at synonymous sites. Under such non-neutrality, substitution rate involves neutral and non-neutral substitutions, and high µ is not necessarily associated with high substitution rate, thus explaining that, unlike high µ, a high substitution rate is associated with gene order rearrangement.

Characterization of five complete Cyrtodactylus mitogenome structures reveals low structural diversity and conservation of repeated sequences in the lineage

Mitochondrial genomes (mitogenomes) of five Cyrtodactylus were determined. Their compositions and structures were similar to most of the available gecko lizard mitogenomes as 13 protein-coding, two rRNA and 22 tRNA genes. The non-coding control region (CR) of almost all Cyrtodactylus mitogenome structures contained a repeated sequence named the 75-bp box family, except for C. auribalteatus which contained the 225-bp box. Sequence similarities indicated that the 225-bp box resulted from the duplication event of 75-bp boxes, followed by homogenization and fixation in C. auribalteatus. The 75-bp box family was found in most gecko lizards with high conservation (55-75% similarities) and could form secondary structures, suggesting that this repeated sequence family played an important role under selective pressure and might involve mitogenome replication and the likelihood of rearrangements in CR. The 75-bp box family was acquired in the common ancestral genome of the gecko lizard, evolving gradually through each lineage by independent nucleotide mutation. Comparison of gecko lizard mitogenomes revealed low structural diversity with at least six types of mitochondrial gene rearrangements. Cyrtodactylus mitogenome structure showed the same gene rearrangement as found in most gecko lizards. Advanced mitogenome information will enable a better understanding of structure evolution mechanisms.

Basal position of two new complete mitochondrial genomes of parasitic Cymothoida (Crustacea: Isopoda) challenges the monophyly of the suborder and phylogeny of the entire order

BACKGROUND

Isopoda is a highly diverse order of crustaceans with more than 10,300 species, many of which are parasitic. Taxonomy and phylogeny within the order, especially those of the suborder Cymothoida Wägele, 1989, are still debated. Mitochondrial (mt) genomes are a useful tool for phylogenetic studies, but their availability for isopods is very limited. To explore these phylogenetic controversies on the mt genomic level and study the mt genome evolution in Isopoda, we sequenced mt genomes of two parasitic isopods, Tachaea chinensis Thielemann, 1910 and Ichthyoxenos japonensis Richardson, 1913, belonging to the suborder Cymothoida, and conducted comparative and phylogenetic mt genomic analyses across Isopoda.

RESULTS

The complete mt genomes of T. chinensis and I. japonensis were 14,616 bp and 15,440 bp in size, respectively, with the A+T content higher than in other isopods (72.7 and 72.8%, respectively). Both genomes code for 13 protein-coding genes, 21 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and possess a control region (CR). Both are missing a gene from the complete tRNA set: T. chinensis lacks trnS1 and I. japonensis lacks trnI. Both possess unique gene orders among isopods. Within the CR of I. japonensis (284 bp), we identified a repetitive region with four tandem repeats. Phylogenetic analysis based on concatenated nucleotide sequences of 13 protein-coding genes showed that the two parasitic cymothoids clustered together and formed a basal clade within Isopoda. However, another parasitic cymothoid, Gyge ovalis Shiino, 1939, formed a sister group with the suborder Limnoriidea Brandt & Poore in Poore, 2002, whereas two free-living cymothoid species were located in the derived part of the phylogram: Bathynomus sp. formed a sister group with the suborder Sphaeromatidea Wägele, 1989, and Eurydice pulchra Leach, 1815 with a clade including Bathynomus sp., Sphaeromatidea and Valvifera G. O. Sars, 1883.

CONCLUSIONS

Our results did not recover the suborders Cymothoida and Oniscidea Latreille, 1802 as monophyletic, with parasitic and free-living cymothoidans forming separate clades. Furthermore, two parasitic cymothoidans formed the sister-clade to all other isopods, separated from Epicaridea Latreille, 1825, which challenges currently prevalent isopod phylogeny. Additional mt genomes of parasitic and free-living isopods might confer a sufficient phylogenetic resolution to enable us to resolve their relationships, and ultimately allow us to better understand the evolutionary history of the entire isopod order.

The first mitochondrial genome from Scopuridae (Insecta: Plecoptera) reveals structural features and phylogenetic implications

This study determined the first complete mitochondrial genome (mitogenome) of a stonefly, Scopura longa (Plecoptera: Scopuridae), and reconstructed a phylogeny based on two datasets of mitogenomes in eighteen available stoneflies to examine the relationships among Plecoptera. The complete mitogenome of S. longa is a circular molecule of 15,798bp in size. It contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). Most PCGs used standard ATN start codons and TAN stop codons. All tRNAs could be folded as typical cloverleaf secondary structures except tRNA^Ser(AGN), with the dihydrouridine (DHU) arm missing. Tandem repeat regions, two potential stem-loop (SL) structures, Poly T structures and G+C-rich regions are detected in the control region. Finally, the phylogenetic relationships among the families within the Arctoperlaria were reconstructed. The topological structures of the two trees were almost identical. The present phylogenetic analysis shows that S. longa belongs to the infraorder Euholognatha. The monophyly of each family is generally well supported based on nucleotide sequences. The Pteronarcyidae is sister to ((Peltoperlidae+Styloperlidae)+(Perlidae+(Perlodidae+Chloroperlidae))).

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies

BACKGROUND

As the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class.

RESULTS

Taxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses.

CONCLUSIONS

Our results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

Genetic homogeneity in South American tomato pinworm, Tuta absoluta: a new invasive pest to oriental region

South American tomato leaf miner, Tuta absoluta (Meyrick, 1917) (Lepidoptera: Gelechiidae), is an important invasive pest of Tomato which invaded India and Nepal in 2014 and 2016, respectively. In the present study, samples from five localities of India and one from Nepal were used for the investigation of genetic diversity of T. absoluta by employing a fragment in the mtDNA gene-encoding cytochrome oxidase I (COI). Based on the partial COI gene, high genetic homogeneity was detected in T. absoluta populations of India and Nepal with rest of the world. Less nucleotide diversity (π 0.00137) was also detected in the populations of T. absoluta from different countries. This is first attempt to analyze molecular data for this new invasive species from India and Nepal.

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

Information from the mitochondrial genomes of two egg parasitoids, Gonatocerus sp. and Telenomus sp., reveals a controversial phylogenetic relationship between Mymaridae and Scelionidae

The taxonomic status and phylogenetic affinities of Mymaridae and Scelionidae are controversial, based on similarities between these families in the characteristics of adults, larvae, and eggs. In this study, we sequenced the mitochondrial (mt) genomes of representatives from these two families and found that the derived secondary structure of tRNA-Arg was the same in each family due to the absence of the D-stem. The segment of "cox1 trnL₂cox2 trnK trnD atp8 atp6 cox3" in Gonatocerus sp. (Mymaridae) is conserved and distinct from those of four other species of Chalcidoidea but similar to that in Proctotrupoidea and Platygastroidea. However, phylogenetic analysis indicated that Gonatocerus sp. was sister group to other species of Chalcidoidea. Comparisons based on complete gene orders may be more useful in a phylogenetic and systematic context, as different branches may exhibit partially homoplastic gene orders.

Actively transcribed and expressed atp8 gene in Mytilus edulis mussels

Background

Animal mitochondrial genomes typically encode 37 genes: 13 proteins, 22 tRNAs and two rRNAs. However, many species represent exceptions to that rule. Bivalvia along with Nematoda and Platyhelminthes are often suspected to fully or partially lack the ATP synthase subunit 8 (atp8) gene. This raises the question as to whether they are really lacking this gene or is this maybe an annotation problem? Among bivalves, Mytilus edulis has been inferred to lack an ATP8 gene since the characterization of its mitochondrial genome in 1992. Even though recent bioinformatic analyses suggested that atp8 is present in Mytilus spp., due to high divergence in predicted amino acid sequences, the existence of a functional atp8 gene in this group remains controversial.

Results

Here we demonstrate that M. edulis mitochondrial open reading frames suggested to be atp8 (in male and female mtDNAs) are actively translated proteins. We also provide evidence that both proteins are an integral part of the ATP synthase complex based on in-gel detection of ATP synthase activity and two-dimensional Blue-Native and SDS polyacrylamide electrophoresis.

Conclusion

Many organisms (e.g., Bivalvia along with Nematoda and Platyhelminthes) are considered to be lacking certain mitochondrial genes often only based on poor similarity between protein coding gene sequences in genetically closed species. In some situations, this may lead to the inference that the ATP8 gene is absent, when it is in fact present, but highly divergent. This shows how important complementary role protein-based approaches, such as those in the present study, can provide to bioinformatic, genomic studies (i.e., ability to confirm the presence of a gene).

Sequencing of complete mitochondrial genomes confirms synonymization of Hyalomma asiaticum asiaticum and kozlovi, and advances phylogenetic hypotheses for the Ixodidae

Phylogeny of hard ticks (Ixodidae) remains unresolved. Mitochondrial genomes (mitogenomes) are increasingly used to resolve phylogenetic controversies, but remain unavailable for the entire large Hyalomma genus. Hyalomma asiaticum is a parasitic tick distributed throughout the Asia. As a result of great morphological variability, two subspecies have been recognised historically; until a morphological data-based synonymization was proposed. However, this hypothesis was never tested using molecular data. Therefore, objectives of this study were to: 1. sequence the first Hyalomma mitogenome; 2. scrutinise the proposed synonymization using molecular data, i.e. complete mitogenomes of both subspecies: H. a. asiaticum and kozlovi; 3. conduct phylogenomic and comparative analyses of all available Ixodidae mitogenomes. Results corroborate the proposed synonymization: the two mitogenomes are almost identical (99.6%). Genomic features of both mitogenomes are standard for Metastriata; which includes the presence of two control regions and all three "Tick-Box" motifs. Gene order and strand distribution are perfectly conserved for the entire Metastriata group. Suspecting compositional biases, we conducted phylogenetic analyses (29 almost complete mitogenomes) using homogeneous and heterogeneous (CAT) models of substitution. The results were congruent, apart from the deep-level topology of prostriate ticks (Ixodes): the homogeneous model produced a monophyletic Ixodes, but the CAT model produced a paraphyletic Ixodes (and thereby Prostriata), divided into Australasian and non-Australasian clades. This topology implies that all metastriate ticks have evolved from the ancestor of the non-Australian branch of prostriate ticks. Metastriata was divided into three clades: 1. Amblyomminae and Rhipicephalinae (Rhipicephalus, Hyalomma, Dermacentor); 2. Haemaphysalinae and Bothriocrotoninae, plus Amblyomma sphenodonti; 3. Amblyomma elaphense, basal to all Metastriata. We conclude that mitogenomes have the potential to resolve the long-standing debate about the evolutionary history of ticks, but heterogeneous evolutionary models should be used to alleviate the effects of compositional heterogeneity on deep-level relationships.

Comparative analyses of the complete mitochondrial genomes of Dosinia clams and their phylogenetic position within Veneridae

Mitochondrial genomes have proved to be a powerful tool in resolving phylogenetic relationship. In order to understand the mitogenome characteristics and phylogenetic position of the genus Dosinia, we sequenced the complete mitochondrial genomes of Dosinia altior and Dosinia troscheli (Bivalvia: Veneridae), compared them with that of Dosinia japonica and established a phylogenetic tree for Veneridae. The mitogenomes of D. altior (17,536 bp) and D. troscheli (17,229 bp) are the two smallest in Veneridae, which include 13 protein-coding genes, 2 ribosomal RNA genes, 22 tRNA genes, and non-coding regions. The mitogenomes of the Dosinia species are similar in size, gene content, AT content, AT- and GC- skews, and gene arrangement. The phylogenetic relationships of family Veneridae were established based on 12 concatenated protein-coding genes using maximum likelihood and Bayesian analyses, which supported that Dosininae and Meretricinae have a closer relationship, with Tapetinae being the sister taxon. The information obtained in this study will contribute to further understanding of the molecular features of bivalve mitogenomes and the evolutionary history of the genus Dosinia.

Extensive mitochondrial gene rearrangements in Ctenophora: insights from benthic Platyctenida

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Phylogeny, evolution and mitochondrial gene order rearrangement in scale worms (Aphroditiformia, Annelida)

Next-generation sequencing (NGS) has become a powerful tool in phylogenetic and evolutionary studies. Here we applied NGS to recover two ribosomal RNA genes (18S and 28S) from 16 species and 15 mitochondrial genomes from 16 species of scale worms representing six families in the suborder Aphroditiformia (Phyllodocida, Annelida), a complex group of polychaetes characterized by the presence of dorsal elytra or scales. The phylogenetic relationship of the several groups of scale worms remains unresolved due to insufficient taxon sampling and low resolution of individual gene markers. Phylogenetic tree topology based on mitochondrial genomes is comparable with that based on concatenated sequences from two mitochondrial genes (cox1 and 16S) and two ribosomal genes (18S and 28S) genes, but has higher statistical support for several clades. Our analyses show that Aphroditiformia is monophyletic, indicating the presence of elytra is an apomorphic trait. Eulepethidae and Aphroditidae together form the sister group to all other families in this suborder, whereas Acoetidae is sister to Iphionidae. Polynoidae is monophyletic, but within this family the deep-sea subfamilies Branchinotogluminae and Macellicephalinae are paraphyletic. Mitochondrial genomes in most scale-worm families have a conserved gene order, but within Polynoidae there are two novel arrangement patterns in the deep-sea clade. Mitochondrial protein-coding genes in polynoids as a whole have evolved under strong purifying selection, but substitution rates in deep-sea species are much higher than those in shallow-water species, indicating that purifying selection is relaxed in deep-sea polynoids. There are positive selected amino acids for some mitochondrial genes of the deep-sea clade, indicating they may involve in the adaption of deep-sea polynoids. Overall, our study (1) provided more evidence for reconstruction of the phylogeny of Aphroditiformia, (2) provided evidence to refute the assumption that mitochondrial gene order in Errantia is conserved, and (3) indicated that the deep-sea extreme environment may have affected the mitochondrial genome evolution rate and gene order arrangement in Polynoidae.

A complete logical approach to resolve the evolution and dynamics of mitochondrial genome in bilaterians

Investigating how recombination might modify gene order during the evolution has become a routine part of mitochondrial genome analysis. A new method of genomic maps analysis based on formal logic is described. The purpose of this method is to 1) use mitochondrial gene order of current taxa as datasets 2) calculate rearrangements between all mitochondrial gene orders and 3) reconstruct phylogenetic relationships according to these calculated rearrangements within a tree under the assumption of maximum parsimony. Unlike existing methods mainly based on the probabilistic approach, the main strength of this new approach is that it calculates all the exact tree solutions with completeness and provides logical consequences as highly robust results. Moreover, this method infers all possible hypothetical ancestors and reconstructs character states for all internal nodes of the trees. We started by testing our method using the deuterostomes as a study case. Then, with sponges as an outgroup, we investigated the evolutionary history of mitochondrial genomes of 47 bilaterian phyla and emphasised the peculiar case of chaetognaths. This pilot work showed that the use of formal logic in a hypothetico-deductive background such as phylogeny (where experimental testing of hypotheses is impossible) is very promising to explore mitochondrial gene order in deuterostomes and should be applied to many other bilaterian clades.

The rearranged mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae), a parasitoid wasp of mantis

The complete sequence of the mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae) is described. The mitogenome was 15,845 bp in size, and contained typical sets of mitochondrial genes. The base composition of the Podagrion sp. mitogenome was also biased toward A + T bases (81.8%). The mitochondrial genome of Podagrion sp. has a weak AT skew (0.07) and a strong GC skew (-0.26). Podagrion sp. exhibits a novel rearrangement compared with the ancestral order, including six protein-coding genes (nad3, cox3, atp6, atp8, cox2 and cox1), which have inverted to the minor strand from the major strand. The A + T-rich region of Podagrion sp., which is located between trnN and trnI, have five tandem repeats. The apomorphic rearrangements, including the conserved block "cox3-atp6-atp8-cox2-cox1-nad5-nad4-nad4l-nad6-cob" and the special locations of trnV and trnA, were mapped onto the phylogeny of Proctotrupomorpha.

The Complete Mitochondrial Genome of the Plant Bug Lygus pratensis Linnaeus (Hemiptera: Miridae)

Lygus pratensis is a phytophagous pest responsible for yield losses in Bt alfalfa and other economic crops in Northwestern China. To better characterize Miridae at the genomic level, the complete mitochondrial (mt) genome of L. pratensis was sequenced and analyzed in this study. The mt genome was amplified via the polymerase chain reaction to generate overlapping fragments. These fragments were then sequenced, spliced, and analyzed to include the examination of nucleotide composition, codon usage, compositional biases, protein-coding genes (PCGs), and RNA secondary structures. Phylogenetic relationships between L. pratensis and other species in different Heteroptera families were also examined. The mt genome was found to be a typical circular genome with a length of 16,591 bp and a total AT content of 75.1%, encoded for 13 PCGs, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (lrRNA and srRNA), and a noncoding control region. The nucleotide composition of the entire mt genome was heavily biased toward A and T. All of the tRNAs were predicted to have classic clover leaf structures, but three of the tRNAs (tRNAAsn, tRNAHis, tRNAHis) were missing the TΨC loop. The control region (2,017 bp), which was found to be located between 12S and tRNAIle, contained three tandem repeat elements. Phylogenetic analyses showed that L. pratensis is closely related to the other three examined Lygus bugs, and that it is a sister group to Apolygus and Adelphocoris. This study confirms the usability of the mt genome in phylogenesis studies pertaining to the Lygus genus, within Miridae.

The Complete Mitochondrial Genome of Coptotermes 'suzhouensis' (syn. Coptotermes formosanus) (Isoptera: Rhinotermitidae) and Molecular Phylogeny Analysis

Coptotermes suzhouensis (Isoptera: Rhinotermitidae) is a significant subterranean termite pest of wooden structures and is widely distributed in southeastern China. The complete mitochondrial DNA sequence of C. suzhouensis was analyzed in this study. The mitogenome was a circular molecule of 15,764 bp in length, which contained 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A+T-rich region with a gene arrangement typical of Isoptera mitogenomes. All PCGs were initiated by ATN codons and terminated by complete termination codons (TAA), except COX2, ND5, and Cytb, which ended with an incomplete termination codon T. All tRNAs displayed a typical clover-leaf structure, except for tRNASer(AGN), which did not contain the stem-loop structure in the DHU arm. The A+T content (69.23%) of the A+T-rich region (949 bp) was higher than that of the entire mitogenome (65.60%), and two different sets of repeat units (A+B) were distributed in this region. Comparison of complete mitogenome sequences with those of Coptotermes formosanus indicated that the two taxa have very high genetic similarity. Forty-one representative termite species were used to construct phylogenetic trees by maximum likelihood, maximum parsimony, and Bayesian inference methods. The phylogenetic analyses also strongly supported (BPP, MLBP, and MPBP = 100%) that all C. suzhouensis and C. formosanus samples gathered into one clade with genetic distances between 0.000 and 0.002. This study provides molecular evidence for a more robust phylogenetic position of C. suzhouensis and inferrs that C. suzhouensis was the synonymy of C. formosanus.

Complete Mitochondrial Genome of Suwallia teleckojensis (Plecoptera: Chloroperlidae) and Implications for the Higher Phylogeny of Stoneflies

Stoneflies comprise an ancient group of insects, but the phylogenetic position of Plecoptera and phylogenetic relations within Plecoptera have long been controversial, and more molecular data is required to reconstruct precise phylogeny. Herein, we present the complete mitogenome of a stonefly, Suwallia teleckojensis, which is 16146 bp in length and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and a control region (CR). Most PCGs initiate with the standard start codon ATN. However, ND5 and ND1 started with GTG and TTG. Typical termination codons TAA and TAG were found in eleven PCGs, and the remaining two PCGs (COII and ND5) have incomplete termination codons. All transfer RNA genes (tRNAs) have the classic cloverleaf secondary structures, with the exception of tRNASer(AGN), which lacks the dihydrouridine (DHU) arm. Secondary structures of the two ribosomal RNAs were shown referring to previous models. A large tandem repeat region, two potential stem-loop (SL) structures, Poly N structure (2 poly-A, 1 poly-T and 1 poly-C), and four conserved sequence blocks (CSBs) were detected in the control region. Finally, both maximum likelihood (ML) and Bayesian inference (BI) analyses suggested that the Capniidae was monophyletic, and the other five stonefly families form a monophyletic group. In this study, S. teleckojensis was closely related to Sweltsa longistyla, and Chloroperlidae and Perlidae were herein supported to be a sister group.