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

• Published in: The journal of clinical endocrinology and metabolism 2020-12, Vol.105 (12), p.1-e4757
• Main Authors: Liu, Yong, Shen, Hui, Greenbaum, Jonathan, Liu, Anqi, Su, Kuan-Jui, Zhang, Li-Shu, Zhang, Lei, Tian, Qing, Hu, Hong-Gang, He, Jin-Sheng, Deng, Hong-Wen
• Format: Article
• Language: English
• Subjects: Adolescent; Adult; Bone mineral density; Child; Child, Preschool; Clinical s; Female; Fractures; Gene Expression; Gene Expression Profiling; Gene Regulatory Networks; Genetic aspects; Genetic Association Studies; Genetic diversity; Genome-wide association studies; Genomes; Humans; Infant; Infant, Newborn; Male; Middle Aged; Osteoporosis; Osteoporosis - genetics; Physiological aspects; Quantitative trait loci; Ribonucleic acid; Risk factors; RNA; RNA splicing; RNA Splicing - genetics; Single-nucleotide polymorphism; Splicing; Transcriptome; Young Adult
• ISSN: 0021-972X; 1945-7197
• DOI: 10.1210/clinem/dgaa572

Abstract Context Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms. Objective Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis. Design, Setting, and Participants Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings. Results In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1. Conclusions Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.