Phylosymbiosis is defined as the relationship in which the microbiome recapitulates the phylogeny of the host and has been demonstrated in a variety of terrestrial organisms, although it has been understudied in fish, the most phylogenetically diverse vertebrate. Given that the species-specificity of fish microbiomes was detected in multiple body parts and differed by body parts, we assumed that the phylogenetic reflection of the microbiome would differ across body parts. Thus, we analyze the difference of phylosymbiotic relationships in the microbial communities found in three body parts (skin, gills, and intestine) of seven wild fish species from four families (Labridae, Sebastidae, Sparidae, and Rajidae) via 16S rRNA gene amplicon sequencing. Fishes were purchased at Docheon port market in Tongyeong City, Korea and were transported to nearby research institutes for aliveness. Mantel tests using dissimilarity values of microbiomes and hosts’ divergence times showed that the differences in microbial communities in all three body parts were related to the hosts’ divergence time. This pattern was the most pronounced in the skin. Furthermore, fishes from the same family showed similar bacterial compositions on their skins and gills, with clear differences depending on the family, with the exception of Labridae. These results suggest that the skin microbiome is particularly vulnerable to evolutionary pressures. We hypothesized that the evolution of the fish immune system and the difference in feeding habits induced the stronger phylosymbiotic signal in the skin. Collectively, this dataset will be useful for understanding the fish microbiome and give insights into phylosymbiosis of aquatic animals across body parts.

Two deep-sea eels collected from the East Mariana Basin in the western Pacific Ocean are described in this study. Based on their morphological features, two eel specimens were assumed to belong to the Gray cutthroat eel family, Synaphobranchidae. Mitochondrial DNA (mtDNA) genes have been widely used as genetic markers to identify fish species. To accurately identify the species of the two eel specimens, we sequenced the mitochondrial cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (16S rRNA) genes from the two eel specimens. The sequences from the specimens were 100% identical. The molecular phylogenetic tree confirmed that the two eel specimens were closely related to Synaphobranchus affinis with a bootstrap value of 100%. This is the first study to report new records of S. affinis from the East Mariana Basin in the western Pacific Ocean.

In this study, using Illumina sequencing, we sequenced first the complete mitochondrial genome (mitogenome) of two deep-sea eels, Histiobranchus bathybius and Simenchelys parasitica, collected from the East Mariana Basin in the Western Pacific Ocean. The complete length of the H. bathybius and S. parasitica mitogenomes were 16,696 and 16,687 bp, respectively, each containing 37 genes (13 protein-coding genes, 22 tRNA genes, and 2 ribosomal RNA genes). To enhance the accuracy of the identification of H. bathybius and S. parasitica, we performed a phylogenetic analysis of multiple deep-sea eels based on the mitochondrial DNA gene (cytochrome c oxidase subunit I [COI]) using the maximum likelihood method. Our phylogenetic tree analysis confirmed that the specimens collected in this study are congeneric species of H. bathybius and S. parasitica reported in previous studies. Based on these results, we report the first complete mitogenomes of H. bathybius and S. parasitica and a new record for the two species in the East Mariana Basin.