Genomic separation of Salish Sea and Pacific outer coast populations of the keystone sea star Pisaster ochraceus

Abstract

Understanding dispersal and population structure in marine species is critical to forecasting ecological responses to climate change. The ochre sea star (Pisaster ochraceus), an intertidal keystone predator along the North American Pacific coast, has long been considered to have little spatial genetic structure due to its extended larval dispersal period and thus high potential for gene flow. Here we characterize the spatial genomic architecture of P. ochraceus using whole-genome sequencing (WGS) data from individuals spanning over 3,000km of coastline from Alaska to southern California. Our analyses of putatively neutral SNPs evince considerable mixing throughout the latitudinal range yet reveal substantial population structure between outer Pacific coast populations and those residing within the semi-enclosed estuary of the Salish Sea, suggesting restricted gene flow and demographic divergence between these regions. Genomic divergence is further supported by evidence of selection, with outlier loci highlighting extended regions of low diversity in the Salish Sea, consistent with recent selective sweeps and potential local adaptation to the distinct estuarine conditions. These findings are consistent with the role of oceanographic barriers or environmental heterogeneity in shaping population structure in P. ochraceus, challenging an earlier notion of range-wide homogeneity. By exploring neutral and adaptive genetic diversity in this keystone marine species, we provide valuable insights for its capacity to persist in a rapidly changing world.