eISSN: 1896-9151
ISSN: 1734-1922
Archives of Medical Science
Current issue Archive Manuscripts accepted About the journal Special issues Editorial board Abstracting and indexing Subscription Contact Instructions for authors Ethical standards and procedures
SCImago Journal & Country Rank
vol. 15
Letter to the Editor

Association of polymorphisms of the TNFRSF11B and TNFSF11 genes with bone mineral density in postmenopausal women from western Mexico

Anahi González-Mercado
Josefina Y. Sánchez-López
Francisco J. Perea-Díaz
Maria T. Magaña-Torres
Mario Salazar-Páramo
Laura González-López
Mirna Gisel González-Mercado
Bertha Ibarra-Cortés

Arch Med Sci 2019; 15 (5): 1352–1356
Online publish date: 2019/08/22
Article file
Get citation
JabRef, Mendeley
Papers, Reference Manager, RefWorks, Zotero

Osteoporosis (OP) is a disease with reduced bone mass and deterioration of the spatial distribution of the trabecular bone structure, leading in consequence to increased bone fragility and risk of fractures [1]. In Mexico its frequency in postmenopausal women is 8.3% to 31% [2]. OP is a multifactorial disease because its development involves physiological, environmental and genetic factors [3]. Genes of the tumor necrosis factor (TNF) family like TNFRSF11B (TNF receptor superfamily member 11b) and TNFSF11 (TNF superfamily member 11) encode for osteoprotegerin (OPG) and receptor activator for nuclear factor κB ligand (RANKL) proteins, respectively, that participate in the signaling pathway OPG/RANK/RANKL, balancing the activity between osteoblasts and osteoclasts to prevent bone loss and to ensure normal bone turnover [4, 5]. Different single nucleotide polymorphisms (SNPs) in genes of this pathway are related to bone mineral density (BMD) or OP [6, 7]. In Mexican postmenopausal women, SNPs of the TNFRSF11B gene were studied but no association with BMD [8, 9] or OP [10] was observed; SNPs of the TNFSF11 gene have not been studied for OP.

The aim of this study was to investigate the association of the polymorphisms 1181 G>C (rs2073618) and 1217 C>T (rs3102734) of TNFRSF11B; as well as -708 T>A, -693 C>G (rs9533155), -657 C>A, -643 C>T (rs9533156) and -290 C>T (rs9525641) of TNFSF11 with BMD.

Five hundred and thirteen postmenopausal Mexican-Mestizo women from western Mexico agreed to participate in the study. They were recruited from Hospital General Regional No. 110 of the Instituto Mexicano del Seguro Social in Guadalajara city. BMD of the lumbar spinal L1-L4 region and of the femoral neck (left and right) was determined by means of DEXA (dual-energy X-ray absorptiometry, Prodigy Advance, GE). A hundred and seventy-four women were matched by age (range was between 42 to 69 years) and body mass index (BMI) (range from 19.7 to 42.1 kg/m2) in two groups by the WHO criteria: women with values of the T-score less than –2.5 SD (OP group, n = 87) and those with values of the T-score above -1 SD (non-OP group, n = 87). A questionnaire that included different biological (age, BMI, age at onset of menarche and menopause, family history of OP, personal and family history of fractures) and lifestyle (smoking, alcoholism, coffee consumption and sedentarism) variables was applied to all participants. This study was approved by the local ethical committee; informed consent was obtained from all women.

Identification of polymorphisms was performed by PCR followed by Sanger DNA sequencing with the BigDye terminator kit v3.1 (Applied Biosystems Foster City). To identify the rs2073618 and rs3102734 variants, a fragment of 274 bp was amplified using the primers forward 5´-GATCAAAGGCAGGCGATACTTCC- 3´ and reverse 5´-CTGGGAGGGAGCGAGTGGAGC- 3´. The PCR mixture (25 μl) contained: genomic DNA (200 ng), 10 pmol each primer, Taq polymerase 1 U, MgCl2 2.0 mM, dNTPs 250 μmol and buffer 1X. The thermocycling conditions were: initial denaturation 95°C/5 min, denaturation 95°C/1 min, annealing 57°C/1 min, extension 72°C/1:30 min, 35 cycles and a final extension 72°C/5 min. The -708 T>A, rs9533155, -657 C>A, rs9533156 and rs9525641 polymorphisms were identified according to Mencej et al. [11].

For categorical variables, the frequencies or percentages were obtained; for continuous variables, the average and standard deviation were calculated. Chi square and Fisher exact tests were used for comparison of genotypic and allelic frequencies of each polymorphism between both groups. The risk for OP was estimated through odds ratio, considering classic models of inheritance. Haplotypes were established inferring phase of family segregation. Linkage disequilibrium (LD) was significant when r 2 ≥ 0.33, using Arlequin v3.01 software. ANOVA and Student t-tests were used to relate quantitative variables with genotypes and alleles, respectively. Other tests used were Spearman and Pearson correlations and linear regression models. Statistical analyses were carried out using SPSS version 20.0 and p < 0.05 was significant.

As expected, BMD values in the group of postmenopausal women with OP (0.834 ±0.052 g/cm2 in lumbar spine L1–L4; 0.788 ±0.086 g/cm2 in left femoral neck; 0.787 ±0.085 g/cm2 in right femoral neck) were lower than in the non-OP group (1.173 ±0.109 g/cm2 in lumbar spine L1-L4; 1.005 ±0.087 g/cm2 in left femoral neck; 1.003 ±0.083 g/cm2 in right femoral neck) (p < 0.001). The other biological and lifestyle variables did not show significant differences between the groups, even when they have been reported as high risk for the development of OP in some populations [12, 13].

The distribution of genotypic and allelic frequencies of the seven SNPs was similar between the groups (p > 0.05). The SNPs –708 T>A and –657 C>A were observed with low frequencies in both groups, and the –708A and –657A alleles were not detected in homozygote state in the study (Table I). The most widely studied polymorphism of TNFRSF11B is 1181 G>C; the distribution of allelic frequencies is variable in women with OP from several populations [1416]. Under classical models of inheritance (dominant, recessive, codominant or additive) neither genotypes nor alleles were associated with risk for OP (p > 0.05). Regarding TNFSF11, in the Egyptian population, the –290T allele was associated with increased risk of developing OP (OR = 1.7, IC: 1.1–2.7, p = 0.019) [17].

Table I

Genotypic and allelic frequencies of the seven polymorphisms studied of the TNFRSF11B and TNFSF11 genes

VariableTNFRSF11B geneTNFSF11 gene
1181 G>C1217 C>T–708 T>A–693 C>G–657 C>A–643 C>T–290 C>T
OP %Non-OP %OP %Non-OP %OP %Non-OP %OP %Non-OP %OP %Non-OP %OP %Non-OP %OP %Non-OP %
Genotypic frequencies:
Allelic frequencies:

[i] OP – osteoporosis, WW – wild homozygote genotype, WM – heterozygote genotype, MM – mutated homozygote genotype, W – wild allele, M – mutated allele. N = 87 OP group, n = 87 non-OP group.

Haplotypes for TNFRSF11B were constructed with two polymorphisms (1181 G>C and 1217 C>T) and three haplotypes were observed; the most frequent was CC (OP 54.6% and non-OP 54.0%). In TNFSF11, the haplotypes were constructed with five SNPs: –708 T>A, –693 C>G, –657 C>A, –643 C>T and –290 C>T, and 13 haplotypes were observed, the most frequent in both groups being TCCTT (OP = 49.4% and non-OP = 54.1%). The distribution of haplotype frequencies was compared between the two groups and no statistically significant differences were found (p > 0.05), with the exception of ACATT haplotype (OP group = 9.8%, non-OP group = 4.0%; p = 0.03). The same combinations of SNPs for the two genes are not reported in the literature.

The SNPs of TNFRSF11B were in linkage equilibrium in the two groups (OP: r 2 = 0.089; non-OP: r 2 = 0.127). For the five sites of TNFSF11, we found that four pairs of loci are in LD in both groups, since r 2 values ranged from 0.6587 to 0.9058 with p < 0.001: 1) –708 T>A/-657 C>A, 2) –693 C>G/–643 C>T, 3) –693 C>G/-290 C>T and 4) –643 C>T/–290 C>T. In women with OP and non-OP from Slovenia, LD was also observed between pairs 2, 3 and 4 [18, 19].

The relationship between genotypes and alleles of the studied polymorphisms with the biological variables in each group was analyzed; only the BMD and BMI presented statistically significant results with some SNPs (Table II). For BMD, no association was found with 1181 G>C and 1217 C>T SNPs of TNFRSF11B; similar results were observed in postmenopausal women from Australia [20], Ireland [14, 21], and Malta [22]; however, the 1181GG genotype was associated with lower BMD in Chinese [7, 16], Slovenian [15, 23], and Koreans [24]; in China, the 1217T allele was related to increased risk of osteoporotic fracture (OR = 1.35; 95% CI: 1.17–1.55; Bonferroni p = 2.6 × 10–4) [25]. Regarding BMI, in the present study, the 1181GG genotype was significantly associated with lower BMI in women with OP (Table II). Although this relationship has not been reported, we can suggest that this polymorphism could have an indirect effect on the variation of the BMD, as was demonstrated by Mendez et al. [26], in which a high BMI is positively correlated with higher BMD. It has been proposed that obesity protects against bone loss through mechanical load, and it has been suggested that the adipose tissue is involved in the homeostasis of the skeleton, possibly through the role of some adipokines in bone remodeling. Also, adipocytes secrete estrogens that help to regulate homeostasis and contribute to the increase of bone mass [26].

Table II

Association between genotypes and alleles of the studied polymorphisms with biological variables

TNFRSF11B1181 G>CBMI [kg/cm2]GG 26.6CC 28.70.05OP
GC 29.30.02
TNFSF11–693 C>GBMD LFN [g/cm2]GG 0.840CC 0.7810.05OP
CG 0.7810.03
BMD RFN [g/cm2]GG 0.846CC 0.7740.02OP
G 0.806C 0.7760.03OP
–643 C>TBMI [kg/cm2]TT 28.0CT 30.20.03Non-OP
BMD LFN [g/cm2]CC 0.833TT 0.7790.05OP
CT 0.7820.04
BMD RFN [g/cm2]CC 0.830TT 0.7740.04OP
C 0.804T 0.7770.04OP
–290 C>TBMI [kg/cm2]TT 27.7CT 30.40.01Non-OP
BMD RFN [g/cm2]CC 0.828CT 0.7780.03OP

[i] BMI – body mass index, OP – osteoporosis, BMD – bone mineral density, LFN – left femoral neck, RFN – right femoral neck. n = 87 OP group, n = 87 non-OP group.

On the other hand, the –693CC, –643TT and –290CT genotypes of TNFSF11 reveal lower BMD in the left and/or right femoral neck than the other genotypes in postmenopausal women with OP, which demonstrates the important role of these SNPs in the variation of BMD (Table II). The linear regression analysis with all studied polymorphisms and BMD of the lumbar spine, and the right and left femoral necks, showed only a significant model between the –693G allele and BMD of the right femoral neck (constant = 0.738, β = 0.028, p = 0.032), indicating that there is a positive relationship since the BMD is increased slightly. These three polymorphisms in Slovenian women with OP were associated with BMD of the lumbar spine and –693 C>G also with BMD of the femoral neck [18]. In the Chinese population –693 C>G or –643 C>T polymorphism was not related to BMD in peri- and postmenopausal women [7]. The –643TT and –290TT genotypes of TNFSF11 in the non-OP group were related to lower BMI, which could suggest a possible influence on the BMD variation as explained above for 1181 G>C [26].

In conclusion, the 1181 G>C and 1217 C>T polymorphisms of TNFRSF11B and –708 T>A, –693 C>G, –657 C>A, –643 C>T and –290 C>T of TNFSF11 were not associated with OP. However, –693CC, –643TT and –290CT genotypes have an effect of lowering BMD in the left and/or right femoral neck. The 1181GG genotype was significantly associated with lower BMI in women with OP as well as –643TT and –290TT genotypes of TNFSF11 in the non-OP group. The –708 T>A and –657 C>A polymorphisms have not been described in the literature.


The work was conducted in División de Genética, CIBO, IMSS, Guadalajara, Jalisco.

The research was supported by a research grant from Consejo Nacional de Ciencia y Tecnología #2003-107.

Conflict of interest

The authors declare no conflict of interest.



Bączyk G, Samborski W, Jaracz K , authors. Evaluation of the quality of life of postmenopausal osteoporotic and osteopenic women with or without fractures. Arch Med Sci. 2016. 12:p. 819–27


González-Mercado A, Sánchez-López JY, Ibarra B , authors. Factores de riesgo para osteoporosis en mujeres posmenopáusicas de Guadalajara, Jalisco. Salud Publica Mex. 2013. 55:p. 627–30


McCormick RK , author. Osteoporosis: integrating biomarkers and other diagnostic correlates into the management of bone fragility. Altern Med Rev. 2007. 12:p. 113–45


Kearns AE, Khosla S, Kostenuik PJ , authors. Receptor activator of nuclear factor kappa b ligand and osteoprotegerin regulation of bone remodeling in health and disease. Endocr Rev. 2008. 29:p. 155–92


Trouvin AP, Goëb V , authors. Receptor activator of nuclear factor-kappaB ligand and osteoprotegerin: maintaining the balance to prevent bone loss. Clin Interv Aging. 2010. 5:p. 345–54


Tu P, Duan P, Zhang RS, et al. , authors. Polymorphisms in genes in the RANKL/RANK/OPG pathway are associated with bone mineral density at different skeletal sites in post-menopausal women. Osteoporos Int. 2015. 26:p. 179–85


Shang M, Lin L, Cui H , authors. Association of genetic polymorphisms of RANK, RANKL and OPG with bone mineral density in Chinese peri- and postmenopausal women. Clin Biochem. 2013. 46:p. 1493–501


Canto-Cetina T, Polanco Reyes L, González Herrera L, et al. , authors. Polymorphism of LRP5, but not of TNFRSF11B, is associated with a decrease in bone mineral density in postmenopausal Maya-Mestizo women. Am J Hum Biol. 2013. 25:p. 713–8


Rojano-Mejía D, Coral-Vázquez RM, Espinosa LC, et al. , authors. TNFRSF11B gene haplotype and its association with bone mineral density variations in postmenopausal Mexican-Mestizo women. Maturitas. 2012. 71:p. 49–54


Zavala-Cerna MG, Moran-Moguel MC, Cornejo-Toledo JA, et al. , authors. Osteoprotegerin polymorphisms in a mexican population with rheumatoid arthritis and generalized osteoporosis: a preliminary report. J Immunol Res. 2015. 2015:p. 376197


Mencej S, Prezelj J, Kocijancic A, Ostanek B, Marc J , authors. Association of TNFSF11 gene promoter polymorphisms with bone mineral density in postmenopausal women. Maturitas. 2006. 55:p. 219–26


Rossini M, Adami S, Bertoldo F, et al. , authors. Guidelines for the diagnosis, prevention and management of osteoporosis. Reumatismo. 2016. 68:p. 1–39


Gronholz MJ , author. Prevention, diagnosis, and management of osteoporosis-related fracture: a multifactoral osteopathic approach. J Am Osteopath Assoc. 2008. 108:p. 575–85


Langdahl BL, Carstens M, Stenkjaer L, Eriksen EF , authors. Polymorphisms in the osteoprotegerin gene are associated with osteoporotic fractures. J Bone Miner Res. 2002. 17:p. 1245–55


Mencej-Bedrač S, Preželj J, Marc J , authors. TNFRSF11B gene polymorphisms 1181G > C and 245T > G as well as haplotype CT influence bone mineral density in postmenopausal women. Maturitas. 2011. 69:p. 263–7


Zhao HY, Liu JM, Ning G, et al. , authors. The influence of Lys3Asn polymorphism in the osteoprotegerin gene on bone mineral density in Chinese postmenopausal women. Osteoporos Int. 2005. 16:p. 1519–24


Mohamed RH, Mohamed RH, El-Shahawy EE , authors. Relationship between RANK and RANKL gene polymorphisms with osteoporosis in rheumatoid arthritis patients. Genet Test Mol Biomarkers. 2016. 20:p. 249–54


Mencej S, Albagha OME, Preželj J, Kocjan T, Marc J , authors. Tumour necrosis factor superfamily member 11 gene promoter polymorphisms modulate promoter activity and influence bone mineral density in postmenopausal women with osteoporosis. J Mol Endocrinol. 2008. 40:p. 273–9


Mencej-Bedrač S, Preželj J, Kocjan T, et al. , authors. The combinations of polymorphisms in vitamin D receptor, osteoprotegerin and tumour necrosis factor superfamily member 11 genes are associated with bone mineral density. J Mol Endocrinol. 2009. 42:p. 239–47


Ueland T, Bollerslev J, Wilson SG, et al. , authors. No associations between OPG gene polymorphisms or serum levels and measures of osteoporosis in elderly Australian women. Bone. 2007. 40:p. 175–81


Wynne F, Drummond F, O’Sullivan K, et al. , authors. Investigation of the genetic influence of the OPG, VDR (Fok1), and COLIA1 Sp1 polymorphisms on BMD in the Irish population. Calcif Tissue Int. 2002. 71:p. 26–35


Vidal C, Brincat M, Xuereb Anastasi A , authors. TNFRSF11B gene variants and bone mineral density in postmenopausal women in Malta. Maturitas. 2006. 53:p. 386–95


Arko B, Preželj J, Kocijancic A, Komel R, Marc J , authors. Association of the osteoprotegerin gene polymorphisms with bone mineral density in postmenopausal women. Maturitas. 2005. 51:p. 270–9


Choi JY, Shin A, Park SK, et al. , authors. Genetic polymorphisms of OPG, RANK, and ESR1 and bone mineral density in Korean postmenopausal women. Calcif Tissue Int. 2005. 77:p. 152–9


Wang C, Zhang Z, Zhang H, et al. , authors. Susceptibility genes for osteoporotic fracture in postmenopausal Chinese women. J Bone Miner Res. 2012. 27:p. 2582–91


Méndez JP, Rojano-Mejía D, Pedraza J, et al. , authors. Bone mineral density in postmenopausal Mexican-Mestizo women with normal body mass index, overweight, or obesity. Menopause. 2013. 20:p. 568–72

Copyright: © 2019 Termedia & Banach. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
Quick links
© 2020 Termedia Sp. z o.o. All rights reserved.
Developed by Bentus.
PayU - płatności internetowe