Connectivity of Caribbean coral populations: complementary insights from empirical and modelled gene flow

Understanding patterns of connectivity among populations of marine organisms is essential for the development of realistic, spatially explicit models of population dynamics. Two approaches, empirical genetic patterns and oceanographic dispersal modelling, have been used to estimate levels of evolutionary connectivity among marine populations but rarely have their potentially complementary insights been combined. Here, a spatially realistic Lagrangian model of larval dispersal and a theoretical genetic model are integrated with the most extensive study of gene flow in a Caribbean marine organism. The 871 genets collected from 26 sites spread over the wider Caribbean subsampled 45.8% of the 1900 potential unique genets in the model. At a coarse scale, significant consensus between modelled estimates of genetic structure and empirical genetic data for populations of the reef-building coral Montastraea annularis is observed. However, modelled and empirical data differ in their estimates of connectivity among northern Mesoamerican reefs indicating that processes other than dispersal may dominate here. Further, the geographic location and porosity of the previously described east-west barrier to gene flow in the Caribbean is refined. A multi-prong approach, integrating genetic data and spatially realistic models of larval dispersal and genetic projection, provides complementary insights into the processes underpinning population connectivity in marine invertebrates on evolutionary timescales.     

Sharp genetic breaks among populations of Haptosquilla pulchella (Stomatopoda) indicate limits to larval transport: patterns, causes, and consequences

To help stem the precipitous decline of coral reef ecosystems world-wide, conservation efforts are focused on establishing interconnected reserve networks to protect threatened populations. Because many coral reef organisms have a planktonic or pelagic larval dispersal phase, it is critical to understand the patterns of ecological connectivity between reserve populations that result from larval dispersal. We used genetics to infer dispersal patterns among 24 Indo-West Pacific populations of the mantis shrimp, Haptosquilla pulchella. Contrary to predictions of high dispersal facilitated by the strong currents of the Indonesian throughflow, mitochondrial DNA sequences from 393 individuals displayed striking patterns of regional genetic differentiation concordant with ocean basins isolated during periods of lowered sea level. Patterns of genetic structuring indicate that although dispersal within geographical regions with semicontiguous coastlines spanning thousands of kilometres may be common, ecologically meaningful connections can be rare among populations separated by as little as 300 km of open ocean. Strong genetic mosaics in a species with high dispersal potential highlight the utility of genetics for identifying regional patterns of genetic connectivity between marine populations and show that the assumption that ocean currents will provide ecological connectivity among marine populations must be empirically tested in the design of marine reserve networks.     

Local retention,dispersal and fluctuating connectivity among populations of a coral reef fish

The persistence and resilience of marine populations in the face of disturbances is directly affected by connectivity among populations. Thus, understanding the magnitude and pattern of connections among populations and the temporal variation in these patterns is critical for the effective management and conservation of marine species. Despite recent advances in our understanding of marine connectivity, few empirical studies have directly measured the magnitude or pattern of connections among populations of marine fishes, and none have explicitly investigated temporal variation in demographic connectivity. We use genetic assignment tests to track the dispersal of 456 individual larval fishes to quantify the extent of connectivity, dispersal, self-recruitment and local retention within and among seven populations of a coral reef fish (Stegastes partitus) over a three-year period. We found that some larvae do disperse long distances (~200 km); however, self-recruitment was a regular phenomenon. Importantly, we found that dispersal distances, self-recruitment, local retention and the pattern of connectivity varied significantly among years. Our data highlight the unpredictable nature of connectivity, and underscore the need for more, temporally replicated, empirical measures of connectivity to inform management decisions.     

Testing the consistency of connectivity patterns for a widely dispersing marine species

Connectivity is widely recognized as an important component in developing effective management and conservation strategies. Although managers are generally most interested in demographic, rather than genetic connectivity, new analytic approaches are able to provide estimates of both demographic and genetic connectivity measures from genetic data. Combining such genetic data with mathematical models represents a powerful approach for accurately determining patterns of population connectivity. Here, we use microsatellite markers to investigate the genetic population structure of the New Zealand Rock Lobster, Jasus edwardsii, which has one of the longest known larval durations of all marine species (>2 years), a very large geographic range (>5500 km), and has been the subject of extensive dispersal modeling. Despite earlier mitochondrial DNA studies finding homogeneous genetic structure, the mathematical model suggests that there are source-sink dynamics for this species. We found evidence of genetic structure in J. edwardsii populations with three distinct genetic groups across New Zealand and a further Australian group; these groups and patterns of gene flow were generally congruent with the earlier mathematical model. Of particular interest was the consistent identification of a self-recruiting population/region from both modeling and genetic approaches. Although there is the potential for selection and harvesting to influence the patterns we observed, we believe oceanographic processes are most likely responsible for the genetic structure observed in J. edwardsii. Our results, using a species at the extreme end of the dispersal spectrum, demonstrate that source-sink population dynamics may still exist for such species.     

Isolation by resistance across a complex coral reef seascape

A detailed understanding of the genetic structure of populations and an accurate interpretation of processes driving contemporary patterns of gene flow are fundamental to successful spatial conservation management. The field of seascape genetics seeks to incorporate environmental variables and processes into analyses of population genetic data to improve our understanding of forces driving genetic divergence in the marine environment. Information about barriers to gene flow (such as ocean currents) is used to define a resistance surface to predict the spatial genetic structure of populations and explain deviations from the widely applied isolation-by-distance model. The majority of seascape approaches to date have been applied to linear coastal systems or at large spatial scales (more than 250 km), with very few applied to complex systems at regional spatial scales (less than 100 km). Here, we apply a seascape genetics approach to a peripheral population of the broadcast-spawning coral Acropora spicifera across the Houtman Abrolhos Islands, a high-latitude complex coral reef system off the central coast of Western Australia. We coupled population genetic data from a panel of microsatellite DNA markers with a biophysical dispersal model to test whether oceanographic processes could explain patterns of genetic divergence. We identified significant variation in allele frequencies over distances of less than 10 km, with significant differentiation occurring between adjacent sites but not between the most geographically distant ones. Recruitment probabilities between sites based on simulated larval dispersal were projected into a measure of resistance to connectivity that was significantly correlated with patterns of genetic divergence, demonstrating that patterns of spatial genetic structure are a function of restrictions to gene flow imposed by oceanographic currents. This study advances our understanding of the role of larval dispersal on the fine-scale genetic structure of coral populations across a complex island system and applies a methodological framework that can be tailored to suit a variety of marine organisms with a range of life-history characteristics.     

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

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

Larval retention and connectivity among populations of corals and reef fishes: history,advances and challenges

The extent of larval dispersal on coral reefs has important implications for the persistence of coral reef metapopulations, their resilience and recovery from an increasing array of threats, and the success of protective measures. This article highlights a recent dramatic increase in research effort and a growing diversity of approaches to the study of larval retention within (self-recruitment) and dispersal among (connectivity) isolated coral reef populations. Historically, researchers were motivated by alternative hypotheses concerning the processes limiting populations and structuring coral reef assemblages, whereas the recent impetus has come largely from the need to incorporate dispersal information into the design of no-take marine protected area (MPA) networks. Although the majority of studies continue to rely on population genetic approaches to make inferences about dispersal, a wide range of techniques are now being employed, from small-scale larval tagging and paternity analyses, to large-scale biophysical circulation models. Multiple approaches are increasingly being applied to cross-validate and provide more realistic estimates of larval dispersal. The vast majority of empirical studies have focused on corals and fishes, where evidence for both extremely local scale patterns of self-recruitment and ecologically significant connectivity among reefs at scales of tens of kilometers (and in some cases hundreds of kilometers) is accumulating. Levels of larval retention and the spatial extent of connectivity in both corals and fishes appear to be largely independent of larval duration or reef size, but may be strongly influenced by geographic setting. It is argued that high levels of both self-recruitment and larval import can contribute to the resilience of reef populations and MPA networks, but these benefits will erode in degrading reef environments.     

Parentage assignment detects frequent and large-scale dispersal in water voles

Estimating the rate and scale of dispersal is essential for predicting the dynamics of fragmented populations, yet empirical estimates are typically imprecise and often negatively biased. We maximized detection of dispersal events between small, subdivided populations of water voles (Arvicola terrestris) using a novel method that combined direct capture-mark-recapture with microsatellite genotyping to identify parents and offspring in different populations and hence infer dispersal. We validated the method using individuals known from trapping data to have dispersed between populations. Local populations were linked by high rates of juvenile dispersal but much lower levels of adult dispersal. In the spring breeding population, 19% of females and 33% of males had left their natal population of the previous year. The average interpopulation dispersal distance was 1.8 km (range 0.3-5.2 km). Overall, patterns of dispersal fitted a negative exponential function. Information from genotyping increased the estimated rate and scale of dispersal by three- and twofold, respectively, and hence represents a powerful tool to provide more realistic estimates of dispersal parameters.     

Genetic population structure of the endemic fourline wrasse (Larabicus quadrilineatus) suggests limited larval dispersal distances in the Red Sea

The connectivity among marine populations is determined by the dispersal capabilities of adults as well as their eggs and larvae. Dispersal distances and directions have a profound effect on gene flow and genetic differentiation within species. Genetic homogeneity over large areas is a common feature of coral reef fishes and can reflect high dispersal capability resulting in high levels of gene flow. If fish larvae return to their parental reef, gene flow would be restricted and genetic differentiation could occur. Larabicus quadrilineatus (Labridae) is considered as an endemic fish species of the Red Sea and Gulf of Aden. The juveniles of this species are cleaner fish that feed on ectoparasites of other fishes. Here, we investigated the genetic population structure and gene flow in L. quadrilineatus among five locations in the Red Sea to infer connectivity among them. To estimate genetic diversity, we analysed 369 bp of 237 mitochondrial DNA control region sequences. Haplotype and nucleotide diversities were higher in the southern than in the northern Red Sea. Analysis of molecular variance (amova) detected the highest significant genetic variation between northern and central/southern populations (Phi(CT) = 0.01; P < 0.001). Migration analysis revealed a several fold higher northward than southward migration, which could be explained by oceanographic conditions and spawning season. Even though the Phi(ST) value of 0.01 is rather low and implies a long larval dispersal distance, estimates based on the isolation-by-distance model show a very low mean larval dispersal distance (0.44-5.1 km) compared to other studies. In order to enable a sustainable ornamental fishery on the fourline wrasse, the results of this study suggest that populations in the northern and southern Red Sea should be managed separately as two different stocks. The rather low larval dispersal distance of about 5 km needs to be considered in the design of marine protected areas to enable connectivity and self-seeding.     

Model-derived dispersal pathways from multiple source populations explain variability of invertebrate larval supply

Background

Predicting the spatial and temporal patterns of marine larval dispersal and supply is a challenging task due to the small size of the larvae and the variability of oceanographic processes. Addressing this problem requires the use of novel approaches capable of capturing the inherent variability in the mechanisms involved.

Methodology/Principal Findings

In this study we test whether dispersal and connectivity patterns generated from a bio-physical model of larval dispersal of the crab Carcinus maenas, along the west coast of the Iberian Peninsula, can predict the highly variable daily pattern of wind-driven larval supply to an estuary observed during the peak reproductive season (March–June) in 2006 and 2007. Cross-correlations between observed and predicted supply were significant (p<0.05) and strong, ranging from 0.34 to 0.81 at time lags of −6 to +5 d. Importantly, the model correctly predicted observed cross-shelf distributions (Pearson r = 0.82, p<0.001, and r = 0.79, p<0.01, in 2006 and 2007) and indicated that all supply events were comprised of larvae that had been retained within the inner shelf; larvae transported to the outer shelf and beyond never recruited. Estimated average dispersal distances ranged from 57 to 198 km and were only marginally affected by mortality.

Conclusions/Significance

The high degree of predicted demographic connectivity over relatively large geographic scales is consistent with the lack of genetic structuring in C. maenas along the Iberian Peninsula. These findings indicate that the dynamic nature of larval dispersal can be captured by mechanistic biophysical models, which can be used to provide meaningful predictions of the patterns and causes of fine-scale variability in larval supply to marine populations.     

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.     

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.     

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

Increasing dispersal duration should result in increasing dispersal distance, facilitating higher gene flow among populations. As such, it has long been predicted that genetic structure (e.g. F(ST) ) among populations of marine species should be strongly correlated with pelagic larval duration (PLD). However, previous studies have repeatedly shown a surprisingly poor correspondence. This result has been frequently interpreted as evidence for larval behaviours or physical oceanographic processes that result in larvae failing to reach their dispersal potential, or error inherent in estimating PLD and F(ST) . This study employed a computer modelling approach to explore the impacts of various uncertainties on the correlation between measures of genetic differentiation such as F(ST) and PLD. Results indicate that variation resulting from PLD estimation error had minor impacts on the correlation between genetic structure and PLD. However, variation in effective population size between species, errors in F(ST) estimation and non-equilibrium F(ST) values all had major impacts, resulting in dramatically weaker correlations between PLD and F(ST) . These results suggest that poor correlations between PLD and F(ST) may result from variation and uncertainty in the terms associated with the calculation of F(ST) values. As such, PLD may be a much stronger determinant of realized larval dispersal than suggested by the weak-to-moderate correlations between PLD and F(ST) reported in empirical studies.     

Population structure is not a simple function of reproductive mode and larval type: insights from tropical corals

1. For a wide range of organisms, heritable variation in life-history characteristics has been shown to be strongly subject to selection, reflecting the impact that variation in characters such as genotypic diversity, duration of larval development and adaptations for dispersal can have on the fitness of offspring and the make-up of populations. Indeed, variation in life-history characteristics, especially reproduction and larval type, have often been used to predict patterns of dispersal and resultant population structures in marine invertebrates. 2. Scleractinian corals are excellent models with which to test this relationship, as they exhibit almost every possible combination of reproductive mode and larval type. Some general patterns are emerging but, contrary to expectations, genetic data suggest that while populations of broadcast spawning species may be genotypically diverse they may be heavily reliant on localized recruitment rather than widespread dispersal of larvae. 3. Here we use microsatellites to test the importance of localized recruitment by comparing the genetic structure of populations of two broadcast spawning corals with contrasting modes of reproduction and larval development; Goniastrea favulus is self-compatible, has sticky, negatively buoyant eggs and larvae and is expected to have restricted dispersal of gametes and larvae. In contrast, Platygyra daedalea is self-incompatibile, spawns positively buoyant egg-sperm bundles and has planktonic development. 4. Surprisingly, spatial-autocorrelation revealed no fine-scale clustering of similar genotypes within sites for G. favulus, but showed a non-random distribution of genotypes in P. daedalea. Both species showed similar levels of genetic subdivision among sites separated by 50-100 m (F(ST) = 0.03), suggesting that larval dispersal may be equivalent in both species. 5. Interestingly, as fragmentation has been considered rare in massive corals, our sample of 284 P. daedalea colonies included 28 replicated genotypes that were each unlikely (P < 0.05) to have been derived independently from sexual reproduction. 6. We conclude that the extreme life history of G. favulus does not produce unusually fine-scale genetic structure and subsequently, that reproductive mode and larval type may not be not good predictors of population structure or dispersal ability.     

Larval connectivity in an effective network of marine protected areas

Acceptance of marine protected areas (MPAs) as fishery and conservation tools has been hampered by lack of direct evidence that MPAs successfully seed unprotected areas with larvae of targeted species. For the first time, we present direct evidence of large-scale population connectivity within an existing and effective network of MPAs. A new parentage analysis identified four parent-offspring pairs from a large, exploited population of the coral-reef fish Zebrasoma flavescens in Hawai'i, revealing larval dispersal distances ranging from 15 to 184 km. In two cases, successful dispersal was from an MPA to unprotected sites. Given high adult abundances, the documentation of any parent-offspring pairs demonstrates that ecologically-relevant larval connectivity between reefs is substantial. All offspring settled at sites to the north of where they were spawned. Satellite altimetry and oceanographic models from relevant time periods indicated a cyclonic eddy that created prevailing northward currents between sites where parents and offspring were found. These findings empirically demonstrate the effectiveness of MPAs as useful conservation and management tools and further highlight the importance of coupling oceanographic, genetic, and ecological data to predict, validate and quantify larval connectivity among marine populations.     

Weak genetic structure indicates strong dispersal limits: a tale of two coral reef fish

The extent of dispersal by pelagic larvae in marine environments, including coral reefs, is central for understanding local population dynamics and designing sustainable marine reserves. We present here the first example of a clear stepping-stone genetic structure throughout the Caribbean basin for a common coral reef species, the French grunt (Haemulon flavolineatum). Analysis of microsatellite DNA markers indicated that French grunt population structure may be characterized by overlapping populations throughout the Caribbean, influenced by independent population dynamics but with no fixed geographical boundaries. In addition, different spatial genetic patterns were found in different oceanographic regions. A second species, the bluehead wrasse (Thalassoma bifasciatum), has a much longer pelagic larval duration than French grunts and showed no explicit spatial pattern of genetic variation. This finding is concordant with the hypothesis of a positive relationship between larval dispersal and duration in the plankton. While the magnitude of the genetic signal of population structure in French grunts was very low (F(ST) approximately 0.003), the pattern of isolation-by-distance throughout the Caribbean indicated considerable population structure with important ecological and conservation significance.     

Rapid recovery of genetic diversity of stomatopod populations on Krakatau: temporal and spatial scales of marine larval dispersal

Although the recovery of terrestrial communities shattered by the massive eruption of Krakatau in 1883 has been well chronicled, the fate of marine populations has been largely ignored. We examined patterns of genetic diversity in populations of two coral reef-dwelling mantis shrimp, Haptosquilla pulchella and Haptosquilla glyptocercus (Stomatopoda: Protosquillidae), on the islands of Anak Krakatau and Rakata. Genetic surveys of mitochondrial cytochrome oxidase c (subunit 1) in these populations revealed remarkably high levels of haplotypic and nucleotide diversity that were comparable with undisturbed populations throughout the Indo-Pacific. Recolonization and rapid recovery of genetic diversity in the Krakatau populations indicates that larval dispersal from multiple and diverse source populations contributes substantially to the demographics of local populations over intermediate temporal (tens to hundreds of years) and spatial scales (tens to hundreds of kilometres). Natural experiments such as Krakatau provide an excellent mechanism to investigate marine larval dispersal and connectivity. Results from stomatopods indicate that marine reserves should be spaced no more than 50-100 km apart to facilitate ecological connectivity via larval dispersal.     

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

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

Close encounters with eddies: oceanographic features increase growth of larval reef fishes during their journey to the reef

Like most benthic marine organisms, coral reef fishes produce larvae that traverse open ocean waters before settling and metamorphosing into juveniles. Where larvae are transported and how they survive is a central question in marine and fisheries ecology. While there is increasing success in modelling potential larval trajectories, our knowledge of the physical and biological processes contributing to larval survivorship during dispersal remains relatively poor. Mesoscale eddies (MEs) are ubiquitous throughout the world''s oceans and their propagation is often accompanied by upwelling and increased productivity. Enhanced production suggests that eddies may serve as important habitat for the larval stages of marine organisms, yet there is a lack of empirical data on the growth rates of larvae associated with these eddies. During three cruises in the Straits of Florida, we sampled larval fishes inside and outside five cyclonic MEs. Otolith microstructure analysis revealed that four of five species of reef fish examined had consistently faster growth inside these eddies. Because increased larval growth often leads to higher survivorship, larvae that encounter MEs during transit are more likely to contribute to reef populations. Successful dispersal in oligotrophic waters may rely on larval encounter with such oceanographic features.