Connectivity and resilience of coral reef metapopulations in marine protected areas: matching empirical efforts to predictive needs

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.

Potential Connectivity of Coldwater Black Coral Communities in the Northern Gulf of Mexico

The black coral Leiopathes glaberrima is a foundation species of deep-sea benthic communities but little is known of the longevity of its larvae and the timing of spawning because it inhabits environments deeper than 50 m that are logistically challenging to observe. Here, the potential connectivity of L. glaberrima in the northern Gulf of Mexico was investigated using a genetic and a physical dispersal model. The genetic analysis focused on data collected at four sites distributed to the east and west of Mississippi Canyon, provided information integrated over many (~10,000) generations and revealed low but detectable realized connectivity. The physical dispersal model simulated the circulation in the northern Gulf at a 1km horizontal resolution with transport-tracking capabilities; virtual larvae were deployed 12 times over the course of 3 years and followed over intervals of 40 days. Connectivity between sites to the east and west of the canyon was hampered by the complex bathymetry, by differences in mean circulation to the east and west of the Mississippi Canyon, and by flow instabilities at scales of a few kilometers. Further, the interannual variability of the flow field surpassed seasonal changes. Together, these results suggest that a) dispersal among sites is limited, b) any recovery in the event of a large perturbation will depend on local larvae produced by surviving individuals, and c) a competency period longer than a month is required for the simulated potential connectivity to match the connectivity from multi-locus genetic data under the hypothesis that connectivity has not changed significantly over the past 10,000 generations.

Spearfishing-induced behavioral changes of an unharvested species inside and outside a marine protected area

By prohibiting fishing, marine protected areas (MPAs) provide a refuge for harvested species. Humans are often perceived as predators by prey and therefore respond fearfully to humans. Thus, fish responses to humans inside and outside of an MPA can provide insights into their perception of humans as a predatory threat. Previous studies have found differences in the distance that harvested species of fish initiate flight (flight initiation distance—FID) from humans inside and outside an MPA, but less is known about unharvested species. We focused on whether the lined bristletooth Ctenochaetus striatus, an unharvested surgeonfish, can discriminate between a snorkeler and a snorkeler with a spear gun inside and outside of a no-take MPA in Mo’orea, French Polynesia. Additionally, we incorporated starting distance (the distance between the person and prey at the start of an experimental approach), a variable that has been found to be important in assessing prey escape decisions in terrestrial species, but that has not been extensively studied in aquatic systems. Lined bristletooth FID was significantly greater in the presence of a spear gun and varied depending on if the spear gun encounter was inside or outside of the MPA. These results imply a degree of sophistication of fish antipredator behavior, generate questions as to how a nontargeted species of fish could acquire fear of humans, and demonstrate that behavioral surveys can provide insights about antipredator behavior.

Reciprocal feedbacks between spatial subsidies and reserve networks in coral reef meta-ecosystems

Top-down processes such as predation and herbivory have been shown to control the dynamics of communities across a range of ecosystems by generating trophic cascades. However, theory is only beginning to describe how these local trophic processes interact with spatial subsidies in the form of material (nutrient, detritus) transport and organismal dispersal to (1) shape the structure of interconnected (meta-) ecosystems and (2) determine their optimal management via reserve networks. Here, we develop a meta-ecosystem model to understand how the reciprocal feedbacks between spatial subsidies and reserve networks modulate the importance of top-down control in a simple herbivorous fish-macroalgae-coral system. We show that in large and isolated reserve networks where connectivity between protected and unprotected areas is limited, spatial subsidies remain largely confined to reserves. This retention of spatial subsidies promotes the top-down control of corals and macroalgae by herbivores inside reserves but reduces it outside reserves. Conversely, in small and aggregated reserves where connectivity between protected and unprotected areas is high, the spillover of spatial subsidies causes a reduction in top-down control of corals and macroalgae by herbivores inside reserves and an increase in the strength of top-down control outside reserves. In addition, we demonstrate that there is a trade-off between local and regional conservation objectives when designing reserve networks: small and aggregated reserves based on the extent of dispersal maximize the abundance of corals and herbivores regionally, whereas large and isolated reserves always maximize the abundance of corals within reserves, regardless of the extent of dispersal. The existence of such "conservation traps," which arise from the fulfillment of population-level objectives within local reserves at the cost of community-level objectives at regional scales, suggests the importance of adopting a more holistic strategy to manage complex and interconnected ecosystems.

No Reef Is an Island: Integrating Coral Reef Connectivity Data into the Design of Regional-Scale Marine Protected Area Networks

We integrated coral reef connectivity data for the Caribbean and Gulf of Mexico into a conservation decision-making framework for designing a regional scale marine protected area (MPA) network that provides insight into ecological and political contexts. We used an ocean circulation model and regional coral reef data to simulate eight spawning events from 2008-2011, applying a maximum 30-day pelagic larval duration and 20% mortality rate. Coral larval dispersal patterns were analyzed between coral reefs across jurisdictional marine zones to identify spatial relationships between larval sources and destinations within countries and territories across the region. We applied our results in Marxan, a conservation planning software tool, to identify a regional coral reef MPA network design that meets conservation goals, minimizes underlying threats, and maintains coral reef connectivity. Our results suggest that approximately 77% of coral reefs identified as having a high regional connectivity value are not included in the existing MPA network. This research is unique because we quantify and report coral larval connectivity data by marine ecoregions and Exclusive Economic Zones (EZZ) and use this information to identify gaps in the current Caribbean-wide MPA network by integrating asymmetric connectivity information in Marxan to design a regional MPA network that includes important reef network connections. The identification of important reef connectivity metrics guides the selection of priority conservation areas and supports resilience at the whole system level into the future.

Patterns, causes, and consequences of marine larval dispersal

Quantifying the probability of larval exchange among marine populations is key to predicting local population dynamics and optimizing networks of marine protected areas. The pattern of connectivity among populations can be described by the measurement of a dispersal kernel. However, a statistically robust, empirical dispersal kernel has been lacking for any marine species. Here, we use genetic parentage analysis to quantify a dispersal kernel for the reef fish Elacatinus lori, demonstrating that dispersal declines exponentially with distance. The spatial scale of dispersal is an order of magnitude less than previous estimates-the median dispersal distance is just 1.7 km and no dispersal events exceed 16.4 km despite intensive sampling out to 30 km from source. Overlaid on this strong pattern is subtle spatial variation, but neither pelagic larval duration nor direction is associated with the probability of successful dispersal. Given the strong relationship between distance and dispersal, we show that distance-driven logistic models have strong power to predict dispersal probabilities. Moreover, connectivity matrices generated from these models are congruent with empirical estimates of spatial genetic structure, suggesting that the pattern of dispersal we uncovered reflects long-term patterns of gene flow. These results challenge assumptions regarding the spatial scale and presumed predictors of marine population connectivity. We conclude that if marine reserve networks aim to connect whole communities of fishes and conserve biodiversity broadly, then reserves that are close in space (<10 km) will accommodate those members of the community that are short-distance dispersers.

Connectivity of Marine Protected Areas and Its Relation with Total Kinetic Energy

The East Continental Shelf (ECS) of Brazil is a hotspot of endemism and biodiversity of reef biota in the South Atlantic, hosting a number of Marine Protected Areas (MPAs). Connectivity of MPAs through larval dispersal influences recruitment, population dynamics, genetic structure and biogeography in coral reef ecosystems. Connectivity of protected reef ecosystem in the ECS was investigated with a hydrodynamic model (ROMS) forcing an Individual Based Model (IBM—Ichthyop), and used groupers (genus Mycteroperca) as functional group. The hydrodynamic output from ROMS was compared with satellite data and showed good agreement with observed surface fields. Eggs were released, in IBM experiments, from April to September along six years (2002–2007) in five MPAs along the ECS. Intrannual variability in recruitment and self-recruitment of grouper larvae was observed, as well as a negative correlation of these population parameters with total Kinetic Energy (KE) used as a metric of the physical environment. Higher KE leads to increased offshore advection of larvae, reduced total recruitment and connectivity of MPAs. Our results indicate high and uni-directional connectivity between MPAs from north to south influenced by the Brazil Current flowing in the same direction. Results also showed that some MPAs act predominantly as “sink” while others are mainly “source” areas.

Fisheries Closed Areas Strengthen Scallop Larval Settlement and Connectivity Among Closed Areas and Across International Open Fishing Grounds: A Model Study

This study examined whether a measured increase in average body size of adult sea scallops inside three fishery closed areas on Georges Bank (GB), United States (US), was sufficient to increase larval supply to closed areas and open fishing areas in both US and Canadian areas of the Bank. The effects of adult scallop density-at-size and fecundity-at-size on egg production were compared among open and closed fishery areas, countries, and time periods before and after the closed areas were established. Estimated egg production was then used to define spawning conditions in a coupled biological-physical larval tracking model that simulated larval development, mortality, and dispersal. Results showed that order of magnitude increases in larval settlement after closure were facilitated by increases in size-dependant egg production inside and dispersal from Closed Areas I and II, but not Nantucket Lightship Closed Area. The distributions of both egg production and larval settlement became more uniform across the Bank, causing the relative contribution of Canadian larvae to US scallop aggregations to decrease after establishment of Closed Areas I and II. Decreases in small and medium-sized scallop density in Canada and decreases in large scallops over the US-Southern Flank after closure caused local declines in egg production but were not sufficient to negatively affect larval settlement at the regional scale. Our model suggests that the establishment of fishery closed areas on GB considerably strengthened larval supply and settlement within and among several adult scallop aggregations.

Sun-Compass Orientation in Mediterranean Fish Larvae

Mortality is very high during the pelagic larval phase of fishes but the factors that determine recruitment success remain unclear and hard to predict. Because of their bipartite life history, larvae of coastal species have to head back to the shore at the end of their pelagic episode, to settle. These settlement-stage larvae are known to display strong sensory and motile abilities, but most work has been focused on tropical, insular environments and on the influence of coast-related cues on orientation. In this study we quantified the in situ orientation behavior of settlement-stage larvae in a temperate region, with a continuous coast and a dominant along-shore current, and inspected both coast-dependent and independent cues. We tested six species: one Pomacentridae, Chromis chromis, and five Sparidae, Boops boops, Diplodus annularis, Oblada melanura, Spicara smaris and Spondyliosoma cantharus. Over 85% of larvae were highly capable of keeping a bearing, which is comparable to the orientation abilities of tropical species. Sun-related cues influenced the precision of bearing-keeping at individual level. Three species, out of the four tested in sufficient numbers, oriented significantly relative to the sun position. These are the first in situ observations demonstrating the use of a sun compass for orientation by wild-caught settlement-stage larvae. This mechanism has potential for large-scale orientation of fish larvae globally.

Variation in responses of fishes across multiple reserves within a network of marine protected areas in temperate waters

Meta-analyses of field studies have shown that biomass, density, species richness, and size of organisms protected by no-take marine reserves generally increase over time. The magnitude and timing of changes in these response variables, however, vary greatly and depend upon the taxonomic groups protected, size and type of reserve, oceanographic regime, and time since the reserve was implemented. We conducted collaborative, fishery-independent surveys of fishes for seven years in and near newly created marine protected areas (MPAs) in central California, USA. Results showed that initially most MPAs contained more and larger fishes than associated reference sites, likely due to differences in habitat quality. The differences between MPAs and reference sites did not greatly change over the seven years of our study, indicating that reserve benefits will be slow to accumulate in California’s temperate eastern boundary current. Fishes in an older reserve that has been closed to fishing since 1973, however, were significantly more abundant and larger than those in associated reference sites. This indicates that reserve benefits are likely to accrue in the California Current ecosystem, but that 20 years or more may be needed to detect significant changes in response variables that are due to MPA implementation. Because of the high spatial and temporal variability of fish recruitment patterns, long-term monitoring is needed to identify positive responses of fishes to protection in the diverse set of habitats in a dynamic eastern boundary current. Qualitative estimates of response variables, such as would be obtained from an expert opinion process, are unlikely to provide an accurate description of MPA performance. Similarly, using one species or one MPA as an indicator is unlikely to provide sufficient resolution to accurately describe the performance of multiple MPAs.

Tracking larvae with molecular markers reveals high relatedness and early seasonal recruitment success in a partially spawning marine bivalve

The partial synchronized spawning strategy adopted by some marine invertebrate broadcast-spawners can lead to the production of many distinct pools of larvae within a single reproductive cycle. Following the fate of these larval groups from birth to settlement with molecular markers might shed light on mechanisms regulating their population recruitment. Larvae and recruits of Mya arenaria, a partially spawning marine bivalve, were monitored and collected over 13 consecutive weeks during an entire reproductive cycle. Each sampled individual (n = 218) was sorted according to size (early veligers, late veligers, post-larval recruits) and genotyped at seven microsatellite loci for comparisons among samples and with adult reference samples (n = 270). While traditional differentiation statistics (e.g., pairwise Δ(ST), allelic richness) suggested the absence of sweepstakes reproductive success, the level of relatedness found within and among larvae and recruit samples suggested otherwise. Four samples out of ten were observed to have higher within-sample relatedness values than randomly expected, including the very first group of early veligers produced in the season (E1) and the last group of post-larvae who survived recruitment (P10). E1 individuals were also found to be more related than randomly expected to individuals of more than 80% of all other samples including the last surviving recruits (P8 and P10). These results suggest that the first larvae produced in the season were the most successful to survive recruitment. Results also show direct evidence for larval retention and demonstrate for the first time larval and post-larval kin aggregation in a marine bivalve.

Delineating biophysical environments of the Sunda Banda Seascape, Indonesia

The Sunda Banda Seascape (SBS), located in the center of the Coral Triangle, is a global center of marine biodiversity and a conservation priority. We proposed the first biophysical environmental delineation of the SBS using globally available satellite remote sensing and model-assimilated data to categorize this area into unique and meaningful biophysical classes. Specifically, the SBS was partitioned into eight biophysical classes characterized by similar sea surface temperature, chlorophyll a concentration, currents, and salinity patterns. Areas within each class were expected to have similar habitat types and ecosystem functions. Our work supplemented prevailing global marine management schemes by focusing in on a regional scale with finer spatial resolution. It also provided a baseline for academic research, ecological assessments and will facilitate marine spatial planning and conservation activities in the area. In addition, the framework and methods of delineating biophysical environments we presented can be expanded throughout the whole Coral Triangle to support research and conservation activities in this important region.

Marine population connectivity: reconciling large-scale dispersal and high self-retention

Predicting connectivity patterns in systems with fluid transport requires descriptions of the spatial distribution of propagules. In contrast to research on terrestrial seed dispersal, where much attention has focused on localized physical factors affecting dispersal, studies of oceanic propagule dispersal have often emphasized the role of large-scale factors. We link these two perspectives by exploring how propagule dispersal in the ocean is influenced by the "coastal boundary layer" (CBL), a region of reduced velocities near the shoreline that might substantially modify local-scale dispersal. We used a simple simulation model to demonstrate that accounting for the CBL markedly alters transport distances, the widths of dispersal distributions, and the fraction of larvae retained near their sites of origin (self-retention). Median dispersal distances were up to 59% shorter in simulations with a CBL than in those without. Self-retention of larvae increased by up to 3 orders of magnitude in the presence of CBLs, but only minor changes arose in the long-distance tails of the distributions, resulting in asymmetric, non-Gaussian kernels analogous to those quantified for terrestrial seed dispersal. Because successfully settling larvae are commonly those that remain close to shore and interact with the CBL, ignoring this pervasive oceanographic feature will substantially alter predictions of population self-persistence, estimates of connectivity, and outcomes of metapopulation analyses.

Lack of population genetic differentiation of a marine ovoviviparous fish Sebastes schlegelii in Northwestern Pacific

Sebastes schlegelii is one of the fishes that aggregate around drifting seaweed during early development. To examine the population genetic structure of S. schlegelii, a 452-bp fragment of the mtDNA control region was sequenced and used to interpret life history characteristics and larval dispersal strategy. Two-hundred and twenty-one individuals from 13 sites across the entire range of S. schlegelii in China, Japan and Korea were analyzed. A neighbor-joining tree and network showed that there were no significant genealogical structures corresponding to sampling locations. AMOVA, pair-wise FST and exact test revealed no significant genetic differentiation among locations. The migration rate among locations was high based on the result of LAMARC. We conclude that larval dispersal with drifting seaweed and the current environmental factors may play an important role in shaping the contemporary phylogeographic pattern and genetic homogeneity of S. schlegelii.

Population connectivity shifts at high frequency within an open-coast marine protected area network

A complete understanding of population connectivity via larval dispersal is of great value to the effective design and management of marine protected areas (MPA). However empirical estimates of larval dispersal distance, self-recruitment, and within season variability of population connectivity patterns and their influence on metapopulation structure remain rare. We used high-resolution otolith microchemistry data from the temperate reef fish Hypsypops rubicundus to explore biweekly, seasonal, and annual connectivity patterns in an open-coast MPA network. The three MPAs, spanning 46 km along the southern California coastline were connected by larval dispersal, but the magnitude and direction of connections reversed between 2008 and 2009. Self-recruitment, i.e. spawning, dispersal, and settlement to the same location, was observed at two locations, one of which is a MPA. Self-recruitment to this MPA ranged from 50–84%; within the entire 60 km study region, self-recruitment accounted for 45% of all individuals settling to study reefs. On biweekly time scales we observed directional variability in alongshore current data and larval dispersal trajectories; if viewed in isolation these data suggest the system behaves as a source-sink metapopulation. However aggregate biweekly data over two years reveal a reef network in which H. rubicundus behaves more like a well-mixed metapopulation. As one of the few empirical studies of population connectivity within a temperate open coast reef network, this work can inform the MPA design process, implementation of ecosystem based management plans, and facilitate conservation decisions.

Kinship analyses identify fish dispersal events on a temperate coastline

Connectivity is crucial for the persistence and resilience of marine species, the establishment of networks of marine protected areas and the delineation of fishery management units. In the marine environment, understanding connectivity is still a major challenge, due to the technical difficulties of tracking larvae. Recently, parentage analysis has provided a means to address this question effectively. To be effective, this method requires limited adult movement and extensive sampling of parents, which is often not possible for marine species. An alternative approach that is less sensitive to constraints in parental movement and sampling could be the reconstruction of sibships. Here, we directly measure connectivity and larval dispersal in a temperate marine ecosystem through both analytical approaches. We use data from 178 single nucleotide polymorphism markers to perform parentage and sibship reconstruction of the black-faced blenny (Tripterygion delaisi) from an open coastline in the Mediterranean Sea. Parentage analysis revealed a decrease in dispersal success in the focal area over 1 km distance and approximately 6.5% of the juveniles were identified as self-recruits. Sibship reconstruction analysis found that, in general, full siblings did not recruit together to the same location, and that the largest distance between recruitment locations was much higher (11.5 km) than found for parent-offspring pairs (1.2 km). Direct measurements of dispersal are essential to understanding connectivity patterns in different marine habitats, and show the degree of self-replenishment and sustainability of populations of marine organisms. We demonstrate that sibship reconstruction allows direct measurements of dispersal and family structure in marine species while being more easily applied in those species for which the collection of the parental population is difficult or unfeasible.

Larval dispersal modeling of pearl oyster Pinctada margaritifera following realistic environmental and biological forcing in Ahe atoll lagoon

Studying the larval dispersal of bottom-dwelling species is necessary to understand their population dynamics and optimize their management. The black-lip pearl oyster (Pinctada margaritifera) is cultured extensively to produce black pearls, especially in French Polynesia's atoll lagoons. This aquaculture relies on spat collection, a process that can be optimized by understanding which factors influence larval dispersal. Here, we investigate the sensitivity of P. margaritifera larval dispersal kernel to both physical and biological factors in the lagoon of Ahe atoll. Specifically, using a validated 3D larval dispersal model, the variability of lagoon-scale connectivity is investigated against wind forcing, depth and location of larval release, destination location, vertical swimming behavior and pelagic larval duration (PLD) factors. The potential connectivity was spatially weighted according to both the natural and cultivated broodstock densities to provide a realistic view of connectivity. We found that the mean pattern of potential connectivity was driven by the southwest and northeast main barotropic circulation structures, with high retention levels in both. Destination locations, spawning sites and PLD were the main drivers of potential connectivity, explaining respectively 26%, 59% and 5% of the variance. Differences between potential and realistic connectivity showed the significant contribution of the pearl oyster broodstock location to its own dynamics. Realistic connectivity showed larger larval supply in the western destination locations, which are preferentially used by farmers for spat collection. In addition, larval supply in the same sectors was enhanced during summer wind conditions. These results provide new cues to understanding the dynamics of bottom-dwelling populations in atoll lagoons, and show how to take advantage of numerical models for pearl oyster management.

Implications of macroalgal isolation by distance for networks of marine protected areas

The global extent of macroalgal forests is declining, greatly affecting marine biodiversity at broad scales through the effects macroalgae have on ecosystem processes, habitat provision, and food web support. Networks of marine protected areas comprise one potential tool that may safeguard gene flow among macroalgal populations in the face of increasing population fragmentation caused by pollution, habitat modification, climate change, algal harvesting, trophic cascades, and other anthropogenic stressors. Optimal design of protected area networks requires knowledge of effective dispersal distances for a range of macroalgae. We conducted a global meta-analysis based on data in the published literature to determine the generality of relation between genetic differentiation and geographic distance among macroalgal populations. We also examined whether spatial genetic variation differed significantly with respect to higher taxon, life history, and habitat characteristics. We found clear evidence of population isolation by distance across a multitude of macroalgal species. Genetic and geographic distance were positively correlated across 49 studies; a modal distance of 50-100 km maintained F(ST) < 0.2. This relation was consistent for all algal divisions, life cycles, habitats, and molecular marker classes investigated. Incorporating knowledge of the spatial scales of gene flow into the design of marine protected area networks will help moderate anthropogenic increases in population isolation and inbreeding and contribute to the resilience of macroalgal forests.

Beyond connectivity: how empirical methods can quantify population persistence to improve marine protected-area design

Demographic connectivity is a fundamental process influencing the dynamics and persistence of spatially structured populations. Consequently, quantifying connectivity is essential for properly designing networks of protected areas so that they achieve their core ecological objective of maintaining population persistence. Recently, many empirical studies in marine systems have provided essential, and historically challenging to obtain, data on patterns of larval dispersal and export from marine protected areas (MPAs). Here, we review the empirical studies that have directly quantified the origins and destinations of individual larvae and assess those studies' relevance to the theory of population persistence and MPA design objectives. We found that empirical studies often do not measure or present quantities that are relevant to assessing population persistence, even though most studies were motivated or contextualized by MPA applications. Persistence of spatial populations, like nonspatial populations, depends on replacement, whether individuals reproduce enough in their lifetime to replace themselves. In spatial populations, one needs to account for the effect of larval dispersal on future recruitment back to the local population through local retention and other connectivity pathways. The most commonly reported descriptor of larval dispersal was the fraction of recruitment from local origin (self-recruitment). Self-recruitment does not inform persistence-based MPA design because it is a fraction of those arriving, not a fraction of those leaving (local retention), so contains no information on replacement. Some studies presented connectivity matrices, which can inform assessments of persistence with additional knowledge of survival and fecundity after recruitment. Some studies collected data in addition to larval dispersal that could inform assessments of population persistence but which were not presented in that way. We describe how three pieces of empirical information are needed to fully describe population persistence in a network of MPAs: (1) lifetime fecundity, (2) the proportion of larvae that are locally retained (or the full connectivity matrix), and (3) survival rate after recruitment. We conclude by linking theory and data to provide detailed guidance to empiricists and practitioners on field sampling design and data presentation that better informs the MPA objective of population persistence.

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

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

Evaluation of rockfish conservation area networks in the United States and Canada relative to the dispersal distance for black rockfish (Sebastes melanops)

Marine reserves networks are implemented as a way to mitigate the impact of fishing on marine ecosystems. Theory suggests that a reserve network will function synergistically when connected by dispersal, but the scale of dispersal is often unknown. On the Pacific coast of the United States and Canada, both countries have recently implemented a number of rockfish conservation areas (RCAs) to protect exploited rockfish species, but no study has evaluated the connectivity within networks in each country or between the two countries. We used isolation-by-distance theory to estimate the scale of dispersal from microsatellite data in the black rockfish, Sebastes melanops, and compared this estimate with the distance between RCAs that would protect this species. Within each country, we found that the distance between RCAs was generally within the confidence intervals of mean dispersal per generation. The distance between these two RCA networks, however, was greater than the average dispersal per generation. The data were also consistent with a genetic break between southern Oregon and central Oregon. We discuss whether additional nearshore RCAs in southern Oregon and Washington would help promote connectivity between RCA's for shallow-water rockfishes.

Recruitment and resilience of a harvested Caribbean octocoral

Disturbance events are an important component of the ecology of coral reefs and increasingly frequent disturbances coupled with a lack of population resilience may contribute to changes in the structure of coral reef communities. The harvest of the Caribbean octocoral Antillogorgia elisabethae provides an opportunity to explore the relationship between adult abundance and recruitment and the manner in which recruitment contributes to the resilience of local populations. Recruitment of A. elisabethae was monitored in 20, 1-m² quadrats at 8 sites along the southern edge of the Little Bahama Bank from 2004 through 2007. A. elisabethae has been harvested in The Bahamas for over fifteen years and all of the sites had been harvested three times, including a harvest during the course of the study. Abundances of adult colonies at those sites as well as a location that had not been harvested were also determined. Recruitment was highly variable, differing between sites, transects within sites, and, depending on the site, between years. Recruitment was best correlated with adult abundance averaged across the surrounding site. Regression analyses suggest abundance on smaller scales had only small effects on recruitment. The effects of the harvesting were site specific ranging from a 38 to 67% reduction in the density of mature colonies. The sites with the most abundant A. elisabethae continued to have the highest abundances after harvesting and there was no significant difference in recruitment before and after harvesting. Population size-structure at 6 of 8 sites that have been harvested multiple times exhibited an overall depletion in small colonies suggesting long term suppression of recruitment and declining populations. Severe depression of adult abundances coupled with local recruitment can create a negative feedback and lead to the decline of local populations. Populations that are dependent on self-recruitment are not resilient to large disturbance events.

Designing effective reserve networks for nonequilibrium metacommunities

The proliferation of efficient fishing practices has promoted the depletion of commercial stocks around the world and caused significant collateral damage to marine habitats. Recent empirical studies have shown that marine reserves can play an important role in reversing these effects. Equilibrium metapopulation models predict that networks of marine reserves can provide similar benefits so long as individual reserves are sufficiently large to achieve self-sustainability, or spaced based on the extent of dispersal of the target species in order to maintain connectivity between neighboring reserves. However, these guidelines have not been tested in nonequilibrium metacommunity models that exhibit the kinds of complex spatiotemporal dynamics typically seen in natural marine communities. Here, we used a spatially explicit predator-prey model whose predictions have been validated in a marine system to show that current guidelines are not optimal for metacommunities. In equilibrium metacommunities, there is a community-level trade-off for designing effective reserves: Networks whose size and spacing are smaller than the extent of dispersal maximize global predator abundance but minimize global prey abundance because of trophic cascades, whereas the converse is true for reserve networks whose size and spacing are larger than the extent of dispersal. In nonequilibrium metacommunities, reserves whose size and spacing match the extent of spatial autocorrelation in adult abundance (i.e., the extent of patchiness) escape this community-level trade-off by maximizing global abundance and persistence of both the prey and the predator. Overall, these results suggest that using the extent of adult patchiness instead of the extent of larval dispersal as the size and spacing of reserve networks is critical for designing community-based management strategies. By emphasizing patchiness over dispersal distance, our results show how the apparent complexity of nonequilibrium communities can actually simplify management guidelines and reduce uncertainty associated with the assessment of dispersal in marine environments.

Contextualising the Last Survivors: Population Structure of Marine Turtles in the Dominican Republic

Nesting by three species of marine turtles persists in the Dominican Republic, despite historic threats and long-term population decline. We conducted a genetic survey of marine turtles in the Dominican Republic in order to link them with other rookeries around the Caribbean. We sequenced a 740bp fragment of the control region of the mitochondrial DNA of 92 samples from three marine turtle species [hawksbill (n = 48), green (n = 2) and leatherback (n = 42)], and incorporated published data from other nesting populations and foraging grounds. The leatherback turtle (Dermochelys coriacea) in the Dominican Republic appeared to be isolated from Awala-Yalimapo, Cayenne, Trinidad and St. Croix but connected with other Caribbean populations. Two distinct nesting populations of hawksbill turtles (Eremochelys imbricata) were detected in the Dominican Republic and exhibited interesting patterns of connectivity with other nesting sites and juvenile and adult male foraging aggregations. The green sea turtle (Chelonia mydas) has almost been extirpated from the Dominican Republic and limited inference could be made from our samples. Finally, results were compared with Lagrangian drifting buoys and published Lagrangian virtual particles that travelled through the Dominican Republic and Caribbean waters. Conservation implications of sink-source effects or genetic isolation derived from these complex inter-connections are discussed for each species and population.

Larval connectivity and the international management of fisheries

Predicting the oceanic dispersal of planktonic larvae that connect scattered marine animal populations is difficult, yet crucial for management of species whose movements transcend international boundaries. Using multi-scale biophysical modeling techniques coupled with empirical estimates of larval behavior and gamete production, we predict and empirically verify spatio-temporal patterns of larval supply and describe the Caribbean-wide pattern of larval connectivity for the Caribbean spiny lobster (Panulirus argus), an iconic coral reef species whose commercial value approaches $1 billion USD annually. Our results provide long sought information needed for international cooperation in the management of marine resources by identifying lobster larval connectivity and dispersal pathways throughout the Caribbean. Moreover, we outline how large-scale fishery management could explicitly recognize metapopulation structure by considering larval transport dynamics and pelagic larval sanctuaries.

Microsatellite analyses of blacktip reef sharks (Carcharhinus melanopterus) in a fragmented environment show structured clusters

The population dynamics of shark species are generally poorly described because highly mobile marine life is challenging to investigate. Here we investigate the genetic population structure of the blacktip reef shark (Carcharhinus melanopterus) in French Polynesia. Five demes were sampled from five islands with different inter-island distances (50–1500 km). Whether dispersal occurs between islands frequently enough to prevent moderate genetic structure is unknown. We used 11 microsatellites loci from 165 individuals and a strong genetic structure was found among demes with both F-statistics and Bayesian approaches. This differentiation is correlated with the geographic distance between islands. It is likely that the genetic structure seen is the result of all or some combination of the following: low gene flow, time since divergence, small effective population sizes, and the standard issues with the extent to which mutation models actually fit reality. We suggest low levels of gene flow as at least a partial explanation of the level of genetic structure seen among the sampled blacktip demes. This explanation is consistent with the ecological traits of blacktip reef sharks, and that the suitable habitat for blacktips in French Polynesia is highly fragmented. Evidence for spatial genetic structure of the blacktip demes we studied highlights that similar species may have populations with as yet undetected or underestimated structure. Shark biology and the market for their fins make them highly vulnerable and many species are in rapid decline. Our results add weight to the case that total bans on shark fishing are a better conservation approach for sharks than marine protected area networks.

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.

Functional endemism: population connectivity, shifting baselines, and the scale of human experience

Quantifying population connectivity is important for visualizing the spatial and temporal scales that conservation measures act upon. Traditionally, migration based on genetic data has been reported in migrants per generation. However, the temporal scales over which this migration may occur do not necessarily accommodate the scales over which human perturbations occur, leaving the potential for a disconnect between population genetic data and conservation action based on those data. Here, we present a new metric called the "Rule of Memory", which helps conservation practitioners to interpret "migrants per generation" in the context both of human modified ecosystems and the cultural memory of those doing the modification. Our rule states that clades should be considered functionally endemic regardless of their actual taxonomic designation if the migration between locations is insufficient to maintain a viable population over the timescales of one human generation (20 years). Since larger animals are more likely to be remembered, we quantify the relationship between migrants per human (N) and body mass of the organism in question (M) with the formula N = 10M(-1). We then use the coral reef fish Pomacentrus moluccensis to demonstrate the taxonomic and spatial scales over which this rule can be applied. Going beyond minimum viable population literature, this metric assesses the probability that a clade's existence will be forgotten by people throughout its range during a period of extirpation. Because conservation plans are predicated on having well-established baselines, a loss of a species over the range of one human generation evokes the likelihood of that species no longer being recognized as a member of an ecosystem, and thus being excluded in restoration or conservation prioritization. [Correction added on 26 December 2012, after first online publication: this formula has been corrected to N=10M(-1)].

Assessing dispersal patterns of fish propagules from an effective mediterranean marine protected area

Successfully enforced marine protected areas (MPAs) have been widely demonstrated to allow, within their boundaries, the recovery of exploited species and beyond their boundaries, the spillover of juvenile and adult fish. Little evidence is available about the so-called ‘recruitment subsidy’, the augmented production of propagules (i.e. eggs and larvae) due to the increased abundance of large-sized spawners hosted within effective MPAs. Once emitted, propagules can be locally retained and/or exported elsewhere. Patterns of propagule retention and/or export from MPAs have been little investigated, especially in the Mediterranean. This study investigated the potential for propagule production and retention/export from a Mediterranean MPA (Torre Guaceto, SW Adriatic Sea) using the white sea bream, Diplodus sargus sargus, as a model species. A multidisciplinary approach was used combining 1) spatial distribution patterns of individuals (post-settlers and adults) assessed through visual census within Torre Guaceto MPA and in northern and southern unprotected areas, 2) Lagrangian simulations of dispersal based on an oceanographic model of the region and data on early life-history traits of the species (spawning date, pelagic larval duration) and 3) a preliminary genetic study using microsatellite loci. Results show that the MPA hosts higher densities of larger-sized spawners than outside areas, potentially guaranteeing higher propagule production. Model simulations and field observation suggest that larval retention within and long-distance dispersal across MPA boundaries allow the replenishment of the MPA and of exploited populations up to 100 km down-current (southward) from the MPA. This pattern partially agrees with the high genetic homogeneity found in the entire study area (no differences in genetic composition and diversity indices), suggesting a high gene flow. By contributing to a better understanding of propagule dispersal patterns, these findings provide crucial information for the design of MPAs and MPA networks effective to replenish fish stocks and enhance fisheries in unprotected areas.

Direct genetic evidence for reproductive philopatry and associated fine-scale migrations in female blacktip reef sharks (Carcharhinus melanopterus) in French Polynesia

Conservation of top predators has been emphasized as essential in an ecosystem due to their role in trophic chain regulation. Optimizing conservation strategies for these endangered marine top predators requires direct estimates of breeding patterns and connectivity as these are essential to understanding the population dynamics. There have been some attempts to investigate breeding patterns of reef sharks from litter reconstruction using molecular analyses. However, direct fine-scale migrations of female sharks for parturition as well as connectivity at a medium scale like between islands remain mostly unknown. We used microsatellite DNA markers and a likelihood-based parentage analysis to determine breeding patterns of female blacktip reef sharks in Moorea (Society Islands, French Polynesia). Most females gave birth at their home island but some migrated to specific nursery areas outside the area they are attached to, sometimes going to another island 50 km away across deep ocean. Our analysis also revealed that females migrated to the same nursery for every birthing event. Many offspring showed a high level of inbreeding indicating an overall reduced population size, restricted movements and dispersal, or specific mating behaviour. Females represent the vectors that transport the genes at nursery grounds, and their fidelity should thus define reproductive units. As females seem to be philopatric, males could be the ones dispersing genes between populations. These results highlight the need to conserve coastal zones where female reef sharks seem to exhibit philopatry during the breeding season.

Persistence of self-recruitment and patterns of larval connectivity in a marine protected area network

The use of marine protected area (MPA) networks to sustain fisheries and conserve biodiversity is predicated on two critical yet rarely tested assumptions. Individual MPAs must produce sufficient larvae that settle within that reserve's boundaries to maintain local populations while simultaneously supplying larvae to other MPA nodes in the network that might otherwise suffer local extinction. Here, we use genetic parentage analysis to demonstrate that patterns of self-recruitment of two reef fishes (Amphiprion percula and Chaetodon vagabundus) in an MPA in Kimbe Bay, Papua New Guinea, were remarkably consistent over several years. However, dispersal from this reserve to two other nodes in an MPA network varied between species and through time. The stability of our estimates of self-recruitment suggests that even small MPAs may be self-sustaining. However, our results caution against applying optimization strategies to MPA network design without accounting for variable connectivity among species and over time.

Reproductive output and duration of the pelagic larval stage determine seascape-wide connectivity of marine populations

Connectivity among marine populations is critical for persistence of metapopulations, coping with climate change, and determining the geographic distribution of species. The influence of pelagic larval duration (PLD) on connectivity has been studied extensively, but relatively little is known about the influence of other biological parameters, such as the survival and behavior of larvae, and the fecundity of adults, on population connectivity. Furthermore, the interaction between the seascape (habitat structure and currents) and these biological parameters is unclear. We explore these interactions using a biophysical model of larval dispersal across the Indo-Pacific. We describe an approach that quantifies geographic patterns of connectivity from demographically relevant to evolutionarily significant levels across a range of species. We predict that at least 95% of larval settlement occurs within 155 km of the source population and within 13 days irrespective of the species' life history, yet long-distant connections remain likely. Self-recruitment is primarily driven by the local oceanography, larval mortality, and the larval precompetency period, whereas broad-scale connectivity is strongly influenced by reproductive output (abundance and fecundity of adults) and the length of PLD. The networks we have created are geographically explicit models of marine connectivity that define dispersal corridors, barriers, and the emergent structure of marine populations. These models provide hypotheses for empirical testing.