Quantifying the Elevation of Mitochondrial DNA Evolutionary Substitution Rates Over Nuclear Rates in the Intertidal Copepod Tigriopus californicus

Absence of heterosis for hypoxia tolerance in F1 hybrids of Tigriopus californicus

Hybridization produces a range of outcomes from advantageous to disadvantageous, and a goal of genetic research is to understand the gene interactions that generate these outcomes. Interactions between cytoplasmic elements, such as mitochondria, and the nucleus may be particularly vulnerable to accruing disadvantageous combinations as a result of their different rates of evolution. Consequently, mitonuclear incompatibilities may play an important role in hybrid outcomes even if their negative impacts could be masked for some fitness measures by heterosis in first-generation (F1) hybrids. We used Tigriopus californicus, a model system for mitonuclear incompatibilities that is also known for exhibiting heterosis in the F1 generation and outbreeding depression in later generations, to test whether heterosis or outbreeding depression would occur when mitonuclear mismatch was paired with a stress that heavily impacts mitochondrial processes-specifically, hypoxia. We generated 284 parental and 436 F1 hybrids from four population crosses (720 total) and compared parental and F1 populations for hypoxia tolerance. We observed that, on average, F1 hybrids were less likely to survive a hypoxia stress test than parental populations, although we did not detect a statistically significant trend (P = 0.246 to 0.614). This suggests that hypoxia may be a particularly intense stressor for mitonuclear coordination and hybridization outcomes vary by trait.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Recovery from hybrid breakdown reveals a complex genetic architecture of mitonuclear incompatibilities

Reproductive isolation is often achieved when genes that are neutral or beneficial in their genomic background become functionally incompatible in a foreign genomic background, causing inviability, sterility or other forms of low fitness in hybrids. Recent studies suggest that mitonuclear interactions are among the initial incompatibilities to evolve at early stages of population divergence across taxa. Yet, the genomic architecture of mitonuclear incompatibilities has rarely been elucidated. We employ an experimental evolution approach starting with low-fitness F₂ interpopulation hybrids of the copepod Tigriopus californicus, in which frequencies of compatible and incompatible nuclear alleles change in response to an alternative mitochondrial background. After about nine generations, we observe a generalized increase in population size and in survivorship, suggesting efficiency of selection against maladaptive phenotypes. Whole genome sequencing of evolved populations showed some consistent allele frequency changes across three replicates of each reciprocal cross, but markedly different patterns between mitochondrial backgrounds. In only a few regions (~6.5% of the genome), the same parental allele was overrepresented irrespective of the mitochondrial background. About 33% of the genome showed allele frequency changes consistent with divergent selection, with the location of these genomic regions strongly differing between mitochondrial backgrounds. In 87% and 89% of these genomic regions, the dominant nuclear allele matched the associated mitochondrial background, consistent with mitonuclear co-adaptation. These results suggest that mitonuclear incompatibilities have a complex polygenic architecture that differs between populations, potentially generating genome-wide barriers to gene flow between closely related taxa.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

Not all mitochondrial DNAs are made equal and the nucleus knows it

The oxidative phosphorylation (OXPHOS) system is the only structure in animal cells with components encoded by two genomes, maternally transmitted mitochondrial DNA (mtDNA), and biparentally transmitted nuclear DNA (nDNA). MtDNA‐encoded genes have to physically assemble with their counterparts encoded in the nucleus to build together the functional respiratory complexes. Therefore, structural and functional matching requirements between the protein subunits of these molecular complexes are rigorous. The crosstalk between nDNA and mtDNA needs to overcome some challenges, as the nuclear‐encoded factors have to be imported into the mitochondria in a correct quantity and match the high number of organelles and genomes per mitochondria that encode and synthesize their own components locally. The cell is able to sense the mito‐nuclear match through changes in the activity of the OXPHOS system, modulation of the mitochondrial biogenesis, or reactive oxygen species production. This implies that a complex signaling cascade should optimize OXPHOS performance to the cellular‐specific requirements, which will depend on cell type, environmental conditions, and life stage. Therefore, the mitochondria would function as a cellular metabolic information hub integrating critical information that would feedback the nucleus for it to respond accordingly. Here, we review the current understanding of the complex interaction between mtDNA and nDNA.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Assessing the fitness consequences of mitonuclear interactions in natural populations

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

First use of molecular evidence to match sexes in the Monstrilloida (Crustacea: Copepoda), and taxonomic implications of the newly recognized and described, partly Maemonstrilla-like females of Monstrillopsis longilobata Lee, Kim & Chang, 2016

Two forms of the monstrilloid copepod, one represented by males, the other by females, were dominant among the monstrilloids collected while sampling zooplankton with a light trap in southern coastal waters of Korea. Morphological examination revealed that the males are conspecific with the previously reported and rather specialized Korean species Monstrillopsis longilobata Lee, Kim & Chang, 2016, hitherto known only from males that have extremely long genital lappets. The females also show several diagnostic features of Monstrillopsis, such as prominent eyes, bilobed fifth legs with the inner lobe unarmed and reduced, and four urosomal somites. In addition, though, these females are extraordinary among all Monstrilloida in that their ovigerous spines are directed ventrally, not posteriorly as in most species or anteriorly as in those assigned to Maemonstrilla. Genetic divergence analyses based on partial mitochondrial cytochrome c oxidase subunit I, complete Internal Transcribed Spacer 1 (ITS1)–5.8S–ITS2, and partial 28S rRNA revealed little or no genetic divergence between the present males and females, thus demonstrating that they are mutually conspecific. The present report (1) provides the first morphological description of female M. longilobata with the proposal of a revised nomenclature for antennular setal armament; (2) presents the molecular evidence for conspecificity of the males and females; and (3) lists several morphological characteristics that are sexually dimorphic in this species, and thus likely also in other monstrilloids. Matters bearing on the validity of the genera Haemocera, Monstrillopsis, and Maemonstrilla are discussed.

Large Variation in the Ratio of Mitochondrial to Nuclear Mutation Rate across Animals: Implications for Genetic Diversity and the Use of Mitochondrial DNA as a Molecular Marker

It is commonly assumed that mitochondrial DNA (mtDNA) evolves at a faster rate than nuclear DNA (nuDNA) in animals. This has contributed to the popularity of mtDNA as a molecular marker in evolutionary studies. Analyzing 121 multilocus data sets and four phylogenomic data sets encompassing 4,676 species of animals, we demonstrate that the ratio of mitochondrial over nuclear mutation rate is highly variable among animal taxa. In nonvertebrates, such as insects and arachnids, the ratio of mtDNA over nuDNA mutation rate varies between 2 and 6, whereas it is above 20, on average, in vertebrates such as scaled reptiles and birds. Interestingly, this variation is sufficient to explain the previous report of a similar level of mitochondrial polymorphism, on average, between vertebrates and nonvertebrates, which was originally interpreted as reflecting the effect of pervasive positive selection. Our analysis rather indicates that the among-phyla homogeneity in within-species mtDNA diversity is due to a negative correlation between mtDNA per-generation mutation rate and effective population size, irrespective of the action of natural selection. Finally, we explore the variation in the absolute per-year mutation rate of both mtDNA and nuDNA using a reduced data set for which fossil calibration is available, and discuss the potential determinants of mutation rate variation across genomes and taxa. This study has important implications regarding DNA-based identification methods in predicting that mtDNA barcoding should be less reliable in nonvertebrates than in vertebrates.

Instant taxonomy: choosing adequate characters for species delimitation and description through congruence between molecular data and quantitative shape analysis

The lack of university funding is one of the major impediments to taxonomy, partly because traditional taxonomic training takes longer than a PhD course. Understanding ranges of phenotypic variability for different morphological structures, and their use as characters for delimitation and description of taxa, is a tedious task. We argue that the advent of molecular barcoding and quantitative shape analysis makes it unnecessary. As an example, we tackle a problematic species-complex of marine copepods from Korea and Japan, approaching it as a starting taxonomist might. Samples were collected from 14 locations and the mitochondrial COI gene was sequenced from 42 specimens. Our phylogenetic analyses reveal four distinct clades in Korea and Japan, and an additional nine belonging to a closely related complex from other parts of the Northern Pacific. Twenty different morphological structures were analysed for one Japanese and two Korean clades using landmark-based two-dimensional geometric morphometrics. Although there is no single morphological character that can distinguish with absolute certainty all three cryptic species, most show statistically significant interspecific differences in shape and size. We use five characters to describe two new species from Korea and to re-describe Tigriopus japonicus Mori, 1938 from near its type locality.

Diversity and evolution of sex determination systems in terrestrial isopods

Sex determination systems are highly variable in many taxa, sometimes even between closely related species. Yet the number and direction of transitions between these systems have seldom been characterized, and the underlying mechanisms are still poorly understood. Here we generated transcriptomes for 19 species of terrestrial isopod crustaceans, many of which are infected by Wolbachia bacterial endosymbionts. Using 88 single-copy orthologous genes, we reconstructed a fully resolved and dated phylogeny of terrestrial isopods. An original approach involving crossings of sex-reversed individuals allowed us to characterize the heterogametic systems of five species (one XY/XX and four ZW/ZZ). Mapping of these and previously known heterogametic systems onto the terrestrial isopod phylogeny revealed between 3 and 13 transitions of sex determination systems during the evolution of these taxa, most frequently from female to male heterogamety. Our results support that WW individuals are viable in many species, suggesting sex chromosomes are at an incipient stage of their evolution. Together, these data are consistent with the hypothesis that nucleo-cytoplasmic conflicts generated by Wolbachia endosymbionts triggered recurrent turnovers of sex determination systems in terrestrial isopods. They further establish terrestrial isopods as a model to study evolutionary transitions in sex determination systems and pave the way to molecularly characterize these systems.

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.

Chromosome-Wide Impacts on the Expression of Incompatibilities in Hybrids of Tigriopus californicus

Chromosome rearrangements such as inversions have been recognized previously as contributing to reproductive isolation by maintaining alleles together that jointly contribute to deleterious genetic interactions and postzygotic reproductive isolation. In this study, an impact of potential incompatibilities merely residing on the same chromosome was found in crosses of populations of the copepod Tigriopus californicus When genetically divergent populations of this copepod are crossed, hybrids show reduced fitness, and deviations from expected genotypic ratios can be used to determine regions of the genome involved in deleterious interactions. In this study, a set of markers was genotyped for a cross of two populations of T. californicus, and these markers show widespread deviations from Mendelian expectations, with entire chromosomes showing marked skew. Despite the importance of mtDNA/nuclear interactions in incompatibilities in this system in previous studies, in these crosses the expected patterns stemming from these interactions are not widely apparent. Females lack recombination in this species, and a striking difference is observed between male and female backcrosses. This suggests that the maintenance of multiple loci on individual chromosomes can enable some incompatibilities, perhaps playing a similar role in the initial rounds of hybridization to chromosomal rearrangements in preserving sets of alleles together that contribute to incompatibilities. Finally, it was observed that candidate pairs of incompatibility regions are not consistently interacting across replicates or subsets of these crosses, despite the repeatability of the deviations at many of the single loci themselves, suggesting that more complicated models of Dobzhansky-Muller incompatibilities may need to be considered.

Co-Speciation of the Ectoparasite Gyrodactylus teuchis (Monogenea, Platyhelminthes) and Its Salmonid Hosts

Co-speciation is a fundamental concept of evolutionary biology and intuitively appealing, yet in practice hard to demonstrate as it is often blurred by other evolutionary processes. We investigate the phylogeographic history of the monogenean ectoparasites Gyrodactylus teuchis and G. truttae on European salmonids of the genus Salmo. Mitochondrial cytochrome oxidase subunit 1 and the nuclear ribosomal internal transcribed spacer 2 were sequenced for 189 Gyrodactylus individuals collected from 50 localities, distributed across most major European river systems, from the Iberian- to the Balkan Peninsula. Despite both anthropogenic and naturally caused admixture of the principal host lineages among major river basins, co-phylogenetic analyses revealed significant global congruence for host and parasite phylogenies, providing firm support for co-speciation of G. teuchis and its salmonid hosts brown trout (S. trutta) and Atlantic salmon (S. salar). The major split within G. teuchis, coinciding with the initial divergence of the hosts was dated to ~1.5 My BP, using a Bayesian framework based on an indirect calibration point obtained from the host phylogeny. The presence of G. teuchis in Europe thus predates some of the major Pleistocene glaciations. In contrast, G. truttae exhibited remarkably low intraspecific genetic diversity. Given the direct life cycle and potentially high transmission potential of gyrodactylids, this finding is interpreted as indication for a recent emergence (<60 ky BP) of G. truttae via a host-switch. Our study thus suggests that instances of two fundamentally different mechanisms of speciation (co-speciation vs. host-switching) may have occurred on the same hosts in Europe within a time span of less than 1.5 My in two gyrodactylid ectoparasite species.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Cytoplasmic genetic variation and extensive cytonuclear interactions influence natural variation in the metabolome

Understanding genome to phenotype linkages has been greatly enabled by genomic sequencing. However, most genome analysis is typically confined to the nuclear genome. We conducted a metabolomic QTL analysis on a reciprocal RIL population structured to examine how variation in the organelle genomes affects phenotypic variation. This showed that the cytoplasmic variation had effects similar to, if not larger than, the largest individual nuclear locus. Inclusion of cytoplasmic variation into the genetic model greatly increased the explained phenotypic variation. Cytoplasmic genetic variation was a central hub in the epistatic network controlling the plant metabolome. This epistatic influence manifested such that the cytoplasmic background could alter or hide pairwise epistasis between nuclear loci. Thus, cytoplasmic genetic variation plays a central role in controlling natural variation in metabolomic networks. This suggests that cytoplasmic genomes must be included in any future analysis of natural variation. DOI: http://dx.doi.org/10.7554/eLife.00776.001.

Postzygotic isolation involves strong mitochondrial and sex-specific effects in Tigriopus californicus, a species lacking heteromorphic sex chromosomes

Detailed studies of the genetics of speciation have focused on a few model systems, particularly Drosophila. The copepod Tigriopus californicus offers an alternative that differs from standard animal models in that it lacks heteromorphic chromosomes (instead, sex determination is polygenic) and has reduced opportunities for sexual conflict, because females mate only once. Quantitative trait loci (QTL) mapping was conducted on reciprocal F2 hybrids between two strongly differentiated populations, using a saturated linkage map spanning all 12 autosomes and the mitochondrion. By comparing sexes, a possible sex ratio distorter was found but no sex chromosomes. Although studies of standard models often find an excess of hybrid male sterility factors, we found no QTL for sterility and multiple QTL for hybrid viability (indicated by non-Mendelian adult ratios) and other characters. Viability problems were found to be stronger in males, but the usual explanations for weaker hybrid males (sex chromosomes, sensitivity of spermatogenesis, sexual selection) cannot fully account for these male viability problems. Instead, higher metabolic rates may amplify deleterious effects in males. Although many studies of standard speciation models find the strongest genetic incompatibilities to be nuclear-nuclear (specifically X chromosome-autosome), we found the strongest deleterious interaction in this system was mito-nuclear. Consistent with the snowball theory of incompatibility accumulation, we found that trigenic interactions in this highly divergent cross were substantially more frequent (>6×) than digenic interactions. This alternative system thus allows important comparisons to studies of the genetics of reproductive isolation in more standard model systems.

Gene conversion yields novel gene combinations in paralogs of GOT1 in the copepod Tigriopus californicus

Background

Gene conversion of duplicated genes can slow the divergence of paralogous copies over time but can also result in other interesting evolutionary patterns. Islands of genetic divergence that persist in the face of gene conversion can point to gene regions undergoing selection for new functions. Novel combinations of genetic variation that differ greatly from the original sequence can result from the transfer of genetic variation between paralogous genes by rare gene conversion events. Genetically divergent populations of the copepod Tigriopus californicus provide an excellent model to look at the patterns of divergence among paralogs across multiple independent evolutionary lineages.

Results

In this study the evolution of a set of paralogous genes encoding putative aspartate transaminase proteins (called GOT1 here) are examined in populations of the copepod T. californicus. One pair of duplicated genes, GOT1p1 and GOT1p2, has regions of high divergence between the copies in the face of apparent on-going gene conversion. The GOT1p2 gene also has unique haplotypes in two populations that appear to have resulted from a transfer of genetic variation via inter-paralog gene conversion. A second pair of duplicated genes GOT1Sr and GOT1Sd also shows evidence of gene conversion, but this gene conversion does not appear to have maintained each as a functional copy in all populations.

Conclusions

The patterns of conservation and sequence divergence across this set of paralogous genes among populations of T. californicus suggest that some interesting evolutionary patterns are occurring at these loci. The results for the GOT1p1/GOT1p2 paralogs illustrate how gene conversion can factor in the creation of a mosaic pattern of regions of high divergence and low divergence. When coupled with rare gene conversion events of divergent regions, this pattern can result in the formation of novel proteins differing substantially from either original protein. The evolutionary patterns across these paralogs show how gene conversion can both constrain and facilitate diversification of genetic sequences.

A disproportionate role for mtDNA in Dobzhansky-Muller incompatibilities?

Evolution in allopatric populations can lead to incompatibilities that result in reduced hybrid fitness and ultimately reproductive isolation upon secondary contact. The Dobzhansky-Muller (DM) model nicely accounts for the evolution of such incompatibilities. Although DM incompatibilities were originally conceived as resulting of interactions between nuclear genes, recent studies have documented cases where incompatibilities have arisen between nuclear and mitochondrial genomes (mtDNA). Although mtDNA comprises only a tiny component (typically <<0.01%) of an organism's genetic material, several features of mtDNA may lead to a disproportionate contribution to the evolution of hybrid incompatibilities: (i) essentially all functions of mtDNA require interaction with nuclear gene products. All mtDNA-encoded proteins are components of the oxidative phosphorylation (OXPHOS) system and all mtDNA-encoded RNAs are part of the mitochondrial protein synthetic machinery; both processes require interaction with nuclear-encoded proteins for function. (ii) Transcription and replication of mtDNA also involve mitonuclear interactions as nuclear-encoded proteins must bind to regulatory motifs in the mtDNA to initiate these processes. (iii) Although features of mtDNA vary amongst taxa, metazoan mtDNA is typically characterized by high nucleotide substitution rates, lack of recombination and reduced effective population sizes that collectively lead to increased chance fixation of mildly deleterious mutations. Combined, these features create an evolutionary dynamic where rapid mtDNA evolution favours compensatory nuclear gene evolution, ultimately leading to co-adaptation of mitochondrial and nuclear genomes. When previously isolated lineages hybridize in nature or in the lab, intergenomic co-adaptation is disrupted and hybrid breakdown is observed; the role of intergenomic co-adaptation in hybrid breakdown and speciation will generally be most pronounced when rates of mtDNA evolution are high or when restricted gene flow results in significant population differentiation.