Asymmetrical directional mutation pressure in the mitochondrial genome of mammals

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Characterization of the mitochondrial genomes of the Mexican endemic bats Corynorhinus mexicanus and Corynorhinus leonpaniaguae (Chiroptera: Vespertilionidae)

BACKGROUND

The genus Corynorhinus is composed of four recognized species: C. rafinesquii, C. townsendii, C. mexicanus, and C. leonpaniaguae, the latter two being endemic to Mexico. According to the IUCN, C. mexicanus is considered "Near Threatened", as its populations are dwindling and habitats are affected by anthropogenic disturbance. Corynorhinus leonpaniaguae has not been assigned to an IUCN Red List risk category due to its recent description.

METHODS AND RESULTS

In this study, the mitochondrial genomes of C. mexicanus and C. leonpaniaguae were assembled and characterized in detail. The mitochondrial genomes (mtDNA) of C. mexicanus and C. leonpaniaguae have lengths of 16,470 and 16,581 bp respectively, with a predominant nucleotide usage of adenine (31.670% and 31.729%, respectively) and thymine (26.15% and 26.18%, respectively). The mtDNA of C. mexicanus and C. leonpaniaguae is composed of 37 coding and non-coding elements: 22 transfer RNAs (tRNA), 13 protein-coding genes (PCGs), two ribosomal RNAs and a non-coding region, the control region, which has a length of 933 bp and 1,149 bp, respectively. All tRNAs exhibited a cloverleaf secondary structure, with the exception of trn-Ser1 which showed a deletion of the dihydrouridine arm in the two species. All PCGs are subjected to purifying selection, with atp8 being the gene showing the highest Ka/Ks value.

CONCLUSIONS

These are the first whole mitogenomic resources developed for C. mexicanus and C. leonpaniaguae and enhance our knowledge of the ecology of these species and aid in their conservation.

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.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Four new complete mitochondrial genomes of Gobioninae fishes (Teleostei: Cyprinidae) and their phylogenetic implications

The subfamily Gobioninae is one of the most diverse fish groups within Cyprinidae. Their taxonomy and phylogenetic relationships are not completely resolved. In this study, the complete mitochondrial genomes (mitogenome) of four Gobioninae species (Microphysogobio elongatus, Microphysogobio chinssuensis, Gobio rivuloides and Rhinogobio nasutus) were sequenced and compared. The mitogenomes of four species ranges from 16603 bp to 16609 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. Most PCGs had significant codon usage bias. Except for the tRNASer (GCT), all the nucleotide substitutions of tRNA loops higher than the stems could fold into a stable secondary structure. The nucleotide compositions of Gobioninae mitogenome were biased toward A/T, and NAD4 was subjected to low purification selection and had a faster evolution rate among 13 PCGs. Bayesian inference and maximum likelihood phylogenetic analyses showed the consistent results. The four sequenced species clustered together with their congener species. However, more samples and mitogenome data are needed to untangle the phylogenetic relationships among genera Microphysogobio, Romanogobio, Hugobio, Biwia and Platysmacheilus.

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

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

In silico prospecting of the mtDNA of Macrobrachium amazonicum from transcriptome data

BACKGROUND

Macrobrachium amazonicum is a freshwater prawn widely distributed in South America that is undergoing speciation, so the denomination "M. amazonicum complex" is used for it. The mitochondrial cytochrome c oxidase subunit I (COI) gene has been used to elucidate this speciation, but heteroplasmies and pseudogenes have been recorded, making separation difficult. Obtaining genes from cDNA (RNA) rather than genomic DNA is an effective tool to mitigate those two types of occurrences. The aim of this study was to assemble in silico the mitochondrial DNA (mtDNA) of the Amazonian coastal population of M. amazonicum inhabiting the state of Pará.

RESULTS

Sequences were obtained from the prawn's transcriptome using the de novo approach. Six libraries of cDNA from the androgen gland, hepatopancreas, and muscle tissue were used. The mtDNA of M. amazonicum was 14,960 bp in length. It contained 13 protein-coding genes, 21 complete transfer RNAs, and the 12S and 16S subunits of ribosomal RNA. All regions were found on the light strand except tRNAGln, which was on the heavy strand. The control region (D-loop) was not recovered, making for a gap of 793 bp. The cladogram showed the formation of the well-defined Macrobrachium clade, with high support value in the established branches (91-100). The three-dimensional spatial conformation of the mtDNA-encoded proteins showed that most of them were mainly composed of major α-helices that typically shows in those proteins inserted in the membrane (mitochondrial).

CONCLUSIONS

It was possible to assemble a large part of the mitochondrial genome of M. amazonicum in silico using data from other genomes deposited in GenBank and to validate it through the similarities between its COI and 16S genes and those from animals of the same region deposited in GenBank. Depositing the M. amazonicum mtDNA sequences in GenBank may help solve the taxonomic problems recorded for the species, in addition to providing complete sequences of candidate coding genes for use as biomarkers in ecological studies.

Mitochondrial quality control in health and cardiovascular diseases

Cardiovascular diseases (CVDs) are one of the primary causes of mortality worldwide. An optimal mitochondrial function is central to supplying tissues with high energy demand, such as the cardiovascular system. In addition to producing ATP as a power source, mitochondria are also heavily involved in adaptation to environmental stress and fine-tuning tissue functions. Mitochondrial quality control (MQC) through fission, fusion, mitophagy, and biogenesis ensures the clearance of dysfunctional mitochondria and preserves mitochondrial homeostasis in cardiovascular tissues. Furthermore, mitochondria generate reactive oxygen species (ROS), which trigger the production of pro-inflammatory cytokines and regulate cell survival. Mitochondrial dysfunction has been implicated in multiple CVDs, including ischemia-reperfusion (I/R), atherosclerosis, heart failure, cardiac hypertrophy, hypertension, diabetic and genetic cardiomyopathies, and Kawasaki Disease (KD). Thus, MQC is pivotal in promoting cardiovascular health. Here, we outline the mechanisms of MQC and discuss the current literature on mitochondrial adaptation in CVDs.

The complete mitochondrial genome of Morishitium polonicum (Trematoda, Cyclocoelidae) and its phylogenetic implications

Trematodes can adversely impact the health and survival of wild animals. The trematode family Cyclocoelidae, which includes large digenean bird parasites, lacks molecular analysis, and reclassifications have not been supported. This study produced the first fully assembled and annotated mitochondrial genome sequence for the trematode Morishitium polonicum. The whole length of the M. polonicum (GenBank accession number: OP930879) mitogenome is 14083 bp, containing 22 transfer ribonucleic acids (tRNAs), 2 ribosomal RNAs (rRNAs, rrnL and rrnS), and a noncoding control section (D-loop) 13777 to 13854 bp in length. The 12 PCG areas have 3269 codons and a total length of 10053 bp, which makes up 71.38% of the mitochondrial genome's overall sequence. Most (10/12) of the PCGs that code for proteins begin with ATG, while the nad4L and nad1 genes have a GTG start codon. Phylogenetic analysis using the concatenated nucleotide sequences of 12 PCGs, and the ML tree analysis results showed that M. polonicum is more closely related to with Echinostomatidae and Fasciolidae, which indicates that the family Cyclocoelidae is more closely associated with Echinochasmidae. This study provides mtDNA information, and analysis of mitogenomic structure and evolution. Moreover, we aimed to understand the phylogenetic relationships of this fluke.

Characterization of the mitochondrial Huso huso genome and new aspects of its organization in the presence of tandem repeats in 12S rRNA

BACKGROUND

The sturgeon group has been economically significant worldwide due to caviar production. Sturgeons consist of 27 species in the world. Mitogenome data could be used to infer genetic diversity and investigate the evolutionary history of sturgeons. A limited number of complete mitogenomes in this family were sequenced. Here, we annotated the mitochondrial Huso huso genome, which revealed new aspects of this species.

RESULTS

In this species, the mitochondrial genome consisted of 13 genes encoding proteins, 22tRNA and 2rRNA, and two non-coding regions that followed other vertebrates. In addition, H. huso had a pseudo-tRNA-Glu between ND6 and Cytb and a 52-nucleotide tandem repeat with two replications in 12S rRNA. This duplication event is probably related to the slipped strand during replication, which could remain in the strand due to mispairing during replication. Furthermore, an 82 bp repeat sequence with three replications was observed in the D-loop control region, which is usually visible in different species. Regulatory elements were also seen in the control region of the mitochondrial genome, which included termination sequences and conserved regulatory blocks. Genomic compounds showed the highest conservation in rRNA and tRNA, while protein-encoded genes and nonencoded regions had the highest divergence. The mitochondrial genome was phylogenetically assayed using 12 protein-encoding genes.

CONCLUSIONS

In H. huso sequencing, we identified a distinct genome organization relative to other species that have never been reported. In recent years, along with the advancement in sequencing identified more genome rearrangements. However, it is an essential aspect of researching the evolution of the mitochondrial genome that needs to be recognized.

The First Complete Mitochondrial Genome of the Genus Pachycondyla (Formicidae, Ponerinae) and Insights into the Phylogeny of Ants

Ants are the standout group among eusocial insects in terms of their exceptional species richness and ecological dominance. The phylogenetic relationships among the group remain elusive. Mitochondrial genome sequences, as a kind of molecular marker, have been widely utilized in the phylogenetic analysis of insects. However, the number of ant mitogenomes published is still very limited. In this study, we utilized next-generation sequencing to determine the complete mitogenome of Pachycondyla annamita (Formicidae, Ponerinae). This is the first mitogenome from the genus Pachycondyla. Two gene rearrangements were identified in the mitogenome, the transposition of trnQ and trnM and the transposition of trnV and rrnS. The secondary structures of tRNAs were predicted. The tRNA genes trnR and trnS1 lacked the dihydrouridine (DHU) arm, and the trnE lacked the TΨC (T) arm. Phylogenetic analyses of the mitochondrial protein-coding genes under maximum likelihood (ML) and Bayesian inference (BI) criteria resulted in conflicting hypotheses. BI analysis using amino acid data with the site-heterogeneous mixture model produced a tree topology congruent with previous studies. The Formicidae was subdivided into two main clades, namely the "poneroid" clade and the "formicoid" clade. A sister group relationship between Myrmicinae and Formicinae was recovered within the "formicoid" clade.

Molecular Phylogeny and Evolution of the Tuerkayana (Decapoda: Brachyura: Gecarcinidae) Genus Based on Whole Mitochondrial Genome Sequences

Tuerkayana is of particular interest because it has been separated, in recent years, from Cardisoma and Discoplax but studies of its taxonomic status, especially from a whole mitochondrial genome perspective, have been lacking. In this study, the mitogenomes of four species (Tuerkayana magnum, Tuerkayana rotundum, Tuerkayana hirtipes, and Tuerkayana celeste) of Tuerkayana are sequenced and contrasted with other species in Brachyura for the first time. The phylogenetic tree of Brachyura, which includes 206 crab species (189 species of Brachyuran and 17 Anomura species) with a complete mitogenome, was constructed to evaluate the phylogenetic position of Tuerkayana and Gecarcinidae within Brachyuran, and explore the monophyly of Gecarcinidae. Furthermore, two single gene trees based on cox1 and 16SrRNA separately within interspecies of Gecarcinidae were reconstructed, providing molecular evidence for Tuerkayana and further clarifying the division of genera in Gecarcinidae. Based on the mitogenome dataset of 206 crabs, the branch-site model was utilized to explore selective pressure in individual codons with CodeML. The strong selective pressure shown in nad6 indicates that it may have played a significant role in the evolution of Gecarcinidae.

Comparative mitochondrial genomics of American nectar-feeding long-nosed bats Leptonycteris spp. with insights into the phylogeny of the family Phyllostomidae

Among leaf-nosed bats (family Phyllostomidae), the genus Leptonycteris (subfamily Glossophaginae), contains three migratory and obligate nectar-dwelling species of great ecological and economic importance; the greater long-nosed bat L. nivalis, the lesser long-nosed bat L. yerbabuenae, and the southern long-nosed bat L. curasoae. According to the IUCN, the three species are categorized as 'vulnerable', 'endangered', and 'near threatened', respectively. In this study, we assembled and characterized in detail the mitochondrial genome of Leptonycteris spp. and examined the phylogenetic position of this genus in the family Phyllostomidae based on protein coding genes (PCGs). The mitogenomes of L. nivalis, L. curasoae, and L. yerbabuenae are 16,708, 16,758, and 16,729 bp in length and each encode 13 PCGs, 22 transfer RNA genes, 2 rRNA genes, and a putative control region (CR). Mitochondrial gene order is identical to that reported before for the family Phyllostomidae. All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that is missing the DHU arm in the three species. All PCGs are exposed to purifying selection with atp8 experiencing the most relaxed purifying selection as the ω ratio was higher for this gene compared to the other PCGs in each species. The CR of each species contains three functional domains: extended termination associated sequence (ETAS), Central, and conserved sequence block domain (CSB). A phylomitogenomic analysis revealed that Leptonycteris is monophyletic and most closely related to the genus Glossophaga. The analysis also supported the monophyly of the family Glossophaginae in the speciose family Phyllostomidae. The mitochondria characterization of these species provides relevant information to develop molecular markers for conservation purposes.

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 Complete Mitochondrial Genome of the Hermit Crab Diogenes edwardsii (Anomura: Diogenidae) and Phylogenetic Relationships within Infraorder Anomura

A complete mitochondrial genome (mitogenome) can provide important information for gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few mitogenomes of hermit crabs (superfamily Paguridae) in the infraorder Anomura have been reported. This study reports the first complete mitogenome of the hermit crab Diogenes edwardsii assembled using high-throughput sequencing. The mitogenome of Diogenes edwardsii is 19,858 bp in length and comprises 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. There are 28 and six genes observed on the heavy and light strands, respectively. The genome composition was highly A + T biased (72.16%), and exhibited a negative AT-skew (-0.110) and positive GC-skew (0.233). Phylogenetic analyses based on the nucleotide dataset of 16 Anomura species indicated that D. edwardsii was closest related to Clibanarius infraspinatus in the same family, Diogenidae. Positive selection analysis showed that two residues located in cox1 and cox2 were identified as positively selected sites with high BEB value (>95%), indicating that these two genes are under positive selection pressure. This is the first complete mitogenome of the genus Diogenes, and this finding helps us to represent a new genomic resource for hermit crab species and provide data for further evolutionary status of Diogenidae in Anomura.

Self-Assembled Peptidyl Aggregates for the Fluorogenic Recognition of Mitochondrial DNA G-Quadruplexes

The selective monitoring of G-quadruplex (G4) structures in living cells is important to elucidate their functions and reveal their value as diagnostic or therapeutic targets. Here we report a fluorogenic probe (CV2) able to selectively light-up parallel G4 DNA over antiparallel topologies. CV2 was constructed by conjugating the excimer-forming CV dye with a peptide sequence (l-Arg-l-Gly-glutaric acid) that specifically recognizes G4s. CV2 forms self-assembled, red excimer-emitting nanoaggregates in aqueous media, but specific binding to G4s triggers its disassembly into rigidified monomeric dyes, leading to a dramatic fluorescence enhancement. Moreover, selective permeation of CV2 stains G4s in mitochondria over the nucleus. CV2 was employed for tracking the folding and unfolding of G4s in living cells, and for monitoring mitochondrial DNA (mtDNA) damage. These properties make CV2 appealing to investigate the possible roles of mtDNA G4s in diseases that involve mitochondrial dysfunction.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Characterization of mitochondrial genome of Indian Ocean blue-spotted maskray, Neotrygon indica and its phylogenetic relationship within Dasyatidae Family

The present study characterized complete mitochondrial genome of Blue-spotted maskray, Neotrygon indica and studied the evolutionary relationship of the species within the Dasyatidae family. The total length of the mitogenome was 17,974 bp including 37 genes and a non-coding control region. The average frequency of nucleotides in protein-coding genes was A: 29.1 %, T: 30.2 %, G: 13.0 % and C: 27.7 % with AT content of 59.3 %. The values of AT and GC skewness were -0.018 and -0.338, respectively. Comparative analyses showed a large number of average synonymous substitutions per synonymous site (Ks) in gene NADH4 (5.07) followed by NADH5 (4.72). High values of average number of non-synonymous substitutions per non-synonymous site (Ka) were observed in genes ATPase8 (0.54) and NADH2 (0.44). Genes NADH4L and NADH2 showed high interspecific genetic distance values of 0.224 ± 0.001 and 0.213 ± 0.002, respectively. Heat map analysis showed variation in codon usage among different species of the Dasyatidae family. The phylogenetic tree showed a sister relationship between the Dasyatinae and the Neotrygoninae subfamilies. Neotrygon indica formed as a sister species to the clade consisting of N. varidens and N. orientalis. Based on the present results, Neotrygon indica could have diverged from the common ancestor of the two latter in the Plio-Pleistocene. The present study showed distinct characteristics of N. indica from its congeners through comparative mitogenomics.

Aonchotheca (Nematoda: Capillariidae) is validated as a separated genus from Capillaria by both mitochondrial and nuclear ribosomal DNA

BACKGROUND

The family Capillariidae is a group of thread-like nematodes of 27 genera and over 300 species that infect a great variety of hosts including humans. Among these, some taxa such as the genus Aonchotheca have remained controversial regarding their systematic status for decades. The aim of the current study was to verify Aonchotheca's systemic status and to further determine whether it is a distinct genus from Capillaria using molecular and phylogenetic analyses.

RESULTS

We sequenced the mitochondrial (mt) genome and nuclear small subunit (18S) rRNA gene of Aonchotheca putorii, a representative species of the genus, and investigated its systematic status in Trichinellida using maximum likelihood and Bayesian inference. The differences in amino acid sequences of 13 protein-coding genes were 12.69-67.35% among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillaria with cox1 (12.69%) and atp8 (67.35%) as the most and the least conserved gene, respectively, and the difference of two mt rRNAs was 18.61-34.15%. Phylogenetic analyses of the complete mt genome and 18S rRNAs unequivocally showed that Aonchotheca was a distinct genus from Capillaria.

CONCLUSIONS

Large difference exists among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillarias. Aonchotheca putorii is the first species in the genus Aonchotheca for which a complete mitogenome has been sequenced. These data are useful for phylogenetics, systematics and the evolution of Capillariidae.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

The Complete Mitogenome of Toxocara vitulorum: Novel In-Sights into the Phylogenetics in Toxocaridae

Toxocara vitulorum (Ascaridida: Nematoda) is one of the most common intestinal nematodes of cattle and buffalos and, therefore, represents a serious threat to their populations worldwide. Despite its significance in veterinary health the epidemiology, population genetics, and molecular ecology of this nematode remain poorly understood. The mitogenome can yield a foundation for studying these areas and assist in the surveillance and control of T. vitulorum. Herein, the first whole mitogenome of T. vitulorum was sequenced utilizing Illumina technology and characterized with bioinformatic pipeline analyses. The entire genome of T. vitulorum was 15,045 bp in length and contained 12 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs). The gene arrangement (GA) of T. vitulorum was similar to those of other Toxocara species under GA3. The whole genome showed significant levels of AT and GC skew. Comparative mitogenomics including sequence identities, Ka/Ks, and sliding window analysis, indicated a purifying selection of 12 PCGs with cox1 and nad6 having the lowest and highest evolutionary rate, respectively. Whole amino acid sequence-based phylogenetic analysis supported a novel sister-species relationship of T. vitulorum with the congeneric species Toxocara canis, Toxocara cati, and Toxocara malaysiensis in the family Toxocaridae. Further, 12 (PCGs) single gene-based phylogenies suggested that nad4 and nad6 genes shared same topological trees with that of the whole genome, suggesting that these genes were suitable as novel genetic markers for phylogenetic and evolutionary studies of Ascaridida species. This complete mitogenome of T. vitulorum refined phylogenetic relationships in Toxocaridae and provided the resource of markers for population genetics, systematics, and epidemiology of this bovine nematode.

Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae

Apodemus is the most common small rodent species in the Palearctic realm and an ideal species for biogeographical research and understanding environmental changes. Elucidating phylogenetic relationships will help us better understand species adaptation and genetic evolution. Due to its stable structure, maternal inheritance, and rapid evolution, the mitogenome has become a hot spot for taxonomic and evolutionary studies. In this research, we determined the mitochondrial genome of Apodemus agrarius ningpoensis and Apodemus draco draco and studied the phylogeny of Muridae using ML and BI trees based on all known complete mitogenomes. The mitochondrial genome of Apodemus agrarius ningpoensis was 16,262 bp, whereas that of Apodemus draco draco was 16,222 bp, and both encoded 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of base composition showed a clear A-T preference. All tRNAs except tRNASer and tRNALys formed a typical trilobal structure. All protein-coding genes contained T- and TAA as stop codons. Phylogeny analysis revealed two main branches in the Muridae family. Apodemus agrarius ningpoensis formed sister species with Apodemus chevrieri, whereas Apodemus draco draco with Apodemus latronum. Our findings provide theoretical basis for future studies focusing on the mitogenome evolution of Apodemus.

Evolutionary inference across eukaryotes identifies universal features shaping organelle gene retention

Mitochondria and plastids power complex life. Why some genes and not others are retained in their organelle DNA (oDNA) genomes remains a debated question. Here, we attempt to identify the properties of genes and associated underlying mechanisms that determine oDNA retention. We harness over 15k oDNA sequences and over 300 whole genome sequences across eukaryotes with tools from structural biology, bioinformatics, machine learning, and Bayesian model selection. Previously hypothesized features, including the hydrophobicity of a protein product, and less well-known features, including binding energy centrality within a protein complex, predict oDNA retention across eukaryotes, with additional influences of nucleic acid and amino acid biochemistry. Notably, the same features predict retention in both organelles, and retention models learned from one organelle type quantitatively predict retention in the other, supporting the universality of these features-which also distinguish gene profiles in more recent, independent endosymbiotic relationships. A record of this paper's transparent peer review process is included in the supplemental information.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Complete mitochondrial genome of Idea leuconoe (Lepidoptera: Danaidae) and related phylogenetic analyses

In the present study, we first sequenced and determined the complete mitochondrial genome (mitogenome) of the giant Danaidae butterfly, Idea leuconoe (Lepidoptera: Danaidae). The mitogenome was a typical closed, circular, double-stranded DNA molecule of 15,278 bp length (GenBank accession number: KR815449), similar to the metazoan mitogenomes containing 37 genes and one A + T-rich region. All the protein-coding genes (PCGs) were initiated with a typical ATN codon. Seven genes (COII, ATP6, COIII, nad4, nad4L, cytb, and nad1) adopted the standard ATG start codon, but the remaining six genes were initiated with ATA. All the 13 PCGs harbored complete termination codons (TAA). The overlap nucleotides ATGATAA were conserved for the ATP8/ATP6 gene. The largest intergenic spacer was located between trnGln and nad2, a common finding in Lepidoptera butterflies. All the transfer RNA genes in the I. leuconoe mitogenome could be folded into typical clover-leaf secondary structures, except for trnSer (AGN) that lacked a dihydrouridine arm. The control region with 94.8% A + T content was 444 bp in length and located between rrnS and trnMet. Finally, the phylogenetic relationships obtained using the maximum likelihood and Bayesian methods confirmed two well-supported phylogenetic trees of Danaidae, Papilionidae, and Nymphalidae from the order Lepidoptera, which were consistent with the traditional morphological classification. Results provided additional information for butterfly phylogenetic analysis and insights into the evolution of genomes.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

The Complete Mitochondrial Genome of Stichopus naso (Aspidochirotida: Stichopodidae: Stichopus) and Its Phylogenetic Position

The mitochondrial genome is widely used to study the molecular evolution of and perform phylogenetic analyses on animals. In this study, the complete mitochondrial genome (mitogenome) of Stichopus naso was sequenced. The mitogenome was 16,239 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The genome composition showed positive AT-skew (0.023) and negative GC-skew (−0.158). The order of the mitochondrial genes was consistent with those from the Stichopus and Isostichopus species, whereas it was different from those of other species of Aspidochirotida. The phylogenetic analysis, based on the nucleotide sequences of 13 PCGs through the methods of Bayesian inference (BI) and maximum likelihood (ML), indicated that S. naso has close relationships with S. horrens and S. monotuberculatus, and belongs to a member of Stichopodidae. Our study provides a reference mitogenome for further molecular evolution studies and phylogenetic research on sea cucumbers.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

Comparative Mitogenome Analysis of Gerbils and the Mitogenome Phylogeny of Gerbillinae (Rodentia: Muridae)

To enrich the mitogenomic database of Gerbillinae (Rodentia: Muridae), mitogenomes of three gerbils from different genera, Meriones tamariscinus (16,393 bp), Brachiones przewalskii (16,357 bp), and Rhombomys opimus (16,352 bp), were elaborated and compared with those of other gerbils in the present study. The three gerbil mitogenomes consisted of 2 ribosomal RNA genes, 13 protein-coding genes (PCGs), 22 transfer RNA genes, and one control region. Here, gerbil mitogenomes have shown unique characteristics in terms of base composition, codon usage, non-coding region, and the replication origin of the light strand. There was no significant correlation between the nucleotide percentage of G + C and the phylogenetic status in gerbils, and between the GC content of PCGs and the leucine count. Phylogenetic relationships of the subfamily Gerbillinae were reconstructed by 7 gerbils that represented four genera based on concatenated mitochondrial DNA data using both Bayesian Inference and Maximum Likelihood. The phylogenetic analysis indicated that M. tamariscinus was phylogenetically distant from the genus Meriones, but has a close relationship with R. opimus. B. przewalskii was closely related to the genus Meriones rather than that of R. opimus.

Comparative phylogeography among eight Neotropical wild cat species: no single evolutionary pattern

The felid species of South America are thought to have arrived on the continent during the Great American Biotic Interchange (GABI) in the Pleistocene. However, molecular and palaeontological data do not agree on how this event affected speciation in felids. Here, we determine both the number of colonization events and the period when felines first migrated from North America to South America. In addition, we evaluate whether similar evolutionary events could have affected the eight Neotropical cat species in their levels of genetic diversity, spatial genetic structure and demographic changes. We analysed four concatenated mitochondrial genes of the jaguar, ocelot, margay, tigrina, pampas cat, Andean cat, puma and jaguarundi. The samples were representative of a wide distribution of these species in Central and South America. Our analysis suggests either three or four colonization events from North America to South America over the past 3 Myr, followed by subsequent speciation events and the attainment of high or very high genetic diversity levels for seven of the species. No unique evolutionary process was detected for any of the current Neotropical cat species.

Mitochondrial DNA Mutagenesis: Feature of and Biomarker for Environmental Exposures and Aging

PURPOSE OF REVIEW

Mitochondrial dysfunction is a hallmark of aging. Mitochondrial genome (mtDNA) instability contributes to mitochondrial dysfunction, and mtDNA mutagenesis may contribute to aging. However, the origin of mtDNA mutations remains somewhat controversial. The goals of this review are to introduce and review recent literature on mtDNA mutagenesis and aging, address recent animal and epidemiological evidence for the effects of chemicals on mtDNA damage and mutagenesis, propose hypotheses regarding the contribution of environmental toxicant exposure to mtDNA mutagenesis in the context of aging, and suggest future directions and approaches for environmental health researchers.

RECENT FINDINGS

Stressors such as pollutants, pharmaceuticals, and ultraviolet radiation can damage the mitochondrial genome or disrupt mtDNA replication, repair, and organelle homeostatic processes, potentially influencing the rate of accumulation of mtDNA mutations. Accelerated mtDNA mutagenesis could contribute to aging, diseases of aging, and sensitize individuals with pathogenic mtDNA variants to stressors. We propose three potential mechanisms of toxicant-induced effects on mtDNA mutagenesis over lifespan: (1) increased de novo mtDNA mutations, (2) altered frequencies of mtDNA mutations, or (3) both. There are remarkably few studies that have investigated the impact of environmental chemical exposures on mtDNA instability and mutagenesis, and even fewer in the context of aging. More studies are warranted because people are exposed to tens of thousands of chemicals, and are living longer. Finally, we suggest that toxicant-induced mtDNA damage and mutational signatures may be a sensitive biomarker for some exposures.

The mtDNA mutation spectrum in the PolG mutator mouse reveals germline and somatic selection

BACKGROUND

Mitochondrial DNA (mtDNA) codes for products necessary for electron transport and mitochondrial gene translation. mtDNA mutations can lead to human disease and influence organismal fitness. The PolG mutator mouse lacks mtDNA proofreading function and rapidly accumulates mtDNA mutations, making it a model for examining the causes and consequences of mitochondrial mutations. Premature aging in PolG mice and their physiology have been examined in depth, but the location, frequency, and diversity of their mtDNA mutations remain understudied. Identifying the locations and spectra of mtDNA mutations in PolG mice can shed light on how selection shapes mtDNA, both within and across organisms.

RESULTS

Here, we characterized somatic and germline mtDNA mutations in brain and liver tissue of PolG mice to quantify mutation count (number of unique mutations) and frequency (mutation prevalence). Overall, mtDNA mutation count and frequency were the lowest in the D-loop, where an mtDNA origin of replication is located, but otherwise uniform across the mitochondrial genome. Somatic mtDNA mutations have a higher mutation count than germline mutations. However, germline mutations maintain a higher frequency and were also more likely to be silent. Cytosine to thymine mutations characteristic of replication errors were the plurality of basepair changes, and missense C to T mutations primarily resulted in increased protein hydrophobicity. Unlike wild type mice, PolG mice do not appear to show strand asymmetry in mtDNA mutations. Indel mutations had a lower count and frequency than point mutations and tended to be short, frameshift deletions.

CONCLUSIONS

Our results provide strong evidence that purifying selection plays a major role in the mtDNA of PolG mice. Missense mutations were less likely to be passed down in the germline, and they were less likely to spread to high frequencies. The D-loop appears to have resistance to mutations, either through selection or as a by-product of replication processes. Missense mutations that decrease hydrophobicity also tend to be selected against, reflecting the membrane-bound nature of mtDNA-encoded proteins. The abundance of mutations from polymerase errors compared with reactive oxygen species (ROS) damage supports previous studies suggesting ROS plays a minimal role in exacerbating the PolG phenotype, but our findings on strand asymmetry provide discussion for the role of polymerase errors in wild type organisms. Our results provide further insight on how selection shapes mtDNA mutations and on the aging mechanisms in PolG mice.

A genomic portrait of Sturnira parvidens: mitochondrial chromosome, repetitive elements, and microsatellite discovery

The yellow-shouldered bat, Sturnira parvidens (Phyllostomidae), is an abundant and widespread species in southern North America and Mesoamerica. Despite its important ecological role, no genomic resources exist for this species. Using low-coverage short Illumina 150 bp pair-end reads sequencing, this study reports the mitochondrial chromosome and nuclear repetitive elements, including microsatellites, in S. parvidens. The mitochondrial genome of S. parvidens is 16,612 bp in length and is comprised of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Repetitive elements constituted ~67% of the nuclear genome while ~33% of the genome represented single- or low-copy sequences. A moderate proportion of repetitive sequences (31% putative families) could not be assigned to known repeat element families. Considering only annotated repetitive elements, the most ubiquitous repetitive elements belonged to Class I-LINE and Satellite DNA, which were considerably more abundant than Class I-LTR elements and Class II-DNA transposons (TcMar-Mariner and hAT-Charlie). A total of 193 microsatellites were identified.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

The population frequency of human mitochondrial DNA variants is highly dependent upon mutational bias

Next-generation sequencing can quickly reveal genetic variation potentially linked to heritable disease. As databases encompassing human variation continue to expand, rare variants have been of high interest, since the frequency of a variant is expected to be low if the genetic change leads to a loss of fitness or fecundity. However, the use of variant frequency when seeking genomic changes linked to disease remains very challenging. Here, I explored the role of selection in controlling human variant frequency using the HelixMT database, which encompasses hundreds of thousands of mitochondrial DNA (mtDNA) samples. I found that a substantial number of synonymous substitutions, which have no effect on protein sequence, were never encountered in this large study, while many other synonymous changes are found at very low frequencies. Further analyses of human and mammalian mtDNA datasets indicate that the population frequency of synonymous variants is predominantly determined by mutational biases rather than by strong selection acting upon nucleotide choice. My work has important implications that extend to the interpretation of variant frequency for non-synonymous substitutions.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

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.

The complete mitochondrial genome of capillariid nematodes (Eucoleus annulatus): A novel gene arrangement and phylogenetic implications

Capillariid nematode is a group of endoparasites of vertebrates with a complex taxonomy, causing significant economic losses to poultry industry. The taxonomic status of the genus Eucoleus remained controversial for several decades. Mitochondrial (mt) DNA provides useful genetic markers for accurate identification of species, but complete mt genome sequences have been lacking for any Capillariid nematodes. In the present study, we decoded the complete mt genome of E. annulatus and examined its phylogenetic relationship with selected members of the class Enoplea nematodes. The circular mt genome of E. annulatus was 14,118 bp, encoded 37 genes with a single non-coding region and showed substantial gene rearrangements (especially tRNA genes) compared to other nematodes studied to date. The complete mt genome of E. annulatus showed a clear A + T bias in nucleotide composition. The number of A (5404) was approximately equal to T (5405) and the GC-skew was negative on average (-0.073). Phylogenetic analyses based on 18S rDNA placed Eucoleus spp. well apart from each other and supported the proposal that Eucoleus and Capillaria are two distinct genera. Similarly, Bayesian inference (BI) and Maximum likelihood (ML) phylogenies based on mtDNA sequences revealed that the family Capillariidae is more closely related to the family Trichuridae than to the family Trichinellidae. This is the first report of the complete mt genome of capillariid nematodes, and it will provide additional genetic markers for studying the molecular epidemiology, population genetics and systematics of capillariid nematodes and should have implications for the molecular diagnosis, prevention, and control of capillariosis in animals.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Potential causes and consequences of rapid mitochondrial genome evolution in thermoacidophilic Galdieria (Rhodophyta)

BACKGROUND

The Cyanidiophyceae is an early-diverged red algal class that thrives in extreme conditions around acidic hot springs. Although this lineage has been highlighted as a model for understanding the biology of extremophilic eukaryotes, little is known about the molecular evolution of their mitochondrial genomes (mitogenomes).

RESULTS

To fill this knowledge gap, we sequenced five mitogenomes from representative clades of Cyanidiophyceae and identified two major groups, here referred to as Galdieria-type (G-type) and Cyanidium-type (C-type). G-type mitogenomes exhibit the following three features: (i) reduction in genome size and gene inventory, (ii) evolution of unique protein properties including charge, hydropathy, stability, amino acid composition, and protein size, and (iii) distinctive GC-content and skewness of nucleotides. Based on GC-skew-associated characteristics, we postulate that unidirectional DNA replication may have resulted in the rapid evolution of G-type mitogenomes.

CONCLUSIONS

The high divergence of G-type mitogenomes was likely driven by natural selection in the multiple extreme environments that Galdieria species inhabit combined with their highly flexible heterotrophic metabolism. We speculate that the interplay between mitogenome divergence and adaptation may help explain the dominance of Galdieria species in diverse extreme habitats.

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.