Analysis of the causal relationship between gut microbiota and bone remodeling growth factor from the gene association

Background: A growing body of evidence indicates a close association between the gut microbiota (GM) and the bone remodeling (BR) process, raising suspicions that the GM may actively participate in BR by modulating the levels of growth factors. However, the precise causal relationship between them remains unclear. Due to many confounding factors, many microorganisms related to BR growth factors have not been identified. We aimed to elucidate the causal relationship between the GM and BR growth factors.

Methods: We evaluated the genome-wide association study (GWAS) summary statistics for GM and five common growth factors associated with BR: namely, bone morphogenetic proteins (BMP), transforming growth factors(TGF), insulin growth factors (IGFs), epidermal growth factors (EGFs), and fibroblast growth factors (FGF). The causal relationship between the GM and BR growth factors was studied by double-sample Mendelian randomized analysis. We used five Mendelian randomization (MR) methods, including inverse variance-weighted (IVW), MR-Egger, simple mode, weighted median, and weighted model methods.

Results: Through MR analysis, a total of 56 bacterial genera were co-identified as associated with BMP, TGF, IGF, EGF, and FGF. Among them, eight genera were found to have a causal relationship with multiple growth factors: Marvinbryantia was causally associated with BMP-6 (P = 0.018, OR = 1.355) and TGF-β2 (P = 0.002, OR = 1.475); Lachnoclostridium, BMP-7 (P = 0.021, OR = 0.73) and IGF-1 (P = 0.046, OR = 0.804); Terrisporobacter, TGF-β (P = 0.02, OR = 1.726) and FGF-23 levels (P = 0.016, OR = 1.76); Ruminiclostridium5, TGF-β levels (P = 0.024, OR = 0.525) and FGFR-2 (P = 0.003, OR = 0.681); Erysipelatoclostridium, TGF-β2 (P = 0.001, OR = 0.739) and EGF and its receptor (EGFR) (P = 0.012, OR = 0.795); Eubacterium_brachy_group, FGFR-2 (P = 0.045, OR = 1.153) and EGF (P = 0.013, OR = 0.7); Prevotella9 with EGFR (P = 0.022, OR = 0.818) and FGFR-2 (P = 0.011, OR = 1.233) and Faecalibacterium with FGF-23 (P = 0.02, OR = 2.053) and IGF-1 (P = 0.005, OR = 0.843).

Conclusion: We confirmed the causal relationship between the GM and growth factors related to BR, which provides a new perspective for the study of BR, through targeted regulation of specific bacteria to prevent and treat diseases and growth factor-mediated BR disorders.