Network-based Transcriptome-wide Expression Study for Postmenopausal Osteoporosis

J Clin Endocrinol Metab. 2020 Aug 1;105(8):2678-2691. doi: 10.1210/clinem/dgaa319.

Abstract

Purpose: Menopause is a crucial physiological transition during a woman's life, and it occurs with growing risks of health issues like osteoporosis. To identify postmenopausal osteoporosis-related genes, we performed transcriptome-wide expression analyses for human peripheral blood monocytes (PBMs) using Affymetrix 1.0 ST arrays in 40 Caucasian postmenopausal women with discordant bone mineral density (BMD) levels.

Methods: We performed multiscale embedded gene coexpression network analysis (MEGENA) to study functionally orchestrating clusters of differentially expressed genes in the form of functional networks. Gene sets net correlations analysis (GSNCA) was applied to assess how the coexpression structure of a predefined gene set differs in high and low BMD groups. Bayesian network (BN) analysis was used to identify important regulation patterns between potential risk genes for osteoporosis. A small interfering ribonucleic acid (siRNA)-based gene silencing in vitro experiment was performed to validate the findings from BN analysis.

Result: MEGENA showed that the "T cell receptor signaling pathway" and the "osteoclast differentiation pathway" were significantly enriched in the identified compact network, which is significantly correlated with BMD variation. GSNCA revealed that the coexpression structure of the "Signaling by TGF-beta receptor complex pathway" is significantly different between the 2 BMD discordant groups; the hub genes in the postmenopausal low and high BMD group are FURIN and SMAD3 respectively. With siRNA in vitro experiments, we confirmed the regulation relationship of TGFBR2-SMAD7 and TGFBR1-SMURF2.

Main conclusion: The present study suggests that biological signals involved in monocyte recruitment, monocyte/macrophage lineage development, osteoclast formation, and osteoclast differentiation might function together in PBMs that contribute to the pathogenesis of postmenopausal osteoporosis.

Keywords: coexpression analysis; osteoclast; postmenopausal osteoporosis; siRNA; transcriptomics.

Publication types

  • Research Support, N.I.H., Extramural

MeSH terms

  • Absorptiometry, Photon
  • Bone Density / genetics*
  • Cohort Studies
  • Computational Biology
  • Female
  • Furin / genetics
  • Gene Expression Profiling
  • Gene Regulatory Networks*
  • Humans
  • Osteoporosis, Postmenopausal / diagnosis
  • Osteoporosis, Postmenopausal / genetics*
  • RNA, Small Interfering / metabolism
  • Receptor, Transforming Growth Factor-beta Type I / genetics
  • Receptor, Transforming Growth Factor-beta Type I / metabolism
  • Receptor, Transforming Growth Factor-beta Type II / genetics
  • Receptor, Transforming Growth Factor-beta Type II / metabolism
  • Signal Transduction / genetics
  • Smad3 Protein / genetics
  • Smad7 Protein / genetics
  • THP-1 Cells
  • Ubiquitin-Protein Ligases / genetics

Substances

  • RNA, Small Interfering
  • SMAD3 protein, human
  • SMAD7 protein, human
  • Smad3 Protein
  • Smad7 Protein
  • SMURF2 protein, human
  • Ubiquitin-Protein Ligases
  • Receptor, Transforming Growth Factor-beta Type I
  • Receptor, Transforming Growth Factor-beta Type II
  • TGFBR1 protein, human
  • TGFBR2 protein, human
  • FURIN protein, human
  • Furin