Delimitation of cryptic species drastically reduces the geographical ranges of marine interstitial ghost-worms (Stygocapitella; Annelida, Sedentaria).

The recognition of cryptic species concealed in traditionally established species may reveal new biogeographical patterns and alter the understanding of how biodiversity is geographically distributed. This is particularly relevant for marine ecosystems where the incidence of cryptic species is high and where species distribution data are often challenging to collect and interpret. Here, we studied specimens of the 'cosmopolitan' interstitial meiofaunal annelid Stygocapitella subterranea Knöllner, 1934 (Parergodrilidae, Orbiniida), obtaining data from four coastlines in the Northern hemisphere. Using phylogenetic tools and several species-delimitation methods (haplotype networks, GMYC, bPTP, maximum likelihood, posterior probability and morphology) we describe eight new Stygocapitella species. With one exception, all species are present along a single coastline, ultimately challenging the idea that Stygocapitella subterranea has a cosmopolitan distribution. We found evidence for several oceanic transitions having occurred in the past as well as a recent translocation, potentially due to human activity. No diagnostic characters were found, and qualitative and quantitative morphological data do not allow an unequivocal differentiation of the identified cryptic species. This suggests that (i) neither traditional diagnostic features nor quantitative morphology suffice to recognise species boundaries in cryptic species complexes, such as the Stygocapitella species complex; and that (ii) the recognition and description of cryptic species is of seminal importance for biodiversity assessments, biogeography and evolutionary biology.

Seascape genetics of the stalked kelp Pterygophora californica and comparative population genetics in the Santa Barbara Channel.

We conducted a population genetic analysis of the stalked kelp, Pterygophora californica, in the Santa Barbara Channel, California, USA. The results were compared with previous work on the genetic differentiation of giant kelp, Macrocystis pyrifera, in the same region. These two sympatric kelps not only share many life history and dispersal characteristics but also differ in that dislodged P. californica does not produce floating rafts with buoyant fertile sporophytes, commonly observed for M. pyrifera. We used a comparative population genetic approach with these two species to test the hypothesis that the ability to produce floating rafts increases the genetic connectivity among kelp patches in the Santa Barbara Channel. We quantified the association of habitat continuity and oceanographic distance with the genetic differentiation observed in stalked kelp, like previously conducted for giant kelp. We compared both overall (across all patches) and pairwise (between patches) genetic differentiation. We found that oceanographic transit time, habitat continuity, and geographic distance were all associated with genetic connectivity in P. californica, supporting similar previous findings for M. pyrifera. Controlling for differences in heterozygosity between kelp species using Jost's DEST , we showed that global differentiation and pairwise differentiation were similar among patches between the two kelp species, indicating that they have similar dispersal capabilities despite their differences in rafting ability. These results suggest that rafting sporophytes do not play a significant role in effective dispersal of M. pyrifera at ecologically relevant spatial and temporal scales.

Isolation by environment in the highly mobile olive ridley turtle (Lepidochelys olivacea) in the eastern Pacific.

Spatial and temporal scales at which processes modulate genetic diversity over the landscape are usually overlooked, impacting the design of conservation management practices for widely distributed species. We examine processes shaping population divergence in highly mobile species by re-assessing the case of panmixia in the iconic olive ridley turtle from the eastern Pacific. We implemented a biophysical model of connectivity and a seascape genetic analysis based on nuclear DNA variation of 634 samples collected from 27 nesting areas. Two genetically distinct populations largely isolated during reproductive migrations and mating were detected, each composed of multiple nesting sites linked by high connectivity. This pattern was strongly associated with a steep environmental gradient and also influenced by ocean currents. These findings relate to meso-scale features of a dynamic oceanographic interface in the eastern tropical Pacific (ETP) region, a scenario that possibly provides different cost-benefit solutions and selective pressures for sea turtles during both the mating and migration periods. We reject panmixia and propose a new paradigm for olive ridley turtles where reproductive isolation due to assortative mating is linked to its environment. Our study demonstrates the relevance of integrative approaches for assessing the role of environmental gradients and oceanographic currents as drivers of genetic differentiation in widely distributed marine species. This is relevant for the conservation management of species of highly mobile behaviour, and assists the planning and development of large-scale conservation strategies for the threatened olive ridley turtles in the ETP.

Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (Parastichopus californicus).

Marine populations are typically characterized by weak genetic differentiation due to the potential for long-distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST  = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2  = .49) than geographic distance (IBD) (R2  = .18), and directional processes played an important role in shaping fine-scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.

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.

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific.

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite-based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large-scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation-by-environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.

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.

Emergent patterns of population genetic structure for a coral reef community.

What shapes variation in genetic structure within a community of codistributed species is a central but difficult question for the field of population genetics. With a focus on the isolated coral reef ecosystem of the Hawaiian Archipelago, we assessed how life history traits influence population genetic structure for 35 reef animals. Despite the archipelago's stepping stone configuration, isolation by distance was the least common type of genetic structure, detected in four species. Regional structuring (i.e. division of sites into genetically and spatially distinct regions) was most common, detected in 20 species and nearly in all endemics and habitat specialists. Seven species displayed chaotic (spatially unordered) structuring, and all were nonendemic generalist species. Chaotic structure also associated with relatively high global FST. Pelagic larval duration (PLD) was not a strong predictor of variation in population structure (R2=0.22), but accounting for higher FST values of chaotic and invertebrate species, compared to regionally structured and fish species, doubled the power of PLD to explain variation in global FST (adjusted R2=0.50). Multivariate correlation of eight species traits to six genetic traits highlighted dispersal ability, taxonomy (i.e. fish vs. invertebrate) and habitat specialization as strongest influences on genetics, but otherwise left much variation in genetic traits unexplained. Considering that the study design controlled for many sampling and geographical factors, the extreme interspecific variation in spatial genetic patterns observed for Hawaìi marine species may be generated by demographic variability due to species-specific abundance and migration patterns and/or seascape and historical factors.

Genetic Structure of the Endangered Coral Cladocora caespitosa Matches the Main Bioregions of the Mediterranean Sea.

Population connectivity studies are a useful tool for species management and conservation planning, particular of highly threatened or endangered species. Here, we evaluated the genetic structure and connectivity pattern of the endangered coral Cladocora caespitosa across its entire distribution range in the Mediterranean Sea. Additionally, we examined the relative importance of sexual and asexual reproduction in the studied populations and their genetic diversity. A total of 541 individuals from 20 localities were sampled and analysed with 19 polymorphic microsatellite markers. Of the genotyped individuals, 482 (89%) had unique multilocus genotypes. Clonality percentages of the populations varied from 0% (in eight populations) to nearly 69% (in one population from Crete). A heterozygosity deficit and a high degree of inbreeding was the general trend in our data set. Population differentiation in C. caespitosa was characterised by significant pairwise F ST values with lower ones observed at an intraregional scale and higher ones, between populations from different biogeographic regions. Genetic structure analyses showed that the populations are divided according to the three main sub-basins of the Mediterranean Sea: the Western (Balearic, Ligurian and Tyrrhenian seas), the Central (Adriatic and Ionian seas) and the Eastern (Levantine and Aegean seas), coinciding with previously described gene flow barriers. However, the three easternmost populations were also clearly separated from one another, and a substructure was observed for the other studied areas. An isolation-by-distance pattern was found among, but not within, the three main population groups. This substructure is mediated mainly by dispersal along the coastline and some resistance to larval movement through the open sea. Despite the low dispersal ability and high self-recruitment rate of C. caespitosa, casual dispersive events between regions seem to be enough to maintain the species' considerable genetic diversity. Understanding the population connectivity and structure of this endangered scleractinian coral allows for more informed conservation decision making.

Fisheries genomics of snapper (Chrysophrys auratus) along the west Australian coast.

The efficacy of fisheries management strategies depends on stock assessment and management actions being carried out at appropriate spatial scales. This requires understanding of spatial and temporal population structure and connectivity, which is challenging in weakly structured and highly connected marine populations. We carried out a population genomics study of the heavily exploited snapper (Chrysophrys auratus) along ~2600 km of the Australian coastline, with a focus on Western Australia (WA). We used 10,903 filtered SNPs in 341 individuals from eight sampling locations to characterize population structure and connectivity in snapper across WA and to assess if current spatial scales of stock assessment and management agree with evidence from population genomics. Our dataset also enabled us to investigate temporal stability in population structure as well as connectivity between WA and its nearest, eastern jurisdictional neighbour. As expected for a species influenced by the extensive ocean boundary current in the region, low genetic differentiation and high connectivity were uncovered across WA. However, we did detect strong isolation by distance and genetic discontinuities in the mid-west and south-east. The discontinuities correlate with boundaries between biogeographic regions, influenced by on-shelf oceanography, and the sites of important spawning aggregations. We also detected temporal instability in genetic structure at one of our sites, possibly due to interannual variability in recruitment in adjacent regions. Our results partly contrast with the current spatial management of snapper in WA, indicating the likely benefits of a review. This study supports the value of population genomic surveys in informing the management of weakly structured and wide-ranging marine fishery resources.

Resolving the origins of invertebrate colonists in the Yangtze River Estuary with molecular markers: Implications for ecological connectivity.

Understanding connectivity over different spatial and temporal scales is fundamental for managing of ecological systems. However, controversy exists for wintertime ecological connectivity between the Yangtze River Estuary (YRE) and inner southwestern Yellow Sea. Here, we investigated ecological connectivity between the YRE and inner southwestern Yellow Sea in wintertime by precisely pinpointing the source of the newly colonized populations of a winter-spawning rocky intertidal invertebrate, Littorina brevicula (Philippi, 1844), on artificial structures along the coast of the Yangtze River Delta (YRD) using mitochondrial ND6 sequences and microsatellite data. Clear phylogeographic and genetic differentiation were detected between natural rocky populations south and north of the YRE, which resulted from the lack of hard substrate for rocky invertebrates in the large YRD coast. For the newly colonized populations on the coast of YRD, most individuals (98%) to the south of ~33.5°N were from natural rocky populations to the south of the YRE and most of those (94%) to the north of ~33.5°N were from the northern natural rocky populations, which demonstrated strong ecological connectivity between the inner southwestern Yellow Sea and the YRE in winter time. We presented the first genetic evidence that demonstrated a northward wintertime coastal current in the inner southwestern Yellow Sea, and precisely illustrated the boundary of the coastal current recently proposed by numerical experiment. These results indicated that the YRE serves as an important source of materials and energy for the inner southwestern Yellow Sea in winter, which can be crucial for the function of the Yellow Sea ecosystem.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea.

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

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.

Large-scale connectivity, cryptic population structure, and relatedness in Eastern Pacific Olive ridley sea turtles (Lepidochelys olivacea).

Endangered species are grouped into genetically discrete populations to direct conservation efforts. Mitochondrial control region (mtCR) haplotypes are used to elucidate deep divergences between populations, as compared to nuclear microsatellites that can detect recent structuring. When prior populations are unknown, it is useful to subject microsatellite data to clustering and/or ordination population inference. Olive ridley sea turtles (Lepidochelys olivacea) are the most abundant sea turtle, yet few studies have characterized olive ridley population structure. Recently, clustering results of olive ridleys in the Eastern Tropical Pacific Ocean suggested weak structuring (F ST = 0.02) between Mexico and Central America. We analyzed mtCR haplotypes, new microsatellite genotypes from Costa Rica, and preexisting microsatellite genotypes from olive ridleys across the Eastern Tropical Pacific, to further explore population structuring in this region. We subjected inferred populations to multiple analyses to explore the mechanisms behind their structuring. We found 10 mtCR haplotypes from 60 turtles nesting at three sites in Costa Rica, but did not detect divergence between Costa Rican sites, or between Central America and Mexico. In Costa Rica, clustering suggested one population with no structuring, but ordination suggested four cryptic clusters with moderate structuring (F ST = 0.08, p < .001). Across the Eastern Tropical Pacific, ordination suggested nine cryptic clusters with moderate structuring (F ST = 0.103, p < .001) that largely corresponded to Mexican and Central American populations. All ordination clusters displayed significant internal relatedness relative to global relatedness (p < .001) and contained numerous sibling pairs. This suggests that broadly dispersed family lineages have proliferated in Eastern Tropical Pacific olive ridleys and corroborates previous work showing basin-wide connectivity and shallow population structure in this region. The existence of broadly dispersed kin in Eastern Tropical Pacific olive ridleys has implications for management of olive ridleys in this region, and adds to our understanding of sea turtle ecology and life history, particularly in light of the natal-homing paradigm.

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.

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.