Abstract:RANKL is a key regulator involved in bone metabolism, and a drug target for osteoporosis. The clinical diagnosis and assessment of osteoporosis are mainly based on bone mineral density (BMD). Previous powerful genomewide association studies (GWASs) have identified multiple intergenic single-nucleotide polymorphisms (SNPs) located over 100kb upstream of RANKL and 65kb downstream of AKAP11 at 13q14.11 for osteoporosis. Whether these SNPs exert their roles on osteoporosis through RANKL is unknown. In this study, we conducted integrative analyses combining expression quantitative trait locus (eQTL), genomic chromatin interaction (high-throughput chromosome conformation capture [Hi-C]), epigenetic annotation, and a series of functional assays. The eQTL analysis identified six potential functional SNPs (rs9533090, rs9594738, r8001611, rs9533094, rs9533095, and rs9594759) exclusively correlated with RANKL gene expression (p<0.001) at 13q14.11. Co-localization analyses suggested that eQTL signal for RANKL and BMD-GWAS signal shared the same causal variants. Hi-C analysis and functional annotation further validated that the first five osteoporosis SNPs are located in a super-enhancer region to regulate the expression of RANKL via long-range chromosomal interaction. Particularly, dual-luciferase assay showed that the region harboring rs9533090 in the super-enhancer has the strongest enhancer activity, and rs9533090 is an allele-specific regulatory SNP. Furthermore, deletion of the region harboring rs9533090 using CRISPR/Cas9 genome editing significantly reduced RANKL expression in both mRNA level and protein level. Finally, we found that the rs9533090-C robustly recruits transcription factor NFIC, which efficiently elevates the enhancer activity and increases the RANKL expression. In summary, we provided a feasible method to identify regulatory noncoding SNPs to distally regulate their target gene underlying the pathogenesis of osteoporosis by using bioinformatics data analyses and experimental validation. Our findings would be a potential and promising therapeutic target for precision medicine in osteoporosis. (c) 2018 American Society for Bone and Mineral Research.

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

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

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

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

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

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

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

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

Two novel pleiotropic loci associated with osteoporosis and abdominal obesity

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

Evidence for Major Pleiotropic Effects on Bone Size Variation from a Principal Component Analysis of 451 Caucasian Families

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies.

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

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

Identification of Novel Pleiotropic SNPs Associated with Osteoporosis and Rheumatoid Arthritis

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

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

Is a gene important for bone resorption a candidate for obesity? An association and linkage study on the RANK (receptor activator of nuclear factor-κB) gene in a large Caucasian sample

A New Discovery of Genetic Polymorphisms of Important Genes in WNT Pathway (LPR5 and AXIN1) Associated with Osteoporosis Susceptibility in Chinese Han Population