Isolation by oceanographic distance explains genetic structure for Macrocystis pyrifera in the Santa Barbara Channel

Ocean currents are expected to be the predominant environmental factor influencing the dispersal of planktonic larvae or spores; yet, their characterization as predictors of marine connectivity has been hindered by a lack of understanding of how best to use oceanographic data. We used a high-resolution oceanographic model output and Lagrangian particle simulations to derive oceanographic distances (hereafter called transport times) between sites studied for Macrocystis pyrifera genetic differentiation. We build upon the classical isolation-by-distance regression model by asking how much additional variability in genetic differentiation is explained when adding transport time as predictor. We explored the extent to which gene flow is dependent upon seasonal changes in ocean circulation. Because oceanographic transport between two sites is inherently asymmetric, we also compare the explanatory power of models using the minimum or the mean transport times. Finally, we compare the direction of connectivity as estimated by the oceanographic model and genetic assignment tests. We show that the minimum transport time had higher explanatory power than the mean transport time, revealing the importance of considering asymmetry in ocean currents when modelling gene flow. Genetic assignment tests were much less effective in determining asymmetry in gene flow. Summer-derived transport times, in particular for the month of June, which had the strongest current speed, greatest asymmetry and highest spore production, resulted in the best-fit model explaining twice the variability in genetic differentiation relative to models that use geographic distance or habitat continuity. The best overall model also included habitat continuity and explained 65% of the variation in genetic differentiation among sites.     

Coalescent and biophysical models of stepping‐stone gene flow in neritid snails

Marine species in the Indo‐Pacific have ranges that can span thousands of kilometres, yet studies increasingly suggest that mean larval dispersal distances are less than historically assumed. Gene flow across these ranges must therefore rely to some extent on larval dispersal among intermediate ‘stepping‐stone’ populations in combination with long‐distance dispersal far beyond the mean of the dispersal kernel. We evaluate the strength of stepping‐stone dynamics by employing a spatially explicit biophysical model of larval dispersal in the tropical Pacific to construct hypotheses for dispersal pathways. We evaluate these hypotheses with coalescent models of gene flow among high‐island archipelagos in four neritid gastropod species. Two of the species live in the marine intertidal, while the other two are amphidromous, living in fresh water but retaining pelagic dispersal. Dispersal pathways predicted by the biophysical model were strongly favoured in 16 of 18 tests against alternate hypotheses. In regions where connectivity among high‐island archipelagos was predicted as direct, there was no difference in gene flow between marine and amphidromous species. In regions where connectivity was predicted through stepping‐stone atolls only accessible to marine species, gene flow estimates between high‐island archipelagos were significantly higher in marine species. Moreover, one of the marine species showed a significant pattern of isolation by distance consistent with stepping‐stone dynamics. While our results support stepping‐stone dynamics in Indo‐Pacific species, we also see evidence for nonequilibrium processes such as range expansions or rare long‐distance dispersal events. This study couples population genetic and biophysical models to help to shed light on larval dispersal pathways.     

Seascape genetics of the stalked kelp Pterygophora californica and comparative population genetics in the Santa Barbara Channel

We conducted a population genetic analysis of the stalked kelp, Pterygophora californica, in the Santa Barbara Channel, California, USA. The results were compared with previous work on the genetic differentiation of giant kelp, Macrocystis pyrifera, in the same region. These two sympatric kelps not only share many life history and dispersal characteristics but also differ in that dislodged P. californica does not produce floating rafts with buoyant fertile sporophytes, commonly observed for M. pyrifera. We used a comparative population genetic approach with these two species to test the hypothesis that the ability to produce floating rafts increases the genetic connectivity among kelp patches in the Santa Barbara Channel. We quantified the association of habitat continuity and oceanographic distance with the genetic differentiation observed in stalked kelp, like previously conducted for giant kelp. We compared both overall (across all patches) and pairwise (between patches) genetic differentiation. We found that oceanographic transit time, habitat continuity, and geographic distance were all associated with genetic connectivity in P. californica, supporting similar previous findings for M. pyrifera. Controlling for differences in heterozygosity between kelp species using Jost's D_EST, we showed that global differentiation and pairwise differentiation were similar among patches between the two kelp species, indicating that they have similar dispersal capabilities despite their differences in rafting ability. These results suggest that rafting sporophytes do not play a significant role in effective dispersal of M. pyrifera at ecologically relevant spatial and temporal scales.     

Matching genetics with oceanography: directional gene flow in a Mediterranean fish species

Genetic connectivity and geographic fragmentation are two opposing mechanisms determining the population structure of species. While the first homogenizes the genetic background across populations the second one allows their differentiation. Therefore, knowledge of processes affecting dispersal of marine organisms is crucial to understand their genetic distribution patterns and for the effective management of their populations. In this study, we use genetic analyses of eleven microsatellites in combination with oceanographic satellite and dispersal simulation data to determine distribution patterns for Serranus cabrilla, a ubiquitous demersal broadcast spawner, in the Mediterranean Sea. Pairwise population F_ST values ranged between ?0.003 and 0.135. Two genetically distinct clusters were identified, with a clear division located between the oceanographic discontinuities at the Ibiza Channel (IC) and the Almeria‐Oran Front (AOF), revealing an admixed population in between. The Balearic Front (BF) also appeared to dictate population structure. Directional gene flow on the Spanish coast was observed as S. cabrilla dispersed from west to east over the AOF, from north to south on the IC and from south of the IC towards the Balearic Islands. Correlations between genetic and oceanographic data were highly significant. Seasonal changes in current patterns and the relationship between ocean circulation patterns and spawning season may also play an important role in population structure around oceanographic fronts.     

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.     

DECOUPLING OF SHORT‐ AND LONG‐DISTANCE DISPERSAL PATHWAYS IN THE ENDEMIC NEW ZEALAND SEAWEED CARPOPHYLLUM MASCHALOCARPUM (PHAEOPHYCEAE,FUCALES)1

The processes that produce and maintain genetic structure in organisms operate at different timescales and on different life‐history stages. In marine macroalgae, gene flow occurs through gamete/zygote dispersal and rafting by adult thalli. Population genetic patterns arise from this contemporary gene flow interacting with historical processes. We analyzed spatial patterns of mitochondrial DNA variation to investigate contemporary and historical dispersal patterns in the New Zealand endemic fucalean brown alga Carpophyllum maschalocarpum (Turner) Grev. Populations bounded by habitat discontinuities were often strongly differentiated from adjoining populations over scales of tens of kilometers and intrapopulation diversity was generally low, except for one region of northeast New Zealand (the Bay of Plenty). There was evidence of strong connectivity between the northern and eastern regions of New Zealand’s North Island and between the North and South Islands of New Zealand and the Chatham Islands (separated by 650 km of open ocean). Moderate haplotypic diversity was found in Chatham Islands populations, while other southern populations showed low diversity consistent with Last Glacial Maximum (LGM) retreat and subsequent recolonization. We suggest that ocean current patterns and prevailing westerly winds facilitate long‐distance dispersal by floating adult thalli, decoupling genetic differentiation of Chatham Island populations from dispersal potential at the gamete/zygote stage. This study highlights the importance of encompassing the entire range of a species when inferring dispersal patterns from genetic differentiation, as realized dispersal distances can be contingent on local or regional oceanographic and historical processes.     

Nonsignificant isolation by distance implies limited dispersal

Understanding the scale of dispersal is an important consideration in the conservation and management of many species. However, in species in which the high‐dispersal stage is characterized by tiny gametes or offspring, it may be difficult to estimate dispersal directly. This is the case for many marine species, whose pelagic larvae are dispersed by ocean currents by several days or weeks before beginning a benthic, more sedentary, adult stage. As consequence of the high‐dispersal larval stage, many marine species have low genetic structure on large spatial scales (Waples 1998 ; Hellberg 2007 ). Despite the high capacity for dispersal, some tagging studies have found that a surprising number of larvae recruit into the population they were released from (self‐recruitment). However, estimates of self‐recruitment are not informative about mean dispersal between subpopulations. To what extent are limited dispersal estimates from tagging studies compatible with high potential for dispersal and low genetic structure? In this issue, a study on five species of coral reef fish used isolation by distance (IBD) between individuals to estimate mean dispersal distances (Puebla et al. 2012 ). They found that mean dispersal was unexpectedly small (<50 km), given relatively low IBD slopes and long pelagic durations. This study demonstrates how low genetic structure is compatible with limited dispersal in marine species. A comprehensive understanding of dispersal in marine species will involve integrating methods that estimate dispersal over different spatial and temporal scales. Genomic data may increase power to resolve these issues but must be applied carefully to this question.     

Population differentiation in the open sea: insights from the pelagic copepod Pleuromamma xiphias

Although a number of recent studies of marine holoplankton have reported significant genetic structure among populations, little is currently known about the biological and oceanographic processes that influence population connectivity in oceanic plankton. In order to examine how depth preferences influence dispersal in oceanic plankton, I characterized the genetic structure of a copepod with diel vertical migration (DVM) (Pleuromamma xiphias), throughout its global distribution, and compared these results to those expected given the interaction of this species' habitat depth with ocean circulation and bathymetry. Mitochondrial COI sequences from 651 individuals from 28 sites in the Indian, Pacific, and Atlantic Oceans revealed highly significant genetic differentiation both within and among ocean basins. Limited dispersal among distinct pelagic provinces seems to have played a major role in population differentiation in this species, with strong genetic breaks observed across known oceanographic fronts or current systems in all three ocean basins. The Indo-West Pacific (IWP) holds a highly distinct genetic population of this species that was sampled in both the western Pacific and eastern Indian Oceans. This suggests that the IWP does not act as a strong barrier to gene flow between basins, as expected, despite the relatively shallow water depth (<200 m) and vertically extensive (>400 m) diel migration of this species. A pattern of isolation by distance was observed in the Indian Ocean with genetic differentiation among samples down to spatial scales of ～800 km, indicating that realized dispersal in P. xiphias occurs over much smaller spatial scales than in previously reported oceanic holoplankton. Given its highly regionalized population genetic structure, P. xiphias may have some capacity to adapt to local oceanographic conditions, and it should not be assumed that populations of this species in distinct pelagic biomes will respond in the same way to shared physical or climatic forcing.     

Small fish,large river: Surprisingly minimal genetic structure in a dispersal‐limited,habitat specialist fish

Genetic connectivity is expected to be lower in species with limited dispersal ability and a high degree of habitat specialization (intrinsic factors). Also, gene flow is predicted to be limited by habitat conditions such as physical barriers and geographic distance (extrinsic factors). We investigated the effects of distance, intervening pools, and rapids on gene flow in a species, the Tuxedo Darter (Etheostoma lemniscatum), a habitat specialist that is presumed to be dispersal‐limited. We predicted that the interplay between these intrinsic and extrinsic factors would limit dispersal and lead to genetic structure even at the small spatial scale of the species range (a 38.6 km river reach). The simple linear distribution of E. lemniscatum allowed for an ideal test of how these factors acted on gene flow and allowed us to test expectations (e.g., isolation‐by‐distance) of linearly distributed species. Using 20 microsatellites from 163 individuals collected from 18 habitat patches, we observed low levels of genetic structure that were related to geographic distance and rapids, though these factors were not barriers to gene flow. Pools separating habitat patches did not contribute to any observed genetic structure. Overall, E. lemniscatum maintains gene flow across its range and is comprised of a single population. Due to the linear distribution of the species, a stepping‐stone model of dispersal best explains the maintenance of gene flow across its small range. In general, our observation of higher‐than‐expected connectivity likely stems from an adaptation to disperse due to temporally unstable and patchy habitat.     

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite‐based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large‐scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation‐by‐environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.     

Asymmetric oceanographic processes mediate connectivity and population genetic structure,as revealed by RADseq,in a highly dispersive marine invertebrate (Parastichopus californicus)

Marine populations are typically characterized by weak genetic differentiation due to the potential for long‐distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (F_ST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R² = .49) than geographic distance (IBD) (R² = .18), and directional processes played an important role in shaping fine‐scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.     

The oceanic concordance of phylogeography and biogeography: a case study in Notochthamalus

Dispersal and adaptation are the two primary mechanisms that set the range distributions for a population or species. As such, understanding how these mechanisms interact in marine organisms in particular – with capacity for long‐range dispersal and a poor understanding of what selective environments species are responding to – can provide useful insights for the exploration of biogeographic patterns. Previously, the barnacle Notochthamalus scabrosus has revealed two evolutionarily distinct lineages with a joint distribution that suggests an association with one of the two major biogeographic boundaries (~30°S) along the coast of Chile. However, spatial and genomic sampling of this system has been limited until now. We hypothesized that given the strong oceanographic and environmental shifts associated with the other major biogeographic boundary (~42°S) for Chilean coastal invertebrates, the southern mitochondrial lineage would dominate or go to fixation in locations further to the south. We also evaluated nuclear polymorphism data from 130 single nucleotide polymorphisms to evaluate the concordance of the signal from the nuclear genome with that of the mitochondrial sample. Through the application of standard population genetic approaches along with a Lagrangian ocean connectivity model, we describe the codistribution of these lineages through a simultaneous evaluation of coastal lineage frequencies, an approximation of larval behavior, and current‐driven dispersal. Our results show that this pattern could not persist without the two lineages having distinct environmental optima. We suggest that a more thorough integration of larval dynamics, explicit dispersal models, and near‐shore environmental analysis can explain much of the coastal biogeography of Chile.     

Seascape genetics of a globally distributed, highly mobile marine mammal: the short-beaked common dolphin (genus Delphinus)

Identifying which factors shape the distribution of intraspecific genetic diversity is central in evolutionary and conservation biology. In the marine realm, the absence of obvious barriers to dispersal can make this task more difficult. Nevertheless, recent studies have provided valuable insights into which factors may be shaping genetic structure in the world's oceans. These studies were, however, generally conducted on marine organisms with larval dispersal. Here, using a seascape genetics approach, we show that marine productivity and sea surface temperature are correlated with genetic structure in a highly mobile, widely distributed marine mammal species, the short-beaked common dolphin. Isolation by distance also appears to influence population divergence over larger geographical scales (i.e. across different ocean basins). We suggest that the relationship between environmental variables and population structure may be caused by prey behaviour, which is believed to determine common dolphins' movement patterns and preferred associations with certain oceanographic conditions. Our study highlights the role of oceanography in shaping genetic structure of a highly mobile and widely distributed top marine predator. Thus, seascape genetic studies can potentially track the biological effects of ongoing climate-change at oceanographic interfaces and also inform marine reserve design in relation to the distribution and genetic connectivity of charismatic and ecologically important megafauna.     

Population genomics meet Lagrangian simulations: Oceanographic patterns and long larval duration ensure connectivity among Paracentrotus lividus populations in the Adriatic and Ionian seas

Connectivity between populations influences both their dynamics and the genetic structuring of species. In this study, we explored connectivity patterns of a marine species with long‐distance dispersal, the edible common sea urchin Paracentrotus lividus, focusing mainly on the Adriatic–Ionian basins (Central Mediterranean). We applied a multidisciplinary approach integrating population genomics, based on 1,122 single nucleotide polymorphisms (SNPs) obtained from 2b‐RAD in 275 samples, with Lagrangian simulations performed with a biophysical model of larval dispersal. We detected genetic homogeneity among eight population samples collected in the focal Adriatic–Ionian area, whereas weak but significant differentiation was found with respect to two samples from the Western Mediterranean (France and Tunisia). This result was not affected by the few putative outlier loci identified in our dataset. Lagrangian simulations found a significant potential for larval exchange among the eight Adriatic–Ionian locations, supporting the hypothesis of connectivity of P. lividus populations in this area. A peculiar pattern emerged from the comparison of our results with those obtained from published P. lividus cytochrome b (cytb) sequences, the latter revealing genetic differentiation in the same geographic area despite a smaller sample size and a lower power to detect differences. The comparison with studies conducted using nuclear markers on other species with similar pelagic larval durations in the same Adriatic–Ionian locations indicates species‐specific differences in genetic connectivity patterns and warns against generalizing single‐species results to the entire community of rocky shore habitats.     

Asymmetric dispersal of southern bull‐kelp (Durvillaea antarctica) adults in coastal New Zealand: testing an oceanographic hypothesis

Coastal populations are often connected by unidirectional current systems, but the biological effects of such asymmetric oceanographic connectivity remain relatively unstudied. We used mtDNA analysis to determine the phylogeographic origins of beach‐cast bull‐kelp (Durvillaea antarctica) adults in the Canterbury Bight, a 180 km coastal region devoid of rocky‐reef habitat in southern New Zealand. A multi‐year, quantitative analysis supports the oceanographically derived hypothesis of asymmetric dispersal mediated by the north‐flowing Southland Current. Specifically, 92% of beach‐cast specimens examined had originated south of the Bight, many drifting north for hundreds of kilometres, and some traversing at least 500 km of ocean from subantarctic sources. In contrast, only 8% of specimens had dispersed south against the prevailing current, and these counter‐current dispersers likely travelled relatively small distances (tens of kilometres). These data show that oceanographic connectivity models can provide robust estimates of passive biological dispersal, even for highly buoyant taxa. The results also indicate that there are no oceanographic barriers to kelp dispersal across the Canterbury Bight, indicating that other ecological factors explain the phylogeographic disjunction across this kelp‐free zone. The large number of long‐distance dispersal events detected suggests drifting macroalgae have potential to facilitate ongoing connectivity between otherwise isolated benthic populations.     

Population structure and dispersal of the coral‐excavating sponge Cliona delitrix

Some excavating sponges of the genus Cliona compete with live reef corals, often killing and bioeroding entire colonies. Important aspects affecting distribution of these species, such as dispersal capability and population structure, remain largely unknown. Thus, the aim of this study was to determine levels of genetic connectivity and dispersal of Cliona delitrix across the Greater Caribbean (Caribbean Sea, Bahamas and Florida), to understand current patterns and possible future trends in their distribution and effects on coral reefs. Using ten species‐specific microsatellite markers, we found high levels of genetic differentiation between six genetically distinct populations: one in the Atlantic (Florida‐Bahamas), one specific to Florida and four in the South Caribbean Sea. In Florida, two independent breeding populations are likely separated by depth. Gene flow and ecological dispersal occur among other populations in the Florida reef tract, and between some Florida locations and the Bahamas. Similarly, gene flow occurs between populations in the South Caribbean Sea, but appears restricted between the Caribbean Sea and the Atlantic (Florida‐Bahamas). Dispersal of C. delitrix was farther than expected for a marine sponge and favoured in areas where currents are strong enough to transport sponge eggs or larvae over longer distances. Our results support the influence of ocean current patterns on genetic connectivity, and constitute a baseline to monitor future C. delitrix trends under climate change.     

Local adaptation and oceanographic connectivity patterns explain genetic differentiation of a marine diatom across the North Sea–Baltic Sea salinity gradient

Drivers of population genetic structure are still poorly understood in marine micro‐organisms. We exploited the North Sea–Baltic Sea transition for investigating the seascape genetics of a marine diatom, Skeletonema marinoi. Eight polymorphic microsatellite loci were analysed in 354 individuals from ten locations to analyse population structure of the species along a 1500‐km‐long salinity gradient ranging from 3 to 30 psu. To test for salinity adaptation, salinity reaction norms were determined for sets of strains originating from three different salinity regimes of the gradient. Modelled oceanographic connectivity was compared to directional relative migration by correlation analyses to examine oceanographic drivers. Population genetic analyses showed distinct genetic divergence of a low‐salinity Baltic Sea population and a high‐salinity North Sea population, coinciding with the most evident physical dispersal barrier in the area, the Danish Straits. Baltic Sea populations displayed reduced genetic diversity compared to North Sea populations. Growth optima of low salinity isolates were significantly lower than those of strains from higher native salinities, indicating local salinity adaptation. Although the North Sea–Baltic Sea transition was identified as a barrier to gene flow, migration between Baltic Sea and North Sea populations occurred. However, the presence of differentiated neutral markers on each side of the transition zone suggests that migrants are maladapted. It is concluded that local salinity adaptation, supported by oceanographic connectivity patterns creating an asymmetric migration pattern between the Baltic Sea and the North Sea, determines genetic differentiation patterns in the transition zone.     

Turtle groups or turtle soup: dispersal patterns of hawksbill turtles in the Caribbean

Despite intense interest in conservation of marine turtles, spatial ecology during the oceanic juvenile phase remains relatively unknown. Here, we used mixed stock analysis and examination of oceanic drift to elucidate movements of hawksbill turtles (Eretmochelys imbricata) and address management implications within the Caribbean. Among samples collected from 92 neritic juvenile hawksbills in the Cayman Islands we detected 11 mtDNA control region haplotypes. To estimate contributions to the aggregation, we performed ‘many‐to‐many’ mixed stock analysis, incorporating published hawksbill genetic and population data. The Cayman Islands aggregation represents a diverse mixed stock: potentially contributing source rookeries spanned the Caribbean basin, delineating a scale of recruitment of 200–2500 km. As hawksbills undergo an extended phase of oceanic dispersal, ocean currents may drive patterns of genetic diversity observed on foraging aggregations. Therefore, using high‐resolution Aviso ocean current data, we modelled movement of particles representing passively drifting oceanic juvenile hawksbills. Putative distribution patterns varied markedly by origin: particles from many rookeries were broadly distributed across the region, while others would appear to become entrained in local gyres. Overall, we detected a significant correlation between genetic profiles of foraging aggregations and patterns of particle distribution produced by a hatchling drift model (Mantel test, r = 0.77, P < 0.001; linear regression, r = 0.83, P < 0.001). Our results indicate that although there is a high degree of mixing across the Caribbean (a ‘turtle soup’), current patterns play a substantial role in determining genetic structure of foraging aggregations (forming turtle groups). Thus, for marine turtles and other widely distributed marine species, integration of genetic and oceanographic data may enhance understanding of population connectivity and management requirements.     

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.     

Reduced genetic diversity and increased reproductive isolation follow population‐level loss of larval dispersal in a marine gastropod

Population‐level consequences of dispersal ability remain poorly understood, especially for marine animals in which dispersal is typically considered a species‐level trait governed by oceanographic transport of microscopic larvae. Transitions from dispersive (planktotrophic) to nondispersive, aplanktonic larvae are predicted to reduce connectivity, genetic diversity within populations, and the spatial scale at which reproductive isolation evolves. However, larval dimorphism within a species is rare, precluding population‐level tests. We show the sea slug Costasiella ocellifera expresses both larval morphs in Florida and the Caribbean, regions with divergent mitochondrial lineages. Planktotrophy predominated at 11 sites, 10 of which formed a highly connected and genetically diverse Caribbean metapopulation. Four populations expressed mainly aplanktonic development and had markedly reduced connectivity, and lower genetic diversity at one mitochondrial and six nuclear loci. Aplanktonic dams showed partial postzygotic isolation in most interpopulation crosses, regardless of genetic or geographic distance to the sire's source, suggesting that outbreeding depression affects fragmented populations. Dams from genetically isolated and neighboring populations also exhibited premating isolation, consistent with reinforcement contingent on historical interaction. By increasing self‐recruitment and genetic drift, the loss of dispersal may thus initiate a feedback loop resulting in the evolution of reproductive isolation over small spatial scales in the sea.