Fine-scale spatial and temporal genomic variation among Dungeness crab Cancer magister larval recruits in the California Current Ecosystem

Dynamic marine environments can shape complex spatial and temporal patterns in the population connectivity of marine species, and this is often exemplified in species with long larval phases. Here, we used a genotyping-by-sequencing (GBS) approach to examine fine-scale spatial and temporal genomic variation among Dungeness crab Cancer magister larval recruits sampled in the California Current Ecosystem. Specifically, we compared samples collected during expected- and late-season recruitment time periods within 2 consecutive years (2017 and 2018) at 2 sites in Oregon, USA (Yaquina Bay and Coos Bay). Evidence was found for high gene flow between the expected- and late-season recruits within each year and at both sites based on 1389 neutral loci. In contrast, strong genetic differentiation was observed between these 2 groups within each year and at both sites based on variation at 2 putatively adaptive loci. Contrary to prediction, the magnitude of genetic differentiation between these 2 seasonal groups was greater in 2017 when the Pacific Decadal Oscillation was stronger, upwelling was weaker, and the spring transition was later. Spatial genetic variation was not observed within 2017 or 2018. Comparing across years, expected- and late-season groups were differentiated at putatively adaptive loci. Interestingly, strong genetic differentiation was also observed between late-season groups across years. We found no evidence for cohesive larval dispersal among recruits based on genetic relatedness estimates. Overall, our findings provide evidence for high connectivity within Dungeness crab, but suggest that selective pressures and ocean conditions influence the genetic composition of larval recruits both intra- and inter-annually.