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.

Abstract 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

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.

Prolonged thermal stress and high levels of solar irradiance can disrupt the coral-algal symbiosis and cause bleaching and lowered overall fitness that lead to the likely death of the cnidarian host. Adaptive bleaching and acclimatization of corals, which posits bleaching as an opportunity for the coral host to switch its currently susceptible endosymbionts to more stress-tolerant taxa, offers hope for survival of reefs amid rapidly warming oceans. In this study, we explored the diversity and distribution of coral-zooxanthellae associations in the context of geospatial patterns of sea surface temperature (SST) and thermal anomalies across the Philippine archipelago. Thermal clusters based on annual sea surface temperature means and each site’s frequency of exposure to heat stress were described using three-decade (1985–2018) remotely sensed data. Haphazard sampling of 628 coral fragments was conducted in 14 reef sites over 3 years (2015–2018). Using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) fingerprinting and sequencing of the zooxanthellae ITS2 region, we characterized endosymbiont diversity within four reef-building coral families across archipelagic thermal regimes. Consistency in dominant Symbiodiniaceae taxon was observed in Acropora spp., Porites spp., and Heliopora coerulea. In contrast, the family Pocilloporidae (Pocillopora spp., Seriatopora spp., and Stylophora pistillata) exhibited biogeographic variability in zooxanthellae composition, concordant with inferred occurrences of sustained thermal stress. Multivariate analyses identify two broad Pocilloporidae clusters that correspond with mean SST ranges and frequency of exposure to bleaching-level thermal stress which are largely supported by ANOSIM. Differences in zooxanthellae assemblages may reflect host-specific responses to ecological or environmental gradients across biogeographic regions. Such patterns of variability provide insight and support for the adaptability and potential resilience of coral communities in geographically and oceanographically complex regions, especially amidst the increasing severity of global and local-scale stressors.

Tropical coral reefs are threatened by local-scale stressors that are exacerbated by global ocean warming and acidification from the post-industrial increase of atmospheric CO2 levels. Despite their observed decline in the past four decades, little is known on how Philippine coral reefs will respond to ocean warming and acidification. This study explored individual and synergistic effects of present-day (pH 8.0, 28°C) and near-future (pH 7.7, 32°C) scenarios of ocean temperature and pH on the adult Favites colemani, a common massive reef-building coral in Bolinao-Anda, Philippines. Changes in seawater temperature drive the physiological responses of F. colemani, whereas changes in pH create an additive effect on survival, growth, and photosynthetic efficiency. Under near-future scenarios, F. colemani showed sustained photosynthetic competency despite the decline in growth rate and zooxanthellae density. F. colemani exhibited specificity with the Cladocopium clade C3u. This coral experienced lower growth rates but survived projected near-future ocean warming and acidification scenarios. Its pH-thermal stress threshold is possibly a consequence of acclimation and adaptation to local environmental conditions and past bleaching events. This research highlights the importance of examining the susceptibility and resilience of Philippine corals to climate-driven stressors for future conservation and restoration efforts in the changing ocean.

The taxonomy of the mottled rabbitfish Siganus fuscescens species complex has long been challenging. In this study, we analyzed microsatellite genotypes, mitochondrial lineages, and morphometric data from 373 S. fuscescens individuals sampled from the northern Philippines and Hong Kong (South China Sea, Philippine Sea and Sulu Sea basins), to examine putative species boundaries in samples comprising three co-occurring mitochondrial lineages previously reported to characterize S. fuscescens (Clade A and Clade B) or S. canaliculatus (Clade C). We report the existence of two cryptic species within S. fuscescens in the northeast region of the South China Sea and northern Philippine Sea, supported by genetic and morphological differences. Individual-based assignment methods recovered concordant groupings of individuals into two nuclear genotype clusters (Cluster 1, Cluster 2) with (1) limited gene flow, if any, between them (FST = 0.241; P < 0.001); (2) low frequency of later-generation hybrids; (3) significant association with mitochondrial Clade A and Clade B, respectively; and (4) subtle yet significant body shape differences as inferred from geometric morphometric analysis. The divergence between mitochondrial Clade C and the two other clades was not matched by genetic differences at microsatellite marker loci. The occurrence of discordant mitonuclear combinations (20.5% of the total number of individuals) is thought to result from mitochondrial introgression, consistent with a scenario of demographic, and presumably spatial, post-Pleistocene expansion of populations from northern regions into a secondary contact zone in the South China Sea. Mitonuclear discordance due to introgression obscures phylogenetic relationships for recently-diverged lineages, and cautions against the use of mitochondrial markers alone for species identification within the mottled rabbitfish species complex in the South China Sea region.