Phylogeography of the crown-of-thorns starfish in the Indian Ocean

Cloning and Selfing Affect Population Genetic Variation in Simulations of Outcrossing, Sexual Sea Stars

AbstractMany sea stars are well known for facultative or obligate asexual reproduction in both the adult and larval life-cycle stages. Some species and lineages are also capable of facultative or obligate hermaphroditic reproduction with self-fertilization. However, models of population genetic variation and empirical analyses of genetic data typically assume only sexual reproduction and outcrossing. A recent reanalysis of previously published empirical data (microsatellite genotypes) from two studies of one of the most well-known sea star species (the crown-of-thorns sea star; Acanthaster sp.) concluded that cloning and self-fertilization in that species are rare and contribute little to patterns of population genetic variation. Here we reconsider that conclusion by simulating the contribution of cloning and selfing to genetic variation in a series of models of sea star demography. Simulated variation in two simple models (analogous to previous analyses of empirical data) was consistent with high rates of cloning or selfing or both. More realistic scenarios that characterize population flux in sea stars of ecological significance, including outbreaks of crown-of-thorns sea stars that devastate coral reefs, invasions by Asterias amurensis, and epizootics of sea star wasting disease that kill Pisaster ochraceus, also showed significant but smaller effects of cloning and selfing on variation within subpopulations and differentiation between subpopulations. Future models or analyses of genetic variation in similar study systems might benefit from simulation modeling to characterize possible contributions of cloning or selfing to genetic variation in population samples or to understand the limits on inferring the effects of cloning or selfing in nature.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

An Investigation into the Genetic History of Japanese Populations of Three Starfish, Acanthaster planci, Linckia laevigata, and Asterias amurensis, Based on Complete Mitochondrial DNA Sequences

Crown-of-thorns starfish, Acanthaster planci (COTS), are common in coral reefs of Indo-Pacific Ocean. Since they are highly fecund predators of corals, periodic outbreaks of COTS cause substantial loss of healthy coral reefs. Using complete mitochondrial DNA sequences, we here examined how COTS outbreaks in the Ryukyu Archipelago, Japan are reflected by the profile of their population genetics. Population genetics of the blue starfish, Linckia laevigata, which lives in the Ryukyu Archipelago, but not break out and the northern Pacific sea star, Asterias amurensis, which lives in colder seawater around the main Islands of Japan, were also examined as controls. Our results showed that As. amurensis has at least two local populations that diverged approximately 4.7 million years ago (MYA), and no genetic exchanges have occurred between the populations since then. Linckia laevigata shows two major populations in the Ryukyu Archipelago that likely diverged ∼6.8 MYA. The two populations, each comprised of individuals collected from coast of the Okinawa Island and those from the Ishigaki Island, suggest the presence of two cryptic species in the Ryukyu Archipelago. On the other hand, population genetics of COTS showed a profile quite different from those of Asterias and Linckia. At least five lineages of COTS have arisen since their divergence ∼0.7 MYA, and each of the lineages is present at the Okinawa Island, Miyako Island, and Ishigaki Island. These results suggest that COTS have experienced repeated genetic bottlenecks that may be associated with or caused by repeated outbreaks.

New insights from nuclear and mitochondrial markers on the genetic diversity and structure of the Indian white shrimp Fenneropenaeus indicus among the marginal seas in the Indian Ocean

Genetic variation in wild stocks of a major commercial shrimp, Fenneropenaeus indicus, from the marginal seas in the Indian Ocean was analysed using polymorphic microsatellite loci and mitochondrial COI gene. The average observed heterozygosity (Ho = 0.44 ± 0.02) and the expected heterozygosity (He = 0.73 ± 0.01) were high across loci and populations indicating high microsatellite variation. Pairwise F_ST and Bayesian clustering indicated the occurrence of four genetically distinct stocks out of the eight sampled populations with implications for specific management approaches. Mantel test for isolation by distance proved that genetic differentiation is not related to geographic distance between populations. Mitochondrial COI sequence analysis showed concordant differentiation pattern as well indicated the relevance of COI in population genetics of shrimps. Pairwise ɸST and phylogenetic and Bayesian analyses revealed four distinct clades, as observed with nuclear markers. Divergence time analysis revealed the origin and initial divergence of F. indicus corresponds to late Miocene and divergence to phylogroups in the Pleistocene. BSP analysis presented a long stable population size with a slight decrease in the late Pleistocene and gradually expanded to the current status. The information here will be useful in commercial shrimp breeding and selection programmes and management of natural stocks of Indian white shrimp.

Seascape genetics of the spiny lobster Panulirus homarus in the Western Indian Ocean: Understanding how oceanographic features shape the genetic structure of species with high larval dispersal potential

This study examines the fine-scale population genetic structure and phylogeography of the spiny lobster Panulirus homarus in the Western Indian Ocean. A seascape genetics approach was used to relate the observed genetic structure based on 21 microsatellite loci to ocean circulation patterns, and to determine the influence of latitude, sea surface temperature (SST), and ocean turbidity (KD490) on population-level processes. At a geospatial level, the genetic clusters recovered corresponded to three putative subspecies, P. h. rubellus from the SW Indian Ocean, P. h. megasculptus from the NW Indian Ocean, and P. h. homarus from the tropical region in-between. Virtual passive Lagrangian particles advected using satellite-derived ocean surface currents were used to simulate larval dispersal. In the SW Indian Ocean, the dispersion of particles tracked over a 4-month period provided insight into a steep genetic gradient observed at the Delagoa Bight, which separates P. h. rubellus and P. h. homarus. South of the contact zone, particles were advected southwestwards by prevailing boundary currents or were retained in nearshore eddies close to release locations. Some particles released in southeast Madagascar dispersed across the Mozambique Channel and reached the African shelf. Dispersal was characterized by high seasonal and inter-annual variability, and a large proportion of particles were dispersed far offshore and presumably lost. In the NW Indian Ocean, particles were retained within the Arabian Sea. Larval retention and self-recruitment in the Arabian Sea could explain the recent genetic divergence between P. h. megasculptus and P. h. homarus. Geographic distance and minimum SST were significantly associated with genetic differentiation in multivariate analysis, suggesting that larval tolerance to SST plays a role in shaping the population structure of P. homarus.

Multiple Facets of Marine Invertebrate Conservation Genomics

Conservation genomics aims to preserve the viability of populations and the biodiversity of living organisms. Invertebrate organisms represent 95% of animal biodiversity; however, few genomic resources currently exist for the group. The subset of marine invertebrates includes the most ancient metazoan lineages and possesses codes for unique gene products and possible keys to adaptation. The benefits of supporting invertebrate conservation genomics research (e.g., likely discovery of novel genes, protein regulatory mechanisms, genomic innovations, and transposable elements) outweigh the various hurdles (rare, small, or polymorphic starting materials). Here we review best conservation genomics practices in the laboratory and in silico when applied to marine invertebrates and also showcase unique features in several case studies of acroporid corals, crown-of-thorns starfish, apple snails, and abalone. Marine conservation genomics should also address how diversity can lead to unique marine innovations, the impact of deleterious variation, and how genomic monitoring and profiling could positively affect broader conservation goals (e.g., value of baseline data for in situ/ex situ genomic stocks).

Contributions of pre- versus post-settlement processes to fluctuating abundance of crown-of-thorns starfishes (Acanthaster spp.)

Numerous hypotheses have been put forward to account for population outbreaks of crown-of-thorns starfishes (CoTS, Acanthaster spp.), which place specific importance on either pre- or post-settlement mechanisms. The purpose of this review is to specifically assess the contributions of pre- versus post-settlement processes in the population dynamics of CoTS. Given the immense reproductive potential of CoTS (>100 million eggs per female), persistent high densities would appear inevitable unless there were significant constraints on larval development, settlement success, and/or early post-settlement growth and survival. In terms of population constraints, pre- and post-settlement processes are both important and have additive effects to suppress densities of juvenile and adult CoTS within reef ecosystems. It is difficult, however, to assess the relative contributions of pre- versus post-settlement processes to population outbreaks, especially given limited data on settlement rates, as well as early post-settlement growth and mortality. Prioritising this research is important to resolve potential effects of anthropogenic activities (e.g., fishing) and habitat degradation on changing population dynamics of CoTS, and will also improve management effectiveness.

Genetic population structure of the convict surgeonfish Acanthurus triostegus: a phylogeographic reassessment across its range

This study investigates the genetic population structure and connectivity of Acanthurus triostegus in five Indo-Pacific biogeographic regions (western and eastern Indian Ocean, western, central and eastern Pacific Ocean), using a mitochondrial DNA marker spanning the ATPase8 and ATPase6 gene regions. In order to assess the phylogeography and genetic population structure of A. triostegus across its range, 35 individuals were sampled from five localities in the western Indian Ocean and complemented with 227 sequences from two previous studies. Results from the overall analysis of molecular variance (AMOVA) without a priori grouping showed evidence of significant differentiation in the Indo-Pacific, with 25 (8.3%) out of 300 pairwise Φ_ST comparisons being significant. However, the hierarchical AMOVA grouping of Indian and Pacific Ocean populations failed to support the vicariance hypothesis, showing a lack of a genetic break between the two ocean basins. Instead, the correlation between pairwise Φ_ST values and geographic distance showed that dispersal of A. triostegus in the Indo-Pacific Ocean follows an isolation-by-distance model. Three haplogroups could be deduced from the haplotype network and phylogenetic tree, with haplogroup 1 and 2 dominating the Indian and the Pacific Ocean, respectively, while haplogroup 3 exclusively occurring in the Hawaiian Archipelago of the central Pacific Ocean.

Monitoring and assessing a 2-year outbreak of the corallivorous seastar Acanthaster planci in Ari Atoll, Republic of Maldives

Outbreaks of the corallivorous crown-of-thorns seastars have received increasing attention due to their negative impacts on coral reefs in the Indo-Pacific Ocean. However, outbreaks in remote and dislocated islands are still poorly understood. This study aims to begin filling informational gaps regarding outbreaks of Acanthaster planci in the remote islands of the central Ari Atoll, Republic of Maldives. The population of A. planci was monitored during three periods over 2 years (2015-2016) to evaluate variations in abundance and to characterise size structure and feeding behaviour. The outbreak appeared to be severe and active throughout the entire study period. The size structure analysis revealed a multimodal distribution dominated by individuals between 20 and 30 cm, suggesting that the outbreak may have resulted from a few nearby mass spawning events. Additionally, the most abundant live coral was Porites, which was also the most consumed genus; however, the electivity index showed a preference for corals of the genera Favites and Pavona. Finally, we also highlighted the need for more geographically extended surveys to better understand local patterns regarding outbreaks of A. planci in the Republic of Maldives.

Sea surface currents and geographic isolation shape the genetic population structure of a coral reef fish in the Indian Ocean

In this contribution, we determine the genetic population structure in the Skunk Clownfish (Amphiprion akallopsisos) across the Indian Ocean, and on a smaller geographic scale in the Western Indian Ocean (WIO). Highly restricted gene flow was discovered between populations on either side of the Indian Ocean using the control region as a mitochondrial marker (mtDNA). We verify this conclusion using 13 microsatellite markers and infer fine scale genetic structuring within the WIO. In total 387 samples from 21 sites were analysed using mtDNA and 13 microsatellite loci. Analysis included estimation of genetic diversity and population differentiation. A haplotype network was inferred using mtDNA. Nuclear markers were used in Bayesian clustering and a principal component analysis. Both markers confirmed strong genetic differentiation between WIO and Eastern Indian Ocean (EIO) populations, and a shallower population structure among Malagasy and East African mainland populations. Limited gene flow across the Mozambique Channel may be explained by its complex oceanography, which could cause local retention of larvae, limiting dispersal between Madagascar and the East African coast. Two other potential current-mediated barriers to larval dispersal suggested in the WIO, the split of the SEC at approximately 10° S and the convergence of the Somali Current with the East African Coast Current at approximately 3° S, were not found to form a barrier to gene flow in this species.

Effects of larvae density and food concentration on Crown-of-Thorns seastar (Acanthaster cf. solaris) development in an automated flow-through system

Coral-eating Crown-of-Thorns Sea stars (Acanthaster spp.) are major contributors to coral reef loss in the Indo-Pacific region. A release from food limitation of their planktotrophic larvae through enhanced pelagic productivity is one of the main hypothesis explaining population outbreaks (‘nutrient limitation hypothesis’). To improve the understanding of these outbreaks we developed an automated flow- through larvae rearing system that maintained food (microalgae) at set levels over the course of four 15d experiments. This resulted in stable food concentrations in experimental tanks. Increased algae concentrations had a significant positive effect on larval development and size at 10 and 15 days post fertilization (dpf). Larvae densities had no effect at 10 dpf. At 15 dpf greater larvae densities were associated with declines in larvae size. Larval development was slowed under higher larvae densities. Thus, the effects of algae concentration and larvae density were additive at 15 dpf, with larvae under low densities at a given algae concentration being further developed than those under higher densities. The development of a flow-through system gives greater insight into the effect of algae and larvae concentrations on Acanthaster development, and the system can be applied to further test the nutrient-limitation hypothesis for present and future outbreaks.

Mechanisms of peripheral phylogeographic divergence in the indo-Pacific: lessons from the spiny lobster Panulirus homarus

Background

There is increasing recognition of the concordance between marine biogeographic and phylogeographic boundaries. However, it is still unclear how population-level divergence translates into species-level divergence, and what are the principal factors that first initiate that divergence, and then maintain reproductive isolation. This study examines the likely forces driving population and lineage divergences in the broadly-distributed Indo-Pacific spiny lobster Panulirus homarus, which has peripheral divergent lineages in the west and east. The study focuses particularly on the West Indian Ocean, which is emerging as a region of unexpected diversity. Mitochondrial control region (mtCR) and COI sequences as well as genotypes of 9 microsatellite loci were examined in 410 individuals from 17 locations grouped into 7 regions from South Africa in the west, and eastward across to Taiwan and the Marquesas Islands. Phylogenetic and population-level analyses were used to test the significance and timing of divergences and describe the genetic relationships among populations.

Results

Analyses of the mtCR revealed high levels of divergence among the seven regions (Ф_ST = 0.594, P < 0.001). Microsatellite analyses also revealed significant divergence among regions, but at a much lower level (F_ST = 0.066, P < 0.001). The results reveal different patterns of mtCR v. nDNA divergence between the two distinct peripheral lineages: a subspecies in South Africa and Madagascar, and a phylogeographically diverged population in the Marquesas. The results also expose a number of other more fine-scale population divergences, particularly in the Indian Ocean.

Conclusions

The divergence of peripheral lineages in the west and east of the species’ range appear to have been initiated and maintained by very different processes. The pattern of mitochondrial and nuclear divergence of the western lineage, implicates processes of parapatric isolation, secondary contact and introgression, and suggests possible maintenance through adaptation and behavioural reproductive isolation. In contrast, the eastern lineage appears to have diverged through a rare colonisation event, maintained through long-term isolation, and matches expectations of the core-periphery hypothesis. The process of active peripheral speciation may be a common force in the Indo-Pacific that helps drive some of the regions’ recognized biogeographic boundaries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1050-8) contains supplementary material, which is available to authorized users.

Genetic Population Structure of the Coral Reef Sea Star Linckia laevigata in the Western Indian Ocean and Indo-West Pacific

The coral reef sea star Linckia laevigata is common on shallow water coral reefs of the Indo-West Pacific. Its large geographic distribution and comprehensive data from previous studies makes it suitable to examine genetic differentiation and connectivity over large geographical scales. Based on partial sequences of the mitochondrial cytochrome oxidase I (COI) gene this study investigates the genetic population structure and connectivity of L. laevigata in the Western Indian Ocean (WIO) and compares it to previous studies in the Indo-Malay-Philippines Archipelago (IMPA). A total of 138 samples were collected from nine locations in the WIO. AMOVA revealed a low but significant ΦST-value of 0.024 for the WIO populations. In the hierarchical AMOVA, the following grouping rejected the hypothesis of panmixia: (1) Kenya (Watamu, Mombasa, Diani) and Tanzanian Island populations (Misali and Jambiani) and (2) the rest of the WIO sites (mainland Tanzania and Madagascar; ΦCT = 0.03). The genetic population structure was stronger and more significant (ΦST = 0.13) in the comparative analysis of WIO and IMPA populations. Three clades were identified in the haplotype network. The strong genetic differentiation (ΦCT = 0.199, P < 0.001) suggests that Indo-West Pacific populations of L. laevigata can be grouped into four biogeographic regions: (1) WIO (2) Eastern Indian Ocean (3) IMPA and (4) Western Pacific. The findings of this study support the existence of a genetic break in the Indo-West Pacific consistent with the effect of lowered sea level during the Pleistocene, which limited gene flow between the Pacific and Indian Ocean.

Low genetic diversity and recent demographic expansion in the red starfish Echinaster sepositus (Retzius 1816)

Understanding the phylogeography and genetic structure of populations and the processes responsible of patterns therein is crucial for evaluating the vulnerability of marine species and developing management strategies. In this study, we explore how past climatic events and ongoing oceanographic and demographic processes have shaped the genetic structure and diversity of the Atlanto-Mediterranean red starfish Echinaster sepositus. The species is relatively abundant in some areas of the Mediterranean Sea, but some populations have dramatically decreased over recent years due to direct extraction for ornamental aquariums and souvenir industries. Analyses across most of the distribution range of the species based on the mitochondrial cytochrome c oxidase subunit I gene and eight microsatellite loci revealed very low intraspecific genetic diversity. The species showed a weak genetic structure within marine basins despite the a priori low dispersal potential of its lecithotrophic larva. Our results also revealed a very recent demographic expansion across the distribution range of the species. The genetic data presented here indicate that the species might be highly vulnerable, due to its low intraspecific genetic diversity.

Genetic structure of the crown-of-thorns seastar in the Pacific Ocean, with focus on Guam

Population outbreaks of the corallivorous crown-of-thorns seastar (COTS), Acanthaster 'planci' L., are among the most important biological disturbances of tropical coral reefs. Over the past 50 years, several devastating outbreaks have been documented around Guam, an island in the western Pacific Ocean. Previous analyses have shown that in the Pacific Ocean, COTS larval dispersal may be geographically restricted to certain regions. Here, we assess the genetic structure of Pacific COTS populations and compared samples from around Guam with a number of distant localities in the Pacific Ocean, and focused on determining the degree of genetic structure among populations previously considered to be isolated. Using microsatellites, we document substantial genetic structure between 14 localities from different geographical regions in the Pacific Ocean. Populations from the 14 locations sampled were found to be structured in three significantly differentiated groups: (1) all locations immediately around Guam, as well as Kingman Reef and Swains Island; (2) Japan, Philippines, GBR and Vanuatu; and (3) Johnston Atoll, which was significantly different from all other localities. The lack of genetic differentiation between Guam and extremely distant populations from Kingman Reef and Swains Island suggests potential long-distance dispersal of COTS in the Pacific.

Putting the ‘Indo’ back into the Indo-Pacific: resolving marine phylogeographic gaps

The Indo-Pacific is an extremely large marine realm that unites two oceans via a restricted Coral Triangle corridor, which was historically subjected to lowered sea levels during global glaciation. Although a strong phylogeographic focus on the Central and West Pacific has produced a large body of research, the Indian Ocean has been largely neglected. This may have serious consequences, because the Indian Ocean hosts a large number of marine centres of endemism, yet a large number of nations rely on its marine resources. We examine reasons for this neglect and review what is known about this region and its connectivity to the Indo-West Pacific. We draw attention to the ‘Leeuwin Effect’, a phenomenon where the southward flow of the Leeuwin Current is responsible for transporting larval propagules from the Coral Triangle region down the coast of Western Australia, resulting in broader Indo-West Pacific rather than Indian Ocean affinities. Given challenges in accessing infrastructure and samples, collaboration will inevitably be key to resolving data gaps. We challenge the assumption that the peak of shallow-water marine biodiversity is solely centred in the Coral Triangle, and raise awareness of a seemingly forgotten hypothesis promoting a secondary peak of biodiversity in the western Indian Ocean.

Reefs and islands of the Chagos Archipelago, Indian Ocean: why it is the world's largest no-take marine protected area

The Chagos Archipelago was designated a no-take marine protected area (MPA) in 2010; it covers 550 000 km2, with more than 60 000 km2 shallow limestone platform and reefs. This has doubled the global cover of such MPAs.It contains 25-50% of the Indian Ocean reef area remaining in excellent condition, as well as the world's largest contiguous undamaged reef area. It has suffered from warming episodes, but after the most severe mortality event of 1998, coral cover was restored after 10 years.Coral reef fishes are orders of magnitude more abundant than in other Indian Ocean locations, regardless of whether the latter are fished or protected.Coral diseases are extremely low, and no invasive marine species are known.Genetically, Chagos marine species are part of the Western Indian Ocean, and Chagos serves as a 'stepping-stone' in the ocean.The no-take MPA extends to the 200 nm boundary, and. includes 86 unfished seamounts and 243 deep knolls as well as encompassing important pelagic species.On the larger islands, native plants, coconut crabs, bird and turtle colonies were largely destroyed in plantation times, but several smaller islands are in relatively undamaged state.There are now 10 'important bird areas', coconut crab density is high and numbers of green and hawksbill turtles are recovering.Diego Garcia atoll contains a military facility; this atoll contains one Ramsar site and several 'strict nature reserves'. Pollutant monitoring shows it to be the least polluted inhabited atoll in the world. Today, strict environmental regulations are enforced.Shoreline erosion is significant in many places. Its economic cost in the inhabited part of Diego Garcia is very high, but all islands are vulnerable.Chagos is ideally situated for several monitoring programmes, and use is increasingly being made of the archipelago for this purpose.