Wingless-type MMTV integration site family, member 16 (wnt16), is a wnt ligand that participates in the regulation of vertebrate skeletal development. Studies have shown that wnt16 can regulate bone metabolism, but its molecular mechanism remains largely undefined. We obtained the wnt16-/- zebrafish model using the CRISPR-Cas9-mediated gene knockout screen with 11 bp deletion in wnt16, which led to the premature termination of amino acid translation and significantly reduced wnt16 expression, thus obtaining the wnt16-/- zebrafish model. The expression of wnt16 in bone-related parts was detected via in situ hybridization. The head, spine, and tail exhibited significant deformities, and the bone mineral density and trabecular bone decreased in wnt16-/- using light microscopy and micro-CT analysis. RNA sequencing was performed to explore the differentially expressed genes (DEGs). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the down-regulated DEGs are mainly concentrated in mTOR, FoxO, and VEGF pathways. Protein–protein interaction (PPI) network analysis was performed with the detected DEGs. Eight down-regulated DEGs including akt1, bnip4, ptena, vegfaa, twsg1b, prkab1a, prkab1b, and pla2g4f.2 were validated by qRT-PCR and the results were consistent with the RNA-seq data. Overall, our work provides key insights into the influence of wnt16 gene on skeletal development.