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.

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

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

Transcription Factor Enrichment Analysis in Enhancers Identifies EZH2 As a Susceptibility Gene for Osteoporosis.

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

An Osteoporosis Risk SNP at 1p36.12 Acts As an Allele-Specific Enhancer to Modulate LINC00339 Expression Via Long-Range Loop Formation

Author Response for "an Osteoporosis Susceptibility Allele at 11p15 Regulates SOX6 Expression by Modulating TCF4 Chromatin Binding"

The Osteoporosis Susceptible SNP Rs4325274 Remotely Regulates the SOX6 Gene Through Enhancers.

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

The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation

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

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

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

LncRNA ZBTB40-IT1 Modulated by Osteoporosis GWAS Risk SNPs Suppresses Osteogenesis

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

Transcriptomic Data Identified Key Transcription Factors for Osteoporosis in Caucasian Women

Genetic Association Study Identified a 20 Kb Regulatory Element in WLS Associated with Osteoporosis and Bone Mineral Density in Han Chinese

USF3 Modulates Osteoporosis Risk by Targeting WNT16, RANKL, RUNX2, and Two GWAS Lead SNPs Rs2908007 and Rs4531631

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

Endochondral Ossification Pathway Genes and Postmenopausal Osteoporosis: Association and Specific Allele Related Serum Bone Sialoprotein Levels in Han Chinese.

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