Large gene overlaps and tRNA processing in the compact mitochondrial genome of the crustacean Armadillidium vulgare

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Novel Gene Rearrangement Pattern in Pachycrepoideus vindemmiae Mitochondrial Genome: New Gene Order in Pteromalidae (Hymenoptera: Chalcidoidea)

The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Mitochondrial phylogenomics provides conclusive evidence that the family Ancyrocephalidae is deeply paraphyletic

BACKGROUND

Unresolved taxonomic classification and paraphyly pervade the flatworm class Monogenea: the class itself may be paraphyletic and split into Polyopisthocotylea and Monopisthocotylea; there are some indications that the monopisthocotylean order Dactylogyridea may also be paraphyletic; single-gene markers and some morphological traits indicate that the family Ancyrocephalidae is paraphyletic and intertwined with the family Dactylogyridae.

METHODS

To attempt to study the relationships of Ancyrocephalidae and Monopisthocotylea using a phylogenetic marker with high resolution, we sequenced mitochondrial genomes of two fish ectoparasites from the family Dactylogyridae: Dactylogyrus simplex and Dactylogyrus tuba. We conducted phylogenetic analyses using three datasets and three methods. Datasets were ITS1 (nuclear) and nucleotide and amino acid sequences of almost complete mitogenomes of almost all available Monopisthocotylea mitogenomes. Methods were maximum likelihood (IQ-TREE), Bayesian inference (MrBayes) and CAT-GTR (PhyloBayes).

RESULTS

Both mitogenomes exhibited the ancestral gene order for Neodermata, and both were compact, with few and small intergenic regions and many and large overlaps. Gene sequences were remarkably divergent for nominally congeneric species, with only trnI exhibiting an identity value > 80%. Both mitogenomes had exceptionally low A + T base content and AT skews. We found evidence of pervasive compositional heterogeneity in the dataset and indications that base composition biases cause phylogenetic artefacts. All six mitogenomic analyses produced unique topologies, but all nine analyses produced topologies that rendered Ancyrocephalidae deeply paraphyletic. Mitogenomic data consistently resolved the order Capsalidea as nested within the Dactylogyridea.

CONCLUSIONS

The analyses indicate that taxonomic revisions are needed for multiple Polyopisthocotylea lineages, from genera to orders. In combination with previous findings, these results offer conclusive evidence that Ancyrocephalidae is a paraphyletic taxon. The most parsimonious solution to resolve this is to create a catch-all Dactylogyridae sensu lato clade comprising the current Ancyrocephalidae, Ancylodiscoididae, Pseudodactylogyridae and Dactylogyridae families, but the revision needs to be confirmed by another marker with a sufficient resolution.

Enhanced Resolution of Evolution and Phylogeny of the Moths Inferred from Nineteen Mitochondrial Genomes

The vast majority (approximately 90%) of Lepidoptera species belong to moths whose phylogeny has been widely discussed and highly controversial. For the further understanding of phylogenetic relationships of moths, nineteen nearly complete mitochondrial genomes (mitogenomes) of moths involved in six major lineages were sequenced and characterized. These mitogenomes ranged from 15,177 bp (Cyclidia fractifasciata) to 15,749 bp (Ophthalmitis albosignaria) in length, comprising of the core 37 mitochondrial genes (13 protein-coding genes (PCGs) + 22 tRNAs + two rRNAs) and an incomplete control region. The order and orientation of genes showed the same pattern and the gene order of trnM-trnI-trnQ showed a typical rearrangement of Lepidoptera compared with the ancestral order of trnI-trnQ-trnM. Among these 13 PCGs, ATP8 exhibited the fastest evolutionary rate, and Drepanidae showed the highest average evolutionary rate among six families involved in 66 species. The phylogenetic analyses based on the dataset of 13 PCGs suggested the relationship of (Notodontidae + (Noctuidae + Erebidae)) + (Geometridae + (Sphingidae + Drepanidae)), which suggested a slightly different pattern from previous studies. Most groups were well defined in the subfamily level except Erebidae, which was not fully consistent across bayesian and maximum likelihood methods. Several formerly unassigned tribes of Geometridae were suggested based on mitogenome sequences despite a not very strong support in partial nodes. The study of mitogenomes of these moths can provide fundamental information of mitogenome architecture, and the phylogenetic position of moths, and contributes to further phylogeographical studies and the biological control of pests.

The complete mitochondrial genomes of Paradiplozoon yarkandense and Paradiplozoon homoion confirm that Diplozoidae evolve at an elevated rate

BACKGROUND

Diplozoidae are monogenean (Monogenea: Polyopisthocotylea) fish parasites characterised by a unique life history: two larvae permanently fuse into an X-shaped "Siamese" organism. Taxonomy and phylogeny of Diplozoidae and Polyopisthocotylea remain unresolved due to the unavailability of molecular markers with sufficiently high resolution. Mitogenomes may be a suitable candidate, but there are currently only 12 available for the Polyopisthocotylea (three for Diplozoidae). The only available study of diplozoid mitogenomes found unique base composition patterns and elevated evolution rates in comparison with other Monogenean mitogenomes.

METHODS

To further explore their evolution and generate molecular data for evolutionary studies, we sequenced the complete mitogenomes of two Diplozoidae species, Paradiplozoon homoion and Paradiplozoon yarkandense, and conducted a number of comparative mitogenomic analyses with other polyopisthocotyleans.

RESULTS

We found further evidence that mitogenomes of Diplozoidae evolve at a unique, elevated rate, which was reflected in their exceptionally long branches, large sizes, unique base composition, skews, and very low gene sequence similarity levels between the two newly sequenced species. They also exhibited remarkably large overlaps between some genes. Phylogenetic analysis of Polyopisthocotylea resolved all major taxa as monophyletic, and Mazocraeidea was split into two major clades: (Diplozoidae) + (all four remaining families: Diclidophoridae, Chauhaneidae, Mazocraeidae and Microcotylidae). It also provided further confirmation that the genus Paradiplozoon is paraphyletic and requires a taxonomic revision, so the two species may have to be renamed Indodiplozoon homoion and Diplozoon yarkandense comb. nov.

CONCLUSIONS

Although our findings indicate that mitogenomes may be a promising tool for resolving the phylogeny of Polyopisthocotylea, elevated evolutionary rates of Diplozoidae may cause phylogenetic artefacts, so future studies should pay caution to this problem. Furthermore, as the reason for their elevated evolution remains unknown, Diplozoidae are a remarkably interesting lineage for other types of evolutionary mitogenomic studies.

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

Structural Features and Phylogenetic Implications of 11 New Mitogenomes of Typhlocybinae (Hemiptera: Cicadellidae)

To explore the characteristics of mitogenomes and discuss the phylogenetic relationships and molecular evolution of the six tribes within Typhlocybinae, 11 complete mitogenomes are newly sequenced and comparatively analyzed. In all of these complete mitogenomes, the number and order of the genes are highly conserved in overall organization. The PCGs initiate with ATN/TTG/GTG and terminate with TAA/TAG/T. Almost all tRNAs are folded into the typical clover-leaf secondary structure. The control region is always variable in length and in numbers of multiple tandem repeat units. The atp8 and nad2 exhibits the highest evolution rate among all the PCGs. Phylogenetic analyses based on whole mitogenome sequences, with three different datasets, using both maximum likelihood and Bayesian methods, indicate the monophyly of Typhlocybinae and its inner tribes, respectively, except for Typhlocybini and Zyginellini that are paraphyletic. Finally, we confirm that Erythroneurini is a subtribe of Dikraneurini.

An Update on Trichoderma Mitogenomes: Complete De Novo Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma harzianum (Hypocreales, Sordariomycetes), an Ex-Neotype Strain CBS 226.95, and Tracing the Evolutionary Divergences of Mitogenomes in Trichoderma

Members of the genus Trichoderma (Hypocreales), widely used as biofungicides, biofertilizers, and as model fungi for the industrial production of CAZymes, have actively been studied for the applications of their biological functions. Recently, the study of the nuclear genomes of Trichoderma has expanded in the directions of adaptation and evolution to gain a better understanding of their ecological traits. However, Trichoderma's mitochondria have received much less attention despite mitochondria being the most necessary element for sustaining cell life. In this study, a mitogenome of the fungus Trichoderma harzianum CBS 226.95 was assembled de novo. A 27,632 bp circular DNA molecule was revealed with specific features, such as the intronless of all core PCGs, one homing endonuclease, and a putative overlapping tRNA, on a closer phylogenetic relationship with T. reesei among hypocrealean fungi. Interestingly, the mitogenome of T. harzianum CBS 226.95 was predicted to have evolved earlier than those of other Trichoderma species and also assumed with a selection pressure in the cox3. Considering the bioavailability, both for the ex-neotype strain of the T. harzianum species complex and the most globally representative commercial fungal biocontrol agent, our results on the T. harzianum CBS 226.95 mitogenome provide crucial information which will be helpful criteria in future studies on Trichoderma.

Two Complete Mitochondrial Genomes of Mileewinae (Hemiptera: Cicadellidae) and a Phylogenetic Analysis

Simple Summary

Mileewinae is a small subfamily of Cicadellidae containing about 160 described species, extensively distributed in the Oriental, Ethiopian and Neotropical regions. Some species are potential pests in agriculture and forestry. The classification of this group has been unstable over the past few decades. Currently, some controversies remain on the monophyly of Mileewinae and phylogenetic relationships of Mileewinae with other subfamilies. To provide further evidence toward answering these questions, two newly completed mitochondrial genomes of Mileewinae species (Mileewa rufivena and Ujna puerana) have been sequenced and analyzed. Results show these two mitochondrial genomes have quite similar structures and features. In phylogenetic analyses, Mileewinae formed a monophyletic group in Cicadellidae in all trees derived from maximum likelihood (ML) and Bayesian inference (BI) methods. In addition, Mileewinae has a closer phylogenetic relationship with Typhlocybinae compared to the Cicadellinae.

Abstract

More studies are using mitochondrial genomes of insects to explore the sequence variability, evolutionary traits, monophyly of groups and phylogenetic relationships. Controversies remain on the classification of the Mileewinae and the phylogenetic relationships between Mileewinae and other subfamilies remain ambiguous. In this study, we present two newly completed mitogenomes of Mileewinae (Mileewa rufivena Cai and Kuoh 1997 and Ujna puerana Yang and Meng 2010) and conduct comparative mitogenomic analyses based on several different factors. These species have quite similar features, including their nucleotide content, codon usage of protein genes and the secondary structure of tRNA. Gene arrangement is identical and conserved, the same as the putative ancestral pattern of insects. All protein-coding genes of U. puerana began with the start codon ATN, while 5 Mileewa species had the abnormal initiation codon TTG in ND5 and ATP8. Moreover, M. rufivena had an intergenic spacer of 17 bp that could not be found in other mileewine species. Phylogenetic analysis based on three datasets (PCG123, PCG12 and AA) with two methods (maximum likelihood and Bayesian inference) recovered the Mileewinae as a monophyletic group with strong support values. All results in our study indicate that Mileewinae has a closer phylogenetic relationship to Typhlocybinae compared to Cicadellinae. Additionally, six species within Mileewini revealed the relationship (U. puerana + (M. ponta + (M. rufivena + M. alara) + (M. albovittata + M. margheritae))) in most of our phylogenetic trees. These results contribute to the study of the taxonomic status and phylogenetic relationships of Mileewinae.

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.

Bioinformatic Prediction of an tRNASec Gene Nested inside an Elongation Factor SelB Gene in Alphaproteobacteria

In bacteria, selenocysteine (Sec) is incorporated into proteins via the recoding of a particular codon, the UGA stop codon in most cases. Sec-tRNASec is delivered to the ribosome by the Sec-dedicated elongation factor SelB that also recognizes a Sec-insertion sequence element following the codon on the mRNA. Since the excess of SelB may lead to sequestration of Sec-tRNASec under selenium deficiency or oxidative stress, the expression levels of SelB and tRNASec should be regulated. In this bioinformatic study, I analyzed the Rhizobiales SelB species because they were annotated to have a non-canonical C-terminal extension. I found that the open reading frame (ORF) of diverse Alphaproteobacteria selB genes includes an entire tRNASec sequence (selC) and overlaps with the start codon of the downstream ORF. A remnant tRNASec sequence was found in the Sinorhizobium melilotiselB genes whose products have a shorter C-terminal extension. Similar overlapping traits were found in Gammaproteobacteria and Nitrospirae. I hypothesized that once the tRNASec moiety is folded and processed, the expression of the full-length SelB may be repressed. This is the first report on a nested tRNA gene inside a protein ORF in bacteria.

Characterization of the Complete Mitochondrial Genomes of Two Species with Preliminary Investigation on Phylogenetic Status of Zyginellini (Hemiptera: Cicadellidae: Typhlocybinae)

To explore the characteristics of mitogenomes and reveal phylogenetic relationships of the tribes of Zyginellini and Typhlocybini in Typhlocybinae, mitogenomes of two species of the Zyginellini, Parazyginella tiani and Limassolla sp., were sequenced. Mitogenomes of both species contain 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a large non-coding region (A + T-rich region). These characteristics are similar to other Membracoidea mitogenomes. All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA/G or with an incomplete T codon. All tRNAs have the typical clover-leaf structure, except trnS1 which has a reduced DHU arm and the acceptor stem of trnR is 5 or 6 bp in some species, an unusual feature here reported for the first time in Typhlocybinae. The A + T-rich region is highly variable in length and in numbers of tandem repeats present. Our analyses indicate that nad6 and atp6 exhibit higher evolutionary rates compared to other PCGs. Phylogenetic analyses by both maximum likelihood and Bayesian methods based on 13 protein-coding genes of 12 species of Typhlocybinae suggest that Zyginellini are paraphyletic with respect to Typhlocybini.

Complete Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma atroviride: Genomic Features, Comparative Analysis and Insight Into the Mitochondrial Evolution in Trichoderma

The improvement of biopesticides for use in the agriculture industry requires an understanding of the biological- and ecological principles underlying their behavior in natural environments. The nuclear genomes of members of the genus Trichoderma, which are representative fungal biocontrol agents, have been actively studied in relation to the unique characteristics of these species as effective producers of CAZymes/secondary metabolites and biopesticides, but their mitochondrial genomes have received much less attention. In this study, the mitochondrial genome of Trichoderma atroviride (Hypocreales, Sordariomycetes), which targets wood-decaying fungal pathogens and has the ability to degrade chemical fungicides, was assembled de novo. A 32,758 bp circular DNA molecule was revealed with specific features, such as a few more protein CDS and trn genes, two homing endonucleases (LAGLIDADG-/GIY-YIG-type), and even a putative overlapping tRNA gene, on a closer phylogenetic relationship with T. gamsii among hypocrealean fungi. Particularly, introns were observed with several footprints likely to be evolutionarily associated with the intron dynamics of the Trichoderma mitochondrial genomes. This study is the first to report the complete de novo mitochondrial genome of T. atroviride, while comparative analyses of Trichoderma mitochondrial genomes were also conducted from the perspective of mitochondrial evolution for the first time.

Architectural instability, inverted skews and mitochondrial phylogenomics of Isopoda: outgroup choice affects the long-branch attraction artefacts

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

Coding palindromes in mitochondrial genes of Nematomorpha

Inverted repeats are common DNA elements, but they rarely overlap with protein-coding sequences due to the ensuing conflict with the structure and function of the encoded protein. We discovered numerous perfect inverted repeats of considerable length (up to 284 bp) embedded within the protein-coding genes in mitochondrial genomes of four Nematomorpha species. Strikingly, both arms of the inverted repeats encode conserved regions of the amino acid sequence. We confirmed enzymatic activity of the respiratory complex I encoded by inverted repeat-containing genes. The nucleotide composition of inverted repeats suggests strong selection at the amino acid level in these regions. We conclude that the inverted repeat-containing genes are transcribed and translated into functional proteins. The survey of available mitochondrial genomes reveals that several other organisms possess similar albeit shorter embedded repeats. Mitochondrial genomes of Nematomorpha demonstrate an extraordinary evolutionary compromise where protein function and stringent secondary structure elements within the coding regions are preserved simultaneously.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Arm-less mitochondrial tRNAs conserved for over 30 millions of years in spiders

BACKGROUND

In recent years, Next Generation Sequencing (NGS) has accelerated the generation of full mitogenomes, providing abundant material for studying different aspects of molecular evolution. Some mitogenomes have been observed to harbor atypical sequences with bizarre secondary structures, which origins and significance could only be fully understood in an evolutionary framework.

RESULTS

Here we report and analyze the mitochondrial sequences and gene arrangements of six closely related spiders in the sister genera Parachtes and Harpactocrates, which belong to the nocturnal, ground dwelling family Dysderidae. Species of both genera have compacted mitogenomes with many overlapping genes and strikingly reduced tRNAs that are among the shortest described within metazoans. Thanks to the conservation of the gene order and the nucleotide identity across close relatives, we were able to predict the secondary structures even on arm-less tRNAs, which would be otherwise unattainable for a single species. They exhibit aberrant secondary structures with the lack of either DHU or TΨC arms and many miss-pairings in the acceptor arm but this degeneracy trend goes even further since at least four tRNAs are arm-less in the six spider species studied.

CONCLUSIONS

The conservation of at least four arm-less tRNA genes in two sister spider genera for about 30 myr suggest that these genes are still encoding fully functional tRNAs though they may be post-transcriptionally edited to be fully functional as previously described in other species. We suggest that the presence of overlapping and truncated tRNA genes may be related and explains why spider mitogenomes are smaller than those of other invertebrates.

The mitochondrial genomes of three skippers: Insights into the evolution of the family Hesperiidae (Lepidoptera)

We sequenced the mitogenomes of Astictopterus jama, Isoteinon lamprospilus and Notocrypta curvifascia to obtain further insight into the mitogenomic architecture evolution and performed phylogenetic reconstruction using 29 Hesperiidae mitogenome sequences. The complete mitogenome sequences of A. jama, I. lamprospilus and N. curvifascia are 15,430, 15,430 and 15,546 bp in size, respectively. All contain 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and an A + T-rich region. Nucleotide composition is A + T biased, and the majority of the protein-coding genes exhibit a negative AT-skew, which is reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region is comprised of nonrepetitive sequences, including the motif ATAGA followed by a poly-T stretch, a microsatellite-like element next to the ATTTA motif, and a poly-A adjacent to tRNAs. Although most genes evolve under a strong purifying selection, the entire nad gene family (especially nad6) exhibits somewhat relaxed purifying selection, and atp8, evolving under a highly relaxed selection, is an outlier in the family Hesperiidae. Several different approaches relatively consistently indicated that nad6, atp8 and nad4 are comparatively fast-evolving genes in this family, which may have implications for future phylogenetic, population genetics and species diagnostics studies. For phylogenetic analyses of Hesperiidae, we tested a few datasets, and found that the one comprising all 37 genes produced the highest node support, indicating that the inclusion of RNAs improves the phylogenetic signal. Results indicate that subfamilies Euschemoninae, Heteropterinae, and Coeliadinae are monophyletic with strong nodal support, but Pyrginae and Eudaminae are paraphyletic. Finally, we confirm that A. jama and I. lamprospilus are close relatives.

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.

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.

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.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.

Polycytidylation of mitochondrial mRNAs in Chlamydomonas reinhardtii

The unicellular photosynthetic organism, Chlamydomonas reinhardtii, represents a powerful model to study mitochondrial gene expression. Here, we show that the 5′- and 3′-extremities of the eight Chlamydomonas mitochondrial mRNAs present two unusual characteristics. First, all mRNAs start primarily at the AUG initiation codon of the coding sequence which is often marked by a cluster of small RNAs. Second, unusual tails are added post-transcriptionally at the 3′-extremity of all mRNAs. The nucleotide composition of the tails is distinct from that described in any other systems and can be partitioned between A/U-rich tails, predominantly composed of Adenosine and Uridine, and C-rich tails composed mostly of Cytidine. Based on 3′ RACE experiments, 22% of mRNAs present C-rich tails, some of them composed of up to 20 consecutive Cs. Polycytidylation is specific to mitochondria and occurs primarily on mRNAs. This unprecedented post-transcriptional modification seems to be a specific feature of the Chlorophyceae class of green algae and points out the existence of novel strategies in mitochondrial gene expression.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Untangling Heteroplasmy, Structure, and Evolution of an Atypical Mitochondrial Genome by PacBio Sequencing

The highly compact mitochondrial (mt) genome of terrestrial isopods (Oniscidae) presents two unusual features. First, several loci can individually encode two tRNAs, thanks to single nucleotide polymorphisms at anticodon sites. Within-individual variation (heteroplasmy) at these loci is thought to have been maintained for millions of years because individuals that do not carry all tRNA genes die, resulting in strong balancing selection. Second, the oniscid mtDNA genome comes in two conformations: a ∼14 kb linear monomer and a ∼28 kb circular dimer comprising two monomer units fused in palindrome. We hypothesized that heteroplasmy actually results from two genome units of the same dimeric molecule carrying different tRNA genes at mirrored loci. This hypothesis, however, contradicts the earlier proposition that dimeric molecules result from the replication of linear monomers-a process that should yield totally identical genome units within a dimer. To solve this contradiction, we used the SMRT (PacBio) technology to sequence mirrored tRNA loci in single dimeric molecules. We show that dimers do present different tRNA genes at mirrored loci; thus covalent linkage, rather than balancing selection, maintains vital variation at anticodons. We also leveraged unique features of the SMRT technology to detect linear monomers closed by hairpins and carrying noncomplementary bases at anticodons. These molecules contain the necessary information to encode two tRNAs at the same locus, and suggest new mechanisms of transition between linear and circular mtDNA. Overall, our analyses clarify the evolution of an atypical mt genome where dimerization counterintuitively enabled further mtDNA compaction.

The mitochondrial genomes of the acoelomorph worms Paratomella rubra, Isodiametra pulchra and Archaphanostoma ylvae

Acoels are small, ubiquitous - but understudied - marine worms with a very simple body plan. Their internal phylogeny is still not fully resolved, and the position of their proposed phylum Xenacoelomorpha remains debated. Here we describe mitochondrial genome sequences from the acoels Paratomella rubra and Isodiametra pulchra, and the complete mitochondrial genome of the acoel Archaphanostoma ylvae. The P. rubra and A. ylvae sequences are typical for metazoans in size and gene content. The larger I. pulchra  mitochondrial genome contains both ribosomal genes, 21 tRNAs, but only 11 protein-coding genes. We find evidence suggesting a duplicated sequence in the I. pulchra mitochondrial genome. The P. rubra, I. pulchra and A. ylvae mitochondria have a unique genome organisation in comparison to other metazoan mitochondrial genomes. We found a large degree of protein-coding gene and tRNA overlap with little non-coding sequence in the compact P. rubra genome. Conversely, the A. ylvae and I. pulchra genomes have many long non-coding sequences between genes, likely driving genome size expansion in the latter. Phylogenetic trees inferred from mitochondrial genes retrieve Xenacoelomorpha as an early branching taxon in the deuterostomes. Sequence divergence analysis between P. rubra sampled in England and Spain indicates cryptic diversity.

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of a Freshwater Mussel Potamilus alatus (Bivalvia: Unionidae)

Doubly uniparental inheritance (DUI) of mitochondrial DNA, found only in some bivalve families and characterized by the existence of gender-associated mtDNA lineages that are inherited through males (M-type) or females (F-type), is one of the very few exceptions to the general rule of strict maternal mtDNA inheritance in animals. M-type sequences are often undetected and hence still underrepresented in the GenBank, which hinders the progress of the understanding of the DUI phenomenon. We have sequenced and analyzed the complete M and F mitogenomes of a freshwater mussel, Potamilus alatus. The M-type was 493 bp longer (M = 16 560, F = 16 067 bp). Gene contents, order and the distribution of genes between L and H strands were typical for unionid mussels. Candidates for the two ORFan genes (forf and morf) were found in respective mitogenomes. Both mitogenomes had a very similar A+T bias: F = 61% and M = 62.2%. The M mitogenome-specific cox2 extension (144 bp) is much shorter than in other sequenced unionid mitogenomes (531-576 bp), which might be characteristic for the Potamilus genus. The overall topology of the phylogenetic tree is in very good agreement with the currently accepted phylogenetic relationships within the Unionidae: both studied sequences were placed within the Ambleminae subfamily clusters in the corresponding M and F clades.