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.

Integrative Analysis of GWAS, Eqtls and Meqtls Data Suggests That Multiple Gene Sets Are Associated with Bone Mineral Density

Integration of Summary Data from GWAS and Eqtl Studies Identified Novel Causal BMD Genes with Functional Predictions.

Twelve New Genomic Loci Associated with Bone Mineral Density

Gene-based Association Analysis Identified Novel Genes Associated with Bone Mineral Density.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Joint Association Analysis Identified 18 New Loci For Bone Mineral Density

Integrative Analysis of GWASs, Human Protein Interaction, and Gene Expression Identified Gene Modules Associated with BMDs.

Integrating Epigenomic Elements and GWASs Identifies BDNF Gene Affecting Bone Mineral Density and Osteoporotic Fracture Risk

Meta-Analysis of Genome-Wide Association Studies Identifies Novel Functional CpG-SNPs Associated with Bone Mineral Density at Lumbar Spine.

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

Integrative analysis of genomic and epigenomic data reveal underlying superenhancer-mediated microRNA regulatory network for human bone mineral density

Bivariate whole genome linkage analysis for femoral neck geometric parameters and total body lean mass.

Pleiotropic Loci Underlying Bone Mineral Density and Bone Size Identified by a Bivariate Genome-Wide Association Analysis

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

A Longitudinal Genome-Wide Association Study of Bone Mineral Density Mean and Variability in the UK Biobank

A Multiethnic Whole Genome Sequencing Study to Identify Novel Loci for Bone Mineral Density.

Meta-analysis of Gene-Based Genome-Wide Association Studies of Bone Mineral Density in Chinese and European Subjects

Contribution Of Genotype And Ethnicity To Bone Mineral Density Variation In Caucasians And Chinese: A Test For Five Candidate Genes For Bone Mass

Integrating Genome- Wide Association Study with Regulatory SNP Annotations Identified Novel Candidate Genes for Osteoporosis

Genome-Wide Association And Follow-Up Replication Studies Identified Adamts18 And Tgfbr3 As Bone Mass Candidate Genes In Different Ethnic Groups

Joint Study of Two Genome-Wide Association Meta-Analyses Identified 20p12.1 and 20q13.33 for Bone Mineral Density.