Genotype-phenotype relationship in a large cohort of osteogenesis imperfecta patients with COL1A1 mutations revealed by a new scoring system

A hemizygous loss-of-function variant in BCORL1 is associated with male infertility and oligoasthenoteratozoospermia

Oligoasthenoteratozoospermia (OAT) is a common type of male infertility; however, its genetic causes remain largely unknown. Some of the genetic determinants of OAT are gene defects affecting spermatogenesis. BCORL1 (BCL6 corepressor like 1) is a transcriptional corepressor that exhibits the OAT phenotype in a knockout mouse model. A hemizygous missense variant of BCORL1 (c.2615T > G:p.Val872Gly) was reported in an infertile male patient with non-obstructive azoospermia (NOA). Nevertheless, the correlation between BCORL1 variants and OAT in humans remains unknown. In this study, we used whole-exome sequencing to identify a novel hemizygous nonsense variant of BCORL1 (c.1564G > T:p.Glu522*) in a male patient with OAT from a Han Chinese family. Functional analysis showed that the variant produced a truncated protein with altered cellular localization and a dysfunctional interaction with SKP1 (S-phase kinase-associated protein 1). Further population screening identified four BCORL1 missense variants in subjects with both OAT (1 of 325, 0.31%) and NOA (4 of 355, 1.13%), but no pathogenic BCORL1 variants among 362 fertile subjects. In conclusion, our findings indicate that BCORL1 is a potential candidate gene in the pathogenesis of OAT and NOA, expanded its disease spectrum and suggested that BCORL1 may play a role in spermatogenesis by interacting with SKP1.

Osteogenesis imperfecta: a cross-sectional study of skeletal and extraskeletal features in a large cohort of Italian patients

Introduction

The present study aims to describe a large cohort of Italian patients affected by osteogenesis imperfecta, providing a picture of the clinical bony and non-bony features and the molecular background to improve knowledge of the disease to inform appropriate management in clinical practice.

Methods

A total of 568 subjects (from 446 unrelated Italian families) affected by osteogenesis imperfecta who received outpatient care at Istituto Ortopedico Rizzoli from 2006 to 2021 were considered in the present study.

Results

Skeletal and extraskeletal features were analyzed showing a lower height (mean z-scores equal to -1.54 for male patients and -1.47 for female patients) compared with the general Italian population. Half of the patient population showed one or more deformities, and most of the patients had suffered a relatively low number of fractures (<10). An alteration in the sclera color was identified in 447 patients. Similarly, several extraskeletal features, like deafness, dental abnormalities, and cardiac problems, were investigated. Additionally, inheritance and genetic background were evaluated, showing that most of the patients have a positive family history and the majority of pathogenic variants detected were on collagen genes, as per literature.

Conclusion

This study supports the definition of a clear picture of the heterogeneous clinical manifestations leading to variable severity in terms of skeletal and extra-skeletal traits and of the genetic background of an Italian population of osteogenesis imperfecta patients. In this perspective, this clearly highlights the crucial role of standardized and structured collection of high-quality data in disease registries particularly in rare disease scenarios, helping clinicians in disease monitoring and follow-up to improve clinical practice.

Scoliosis in osteogenesis imperfecta: identifying the genetic and non-genetic factors affecting severity and progression from longitudinal data of 290 patients

BACKGROUND

Scoliosis is widely prevalent among osteogenesis imperfecta (OI) patients, and is progressive with age. However, factors affecting scoliosis in OI are not well known.

METHODS

We retrospectively retrieved longitudinal radiographic and clinical records of consecutive OI patients seeking treatments at our hospital from 2014 to 2022, graded their pre-operative spinal conditions into four outcome groups, estimated their progression rates, and descriptively and inferentially analyzed the genetic and non-genetic factors that may affect the outcomes and progression rates.

RESULTS

In all, 290 OI patients met the inclusion criteria, where 221 had genetic records. Of these 221, about 2/3 had mutations in COL1A1 or COL1A2, followed by mutations in WNT1 (9.0%), IFITM5 (9.0%) and other OI risk genes. With an average age of 12.0 years (interquartile range [IQR] 6.9-16.1), 70.7% of the cohort had scoliosis (Cobb angle > 10°), including 106 (36.5%) mild (10°-25°), 40 (13.8%) moderate (25°-50°), and 59 (20.3%) severe (> 50°) scoliosis patients. Patients with either COL1A1 and COL1A2 were strongly biased toward having mild or no scoliosis, whereas patients with mutations in IFITM5, WNT1 and other recessive genes were more evenly distributed among the four outcome grades. Lower-limb discrepancy, bone mineral density (BMD) and age of first drug used were all significantly correlated with severity outcomes. Using multivariate logistic regression, we estimated that each year older adds an odds ratio of 1.13 (95% confidence interval [CI] 1.07-1.2) in progression into advanced stages of scoliosis. We estimated a cohort-wide progression rate of 2.7 degrees per year (95% CI 2.4-3.0). Early-onset patients experienced fast progressions during both infantile and adolescent stages. Twenty-five of the 59 (42.8%) patients with severe scoliosis underwent spinal surgeries, enjoying an average Cobb angle reduction of 33° (IQR 23-40) postoperatively.

CONCLUSION

The severity and progression of scoliosis in osteogenesis imperfecta were affected by genetic factors including genotypes and mutation types, and non-genetic factors including age and BMD. As compared with COL1A1, mutations in COL1A2 were less damaging while those on IFITM5 and other recessive genes conferred damaging effects. Progression rates were the fastest in the adolescent adult age-group.

The patient clinical journey and socioeconomic impact of osteogenesis imperfecta: a systematic scoping review

BACKGROUND

Osteogenesis imperfecta (OI) is a rare heritable connective tissue disorder primarily characterised by skeletal deformity and fragility, and an array of secondary features. The purpose of this review was to capture and quantify the published evidence relating specifically to the clinical, humanistic, and economic impact of OI on individuals, their families, and wider society.

METHODS

A systematic scoping review of 11 databases (MEDLINE, MEDLINE in-progress, EMBASE, CENTRAL, PsycINFO, NHS EED, CEA Registry, PEDE, ScHARRHUd, Orphanet and Google Scholar), supplemented by hand searches of grey literature, was conducted to identify OI literature published 1st January 1995-18th December 2021. Searches were restricted to English language but without geographical limitations. The quality of included records was assessed using the AGREE II checklist and an adapted version of the JBI cross-sectional study checklist.

RESULTS

Of the identified 7,850 records, 271 records of 245 unique studies met the inclusion criteria; overall, 168 included records examined clinical aspects of OI, 67 provided humanistic data, 6 reported on the economic impact of OI, and 30 provided data on mixed outcomes. Bone conditions, anthropometric measurements, oral conditions, diagnostic techniques, use of pharmacotherapy, and physical functioning of adults and children with OI were well described. However, few records included current care practice, diagnosis and monitoring, interactions with the healthcare system, or transition of care across life stages. Limited data on wider health concerns beyond bone health, how these concerns may impact health-related quality of life, in particular that of adult men and other family members, were identified. Few records described fatigue in children or adults. Markedly few records provided data on the socioeconomic impact of OI on patients and their caregivers, and associated costs to healthcare systems, and wider society. Most included records had qualitative limitations.

CONCLUSION

Despite the rarity of OI, the volume of recently published literature highlights the breadth of interest in the OI field from the research community. However, significant data gaps describing the experience of OI for individuals, their families, and wider society warrant further research to capture and quantify the full impact of OI.

Phenotypic Spectrum and Molecular Basis in a Chinese Cohort of Osteogenesis Imperfecta With Mutations in Type I Collagen

Osteogenesis imperfecta (OI) is a rare inherited connective tissue dysplasia characterized with skeletal fragility, recurrent fractures and bone deformity, predominantly caused by mutations in the genes COL1A1 or COL1A2 that encode the chains of type I collagen. In the present study, clinical manifestations and genetic variants were analysed from 187 Chinese OI patients, majority of whom are of southern Chinese origin. By targeted sequencing, 63 and 58 OI patients were found carrying mutations in COL1A1 and COL1A2 respectively, including 8 novel COL1A1 and 7 novel COL1A2 variants. We validated a novel splicing mutation in COL1A1. A diverse mutational and phenotypic spectrum was observed, coupling with the heterogeneity observed in the transcriptomic data derived from osteoblasts of six patients from our cohort. Missense mutations were significantly associated (χ²p = 0.0096) with a cluster of patients with more severe clinical phenotypes. Additionally, the severity of OI was more correlated with the quality of bones, rather than the bone mineral density. Bone density is most responsive to bisphosphonate treatment during the juvenile stage (10–15 years old). In contrast, height is not responsive to bisphosphonate treatment. Our findings expand the mutational spectrum of type I collagen genes and the genotype-phenotype correlation in Chinese OI patients. The observation of effective bisphosphonate treatment in an age-specific manner may help to improve OI patient management.

Cardiovascular abnormalities and its correlation with genotypes of children with osteogenesis imperfecta

BACKGROUND AND OBJECTIVES

Osteogenesis imperfecta (OI) is a rare disorder of abnormal production or modification of type I collagen, which is caused by mutations in COL1A1, COL1A2 or other genes. We investigate the cardiac abnormalities and its correlation with pathogenic mutations in OI children.

METHODS

A cross-sectional comparative study was completed in a relatively large sample of OI children, who were matched in body surface area (BSA) with healthy controls. All echocardiography was performed by experienced cardiologists using Vivid 7 equipment (GE Medical Systems, Horton, Norway). The resting standard 12-lead electrocardiogram (ECG) were obtained in OI patients by FX-8600 machine. Skeletal phenotypes of OI patients were evaluated, including information of bone fractures, deformities, motility, and bone mineral density (BMD). Pathogenic mutations of OI were detected by a next-generation sequencing panel and confirmed by Sanger sequencing.

RESULTS

A total of 69 OI children and 42 healthy children matched in BSA were enrolled. Abnormalities of echocardiography were found in 6 OI children, including enlarged left atrium (n=5), increased internal diameter of the left ventricle (n=1), who all carried the COL1A1 mutation. Mild regurgitation of mitral or tricuspid valves was observed in 26 OI patients. Abnormal ECG manifestations were found in 8 OI children, including deep Q wave, T wave change, premature ventricular complexes, short P-R interval, incomplete bundle branch block and high voltage of left ventricular. Compared with healthy controls, OI children had significant larger values in the main pulmonary artery (1.84 vs 1.60 cm, P < 0.01), left atrial diameter (2.58 vs 2.11 cm, P < 0.001), left ventricular internal dimension at end-diastolic (LVEDd) (3.85 vs 3.50 cm, P < 0.05) and lower left ventricular ejection fraction (LVEF) (68.40% vs 71.74%, P < 0.01). Moreover, OI patients with COL1A1 mutation tended to have greater main pulmonary artery, larger diameters of left atrial and LVEDd, and lower LVEF than healthy controls. COL1A1 mutation was correlated to dilated MPA (β = 1.557, P < 0.01), LAD (β = 3.915, P < 0.001), and LVEDd (β = 2.714, P < 0.01), and decreased LVEF (β = -3.249, P < 0.01).

CONCLUSIONS

Cardiovascular alterations were identified in OI children, including increased dimensions of the main pulmonary artery and left chamber, and low LVEF. The cardiovascular abnormalities seemed to be correlated to COL1A1 mutation and defects of type I collagen, which expanded our understandings of the cardiac phenotypes of OI children.

Comparing Clinical and Genetic Characteristics of De Novo and Inherited COL1A1/COL1A2 Variants in a Large Chinese Cohort of Osteogenesis Imperfecta

PURPOSE

Nearly 85%-90% of osteogenesis imperfecta (OI) cases are caused by autosome dominant mutations of COL1A1 and COL1A2 genes, of which de novo mutations cover a large proportion, whereas their characteristics remain to be elucidated. This study aims to compare the differences in clinical and genetic characteristics of de novo and inherited COL1A1/COL1A2 mutations of OI, assess the average paternal and maternal age at conception in de novo mutations, and research the rate of nonpenetrance in inherited mutations.

MATERIALS AND METHODS

A retrospective comparison between de novo and inherited mutations was performed among 135 OI probands with COL1A1/COL1A2 mutations. Mutational analyses of all probands and their family members were completed by Sanger sequencing. A new clinical scoring system was developed to assess the clinical severity of OI quantitatively.

RESULTS

A total of 51 probands (37.78%) with de novo mutations and 84 probands (62.22%) with inherited mutations were grouped by the results of the parental gene verification. The proportion of clinical type III (P<0.001) and clinical scores (P<0.001) were significantly higher in de novo mutations. Missense mutations covered a slightly higher proportion of de novo COL1A1 mutations (46.34%) compared with inherited COL1A1 mutations (33.33%), however, lacking a significant difference (P=0.1923). The mean BMD Z/T-score at the lumbar spine in de novo mutations was -2.3 ± 1.5, lower than inherited mutations (-1.7 ± 1.8), but lacking statistical significance (P=0.0742). There was no significant difference between the two groups in OI-related phenotypes (like fracture frequency, blue sclera, and hearing loss) and biochemical indexes. In de novo mutations, the average paternal and maternal age at conception was 29.2 (P<0.05) and 26.8 (P<0.0001), respectively, which were significantly younger than the average gestational age of the population. Additionally, 98.04% of pedigrees (50/51) with de novo mutations were spontaneous conception. The rate of nonpenetrance of parents with pathogenic variants in the inherited mutation group was 25.64% (20/78).

CONCLUSIONS

Our data revealed that the proportion of clinical type III and clinical scores were significantly higher in de novo mutations than in inherited mutations, demonstrating that de novo mutations are more damaging because they have not undergone purifying selection.

Total Flavonoids of Rhizoma Drynariae Promotes Differentiation of Osteoblasts and Growth of Bone Graft in Induced Membrane Partly by Activating Wnt/β-Catenin Signaling Pathway

Total flavonoids of Rhizoma drynariae (TFRD), a Chinese medicine, is widely used in the treatment of fracture, bone defect, osteoporosis and other orthopedic diseases, and has achieved good effects. Purpose of this trial was to explore efficacy of TFRD on bone graft’s mineralization and osteoblasts’ differentiation in Masquelet induced membrane technique in rats. Forty male Sprague-Dawley rats were randomly divided into high dose group (H-TFRD), middle dose group (M-TFRD), low dose group (L-TFRD) and control group (control). The critical size bone defect model of rats was established with 10 rats in each group. Polymethyl methacrylate (PMMA) spacer was implanted into the defect of right femur in rats. After the formation of the induced membrane, autogenous bone was implanted into the induced membrane. After 12 weeks of bone graft, bone tissues in the area of bone graft were examined by X-ray, Micro-CT, hematoxylin-eosin (HE) and Masson trichrome staining to evaluate the growth of the bone graft. The β-catenin, c-myc, COL1A1, BMP-2 and OPN in bone graft were quantitatively analyzed by Western blot and Immunohistostaining. Osteoblasts were cultured in the medium containing TFRD. Cell Counting Kit-8 (CCK-8) method, Alkaline phosphatase (ALP) and Alizarin Red S (ARS) staining, Western blot, RT-PCR and other methods were used to detect the effects of TFRD on the proliferation of osteoblasts and the regulation of Wnt/β-catenin signaling pathway. In vivo experiments showed that the growth and mineralization of bone graft in TFRD group was better. Moreover, the expression of Wnt/β-catenin and osteogenesis-related proteins in bone tissue of TFRD group was more than that in other groups. In vitro experiments indicated that osteoblasts proliferated faster, activity of ALP was higher, number of mineralized nodules and proteins related to osteogenesis were more in TFRD group. But blocking Wnt/β-catenin signaling pathway could limit these effects. Therefore, TFRD could promote mineralization of bone graft and differentiation of osteoblasts in a dose-dependent manner during growing period of the bone graft of induced membrane technique, which is partly related to the activation of Wnt/β-catenin signaling pathway.

Genetic analysis in Japanese patients with osteogenesis imperfecta: Genotype and phenotype spectra in 96 probands

BACKGROUND

Osteogenesis imperfecta (OI) is a rare connective-tissue disorder characterized by bone fragility. Approximately 90% of all OI cases are caused by variants in COL1A1 or COL1A2. Additionally, IFITM5 variants are responsible for the unique OI type 5. We previously analyzed COL1A1/2 variants in 22 Japanese families with OI through denaturing high-performance liquid chromatography screening, but our detection rate was low (41%).

METHODS

To expand the genotype-phenotype correlations, we performed a genetic analysis of COL1A1/2 and IFITM5 in 96 non-consanguineous Japanese OI probands by Sanger sequencing.

RESULTS

Of these individuals, 54, 41, and 1 had type 1 (mild), type 2-4 (moderate-to-severe), and type 5 phenotypes, respectively. In the mild group, COL1A1 nonsense and splice-site variants were prevalent (n = 30 and 20, respectively), but there were also COL1A1 and COL1A2 triple-helical glycine substitutions (n = 2 and 1, respectively). In the moderate-to-severe group, although COL1A1 and COL1A2 glycine substitutions were common (n = 14 and 18, respectively), other variants were also detected. The single case of type 5 had the characteristic c.-14C>T variant in IFITM5.

CONCLUSION

These results increase our previous detection rate for COL1A1/2 variants to 99% and provide insight into the genotype-phenotype correlations in OI.

Clinical and genetic analysis in 185 Chinese probands of osteogenesis imperfecta

INTRODUCTION

Osteogenesis imperfecta (OI) is a well-known heritable disorder of connective tissue characterized by skeletal fragility and low bone mass. Nearly 90% of patients with OI have disease variants in COL1A1 and COL1A2 that encode for the α1 and α2 chains of type I collagen.

MATERIALS AND METHODS

A retrospective analysis of 185 probands who were diagnosed with OI in Shanghai Jiao Tong University Affiliated Sixth People's Hospital from March 2005 to December 2019 was performed.

RESULTS

A total of 140 mutations in COL1A1 and 45 mutations in COL1A2 were identified, of which 18 variations were novel. In the phenotype analysis, there were more sporadic cases than familial OI cases in China (54.6% vs. 45.4%, P < 0.001). A total of 98.9% of patients presented with a fracture history. The most common fracture sites were extremity long bones (femur, tibia-fibula and radius-ulna accounted for 36.6%, 17.1% and 11.7%, respectively). Patients with OI types III and IV, especially type III, had a higher proportion of dentinogenesis imperfecta (DI) than patients with OI type I (55% vs. 28%, P < 0.001). Interestingly, G767S and D1219N in COL1A1 and G337S in COL1A2 were the most frequent (3.52%, 2.11% and 8.89%, respectively), which seem to be hotspot mutations in the COL1A1 and COL1A2 genes in Chinese patients.

CONCLUSIONS

This study describes the mutations in the main pathogenic genes, COL1A1 and COL1A2, and the clinical characteristics of osteogenesis imperfecta in China. Furthermore, these findings help reveal the genetic basis of Asian OI patients and contribute to genetic counselling.

Histopathology of osteogenesis imperfecta bone. Supramolecular assessment of cells and matrices in the context of woven and lamellar bone formation using light, polarization and ultrastructural microscopy

Diaphyseal long bone cortical tissue from 30 patients with lethal perinatal Sillence II and progressively deforming Sillence III osteogenesis imperfecta (OI) has been studied at multiple levels of structural resolution. Interpretation in the context of woven to lamellar bone formation by mesenchymal osteoblasts (MOBLs) and surface osteoblasts (SOBLs) respectively demonstrates lamellar on woven bone synthesis as an obligate self-assembly mechanism and bone synthesis following the normal developmental pattern but showing variable delay in maturation caused by structurally abnormal or insufficient amounts of collagen matrix. The more severe the variant of OI is, the greater the persistence of woven bone and the more immature the structural pattern; the pattern shifts to a structurally stronger lamellar arrangement once a threshold accumulation for an adequate scaffold of woven bone has been reached. Woven bone alone characterizes lethal perinatal variants; variable amounts of woven and lamellar bone occur in progressively deforming variants; and lamellar bone increasingly forms rudimentary and then partially compacted osteons not reaching full compaction. At differing levels of microscopic resolution: lamellar bone is characterized by short, obliquely oriented lamellae with a mosaic appearance in progressively deforming forms; polarization defines tissue conformations and localizes initiation of lamellar formation; ultrastructure of bone forming cells shows markedly dilated rough endoplasmic reticulum (RER) and prominent Golgi bodies with disorganized cisternae and swollen dispersed tubules and vesicles, structural indications of storage disorder/stress responses and mitochondrial swelling in cells with massively dilated RER indicating apoptosis; ultrastructural matrix assessments in woven bone show randomly oriented individual fibrils but also short pericellular bundles of parallel oriented fibrils positioned obliquely and oriented randomly to one another and in lamellar bone show unidirectional fibrils that deviate at slight angles to adjacent bundles and obliquely oriented fibril groups consistent with twisted plywood fibril organization. Histomorphometric indices, designed specifically to document woven and lamellar conformations in normal and OI bone, establish ratios for: i) cell area/total area X 100 indicating the percentage of an area occupied by cells (cellularity index) and ii) total area/number of cells (pericellular matrix domains). Woven bone is more cellular than lamellar bone and OI bone is more cellular than normal bone, but these findings occur in a highly specific fashion with values (high to low) encompassing OI woven, normal woven, OI lamellar and normal lamellar conformations. Conversely, for the total area/number of cells ratio, pericellular matrix accumulations in OI woven are smallest and normal lamellar largest. Since genotype-phenotype correlation is not definitive, interposing histologic/structural analysis allowing for a genotype-histopathologic-phenotype correlation will greatly enhance understanding and clinical management of OI.

4-Phenylbutyric acid enhances the mineralization of osteogenesis imperfecta iPSC-derived osteoblasts

Osteogenesis imperfecta (OI) is a heritable brittle bone disease mainly caused by mutations in the two type I collagen genes. Collagen synthesis is a complex process including trimer formation, glycosylation, secretion, extracellular matrix (ECM) formation, and mineralization. Using OI patient-derived fibroblasts and induced pluripotent stem cells (iPSCs), we investigated the effect of 4-phenylbutyric acid (4-PBA) on collagen synthesis to test its potential as a new treatment for OI. Endoplasmic reticulum (ER) retention of type I collagen was observed by immunofluorescence staining in OI patient-derived fibroblasts with glycine substitution and exon skipping mutations. Liquid chromatography-mass spectrometry analysis revealed excessive glycosylation of secreted type I collagen at the specific sites in OI cells. The misfolding of the type I collagen triple helix in the ECM was demonstrated by the incorporation of heat-dissociated collagen hybridizing peptide in OI cells. Type I collagen was produced excessively by OI fibroblasts with a glycine mutation, but this excessive production was normalized when OI fibroblasts were cultured on control fibroblast-derived ECM. We also found that mineralization was impaired in osteoblasts differentiated from OI iPSCs. In summary, treatment with 4-PBA normalizes the excessive production of type I collagen, reduces ER retention, partially improves misfolding of the type I collagen helix in ECM, and improves osteoblast mineralization. Thus, 4-PBA may improve not only ER retention, but also type I collagen synthesis and mineralization in human cells from OI patients.

Fermented Oyster Extract Promotes Osteoblast Differentiation by Activating the Wnt/β-Catenin Signaling Pathway, Leading to Bone Formation

The Pacific oyster, Crassostrea gigas, is well-known as a nutritious food. Recently, we revealed that fermented extract of C. gigas (FO) inhibited ovariectomy-induced osteoporosis, resulting from suppression of osteoclastogenesis. However, since the beneficial effect of FO on osteogenesis is poorly understood, it was examined in mouse preosteoblast MC3T3-E1 cells, human osteosarcoma MG-63 osteoblast-like cells, and zebrafish larvae in this study. We found that FO increased mitochondrial activity from days 1 to 7; however, total cell number of MC3T3-E1 cells gradually decreased without any change in cell viability, which suggests that FO stimulates the differentiation of MC3T3-E1 cells. FO also promoted the expression of osteoblast marker genes, including runt-related transcription factor 2 (mRUNX2), alkaline phosphatase (mALP), collagen type I α1 (mCol1α1), osteocalcin (mOCN), osterix (mOSX), bone morphogenetic protein 2 (mBMP2), and mBMP4 in MC3T3-E1 cells accompanied by a significant increase in ALP activity. FO also increased nuclear translocation of RUNX2 and OSX transcription factors, ALP activity, and calcification in vitro along with the upregulated expression of osteoblast-specific marker proteins such as RUNX2, ALP, Col1α1, OCN, OSX, and BMP4. Additionally, FO enhanced bone mineralization (calcein intensity) in zebrafish larvae at 9 days post-fertilization comparable to that in the β-glycerophosphate (GP)-treated group. All the tested osteoblast marker genes, including zRUNX2a, zRUNX2b, zALP, zCol1a1, zOCN, zBMP2, and zBMP4, were also remarkably upregulated in the zebrafish larvae in response to FO. It also promoted tail fin regeneration in adult zebrafish as same as the GP-treated groups. Furthermore, not only FO positively regulate β-catenin expression and Wnt/β-catenin luciferase activity, but pretreatment with a Wnt/β-catenin inhibitor (FH535) also significantly decreased FO-mediated bone mineralization in zebrafish larvae, which indicates that FO-induced osteogenesis depends on the Wnt/β-catenin pathway. Altogether, the current study suggests that the supplemental intake of FO has a beneficial effect on osteogenesis.