Whole genome sequencing reveals fine-scale climate associated adaptive divergence near the range limits of a temperate reef fish.

Adaptation to ocean climate is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus), is a temperate reef fish that displays both pelagic early life history stages and strong site-associated homing as adults; the species is also presently of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. Here we produce a chromosome-level genome assembly for cunner and characterize spatial population structure throughout Atlantic Canada using whole genome resequencing. The genome assembly spanned 0.72 Gbp and resolved 24 chromosomes; whole genome resequencing of 803 individuals from 20 locations spanning from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four distinct regional groups in Atlantic Canada, including three near the range edge in Newfoundland. Pairwise FST and selection scans revealed consistent signals of differentiation and selection at discrete genomic regions including adjacent peaks on chromosome 10 recurring across multiple pairwise comparisons (i.e., FST 0.5-0.75). Redundancy analysis suggested significant association of environmental variables related to benthic temperature and oxygen range with genomic structure, again highlighting the previously identified region on chromosome 10. Our results suggest that climate associated adaptation in this temperate reef fish drives regional diversity despite high early life history dispersal potential.