Search

Abstract Morphologically cryptic and pseudo-cryptic species pose a challenge to taxonomic identification and assessments of species diversity and distributions. Such is the case for the sea cucumber Stichopus horrens, commonly confused with Stichopus monotuberculatus. Here, we used mitochondrial cytochrome oxidase subunit I (COI) and microsatellite markers to examine genetic diversity in Stichopus cf. horrens throughout the Philippine archipelago, to aid species identification and clarify species boundaries. Phylogenetic analysis reveals two recently diverged COI lineages (Clade A and Clade B; c. 1.35–2.54 Mya) corresponding to sequence records for specimens identified as S. monotuberculatus and S. horrens, respectively. Microsatellite markers reveal two significantly differentiated genotype clusters broadly concordant with COI lineages (Cluster 1, Cluster 2). A small proportion of individuals were identified as later-generation hybrids indicating limited contemporary gene flow between genotype clusters, thus confirming species boundaries. Morphological differences in papillae distribution and form are observed for the two species, however tack-like spicules from the dorsal papillae are not a reliable diagnostic character. An additional putative cryptic species was detected within Clade B-Cluster 2 specimens warranting further examination. We propose that these lineages revealed by COI and genotype data be referred to as Stichopus cf. horrens species complex.

The sandfish Holothuria (Metriatyla) scabra, 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 H. scabra 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 Fst=0.016), with greater separation evident between geographically distant sites (Fst 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.

Recurring outbreaks of crown-of-thorns starfish (COTS) severely damage healthy corals, especially in the Western Pacific Ocean. To obtain a better understanding of population genetics of COTS and historical colonization across the Pacific Ocean, complete mitochondrial genomes were sequenced from 243 individuals collected in 11 reef regions. Our results indicate that Pacific COTS (Acanthaster cf. solaris) comprise two major clades, an East-Central Pacific (ECP) clade and a Pan-Pacific (PP) clade, separation of which was supported by high bootstrap value. The ECP clade consists of COTS from French Polynesia, Fiji, Vanuatu and the Great Barrier Reef (GBR). The Hawaii population is unique within this clade, while California COTS are included in EPC clade. On the other hand, the PP clade comprises multiple lineages that contain COTS from Vietnam, the Philippines, Japan, Papua New Guinea, Micronesia, the Marshall Islands, GBR, Vanuatu, Fiji and French Polynesia. For example, a lineage of the PP clade, which has the largest geographic distribution, includes COTS from all of these locations. These results suggest two alternative histories of current geographic distributions of COTS in the Pacific Ocean, an ECP clade ancestry or Western Pacific clade ancestry. Although further questions remain to be explored, this discovery provides an evolutionary context for the interpretation of COTS population structure which will aid future coral reef research in the Pacific Ocean, and ultimately improve reef management of COTS.

Stock identification and delineation are important in the management and conservation of marine resources. These were highlighted as priority research areas for Bali sardinella (Sardinella lemuru) which is among the most commercially important fishery resources in the Philippines. Previous studies have already assessed the stocks of S. lemuru between Northern Mindanao Region (NMR) and Northern Zamboanga Peninsula (NZP), yielding conflicting results. Phenotypic variation suggests distinct stocks between the two regions, while mitochondrial DNA did not detect evidence of genetic differentiation for this high gene flow species. This paper tested the hypothesis of regional structuring using genome-wide single nucleotide polymorphisms (SNPs) acquired through restriction site-associated DNA sequencing (RADseq). We examined patterns of population genomic structure using a full panel of 3,573 loci, which was then partitioned into a neutral panel of 3,348 loci and an outlier panel of 31 loci. Similar inferences were obtained from the full and neutral panels, which were contrary to the inferences from the outlier panel. While the full and neutral panels suggested a panmictic population (global FST ∼ 0, p > 0.05), the outlier panel revealed genetic differentiation between the two regions (global FST = 0.161, p = 0.001; FCT = 0.263, p < 0.05). This indicated that while gene flow is apparent, selective forces due to environmental heterogeneity between the two regions play a role in maintaining adaptive variation. Annotation of the outlier loci returned five genes that were mostly involved in organismal development. Meanwhile, three unannotated loci had allele frequencies that correlated with sea surface temperature. Overall, our results provided support for local adaptation despite high levels of gene flow in S. lemuru. Management therefore should not only focus on demographic parameters (e.g., stock size and catch volume), but also consider the preservation of adaptive variation.