Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

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.

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.

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.

Mitochondrial genetic variation as a potential mediator of intraspecific behavioural diversity

Mitochondrial genes play an essential role in energy metabolism. Variation in the mitochondrial DNA (mtDNA) sequence often exists within species, and this variation can have consequences for energy production and organismal life history. Yet, despite potential links between energy metabolism and the expression of animal behaviour, mtDNA variation has been largely neglected to date in studies investigating intraspecific behavioural diversity. We outline how mtDNA variation and interactions between mitochondrial and nuclear genotypes may contribute to the expression of individual-to-individual behavioural differences within populations, and why such effects may lead to sex differences in behaviour. We contend that integration of the mitochondrial genome into behavioural ecology research may be key to fully understanding the evolutionary genetics of animal behaviour.

Mitonuclear interactions modulate nutritional preference

In nature, organisms are faced with constant nutritional options which fuel key life-history traits. Studies have shown that species can actively make nutritional decisions based on internal and external cues. Metabolism itself is underpinned by complex genomic interactions involving components from both nuclear and mitochondrial genomes. Products from these two genomes must coordinate how nutrients are extracted, used and recycled. Given the complicated nature of metabolism, it is not well understood how nutritional choices are affected by mitonuclear interactions. This is under the rationale that changes in genomic interactions will affect metabolic flux and change physiological requirements. To this end we used a large Drosophila mitonuclear genetic panel, comprising nine isogenic nuclear genomes coupled to nine mitochondrial haplotypes, giving a total of 81 different mitonuclear genotypes. We use a capillary-based feeding assay to screen this panel for dietary preference between carbohydrate and protein. We find significant mitonuclear interactions modulating nutritional choices, with these epistatic interactions also being dependent on sex. Our findings support the notion that complex genomic interactions can place a constraint on metabolic flux. This work gives us deeper insights into how key metabolic interactions can have broad implications on behaviour.

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.

Ecology, the pace-of-life, epistatic selection and the maintenance of genetic variation in life-history genes

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

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.

Mapping mitonuclear epistasis using a novel recombinant yeast population

Genetic variation in mitochondrial and nuclear genomes can perturb mitonuclear interactions and lead to phenotypic differences between individuals and populations. Despite their importance to most complex traits, it has been difficult to identify the interacting mitonuclear loci. Here, we present a novel advanced intercrossed population of Saccharomyces cerevisiae yeasts, called the Mitonuclear Recombinant Collection (MNRC), designed explicitly for detecting mitonuclear loci contributing to complex traits. For validation, we focused on mapping genes that contribute to the spontaneous loss of mitochondrial DNA (mtDNA) that leads to the petite phenotype in yeast. We found that rates of petite formation in natural populations are variable and influenced by genetic variation in nuclear DNA, mtDNA and mitonuclear interactions. We mapped nuclear and mitonuclear alleles contributing to mtDNA stability using the MNRC by integrating a term for mitonuclear epistasis into a genome-wide association model. We found that the associated mitonuclear loci play roles in mitotic growth most likely responding to retrograde signals from mitochondria, while the associated nuclear loci with main effects are involved in genome replication. We observed a positive correlation between growth rates and petite frequencies, suggesting a fitness tradeoff between mitotic growth and mtDNA stability. We also found that mtDNA stability was correlated with a mobile mitochondrial GC-cluster that is present in certain populations of yeast and that selection for nuclear alleles that stabilize mtDNA may be rapidly occurring. The MNRC provides a powerful tool for identifying mitonuclear interacting loci that will help us to better understand genotype-phenotype relationships and coevolutionary trajectories.

Population admixtures in medaka inferred by multiple arbitrary amplicon sequencing

Cost-effective genotyping can be achieved by sequencing PCR amplicons. Short 3-10 base primers can arbitrarily amplify thousands of loci using only a few primers. To improve the sequencing efficiency of the multiple arbitrary amplicon sequencing (MAAS) approach, we designed new primers and examined their efficiency in sequencing and genotyping. To demonstrate the effectiveness of our method, we applied it to examining the population structure of the small freshwater fish, medaka (Oryzias latipes). We obtained 2987 informative SNVs with no missing genotype calls for 67 individuals from 15 wild populations and three artificial strains. The estimated phylogenic and population genetic structures of the wild populations were consistent with previous studies, corroborating the accuracy of our genotyping method. We also attempted to reconstruct the genetic backgrounds of a commercial orange mutant strain, Himedaka, which has caused a genetic disturbance in wild populations. Our admixture analysis focusing on Himedaka showed that at least two wild populations had genetically been contributed to the nuclear genome of this mutant strain. Our genotyping methods and results will be useful in quantitative assessments of genetic disturbance by this commercially available strain.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer’s, Parkinson’s, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Interindividual variation in maximum aerobic metabolism varies with gill morphology and myocardial bioenergetics

This study asked whether interindividual variation in maximum and standard aerobic metabolic rates of the Gulf killifish, Fundulus grandis, correlate with gill morphology and cardiac mitochondrial bioenergetics, traits reflecting critical steps in the O2 transport cascade from the environment to the tissues. Maximum metabolic rate (MMR) was positively related to body mass, total gill filament length, and myocardial oxygen consumption during maximum oxidative phosphorylation (multiple R2=0.836). Standard metabolic rate (SMR) was positively related to body mass, total gill filament length, and myocardial oxygen consumption during maximum electron transport system activity (multiple R2=0.717). After controlling for body mass, individuals with longer gill filaments, summed over all gill arches, or greater cardiac respiratory capacity had higher whole-animal metabolic rates. The overall model fit and the explanatory power of individual predictor variables were better for MMR than for SMR, suggesting that gill morphology and myocardial bioenergetics are more important in determining active rather than resting metabolism. After accounting for body mass, heart ventricle mass was not related to variation in MMR or SMR, indicating that the quality of the heart (i.e., the capacity for mitochondrial metabolism) was more influential than heart size. Finally, the myocardial oxygen consumption required to offset the dissipation of the transmembrane proton gradient in the absence of ATP synthesis was not correlated with either MMR or SMR. The results support the idea that interindividual variation in aerobic metabolism, particularly maximum metabolic rate, is associated with variation in specific steps in the O2 transport cascade.

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

Exploring Thermal Sensitivities and Adaptations of Oxidative Phosphorylation Pathways

Temperature shifts are a major challenge to animals; they drive adaptations in organisms and species, and affect all physiological functions in ectothermic organisms. Understanding the origin and mechanisms of these adaptations is critical for determining whether ectothermic organisms will be able to survive when faced with global climate change. Mitochondrial oxidative phosphorylation is thought to be an important metabolic player in this regard, since the capacity of the mitochondria to produce energy greatly varies according to temperature. However, organism survival and fitness depend not only on how much energy is produced, but, more precisely, on how oxidative phosphorylation is affected and which step of the process dictates thermal sensitivity. These questions need to be addressed from a new perspective involving a complex view of mitochondrial oxidative phosphorylation and its related pathways. In this review, we examine the effect of temperature on the commonly measured pathways, but mainly focus on the potential impact of lesser-studied pathways and related steps, including the electron-transferring flavoprotein pathway, glycerophosphate dehydrogenase, dihydroorotate dehydrogenase, choline dehydrogenase, proline dehydrogenase, and sulfide:quinone oxidoreductase. Our objective is to reveal new avenues of research that can address the impact of temperature on oxidative phosphorylation in all its complexity to better portray the limitations and the potential adaptations of aerobic metabolism.

Advancing human disease research with fish evolutionary mutant models

Model organism research is essential to understand disease mechanisms. However, laboratory-induced genetic models can lack genetic variation and often fail to mimic the spectrum of disease severity. Evolutionary mutant models (EMMs) are species with evolved phenotypes that mimic human disease. EMMs complement traditional laboratory models by providing unique avenues to study gene-by-environment interactions, modular mutations in noncoding regions, and their evolved compensations. EMMs have improved our understanding of complex diseases, including cancer, diabetes, and aging, and illuminated mechanisms in many organs. Rapid advancements of sequencing and genome-editing technologies have catapulted the utility of EMMs, particularly in fish. Fish are the most diverse group of vertebrates, exhibiting a kaleidoscope of specialized phenotypes, many that would be pathogenic in humans but are adaptive in the species' specialized habitat. Importantly, evolved compensations can suggest avenues for novel disease therapies. This review summarizes current research using fish EMMs to advance our understanding of human disease.

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.

Molecular characterization of ubiquitin-conjugating enzyme gene ube2h and siRNA-mediated regulation on targeting p53 in turbot, Scophthalmus maximus

Ubiquitin-conjugating enzymes are key factors in the ubiquitin proteasome pathway (UPP), which play key roles in ubiquitination. These enzymes affect the efficiency of UPP during stress conditions. P53 has important control of cell cycle arrest and apoptosis in response to cellular stress; these modifications are critical for the stability and transcriptional activity of p53 as the protein activates downstream target genes that dictate the cellular response. However, few studies have investigated the effects of thermal stress in turbot (Scophthalmus maximus), specifically the UPP signaling pathway, and the crosstalk between the ube2h and p53. In this study, the rapid amplification of cDNA ends was used to obtain a full-length cDNA of the turbot UBE2H gene (Sm-ube2h) and perform bioinformatics analysis. Our results showed that the cDNA of the Sm-ube2h was 718 bp in length, encoding a 189 amino acid protein, with a theoretical isoelectric point of 4.77. It also contained a catalytic (UBCc) domain. Expression of Sm-ube2h in different tissues was detected and quantified by qPCR, which was highest in the spleen and lowest in the liver. We also investigated the Sm-ube2h expression profiles in the liver and heart after thermal stress, and changes in Sm-ube2h and p53 under thermal stress, upon RNA interference. Our data speculated that Sm-ube2h and p53 exhibited antagonistic effects under normal temperature conditions after ube2h interference, but displayed synergistic effects under thermal stress, suggesting the crosstalk between UPP and p53 signaling pathway. Our results improved our understanding of the underlying molecular mechanism of thermal tolerance in turbot.

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.

Interindividual plasticity in metabolic and thermal tolerance traits from populations subjected to recent anthropogenic heating

To better understand temperature's role in the interaction between local evolutionary adaptation and physiological plasticity, we investigated acclimation effects on metabolic performance and thermal tolerance among natural Fundulus heteroclitus (small estuarine fish) populations from different thermal environments. Fundulus heteroclitus populations experience large daily and seasonal temperature variations, as well as local mean temperature differences across their large geographical cline. In this study, we use three populations: one locally heated (32°C) by thermal effluence (TE) from the Oyster Creek Nuclear Generating Station, NJ, and two nearby reference populations that do not experience local heating (28°C). After acclimation to 12 or 28°C, we quantified whole-animal metabolic (WAM) rate, critical thermal maximum (CTmax) and substrate-specific cardiac metabolic rate (CaM, substrates: glucose, fatty acids, lactate plus ketones plus ethanol, and endogenous (i.e. no added substrates)) in approximately 160 individuals from these three populations. Populations showed few significant differences due to large interindividual variation within populations. In general, for WAM and CTmax, the interindividual variation in acclimation response (log2 ratio 28/12°C) was a function of performance at 12°C and order of acclimation (12-28°C versus 28-12°C). CTmax and WAM were greater at 28°C than 12°C, although WAM had a small change (2.32-fold) compared with the expectation for a 16°C increase in temperature (expect 3- to 4.4-fold). By contrast, for CaM, the rates when acclimatized and assayed at 12 or 28°C were nearly identical. The small differences in CaM between 12 and 28°C temperature were partially explained by cardiac remodeling where individuals acclimatized to 12°C had larger hearts than individuals acclimatized to 28°C. Correlation among physiological traits was dependent on acclimation temperature. For example, WAM was negatively correlated with CTmax at 12°C but positively correlated at 28°C. Additionally, glucose substrate supported higher CaM than fatty acid, and fatty acid supported higher CaM than lactate, ketones and alcohol (LKA) or endogenous. However, these responses were highly variable with some individuals using much more FA than glucose. These findings suggest interindividual variation in physiological responses to temperature acclimation and indicate that additional research investigating interindividual may be relevant for global climate change responses in many species.

Mitochondrial physiology and responses to elevated hydrogen sulphide in two isogenic lineages of an amphibious mangrove fish

Hydrogen sulphide (H2S) is toxic and can act as a selective pressure on aquatic organisms, facilitating a wide range of adaptations for life in sulphidic environments. Mangrove rivulus (Kryptolebias marmoratus) inhabit mangrove swamps and have developed high tolerance to environmental H2S. They are hermaphroditic and can self-fertilize, producing distinct isogenic lineages with different sensitivity to H2S. Here, we tested the hypothesis that observed differences in responses to H2S are the result of differences in mitochondrial functions. For this purpose, we performed two experimental series, testing (1) the overall mitochondrial oxidizing capacities and (2) the kinetics of apparent H2S mitochondrial oxidation and inhibition in two distinct lineages of mangrove rivulus, originally collected from Belize and Honduras. We used permeabilized livers from both lineages, measured mitochondrial oxidation, and monitored changes during gradual increases of sulphide. Ultimately, we determined that each lineage has a distinct strategy for coping with elevated H2S, indicating divergences in mitochondrial function and metabolism. The Honduras lineage has higher anaerobic capacity substantiated by higher lactate dehydrogenase activity and higher apparent H2S oxidation rates, likely enabling them to tolerate H2S by escaping aquatic H2S in a terrestrial environment. However, Belize fish have increased cytochrome c oxidase and citrate synthase activities as well as increased succinate contribution to mitochondrial respiration, allowing them to tolerate higher levels of aquatic H2S without inhibition of mitochondrial oxygen consumption. Our study reveals distinct physiological strategies in genetic lineages of a single species, indicating possible genetic and/or functional adaptations to sulphidic environments at the mitochondrial level.

The Role of Mitonuclear Incompatibility in Bipolar Disorder Susceptibility and Resilience Against Environmental Stressors

It has been postulated that mitochondrial dysfunction has a significant role in the underlying pathophysiology of bipolar disorder (BD). Mitochondrial functioning plays an important role in regulating synaptic transmission, brain function, and cognition. Neuronal activity is energy dependent and neurons are particularly sensitive to changes in bioenergetic fluctuations, suggesting that mitochondria regulate fundamental aspects of brain function. Vigorous evidence supports the role of mitochondrial dysfunction in the etiology of BD, including dysregulated oxidative phosphorylation, general decrease of energy, altered brain bioenergetics, co-morbidity with mitochondrial disorders, and association with genetic variants in mitochondrial DNA (mtDNA) or nuclear-encoded mitochondrial genes. Despite these advances, the underlying etiology of mitochondrial dysfunction in BD is unclear. A plausible evolutionary explanation is that mitochondrial-nuclear (mitonuclear) incompatibility leads to a desynchronization of machinery required for efficient electron transport and cellular energy production. Approximately 1,200 genes, encoded from both nuclear and mitochondrial genomes, are essential for mitochondrial function. Studies suggest that mitochondrial and nuclear genomes co-evolve, and the coordinated expression of these interacting gene products are essential for optimal organism function. Incompatibilities between mtDNA and nuclear-encoded mitochondrial genes results in inefficiency in electron flow down the respiratory chain, differential oxidative phosphorylation efficiency, increased release of free radicals, altered intracellular Ca²⁺ signaling, and reduction of catalytic sites and ATP production. This review explores the role of mitonuclear incompatibility in BD susceptibility and resilience against environmental stressors.

Comparative Mitogenomic Analysis Reveals Dynamics of Intron Within and Between Tricholoma Species and Phylogeny of Basidiomycota

The genus of Tricholoma is a group of important ectomycorrhizal fungi. The overlapping of morphological characteristics often leads to the confusion of Tricholoma species classification. In this study, the mitogenomes of five Tricholoma species were sequenced based on the next-generation sequencing technology, including T. matsutake SCYJ1, T. bakamatsutake, T. terreum, T. flavovirens, and T. saponaceum. These five mitogenomes were all composed of circular DNA molecules, with sizes ranging from 49,480 to 103,090 bp. Intergenic sequences were considered to be the main factor contributing to size variations of Tricholoma mitogenomes. Comparative mitogenomic analysis showed that the introns of the Agaricales mitogenome experienced frequent loss/gain events. In addition, potential gene transfer was detected between the mitochondrial and nuclear genomes of the five species of Tricholoma. Evolutionary analysis showed that the rps3 gene of the Tricholoma species was under positive selection or relaxed selection in the evolutionary process. In addition, large-scale gene rearrangements were detected between some Tricholoma species. Phylogenetic analysis using the Bayesian inference and maximum likelihood methods based on a combined mitochondrial gene set yielded identical and well-supported tree topologies. This study promoted the understanding of the genetics, evolution, and phylogeny of the Tricholoma genus and related species.

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747-Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA-sequencing, RAD-sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA-sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO-lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage-specific selection on mitochondria-related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.

Mitochondrial DNA lesions and copy number are strain dependent in endurance-trained mice

In this pilot work, we selected two inbred strains that respond well to endurance training (ET) (FVB/NJ, and SJL/J strains), and two strains that respond poorly (BALB/cByJ and NZW/LacJ), to determine the effect of a standardized ET treadmill program on mitochondrial and nuclear DNA (nucDNA) integrity, and mitochondrial DNA (mtDNA) copy number. DNA was isolated from plantaris muscles (n = 37) and a gene-specific quantitative PCR-based assay was used to measure DNA lesions and mtDNA copy number. Mean mtDNA lesions were not different within strains in the sedentary or exercise-trained states. However, mtDNA lesions were significantly higher in trained low-responding NZW/LacJ mice (0.24 ± 0.06 mtDNA lesions/10 Kb) compared to high-responding strains (mtDNA lesions/10 Kb: FVB/NJ = 0.11 ± 0.01, p = .049; SJL/J = 0.04 ± 0.02; p = .003). ET did not alter mean mtDNA copy numbers for any strain, although both sedentary and trained FVB/NJ mice had significantly higher mtDNA copies (99,890 ± 4,884 mtDNA copies) compared to low-responding strains (mtDNA copies: BALB/cByJ = 69,744 ± 4,675; NZW/LacJ = 65,687 ± 5,180; p < .001). ET did not change nucDNA lesions for any strain, however, SJL/J had the lowest mean nucDNA lesions (3.5 ± 0.14 nucDNA lesions/6.5 Kb) compared to all other strains (nucDNA lesions/6.5 Kb: FVB/NJ = 4.4 ± 0.11; BALB/cByJ = 4.7 ± 0.09; NZW/LacJ = 4.4 ± 0.11; p < .0001). Our results demonstrate strain differences in plantaris muscle mtDNA lesions in ET mice and, independent of condition, differences in mean mtDNA copy and nucDNA lesions between strains.

Adjustments of cardiac mitochondrial phenotype in a warmer thermal habitat is associated with oxidative stress in European perch, Perca fluviatilis

Mitochondria are playing key roles in setting the thermal limits of fish, but how these organelles participate in selection mechanisms during extreme thermal events associated with climate warming in natural populations is unclear. Here, we investigated the thermal effects on mitochondrial metabolism, oxidative stress, and mitochondrial gene expression in cardiac tissues of European perch (Perca fluviatilis) collected from an artificially heated ecosystem, the "Biotest enclosure", and an adjacent reference area in the Baltic sea with normal temperatures (~ 23 °C and ~ 16 °C, respectively, at the time of capture in summer). Fish were sampled one month after a heat wave that caused the Biotest temperatures to peak at ~ 31.5 °C, causing significant mortality. When assayed at 23 °C, Biotest perch maintained high mitochondrial capacities, while reference perch displayed depressed mitochondrial functions relative to measurements at 16 °C. Moreover, mitochondrial gene expression of nd4 (mitochondrial subunit of complex I) was higher in Biotest fish, likely explaining the increased respiration rates observed in this population. Nonetheless, cardiac tissue from Biotest perch displayed higher levels of oxidative damage, which may have resulted from their chronically warm habitat, as well as the extreme temperatures encountered during the preceding summer heat wave. We conclude that eurythermal fish such as perch are able to adjust and maintain mitochondrial capacities of highly aerobic organs such as the heart when exposed to a warming environment as predicted with climate change. However, this might come at the expense of exacerbated oxidative stress, potentially threatening performance in nature.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

Clues of in vivo nuclear gene regulation by mitochondrial short non-coding RNAs

Gene expression involves multiple processes, from transcription to translation to the mature, functional peptide, and it is regulated at multiple levels. Small RNA molecules are known to bind RNA messengers affecting their fate in the cytoplasm (a process generically termed ‘RNA interference’). Such small regulatory RNAs are well-known to be originated from the nuclear genome, while the role of mitochondrial genome in RNA interference was largely overlooked. However, evidence is growing that mitochondrial DNA does provide the cell a source of interfering RNAs. Small mitochondrial highly transcribed RNAs (smithRNAs) have been proposed to be transcribed from the mitochondrion and predicted to regulate nuclear genes. Here, for the first time, we show in vivo clues of the activity of two smithRNAs in the Manila clam, Ruditapes philippinarum. Moreover, we show that smithRNAs are present and can be annotated in representatives of the three main bilaterian lineages; in some cases, they were already described and assigned to a small RNA category (e.g., piRNAs) given their biogenesis, while in other cases their biogenesis remains unclear. If mitochondria may affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for them to interact with the nucleus and makes metazoan mitochondrial DNA a much more complex genome than previously thought.

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.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

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.

The Genomic Origins of Small Mitochondrial RNAs: Are They Transcribed by the Mitochondrial DNA or by Mitochondrial Pseudogenes within the Nucleus (NUMTs)?

Several studies have linked mitochondrial genetic variation to phenotypic modifications; albeit the identity of the mitochondrial polymorphisms involved remains elusive. The search for these polymorphisms led to the discovery of small noncoding RNAs, which appear to be transcribed by the mitochondrial DNA ("small mitochondrial RNAs"). This contention is, however, controversial because the nuclear genome of most animals harbors mitochondrial pseudogenes (NUMTs) of identical sequence to regions of mtDNA, which could alternatively represent the source of these RNAs. To discern the likely contributions of the mitochondrial and nuclear genome to transcribing these small mitochondrial RNAs, we leverage data from six vertebrate species exhibiting markedly different levels of NUMT sequence. We explore whether abundances of small mitochondrial RNAs are associated with levels of NUMT sequence across species, or differences in tissue-specific mtDNA content within species. Evidence for the former would support the hypothesis these RNAs are primarily transcribed by NUMT sequence, whereas evidence for the latter would provide strong evidence for the counter hypothesis that these RNAs are transcribed directly by the mtDNA. No association exists between the abundance of small mitochondrial RNAs and NUMT levels across species. Moreover, a sizable proportion of transcripts map exclusively to the mtDNA sequence, even in species with highest NUMT levels. Conversely, tissue-specific abundances of small mitochondrial RNAs are strongly associated with the mtDNA content. These results support the hypothesis that small mitochondrial RNAs are primarily transcribed by the mitochondrial genome and that this capacity is conserved across Amniota and, most likely, across most metazoan lineages.

Thermal tolerance and fish heart integrity: fatty acids profiles as predictors of species resilience

The cardiovascular system is a major limiting system in thermal adaptation, but the exact physiological mechanisms underlying responses to thermal stress are still not completely understood. Recent studies have uncovered the possible role of reactive oxygen species production rates of heart mitochondria in determining species' upper thermal limits. The present study examines the relationship between individual response to a thermal challenge test (CT_max), susceptibility to peroxidation of membrane lipids, heart fatty acid profiles and cardiac antioxidant enzyme activities in two salmonid species from different thermal habitats (Salvelinus alpinus, Salvelinus fontinalis) and their hybrids. The susceptibility to peroxidation of membranes in the heart was negatively correlated with individual thermal tolerance. The same relationship was found for arachidonic and eicosapentaenoic acid. Total H₂O₂ buffering activity of the heart muscle was higher for the group with high thermal resistance. These findings underline a potential general causative relationship between sensitivity to oxidative stress, specific fatty acids, antioxidant activity in the cardiac muscle and thermal tolerance in fish and likely other ectotherms. Heart fatty acid profile could be indicative of species resilience to global change, and more importantly the plasticity of this trait could predict the adaptability of fish species or populations to changes in environmental temperature.

Controlling for body size leads to inferential biases in the biological sciences

Many traits correlate with body size. Studies that seek to uncover the ecological factors that drive evolutionary responses in traits typically examine these responses relative to associated changes in body size using multiple regression analysis. However, it is not well appreciated that in the presence of strongly correlated variables, the partial (i.e., relative) regression coefficients often change sign compared to the original coefficients. Such sign reversals are difficult to interpret in a biologically meaningful way, and could lead to erroneous evolutionary inferences if the true mechanism underlying the sign reversal differed from the proposed mechanism. Here, we use simulations to demonstrate that sign reversal occurs over a wide range of parameter values common in the biological sciences. Further, as a case‐in‐point, we review the literature on brain size evolution; a field that explores how ecological traits relate to the evolution of relative brain size (brain size relative to body size). We find that most studies show sign reversals and thus that the inferences of many studies in this field may be inconclusive. Finally, we propose some approaches to mitigating this issue.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Mitochondria, sex and variation in routine metabolic rate

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (ṀO_{2 routine} ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with ṀO_{2 routine} , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with ṀO_{2 routine} involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in ṀO_{2 routine} . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on ṀO_{2 routine} are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

Differences in mitochondrial efficiency explain individual variation in growth performance

The physiological causes of intraspecific differences in fitness components such as growth rate are currently a source of debate. It has been suggested that differences in energy metabolism may drive variation in growth, but it remains unclear whether covariation between growth rates and energy metabolism is: (i) a result of certain individuals acquiring and consequently allocating more resources to growth, and/or is (ii) determined by variation in the efficiency with which those resources are transformed into growth. Studies of individually housed animals under standardized nutritional conditions can help shed light on this debate. Here we quantify individual variation in metabolic efficiency in terms of the amount of adenosine triphosphate (ATP) generated per molecule of oxygen consumed by liver and muscle mitochondria and examine its effects, both on the rate of protein synthesis within these tissues and on the rate of whole-body growth of individually fed juvenile brown trout (Salmo trutta) receiving either a high or low food ration. As expected, fish on the high ration on average gained more in body mass and protein content than those maintained on the low ration. Yet, growth performance varied more than 10-fold among individuals on the same ration, resulting in some fish on low rations growing faster than others on the high ration. This variation in growth for a given ration was related to individual differences in mitochondrial properties: a high whole-body growth performance was associated with high mitochondrial efficiency of ATP production in the liver. Our results show for the first time, to our knowledge, that among-individual variation in the efficiency with which substrates are converted into ATP can help explain marked variation in growth performance, independent of food intake. This study highlights the existence of inter-individual differences in mitochondrial efficiency and its potential importance in explaining intraspecific variation in whole-animal performance.

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.

Winter temperature correlates with mtDNA genetic structure of yellow-necked mouse population in NE Poland

We analysed a fragment (247 bp) of cytochrome b of mitochondrial DNA sequenced using 353 samples of yellow-necked mice Apodemus flavicollis trapped in seven forests and along three woodlot transects in north-eastern Poland. Our aims were to identify the phylogeographic pattern and mtDNA structure of the population and to evaluate the role of environmental conditions in shaping the spatial pattern of mtDNA diversity. We found out that three European haplogroups occurred sympatrically in north-eastern Poland. Inferences based on mtDNA haplotype distribution and frequency defined five subpopulations. The mtDNA-based structure of mice significantly correlated with winter temperature: frequency of Haplogroup 1 was positively, and that of Haplogroup 3 negatively correlated to mean temperature of January in the year of trapping. Synthesis of the published pan-European data on the species phylogeography also showed that the possibly 'thermophilous' Haplogroup 1 has the westernmost occurrence, whereas the more 'cold-resistant' Haplogroup 3 occurs much further to north-east than the other haplogroups. The observed patter may be a byproduct of the tight coevolution with nuclear genes, as we have earlier found that - in mice population in NE Poland - the spatial pattern of nuclear DNA was best explained by January temperature. Alternatively, the observed association of mitochondrial genetic variation with temperature is possible to be adaptive as cytochrome b is involved in the process of ATP production via oxidative phosphorylation.

Sexual conflict through mother's curse and father's curse

In contrast with autosomes, lineages of sex chromosomes reside for different amounts of time in males and females, and this transmission asymmetry makes them hotspots for sexual conflict. Similarly, the maternal inheritance of the mitochondrial genome (mtDNA) means that mutations that are beneficial in females can spread in a population even if they are deleterious in males, a form of sexual conflict known as Mother's Curse. While both Mother's Curse and sex chromosome induced sexual conflict have been well studied on their own, the interaction between mitochondrial genes and genes on sex chromosomes is poorly understood. Here, we use analytical models and computer simulations to perform a comprehensive examination of how transmission asymmetries of nuclear, mitochondrial, and sex chromosome-linked genes may both cause and resolve sexual conflicts. For example, the accumulation of male-biased Mother's Curse mtDNA mutations will lead to selection in males for compensatory nuclear modifier loci that alleviate the effect. We show how the Y chromosome, being strictly paternally transmitted provides a particularly safe harbor for such modifiers. This analytical framework also allows us to discover a novel kind of sexual conflict, by which Y chromosome-autosome epistasis may result in the spread of male beneficial but female deleterious mutations in a population. We christen this phenomenon Father's Curse. Extending this analytical framework to ZW sex chromosome systems, where males are the heterogametic sex, we also show how W-autosome epistasis can lead to a novel kind of nuclear Mother's Curse. Overall, this study provides a comprehensive framework to understand how genetic transmission asymmetries may both cause and resolve sexual conflicts.

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

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

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.

Genotype to phenotype: Diet-by-mitochondrial DNA haplotype interactions drive metabolic flexibility and organismal fitness

Diet may be modified seasonally or by biogeographic, demographic or cultural shifts. It can differentially influence mitochondrial bioenergetics, retrograde signalling to the nuclear genome, and anterograde signalling to mitochondria. All these interactions have the potential to alter the frequencies of mtDNA haplotypes (mitotypes) in nature and may impact human health. In a model laboratory system, we fed four diets varying in Protein: Carbohydrate (P:C) ratio (1:2, 1:4, 1:8 and 1:16 P:C) to four homoplasmic Drosophila melanogaster mitotypes (nuclear genome standardised) and assayed their frequency in population cages. When fed a high protein 1:2 P:C diet, the frequency of flies harbouring Alstonville mtDNA increased. In contrast, when fed the high carbohydrate 1:16 P:C food the incidence of flies harbouring Dahomey mtDNA increased. This result, driven by differences in larval development, was generalisable to the replacement of the laboratory diet with fruits having high and low P:C ratios, perturbation of the nuclear genome and changes to the microbiome. Structural modelling and cellular assays suggested a V161L mutation in the ND4 subunit of complex I of Dahomey mtDNA was mildly deleterious, reduced mitochondrial functions, increased oxidative stress and resulted in an increase in larval development time on the 1:2 P:C diet. The 1:16 P:C diet triggered a cascade of changes in both mitotypes. In Dahomey larvae, increased feeding fuelled increased β-oxidation and the partial bypass of the complex I mutation. Conversely, Alstonville larvae upregulated genes involved with oxidative phosphorylation, increased glycogen metabolism and they were more physically active. We hypothesise that the increased physical activity diverted energy from growth and cell division and thereby slowed development. These data further question the use of mtDNA as an assumed neutral marker in evolutionary and population genetic studies. Moreover, if humans respond similarly, we posit that individuals with specific mtDNA variations may differentially metabolise carbohydrates, which has implications for a variety of diseases including cardiovascular disease, obesity, and perhaps Parkinson's Disease.

Inference of epistatic effects in a key mitochondrial protein

We use Potts model inference to predict pair epistatic effects in a key mitochondrial protein-cytochrome c oxidase subunit 2-for ray-finned fishes. We examine the effect of phylogenetic correlations on our predictions using a simple exact fitness model, and we find that, although epistatic effects are underpredicted, they maintain a roughly linear relationship to their true (model) values. After accounting for this correction, epistatic effects in the protein are still relatively weak, leading to fitness valleys of depth 2Ns≃-5 in compensatory double mutants. Interestingly, positive epistasis is more pronounced than negative epistasis, and the strongest positive effects capture nearly all sites subject to positive selection in fishes, similar to virus proteins evolving under selection pressure in the context of drug therapy.

The Cytoplasm Affects the Epigenome in Drosophila melanogaster

Cytoplasmic components and their interactions with the nuclear genome may mediate patterns of phenotypic expression to form a joint inheritance system. However, proximate mechanisms underpinning these interactions remain elusive. To independently assess nuclear genetic and epigenetic cytoplasmic effects, we created a full-factorial design in which representative cytoplasms and nuclear backgrounds from each of two geographically disjunct populations of Drosophila melanogaster were matched together in all four possible combinations. To capture slowly-accumulating epimutations in addition to immediately occurring ones, these constructed populations were examined one year later. We found the K4 methylation of histone H3, H3K4me3, an epigenetic marker associated with transcription start-sites had diverged across different cytoplasms. The loci concerned mainly related to metabolism, mitochondrial function, and reproduction. We found little overlap (<8%) in sites that varied genetically and epigenetically, suggesting that epigenetic changes have diverged independently from any cis-regulatory sequence changes. These results are the first to show cytoplasm-specific effects on patterns of nuclear histone methylation. Our results highlight that experimental nuclear-cytoplasm mismatch may be used to provide a platform to identify epigenetic candidate loci to study the molecular mechanisms of cyto-nuclear interactions.