Bone mineral density (BMD) is a key indicator for diagnosis and treatment for osteoporosis; the reduction of BMD could increase the risk of osteoporotic fracture. It was very recently found that Piezo1 mediated mechanically evoked responses in bone and further participated in bone formation in mice. Here, we performed cross phenotype meta-analysis for human BMD at lumbar spine (LS), femoral neck (FN), distal radius/forearm (FA) and heel and screened out 14 top SNPs for PIEZO1, these SNPs were overlapped with putative enhancers, DNase-I hypersensitive sites and active promoter flanking regions. We found that the signal of the best SNP rs62048221 was mainly from heel ultrasound estimated BMD (-0.02 SD per T allele, P = 8.50E-09), where calcaneus supported most of the mechanical force of body when standing, walking and doing physical exercises. Each copy of the effect allele T of SNP rs62048221 was associated with a decrease of 0.0035 g/cm BMD (P = 4.6E-27, SE = 0.0003) in UK Biobank data within 477,760 samples. SNP rs62048221 was located at the enhancer region (HEDD enhancer ID 2331049) of gene PIEZO1, site-directed ChIP assays in human mesenchymal stem cells (hMSCs) showed significant enrichment of H3K4me1 and H3K27ac in this region, luciferase assays showed that rs62048221 could significantly affect the activity of the enhancer where it resides. Our results first suggested that SNP rs62048221 might mediate the PIEZO1 expression level via modulating the activity of cis-regulatory elements and then further affect the BMD.Copyright © 2018. Published by Elsevier Inc.