ABSTRACT
The large geographic distribution of the eastern oyster, Crassostrea virginica, makes it an ideal species to test how populations have adapted to latitudinal gradients in temperature. Despite inhabiting distinct thermal regimes, populations of C. virginica near the species' southern and northern geographic range show no population differences in their physiological response to temperature. In this study, we used comparative transcriptomics to understand how oysters from either end of the species' range maintain enantiostasis across three acclimation temperatures (10, 20, and 30°C). With this approach, we identified genes that were differentially expressed in response to temperature between individuals of C. virginica collected from New Brunswick, Canada and Louisiana, USA. We observed a core set of genes whose expression responded to temperature in both populations, but also an even larger set of genes with expression patterns that were unique to each population. Intriguingly, the genes with population-specific responses to temperature had elevated FST and Ka/Ks ratios compared to the genome-wide average. In contrast, genes showing only a response to temperature were found to only have elevated FST values suggesting that divergent FST may be due to selection on linked regulatory regions rather than positive selection on protein coding regions. Taken together, our results suggest that, despite coarse-scale physiological similarities, natural selection has shaped divergent gene expression responses to temperature in geographically separated populations of this broadly eurythermal marine invertebrate.