Acclimation and acute temperature effects on population differences in oxidative phosphorylation

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.

An updated review of cold shock and cold stress in fish

Temperature is critical in regulating virtually all biological functions in fish. Low temperature stress (cold shock/stress) is an often-overlooked challenge that many fish face as a result of both natural events and anthropogenic activities. In this study, we present an updated review of the cold shock literature based on a comprehensive literature search, following an initial review on the subject by M.R. Donaldson and colleagues, published in a 2008 volume of this journal. We focus on how knowledge on cold shock and fish has evolved over the past decade, describing advances in the understanding of the generalized stress response in fish under cold stress, what metrics may be used to quantify cold stress and what knowledge gaps remain to be addressed in future research. We also describe the relevance of cold shock as it pertains to environmental managers, policymakers and industry professionals, including practical applications of cold shock. Although substantial progress has been made in addressing some of the knowledge gaps identified a decade ago, other topics (e.g., population-level effects and interactions between primary, secondary and tertiary stress responses) have received little or no attention despite their significance to fish biology and thermal stress. Approaches using combinations of primary, secondary and tertiary stress responses are crucial as a research priority to better understand the mechanisms underlying cold shock responses, from short-term physiological changes to individual- and population-level effects, thereby providing researchers with better means of quantifying cold shock in laboratory and field settings.

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.

Acclimation to warm temperatures has important implications for mitochondrial function in Atlantic salmon (Salmo salar)

In fish, the capacity of thermal acclimation to preserve cardiac mitochondrial function under future warming scenarios is important to understand given the central roles that cardiac energy metabolism and performance play in this taxa's thermal tolerance. We acclimated Atlantic salmon to 12 and 20°C (for >2 months), and investigated the effects of acute and chronic warming on cardiac mitochondrial respiration and reactive oxygen species (ROS) production (release rate) using high-resolution fluorespirometry. Further, we compared the sensitivity of mitochondrial respiration to nitric oxide (i.e. the NO IC50), and assessed the mitochondrial response to anoxia-reoxygenation (AR). Acute exposure to 20°C increased maximal mitochondrial respiration by ∼55%; however, the mitochondria's complex I respiratory control ratio was 17% lower and ROS production was increased by ≥60%. Acclimation to 20°C: (1) preserved mitochondrial coupling and aerobic capacity; (2) decreased the mitochondria's ROS production by ∼30%; (3) increased the mitochondria's NO IC50 by ∼23%; and (4) improved mitochondrial membrane integrity at 20°C. AR did not affect mitochondrial function at 12°C, but acute exposure to 20°C and AR depressed maximal mitochondrial respiration (by ∼9%) and coupling (by ∼16%) without impacting ROS production. Finally, warm acclimation did not improve the capacity of mitochondria to recover from AR, indicating that there was no 'cross-tolerance' between these challenges. Our findings provide compelling evidence that thermal plasticity of cardiac mitochondrial function contributes to the Atlantic salmon's capability to survive at ≥20°C for prolonged periods, but call into question whether this plasticity may allow them to withstand high temperatures when combined with other stressors.

Improved mitochondrial function in salmon (Salmo salar) following high temperature acclimation suggests that there are cracks in the proverbial 'ceiling'

Mitochondrial function can provide key insights into how fish will respond to climate change, due to its important role in heart performance, energy metabolism and oxidative stress. However, whether warm acclimation can maintain or improve the energetic status of the fish heart when exposed to short-term heat stress is not well understood. We acclimated Atlantic salmon, a highly aerobic eurythermal species, to 12 and 20 °C, then measured cardiac mitochondrial functionality and integrity at 20 °C and at 24, 26 and 28 °C (this species' critical thermal maximum ± 2 °C). Acclimation to 20 °C vs. 12 °C enhanced many aspects of mitochondrial respiratory capacity and efficiency up to 24 °C, and preserved outer mitochondrial membrane integrity up to 26 °C. Further, reactive oxygen species (ROS) production was dramatically decreased at all temperatures. These data suggest that salmon acclimated to 'normal' maximum summer temperatures are capable of surviving all but the most extreme ocean heat waves, and that there is no 'tradeoff' in heart mitochondrial function when Atlantic salmon are acclimated to high temperatures (i.e., increased oxidative phosphorylation does not result in heightened ROS production). This study suggests that fish species may show quite different acclimatory responses when exposed to prolonged high temperatures, and thus, susceptibility to climate warming.

What do warming waters mean for fish physiology and fisheries?

Environmental signals act primarily on physiological systems, which then influence higher-level functions such as movement patterns and population dynamics. Increases in average temperature and temperature variability associated with global climate change are likely to have strong effects on fish physiology and thereby on populations and fisheries. Here we review the principal mechanisms that transduce temperature signals and the physiological responses to those signals in fish. Temperature has a direct, thermodynamic effect on biochemical reaction rates. Nonetheless, plastic responses to longer-term thermal signals mean that fishes can modulate their acute thermal responses to compensate at least partially for thermodynamic effects. Energetics are particularly relevant for growth and movement, and therefore for fisheries, and temperature can have pronounced effects on energy metabolism. All energy (ATP) production is ultimately linked to mitochondria, and temperature has pronounced effects on mitochondrial efficiency and maximal capacities. Mitochondria are dependent on oxygen as the ultimate electron acceptor so that cardiovascular function and oxygen delivery link environmental inputs with energy metabolism. Growth efficiency, that is the conversion of food into tissue, changes with temperature, and there are indications that warmer water leads to decreased conversion efficiencies. Moreover, movement and migration of fish relies on muscle function, which is partially dependent on ATP production but also on intracellular calcium cycling within the myocyte. Neuroendocrine processes link environmental signals to regulated responses at the level of different tissues, including muscle. These physiological processes within individuals can scale up to population responses to climate change. A mechanistic understanding of thermal responses is essential to predict the vulnerability of species and populations to climate change.

Effects of mitonuclear combination and thermal acclimation on the energetic phenotype

Activity of the oxidative phosphorylation complexes rely on intimately associated subunits encoded by the mitochondrial and nuclear genomes. Given the key role of this system in adenosine triphosphate production, genes from both genomes must coevolve. A combination of northern redbelly dace (Chrosomus eos) or finescale dace (C. neogaeus) mitochondrial genome with a C. eos nuclear genome allows for a close examination of a naturally occurring disruption of mitonuclear coevolution. We, therefore, investigated the combined effect of mitonuclear genotypes, acclimation, and temperature on the activity of enzymes linked with the energy metabolism in a sympatric population of wild type and cybrid. As expected, the activity of the nuclear-encoded citrate synthase was only influenced by temperature while the cytochrome c oxidase (composed of nuclear and mitochondrial subunits from wild type and cybrid individuals) responded differently to temperature. This study provides clear evidence of the extent by which mitonuclear coadaptation could influence aerobic metabolism.

Powerhouses in the cold: mitochondrial function during thermal acclimation in montane mayflies

Mitochondria provide the vast majority of cellular energy available to eukaryotes. Therefore, adjustments in mitochondrial function through genetic changes in mitochondrial or nuclear-encoded genes might underlie environmental adaptation. Environmentally induced plasticity in mitochondrial function is also common, especially in response to thermal acclimation in aquatic systems. Here, we examined mitochondrial function in mayfly larvae (Baetis and Drunella spp.) from high and low elevation mountain streams during thermal acclimation to ecologically relevant temperatures. A multi-substrate titration protocol was used to evaluate different respiratory states in isolated mitochondria, along with cytochrome oxidase and citrate synthase activities. In general, maximal mitochondrial respiratory capacity and oxidative phosphorylation coupling efficiency decreased during acclimation to higher temperatures, suggesting montane insects may be especially vulnerable to rapid climate change. Consistent with predictions of the climate variability hypothesis, mitochondria from Baetis collected at a low elevation site with highly variable daily and seasonal temperatures exhibited greater thermal tolerance than Baetis from a high elevation site with comparatively stable temperatures. However, mitochondrial phenotypes were more resilient than whole-organism phenotypes in the face of thermal stress. These results highlight the complex relationships between mitochondrial and organismal genotypes, phenotypes and environmental adaptation. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

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

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

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.

Evolutionary Toxicogenomics of the Striped Killifish (Fundulus majalis) in the New Bedford Harbor (Massachusetts, USA)

In this paper, we used a Genotyping-by-Sequencing (GBS) approach to find and genotype more than 4000 genome-wide SNPs (Single Nucleotide Polymorphisms) from striped killifish exposed to a variety of polychlorinated biphenyls (PCBs) and other aromatic pollutants in New Bedford Harbor (NBH, Massachusetts, USA). The aims of this study were to identify the genetic consequences of exposure to aquatic pollutants and detect genes that may be under selection. Low genetic diversity (H_E and π) was found in the site exposed to the highest pollution level, but the pattern of genetic diversity did not match the pollution levels. Extensive connectivity was detected among sampling sites, which suggests that balanced gene flow may explain the lack of genetic variation in response to pollution levels. Tests for selection identified 539 candidate outliers, but many of the candidate outliers were not shared among tests. Differences among test results likely reflect different test assumptions and the complex pollutant mixture. Potentially, selectively important loci are associated with 151 SNPs, and enrichment analysis suggests a likely involvement of these genes with pollutants that occur in NBH. This result suggests that selective processes at genes targeted by pollutants may be occurring, even at a small geographical scale, and may allow the local striped killifish to resist the high pollution levels.

Patterns of mitochondrial membrane remodeling parallel functional adaptations to thermal stress

The effect of temperature on mitochondrial performance is thought to be partly due to its effect on mitochondrial membranes. Numerous studies have shown that thermal acclimation and adaptation can alter the amount of inner-mitochondrial membrane (IMM), but little is known about the capacity of organisms to modulate mitochondrial membrane composition. Using northern and southern subspecies of Atlantic killifish (Fundulus heteroclitus) that are locally adapted to different environmental temperatures, we assessed whether thermal acclimation altered liver mitochondrial respiratory capacity or the composition and amount of IMM. We measured changes in phospholipid headgroups and headgroup-specific fatty acid (FA) remodeling, and used respirometry to assess mitochondrial respiratory capacity. Acclimation to 5°C and 33°C altered mitochondrial respiratory capacity in both subspecies. Northern F. heteroclitus exhibited greater mitochondrial respiratory capacity across acclimation temperatures, consistent with previously observed subspecies differences in whole-organism aerobic metabolism. Mitochondrial phospholipids were altered following thermal acclimation, and the direction of these changes was largely consistent between subspecies. These effects were primarily driven by remodeling of specific phospholipid classes and were associated with shifts in metabolic phenotypes. There were also differences in membrane composition between subspecies that were driven largely by differences in phospholipid classes. Changes in respiratory capacity between subspecies and with acclimation were largely but not completely accounted for by alterations in the amount of IMM. Taken together, these results support a role for changes in liver mitochondrial function in the ectothermic response to thermal stress during both acclimation and adaptation, and implicate lipid remodeling as a mechanism contributing to these changes.

Swim Training Modulates Skeletal Muscle Energy Metabolism, Oxidative Stress, and Mitochondrial Cholesterol Content in Amyotrophic Lateral Sclerosis Mice

Recently, in terms of amyotrophic lateral sclerosis (ALS), much attention has been paid to the cell structures formed by the mitochondria and the endoplasmic reticulum membranes (MAMs) that are involved in the regulation of Ca²⁺ signaling, mitochondrial bioenergetics, apoptosis, and oxidative stress. We assumed that remodeling of these structures via swim training may accompany the prolongation of the ALS lifespan. In the present study, we used transgenic mice with the G93A hmSOD1 gene mutation. We examined muscle energy metabolism, oxidative stress parameters, and markers of MAMs (Caveolin-1 protein level and cholesterol content in crude mitochondrial fraction) in groups of mice divided according to disease progression and training status. The progression of ALS was related to the lowering of Caveolin-1 protein levels and the accumulation of cholesterol in a crude mitochondrial fraction. These changes were associated with aerobic and anaerobic energy metabolism dysfunction and higher oxidative stress. Our data indicated that swim training prolonged the lifespan of ALS mice with accompanying changes in MAM components. Swim training also maintained mitochondrial function and lowered oxidative stress. These data suggest that modification of MAMs might play a crucial role in the exercise-induced deceleration of ALS development.

Intraspecific variation and plasticity in mitochondrial oxygen binding affinity as a response to environmental temperature

Mitochondrial function has been suggested to underlie constraints on whole-organism aerobic performance and associated hypoxia and thermal tolerance limits, but most studies have focused on measures of maximum mitochondrial capacity. Here we investigated whether variation in mitochondrial oxygen kinetics could contribute to local adaptation and plasticity in response to temperature using two subspecies of the Atlantic killifish (Fundulus heteroclitus) acclimated to a range of temperatures (5, 15, and 33 °C). The southern subspecies of F. heteroclitus, which has superior thermal and hypoxia tolerances compared to the northern subspecies, exhibited lower mitochondrial O₂ P50 (higher O₂ affinity). Acclimation to thermal extremes (5 or 33 °C) altered mitochondrial O₂ P50 in both subspecies consistent with the effects of thermal acclimation on whole-organism thermal tolerance limits. We also examined differences between subspecies and thermal acclimation effects on whole-blood Hb O₂-P50 to assess whether variation in oxygen delivery is involved in these responses. In contrast to the clear differences between subspecies in mitochondrial O₂-P50 there were no differences in whole-blood Hb-O₂ P50 between subspecies. Taken together these findings support a general role for mitochondrial oxygen kinetics in differentiating whole-organism aerobic performance and thus in influencing species responses to environmental change.

Remodeling pathway control of mitochondrial respiratory capacity by temperature in mouse heart: electron flow through the Q-junction in permeabilized fibers

Fuel substrate supply and oxidative phosphorylation are key determinants of muscle performance. Numerous studies of mammalian mitochondria are carried out (i) with substrate supply that limits electron flow, and (ii) far below physiological temperature. To analyze potentially implicated biases, we studied mitochondrial respiratory control in permeabilized mouse myocardial fibers using high-resolution respirometry. The capacity of oxidative phosphorylation at 37 °C was nearly two-fold higher when fueled by physiological substrate combinations reconstituting tricarboxylic acid cycle function, compared with electron flow measured separately through NADH to Complex I or succinate to Complex II. The relative contribution of the NADH pathway to physiological respiratory capacity increased with a decrease in temperature from 37 to 25 °C. The apparent excess capacity of cytochrome c oxidase above physiological pathway capacity increased sharply under hypothermia due to limitation by NADH-linked dehydrogenases. This mechanism of mitochondrial respiratory control in the hypothermic mammalian heart is comparable to the pattern in ectotherm species, pointing towards NADH-linked mt-matrix dehydrogenases and the phosphorylation system rather than electron transfer complexes as the primary drivers of thermal sensitivity at low temperature. Delineating the link between stress and remodeling of oxidative phosphorylation is important for understanding metabolic perturbations in disease evolution and cardiac protection.

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.

Thermal acclimation and subspecies-specific effects on heart and brain mitochondrial performance in a eurythermal teleost (Fundulus heteroclitus)

Mitochondrial performance may play a role in setting whole-animal thermal tolerance limits and their plasticity, but the relative roles of adjustments in mitochondrial performance across different highly aerobic tissues remain poorly understood. We compared heart and brain mitochondrial responses to acute thermal challenges and to thermal acclimation using high-resolution respirometry in two locally adapted subspecies of Atlantic killifish (Fundulus heteroclitus). We predicted that 5°C acclimation to would result in compensatory increases in mitochondrial performance, while 33°C acclimation would cause suppression of mitochondrial function to minimize the effects of high temperature on mitochondrial metabolism. In contrast, acclimation to both 33 and 5°C decreased mitochondrial performance compared to fish acclimated to 15°C. These adjustments could represent an energetic cost saving mechanism at temperature extremes. Acclimation responses were similar in both heart and brain; however, this effect was smaller in the heart which might indicate its importance in maintaining whole-animal thermal performance. Alternatively, larger acclimation effects in the brain might indicate greater thermal sensitivity compared to the heart. We detected only modest differences between subspecies that were dependent on the tissue assayed. These data demonstrate extensive plasticity in mitochondrial performance following thermal acclimation in killifish, and indicate that the extent of these responses differs between tissues, highlighting the importance and complexity of mitochondrial regulation in thermal acclimation in eurytherms.

Metabolic and regulatory responses involved in cold acclimation in Atlantic killifish, Fundulus heteroclitus

Ectotherms often respond to prolonged cold exposure by increasing mitochondrial capacity via elevated mitochondrial volume density [V _V(mit,f)]. In fish, higher V _V(mit,f) is typically associated with increased expression of nuclear respiratory factor 1 (Nrf1), a transcription factor that induces expression of nuclear-encoded respiratory genes. To examine if nrf1 expression or the expression of other genes that regulate mitochondrial biogenesis contribute to changes in whole-organism metabolic rate during cold acclimation, we examined the time course of changes in the expression of these genes and in metabolic rate in Atlantic killifish, Fundulus heteroclitus. Cold acclimation rapidly decreased metabolic rate, but increased the expression of nrf1 more gradually, with a time course that depended on how rapidly the fish were transitioned to low temperature. Cold-induced nrf1 expression was not associated with increases in biochemical indicators of mitochondrial respiratory capacity, suggesting that cold-induced mitochondrial biogenesis may occur without increases in oxidative capacity in this species. These observations imply that changes in nrf1 expression and metabolic rate due to cold acclimation occur through different physiological mechanisms, and that increases in V _V(mit,f) are likely not directly related to changes in metabolic rate with cold acclimation in this species. However, nrf1 expression differed between northern and southern killifish subspecies regardless of acclimation temperature, consistent with observed differences in metabolic rate and V _V(mit,f) at 5 °C between these subspecies. Taken together, these results reveal substantial complexity in the regulation of V _V(mit,f) and mitochondrial capacity with temperature in fish and the relationship of these parameters to metabolic rate.

A Cost-Effective Approach to Sequence Hundreds of Complete Mitochondrial Genomes

We present a cost-effective approach to sequence whole mitochondrial genomes for hundreds of individuals. Our approach uses small reaction volumes and unmodified (non-phosphorylated) barcoded adaptors to minimize reagent costs. We demonstrate our approach by sequencing 383 Fundulus sp. mitochondrial genomes (192 F. heteroclitus and 191 F. majalis). Prior to sequencing, we amplified the mitochondrial genomes using 4–5 custom-made, overlapping primer pairs, and sequencing was performed on an Illumina HiSeq 2500 platform. After removing low quality and short sequences, 2.9 million and 2.8 million reads were generated for F. heteroclitus and F. majalis respectively. Individual genomes were assembled for each species by mapping barcoded reads to a reference genome. For F. majalis, the reference genome was built de novo. On average, individual consensus sequences had high coverage: 61-fold for F. heteroclitus and 57-fold for F. majalis. The approach discussed in this paper is optimized for sequencing mitochondrial genomes on an Illumina platform. However, with the proper modifications, this approach could be easily applied to other small genomes and sequencing platforms.

Steep, coincident, and concordant clines in mitochondrial and nuclear-encoded genes in a hybrid zone between subspecies of Atlantic killifish, Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus-specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non-neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non-neutral nuclear clines, SNPs in two nuclear-encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin-binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non-unimodal, suggesting selection against advanced-generation hybrids, possibly due to effects on processes involved in oxidative metabolism.

Gene by environmental interactions affecting oxidative phosphorylation and thermal sensitivity

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only metabolic pathway with proteins encoded by both nuclear and mitochondrial genomes. In studies examining mitonuclear interactions among distant populations within a species or across species, the interactions between these two genomes can affect metabolism, growth, and fitness, depending on the environment. However, there is little data on whether these interactions impact natural populations within a single species. In an admixed Fundulus heteroclitus population with northern and southern mitochondrial haplotypes, there are significant differences in allele frequencies associated with mitochondrial haplotype. In this study, we investigate how mitochondrial haplotype and any associated nuclear differences affect six OxPhos parameters within a population. The data demonstrate significant OxPhos functional differences between the two mitochondrial genotypes. These differences are most apparent when individuals are acclimated to high temperatures with the southern mitochondrial genotype having a large acute response and the northern mitochondrial genotype having little, if any acute response. Furthermore, acute temperature effects and the relative contribution of Complex I and II depend on acclimation temperature: when individuals are acclimated to 12°C, the relative contribution of Complex I increases with higher acute temperatures, whereas at 28°C acclimation, the relative contribution of Complex I is unaffected by acute temperature change. These data demonstrate a complex gene by environmental interaction affecting the OxPhos pathway.