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.

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.

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.

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.

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.

Ecological population genomics in the marine environment.

Marine species live in a wide diversity of environments and yet, because of their pelagic life stages, are thought to be well-connected: they have high migration rates that inhibit significant population structure. Recent innovations in sequencing technologies now provide information on nucleotide polymorphisms at thousands to tens of thousands of loci based on whole genomes, reduced representative portions of genomes (0.1-1%) or a majority of expressed mRNAs. Data from these genomic approaches are used to define and quantify single-nucleotide polymorphisms (SNPs). These SNP data tend to agree with data from older technologies (allozymes or microsatellites), which support well-connected populations with few genetic differences among populations. However, these studies also find few percentages of SNPs (1-5%) that readily distinguish genetic differences among populations on relatively small geographic scales. The magnitudes of the genetic differences (FST values) suggest that hundreds of loci with significant differences are due to positive selective pressures. Thus, these data suggest that natural selection is effectively altering allele frequencies at 100s of loci in marine populations. In this manuscript, we provide examples of these studies, the strengths and weaknesses of different genomic approaches as well as important technical aspects associated with genomic approaches.

Acclimation and acute temperature effects on population differences in oxidative phosphorylation.

Temperature changes affect metabolism on acute, acclamatory, and evolutionary time scales. To better understand temperature's affect on metabolism at these different time scales, we quantified cardiac oxidative phosphorylation (OxPhos) in three Fundulus taxa acclimated to 12 and 28°C and measured at three acute temperatures (12, 20, and 28°C). The Fundulus taxa (northern Maine and southern Georgia F. heteroclitus, and a sister taxa, F. grandis) were used to identify evolved changes in OxPhos. Cardiac OxPhos metabolism was quantified by measuring six traits: state 3 (ADP and substrate-dependent mitochondrial respiration); E state (uncoupled mitochondrial activity); complex I, II, and IV activities; and LEAK ratio. Acute temperature affected all OxPhos traits. Acclimation only significantly affected state 3 and LEAK ratio. Populations were significantly different for state 3. In addition to direct effects, there were significant interactions between acclimation and population for complex I and between population and acute temperature for state 3. Further analyses suggest that acclimation alters the acute temperature response for state 3, E state, and complexes I and II: at the low acclimation temperature, the acute response was dampened at low assay temperatures, and at the high acclimation temperature, the acute response was dampened at high assay temperatures. Closer examination of the data also suggests that differences in state 3 respiration and complex I activity between populations were greatest between fish acclimated to low temperatures when assayed at high temperatures, suggesting that differences between the populations become more apparent at the edges of their thermal range.

Effects of Anthropogenic Pollution on the Oxidative Phosphorylation Pathway of Hepatocytes from Natural Populations of Fundulus heteroclitus.

Persistent organic pollutants (POPs), including polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs), potentially target mitochondria and cause toxicity. We compared the effects of POPs on mitochondrial respiration by measuring oxidative phosphorylation (OxPhos) metabolism in hepatocytes isolated from lab-depurated Fundulus heteroclitus from a Superfund site contaminated with PAHs (Elizabeth River VA, USA) relative to OxPhos metabolism in individuals from a relatively clean, reference population (King's Creek VA, USA). In individuals from the polluted Elizabeth River population, OxPhos metabolism displayed lower LEAK and lower activities in complex III, complex IV, and E State, but higher activity in complex I compared to individuals from the reference King's Creek population. To test the supposition that these differences were due to or related to the chronic PAH contamination history of the Elizabeth River population, we compared the OxPhos functions of undosed individuals from the polluted and reference populations to individuals from these populations dosed with a PAH {benzo [α] pyrene (BaP)} or a PCB {PCB126 (3,3',4,4',5-pentachlorobiphenyl)}, respectively. Exposure to PAH or PCB affected OxPhos in the reference King's Creek population but had no detectable effects on the polluted Elizabeth River population. Thus, PAH exposure significantly increased LEAK, and exposure to PCB126 significantly decreased State 3, E state and complex I activity in the reference King's Creek population. These data strongly implicate an evolved tolerance in the Elizabeth River fish where dosed fish are not affected by PAH exposure and undosed fish show decreased LEAK and increased State 3 and E state.

Phenotypic plasticity in gene expression contributes to divergence of locally adapted populations of Fundulus heteroclitus.

We examine the interaction between phenotypic plasticity and evolutionary adaptation using muscle gene expression levels among populations of the fish Fundulus heteroclitus acclimated to three temperatures. Our analysis reveals shared patterns of phenotypic plasticity due to thermal acclimation as well as non-neutral patterns of variation among populations adapted to different thermal environments. For the majority of significant differences in gene expression levels, phenotypic plasticity and adaptation operate on different suites of genes. The subset of genes that demonstrate both adaptive differences and phenotypic plasticity, however, exhibit countergradient variation of expression. Thus, expression differences among populations counteract environmental effects, reducing the phenotypic differentiation between populations. Finally, gene-by-environment interactions among genes with non-neutral patterns of expression suggest that the penetrance of adaptive variation depends on the environmental conditions experienced by the individual.