Understanding the spatial scale of local adaptation and the factors associated with adaptive diversity are important objectives for ecology and evolutionary biology, and have significant implications for effective conservation and management of wild populations and natural resources. In this study, we used an environmental association analysis (EAA) to identify important bioclimatic variables correlated with putatively adaptive genetic variation in a benthic marine invertebrate – the giant California sea cucumber (Parastichopus californicus) – spanning coastal British Columbia and southeastern Alaska. We used a redundancy analysis (RDA) with 3,699 SNPs obtained using RAD sequencing to detect candidate markers associated with 11 bioclimatic variables, including sea bottom and surface conditions, across two spatial scales (entire study area and within sub-regions). At the broadest scale, RDA revealed 59 candidate SNPs, 86% of which were associated with mean bottom temperature. Similar patterns were identified when population structure was accounted for. Additive polygenic scores, which provide a measure of the cumulative signal across all candidate SNPs, were strongly correlated with mean bottom temperature, consistent with spatially varying selection across a thermal gradient. At a finer scale, 23 candidate SNPs were detected, primarily associated with surface salinity (26%) and bottom current velocity (17%). Our findings suggest that environmental variables may play a role as drivers of spatially varying selection for P. californicus. These results provide context for future studies to evaluate the genetic basis of local adaptation in P. californicus and help inform the relevant scales and environmental variables for in situ field studies of putative adaptive variation in marine invertebrates.
Usage Notes:
01-RDA_candidate_snps
script to perform RDA for detecting candidate SNPs
02-compute_polygenic_scores
script to compute polygenic scores across candidate SNPs
3699snps_freqs
allele frequencies for RDA
candidates.frq
candidate frequencies for polygenic scores
meanbottomtemp_atsites
bottom temperature at sites for polygenic scores