Abstract:In this study, we integrated large-scale GWAS summary data and used the predicted transcriptome-wide association study method to discover novel genes associated with osteoporosis. We identified 204 candidate genes, which provide novel clues for understanding the genetic mechanism of osteoporosis and indicate potential therapeutic targets. Introduction: Osteoporosis is a highly polygenetic disease characterized by low bone mass and deterioration of the bone microarchitecture. Our objective was to discover novel candidate genes associated with osteoporosis. Methods: To identify potential causal genes of the associated loci, we investigated trait-gene expression associations using the transcriptome-wide association study (TWAS) method. This method directly imputes gene expression effects from genome-wide association study (GWAS) data using a statistical prediction model trained on GTEx reference transcriptome data. We then performed a colocalization analysis to evaluate the posterior probability of biological patterns: associations characterized by a single causal variant or multiple distinct causal variants. Finally, a functional enrichment analysis of gene sets was performed using the VarElect and CluePedia tools, which assess the causal relationships between genes and a disease and search for potential gene's functional pathways. The osteoporosis-associated genes were further confirmed based on the differentially expressed genes profiled from mRNA expression data of bone tissue. Results: Our analysis identified 204 candidate genes, including 154 genes that have been previously associated with osteoporosis, 50 genes that have not been previously discovered. A biological function analysis found that 20 of the candidate genes were directly associated with osteoporosis. Further analysis of multiple gene expression profiles showed that 15 genes were differentially expressed in patients with osteoporosis. Among these, SLC11A2, MAP2K5, NFATC4, and HSP90B1 were enriched in four pathways, namely, mineral absorption pathway, MAPK signaling pathway, Wnt signaling pathway, and PI3K-Akt signaling pathway, which indicates a causal relationship with the occurrence of osteoporosis. Conclusions: We demonstrated that transcriptome fine-mapping identifies more osteoporosis-related genes and provides key insight into the development of novel targeted therapeutics for the treatment of osteoporosis.

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

Integrating Transcriptome-Wide Association Study and Mrna Expression Profiling Identifies Novel Genes Associated with Bone Mineral Density

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

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

Integrative Genomic Analysis Predicts Novel Functional Enhancer-Snps for Bone Mineral Density

Integrating genome-wide association and transcriptome prediction model identifies novel target genes for osteoporosis

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

Are Bone Mineral Density Loci Associated with Hip Osteoporotic Fractures? A Validation Study on Previously Reported Genome-Wide Association Loci in a Chinese Population

Targeted Next Generation Sequencing of Nine Osteoporosis-Related Genes in the Wnt Signaling Pathway among Chinese Postmenopausal Women

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

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

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

Replication Study of Candidate Genes/loci Associated with Osteoporosis Based on Genome-Wide Screening

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

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

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

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

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

Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis.

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