# Gene-gene and gene-environment interactions of CYP19A1, ESR1, IL6, IL6R, IL1β, RANK, and RANKL variants in relation to osteoporosis and hip fracture risk in Mexican women

## TL;DR

This study explores how genetic variants and their interactions with each other and the environment influence osteoporosis and hip fracture risk in Mexican women.

## Contribution

The study identifies novel gene-gene and gene-environment interactions in Mexican women related to osteoporosis and hip fracture risk.

## Key findings

- RANK rs3018362 GA and AA genotypes increase osteoporosis risk, while ESR1 rs2234693 CC genotype reduces it.
- A gene-gene interaction between ESR1 and RANK increases osteoporosis risk, and another involving ESR1, IL6R, IL1β, and RANKL is linked to hip fracture.
- Gene-environment interactions involving CYP19A1, IL1β, and fracture probability are significant in osteoporosis susceptibility.

## Abstract

Osteoporosis is a complex disease influenced by genetic variants, environmental factors, and comorbidities. While individual single-nucleotide variants (SNVs) have been associated with disease risk, limited data are available on how gene–gene and gene–environment interactions contribute to osteoporosis and fracture susceptibility in Mexican women.

In this case–control study, we evaluated the association of SNVs in estrogen receptor alpha (ESR1), aromatase (CYP19A1), interleukin 6 (IL6) and its receptor (IL6R), interleukin 1 beta (IL1β), Receptor activator of nuclear factor κ B (RANK) and its ligand (RANKL) genes, with the risk of osteoporosis and hip fracture, as well as their gene-gene and gene-environment interactions, in 609 Mexican women (169 with osteoporosis, 205 with hip fracture, and 235 controls), by real time PCR with TaqMan probes. In addition, multifactor dimensionality reduction (MDR) software was used to detect gene–gene interactions and gene-environment interactions.

The GA and AA genotypes of RANK rs3018362 were significantly associated with increased osteoporosis risk (OR = 2.08 [1.08–3.98] and 2.76 [1.21–6.30], respectively) while the CC genotype of ESR1 rs2234693 was associated with reduced risk (OR = 0.28 [0.11–0.69]). For hip fracture, ESR1 rs2234693 (CC genotype) was protective (OR = 0.30 [0.12–0.75]), whereas RANK rs3018362 (AA genotype) increased risk (OR = 2.4 [1.01–6.06]). A significant gene–gene interaction between ESR1 (rs2228480) and RANK (rs3018362) increased osteoporosis risk (OR = 2.1 [1.4–3.2], CVC = 10/10), and a gene–gene model involving ESR1, IL6R, IL1β, and RANKL was identified for hip fracture (CVC = 8/10). In osteoporosis, a gene–environment interaction was observed between CYP19A1 SNVs (rs700518, rs1062033, rs4775936, rs767199), IL1β rs16944, and 10-year probability of major fracture (CVC = 10/10).

Our findings suggest that RANK and ESR1 variants are independently and interactively associated with osteoporosis and fracture risk, and that gene–gene and gene–environment interactions play a critical role in disease susceptibility among Mexican women.