Strong spatial genetic structure in a Baltic Sea herbivore due to recent range expansion, multiple bottlenecks and low connectivity

In the Baltic Sea, recent range expansions following the opening of the Danish straits have resulted in a low-diversity ecosystem, both among and within species. However, relatively little is known about population genetic patterns within the basin, except for in a few commercially caught species and some primary producers thought to be ecosystem engineers. Here, we investigate the population genetic structure of the ecologically important crustaceanIdotea balthicathroughout the Baltic Sea using an array of 33,774 genome-wide SNP markers derived from 2b-RAD sequencing. We also generate a biophysical connectivity matrix, with which we compare the genomic data. We find strong population structure on small scales across the Baltic Sea, and that genomic patterns in most cases closely match biophysical connectivity, suggesting that current patterns are important for dispersal of this species. We also find a strong signal of multiple bottlenecks during the initial range expansion, in the form of reduced heterozygosity along the historical expansion front. The lack of gene flow among sampling sites in the Baltic Sea environmental gradient potentiates local adaptation, while at the same time also increasing genetic drift in low-diversity areas.