Abstract:Aiming to identify more genomic loci associated with bone mineral density (BMD), we conducted a joint association analysis of 2 genome-wide association study (GWAS) by the integrative association method multi-trait analysis of GWAS (MTAG). The first one is the single GWAS of estimated heel BMD (eBMD) in the UK biobank (UKB) cohort (N = 426,824), and the second one is the GWAS meta-analysis of total body BMD (TB-BMD) in 66,628 participants from 30 studies. Approximate conditional association analysis was performed in the identified novel loci to identify secondary association signal. Statistical fine-mapping was conducted to prioritize plausible credible risk variants (CRVs). Candidate genes were prioritized based on the analyses of cis- expression quantitative trait locus (cis-eQTL) and cis-protein QTL (cis-pQTL) information as well as the functional category of the SNP. By integrating the information carried in over 490,000 participants, this largest joint analysis of BMD GWAS identified 12 novel genomic loci at the genome-wide significance level (GWS, p = 5.0 × 10−8), nine of which were for eBMD and four were for TB-BMD, explaining an additional 0.11 and 0.23% heritability for the two traits, respectively. These loci include 1p33 (lead SNP rs10493130, peBMD = 3.19 × 10−8), 5q13.2 (rs4703589, peBMD = 4.78 × 10−8), 5q31.3 (rs9324887, pTB−BMD = 1.36 × 10−9), 6p21.32 (rs6905837, peBMD = 3.32 × 10−8), 6q14.1 (rs10806234, peBMD = 2.63 × 10−8), 7q21.11 (rs10806234, pTB−BMD = 3.37 × 10−8), 8q24.12 (rs11995866, peBMD = 6.72 × 10−9), 12p13.31 (rs1639122, peBMD = 4.43 × 10−8), 12p12.1 (rs58489179, peBMD = 4.74 × 10−8), 12q24.23 (rs75499226, peBMD = 1.44 × 10−8), 19q13.31 (rs7255083, pTB−BMD = 2.18 × 10−8) and 22q11.23 (rs13056137, pTB−BMD = 2.54 × 10−8). All lead SNPs in these 12 loci are nominally significant in both original studies as well as consistent in effect direction between them, providing solid evidence of replication. Approximate conditional analysis identified one secondary signal in 5q13.2 (rs11738874, pconditional = 5.06 × 10−9). Statistical fine-mapping analysis prioritized 269 CRVs. A total of 65 candidate genes were prioritized, including those with known biological function to bone development (such as FGF1, COL11A2 and DEPTOR). Our findings provide novel insights into a better understanding of the genetic mechanism underlying bone development as well as candidate genes for future functional investigation.

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

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

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

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

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

Twelve New Genomic Loci Associated with Bone Mineral Density

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

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

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

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

Genomic structural variations link multiple genes to bone mineral density in a multi-ethnic cohort study: Louisiana osteoporosis study

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

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

Identifying Pleiotropic SNPs Associated With Femoral Neck and Heel Bone Mineral Density

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

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

Genome-wide Association Study of Lncrna Polymorphisms with Bone Mineral Density

Genome-wide identification of m 6 A-associated SNPs as potential functional variants for bone mineral density

Joint Association Analysis Identified 18 New Loci For Bone Mineral Density

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

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