Unveiling the transcriptomic landscape and the potential antagonist feedback mechanisms of TGF-β superfamily signaling module in bone and osteoporosis

Abstract

Abstract Background TGF-β superfamily signaling is indispensable for bone homeostasis. However, the global expression profiles of all the genes that make up this signaling module in bone and bone-related diseases have not yet been well characterized. Methods Transcriptomic datasets from human bone marrows, bone marrow-derived mesenchymal stem cells (MSCs) and MSCs of primary osteoporotic patients were used for expression profile analyses. Protein treatments, gene quantification, reporter assay and signaling dissection in MSC lines were used to clarify the interactive regulations and feedback mechanisms between TGF-β superfamily ligands and antagonists. Ingenuity Pathway Analysis was used for network construction. Results We identified TGFB1 in the ligand group that carries out SMAD2/3 signaling and BMP8A, BMP8B and BMP2 in the ligand group that conducts SMAD1/5/8 signaling have relatively high expression levels in normal bone marrows and MSCs. Among 16 antagonist genes, the dominantly expressed TGF-β superfamily ligands induced only NOG, GREM1 and GREM2 via different SMAD pathways in MSCs. These induced antagonist proteins further showed distinct antagonisms to the treated ligands and thus would make up complicated negative feedback networks in bone. We further identified TGF-β superfamily signaling is enriched in MSCs of primary osteoporosis. Enhanced expression of the genes mediating TGF-β-mediated SMAD3 signaling and the genes encoding TGF-β superfamily antagonists served as significant features to osteoporosis. Conclusion Our data for the first time unveiled the transcription landscape of all the genes that make up TGF-β superfamily signaling module in bone. The feedback mechanisms and regulatory network prediction of antagonists provided novel hints to treat osteoporosis.