Distribution and population structure in the naked goby Gobiosoma bosc (Perciformes: Gobiidae) along a salinity gradient in two western Atlantic estuaries

Many species of fish produce larvae that undergo a prolonged dispersal phase. However, evidence from a number of recent studies on demersal fishes suggests that the dispersal of propagules may not be strongly correlated with gene flow. Instead, other factors like larval behavior and the availability of preferred settlement habitat may be more important to maintaining population structure. We used an ecologically important benthic fish species, Gobiosoma bosc (naked goby), to investigate local and regional scale population structure and gene flow along a salinity gradient (∼3 ppt to ∼18 ppt) in two North Carolina estuaries. G. bosc is an abundant and geographically widespread species that requires complex but patchy microhabitat (e.g. oyster reefs, rubble, woody debris) for reproduction and refuge. We sequenced 155 fish from 10 sites, using a common barcoding gene (COI). We also included recent sequence data from GenBank to determine how North Carolina populations fit into the larger biogeographic understanding of this species. In North Carolina, we found a significant amount of gene flow within and between estuaries. Our analysis also showed high predicted genetic diversity based upon a large number of rare haplotypes found within many of our sampled populations. Moreover, we detected a number of new haplotypes in North Carolina that had not yet been observed in prior work. Sampling along a salinity gradient did not reveal any significant positive or negative correlations between salinity and genetic diversity, nor the proportion of singleton haplotypes, with the exception of a positive correlation between salinity standard deviation and genetic diversity. We also found evidence that an introduced European population of naked gobies may have originated from an Atlantic source population. Altogether, this system offers a compelling way to evaluate whether factors other than dispersal per se mediate recruitment in an estuarine-dependent species of fish with a larval dispersal phase. It also demonstrates the importance of exploring both smaller and larger scale population structure in marine organisms to better understand local and regional patterns of population connectivity and gene flow.