Molecular mechanisms of intermuscular bone development in fish: a review

Deletion of Asb15b gene can lead to a significant decrease in zebrafish intermuscular bone

Intermuscular bones, which are present in numerous economically significant fish species, have a negative impact on the development of aquaculture. The Asb15b gene, primarily expressed in skeletal muscle, plays a crucial role in regulating protein turnover and the development of muscle fibers. It stimulates protein synthesis and controls the differentiation of muscle fibers. In this study, we employed CRISPR/Cas9 technology to generate homozygous zebrafish strains with 7 bp and 49 bp deletions in the Asb15b gene. Subsequent analyses using skeleton staining demonstrated a substantial reduction in the number of intermuscular bones in adult Asb15b-/- -7 bp and Asb15b-/- -49 bp mutants compared to the wild-type zebrafish, with decreases of 30 % (P < 0.001) and 40 % (P < 0.0001), respectively. Histological experiments further revealed that the diameter and number of muscle fibers in adult Asb15b-/- mutants did not exhibit significant changes when compared to wild-type zebrafish. Moreover, qRT-PCR experiments demonstrated significant differences in the expression of bmp6 and runx2b genes, which are key regulators of intermuscular bone development, during different stages of intermuscular bone development in Asb15b-/- mutants. This study strongly suggests that the Asb15b gene plays a crucial role in regulating intermuscular bone development in fish and lays the groundwork for further exploration of the role of the Asb15b gene in zebrafish intermuscular bone development.

Genetic regulation of injury-induced heterotopic ossification in adult zebrafish

Heterotopic ossification is the inappropriate formation of bone in soft tissues of the body. It can manifest spontaneously in rare genetic conditions or as a response to injury, known as acquired heterotopic ossification. There are several experimental models for studying acquired heterotopic ossification from different sources of damage. However, their tenuous mechanistic relevance to the human condition, invasive and laborious nature and/or lack of amenability to chemical and genetic screens, limit their utility. To address these limitations, we developed a simple zebrafish injury model that manifests heterotopic ossification with high penetrance in response to clinically emulating injuries, as observed in human myositis ossificans traumatica. Using this model, we defined the transcriptional response to trauma, identifying differentially regulated genes. Mutant analyses revealed that an increase in the activity of the potassium channel Kcnk5b potentiates injury response, whereas loss of function of the interleukin 11 receptor paralogue (Il11ra) resulted in a drastically reduced ossification response. Based on these findings, we postulate that enhanced ionic signalling, specifically through Kcnk5b, regulates the intensity of the skeletogenic injury response, which, in part, requires immune response regulated by Il11ra.

Effect of runx2b deficiency in intermuscular bones on the regulatory network of lncRNA-miRNA-mRNA

Intermuscular bones (IBs) are mineralized spicules that negatively impact the quality and value of fish products. Runx2b is a crucial modulator in promoting bone formation through regulating osteoblast differentiation. Previous studies suggested that loss of runx2b gene completely inhibited IBs formation in zebrafish. However, how the whole transcriptome, including mRNA and non-coding RNA (ncRNA), affects the IBs development in runx2b-/- zebrafish are not known. The aim of this study was to identify the regulatory networks of differentially expressed (DE) lncRNAs, miRNAs, and mRNAs in zebrafish with and without IBs (runx2b+/+ fish and runx2b-/- fish) utilizing high-throughput sequencing techniques. All together there are 1051 mRNAs, 456 lncRNAs, and 18 miRNAs differentially expressed were found between these two strains. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) has highlighted significant pathways linked to the development of IBs, specifically the TGF-beta and Wnt signaling pathways, and a number of genes concentrated on these two signaling pathways related to the formation of IBs. Further, 1989 competing endogenous RNA (ceRNA) networks were created according to the correlation among mRNAs, miRNAs and lncRNAs. The ceRNA networks results revealed 52 ceRNA pairs related to the IBs formation, consisting of 52 mRNAs, 37 lncRNAs, and 6 miRNAs. Of these, we found that dre-miR-2189 was the key element of ceRNA pairs, interacting with 19 mRNAs and 11 lncRNAs, and MSTRG.13175.1 could regulate sp7 expression by interacting with dre-miR-2189 to function in osteogenic differentiation. Subsequent experiments at the cellular level also revealed the interaction mechanism. The outcomes indicated a crucial role of miRNAs and lncRNAs in the development of fish IBs, which offer new views into the functions of ncRNAs involved in IBs formation.

A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis)

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

Role of p53 in promoting BMP9‑induced osteogenic differentiation of mesenchymal stem cells through TGF‑β1

Known as a tumour suppressor gene, p53 also plays a key role in controlling the differentiation of mesenchymal stem cells (MSCs). Bone morphogenetic protein 9 (BMP9) has been identified as a potent factor in inducing osteogenic differentiation of MSCs, but its relationship with p53 remains unclear. The present study revealed that TP53 was expressed at higher levels in MSCs from patients with osteoporosis and was associated with the top 10 core central genes found in the current osteoporosis genetic screen. p53 was expressed in C2C12, C3H10T1/2, 3T3-L1, MEFs, and MG-63 cell lines, and could be upregulated by BMP9, as measured by western blotting and reverse-transcription quantitative PCR (RT-qPCR). Furthermore, overexpression of p53 increased the mRNA and protein levels of osteogenic marker Runx2 and osteopontin, as evaluated by western blotting and RT-qPCR in BMP9-induced MSCs, whereas the p53 inhibitor pifithrin (PFT)-α attenuated these effects. The same trend was found in alkaline phosphatase activities and matrix mineralization, as measured by alkaline phosphatase staining and alizarin red S staining. Moreover, p53 overexpression reduced adipo-differentiation markers of PPARγ and lipid droplet formation, as measured by western blotting, RT-qPCR and oil red O staining, respectively, whereas PFT-α facilitated adipo-differentiation in MSCs. In addition, p53 promoted TGF-β1 expression and inhibition of TGF-β1 by LY364947 partially attenuated the effects of p53 on promoting BMP9-induced MSC osteo-differentiation and inhibiting adipo-differentiation. The inhibitory effect of PFT-α on osteogenic markers and the promoting effect on adipogenic markers can be reversed when combined with TGF-β1. TGF-β1 may enhance the promotion of osteo-differentiation of MSCs by p53 through inhibition of adipo-differentiation. Collectively, by promoting BMP9-induced MSCs bone differentiation and inhibiting adipose differentiation, p53 may be a novel therapeutic target for bone-related diseases.

Transcriptomic analysis of intermuscular bone development in barbel steed (Hemibarbus labeo)

Intermuscular bones (IBs), which are little, bony spicules in muscle, are embedded in lower teleosts' myosepta. Despite the importance of studying IB development in freshwater aquaculture species, the genes associated with IB development need to be further explored. In the present study, we identified four stages of IB development in barbel steed (Hemibarbus labeo), namely stage 1: IBs have not emerged, stage 2: a few small IBs have emerged in the tail, stage 3: longer IBs gradually emerged in the tail and stage 4: all of the IBs in the tail are mature and long, via Alizarin red staining. Subsequently, we used the HiseqXTen platform to sequence and de novo assemble the transcriptome of epaxial muscle (between 35th and 40th myomere) of barbel steed at 29 days (stage 1) and 42 days (stage 3) after hatching. A total of 190,814 unigenes were obtained with an average length and N50 of 648 bp and 1027 bp, respectively. We found 2174 differentially expressed genes (DEGs) between stages 1 and 3, of which 378 and 1796 were up- and down-regulated, respectively. Functional enrichment analysis showed that several DEGs functioned in ossification, positive regulation of osteoblast differentiation, osteoblast differentiation, and BMP signaling pathway, and were further enriched in signal pathway, including osteoclast differentiation, TGF-β signaling pathway, cytokine-cytokine receptor interaction, Jak-STAT signaling pathway, and other KEEG pathways. In conclusion, we identified genes that may be related to IB development, such as kazal type serine peptidase inhibitor domain 1 (KAZALD1), extracellular matrix protein 1 (ECM1), tetranectin, bone morphogenetic protein 1 (bmp1), acid phosphatase 5 (ACP5), collagen type XI alpha 1 chain (COL11A1), matrix metallopeptidase 9 (MMP9), pannexin-3 (PANX3), sp7 transcription factor (Sp7), and c-x-c motif chemokine ligand 8 (CXCL8), by comparing the transcriptomes of epaxial muscle before and after IB ossification. This study provided a theoretical basis for identifying the molecular mechanisms underlying IB development in fish.

The inclusion of extract from aerial part of Scutellaria baicalensis in feeding of pearl gentian grouper (Epinephelus fuscoguttatus♀ × Epinephelus lanceo-latus♂) promotes growth and immunity

The root of Scutellaria baicalensis (Scutellaria Radix) has been used as herbal medicine for years in China; however, its stem and leaf (aerial part) are considered as waste. The water extract of aerial part of S. baicalensis, named as SBA, having anti-microbial property has been applied in fish aquaculture. To extend the usage of SBA in fish feeding, SBA was employed to feed pearl gentian grouper (a hybrid of Epinephelus fuscoguttatus♀ × Epinephelus lanceolatus♂), and subsequently the total fish output, the levels of digestive enzymes and inflammatory cytokines were determined. Feeding the fish with different doses of SBA for two months, the body length and weight were significantly increased by 5%-10%. In parallel, the expressions of alkaline phosphatase and growth-related factors in bone, liver and muscle of SBA-fed fish were doubled, which could account the growth promoting effect of SBA. Besides, the activity of digestive enzyme, lipase, and the expressions of anti-inflammatory cytokines were markedly stimulated by 2-3 times under the feeding of 3% SBA-containing diet. The results indicated the growth promoting activity of SBA in culture of pearl gentian grouper, as well as the effect of SBA in strengthening the immunity. These beneficial effects of SBA feeding can increase the total yield of pearl gentian grouper in aquaculture. Thus, the re-cycle of waste products during the farming of S. baicalensis herb in serving as fish feeding should be encouraged.

Transcriptomic Analysis Provides Insights to Reveal the bmp6 Function Related to the Development of Intermuscular Bones in Zebrafish

Intermuscular bones (IBs) are small, hard-boned spicules located in the muscle tissue that mainly exist in the myosepta of lower teleosts, which hurt the edibleness and economic value of fish. The study of the development of IBs is very important for freshwater aquaculture fish, but the molecular mechanism of its formation and the key regulatory genes remain unclear. In this study, we first constructed two types of zebrafish mutants (the mutants losing IBs and the mutants with partial deletion of IBs) by knocking out bmp6. We then carried out a transcriptomic analysis to reveal the role of bmp6 in the developmental mechanism of IBs; we used the caudal musculoskeletal tissues of these mutants and wild-type zebrafish at three development stages (20, 45, and 60 dph) to perform transcriptomic analysis. The results showed that the deficiency of bmp6 upregulated sik1 and activated the TNF-A signaling via the NF-KB pathway, which inhibited the development of osteoblasts and promoted osteoclast formation, thereby inhibiting the formation of IBs. These results provided insights to understand the role of bmp6 in the development of IBs in zebrafish and are useful for selective breeding of IBs in cyprinids.

Bulked Segregant Analysis and Association Analysis Identified the Polymorphisms Related to the Intermuscular Bones in Common Carp (Cyprinus carpio)

Simple Summary

Many widely cultured freshwater fish species, such as common carp, belong to the Cyprinidae family. However, most cyprinids have numerous and complex intermuscular bones (IBs), resulting in an adverse effect on cyprinid fish meat processing and consumption. Numerous studies have been trying to understand the development mechanism of IBs and to identify the SNPs associated with the total IB number. However, the SNPs associated with different forms of IBs have been studied less thoroughly. The joint effects of the SNPs on IB development also remain poorly understood. The common carp has numerous geographical populations and domesticated strains, diversifying its phenotypes. The question of whether consensus IB-related SNPs or genes exist among multiple strains of common carp has also not yet been answered. Selective breeding of IB-reduced common carp has been hindered due to a lack of effective molecular markers. To answer these questions, we performed bulked segregant analysis (BSA) to detect the consensus SNPs in three strains. The consensus BSA-SNPs and the other SNPs in their flanking regions were validated in additional individuals. The SNPs associated with the frequency of different IB types were identified. We examined the joint effects of significant SNPs on the numbers of different types of IBs. The identified genetic markers may benefit future selective breeding and reduce the IB number in common carp.

Abstract

The allotetraploid common carp is one of the most important freshwater food fish. However, the IBs found in allotetraploid common carp increase the difficulty in fish meat processing and consumption. Although candidate genes associated with the total IB number have been identified, the SNPs associated with the numbers of the total IBs and different forms of IBs have not yet been identified, hindering the breeding of IB-reduced common carp. Herein, the numbers of different types of IBs in three common carp strains were measured. Using whole-genome resequencing and bulked segregant analysis in three pairs of IB-more and IB-less groups, we identified the consensus nonsynonymous SNPs in three strains of common carp. Screening the flanking regions of these SNPs led to the detection of other SNPs. Association study detected 21 SNPs significantly associated with the number of total IBs, epineural-IBs, and ten detailed types of IBs. We observed the joint effects of multiple SNPs on each associated IB number with an improved explained percentage of phenotypic variation. The resulting dataset provides a resource to understand the molecular mechanisms of IB development in different common carp strains. These SNPs are potential markers for future selection to generate IB-reduced common carp.