Bone-related complications in Gaucher disease are considered to be poorly responsive to specific enzyme replacement therapy. Polymorphisms of candidate genes associated with low bone density were investigated to see whether they are correlated with bone mineral density (BMD) and bone involvement in Gaucher disease. Genotyping for polymorphisms in candidate genes (interleukins 1alpha and 1beta, interleukin-1 receptor antagonist; cytochrome P450; collagen 1A1; low-density Lipoprotein Receptor; bone morphogenic protein 4; vitamin D receptor; and estrogen receptor 2beta) were performed using standard methodologies. BMD was measured by dual energy X-ray absorptiometry (DXA). One hundred and ninety-four patients and 100 controls were genotyped for the above polymorphisms. Thirteen haplotypes were obtained, with several correlations with BMD in patients; also, a haplotype (T889-T3954-C511-240VNTR of IL1) was significantly correlated with T-scores and Z-score for femur neck and lumbar spine (p = 0.01) in patients. Haplotypes of bone-specific candidate genes associated with BMD may predict severity of these features in Gaucher disease.