Strong habitat-specific phenotypic plasticity but no genome-wide differentiation across a rainforest gradient in an African butterfly.

Habitat-specific thermal responses are well documented in various organisms and likely determine the vulnerability of populations to climate change. However, the underlying roles of genetics and plasticity that shape such habitat-specific patterns are rarely investigated together. Here we examined the thermal plasticity of the butterfly Bicyclus dorothea originating from rainforest and ecotone habitats in Cameroon under common garden conditions. We also sampled wild-caught butterflies from forest and ecotone sites and used RADseq to explore genome-wide population differentiation. We found differences in the level of phenotypic plasticity across habitats. Specifically, ecotone populations exhibited greater sensitivity in wing eyespot features with variable development temperatures relative to rainforest populations. Known adaptive roles of wing eyespots in Bicyclus species suggest that this morphological plasticity is likely under divergent selection across environmental gradients. However, we found no distinct population structure of genome-wide variation between habitats, suggesting high level of ongoing gene flow between habitats is homogenizing most parts of the genome.