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.

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

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

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

SNP Rs17079281 Decreases Lung Cancer Risk Through Creating an YY1-Binding Site to Suppress DCBLD1 Expression

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

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

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

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

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

An Intronic Risk SNP rs12454712 for Central Obesity Acts As an Allele-Specific Enhancer To Regulate BCL2 Expression

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

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

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

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

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

GATA4 and estrogen receptor alpha bind at SNPs rs9921222 and rs10794639 to regulate AXIN1 expression in osteoblasts

Breast cancer-associated SNP rs72755295 is a cis-regulatory variation for human EXO1

Two novel pleiotropic loci associated with osteoporosis and abdominal obesity

The Genetic Polymorphisms of Key Genes in WNT Pathway (LRP5 and AXIN1) Was Associated with Osteoporosis Susceptibility in Chinese Han Population

Identification of rs11615992 as a novel regulatory SNP for human P2RX7 by allele-specific expression

Genome-wide association and functional interrogation identified a variant at 3p26.1 modulating ovarian cancer survival among Chinese women