Abstract
Reef stonefish (Synanceia verrucosa) is one of the most venomous fishes, but its biomedical study has been restricted to molecular cloning and purification of its toxins, instead of high-throughput genetic research on related toxin genes. In this study, we constructed a chromosome-level haplotypic genome assembly for the reef stonefish. The genome was assembled into 24 pseudo-chromosomes, and the length totaled 689.74 Mb, reaching a contig N50 of 11.97 Mb and containing 97.8% of complete BUSCOs. A total of 24,050 protein-coding genes were annotated, of which metalloproteinases, C-type lectins, and stonustoxins (sntx) were the most abundant putative toxin genes. Multitissue transcriptomic and venom proteomic data showed that sntx genes, especially those clustered within a 50-kb region on the chromosome 2, had higher transcription levels than other types of toxins as well as those sntx genes scatteringly distributed on other chromosomes. Further comparative genomic analysis predicted an expansion of sntx-like genes in the Percomorpha lineage including nonvenomous fishes, but Scorpaenoidei species experienced extra independent sntx duplication events, marking the clear-cut origin of authentic toxic stonustoxins. In summary, this high-quality genome assembly and related comparative analysis of toxin genes highlight valuable genetic differences for potential involvement in the evolution of venoms among Scorpaeniformes fishes.