PATTERN BENEATH THE CHAOS: THE EFFECT OF RECRUITMENT ON GENETIC PATCHINESS IN AN INTERTIDAL LIMPET

Species' attributes predict the relative magnitude of ecological and genetic recovery following mass mortality

Theoretically, species' characteristics should allow estimation of dispersal potential and, in turn, explain levels of population genetic differentiation. However, a mismatch between traits and genetic patterns is often reported for marine species, and interpreted as evidence that life-history traits do not influence dispersal. Here, we couple ecological and genomic methods to test the hypothesis that species with attributes favouring greater dispersal potential-e.g., longer pelagic duration, higher fecundity and larger population size-have greater realized dispersal overall. We used a natural experiment created by a large-scale and multispecies mortality event which created a "clean slate" on which to study recruitment dynamics, thus simplifying a usually complex problem. We surveyed four species of differing dispersal potential to quantify the abundance and distribution of recruits and to genetically assign these recruits to probable parental sources. Species with higher dispersal potential recolonized a broader extent of the impacted range, did so more quickly and recovered more genetic diversity than species with lower dispersal potential. Moreover, populations of taxa with higher dispersal potential exhibited more immigration (71%-92% of recruits) than taxa with lower dispersal potential (17%-44% of recruits). By linking ecological with genomic perspectives, we demonstrate that a suite of interacting life-history and demographic attributes do influence species' realized dispersal and genetic neighbourhoods. To better understand species' resilience and recovery in this time of global change, integrative eco-evolutionary approaches are needed to more rigorously evaluate the effect of dispersal-linked attributes on realized dispersal and population genetic differentiation.

High gene flow in the silverlip pearl oyster Pinctada maxima between inshore and offshore sites near Eighty Mile Beach in Western Australia

An understanding of stock recruitment dynamics in fisheries is fundamental to successful management. Pinctada maxima is a bivalve mollusc widely distributed throughout the Indo-Pacific and is the main species targeted for cultured pearl and pearl shell production in Australia. Pearl production in Australia relies heavily on wild-caught individuals, the majority of which come from the Eighty Mile Beach region near Broome in Western Australia. In this study, we used a genotyping by sequencing approach to explore fine-scale patterns of genetic connectivity among inshore shallow and offshore deep populations of P. maxima near Eighty Mile Beach. Our results revealed high-levels of gene flow among inshore and offshore sites and no differences in genetic diversity between depths. Global estimates of genetic differentiation were low (F ST = 0.006) but significantly different from zero, and pairwise estimates of genetic differentiation among sites were significant in only 3% of comparisons. Moreover, Bayesian clustering detected no separation of inshore and offshore sample sites, and instead showed all samples to be admixed among sites, locations and depths. Despite an absence of any clear spatial clustering among sites, we identified a significant pattern of isolation by distance. In a dynamic environment like Eighty Mile Beach, genetic structure can change from year-to-year and successive dispersal and recruitment events over generations likely act to homogenize the population. Although we cannot rule out the null hypothesis of panmixia, our data indicate high levels of dispersal and connectivity among inshore and offshore fishing grounds.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Sweepstake reproductive success and collective dispersal produce chaotic genetic patchiness in a broadcast spawner

A long-standing paradox of marine populations is chaotic genetic patchiness (CGP), temporally unstable patterns of genetic differentiation that occur below the geographic scale of effective dispersal. Several mechanisms are hypothesized to explain CGP including natural selection, spatiotemporal fluctuations in larval source populations, self-recruitment, and sweepstake reproduction. Discriminating among them is extremely difficult but is fundamental to understanding how marine organisms reproduce and disperse. Here, we report a notable example of CGP in the Antarctic limpet, an unusually tractable system where multiple confounding explanations can be discounted. Using population genomics, temporally replicated sampling, surface drifters, and forward genetic simulations, we show that CGP likely arises from an extreme sweepstake event together with collective larval dispersal, while selection appears to be unimportant. Our results illustrate the importance of neutral demographic forces in natural populations and have important implications for understanding the recruitment dynamics, population connectivity, local adaptation, and resilience of marine populations.

Investigation of mechanisms underlying chaotic genetic patchiness in the intertidal marbled crab Pachygrapsus marmoratus (Brachyura: Grapsidae) across the Ligurian Sea

Background

Studies on marine community dynamics and population structures are limited by the lack of exhaustive knowledge on the larval dispersal component of connectivity. Genetic data represents a powerful tool in understanding such processes in the marine realm. When dealing with dispersion and connectivity in marine ecosystems, many evidences show patterns of genetic structure that cannot be explained by any clear geographic trend and may show temporal instability. This scenario is usually referred to as chaotic genetic patchiness, whose driving mechanisms are recognized to be selection, temporal shifts in local population dynamics, sweepstakes reproductive success and collective dispersal.

In this study we focused on the marbled crab Pachygrapsus marmoratus that inhabits the rocky shores of the Mediterranean Sea, Black Sea and East Atlantic Ocean, and disperses through planktonic larvae for about 1 month. P. marmoratus exhibits unexpectedly low connectivity levels at local scale, although well-defined phylogeographic patterns across the species’ distribution range were described. This has been explained as an effect of subtle geographic barriers or due to sweepstake reproductive success. In order to verify a chaotic genetic patchiness scenario, and to explore mechanisms underlying it, we planned our investigation within the Ligurian Sea, an isolated basin of the western Mediterranean Sea, and we genotyped 321 individuals at 11 microsatellite loci.

Results

We recorded genetic heterogeneity among our Ligurian Sea samples with the occurrence of genetic clusters not matching the original populations and a slight inter-population divergence, with the geographically most distant populations being the genetically most similar ones. Moreover, individuals from each site were assigned to all the genetic clusters. We also recorded evidences of self-recruitment and a higher than expected within-site kinship.

Conclusions

Overall, our results suggest that the chaotic genetic patchiness we found in P. marmoratus Ligurian Sea populations is the result of a combination of differences in reproductive success, en masse larval dispersion and local larval retention. This study defines P. marmoratus as an example of marine spawner whose genetic pool is not homogenous at population level, but rather split in a chaotic mosaic of slightly differentiated genetic patches derived from complex and dynamic ecological processes.

Genetic structure is stronger across human-impacted habitats than among islands in the coral Porites lobata

We examined genetic structure in the lobe coral Porites lobata among pairs of highly variable and high-stress nearshore sites and adjacent less variable and less impacted offshore sites on the islands of Oahu and Maui, Hawaii. Using an analysis of molecular variance framework, we tested whether populations were more structured by geographic distance or environmental extremes. The genetic patterns we observed followed isolation by environment, where nearshore and adjacent offshore populations showed significant genetic structure at both locations (AMOVA F ST = 0.04∼0.19, P < 0.001), but no significant isolation by distance between islands. Strikingly, corals from the two nearshore sites with higher levels of environmental stressors on different islands over 100 km apart with similar environmentally stressful conditions were genetically closer (FST = 0.0, P = 0.73) than those within a single location less than 2 km apart (FST = 0.04∼0.08, P < 0.01). In contrast, a third site with a less impacted nearshore site (i.e., less pronounced environmental gradient) showed no significant structure from the offshore comparison. Our results show much stronger support for environment than distance separating these populations. Our finding suggests that ecological boundaries from human impacts may play a role in forming genetic structure in the coastal environment, and that genetic divergence in the absence of geographical barriers to gene flow might be explained by selective pressure across contrasting habitats.

Sweepstakes reproductive success is absent in a New Zealand snapper (Chrysophrus auratus) population protected from fishing despite "tiny" Ne /N ratios elsewhere

A landmark study published in 2002 estimated a very small N_e /N ratio (around 10^-5 ) in a population of pink snapper (Chrysophrys auratus, Forster, 1801) in the Hauraki Gulf in New Zealand. It epitomized the tiny N_e /N ratios (<10^-3 ) reported in marine species due to the hypothesized operation of sweepstakes reproductive success (SRS). Here we re-evaluate the occurrence of SRS in marine species and the potential effect of fishing on the N_e /N ratio by studying the same species in the same region, but in a population that has been protected from fishing since 1975. We combine empirical, simulation and model-based approaches to estimate N_e (and N_b ) from genotypes of 1,044 adult fish and estimate N using recapture-probabilities. The estimated N_e /N ratio was much larger (0.33, SE: 0.14) than expected. The magnitude of estimates of population-wide variance in individual lifetime reproductive success (10-18) suggested that the sweepstakes effect was negligible in the study population. After evaluating factors that could explain the contrast between studies - experimental design, life history differences, environmental effects and the influence of exploitation on the N_e /N ratio - we conclude that the low N_e of the Hauraki Gulf population is associated with demographic instability in the harvested compared to the protected population despite circumstantial evidence that the 2002 study may have underestimated N_e . This study has broad implications for the prevailing view that reproductive success in the sea is largely driven by chance, and for genetic monitoring of populations using the N_e /N ratio and N_b .

Genetic homogeneity of the Pacific thread herring (Opisthonema libertate) (Günther, 1867) in the Eastern Pacific, inferred from mtDNA sequences

In the present study, the population genetic structure of the Pacific thread herring (Opisthonema libertate) was analyzed through mitochondrial DNA (mtDNA) control region sequences. Organisms were collected from June 2015 to July 2015 from four commercial landing sites (Bahia Magdalena, Guaymas and Mazatlan, Mexico, and Puntarenas, Costa Rica) and one artisanal sampling site (Puerto La Libertad, El Salvador). A total of 125 sequences were analyzed. High levels of haplotype (h = 0.990) and nucleotide (π = 0.030) diversity were found. Pairwise Φst comparisons indicated differences attributed mainly to the organisms from El Salvador. However, Bayesian inferences did not support the existence of different populations. The haplotype distribution between locations did not show a clear phylogeographic pattern. Mismatched distribution showed a unimodal pattern for the five sampled areas, indicative of sudden demographic expansion. These results were supported by Bayesian skyline plot. Our results do not support the hypothesis that the Pacific thread herring presents a population genetic structure. Future genetic comparisons should include a larger number of samples as well as more polymorphic molecular markers to further support our results.

Sand Dollar Larvae Show Within-Population Variation in Their Settlement Induction by Turbulence

Settlement-the generally irreversible transition from a planktonic phase to a benthic phase-is a critical stage in the life history of many shoreline organisms. It is reasonable to expect that larvae are under intense selection pressure to identify appropriate settlement habitat. Several decades of studies have focused mainly on local indicators that larvae use to identify suitable habitat, such as olfactory cues that indicate the presence of conspecifics or a favored food source. Our recent work has shown that the larvae of seashore-dwelling echinoids (sea urchins, sand dollars, and kin) can be primed to settle following a brief exposure to a broader-scale indicator of their approach to shore: an increase in fluid turbulence. Here we demonstrate that this priming shows within-population variation: the offspring of certain Pacific sand dollar (Dendraster excentricus) parents-both specific fathers and specific mothers, regardless of the other parent-are more responsive to turbulence than others. In particular, the observation of the effect correlating, in some cases, with specific fathers leads us to conclude that these behavioral differences are likely genetic and thus heritable. We also report that turbulence exposure causes larvae to temporarily sink to the bottom of a container of seawater and that larvae that respond in this way are also more likely to subsequently settle. We hypothesize a two-step scenario for the evolution of turbulence responsiveness at settlement and suggest that the evolutionary origin of these behaviors could be a driving force for population differentiation and speciation.

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.

INTERTIDAL MICROHABITAT AND SELECTION AT MPI: INTERLOCUS CONTRASTS IN THE NORTHERN ACORN BARNACLE, SEMIBALANUS BALANOIDES

Barnacles were sampled from various microhabitats in the rocky intertidal at multiple sites in two years. At sites in which there were large differences among microhabitats in temperature profiles, Mpi genotype frequencies were consistently and significantly different. Genotype frequencies for another allozyme locus (Gpi) as well as a DNA marker shown to be neutral (the mtDNA control region) were statistically homogeneous among thermal microhabitats at all sites in both years. The data indicate that temperature and/or desiccation mediated selection is operating at Mpi or a linked locus and that Mpi genotypes experience differential mortality in the various habitat types. If the relative fitness of genotypes is dependent on habitat type, the Mpi polymorphism may be actively maintained by a Levene model of balancing selection (Levene 1953). Because barnacle larvae are produced in abundance each year and spend five to eight weeks dispersing in the water column, there is little opportunity for the accumulation of adaptive divergence over the environmental grain size relevant in intertidal habitats. The Mpi polymorphism may be an important component of a suite of traits involved in the adaptation of barnacles to heterogeneous environments. Due to the relatively high concentration of mannose in a variety of algal groups, the metabolism of mannose may substantially affect individual performance and fitness in this and other species that feed on algae and phytoplankton. Because the Mpi locus is one of the most strongly polymorphic in marine organisms, these findings may be relevant for a diversity of other such species.

DEPENDENCE OF GENE FLOW ON GEOGRAPHIC DISTANCE IN TWO SOLITARY CORALS WITH DIFFERENT LARVAL DISPERSAL CAPABILITIES

When the level of gene flow among populations depends upon the geographic distance separating them, genetic differentiation is relatively enhanced. Although the larval dispersal capabilities of marine organisms generally correlate with inferred levels of average gene flow, the effect of different modes of larval development on the association between gene flow and geographic distance remains unknown. In this paper, I examined the relationship between gene flow and distance in two co-occurring solitary corals. Balanophyllia elegans broods large, nonfeeding planulae that generally crawl only short distances from their place of birth before settling. In contrast, Paracyathus stearnsii free-spawns and produces small planktonic larvae presumably capable of broad dispersal by oceanic currents. I calculated F-statistics using genetic variation at six (P. stearnsii) or seven (B. elegans) polymorphic allozyme loci revealed by starch gel electrophoresis, and used these F-statistics to infer levels of gene flow. Average levels of gene flow among twelve Californian localities agreed with previous studies: the species with planktonic, feeding larvae was less genetically subdivided than the brooding species. In addition, geographic isolation between populations appeared to affect gene flow between populations in very different ways in the two species. In the brooding B. elegans, gene flow declined with increasing separation, and distance explained 31% of the variation in gene flow. In the planktonically dispersed P. stearnsii distance of separation between populations at the scale studied (10-1000 km) explained only 1% of the variation in gene flow between populations. The mechanisms generating geographic genetic differentiation in species with different modes of larval development should vary fundamentally as a result of these qualitative differences in the dependence of gene flow on distance.

Current hypotheses to explain genetic chaos under the sea

Chaotic genetic patchiness (CGP) refers to surprising patterns of spatial and temporal genetic structure observed in some marine species at a scale where genetic variation should be efficiently homogenized by gene flow via larval dispersal. Here we review and discuss 4 mechanisms that could generate such unexpected patterns: selection, sweepstakes reproductive success, collective dispersal, and temporal shifts in local population dynamics. First, we review examples where genetic differentiation at specific loci was driven by diversifying selection, which was historically the first process invoked to explain CGP. Second, we turn to neutral demographic processes that may drive genome-wide effects, and whose effects on CGP may be enhanced when they act together. We discuss how sweepstakes reproductive success accelerates genetic drift and can thus generate genetic structure, provided that gene flow is not too strong. Collective dispersal is another mechanism whereby genetic structure can be maintained regardless of dispersal intensity, because it may prevent larval cohorts from becoming entirely mixed. Theoretical analyses of both the sweepstakes and the collective dispersal ideas are presented. Finally, we discuss an idea that has received less attention than the other ones just mentioned, namely temporal shifts in local population dynamics.

Patterns of Fish Connectivity between a Marine Protected Area and Surrounding Fished Areas

Patterns of connectivity and self-recruitment are recognized as key factors shaping the dynamics of marine populations. Connectivity is also essential for maintaining and restoring natural ecological processes with genetic diversity contributing to the adaptation and persistence of any species in the face of global disturbances. Estimates of connectivity are crucial to inform the design of both marine protected areas (MPAs) and MPA networks. Among several approaches, genetic structure is frequently used as a proxy for patterns of connectivity. Using 8 microsatellite loci, we investigated genetic structure of the two-banded sea bream Diplodus vulgaris, a coastal fish that is both commercially and ecologically important. Adults were sampled in 7 locations (stretches of coastline approximately 8 km long) and juveniles in 14 sites (~100 to 200 m of coastline) along 200 km of the Apulian Adriatic coast (SW Adriatic Sea), within and outside an MPA (Torre Guaceto MPA, Italy). Our study found similar genetic diversity indices for both the MPA and the surrounding fished areas. An overall lack of genetic structure among samples suggests high gene flow (i.e. connectivity) across a scale of at least 200 km. However, some local genetic divergences found in two locations demonstrate some heterogeneity in processes renewing the population along the Apulian Adriatic coast. Furthermore, two sites appeared genetically divergent, reinforcing our observations within the genetic makeup of adults and confirming heterogeneity in early stage genetics that can come from either different supply populations or from chaotic genetic patchiness occurring under temporal variation in recruitment and in the reproductive success. While the specific role of the MPA is not entirely known in this case, these results confirm the presence of regional processes and the key role of connectivity in maintaining the local population supply.

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.

Highly Connected Populations and Temporal Stability in Allelic Frequencies of a Harvested Crab from the Southern Pacific Coast

For marine invertebrates with a benthic adult form and a planktonic larva phase, the connectivity among populations is mainly based on larval dispersal. While an extended larval phase will promote gene flow, other factors such as an intensive fishery and geographical barriers could lead to changes in genetic variability. In this study, the population genetic structure of the commercial crab Metacarcinus edwardsii was analyzed along 700 km of the Chilean coast. The analysis, based on eight microsatellite loci genotyped from megalopae and adult crabs, considered temporal and spatial patterns of genetic variation. The results showed no evidence of spatial patterns in genetic structure, suggesting high connectivity among the sampling sites. The temporal analysis showed no evidence of changes in allele frequencies and no evidence of a recent bottleneck. The lack of spatial structure and allele variation over time could be explained by the interaction of factors such as i) low reproductive variance due to the capability of females to store sperm in the seminal receptacle, which can be used for successive broods, ii) high larval dispersal and iii) high individual reproductive output. Using our data as priors, a genetic modelling approach coincided, predicting this temporal and spatial stability. The same analysis showed that a reduction in population size leads to the loss of genetic variability in populations, as well as of the genetic cohesiveness between populations, pointing out the importance management for species under exploitation, such as M. edwardsii.

Moderate genetic drift is driven by extreme recruitment events in the invasive mollusk Crepidula fornicata

Effective population size (Ne) is a measure of genetic drift and is thus a central parameter in evolution, conservation genetics and invasion biology. Interestingly, in native marine species, Ne is typically several orders of magnitude lower than the census size. This pattern has often been explained by high fecundity, variation in reproductive success and pronounced early mortality, resulting in genetic drift across generations. Data documenting genetic drift and/or Ne in marine invasive species are, however, still scarce. We examined the importance of genetic drift in the invasive species Crepidula fornicata by genotyping 681 juveniles sampled during each annual recruitment peak over nine consecutive years in the Bay of Morlaix (Brittany, France). Observed variations in genetic diversity were partially explained by variation in recruitment intensity. In addition, substantial temporal genetic differentiation was documented (that is, genetic drift), and was attributed to nonrandom variance in the reproductive success of different breeding groups across years in the study species. Using a set of single-sample and temporal estimators for Ne, we estimated Ne to be three or four orders of magnitude smaller than the census size (Nc). On one hand, this reduction in Ne relative to Nc appeared congruent with, although slight higher than, values commonly observed in native marine species. Particular life-history traits of this invasive species may play an important role in buffering genetic drift. On the other hand, Ne still remained far below Nc, hence, possibly reducing the efficiency of selection effects.

Understanding the Spatial Scale of Genetic Connectivity at Sea: Unique Insights from a Land Fish and a Meta-Analysis

Quantifying the spatial scale of population connectivity is important for understanding the evolutionary potential of ecologically divergent populations and for designing conservation strategies to preserve those populations. For marine organisms like fish, the spatial scale of connectivity is generally set by a pelagic larval phase. This has complicated past estimates of connectivity because detailed information on larval movements are difficult to obtain. Genetic approaches provide a tractable alternative and have the added benefit of estimating directly the reproductive isolation of populations. In this study, we leveraged empirical estimates of genetic differentiation among populations with simulations and a meta-analysis to provide a general estimate of the spatial scale of genetic connectivity in marine environments. We used neutral genetic markers to first quantify the genetic differentiation of ecologically-isolated adult populations of a land dwelling fish, the Pacific leaping blenny (Alticus arnoldorum), where marine larval dispersal is the only probable means of connectivity among populations. We then compared these estimates to simulations of a range of marine dispersal scenarios and to collated FST and distance data from the literature for marine fish across diverse spatial scales. We found genetic connectivity at sea was extensive among marine populations and in the case of A. arnoldorum, apparently little affected by the presence of ecological barriers. We estimated that ~5000 km (with broad confidence intervals ranging from 810-11,692 km) was the spatial scale at which evolutionarily meaningful barriers to gene flow start to occur at sea, although substantially shorter distances are also possible for some taxa. In general, however, such a large estimate of connectivity has important implications for the evolutionary and conservation potential of many marine fish communities.

Kin-Aggregations Explain Chaotic Genetic Patchiness, a Commonly Observed Genetic Pattern, in a Marine Fish

The phenomenon of chaotic genetic patchiness is a pattern commonly seen in marine organisms, particularly those with demersal adults and pelagic larvae. This pattern is usually associated with sweepstakes recruitment and variable reproductive success. Here we investigate the biological underpinnings of this pattern in a species of marine goby Coryphopterus personatus. We find that populations of this species show tell-tale signs of chaotic genetic patchiness including: small, but significant, differences in genetic structure over short distances; a non-equilibrium or "chaotic" pattern of differentiation among locations in space; and within locus, within population deviations from the expectations of Hardy-Weinberg equilibrium (HWE). We show that despite having a pelagic larval stage, and a wide distribution across Caribbean coral reefs, this species forms groups of highly related individuals at small spatial scales (<10 metres). These spatially clustered family groups cause the observed deviations from HWE and local population differentiation, a finding that is rarely demonstrated, but could be more common than previously thought.

Temporal Stability of Genetic Structure in a Mesopelagic Copepod

Although stochasticity in oceanographic conditions is known to be an important driver of temporal genetic change in many marine species, little is known about whether genetically distinct plankton populations can persist in open ocean habitats. A prior study demonstrated significant population genetic structure among oceanic gyres in the mesopelagic copepod Haloptilus longicornis in both the Atlantic and Pacific Oceans, and we hypothesized that populations within each gyre represent distinct gene pools that persist over time. We tested this expectation through basin-scale sampling across the Atlantic Ocean in 2010 and 2012. Using both mitochondrial (mtCOII) and microsatellite markers (7 loci), we show that the genetic composition of populations was stable across two years in both the northern and southern subtropical gyres. Genetic variation in this species was partitioned among ocean gyres (F_CT = 0.285, P < 0.0001 for mtCOII, F_CT = 0.013, P < 0.0001 for microsatellites), suggesting strong spatial population structure, but no significant partitioning was found among sampling years. This temporal persistence of population structure across a large geographic scale was coupled with chaotic genetic patchiness at smaller spatial scales, but the magnitude of genetic differentiation was an order of magnitude lower at these smaller scales. Our results demonstrate that genetically distinct plankton populations persist over time in highly-dispersive open ocean habitats, and this is the first study to rigorously test for temporal stability of large scale population structure in the plankton.

Does natural selection explain the fine scale genetic structure at the nuclear exon Glu-5' in blue mussels from Kerguelen?

The Kerguelen archipelago, isolated in the Southern Ocean, shelters a blue mussel Mytilus metapopulation far from any influence of continental populations or any known hybrid zone. The finely carved coast leads to a highly heterogeneous habitat. We investigated the impact of the environment on the genetic structure in those Kerguelen blue mussels by relating allele frequencies to habitat descriptors. A total sample comprising up to 2248 individuals from 35 locations was characterized using two nuclear markers, mac-1 and Glu-5', and a mitochondrial marker (COI). The frequency data from 9 allozyme loci in 9 of these locations were also reanalyzed. Two other nuclear markers (EFbis and EFprem's) were monomorphic. Compared to Northern Hemisphere populations, polymorphism in Kerguelen blue mussels was lower for all markers except for the exon Glu-5'. At Glu-5', genetic differences were observed between samples from distinct regions (F CT = 0.077), as well as within two regions, including between samples separated by <500 m. No significant differentiation was observed in the AMOVA analyses at the two other markers (mac-1 and COI). Like mac-1, all allozyme loci genotyped in a previous publication, displayed lower differentiation (Jost's D) and F ST values than Glu-5'. Power simulations and confidence intervals support that Glu-5' displays significantly higher differentiation than the other loci (except a single allozyme for which confidence intervals overlap). AMOVA analyses revealed significant effects of the giant kelp Macrocystis and wave exposure on this marker. We discuss the influence of hydrological conditions on the genetic differentiation among regions. In marine organisms with high fecundity and high dispersal potential, gene flow tends to erase differentiation, but this study showed significant differentiation at very small distance. This may be explained by the particular hydrology and the carved coastline of the Kerguelen archipelago, together with spatially variable selection at Glu-5'.

Depth as an organizing force in Pocillopora damicornis: intra-reef genetic architecture

Relative to terrestrial plants, and despite similarities in life history characteristics, the potential for corals to exhibit intra-reef local adaptation in the form of genetic differentiation along an environmental gradient has received little attention. The potential for natural selection to act on such small scales is likely increased by the ability of coral larval dispersal and settlement to be influenced by environmental cues. Here, we combine genetic, spatial, and environmental data for a single patch reef in Kāne‘ohe Bay, O‘ahu, Hawai‘i, USA in a landscape genetics framework to uncover environmental drivers of intra-reef genetic structuring. The genetic dataset consists of near-exhaustive sampling (n = 2352) of the coral, Pocillopora damicornis at our study site and six microsatellite genotypes. In addition, three environmental parameters – depth and two depth-independent temperature indices – were collected on a 4 m grid across 85 locations throughout the reef. We use ordinary kriging to spatially interpolate our environmental data and estimate the three environmental parameters for each colony. Partial Mantel tests indicate a significant correlation between genetic relatedness and depth while controlling for space. These results are also supported by multi-model inference. Furthermore, spatial Principle Component Analysis indicates a statistically significant genetic cline along a depth gradient. Binning the genetic dataset based on size-class revealed that the correlation between genetic relatedness and depth was significant for new recruits and increased for larger size classes, suggesting a possible role of larval habitat selection as well as selective mortality in structuring intra-reef genetic diversity. That both pre- and post-recruitment processes may be involved points to the adaptive role of larval habitat selection in increasing adult survival. The conservation importance of uncovering intra-reef patterns of genetic diversity is discussed.

First evidence of inbreeding, relatedness and chaotic genetic patchiness in the holoplanktonic jellyfish Pelagia noctiluca (Scyphozoa, Cnidaria)

Genetic drift and non-random mating seldom influence species with large breeding populations and high dispersal potential, characterized by unstructured gene pool and panmixia at a scale lower than the minimum dispersal range of individuals. In the present study, a set of nine microsatellite markers was developed and used to investigate the spatio-temporal genetic patterns of the holoplanktonic jellyfish Pelagia noctiluca (Scyphozoa) in the Southern Tyrrhenian Sea. Homozygote excess was detected at eight loci, and individuals exhibited intra-population relatedness higher than expected by chance in at least three samples. This result was supported by the presence of siblings in at least 5 out 8 samples, 4 of which contained full-sib in addition to half-sib dyads. Having tested and ruled out alternative explanations as null alleles, our results suggest the influence of reproductive and behavioural features in shaping the genetic structure of P. noctiluca, as outcomes of population genetics analyses pointed out. Indeed, the genetic differentiation among populations was globally small but highlighted: a) a spatial genetic patchiness uncorrelated with distance between sampling locations, and b) a significant genetic heterogeneity between samples collected in the same locations in different years. Therefore, despite its extreme dispersal potential, P. noctiluca does not maintain a single homogenous population, but rather these jellyfish appear to have intra-bloom localized recruitment and/or individual cohesiveness, whereby siblings more likely swarm together as a single group and remain close after spawning events. These findings provide the first evidence of family structures and consequent genetic patchiness in a species with highly dispersive potential throughout its whole life cycle, contributing to understanding the patterns of dispersal and connectivity in marine environments.

Temporal genetic structure in a poecilogonous polychaete: the interplay of developmental mode and environmental stochasticity

Background

Temporal variation in the genetic structure of populations can be caused by multiple factors, including natural selection, stochastic environmental variation, migration, or genetic drift. In benthic marine species, the developmental mode of larvae may indicate a possibility for temporal genetic variation: species with dispersive planktonic larvae are expected to be more likely to show temporal genetic variation than species with benthic or brooded non-dispersive larvae, due to differences in larval mortality and dispersal ability. We examined temporal genetic structure in populations of Pygospio elegans, a poecilogonous polychaete with within-species variation in developmental mode. P. elegans produces either planktonic, benthic, or intermediate larvae, varying both among and within populations, providing a within-species test of the generality of a relationship between temporal genetic variation and larval developmental mode.

Results

In contrast to our expectations, our microsatellite analyses of P. elegans revealed temporal genetic stability in the UK population with planktonic larvae, whereas there was variation indicative of drift in temporal samples of the populations from the Baltic Sea, which have predominantly benthic and intermediate larvae. We also detected temporal variation in relatedness within these populations. A large temporal shift in genetic structure was detected in a population from the Netherlands, having multiple developmental modes. This shift could have been caused by local extiction due to extreme environmental conditions and (re)colonization by planktonic larvae from neighboring populations.

Conclusions

In our study of P. elegans, temporal genetic variation appears to be due to not only larval developmental mode, but also the stochastic environment of adults. Large temporal genetic shifts may be more likely in marine intertidal habitats (e.g. North Sea and Wadden Sea) which are more prone to environmental stochasticity than the sub-tidal Baltic habitats. Sub-tidal and/or brackish (less saline) habitats may support smaller P. elegans populations and these may be more susceptible to the effects of random genetic drift. Moreover, higher frequencies of asexual reproduction and the benthic larval developmental mode in these populations leads to higher relatedness and contributes to drift. Our results indicate that a general relationship between larval developmental mode and temporal genetic variation may not exist.

Restriction to large-scale gene flow vs. regional panmixia among cold seep Escarpia spp. (Polychaeta, Siboglinidae)

The history of colonization and dispersal in fauna distributed among deep-sea chemosynthetic ecosystems remains enigmatic and poorly understood because of an inability to mark and track individuals. A combination of molecular, morphological and environmental data improves understanding of spatial and temporal scales at which panmixia, disruption of gene flow or even speciation may occur. Vestimentiferan tubeworms of the genus Escarpia are important components of deep -sea cold seep ecosystems, as they provide long-term habitat for many other taxa. Three species of Escarpia, Escarpia spicata [Gulf of California (GoC)], Escarpia laminata [Gulf of Mexico (GoM)] and Escarpia southwardae (West African Cold Seeps), have been described based on morphology, but are not discriminated through the use of mitochondrial markers (cytochrome oxidase subunit 1; large ribosomal subunit rDNA, 16S; cytochrome b). Here, we also sequenced the exon-primed intron-crossing Haemoglobin subunit B2 intron and genotyped 28 microsatellites to (i) determine the level of genetic differentiation, if any, among the three geographically separated entities and (ii) identify possible population structure at the regional scale within the GoM and West Africa. Results at the global scale support the occurrence of three genetically distinct groups. At the regional scale among eight sampling sites of E. laminata (n = 129) and among three sampling sites of E. southwardae (n = 80), no population structure was detected. These findings suggest that despite the patchiness and isolation of seep habitats, connectivity is high on regional scales.

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.

AFLPs reveal different population genetic structure under contrasting environments in the marine snail Nucella lapillus L

Dispersal has received growing attention in marine ecology, particularly since evidence obtained with up-to-date techniques challenged the traditional view. The dogwhelk Nucella lapillus L., a sedentary gastropod with direct development, is a good example: dispersal was traditionally assumed to be limited until studies with microsatellites disputed this idea. To shed some light on this controversy, the genetic structure of dogwhelk populations in northwest Spain was investigated with highly polymorphic AFLP markers giving special attention to the influence of hydrodynamic stress. In agreement with the expectations for a poor disperser, our results show a significant genetic structure at regional (<200 km) and areal scales (<15 km). However, the spatial genetic structure varied with wave-exposure in the present case study: IBD was evident under sheltered conditions but absent from the exposed area where genetic differentiation was stronger. Our results provide evidence that differences in wave-exposure can exert a detectable influence on the genetic structure of coastal organisms, even in species without a planktonic larva.

Variance in the reproductive success of flat oyster Ostrea edulis L. assessed by parentage analyses in natural and experimental conditions

In order to document further the phenomena of variance in reproductive success in natural populations of the European flat oyster Ostrea edulis, two complementary studies based on natural and experimental populations were conducted. The first part of this work was focused on paternity analyses using a set of four microsatellite markers for larvae collected from 13 brooding females sampled in Quiberon Bay (Brittany, France). The number of individuals contributing as the male parent to each progeny assay was highly variable, ranging from 2 to more than 40. Moreover, paternal contributions showed a much skewed distribution, with some males contributing to 50-100% of the progeny assay. The second part of this work consisted of the analysis of six successive cohorts experimentally produced from an acclimated broodstock (62 wild oysters sampled in the Quiberon Bay). Allelic richness was significantly higher in the adult population than in the temporal cohorts collected. Genetic differentiation (F(st) estimates) was computed for each pair of samples and all significant values ranged from 0.7 to 11.9%. A limited effective number of breeders (generally below 25) was estimated in the six temporal cohorts. The study gives first indications of the high variance in reproductive success as well as a reduced effective size, not only under experimental conditions but also in the wild. Surprisingly, the pool of the successive cohorts, based on the low number of loci used, appeared to depict a random and representative set of alleles of the progenitor population, indicating that the detection of patterns of temporal genetic differentiation at a local scale most likely depends on the sampling window.

The concept of population in clonal organisms: mosaics of temporally colonized patches are forming highly diverse meadows of Zostera marina in Brittany

Seagrasses structure some of the world's key coastal ecosystems presently in decline due to human activities and global change. The ability to cope with environmental changes and the possibilities for shifts in distribution range depend largely on their evolvability and dispersal potential. As large-scale data usually show strong genetic structure for seagrasses, finer-grained work is needed to understand the local processes of dispersal, recruitment and colonization that could explain the apparent lack of exchange across large distances. We aimed to assess the fine-grained genetic structure of one of the most important and widely distributed seagrasses, Zostera marina, from seven meadows in Brittany, France. Both classic population genetics and network analysis confirmed a pattern of spatial segregation of polymorphism at both regional and local scales. One location exhibiting exclusively the variety 'angustifolia' did not appear more differentiated than the others, but instead showed a central position in the network analysis, confirming the status of this variety as an ecotype. This phenotypic diversity and the high allelic richness at nine microsatellites (2.33-9.67 alleles/locus) compared to levels previously reported across the distribution range, points to Brittany as a centre of diversity for Z. marina at both genetic and phenotypic levels. Despite dispersal potential of several 100 m, a significant pattern of genetic differentiation, even at fine-grained scale, revealed 'genetic patchiness'. Meadows seem to be composed of a mosaic of clones with distinct origins in space and time, a result that calls into question the accuracy of the concept of populations for such partially clonal species.

Genetic structure of a recent climate change-driven range extension

The life-history strategies of some species make them strong candidates for rapid exploitation of novel habitat under new climate regimes. Some early-responding species may be considered invasive, and negatively impact on 'naïve' ecosystems. The barrens-forming sea urchin Centrostephanus rodgersii is one such species, having a high dispersal capability and a high-latitude range margin limited only by a developmental temperature threshold. Within this species' range in eastern Australian waters, sea temperatures have increased at greater than double the global average rate. The coinciding poleward range extension of C. rodgersii has caused major ecological changes, threatening reef biodiversity and fisheries productivity. We investigated microsatellite diversity and population structure associated with range expansion by this species. Generalized linear model analyses revealed no reduction in genetic diversity in the newly colonized region. A 'seascape genetics' analysis of genetic distances found no spatial genetic structure associated with the range extension. The distinctive genetic characteristic of the extension zone populations was reduced population-specific F(ST), consistent with very rapid population expansion. Demographic and genetic simulations support our inference of high connectivity between pre- and post-extension zones. Thus, the range shift appears to be a poleward extension of the highly-connected rangewide population of C. rodgersii. This is consistent with advection of larvae by the intensified warm water East Australian current, which has also increased Tasmanian Sea temperatures above the species' lower developmental threshold. Thus, ocean circulation changes have improved the climatic suitability of novel habitat for C. rodgersii and provided the supply of recruits necessary for colonization.

Microsatellite markers reveal shallow genetic differentiation between cohorts of the common sea urchin Paracentrotus lividus (Lamarck) in northwest Mediterranean

Temporal variability was studied in the common sea urchin Paracentrotus lividus through the analysis of the genetic composition of three yearly cohorts sampled over two consecutive springs in a locality in northwestern Mediterranean. Individuals were aged using growth ring patterns observed in tests and samples were genotyped for five microsatellite loci. No reduction of genetic diversity was observed relative to a sample of the adult population from the same location or within cohorts across years. F(ST )and amova results indicated that the differentiation between cohorts is rather shallow and not significant, as most variability is found within cohorts and within individuals. This mild differentiation translated into estimates of effective population size of 90-100 individuals. When the observed excess of homozygotes was taken into account, the estimate of the average number of breeders increased to c. 300 individuals. Given our restricted sampling area and the known small-scale heterogeneity in recruitment in this species, our results suggest that at stretches of a few kilometres of shoreline, large numbers of progenitors are likely to contribute to the larval pool at each reproduction event. Intercohort variation in our samples is six times smaller than spatial variation between adults of four localities in the western Mediterranean. Our results indicate that, notwithstanding the stochastic events that take place during the long planktonic phase and during the settlement and recruitment processes, reproductive success in this species is high enough to produce cohorts genetically diverse and with little differentiation between them. Further research is needed before the link between genetic structure and underlying physical and biological processes can be well established.

Spatial and temporal population genetic structure of four northeastern Pacific littorinid gastropods: the effect of mode of larval development on variation at one mitochondrial and two nuclear DNA markers

We investigated the effect of development mode on the spatial and temporal population genetic structure of four littorinid gastropod species. Snails were collected from the same three sites on the west coast of Vancouver Island, Canada in 1997 and again in 2007. DNA sequences were obtained for one mitochondrial gene, cytochrome b (Cyt b), and for up to two nuclear genes, heat shock cognate 70 (HSC70) and aminopeptidase N intron (APN54). We found that the mean level of genetic diversity and long-term effective population sizes (N(e)) were significantly greater for two species, Littorina scutulata and L. plena, that had a planktotrophic larval stage than for two species, Littorina sitkana and L. subrotundata, that laid benthic egg masses which hatched directly into crawl-away juveniles. Predictably, two poorly dispersing species, L. sitkana and L. subrotundata, showed significant spatial genetic structure at an 11- to 65-km geographical scale that was not observed in the two planktotrophic species. Conversely, the two planktotrophic species had more temporal genetic structure over a 10-year interval than did the two direct-developing species and showed highly significant temporal structure for spatially pooled samples. The greater temporal genetic variation of the two planktotrophic species may have been caused by their high fecundity, high larval dispersal, and low but spatially correlated early survivorship. The sweepstakes-like reproductive success of the planktotrophic species could allow a few related females to populate hundreds of kilometres of coastline and may explain their substantially larger temporal genetic variance but lower spatial genetic variance relative to the direct-developing species.

Coalescence times and FST under a skewed offspring distribution among individuals in a population

Estimates of gene flow between subpopulations based on F(ST) (or N(ST)) are shown to be confounded by the reproduction parameters of a model of skewed offspring distribution. Genetic evidence of population subdivision can be observed even when gene flow is very high, if the offspring distribution is skewed. A skewed offspring distribution arises when individuals can have very many offspring with some probability. This leads to high probability of identity by descent within subpopulations and results in genetic heterogeneity between subpopulations even when Nm is very large. Thus, we consider a limiting model in which the rates of coalescence and migration can be much higher than for a Wright-Fisher population. We derive the densities of pairwise coalescence times and expressions for F(ST) and other statistics under both the finite island model and a many-demes limit model. The results can explain the observed genetic heterogeneity among subpopulations of certain marine organisms despite substantial gene flow.

Fine-scale spatial genetic structure and dispersal among spotted salamander (Ambystoma maculatum) breeding populations

We examined fine-scale genetic variation among breeding aggregations of the spotted salamander (Ambystoma maculatum) to quantify dispersal, interpopulation connectivity and population genetic structure. Spotted salamanders rely on temporary ponds or wetlands for aggregate breeding. Adequate breeding sites are relatively isolated from one another and field studies suggest considerable adult site fidelity; therefore, we expected to find population structure and differentiation at small spatial scales. We used microsatellites to estimate population structure and dispersal among 29 breeding aggregations in Tompkins County, New York, USA, an area encompassing 1272 km(2). Bayesian and frequency-based analyses revealed fine-scale genetic structure with two genetically defined demes: the North deme included seven breeding ponds, and the South deme included 13 ponds. Nine ponds showed evidence of admixture between these two genetic pools. Bayesian assignment tests for detection of interpopulation dispersal indicate that immigration among ponds is common within demes, and that certain populations serve as sources of immigrants to neighbouring ponds. Likewise, spatial genetic correlation analyses showed that populations < or = 4.8 km distant from each other show significant genetic correlation that is not evident at higher scales. Within-population levels of relatedness are consistently larger than expected if mating were completely random across ponds, and in the case of a few ponds, within-population processes such as inbreeding or reproductive skew contribute significantly to differentiation from neighbouring ponds. Our data underscore the importance of these within-population processes as a source of genetic diversity across the landscape, despite considerable population connectivity. Our data further suggest that spotted salamander breeding groups behave as metapopulations, with population clusters as functional units, but sufficient migration among demes to allow for potential rescue and recolonization. Amphibian habitats are becoming increasingly fragmented and a clear understanding of dispersal and patterns of population connectivity for taxa with different ecologies and life histories is crucial for their conservation.