Mitochondrial genotype and phenotypic plasticity of gene expression in response to cold acclimation in killifish

Differential Mitochondrial Genome Expression of Three Sympatric Lizards in Response to Low-Temperature Stress

Ecological factors related to climate extremes have a significant influence on the adaptability of organisms, especially for ectotherms such as reptiles that are sensitive to temperature change. Climate extremes can seriously affect the survival and internal physiology of lizards, sometimes even resulting in the loss of local populations or even complete extinction. Indeed, studies have shown that the expression levels of the nuclear genes and mitochondrial genomes of reptiles change under low-temperature stress. At present, the temperature adaptability of reptiles has rarely been studied at the mitochondrial genome level. In the present study, the mitochondrial genomes of three species of lizards, Calotes versicolor, Ateuchosaurus chinensis, and Hemidactylus bowringii, which live in regions of sympatry, were sequenced. We used RT-qPCR to explore the level of mitochondrial gene expression under low-temperature stress, as compared to a control temperature. Among the 13 protein-coding genes (PCGs), the steady-state transcript levels of ND4L, ND1, ATP6, and COII were reduced to levels of 0.61 ± 0.06, 0.50 ± 0.08, 0.44 ± 0.16, and 0.41 ± 0.09 in C. versicolor, respectively, compared with controls. The transcript levels of the ND3 and ND6 genes fell to levels of just 0.72 ± 0.05 and 0.67 ± 0.05 in H. bowringii, compared with controls. However, the transcript levels of ND3, ND5, ND6, ATP6, ATP8, Cytb, and COIII genes increased to 1.97 ± 0.15, 2.94 ± 0.43, 1.66 ± 0.07, 1.59 ± 0.17, 1.46 ± 0.04, 1.70 ± 0.16, and 1.83 ± 0.07 in A. chinensis. Therefore, the differences in mitochondrial gene expression may be internally related to the adaptative strategy of the three species under low-temperature stress, indicating that low-temperature environments have a greater impact on A. chinensis, with a small distribution area. In extreme environments, the regulatory trend of mitochondrial gene expression in reptiles is associated with their ability to adapt to extreme climates, which means differential mitochondrial genome expression can be used to explore the response of different lizards in the same region to low temperatures. Our experiment aims to provide one new research method to evaluate the potential extinction of reptile species in warm winter climates.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Three steps down: Metabolic depression in winter-acclimatized crucian carp (Carassius carassius L.)

Acclimatization of certain ectothermic vertebrates to winter conditions is associated with reduced energy consumption (winter dormancy). Principally, this may be achieved by reducing movement activity, depression of basal cellular functions, or by switching from aerobic to anaerobic energy production to sustain low energy consumption during anoxia. Therefore, we determined standard (SMR), routine (RMR) and anoxic (AMR) metabolic rates in summer- (SumA; 20 °C) and winter-acclimatized (WinA; 2 °C) crucian carp (Carassius carassius), an anoxia-tolerant teleost fish. At 20 °C, RMR was 39% lower in WinA than SumA fish (p < 0.05), indicating reduced movements in winter. SMR, measured for the first time in crucian carp, was 45% lower in WinA than SumA fish at 10 °C (p < 0.05), indicating significant reduction of energy consumption in vital cellular processes. At 2 °C, AMR - measured from ethanol production- was 78% and 97% of the SMR and RMR levels at the same temperature, respectively. The current study revealed that the winter dormancy in anoxia-tolerant crucian carp is achieved in 3 different steps: (1) by active reduction of SMR, possibly in anticipation of seasonal anoxia, (2) acute Arrhenius Q10 effect that slows SMR as well as RMR, and (3) direct response to the absence of oxygen (AMR). Furthermore, the anoxic energy production was strongly dependent on body mass with scaling exponents of -0.335 and - 0.421 for WinA and SumA fish, respectively.

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.

Extensive alternative splicing triggered by mitonuclear mismatch in naturally introgressed Rhinolophus bats

Mitochondrial function needs strong interactions of mitochondrial and nuclear (mitonuclear) genomes, which can be disrupted by mitonuclear mismatch due to mitochondrial DNA (mtDNA) introgression between two formerly isolated populations or taxa. This mitonuclear disruption may cause severe cellular stress in mismatched individuals. Gene expression changes and alternative splicing (AS) are two important transcriptional regulations to respond to environmental or cellular stresses. We previously identified a naturally introgressed population in the intermediate horseshoe bat (Rhinolophus affinis). Individuals from this population belong to R. a. himalayanus and share almost identical nuclear genetic background; however, some of them had mtDNA from another subspecies (R. a. macrurus). With this unique natural system, we examined gene expression changes in six tissues between five mitonuclear mismatched and five matched individuals. A small number of differentially expressed genes (DEGs) were identified, and functional enrichment analysis revealed that most DEGs were related to immune response although some may be involved in response to oxidative stress. In contrast, we identified extensive AS events and alternatively spliced genes (ASGs) between mismatched and matched individuals. Functional enrichment analysis revealed that multiple ASGs were directly or indirectly associated with energy production in mitochondria which is vital for survival of organism. To our knowledge, this is the first study to examine the role of AS in responding to cellular stress caused by mitonuclear mismatch in natural populations. Our results suggest that AS may play a more important role than gene expression regulation in responding to severe environmental or cellular stresses.

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.

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.

Transcriptional Programs Underlying Cold Acclimation of Common Carp (Cyprinus carpio L.)

Properly regulated transcriptional responses to environmental perturbations are critical for the fitness of fish. Although gene expression profiles in the tissues of common carp upon cold stress were previously characterized, the transcriptional programs underlying cold acclimation are still not well known. In this study, the ability of three common carp strains including Hebao red carp (HB), Songpu mirror carp (SPM) and Yellow river carp (YR) to establish cold resistance after acclimation to a mild hypothermia stress at 18°C for 24 h was confirmed by measurements of the critical thermal minimums (CTMin). The gene expression profiles of the brain and the heart from these strains under both control and cold-acclimated conditions were characterized with RNA-sequencing. The data of the three common carp strains with different genetic background were combined in the differential gene expression analyses to balance the effects of genetic diversity on gene expression. Marked effects of tissue origins on the cold-induced transcriptional responses were revealed by comparing the differentially expressed gene (DEG) lists of the two tissues. Functional categories including spliceosome and RNA splicing were highly enriched in the DEGs of both tissues. However, steroid biosynthesis was specifically enriched in DEGs of the brain and response to unfolded protein was solely enriched in DEGs of the heart. Consistent with the up-regulation of the genes involved in cholesterol biosynthesis, total cholesterol content of the brain was significantly increased upon cold stress. Moreover, cold-induced alternative splicing (AS) events were explored and AS of the rbmx (RNA-binding motif protein, X chromosome) gene was confirmed by real-time quantitative PCR. Finally, a core set of cold responsive genes (CRGs) were defined by comparative transcriptomic analyses. Our data provide insights into the transcriptional programs underlying cold acclimation of common carp and offer valuable clues for further investigating the genetic determinants for cold resistance of farmed fish.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

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.

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.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Patterns of alternative splicing in response to cold acclimation in fish

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Sensory Ecology of Predator-Induced Phenotypic Plasticity

Ecological communities are organized in trophic levels that share manifold interactions forming complex food webs. Infochemicals can further modify these interactions, e.g., by inducing defenses in prey. The micro-crustacean Daphnia is able to respond to predator-specific chemical cues indicating an increased predation risk. Daphnia shows plastic responses by adapting its morphology, behavior, and physiology, increasing organism, and population fitness. This stabilizes community structures. This review will describe the progress that has been made in understanding the high degree of plasticity observed in the model crustacean Daphnia. I summarize current knowledge on the processes of predator detection, ranging from the nature of biologically active chemical cues to the underlying neurophysiological mechanisms. With this, I aim to provide a comprehensive overview on the molecular mechanisms of ad hoc environmental phenotypic adaptation. In times of climate change and pollution understanding information transfer in aquatic systems is valuable as it will allow us to predict whether and how community structures are being affected.

Gene flow prevents mitonuclear co-adaptation: A comparative portrait of sympatric wild types and cybrids in the fish Chrosomus eos

Allospecific mtDNA can occasionally be beneficial for the fitness of populations. It is, however, difficult to assess the effect of mtDNA in natural conditions due to genetic and/or environmental interactions. In the fish Chrosomus eos, the transfer of C. neogaeus mitochondria occurs in a single generation and results in natural cybrids. For a few lakes in Quebec, C. eos can harbor either a C. eos mtDNA (wild types) or a C. neogaeus mtDNA (cybrids). Moreover, mtDNA of cybrids originated either from Mississippian or Atlantic glacial refuges. Such diversity provides a useful system for in situ assessment of allospecific mtDNA effects. We determined genetic, epigenetic and transcriptomic variation as well as mitochondrial enzymatic activity (complex IV) changes among wild types and cybrids either in sympatry or allopatry. Wild types and cybrids did not segregate spatially within a lake. Moreover, no significant genetic differentiation was detected among wild types and cybrids indicating sustained gene flow. Mitochondrial complex IV activity was higher for cybrids in both sympatry and allopatry while no difference was detected among cybrid haplotypes. Epigenetic and transcriptomic analyses revealed only subtle differences between sympatric wild types and cybrids compared to differences between sites. Altogether, these results indicate a limited influence of allospecific mtDNA in nuclear gene expression when controlling for genetic and environmental effects. The absence of a reproductive barrier between wild types and cybrids results in random association of either C. eos or C. neogaeus mtDNA with C. eos nDNA at each generation, and prevents mitonuclear co-adaptation in sympatry.

Similarities in temperature-dependent gene expression plasticity across timescales in threespine stickleback (Gasterosteus aculeatus)

Phenotypic plasticity occurs at a variety of timescales, but little is known about the degree to which plastic responses at different timescales are associated with similar underlying molecular processes, which is critical for assessing the effects of plasticity on evolutionary trajectories. To address this issue, we identified differential gene expression in response to developmental temperature in the muscle transcriptome of adult threespine stickleback (Gasterosteus aculeatus) exposed to 12, 18 and 24°C until hatch and then held at 18°C for 9 months and compared these results to differential gene expression in response to adult thermal acclimation in stickleback developed at 18°C and then acclimated to 5 and 25°C as adults. Adult thermal acclimation affected the expression of 7,940 and 7,015 genes in response to cold and warm acclimation, respectively, and 4,851 of these genes responded in both treatments. In contrast, the expression of only 33 and 29 genes was affected by cold and warm development, respectively. The majority of the genes affected by developmental temperature were also affected by adult acclimation temperature. Many genes that were differentially expressed as a result of adult acclimation were associated with previously identified temperature-dependent effects on DNA methylation patterns, suggesting a role of epigenetic mechanisms in regulating gene expression plasticity during acclimation. Taken together, these results demonstrate similarities between the persistent effects of developmental plasticity on gene expression and the effects of adult thermal acclimation, emphasizing the potential for mechanistic links between plasticity acting at these different life stages.

Subspecies differences in thermal acclimation of mitochondrial function and the role of uncoupling proteins in killifish

Thermal effects on mitochondrial efficiency and ATP production can influence whole-animal thermal tolerance and performance. Thus, organisms may have the capacity to alter mitochondrial processes through acclimation or adaptation to mitigate these effects. One possible mechanism is through the action of uncoupling proteins (UCPs) which can decrease the proton motive force independent of the production of ATP. To test this hypothesis, we examined the mRNA expression patterns of UCP isoforms and characterized the effects of thermal acclimation and putative local thermal adaptation on mitochondrial capacity, proton leak, and P/O ratios in two subspecies of Atlantic killifish (Fundulus heteroclitus). Ucp1 was the dominant isoform in liver and was more highly expressed in northern killifish. We found that cold acclimation increased mitochondrial capacity (state III and maximum substrate oxidation capacity), state II membrane potential, proton leak, and P/O ratios in northern, but not southern killifish liver mitochondria. Palmitate-induced mitochondrial uncoupling was detected in northern, but not southern, killifish liver mitochondria, consistent with the differences in mRNA expression between the subspecies. Taken together, our data suggest that mitochondrial function is more plastic in response to thermal acclimation in northern killifish than southern killifish and that UCP1 may play a role in regulating the proton motive force in northern, but not southern killifish in response to thermal acclimation. These data demonstrate the potential for adaptive variation in mitochondrial plasticity in response to cold.