Commensal Gut Microbiota Immunomodulatory Actions in Bone Marrow and Liver have Catabolic Effects on Skeletal Homeostasis in Health

Novince, Chad M.; Whittow, Carolyn R.; Aartun, Johannes D.; Hathaway, Jessica D.; Poulides, Nicole; Chavez, Michael B.; Steinkamp, Heidi M.; Kirkwood, Kaeleigh A.; Huang, Emily; Westwater, Caroline; Kirkwood, Keith L.

doi:10.1038/s41598-017-06126-x

Download PDF

Article
Open access
Published: 18 July 2017

Commensal Gut Microbiota Immunomodulatory Actions in Bone Marrow and Liver have Catabolic Effects on Skeletal Homeostasis in Health

Chad M. Novince¹,
Carolyn R. Whittow¹,
Johannes D. Aartun¹,
Jessica D. Hathaway¹,
Nicole Poulides¹,
Michael B. Chavez¹,
Heidi M. Steinkamp¹,
Kaeleigh A. Kirkwood¹,
Emily Huang¹,
Caroline Westwater^1,2 &
…
Keith L. Kirkwood^1,2

Scientific Reports volume 7, Article number: 5747 (2017) Cite this article

4580 Accesses
80 Citations
12 Altmetric
Metrics details

Subjects

Abstract

Despite knowledge the gut microbiota regulates bone mass, mechanisms governing the normal gut microbiota’s osteoimmunomodulatory effects on skeletal remodeling and homeostasis are unclear in the healthy adult skeleton. Young adult specific-pathogen-free and germ-free mice were used to delineate the commensal microbiota’s immunoregulatory effects on osteoblastogenesis, osteoclastogenesis, marrow T-cell hematopoiesis, and extra-skeletal endocrine organ function. We report the commensal microbiota has anti-anabolic effects suppressing osteoblastogenesis and pro-catabolic effects enhancing osteoclastogenesis, which drive bone loss in health. Suppression of Sp7(Osterix) and Igf1 in bone, and serum IGF1, in specific-pathogen-free mice suggest the commensal microbiota’s anti-osteoblastic actions are mediated via local disruption of IGF1-signaling. Differences in the RANKL/OPG Axis in vivo, and RANKL-induced maturation of osteoclast-precursors in vitro, indicate the commensal microbiota induces sustained changes in RANKL-mediated osteoclastogenesis. Candidate mechanisms mediating commensal microbiota’s pro-osteoclastic actions include altered marrow effector CD4⁺T-cells and a novel Gut-Liver-Bone Axis. The previously unidentified Gut-Liver-Bone Axis intriguingly implies the normal gut microbiota’s osteoimmunomodulatory actions are partly mediated via immunostimulatory effects in the liver. The molecular underpinnings defining commensal gut microbiota immunomodulatory actions on physiologic bone remodeling are highly relevant in advancing the understanding of normal osteoimmunological processes, having implications for the prevention of skeletal deterioration in health and disease.

Metabolic rewiring promotes anti-inflammatory effects of glucocorticoids

Article 10 April 2024

Human skeletal muscle aging atlas

Article Open access 15 April 2024

Multimodal cell atlas of the ageing human skeletal muscle

Article Open access 22 April 2024

Introduction

Gut microbiota interactions with the host modulates gastrointestinal processes, metabolism and immunity^1,2,3,4,5, having implications for the development and homeostasis of host tissues^{6, 7}. Extensive research has focused on the commensal gut microbiota immunoregulatory effects in the context of resistance to pathogenic microbes and intestinal homeostasis^{8, 9}, and more recently investigations have begun to define the normal gut microbiota’s role in the pathophysiology of metabolic and autoimmune disease states^{4, 6, 9, 10}. Central to this investigation, the commensal gut microbiota’s influence on physiologic tissue remodeling and homeostasis at extra-gastrointestinal sites is largely unknown¹¹.

The study of osteoimmunology has elucidated that innate-immunity, marrow effector T-cells, and diverse endocrine organs regulate osteoclast-osteoblast mediated bone remodeling, both in health and disease^{12,13,14,15,16,17}. Bone remodeling is a continuous dynamic skeletal renewal process in which monocyte-myeloid derived osteoclast cells resorb old bone matrix, and mesenchymal derived osteoblast cells subsequently form new bone matrix. Skeletal homeostasis occurs when balanced osteoclast-osteoblast actions remodel the adult skeleton, and there is no net gain/loss of osseous tissue^{12,13,14, 18}. Understanding that resident gut microbes significantly modulates host immunity^1,2,3,4,5 and metabolic-endocrine processes^{4, 6, 7} in health, the commensal gut microbiota is a strong candidate immunoregulator of physiologic bone remodeling processes in the adult skeleton.

The impact of resident gut microbes on skeletal physiology was introduced by two early life antibiotic treatment studies^{19, 20}, and a growing Conventional (Conv) vs. germ-free (GF) mouse model investigation²¹, which revealed the normal gut microbiota has inhibitory effects on the accrual of bone mass in the post-natal developing skeleton of C57BL/6 mice. Subsequent GF mouse model investigations demonstrating that the normal gut microbiota supports skeletal growth and bone formation in BALB/c²² and CB6F1²³ mice, notably imply that the commensal microbiota impact on skeletal physiology is regulated by mouse strain genetic determinants. While resident gut microbes have been reported to have pro-osteoclastic immunomodulatory effects impairing the accrual of bone mass in growing C57BL/6 mice²¹, and pro-anabolic actions enhancing liver IGF1 mediated skeletal growth in BALB/c²² and CB6F1²³ mice, the normal gut microbiota’s osteoimmunomodulatory effects on skeletal physiology are largely unknown in the healthy adult skeleton.

This osteoimmunology study in young adult C57BL/6 SPF vs. GF mice introduces the commensal gut microbiota as a potent immunoregulator of osteoclast-osteoblast mediated bone remodeling processes in the healthy adult skeleton. 3 month-old mice were investigated since this is the age when bone modeling (growth) is considered principally complete, and physiologic bone remodeling (turnover) is highly metabolically active in the C57BL/6 strain^{24, 25}. Appreciating that the adult skeleton continues to accrue peak bone mineral content until about age 30 years in humans^{26, 27}, optimizing skeletal remodeling processes during young adulthood has significant implications for the prevention of osteoporosis and other aging associated skeletal deterioration states.

While experimental disease studies have elucidated that the commensal gut microbiota drives skeletal deterioration in pathophysiologic states^28,29,30, this report reveals the commensal gut microbiota has pro-osteoclastic and anti-osteoblastic actions which drive catabolic effects on skeletal tissue homeostasis in health. Identification of a novel Gut-Liver-Bone Axis intriguingly implies the commensal gut microbiota’s catabolic actions on physiologic skeletal homeostasis are in part mediated by pro-inflammatory immunostimulatory effects in liver. Having broad systemic implications, this investigation highlights the commensal gut microbiota as a critical regulator of physiologic tissue remodeling and homeostasis at extra-gastrointestinal sites¹¹, a biological phenomenon which is currently unappreciated and not well understood.

Results

Commensal microbiota has anti-anabolic effects on trabecular bone remodeling

Histomorphometric and micro-CT studies were carried out in the trabecular bone compartment, based on the anatomical region being the most metabolically active site of bone remodeling in the young adult bone organ^{31, 32}. Histomorphometric analysis revealed decreased trabecular bone area (B.Ar/T.Ar) in the distal femur of SPF vs. GF mice (Fig. 1b,c). Micro-CT analysis of the proximal tibia demonstrated an osteopenic trabecular bone phenotype in SPF vs. GF mice (Fig. 1e). Consistent with the reduced trabecular bone quantity found in the distal femur via histomorphometry (Fig. 1b), micro-CT analysis revealed decreased trabecular bone volume (BV/TV) in the proximal tibia of SPF mice (Fig. 1f). Inferior trabecular bone micro-architecture properties in the proximal tibia of SPF mice were characterized by reduced trabecular number (Tb.N) (Fig. 1g), and a trend towards increased trabecular separation (Tb.Sp) (Fig. 1i). Micro-CT analysis of cortical bone parameters revealed no differences in femur length, and marginally reduced cortical bone (Ct.Ar/Tt.Ar) in the femur mid-diaphysis of SPF mice (Supplementary Fig. S1).

Recognizing that the normal gut microbiota blunts the accrual of bone mass during bone modeling (growth) processes in the developing C57BL/6 skeleton^19,20,21, dynamic bone formation indices were analyzed (Fig. 1j–l) to validate that alterations in bone remodeling are contributing to the osteopenic trabecular skeletal phenotype found in young adult SPF vs. GF mice. Mineral apposition rate (MAR) (Fig. 1k) and bone formation rate (BFR) (Fig. 1l) were suppressed in the distal femur trabecular bone of SPF mice, which demonstrates that the commensal gut microbiota has anti-anabolic effects on trabecular bone remodeling in the adult skeleton.

Commensal microbiota suppresses osteoblast differentiation and function

Considering the blunted dynamic bone formation indices found in the remodeling trabecular bone of SPF mice (Fig. 1j–l), bone marrow stromal cell (BMSC) osteoblast-progenitors were isolated to elucidate commensal microbiota induced alterations in osteoblast potential (Fig. 2a–g). BMSC expansion over time was blunted in SPF vs. GF BMSC cultures (Fig. 2a), and gene expression studies in untreated day-4 BMSC cultures revealed alterations in intrinsic mesenchymal-stromal cell differentiation potential (Fig. 2b–e). Col2a1 (Fig. 2c), an early marker for chondrogenesis, and Runx2 (Fig. 2d) and Sp7 (Fig. 2e), transcription factors critical for commitment to and maturation in the osteoblastic lineage³³, were significantly decreased in day-4 BMSC cultures from SPF vs. GF mice. Consistent with the reduced osteoblastic differentiation potential observed in untreated day-4 BMSC cultures from SPF mice (Fig. 2d,e), 21-days mineralization treatment (von Kossa assay) induced less mineralization in SPF vs. GF BMSC cultures (Fig. 2f,g).

SPF mice have decreased Bglap and Igf1 in bone, and lower serum IGF1

In light of the blunted mineralization in SPF trabecular bone (Fig. 1k) and BMSC cultures (Fig. 2f,g), Bglap(Osteocalcin) was assessed as a marker for mature osteoblast function in bone marrow and calvaria (Fig. 2h). While bone marrow and calvaria are both heterogeneous in their cellular composition, the rationale for including calvaria is based on its more homogeneous stromal-osteoblastic cellular composition better reflecting osteoblast specific gene expression. Corroborating the suppressed osteoblast function detected in vivo (Fig. 1k) and in vitro (Fig. 2f,g) in SPF vs. GF mice, Bglap was marginally down-regulated in SPF bone marrow, and significantly decreased in SPF calvaria (Fig. 2h).

Considering Sp7(Osterix) was decreased in BMSC cultures from SPF mice (Fig. 2e), and that IGF1 signaling upregulates Sp7 mediated osteoblast maturation and function^{33, 34}, Igf1 expression and serum IGF1 were evaluated (Fig. 2i,j). Igf1 was decreased in SPF vs. GF bone marrow and calvaria (Fig. 2i), which suggest the commensal gut microbiota’s inhibitory effects on osteoblastogenesis are potentially mediated through blunted IGF1 signaling in osteoblastic cells³⁵. Indirect evidence from a transgenic mouse model locally deficient in IGF1 signaling within osteoblastic cells, which has a skeletal phenotype (reduced BV/TV, decreased Tb.N, increased Tb.Sp)³⁵ strikingly similar to SPF vs. GF mice (Fig. 1), indicates the commensal gut microbiota impact on trabecular bone morphology may occur in part by impaired IGF1 signaling within skeletal tissue. Consistent with the reduced Igf1 expression in SPF vs. GF bone tissues (Fig. 2i), serum IGF1 was 18.3% lower in SPF mice (Fig. 2j). Appreciating that liver-derived IGF1 constitutes 70% of circulating IGF1³⁶, differences were ruled out in liver Igf1 expression (Fig. 2i). While the suppressed osteoblastogenesis phenotype in SPF vs GF mice (Fig. 1j–l; Fig. 2) delineates the commensal microbiota’s catabolic effects on skeletal remodeling are mediated in part by blunted osteoblast bone formation, recognizing that bone remodeling occurs through dual osteoclast-osteoblast actions, investigations were carried out to elucidate the commensal microbiota impact on osteoclastogenesis.

Commensal microbiota enhances osteoclast size and eroded bone perimeter

Histomorphometric analysis of tartrate-resistant acid phosphatase (TRAP) stained distal femur sections was performed to investigate the commensal microbiota’s effects on in vivo osteoclastogenesis (Fig. 3a–g). SPF vs. GF mice had similar numbers of osteoclasts lining the trabecular bone perimeter (N.Oc/B.Pm) (Fig. 3b), which suggest the commensal microbiota does not alter the commitment of monocyte-myeloid cells to the osteoclast lineage. The average osteoclast cell size (Oc.Ar/Oc) was 2.5× larger in SPF mice (Fig. 3c), which resulted in a 2× greater osteoclast perimeter per bone perimeter (Oc.Pm/B.Pm) in SPF vs. GF mice (Fig. 3d). The substantially increased Oc.Ar/Oc (Fig. 3c) in SPF mice, implies the commensal microbiota enhances osteoclast maturation.

Eroded bone perimeter analysis (Fig. 3e–g) was performed to assess alterations in osteoclast function. SPF vs. GF mice had an increased eroded perimeter per bone perimeter (E.Pm/B.Pm) (Fig. 3e), demonstrating that the commensal microbiota upregulates bone resorption in health. Findings that the increased osteoclast-positive eroded perimeter per bone perimeter (Oc + E.Pm/B.Pm) (Fig. 3f), paralleled the 2× greater Oc.Pm/B.Pm (Fig. 3d) in SPF vs. GF mice, indicates the commensal microbiota’s pro-resorptive actions are attributed to mechanisms enhancing osteoclast size/maturation (Fig. 3c).

Commensal microbiota modulates the RANKL/OPG Axis

Based on the in vivo histomorphometry findings revealing enhanced osteoclast size/maturation in SPF mice (Fig. 3c), the Tnfsf11(Rankl)/Tnfrsf11b(Opg) Axis was assessed (Fig. 3h–j) to evaluate alterations in critical osteoclastic signaling factors^{13, 15, 18}. RANKL, which signals at the RANK receptor on pre-osteoclast/osteoclast cells, is required for osteoclast differentiation/maturation. Due to OPG functioning as the RANK decoy receptor, it is imperative to assess the RANKL/OPG ratio when evaluating RANKL levels. Gene expression analysis revealed a trend towards a higher Tnfsf11:Tnfrsf11b ratio in marrow, and a significantly upregulated Tnfsf11:Tnfrsf11b ratio in calvaria of SPF vs. GF mice (Fig. 3j). The increased Tnfsf11:Tnfrsf11b ratio findings in SPF mice (Fig. 3j) were interestingly attributed to enhanced Tnfrsf11b expression (Fig. 3i), not alterations in Tnfsf11 expression (Fig. 3h). Recognizing that Tnfrsf11b is primarily expressed by stromal-osteoblastic cells in the bone environment, the findings that Tnfrsf11b upregulation is a trend in bone marrow and significant in calvaria (Fig. 3i) are in-line with calvaria vs. marrow having a more homogeneous stromal-osteoblastic cellular composition. Appreciating that an increased Tnfsf11:Tnfrsf11b ratio indicates higher levels of unbound RANKL available in the bone environment to activate RANK signaling, this notably implies that differences in in vivo RANKL signaling contribute to the superior osteoclast size/maturation phenotype in SPF vs. GF mice (Fig. 3c).

Commensal microbiota dynamically regulates osteoclast-precursor cell maturation

Osteoclast-precursor (OCP) in vitro differentiation assays (Fig. 4a–p) were utilized to further define the commensal microbiota’s regulatory effects on osteoclast maturation, including the role of RANKL-signaling. Isolated marrow hematopoietic progenitor cells (HPCs) were labeled with CD11b MicroBeads, magnetic cell sorting was applied to separate CD11b^neg HPCs, and cells were stimulated in culture (primed with CSF1) to enrich for CD11b^neg osteoclast-precursor (OCP) cells having high osteoclastic potential^37,38,39. CD11b^neg OCP cultures were then subjected to stimulation with control (CSF1 alone) or treatment (CSF1 & RANKL) media for 3, 5 and 7 days^37,38,39. Cultures stimulated for 3, 5 and 7 days were TRAP stained for cytomorphometric analysis of cellular differentiation endpoints (Fig. 4; Supplementary Fig. S2), to evaluate cell level alterations in RANKL-induced osteoclast differentiation. Appreciating that day-3 of the culture system is when OCPs begin fusing to form multi-nucleated osteoclasts^37,38,39, gene expression studies were carried out in day-3 OCP cultures (Fig. 4e–h) to detect early transcription level alterations in RANKL-stimulated osteoclast differentiation.

Cytomorphometric analysis of day-3 OCP cultures showed no differences in SPF vs. GF OCP cellular endpoints (Fig. 4a–d), which implies commensal microbiota immunomodulatory effects do not alter RANKL-induced early commitment of pre-osteoclastic cells to the osteoclast lineage. Day-3 OCP culture gene expression analysis revealed alterations in RANKL-stimulated osteoclastic genes, which unexpectedly indicated enhanced cell level differentiation potential outcomes in GF vs. SPF OCPs at later time points in the culture system. Nfatc1 (Fig. 4e), the master transcription factor for osteoclastogenesis, and Tnfrsf11a (Fig. 4f), the RANKL receptor and a surrogate marker for in vitro osteoclast maturation, were decreased in SPF vs. GF CD11b^neg OCPs. Appreciating that signaling at both the CSF1R and RANK receptor are critical and necessary for osteoclast differentiation, differences in Csf1r (Fig. 4g) were ruled out to delineate that alterations in SPF vs. GF OCP outcomes are mediated by RANKL-stimulation.

Day-5 OCP culture cytomorphometric analysis elucidated cell level alterations in OCP differentiation potential (Fig. 4i–l), which notably paralleled the Nfatc1 (Fig. 4e) and Tnfrsf11a (Fig. 4f) gene expression level differences found in day-3 cultures. Day-5 OCP culture osteoclast size (Oc.Ar/Oc) (Fig. 4k) and numbers of nuclei per osteoclast (N.Nc/Oc) (Fig. 4l) were decreased in SPF vs. GF CD11b^neg OCP cultures. Cytomorphometric analysis of day-7 OCP cultures (Fig. 4m–p), which reflects terminal osteoclast differentiation potential in the OCP culture system, validated the observed superior in vivo osteoclast maturation phenotype noted in distal femur trabecular bone of SPF mice (Fig. 3c). The 2.5× larger osteoclast size (Oc.Ar/Oc) lining the trabecular bone in SPF mice (Fig. 3c), was remarkably paralleled by 2.5× larger Oc.Ar/Oc in SPF vs. GF OCP cultures (Fig. 4o). Number of nuclei per osteoclast (N.Nc/Oc), another indicator of osteoclast maturation, was also greater in day-7 SPF OCP cultures (Fig. 4p).

The study of an enriched OCP population (CD11b^neg) facilitated elucidating subtle alterations in the osteoclastogenic capacity of SPF vs. GF mice, which appears to be secondary to commensal microbiota immunoregulatory effects modulating RANKL-induced early osteoclast fusion genes. Dcstamp, a transmembrane protein critical for osteoclast fusion which is upregulated by RANKL-signaling⁴⁰, was profoundly increased in day-3 GF OCP cultures (Fig. 4h). Appreciating that Dcstamp is essential for RANKL-induced osteoclastogenesis⁴⁰, and considering Dcstamp was substantially increased in CD11b^neg OCPs from GF vs. SPF mice (Fig. 4h), it appears the absence of the commensal microbiota in vivo immunostimulation sensitizes CD11b^neg OCPs to RANKL-stimulated osteoclast differentiation/maturation. The realization that the enhanced osteoclast differentiation endpoints found in the GF vs. SPF OCP cultures at day-5 (Fig. 4k,l) were lost at day-7 (Fig. 4o,p), suggests that the sensitized RANKL-signaling in GF CD11b^neg OCPs was transient.

The altered susceptibility of SPF vs. GF OCPs to RANKL stimulation in the OCP culture system provides novel insight into the commensal microbiota immunomodulatory impact on physiologic osteoclastogenesis. Osteoclasts are critical for normal skeletal turnover and repair processes, which is evidenced by anti-resorptive medications’ deleterious effects (osteonecrosis of the jaw, atypical femoral fractures). Considering that RANKL is locally upregulated at osseous micro-damage⁴¹ and early fracture sites⁴² to transiently stimulate osteoclastic bone resorption necessary for normal osseous repair processes, based on the finding that RANKL-stimulation induced a more rapid onset and shorter duration pro-osteoclastogenic effects in GF vs. SPF OCPs (Fig. 4), the absence of the commensal microbiota may be advantageous for acute skeletal repair. Recognizing that increased chronic systemic inflammation has catabolic effects on systemic skeletal remodeling¹² and disrupts local osseous tissue repair⁴³, the enhanced terminal osteoclast maturation phenotype shown in SPF mice (Fig. 3c; Fig. 4m–p) implies that commensal microbiota pro-osteoclastic immunostimulatory effects are detrimental to skeletal tissue homeostasis in health.

SPF mice have increased pro-inflammatory cytokines in marrow & liver, and higher serum TNF

In order to further elucidate osteoimmunoregulatory mechanisms mediating the commensal gut microbiota’s actions supporting terminal osteoclast maturation (Fig. 3c; Fig. 4o,p), the expression of candidate pro-inflammatory innate-immune cytokines were assessed in bone marrow (Fig. 5a) and at hepato-gastrointestinal sites (Fig. 5b–d). Tnf, a potent pro-inflammatory cytokine which enhances RANKL-mediated osteoclastogenesis and inhibits osteoblastogenesis^{18, 44}, was increased in SPF vs. GF marrow (Fig. 5a).

Due to the proximity of the resident gut microbiota, upregulated cytokine expression was anticipated in the ileum of SPF mice. There was no difference in Tnf, Csf1, Ccl2 or Cxcl1 and surprisingly Il6 was elevated in the GF ileum (Fig. 5b). The increased Il6 expression in the ileum of GF vs. SPF mice (Fig. 5b) is consistent with a prior report in the colon of GF vs. Conv mice²¹, which implies that in absence of exogenous gut microbiota immuno-stimulation, Il6 is endogenously upregulated to support gut homeostasis/function.

Unexpectedly, Tnf, Il6, Csf1, Ccl2 and Cxcl1 were all elevated in the liver (Fig. 5c) of SPF vs. GF mice. These findings are in line with recent seminal reports demonstrating that commensal gut microbiota derived ligands translocate into the circulation under physiologic conditions to directly stimulate the liver innate-immune response in health^{45, 46}. Corroborating the increased innate-immune cytokine expression findings in liver, but not ileum of SPF mice (Fig. 5b,c), the frequency of CD11b⁺LY6G⁻F4/80⁺LY6C^hi (inflammatory monocyte) cells^{47, 48} was enhanced in the draining liver (celiac, portal) lymph nodes (LLNs), but not the mesenteric lymph nodes (MLNs) of SPF mice (Fig. 5e). The upregulated innate-immune cytokine expression in liver and inflammatory monocyte cell frequency in draining LLNs, but no differences in ileum and MLNs of SPF mice, led to the postulation that the commensal gut microbiota has osteoimmunomodulatory actions mediated through a previously unidentified Gut-Liver-Bone Axis (Fig. 5f).

As a means to substantiate the authors’ novel Gut-Liver-Bone axis theory, the expression of pattern-recognition receptors (toll-like receptors, NOD-like receptors) and critical downstream signal transduction factors were evaluated in SPF vs. GF livers (Supplementary Fig. S3). Toll-like receptor (Tlr)2, which primarily recognizes extracellular microbial cell wall components, and Tlr3, which detects viral nucleic acids, were upregulated in the liver of SPF vs. GF mice (Supplementary Fig. S3a). Findings that critical regulators of TLR2 mediated signal transduction (Myd88, Tirap(Mal), Irak4) were increased in SPF livers, while there was no difference in indispensable regulators of TLR3 mediated signal transduction (Ticam1(Trif)) (Supplementary Fig. S3b), indicate that the commensal gut microbiota’s pro-inflammatory actions in the liver are primarily mediated through TLR2 signaling. Notably, this investigation supports findings from prior reports demonstrating that normal gut microbiota derived ligands translocate into the circulation under physiologic conditions to directly stimulate the liver innate-immune response in health^{45, 46}. While pattern-recognition receptor signaling in the liver has been extensively investigated under pathophysiological states, the current report delineates that resident gut microbes substantially upregulate TLR2 signaling in the liver in health.

Congruent with the enhanced Tnf expression in marrow and liver (Fig. 5a,c), serum TNF levels were significantly elevated in SPF mice (Fig. 5g). Serum CSF1 was similar in SPF vs. GF mice (Fig. 5h), which is consistent with knowledge that hepatic/splenic macrophages tightly regulate circulating CSF1 levels via receptor-mediated endocytosis⁴⁹.

T_H17, CD4⁺IL17⁺ and CD4⁺IFNγ⁺ cells are increased in SPF marrow

Despite knowledge that the commensal gut microbiota directs T-lymphocyte mediated immunity, the commensal gut microbiota immuno-regulatory effects on marrow CD4⁺/CD8⁺ T-cell hematopoiesis in the healthy adult skeleton is unclear. Recognizing the field of osteoimmunology has shown that specific bone marrow T-lymphocytic cells regulate osteoclastogenesis and bone remodeling in the adult skeleton^12,13,14,15, flow cytometric analysis was employed to delineate alterations in these CD4⁺/CD8⁺ T-lymphocytic cell populations. While the overall frequency of marrow CD3⁺CD4⁺CD8⁻ (helper) T-cells and CD3⁺CD4⁻CD8⁺ (cytotoxic) T-cells were similar in SPF vs. GF mice (Supplementary Fig. S4), transcription factor expression analysis (Fig. 6a–c) and intracellular cytokine expression analysis (Fig. 6d–g) elucidated differences in the frequency of specific helper CD4⁺T-cell subsets and their characteristic cytokines.

Helper T-cell subsets of interest were based on the expression of characteristic cytokines known to regulate osteoclastogenesis: T_REG cells – IL10 (anti-osteoclastogenic); T_H17 cells – IL17a (pro-osteoclastogenic); T_H1 cells – IFNγ (anti-osteoclastogenic)¹⁵. Whereas there were no differences in the frequency of marrow CD4⁺FOXP3⁺ (T_REG) cells (Fig. 6a) or CD4⁺T-bet⁺ (T_H1) cells (Fig. 6c), the frequency of marrow CD4⁺RORγt⁺ (T_H17) cells (Fig. 6b) was increased SPF vs. GF mice. Consistent with a greater frequency of T_H17 cells (Fig. 6b), intracellular cytokine expression analysis revealed a higher frequency of CD3⁺CD4⁺CD8⁻IL10⁻IL17a⁺/CD3⁺CD4⁺CD8⁻IFNγ⁻IL17a⁺ (CD4⁺IL17a⁺) cells (Fig. 6e,f) in SPF marrow. Considering the lack of alterations in marrow T_H1 cells in SPF vs. GF mice (Fig. 6c), the upregulated frequency of CD3⁺CD4⁺CD8⁻IFNγ⁺IL17a⁻ (CD4⁺IFNγ⁺) cells (Fig. 6g) in SPF marrow is likely attributed to IFNγ expressing non-T_H1 helper T-cells. Intracellular cytokine expression analysis was also carried out in marrow cytotoxic CD8⁺T-cells to rule out alterations in T_C17 and T_C1 cell derived IL17a (CD8⁺IL17a⁺ cells) and IFNγ (CD8⁺IFNγ⁺ cells) (Supplementary Fig. S5).

While a recent experimental postmenopausal osteoporosis study in C57BL/6 SPF vs. GF mice demonstrated the commensal gut microbiota induces the upregulation of marrow TNF, IL17a, and IFNγ in sex steroid deprivation states²⁸, the current study reveals the commensal gut microbiota enhances marrow T-lymphocyte expression of osteoclastogenic cytokines in health. The flow cytometry findings demonstrate the complex nature of lymphocytic cell interactions with bone cells in the marrow environment, providing novel insight into commensal microbiota in vivo immuno-stimulatory effects which could have altered the osteoclast maturation phenotype found in the in vitro OCP system (Fig. 4). Considering that IFNγ suppresses osteoclastogenesis through inhibition of RANKL-signaling via direct targeting of maturing osteoclast precursors^{50, 51}, the suppressed sensitivity to RANKL-stimulation observed in day-5 SPF CD11b^neg OCP cultures (Fig. 4i-l) may be attributed to the enhanced frequency of marrow CD4⁺IFNγ⁺ cells in SPF vs. GF mice (Fig. 6g). Recognizing that IL17a is a pro-resorptive cytokine having potent synergistic effects on TNF pro-osteoclastic actions^{14, 52}, the higher marrow Tnf (Fig. 5a)/circulating TNF (Fig. 5g) levels and increased frequency of marrow CD4⁺IL17a⁺ cells (Fig. 6e,f) in SPF vs. GF mice may mediate the enhanced terminal osteoclast maturation potential exhibited in day-7 SPF OCP cultures (Fig. 4m–p).

Pro-inflammatory effector helper T-cells are increased in LLNs, but not MLNs of SPF mice

Microbial-associated molecular patterns (MAMPs) derived from normal gut microbiota structural components (lipopolysaccharide, peptidoglycan, nucleic acids, etc.) are recognized by host cells via pattern-recognition receptors (toll-like receptors, NOD-like receptors, etc.). Central to the authors’ postulated novel Gut-Bone-Liver Axis (Fig. 5f), normal gut microbiota derived (MAMP) byproducts circulating through liver sinusoidal vessels have been shown to induce innate-immune signaling at liver sinusoid endothelial cells and closely residing tissue resident Kupffer (liver tissue macrophage) and dendritic cells^{45, 46}. Considering tissue resident dendritic cells in gut and liver present antigens in draining lymph nodes to regulate the commitment of naïve T-cells to effector T-cell subsets, T-cell hematopoiesis was assessed in LLNs and MLNs (Fig. 7) to further explore the commensal gut microbiota osteoimmunomodulatory effects mediated through hepato-gastrointestinal tissues.

Consistent with prior reports that lymphocytic cells are reduced in the gastrointestinal tissues of GF vs. SPF mice^{1, 53, 54}, the frequency of overall CD3⁺CD4⁺CD8⁻ (helper) T-cells, CD3⁺CD4⁻CD8⁺ (cytotoxic) T-cells, and CD3⁺CD4⁻CD8^-TCRγδ⁺ (gamma-delta) T-cells were decreased in the MLNs of GF mice (Fig. 7a–c). Transcription factor expression analysis of lymph node cells (Fig. 7d–f) was supplemented with cell surface activation markers, due to experimental design not including intracellular cytokine expression analysis. In line with a previous investigation in SPF vs GF mice⁵⁴, the frequency of CD4⁺CD25⁺FOXP3⁺ (T_REG) cells were upregulated in MLNs from GF mice (Fig. 7d), which may be attributed to the recognition of self-antigens or a hypersensitivity to food antigens in the absence of commensal gut microbiota immuno-stimulation. While there were no differences in the MLNs, the frequencies of CD4⁺CD196⁺RORγt⁺ (T_H17) cells and CD4⁺183⁺T-bet⁺ (T_H1) cells were marginally (Fig. 7e) and significantly (Fig. 7f) increased in the LLNs of SPF vs. GF mice.

The lack of alterations in T_H1/T_H17 cells in SPF MLNs is consistent with current knowledge the MLNs function as the primary site for oral tolerance induction to commensal gut microbiota antigens^{55, 56}. In view of indirect evidence that supra-physiologic T_H1/T_H17 cell priming in MLNs drives pathophysiologic inflammation in gastrointestinal disease states^{56, 57}, the finding that SPF mice have similarly increased effector helper T-cell populations in LLNs suggests that commensal gut microbiota pro-inflammatory actions in health are mediated through the liver. The realization that resident gut microbes induce sustained alterations in the adaptive-immune response in draining LLNs, but not MLNs of SPF mice, supports the authors’ postulate that the commensal gut microbiota has osteoimmunomodulatory actions mediated through a previously unidentified Gut-Liver-Bone Axis.

GF mice colonization with SPF mice gut microbiota normalizes osteoimmunomodulatory outcomes

A conventionalization study was executed to validate that the commensal gut microbiota mediates the osteoimmunomodulatory effects found in the SPF vs. GF mouse model (Figs 1–7). Conventionalized (ConvD) mice were generated by microbially associating 8 week-old GF C57BL/6 littermate mice with gut microbiota from age matched SPF C57BL/6 littermate mice (Fig. 8a). Four weeks following conventionalization there were no differences in trabecular bone (Fig. 8b–f), critical osteoblastic signaling factors (Fig. 8g,h), or critical osteoclastic signaling factors (Fig. 8i,j) in ConvD vs SPF mice. Pro-inflammatory innate-immune cytokines normalized in the marrow (Fig. 8k) and liver (Fig. 8m) of ConvD mice. Findings that Tnf and Cxcl1 were upregulated in the ileum of ConvD vs SPF mice (Fig. 8l) is in line with an acute innate immune response to the recently associated microbiota colonizing the proximal gut tissues, which is likely secondary to a transient increase in neutrophils expressing high levels of TNF and CXCL1.

Discussion

This osteoimmunology study in young adult mice introduces the commensal gut microbiota as a critical immunoregulator of osteoclast-osteoblast mediated bone remodeling in the healthy adult skeleton. As opposed to prior investigations of the normal gut microbiota’s impact on physiological skeletal growth^{19,20,21,22,23} and pathophysiological skeletal deterioration^28,29,30, this report is the first to delineate the commensal gut microbiota’s osteoimmunomodulatory effects on both osteoclastogenesis and osteoblastogenesis in the remodeling healthy adult skeleton.

Earlier investigations of the normal gut microbiota’s influence on skeletal physiology in growing GF mouse models are essentially different than the current investigation of the normal gut microbiota’s impact on bone remodeling in young adult GF mice. To properly appreciate our study outcomes, it is necessary to understand that bone modeling (growth) in the developing skeleton is a fundamentally different skeletal metabolic process than bone remodeling (turnover) in the adult skeleton. Sjogren et al.²¹ reported the normal gut microbiota has inhibitory effects on bone mass accrual in the growing skeleton of C57BL/6 Conv vs. GF mice, which the authors attributed to increased osteoclast numbers secondary to immunostimulatory effects in the gut. In contrast, more recent work showing that the normal gut microbiota supports skeletal growth in BALB/c Conv vs. GF mice²², highlights the potential influence of mouse strain genetic background on gut microbiota osteoimmunoregulatory effects. Investigations in BALB/c Conv vs. GF mice²² and CB6F1 (mixed BALB/c and C57BL/6 background) ConvD vs. GF mice²³ delineated that the normal gut microbiota enhances skeletal growth and bone formation through pro-anabolic actions upregulating liver and serum IGF1. While the current report and prior osteoimmunology investigations in the GF murine model focus on commensal gut microbiota induced host immune response effects, it should be emphasized that confounding study variables (e.g., maternal care effects, environmental factors, diet, dissimilar microbiota across SPF housing facilities, etc.) potentially impact outcomes reported within/across studies.

The current investigation reveals that the normal gut microbiota has catabolic effects on skeletal homeostasis in the healthy adult skeleton, secondary to immunomodulatory actions suppressing osteoblastogenesis and enhancing osteoclastogenesis (Supplementary Fig. S6). SPF mice exhibited blunted bone formation in vivo, and reduced osteoblastic cell differentiation and mineralization in vitro, demonstrating the commensal microbiota has anti-anabolic effects on bone remodeling. Decreased Sp7(Osterix) in and Igf1 in bone tissues and IGF1 in serum from SPF mice, indicate the commensal microbiota’s anti-osteoblastic actions are potentially mediated via suppression of local IGF1-signaling in skeletal tissues. Considering that prior investigations of mice on the BALB/c background have reported the gut microbiota enhances skeletal growth via upregulating liver derived IGF1^{22, 23}, current study findings demonstrating that resident gut microbes suppress skeletal tissue derived IGF1 in C57BL/6 mice highlights IGF1 as a potential genetic determinant modulating the gut microbiota’s osteoimmunomodulatory effects. These findings are in line with prior investigations delineating that mouse strain background significantly impacts skeletal physiology^{58, 59}. Knowledge that BALB/c and C57BL/6 mice differ in polyreactive IgA abundance which impacts the generation of antigen-specific IgA and microbiota diversity⁶⁰, provides indirect evidence that genetic influence on host – microbiota interactions may have implications for the gut microbiota impact on osteoimmunology.

While Sjogren et al.²¹ reported the normal gut microbiota increases osteoclast numbers in growing Conv vs. GF mice, the current study has elucidated the commensal gut microbiota pro-catabolic effects on bone remodeling in the adult skeleton are secondary to immunomodulatory actions enhancing osteoclast size/maturation. Upregulated Tnfsf11:Tnfrsf11b (Rankl:Opg) ratio in vivo, and superior RANKL-stimulated terminal OCP maturation outcomes in vitro, in SPF vs. GF mice delineated that commensal gut microbiota pro-osteoclastic actions are related to sustained alterations in RANKL-signaling. Marrow Tnf, T_H17 cells, and CD4⁺IL17a⁺ T-cells were increased in SPF mice, implying the commensal microbiota’s pro-osteoclastic actions are secondary to immunomodulatory effects directing marrow effector CD4⁺ T-cell hematopoiesis. Unexpected findings that SPF mice had upregulated pro-inflammatory innate-immune cytokines in liver, and increased pro-inflammatory innate- and adaptive-immune cells in LLNs, suggests the commensal gut microbiota’s osteoimmunomodulatory actions are partly mediated by a previously unidentified Gut-Liver-Bone Axis (Supplementary Fig. S6).

Appreciating that all of the toll-like receptors (TLRs) are expressed in the liver, and that TLR2, TLR4, TLR5 and TLR9 recognize the majority of gut microbiota derived ligands in the portal blood^{45, 61}, future research is necessary to determine how signaling at specific TLRs in the liver regulates the normal gut microbiota’s osteoimmunomodulatory effects. Recognizing that distinct commensal gut bacteria disproportionately modulate the host immune response^1,2,3,4,5, ongoing research is indicated to discern whether specific gut commensals drive the upregulated T_H1/T_H17 phenotype found in the LLNs and marrow of SPF mice. Considering indirect evidence that probiotic administration protects against estrogen depletion induced bone loss in experimental murine osteoporosis models^{28, 62, 63}, non-invasive interventions (dietary modulation, pre/probiotics administration) in the gut microbiome need to be researched as a means to optimize skeletal tissue remodeling and homeostasis in health.

Findings from this investigation of the impact of the commensal gut microbiota vs. studies assessing the impact of pathophysiologic gastrointestinal conditions on skeletal metabolism underscore the necessity of understanding health when discerning mechanisms causing disease. Mechanisms that have been extensively reported to drive pathophysiologic bone loss in inflammatory hepato-gastrointestinal disease states, including increased liver CXCL1, CCL2, CSF1, IL6 and TNF, upregulated circulating TNF, decreased circulating IGF1, increased marrow T_H17 cells/IL17a, unbalanced RANKL:OPG ratio, enhanced osteoclastogenesis, and blunted osteoblastogenesis^{44, 64,65,66,67,68}, are astonishingly the same immunomodulatory mechanisms which appear to mediate the commensal gut microbiota’s catabolic effects on skeletal homeostasis in health.

This research defining the commensal gut microbiota immunomodulatory effects on osteoclast-osteoblast mediated skeletal remodeling in health calls into question what is “normal” bone remodeling, a phenomena that is currently poorly understood. Based on the generally accepted empirical theory governing bone remodeling, which dictates that osteoclastic-osteoblastic mediated actions are “coupled” (balanced) in health and “uncoupled” (unbalanced) in disease^{12,13,14, 18}, findings reported here alarmingly imply the commensal gut microbiota induces a pathologic skeletal dys-homeostasis in health. The realization that the commensal gut microbiota has catabolic effects on endogenously programmed bone remodeling infers that bone remodeling is a not a simple “coupled” vs. “uncoupled” process, but rather a physiological gradient. The introduction of the commensal gut microbiota as a critical immunoregulator of normal osteoclast-osteoblast mediated bone remodeling processes in the healthy adult skeleton advances our understanding of skeletal physiology, having significant implications for the prevention of skeletal deterioration in health and disease.

Methods

Animals

Germ-free (GF) C57BL/6 mice were obtained from GF & Gnotobiotic Mouse Facilities at University of Michigan, and bred and maintained in sterile isolators at Medical University of South Carolina (MUSC) Gnotobiotic Animal Core. Specific-pathogen-free (SPF) C57BL/6 mice were purchased from Taconic, and maintained in ventilated cages in a SPF vivarium at MUSC. Mice were aged to 11–12 weeks for experiments. Animal procedures were approved by the MUSC Institutional Animal Care and Use Committee, and carried out in accordance with the approved guidelines.

Conventionalization Study

Conventionalized (ConvD) mice were generated by transferring four 8 week-old GF C57BL/6 littermate mice from the MUSC Gnotobiotic Animal Core to the MUSC SPF vivarium facility. ConvD mice microbial association was performed via fecal inoculum derived from pooled fresh feces of four age matched SPF C57BL/6 littermate mice. Fresh feces from SPF mice was homogenized in PBS, and a fecal inoculum suspension was prepared at a concentration of 0.1 g/ml. At age 8 weeks, and weekly intervals thereafter until age 11 weeks, ConvD mice were administered fresh fecal inoculum (derived from age matched SPF mice) via oral gavage at an approximate dose of 200 uL per mouse; residual fecal inoculum suspension was deposited on the ConvD mice fur. Dirty bedding was transferred from SPF mice caging to ConvD mice caging during weekly inoculation treatments. ConvD and SPF mice were housed in the same SPF vivarium; cages were maintained side by side on the same rack, and on the same shelf during the entire procedure. Further precautions were taken to minimize cage effects, including that animal handling and cage changing/maintenance was carried out by designated personnel trained in gnotobiotic animal husbandry. One week following the final inoculation treatment, ConvD and SPF littermate mice were sacrificed at age 12 weeks.

Histomorphometry

Femurs were fixed in 10% phosphate-buffered-formalin, dehydrated in graded EtOH and xylene, and embedded undecalcified in modified-methylmethacrylate⁶⁹. Serial para-sagittal sections were cut through distal femur for trabecular bone analyses. 4 um tartrate-resistant acid phosphatase (TRAP) stained sections were used for quantifying osteoclast cellular endpoints and assessing eroded bone perimeter. 8 um sections were stained with toluidine blue for bone area analysis. 8 um unstained sections were used for dynamic indices of bone formation. 20 mg/kg calcein was administered via intraperitoneal injection 5 and 2 days prior to euthanasia⁶⁹. Analyses were limited to the secondary spongiosa, beginning 250 µm proximal to the growth plate and extended 1000 µm proximally (50 µm from endocortical surfaces). Data were collected semi-automatically via an Olympus BX61 microscope and Visiopharm software. Data are reported in accordance with standardized nomenclature⁷⁰.

Micro-CT

Tibiae and femurs were fixed in 10% phosphate-buffered-formalin, and stored in 70% EtOH. Specimens were scanned with Scanco Medical µCT 40 Scanner, using the following acquisition parameters: X-ray tube potential = 55 kVp; X-ray intensity = 145 µA; Integration time = 200 ms; Isotropic voxel size = 6 µm³. Calibrated three-dimensional images were reconstructed. Tibia trabecular bone morphology was analyzed using Analyze 12.0 Bone Microarchitecture Analysis software (Analyze Direct). For trabecular analysis, transverse CT slices were analyzed beginning 250 µm distal to the proximal growth plate and extending 1200 µm distally; fixed threshold of 1750 Hounsfield Units was used to discriminate mineralized tissue. Femur length and cortical bone morphology were analyzed using Analyze 12.0 Bone Microarchitecture Analysis software (Analyze Direct). For cortical analysis, transverse CT slices were analyzed in a 1000 µm segment of the mid-diaphysis; fixed threshold of 2500 Hounsfield Units was used to discriminate mineralized tissue. Data are reported in accordance with standardized nomenclature⁷¹.

Bone Marrow Cultures

Femurs and tibiae were harvested and marrow was flushed with α-MEM media, 10% FBS (Hyclone), 1% PSG (2 mM glutamine, 100 U/ml penicillin, 100 mg/ml streptomycin). For each animal, marrow cells were disassociated, counted, and plated at 3 × 10⁶ cells/cm² in a 60 mm dish. 24-hours following plating, hematopoietic progenitor cells (HPCs) were isolated for osteoclast-precursor (OCP) assays via decanting off the non-adherent cells and culture media. Fresh media was added back to the bone marrow cultures, and 48-hours later adherent cells were resourcefully harvested for bone marrow stromal cell (BMSC) assays. Notably, marrow cells were not combined from animals for initial bone marrow cultures or subsequent OCP/BMSC assays; n-values reported for in vitro assays represent biological replicates.

Osteoclast-Precursor (OCP) Assays

First passage HPCs (non-adherent cells isolated from bone marrow cultures) were washed, and incubated with CD11b MicroBeads (BD Biosciences). AutoMACS Sorter (Miltenyi Biotec) was applied to separate CD11b^neg HPCs, which were then stimulated in culture (primed with CSF1) to enrich for CD11b^neg osteoclast-precursor (OCP) cells having high osteoclastic potential^37,38,39.

An initial sort of the HPCs was performed by CD11b positive selection using the Possel separation, which separated the HPCs into CD11b^hi and CD11b^low/neg populations. A subsequent sort of the CD11b^low/neg population was performed by CD11b positive selection using the Possel-s separation, which separated the CD11b^low/neg population into CD11b^low and CD11b^neg populations. CD11b^neg cells were washed, counted, and plated for assays in α-MEM media, 10% FBS (Hyclone), 1% PSG. Gene expression assay: CD11b^neg cells were plated at 1.0 × 10⁵ cells/cm² in 12-well plates, and primed for 36 hours with 10 ng/ml CSF1 (R&D Systems) to drive the cells down the pre-osteoclastic lineage. CD11b^neg OCP cultures were then stimulated with fresh control (25 ng/ml CSF1; R&D Systems) or treatment (25 ng/ml CSF1 and 50 ng/ml RANKL; R&D Systems) media for 3 days; media was refreshed every other day. Day-3 control cultures (25 ng/ml CSF1) and treatment cultures (25 ng/ml CSF1 and 50 ng/ml RANKL) were harvested for quantitative real-time PCR (qRT-PCR) mRNA analysis. Gene expression assay was performed in duplicate (technical replicate) cultures. TRAP stain assay: CD11b^neg cells were plated at 1.5 × 10⁵ cells/cm² in 96-well plates, and primed for 36 hours with 10 ng/ml CSF1 (R&D Systems) to drive the cells down the pre-osteoclastic lineage. CD11b^neg OCP cultures were then stimulated with fresh control (25 ng/ml CSF1; R&D Systems) or treatment (25 ng/ml CSF1 and 50 ng/ml RANKL; R&D Systems) media for 3, 5 and 7 days; media was refreshed every other day. Day-3/5/7 control cultures (25 ng/ml CSF1) and treatment cultures (25 ng/ml CSF1 and 50 ng/ml RANKL) were stained via the TRAP method, as previously reported^{38, 39}. TRAP stain assays was performed in triplicate (technical replicate) cultures; three fields of view analyzed per technical replicate culture. TRAP + cell with > 3 nuclei considered an osteoclast.

Bone Marrow Stromal Cell (BMSC) Assays

First passage BMSCs were isolated from bone marrow cultures; adherent cells were washed, trypsinized, counted, and plated for assays in α-MEM media, 10% FBS (Atlanta Biologicals), 1% PSG. Assays were performed in duplicate (technical replicate) cultures; medium changed every other day. Cell expansion assay: BMSCs plated at 2.0 × 10⁴ cells/cm² in 24-well plates; untreated cultures harvested at day-1/3/5/7/9 for cell counts. Differentiation potential assay: BMSCs plated at 2.0 × 10⁴ cells/cm² in 12-well plates; untreated pre-confluent (day-4) cultures were harvested for qRT-PCR mRNA analysis. Recognizing that cell confluency in the BMSC culture system drives the BMSCs down the osteoblastogenic lineage, cultures were harvested at a pre-confluent state to facilitate evaluating BMSC multipotent differentiation potential down the mesenchymal (osteoblastogenic, adipogenic, chrondrogenic) lineages. Mineralization assay: BMSCs plated at 1.0 × 10⁵ cells/cm² in 12-well plates. Confluent (day-4) cultures were treated with mineralization medium (50 mg/ml ascorbic acid, 10 mM β-glycerophosphate) for 21 days. Mineralization detected by von Kossa method, as reported previously⁶⁹.

Quantitative Real-Time PCR (qRT-PCR) and NanoString

Whole calvaria, spleen, liver, and Ileum were flash frozen, pulverized, and homogenized in TRIzol Reagent (Invitrogen). Femur and tibia bone marrow were flushed with TRIzol. Cultures were washed twice with 1× PBS, and TRIzol was directly applied. RNA was isolated via the TRIzol method, following manufacturer’s protocol. Total RNA was quantified via NanoDrop 1000 (Thermo Scientific). qRT-PCR gene expression analysis: cDNA was synthesized from in vitro cell isolates and in vivo tissue isolates using Taqman Random Hexamers and Reverse Transcription Reagents, according to the manufacturer’s protocol. cDNA was amplified via the StepOnePlus System (Applied Biosystems) protocol, using TaqMan Universal PCR Master Mix and TaqMan gene expression primers/probes. Relative quantification of mRNA was performed via the comparative C _T method (ΔΔCT); Gapdh was utilized as an internal control gene⁷². NanoString gene expression analysis: nCounter PanCancer Immune Profiling for Mouse gene expression panel (NanoString Technologies) was applied to assess pattern-recognition receptor (PRR) signaling in liver. Hybridization was carried out and products were run on the nCounter preparation station, according to the manufacturer’s protocol. Data collected by the nCounter digital analyzer was evaluated via nSolver Analysis Software v2.6 (NanoString Technologies). Data were normalized to the geometric means of spiked-in positive controls and internal control genes. Absolute quantification of RNA reported as normalized RNA counts.

Flow Cytometry

Live Cell Analysis

Mesenteric and liver (celiac, portal) lymph nodes cells were isolated, washed, and counted. Cells were treated with FcR-block (Miltenyi Biotec) and stains were performed. Inflammatory monocytes: anti-CD11b-APC (Miltenyi Biotec, clone M1/70.15.11.5), anti-Ly6g-PacB (Biolegend, clone 18A), anti-F4/80-PE (eBioscience, clone BM8), anti-Ly6c-FITC (Novus Biologicals, clone HK1.4). T-cell hematopoiesis: anti-CD3e-APC (Miltenyi Biotec, clone 145-2C11), anti-TCRγδ-e450 (eBioscience, clone eBioGL3), anti-CD8a-PE (Miltenyi Biotec, clone 53-6.7), anti-CD4-FITC (Miltenyi Biotec, clone GK1.5). Dead cells were excluded via propidium iodide viability dye (Miltenyi Biotec). Data was acquired by the MACSQuant System (Miltenyi Biotec). Analyses were performed via FlowJo VX software (TreeStar).

Transcription Factor Analysis

Femoral bone marrow, mesenteric and liver (celiac, portal) lymph nodes cells were isolated, washed, and counted. Cells were treated with FcR-block (Miltenyi Biotec), surface stains were performed and intracellular stains were carried out following the fixation-permeabilization buffer manufacturer’s protocol (eBioscience). T _REG cells: anti-CD4-FITC (Miltenyi Biotec, clone GK1.5), anti-CD25-PE (eBioscience, clone PC61.5), anti-FOXP3-APC (eBioscience, clone FJK-16s). T _H 17 cells: anti-CD4-FITC (Miltenyi Biotec, clone GK1.5), anti-CD196-PE (Miltenyi Biotec, clone REA277), anti- RORγt-APC (Miltenyi Biotec, clone REA278). T _H 1 cells: anti-CD4-FITC (Miltenyi Biotec, clone GK1.5), anti-CD183-PE (Miltenyi Biotec, clone CXCR3-173), anti-T-bet-APC (Miltenyi Biotec, clone REA102). Dead cells were excluded via e450 viability dye (eBioscience). Data was acquired by the MACSQuant System. Analyses were performed via FlowJo VX software.

Intracellular Cytokine Analysis

Femoral bone marrow cells were isolated and plated, consistent with methods reported for In Vitro Assays. Whole marrow cultures were maintained undisturbed for 4-hours; cultures were subsequently stimulated for 1-hour with 1 µl/ml PMA and 2 µl/ml Ionomycin, followed by 4-hours stimulation with 1 µl/ml Monensin. Cells were harvested, washed, and counted. Surface stains were performed and intracellular stains were carried out following the fixation-permeabilization buffer manufacturer’s protocol (BioLegend): anti-CD3e-FITC (Tonbo Biosciences, clone 145-2C11), anti-CD4-APC/Cy7 (BD Pharmingen, clone GK1.5), anti-CD8a-PerCP/Cy5.5 (BioLegend, clone 53-6.7), anti-IL10-APC (BioLegend, clone JES5-16E3), anti-IL17a-PE (eBioscience, clone eBio17B7), anti-IFNγ-BV421 (BioLegend, clone XMG1.2). Data was acquired by the FACSVerse System (BD Biosciences). Analyses were performed via FlowJo VX software.

Serum Biochemical Assays

Whole blood was collected via cardiac puncture at euthanasia; serum was isolated, and stored at −80 °C. Mouse Quantikine ELISA Immunoassays were performed following manufacturer’s protocol (R&D Systems).

Statistical Analysis

Unpaired t tests were performed using GraphPad Prism 6.0. Data are presented as mean ± SEM; significance is p < 0.05.

References

Round, J. L. & Mazmanian, S. K. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol 9, 313–323 (2009).
Article CAS PubMed PubMed Central Google Scholar
Ivanov, I. I. & Honda, K. Intestinal commensal microbes as immune modulators. Cell Host Microbe 12, 496–508 (2012).
Article CAS PubMed PubMed Central Google Scholar
Belkaid, Y., Bouladoux, N. & Hand, T. W. Effector and memory T cell responses to commensal bacteria. Trends Immunol 34, 299–306 (2013).
Article CAS PubMed PubMed Central Google Scholar
Brestoff, J. R. & Artis, D. Commensal bacteria at the interface of host metabolism and the immune system. Nat Immunol 14, 676–684 (2013).
Article CAS PubMed PubMed Central Google Scholar
Rolig, A. S., Parthasarathy, R., Burns, A. R., Bohannan, B. J. & Guillemin, K. Individual Members of the Microbiota Disproportion-ately Modulate Host Innate Immune Responses. Cell Host Microbe 18, 613–620 (2015).
Article CAS PubMed PubMed Central Google Scholar
Sommer, F. & Backhed, F. The gut microbiota–masters of host development and physiology. Nat Rev Microbiol 11, 227–238 (2013).
Article CAS PubMed Google Scholar
Clemente, J. C., Ursell, L. K., Parfrey, L. W. & Knight, R. The impact of the gut microbiota on human health: an integrative view. Cell 148, 1258–1270 (2012).
Article CAS PubMed PubMed Central Google Scholar
Kamada, N., Chen, G. Y., Inohara, N. & Nunez, G. Control of pathogens and pathobionts by the gut microbiota. Nat Immunol 14, 685–690 (2013).
Article CAS PubMed PubMed Central Google Scholar
Belkaid, Y. & Hand, T. W. Role of the microbiota in immunity and inflammation. Cell 157, 121–141 (2014).
Article CAS PubMed PubMed Central Google Scholar
Kamada, N., Seo, S. U., Chen, G. Y. & Nunez, G. Role of the gut microbiota in immunity and inflammatory disease. Nat Rev Immunol 13, 321–335 (2013).
Article CAS PubMed Google Scholar
Schroeder, B. O. & Backhed, F. Signals from the gut microbiota to distant organs in physiology and disease. Nat Med 22, 1079–1089 (2016).
Article CAS PubMed Google Scholar
Zaidi, M. Skeletal remodeling in health and disease. Nat Med 13, 791–801 (2007).
Article CAS PubMed Google Scholar
Lorenzo, J., Horowitz, M. & Choi, Y. Osteoimmunology: interactions of the bone and immune system. Endocr Rev 29, 403–440 (2008).
Article CAS PubMed PubMed Central Google Scholar
Pacifici, R. T cells: critical bone regulators in health and disease. Bone 47, 461–471 (2010).
Article CAS PubMed PubMed Central Google Scholar
Takayanagi, H. Osteoimmunology and the effects of the immune system on bone. Nat Rev Rheumatol 5, 667–676 (2009).
Article CAS PubMed Google Scholar
Bouillon, R. et al. The past 10 years-new hormones, new functions, new endocrine organs. Nat Rev Endocrinol 11, 681–686 (2015).
Article CAS PubMed Google Scholar
Zhang, Q., Riddle, R. C. & Clemens, T. L. Bone and the regulation of global energy balance. J Intern Med 277, 681–689 (2015).
Article CAS PubMed PubMed Central Google Scholar
Boyce, B. F. & Xing, L. Functions of RANKL/RANK/OPG in bone modeling and remodeling. Arch Biochem Biophys 473, 139–146 (2008).
Article CAS PubMed PubMed Central Google Scholar
Cho, I. et al. Antibiotics in early life alter the murine colonic microbiome and adiposity. Nature 488, 621–626 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
Nobel, Y. R. et al. Metabolic and metagenomic outcomes from early-life pulsed antibiotic treatment. Nat Commun 6, 7486 (2015).
Article PubMed PubMed Central Google Scholar
Sjogren, K. et al. The gut microbiota regulates bone mass in mice. J Bone Miner Res 27, 1357–1367 (2012).
Article PubMed PubMed Central Google Scholar
Schwarzer, M. et al. Lactobacillus plantarum strain maintains growth of infant mice during chronic undernutrition. Science 351, 854–857 (2016).
Article ADS CAS PubMed Google Scholar
Yan, J. et al. Gut microbiota induce IGF-1 and promote bone formation and growth. Proc Natl Acad Sci USA 113, E7554–e7563 (2016).
Article CAS PubMed PubMed Central Google Scholar
Ferguson, V. L., Ayers, R. A., Bateman, T. A. & Simske, S. J. Bone development and age-related bone loss in male C57BL/6J mice. Bone 33, 387–398 (2003).
Article PubMed Google Scholar
Glatt, V., Canalis, E., Stadmeyer, L. & Bouxsein, M. L. Age-related changes in trabecular architecture differ in female and male C57BL/6J mice. J Bone Miner Res 22, 1197–1207 (2007).
Article PubMed Google Scholar
Baxter-Jones, A. D., Faulkner, R. A., Forwood, M. R., Mirwald, R. L. & Bailey, D. A. Bone mineral accrual from 8 to 30 years of age: an estimation of peak bone mass. J Bone Miner Res 26, 1729–1739 (2011).
Article PubMed Google Scholar
Recker, R. R. et al. Bone gain in young adult women. Jama 268, 2403–2408 (1992).
Article CAS PubMed Google Scholar
Li, J. Y. et al. Sex steroid deficiency-associated bone loss is microbiota dependent and prevented by probiotics. J Clin Invest 126, 2049–2063 (2016).
Article PubMed PubMed Central Google Scholar
Tlaskalova-Hogenova, H. et al. Commensal bacteria (normal microflora), mucosal immunity and chronic inflammatory and autoimmune diseases. Immunol Lett 93, 97–108 (2004).
Article CAS PubMed Google Scholar
Wu, H. J. et al. Gut-residing segmented filamentous bacteria drive autoimmune arthritis via T helper 17 cells. Immunity 32, 815–827 (2010).
Article CAS PubMed PubMed Central Google Scholar
Clarke, B. Normal bone anatomy and physiology. Clin J Am Soc Nephrol 3(Suppl 3), S131–139 (2008).
Article CAS PubMed PubMed Central Google Scholar
Farr, J. N. & Khosla, S. Skeletal changes through the lifespan–from growth to senescence. Nat Rev Endocrinol 11, 513–521 (2015).
PubMed PubMed Central Google Scholar
Nakashima, K. et al. The novel zinc finger-containing transcription factor osterix is required for osteoblast differentiation and bone formation. Cell 108, 17–29 (2002).
Article CAS PubMed Google Scholar
Celil, A. B. & Campbell, P. G. BMP-2 and insulin-like growth factor-I mediate Osterix (Osx) expression in human mesenchymal stem cells via the MAPK and protein kinase D signaling pathways. J Biol Chem 280, 31353–31359 (2005).
Article CAS PubMed Google Scholar
Zhang, M. et al. Osteoblast-specific knockout of the insulin-like growth factor (IGF) receptor gene reveals an essential role of IGF signaling in bone matrix mineralization. J Biol Chem 277, 44005–44012 (2002).
Article CAS PubMed Google Scholar
Sjogren, K. et al. Effects of liver-derived insulin-like growth factor I on bone metabolism in mice. J Bone Miner Res 17, 1977–1987 (2002).
Article CAS PubMed Google Scholar
Jacquin, C., Gran, D. E., Lee, S. K., Lorenzo, J. A. & Aguila, H. L. Identification of multiple osteoclast precursor populations in murine bone marrow. J Bone Miner Res 21, 67–77 (2006).
Article PubMed Google Scholar
Valerio, M. S. et al. MKP-1 signaling events are required for early osteoclastogenesis in lineage defined progenitor populations by disrupting RANKL-induced NFATc1 nuclear translocation. Bone 60, 16–25 (2014).
Article CAS PubMed Google Scholar
Herbert, B. A., Valerio, M. S., Gaestel, M. & Kirkwood, K. L. Sexual Dimorphism in MAPK-Activated Protein Kinase-2 (MK2) Regulation of RANKL-Induced Osteoclastogenesis in Osteoclast Progenitor Subpopulations. PLoS One 10, e0125387 (2015).
Article PubMed PubMed Central Google Scholar
Kukita, T. et al. RANKL-induced DC-STAMP is essential for osteoclastogenesis. J Exp Med 200, 941–946 (2004).
Article CAS PubMed PubMed Central Google Scholar
Kennedy, O. D. et al. Activation of resorption in fatigue-loaded bone involves both apoptosis and active pro-osteoclastogenic signaling by distinct osteocyte populations. Bone 50, 1115–1122 (2012).
Article ADS PubMed PubMed Central Google Scholar
Kon, T. et al. Expression of osteoprotegerin, receptor activator of NF-kappaB ligand (osteoprotegerin ligand) and related proinflammatory cytokines during fracture healing. J Bone Miner Res 16, 1004–1014 (2001).
Article CAS PubMed Google Scholar
Claes, L., Recknagel, S. & Ignatius, A. Fracture healing under healthy and inflammatory conditions. Nat Rev Rheumatol 8, 133–143 (2012).
Article CAS PubMed Google Scholar
Redlich, K. & Smolen, J. S. Inflammatory bone loss: pathogenesis and therapeutic intervention. Nat Rev Drug Discov 11, 234–250 (2012).
Article CAS PubMed Google Scholar
Lunz, J. G. 3rd, Specht, S. M., Murase, N., Isse, K. & Demetris, A. J. Gut-derived commensal bacterial products inhibit liver dendritic cell maturation by stimulating hepatic interleukin-6/signal transducer and activator of transcription 3 activity. Hepatology 46, 1946–1959 (2007).
Article CAS PubMed Google Scholar
Corbitt, N. et al. Gut bacteria drive Kupffer cell expansion via MAMP-mediated ICAM-1 induction on sinusoidal endothelium and influence preservation-reperfusion injury after orthotopic liver transplantation. Am J Pathol 182, 180–191 (2013).
Article CAS PubMed PubMed Central Google Scholar
Shi, C. & Pamer, E. G. Monocyte recruitment during infection and inflammation. Nat Rev Immunol 11, 762–774 (2011).
Article CAS PubMed PubMed Central Google Scholar
Grainger, J. R. et al. Inflammatory monocytes regulate pathologic responses to commensals during acute gastrointestinal infection. Nat Med 19, 713–721 (2013).
Article CAS PubMed PubMed Central Google Scholar
Bartocci, A. et al. Macrophages specifically regulate the concentration of their own growth factor in the circulation. Proc Natl Acad Sci USA 84, 6179–6183 (1987).
Article ADS CAS PubMed PubMed Central Google Scholar
Takayanagi, H. et al. T-cell-mediated regulation of osteoclastogenesis by signalling cross-talk between RANKL and IFN-gamma. Nature 408, 600–605 (2000).
Article ADS CAS PubMed Google Scholar
Gao, Y. et al. IFN-gamma stimulates osteoclast formation and bone loss in vivo via antigen-driven T cell activation. J Clin Invest 117, 122–132 (2007).
Article CAS PubMed Google Scholar
Shen, F., Ruddy, M. J., Plamondon, P. & Gaffen, S. L. Cytokines link osteoblasts and inflammation: microarray analysis of interleukin-17- and TNF-alpha-induced genes in bone cells. J Leukoc Biol 77, 388–399 (2005).
Article CAS PubMed Google Scholar
Dobber, R., Hertogh-Huijbregts, A., Rozing, J., Bottomly, K. & Nagelkerken, L. The involvement of the intestinal microflora in the expansion of CD4+ T cells with a naive phenotype in the periphery. Dev Immunol 2, 141–150 (1992).
Article CAS PubMed PubMed Central Google Scholar
Ostman, S., Rask, C., Wold, A. E., Hultkrantz, S. & Telemo, E. Impaired regulatory T cell function in germ-free mice. Eur J Immunol 36, 2336–2346 (2006).
Article PubMed Google Scholar
Worbs, T. et al. Oral tolerance originates in the intestinal immune system and relies on antigen carriage by dendritic cells. J Exp Med 203, 519–527 (2006).
Article CAS PubMed PubMed Central Google Scholar
Koboziev, I., Karlsson, F. & Grisham, M. B. Gut-associated lymphoid tissue, T cell trafficking, and chronic intestinal inflammation. Ann N Y Acad Sci 1207(Suppl 1), E86–93 (2010).
Article ADS PubMed PubMed Central Google Scholar
Sakuraba, A. et al. Th1/Th17 immune response is induced by mesenteric lymph node dendritic cells in Crohn’s disease. Gastroenterology 137, 1736–1745 (2009).
Article CAS PubMed Google Scholar
Rosen, C. J. et al. Mapping quantitative trait loci for serum insulin-like growth factor-1 levels in mice. Bone 27, 521–528 (2000).
Article CAS PubMed Google Scholar
Judex, S., Garman, R., Squire, M., Donahue, L. R. & Rubin, C. Genetically based influences on the site-specific regulation of trabecular and cortical bone morphology. J Bone Miner Res 19, 600–606 (2004).
Article PubMed Google Scholar
Fransen, F. et al. BALB/c and C57BL/6 Mice Differ in Polyreactive IgA Abundance, which Impacts the Generation of Antigen-Specific IgA and Microbiota Diversity. Immunity 43, 527–540 (2015).
Article CAS PubMed Google Scholar
Bigorgne, A.E. & Crispe, I.N. TLRs in Hepatic Cellular Crosstalk. Gastroenterol Res Pract 2010 (2010).
Britton, R. A. et al. Probiotic L. reuteri treatment prevents bone loss in a menopausal ovariectomized mouse model. J Cell Physiol 229, 1822–1830 (2014).
Article CAS PubMed PubMed Central Google Scholar
Ohlsson, C. et al. Probiotics protect mice from ovariectomy-induced cortical bone loss. PLoS One 9, e92368 (2014).
Article ADS PubMed PubMed Central Google Scholar
Tilg, H., Moschen, A. R., Kaser, A., Pines, A. & Dotan, I. Gut, inflammation and osteoporosis: basic and clinical concepts. Gut 57, 684–694 (2008).
Article CAS PubMed Google Scholar
Nakchbandi, I. A. & van der Merwe, S. W. Current understanding of osteoporosis associated with liver disease. Nat Rev Gastroenterol Hepatol 6, 660–670 (2009).
Article PubMed Google Scholar
Ciucci, T. et al. Bone marrow Th17 TNFalpha cells induce osteoclast differentiation, and link bone destruction to IBD. Gut 64, 1072–1081 (2015).
Article CAS PubMed Google Scholar
Heymann, F. & Tacke, F. Immunology in the liver - from homeostasis to disease. Nat Rev Gastroenterol Hepatol 13, 88–110 (2016).
Article CAS PubMed Google Scholar
Harris, L., Senagore, P., Young, V. B. & McCabe, L. R. Inflammatory bowel disease causes reversible suppression of osteoblast and chondrocyte function in mice. Am J Physiol Gastrointest Liver Physiol 296, G1020–1029 (2009).
Article CAS PubMed PubMed Central Google Scholar
Novince, C. M. et al. Proteoglycan 4: a dynamic regulator of skeletogenesis and parathyroid hormone skeletal anabolism. J Bone Miner Res 27, 11–25 (2012).
Article CAS PubMed PubMed Central Google Scholar
Dempster, D. W. et al. Standardized nomenclature, symbols, and units for bone histomorphometry: a 2012 update of the report of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res 28, 2–17 (2013).
Article PubMed PubMed Central Google Scholar
Bouxsein, M. L. et al. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J Bone Miner Res 25, 1468–1486 (2010).
Article PubMed Google Scholar
Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 3, 1101–1108 (2008).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

This work was supported by the National Institutes of Health: K08DE025337 (CMN), R01DE021423 (KLK), P30GM103331 (KLK), T32DE017551 (KLK), R25DE022677 (KLK); UL1TR000062 (CMN). The authors thank Abigail Lauer from the Biostatistical Unit of the MUSC Center for Oral Health Research Clinical Core for statistical consultation. The authors thank Courtney LaCosta for proof-reading assistance.

Author information

Authors and Affiliations

Department of Oral Health Sciences and Center for Oral Health Research, College of Dental Medicine, Medical University of South Carolina, Charleston, South Carolina, 29425, USA
Chad M. Novince, Carolyn R. Whittow, Johannes D. Aartun, Jessica D. Hathaway, Nicole Poulides, Michael B. Chavez, Heidi M. Steinkamp, Kaeleigh A. Kirkwood, Emily Huang, Caroline Westwater & Keith L. Kirkwood
Department of Microbiology and Immunology, Hollings Cancer Center, Medical University of South Carolina, Charleston, South Carolina, 29425, USA
Caroline Westwater & Keith L. Kirkwood

Authors

Chad M. Novince
View author publications
You can also search for this author in PubMed Google Scholar
Carolyn R. Whittow
View author publications
You can also search for this author in PubMed Google Scholar
Johannes D. Aartun
View author publications
You can also search for this author in PubMed Google Scholar
Jessica D. Hathaway
View author publications
You can also search for this author in PubMed Google Scholar
Nicole Poulides
View author publications
You can also search for this author in PubMed Google Scholar
Michael B. Chavez
View author publications
You can also search for this author in PubMed Google Scholar
Heidi M. Steinkamp
View author publications
You can also search for this author in PubMed Google Scholar
Kaeleigh A. Kirkwood
View author publications
You can also search for this author in PubMed Google Scholar
Emily Huang
View author publications
You can also search for this author in PubMed Google Scholar
Caroline Westwater
View author publications
You can also search for this author in PubMed Google Scholar
Keith L. Kirkwood
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

C.M.N., C.W. and K.L.K. designed research; C.M.N., C.R.W., J.D.A., J.D.H., N.P., M.B.C., H.M.S., K.A.K., E.H. and C.W. performed research; C.M.N., C.R.W., J.D.A., J.D.H., N.P., K.A.K. and E.H. analyzed data; C.M.N. and K.L.K. wrote manuscript; all authors reviewed the manuscript.

Corresponding author

Correspondence to Chad M. Novince.

Ethics declarations

Competing Interests

The authors declare that they have no competing interests.

Additional information

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Novince, C.M., Whittow, C.R., Aartun, J.D. et al. Commensal Gut Microbiota Immunomodulatory Actions in Bone Marrow and Liver have Catabolic Effects on Skeletal Homeostasis in Health. Sci Rep 7, 5747 (2017). https://doi.org/10.1038/s41598-017-06126-x

Download citation

Received: 04 October 2016
Accepted: 08 June 2017
Published: 18 July 2017
DOI: https://doi.org/10.1038/s41598-017-06126-x

This article is cited by

Modulation of bone remodeling by the gut microbiota: a new therapy for osteoporosis
- Zhengtian Lyu
- Yongfei Hu
- Dan Liu
Bone Research (2023)
Commensal gut bacterium critically regulates alveolar bone homeostasis
- Jessica D. Hathaway-Schrader
- Matthew D. Carson
- Chad M. Novince
Laboratory Investigation (2022)
Phages from Ganges River curtail in vitro biofilms and planktonic growth of drug resistant Klebsiella pneumoniae in a zebrafish infection model
- Niranjana Sri Sundaramoorthy
- Subramaniam Thothathri
- Saisubramanian Nagarajan
AMB Express (2021)
The Gut Microbiome: a New Frontier in Musculoskeletal Research
- R. Li
- C. G. Boer
- Carolina Medina-Gomez
Current Osteoporosis Reports (2021)
Gut Microbiota in Bone Health and Diabetes
- Julie Kristine Knudsen
- Peter Leutscher
- Suzette Sørensen
Current Osteoporosis Reports (2021)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.