Theoretical limits to the correlation between pelagic larval duration and population genetic structure

Daily age, growth rate, and pelagic larval duration of commercially important snapper species in Abrolhos National Marine Park

The age and daily growth of fish are registered through the deposition of increments in their otoliths, which are concretions formed by the precipitation of substances present in the endolymphatic fluid, mainly calcium carbonate (CaCO3). Faced with the need to fill some of the gaps in the knowledge on the occurrence and duration of the initial stages of snapper species' life cycles in the Abrolhos Bank, this study aimed to describe the growth rates, age, and period of pelagic larval duration (PLD) of three snapper species during the larval pre-settlement phase, in the Abrolhos Bank region. The post-larvae were captured using light traps. Otoliths were removed from 117 samples of snapper species; however, only 69 were viable for age estimation, of which 15 were Lutjanus analis, 25 were Lutjanus jocu, and 29 were Lutjanus synagris. Together, the samples presented individuals with total lengths ranging from 16.14 to 24.76 mm and ages from 21 to 39 days. Settlement marks were found for all three species, and the average PLD was ~25 days. The somatic growth of the snapper species was positively correlated with otolith growth. L. jocu presented the greatest daily growth compared to the other species. The three species use the Abrolhos Bank as a larval settlement site, demonstrating plasticity by using different habitats throughout their lives.

Demographic changes and life-history strategies predict the genetic diversity in crabs

Uncovering what predicts genetic diversity (GD) within species can help us access the status of populations and their evolutionary potential. Traits related to effective population size show a proportional association to GD, but evidence supports life-history strategies and habitat as the drivers of GD variation. Instead of investigating highly divergent taxa, focusing on one group could help to elucidate the factors influencing the GD. Additionally, most empirical data is based on vertebrate taxa; therefore, we might be missing novel patterns of GD found in neglected invertebrate groups. Here, we investigated the predictors of the GD in crabs (Brachyura) by compiling the most comprehensive cytochrome c oxidase subunit I (COI) available. Eight predictor variables were analysed across 150 species (16 992 sequences) using linear models (multiple linear regression) and comparative methods (PGLS). Our results indicate that population size fluctuation represents the most critical trait predicting GD, with species that have undergone bottlenecks followed by population expansion showing lower GD. Egg size, pelagic larval duration and habitat might play a role probably because of their association with how species respond to disturbances. Ultimately, K-strategists that have undergone bottlenecks are the species showing lower GD. Some variables do not show an association with GD as expected, most likely due to the taxon-specific role of some predictors, which should be considered in further investigations and generalizations. This work highlights the complexity underlying the predictors of GD and adds results from a marine invertebrate group to the current understanding of this topic.

High-throughput SNPs dataset reveal restricted population connectivity of marine gastropod within the narrow distribution range of peripheral oceanic islands

Molecular studies based on the high resolution genetic markers help us to grasp the factor shaping the genetic structure of marine organisms. Ecological factors linking to life history traits have often explained the process of genetic structuring in open and connectable oceanic environments. Besides, population genetic divergence can be affected by fragmented habitat, oceanic current, and past geographical events. In the present study, we demonstrated the genetic differentiation of marine gastropod Monodonta sp. within a narrow range of peripheral oceanic islands, the Ogasawara Islands. Genetic analyses were performed not only with a mitochondrial DNA marker but also with a high-throughput SNPs dataset obtained by ddRAD-seq. The results of the mtDNA analyses did not show genetic divergence among populations, while the SNPs dataset detected population genetic differentiation. Population demographic analyses and gene flow estimation suggested that the genetic structure was formed by sea level fluctuation associated with the past climatic change and regulated by temporal oceanographic conditions. These findings provide important insights into population genetic patterns in open and connectable environments.

The role of wind in controlling the connectivity of blue mussels (Mytilus edulis L.) populations

BACKGROUND

Larval connectivity between distinct benthic populations is essential for their persistence. Although connectivity is difficult to measure in situ, it can be predicted via models that simulate biophysical interactions between larval behaviour and ocean currents. The blue mussel (Mytilus Edulis L.) is widespread throughout the northern hemisphere and extensively commercialised worldwide. In the Irish Sea, this industry represents ~ 50% of Welsh shellfisheries, where cultivation is mainly based on wild spat. However, the main sources and amount of spat varied interannually (1100 tonnes harvest in 2014 against zero in 2018). The aim of this study is to characterise the structure and dynamics of the blue mussel metapopulation within the northern part of the Irish Sea.

METHODS

We develop a Lagrangian particle tracking model, driven by a high-resolution (from 30 to 5000 m) validated unstructured coastal hydrodynamic model of the Irish Sea, to simulate spatial and temporal variability of larval dispersal and connectivity between distinct mussel populations and potential settlement areas.

RESULTS

Our results showed that: (1) larvae positioned near the surface were strongly influenced by wind-driven currents suggesting that connectivity networks had the potential to span hundreds of kilometres; (2) in contrast, larvae positioned deeper in the water column were driven by tidal currents, producing intricate spatial patterns of connectivity between mussel beds over tens of kilometres that were consistent over time.

CONCLUSIONS

Dispersal of mussel larvae in the tidally energetic Irish Sea during the April-May spawning season is potentially driven by wind-driven surface currents, as confirmed by fisherman observations of inter-annual variability in wild spat collection. These results have important implications for metapopulation dynamics within the context of climate change and sustainable shellfisheries management (i.e. gain and loss of populations and harvest areas according to wind conditions).

Dongsha Atoll is an important stepping-stone that promotes regional genetic connectivity in the South China Sea

Background

Understanding region-wide patterns of larval connectivity and gene flow is crucial for managing and conserving marine biodiversity. Dongsha Atoll National Park (DANP), located in the northern South China Sea (SCS), was established in 2007 to study and conserve this diverse and remote coral atoll. However, the role of Dongsha Atoll in connectivity throughout the SCS is seldom studied. In this study, we aim to evaluate the role of DANP in conserving regional marine biodiversity.

Methods

In total, 206 samples across nine marine species were collected and sequenced from Dongsha Atoll, and these data were combined with available sequence data from each of these nine species archived in the Genomic Observatories Metadatabase (GEOME). Together, these data provide the most extensive population genetic analysis of a single marine protected area. We evaluate metapopulation structure for each species by using a coalescent sampler, selecting among panmixia, stepping-stone, and island models of connectivity in a likelihood-based framework. We then completed a heuristic graph theoretical analysis based on maximum dispersal distance to get a sense of Dongsha’s centrality within the SCS.

Results

Our dataset yielded 111 unique haplotypes across all taxa at DANP, 58% of which were not sampled elsewhere. Analysis of metapopulation structure showed that five out of nine species have strong regional connectivity across the SCS such that their gene pools are effectively panmictic (mean pelagic larval duration (PLD) = 78 days, sd = 60 days); while four species have stepping-stone metapopulation structure, indicating that larvae are exchanged primarily between nearby populations (mean PLD = 37 days, sd = 15 days). For all but one species, Dongsha was ranked within the top 15 out of 115 large reefs in the South China Sea for betweenness centrality. Thus, for most species, Dongsha Atoll provides an essential link for maintaining stepping-stone gene flow across the SCS.

Conclusions

This multispecies study provides the most comprehensive examination of the role of Dongsha Atoll in marine connectivity in the South China Sea to date. Combining new and existing population genetic data for nine coral reef species in the region with a graph theoretical analysis, this study provides evidence that Dongsha Atoll is an important hub for sustaining connectivity for the majority of coral-reef species in the region.

Weak biodiversity connectivity in the European network of no-take marine protected areas

The need for international cooperation in marine resource management and conservation has been reflected in the increasing number of agreements aiming for effective and well-connected networks of Marine Protected Areas (MPAs). However, the extent to which individual MPAs are connected remains mostly unknown. Here, we use a biophysical model tuned with empirical data on species dispersal ecology to predict connectivity of a vast spectrum of biodiversity in the European network of marine reserves (i.e., no-take MPAs). Our results highlight the correlation between empirical propagule duration data and connectivity potential and show weak network connectivity and strong isolation for major ecological groups, resulting from the lack of direct connectivity corridors between reserves over vast regions. The particularly high isolation predicted for ecosystem structuring species (e.g., corals, sponges, macroalgae and seagrass) might potentially undermine biodiversity conservation efforts if local retention is insufficient and unmanaged populations are at risk. Isolation might also be problematic for populations' persistence in the light of climate change and expected species range shifts. Our findings provide novel insights for management directives, highlighting the location of regions requiring additional marine reserves to function as stepping-stone connectivity corridors.

Life history, climate and biogeography interactively affect worldwide genetic diversity of plant and animal populations

Understanding how biological and environmental factors interactively shape the global distribution of plant and animal genetic diversity is fundamental to biodiversity conservation. Genetic diversity measured in local populations (GDP) is correspondingly assumed representative for population fitness and eco-evolutionary dynamics. For 8356 populations across the globe, we report that plants systematically display much lower GDP than animals, and that life history traits shape GDP patterns both directly (animal longevity and size), and indirectly by mediating core-periphery patterns (animal fecundity and plant dispersal). Particularly in some plant groups, peripheral populations can sustain similar GDP as core populations, emphasizing their potential conservation value. We further find surprisingly weak support for general latitudinal GDP trends. Finally, contemporary rather than past climate contributes to the spatial distribution of GDP, suggesting that contemporary environmental changes affect global patterns of GDP. Our findings generate new perspectives for the conservation of genetic resources at worldwide and taxonomic-wide scales.

Range-Wide Population Structure of 3 Deepwater Eteline Snappers Across the Indo-Pacific Basin

Deep-sea habitats may drive unique dispersal and demographic patterns for fishes, but population genetic analyses to address these questions have rarely been conducted for fishes in these environments. This study investigates the population structure of 3 tropical deepwater snappers of the genus Etelis that reside at 100-400 m depth, with broad and overlapping distributions in the Indo-Pacific. Previous studies showed little population structure within the Hawaiian Archipelago for 2 of these species: Etelis coruscans and E. carbunculus. Here we extend sampling to the entire geographic range of each species to resolve the population genetic architecture for these 2 species, as well as a recently exposed cryptic species (Etelis sp.). One goal was to determine whether deepwater snappers are more dispersive than shallow-water fishes. A second goal was to determine whether submesophotic fishes have older, more stable populations than shallow reef denizens that are subject to glacial sea-level fluctuations. Both goals are pertinent to the management of these valuable food fishes. A total of 1153 specimens of E. coruscans from 15 geographic regions were analyzed, along with 1064 specimens of E. carbunculus from 11 regions, and 590 specimens of E. sp. from 16 regions. The first 2 species were analyzed with mtDNA and 9-11 microsatellite loci, while E. sp. was analyzed with mtDNA only. Etelis coruscans had a non-significant microsatellite global FST, but significant global mtDNA Ф ST = 0.010 (P = 0.0007), with the isolation of Seychelles in the western Indian Ocean, and intermittent signals of isolation for the Hawaiian Archipelago. Etelis carbunculus had a non-significant microsatellite global FST, and significant global mtDNA Ф ST = 0.021 (P = 0.0001), with low but significant levels of isolation for Hawai'i, and divergence between Tonga and Fiji. Etelis sp. had mtDNA Ф ST = 0.018 (P = 0.0005), with a strong pattern of isolation for both Seychelles and Tonga. Overall, we observed low population structure, shallow mtDNA coalescence (similar to near-shore species), and isolation at the fringes of the Indo-Pacific basin in Hawai'i and the western Indian Ocean. While most shallow-water species have population structure on the scale of biogeographic provinces, deepwater snapper populations are structured on the wider scale of ocean basins, more similar to pelagic fishes than to shallow-water species. This population structure indicates the capacity for widespread dispersal throughout the Indo-Pacific region.

Against all odds: a tale of marine range expansion with maintenance of extremely high genetic diversity

The displacement of species from equatorial latitudes to temperate locations following the increase in sea surface temperatures is among the significant reported consequences of climate change. Shifts in the distributional ranges of species result in fish communities tropicalisation, i.e., high latitude colonisations by typically low latitude distribution species. These movements create new interactions between species and new trophic assemblages. The Senegal seabream, Diplodus bellottii, may be used as a model to understand the population genetics of these invasions. In the last decades, this species has undergone an outstanding range expansion from its African area of origin to the Atlantic coast of the Iberian Peninsula, where now occurs abundantly. Mitochondrial and nuclear markers revealed a striking high haplotypic nucleotide and genetic diversity values, along with significant population differentiation throughout the present-day geographical range of the Senegal seabream. These results are not consistent with the central-marginal hypothesis, nor with the expectations of a leptokurtic distribution of individuals, as D. bellottii seems to be able to retain exceptional levels of diversity in marginal and recently colonised areas. We discuss possible causes for hyperdiversity and lack of geographical structure and subsequent implications for fisheries.

Highly regional population structure of Spondyliosoma cantharus depicted by nuclear and mitochondrial DNA data

Resolution of population structure represents an effective way to define biological stocks and inform efficient fisheries management. In the present study, the phylogeography of the protogynous sparid Spondyliosoma cantharus, in the East Atlantic and Mediterranean Sea, was investigated with nuclear (S7) and mitochondrial (cytochrome b) DNA markers. Significant divergence of four regional genetic groups was observed: North Eastern Atlantic, Mediterranean Sea, Western African Transition (Cape Verde) and Gulf of Guinea (Angola). The two southern populations (Cape Verde and Angola) each comprised reciprocally monophyletic mtDNA lineages, revealed low levels of diversity in Cape Verde and high diversity for Angola despite being represented by only 14 individuals. A complete divergence between North Atlantic and Mediterranean populations was depicted by the mitochondrial marker, but a highly shared nuclear haplotype revealed an incomplete lineage sorting between these regions. Bayesian skyline plots and associated statistics revealed different dynamics among the four regions. Cape Verde showed no expansion and the expansion time estimated for Angola was much older than for the other regions. Mediterranean region seems to have experienced an early population growth but has remained with a stable population size for the last 30000 years while the North Atlantic population has been steadily growing. The lack of genetic structuring within these regions should not be taken as evidence of demographic panmixia in light of potential resolution thresholds and previous evidence of intra-regional phenotypic heterogeneity.

Dwarf brooder versus giant broadcaster: combining genetic and reproductive data to unravel cryptic diversity in an Antarctic brittle star

Poecilogony, or multiple developmental modes in a single species, is exceedingly rare. Several species described as poecilogenous were later demonstrated to be multiple (cryptic) species with a different developmental mode. The Southern Ocean is known to harbor a high proportion of brooders (Thorson's Rule) but with an increasing number of counter examples over recent years. Here we evaluated poecilogony vs. crypticism in the brittle star Astrotoma agassizii across the Southern Ocean. This species was initially described from South America as a brooder before some pelagic stages were identified in Antarctica. Reproductive and mitochondrial data were combined to unravel geographic and genetic variation of developmental modes. Our results indicate that A. agassizii is composed of seven well-supported and deeply divergent clades (I: Antarctica and South Georgia; II: South Georgia and Sub-Antarctic locations including Kerguelen, Patagonian shelf, and New Zealand; III-VI-VII: Patagonian shelf, IV-V: South Georgia). Two of these clades demonstrated strong size dimorphism when in sympatry and can be linked to differing developmental modes (Clade V: dwarf brooder vs. Clade I: giant broadcaster). Based on their restricted geographic distributions and on previous studies, it is likely that Clades III-VI-VII are brooders. Clade II is composed of different morphological species, A. agassizii and A. drachi, the latter originally used as the outgroup. By integrating morphology, reproductive, and molecular data we conclude that the variation identified in A. agassizii is best described as crypticism rather than poecilogony.

Do differences in mating behaviour lead to differences in connectivity patterns of reef fishes? Insights from two sympatric surgeonfish species in the Indian Ocean

Disentangling the contribution of biotic and abiotic factors in the structuring of the genetic diversity of reef species is critical to illuminate the diversification of evolutionary lineages in marine environment. Howevr, previous studies have mainly focused on determining the influence of pelagic larval duration on the connectivity and demographic history of reef fishes, whereas few studies have examined the effects of other biotic factors, such as mating behaviour and habitat preference. Here, we use mitochondrial DNA (ATPase 6/8) and ten microsatellite loci to compare the population genetic structure and demographic history of the spawning aggregating Acanthurus triostegus with the monogamous spawning Acanthurus leucosternon. Pairwise comparisons and discriminant analysis of principal components showed that the genetic structuring patterns of the two species are not consistent with the influence of mating behaviour, suggesting the possible role of other biotic and abiotic factors. However, demographic history estimates revealed that these species may have responded differently to sea level fluctuations during the glacial maxima.

A coalescent sampler successfully detects biologically meaningful population structure overlooked by F-statistics

Assessing the geographic structure of populations has relied heavily on Sewell Wright's F-statistics and their numerous analogues for many decades. However, it is well appreciated that, due to their nonlinear relationship with gene flow, F-statistics frequently fail to reject the null model of panmixia in species with relatively high levels of gene flow and large population sizes. Coalescent genealogy samplers instead allow a model-selection approach to the characterization of population structure, thereby providing the opportunity for stronger inference. Here, we validate the use of coalescent samplers in a high gene flow context using simulations of a stepping-stone model. In an example case study, we then re-analyze genetic datasets from 41 marine species sampled from throughout the Hawaiian archipelago using coalescent model selection. Due to the archipelago's linear nature, it is expected that most species will conform to some sort of stepping-stone model (leading to an expected pattern of isolation by distance), but F-statistics have only supported this inference in ~10% of these datasets. Our simulation analysis shows that a coalescent sampler can make a correct inference of stepping-stone gene flow in nearly 100% of cases where gene flow is ≤100 migrants per generation (equivalent to F ST = 0.002), while F-statistics had mixed results. Our re-analysis of empirical datasets found that nearly 70% of datasets with an unambiguous result fit a stepping-stone model with varying population sizes and rates of gene flow, although 37% of datasets yielded ambiguous results. Together, our results demonstrate that coalescent samplers hold great promise for detecting weak but meaningful population structure, and defining appropriate management units.

Hitchhiking the high seas: Global genomics of rafting crabs

Population differentiation and diversification depend in large part on the ability and propensity of organisms to successfully disperse. However, our understanding of these processes in organisms with high dispersal ability is biased by the limited genetic resolution offered by traditional genotypic markers. Many neustonic animals disperse not only as pelagic larvae, but also as juveniles and adults while drifting or rafting at the surface of the open ocean. In theory, the heightened dispersal ability of these animals should limit opportunities for species diversification and population differentiation. To test these predictions, we used next-generation sequencing of genomewide restriction-site-associated DNA tags (RADseq) and traditional mitochondrial DNA sequencing, to investigate the species-level relationships and global population structure of Planes crabs collected from oceanic flotsam and sea turtles. Our results indicate that species diversity in this clade is low-likely three closely related species-with no evidence of cryptic or undescribed species. Moreover, our results indicate weak population differentiation among widely separated aggregations with genetic indices showing only subtle genetic discontinuities across all oceans of the world (RADseq F ST = 0.08-0.16). The results of this study provide unprecedented resolution of the systematics and global biogeography of this group and contribute valuable information to our understanding of how theoretical dispersal potential relates to actual population differentiation and diversification among marine organisms. Moreover, these results demonstrate the limitations of single gene analyses and the value of genomic-level resolution for estimating contemporary population structure in organisms with large, highly connected populations.

Stepping stones to isolation: Impacts of a changing climate on the connectivity of fragmented fish populations

In the marine environment, understanding the biophysical mechanisms that drive variability in larval dispersal and population connectivity is essential for estimating the potential impacts of climate change on the resilience and genetic structure of populations. Species whose populations are small, isolated and discontinuous in distribution will differ fundamentally in their response and resilience to environmental stress, compared with species that are broadly distributed, abundant and frequently exchange conspecifics. Here, we use an individual-based modelling approach, combined with a population genetics projection model, to consider the impacts of a warming climate on the population connectivity of two contrasting Antarctic fish species, Notothenia rossii and Champsocephalus gunnari. Focussing on the Scotia Sea region, sea surface temperatures are predicted to increase significantly by the end of the 21st century, resulting in reduced planktonic duration and increased egg and larval mortality. With shorter planktonic durations, the results of our study predict reduced dispersal of both species across the Scotia Sea, from Antarctic Peninsula sites to islands in the north and east, and increased dispersal among neighbouring sites, such as around the Antarctic Peninsula. Increased mortality modified the magnitude of population connectivity but had little effect on the overall patterns. Whilst the predicted changes in connectivity had little impact on the projected regional population genetic structure of N. rossii, which remained broadly genetically homogeneous within distances of ~1,500 km, the genetic isolation of C. gunnari populations in the northern Scotia Sea was predicted to increase with rising sea temperatures. Our study highlights the potential for increased isolation of island populations in a warming world, with implications for the resilience of populations and their ability to adapt to ongoing environmental change, a matter of high relevance to fisheries and ecosystem-level management.

Pelagic larval duration, growth rate, and population genetic structure of the tidepool snake moray Uropterygius micropterus around the southern Ryukyu Islands, Taiwan, and the central Philippines

The relationships between pelagic larval duration (PLD) and geographic distribution patterns or population genetic structures of fishes remain obscure and highly variable among species. To further understand the early life history of the tidepool snake moray Uropterygius micropterus and the potential relationship between PLD and population genetic structure of this species, otolith microstructure and population genetics based on concatenated mtDNA sequence (cytochrome b and cytochrome oxidase subunit I, 1,336 bp) were analyzed for 195 specimens collected from eight locations around the southern Ryukyu Islands, Taiwan, and the central Philippines. Eels with longer PLD and lower otolith growth rates were observed at relatively higher latitudes with lower water temperatures (54.6 ± 7.7 days and 1.28 ± 0.16 µm day^-1 on Ishigaki Island, Japan, vs. 43.9 ± 4.9 days and 1.60 ± 0.19 µm day^-1 on Badian, the Philippines), suggesting that leptocephali grew faster and had shortened pelagic periods in warmer waters. Meanwhile, the eels along the southwest coast of Taiwan had relatively longer PLD (57.9 ± 10.5 days), which might be associated with the more complex ocean current systems compared to their counterparts collected along the east coast of Taiwan (52.6 ± 8.0 days). However, the southwestern and eastern Taiwan groups had similar otolith growth rates (1.33 ± 0.19 µm day^-1 vs. 1.36 ± 0.16 µm day^-1). Despite the intergroup variation in PLD, genetic analysis revealed fluent gene flow among the tidepool snake morays in the study regions, implying that intraspecies PLD variation had a weak effect on genetic structure. The leptocephalus stage might have ensured the widespread gene flow among the study areas and leptocephalus growth was likely influenced by regional water temperature.

Chasing genetic structure in coralligenous reef invertebrates: patterns, criticalities and conservation issues

Conservation of coastal habitats is a global issue, yet biogenic reefs in temperate regions have received very little attention. They have a broad geographic distribution and are a key habitat in marine ecosystems impacted by human activities. In the Mediterranean Sea coralligenous reefs are biodiversity hot spots and are classified as sensitive habitats deserving conservation. Genetic diversity and structure influence demographic, ecological and evolutionary processes in populations and play a crucial role in conservation strategies. Nevertheless, a comprehensive view of population genetic structure of coralligenous species is lacking. Here, we reviewed the literature on the genetic structure of sessile and sedentary invertebrates of the Mediterranean coralligenous reefs. Linear regression models and meta-analytic approaches are used to assess the contributions of genetic markers, phylum, pelagic larval duration (PLD) and geographical distance to the population genetic structure. Our quantitative approach highlight that 1) most species show a significant genetic structure, 2) structuring differs between phyla, and 3) PLD does not appear to be a major driver of the structuring. We discuss the implication of these finding for the management and conservation, suggesting research areas that deserve attention, and providing recommendations for broad assessment and monitoring of genetic diversity in biogenic reefs species.

Strong population structure but no equilibrium yet: Genetic connectivity and phylogeography in the kelp Saccharina latissima (Laminariales, Phaeophyta)

Kelp aquaculture is globally developing steadily as human food source, along with other applications. One of the newer crop species is Saccharina latissima, a northern hemisphere kelp inhabiting temperate to arctic rocky shores. To protect and document its natural genetic variation at the onset of this novel aquaculture, as well as increase knowledge on its taxonomy and phylogeography, we collected new genetic data, both nuclear and mitochondrial, and combined it with previous knowledge to estimate genetic connectivity and infer colonization history. Isolation‐with‐migration coalescent analyses demonstrate that gene flow among the sampled locations is virtually nonexistent. An updated scenario for the origin and colonization history of S. latissima is developed as follows: We propose that the species (or species complex) originated in the northwest Pacific, crossed to the northeast Pacific in the Miocene, and then crossed the Bering Strait after its opening ~5.5 Ma into the Arctic and northeast Atlantic. It subsequently crossed the Atlantic from east to west. During the Pleistocene, it was compressed in the south with evidence for northern refugia in Europe. Postglacial recolonization led to secondary contact in the Canadian Arctic. Saccharina cichorioides is shown to probably belong to the S. latissima species complex and to derive from ancestral populations in the Asian North Pacific. Our novel approach of comparing inferred gene flow based on coalescent analysis versus Wright's island model suggests that equilibrium levels of differentiation have not yet been reached in Europe and, hence, that genetic differentiation is expected to increase further if populations are left undisturbed.

Three-year monitoring of genetic diversity reveals a micro-connectivity pattern and local recruitment in the broadcast marine species Paracentrotus lividus

In conservation and management of marine biological resources, a knowledge of connectivity is necessary to understand how local populations are naturally replenished by the arrival of new recruits from source populations. At small geographical scales, species experiencing moderate to long pelagic larval phases are mostly genetically homogeneous, which hinders inferences about local connectivity. Recent studies demonstrated that assessing genetic relatedness and kinship could provide information about local connectivity in populations with high levels of gene flow. Here, we were interested in deciphering the structure and connectivity of populations of the sea urchin Paracentrotus lividus, by monitoring populations at 11 localities distributed along a 225-km coast-line in the south-eastern French Mediterranean Sea. Using 12 microsatellite loci, we found a weak but significant genetic differentiation and observed a transient genetic differentiation among locations within temporal cohorts, without any correlation with the distance between locations, interpreted as unexplainable chaotic genetic patchiness. Among temporal cohorts, the more related individuals were mainly found within locations and the observed local differentiation (F_ST) correlated with the proportion of kin within locations, suggesting that larvae dispersed cohesively. Specifically, we could also reveal that populations flanking Cape Sicié were influenced by eastern populations and that local recruitment was a frequent occurrence. Overall, our results contribute to the growing number of studies showing that connectivity can be reliably assessed at a fine spatial scale even in genetically homogenous populations.

A multispecies approach reveals hot spots and cold spots of diversity and connectivity in invertebrate species with contrasting dispersal modes

Genetic diversity is crucial for species' maintenance and persistence, yet is often overlooked in conservation studies. Species diversity is more often reported due to practical constraints, but it is unknown if these measures of diversity are correlated. In marine invertebrates, adults are often sessile or sedentary and populations exchange genes via dispersal of gametes and larvae. Species with a larval period are expected to have more connected populations than those without larval dispersal. We assessed the relationship between measures of species and genetic diversity, and between dispersal ability and connectivity. We compiled data on genetic patterns and life history traits in nine species across five phyla. Sampling sites spanned 600 km in the northwest Mediterranean Sea and focused on a 50-km area near Marseilles, France. Comparative population genetic approaches yielded three main results. (i) Species without larvae showed higher levels of genetic structure than species with free-living larvae, but the role of larval type (lecithotrophic or planktotrophic) was negligible. (ii) A narrow area around Marseilles, subject to offshore advection, limited genetic connectivity in most species. (iii) We identified sites with significant positive contributions to overall genetic diversity across all species, corresponding with areas near low human population densities. In contrast, high levels of human activity corresponded with a negative contribution to overall genetic diversity. Genetic diversity within species was positively and significantly linearly related to local species diversity. Our study suggests that local contribution to overall genetic diversity should be taken into account for future conservation strategies.

Phylogeography and population structure of - grypotus (Richardson, 1846) as revealed by mitochondrial control region sequences

The 137 individuals of Johniusgrypotus were collected from seven localities from the Bohai Sea to the East China Sea. A 549 base pair (bp) fragment of the hypervariable region of the mtDNA control region was sequenced to examine genetic diversity and population structure. The populations of J.grypotus showed high haplotype diversity (h) with a range from 0.7500 to 0.9740 and low nucleotide diversity (π) with a range from 0.0024 to 0.0067. Low and non-significant genetic differentiation was estimated among populations except for North Yellow Sea population, which has a significant genetic difference with other populations. The demographic history examined by mismatch distribution analyses and Bayesian skyline plot (BSP) analyses revealed that a sudden population expansion occurred almost 20 to 40 thousand years before. Relatively recent population expansion in the last glacial period, large dispersal of eggs or larvae carried by coastal current, and the homogeneity of living environment may have an important influence on the population genetic pattern.

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).

Community assembly of coral reef fishes along the Melanesian biodiversity gradient

The Indo-Pacific is home to Earth’s most biodiverse coral reefs. Diversity on these reefs decreases from the Coral Triangle east through the islands of Melanesia. Despite this pattern having been identified during the early 20^th century, our knowledge about the interaction between pattern and process remains incomplete. To evaluate the structure of coral reef fish communities across Melanesia, we obtained distributional records for 396 reef fish species in five taxa across seven countries. We used hierarchical clustering, nestedness, and multiple linear regression analyses to evaluate the community structure. We also compiled data on life history traits (pelagic larval duration, body size and schooling behavior) to help elucidate the ecological mechanisms behind community structure. Species richness for these taxa along the gradient was significantly related to longitude but not habitat area. Communities are significantly nested, indicating that species-poor communities are largely composed of subsets of the species found on species rich reefs. These trends are robust across taxonomic groups except for the Pomacentridae, which exhibit an anti-nested pattern, perhaps due to a large number of endemic species. Correlations between life history traits and the number of reefs on which species occurred indicate that dispersal and survival ability contribute to determining community structure. We conclude that distance from the Coral Triangle dominates community structure in reef fish; however, conservation of the most species-rich areas will not be sufficient alone to conserve the vivid splendor of this region.

Integrative taxonomy of the ornamental 'peppermint' shrimp public market and population genetics of Lysmata boggessi, the most heavily traded species worldwide

The ornamental trade is a worldwide industry worth >15 billion USD with a problem of rampant product misidentification. Minimizing misidentification is critical in the face of overexploitation of species in the trade. We surveyed the peppermint shrimp ornamental marketplace in the southeastern USA, the most intense market for peppermint shrimps worldwide, to characterize the composition of species in the trade, reveal the extent of misidentification, and describe the population genetics of the true target species. Shrimps were bought from aquarium shops in FL, GA, SC, and NC. We demonstrated, contrary to popular belief (information from dealers), that the most heavily traded species in the market was Lysmata boggessi, an endemic species to the eastern Gulf of Mexico, and not Lysmata wurdemanni. Importantly, only when color pattern or genetic markers in conjunction with morphological traits were employed, was it was possible to unequivocally identify L. boggessi as the only species in the trade. The intensity of the market for peppermint shrimps in the USA has led to L. boggessi being the most traded species worldwide. Misidentification in the shrimp aquarium trade is accidental and involuntary, and is explained by remarkable similarity among congeneric species. Using sequences of the 16S-mt-DNA marker, we found no indication of population genetic structure in the endemic L. boggessi across  550 km of linear coast. Therefore, this species can be considered genetically homogeneous and a single fished stock. Still, we argue in favor of additional studies using more powerful markers (e.g., SNPs) capable of revealing genetic structure at a finer spatial-scale. Our results will help advance management and conservation policies in this lucrative yet understudied fishery. Future studies of other ornamental fisheries will benefit from using an integrative taxonomic approach, as we demonstrate here.

Population differentiation or species formation across the Indian and the Pacific Oceans? An example from the brooding marine hydrozoan Macrorhynchia phoenicea

Assessing population connectivity is necessary to construct effective marine protected areas. This connectivity depends, among other parameters, inherently on species dispersal capacities. Isolation by distance (IBD) is one of the main modes of differentiation in marine species, above all in species presenting low dispersal abilities. This study reports the genetic structuring in the tropical hydrozoan Macrorhynchia phoenicea α (sensu Postaire et al., 2016a), a brooding species, from 30 sampling sites in the Western Indian Ocean and the Tropical Southwestern Pacific, using 15 microsatellite loci. At the local scale, genet dispersal relied on asexual propagation at short distance, which was not found at larger scales. Considering one representative per clone, significant positive F_IS values (from −0.327*** to 0.411***) were found within almost all sites. Gene flow was extremely low at all spatial scales, among sites within islands (<10 km distance) and among islands (100 to >11,000 km distance), with significant pairwise F_ST values (from 0.035*** to 0.645***). A general pattern of IBD was found at the Indo‐Pacific scale, but also within ecoregions in the Western Indian Ocean province. Clustering and network analyses identified each island as a potential independent population, while analysis of molecular variance indicated that population genetic differentiation was significant at small (within island) and intermediate (among islands within province) spatial scales. As shown by this species, a brooding life cycle might be corollary of the high population differentiation found in some coastal marine species, thwarting regular dispersal at distances more than a few kilometers and probably leading to high cryptic diversity, each island housing independent evolutionary lineages.

Invertebrate population genetics across Earth's largest habitat: The deep-sea floor

Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non-chemosynthetic ecosystems on the deep-sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow-water species. Generally, populations at similar depths were well connected over 100s-1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s-1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean-wide (under 4%), and 48% were Atlantic-focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, "ecosystem engineers" and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single-locus mitochondrial genes revealing a common pattern of non-neutrality, consistent with demographic instability or selective sweeps; similar to deep-sea hydrothermal vent fauna. The absence of a clear difference between vent and non-vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single-locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta-analyses where broad inferences about deep-sea ecology could be made.

Historical gene flow constraints in a northeastern Atlantic fish: phylogeography of the ballan wrasse Labrus bergylta across its distribution range

The distribution and demographic patterns of marine organisms in the north Atlantic were largely shaped by climatic changes during the Pleistocene, when recurrent glacial maxima forced them to move south or to survive in northern peri-glacial refugia. These patterns were also influenced by biological and ecological factors intrinsic to each species, namely their dispersion ability. The ballan wrasse (Labrus bergylta), the largest labrid fish along Europe's continental margins, is a target for fisheries and aquaculture industry. The phylogeographic pattern, population structure, potential glacial refugia and recolonization routes for this species were assessed across its full distribution range, using mitochondrial and nuclear markers. The existence of a marked population structure can reflect both recolonization from three distinct glacial refugia and current and past oceanographic circulation patterns. Although isolated in present times, shared haplotypes between continental and Azores populations and historical exchange of migrants in both directions point to a common origin of L. bergylta. This situation is likely to be maintained and/or accentuated by current circulation patterns in the north Atlantic, and may lead to incipient speciation in the already distinct Azorean population. Future monitoring of this species is crucial to evaluate how this species is coping with current environmental changes.

Navigating the currents of seascape genomics: how spatial analyses can augment population genomic studies

Population genomic approaches are making rapid inroads in the study of non-model organisms, including marine taxa. To date, these marine studies have predominantly focused on rudimentary metrics describing the spatial and environmental context of their study region (e.g., geographical distance, average sea surface temperature, average salinity). We contend that a more nuanced and considered approach to quantifying seascape dynamics and patterns can strengthen population genomic investigations and help identify spatial, temporal, and environmental factors associated with differing selective regimes or demographic histories. Nevertheless, approaches for quantifying marine landscapes are complicated. Characteristic features of the marine environment, including pelagic living in flowing water (experienced by most marine taxa at some point in their life cycle), require a well-designed spatial-temporal sampling strategy and analysis. Many genetic summary statistics used to describe populations may be inappropriate for marine species with large population sizes, large species ranges, stochastic recruitment, and asymmetrical gene flow. Finally, statistical approaches for testing associations between seascapes and population genomic patterns are still maturing with no single approach able to capture all relevant considerations. None of these issues are completely unique to marine systems and therefore similar issues and solutions will be shared for many organisms regardless of habitat. Here, we outline goals and spatial approaches for landscape genomics with an emphasis on marine systems and review the growing empirical literature on seascape genomics. We review established tools and approaches and highlight promising new strategies to overcome select issues including a strategy to spatially optimize sampling. Despite the many challenges, we argue that marine systems may be especially well suited for identifying candidate genomic regions under environmentally mediated selection and that seascape genomic approaches are especially useful for identifying robust locus-by-environment associations.

Beyond the Coral Triangle: high genetic diversity and near panmixia in Singapore's populations of the broadcast spawning sea star Protoreaster nodosus

The Coral Triangle is widely considered the most important centre of marine biodiversity in Asia while areas on its periphery such as the South China Sea, have received much less interest. Here, we demonstrate that a small population of the knobbly sea star Protoreaster nodosus in Singapore has similarly high levels of genetic diversity as comparable Indonesian populations from the Coral Triangle. The high genetic diversity of this population is remarkable because it is maintained despite decades of continued anthropogenic disturbance. We postulate that it is probably due to broadcast spawning which is likely to maintain high levels of population connectivity. To test this, we analysed 6140 genome-wide single nucleotide polymorphism (SNP) loci for Singapore's populations and demonstrate a pattern of near panmixia. We here document a second case of high genetic diversity and low genetic structure for a broadcast spawner in Singapore, which suggests that such species have high resilience against anthropogenic disturbances. The study demonstrates the feasibility and power of using genome-wide SNPs for connectivity studies of marine invertebrates without a sequenced genome.

Connectivity in the cold: the comparative population genetics of vent-endemic fauna in the Scotia Sea, Southern Ocean

We report the first comparative population genetics study for vent fauna in the Southern Ocean using cytochrome C oxidase I and microsatellite markers. Three species are examined: the kiwaid squat lobster, Kiwa tyleri, the peltospirid gastropod, Gigantopelta chessoia, and a lepetodrilid limpet, Lepetodrilus sp., collected from vent fields 440 km apart on the East Scotia Ridge (ESR) and from the Kemp Caldera on the South Sandwich Island Arc, ~95 km eastwards. We report no differentiation for all species across the ESR, consistent with panmixia or recent range expansions. A lack of differentiation is notable for Kiwa tyleri, which exhibits extremely abbreviated lecithotrophic larval development, suggestive of a very limited dispersal range. Larval lifespans may, however, be extended by low temperature-induced metabolic rate reduction in the Southern Ocean, muting the impact of dispersal strategy on patterns of population structure. COI diversity patterns suggest all species experienced demographic bottlenecks or selective sweeps in the past million years and possibly at different times. ESR and Kemp limpets are divergent, although with evidence of very recent ESR-Kemp immigration. Their divergence, possibility indicative of incipient speciation, along with the absence of the other two species at Kemp, may be the consequence of differing dispersal capabilities across a ~1000 m depth range and/or different selective regimes between the two areas. Estimates of historic and recent limpet gene flow between the ESR and Kemp are consistent with predominantly easterly currents and potentially therefore, cross-axis currents on the ESR, with biogeographic implications for the region.

A synthesis of genetic connectivity in deep-sea fauna and implications for marine reserve design

With anthropogenic impacts rapidly advancing into deeper waters, there is growing interest in establishing deep-sea marine protected areas (MPAs) or reserves. Reserve design depends on estimates of connectivity and scales of dispersal for the taxa of interest. Deep-sea taxa are hypothesized to disperse greater distances than shallow-water taxa, which implies that reserves would need to be larger in size and networks could be more widely spaced; however, this paradigm has not been tested. We compiled population genetic studies of deep-sea fauna and estimated dispersal distances for 51 studies using a method based on isolation-by-distance slopes. Estimates of dispersal distance ranged from 0.24 km to 2028 km with a geometric mean of 33.2 km and differed in relation to taxonomic and life-history factors as well as several study parameters. Dispersal distances were generally greater for fishes than invertebrates with the Mollusca being the least dispersive sampled phylum. Species that are pelagic as adults were more dispersive than those with sessile or sedentary lifestyles. Benthic species from soft-substrate habitats were generally less dispersive than species from hard substrate, demersal or pelagic habitats. As expected, species with pelagic and/or feeding (planktotrophic) larvae were more dispersive than other larval types. Many of these comparisons were confounded by taxonomic or other life-history differences (e.g. fishes being more dispersive than invertebrates) making any simple interpretation difficult. Our results provide the first rough estimate of the range of dispersal distances in the deep sea and allow comparisons to shallow-water assemblages. Overall, dispersal distances were greater for deeper taxa, although the differences were not large (0.3-0.6 orders of magnitude between means), and imbalanced sampling of shallow and deep taxa complicates any simple interpretation. Our analyses suggest the scales of dispersal and connectivity for reserve design in the deep sea might be comparable to or slightly larger than those in shallow water. Deep-sea reserve design will need to consider the enormous variety of taxa, life histories, hydrodynamics, spatial configuration of habitats and patterns of species distributions. The many caveats of our analyses provide a strong impetus for substantial future efforts to assess connectivity of deep-sea species from a variety of habitats, taxonomic groups and depth zones.

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.

Oceanography and life history predict contrasting genetic population structure in two Antarctic fish species

Understanding the key drivers of population connectivity in the marine environment is essential for the effective management of natural resources. Although several different approaches to evaluating connectivity have been used, they are rarely integrated quantitatively. Here, we use a 'seascape genetics' approach, by combining oceanographic modelling and microsatellite analyses, to understand the dominant influences on the population genetic structure of two Antarctic fishes with contrasting life histories, Champsocephalus gunnari and Notothenia rossii. The close accord between the model projections and empirical genetic structure demonstrated that passive dispersal during the planktonic early life stages is the dominant influence on patterns and extent of genetic structuring in both species. The shorter planktonic phase of C. gunnari restricts direct transport of larvae between distant populations, leading to stronger regional differentiation. By contrast, geographic distance did not affect differentiation in N. rossii, whose longer larval period promotes long-distance dispersal. Interannual variability in oceanographic flows strongly influenced the projected genetic structure, suggesting that shifts in circulation patterns due to climate change are likely to impact future genetic connectivity and opportunities for local adaptation, resilience and recovery from perturbations. Further development of realistic climate models is required to fully assess such potential impacts.

Phylogeography of the California sheephead, Semicossyphus pulcher: the role of deep reefs as stepping stones and pathways to antitropicality

In the past decade, the study of dispersal of marine organisms has shifted from focusing predominantly on the larval stage to a recent interest in adult movement. Antitropical distributions provide a unique system to assess vagility and dispersal. In this study, we have focused on an antitropical wrasse genus, Semicossyphus, which includes the California sheephead, S. pulcher, and Darwin's sheephead, S. darwini. Using a phylogenetic approach based on mitochondrial and nuclear markers, and a population genetic approach based on mitochondrial control region sequences and 10 microsatellite loci, we compared the phylogenetic relationships of these two species, as well as the population genetic characteristics within S. pulcher. While S. pulcher and S. darwini are found in the temperate eastern Pacific regions of the northern and southern hemispheres, respectively, their genetic divergence was very small (estimated to have occurred between 200 and 600 kya). Within S. pulcher, genetic structuring was generally weak, especially along mainland California, but showed weak differentiation between Sea of Cortez and California, and between mainland California and Channel Islands. We highlight the congruence of weak genetic differentiation both within and between species and discuss possible causes for maintenance of high gene flow. In particular, we argue that deep and cooler water refugia are used as stepping stones to connect distant populations, resulting in low levels of genetic differentiation.

Speciation in fishes

The field of speciation has seen much renewed interest in the past few years, with theoretical and empirical advances that have moved it from a descriptive field to a predictive and testable one. The goal of this review is to provide a general background on research on speciation as it pertains to fishes. Three major components to the question are first discussed: the spatial, ecological and sexual factors that influence speciation mechanisms. We then move to the latest developments in the field of speciation genomics. Affordable and rapidly available, massively parallel sequencing data allow speciation studies to converge into a single comprehensive line of investigation, where the focus has shifted to the search for speciation genes and genomic islands of speciation. We argue that fish present a very diverse array of scenarios, making them an ideal model to study speciation processes.

Combined analyses of kinship and FST suggest potential drivers of chaotic genetic patchiness in high gene-flow populations

We combine kinship estimates with traditional F-statistics to explain contemporary drivers of population genetic differentiation despite high gene flow. We investigate range-wide population genetic structure of the California spiny (or red rock) lobster (Panulirus interruptus) and find slight, but significant global population differentiation in mtDNA (ΦST = 0.006, P = 0.001; D(est_Chao) = 0.025) and seven nuclear microsatellites (F(ST) = 0.004, P < 0.001; D(est_Chao) = 0.03), despite the species' 240- to 330-day pelagic larval duration. Significant population structure does not correlate with distance between sampling locations, and pairwise FST between adjacent sites often exceeds that among geographically distant locations. This result would typically be interpreted as unexplainable, chaotic genetic patchiness. However, kinship levels differ significantly among sites (pseudo-F(16,988) = 1.39, P = 0.001), and ten of 17 sample sites have significantly greater numbers of kin than expected by chance (P < 0.05). Moreover, a higher proportion of kin within sites strongly correlates with greater genetic differentiation among sites (D(est_Chao), R(2) = 0.66, P < 0.005). Sites with elevated mean kinship were geographically proximate to regions of high upwelling intensity (R(2) = 0.41, P = 0.0009). These results indicate that P. interruptus does not maintain a single homogenous population, despite extreme dispersal potential. Instead, these lobsters appear to either have substantial localized recruitment or maintain planktonic larval cohesiveness whereby siblings more likely settle together than disperse across sites. More broadly, our results contribute to a growing number of studies showing that low F(ST) and high family structure across populations can coexist, illuminating the foundations of cryptic genetic patterns and the nature of marine dispersal.

Hierarchical population genetic structure in a direct developing antarctic marine invertebrate

Understanding the relationship between life-history variation and population structure in marine invertebrates is not straightforward. This is particularly true of polar species due to the difficulty of obtaining samples and a paucity of genomic resources from which to develop nuclear genetic markers. Such knowledge, however, is essential for understanding how different taxa may respond to climate change in the most rapidly warming regions of the planet. We therefore used over two hundred polymorphic Amplified Fragment Length Polymorphisms (AFLPs) to explore population connectivity at three hierachical spatial scales in the direct developing Antarctic topshell Margarella antarctica. To previously published data from five populations spanning a 1500 km transect along the length of the Western Antarctic Peninsula, we added new AFLP data for four populations separated by up to 6 km within Ryder Bay, Adelaide Island. Overall, we found a nonlinear isolation-by-distance pattern, suggestive of weaker population structure within Ryder Bay than is present over larger spatial scales. Nevertheless, significantly positive F_st values were obtained in all but two of ten pairwise population comparisons within the bay following Bonferroni correction for multiple tests. This is in contrast to a previous study of the broadcast spawner Nacella concinna that found no significant genetic differences among several of the same sites. By implication, the topshell's direct-developing lifestyle may constrain its ability to disperse even over relatively small geographic scales.

Population structure in the native range predicts the spread of introduced marine species

Forecasting invasion success remains a fundamental challenge in invasion biology. The effort to identify universal characteristics that predict which species become invasive has faltered in part because of the diversity of taxa and systems considered. Here, we use an alternative approach focused on the spread stage of invasions. FST, a measure of alternative fixation of alleles, is a common proxy for realized dispersal among natural populations, summarizing the combined influences of life history, behaviour, habitat requirements, population size, history and ecology. We test the hypothesis that population structure in the native range (FST) is negatively correlated with the geographical extent of spread of marine species in an introduced range. An analysis of the available data (29 species, nine phyla) revealed a significant negative correlation (R(2) = 0.245-0.464) between FST and the extent of spread of non-native species. Mode FST among pairwise comparisons between populations in the native range demonstrated the highest predictive power (R(2) = 0.464, p < 0.001). There was significant improvement when marker type was considered, with mtDNA datasets providing the strongest relationship (n = 21, R(2) = 0.333-0.516). This study shows that FST can be used to make qualitative predictions concerning the geographical extent to which a non-native marine species will spread once established in a new area.

Going where traditional markers have not gone before: utility of and promise for RAD sequencing in marine invertebrate phylogeography and population genomics

Characterization of large numbers of single-nucleotide polymorphisms (SNPs) throughout a genome has the power to refine the understanding of population demographic history and to identify genomic regions under selection in natural populations. To this end, population genomic approaches that harness the power of next-generation sequencing to understand the ecology and evolution of marine invertebrates represent a boon to test long-standing questions in marine biology and conservation. We employed restriction-site-associated DNA sequencing (RAD-seq) to identify SNPs in natural populations of the sea anemone Nematostella vectensis, an emerging cnidarian model with a broad geographic range in estuarine habitats in North and South America, and portions of England. We identified hundreds of SNP-containing tags in thousands of RAD loci from 30 barcoded individuals inhabiting four locations from Nova Scotia to South Carolina. Population genomic analyses using high-confidence SNPs resulted in a highly-resolved phylogeography, a result not achieved in previous studies using traditional markers. Plots of locus-specific FST against heterozygosity suggest that a majority of polymorphic sites are neutral, with a smaller proportion suggesting evidence for balancing selection. Loci inferred to be under balancing selection were mapped to the genome, where 90% were located in gene bodies, indicating potential targets of selection. The results from analyses with and without a reference genome supported similar conclusions, further highlighting RAD-seq as a method that can be efficiently applied to species lacking existing genomic resources. We discuss the utility of RAD-seq approaches in burgeoning Nematostella research as well as in other cnidarian species, particularly corals and jellyfishes, to determine phylogeographic relationships of populations and identify regions of the genome undergoing selection.

Consequences of a poecilogonous life history for genetic structure in coastal populations of the polychaete Streblospio benedicti

In many species, alternative developmental pathways lead to the production of two distinct phenotypes, promoting the evolution of morphological novelty and diversification. Offspring type in marine invertebrates influences transport time by ocean currents, which dictate dispersal potential and gene flow, and thus has sweeping evolutionary effects on the potential for local adaptation and on rates of speciation, extinction and molecular evolution. Here, we use the polychaete Streblospio benedicti to investigate the effects of dimorphic offspring type on gene flow and genetic structure in coastal populations. We use 84 single nucleotide polymorphism (SNP) markers for this species to assay populations on the East and West Coasts of the United States. Using these markers, we found that in their native East Coast distribution, populations of S. benedicti have high-population genetic structure, but this structure is associated primarily with geographic separation rather than developmental differences. Interestingly, very little genetic differentiation is recovered between individuals of different development types when they occur in the same or nearby populations, further supporting that this is a true case of poecilogony. In addition, we were able to demonstrate that the recently introduced (~100 ya) West Coast populations probably originated from a lecithotrophic population near Delaware.

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.