Asymmetrical directional mutation pressure in the mitochondrial genome of mammals

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mitonuclear Interactions in the Maintenance of Mitochondrial Integrity

In eukaryotic cells, mitochondria originated in an α-proteobacterial endosymbiont. Although these organelles harbor their own genome, the large majority of genes, originally encoded in the endosymbiont, were either lost or transferred to the nucleus. As a consequence, mitochondria have become semi-autonomous and most of their processes require the import of nuclear-encoded components to be functional. Therefore, the mitochondrial-specific translation has evolved to be coordinated by mitonuclear interactions to respond to the energetic demands of the cell, acquiring unique and mosaic features. However, mitochondrial-DNA-encoded genes are essential for the assembly of the respiratory chain complexes. Impaired mitochondrial function due to oxidative damage and mutations has been associated with numerous human pathologies, the aging process, and cancer. In this review, we highlight the unique features of mitochondrial protein synthesis and provide a comprehensive insight into the mitonuclear crosstalk and its co-evolution, as well as the vulnerabilities of the animal mitochondrial genome.

Complete Mitochondrial Genomes of Two Toxin-Accumulated Nassariids (Neogastropoda: Nassariidae: Nassarius) and Their Implication for Phylogeny

The Indo-Pacific nassariids (genus Nassarius) possesses the highest diversity within the family Nassariidae. However, the previous shell or radula-based classification of Nassarius is quite confusing due to the homoplasy of certain morphological characteristics. The toxin accumulators Nassarius glans and Nassarius siquijorensis are widely distributed in the subtidal regions of the Indo-Pacific Ocean. In spite of their biological significance, the phylogenetic positions of N. glans and N. siquijorensis are still undetermined. In the present study, the complete mitochondrial genomes of N. glans and N. siquijorensis were sequenced. The present mitochondrial genomes were 15,296 and 15,337 bp in length, respectively, showing negative AT skews and positive GC skews as well as a bias of AT rich on the heavy strand. They contained 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and several noncoding regions, and their gene order was identical to most caenogastropods. Based on the nucleotide sequences combining 13 protein coding genes and two rRNA genes, a well-supported phylogeny of Nassarius was reconstructed, and several morphological synapomorphies were observed corresponding to the phylogenetic framework. In addition, the sister group relationship between N. variciferus and the remaining toxin-accumulated nassariids was determined, suggesting that the phylogeny might be related to their diet. The divergence time estimation analysis revealed a correlation between speciation events of nassariids and glacial cycles during the Pliocene-Pleistocene epoch.

The complete mitogenome of Arion vulgaris Moquin-Tandon, 1855 (Gastropoda: Stylommatophora): mitochondrial genome architecture, evolution and phylogenetic considerations within Stylommatophora

Stylommatophora is one of the most speciose orders of Gastropoda, including terrestrial snails and slugs, some of which are economically important as human food, agricultural pests, vectors of parasites or due to invasiveness. Despite their great diversity and relevance, the internal phylogeny of Stylommatophora has been debated. To date, only 34 stylommatophoran mitogenomes were sequenced. Here, the complete mitogenome of an invasive pest slug, Arion vulgaris Moquin-Tandon, 1855 (Stylommatophora: Arionidae), was sequenced using next generation sequencing, analysed and compared with other stylommatophorans. The mitogenome of A. vulgaris measures 14,547 bp and contains 13 protein-coding, two rRNA, 22 tRNA genes, and one control region, with an A + T content of 70.20%. All protein coding genes (PCGs) are initiated with ATN codons except for COX1, ND5 and ATP8 and all are ended with TAR or T-stop codons. All tRNAs were folded into a clover-leaf secondary structure except for trnC and trnS1 (AGN). Phylogenetic analyses confirmed the position of A. vulgaris within the superfamily Arionoidea, recovered a sister group relationship between Arionoidea and Orthalicoidea, and supported monophyly of all currently recognized superfamilies within Stylommatophora except for the superfamily Helicoidea. Initial diversification time of the Stylommatophora was estimated as 138.55 million years ago corresponding to Early Cretaceous. The divergence time of A. vulgaris and Arion rufus (Linnaeus, 1758) was estimated as 15.24 million years ago corresponding to one of Earth's most recent, global warming events, the Mid-Miocene Climatic Optimum. Furthermore, selection analyses were performed to investigate the role of different selective forces shaping stylommatophoran mitogenomes. Although purifying selection is the predominant selective force shaping stylommatophoran mitogenomes, six genes (ATP8, COX1, COX3, ND3, ND4 and ND6) detected by the branch-specific aBSREL approach and three genes (ATP8, CYTB and ND4L) detected by codon-based BEB, FUBAR and MEME approaches were exposed to diversifying selection. The positively selected substitutions at the mitochondrial PCGs of stylommatophoran species seems to be adaptive to environmental conditions and affecting mitochondrial ATP production or protection from reactive oxygen species effects. Comparative analysis of stylommatophoran mitogenome rearrangements using MLGO revealed conservatism in Stylommatophora; exceptions refer to potential apomorphies for several clades including rearranged orders of trnW-trnY and of trnE-trnQ-rrnS-trnM-trnL2-ATP8-trnN-ATP6-trnR clusters for the genus Arion. Generally, tRNA genes tend to be rearranged and tandem duplication random loss, transitions and inversions are the most basic mechanisms shaping stylommatophoran mitogenomes.

Deamination gradients within codons after 1<->2 position swap predict amino acid hydrophobicity and parallel β-sheet conformational preference

Deaminations C->T and A->G are frequent mutations producing nucleotide content gradients across genomes proportional to singlestrandedness during replication/transcription. Hence, within single codons, deamination risks increase from first to third codon positions, while second codon positions are functionally most crucial. Here genetic codes are analyzed assuming that after anticodons protected codons from deaminations, first and second codon positions swapped (N2N1N3->N1N2N3), with lowest deamination risks for N2 in presumed primitive N2N1N3 codons. N2N1N3, not standard N1N2N3, codon structure minimizes deaminations inversely proportionally to cognate amino acid hydrophobicity and parallel betasheet conformational preference. For N1N2N3, deamination minimization increases with genetic code integration order of cognate amino acids: during the presumed N2N1N3->N1N2N3 codon structure transition, protein synthesis combined direct codon-amino acid interactions for late amino acids and tRNA-based translation for early amino acids. Hence N2N1N3 codons would correspond to tRNA-free translation by spontaneous codon-amino acid affinities, and tRNA-mediated translation presumably caused N2N1N3->N1N2N3 swaps. Results show that rational, not arbitrary rules link codon and amino acid structures. Some analyses detect mitochondrial RNAs and peptides in public data corresponding to systematic position swaps, suggesting occasional swapping polymerase activity.

RNA Rings Strengthen Hairpin Accretion Hypotheses for tRNA Evolution: A Reply to Commentaries by Z.F. Burton and M. Di Giulio

Theoretical minimal RNA ring design ensures coding over the shortest length once for each coding signal (start and stop codons, and each amino acid) and their hairpin configuration. These constraints define 25 RNA rings which surprisingly resemble ancestral tRNA loops, suggesting commonalities between RNA ring design and proto-tRNAs. RNA rings share several other properties with tRNAs, suggesting that primordial RNAs were multifunctional peptide coding sequences and structural RNAs. Two hypotheses, respectively, by M. Di Giulio and Z.F. Burton, derived from cloverleaf structural symmetries suggest that two and three, respectively, stem-loop hairpins agglutinated into tRNAs. Their authors commented that their respective structure-based hypotheses reflect better tRNA structure than RNA rings. Unlike these hypotheses, RNA ring design uses no tRNA-derived information, rendering model predictive power comparisons senseless. Some analyses of RNA ring primary and secondary structures stress RNA ring splicing in their predicted anticodon's midst, indicating ancestrality of split tRNAs, as the two-piece model predicts. Advancement of knowledge, rather than of specific hypotheses, gains foremost by examining independent hypotheses for commonalities, and only secondarily for discordances. RNA rings mimick ancestral biomolecules including tRNAs, and their evolution, and constitute an interesting synthetic system for early prebiotic evolution tests/simulations.

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.

The primordial tRNA acceptor stem code from theoretical minimal RNA ring clusters

BACKGROUND

Theoretical minimal RNA rings code by design over the shortest length once for each of the 20 amino acids, a start and a stop codon, and form stem-loop hairpins. This defines at most 25 RNA rings of 22 nucleotides. As a group, RNA rings mimick numerous prebiotic and early life biomolecular properties: tRNAs, deamination gradients and replication origins, emergence of codon preferences for the natural circular code, and contents of early protein coding genes. These properties result from the RNA ring's in silico design, based mainly on coding nonredundancy among overlapping translation frames, as the genetic code's codon-amino acid assignments determine. RNA rings resemble ancestral tRNAs, defining RNA ring anticodons and corresponding cognate amino acids. Surprisingly, all examined RNA ring properties coevolve with genetic code integration ranks of RNA ring cognates, as if RNA rings mimick prebiotic and early life evolution.

METHODS

Distances between RNA rings were calculated using different evolutionary models. Associations between these distances and genetic code evolutionary hypotheses detect evolutionary models best describing RNA ring diversification.

RESULTS

Here pseudo-phylogenetic analyses of RNA rings produce clusters corresponding to the primordial code in tRNA acceptor stems, more so when substitution matrices from neutrally evolving pseudogenes are used rather than from functional protein coding genes reflecting selection for conserving amino acid properties.

CONCLUSIONS

Results indicate RNA rings with recent cognates evolved from those with early cognates. Hence RNA rings, as designed by the genetic code's structure, simulate tRNA stem evolution and prebiotic history along neutral chemistry-driven mutation regimes.

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.

The complete mitochondrial genome of Harpago chiragra and Lambis lambis (Gastropoda: Stromboidea): implications on the Littorinimorpha phylogeny

The complete mitochondrial genomes of Harpago chiragra and Lambis lambis (Strombidae) were determined with the size of 15,460 bp and 15,481 bp, respectively, and both sequences contained 13 protein-coding genes, 22 tRNAs, and two rRNAs. H. chiragra and L. lambis have similar mitochondrial features, corresponding to typical gastropod mitochondrial genomes, such as the conserved gene order, a high A + T content (66.22% for H. chiragra and 66.10% for L. lambis), and preference for A + T-rich codons. The start or termination codon of same protein-coding gene in H. chiragra was consistent with that in L. lambis, except for the termination codon of cox1 gene (TAG for H. chiragra and TAA for L. lambis) and the start codon of nad4 (GTG for H. chiragra and ATG for L. lambis). Pairwise sequence alignments detected different degrees of variations in H. chiragra and L. lambis mitochondrial genomes; and the two species had lower levels of genetic distance (0.202 for nucleotide sequence) and closest relationships as compared to Strombus gigas and Oncomelania hupensis. The 13 partitioned nucleotide sequences of protein coding genes of H. chiragra and L. lambis were aligned with representatives of the main lineages of gastropods and their phylogenetic relationships were inferred. H. chiragra and L. lambis share the same gene order as Littorinimorpha species, except Vermetoidea, which demonstrate a gene rearrangement in species. The reconstructed phylogeny supports three major clades within Littorinimorpha: 1) Stromboidea, Tonnoidea, Littorinoidea, and Naticoidea, 2) Rissooidea and Truncatelloidea, and 3) Vermetoidea. In addition, a relaxed molecular clock calibrated with fossils dated the diversification of Strombidae near 112 (44–206) Mya and a possible radiation is detected to occur between 45–75 Mya, providing implications to understand the Cenozoic replacement event (65–135 Mya) of Aporrhaidae by Strombidae.

Two novel cricetine mitogenomes: Insight into the mitogenomic characteristics and phylogeny in Cricetinae (Rodentia: Cricetidae)

Both Cricetus cricetus and Phodopus sungorus mitochondrial genomes (mitogenomes) were sequenced and elaborated for the first time in the present study. Their mitogenomes contained 37 genes and showed typical characteristics of the vertebrate mitogenome. Comparative analysis of 10 cricetine mitogenomes indicated that they shared similar characteristics with those of other cricetines in terms of genes arrangement, nucleotide composition, codon usage, tRNA structure, nucleotide skew and the origin of replication of light strand. Phylogenetic relationship of the subfamily Cricetinae was reconstructed using mitogenomes data with the methods of Bayesian Inference and Maximum Likelihood. Phylogenetic analysis indicated that Cricetulus kamensis was at basal position and phylogenetically distant from all other Cricetulus species but had a close relationship with the group of Phodopus, and supported that the genus Urocricetus deserved as a separate genus rank. The phylogenetic status of Tscherskia triton represented a separate clade corresponding to a diversified cricetine lineage (Cricetulus, Allocricetulus, and Cricetus).

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

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

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 Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Localized Context-Dependent Effects of the "Ambush" Hypothesis: More Off-Frame Stop Codons Downstream of Shifty Codons

The ambush hypothesis speculates that off-frame stop codons increase translational efficiency after ribosomal frameshifts by stopping early frameshifted translation. Some evidences fit this hypothesis: (1) synonymous codon usages increase with their potential contribution to off-frame stops; (2) the genetic code assigns frequent amino acids to codon families contributing to off-frame stops; (3) positive biases for off-frame stops (AT rich) occur despite adverse nucleotide (GC) biases; and (4) mitochondrial off-frame stop codon densities increase with ribosomal structural instability, potential proxy of frameshift frequencies. In this study, analyses of vertebrate mitogenes and tRNA synthetase genes from all superkingdoms and viruses test a new prediction of the ambush hypothesis: sequences immediately downstream of frameshift-inducing homopolymer codons (AAA, CCC, GGG, and TTT) are off-frame stop rich. Codons immediately downstream of homopolymer codons form more than average off-frame stops, biases are stronger than for corresponding upstream distances and for any other group of synonymous codons. Sequences downstream of that high-density region are off-frame stop depleted. This decrease suggests that off-frame stops, combined with suppressor tRNAs regulate translation of overlapping coding sequences. Results show the predictive power of the ambush hypothesis, from macroevolutionary (genetic code structure) to detailed gene sequence anatomy.

Mutation signatures in germline mitochondrial genome provide insights into human mitochondrial evolution and disease

Variations in mitochondrial DNA (mtDNA) have been fundamental for understanding human evolution and are causative for a plethora of inherited mitochondrial diseases, but the mutation signatures of germline mtDNA and their value in understanding mitochondrial pathogenicity remain unknown. Here, we carried out a systematic analysis of mutation patterns in germline mtDNA based on 97,566 mtDNA variants from 45,494 full-length sequences and revealed a highly non-stochastic and replication-coupled mutation signature characterized by nucleotide-specific mutation pressure (G > T>A > C) and position-specific selection pressure, suggesting the existence of an intensive mutation-selection interplay in germline mtDNA. We provide evidence that this mutation-selection interplay has strongly shaped the mtDNA sequence during evolution, which not only manifests as an oriented alteration of amino acid compositions of mitochondrial encoded proteins, but also explains the long-lasting mystery of CpG depletion in mitochondrial genome. Finally, we demonstrated that these insights may be integrated to better understand the pathogenicity of disease-implicated mitochondrial variants.

G-quadruplex dynamics contribute to regulation of mitochondrial gene expression

Single-stranded DNA or RNA sequences rich in guanine (G) can adopt non-canonical structures known as G-quadruplexes (G4). Mitochondrial DNA (mtDNA) sequences that are predicted to form G4 are enriched on the heavy-strand and have been associated with formation of deletion breakpoints. Increasing evidence supports the ability of mtDNA to form G4 in cancer cells; however, the functional roles of G4 structures in regulating mitochondrial nucleic acid homeostasis in non-cancerous cells remain unclear. Here, we demonstrate by live cell imaging that the G4-ligand RHPS4 localizes primarily to mitochondria at low doses. We find that low doses of RHPS4 do not induce a nuclear DNA damage response but do cause an acute inhibition of mitochondrial transcript elongation, leading to respiratory complex depletion. We also observe that RHPS4 interferes with mtDNA levels or synthesis both in cells and isolated mitochondria. Importantly, a mtDNA variant that increases G4 stability and anti-parallel G4-forming character shows a stronger respiratory defect in response to RHPS4, supporting the conclusion that mitochondrial sensitivity to RHPS4 is G4-mediated. Taken together, our results indicate a direct role for G4 perturbation in mitochondrial genome replication, transcription processivity, and respiratory function in normal cells.

Phenomenological Load on Model Parameters Can Lead to False Biological Conclusions

When a substitution model is fitted to an alignment using maximum likelihood, its parameters are adjusted to account for as much site-pattern variation as possible. A parameter might therefore absorb a substantial quantity of the total variance in an alignment (or more formally, bring about a substantial reduction in the deviance of the fitted model) even if the process it represents played no role in the generation of the data. When this occurs, we say that the parameter estimate carries phenomenological load (PL). Large PL in a parameter estimate is a concern because it not only invalidates its mechanistic interpretation (if it has one) but also increases the likelihood that it will be found to be statistically significant. The problem of PL was not identified in the past because most off-the-shelf substitution models make simplifying assumptions that preclude the generation of realistic levels of variation. In this study, we use the more realistic mutation-selection framework as the basis of a generating model formulated to produce data that mimic an alignment of mammalian mitochondrial DNA. We show that a parameter estimate can carry PL when 1) the substitution model is underspecified and 2) the parameter represents a process that is confounded with other processes represented in the data-generating model. We then provide a method that can be used to identify signal for the process that a given parameter represents despite the existence of PL.

Mitochondrial DNA Integrity: Role in Health and Disease

As the primary cellular location for respiration and energy production, mitochondria serve in a critical capacity to the cell. Yet, by virtue of this very function of respiration, mitochondria are subject to constant oxidative stress that can damage one of the unique features of this organelle, its distinct genome. Damage to mitochondrial DNA (mtDNA) and loss of mitochondrial genome integrity is increasingly understood to play a role in the development of both severe early-onset maladies and chronic age-related diseases. In this article, we review the processes by which mtDNA integrity is maintained, with an emphasis on the repair of oxidative DNA lesions, and the cellular consequences of diminished mitochondrial genome stability.

Potential Roles for G-Quadruplexes in Mitochondria

Some DNA or RNA sequences rich in guanine (G) nucleotides can adopt noncanonical conformations known as G-quadruplexes (G4). In the nuclear genome, G4 motifs have been associated with genome instability and gene expression defects, but they are increasingly recognized to be regulatory structures. Recent studies have revealed that G4 structures can form in the mitochondrial genome (mtDNA) and potential G4 forming sequences are associated with the origin of mtDNA deletions. However, little is known about the regulatory role of G4 structures in mitochondria. In this short review, we will explore the potential for G4 structures to regulate mitochondrial function, based on evidence from the nucleus.

Characterization of the mitochondrial genome of Arge bella Wei & Du sp. nov. (Hymenoptera: Argidae)

We describe Arge bella Wei & Du sp. nov., a large and beautiful species of Argidae from south China, and report its mitochondrial genome based on high-throughput sequencing data. We present the gene order, nucleotide composition of protein-coding genes (PCGs), and the secondary structures of RNA genes. The nearly complete mitochondrial genome of A. bella has a length of 15,576 bp and a typical set of 37 genes (22 tRNAs, 13 PCGs, and 2 rRNAs). Three tRNAs are rearranged in the A. bella mitochondrial genome as compared to the ancestral type in insects: trnM and trnQ are shuffled, while trnW is translocated from the trnW-trnC-trnY cluster to a location downstream of trnI. All PCGs are initiated by ATN codons, and terminated with TAA, TA or T as stop codons. All tRNAs have a typical cloverleaf secondary structure, except for trnS1. H821 of rrnS and H976 of rrnL are redundant. A phylogenetic analysis based on mitochondrial genome sequences of A. bella, 21 other symphytan species, two apocritan representatives, and four outgroup taxa supports the placement of Argidae as sister to the Pergidae within the symphytan superfamily Tenthredinoidea.

New Insight into Parrots' Mitogenomes Indicates That Their Ancestor Contained a Duplicated Region

Mitochondrial genomes of vertebrates are generally thought to evolve under strong selection for size reduction and gene order conservation. Therefore, a growing number of mitogenomes with duplicated regions changes our view on the genome evolution. Among Aves, order Psittaciformes (parrots) is especially noteworthy because of its large morphological, ecological, and taxonomical diversity, which offers an opportunity to study genome evolution in various aspects. Former analyses showed that tandem duplications comprising the control region with adjacent genes are restricted to several lineages in which the duplication occurred independently. However, using an appropriate polymerase chain reaction strategy, we demonstrate that early diverged parrot groups contain mitogenomes with the duplicated region. These findings together with mapping duplication data from other mitogenomes onto parrot phylogeny indicate that the duplication was an ancestral state for Psittaciformes. The state was inherited by main parrot groups and was lost several times in some lineages. The duplicated regions were subjected to concerted evolution with a frequency higher than the rate of speciation. The duplicated control regions may provide a selective advantage due to a more efficient initiation of replication or transcription and a larger number of replicating genomes per organelle, which may lead to a more effective energy production by mitochondria. The mitogenomic duplications were associated with phenotypic features and parrots with the duplicated region can live longer, show larger body mass as well as predispositions to a more active flight. The results have wider implications on the presence of duplications and their evolution in mitogenomes of other avian groups.

Accelerated Laboratory Evolution Reveals the Influence of Replication on the GC Skew in Escherichia coli

Most bacterial genomes display contrasting strand asymmetry in a variety of features, such as nucleotide composition and gene orientation, of the two replichores separated by the replication origin and terminus. The cause for the polarization is often attributed to mutations arising from the asymmetric replication machinery. Notably, a base compositional bias known as a GC skew is focused on as a footprint of the bacterial genome evolution driven by DNA replication. Previously, although a replication driven mutation pattern responsible for the GC skew formation or the related mathematical models have been well reported, an exact impact of the replication-related elements on the genomic structure is yet actively debated, and not confirmed experimentally. However, the GC skew formation is very time consuming and challenging in the laboratory. We, therefore, used cytosine deaminase as a DNA mutator, and by monitoring the mutations during an accelerated laboratory evolution procedure with Illumina sequencing, we enabled the trial and error of the GC skew formation in high resolution. Using this technology, we succeeded in reconfirming the influence of bacterial replication machinery on the genomic structure at high resolution.

The mitochondrial genome of Greenidea psidii van der Goot (Hemiptera: Aphididae: Greenideinae) and comparisons with other Aphididae aphids

The complete mitochondrial genome of Greenidea psidii was sequenced and compared with the genomes of other aphid species. The G. psidii mitogenome is 16,202 bp long with an A + T content of 85.4%, comprising 37 genes arranged in the same order as the inferred insect ancestral arrangement, a control region and a special repeat region between trnE and trnF. G. psidii is the first representative possessing such a repeat region from the subfamily Greenideinae. These repeat motifs found in Aphidinae, Eriosomatinae and Greenideinae are species-specific, differing in nucleotide composition, length and copy number across different species. All reported complete aphid mitogenomes are A- and C-skewed in the whole genomes and show opposite skewness for C and G between the majority and minority strands. The ratios of nonsynonymous to synonymous substitution rates indicated that the evolution of the aphid mitogenomes has been dominated by purifying selection. Based on whole mitogenome sequences, the phylogenetic relationships among all aphid species with complete mitogenomes were investigated. G. psidii was robustly clustered with Cervaphis quercus, confirming its affiliation with the subfamily Greenideinae.

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.

In Cold Blood: Compositional Bias and Positive Selection Drive the High Evolutionary Rate of Vampire Bats Mitochondrial Genomes

Mitochondrial genomes of animals have long been considered to evolve under the action of purifying selection. Nevertheless, there is increasing evidence that they can also undergo episodes of positive selection in response to shifts in physiological or environmental demands. Vampire bats experienced such a shift, as they are the only mammals feeding exclusively on blood and possessing anatomical adaptations to deal with the associated physiological requirements (e.g., ingestion of high amounts of liquid water and iron). We sequenced eight new chiropteran mitogenomes including two species of vampire bats, five representatives of other lineages of phyllostomids and one close outgroup. Conducting detailed comparative mitogenomic analyses, we found evidence for accelerated evolutionary rates at the nucleotide and amino acid levels in vampires. Moreover, the mitogenomes of vampire bats are characterized by an increased cytosine (C) content mirrored by a decrease in thymine (T) compared with other chiropterans. Proteins encoded by the vampire bat mitogenomes also exhibit a significant increase in threonine (Thr) and slight reductions in frequency of the hydrophobic residues isoleucine (Ile), valine (Val), methionine (Met), and phenylalanine (Phe). We show that these peculiar substitution patterns can be explained by the co-occurrence of both neutral (mutational bias) and adaptive (positive selection) processes. We propose that vampire bat mitogenomes may have been impacted by selection on mitochondrial proteins to accommodate the metabolism and nutritional qualities of blood meals.

Complete mitochondrial genome of bamboo grasshopper, Ceracris fasciata, and the phylogenetic analyses and divergence time estimation of Caelifera (Orthoptera)

The bamboo grasshopper Ceracris fasciata is regarded as a major pest species because of the damage it causes to bamboo, and its classification within the families and subfamilies of the suborder Caelifera remains unclear. Thus, we attempted to resolve these questions using molecular biology methods and analyses. Our results are as follows: (1) the complete mitochondrial genome of C. fasciata is 15,569 bp in length. The mitochondrial genome contains a standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region in the same order as those of the other analysed Caeliferan species. The putative start codon for the COX1 gene in C. fasciata is ACC, although it is not defined in other genes. The presence of tandem repeats of different sizes in the A + T-rich region may lead to size differences in other mitochondrial genomes. The mitochondrial genome of C. fasciata harbours the typical 37 genes and an A + T-rich region, and it shows similar characteristics to those of other grasshopper species. Characterization of the mitochondrial genome has enriched our knowledge of the mitochondrial genomes of Orthoptera around the world. Therefore, the phylogenetic relationships in Orthoptera can be re-examined. (2) In phylogenetic analyses, the monophyly of Orthoptera and its two suborders (Caelifera and Ensifera) has been consistently recovered based on most of the datasets selected, regardless of the optimal criteria. Our results do not support the monophyly of the subfamily Oedipodinae of Caelifera. We found that Phlaeoba albonema of the Acridinae is sorted into a clade with Ceracris in all our phylogenetic trees, and field experiments show that Phlaeoba always lives with Ceracris in the same ecotopes. Therefore, we suggest that Phlaeoba should be classified as a member of the Oedipodinae. We found that C. fasciata always clustered with Ceracris kiangsu, and both were sisters to Ceracris versicolor. Therefore, the genetic relationship between C. fasciata and C. kiangsu is closer than that between C. fasciata and C. versicolor. (3) The oldest estimated time of divergence of Ensifera in this context was determined to be 146.16 million years ago (Mya), or around the late Jurassic or early Cretaceous. We estimated that katydids (Grylloidea) likely diverged from other groups in the early Cretaceous. According to our divergence time analyses, we concluded that the ancestral Acrididae probably originated in the early Paleogene, and it is likely that the major diversification events happened at the middle Paleogene, well into the next geologic time. We estimated that crickets (Tettigoniidae) likely diverged from other groups in the early Cretaceous. Acrididae and Romaleinae group, Pyrgacrididae and Ommexechidae group, the youngest two clades we observed, were estimated to have diverged 58.79 Mya, between the middle and early Paleogene. C. versicolor is a sister to the group containing C. kiangsu and C. fasciata. First, C. versicolor diverged from the sister group (C. kiangsu + C. fasciata) around 44.81 Mya, and then the C. kiangsu and C. fasciata group separated at 43.04 Mya.

Protein Sequences Recapitulate Genetic Code Evolution

Several hypotheses predict ranks of amino acid assignments to genetic code's codons. Analyses here show that average positions of amino acid species in proteins correspond to assignment ranks, in particular as predicted by Juke's neutral mutation hypothesis for codon assignments. In all tested protein groups, including co- and post-translationally folding proteins, 'recent' amino acids are on average closer to gene 5' extremities than 'ancient' ones. Analyses of pairwise residue contact energies matrices suggest that early amino acids stereochemically selected late ones that stablilize residue interactions within protein cores, presumably producing 5'-late-to-3'-early amino acid protein sequence gradients. The gradient might reduce protein misfolding, also after mutations, extending principles of neutral mutations to protein folding. Presumably, in self-perpetuating and self-correcting systems like the genetic code, initial conditions produce similarities between evolution of the process (the genetic code) and 'ontogeny' of resulting structures (here proteins), producing apparent teleonomy between process and product.

Intraspecific mitochondrial genome comparison identified CYTB as a high-resolution population marker in a new pest Athetis lepigone

A new outbreak pest, Athetis lepigone (Möschler) (Lepidoptera: Noctuidae), has caused severe economic loss in maize crops in China. In order to conduct population genetics study with a more polymorphic and scientific mitochondrial marker, we sequenced the complete mitochondrial genomes of 13 different A. lepigone individuals. Intraspecific comparison of all PCGs showed that the NADH dehydrogenase and cytochrome b genes had the highest nucleotide diversity. We also found evidence of episodic positive selection on two amino acids, which are encoded by NADH dehydrogenase genes (ND3 and ND4L), against a background of widespread neutral selection of all other mitochondrial PCGs. The genetic divergence observed in this study indicated that the cytochrome b gene (CYTB) is better than COI at recovering population structure. The preliminary population genetic analysis illustrated strong gene flow among A. lepigone populations in China. Our study provides basic information for further research on population genetics of A. lepigone.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

Characterisation of the complete mitochondrial genome of Helice wuana (Grapsoidea: Varunidae) and comparison with other Brachyuran crabs

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Herein, we sequenced, annotated, and characterised the mitogenome of the crab Helice wuana to better understand its molecular evolution and phylogeny. The 16,359bp mitogenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region. The genome composition is highly A+T biased 68.42%, and exhibits a negative AT-skew (-0.036) and GC-skew (-0.269) among Brachyura crabs. Gene rearrangements were detected, as was tandem duplication followed by random loss, which explains the translocation of mitochondrial genes. Phylogenetic analysis showed that H. wuana and H. tientsinensis clustered on one branch with high nodal support values. These results confirm that the placement of H. wuana within the Varunidae family of Thoracotrematan crabs. This study will provided a better understanding for gene rearrangements and crab evolution in the further.

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

The complete mitochondrial DNA sequence of the pantropical earthworm Pontoscolex corethrurus (Rhinodrilidae, Clitellata): Mitogenome characterization and phylogenetic positioning

Pontoscolexcorethrurus (Müller, 1857) plays an important role in tropical soil ecosystems and has been widely used as an animal model for a large variety of ecological studies in particular due to its common presence and generally high abundance in human-disturbed tropical soils. In this study we describe the complete mitochondrial genome of the peregrine earthworm P.corethrurus. This is the first record of a mitochondrial genome within the Rhinodrilidae family. Its mitochondrial genome is 14 835 bp in length containing 37 genes (13 protein-coding genes (PCG) 2 rRNA genes and 22 tRNA genes). It has the same gene content and structure as in other sequenced earthworms but unusual among invertebrates it hasseveral overlapping open reading frames. All genes are encoded on the same strand. Most of the PCGs use ATG as the start codon except for ND3 which uses GTG as the start codon. The A+T content of the mitochondrial genome is 59.9% (31.8% A 28.1% T 14.6% G and 25.6% for C). The annotated genome sequence has been deposited in GenBank under the accession number KT988053.

The complete mitochondrial genome of Euproctis similis (Lepidoptera: Noctuoidea: Erebidae) and phylogenetic analysis

The mitochondrial genome (mitogenome) can provide information for phylogenetic analyses and evolutionary biology. We sequenced, annotated, and characterized the mitogenome of Euproctis similis. The complete mitogenome is 15,437bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a control region (A+T-rich region). The A+T content in the mitogenome was 80.16%. All PCGs use standard ATN as a start codon, with the exception of cytochrome c coxidase 1 (cox1) with CGA. A gene rearrangement (trnM) was found. All transfer RNA (tRNA) genes have a typical clover-leaf structure except for trnS1 (AGN). Phylogenetic analysis was performed using Bayesian Inference and Maximum Likelihood based on the amino acid and nucleotide sequences of 13 mitochondrial PCGs. The well-supported phylogenetic relationships can be generally described as: Notodontidae+(Erebidae+(Nolidae+(Euteliidae+Noctuidae))). The tree support that E. similis shares a close ancestry with Erebidae insects. Our results indicate that Erebidae is a sister group to the other families (Euteliidae, Nolidae, and Noctuidae).

Optimization of amino acid replacement costs by mutational pressure in bacterial genomes

Mutations are considered a spontaneous and random process, which is important component of evolution because it generates genetic variation. On the other hand, mutations are deleterious leading to non-functional genes and energetically costly repairs. Therefore, one can expect that the mutational pressure is optimized to simultaneously generate genetic diversity and preserve genetic information. To check if empirical mutational pressures are optimized in these ways, we compared matrices of nucleotide mutation rates derived from bacterial genomes with their best possible alternatives that minimized or maximized costs of amino acid replacements associated with differences in their physicochemical properties (e.g. hydropathy and polarity). It should be noted that the studied empirical nucleotide substitution matrices and the costs of amino acid replacements are independent because these matrices were derived from sites free of selection on amino acid properties and the amino acid costs assumed only amino acid physicochemical properties without any information about mutation at the nucleotide level. Obtained results indicate that the empirical mutational matrices show a tendency to minimize costs of amino acid replacements. It implies that bacterial mutational pressures can evolve to decrease consequences of amino acid substitutions. However, the optimization is not full, which enables generation of some genetic variability.

Complete mitochondrial genomes of Trisidos kiyoni and Potiarca pilula: Varied mitochondrial genome size and highly rearranged gene order in Arcidae

We present the complete mitochondrial genomes (mitogenomes) of Trisidos kiyoni and Potiarca pilula, both important species from the family Arcidae (Arcoida: Arcacea). Typical bivalve mtDNA features were described, such as the relatively conserved gene number (36 and 37), a high A + T content (62.73% and 61.16%), the preference for A + T-rich codons, and the evidence of non-optimal codon usage. The mitogenomes of Arcidae species are exceptional for their extraordinarily large and variable sizes and substantial gene rearrangements. The mitogenome of T. kiyoni (19,614 bp) and P. pilula (28,470 bp) are the two smallest Arcidae mitogenomes. The compact mitogenomes are weakly associated with gene number and primarily reflect shrinkage of the non-coding regions. The varied size in Arcidae mitogenomes reflect a dynamic history of expansion. A significant positive correlation is observed between mitogenome size and the combined length of cox1-3, the lengths of Cytb, and the combined length of rRNAs (rrnS and rrnL) (P < 0.001). Both protein coding genes (PCGs) and tRNA rearrangements is observed in P. pilula and T. kiyoni mitogenomes. This analysis imply that the complicated gene rearrangement in mitochondrial genome could be considered as one of key characters in inferring higher-level phylogenetic relationship of Arcidae.

Body mass-corrected molecular rate for bird mitochondrial DNA

Mitochondrial DNA remains one of the most widely used molecular markers to reconstruct the phylogeny and phylogeography of closely related birds. It has been proposed that bird mitochondrial genomes evolve at a constant rate of ~0.01 substitution per site per million years, that is that they evolve according to a strict molecular clock. This molecular clock is often used in studies of bird mitochondrial phylogeny and molecular dating. However, rates of mitochondrial genome evolution vary among bird species and correlate with life history traits such as body mass and generation time. These correlations could cause systematic biases in molecular dating studies that assume a strict molecular clock. In this study, we overcome this issue by estimating corrected molecular rates for birds. Using complete or nearly complete mitochondrial genomes of 475 species, we show that there are strong relationships between body mass and substitution rates across birds. We use this information to build models that use bird species' body mass to estimate their substitution rates across a wide range of common mitochondrial markers. We demonstrate the use of these corrected molecular rates on two recently published data sets. In one case, we obtained molecular dates that are twice as old as the estimates obtained using the strict molecular clock. We hope that this method to estimate molecular rates will increase the accuracy of future molecular dating studies in birds.

Comparative Large-Scale Mitogenomics Evidences Clade-Specific Evolutionary Trends in Mitochondrial DNAs of Bivalvia

Despite the figure of complete bivalve mitochondrial genomes keeps growing, an assessment of the general features of these genomes in a phylogenetic framework is still lacking, despite the fact that bivalve mitochondrial genomes are unusual under different aspects. In this work, we constructed a dataset of one hundred mitochondrial genomes of bivalves to perform the first systematic comparative mitogenomic analysis, developing a phylogenetic background to scaffold the evolutionary history of the class' mitochondrial genomes. Highly conserved domains were identified in all protein coding genes; however, four genes (namely, atp6, nad2, nad4L, and nad6) were found to be very divergent for many respects, notwithstanding the overall purifying selection working on those genomes. Moreover, the atp8 gene was newly annotated in 20 mitochondrial genomes, where it was previously declared as lacking or only signaled. Supernumerary mitochondrial proteins were compared, but it was possible to find homologies only among strictly related species. The rearrangement rate on the molecule is too high to be used as a phylogenetic marker, but here we demonstrate for the first time in mollusks that there is correlation between rearrangement rates and evolutionary rates. We also developed a new index (HERMES) to estimate the amount of mitochondrial evolution. Many genomic features are phylogenetically congruent and this allowed us to highlight three main phases in bivalve history: the origin, the branching of palaeoheterodonts, and the second radiation leading to the present-day biodiversity.

Decreased Mitochondrial Mutagenesis during Transformation of Human Breast Stem Cells into Tumorigenic Cells

Rare stochastic mutations may accumulate during dormancy of stem-like cells, but technical limitations in DNA sequencing have limited exploring this possibility. In this study, we employed a recently established deep-sequencing method termed Duplex Sequencing to conduct a genome-wide analysis of mitochondrial (mt) DNA mutations in a human breast stem cell model that recapitulates the sequential stages of breast carcinogenesis. Using this method, we found significant differences in mtDNA among normal stem cells, immortal/preneoplastic cells, and tumorigenic cells. Putative cancer stem-like cell (CSC) populations and mtDNA copy numbers increased as normal stem cells become tumorigenic cells. Transformed cells exhibited lower rare mutation frequencies of whole mtDNA than did normal stem cells. The predicted mtDNA rare mutation pathogenicity was significantly lower in tumorigenic cells than normal stem cells. Major rare mutation types in normal stem cells are C>T/G>A and T>C/A>G transitions, while only C>T/G>A are major types in transformed cells. We detected a total of 1,220 rare point mutations, 678 of which were unreported previously. With only one possible exception (m10342T>C), we did not find specific mutations characterizing mtDNA in human breast CSCs; rather, the mitochondrial genome of CSCs displayed an overall decrease in rare mutations. On the basis of our work, we suggest that this decrease (in particular T>C/A>G transitions), rather than the presence of specific mitochondrial mutations, may constitute an early biomarker for breast cancer detection. Our findings support the hypothesis that the mitochondrial genome is altered greatly as a result of the transformation of normal stem cells to CSCs, and that mtDNA mutation signatures may aid in delineating normal stem cells from CSCs. Cancer Res; 76(15); 4569-78. ©2016 AACR.