Functional conservatism in mitochondrial evolution: insight from hybridization of arctic and brook charrs

Phylogeography using mitogenomes: A rare Dipodidae, Sicista betulina , in North‐western Europe

Repeated climatic and vegetation changes during the Pleistocene have shaped biodiversity in Northern Europe including Denmark. The Northern Birch Mouse (Sicista betulina) was one of the first small rodent species to colonize Denmark after the Late Glacial Maximum. This study analyses complete mitochondrial genomes and two nuclear genes of the Northern Birch Mouse to investigate the phylogeographical pattern in North‐western Europe and test whether the species colonized Denmark through several colonization events. The latter was prompt by (i) the present‐day distinct northern and southern Danish distribution and (ii) the subfossil record of Northern Birch Mouse, supporting early Weichselian colonization. Samples from Denmark, Norway, Sweden, Russia, Latvia, Estonia, and Slovakia were included. Mitogenomes were obtained from 54 individuals, all representing unique mitogenomes supporting high genetic variation. Bayesian phylogenetic analysis identified two distinct evolutionary linages in Northern Europe diverging within the Elster glaciation period. The results of the two nuclear genomes showed lower genetic differentiation but supported the same evolutionary history. This suggests an allopatric origin of the clades followed by secondary contact. Individuals from southern Denmark were only found in one clade, while individuals from other areas, including northern Denmark, were represented in both clades. Nevertheless, we found no evidence for repeated colonization's explaining the observed fragmented distribution of the species today. The results indicated that the mitogenome pattern of the Northern Birch Mouse population in southern Denmark was either (i) due to the population being founded from northern Denmark, (ii) a result of climatic and anthropogenic effects reducing population size increasing genetic drift or (iii) caused by sampling bias.

Gene flow prevents mitonuclear co-adaptation: A comparative portrait of sympatric wild types and cybrids in the fish Chrosomus eos

Allospecific mtDNA can occasionally be beneficial for the fitness of populations. It is, however, difficult to assess the effect of mtDNA in natural conditions due to genetic and/or environmental interactions. In the fish Chrosomus eos, the transfer of C. neogaeus mitochondria occurs in a single generation and results in natural cybrids. For a few lakes in Quebec, C. eos can harbor either a C. eos mtDNA (wild types) or a C. neogaeus mtDNA (cybrids). Moreover, mtDNA of cybrids originated either from Mississippian or Atlantic glacial refuges. Such diversity provides a useful system for in situ assessment of allospecific mtDNA effects. We determined genetic, epigenetic and transcriptomic variation as well as mitochondrial enzymatic activity (complex IV) changes among wild types and cybrids either in sympatry or allopatry. Wild types and cybrids did not segregate spatially within a lake. Moreover, no significant genetic differentiation was detected among wild types and cybrids indicating sustained gene flow. Mitochondrial complex IV activity was higher for cybrids in both sympatry and allopatry while no difference was detected among cybrid haplotypes. Epigenetic and transcriptomic analyses revealed only subtle differences between sympatric wild types and cybrids compared to differences between sites. Altogether, these results indicate a limited influence of allospecific mtDNA in nuclear gene expression when controlling for genetic and environmental effects. The absence of a reproductive barrier between wild types and cybrids results in random association of either C. eos or C. neogaeus mtDNA with C. eos nDNA at each generation, and prevents mitonuclear co-adaptation in sympatry.

The effects of allospecific mitochondrial genome on the fitness of northern redbelly dace (Chrosomus eos)

Instantaneous mitochondrial introgression events allow the disentangling of the effects of hybridization from those of allospecific mtDNA. Such process frequently occurred in the fish Chrosomus eos, resulting in cybrid individuals composed of a C. eos nuclear genome but with a C. neogaeus mtDNA. This provides a valuable model to address the fundamental question: How well do introgressed individuals perform in their native environment? We infer where de novo production of cybrids occurred to discriminate native environments from those colonized by cybrids in 25 sites from two regions (West-Qc and East-Qc) in Quebec (Canada). We then compared the relative abundance of wild types and cybrids as a measure integrating both fitness and de novo production of cybrids. According to mtDNA variation, 12 introgression events are required to explain the diversity of cybrids. Five cybrid lineages could not be associated with in situ introgression events. This includes one haplotype carried by 93% of the cybrids expected to have colonized West-Qc. These cybrids also displayed a nearly complete allopatric distribution with wild types. We still inferred de novo production of cybrids at seven sites, that accounted for 70% of the cybrids in East-Qc. Wild-type and cybrid individuals coexist in all East-Qc sites while cybrids were less abundant. Allopatry of cybrids restricted to the postglacial expansion suggests the existence of higher fitness for cybrids in specific conditions, allowing for the colonization of different environments and expanding the species' range. However, allospecific mtDNA does not provide a higher fitness to cybrids in their native environment compared to wild types, making the success of an introgressed lineage uncertain.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

The on-again, off-again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein-coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species-specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein-coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four-fold higher dN /dS ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well-preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.

The elusive nature of adaptive mitochondrial DNA evolution of an arctic lineage prone to frequent introgression

Mitochondria play a fundamental role in cellular metabolism, being responsible for most of the energy production of the cell in the oxidative phosphorylation (OXPHOS) pathway. Mitochondrial DNA (mtDNA) encodes for key components of this process, but its direct role in adaptation remains far from understood. Hares (Lepus spp.) are privileged models to study the impact of natural selection on mitogenomic evolution because 1) species are adapted to contrasting environments, including arctic, with different metabolic pressures, and 2) mtDNA introgression from arctic into temperate species is widespread. Here, we analyzed the sequences of 11 complete mitogenomes (ten newly obtained) of hares of temperate and arctic origins (including two of arctic origin introgressed into temperate species). The analysis of patterns of codon substitutions along the reconstructed phylogeny showed evidence for positive selection in several codons in genes of the OXPHOS complexes, most notably affecting the arctic lineage. However, using theoretical models, no predictable effect of these differences was found on the structure and physicochemical properties of the encoded proteins, suggesting that the focus of selection may lie on complex interactions with nuclear encoded peptides. Also, a cloverleaf structure was detected in the control region only from the arctic mtDNA lineage, which may influence mtDNA replication and transcription. These results suggest that adaptation impacted the evolution of hare mtDNA and may have influenced the occurrence and consequences of the many reported cases of massive mtDNA introgression. However, the origin of adaptation remains elusive.

The complete mitochondrial genome of the land snail Cornu aspersum (Helicidae: Mollusca): intra-specific divergence of protein-coding genes and phylogenetic considerations within Euthyneura

The complete sequences of three mitochondrial genomes from the land snail Cornu aspersum were determined. The mitogenome has a length of 14050 bp, and it encodes 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. It also includes nine small intergene spacers, and a large AT-rich intergenic spacer. The intra-specific divergence analysis revealed that COX1 has the lower genetic differentiation, while the most divergent genes were NADH1, NADH3 and NADH4. With the exception of Euhadra herklotsi, the structural comparisons showed the same gene order within the family Helicidae, and nearly identical gene organization to that found in order Pulmonata. Phylogenetic reconstruction recovered Basommatophora as polyphyletic group, whereas Eupulmonata and Pulmonata as paraphyletic groups. Bayesian and Maximum Likelihood analyses showed that C. aspersum is a close relative of Cepaea nemoralis, and with the other Helicidae species form a sister group of Albinaria caerulea, supporting the monophyly of the Stylommatophora clade.

Methods for assessing mitochondrial function in diabetes

A growing body of research is investigating the potential contribution of mitochondrial function to the etiology of type 2 diabetes. Numerous in vitro, in situ, and in vivo methodologies are available to examine various aspects of mitochondrial function, each requiring an understanding of their principles, advantages, and limitations. This review provides investigators with a critical overview of the strengths, limitations and critical experimental parameters to consider when selecting and conducting studies on mitochondrial function. In vitro (isolated mitochondria) and in situ (permeabilized cells/tissue) approaches provide direct access to the mitochondria, allowing for study of mitochondrial bioenergetics and redox function under defined substrate conditions. Several experimental parameters must be tightly controlled, including assay media, temperature, oxygen concentration, and in the case of permeabilized skeletal muscle, the contractile state of the fibers. Recently developed technology now offers the opportunity to measure oxygen consumption in intact cultured cells. Magnetic resonance spectroscopy provides the most direct way of assessing mitochondrial function in vivo with interpretations based on specific modeling approaches. The continuing rapid evolution of these technologies offers new and exciting opportunities for deciphering the potential role of mitochondrial function in the etiology and treatment of diabetes.

Mitochondrial genome organization and divergence in hybridizing central European waterfrogs of the Pelophylax esculentus complex (Anura, Ranidae)

Natural transfer of mitochondrial DNA has occurred between three western Palaearctic waterfrog taxa: Pelophylax lessonae, Pelophylax ridibundus and their hybridogenetic hybrid, Pelophylax kl. esculentus. The transfer is asymmetric with most P. kl. esculentus and approximately one third of all central European P. ridibundus having mtDNA derived from P. lessonae (L-mtDNA). We obtained complete nucleotide sequences of multiple mitochondrial genomes (15,376-78 bp without control regions) from all 3 taxa, including a P. ridibundus frog with introgressed L-mtDNA. The gene content and organization of the mitogenomes correspond to those typical of neobatrachians. Divergence between the mtDNAs of P. lessonae and P. ridibundus is high with an uncorrected p-distance of 11.9% across the entire mitogenome. However, the rate of nucleotide substitution depends on the degree of functional constraint with up to 30-fold differences in levels of divergence. In general, mitochondrial genes encoding the translational machinery evolve very slowly, whereas genes encoding polypeptides of the electron transport system, especially the ND genes, evolve rapidly. Only 25 of 211-213 observed amino acid replacements could be classified as radical and are therefore more likely to be exposed to selection. A disproportionately high number of amino acid substitutions has occurred in the ND4, ND4L and cytb genes of the P. lessonae lineage (including 36% of all radical changes). In contrast to the interspecific divergence, nucleotide polymorphism within L- and R-mtDNA is very low: L-mtDNA haplotypes differed on average by only 19 nucleotides, while there was no variation within two mtDNAs derived from P. ridibundus. This is an expected finding considering that we have sampled a post-glacial expansion area. Moreover, the introgressed L-mtDNA on a P. ridibundus background differed from other L-mtDNAs by only a few substitutions, indicative of a very recent introgression event. We discuss our findings in the context of natural selection acting on L-mtDNA and its potential significance in cytonuclear epistasis.

Mitochondrial phylogenomics of the Bivalvia (Mollusca): searching for the origin and mitogenomic correlates of doubly uniparental inheritance of mtDNA

Background

Doubly uniparental inheritance (DUI) is an atypical system of animal mtDNA inheritance found only in some bivalves. Under DUI, maternally (F genome) and paternally (M genome) transmitted mtDNAs yield two distinct gender-associated mtDNA lineages. The oldest distinct M and F genomes are found in freshwater mussels (order Unionoida). Comparative analyses of unionoid mitochondrial genomes and a robust phylogenetic framework are necessary to elucidate the origin, function and molecular evolutionary consequences of DUI. Herein, F and M genomes from three unionoid species, Venustaconcha ellipsiformis, Pyganodon grandis and Quadrula quadrula have been sequenced. Comparative genomic analyses were carried out on these six genomes along with two F and one M unionoid genomes from GenBank (F and M genomes of Inversidens japanensis and F genome of Lampsilis ornata).

Results

Compared to their unionoid F counterparts, the M genomes contain some unique features including a novel localization of the trnH gene, an inversion of the atp8-trnD genes and a unique 3'coding extension of the cytochrome c oxidase subunit II gene. One or more of these unique M genome features could be causally associated with paternal transmission. Unionoid bivalves are characterized by extreme intraspecific sequence divergences between gender-associated mtDNAs with an average of 50% for V. ellipsiformis, 50% for I. japanensis, 51% for P. grandis and 52% for Q. quadrula (uncorrected amino acid p-distances). Phylogenetic analyses of 12 protein-coding genes from 29 bivalve and five outgroup mt genomes robustly indicate bivalve monophyly and the following branching order within the autolamellibranch bivalves: ((Pteriomorphia, Veneroida) Unionoida).

Conclusion

The basal nature of the Unionoida within the autolamellibranch bivalves and the previously hypothesized single origin of DUI suggest that (1) DUI arose in the ancestral autolamellibranch bivalve lineage and was subsequently lost in multiple descendant lineages and (2) the mitochondrial genome characteristics observed in unionoid bivalves could more closely resemble the DUI ancestral condition. Descriptions and comparisons presented in this paper are fundamental to a more complete understanding regarding the origins and consequences of DUI.

Comparative genomics of Drosophila mtDNA: Novel features of conservation and change across functional domains and lineages

To gain insight on mitochondrial DNA (mtDNA) evolution, we assembled and analyzed the mitochondrial genomes of Drosophila erecta, D. ananassae, D. persimilis, D. willistoni, D. mojavensis, D. virilis and D. grimshawi together with the sequenced mtDNAs of the melanogaster subgroup. Genomic comparisons across the well-defined Drosophila phylogeny impart power for detecting conserved mtDNA regions that maintain metabolic function and regions that evolve uniquely on lineages. Evolutionary rate varies across intergenic regions of the mtDNA. Rapidly evolving intergenic regions harbor the majority of mitochondrial indel divergence. In contrast, patterns of nearly perfect conservation within intergenic regions reveal a refined set of nucleotides underlying the binding of transcription termination factors. Sequencing of 5' cDNA ends indicates that cytochrome C oxidase I (CoI) has a novel (T/C)CG start codon and that perfectly conserved regions upstream of two NADH dehydrogenase (ND) genes are transcribed and likely extend these protein sequences. Substitutions at synonymous sites in the Drosophila mitochondrial proteomes reflect a mutation process that is biased toward A and T nucleotides and differs between mtDNA strands. Differences in codon usage bias across genes reveal that weak selection at silent sites may offset the mutation bias. The mutation-selection balance at synonymous sites has also diverged between the Drosophila and Sophophora lineages. Rates of evolution are highly heterogeneous across the mitochondrial proteome, with ND accumulating many more amino acid substitutions than CO. These oxidative phosphorylation complex-specific rates of evolution vary across lineages and may reflect physiological and ecological change across the Drosophila phylogeny.

Near neutrality: leading edge of the neutral theory of molecular evolution

The nearly neutral theory represents a development of Kimura's neutral theory of molecular evolution that makes testable predictions that go beyond a mere null model. Recent evidence has strongly supported several of these predictions, including the prediction that slightly deleterious variants will accumulate in a species that has undergone a severe bottleneck or in cases where recombination is reduced or absent. Because bottlenecks often occur in speciation and slightly deleterious mutations in coding regions will usually be nonsynonymous, we should expect that the ratio of nonsynonymous to synonymous fixed differences between species should often exceed the ratio of nonsynonymous to synonymous polymorphisms within species. Many data support this prediction, although they have often been wrongly interpreted as evidence for positive Darwinian selection. The use of conceptually flawed tests for positive selection has become widespread in recent years, seriously harming the quest for an understanding of genome evolution. When properly analyzed, many (probably most) claimed cases of positive selection will turn out to involve the fixation of slightly deleterious mutations by genetic drift in bottlenecked populations. Slightly deleterious variants are a transient feature of evolution in the long term, but they have substantially affected contemporary species, including our own.

The unusual system of doubly uniparental inheritance of mtDNA: isn't one enough?

Mitochondria possess their own genetic material (mitochondrial DNA or mtDNA), whose gene products are involved in mitochondrial respiration and oxidative phosphorylation, transcription, and translation. In animals, mitochondrial DNA is typically transmitted to offspring by the mother alone. The discovery of 'doubly uniparental inheritance' (DUI) of mtDNA in some bivalves has challenged the paradigm of strict maternal inheritance (SMI). In this review, we survey recent advances in our understanding of DUI, which is a peculiar system of cytoplasmic DNA inheritance that involves distinct maternal and paternal routes of mtDNA transmission, a novel extension of a mitochondrial gene (cox2), recombination, and periodic 'role-reversals' of the normally male and female-transmitted mitochondrial genomes. DUI provides a unique opportunity for studying nuclear-cytoplasmic genome interactions and the evolutionary significance of different modes of mitochondrial inheritance.