Complex patterns of genetic structure in the sea cucumber Holothuria (Metriatyla) scabra from the Philippines: implications for aquaculture and fishery management

The sandfish <i>Holothuria (Metriatyla) scabra</i>, is a high-value tropical sea cucumber harvested from wild stocks for over four centuries in multi-species fisheries across its Indo-Pacific distribution, for the global bêche-de-mer (BDM) trade. Within Southeast Asia, the Philippines is an important centre of the BDM trade, however overharvesting and largely open fishery management have resulted in declining catch volumes. Sandfish mariculture has been developed to supplement BDM supply and assist restocking efforts; however, it is heavily reliant on wild populations for broodstock supply. Consequently, to inform fishery, mariculture, germplasm and translocation management policies for both wild and captive resources, a high-resolution genomic audit of 16 wild sandfish populations was conducted, employing a proven genotyping-by-sequencing approach for this species (DArTseq). Genomic data (8,266 selectively-neutral and 117 putatively-adaptive SNPs) were used to assess fine-scale genetic structure, diversity, relatedness, population connectivity and local adaptation at both broad (biogeographic region) and local (within-biogeographic region) scales. An independent hydrodynamic particle dispersal model was also used to assess population connectivity. The overall pattern of population differentiation at the country level for <i>H. scabra</i> in the Philippines is complex, with nine genetic stocks and respective management units delineated across 5 biogeographic regions: (1) Celebes Sea, (2) North and (3) South Philippine Seas, (4) South China and Internal Seas and (5) Sulu Sea. Genetic connectivity is highest within proximate marine biogeographic regions (mean <i>F</i>_st=0.016), with greater separation evident between geographically distant sites (<i>F</i>_st range=0.041–0.045). Signatures of local adaptation were detected among six biogeographic regions, with genetic bottlenecks at 5 sites, particularly within historically heavily-exploited locations in the western and central Philippines. Genetic structure is influenced by geographic distance, larval dispersal capacity, species-specific larval development and settlement attributes, variable ocean current-mediated gene flow, source and sink location geography and habitat heterogeneity across the archipelago. Data reported here will inform accurate and sustainable fishery regulation, conservation of genetic diversity, direct broodstock sourcing for mariculture and guide restocking interventions across the Philippines.