Evidence for compensatory evolution of ribosomal proteins in response to rapid divergence of mitochondrial rRNA

Selection Increases Mitonuclear DNA Discordance but Reconciles Incompatibility in African Cattle

Mitochondrial function relies on the coordinated interactions between genes in the mitochondrial DNA and nuclear genomes. Imperfect interactions following mitonuclear incompatibility may lead to reduced fitness. Mitochondrial DNA introgressions across species and populations are common and well documented. Various strategies may be expected to reconcile mitonuclear incompatibility in hybrids or admixed individuals. African admixed cattle (Bos taurus × B. indicus) show sex-biased admixture, with taurine (B. taurus) mitochondrial DNA and a nuclear genome predominantly of humped zebu (B. indicus). Here, we leveraged local ancestry inference approaches to identify the ancestry and distribution patterns of nuclear functional genes associated with the mitochondrial oxidative phosphorylation process in the genomes of African admixed cattle. We show that most of the nuclear genes involved in mitonuclear interactions are under selection and of humped zebu ancestry. Variations in mitochondrial DNA copy number may have contributed to the recovery of optimal mitochondrial function following admixture with the regulation of gene expression, alleviating or nullifying mitochondrial dysfunction. Interestingly, some nuclear mitochondrial genes with enrichment in taurine ancestry may have originated from ancient African aurochs (B. primigenius africanus) introgression. They may have contributed to the local adaptation of African cattle to pathogen burdens. Our study provides further support and new evidence showing that the successful settlement of cattle across the continent was a complex mechanism involving adaptive introgression, mitochondrial DNA copy number variation, regulation of gene expression, and selection of ancestral mitochondria-related genes.

Rapamycin's lifespan effect is modulated by mito-nuclear epistasis in Drosophila

The macrolide drug rapamycin is a benchmark anti-ageing drug, which robustly extends lifespan of diverse organisms. For any health intervention, it is paramount to establish whether benefits are distributed equitably among individuals and populations, and ideally to match intervention to recipients' needs. However, how responses to rapamycin vary is surprisingly understudied. Here we investigate how among-population variation in both mitochondrial and nuclear genetics shapes rapamycin's effects on lifespan. We show that epistatic "mito-nuclear" interactions, between mitochondria and nuclei, modulate the response to rapamycin treatment. Differences manifest as differential demographic effects of rapamycin, with altered age-specific mortality rate. However, a fitness cost of rapamycin early in life does not show a correlated response, suggesting that mito-nuclear epistasis can decouple costs and benefits of treatment. These findings suggest that a deeper understanding of how variation in mitochondrial and nuclear genomes shapes physiology may facilitate tailoring of anti-ageing therapy to individual need.

Nuclear compensatory evolution driven by mito-nuclear incompatibilities

Mitochondrial function relies on the coordinated expression of mitochondrial and nuclear genes, exhibiting remarkable resilience despite high mitochondrial mutation rates. The nuclear compensation mechanism suggests deleterious mitochondrial alleles drive compensatory nuclear mutations to preserve mito-nuclear compatibility. However, prevalence and factors conditioning this phenomenon remain debated due to its conflicting evidence. Here, we investigate how mito-nuclear incompatibilities impact substitutions in a model for species radiation. Mating success depends on genetic compatibility (nuclear DNA) and spatial proximity. Populations evolve from partially compatible mito-nuclear states, simulating mitochondrial DNA (mtDNA) introgression. Mutations do not confer advantages nor disadvantages, but individual fecundity declines with increasing incompatibilities, selecting for mito-nuclear coordination. We find that selection for mito-nuclear compatibility affects each genome differently based on their initial state. In compatible gene pairs, selection reduces substitutions in both genomes, while in incompatible nuclear genes, it consistently promotes compensation, facilitated by more mismatches. Interestingly, high mitochondrial mutation rates can reduce nuclear compensation by increasing mtDNA rectification, while substitutions in initially compatible nuclear gene are boosted. Finally, the presence of incompatibilities accelerates species radiation, but equilibrium richness is not directly correlated to substitution rates, revealing the complex dynamics triggered by mitochondrial introgression and mito-nuclear coevolution. Our study provides a perspective on nuclear compensation and the role of mito-nuclear incompatibilities in speciation by exploring extreme scenarios and identifying trends that empirical data alone cannot reveal. We emphasize the challenges in detecting these dynamics and propose analyzing specific genomic signatures could shed light on this evolutionary process.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BACKGROUND

Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.

RESULTS

We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.

CONCLUSIONS

Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Mitonuclear interactions shape both direct and parental effects of diet on fitness and involve a SNP in mitoribosomal 16s rRNA

Nutrition is a primary determinant of health, but responses to nutrition vary with genotype. Epistasis between mitochondrial and nuclear genomes may cause some of this variation, but which mitochondrial loci and nutrients participate in complex gene-by-gene-by-diet interactions? Furthermore, it remains unknown whether mitonuclear epistasis is involved only in the immediate responses to changes in diet, or whether mitonuclear genotype might modulate sensitivity to variation in parental nutrition, to shape intergenerational fitness responses. Here, in Drosophila melanogaster, we show that mitonuclear epistasis shapes fitness responses to variation in dietary lipids and amino acids. We also show that mitonuclear genotype modulates the parental effect of dietary lipid and amino acid variation on offspring fitness. Effect sizes for the interactions between diet, mitogenotype, and nucleogenotype were equal to or greater than the main effect of diet for some traits, suggesting that dietary impacts cannot be understood without first accounting for these interactions. Associating phenotype to mtDNA variation in a subset of populations implicated a C/T polymorphism in mt:lrRNA, which encodes the 16S rRNA of the mitochondrial ribosome. This association suggests that directionally different responses to dietary changes can result from variants on mtDNA that do not change protein coding sequence, dependent on epistatic interactions with variation in the nuclear genome.

Mutation pressure, drift, and the pace of molecular coevolution

Most aspects of the molecular biology of cells involve tightly coordinated intermolecular interactions requiring specific recognition at the nucleotide and/or amino acid levels. This has led to long-standing interest in the degree to which constraints on interacting molecules result in conserved vs. accelerated rates of sequence evolution, with arguments commonly being made that molecular coevolution can proceed at rates exceeding the neutral expectation. Here, a fairly general model is introduced to evaluate the degree to which the rate of evolution at functionally interacting sites is influenced by effective population sizes (Ne), mutation rates, strength of selection, and the magnitude of recombination between sites. This theory is of particular relevance to matters associated with interactions between organelle- and nuclear-encoded proteins, as the two genomic environments often exhibit dramatic differences in the power of mutation and drift. Although genes within low Ne environments can drive the rate of evolution of partner genes experiencing higher Ne, rates exceeding the neutral expectation require that the former also have an elevated mutation rate. Testable predictions, some counterintuitive, are presented on how patterns of coevolutionary rates should depend on the relative intensities of drift, selection, and mutation.

Evolution of loss of heterozygosity patterns in hybrid genomes of Candida yeast pathogens

BACKGROUND

Hybrids are chimeric organisms with highly plastic heterozygous genomes that may confer unique traits enabling the adaptation to new environments. However, most evolutionary theory frameworks predict that the high levels of genetic heterozygosity present in hybrids from divergent parents are likely to result in numerous deleterious epistatic interactions. Under this scenario, selection is expected to favor recombination events resulting in loss of heterozygosity (LOH) affecting genes involved in such negative interactions. Nevertheless, it is so far unknown whether this phenomenon actually drives genomic evolution in natural populations of hybrids. To determine the balance between selection and drift in the evolution of LOH patterns in natural yeast hybrids, we analyzed the genomic sequences from fifty-five hybrid strains of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis, which derived from at least six distinct natural hybridization events.

RESULTS

We found that, although LOH patterns in independent hybrid clades share some level of convergence that would not be expected from random occurrence, there is an apparent lack of strong functional selection. Moreover, while mitosis is associated with a limited number of inter-homeologous chromosome recombinations in these genomes, induced DNA breaks seem to increase the LOH rate. We also found that LOH does not accumulate linearly with time in these hybrids. Furthermore, some C. orthopsilosis hybrids present LOH patterns compatible with footprints of meiotic recombination. These meiotic-like patterns are at odds with a lack of evidence of sexual recombination and with our inability to experimentally induce sporulation in these hybrids.

CONCLUSIONS

Our results suggest that genetic drift is the prevailing force shaping LOH patterns in these hybrid genomes. Moreover, the observed LOH patterns suggest that these are likely not the result of continuous accumulation of sporadic events-as expected by mitotic repair of rare chromosomal breaks-but rather of acute episodes involving many LOH events in a short period of time.

Adaptive evolution characteristics of mitochondrial genomes in genus Aparapotamon (Brachyura, Potamidae) of freshwater crabs

BACKGROUND

Aparapotamon, a freshwater crab genus endemic to China, includes 13 species. The distribution of Aparapotamon spans the first and second tiers of China's terrain ladder, showing great altitudinal differences. To study the molecular mechanisms of adaptive evolution in Aparapotamon, we performed evolutionary analyses, including morphological, geographical, and phylogenetic analyses and divergence time estimation. We sequenced the mitogenomes of Aparapotamon binchuanense and Aparapotamon huizeense for the first time and resequenced three other mitogenomes of Aparapotamon grahami and Aparapotamon gracilipedum. These sequences were combined with NCBI sequences to perform comparative mitogenome analysis of all 13 Aparapotamon species, revealing mitogenome arrangement and the characteristics of protein-coding and tRNA genes.

RESULTS

A new species classification scheme of the genus Aparapotamon has been detected and verified by different aspects, including geographical, morphological, phylogenetics and comparative mitogenome analyses. Imprints from adaptive evolution were discovered in the mitochondrial genomes of group A, including the same codon loss at position 416 of the ND6 gene and the unique arrangement pattern of the tRNA-Ile gene. Multiple tRNA genes conserved or involved in adaptive evolution were detected. Two genes associated with altitudinal adaptation, ATP8 and ND6, which experienced positive selection, were identified for the first time in freshwater crabs.

CONCLUSIONS

Geological movements of the Qinghai-Tibet Plateau and Hengduan Mountains likely strongly impacted the speciation and differentiation of the four Aparapotamon groups. After some group A species dispersed from the Hengduan Mountain Range, new evolutionary characteristics emerged in their mitochondrial genomes, facilitating adaptation to the low-altitude environment of China's second terrain tier. Ultimately, group A species spread to high latitudes along the upper reaches of the Yangtze River, showing faster evolutionary rates, higher species diversity and the widest distribution.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Mitochondrial RNA processing defect caused by a SUPV3L1 mutation in two siblings with a novel neurodegenerative syndrome

SUPV3L1 encodes a helicase that is mainly localized in the mitochondria. It has been shown in vitro to possess both double-stranded RNA and DNA unwinding activity that is ATP-dependent. Here we report the first two patients for this gene who presented with a homozygous preliminary stop codon resulting in a C-terminal truncation of the SUPV3L1 protein. They presented with a characteristic phenotype of neurodegenerative nature with progressive spastic paraparesis, growth restriction, hypopigmentation, and predisposition to autoimmune disease. Ophthalmological examination showed severe photophobia with corneal erosions, optic atrophy, and pigmentary retinopathy, while neuroimaging showed atrophy of the optic chiasm and the pons with calcification of putamina, with intermittent and mild elevation of lactate. We show that the amino acids that are eliminated by the preliminary stop codon are highly conserved and are predicted to form an amphipathic helix. To investigate if the mutation causes mitochondrial dysfunction, we examined fibroblasts of the proband. We observed very low expression of the truncated protein, a reduction in the mature ND6 mRNA species as well as the accumulation of double-stranded RNA. Lentiviral complementation with the full-length SUPV3L1 cDNA partly restored the observed RNA phenotypes, supporting that the SUPV3L1 mutation in these patients is pathogenic and the cause of the disease.

The role of mitonuclear incompatibilities in allopatric speciation

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Temperature-Specific and Sex-Specific Fitness Effects of Sympatric Mitochondrial and Mito-Nuclear Variation in Drosophila obscura

Simple Summary

Does variation in the mitochondrial DNA sequence influence the survival and reproduction of an individual? What is the purpose of genetic variation of the mitochondrial DNA between individuals from the same population? As a simple laboratory model, Drosophila species can give us the answer to this question. Creating experimental lines with different combinations of mitochondrial and nuclear genomic DNA and testing how successful these lines were in surviving in different experimental set-ups enables us to deduce the effect that both genomes have on fitness. This study on D. obscura experimentally validates theoretical models that explain the persistence of mitochondrial DNA variation within populations. Our results shed light on the various mechanisms that maintain this type of variation. Finally, by conducting the experiments on two experimental temperatures, we have shown that environmental variations can support mitochondrial DNA variation within populations.

Abstract

The adaptive significance of sympatric mitochondrial (mtDNA) variation and the role of selective mechanisms that maintain it are debated to this day. Isofemale lines of Drosophila obscura collected from four populations were backcrossed within populations to construct experimental lines, with all combinations of mtDNA Cyt b haplotypes and nuclear genetic backgrounds (nuDNA). Individuals of both sexes from these lines were then subjected to four fitness assays (desiccation resistance, developmental time, egg-to-adult viability and sex ratio) on two experimental temperatures to examine the role of temperature fluctuations and sex-specific selection, as well as the part that interactions between the two genomes play in shaping mtDNA variation. The results varied across populations and fitness components. In the majority of comparisons, they show that sympatric mitochondrial variants affect fitness. However, their effect should be examined in light of interactions with nuDNA, as mito-nuclear genotype was even more influential on fitness across all components. We found both sex-specific and temperature-specific differences in mitochondrial and mito-nuclear genotype ranks in all fitness components. The effect of temperature-specific selection was found to be more prominent, especially in desiccation resistance. From the results of different components tested, we can also infer that temperature-specific mito-nuclear interactions rather than sex-specific selection on mito-nuclear genotypes have a more substantial role in preserving mtDNA variation in this model species.

Mitonuclear Interactions Produce Diverging Responses to Mild Stress in Drosophila Larvae

Mitochondrial function depends on direct interactions between respiratory proteins encoded by genes in two genomes, mitochondrial and nuclear, which evolve in very different ways. Serious incompatibilities between these genomes can have severe effects on development, fitness and viability. The effect of subtle mitonuclear mismatches has received less attention, especially when subject to mild physiological stress. Here, we investigate how two distinct physiological stresses, metabolic stress (high-protein diet) and redox stress [the glutathione precursor N-acetyl cysteine (NAC)], affect development time, egg-to-adult viability, and the mitochondrial physiology of Drosophila larvae with an isogenic nuclear background set against three mitochondrial DNA (mtDNA) haplotypes: one coevolved (WT) and two slightly mismatched (COX and BAR). Larvae fed the high-protein diet developed faster and had greater viability in all haplotypes. The opposite was true of NAC-fed flies, especially those with the COX haplotype. Unexpectedly, the slightly mismatched BAR larvae developed fastest and were the most viable on both treatments, as well as control diets. These changes in larval development were linked to a shift to complex I-driven mitochondrial respiration in all haplotypes on the high-protein diet. In contrast, NAC increased respiration in COX larvae but drove a shift toward oxidation of proline and succinate. The flux of reactive oxygen species was increased in COX larvae treated with NAC and was associated with an increase in mtDNA copy number. Our results support the notion that subtle mitonuclear mismatches can lead to diverging responses to mild physiological stress, undermining fitness in some cases, but surprisingly improving outcomes in other ostensibly mismatched fly lines.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Phenotypic Variation in Mitochondria-Related Performance Traits Across New Zealand Snail Populations

Mitochondrial function is critical for energy homeostasis and should shape how genetic variation in metabolism is transmitted through levels of biological organization to generate stability in organismal performance. Mitochondrial function is encoded by genes in two distinct and separately inherited genomes-the mitochondrial genome and the nuclear genome-and selection is expected to maintain functional mito-nuclear interactions. The documented high levels of polymorphism in genes involved in these mito-nuclear interactions and wide variation for mitochondrial function demands an explanation for how and why variability in such a fundamental trait is maintained. Potamopyrgus antipodarum is a New Zealand freshwater snail with coexisting sexual and asexual individuals and, accordingly, contrasting systems of separate vs. co-inheritance of nuclear and mitochondrial genomes. As such, this snail provides a powerful means to dissect the evolutionary and functional consequences of mito-nuclear variation. The lakes inhabited by P. antipodarum span wide environmental gradients, with substantial across-lake genetic structure and mito-nuclear discordance. This situation allows us to use comparisons across reproductive modes and lakes to partition variation in cellular respiration across genetic and environmental axes. Here, we integrated cellular, physiological, and behavioral approaches to quantify variation in mitochondrial function across a diverse set of wild P. antipodarum lineages. We found extensive across-lake variation in organismal oxygen consumption and behavioral response to heat stress and differences across sexes in mitochondrial membrane potential but few global effects of reproductive mode. Taken together, our data set the stage for applying this important model system for sexual reproduction and polyploidy to dissecting the complex relationships between mito-nuclear variation, performance, plasticity, and fitness in natural populations.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

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.

Mitochondrial Short-Term Plastic Responses and Long-Term Evolutionary Dynamics in Animal Species

How do species respond or adapt to environmental changes? The answer to this depends partly on mitochondrial epigenetics and genetics, new players in promoting adaptation to both short- and long-term environmental changes. In this review, we explore how mitochondrial epigenetics and genetics mechanisms, such as mtDNA methylation, mtDNA-derived noncoding RNAs, micropeptides, mtDNA mutations, and adaptations, can contribute to animal plasticity and adaptation. We also briefly discuss the challenges in assessing mtDNA adaptive evolution. In sum, this review covers new advances in the field of mitochondrial genomics, many of which are still controversial, and discusses processes still somewhat obscure, and some of which are still quite speculative and require further robust experimentation.

Correlated Evolution of Large DNA Fragments in the 3D Genome of Arabidopsis thaliana

In eukaryotes, the three-dimensional (3D) conformation of the genome is far from random, and this nonrandom chromatin organization is strongly correlated with gene expression and protein function, which are two critical determinants of the selective constraints and evolutionary rates of genes. However, whether genes and other elements that are located close to each other in the 3D genome evolve in a coordinated way has not been investigated in any organism. To address this question, we constructed chromatin interaction networks (CINs) in Arabidopsis thaliana based on high-throughput chromosome conformation capture data and demonstrated that adjacent large DNA fragments in the CIN indeed exhibit more similar levels of polymorphism and evolutionary rates than random fragment pairs. Using simulations that account for the linear distance between fragments, we proved that the 3D chromosomal organization plays a role in the observed correlated evolution. Spatially interacting fragments also exhibit more similar mutation rates and functional constraints in both coding and noncoding regions than the random expectations, indicating that the correlated evolution between 3D neighbors is a result of combined evolutionary forces. A collection of 39 genomic and epigenomic features can explain much of the variance in genetic diversity and evolutionary rates across the genome. Moreover, features that have a greater effect on the evolution of regional sequences tend to show higher similarity between neighboring fragments in the CIN, suggesting a pivotal role of epigenetic modifications and chromatin organization in determining the correlated evolution of large DNA fragments in the 3D genome.

Mitochondrial-nuclear coadaptation revealed through mtDNA replacements in Saccharomyces cerevisiae

BACKGROUND

Mitochondrial function requires numerous genetic interactions between mitochondrial- and nuclear- encoded genes. While selection for optimal mitonuclear interactions should result in coevolution between both genomes, evidence for mitonuclear coadaptation is challenging to document. Genetic models where mitonuclear interactions can be explored are needed.

RESULTS

We systematically exchanged mtDNAs between 15 Saccharomyces cerevisiae isolates from a variety of ecological niches to create 225 unique mitochondrial-nuclear genotypes. Analysis of phenotypic profiles confirmed that environmentally-sensitive interactions between mitochondrial and nuclear genotype contributed to growth differences. Exchanges of mtDNAs between strains of the same or different clades were just as likely to demonstrate mitonuclear epistasis although epistatic effect sizes increased with genetic distances. Strains with their original mtDNAs were more fit than strains with synthetic mitonuclear combinations when grown in media that resembled isolation habitats.

CONCLUSIONS

This study shows that natural variation in mitonuclear interactions contributes to fitness landscapes. Multiple examples of coadapted mitochondrial-nuclear genotypes suggest that selection for mitonuclear interactions may play a role in helping yeasts adapt to novel environments and promote coevolution.

Genetic hitchhiking, mitonuclear coadaptation, and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

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.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Mitonuclear Compensatory Coevolution

In bilaterian animals, the mitochondrial genome is small, haploid, does not typically recombine, and is subject to accumulation of deleterious alleles via Muller's ratchet. These basic features of the genomic architecture present a paradox: mutational erosion of these genomes should lead to decline in mitochondrial function over time, yet no such decline is observed. Compensatory coevolution, whereby the nuclear genome evolves to compensate for the deleterious alleles in the mitochondrial genome, presents a potential solution to the paradox of Muller's ratchet without loss of function. Here, I review different proposed forms of mitonuclear compensatory coevolution. Empirical evidence from diverse eukaryotic taxa supports the mitonuclear compensatory coevolution hypothesis, but the ubiquity and importance of such compensatory coevolution remains a topic of debate.

Strong selective effects of mitochondrial DNA on the nuclear genome

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F₂ hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Sexually Antagonistic Mitonuclear Coevolution in Duplicate Oxidative Phosphorylation Genes

Mitochondrial function is critical in eukaryotes. To maintain an adequate supply of energy, precise interactions must be maintained between nuclear- and mitochondrial-encoded gene products. Such interactions are paramount in chimeric enzymes such as the oxidative phosphorylation (OXPHOS) complexes. Mutualistic coevolution between the two genomes has therefore been suggested to be a critical, ubiquitous feature of eukaryotes that acts to maintain cellular function. However, mitochondrial genomes can also act selfishly and increase their own transmission at the expense of organismal function. For example, male-harming mutations are predisposed to accumulate in mitochondrial genomes due to their maternal inheritance ("mother's curse"). Here, we investigate sexually antagonistic mitonuclear coevolution in nuclear-encoded OXPHOS paralogs from mammals and Drosophila. These duplicate genes are highly divergent but must interact with the same set of mitochondrial-encoded genes. Many such paralogs show testis-specific expression, prompting previous hypotheses suggesting they may have evolved under selection to counteract male-harming mitochondrial mutations. We found increased rates of evolution in OXPHOS paralogs with testis-specific expression in mammals and Drosophila, supporting this hypothesis. However, further analyses suggested such patterns may be due to relaxed, not positive selection, especially in Drosophila. Structural data also suggest that mitonuclear interactions do not play a major role in the evolution of many OXPHOS paralogs in a consistent way. In conclusion, no single OXPHOS paralog met all our criteria for being under selection to counteract male-harming mitochondrial mutations. We discuss alternative explanations for the drastic patterns of evolution in these genes, including mutualistic mitonuclear coevolution, adaptive subfunctionalization after gene duplication, and relaxed selection on OXPHOS in male tissues.

Impact of mitonuclear interactions on life-history responses to diet

Mitochondria are central to both energy metabolism and biosynthesis. Mitochondrial function could therefore influence resource allocation. Critically, mitochondrial function depends on interactions between proteins encoded by the mitochondrial and nuclear genomes. Severe incompatibilities between these genomes can have pervasive effects on both fitness and longevity. How milder deficits in mitochondrial function affect life-history trade-offs is less well understood. Here, we analyse how mitonuclear interactions affect the trade-off between fecundity and longevity in Drosophila melanogaster. We consider a panel of 10 different mitochondrial DNA haplotypes against two contrasting nuclear backgrounds (w¹¹¹⁸ (WE) and Zim53 (ZIM)) in response to high-protein versus standard diet. We report strikingly different responses between the two nuclear backgrounds. WE females have higher fecundity and decreased longevity on high protein. ZIM females have much greater fecundity and shorter lifespan than WE flies on standard diet. High protein doubled their fecundity with no effect on longevity. Mitochondrial haplotype reflected nuclear life-history trade-offs, with a negative correlation between longevity and fecundity in WE flies and no correlation in ZIM flies. Mitonuclear interactions had substantial effects but did not reflect genetic distance between mitochondrial haplotypes. We conclude that mitonuclear interactions can have significant impact on life-history trade-offs, but their effects are not predictable by relatedness. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

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'.

Molecular function limits divergent protein evolution on planetary timescales

Functional conservation is known to constrain protein evolution. Nevertheless, the long-term divergence patterns of proteins maintaining the same molecular function and the possible limits of this divergence have not been explored in detail. We investigate these fundamental questions by characterizing the divergence between ancient protein orthologs with conserved molecular function. Our results demonstrate that the decline of sequence and structural similarities between such orthologs significantly slows down after ~1-2 billion years of independent evolution. As a result, the sequence and structural similarities between ancient orthologs have not substantially decreased for the past billion years. The effective divergence limit (>25% sequence identity) is not primarily due to protein sites universally conserved in all linages. Instead, less than four amino acid types are accepted, on average, per site across orthologous protein sequences. Our analysis also reveals different divergence patterns for protein sites with experimentally determined small and large fitness effects of mutations.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

CyMIRA: The Cytonuclear Molecular Interactions Reference for Arabidopsis

The function and evolution of eukaryotic cells depend upon direct molecular interactions between gene products encoded in nuclear and cytoplasmic genomes. Understanding how these cytonuclear interactions drive molecular evolution and generate genetic incompatibilities between isolated populations and species is of central importance to eukaryotic biology. Plants are an outstanding system to investigate such effects because of their two different genomic compartments present in the cytoplasm (mitochondria and plastids) and the extensive resources detailing subcellular targeting of nuclear-encoded proteins. However, the field lacks a consistent classification scheme for mitochondrial- and plastid-targeted proteins based on their molecular interactions with cytoplasmic genomes and gene products, which hinders efforts to standardize and compare results across studies. Here, we take advantage of detailed knowledge about the model angiosperm Arabidopsis thaliana to provide a curated database of plant cytonuclear interactions at the molecular level. CyMIRA (Cytonuclear Molecular Interactions Reference for Arabidopsis) is available at http://cymira.colostate.edu/ and https://github.com/dbsloan/cymira and will serve as a resource to aid researchers in partitioning evolutionary genomic data into functional gene classes based on organelle targeting and direct molecular interaction with cytoplasmic genomes and gene products. It includes 11 categories (and 27 subcategories) of different cytonuclear complexes and types of molecular interactions, and it reports residue-level information for cytonuclear contact sites. We hope that this framework will make it easier to standardize, interpret, and compare studies testing the functional and evolutionary consequences of cytonuclear interactions.

Evolutionary Rate Correlation between Mitochondrial-Encoded and Mitochondria-Associated Nuclear-Encoded Proteins in Insects

The mitochondrion is a pivotal organelle for energy production, and includes components encoded by both the mitochondrial and nuclear genomes. Functional and evolutionary interactions are expected between the nuclear- and mitochondrial-encoded components. The topic is of broad interest in biology, with implications to genetics, evolution, and medicine. Here, we compare the evolutionary rates of mitochondrial proteins and ribosomal RNAs to rates of mitochondria-associated nuclear-encoded proteins, across the major orders of holometabolous insects. There are significant evolutionary rate correlations (ERCs) between mitochondrial-encoded and mitochondria-associated nuclear-encoded proteins, which are likely driven by different rates of mitochondrial sequence evolution and correlated changes in the interacting nuclear-encoded proteins. The pattern holds after correction for phylogenetic relationships and considering protein conservation levels. Correlations are stronger for both nuclear-encoded OXPHOS proteins that are in contact with mitochondrial OXPHOS proteins and for nuclear-encoded mitochondrial ribosomal amino acids directly contacting the mitochondrial rRNAs. We find that ERC between mitochondrial- and nuclear-encoded proteins is a strong predictor of nuclear-encoded proteins known to interact with mitochondria, and ERC shows promise for identifying new candidate proteins with mitochondrial function. Twenty-three additional candidate nuclear-encoded proteins warrant further study for mitochondrial function based on this approach, including proteins in the minichromosome maintenance helicase complex.

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.

A systematic review and meta-analysis reveals pervasive effects of germline mitochondrial replacement on components of health

BACKGROUND

Mitochondrial replacement, a form of nuclear transfer, has been proposed as a germline therapy to prevent the transmission of mitochondrial diseases. Mitochondrial replacement therapy has been licensed for clinical application in the UK, and already carried out in other countries, but little is known about negative or unintended effects on the health of offspring born using this technique.

OBJECTIVE AND RATIONALE

Studies in invertebrate models have used techniques that achieve mitochondrial replacement to create offspring with novel combinations of mitochondrial and nuclear genotype. These have demonstrated that the creation of novel mitochondrial-nuclear interactions can lead to alterations in offspring characteristics, such as development rates, fertility and longevity. However, it is currently unclear whether such interactions could similarly affect the outcomes of vertebrate biomedical studies, which have sought to assess the efficacy of the replacement therapy.

SEARCH METHODS

This systematic review addresses whether the effects of mitochondrial replacement on offspring characteristics differ in magnitude between biological (conducted on invertebrate models, with an ecological or evolutionary focus) and biomedical studies (conducted on vertebrate models, with a clinical focus). Studies were selected based on a key-word search in 'Web of Science', complemented by backward searches of reviews on the topic of mitochondrial-nuclear (mito-nuclear) interactions. In total, 43 of the resulting 116 publications identified in the search contained reliable data to estimate effect sizes of mitochondrial replacement. We found no evidence of publication bias when examining effect-size estimates across sample sizes.

OUTCOMES

Mitochondrial replacement consistently altered the phenotype, with significant effects at several levels of organismal performance and health, including gene expression, anatomy, metabolism and life-history. Biomedical and biological studies, while differing in the methods used to achieve mitochondrial replacement, showed only marginally significant differences in effect-size estimates (-0.233 [CI: -0.495 to -0.011]), with larger effect-size estimates in biomedical studies (0.697 [CI: 0.450-0.956]) than biological studies (0.462 [CI: 0.287-0.688]). Humans showed stronger effects than other species. Effects of mitochondrial replacement were also stronger in species with a higher basal metabolic rate. Based on our results, we conducted the first formal risk analysis of mitochondrial replacement, and conservatively estimate negative effects in at least one in every 130 resulting offspring born to the therapy.

WIDER IMPLICATIONS

Our findings suggest that mitochondrial replacement may routinely affect offspring characteristics across a wide array of animal species, and that such effects are likely to extend to humans. Studies in invertebrate models have confirmed mito-nuclear interactions as the underpinning cause of organismal effects following mitochondrial replacement. This therefore suggests that mito-nuclear interactions are also likely to be contributing to effects seen in biomedical studies, on vertebrate models, whose effect sizes exceeded those of biological studies. Our results advocate the use of safeguards that could offset any negative effects (defining any unintended effect as being negative) mediated by mito-nuclear interactions following mitochondrial replacement in humans, such as mitochondrial genetic matching between donor and recipient. Our results also suggest that further research into the molecular nature of mito-nuclear interactions would be beneficial in refining the clinical application of mitochondrial replacement, and in establishing what degree of variation between donor and patient mitochondrial DNA haplotypes is acceptable to ensure 'haplotype matching'.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

Homologous recombination changes the context of Cytochrome b transcription in the mitochondrial genome of Silene vulgaris KRA

Background

Silene vulgaris (bladder campion) is a gynodioecious species existing as two genders – male-sterile females and hermaphrodites. Cytoplasmic male sterility (CMS) is generally encoded by mitochondrial genes, which interact with nuclear fertility restorer genes. Mitochondrial genomes of this species vary in DNA sequence, gene order and gene content. Multiple CMS genes are expected to exist in S. vulgaris, but little is known about their molecular identity.

Results

We assembled the complete mitochondrial genome from the haplotype KRA of S. vulgaris. It consists of five chromosomes, two of which recombine with each other. Two small non-recombining chromosomes exist in linear, supercoiled and relaxed circle forms. We compared the mitochondrial transcriptomes from females and hermaphrodites and confirmed the differentially expressed chimeric gene bobt as the strongest CMS candidate gene in S. vulgaris KRA. The chimeric gene bobt is co-transcribed with the Cytochrome b (cob) gene in some genomic configurations. The co-transcription of a CMS factor with an essential gene may constrain transcription inhibition as a mechanism for fertility restoration because of the need to maintain appropriate production of the necessary protein. Homologous recombination places the gene cob outside the control of bobt, which allows for the suppression of the CMS gene by the fertility restorer genes. We found the loss of three editing sites in the KRA mitochondrial genome and identified four sites with highly distinct editing rates between KRA and another S. vulgaris haplotypes (KOV). Three of these highly differentially edited sites were located in the transport membrane protein B (mttB) gene. They resulted in differences in MttB protein sequences between haplotypes.

Conclusions

Frequent homologous recombination events that are widespread in plant mitochondrial genomes may change chromosomal configurations and also the control of gene transcription including CMS gene expression. Posttranscriptional processes, e.g. RNA editing shall be evaluated in evolutionary and co-evolutionary studies of mitochondrial genes, because they may change protein composition despite the sequence identity of the respective genes. The investigation of natural populations of wild species such as S. vulgaris are necessary to reveal important aspects of CMS missed in domesticated crops, the traditional focus of the CMS studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5254-0) contains supplementary material, which is available to authorized users.

Does mitochondrial DNA evolution in metazoa drive the origin of new mitochondrial proteins?

Most eukaryotic cells contain mitochondria with a genome that evolved from their α-proteobacterial ancestor. In the course of eukaryotic evolution, the mitochondrial genome underwent a dramatic reduction in size, caused by the loss and translocation of genes. This required adjustments in mitochondrial gene expression mechanisms and resulted in a complex collaborative system of mitochondrially encoded transfer RNAs and ribosomal RNAs with nuclear encoded proteins to express the mitochondrial encoded oxidative phosphorylation (OXPHOS) proteins. In this review, we examine mitochondrial gene expression from an evolutionary point of view: to what extent can we correlate changes in the mitochondrial genome in the evolutionary lineage leading to human with the origin of new nuclear encoded proteins. We dated the evolutionary origin of mitochondrial proteins that interact with mitochondrial DNA or its RNA and/or protein products in a systematic manner and compared them with documented changes in the mitochondrial DNA. We find anecdotal but accumulating evidence that metazoan RNA-interacting proteins arose in conjunction with changes of the mitochondrial DNA. We find no substantial evidence for such compensatory evolution in new OXPHOS proteins, which appear to be constrained by the ability to form supercomplexes. © 2018 IUBMB Life, 70(12):1240-1250, 2018.

Cytonuclear integration and co-evolution

The partitioning of genetic material between the nucleus and cytoplasmic (mitochondrial and plastid) genomes within eukaryotic cells necessitates coordinated integration between these genomic compartments, with important evolutionary and biomedical implications. Classic questions persist about the pervasive reduction of cytoplasmic genomes via a combination of gene loss, transfer and functional replacement - and yet why they are almost always retained in some minimal form. One striking consequence of cytonuclear integration is the existence of 'chimeric' enzyme complexes composed of subunits encoded in two different genomes. Advances in structural biology and comparative genomics are yielding important insights into the evolution of such complexes, including correlated sequence changes and recruitment of novel subunits. Thus, chimeric cytonuclear complexes provide a powerful window into the mechanisms of molecular co-evolution.