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

Background

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

Methods

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

Results

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

Conclusions

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

Population genetics and historical demographic inferences of the blue crab Callinectes sapidus in the US based on microsatellites

The native range of the blue crab Callinectes sapidus spans Nova Scotia to northern Argentina. In the US, it constitutes a keystone species in estuarine habitats of the Atlantic coast and Gulf of Mexico (GOM), serving as both predator and prey to other species, and also has historically represented a multi-billion dollar fishery. Knowledge relevant to effective management and monitoring of this ecologically and economically important species, such as levels of population genetic differentiation and genetic diversity, is necessary. Although several population genetics studies have attempted to address these questions in one or more parts of its distribution, conflicting results and potential problems with the markers used, as well as other issues, have obscured our understanding on them. In this study, we examined large-scale genetic connectivity of the blue crab in the US, using 16 microsatellites, and genotyped individuals from Chesapeake Bay, in the US Atlantic, and from nine localities along the US GOM coast. Consistent with the high long-distance dispersal potential of this species, very low levels of genetic differentiation were detected for the blue crab among the ten US localities examined, suggesting it constitutes a large panmictic population within this region. Estimations of genetic diversity for the blue crab appear to be high in the US, and provide a baseline for monitoring temporal changes in this species. Demographic analyses indicate a recent range expansion of the US population, probably during the Holocene. In addition, capitalizing on published microsatellite data from southern Brazil, our analyses detected high genetic differentiation between localities in the US and Brazil. These results point to the need for examination of genetic diversity and differentiation along the area spanning the US to southern Brazil.

A multidisciplinary analytical framework to delineate spawning areas and quantify larval dispersal in coastal fish

Assessing larval dispersal is essential to understand the structure and dynamics of marine populations. However, knowledge about early-life dispersal is sparse, and so is our understanding of the spawning process, perhaps the most obscure component of biphasic life cycles. Indeed, poorly known species-specific spawning modality and species-specific early-life traits, as well as the high spatio-temporal variability of the oceanic circulation experienced during larval drift, hamper our ability to appraise the realized connectivity of coastal fishes. Here, we propose an analytical framework which combines Lagrangian modelling, network theory, otolith analyses and biogeographical information to pinpoint and characterize larval sources which are then grouped into discrete spawning areas. Such well-delineated larval sources allow improving the quantitative evaluations of both dispersal scales and connectivity patterns. To illustrate its added value, our approach is applied to two case-studies focusing on Diplodus sargus and Diplodus vulgaris in the Adriatic sea. We evidence robust correlations between otolith geochemistry and modelled spawning areas to assess their relative importance for the larval replenishment of the Apulian coast. Our results show that, contrary to D. sargus, D. vulgaris larvae originate from both eastern and western Adriatic shorelines. Our findings also suggest that dispersal distances and dispersal surfaces scale differently with the pelagic larval duration. Furthermore, 30.8% of D. sargus larvae and 23.6% of D. vulgaris larvae of the Apulian populations originate from Marine protected area (MPA), exemplifying larval export from MPAs to surrounding unprotected areas. This flexible multidisciplinary framework, which can be adjusted to any coastal fish and oceanic system, exploits the explanatory power of a dispersal model, fine-tuned and backed-up by observations, to provide more reliable scientific basis for the management and conservation of marine ecosystems.

Modeled larval connectivity of a multi-species reef fish and invertebrate assemblage off the coast of Moloka'i, Hawai'i

We use a novel individual-based model (IBM) to simulate larval dispersal around the island of Moloka‘i in the Hawaiian Archipelago. Our model uses ocean current output from the Massachusetts Institute of Technology general circulation model (MITgcm) as well as biological data on four invertebrate and seven fish species of management relevance to produce connectivity maps among sites around the island of Moloka‘i. These 11 species span the range of life history characteristics of Hawaiian coral reef species and show different spatial and temporal patterns of connectivity as a result. As expected, the longer the pelagic larval duration (PLD), the greater the proportion of larvae that disperse longer distances, but regardless of PLD (3–270 d) most successful dispersal occurs either over short distances within an island (<30 km) or to adjacent islands (50–125 km). Again, regardless of PLD, around the island of Moloka‘i, connectivity tends to be greatest among sites along the same coastline and exchange between northward, southward, eastward and westward-facing shores is limited. Using a graph-theoretic approach to visualize the data, we highlight that the eastern side of the island tends to show the greatest out-degree and betweenness centrality, which indicate important larval sources and connectivity pathways for the rest of the island. The marine protected area surrounding Kalaupapa National Historical Park emerges as a potential source for between-island larval connections, and the west coast of the Park is one of the few regions on Moloka‘i that acts as a net larval source across all species. Using this IBM and visualization approach reveals patterns of exchange between habitat regions and highlights critical larval sources and multi-generational pathways to indicate priority areas for marine resource managers.

Life histories predict genetic diversity and population structure within three species of octopus targeted by small-scale fisheries in Northwest Mexico

The fishery for octopus in Northwest Mexico has increased to over 2,000 tons annually, but to date the specific composition of the catch has been ignored. With at least three main species targeted by artisanal fisheries in the region with distinct life histories, the lack of basic biological information about the distribution, metapopulation size and structure of each species could impede effective fisheries management to avoid overexploitation. We tested if different life histories of three species of octopus could help predict observed patterns of genetic diversity, population dynamics, structure and connectivity and how this information could be relevant to the sustainable management of the fishery. We sequenced two mitochondrial genes and genotyped seven nuclear microsatellite loci to identify the distribution of each species in 20 locations from the Gulf of California and the west coast of the Baja California peninsula. We tested five hypotheses derived from population genetic theory based on differences in the fecundity and dispersal potential for each species. We discovered that Octopus bimaculoides with low fecundity and direct development (without a planktonic phase) had lower average effective population size and genetic diversity, but higher levels of kinship, population structure, and richness of private alleles, than the other two species. These features indicated limited dispersal and high local recruitment. In contrast, O. bimaculatus and O. hubbsorum with higher fecundity and planktonic phase as paralarvae had higher effective population size and genetic diversity, and overall lower kinship and population structure than O. bimaculoides. These observations supported higher levels of gene flow over a larger geographical scale. O. bimaculatus with the longest planktonic paralarval duration and therefore larger dispersal potential had differences in the calculated parameters possibly associated with increased connectivity. We propose O. bimaculoides is more susceptible to over exploitation of small, isolated populations and could have longer recovery times than the other two species. This species may benefit from distinct fishery management within each local population. O. bimaculatus and O. hubbsorum may benefit from fishery management that takes into account metapopulation structure over larger geographic scales and the directionality and magnitude of larval dispersal driven by ocean currents and population connectivity among individuals of each locality. The distribution of each species and variations in their reproductive phenology is also important to consider when establishing marine reserves or seasonal fishing closures.

Characterization of 32 microsatellite loci for the Pacific red snapper, Lutjanus peru, through next generation sequencing

We developed a set of hypervariable microsatellite markers for the Pacific red snapper (Lutjanus peru), an economically important marine fish for small-scale fisheries in the west coast of Mexico. We performed shotgun genome sequencing with the 454 XL titanium chemistry and used bioinformatic tools to search for perfect microsatellite loci. We selected 66 primer pairs that were synthesized and genotyped in an ABI PRISM 3730XL DNA sequencer in 32 individuals from the Gulf of California. We estimated levels of genetic diversity, deviations from linkage and Hardy-Weinberg equilibrium, estimated the frequency of null alleles and the probability of individual identity for the new markers. We reanalyzed 16 loci in 16 individuals to estimate genotyping error rates. Eighteen loci failed to amplify, 16 loci were discarded due to unspecific amplifications and 32 loci (14 tetranucleotide and 18 dinucleotide) were successfully scored. The average number of alleles per locus was 21 (±6.87, SD) and ranged from 8 to 34. The average observed and expected heterozygosities were 0.787 (±0.144 SD, range 0.250-0.935) and 0.909 (±0.122 SD, range 0.381-0.965), respectively. No significant linkage was detected. Eight loci showed deviations from Hardy-Weinberg equilibrium, and from these, four loci showed moderate null allele frequencies (0.104-0.220). The probability of individual identity for the new loci was 1.46^-62. Genotyping error rates averaged 9.58%. The new markers will be useful to investigate patterns of larval dispersal, metapopulation dynamics, fine-scale genetic structure and diversity aimed to inform the implementation of spatially explicit fisheries management strategies in the Gulf of California.

Assessing larval connectivity for marine spatial planning in the Adriatic

There are plans to start building offshore marine renewable energy devices throughout the Mediterranean and the Adriatic has been identified as a key location for wind farm developments. The development of offshore wind farms in the area would provide hard substrata for the settlement of sessile benthos. Since the seafloor of the Adriatic is predominantly sedimentary this may alter the larval connectivity of benthic populations in the region. Here, we simulated the release of larvae from benthic populations along the coasts of the Adriatic Sea using coupled bio-physical models and investigated the effect of pelagic larval duration on dispersal. Our model simulations show that currents typically carry particles from east to west across the Adriatic, whereas particles released along western coasts tend to remain there with the Puglia coast of Italy acting as a sink for larvae from benthic populations. We identify areas of high connectivity, as well as areas that are much more isolated, and discuss how these results can be used to inform marine spatial planning and the licensing of offshore marine renewable energy developments.

Asymmetric connectivity of spawning aggregations of a commercially important marine fish using a multidisciplinary approach

Understanding patterns of larval dispersal is key in determining whether no-take marine reserves are self-sustaining, what will be protected inside reserves and where the benefits of reserves will be observed. We followed a multidisciplinary approach that merged detailed descriptions of fishing zones and spawning time at 17 sites distributed in the Midriff Island region of the Gulf of California with a biophysical oceanographic model that simulated larval transport at Pelagic Larval Duration (PLD) 14, 21 and 28 days for the most common and targeted predatory reef fish, (leopard grouper Mycteroperca rosacea). We tested the hypothesis that source–sink larval metapopulation dynamics describing the direction and frequency of larval dispersal according to an oceanographic model can help to explain empirical genetic data. We described modeled metapopulation dynamics using graph theory and employed empirical sequence data from a subset of 11 sites at two mitochondrial genes to verify the model predictions based on patterns of genetic diversity within sites and genetic structure between sites. We employed a population graph describing a network of genetic relationships among sites and contrasted it against modeled networks. While our results failed to explain genetic diversity within sites, they confirmed that ocean models summarized via graph and adjacency distances over modeled networks can explain seemingly chaotic patterns of genetic structure between sites. Empirical and modeled networks showed significant similarities in the clustering coefficients of each site and adjacency matrices between sites. Most of the connectivity patterns observed towards downstream sites (Sonora coast) were strictly asymmetric, while those between upstream sites (Baja and the Midriffs) were symmetric. The best-supported gene flow model and analyses of modularity of the modeled networks confirmed a pulse of larvae from the Baja Peninsula, across the Midriff Island region and towards the Sonoran coastline that acts like a larval sink, in agreement with the cyclonic gyre (anti-clockwise) present at the peak of spawning (May–June). Our approach provided a mechanistic explanation of the location of fishing zones: most of the largest areas where fishing takes place seem to be sustained simultaneously by high levels of local retention, contribution of larvae from upstream sites and oceanographic patterns that concentrate larval density from all over the region. The general asymmetry in marine connectivity observed highlights that benefits from reserves are biased towards particular directions, that no-take areas need to be located upstream of targeted fishing zones, and that some fishing localities might not directly benefit from avoiding fishing within reserves located adjacent to their communities. We discuss the implications of marine connectivity for the current network of marine protected areas and no-take zones, and identify ways of improving it.