Mitochondrial Involvement in Vertebrate Speciation? The Case of Mito-nuclear Genetic Divergence in Chameleons

Mitochondrial-related genome-wide Mendelian randomization identifies putatively causal genes in the pathogenesis of sepsis

BACKGROUND

The dysfunction of mitochondria has been associated with the development of sepsis, but the specific mitochondrial-related genes and their roles in sepsis have not been fully elucidated. We employed Mendelian randomization and colocalization analysis to investigate the association between mitochondrial-related genes and sepsis by integrating multi-omics data.

METHODS

Summary-level data on mitochondrial gene methylation, expression, and protein abundance levels were obtained from corresponding studies on methylation, expression, and protein quantitative trait loci, respectively. Genetic associations with sepsis were obtained from the genome-wide association studies catalog database. We used the MitoCarta3.0 database, which contains an updated list of 1,136 human mitochondrial genes, to identify mitochondrial genes. To assess the associations between mitochondrial gene-related molecular features and sepsis, we conducted summary data-based Mendelian randomization analysis. In addition, we performed colocalization analysis to determine whether the identified signal pairs shared a causal genetic variant.

RESULTS

After integrating the multi-omics data between methylation quantitative trait loci- expression quantitative trait loci and expression quantitative trait loci-protein quantitative trait loci, we identified FIS1 as having tier 1 evidence for its association with sepsis. Methylation of cg01299997 in FIS1 was found to be associated with lower expression of FIS1, an increased risk of sepsis, and a positive role in cg01299997 methylation. Furthermore, NUDT2, IMMP2L, LYRM4, MRPL10, MRPL17, MTIF3, and TFAM genes were associated with sepsis risk with tier 2 evidence. Both gene expression and protein abundance levels of NUDT2 were observed to be associated with an increased risk of sepsis. In addition, the ATP5MC1 and VWA8 genes were associated with sepsis risk with tier 3 evidence. Among these tertiary genes, ATP5MC1 gene expression level showed a negative correlation (posterior probability of H4 = 0.9242), whereas the gene expression level of VWA8 exhibited a positive correlation (posterior probability of H4 = 0.7270).

CONCLUSION

We found that the mitochondrial FIS1, NUDT2, IMMP2L, LYRM4, MRPL10, MRPL17, MTIF3, TFAM, ATP5MC1, and VWA8 genes were putatively associated with sepsis risk with evidence from multi-omics levels. This study identified mitochondrial genes in relation to sepsis, which may enhance the understanding of the pathogenic mechanisms of sepsis development.

Vertebrates show coordinated elevated expression of mitochondrial and nuclear genes after birth

Interactions between mitochondrial and nuclear factors are essential to life. Nevertheless, the importance of coordinated regulation of mitochondrial-nuclear gene expression (CMNGE) to changing physiological conditions is poorly understood and is limited to certain tissues and organisms. We hypothesized that CMNGE is important for development across vertebrates and, hence, should be conserved. As a first step, we analyzed more than 1400 RNA-seq experiments performed during prenatal development, in neonates, and in adults across vertebrate evolution. We find conserved sharp elevation of CMNGE after birth, including oxidative phosphorylation (OXPHOS) and mitochondrial ribosome genes, in the heart, hindbrain, forebrain, and kidney across mammals, as well as in Gallus gallus and in the lizard Anolis carolinensis This is accompanied by elevated expression of TCA cycle enzymes and reduction in hypoxia response genes, suggesting a conserved cross-tissue metabolic switch after birth/hatching. Analysis of about 70 known regulators of mitochondrial gene expression reveals consistently elevated expression of PPARGC1A (also known as Pgc-1alpha) and CEBPB after birth/hatching across organisms and tissues, thus highlighting them as candidate regulators of CMNGE upon transition to the neonate. Analyses of Danio rerio, Xenopus tropicalis, Caenorhabditis elegans, and Drosophila melanogaster reveal elevated CMNGE prior to hatching in X. tropicalis and in D. melanogaster, which is associated with the emergence of muscle activity. Lack of such an ancient pattern in mammals and in chickens suggests that it was lost during radiation of terrestrial vertebrates. Taken together, our results suggest that regulated CMNGE after birth reflects an essential metabolic switch that is under strong selective constraints.

Genomic data support reticulate evolution in whiptail lizards from the Brazilian Caatinga

Species relationships have traditionally been represented by phylogenetic trees, but not all evolutionary histories fit into bifurcating divergence models. Introgressive hybridization challenges this assumption by sometimes [or maybe often] leading to mitochondrial introgression, wherein one species' mitochondrial genome is entirely replaced by another's (mitochondrial capture). Such processes result in mitonuclear discrepancies, complicating species delimitation and phylogenetic inference. In our study, we used ultraconserved elements (UCE) and mitogenomic data to investigate the evolutionary history of the Ameivula ocellifera complex, a group of South American whiptail lizards widely distributed in semiarid environments of the Caatinga Domain in Brazil. We examine mitonuclear discordances, assessing reticulate evolution, evaluating species limits, and testing for adaptive mitochondrial capture that could explain higher introgression in the mitochondrial genome compared to nuclear DNA. Our findings support the occurrence of an ancient reticulation event during the diversification of these lizards, driven by introgressive hybridization, leading to mitochondrial capture, and explaining mitonuclear discrepancies. Overall, we did not find clear evidence of positive selection across mitochondrial protein-coding genes suggesting adaptive mitochondrial capture of individuals with introgressed mtDNA. Thus, the genetic diversification and mitogenome evolution could be neutral, with selection against hybridization in the autosomal loci only, or even mediated by mitonuclear incompatibilities. Analyses of mtDNA genomes alongside network and species delimitation methods were crucial for identifying and validating individuals with introgressed mtDNA as a distinct species, demonstrating the potential of genome sampling, and using innovative analytical techniques for elucidating speciation processes in the presence of introgressive hybridization.

Random Mating in a Hybrid Zone Between Two Putative Climate-Adapted Bird Lineages With Predicted Mitonuclear Incompatibilities

Biochemical and evolutionary interactions between mitochondrial and nuclear genomes ('mitonuclear interactions') are proposed to underpin fundamental aspects of biology including evolution of sexual reproduction, adaptation and speciation. We investigated the role of pre-mating isolation in maintaining functional mitonuclear interactions in wild populations bearing diverged, putatively co-adapted mitonuclear genotypes. Two lineages of eastern yellow robin Eopsaltria australis-putatively climate-adapted to 'inland' and 'coastal' climates-differ by ~7% of mitogenome nucleotides, whereas nuclear genome differences are concentrated into a sex-linked region enriched with mitochondrial functions. Female-specific selection and male-mediated gene flow across the hybrid zone where the lineages coexist and interbreed can explain this pattern. It remains unknown whether lineage divergence is driven by intrinsic incompatibilities (particularly in females; Haldane's rule), extrinsic selection, both, or other drivers. We tested whether lineage divergence could be facilitated by non-random mate-pairing with respect to partners' mitolineage and/or mitonuclear genes encoded by the Z sex-chromosome, which differ between the lineages. We used field-, Z-linked- and mitolineage data from two locations where lineages hybridise to test whether females mate disproportionately with (1) males of their own mitolineage and/or bearing similar Z-linked variation, as might be expected if hybrids experience intrinsic incompatibilities, or (2) putatively locally-adapted males, as might be expected under environmental selection. Comparing field observations with simulations provided no evidence of non-random mating, thus drivers of observed population genetic patterns are consistent with reduced female gene flow likely acting post-mating. Future tests of female-biased mortality at different life stages and habitat selection may clarify mechanisms of selection.

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Association Analysis Provides Insights into Plant Mitonuclear Interactions

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

Diet and mitonuclear haplotype interactions affect growth rate in a slime mould

Trait expression in metazoans is strongly influenced by the balance of macronutrients (i.e. protein, carbohydrate and fat) in the diet. At the same time, an individual's genetic background seems to regulate the magnitude of phenotypic response to a particular diet. It needs to be better understood whether interactions between diet, genetic background and trait expression are found in unicellular eukaryotes. A protist-the slime mould, Physarum polycephalum can choose diets based on protein-to-carbohydrate (P:C) content to support optimal growth rate. Yet, the role of genetic background (variation in the mitochondrial and nuclear DNAs) in mediating growth rate response to dietary P:C ratios in the slime mould is unknown. Here, we studied the effects of interactions between mitochondrial and nuclear DNA haplotypes and diet (i.e. G × G × E interactions) on the growth rate of P. polycephalum. A genetic panel of six distinct strains of P. polycephalum that differ in their mitochondrial and nuclear DNA haplotypes was used to measure growth rate across five diets that varied in their P:C ratio and total calories. We first determined the strains' growth rate (total biomass and surface area) when grown on a set menu with access to a particular diet. We then assessed whether the growth rate of strains increased on a buffet menu with access to all diets. Our findings show that the growth rate of P. polycephalum is generally higher on diets containing more carbohydrates than protein and that total calories negatively affect the growth rate. Three-way interactions between mitochondrial, nuclear haplotypes and dietary P:C ratios affected the strains' surface area of growth but not biomass. Intriguingly, strains did not increase their surface area and biomass when they had access to all diets on the buffet menu. Our findings have broad implications for our understanding of the effect of mitonuclear interactions on trait expression across diverse eukaryotic lineages.

Investigating the Impact of a Curse: Diseases, Population Isolation, Evolution and the Mother's Curse

The mitochondrion was characterized for years as the energy factory of the cell, but now its role in many more cellular processes is recognized. The mitochondrion and mitochondrial DNA (mtDNA) also possess a set of distinct properties, including maternal inheritance, that creates the Mother's Curse phenomenon. As mtDNA is inherited from females to all offspring, mutations that are harmful to males tend to accumulate more easily. The Mother's Curse is associated with various diseases, and has a significant effect on males, in many cases even affecting their reproductive ability. Sometimes, it even leads to reproductive isolation, as in crosses between different populations, the mitochondrial genome cannot cooperate effectively with the nuclear one resulting in a mito-nuclear incompatibility and reduce the fitness of the hybrids. This phenomenon is observed both in the laboratory and in natural populations, and have the potential to influence their evolution and speciation. Therefore, it turns out that the study of mitochondria is an exciting field that finds many applications, including pest control, and it can shed light on the molecular mechanism of several diseases, improving successful diagnosis and therapeutics. Finally, mito-nuclear co-adaptation, paternal leakage, and kin selection are some mechanisms that can mitigate the impact of the Mother's Curse.

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.

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.

OUP accepted manuscript

Recognition of mammalian mitochondrial promoters requires the concerted action of mitochondrial RNA polymerase (mtRNAP) and transcription initiation factors TFAM and TFB2M. In this work, we found that transcript slippage results in heterogeneity of the human mitochondrial transcripts in vivo and in vitro. This allowed us to correctly interpret the RNAseq data, identify the bona fide transcription start sites (TSS), and assign mitochondrial promoters for > 50% of mammalian species and some other vertebrates. The divergent structure of the mammalian promoters reveals previously unappreciated aspects of mtDNA evolution. The correct assignment of TSS also enabled us to establish the precise register of the DNA in the initiation complex and permitted investigation of the sequence-specific protein-DNA interactions. We determined the molecular basis of promoter recognition by mtRNAP and TFB2M, which cooperatively recognize bases near TSS in a species-specific manner. Our findings reveal a role of mitochondrial transcription machinery in mitonuclear coevolution and speciation.

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.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Signatures of mitonuclear coevolution in a warbler species complex

Divergent mitonuclear coadaptation could facilitate speciation. We investigate this possibility in two hybridizing species of warblers, Setophaga occidentalis and S. townsendi, in western North America. Inland S. townsendi harbor distinct mitochondrial DNA haplotypes from those of S. occidentalis. These populations also differ in several nuclear DNA regions. Coastal S. townsendi demonstrate mixed mitonuclear ancestry from S. occidentalis and inland S. townsendi. Of the few highly-differentiated chromosomal regions between inland S. townsendi and S. occidentalis, a 1.2 Mb gene block on chromosome 5 is also differentiated between coastal and inland S. townsendi. Genes in this block are associated with fatty acid oxidation and energy-related signaling transduction, thus linked to mitochondrial functions. Genetic variation within this candidate gene block covaries with mitochondrial DNA and shows signatures of divergent selection. Spatial variation in mitonuclear ancestries is correlated with climatic conditions. Together, these observations suggest divergent mitonuclear coadaptation underpins cryptic differentiation in this species complex.

Pervasive Mitonuclear Coadaptation Underlies Fast Development in Interpopulation Hybrids of a Marine Crustacean

Cellular energy production requires coordinated interactions between genetic components from the nuclear and mitochondrial genomes. This coordination results in coadaptation of interacting elements within populations. Interbreeding between divergent gene pools can disrupt coadapted loci and result in hybrid fitness breakdown. While specific incompatible loci have been detected in multiple eukaryotic taxa, the extent of the nuclear genome that is influenced by mitonuclear coadaptation is not clear in any species. Here, we used F₂ hybrids between two divergent populations of the copepod Tigriopus californicus to examine mitonuclear coadaptation across the nuclear genome. Using developmental rate as a measure of fitness, we found that fast-developing copepods had higher ATP synthesis capacity than slow developers, suggesting variation in developmental rates is at least partly associated with mitochondrial dysfunction. Using Pool-seq, we detected strong biases for maternal alleles across 7 (of 12) chromosomes in both reciprocal crosses in high-fitness hybrids, whereas low-fitness hybrids showed shifts toward the paternal population. Comparison with previous results on a different hybrid cross revealed largely different patterns of strong mitonuclear coadaptation associated with developmental rate. Our findings suggest that functional coadaptation between interacting nuclear and mitochondrial components is reflected in strong polygenic effects on this life-history phenotype, and reveal that molecular coadaptation follows independent evolutionary trajectories among isolated populations.

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.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

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.

Extreme mito-nuclear discordance in a peninsular lizard: the role of drift, selection, and climate

Nuclear and mitochondrial genomes coexist within cells but are subject to different tempos and modes of evolution. Evolutionary forces such as drift, mutation, selection, and migration are expected to play fundamental roles in the origin and maintenance of diverged populations; however, divergence may lag between genomes subject to different modes of inheritance and functional specialization. Herein, we explore whole mitochondrial genome data and thousands of nuclear single nucleotide polymorphisms to evidence extreme mito-nuclear discordance in the small black-tailed brush lizard, Urosaurus nigricaudus, of the Peninsula of Baja California, Mexico and southern California, USA, and discuss potential drivers. Results show three deeply divergent mitochondrial lineages dating back to the later Miocene (ca. 5.5 Ma) and Pliocene (ca. 2.8 Ma) that likely followed geographic isolation due to trans-peninsular seaways. This contrasts with very low levels of genetic differentiation in nuclear loci (F_ST < 0.028) between mtDNA lineages. Analyses of protein-coding genes reveal substantial fixed variation between mitochondrial lineages, of which a significant portion comes from non-synonymous mutations. A mixture of drift and selection is likely responsible for the rise of these mtDNA groups, albeit with little evidence of marked differences in climatic niche space between them. Finally, future investigations can look further into the role that mito-nuclear incompatibilities and mating systems play in explaining contrasting nuclear gene flow.

Investigating mitonuclear interactions in human admixed populations

To function properly, mitochondria utilize products of 37 mitochondrial and >1,000 nuclear genes, which should be compatible with each other. Discordance between mitochondrial and nuclear genetic ancestry could contribute to phenotypic variation in admixed populations. Here, we explored potential mitonuclear incompatibility in six admixed human populations from the Americas: African Americans, African Caribbeans, Colombians, Mexicans, Peruvians and Puerto Ricans. By comparing nuclear versus mitochondrial ancestry in these populations, we first show that mitochondrial DNA (mtDNA) copy number decreases with increasing discordance between nuclear and mtDNA ancestry. The direction of this effect is consistent across mtDNA haplogroups of different geographic origins. This observation indicates suboptimal regulation of mtDNA replication when its components are encoded by nuclear and mtDNA genes with different ancestry. Second, while most populations analysed exhibit no such trend, in African Americans and Puerto Ricans, we find a significant enrichment of ancestry at nuclear-encoded mitochondrial genes towards the source populations contributing the most prevalent mtDNA haplogroups (African and Native American, respectively). This possibly reflects compensatory effects of selection in recovering mitonuclear interactions optimized in the source populations. Our results provide evidence of mitonuclear interactions in human admixed populations and we discuss their implications for human health and disease.

A common pattern of DNase I footprinting throughout the human mtDNA unveils clues for a chromatin-like organization

Human mitochondrial DNA (mtDNA) is believed to lack chromatin and histones. Instead, it is coated solely by the transcription factor TFAM. We asked whether mtDNA packaging is more regulated than once thought. To address this, we analyzed DNase-seq experiments in 324 human cell types and found, for the first time, a pattern of 29 mtDNA Genomic footprinting (mt-DGF) sites shared by ∼90% of the samples. Their syntenic conservation in mouse DNase-seq experiments reflect selective constraints. Colocalization with known mtDNA regulatory elements, with G-quadruplex structures, in TFAM-poor sites (in HeLa cells) and with transcription pausing sites, suggest a functional regulatory role for such mt-DGFs. Altered mt-DGF pattern in interleukin 3-treated CD34⁺ cells, certain tissue differences, and significant prevalence change in fetal versus nonfetal samples, offer first clues to their physiological importance. Taken together, human mtDNA has a conserved protein-DNA organization, which is likely involved in mtDNA regulation.

The Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

The First Mitochondrial Genomics and Evolution SMBE-Satellite Meeting: A New Scientific Symbiosis

The central role of the mitochondrion for cellular and organismal metabolism is well known, yet its functional role in evolution has rarely been featured in leading international conferences. Moreover, the contribution of mitochondrial genetics to complex disease phenotypes is particularly important, and although major advances have been made in the field of genomics, mitochondrial genomic data have in many cases been overlooked. Accumulating data and new knowledge support a major contribution of this maternally inherited genome, and its interactions with the nucleus, to both major evolutionary processes and diverse disease phenotypes. These advances encouraged us to assemble the first Mitochondrial Genomics and Evolution (MGE) meeting-an SMBE satellite and Israeli Science foundation international conference (Israel, September 2017). Here, we report the content and outcome of the MGE meeting (https://www.mge2017.com/; last accessed November 5, 2017).

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.

Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only pathway with both nuclear and mitochondrial encoded proteins. The importance of the interactions between these two genomes has recently received more attention because of their potential evolutionary effects and how they may affect human health and disease. In many different organisms, healthy nuclear and mitochondrial genome hybrids between species or among distant populations within a species affect fitness and OxPhos functions. However, what is less understood is whether these interactions impact individuals within a single natural population. The significance of this impact depends on the strength of selection for mito-nuclear interactions. We examined whether mito-nuclear interactions alter allele frequencies for ~11,000 nuclear SNPs within a single, natural Fundulus heteroclitus population containing two divergent mitochondrial haplotypes (mt-haplotypes). Between the two mt-haplotypes, there are significant nuclear allele frequency differences for 349 SNPs with a p-value of 1% (236 with 10% FDR). Unlike the rest of the genome, these 349 outlier SNPs form two groups associated with each mt-haplotype, with a minority of individuals having mixed ancestry. We use this mixed ancestry in combination with mt-haplotype as a polygenic factor to explain a significant fraction of the individual OxPhos variation. These data suggest that mito-nuclear interactions affect cardiac OxPhos function. The 349 outlier SNPs occur in genes involved in regulating metabolic processes but are not directly associated with the 79 nuclear OxPhos proteins. Therefore, we postulate that the evolution of mito-nuclear interactions affects OxPhos function by acting upstream of OxPhos.

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

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

Mitochondrial-nuclear co-evolution leads to hybrid incompatibility through pentatricopeptide repeat proteins

Mitochondrial-nuclear incompatibility has a major role in reproductive isolation between species. However, the underlying mechanism and driving force of mitochondrial-nuclear incompatibility remain elusive. Here, we report a pentatricopeptide repeat-containing (PPR) protein, Ccm1, and its interacting partner, 15S rRNA, to be involved in hybrid incompatibility between two yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus S. bayanus-Ccm1 has reduced binding affinity for S. cerevisiae-15S rRNA, leading to respiratory defects in hybrid cells. This incompatibility can be rescued by single mutations on several individual PPR motifs, demonstrating the highly evolvable nature of PPR proteins. When we examined other PPR proteins in the closely related Saccharomyces sensu stricto yeasts, about two-thirds of them showed detectable incompatibility. Our results suggest that fast co-evolution between flexible PPR proteins and their mitochondrial RNA substrates may be a common driving force in the development of mitochondrial-nuclear hybrid incompatibility.

Divergent evolution of life span associated with mitochondrial DNA evolution

Mitochondria play a key role in ageing. The pursuit of genes that regulate variation in life span and ageing have shown that several nuclear-encoded mitochondrial genes are important. However, the role of mitochondrial encoded genes (mtDNA) is more controversial and our appreciation of the role of mtDNA for the evolution of life span is limited. We use replicated lines of seed beetles that have been artificially selected for long or short life for >190 generations, now showing dramatic phenotypic differences, to test for a possible role of mtDNA in the divergent evolution of ageing and life span. We show that these divergent selection regimes led to the evolution of significantly different mtDNA haplotype frequencies. Selection for a long life and late reproduction generated positive selection for one specific haplotype, which was fixed in most such lines. In contrast, selection for reproduction early in life led to both positive selection as well as negative frequency-dependent selection on two different haplotypes, which were both present in all such lines. Our findings suggest that the evolution of life span was in part mediated by mtDNA, providing support for the emerging general tenet that adaptive evolution of life-history syndromes may involve mtDNA.

Mitochondrial 16S rRNA Is Methylated by tRNA Methyltransferase TRMT61B in All Vertebrates

The mitochondrial ribosome, which translates all mitochondrial DNA (mtDNA)-encoded proteins, should be tightly regulated pre- and post-transcriptionally. Recently, we found RNA-DNA differences (RDDs) at human mitochondrial 16S (large) rRNA position 947 that were indicative of post-transcriptional modification. Here, we show that these 16S rRNA RDDs result from a 1-methyladenosine (m¹A) modification introduced by TRMT61B, thus being the first vertebrate methyltransferase that modifies both tRNA and rRNAs. m¹A947 is conserved in humans and all vertebrates having adenine at the corresponding mtDNA position (90% of vertebrates). However, this mtDNA base is a thymine in 10% of the vertebrates and a guanine in the 23S rRNA of 95% of bacteria, suggesting alternative evolutionary solutions. m¹A, uridine, or guanine may stabilize the local structure of mitochondrial and bacterial ribosomes. Experimental assessment of genome-edited Escherichia coli showed that unmodified adenine caused impaired protein synthesis and growth. Our findings revealed a conserved mechanism of rRNA modification that has been selected instead of DNA mutations to enable proper mitochondrial ribosome function.

Two solutions were selected during evolution to allow proper function of the vertebrate mitochondrial 16S ribosomal RNAeither RNA methylation by a tRNA methyltransferase or ancient evolutionary mutation.

RNA modifications constitute an important layer of information, with functional implications that are not written in the underlying DNA sequence. Recently, we observed an apparent RNA-DNA difference (RDD) at position 947 of the human mitochondrial 16S ribosomal RNA (rRNA), but its nature and mechanism were unclear. Here we show that this disparity reflects an m¹A modification (methylation at position 1 of the adenine moiety), and demonstrated by a combination of knock-down experiments in cells and in vitro methylation assays that the tRNA methyltransferase TRMT61B is the best candidate enzyme to introduce this modification. We also show that this modification is present in most of the 16S rRNA molecules in isolated mitochondrial ribosomes, and that it occurs in all vertebrates with an adenine (90% of the vertebrates), but not in those with a thymidine at this 16S rRNA position. Finally, as the first step towards understanding the functional importance of this rRNA modification, we used a genome-edited bacterial system to demonstrate that an unmodified adenine reduced the growth and translation rates of the bacteria as compared to both wild-type bacteria and mutant bacteria with a thymidine in the relevant position. Hence, three solutions were selected during evolution to allow proper function of the mitochondrial 16S rRNA—either RNA modification or two alternative ancient evolutionary DNA mutations.

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population-specific mitonuclear coadaptations create between-population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.