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

Environmental variation 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 with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise F_ST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. F_ST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.