Intermuscular bones (IBs) are mineralized spicules, present in the myosepta of many, but not all, teleost species. IBs are often small and sharp, and they consequently limit how the fish can be processed; the IBs may cause injury or trauma if lodged in consumers’ throats or mouths, and therefore affect the appeal of the fish to many consumers. The development of IBs in teleosts is still not fully understood and the molecular basis of IB development remains to be established. Here, the characteristics of IB tissue are evaluated based on single-cell transcriptomics in wild-type zebrafish. The analysis defined 18 distinct cell types. Differentiation trajectories showed that IBs are derived from tendons and that a core tendon-osteoblast cell lineage is related to IB formation. In particular, the functions of 10 candidate genes were evaluated via CRISPR-Cas9 mutants. Among those, runx2b^−/− mutants completely lost IBs, while swimming performance, growth and bone mineral density were not significantly different from runx2b^+/+ zebrafish. Comparative single-cell RNA sequencing (scRNA-seq) analysis in runx2b^−/− and runx2b^+/+ zebrafish revealed the role of osteoblasts in IB formation. In addition, differentially expressed genes were enriched in the transforming growth factor β/bone morphogenetic protein (TGF-β/BMP) pathway after runx2b deletion. This study provides evidence for the crucial role of runx2b regulation in IB formation. Genetic breeding can target runx2b regulation and generate strains of commercial fish species without IBs, which can improve the safe consumption and economic value of many farmed fish species.

Bone morphogenetic protein 7 (BMP7) belongs to the transforming growth factor β (TGF-β) family, which not only induces cartilage and bone formation, but also regulates eye development and melanoma tumorigenesis in mammals. In teleosts, BMP7 differentiates into two subtypes, bmp7a and bmp7b, which have clearly differentiated structures. To fully understand the functional differentiation of bmp7a and bmp7b in fish species, we successfully constructed bmp7a and bmp7b gene deletion mutants in zebrafish using CRISPR/Cas9-mediated gene editing technology. Our results showed that bmp7a mutation caused abnormal development of the embryo’s dorsal-ventral pattern that led to death; bmp7b mutation induced growth inhibition and increased melanin production in the skin and eye of mutants. Histological analysis revealed that melanin in the retina of the eyes in bmp7b mutants increased, and behavioral observation showed that the vision and sensitivity to food of the mutants were reduced. Transcriptome analysis of the skin and eye tissues showed that the expression changes of wnt7ba and gna14 in bmp7b mutants might promote the increase of melanin. Additionally, the eye transcriptome analysis indicated that changes in the structure of the eyes in bmp7b mutants led to defects in phototransduction, and seven DEGs (rgs9a, rgs9b, rcvrn2, guca1d, grk1b, opn1mw4, and gc2) were identified as key candidate genes that affected the photonic response of the eyes. The study revealed the functional differentiation of bmp7a and bmp7b in teleosts and the first report about the inhibitory effect of bmp7b on melanogenesis may provide useful information for the future research on human melanoma-related diseases.

Single nucleotide polymorphisms (SNPs) within microRNAs (miRNAs) and their target binding sites can influence miRNA biogenesis and target regulation, thereby participating in a variety of diseases and biological processes. Current miRNA-related SNP databases are often species-limited or based on outdated data. Therefore, we updated our miRNASNP database to version 4 by updating data, expanding the species from Homo sapiens to 17 species, and introducing several new features. In miRNASNP-v4, 82 580 SNPs in miRNAs and 24 836 179 SNPs in 3′UTRs of genes across 17 species were identified and their potential effects on miRNA secondary structure and target binding were characterized. In addition, compared to the last release, miRNASNP-v4 includes the following improvements: (i) gene enrichment analysis for gained or lost miRNA target genes; (ii) identification of miRNA-related SNPs associated with drug response and immune infiltration in human cancers; (iii) inclusion of experimentally supported immune-related miRNAs and (iv) online prediction tools for 17 animal species. With the extensive data and user-friendly web interface, miRNASNP-v4 will serve as an invaluable resource for functional studies of SNPs and miRNAs in multiple species. The database is freely accessible at http://gong_lab.hzau.edu.cn/miRNASNP/.

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.