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.

BackgroundBoth obesity (OB) and periodontitis (PD) are chronic non-communicable diseases, and numerous epidemiological studies have demonstrated the association between these two diseases. However, the molecular mechanisms that could explain the association between OB and PD are largely unclear. This study aims to investigate the common gene signatures and biological pathways in OB and PD through bioinformatics analysis of publicly available transcriptome datasets.MethodsThe RNA expression profile datasets of OB (GSE104815) and PD (GSE106090) were used as training data, and GSE152991 and GSE16134 as validation data. After screening for differentially expressed genes (DEGs) shared by OB and PD, gene enrichment analysis, protein-protein interaction (PPI) network construction, GeneMANIA analysis, immune infiltration analysis and gene set enrichment analysis (GSEA) were performed. In addition, receiver operating characteristic (ROC) curves were used to assess the predictive accuracy of the hub gene. Finally, we constructed the hub gene-associated TF-miRNA-mRNA regulatory network.ResultsWe identified a total of 147 DEGs shared by OB and PD (38 down-regulated and 109 up-regulated). Functional analysis showed that these genes were mainly enriched in immune-related pathways such as B cell receptor signalling, leukocyte migration and cellular defence responses. 14 hub genes (FGR, MNDA, NCF2, FYB1, EVI2B, LY86, IGSF6, CTSS, CXCR4, LCK, FCN1, CXCL2, P2RY13, MMP7) showed high sensitivity and specificity in the ROC curve analysis. The results of immune infiltration analysis showed that immune cells such as macrophages, activated CD4 T cells and immune B cells were present at high infiltration levels in both OB and PD samples.The results of GeneMANIA analysis and GSEA analysis suggested that five key genes (FGR, LCK, FYB1, LY86 and P2RY13) may be strongly associated with macrophages. Finally, we constructed a TF-miRNA-mRNA regulatory network consisting of 233 transcription factors (TFs), 8 miRNAs and 14 mRNAs based on the validated information obtained from the database.ConclusionsFive key genes (FGR, LCK, FYB1, LY86, P2RY13) may be important biomarkers of OB and PD. These genes may play an important role in the pathogenesis of OB and PD by affecting macrophage activity and participating in immune regulation and inflammatory responses.

Osteoporosis is a metabolic bone disease that occurs during aging, characterized by low bone mineral density (BMD) and a high risk of trauma fracture. While current pharmacological interventions provide symptomatic benefits, they are unsatisfactory and have major side effects. In this study, we used multi-omics data and drug similarity to construct osteoporosis driver signaling networks (ODSN) and drug functional networks (DFN), respectively. By integrating ODSN and DFN with treatment transcriptional responses, we observed 8 drugs that demonstrated strong targeting effects on ODSN. Mendelian Randomization analysis determines the causal effect on BMD using cis-eQTLs of the drug targets and BMD GWAS data. The findings suggested Acebutolol and Amiloride may increase BMD, while Acenocoumarol, Aminocaproic acid and Armodafinil may enhance bone loss. Zebrafish experiments experimentally showed Acebutolol hydrochloride and Amiloride hydrochloride had significant protective effects on osteoporosis in zebrafish embryos induced by Dexamethasone. Also, Acenocoumarol reduced bone mineralization compared with the control group. The findings suggest that the hypertension drugs Acebutolol and Amiloride warrant further investigation in functional mechanistic experiments to evaluate their effectiveness for osteoporosis treatments.

Background: Lipid pathways have been implicated in the pathogenesis of osteoporosis (OP). Lipid-lowering drugs may be used to prevent and treat OP. However, the causal interpretation of results from traditional observational designs is controversial by confounding. We aimed to investigate the causal association between genetically proxied lipid-lowering drugs and OP risk.Methods: We conducted two-step Mendelian randomization (MR) analyses to investigate the causal association of genetically proxied lipid-lowering drugs on the risk of OP. The first step MR was used to estimate the associations of drug target genes expression with low-density lipoprotein cholesterol (LDL-C) levels. The significant SNPs in the first step MR were used as instrumental variables in the second step MR to estimate the associations of LDL-C levels with forearm bone mineral density (FA-BMD), femoral neck BMD (FN-BMD), lumbar spine BMD (LS-BMD) and fracture. The significant lipid-lowering drugs after MR analyses were further evaluated for their effects on bone mineralization using a dexamethasone-induced OP zebrafish model.Results: The first step MR analysis found that the higher expression of four genes (HMGCR, NPC1L1, PCSK9 and PPARG) was significantly associated with a lower LDL-C level. The genetically decreased LDL-C level mediated by the PPARG was significantly associated with increased FN-BMD (BETA = −1.38, p = 0.001) and LS-BMD (BETA = −2.07, p = 3.35 × 10−5) and was marginally significantly associated with FA-BMD (BETA = −2.36, p = 0.008) and reduced fracture risk (OR = 3.47, p = 0.008). Bezafibrate (BZF) and Fenofibric acid (FBA) act as PPARG agonists. Therefore genetically proxied BZF and FBA had significant protective effects on OP. The dexamethasone-induced OP zebrafish treated with BZF and FBA showed increased bone mineralization area and integrated optical density (IOD) with alizarin red staining.Conclusion: The present study provided evidence that BZF and FBA can increase BMD, suggesting their potential effects in preventing and treating OP. These findings potentially pave the way for future studies that may allow personalized selection of lipid-lowering drugs for those at risk of OP.