Local and regional scale habitat heterogeneity contribute to genetic adaptation in a commercially important marine mollusc (Haliotis rubra) from southeastern Australia

Characterising adaptive genetic divergence among conspecific populations is often achieved by studying genetic variation across defined environmental gradients. In marine systems this is challenging due to a paucity of information on habitat heterogeneity at local and regional scales and a dependency on sampling regimes that are typically limited to broad longitudinal and latitudinal environmental gradients. As a result, the spatial scales at which selection processes operate and the environmental factors that contribute to genetic adaptation in marine systems are likely to be unclear. In this study we explore patterns of adaptive genetic structuring in a commercially‐ harvested abalone species (Haliotis rubra) from southeastern Australia, using a panel of genome‐wide SNP markers (5,239 SNPs), and a sampling regime informed by marine LiDAR bathymetric imagery and 20‐year hindcasted oceanographic models. Despite a lack of overall genetic structure across the sampling distribution, significant genotype associations with heterogeneous habitat features were observed at local and regional spatial scales, including associations with wave energy, ocean current, sea surface temperature, and geology. These findings provide insights into the potential resilience of the species to changing marine climates and the role of migration and selection on recruitment processes, with implications for conservation and fisheries management. This study points to the spatial scales at which selection processes operate in marine systems and highlights the benefits of geospatially‐informed sampling regimes for overcoming limitations associated with marine population genomic research.     

Divergence genetics analysis reveals historical population genetic processes leading to contrasting phylogeographic patterns in co-distributed species

Coalescent samplers are computational time machines for inferring the historical demographic genetic processes that have given rise to observable patterns of spatial genetic variation among contemporary populations. We have used traditional characterizations of population structure and coalescent‐based inferences about demographic processes to reconstruct the population histories of two co‐distributed marine species, the frilled dog whelk, Nucella lamellosa, and the bat star, Patiria miniata. Analyses of population structure were consistent with previous work in both species except that additional samples of N. lamellosa showed a larger regional genetic break on Vancouver Island (VI) rather than between the southern Alexander Archipelago as in P. miniata. Our understanding of the causes, rather than just the patterns, of spatial genetic variation was dramatically improved by coalescent analyses that emphasized variation in population divergence times. Overall, gene flow was greater in bat stars (planktonic development) than snails (benthic development) but spatially homogeneous within species. In both species, these large phylogeographic breaks corresponded to relatively ancient divergence times between populations rather than regionally restricted gene flow. Although only N. lamellosa shows a large break on VI, population separation times on VI are congruent between species, suggesting a similar response to late Pleistocene ice sheet expansion. The absence of a phylogeographic break in P. miniata on VI can be attributed to greater gene flow and larger effective population size in this species. Such insights put the relative significance of gene flow into a more comprehensive historical biogeographic context and have important implications for conservation and landscape genetic studies that emphasize the role of contemporary gene flow and connectivity in shaping patterns of population differentiation.     

Polymorphism in developmental mode and its effect on population genetic structure of a spionid polychaete, Pygospio elegans

Population genetic structure of sedentary marine species is expected to be shaped mainly by the dispersal ability of their larvae. Long-lived planktonic larvae can connect populations through migration and gene flow, whereas species with nondispersive benthic or direct-developing larvae are expected to have genetically differentiated populations. Poecilogonous species producing different larval types are ideal when studying the effect of developmental mode on population genetic structure and connectivity. In the spionid polychaete Pygospio elegans, different larval types have been observed between, and sometimes also within, populations. We used microsatellite markers to study population structure of European P. elegans from the Baltic Sea (BS) and North Sea (NS). We found that populations with planktonic larvae had higher genetic diversity than did populations with benthic larvae. However, this pattern may not be related to developmental mode, since in P. elegans, developmental mode may be associated with geography. Benthic larvae were more commonly seen in the brackish BS and planktonic larvae were predominant in the NS, although both larval types also are found from both areas. Significant isolation-by-distance (IBD) was found overall and within regions. Most of the pair-wise F(ST) comparisons among populations were significant, although some geographically close populations with planktonic larvae were found to be genetically similar. However, these results, together with the pattern of IBD, autocorrelation within populations, as well as high estimated local recruitment, suggest that dispersal is limited in populations with planktonic larvae as well as in those with benthic larvae. The decrease in salinity between the NS and BS causes a barrier to gene flow in many marine species. In P. elegans, low, but significant, differentiation was detected between the NS and BS (3.34% in AMOVA), but no clear transition zone was observed, indicating that larvae are not hampered by the change in salinity.     

Highly contrasted population genetic structures in a host–parasite pair in the Caribbean Sea

Evolution and population genetic structure of marine species across the Caribbean Sea are shaped by two complex factors: the geological history and the present pattern of marine currents. Characterizing and comparing the genetic structures of codistributed species, such as host–parasite associations, allow discriminating the relative importance of environmental factors and life history traits that influenced gene flow and demographic events. Using microsatellite and Cytochrome Oxidase I markers, we investigated if a host–parasite pair (the heart urchin Meoma ventricosa and its parasitic pea crab Dissodactylus primitivus) exhibits comparable population genetic structures in the Caribbean Sea and how the observed patterns match connectivity regions from predictive models and other taxa. Highly contrasting patterns were found: the host showed genetic homogeneity across the whole studied area, whereas the parasite displayed significant differentiation at regional and local scales. The genetic diversity of the parasitic crabs (both in microsatellites and COI) was distributed in two main groups, Panama–Jamaica–St Croix on the one hand, and the South‐Eastern Caribbean on the other. At a smaller geographical scale, Panamanian and Jamaican parasite populations were genetically more similar, while more genetic differentiation was found within the Lesser Antilles. Both species showed a signature of population expansion during the Quaternary. Some results match predictive models or data from previous studies (e.g., the Western‐Eastern dichotomy in the parasite) while others do not (e.g., genetic differentiation within the Lesser Antilles). The sharp dissimilarity of genetic structure of these codistributed species outlines the importance of population expansion events and/or contrasted patterns of gene flow. This might be linked to differences in several life history traits such as fecundity (higher for the host), swimming capacity of larval stages (higher for the parasite), and habitat availability (higher for the host).     

Gene flow in the anemone Anthopleura elegantissima limits signatures of local adaptation across an extensive geographic range

Species inhabiting marine environments face a wide range of environmental conditions that vary spatially across several orders of magnitude. The selective pressures that these conditions impose on marine organisms, in combination with potentially high rates of gene flow between distant populations, make it difficult to predict the extent to which these populations can locally adapt. Here, I identify how selection and gene flow influence the population genetic structure of the anemone Anthopleura elegantissima along the Pacific coast of North America. Isolation by distance is the dominant pattern across the range of this species, with a genetic break near Pt. Conception, CA. Furthermore, demographic modelling suggests that this species was historically confined to southerly latitudes before expanding northward. Outlier analyses identify 24 loci under selection (out of ~1,100), but the same analysis on simulated genetic data generated using the most likely demographic model erroneously identified the same number of loci under selection, if not more. Taken together, these results suggest that demographic processes are the dominant force shaping population genetic patterns in A. elegantissima along the Pacific coast of North America. I discuss these patterns in terms of the evolutionary history of A. elegantissima, the potential for local adaptation, and their consequences with respect to interactions with the endosymbiont Breviolum muscatinei across their geographic range.     

From global to local genetic structuring in the red gorgonian Paramuricea clavata: the interplay between oceanographic conditions and limited larval dispersal

Defining the scale of connectivity among marine populations and identifying the barriers to gene flow are tasks of fundamental importance for understanding the genetic structure of populations and for the design of marine reserves. Here, we investigated the population genetic structure at three spatial scales of the red gorgonian Paramuricea clavata (Cnidaria, Octocorallia), a key species dwelling in the coralligenous assemblages of the Mediterranean Sea. Colonies of P. clavata were collected from 39 locations across the Mediterranean Sea from Morocco to Turkey and analysed using microsatellite loci. Within three regions (Medes, Marseille and North Corsica), sampling was obtained from multiple locations and at different depths. Three different approaches (measures of genetic differentiation, Bayesian clustering and spatially explicit maximum‐difference algorithm) were used to determine the pattern of genetic structure. We identified genetic breaks in the spatial distribution of genetic diversity, which were concordant with oceanographic conditions in the Mediterranean Sea. We revealed a high level of genetic differentiation among populations and a pattern of isolation by distance across the studied area and within the three regions, underlining short effective larval dispersal in this species. We observed genetic differentiation among populations in the same locality dwelling at different depths, which may be explained by local oceanographic conditions and which may allow a process of local adaptation of the populations to their environment. We discuss the implications of our results for the conservation of the species, which is exposed to various threats.     

Ocean circulation model predicts high genetic structure observed in a long‐lived pelagic developer

Understanding the movement of genes and individuals across marine seascapes is a long‐standing challenge in marine ecology and can inform our understanding of local adaptation, the persistence and movement of populations, and the spatial scale of effective management. Patterns of gene flow in the ocean are often inferred based on population genetic analyses coupled with knowledge of species' dispersive life histories. However, genetic structure is the result of time‐integrated processes and may not capture present‐day connectivity between populations. Here, we use a high‐resolution oceanographic circulation model to predict larval dispersal along the complex coastline of western Canada that includes the transition between two well‐studied zoogeographic provinces. We simulate dispersal in a benthic sea star with a 6–10 week pelagic larval phase and test predictions of this model against previously observed genetic structure including a strong phylogeographic break within the zoogeographical transition zone. We also test predictions with new genetic sampling in a site within the phylogeographic break. We find that the coupled genetic and circulation model predicts the high degree of genetic structure observed in this species, despite its long pelagic duration. High genetic structure on this complex coastline can thus be explained through ocean circulation patterns, which tend to retain passive larvae within 20–50 km of their parents, suggesting a necessity for close‐knit design of Marine Protected Area networks.     

Genetic effects of habitat fragmentation and population isolation on Etheostoma raneyi (Percidae)

The use of genetic methods to quantify the effects of anthropogenic habitat fragmentation on population structure has become increasingly common. However, in today’s highly fragmented habitats, researchers have sometimes concluded that populations are currently genetically isolated due to habitat fragmentation without testing the possibility that populations were genetically isolated before European settlement. Etheostoma raneyi is a benthic headwater fish restricted to river drainages in northern Mississippi, USA, that has a suite of adaptive traits that correlate with poor dispersal ability. Aquatic habitat within this area has been extensively modified, primarily by flood-control projects, and populations in headwater streams have possibly become genetically isolated from one another. We used microsatellite markers to quantify genetic structure as well as contemporary and historical gene flow across the range of the species. Results indicated that genetically distinct populations exist in each headwater stream analyzed, current gene flow rates are lower than historical rates, most genetic variation is partitioned among populations, and populations in the Yocona River drainage show lower levels of genetic diversity than populations in the Tallahatchie River drainage and other Etheostoma species. All populations have negative F_IS scores, of which roughly half are significant relative to Hardy–Weinberg expectations, perhaps due to small population sizes. We conclude that anthropogenic habitat alteration and fragmentation has had a profoundly negative impact on the species by isolating E. raneyi within headwater stream reaches. Further research is needed to inform conservation strategies, but populations in the Yocona River drainage are in dire need of management action. Carefully planned human-mediated dispersal and habitat restoration should be explored as management options across the range of the species.     

Habitat continuity and stepping‐stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping‐stone oceanographic transport and habitat continuity, using as model an ecosystem‐structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping‐stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life‐history traits. Our results highlight the importance of spatially explicit modelling of stepping‐stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.     

Niche shifts and range expansions along cordilleras drove diversification in a high‐elevation endemic plant genus in the tropical Andes

The tropical Andes represent one of the world's biodiversity hot spots, but the evolutionary drivers generating their striking species diversity still remain poorly understood. In the treeless high‐elevation Andean environments, Pleistocene glacial oscillations and niche differentiation are frequently hypothesized diversification mechanisms; however, sufficiently densely sampled population genetic data supporting this are still lacking. Here, we reconstruct the evolutionary history of Loricaria (Asteraceae), a plant genus endemic to the Andean treeless alpine zone, based on comprehensive population‐level sampling of 289 individuals from 67 populations across the entire distribution ranges of its northern Andean species. Partly incongruent AFLP and plastid DNA markers reveal that the distinct genetic structure was shaped by a complex interplay of biogeography (spread along and across the cordilleras), history (Pleistocene glacial oscillations) and local ecological conditions. While plastid variation documents an early split or colonization of the northern Andes by at least two lineages, one of which further diversified, a major split in the AFLP data correlate with altitudinal ecological differentiation. This suggests that niche shifts may be important drivers of Andean diversification not only in forest–alpine transitions, but also within the treeless alpine zone itself. The patterns of genetic differentiation at the intraspecific level reject the hypothesized separation in spatially isolated cordilleras and instead suggest extensive gene flow among populations from distinct mountain chains. Our study highlights that leveraging highly variable markers against extensive population‐level sampling is a promising approach to address mechanisms of rapid species diversifications.     

Range‐wide population genetic structure of the Caribbean sea fan coral,Gorgonia ventalina

The population structure of benthic marine organisms is of central relevance to the conservation and management of these often threatened species, as well as to the accurate understanding of their ecological and evolutionary dynamics. A growing body of evidence suggests that marine populations can be structured over short distances despite theoretically high dispersal potential. Yet the proposed mechanisms governing this structure vary, and existing empirical population genetic evidence is of insufficient taxonomic and geographic scope to allow for strong general inferences. Here, we describe the range‐wide population genetic structure of an ecologically important Caribbean octocoral, Gorgonia ventalina. Genetic differentiation was positively correlated with geographic distance and negatively correlated with oceanographically modelled dispersal probability throughout the range. Although we observed admixture across hundreds of kilometres, estimated dispersal was low, and populations were differentiated across distances <2 km. These results suggest that populations of G. ventalina may be evolutionarily coupled via gene flow but are largely demographically independent. Observed patterns of differentiation corroborate biogeographic breaks found in other taxa (e.g. an east/west divide near Puerto Rico), and also identify population divides not discussed in previous studies (e.g. the Yucatan Channel). High genotypic diversity and absence of clonemates indicate that sex is the primary reproductive mode for G. ventalina. A comparative analysis of the population structure of G. ventalina and its dinoflagellate symbiont, Symbiodinium, indicates that the dispersal of these symbiotic partners is not coupled, and symbiont transmission occurs horizontally.     

Self‐recruitment and sweepstakes reproduction amid extensive gene flow in a coral‐reef fish

Identifying patterns of larval dispersal within marine metapopulations is vital for effective fisheries management, appropriate marine reserve design, and conservation efforts. We employed genetic markers (microsatellites) to determine dispersal patterns in bicolour damselfish (Pomacentridae: Stegastes partitus). Tissue samples of 751 fish were collected in 2004 and 2005 from 11 sites encompassing the Exuma Sound, Bahamas. Bayesian parentage analysis identified two parent–offspring pairs, which is remarkable given the large population sizes and 28 day pelagic larval duration of bicolour damselfish. The two parent–offspring pairs directly documented self‐recruitment at the two northern‐most sites, one of which is a long‐established marine reserve. Principal coordinates analyses of pair‐wise relatedness values further indicated that self‐recruitment was common in all sampled populations. Nevertheless, measures of genetic differentiation (F_ST) and results from assignment methods suggested high levels of gene flow among populations. Comparisons of heterozygosity and relatedness among samples of adults and recruits indicated spatially and temporally independent sweepstakes events, whereby only a subset of adults successfully contribute to subsequent generations. These results indicate that self‐recruitment and sweepstakes reproduction are the predominant, ecologically‐relevant processes that shape patterns of larval dispersal in this system.     

HETEROGENEOUS GENOMIC DIFFERENTIATION IN MARINE THREESPINE STICKLEBACKS: ADAPTATION ALONG AN ENVIRONMENTAL GRADIENT

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on‐going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F_ST ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F_ST ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers—45% of which were identified as outliers—correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.     

A comparative study of asymmetric migration events across a marine biogeographic boundary

In many nonclonal, benthic marine species, geographic distribution is mediated by the dispersal of their larvae. The dispersal and recruitment of marine larvae may be limited by temperature gradients that can affect mortality or by ocean currents that can directly affect the movements of pelagic larvae. We focus on Point Conception, a well-known biogeographic boundary between the Californian and Oregonian biogeographic provinces, to investigate whether ocean currents affect patterns of gene flow in intertidal marine invertebrates. The predominance of pelagically dispersing species with northern range limits at Point Conception suggests that ocean currents can affect species distributions by erecting barriers to the dispersal of planktonic larvae. In this paper, we investigate whether the predominantly southward currents have left a recognizable genetic signature in species with pelagically dispersing larvae whose ranges span Point Conception. We use patterns of genetic diversity and a new method for inferring cladistic migration events to test the hypothesis that southward currents increase southward gene flow for species with pelagically dispersing larvae. We collected mitochondrial DNA (mtDNA) sequence data for the barnacles Balanus glandula and Chthamalus fissus and also reanalyzed a previously published mtDNA dataset (Strongylocentrotus purpuratus, Edmands et al. 1996). For all three species, our cladistic approach identified an excess of southward migration events across Point Conception. In data from a fourth species with nondispersing larvae (Nucella emarginata, Marko 1998), our method suggests that ocean currents have not played a role in generating genetic structure.     

Historical influences dominate the population genetic structure of a sedentary marine fish,Atlantic wolffish (Anarhichas lupus), across the North Atlantic Ocean

Genetic variation was assessed in Atlantic wolffish, Anarhichas lupus, across the North Atlantic Ocean using microsatellite and amplified fragment length polymorphism (AFLP) markers. Despite unusual life history attributes such as large benthic eggs, large larvae, a limited pelagic stage and relatively sedentary adults, which suggest potential for strong population structure, range‐wide F_ST values were comparable to other marine fishes (≤0.035). Nevertheless, both significant genetic differentiation among regions and isolation by distance were observed, suggesting limited dispersal in this species. AFLP loci, evaluated on a subset of samples, revealed slightly higher F_ST values, but similar patterns of differentiation and isolation‐by‐distance estimates, compared to microsatellites. The genetic structure of Atlantic wolffish has likely been shaped by its post‐glacial history of recolonization, North Atlantic current patterns and continuity of habitat on continental shelves.     

Using hierarchical sampling to understand scales of spatial variation in early coral recruitment

Ecological patterns are created by processes acting over multiple spatial and temporal scales. By combining spatially explicit sampling with variance components models, the relative importance of spatial scale to overall variability can be determined. We used a spatially structured experimental design in the Mombasa Marine National Park in Kenya to quantify variation in coral recruitment across four spatial scales (~1–1,000 m) and to generate hypotheses about processes affecting recruitment and potential sources of post-settlement mortality during early life history. For the dominant recruiting corals (Pocillopora spp.), variation in recruitment on surfaces protected from fish grazing was greatest at the largest spatial scale examined (1,000 m). We hypothesize that recruitment on protected surfaces varies mainly with larval delivery due to different lagoonal circulation and water flow between sites. Conversely, variation on surfaces exposed to fishes was greatest at the smallest spatial scale (1 m). We hypothesize that recruitment on exposed surfaces mainly reflects local differences in the scale and intensity of fish grazing, which may obscure larval delivery patterns. Spatial variation in recruitment can affect many ecological processes and factors, including growth, survival to maturity, the distribution of habitat, and variation in species interaction strengths. This study demonstrates how spatially explicit sampling, followed by variance components modeling to partition variance across scales, can help to identify potential drivers of patterns at each relevant scale.     

Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea

The relatively recent fields of terrestrial landscape and marine seascape genetics seek to identify the influence of biophysical habitat features on the spatial genetic structure of populations or individuals. Over the last few years, there has been accumulating evidence for the effect of environmental heterogeneity on patterns of gene flow and connectivity in marine systems. Here, we investigate the population genetic patterns of an anemonefish, Amphiprion bicinctus, along the Saudi Arabian coast of the Red Sea. We collected nearly one thousand samples from 19 locations, spanning approximately 1500 km, and genotyped them at 38 microsatellite loci. Patterns of gene flow appeared to follow a stepping‐stone model along the northern and central Red Sea, which was disrupted by a distinct genetic break at a latitude of approximately 19°N. The Red Sea is characterized by pronounced environmental gradients along its axis, roughly separating the northern and central from the southern basin. Using mean chlorophyll‐a concentrations as a proxy for this gradient, we ran tests of isolation by distance (IBD, R² = 0.52) and isolation by environment (IBE, R² = 0.64), as well as combined models using partial Mantel tests and multiple matrix regression with randomization (MMRR). We found that genetic structure across our sampling sites may be best explained by a combined model of IBD and IBE (Mantel: R² = 0.71, MMRR: R² = 0.86). Our results highlight the potential key role of environmental patchiness in shaping patterns of gene flow in species with pelagic larval dispersal. We support growing calls for the integration of biophysical habitat characteristics into future studies of population genetic structure.     

A modified stepping-stone model of population structure in red drum, Sciaenops ocellatus (Sciaenidae), from the northern Gulf of Mexico

Genetic studies of population or stock structure in exploited marine fishes typically are designed to determine whether geographic boundaries useful for conservation and management planning are identifiable. Implicit in many such studies is the notion that subpopulations or stocks, if they exist, have fixed territories with little or no gene exchange between them. Herein, we review our long-term genetic studies of red drum (Sciaenops ocellatus), an estuarine-dependent sciaenid fish in the Gulf of Mexico and western Atlantic Ocean. Significant differences in frequencies of mitochondrial DNA haplotypes and of alleles at nuclear-encoded microsatellites occur among red drum sampled across the northern Gulf of Mexico. The spatial distribution of the genetic variation, however, follows a pattern of isolation-by-distance consistent with the hypothesis that gene flow occurs among subpopulations and is an inverse (and continuous) function of geographic distance. However, successful reproduction and recruitment of red drum depend on estuarine habitats that have geographically discrete boundaries. We hypothesize that population structure in red drum follows a modified one-dimensional, linear stepping-stone model where gene exchange occurs primarily (but not exclusively) between adjacent bays and estuaries distributed linearly along the coastline. Gene flow does occur among estuaries that are not adjacent but probabilities of gene exchange decrease as a function of geographic distance. Implications of our hypothesis are discussed in terms of inferences drawn from patterns of isolation-by-distance and relative to conservation and management of estuarine-dependent species like red drum. Based on estimates of the ratio of genetic effective population size and census size in red drum, observed patterns of gene flow in red drum may play a significant role in recruitment.     

New view of population genetics of zooplankton: RAD‐seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium‐to‐small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next‐generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site‐associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High‐density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.     

Using dense locality sampling resolves the subtle genetic population structure of the dispersive fish species Plecoglossus altivelis

In dispersive species with continuous distributions, genetic differentiation between local populations is often absent or subtle and thus difficult to detect. To incorporate such subtle differentiation into management plans, it may be essential to analyse many samples from many localities using adequate numbers of high‐resolution genetic markers. Here, we evaluated the usefulness of dense locality sampling in resolving genetic population structure in the ayu (Plecoglossus altivelis), a dispersive fish important in Japanese inland fisheries. Genetic variability in, and differentiation between, ayu populations around the Japan–Ryukyu Archipelago were investigated in 4746 individuals collected from 120 localities by genotyping 12 microsatellite markers. These individuals represented the two subspecies of ayu, namely the Ryukyuan subspecies (Plecoglossus altivelis ryukyuensis) and both amphidromous and landlocked forms of the nominotypical subspecies (P. a. altivelis) along the archipelago. We successfully detected an absence of genetic differentiation within the landlocked form and subtle but significant differentiation and clear geographic patterns of genetic variation among populations of the amphidromous form, which had been considered genetically homogeneous. This suggests that dense locality sampling effectively resolves subtle differences in genetic population structure, reducing stochastic deviation in the detection of genetic differentiation and geographic patterns in local populations of this dispersive species. Resampling analyses based on empirical data sets clearly demonstrate the effectiveness of increasing the number of locality samples for stable and reliable estimations of genetic fixation indices. The genetic population structure observed within the amphidromous form provides useful information for identifying management or conservation units in ayu.