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

Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome-wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined superenhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range ‘mSE-miRNA’ approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).