The skeleton has the ability to alter its mass, geometry, and strength in response to mechanical stress. In order to elucidate the molecular mechanisms underlying this phenomenon, differential display reverse transcriptase-polymerase chain reaction (DDRT-PCR) was used to analyze gene expression in endocortical bone of mature female rats. Female Sprague-Dawley rats, approximately 8 months old, received either a sham or bending load using a four-point loading apparatus on the right tibia. RNA was collected at 1 h and 24 h after load was applied, reverse-transcribed into cDNA, and used in DDRT-PCR. Parallel display of samples from sham and loaded bones on a sequencing gel showed several regulated bands. Further analysis of seven of these bands allowed us to isolate two genes that are regulated in response to a loading stimulus. Nucleotide analysis showed that one of the differentially expressed bands shares 99% sequence identity with rat osteopontin (OPN), a noncollagenous bone matrix protein. Northern blot analysis confirms that OPN mRNA expression is increased by nearly 4-fold, at 6 h and 24 h after loading. The second band shares 90% homology with mouse myeloperoxidase (MPO), a bactericidal enzyme found primarily in neutrophils and monocytes. Semiquantitative PCR confirms that MPO expression is decreased 4- to 10-fold, at 1 h and 24 h after loading. Tissue distribution analysis confirmed MPO expression in bone but not in other tissues examined. In vitro analysis showed that MPO expression was not detectable in total RNA from UMR 106 osteoblastic cells or in confluent primary cultures of osteoblasts derived from either rat primary spongiosa or diaphyseal marrow. Database analysis suggests that MPO is expressed by osteocytes. These findings reinforce the association of OPN expression to bone turnover and describes for the first time, decreased expression of MPO during load-induced bone formation. These results suggest a role for both OPN and MPO expression in bone cell function.