Trans-Acting Genotypes Associated with mRNA Expression Affect Metabolic and Thermal Tolerance Traits.

Evolutionary processes driving physiological trait variation depend on the underlying genomic mechanisms. Evolution of these mechanisms depends on the genetic complexity (involving many genes) and how gene expression impacting the traits is converted to phenotype. Yet, genomic mechanisms that impact physiological traits are diverse and context dependent (e.g., vary by environment and tissues), making them difficult to discern. We examine the relationships between genotype, mRNA expression, and physiological traits to discern the genetic complexity and whether the gene expression affecting the physiological traits is primarily cis- or trans-acting. We use low-coverage whole genome sequencing and heart- or brain-specific mRNA expression to identify polymorphisms directly associated with physiological traits and expressed quantitative trait loci (eQTL) indirectly associated with variation in six temperature specific physiological traits (standard metabolic rate, thermal tolerance, and four substrate specific cardiac metabolic rates). Focusing on a select set of mRNAs belonging to co-expression modules that explain up to 82% of temperature specific traits, we identified hundreds of significant eQTL for mRNA whose expression affects physiological traits. Surprisingly, most eQTL (97.4% for heart and 96.7% for brain) were trans-acting. This could be due to higher effect size of trans- versus cis-acting eQTL for mRNAs that are central to co-expression modules. That is, we may have enhanced the identification of trans-acting factors by looking for single nucleotide polymorphisms associated with mRNAs in co-expression modules that broadly influence gene expression patterns. Overall, these data indicate that the genomic mechanism driving physiological variation across environments is driven by trans-acting heart- or brain-specific mRNA expression.

Convergence or redundancy: alternative views about the evolutionary genomics of character displacement.

An evolutionary debate contrasts the importance of genetic convergence versus genetic redundancy. In genetic convergence, the same adaptive trait evolves because of similar genetic changes. In genetic redundancy, the adaptive trait evolves using different genetic combinations, and populations might not share the same genetic changes. Here we address this debate by examining single nucleotide polymorphisms (SNPs) associated with the rapid evolution of character displacement in Anolis carolinensis populations inhabiting replicate islands with and without a competitor species (1Spp and 2Spp islands, respectively). We identify 215-outliers SNPs that have improbably large FST values, low nucleotide variation, greater linkage than expected and that are enriched for genes underlying animal movement. The pattern of SNP divergence between 1Spp and 2Spp populations supports both genetic convergence and genetic redundancy for character displacement. In support of genetic convergence: all 215-outliers SNPs are shared among at least three of the five 2Spp island populations, and 23% of outlier SNPS are shared among all five 2Spp island populations. In contrast, in support of genetic redundancy: many outlier SNPs only have meaningful allele frequency differences between 1Spp and 2Spp islands on a few 2Spp islands. That is, on at least one of the 2Spp islands, 77% of outlier SNPs have allele frequencies more similar to those on 1Spp islands than to those on 2Spp islands. Focusing on genetic convergence is scientifically rigorous because it relies on replication. Yet, this focus distracts from the possibility that there are multiple, redundant genetic solutions that enhance the rate and stability of adaptive change.

Feeling the heat: variation in thermal sensitivity within and among populations.

Physiology defines individual responses to global climate change and species distributions across environments. Physiological responses are driven by temperature on three time scales: acute, acclimatory and evolutionary. Acutely, passive temperature effects often dictate an expected 2-fold increase in metabolic processes for every 10°C change in temperature (Q10). Yet, these acute responses often are mitigated through acclimation within an individual or evolutionary adaptation within populations over time. Natural selection can influence both responses and often reduces interindividual variation towards an optimum. However, this interindividual physiological variation is not well characterized. Here, we quantified responses to a 16°C temperature difference in six physiological traits across nine thermally distinct Fundulus heteroclitus populations. These traits included whole-animal metabolism (WAM), critical thermal maximum (CTmax) and substrate-specific cardiac metabolism measured in approximately 350 individuals. These traits exhibited high variation among both individuals and populations. Thermal sensitivity (Q10) was determined, specifically as the acclimated Q10, in which individuals were both acclimated and assayed at each temperature. The interindividual variation in Q10 was unexpectedly large: ranging from 0.6 to 5.4 for WAM. Thus, with a 16°C difference, metabolic rates were unchanged in some individuals, while in others they were 15-fold higher. Furthermore, a significant portion of variation was related to habitat temperature. Warmer populations had a significantly lower Q10 for WAM and CTmax after acclimation. These data suggest that individual variation in thermal sensitivity reflects different physiological strategies to respond to temperature variation, providing many different adaptive responses to changing environments.

Sequencing Bait: Nuclear and Mitogenome Assembly of an Abundant Coastal Tropical and Subtropical Fish, Atherinomorus stipes.

Genetic data from nonmodel species can inform ecology and physiology, giving insight into a species' distribution and abundance as well as their responses to changing environments, all of which are important for species conservation and management. Moreover, reduced sequencing costs and improved long-read sequencing technology allows researchers to readily generate genomic resources for nonmodel species. Here, we apply Oxford Nanopore long-read sequencing and low-coverage (∼1x) whole genome short-read sequencing technology (Illumina) to assemble a genome and examine population genetics of an abundant tropical and subtropical fish, the hardhead silverside (Atherinomorus stipes). These fish are found in shallow coastal waters and are frequently included in ecological models because they serve as abundant prey for commercially and ecologically important species. Despite their importance in sub-tropical and tropical ecosystems, little is known about their population connectivity and genetic diversity. Our A. stipes genome assembly is about 1.2 Gb with comparable repetitive element content (∼47%), number of protein duplication events, and DNA methylation patterns to other teleost fish species. Among five sampled populations spanning 43 km of South Florida and the Florida Keys, we find little population structure suggesting high population connectivity.

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 (CTmax), 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.

Cardiac physiology and metabolic gene expression during late organogenesis among F. heteroclitus embryo families from crosses between pollution-sensitive and -resistant parents.

BACKGROUND: The teleost fish Fundulus heteroclitus inhabit estuaries heavily polluted with persistent and bioaccumulative chemicals. While embryos of parents from polluted sites are remarkably resistant to toxic sediment and develop normally, embryos of parents from relatively clean estuaries, when treated with polluted sediment extracts, are developmentally delayed, displaying deformities characteristic of pollution-induced embryotoxicity. To gain insight into parental effects on sensitive and resistant phenotypes during late organogenesis, we established sensitive, resistant, and crossed embryo families using five female and five male parents from relatively clean and predominantly PAH-polluted estuaries each, measured heart rates, and quantified individual embryo expression of 179 metabolic genes. RESULTS: Pollution-induced embryotoxicity manifested as morphological deformities, significant developmental delays, and altered cardiac physiology was evident among sensitive embryos resulting from crosses between females and males from relatively clean estuaries. Significantly different heart rates among several geographically unrelated populations of sensitive, resistant, and crossed embryo families during late organogenesis and pre-hatching suggest site-specific adaptive cardiac physiology phenotypes relative to pollution exposure. Metabolic gene expression patterns (32 genes, 17.9%, at p < 0.05; 11 genes, 6.1%, at p < 0.01) among the embryo families indicate maternal pollutant deposition in the eggs and parental effects on gene expression and metabolic alterations. CONCLUSION: Heart rate differences among sensitive, resistant, and crossed embryos is a reliable phenotype for further explorations of adaptive mechanisms. While metabolic gene expression patterns among embryo families are suggestive of parental effects on several differentially expressed genes, a definitive adaptive signature and metabolic cost of resistant phenotypes is unclear and shows unexpected sensitive-resistant crossed embryo expression profiles. Our study highlights physiological and metabolic gene expression differences during a critical embryonic stage among pollution sensitive, resistant, and crossed embryo families, which may contribute to underlying resistance mechanisms observed in natural F. heteroclitus populations living in heavily contaminated estuaries.

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.

PAH-pollution effects on sensitive and resistant embryos: Integrating structure and function with gene expression.

Polycyclic aromatic hydrocarbons (PAHs) are among the most widespread natural and anthropogenic pollutants, and some PAHs are proven developmental toxicants. We chemically characterized clean and heavily polluted sites and exposed fish embryos to PAH polluted sediment extracts during four critical developmental stages. Embryos were collected from Fundulus heteroclitus populations inhabiting the clean and heavily polluted Superfund estuary. Embryos of parents from the clean sites are sensitive to PAH pollutants while those of parents from the heavily polluted site are resistant. Chemical analysis of embryos suggests PAH accumulation and pollution-induced toxicity among sensitive embryos during development that ultimately kills all sensitive embryos before hatching, while remarkably, the resistant embryos develop normally. The adverse effects on sensitive embryos are manifested as developmental delays, reduced heart rates, and severe heart, liver, and kidney morphological abnormalities. Gene expression analysis of early somitogenesis, heartbeat initiation, late organogenesis, and pre-hatching developmental stages reveals genes whose expression significantly differs between sensitive and resistant embryo populations and helps to explain mechanisms of sensitivity and resistance to polluted environments during vertebrate animal development.

Polygenic Selection within a Single Generation Leads to Subtle Divergence among Ecological NichesINc.

Selection on standing genetic variation may be effective enough to allow for adaptation to distinct niche environments within a single generation. Minor allele frequency changes at multiple, redundant loci of small effect can produce remarkable phenotypic shifts. Yet, demonstrating rapid adaptation via polygenic selection in the wild remains challenging. Here we harness natural replicate populations that experience similar selection pressures and harbor high within-, yet negligible among-population genetic variation. Such populations can be found among the teleost Fundulus heteroclitus that inhabits marine estuaries characterized by high environmental heterogeneity. We identify 10,861 single nucleotide polymorphisms in F. heteroclitus that belong to a single, panmictic population yet reside in environmentally distinct niches (one coastal basin and three replicate tidal ponds). By sampling at two time points within a single generation, we quantify both allele frequency change within as well as spatial divergence among niche subpopulations. We observe few individually significant allele frequency changes yet find that the "number" of moderate changes exceeds the neutral expectation by 10-100%. We find allele frequency changes to be significantly concordant in both direction and magnitude among all niche subpopulations, suggestive of parallel selection. In addition, within-generation allele frequency changes generate subtle but significant divergence among niches, indicative of local adaptation. Although we cannot distinguish between selection and genotype-dependent migration as drivers of within-generation allele frequency changes, the trait/s determining fitness and/or migration likelihood appear to be polygenic. In heterogeneous environments, polygenic selection and polygenic, genotype-dependent migration offer conceivable mechanisms for within-generation, local adaptation to distinct niches.

Evolutionary Physiology and Genomics in the Highly Adaptable Killifish (Fundulus heteroclitus).

By investigating evolutionary adaptations that change physiological functions, we can enhance our understanding of how organisms work, the importance of physiological traits, and the genes that influence these traits. This approach of investigating the evolution of physiological adaptation has been used with the teleost fish Fundulus heteroclitus and has produced insights into (i) how protein polymorphisms enhance swimming and development; (ii) the role of equilibrium enzymes in modulating metabolic flux; (iii) how variation in DNA sequences and mRNA expression patterns mitigate changes in temperature, pollution, and salinity; and (iv) the importance of nuclear-mitochondrial genome interactions for energy metabolism. Fundulus heteroclitus provides so many examples of adaptive evolution because their local population sizes are large, they have significant standing genetic variation, and they experience large ranges of environmental conditions that enhance the likelihood that adaptive evolution will occur. Thus, F. heteroclitus research takes advantage of evolutionary changes associated with exposure to diverse environments, both across the North American Atlantic coast and within local habitats, to contrast neutral versus adaptive divergence. Based on evolutionary analyses contrasting neutral and adaptive evolution in F. heteroclitus populations, we conclude that adaptive evolution can occur readily and rapidly, at least in part because it depends on large amounts of standing genetic variation among many genes that can alter physiological traits. These observations of polygenic adaptation enhance our understanding of how evolution and physiological adaptation progresses, thus informing both biological and medical scientists about genotype-phenotype relationships. © 2020 American Physiological Society. Compr Physiol 10:637-671, 2020.

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 (HE 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.

Population genomics of rapid evolution in natural populations: polygenic selection in response to power station thermal effluents.

BACKGROUND: Examples of rapid evolution are common in nature but difficult to account for with the standard population genetic model of adaptation. Instead, selection from the standing genetic variation permits rapid adaptation via soft sweeps or polygenic adaptation. Empirical evidence of this process in nature is currently limited but accumulating. RESULTS: We provide genome-wide analyses of rapid evolution in Fundulus heteroclitus populations subjected to recently elevated temperatures due to coastal power station thermal effluents using 5449 SNPs across two effluent-affected and four reference populations. Bayesian and multivariate analyses of population genomic structure reveal a substantial portion of genetic variation that is most parsimoniously explained by selection at the site of thermal effluents. An FST outlier approach in conjunction with additional conservative requirements identify significant allele frequency differentiation that exceeds neutral expectations among exposed and closely related reference populations. Genomic variation patterns near these candidate loci reveal that individuals living near thermal effluents have rapidly evolved from the standing genetic variation through small allele frequency changes at many loci in a pattern consistent with polygenic selection on the standing genetic variation. CONCLUSIONS: While the ultimate trajectory of selection in these populations is unknown and we survey only a minority of genomic loci, our findings suggest that polygenic models of adaptation may play important roles in large, natural populations experiencing recent selection due to environmental changes that cause broad physiological impacts.

Stable genetic structure and connectivity in pollution-adapted and nearby pollution-sensitive populations of Fundulus heteroclitus.

Populations of the non-migratory estuarine fish Fundulus heteroclitus inhabiting the heavily polluted New Bedford Harbour (NBH) estuary have shown inherited tolerance to local pollutants introduced to their habitats in the past 100 years. Here we examine two questions: (i) Is there pollution-driven selection on the mitochondrial genome across a fine geographical scale? and (ii) What is the pattern of migration among sites spanning a strong pollution gradient? Whole mitochondrial genomes were analysed for 133 F. heteroclitus from seven nearby collection sites: four sites along the NBH pollution cline (approx. 5 km distance), which had pollution-adapted fish, as well as one site adjacent to the pollution cline and two relatively unpolluted sites about 30 km away, which had pollution-sensitive fish. Additionally, we used microsatellite analyses to quantify genetic variation over three F. heteroclitus generations in both pollution-adapted and sensitive individuals collected from two sites at two different time points (1999/2000 and 2007/2008). Our results show no evidence for a selective sweep of mtDNA in the polluted sites. Moreover, mtDNA analyses revealed that both pollution-adapted and sensitive populations harbour similar levels of genetic diversity. We observed a high level of non-synonymous mutations in the most polluted site. This is probably associated with a reduction in Ne and concomitant weakening of purifying selection, a demographic expansion following a pollution-related bottleneck or increased mutation rates. Our demographic analyses suggest that isolation by distance influences the distribution of mtDNA genetic variation between the pollution cline and the clean populations at broad spatial scales. At finer scales, population structure is patchy, and neither spatial distance, pollution concentration or pollution tolerance is a good predictor of mtDNA variation. Lastly, microsatellite analyses revealed stable population structure over the last decade.

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.

The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish.

Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor-based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.

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.

Heritable oxidative phosphorylation differences in a pollutant resistant Fundulus heteroclitus population.

Populations can adapt to stress including recent anthropogenic pollution. Our published data suggests heritable differences in hepatocyte oxidative phosphorylation (OxPhos) metabolism in field-caught killifish (Fundulus heteroclitus) from the highly polluted Elizabeth River, VA, USA, relative to fish from a nearby, relatively unpolluted reference site in King's Creek VA. Consistent with other studies showing that Elizabeth River killifish are resistant to some of the toxic effects of certain contaminants, OxPhos measurements in hepatocytes from field-caught King's Creek but not field-caught Elizabeth River killifish were altered by acute benzo [a] pyrene exposures. To more definitively test whether the enhanced OxPhos metabolism and toxicity resistance are heritable, we measured OxPhos metabolism in a laboratory-reared F3 generation from the Elizabeth River population versus a laboratory-reared F1 generation from the King's Creek population and compared these results to previous data from the field-caught fish. The F3 Elizabeth River fish compared to F1 King's Creek fish had significantly higher State 3 respiration (routine metabolism) and complex II activity, and significantly lower complex I activity. The consistently higher routine metabolism in the F3 and field-caught Elizabeth River fish versus F1 and field-caught King's Creek fish implies a heritable change in OxPhos function. The observation that LEAK, E-State, Complex I and Complex II were different in laboratory bred versus field-caught fish suggests that different physiological mechanisms produce the enhanced OxPhos differences. Finally, similar to field-caught Elizabeth River fish, acute benzo [a] pyrene exposure did not affect OxPhos function of the laboratory-reared F3 generation, supporting the heritability of the toxicity resistance. Overall, these results suggest that the Elizabeth River population has evolved genetic changes in physiological homeostasis that enhance routine metabolism, and we speculate that these genetic changes interact with environmental factors altering the physiological mechanisms (e.g., alter LEAK, Complex I, and electron transfer system capacity) used to achieve this enhanced 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.