The West Pacific diversity hotspot as a source or sink for new species? Population genetic insights from the Indo-Pacific parrotfish Scarus rubroviolaceus

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

The origin of the parrotfish species Scarus compressus in the Tropical Eastern Pacific: region-wide hybridization between ancient species pairs

Background

In the Tropical Eastern Pacific (TEP), four species of parrotfishes with complex phylogeographic histories co-occur in sympatry on rocky reefs from Baja California to Ecuador: Scarus compressus, S. ghobban, S. perrico, and S. rubroviolaceus. The most divergent, S. perrico, separated from a Central Indo-Pacific ancestor in the late Miocene (6.6 Ma). We tested the hypothesis that S. compressus was the result of ongoing hybridization among the other three species by sequencing four nuclear markers and a mitochondrial locus in samples spanning 2/3 of the latitudinal extent of the TEP.

Results

A Structure model indicated that K = 3 fit the nuclear data and that S. compressus individuals had admixed genomes. Our data could correctly detect and assign pure adults and F1 hybrids with > 0.90 probability, and correct assignment of F2s was also high in some cases. NewHybrids models revealed that 89.8% (n = 59) of the S. compressus samples were F1 hybrids between either S. perrico × S. ghobban or S. perrico × S. rubroviolaceus. Similarly, the most recently diverged S. ghobban and S. rubroviolaceus were hybridizing in small numbers, with half of the admixed individuals assigned to F1 hybrids and the remainder likely > F1 hybrids. We observed strong mito-nuclear discordance in all hybrid pairs. Migrate models favored gene flow between S. perrico and S. ghobban, but not other species pairs.

Conclusions

Mating between divergent species is giving rise to a region-wide, multispecies hybrid complex, characterized by a high frequency of parental and F1 genotypes but a low frequency of > F1 hybrids. Trimodal structure, and evidence for fertility of both male and female F1 hybrids, suggest that fitness declines sharply in later generation hybrids. In contrast, the hybrid population of the two more recently diverged species had similar frequencies of F1 and > F1 hybrids, suggesting accelerating post-mating incompatibility with time. Mitochondrial genotypes in hybrids suggest that indiscriminate mating by male S. perrico is driving pre-zygotic breakdown, which may reflect isolation of this endemic species for millions of years resulting in weak selection for conspecific mate recognition. Despite overlapping habitat use and high rates of hybridization, species boundaries are maintained by a combination of pre- and post-mating processes in this complex.

Cryptic ecological and geographic diversification in coral-associated nudibranchs

Coral reefs are among the most biologically diverse ecosystems of the world, yet little is known about the processes creating and maintaining their diversity. Ecologically, corallivory in nudibranchs resembles phytophagy in insects- a process that for decades has served as a model for ecological speciation via host shifting. This study uses extensive field collections, DNA sequencing, and phylogenetic analyses to reconstruct the evolutionary history of coral-associated nudibranchs and assess the relative roles that host shifting and geography may have played in their diversification. We find that the number of species is three times higher than the number previously known to science, with evidence for both allopatric and ecological divergence through host shifting and host specialization. Results contribute to growing support for the importance of ecological diversification in marine environments and provide evidence for new species in the genus Tenellia.

Contrasting population genetic structure in three aggregating groupers (Percoidei: Epinephelidae) in the Indo-West Pacific: the importance of reproductive mode

Background

Understanding the factors shaping population genetic structure is important for evolutionary considerations as well as for management and conservation. While studies have revealed the importance of palaeogeographic changes in shaping phylogeographic patterns in multiple marine fauna, the role of reproductive behaviour is rarely considered in reef fishes. We investigated the population genetics of three commercially important aggregating grouper species in the Indo-West Pacific, namely the camouflage grouper Epinephelus polyphekadion, the squaretail coral grouper Plectropomus areolatus, and the common coral trout P. leopardus, with similar life histories but distinct spatio-temporal characteristics in their patterns of forming spawning aggregations.

Results

By examining their mitochondrial control region and 9–11 microsatellite markers, we found an overarching influence of palaeogeographic events in the population structure of all species, with genetic breaks largely coinciding with major biogeographic barriers. The divergence time of major lineages in these species coincide with the Pleistocene glaciations. Higher connectivity is evident in E. polyphekadion and P. areolatus that assemble in larger numbers at fewer spawning aggregations and in distinctive offshore locations than in P. leopardus which has multiple small, shelf platform aggregations.

Conclusions

While palaeogeographic events played an important role in shaping the population structure of the target species, the disparity in population connectivity detected may be partly attributable to differences in their reproductive behaviour, highlighting the need for more investigations on this characteristic and the need to consider reproductive mode in studies of connectivity and population genetics.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1284-0) contains supplementary material, which is available to authorized users.

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.

Growth and maturation of the redlip parrotfish Scarus rubroviolaceus

This study presents age-based life-history information for the red lip parrotfish Scarus rubroviolaceus based on a 5 year sampling programme from the commercial fishery of American Samoa. Females reached sexual maturity at 31·9 cm fork length (L_F ) and 2·6 years and sex change occurred at 42·3 cm L_F , although not all females change sex through their ontogeny. The maximum observed age was 14 years and c. 65% of the fishery harvest was above the median L_F at sex change.

The role of peripheral endemism and habitat associations in the evolution of the Indo-West Pacific tuskfishes (Labridae: Choerodon)

The unrivalled level of biodiversity across the tropical Indo-Australian Archipelago (IAA) has been the subject of wide debate. Attempts to understand its origins have focussed on the timing of speciation, rates of diversification and the directionality of colonisation across geographical and climatic gradients in an array of marine groups. We investigate origins and evolution in the Choerodon tuskfishes, a group of labrids whose centre of diversity coincides with this region. Mitochondrial (COI, 16S) and nuclear (RAG2, Tmo4c4) molecular phylogenies and biogeographic analyses, coupled with molecular clock dating, were inferred from 19 of the 23 valid Choerodon species. Two additional, undescribed Choerodon species were also included, showing reciprocal monophyly in both genomes, confirming their species level status. Choerodon diverged from their ancestral sister group, the Odacines, at the onset of the Miocene, coinciding with the collision of the Australian and Eurasian Plates when extensive areas of shallow-water habitat formed. Despite subsequent evolutionary patterns being partially obscured by overlapping distribution ranges between many species and a lack of clear evidence for climatically driven lineage divergences, our data support an evolutionary scenario of peripheral endemics budding from once widespread populations across this biodiversity hotspot. Interestingly, these peripheral endemics tend to occupy more specialised reef or non-reef habitats whereas widespread groups appear to generally take advantage of both reef and non-reef environments. Our results are discussed in light of the most accredited hypotheses proposed to explain species richness in the IAA, with some support for processes such as centrifugal speciation.

Multiple colonizations lead to cryptic biodiversity in an island ecosystem: comparative phylogeography of anchialine shrimp species in the Ryukyu Archipelago, Japan

Archipelagos of the Indo-West Pacific are considered to be among the richest in the world in biodiversity, and phylogeographic studies generally support either the center of origin or the center of accumulation hypothesis to explain this pattern. To differentiate between these competing hypotheses for organisms from the Indo-West Pacific anchialine ecosystem, defined as coastal bodies of mixohaline water fluctuating with the tides but having no direct oceanic connections, we investigated the genetic variation, population structure, and evolutionary history of three caridean shrimp species (Antecaridina lauensis, Halocaridinides trigonophthalma, and Metabetaeus minutus) in the Ryukyu Archipelago, Japan. We used two mitochondrial genes--cytochrome c oxidase subunit I (COI) and large ribosomal subunit (16S-rDNA)--complemented with genetic examination of available specimens from the same or closely related species from the Indian and Pacific Oceans. In the Ryukyus, each species encompassed 2-3 divergent (9.52%-19.2% COI p-distance) lineages, each having significant population structure and varying geographic distributions. Phylogenetically, the A. lauensis and M. minutus lineages in the Ryukyus were more closely related to ones from outside the archipelago than to one another. These results, when interpreted in the context of Pacific oceanographic currents and geologic history of the Ryukyus, imply multiple colonizations of the archipelago by the three species, consistent with the center of accumulation hypothesis. While this study contributes toward understanding the biodiversity, ecology, and evolution of organisms in the Ryukyus and the Indo-West Pacific, it also has potential utility in establishing conservation strategies for anchialine fauna of the Pacific Basin in general.

Phylogeography of the Indo-West Pacific maskrays (Dasyatidae, Neotrygon): a complex example of chondrichthyan radiation in the Cenozoic

Maskrays of the genus Neotrygon (Dasyatidae) have dispersed widely in the Indo-West Pacific being represented largely by an assemblage of narrow-ranging coastal endemics. Phylogenetic reconstruction methods reproduced nearly identical and statistically robust topologies supporting the monophyly of the genus Neotrygon within the family Dasyatidae, the genus Taeniura being consistently basal to Neotrygon, and Dasyatis being polyphyletic to the genera Taeniurops and Pteroplatytrygon. The Neotrygon kuhlii complex, once considered to be an assemblage of color variants of the same biological species, is the most derived and widely dispersed subgroup of the genus. Mitochondrial (COI, 16S) and nuclear (RAG1) phylogenies used in synergy with molecular dating identified paleoclimatic fluctuations responsible for periods of vicariance and dispersal promoting population fragmentation and speciation in Neotrygon. Signatures of population differentiation exist in N. ningalooensis and N. annotata, yet a large-scale geological event, such as the collision between the Australian and Eurasian Plates, coupled with subsequent sea-level falls, appears to have separated a once homogeneous population of the ancestral form of N. kuhlii into southern Indian Ocean and northern Pacific taxa some 4–16 million years ago. Repeated climatic oscillations, and the subsequent establishment of land and shallow sea connections within and between Australia and parts of the Indo-Malay Archipelago, have both promoted speciation and established zones of secondary contact within the Indian and Pacific Ocean basins.

Population differentiation in the open sea: insights from the pelagic copepod Pleuromamma xiphias

Although a number of recent studies of marine holoplankton have reported significant genetic structure among populations, little is currently known about the biological and oceanographic processes that influence population connectivity in oceanic plankton. In order to examine how depth preferences influence dispersal in oceanic plankton, I characterized the genetic structure of a copepod with diel vertical migration (DVM) (Pleuromamma xiphias), throughout its global distribution, and compared these results to those expected given the interaction of this species' habitat depth with ocean circulation and bathymetry. Mitochondrial COI sequences from 651 individuals from 28 sites in the Indian, Pacific, and Atlantic Oceans revealed highly significant genetic differentiation both within and among ocean basins. Limited dispersal among distinct pelagic provinces seems to have played a major role in population differentiation in this species, with strong genetic breaks observed across known oceanographic fronts or current systems in all three ocean basins. The Indo-West Pacific (IWP) holds a highly distinct genetic population of this species that was sampled in both the western Pacific and eastern Indian Oceans. This suggests that the IWP does not act as a strong barrier to gene flow between basins, as expected, despite the relatively shallow water depth (<200 m) and vertically extensive (>400 m) diel migration of this species. A pattern of isolation by distance was observed in the Indian Ocean with genetic differentiation among samples down to spatial scales of ∼800 km, indicating that realized dispersal in P. xiphias occurs over much smaller spatial scales than in previously reported oceanic holoplankton. Given its highly regionalized population genetic structure, P. xiphias may have some capacity to adapt to local oceanographic conditions, and it should not be assumed that populations of this species in distinct pelagic biomes will respond in the same way to shared physical or climatic forcing.

High connectivity in the deepwater snapper Pristipomoides filamentosus (Lutjanidae) across the Indo-Pacific with isolation of the Hawaiian archipelago

In the tropical Indo-Pacific, most phylogeographic studies have focused on the shallow-water taxa that inhabit reefs to approximately 30 m depth. Little is known about the large predatory fishes, primarily snappers (subfamily Etelinae) and groupers (subfamily Epinephelinae) that occur at 100-400 m. These long-lived, slow-growing species support fisheries across the Indo-Pacific, yet no comprehensive genetic surveys within this group have been conducted. Here we contribute the first range-wide survey of a deepwater Indo-Pacific snapper, Pristipomoides filamentosus, with special focus on Hawai'i. We applied mtDNA cytochrome b and 11 microsatellite loci to 26 samples (N=1,222) collected across 17,000 km from Hawai'i to the western Indian Ocean. Results indicate that P. filamentosus is a highly dispersive species with low but significant population structure (mtDNA Φ(ST)=0.029, microsatellite F(ST)=0.029) due entirely to the isolation of Hawai'i. No population structure was detected across 14,000 km of the Indo-Pacific from Tonga in the Central Pacific to the Seychelles in the western Indian Ocean, a pattern rarely observed in reef species. Despite a long pelagic phase (60-180 days), interisland dispersal as adults, and extensive gene flow across the Indo-Pacific, P. filamentosus is unable to maintain population connectivity with Hawai'i. Coalescent analyses indicate that P. filamentosus may have colonized Hawai'i 26 K-52 K y ago against prevailing currents, with dispersal away from Hawai'i dominating migration estimates. P. filamentosus harbors low genetic diversity in Hawai'i, a common pattern in marine fishes, and our data indicate a single archipelago-wide stock. However, like the Hawaiian Grouper, Hyporthodus quernus, this snapper had several significant pairwise comparisons (F(ST)) clustered around the middle of the archipelago (St. Rogatien, Brooks Banks, Gardner) indicating that this region may be isolated or (more likely) receives input from Johnston Atoll to the south.

Phylogeography of the reef fish Cephalopholis argus (Epinephelidae) indicates Pleistocene isolation across the Indo-Pacific Barrier with contemporary overlap in The Coral Triangle

BACKGROUND

The Coral Triangle (CT), bounded by the Philippines, the Malay Peninsula, and New Guinea, is the epicenter of marine biodiversity. Hypotheses that explain the source of this rich biodiversity include 1) the center of origin, 2) the center of accumulation, and 3) the region of overlap. Here we contribute to the debate with a phylogeographic survey of a widely distributed reef fish, the Peacock Grouper (Cephalopholis argus; Epinephelidae) at 21 locations (N = 550) using DNA sequence data from mtDNA cytochrome b and two nuclear introns (gonadotropin-releasing hormone and S7 ribosomal protein).

RESULTS

Population structure was significant (ΦST = 0.297, P < 0.001; FST = 0.078, P < 0.001; FST = 0.099, P < 0.001 for the three loci, respectively) among five regions: French Polynesia, the central-west Pacific (Line Islands to northeastern Australia), Indo-Pacific boundary (Bali and Rowley Shoals), eastern Indian Ocean (Cocos/Keeling and Christmas Island), and western Indian Ocean (Diego Garcia, Oman, and Seychelles). A strong signal of isolation by distance was detected in both mtDNA (r = 0.749, P = 0.001) and the combined nuclear loci (r = 0.715, P < 0.001). We detected evidence of population expansion with migration toward the CT. Two clusters of haplotypes were detected in the mtDNA data (d = 0.008), corresponding to the Pacific and Indian Oceans, with a low level of introgression observed outside a mixing zone at the Pacific-Indian boundary.

CONCLUSIONS

We conclude that the Indo-Pacific Barrier, operating during low sea level associated with glaciation, defines the primary phylogeographic pattern in this species. These data support a scenario of isolation on the scale of 105 year glacial cycles, followed by population expansion toward the CT, and overlap of divergent lineages at the Pacific-Indian boundary. This pattern of isolation, divergence, and subsequent overlap likely contributes to species richness at the adjacent CT and is consistent with the region of overlap hypothesis.

Escaping paradise: Larval export from Hawaii in an Indo-Pacific reef fish, the Yellow Tang (Zebrasoma flavescens)

The depauperate marine ecosystems of the Hawaiian Archipelago share a high proportion of species with the southern and western Pacific, indicating historical and/or ongoing connections across the large oceanic expanse separating Hawaii from its nearest neighbors. The rate and direction of these interactions are, however, unknown. While previous biogeographic studies have consistently described Hawaii as a diversity sink, prevailing currents likely offer opportunities for larval export. To assess interactions between the remote reefs of the Hawaiian Archipelago and the species rich communities of the Central and West Pacific, we surveyed 14 nuclear microsatellite loci (nDNA; n = 857) and a 614 bp segment of mitochondrial cytochrome b (mtDNA; n = 654) in the Yellow Tang (Zebrasoma flavescens). Concordant frequency shifts in both nDNA and mtDNA reveal significant population differentiation among three West Pacific sites and Hawaii (nDNA F' _CT = 0.116, mtDNA ϕ _CT = 0.098, P < 0.001). SAMOVA analyses of microsatellite data additionally indicate fine scale differentiation within the 2600 km Hawaiian Archipelago (F' _SC = 0.026; P < 0.001), with implications for management of this heavily-exploited aquarium fish. Mismatch analyses indicate the oldest contemporary populations are in the Hawaiian Archipelago (circa 318,000 y), with younger populations in the West Pacific (91,000 - 175,000 y). Estimates of Yellow Tang historical demography contradict expectations of Hawaii as a population sink, and instead indicate asymmetrical gene flow, with Hawaii exporting rather than importing Yellow Tang larvae.