In this large population-based cohort study, variants in ESR2 were associated with increased risk of vertebral and incident fragility fracture in postmenopausal women. Interaction of ESR2 with ESR1 and IGF1 was determined and revealed a deleterious genetic combination that enhances the risk of osteoporotic fracture.

Introduction: Osteoporosis is a complex disease with strong genetic influence, but the genes involved are ill-defined. We examined estrogen receptor beta (ESR2) polymorphisms in interaction with estrogen receptor alpha (ESR1) and insulin-like growth factor I (IGF1) variants in relation to the risk of osteoporotic fracture, BMD, and bone geometry.

Materials and methods: In the Rotterdam study, a prospective population-based cohort of elderly white individuals, we studied six single nucleotide polymorphisms (SNPs) in ESR2 (n = 6343, 60% women). We analyzed the genetic variants in the form of haplotypes reconstructed by a statistical method. Results refer to the most frequent ESR2 haplotype 1 estimated from two SNPs in intron 2 and the 3'-untranslated region (UTR). Outcomes included vertebral and incident nonvertebral fractures, BMD, and hip structural analysis (HSA). We also studied the interaction with (the most frequent) ESR1 haplotype 1 estimated from the PvuII and XbaI polymorphisms and an IGF1 promoter CA-repeat.

Results: Compared with ESR2 haplotype 1 noncarriers, female homozygous carriers had a 1.8- and 1.4-fold increased risk of vertebral and fragility fractures. HSA showed that ESR2 haplotype 1 homozygote women had 2.6% thinner cortices, 1.0% increased neck width, and 4.3% higher bone instability (buckling ratios). For testing the gene interaction, we assumed a recessive model of ESR2 haplotype 1. Female homozygous carriers of ESR2 haplotype 1 and noncarriers of ESR1 haplotype 1 had a 3.5- and 1.8-fold increased risk of vertebral and fragility fractures (p(interaction) = 0.10). Such effects and interactions were stronger in women homozygous for the IGF1 192-bp allele, with 9.3-fold increased risk (p(interaction) = 0.002) for vertebral and 4.0-fold increased risk (p(interaction) = 0.01) for fragility fractures. Multilocus interaction analyses of fracture endured correction for multiple testing using Monte-Carlo simulations (p(interaction) = 0.02 for vertebral and p(interaction) = 0.03 for fragility fractures). Similar patterns of interaction were observed for BMD, cortical thickness, bone strength (section modulus), and instability (buckling ratio). In men, no such effects were observed.

Conclusions: Variants of ESR2 alone and in interaction with ESR1 and IGF1 influence the risk of fracture in postmenopausal women. These findings reinforce the polygenic and complex character of osteoporosis.