The variation in several of the risk factors for osteoporotic fracture, including bone mineral density (BMD), has been shown to be strongly influenced by genetic differences. However, the genetic architecture of BMD is complex in both humans and in model organisms. We previously reported quantitative trait locus (QTL) results for BMD from a genome screen of 828 F2 progeny of Copenhagen and dark agouti rats. These progeny also provide an excellent opportunity to search for epistatic effects, or interaction between genetic loci, that contribute to fracture risk. Microsatellite marker data from a 20-cM genome screen was analyzed along with weight-adjusted bone density (DXA and pQCT) phenotypic data using the R/qtl software package. Genotype and phenotype data were permuted to determine genome-wide significance thresholds for the full model and epistasis (interaction) LOD scores corresponding to an alpha level of 0.01. A novel locus on chromosome 15 and a previously reported chromosome 14 QTL demonstrated a strong epistatic effect on BMD at the femur by DXA (LOD = 5.4). Two novel QTLs on chromosomes 2 and 12 were found to interact to affect total BMD at the femur midshaft by pQCT (LOD = 5.0). These results provide new information regarding the mode of action of previously identified QTL in the rat, as well as identifying novel loci that act in combination with known QTL or with other novel loci to contribute to BMD variation.