Near-Chromosomal-Level Genome Assembly of the Sea Urchin Echinometra lucunter, a Model for Speciation in the Sea

Echinometra lucunter, the rock-boring sea urchin, is a widely distributed echinoid and a model for ecological studies of reproduction, responses to climate change, and speciation. We present a near chromosome-level genome assembly of E. lucunter, including 21 scaffolds larger than 10 Mb predicted to represent each of the chromosomes of the species. The 760.4 Mb assembly includes a scaffold N50 of 30.0 Mb and BUSCO (benchmarking universal single-copy orthologue) single copy and a duplicated score of 95.8% and 1.4%, respectively. Ab-initio gene model prediction and annotation with transcriptomic data constructed 33,989 gene models composing 50.4% of the assembly, including 37,036 transcripts. Repetitive elements make up approximately 39.6% of the assembly, and unresolved gap sequences are estimated to be 0.65%. Whole genome alignment with Echinometra sp. EZ revealed high synteny and conservation between the two species, further bolstering Echinometra as an emerging genus for comparative genomics studies. This genome assembly represents a high-quality genomic resource for future evolutionary and developmental studies of this species and more broadly of echinoderms.

Accumulation of gene copy number variations during the early phase of free-spawning abalone speciation

The genetic basis of speciation in free-spawning marine invertebrates is poorly understood. Although gene copy number variations (GCNVs) and nucleotide variations possibly trigger the speciation of these organisms, empirical evidence for such a hypothesis is limited. In this study, we searched for genomic signatures of GCNVs that may contribute to the speciation of Western Pacific abalone species. Whole-genome sequencing data suggested the existence of significant amounts of GCNVs in closely related abalones, Haliotis discus and H. madaka, in the early phase of speciation. In addition, the degree of interspecies genetic differentiation in the genes where GCNVs were estimated was higher than that in other genes, suggesting that nucleotide divergence also accumulates in the genes with GCNVs. GCNVs in some genes were also detected in other related abalone species, suggesting that these GCNVs are derived from both ancestral and de novo mutations. Our findings suggest that GCNVs have been accumulated in the early phase of free-spawning abalone speciation.

Ten years of demographic modelling of divergence and speciation in the sea

Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the 'grey zone of speciation' for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Phylogeography of Lionfishes (Pterois) Indicate Taxonomic Over Splitting and Hybrid Origin of the Invasive Pterois volitans

The lionfish is an iconic marine fish, and recently renowned for a disastrous introduction into the West Atlantic. Genetic surveys of the putative invaders (Pterois volitans and Pterois miles) in their natural Indo-Pacific range can illuminate both topics. Previous research indicated that P. volitans and P. miles are sister species that hybridize in the invasive range, but hybridization in the native range is unknown. Here, we apply mtDNA COI and 2 nuclear introns (S7 RP1 and Gpd2) from 229 lionfish including the 2 invaders and 2 closely-related taxa (44 P. miles, 91 P. volitans, 31 Pterois lunulata, and 63 Pterois russelii) from 10 locations in their native ranges. Genetic data are supplemented with key morphological characters: dorsal, anal, and pectoral fin ray counts. We observed 2 lineages (d = 4.07%, 0.89%, and 2.75% at COI, S7 RP1, and Gpd2, respectively) among the 4 putative species: an Indian Ocean lineage represented by P. miles, and a Pacific Ocean lineage represented by P. lunulata and P. russelii. All specimens of the invasive P. volitans appear to be hybrids between the Indian Ocean P. miles and a Pacific lineage encompassing P. lunulata/russelii, a conclusion supported by both genetics and morphology. The divergences between Indian and Pacific forms are within the range of species-level partitions in fishes, and we recommend retention of the names P. miles and P. russelii for Indian and Pacific forms. The hybrid origin of the Atlantic invasion invokes the possibility of heterosis as a contributing factor to invasion success.

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.