Molecular disease map of bone characterizing the postmenopausal osteoporosis phenotype

Genome‐wide gene expressions in bone biopsies from patients with postmenopausal osteoporosis and healthy controls were profiled, to identify osteoporosis candidate genes. All osteoporotic patients (n = 27) in an unbiased cohort of Norwegian women presented with bone mineral density (BMD) T‐scores of less than −2.5 SD and one or more confirmed low‐energy fracture(s). A validation group (n = 18) had clinical and laboratory parameters intermediate to the control (n = 39) and osteoporosis groups. RNA from iliac crest bone biopsies were analyzed by Affymetrix microarrays and real‐time reverse‐transcriptase polymerase chain reaction (RT‐PCR). Differentially expressed genes in osteoporosis versus control groups were identified using the Bayesian ANOVA for microarrays (BAMarray) method, whereas the R‐package Limma (Linear Models for Microarray Data) was used to determine whether these transcripts were explained by disease, age, body mass index (BMI), or combinations thereof. Laboratory tests showed normal ranges for the cohort. A total of 609 transcripts were differentially expressed in osteoporotic patients relative to controls; 256 transcripts were confirmed for disease when controlling for age or BMI. Most of the osteoporosis susceptibility genes (80%) also were confirmed to be regulated in the same direction in the validation group. Furthermore, 217 of 256 transcripts were correlated with BMD (adjusted for age and BMI) at various skeletal sites (|r| > 0.2, p