# High levels of connectivity over large distances in the diadematid sea urchin Centrostephanus sylviae

David Veliz, Noemi Rojas-Hernandez, Pablo Fibla, Boris Dewitte, Sebastian

Cornejo-Guzman, C. Parada

► To cite this version:

David Veliz, Noemi Rojas-Hernandez, Pablo Fibla, Boris Dewitte, Sebastian Cornejo-Guzman, et al. High levels of connectivity over large distances in the diadematid sea urchin Centrostephanus sylviae. PLoS ONE, 2021, 16 (11), pp.e0259595. 10.1371/journal.pone.0259595. hal-04739153

HAL Id: hal-04739153

https://cnrs.hal.science/hal-04739153v1

Submitted on 16 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# High levels of connectivity over large distances in the diadematid sea urchin *Centrostephanus sylviae*

David Veliz<sup>1,2,3*</sup>, Noemi Rojas-Hernández<sup>1</sup>, Pablo Fibla<sup>1</sup>, Boris Dewitte<sup>2,4,5</sup>, Sebastián Cornejo-Guzmán<sup>2,6</sup>, Carolina Parada<sup>2,6</sup>

1. Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile,

2. Núcleo Milenio de Ecología y Manejo Sustentable (ESMOI), Coquimbo, Chile, 3. Instituto de Ecología y Biodiversidad (IEB), Santiago, Chile, 4. Centro de Estudios Avanzados en Zonas Áridas (CEAZA), Coquimbo, Chile, 5. Université de Toulouse, CERFACS/CNRS, Toulouse, France, 6. Departamento de Geofísica, Universidad de Concepción, Concepción, Chile

* dveliz@uchile.cl

<figure><img src="image_2.png" /></figure>

## Abstract

Most benthic marine invertebrates with sedentary benthic adult phases have planktonic larvae that permit connectivity between geographically isolated populations. Planktonic larval duration and oceanographic processes are vital to connecting populations of species inhabiting remote and distant islands. In the present study, we analyzed the population genetic structure of the sea urchin *Centrostephanus sylviae*, which inhabits only the Juan Fernández Archipelago and the Desventuradas islands, separated by more than 800 km. For 92 individuals collected from Robinson Crusoe and Selkirk Islands (Juan Fernández Archipelago) and San Ambrosio Island (Desventuradas Islands), 7,067 single nucleotide polymorphisms (SNPs) were obtained. The results did not show a spatial genetic structure for *C. sylviae*; relative high migration rates were revealed between the islands. An analysis of the water circulation pattern in the area described a predominant northward water flow with periods of inverted flow, suggesting that larvae could move in both directions. Overall, this evidence suggests that *C. sylviae* comprises a single large population composed of individuals separated by more than 800 km.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Veliz D, Rojas-Hernández N, Fibla P, Dewitte B, Cornejo-Guzmán S, Parada C (2021) High levels of connectivity over large distances in the diadematid sea urchin *Centrostephanus sylviae*. PLoS ONE 16(11): e0259595. https://doi.org/10.1371/journal.pone.0259595

**Editor:** Tzen-Yuh Chiang, National Cheng Kung University, TAIWAN

**Received:** June 28, 2021

**Accepted:** October 21, 2021

**Published:** November 4, 2021

**Copyright:** © 2021 Veliz et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The datasets presented in this study can be found in: https://doi.org/10.34691/FK2/KWJDDA.

**Funding:** Dv, CP, BD: grant ANID FONDECYT 1191606. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. DV: FONDECYT 1200589 and Grant ANID ACE210006. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. BD: CLAP project, ANID (R20F0008)

## Introduction

Knowledge about the genetic structure of populations and gene flow between them is a central issue in population ecology. This information allows us to determine the importance of each population in a metapopulation context and the symmetry of reciprocal migration; key information for species management and conservation plans [1].

Gene flow among populations of benthic marine organisms occurs mainly by the interpopulation movement of larvae [2]. The life cycle of most of these organisms includes a low-mobility or sedentary benthic adult phase and a planktonic larval phase [3], the latter of vital importance in the gene flow in fragmented populations [4]. Larval development duration, larval behavior, and oceanographic processes have been described as the main factors affecting larval dispersal of these benthic organisms [5, 6]. It is known that species with long planktonic
and FutureMares European Union's Horizon 2020 research and innovation programme (grant agreement No 869300). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

larvae stages are capable of the largest potential connectivity distance, thereby promoting low population genetic differentiation [7, 8]. Furthermore, oceanographic processes are important in determining the direction of larval movement, in some cases promoting asymmetric gene flow between populations [9–11].

Genetic population variability has been used to study the genetic structure and gene flow among geographically isolated populations. Specifically, the use of thousands of single nucleotide polymorphisms (SNPs) obtained via next-generation sequencing represents a cost-effective alternative to obtain genetic data in non-model organisms [12].

The fragmented nature of the island systems, the span of the larval period and the type and scale of predominant oceanographic processes are all relevant in gene flow of the species inhabiting these islands; larval displacement is fundamental in determining the geographic extension of these populations [13, 14]. Two island systems of Chile separated by 800 km share most of their marine species [15]: the Juan Fernández Archipelago, encompassed by Robinson Crusoe-Santa Clara and Alejandro Selkirk islands, and the Desventuradas Islands, comprised of San Félix and San Ambrosio islands. Due to their high degree of geographic isolation from other areas in the Pacific Ocean, both island systems have a high degree of endemism in both terrestrial and marine biota. The Humboldt Current creates a strong biogeographic barrier between these island systems and the mainland [15]. Porobić et al. (2012) [16] described both mesoscale eddies (connectivity between the islands of the Juan Fernández Archipelago) and meanders (northward flow from the Juan Fernández Archipelago to the Desventuradas Islands) are important oceanographic features in this area, providing physical connectivity and establishing significant particle flow among the different island groups. Previous studies of population connectivity in the Juan Fernández Archipelago and the Desventuradas Islands have only been carried out for the lobster *Jasus frontalis*, which has a larval development of one year [17]. Genetic data and biophysical modeling revealed high connectivity between these island systems for this lobster species [16, 18], suggesting that both systems exhibit high gene flow for a species with a planktonic larval duration of one year.

The sea urchin *Centrostephanus sylviae* Fell 1975 [19] is a species endemic to the Juan Fernández Archipelago and the Desventuradas Islands. There is little published knowledge about the biology and ecology of this species. Friedlander et al. (2016) [15] found higher abundances in the Desventuradas Islands than the Juan Fernández Archipelago. Although the duration of the planktonic development of *C. sylviae* is unknown, a closely related species, *C. rodgersii*, has a larval development of four months [20]. It has been suggested that this period of larval development allows long distance larval dispersal for *C. rodgersii* between Australia and New Zealand [21].

The objective of this study was to estimate the degree of connectivity and population differentiation of the sea urchin *Centrostephanus sylviae* in the Juan Fernández Archipelago and the Desventuradas Islands. Using thousands of SNPs, population genetic structure and the degree of connectivity between populations (migrants between islands) of this species were described. Furthermore, a high-resolution reanalysis model analyzed the climatological currents between the islands to complement the genomic analysis. Considering the previous evidence for *C. rodgersii* in Australia and New Zealand, it was expected that the sea urchin *C. sylviae* would not show genetic differences between populations in the Juan Fernández Archipelago and the Desventuradas Islands.

# Materials and methods

## Sample collection

A total of 92 specimens of the sea urchin *Centrostephanus sylviae* were collected from three oceanic islands (Fig 1), 61 of which were collected from the Juan Fernández Archipelago. Of
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Sampling sites of *Centrostephanus sylviae* in the Juan Fernández Archipelago and Desventuradas Island. Red points represent the sampling areas on each island: Robinson Crusoe (33°36'S, 78°48'W), Alejandro Selkirk (33°48'S, 80°40'W) and San Ambrosio (26°20'S, 79°53'W). Maps drawn using library maps [22] and mapdata [23] implemented in R software [24].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g001

these 61 samples, 29 were collected from Robinson Island (October 2019) and 32 from Selkirk Island (December 2019). In addition, 31 individuals were collected from San Ambrosio Island of the Desventuradas Islands in January 2020. Samples were obtained by scuba diving in shallow water. The Aristotle's lantern of each specimen was preserved in 95% alcohol for genetic analysis. Sampling was performed under Fondecyt grant 1191606 and the Monitoring Program for Crustacean Fisheries of the Juan Fernández Archipelago (IFOP grant 2019) and sampling permission was granted by the Port Captain at Juan Fernández Archipelago, Chile.

## Genetic analysis

A small piece of muscle collected from each Aristotle's lantern was placed in 96-well PCR plates with 80 µL of 99% ethanol for genetic analysis. DNA extraction and massive sequencing were performed at Dart Diversity Arrays Technology Pty Ltd., Canberra, Australia (www.diversityarrays.com). DNA was digested with the restriction enzymes Sbfl and PstI following Kilian et al. (2012) [25]. Fragments > 200 bp were ligated with an 8 base pair barcode and amplified by PCR. The PCR products were standardized and sequenced on a HiSeq 2500 (Illumina Inc, San Diego, USA).
Quality control and initial SNP calling

Dart Diversity's bioinformatics service performed demultiplexing and removal of the DNA barcodes. More information about the detection of SNPs is described in Kilian et al. (2012). Raw SNP data from Dart were filtered using the dartR library [26] implemented in the R program [24]. This filter eliminated i) only one SNPs was retained in the reads containing two or more SNPs, ii) loci with a read depth below five or above 200, iii) loci with < 99% reproducibility, iv) monomorphic loci, v) loci with > 10% missing data, vi) individuals with > 5% missing data and vii) all SNPs with a minimum allele frequency (MAF) < 1%.

To avoid bias in estimating the differences between populations, the loci identified as under selection were eliminated. Four different approaches were used for this: i) the method based on likelihood implemented in the outflank function of the dartR library, ii) the method based on principal components, implemented in the PCAadapt library [27], also implemented in the R software, iii) the Bayesian method implemented in the BayeScan program [28], and iv) the method based on the relationship between FST and heterozygosity implemented in the Fsthet library [29] of the R software. In addition, all loci that showed evidence of selection in at least 2/4 methods implemented were removed. Finally, loci with significant departures from Hardy-Weinberg equilibrium in all sampling sites were removed using the dartR library, and loci showing linkage disequilibrium > 0.5 in all sampling sites were filtered with PLINK 2.0 software [30].

Genetic diversity

Expected heterozygosity (Hexp.), observed heterozygosity (Hobs.), and the inbreeding coefficient (FIS) were estimated using GENETIX v 4.05 [31]. Allelic richness (AR) was estimated with the divBasic function of the diveRsity library [32].

Population genetic structure

Population genetic structure was estimated using the SNP database for the three sampling sites of *C. sylviae*. Three methods were used for this: i) A principal coordinates analysis (PCoA) was used to observe the distribution of individuals in a multivariate space qualitatively, using the dartR library implemented in the software R; ii) the paired FST index was assessed with the gl. fst.pop function implemented in the dartR library; iii) the most probable number of genetic clusters (K) was estimated using the Bayesian approach implemented in the Structure software [33]. The admixture model and correlation of allele frequencies were used as input. The procedure was performed three times for each K (from K = 1 to K = 4) with a burn-in of 100,000 iterations and an after-burn-in of 200,000 iterations. The most likely number of clusters was estimated using the output in the Structure Harvester program [34]. The probability of each K was estimated using the method described in the Structure manual [35].

Finally, the direction and magnitude of the gene flow of *C. sylviae* between the islands were estimated using the divMigrate function of the diveRsity library [32] implemented in the R software. The effective number of migrants (Nm) was used as a distance measurement, and the asymmetry of the flow between the islands was tested using a bootstrap of 1000 iterations.

Oceanographic circulation in the area

GLOBAGLOBAL_REANALYSIS PHY_001_030 (CGLOPHY030) and
CMEMSGLOBALengerSIS FORECAST PHY_001_024 (CGLOPHY024) products dis-
tributed by CMEMS (http://marine.copernicus.eu/) were used to understand circulation pat-
terns between the Juan Fernandez Archipelago and the Desventuradas Islands. Briefly,
CGLOPHY030 is a daily global reanalysis model [36] providing resolution and depths, with output spanning between 1993 to 2018, and GLOPHY024 is a daily global forecast model [37] of a high-resolution regular grid 1/12° (~8 km) and 50 depth levels (from surface to 5500 m). The model output spanned from January 2016 to the present. It is based on the hydrodynamic model of the ocean Nucleus for European Modelling of the Ocean (NEMO) [38] that assimilates altimetry, sea surface temperature satellite data, vertical temperature, and salinity. NEMO v3.1 was used in our study.

Both model outputs were used to characterize the meridional and longitudinal velocities and transport patterns between the Juan Fernandez Archipelago and the Desventuradas Islands from 1993 to 2020. One-hundred meters of averaged U and V velocity components were extracted in the region between 81°–78°W and 35°–25.5°S that includes both island groups. Seasonal spatial climatology of the meridional and longitudinal velocities was also estimated. A time series of meridional and longitudinal velocities and transport between regions, averaged for the first 100 m and between 1000 and 100 m, was estimated for the same period (1993–2020).

## Results

In the 92 individuals analyzed, 85,933 raw SNPs were obtained. After applying filters and removing outlier loci, 7,067 SNPs and 89 individuals were retained: 28 from Robinson, 31 from Selkirk, and 30 from San Ambrosio. The allelic richness and heterozygosity of these filtered data showed similar values among the three islands studied (Table 1).

Population genetic structure analysis did not detect significant differences among the studied sites. The PCoA showed an overlap of the individuals in the multivariate space, suggesting no segregation of the populations (Fig 2). Significant differences in FST were not detected between the pairs of sites compared (Table 2). The Bayesian analysis implemented in the Structure software showed that K = 1 had the lowest value of ln (K) (mean LNP [K] = -365659.3, p = 0.999). However, Evanno's method indicated K = 3 as the most likely number of clusters. The graphs from K = 1 to K = 4 (Fig 3) suggest no clear genetic separation between the islands; thus, K = 1 is the most biologically viable explanation. Finally, the analysis performed with Div-Migrate suggests strong gene flow between the islands (Fig 4); the bootstrap did not detect significant asymmetric gene flow between the islands (p > 0.05).

No seasonal changes were observed at the meridional velocities, with a predominance of northward velocities between the region of the Juan Fernandez Archipelago and the Desventuradas Islands (Fig 5). Southward velocities were not significant. Within the Juan Fernandez Archipelago, an eastward flow (light blue) between Selkirk and Robinson islands characterized the velocities during all seasons (see the flow inside the black squares) except for fall, when the flow associated to the southward region of the islands reversed to westward velocities (oranges) connecting Robinson towards Selkirk islands (Fig 6).

Monthly meridional velocity variability time series showed a northward predominance between Juan Fernandez and Desventuradas in the studied years (1993 to 2020) in the surface depth (mean 0–100 m) and the deep strata (100–1000 m), with brief periods of inverted flows

<table><caption>Table 1. Summary of data used in the analysis of the sea urchin *C. sylviae* including sample size (n), allelic richness (AR) observed heterozygosity (Hobs.), expected heterozygosity (Hexp.), non-biased heterozygosity (Hn.b.) and FIS at each study site.</caption><thead><tr><th>Site</th><th>n</th><th>AR</th><th>Hobs.</th><th>Hexp.</th><th>Hn.b.</th><th>FIS</th></tr></thead><tbody><tr><td>Robinson</td><td>28</td><td>1.71</td><td>0.1095</td><td>0.1466</td><td>0.1494</td><td>0.2709</td></tr><tr><td>Selkirk</td><td>31</td><td>1.72</td><td>0.1107</td><td>0.1470</td><td>0.1495</td><td>0.26288</td></tr><tr><td>San Ambrosio</td><td>30</td><td>1.7</td><td>0.1024</td><td>0.1431</td><td>0.1456</td><td>0.30043</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0259595.t001
<figure><img src="image_5.png" /><figcaption>Fig 2. Principal coordinate analysis (PCoA) performed with sea urchin *C. sylviae* data obtained from Robinson Island (ROB), San Ambrosio Island (SAM), and Selkirk Island (SK).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g002

(southward). Monthly longitudinal velocity time series showed eastward and westward flow connecting Selkirk and Robinson islands (Fig 7), with a similar current pattern in the surface and deep strata.

Discussion

As expected, our results suggest that the sea urchin *Centrostephanus sylviae* inhabiting the Desventuradas Islands and the Juan Fernández Archipelago are a single panmictic genetic population with strong gene flow between these two systems, even though they are separated by more than 800 km. Since adult echinoids have low mobility, migration and gene flow occur mainly via planktonic larvae.

Previous evidence showed high dispersal of larvae in sea urchins, which allows their populations to cover large geographic areas. For example, *C. rodgersii* shows weak genetic differentiation and no isolation over thousands of kilometers between eastern Australia and northern New Zealand [21]. Other sea urchins have similarly shown no large-scale genetic differences: *Strongylocentrotus droebachiensis* inhabiting Atlantic Canada [39], *Paracentrotus lividus* in the Mediterranean [40], *Tripneustes gratilla* in the Philippines [41], and *S. franciscanus* in California [42], among others.

<table><caption>Table 2. Pairwise FST (above diagonal) and P-value (below diagonal) values for individuals collected in the islands studied.</caption><thead><tr><th></th><th>Robinson</th><th>Selkirk</th><th>San Ambrosio</th></tr></thead><tbody><tr><th>Robinson</th><td></td><td></td><td></td></tr><tr><th>Selkirk</th><td>0.609</td><td>-0.00011</td><td>0.00056</td></tr><tr><th>San Ambrosio</th><td>0.139</td><td>0.476</td><td>0.00003</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0259595.t002
Other benthic invertebrate taxa (e.g., marine crabs) with similar planktonic larval duration have large geographic populations, suggesting high dispersal of the larvae. For example, *Cancer*

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Population structure of the sea urchin *C. sylviae* inferred using the software STRUCTURE for K = 1 to K = 4 of 89 individuals from the three studied islands. A vertical line represents each individual, and each color represents the probability of belonging to one of the genetic clusters. ROB = Robinson Island, SAM = San Ambrosio Island, and SK = Selkirk Island.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g003

*pagurus* have a planktonic larval duration of 3 months [43] and did not show genetic differ-
ences among populations sampled in the Swedish Sea [44]. Moreover, Rojas-Hernández et al.
(2016) [45] did not detect differences along 700 km of coastline in Chile in the commercial
crab *Metacarcinus edwardsii*, which has a planktonic development of 60 days (Quintana, 1989)
[46]. In general, the evidence of sea urchins and crabs suggests that a planktonic development
of at least 60 days allows high dispersion between geographically remote locations.

In the same studied islands, Porobic et al. (2013) [18] described the population genetic
structure of the lobster *Jasus frontalis*. In this study, the individuals sampled in the Juan Fer-
nández Archipelago and Desventuradas Islands systems showed a single panmictic population
and high gene flow between islands, revealed by COI gene variability. This evidence is consis-
tent with the pattern detected in the current study. Furthermore, Porobic et al. (2012) [16]
used a biophysical model to study the connectivity of *J. frontalis* populations between the Juan
Fernández Archipelago and the Desventuradas Islands. The authors detected an asymmetric
connectivity pattern between these islands, consisting of significantly greater larvae dispersal
from Juan Fernández to the Desventuradas Islands. Our analysis using Nemo corroborates the
analysis performed by Porobic et al. (2012) [16] that used the hydrodynamic model OFES
(Ocean Model for the Earth Simulator) described by Masumoto et al. (2004) [47]. The velocity
field and time-series analyses in the present study support this asymmetric relationship in the
circulation pattern, revealing a predominance of northward velocities connecting the Juan Fer-
nandez system to the Desventuradas Islands. The genetic analysis of *C. sylviae* did not show an
asymmetric gene flow pattern as expected based on the biophysical model, suggesting that
other factors could be involved in promoting the high gene flow observed in the current study.
Relative migration network

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Migration network for the sampling sites of *C. sylviae* on Robinson Island (ROB), Selkirk Island (SK), and Sam Ambrosio Island (SAM). This figure was generated with the divMigrate function implemented in diveRsity (Keenan et al., 2013 [32]). Each vertex represents a sampling site, and each arrow represents the magnitude and direction of the relative migration.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g004

Therefore, it is necessary to continue studying dispersal patterns in this area to understand how the larval behavior, circulation (both mean and turbulent flow) and the biogeochemical environment relate to species connectivity between these two island systems, since it is known that the crustacean (e.g., [48]) and sea urchin larvae [49] have larval behavior of vertical migration. It was described that mean eddy activity corresponding to transient small scale (~10 to 100 km) recirculation patterns (eddies) is at a relative maximum between the two island systems, which could selectively affect the dispersion of species with planktonic larvae. Furthermore, Desventuradas and Juan Fernandez are in the neighborhood of an intense oxygen minimum zone [50], so that sharp oxygen gradients could also influence the species in a distinctive way.

Finally, *C. sylvae* larvae could move from these islands separated by more than 800 km and continue to connect these geographic groups present on the different islands studied. This evidence is crucial to understand how species are connected in island systems that have marine protected areas. In 2015, Chile created the Nazca-Desventuradas Marine Park of approximately 300,000 km² [51] which includes the Desventuradas Islands (San Ambrosio and San Felix islands), with the main goal being the conservation of the biological processes in the area.
<figure><img src="image_8.png" /><figcaption>Fig 5. Seasonal (JFM, summer; AMJ, fall; JAS, winter; OND, spring) meridional (positive, northward; negative, southward) velocities averaged over the first 100 m between Juan Fernandez Archipelago and the Desventuradas Islands, 25.5–35°S and 81.5–78°W. Data was obtained with Copernicus model products freely distributed by CMEMS (http://marine.copernicus.eu/).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g005

<figure><img src="image_9.png" /><figcaption>Fig 6. Seasonal (JFM, summer; AMJ, fall; JAS, winter; OND, spring) longitudinal (positive, westward; negative, eastward) velocities averaged over the first 100 m in the Juan Fernandez Archipelago and the Desventuradas Islands, 25.5–35°S and 81.5–78°W. Data was obtained with Copernicus model products freely distributed by CMEMS (http://marine.copernicus.eu/).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g006
<figure><img src="image_10.png" /><figcaption>**Fig 7.** a) Monthly meridional mean velocity (blue line, mean 0–100 m; red line, mean 100–1000 m) time series for the region 25.5–35°S and 81.5–78°W (positive, northward; negative, southward) explaining connectivity between Juan Fernandez Archilepago (JFA) and Desventuradas Islands (DA); b) Monthly longitudinal mean velocity (blue line, mean 0–100 m; red line, mean 100–1000 m) time series (positive, eastward; negative, westward) explaining connectivity between Robinson and Selkirk islands. Data was obtained with Copernicus model products freely distributed by CMEMS (http://marine.copernicus.eu/).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0259595.g007

It is important to note that the Desventuradas Islands and Robinson Crusoe Island (Juan Fernandez Archipelago) are uninhabited and inhabited by humans, respectively. Thus, high connectivity suggests that species exploited on the inhabited island (Robinson Crusoe Island) could be currently subsidized from this Marine Park. New investigations of species with different planktonic larval durations may help determine the level of connectivity of the benthic marine communities shared between these two island systems and, incidentally, may elucidate the role of the marine protected area in the face of exploitation pressure in the area.

**Acknowledgments**

We would like to thank Pablo Manriquez and the fishermen from Juan Fernández for sampling in Selkirk, Robinson and Sam Ambrosio islands.

**Author Contributions**

**Conceptualization:** David Veliz, Boris Dewitte, Carolina Parada.

**Data curation:** David Veliz, Noemi Rojas-Hernández, Boris Dewitte, Carolina Parada.
**Formal analysis:** David Veliz, Noemi Rojas-Hernández, Pablo Fibla, Sebastián Cornejo-Guzmán, Carolina Parada.

**Funding acquisition:** David Veliz, Boris Dewitte, Carolina Parada.

**Investigation:** David Veliz, Noemi Rojas-Hernández, Boris Dewitte, Carolina Parada.

**Methodology:** David Veliz, Noemi Rojas-Hernández, Pablo Fibla, Boris Dewitte, Sebastián Cornejo-Guzmán, Carolina Parada.

**Project administration:** David Veliz, Noemi Rojas-Hernández, Boris Dewitte, Sebastián Cornejo-Guzmán, Carolina Parada.

**Resources:** David Veliz, Carolina Parada.

**Software:** Sebastián Cornejo-Guzmán, Carolina Parada.

**Supervision:** David Veliz, Noemi Rojas-Hernández, Carolina Parada.

**Validation:** David Veliz, Noemi Rojas-Hernández, Sebastián Cornejo-Guzmán, Carolina Parada.

**Visualization:** David Veliz, Noemi Rojas-Hernández, Sebastián Cornejo-Guzmán, Carolina Parada.

**Writing – original draft:** David Veliz, Noemi Rojas-Hernández, Pablo Fibla, Boris Dewitte, Sebastián Cornejo-Guzmán, Carolina Parada.

**Writing – review & editing:** David Veliz, Noemi Rojas-Hernández, Pablo Fibla, Boris Dewitte, Sebastián Cornejo-Guzmán, Carolina Parada.

## References

1. Segelbacher G, Cushman SA, Epperson BK, Fortin M-J, Francois O, Hardy OJ, et al. Applications of landscape genetics in conservation biology: concepts and challenges. Conserv Genet. 2010; 11: 375–385.

2. Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Annu Rev Mar Sci. 2009; 1:443–466. https://doi.org/10.1146/annurev.marine.010908.163757 PMID: 21141044

3. Thorson G. Reproductive and larval ecology of marine bottom invertebrates. Biol Rev. 1950; 25:1–45. https://doi.org/10.1111/j.1469-185x.1950.tb00585.x PMID: 24537188

4. Jones GP, Almany GR, Russ GR, Sale PF, Steneck RS, van Oppen MJH, et al. Larval retention and connectivity among populations of corals and reef fishes: history, advances and challenges. Coral Reefs. 2009; 28:307–325.

5. Jenkins SR, Hawkins SJ. Barnacle larval supply to sheltered rocky shores: a limiting factor? Hydrobiologia. 2003; 503:143–151.

6. Pineda J, Hare JA, Sponaugle S. Consequences for population connectivity. Oceanography. 2007; 20:22–39.

7. Doherty PJ, Planes S, Mather P. Gene flow and larval duration in seven species of fish from the great barrier reef. Ecology. 1995; 76:2373–2391.

8. Tremi EA, Roberts JJ, Chao Y, Halpin PN, Possingham HP, Riginos C. Reproductive output and duration of the pelagic larval stage determine seascape-wide connectivity of marine populations. Integr Comp Biol. 2012; 52: 525–537. https://doi.org/10.1093/icb/ics101 PMID: 22821585

9. Xuereb A, Benestan L, Normandéau E, Daigle RM, Curtis JMR, Bernatchez L, et al. Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (*Parastichopus californicus*). Mol Ecol. 2018; 27: 2347–2364. https://doi.org/10.1111/mec.14589 PMID: 29654703

10. Meerhoff E, Yannicelli B, Dewitte B, Díaz-Cabrera E, Vega-Retter C, Ramos M, et al. Asymmetric connectivity of the lobster *Panulirus pascuensis* in remote islands of the Southern Pacific: importance for its management and conservation. Bulletin of Marine Science. 2018; 94(3):753–774.
11. Riginos C, Hock K, Matias AM, Mumby PJ, van Oppen MJH, Lukoschek V. Asymmetric dispersal is a critical element of concordance between biophysical dispersal models and spatial genetic structure in Great Barrier Reef corals. Diversity Distrib. 2019; 25:1684–1696.

12. Helyar SJ, Hemmer-Hansen J, Bekkevold D, Taylor MI, Ogden R, Limborg MT, et al. Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges. Mol Ecol Resour. 2011; 11:123–136. https://doi.org/10.1111/j.1755-0998.2010.02943.x PMID: 21429169

13. Hellberg ME. Dependence of gene flow on geographic distance in two solitary corals with different larval dispersal capabilities. Evolution. 1996; 50:1167–1175. https://doi.org/10.1111/j.1558-5646.1996.tb02357.x PMID: 28565289

14. Bell JJ. Connectivity between island Marine Protected Areas and the mainland. Biol Conserv. 2008; 141:2807–2820.

15. Friedlander AM, Ballesteros E, Caselle JE, Gaymer CF, Palma AT, Petit I, et al. Marine Biodiversity in Juan Fernández and Desventuradas Islands, Chile: Global Endemism Hotspots. PLoS One. 2016; 11: e0145059. https://doi.org/10.1371/journal.pone.0145059 PMID: 26734732

16. Porobić J, Parada C, Ernst B, Hormazábal SE, Combes V. Modelación de la conectividad de las subpopaciones de la langosta de Juan Fernández (*Jasus frontalis*), a través de un modelo biofísico. Lat Am J Aquat Res. 2012; 40:613–632.

17. Arana P. Perspectivas históricas y proyecciones de la actividad pesquera realizada en el archipiélago de Juan Fernández, Chile. In: Islas Oceánicas Chilenas: Conocimiento Científico y Necesidades de Investigaciones (Castilla JC ed.). Ediciones Universidad Católica de Chile. 1987; pp. 319–353.

18. Porobić J, Canales-Aguirre CB, Ernst B, Galleguillos R, Hernández CE. Biogeography and historical demography of the Juan Fernández rock lobster, *Jasus frontalis* (Milne Edwards, 1837). J Hered. 2013; 104:223–233. https://doi.org/10.1093/jhered/ess141 PMID: 23341535

19. Fell FJ. The echinoid genus *Centrostephanus* in the South Pacific Ocean with a description of a new species. J R Soc N Z. 1975; 5: 179–193.

20. Huggett MJ, King CK, Williamson JE, Steinberg PD. Larval development and metamorphosis of the Australian diadematid sea urchin *Centrostephanus rodgersii*. Invertebr Reprod Dev. 2005; 47:197–204.

21. Banks SC, Piggott LMP, Williamson JE, Bové U, Holbrook NJ, Beheregaray LB. Oceanic variability and coastal topography shape genetic structure in a long-dispersing sea urchin. Ecology. 2007; 88: 3055–3064. https://doi.org/10.1890/07-0091.1 PMID: 18229840

22. Becker RA, Wilks AR. R version by Ray Brownrigg. Enhancements by TP Minka and A Deckmyn. maps: Draw Geographical Maps. R package version 3.3.0. 2018a. https://CRAN.R-project.org/package=maps.

23. Becker RA, Wilks AR. R version by Ray Brownrigg. mapdata: Extra Map Databases. 2018b. R package version 2.3.0. https://CRAN.R-project.org/package=mapdata.

24. R Core Team. R: A language and environment for statistical computing. 2021; R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.

25. Kilian A, Wenzl P, Huttner E, Carling J, Xia L, Blois H, et al. Diversity arrays technology: a generic genome profiling technology on open platforms. Methods Mol Biol. 2012; 888:67–89. https://doi.org/10.1007/978-1-61779-870-2_5 PMID: 22665276

26. Gruber B, Unmack PJ, Berry OF, Georges A. dartr: An R package to facilitate analysis of SNP data generated from reduced representation genome sequencing. Mol Ecol Resour. 2018; 18:691–699. https://doi.org/10.1111/1755-0998.12745 PMID: 29266847

27. Luu K, Bazin E, Blum MGB. pcadapt: an R package to perform genome scans for selection based on principal component analysis. Mol Ecol Resour. 2017; 17:67–77. https://doi.org/10.1111/1755-0998.12592 PMID: 27601374

28. Foll M, Gaggiotti OE. A genome scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics. 2008; 180:977–993. https://doi.org/10.1534/genesics.108.092221 PMID: 18780740

29. Flanagan SP, Jones AG. Constraints on the Fst-heterozygosity outlier approach. J Hered. 2017; 108:561–573. https://doi.org/10.1093/jhered/esx048 PMID: 28486592

30. Chang, CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015; 4. URL: www.cog-genomics.org/plink/2.0/.

31. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, Logiciel sous Windows pour la Genetique des Populations. Laboratoire Genome, populations, interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier, France. 1996. French.
32. Keenan K, McGinnity P, Cross TF, Crozier WW, Prodöhl PA. diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol Evol. 2013; 4:782-788.

33. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155:945-959. https://doi.org/10.1093/genetics/155.2.945 PMID: 10835412

34. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012; 4:359-361.

35. Pritchard JK, Wen X, Falush D. Documentation for structure software: Version 2.3. 2010; University of Oxford http://pritch.bsd.uchicago.edu/structure.html.

36. Fernandez E, Lellouche JM. Product User Manual: For the GLOBAL Ocean Sea Physical Reanalysis Product GLOBAL_REANALYSIS PHY 001_030. Report CMEMS-GLO-PUM, (001-030), 15. 2018; https://resources.marine.copernicus.eu/documents/PUM/CMEMS-GLO-PUM-001-030.pdf.

37. von Schuckmann K, Cazenave A, Chambers D, Hansen J, Josey S, Kosaka Y, et al. An imperative to monitor Earth's energy imbalance. Nat Clim Chang. 2016; 6: 138-144.

38. Gurvan M, Bourdallé-Badie R,nut J, Clementi E, Coward A, Ethé C, et al. Notes Du Pôle De Modélisat ion De L'institut Pierre-Simon Laplace (IPSL). Zenodo. 2019. http://doi.org/10.5281/zenodo.3878122.

39. Addison JA, Hart MW (2004). Analysis of population structure of a green sea urchin (*Strongylocentrotus droebanchiensis*) using microsatellites. Mar Biol. 2004; 144: 243-251.

40. Duran S, Palacín C, Becerro MA, Turon X, Giribet G. Genetic diversity and population structure of the commercially harvested sea urchin *Paracentrotus lividus* (Echinodermata, Echinoidea). Mol Ecol. 2004; 13:3317-3328. https://doi.org/10.1111/j.1365-294X.2004.02338.x PMID: 15487992

41. Casilagan ILN, Juinio-Meñez MA, Crandall ED. Genetic diversity, population structure, and demographic history of exploited sea urchin populations (*Tripneustes gratilla*) in the Philippines. J Exp Mar Biol Ecol. 2013; 449:284-293.

42. Benham CE, Supernault KJ, Burton RS. Genetic assessment of the population connectivity of the red urchin (*Strongylocentrotus franciscanus*). J Exp Mar Biol Ecol. 2012; 432-433: 47-54.

43. Weiss M, Thatje S, Heilmayer O, Anger K, Brey T, Keller M. Influence of temperature on the larval development of the edible crab, *Cancer pagurus*. J Mar Biol Assoc UK. 2009; 89: 753-759.

44. Ungfors A, McKeown NJ, Shaw PW, Andre C. Lack of spatial genetic variation in the edible crab (*Cancer pagurus*) in the Kattegat Skagerrak area. ICES J Mar Sci. 2009; 66: 462-469.

45. Rojas-Hernandez N, Veliz D, Riveros MP, Fuentes JP, Pardo LM. Highly connected populations and temporal stability in allelic frequencies of a harvested crab from the southern Pacific Coast. PlosOne. 2016; 11(11): e0166029. https://doi.org/10.1371/journal.pone.0166029 PMID: 27814382

46. Quintana R. Larval development of the edible crab, *Cancer edwardsi* Bell, 1835 under laboratory conditions (Decapoda, Brachyura). Rep Usa Mar Biol Inst Kochi Univ. 1983; 5:1-19.

47. Masumoto Y, Sasaki H, Kagimoto T, Komori N, Ishida A, Sasai Y, et al. A fifty-year eddy-resolving simulation of the world ocean-Preliminary outcomes of OFES (OGCM for the Earth Simulator). J Earth Simul. 2004: 1:35-56.

48. Marta-Almeida M, Dubert J, Peliz A, Queiroga H. Influence of vertical migration pattern on retention of crab larvae in a seasonal upwelling system. Mar Ecol Progr Ser. 2006; 307:1-19.

49. Pennington JT, Emlet RB. Ontogenetic and diel vertical migration of a planktonic echinoid larva, *Dendraster excentricus* (Eschscholtz): Occurrence, causes, and probable consequences. J Exp Mar Biol Ecol. 1986; 104:69-95.

50. Dewitte B, Conejero C, Ramos M, Bravo L, Garçon V, Parada C, et al. Understanding the impact of climate change in the oceanic circulation in the Chilean oceanic island ecoregions. Aquat Conserv. 2021; 31, 232-252.

51. White TD, Ong T, Ferretti F, Block BA, McCauley DJ, Micheli F, et al. Tracking the response of industrial fishing fleets to large marine protected areas in the Pacific Ocean. Conserv Biol. 2020; 34: 1571-1578. https://doi.org/10.1111/cobi.13584 PMID: 33031635