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.

Genome-wide identification of m 6 A-associated SNPs as potential functional variants 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 Genomic Analysis Predicts Novel Functional Enhancer-Snps for Bone Mineral Density

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

Genome-wide Enrichment of M6a-Associated Single-Nucleotide Polymorphisms in the Lipid Loci

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

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

Examination of the Associations Between M6a-Associated Single-Nucleotide Polymorphisms and Blood Pressure

Association tests of interleukin-6 ( IL-6 ) and type II tumor necrosis factor receptor ( TNFR2 ) genes with bone mineral density in Caucasians using a re-sampling approach

Twelve New Genomic Loci Associated with Bone Mineral Density

Detection of M(6)a-Associated SNPs As Potential Functional Variants for Coronary Artery Disease

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

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

Genetic variants in the SOX6 gene are associated with bone mineral density in both Caucasian and Chinese populations

SNP rs6265 Regulates Protein Phosphorylation and Osteoblast Differentiation and Influences BMD in Humans

Identification of IDUA and WNT16 Phosphorylation-Related Non-Synonymous Polymorphisms for Bone Mineral Density in Meta-Analyses of Genome-Wide Association Studies.

A Functional SNP Regulated by Mir-196A-3p in the 3'UTR of FGF2 is Associated with Bone Mineral Density in the Chinese Population.

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

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

Joint Association Analysis Identified 18 New Loci For Bone Mineral Density

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