Integrated single-cell transcriptome and comparative genome analysis reveals the origin of intermuscular bones in zebrafish

The evolutionary process of intermuscular bones (IBs) is complex, the molecular regulatory mechanisms of their development are not clear, and even the genes involved in the evolution and development of IBs are poorly understood. In this study, comparative genomic analysis of four fish species with IBs and eleven fish species without IBs identified 106 genes that are more conservatively evolved in fish species with IBs, but highly variable in fish species without IBs. These genes are mainly involved in swimming behavior and BMP signaling pathways. We performed single-cell transcriptome sequencing of IBs origin tissues in zebrafish before and after IBs formation and found that osteoblasts and mesenchymal stem cells (MSCs) increased significantly after IBs formation. RNA velocity analysis showed that osteoblasts in IBs differentiate from MSCs, and the differentiation trajectory of MSCs into osteoblasts was successfully constructed by pseudo-time analysis. Combined with the results of multi-omics analysis, seven candidate genes associated with IBs development were screened and knocked out in zebrafish. It was found that foxn3 mutation resulted in a delay in IB development, whereas bmp6 mutation resulted in a total loss of IB. By comparing the transcriptome of IBs tissues between bmp6^+/+ zebrafish and bmp6^-/- zebrafish, we found that bmp6 deletion may inhibit the differentiation of MSCs into osteoblasts while promoting the formation of osteoclasts and ultimately inhibiting the formation of IBs. This study provides new insights into the molecular regulatory mechanisms and evolutionary processes of IB development.