The osteoporosis susceptible SNP rs4325274 remotely regulates the SOX6 gene through enhancers.

Osteoporosis is a typical polygenic disease, and its heritability is as high as 85%. The incidence of osteoporosis has jumped to the fifth among the common diseases. Although a large number of osteoporosis-susceptible SNPs have been identified, most of them are in the non-coding regions of the genome and the functional mechanisms are unknown. The purpose of this study was to explore the function of non-coding osteoporosis-susceptible SNP rs4325274 and dissect the molecular regulatory mechanisms through integrating bioinformatics analysis and functional experiments. Firstly, we found the SNP rs4325274 resided in a putative enhancer element through functional annotation. eQTL and Hi-C analysis found that the SOX6 gene might be a potential distal target of rs4325274. We conducted the motif prediction using multiple databases and verified the result using ChIP-seq data from GEO database. The result showed that the transcription factor HNF1A could preferentially bind to SNP rs4325274-G allele. We further demonstrated that SNP rs4325274 acted as an enhancer regulating SOX6 gene expression by using dual-luciferase reporter assays. Knockdown of HNF1A decreased the SOX6 gene expression. Taken together, our results uncovered a new mechanism of a non-coding functional SNP rs4325274 as a distal enhancer to modulate SOX6 expression, which provides new insights into deciphering molecular regulatory mechanisms underlying non-coding susceptibility SNPs on complex diseases.