In absence of local adaptation, plasticity and spatially varying selection rule: a view from genomic reaction norms in a panmictic species (Anguilla rostrata)

Abstract

Background: American eel (Anguilla rostrata) is one of the few species for which panmixia has been demonstrated at the scale of the entire species. As such, the development of long term local adaptation is impossible. However, both plasticity and spatially varying selection have been invoked in explaining how American eel may cope with an unusual broad scope of environmental conditions. Here, we address this question through transcriptomic analyses and genomic reaction norms of eels from two geographic origins reared in controlled environments.

Results: The null hypothesis of no difference in gene expression between eels from the two origins was rejected. Many unique transcripts and two out of seven gene clusters showed significant difference in expression, both at time of capture and after three months of common rearing. Differences in expression were observed at numerous genes representing many functional groups when comparing eels from a same origin reared under different salinity conditions. Plastic response to different rearing conditions varied among gene clusters with three clusters showing significant origin-environment interactions translating into differential genomic norms of reaction. Most genes and functional categories showing differences between origins were previously shown to be differentially expressed in a study comparing transcription profiles between adult European eels acclimated to different salinities.

Conclusions: These results emphasize that while plasticity in expression may be important, there is also a role for local genetic (and/or epigenetic) differences in explaining differences in gene expression between eels from different geographic origins. Such differences match those reported in genetically distinct populations in other fishes, both in terms of the proportion of genes that are differentially expressed and the diversity of biological functions involved. We thus propose that genetic differences between glass eels of different origins caused by spatially varying selection due to local environmental conditions translates into transcriptomic differences (including different genomic norms of reaction) which may in turn explain part of the phenotypic variance observed between different habitats colonized by eels.

Figure 1

Cluster mean fold change in Log-2 (and 95% confidence limits). From left to right: (a) between glass eels at the river mouth (T0: MR/GRB), (b) between elvers of both origins after three months rearing in BW and (c) in FW, (d) between environments for elvers from MR, and (e) GRB. GLM procedure was done, followed by K-means analysis where clusters with the same letter are not significantly different (α = 0.05).

Figure 2

Clusters genomic reaction norm, where MR and GRB absolute fold changes are illustrated on the same standardised scale. K-mean categories defined in Figure 1 are given for each cluster.