Abstract:BackgroundExisting microarray studies of bone mineral density (BMD) have been critical for understanding the pathophysiology of osteoporosis, and have identified a number of candidate genes. However, these studies were limited by their relatively small sample sizes and were usually analyzed individually. Here, we propose a novel network-based meta-analysis approach that combines data across six microarray studies to identify functional modules from human protein-protein interaction (PPI) data, and highlight several differentially expressed genes (DEGs) and a functional module that may play an important role in BMD regulation in women.MethodsExpression profiling studies were identified by searching PubMed, Gene Expression Omnibus (GEO) and ArrayExpress. Two meta-analysis methods were applied across different gene expression profiling studies. The first, a nonparametric Fisher's method, combined p-values from individual experiments to identify genes with large effect sizes. The second method combined effect sizes from individual datasets into a meta-effect size to gain a higher precision of effect size estimation across all datasets. Genes with Q test's p-values <0.05 or I-2 values >50% were assessed by a random effects model and the remainder by a fixed effects model. Using Fisher's combined p-values, functional modules were identified through an integrated analysis of microarray data in the context of large protein-protein interaction (PPI) networks. Two previously published meta-analysis studies of genomewide association (GWA) datasets were used to determine whether these module genes were genetically associated with BMD. Pathway enrichment analysis was performed with a hypergeometric test.ResultsSix gene expression datasets were identified, which included a total of 249 (129 high BMD and 120 low BMD) female subjects. Using a network-based meta-analysis, a consensus module containing 58 genes (nodes) and 83 edges was detected. Pathway enrichment analysis of the 58 module genes revealed that these genes were enriched in several important KEGG pathways including Osteoclast differentiation, B cell receptor signaling pathway, MAPK signaling pathway, Chemokine signaling pathway and Insulin signaling pathway. The importance of module genes was replicated by demonstrating that most module genes were genetically associated with BMD in the GWAS data sets. Meta-analyses were performed at the individual gene level by combining p-values and effect sizes. Five candidate genes (ESR1, MAP3K3, PYGM, RAC1 and SYK) were identified based on gene expression meta-analysis, and their associations with BMD were also replicated by two BMD metaanalysis studies.ConclusionsIn summary, our network-based meta-analysis not only identified important differentially expressed genes but also discovered biologically meaningful functional modules for BMD determination. Our study may provide novel therapeutic targets for osteoporosis in women.

A novel approach for correction of crosstalk effects in pathway analysis and its application in osteoporosis research

Integrative Analysis of Transcriptomic and Epigenomic Data to Reveal Regulation Patterns for Bmd Variation

Attenuated monocyte apoptosis, a new mechanism for osteoporosis suggested by a transcriptome-wide expression study of monocytes.

Screening pathways and hub genes involved in osteoclastogenesis by gene expression analysis of circulating monocytes based on Gibbs sampling

Monocytes Affect Bone Mineral Density in Pre- and Postmenopausal Women Through Ribonucleoprotein Complex Biogenesis by Integrative Bioinformatics Analysis.

A Bootstrap-Based Regression Method for Comprehensive Discovery of Differential Gene Expressions: an Application to the Osteoporosis Study

Multi-omics Data Integration for Identifying Osteoporosis Biomarkers and Their Biological Interaction and Causal Mechanisms

Analysis and Validation of Hub Genes in Blood Monocytes of Postmenopausal Osteoporosis Patients

Network-based Transcriptome-wide Expression Study for Postmenopausal Osteoporosis

Transcriptomic Data Identified Key Transcription Factors for Osteoporosis in Caucasian Women

In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis.

Study on Key Genes and Regulatory Networks Associated with Osteoporosis by Microarray Technology

Unraveling the Unexplored Complexity of Osteoporosis: Deciphering Targeted Therapeutic Strategies Through Integrated Bioinformatics Analysis, Mendelian Randomization, and Drug Interactions

Integrative Analysis of Genomics and Transcriptome Data to Identify Regulation Networks in Female Osteoporosis

Network-Based Meta-Analyses of Associations of Multiple Gene Expression Profiles with Bone Mineral Density Variations in Women

Screening of Important Markers in Peripheral Blood Mononuclear Cells to Predict Female Osteoporosis Risk Using LASSO Regression Algorithm and SVM Method

Identification of differentially expressed genes, signaling pathways and immune infiltration in postmenopausal osteoporosis by integrated bioinformatics analysis

Investigating the impact of Wnt pathway-related genes on biomarker and diagnostic model development for osteoporosis in postmenopausal females

RETRACTED ARTICLE: Key Genes Associated with Osteoporosis Revealed by Genome Wide Gene Expression Analysis

Integrative analyses of genetic mechanisms responsible for bone–fat imbalance in osteoporosis

Network-Based Proteomic Analysis for Postmenopausal Osteoporosis in Caucasian Females