Osteoprotegerin Polymorphisms in a Mexican Population with Rheumatoid Arthritis and Generalized Osteoporosis: A Preliminary Report

Abstract

Bone disease in rheumatoid arthritis (RA) is a complex phenomenon where genetic risk factors have been partially evaluated. The system formed by receptor activator for nuclear factor-κB (RANK), receptor activator for nuclear factor-κB ligand (RANKL), and osteoprotegerin (OPG): RANK/RANKL/OPG is a crucial molecular pathway for coupling between osteoblasts and osteoclasts, since OPG is able to inhibit osteoclast differentiation and activation. We aim to evaluate the association between SNPs C950T (rs2073617), C209T (rs3134069), T245G (rs3134070) in the TNFRSF11B (OPG) gene, and osteoporosis in RA. We included 81 women with RA and 52 healthy subjects in a cross-sectional study, genotyped them, and measured bone mineral density (BMD) at the lumbar spine and the femoral neck. Mean age in RA was 50±12 with disease duration of 12±8 years. According to BMD results, 23 (33.3%) were normal and 46 (66.7%) had osteopenia/osteoporosis. We found a higher prevalence of C allele for C950T SNP in RA. Polymorphisms C209T and T245G did not reach statistical significance in allele distribution. Further studies including patients from other regions of Latin America with a multicenter design to increase the sample size are required to confirm our findings and elucidate if C950T SNP could be associated with osteoporosis in RA.

Maria Guadalupe Zavala-Cerna 1 and Maria Cristina Moran-Moguel 2 and Jesus Alejandro Cornejo-Toledo 2 and Norma Guadalupe Gonzalez-Montoya 3 and Jose Sanchez-Corona 2 and Mario Salazar-Paramo 4 and Arnulfo Hernan Nava-Zavala 1, 5, 6 and Erika Anita Aguilar-Chavez 7 and Miriam Fabiola Alcaraz-Lopez 8 and Alicia Guadalupe Gonzalez-Sanchez 3 and Laura Gonzalez-Lopez 9 and Jorge Ivan Gamez-Nava 10

Academic Editor:Giacomina Brunetti

1, Programa Internacional, Facultad de Medicina, Universidad Autonoma de Guadalajara, 44100 Zapopan, JAL, Mexico
2, Division de Medicina Molecular del Centro de Investigacion Biomedica de Occidente, Instituto Mexicano del Seguro Social, 44340 Guadalajara, JAL, Mexico
3, Programa de Becarios en Investigacion del Instituto Mexicano del Seguro Social, Programa de Doctorado en Farmacologia, Centro Universitario de Ciencias de la Salud, Universidad de Guadalajara, 44340 Guadalajara, JAL, Mexico
4, Division de Investigacion en Salud, UMAE, Hospital de Especialidades, CMNO, IMSS and Departamento de Fisiologia, CUCS, Universidad de Guadalajara, 44340 Guadalajara, JAL, Mexico
5, Unidad de Investigacion en Epidemiologia Clinica, Hospital de Especialidades, Centro Medico Nacional de Occidente del Instituto Mexicano del Seguro Social, 44340 Guadalajara, JAL, Mexico
6, Servicio de Medicina Interna, Inmunologia y Reumatologia, Hospital General de Occidente, Secretaria de Salud Jalisco, 45170 Guadalajara, JAL, Mexico
7, Unidad de Medicina Familiar 2, Instituto Mexicano del Seguro Social, 44340 Guadalajara, JAL, Mexico
8, Departamento de Medicina Interna-Reumatologia, Hospital General de Zona 14, Instituto Mexicano del Seguro Social, 44860 Guadalajara, JAL, Mexico
9, Departamento de Medicina Interna-Reumatologia del Hospital General Regional No. 110, Instituto Mexicano del Seguro Social y Programa de Doctorado en Salud Publica, CUCS, Universidad de Guadalajara, 44340 Guadalajara, JAL, Mexico
10, Unidad de Investigacion en Epidemiologia Clinica, Hospital de Especialidades, Centro Medico Nacional de Occidente del Instituto Mexicano del Seguro Social y Programa de Doctorado en Salud Publica, CUCS, Universidad de Guadalajara, 44340 Guadalajara, JAL, Mexico

Received 17 July 2014; Revised 12 September 2014; Accepted 15 September 2014; 3 May 2015

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Rheumatoid arthritis (RA) is a chronic inflammatory disease that affects synovial joints [1]. Its prevalence in Mexico is considerably elevated ranging from 0.3 to 1% [2]. Currently, the diagnosis of osteoporosis is made on the basis of a measurement by dual-energy x-ray absorptiometry (DXA) at the spine and hip showing a decrement in bone mineral density (BMD) lower than 2.5 standard deviations in comparison with healthy young adults [3]. However, bone involvement in RA is complex; an increased proportion of patients may suffer from osteoporosis; nevertheless causal effects are considered to be multifactorial. Previous studies performed in our center showed a prevalence of axial osteoporosis in Mexican women with RA of 31% [4].

Genetic risk factors constitute a variable that contributes to influencing the development of osteoporosis in RA, in this context; some genetic polymorphisms may influence the rate of presentation of osteoporosis. Nevertheless such genetic changes do not explain entirely osteoporosis presence or fracture related events. Strong evidence supports the fact that osteoclasts have an important role in osteoporosis and bone erosion formation [5, 6].

The system conformed by nuclear factor-κ B (RANK)/receptor activator for nuclear factor-κ B ligand (RANKL) and osteoprotegerin (OPG) are crucial, for the bone remodeling process. Furthermore, it appeared that not only osteoblasts but also activated [figure omitted; refer to PDF] lymphocytes play a crucial role in the pathogenesis of rheumatoid arthritis (RA) and osteoporosis, since they can produce RANKL, which stimulates the differentiation and activation of osteoclasts [7]. These three proteins belong to the superfamily of the TNF- [figure omitted; refer to PDF] and have different patterns of expression; RANKL is expressed in osteoblasts and its function is to participate in osteoclasts differentiation [7, 8]. Osteoprotegerin (OPG) lacks a transmembrane domain, and hence it is a receptor that can be secreted. OPG recognizes and is able to bind RANKL, blocking its interaction with RANK inhibiting osteoclastic differentiation and activation [9-11]. Since OPG is homologous to RANK it acts as a false receptor for RANKL [12-14]. OPG is secreted mainly by stromal cells but can also be secreted by other types of cells [15]. Different studies have been able to identify single nucleotide polymorphisms (SNPs) in TNFSF11 (RANKL) and TNFRSF11B (OPG) genes in association with osteoporosis or low bone density [12, 16]. It has been documented that OPG overexpression induces an increase in bone density and protects from the development to osteoporosis. In contrast mice that are OPG-deficient develop osteoporosis very rapidly [14, 15, 17]. OPG overexpression blocks osteoclasts maturation; contrariwise, an underexpression of this molecule results in an excessive increase of bone resorption and osteoporosis [18]. Mutations in the promoter region of TNFRSF11B could have influence on the transcription and translation rate of OPG. A study in European population (Eslovenia) showed that two SNPs named rs3134069 and rs3134070 in the TNFRSF11B could form a haplotype with susceptibility to osteoporosis [16, 19]. To date, there is a lack of information evaluating the relationship between polymorphisms of OPG gene and bone mineral density.

Therefore, the aim of the present study was to assess if there is an association between SNPs C950T (rs2073617), C2097 (rs3134069), and T245G (rs3134070) in the TNFRSF11B (OPG) gene and osteoporosis in rheumatoid arthritis patients.

2. Patients and Methods

2.1. Study Subjects

We carried out a cross-sectional study that included 81 consecutive unrelated women who were diagnosed with RA from March 2007 to March 2009. These patients were referred from an outpatient rheumatology clinic in a secondary care center in Guadalajara, Mexico (Hospital General Regional 110, IMSS). Inclusion criteria were as follows: age ≥ 18 years at entry; self-identified race Mexican Mestizo (defined as individuals who were born in Mexico and were descendants of the original autochthonous inhabitants of the region and of individuals mainly Spaniards); fulfillment of the 1987 American College of Rheumatology (ACR) classification criteria [20]; and no previous BMD measurement. Patients were excluded if pregnancy was present and if they had an overlap syndrome, were receiving bisphosphonates or parathyroid hormone therapy, or had a comorbidity associated with low BMD, such as diabetes mellitus, thyroid disease, or chronic renal failure; patients were not related to each other.

A healthy control group matched by age and sex was obtained from blood donors who attended to "Centro Medico Nacional de Occidente" blood bank; all blood donors were of Mexican Mestizo origin and denied at the time of the study having any disease (acute or chronic).

2.2. Clinical Assessment

Each patient was interviewed using a structured questionnaire to record demographic information, general risk factors for osteoporosis (i.e., age, weight, and height), clinical characteristics, and RA treatment. At the time of the evaluation two trained researchers assessed RA disease activity by systematical evaluation of the DAS-28 index [21]; to determine patient disability the Spanish modified version of the Health Assessment Questionnaire-Disability Index (HAQ-DI) [22] was used. Global functional status was evaluated according to the Steinbrocker classification [23]. Rheumatoid factor and C-reactive protein (CRP) were measured by nephelometry using commercially available kits. Erythrocyte sedimentation rate (ESR, mm/h) was measured using the Wintrobe method.

2.3. BMD Measurement

BMD was measured (g/cm² ) by DEXA using a Lunar Prodigy densitometer (GE Medical Systems Lunar, Madison, WI, software V 8.8; GE Medical Systems); the anatomical regions assessed were lumbar spine in the posterior-anterior projection (L1-L4) and femoral neck. The coefficient of variation during the measurement of a standard phantom in our laboratory is 0.7%. The coefficient of variation was 2.4% at the lumbar spine and 1.6% at the femoral neck. All scans were performed by the same experienced technician, who was blinded to the characteristics of patients. Each patient was classified into one of the following categories proposed by WHO: normal BMD, defined as having a measurement of BMD within 1 SD of the BMD of normal young adult women (equivalent to [figure omitted; refer to PDF] score > -1); osteopenia, defined as having a [figure omitted; refer to PDF] score in BMD between -1 and -2.4 SD; and osteoporosis, defined as having a [figure omitted; refer to PDF] score in BMD <= -2.5 SD. Patients with a [figure omitted; refer to PDF] score < -1.0 were considered to have low BMD in consequence; for analysis purposes in our study we decided to include osteopenia and osteoporosis in one single group denominated as "low BMD".

2.4. Genotyping

DNA from 81 patients with rheumatoid arthritis was extracted from blood samples using conventional methods [24] and stored frozen at -80°C.

We chose three SNPs that are present in the promoter region (5[variant prime]UTR) of the TNFRSF11B (OPG) gene; such polymorphisms were previously reported in association with osteoporosis in postmenopausal women or related conditions [19, 25-28]. We performed 3 polymerase chain reaction-restriction fragment length polymorphisms (PCR_RFLPs) protocols as previously described [26, 29, 30]. Table 1 shows primer sequence, products obtained after PCR, restriction enzymes used, their recognition site, and end products after restriction.

Table 1: Genotyping strategies for TNFRSF11B polymorphism variants detection.

SNP	Primers	Band size (bp)	Restriction enzyme	Recognized sequence	Band size after digestion (bp)
C950Trs2073617	5[variant prime]-TGCGTCCGGATCTTGGCTGGATCGG-3[variant prime]	548	Hinc II	GTY RAC	C 248 and 287T 548
5[variant prime]-ACTTACCACGAGCGCGCAGCACAGCAA-3[variant prime]
C209Trs3134070	5[variant prime]-CGAACCCTAGAGCAAAGTGC-3[variant prime]	271	Taq I	T CGA	C 212, 31 and 28T 212 and 59
5[variant prime]-TGTCTGATTGGCCCTAAAGC-3[variant prime]
T245Grs3134069	5[variant prime]-CGAACCCTAGAGCAAAGTGC-3[variant prime]	271	Hinf I	G ANTC	A 195, and 176C 127, 76 and 69
5[variant prime]-TGTCTGATTGGCCCTAAAGC-3[variant prime]

SNP: single nucleotide polymorphism and bp: base pair.

Each PCR reaction was carried out in 10 μ L final volume containing (final concentrations) 1X buffer (200 mM Tris-HCl pH 8.4, 500 mM KCl, and 4 mM MgCl₂ ); 5 pmol/mL each of the pair primers according to polymorphism (Table 1); 10 mM each of the four deoxyribonucleoside triphosphates; 1 U of taq DNA polymerase (Invitrogen, Carlsbad, CA, USA); and 200 to 300 ng DNA template. The PCR products were visualized by electrophoresis in 10% polyacrylamide gels at 150 V for 90 min, followed by silver staining.

2.5. Statistical Analysis

Allelic and genotypic frequencies for each SNP were determined by gene count. Hardy-Weinberg equilibrium was tested, using [figure omitted; refer to PDF] test. Genotypic differences between RA patients with and without osteoporosis were evaluated by Mantel-Haenszel test using the EPI INFO (version 6.04d) statistical program. For allelic differences, we used Fisher's exact test. Odds ratios (OR) and 95% confidence interval (95% CI) were computed to evaluate the risk for osteoporosis conferred by presence of the risk alleles. All analyses used two-sided tails with a [figure omitted; refer to PDF] value of <=0.05 used as significance criterion.

3. Results

3.1. Demographics and Clinical

We included 81 patients, classified as RA according to ACR classification criteria [20], and 52 healthy subjects. For the RA group, mean age was [figure omitted; refer to PDF] years, 100% of the studied subjects were women, mean duration disease was [figure omitted; refer to PDF] years, and all of them were receiving treatment, some of them with monotherapy and others combined therapy; the number and percentages of patients taking a specific drug are listed in Table 2.

Table 2: Clinical characteristics of study patients.

Characteristic	[figure omitted; refer to PDF] = 81
Age (years), median (range)	49.9 (21-79)
Menopausal, [figure omitted; refer to PDF] (%)	43 (53.1)
Weight (kg), median (range)	65.8 (43-94)
BMI, median (range)	26.87 (18-38)
Oral contraceptives use, [figure omitted; refer to PDF] (%)	41 (50.6)
Current smoker, [figure omitted; refer to PDF] (%)	17 (21)
Family history of fractures [figure omitted; refer to PDF] (%)	6 (7.4)
Alcohol consumption 3 or more units/day, [figure omitted; refer to PDF] (%)	3 (3.37)
Duration of RA (years), median (range)	11.8 (1.5-35)
DAS-28 score, median (range)	3.8 (0-6.7)
Global functional status III-IV, [figure omitted; refer to PDF] (%)	7 (7.8)
HAQ-DI score, median (range)	0.7 (0-2.8)
Medications	[figure omitted; refer to PDF] (%)
Glucocorticoid	72 (88.9)
Prednisone dosage (mg/day), median (range)	4.4 (0-10)
Methotrexate	56 (70.9)
Sulfasalazine	39 (49.4)
Chloroquine	27 (34.2)
Azathioprine	25 (31.6)
D-Penicillamine	10 (12.7)
Biologic agents	8 (10.1)
BMD (WHO classification)	[figure omitted; refer to PDF] (%)
Normal	23 (33.3)
Osteopenia/osteoporosis	46 (66.7)

BMI: body mass index, RA: rheumatoid arthritis, DAS-28: disease activity score, HAQ-DI: health assessment questionnaire-disability index, BMD: bone mineral density, and WHO: World Health Organization.

RA patients were evaluated for BMD according to the classification criteria of the WHO for osteoporosis. According to the results of bone densitometry 23 patients (33.3%) had normal BMD, 31 (44.9%) had osteopenia (low bone density), and 15 (21.7%) had osteoporosis. For study purposes we grouped osteopenia with osteoporosis in one group constituted by 46 patients (66.7%) (see Table 2). From these, 23 patients (32.9%) had osteopenia and 13 (18.6%) had osteoporosis in the lumbar spine, whereas 35 patients (50.7%) had osteopenia and 7 (10.1%) had osteoporosis of the femoral neck. Clinical risk factors and factors related to RA severity for low BMD were registered and are listed in Table 2.

3.2. Molecular Analysis

We performed at least three PCR-RFLP tests for every subject included in the present study but were able to obtain complete genotyping results in 77 DNA samples from RA patients and 52 healthy subjects. Healthy subjects group was in Hardy-Weinberg equilibrium (data not shown).

In Table 3 genotype and allele frequencies of three SNPs in the OPG gene in healthy subjects and rheumatoid arthritis patients are shown; we found significant difference in the distribution of the C950T (rs2073617) alleles, the C allele being more prevalent in the RA group. Polymorphisms C209T and T245G did not reach statistical significance.

Table 3: Genotype and allele frequencies of polymorphisms in the TNFRSF11B (OPG) gene in healthy subjects and rheumatoid arthritis (RA) patients.

SNP	Genotypes	HS	RA	Allele	HS	RA	[figure omitted; refer to PDF]	OR
(95% CI)
C950Trs2073617	CC	9 (17.6)	28 (43.1)	C	47 (46.1)	82 (63.1)	0.01	0.50 (0.30-0.85)
CT	29 (56.9)	26 (40)	T	55 (53.9)	48 (36.9)
TT	13 (25.5)	11 (16.9)
C209Trs3134070	CC	39 (75.0)	52 (72.2)	C	89 (85.6)	122 (84.7)	1.0	1.1 (0.53-2.2)
CT	11 (21.2)	18 (25.0)	T	15 (14.4)	22 (15.3)
TT	2 (3.8)	2 (2.8)
T245Grs3134069	AA	41 (78.8)	66 (85.7)	A	93 (89.4)	142 (92.2)	0.51	0.71 (0.30-1.69)
AC	11 (21.2)	10 (13)	C	11 (10.6)	12 (7.8)
CC	0 (0)	1 (1.3)

OPG: osteoprotegerin, SNP: single nucleotide polymorphism, RA: rheumatoid arthritis, HS: healthy subjects, OR: odds ratios, and 95% CI: 95% confidence intervals. In order to compute OR (95% CI) the following alleles were used as reference: in correspondence with the allele of risk: rs2073617 T, rs3134070 T, and rs3134069 C.

We further analyzed genetic variants in RA patients and divided this group with respect to bone mass index status according to the WHO classification in normal or osteopenia/osteoporosis; there was no significant difference in allele distribution between RA patients with and without osteoporosis; Table 4 shows genotype and allele frequencies in these subgroups.

Table 4: Genotype and allele frequencies of polymorphisms in the TNFRSF11B (OPG) gene in rheumatoid arthritis (RA) patients with normal and low bone mineral density (BMD).

SNP	Genotypes	RA with normal BMD	RA with low BMD	Allele	RA with normal BMD	RA with low BMD	P	OR
(95% CI)
C950Trs2073617	CC	9 (50.0)	16 (44.4)	C	25 (69.4)	46 (63.9)	0.67	1.3 (0.55-3.03)
CT	7 (38.9)	14 (38.9)	T	11 (30.6)	26 (36.1)
TT	2 (11.1)	6 (16.7)
C209Trs3134070	CC	17 (81.0)	29 (69.0)	C	37 (88.1)	70 (83.3)	0.60	1.5 (0.49-4.43)
CT	3 (14.3)	12 (28.6)	T	5 (11.9)	14 (16.7)
TT	1 (4.8)	1 (2.4)
T245Grs3134069	AA	20 (90.9)	39 (88.6)	A	42 (95.5)	82 (93.2)	0.72	1.5 (0.29-7.95)
AC	2 (9.1)	4 (9.1)	C	2 (4.5)	6 (6.8)
CC	0 (0)	1 (2.3)

OPG: osteoprotegerin, SNP: single nucleotide polymorphism, RA: rheumatoid arthritis, BMD: bone mineral density, OR: odds ratios, and 95% CI: 95% confidence intervals. In order to compute OR (95% CI) the following alleles were used as reference: in correspondence with the allele of risk: rs2073617 T, rs3134070 T, and rs3134069 C.

4. Discussion

Osteoporosis represents a major health problem throughout the world and the most serious consequence of osteoporosis is hip fracture, which has a high associated morbidity and mortality [31]. Systemic effects of RA in bone remodeling include loss of axial and appendicular bone mass, associated with an increase in fracture risk (9.6 per 1,000 person-years) [32]. Clinical patterns of bone mass loss include generalized osteoporosis (axial and peripheral), juxta-articular osteoporosis which is adjacent to synovial membrane, and marginal and subchondral erosions which are directly associated with inflamed synovial tissue [33]. A high proportion of patients will suffer from hip and spine fractures causing limitations in physical abilities and social capabilities. Osteoporosis presence in RA patients involves multiple factors, such as disease activity, the use of some prescribed drugs particularly DMARDs and corticosteroids, female sex, older age, and menopause [34].

The development of erosive lesions as well as generalized osteoporosis in RA results from an imbalance between bone mass formation and bone resorption process. In RA the unsteadiness between osteoblastic bone formation (mesenchymal origin) and osteoclastic bone resorption (derived from monocytes/macrophages) promotes bone resorption by osteoclasts and as a consequence loss in bone strength, which may lead to the development of osteoporosis and fractures.

Multiple interactions between bone and the immune system have been described; bone is a major target of chronic inflammation, since it increases bone resorption and results in suppressed local bone formation, causing a wide spectrum of bone involvement in RA [35]. Lymphocyte or macrophage-derived cytokines are among the most potent mediators of inflammation, involved in upregulation of osteoclasts formation, activation, and signaling pathways [36].

Osteoclasts progenitors are recruited under physiological conditions from the haematopoyetic system and can get stimulated by some cytokines and hormones to form precursors that will eventually express cell surface markers of a well-differentiated osteoclast; furthermore, some inflammatory synovial cells, present in RA, particularly mononuclear cells and giant multinucleated cells can be differentiated into osteoclasts [5], under proinflammatory cytokines stimuli, including tumor necrosis factor alpha (TNF-α ), interleukin 1 (IL-1), and interleukin 6 (IL-6) [6].

Recently Pathak et al. described that sera from RA active patients contribute to bone loss by two mechanisms: (1) directly inhibiting osteoblasts proliferation and differentiation and (2) enhancing osteoblast-mediated osteoclastogenesis via RANKL and IL-6 [37].

Another key component that has influence over the osteoclastic differentiation is RANK, a receptor for TNF, attached to the membrane of osteoclast precursor cells that can bind to RNAKL through interaction with osteoblastic stromal cells [9, 38, 39]. RANK is essential for signal transduction that will lead to ostoclastic differentiation. Osteoprotegerin (OPG) also known as inhibition factor for osteoclastogenesis is a member of the TNFRSF, located in chromosome 8q24 [17] with a length of 29 Kb [13]; it is synthesized originally as a propeptide of 400 aa and is secreted as a soluble protein after the cleavage of 21 aa that corresponds to the transmembrane and cytoplasmic domains [17]. OPG mRNA is expressed in different human tissues (lung, heart, kidney, liver, intestines, stomach, cerebrum, thyroid gland, and bone marrow) apart from the bone, where its primary function appears to be inhibition of osteoclasts maturation and activation. In this study, we evaluated the association of three TNFRSF11B promoter polymorphisms with presence of osteoporosis in RA. Our study demonstrated that the C allele of the C950T was more frequent in RA compared with HS; thereafter we analyzed the RA group to seek for differences in allele distribution between patients with and without osteoporosis; nevertheless we did not observe a significant difference. The other two polymorphisms did not reach statistical significance.

The observed association in the present study could be caused by direct functional effect of the SNP or to linkage disequilibrium (LD) with another functional variant. According to the HapMap data, there is a certain amount of LD between these SNPs, indicating that although it is possible that different causative variants exist in TNFRSF11B , it is also possible that only one of them turns out to have functional consequences. The C950T SNP is located in the promoter region, near the TGF-β response area and 129 bp upstream the TATA box [40]. A functional study of TNFRSF11B polymorphisms through luciferase reporter assays revealed that gene expression was affected by the C950T polymorphism, since it was significantly higher in the presence of the C allele [27].

Previous association studies that include polymorphisms in the TNFRSF11B gene were performed in osteoporotic patients. In 2002 Ohmori et al. showed a significant association between C950T SNP and lower BMD [25]; in contrast with such results in 2004 Brandström et al. showed that C950T was not associated with bone mineral density in elderly Swedish women [41], and later Arko et al. in 2005 tested the three polymorphisms included in the present study and others searching for an association with bone mineral density in postmenopausal women; the linkage disequilibrium was confirmed but no associations were found [26]. For SNP T245G, similar frequencies as obtained in this study have been described previously by Zajickova et al. in postmenopausal women [42]. Only a few studies have looked for an association between these SNPs and RA; in 2010 Assmann et al. performed a case-control study and found an association between SNPs in the RANK and RANKL gene but not in the OPG gene; nevertheless none of the SNPs included in our study were genotyped by them [43]. In a recent study the SNP A163G in the osteoprotegerin gene was associated with osteoporosis in RA [44]. This polymorphism was not included in the present study but a lack of association with this particular polymorphism and osteoporosis in RA was previously confirmed by our study group [45].

Because this work constitutes a preliminary report, a limitation in our study is a small sample size that may limit our statistical power to detect small differences between groups although we consider that our sample would be sufficient to identify polymorphisms that may confer a high risk for low BMD in RA and did not show any trend to suspect that increasing the sample size a statistical significance would be achieved. We recognize that future multicenter studies should be performed, increasing the sample size, with our results being useful as a referent for Mexican patients with RA to identify that these polymorphisms probably do not confer a high risk for low BMD in our population with RA. Further studies in other populations are required to confirm our findings.

In conclusion, we confirmed the association of C allele of the C950T SNP in the TNFRSF11B gene with RA patients, but we were not able to support our main hypothesis that this polymorphism in the OPG gene is associated with osteoporosis in RA. Further studies including other regions of Mexico and Latin America with an increased sample size are required to confirm our findings and further elucidate if C950T SNP in the OPG gene could be associated with osteoporosis and not only with RA presence.

Acknowledgment

This work was supported by Grant from Instituto Mexicano del Seguro Social through FIS/IMSS/PROT/C2007/025.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] M. Akil, R. S. Amos, "Rheumatoid arthritis, I: clinical features and diagnosis," British Medical Journal , vol. 310, no. 6979, pp. 587-590, 1995.

[2] M. H. Cardiel, J. Rojas-Serrano, "Community based study to estimate prevalence, burden of illness and help seeking behavior in rheumatic diseases in Mexico City. A COPCORD study," Clinical and Experimental Rheumatology , vol. 20, no. 5, pp. 617-624, 2002.

[3] J. A. Kanis, L. J. Melton III, C. Christiansen, C. C. Johnston, N. Khaltaev, "The diagnosis of osteoporosis," Journal of Bone and Mineral Research , vol. 9, no. 8, pp. 1137-1141, 1994.

[4] F. Alcaraz-Lopez Miriam, G.-L. L. Aguilar-Chavez Erika, A. Lopez-Olivo Maria, C. Loaiza-Cardenas, I. Gamez-Nava Jorge, "Performance of albrand index for identifying low bone density in rheumatoid arthritis," The Journal of Rheumatology , vol. 33, no. 2, pp. 408, 2006.

[5] E. M. Gravallese, "Bone destruction in arthritis,", supplement 2 Annals of the Rheumatic Diseases , vol. 61, pp. ii84-ii86, 2002.

[6] S. Roux, P. Orcel, "Bone loss: factors that regulate osteoclast differentiation: an update," Arthritis Research , vol. 2, no. 6, pp. 451-456, 2000.

[7] D. L. Lacey, E. Timms, H.-L. Tan, M. J. Kelley, C. R. Dunstan, T. Burgess, R. Elliott, A. Colombero, G. Elliott, S. Scully, H. Hsu, J. Sullivan, N. Hawkins, E. Davy, C. Capparelli, A. Eli, Y.-X. Qian, S. Kaufman, I. Sarosi, V. Shalhoub, G. Senaldi, J. Guo, J. Delaney, W. J. Boyle, "Osteoprotegerin ligand is a cytokine that regulates osteoclast differentiation and activation," Cell , vol. 93, no. 2, pp. 165-176, 1998.

[8] N. Udagawa, N. Takahashi, H. Yasuda, A. Mizuno, K. Itoh, Y. Ueno, T. Shinki, M. T. Gillespie, T. J. Martin, K. Higashio, T. Suda, "Osteoprotegerin produced by osteoblasts is an important regulator in osteoclast development and function," Endocrinology , vol. 141, no. 9, pp. 3478-3484, 2000.

[9] H. Hsu, D. L. Lacey, C. R. Dunstan, I. Solovyev, A. Colombero, E. Timms, H.-L. Tan, G. Elliott, M. J. Kelley, I. Sarosi, L. Wang, X.-Z. Xia, R. Elliott, L. Chiu, T. Black, S. Scully, C. Capparelli, S. Morony, G. Shimamoto, M. B. Bass, W. J. Boyle, "Tumor necrosis factor receptor family member RANK mediates osteoclast differentiation and activation induced by osteoprotegerin ligand," Proceedings of the National Academy of Sciences of the United States of America , vol. 96, no. 7, pp. 3540-3545, 1999.

[10] M. C. Bezerra, J. F. Carvalho, A. S. Prokopowitsch, R. M. R. Pereira, "RANK, RANKL and osteoprotegerin in arthritic bone loss," Brazilian Journal of Medical and Biological Research , vol. 38, no. 2, pp. 161-170, 2005.

[11] D. H. Jones, Y.-Y. Kong, J. M. Penninger, "Role of RANKL and RANK in bone loss and arthritis,", supplement 2 Annals of the Rheumatic Diseases , vol. 61, pp. ii32-ii39, 2002.

[12] S.-S. Dong, X.-G. Liu, Y. Chen, Y. Guo, L. Wang, J. Zhao, D.-H. Xiong, X.-H. Xu, R. R. Recker, H.-W. Deng, "Association analyses of RANKL/RANK/OPG gene polymorphisms with femoral neck compression strength index variation in caucasians," Calcified Tissue International , vol. 85, no. 2, pp. 104-112, 2009.

[13] Y.-H. Hsu, T. Niu, H. A. Terwedow, X. Xu, Y. Feng, Z. Li, J. D. Brain, C. J. Rosen, N. Laird, X. Xu, "Variation in genes involved in the RANKL/RANK/OPG bone remodeling pathway are associated with bone mineral density at different skeletal sites in men," Human Genetics , vol. 118, no. 5, pp. 568-577, 2006.

[14] C. Liu, T. S. Walter, P. Huang, S. Zhang, X. Zhu, Y. Wu, L. R. Wedderburn, P. Tang, R. J. Owens, D. I. Stuart, J. Ren, B. Gao, "Structural and functional insights of RANKL-RANK interaction and signaling," Journal of Immunology , vol. 184, no. 12, pp. 6910-6919, 2010.

[15] A.-P. Trouvin, V. Goëb, "Receptor activator of nuclear factor-κ B ligand and osteoprotegerin: maintaining the balance to prevent bone loss," Clinical Interventions in Aging , vol. 5, pp. 345-354, 2010.

[16] J. B. Richards, F. K. Kavvoura, F. Rivadeneira, "Collaborative meta-analysis: associations of 150 candidate genes with osteoporosis and osteoporotic fracture," Annals of Internal Medicine , vol. 151, no. 8, pp. 528-537, 2009.

[17] P. Narducci, R. Bareggi, V. Nicolin, "Receptor Activator for Nuclear Factor kappa B Ligand (RANKL) as an osteoimmune key regulator in bone physiology and pathology," Acta Histochemica , vol. 113, no. 2, pp. 73-81, 2011.

[18] W. J. Boyle, W. S. Simonet, D. L. Lacey, "Osteoclast differentiation and activation," Nature , vol. 423, no. 6937, pp. 337-342, 2003.

[19] B. Arko, J. Preelj, R. Komel, A. Kocijani, P. Hudler, J. Marc, "Sequence variations in the osteoprotegerin gene promoter in patients with postmenopausal osteoporosis," Journal of Clinical Endocrinology and Metabolism , vol. 87, no. 9, pp. 4080-4084, 2002.

[20] F. C. Arnett, S. M. Edworthy, D. A. Bloch, D. J. McShane, J. F. Fries, N. S. Cooper, L. A. Healey, S. R. Kaplan, M. H. Liang, H. S. Luthra, T. A. Medsger Jr., D. M. Mitchell, D. H. Neustadt, R. S. Pinals, J. G. Schaller, J. T. Sharp, R. L. Wilder, G. G. Hunder, "The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis," Arthritis and Rheumatism , vol. 31, no. 3, pp. 315-324, 1988.

[21] M. L. L. Prevoo, M. A. van'T Hof, H. H. Kuper, M. A. van Leeuwen, L. B. A. van de Putte, P. L. C. M. van Riel, "Modified disease activity scores that include twenty-eight-joint counts development and validation in a prospective longitudinal study of patients with rheumatoid arthritis," Arthritis & Rheumatism , vol. 38, no. 1, pp. 44-48, 1995.

[22] M. H. Cardiel, M. Abello-Banfi, R. Ruiz-Mercado, D. Alarcon-Segovia, "How to measure health status in rheumatoid arthritis in non-English speaking patients: validation of a Spanish version of the Health Assessment Questionnaire Disability Index (Spanish HAQ-DI)," Clinical and Experimental Rheumatology , vol. 11, no. 2, pp. 117-121, 1993.

[23] M. C. Hochberg, R. W. Chang, I. Dwosh, S. Lindsey, T. Pincus, F. Wolfe, "The American College of Rheumatology 1991 revised criteria for the classification of global functional status in rheumatoid arthritis," Arthritis and Rheumatism , vol. 35, no. 5, pp. 498-502, 1992.

[24] S. Gustincich, G. Manfioletti, G. del Sal, C. Schneider, P. Carninci, "A fast method for high-quality genomic DNA extraction from whole human blood," BioTechniques , vol. 11, no. 3, pp. 298-302, 1991.

[25] H. Ohmori, Y. Makita, M. Funamizu, K. Hirooka, T. Hosoi, H. Orimo, T. Suzuki, K. Ikari, T. Nakajima, I. Inoue, A. Hata, "Linkage and association analyses of the osteoprotegerin gene locus with human osteoporosis," Journal of Human Genetics , vol. 47, no. 8, pp. 400-406, 2002.

[26] B. Arko, J. Prezelj, A. Kocijancic, R. Komel, J. Marc, "Association of the osteoprotegerin gene polymorphisms with bone mineral density in postmenopausal women," Maturitas , vol. 51, no. 3, pp. 270-279, 2005.

[27] C. Vidal, R. Formosa, A. Xuereb-Anastasi, "Functional polymorphisms within the TNFRSF11B (osteoprotegerin) gene increase the risk for low bone mineral density," Journal of Molecular Endocrinology , vol. 47, no. 3, pp. 327-333, 2011.

[28] G. Beyens, A. Daroszewska, F. de Freitas, E. Fransen, F. Vanhoenacker, L. Verbruggen, H.-G. Zmierczak, R. Westhovens, J. van Offel, S. H. Ralston, J.-P. Devogelaer, W. van Hul, "Identification of sex-specific associations between polymorphisms of the Osteoprotegerin gene, TNFRSF11B, and Paget's disease of bone," Journal of Bone and Mineral Research , vol. 22, no. 7, pp. 1062-1071, 2007.

[29] M. Soufi, M. Schoppet, A. M. Sattler, M. Herzum, B. Maisch, L. C. Hofbauer, J. R. Schaefer, "Osteoprotegerin gene polymorphisms in men with coronary artery disease," Journal of Clinical Endocrinology and Metabolism , vol. 89, no. 8, pp. 3764-3768, 2004.

[30] D. Pitocco, G. Zelano, G. Gioffrè, E. Di Stasio, F. Zaccardi, F. Martini, T. Musella, G. Scavone, M. Galli, S. Caputo, L. Mancini, G. Ghirlanda, "Association between osteoprotegerin G1181C and T245G polymorphisms and diabetic charcot neuroarthropathy: a case-control study," Diabetes Care , vol. 32, no. 9, pp. 1694-1697, 2009.

[31] C. Cooper, G. Campion, L. J. Melton III, "Hip fractures in the elderly: a world-wide projection," Osteoporosis International , vol. 2, no. 6, pp. 285-289, 1992.

[32] S. Y. Kim, S. Schneeweiss, J. Liu, G. W. Daniel, C.-L. Chang, K. Garneau, D. H. Solomon, "Risk of osteoporotic fracture in a large population-based cohort of patients with rheumatoid arthritis," Arthritis Research and Therapy , vol. 12, no. 4, article R154, 2010.

[33] S. R. Goldring, "The final pathogenetic steps in focal bone erosions in rheumatoid arthritis," Annals of the Rheumatic Diseases , vol. 59, no. 1, pp. i72-i74, 2000.

[34] R. M. R. Pereira, J. F. de Carvalho, E. Canalis, "Glucocorticoid-induced osteoporosis in rheumatic diseases," Clinics , vol. 65, no. 11, pp. 1197-1205, 2010.

[35] P. Geusens, W. F. Lems, "Osteoimmunology and osteoporosis," Arthritis Research and Therapy , vol. 13, no. 5, article 242, 2011.

[36] G. Mori, P. D'Amelio, R. Faccio, G. Brunetti, "The interplay between the bone and the immune system," Clinical and Developmental Immunology , vol. 2013, 2013.

[37] J. L. Pathak, N. Bravenboer, P. Verschueren, "Inflammatory factors in the circulation of patients with active rheumatoid arthritis stimulate osteoclastogenesis via endogenous cytokine production by osteoblasts," Osteoporosis International , vol. 25, no. 10, pp. 2453-2463, 2014.

[38] J. Li, I. Sarosi, X.-Q. Yan, "RANK is the intrinsic hematopoietic cell surface receptor that controls osteoclastogenesis and regulation of bone mass and calcium metabolism," Proceedings of the National Academy of Sciences of the United States of America , vol. 97, no. 4, pp. 1566-1571, 2000.

[39] T. L. Burgess, Y.-X. Qian, S. Kaufman, B. D. Ring, G. Van, C. Capparelli, M. Kelley, H. Hsu, W. J. Boyle, C. R. Dunstan, S. Hu, D. L. Lacey, "The ligand for osteoprotegerin (OPGL) directly activates mature osteoclasts," Journal of Cell Biology , vol. 145, no. 3, pp. 527-538, 1999.

[40] K. Thirunavukkarasu, R. R. Miles, D. L. Halladay, X. Yang, R. J. S. Galvin, S. Chandrasekhar, T. J. Martin, J. E. Onyia, "Stimulation of osteoprotegerin (OPG) gene expression by transforming growth factor-β (TGF-β ). Mapping of the OPG promoter region that mediates TGF-β effects," The Journal of Biological Chemistry , vol. 276, no. 39, pp. 36241-36250, 2001.

[41] H. Brandström, P. Gerdhem, F. Stiger, "Single nucleotide polymorphisms in the human gene for osteoprotegerin are not related to bone mineral density or fracture in elderly women," Calcified Tissue International , vol. 74, no. 1, pp. 18-24, 2004.

[42] K. Zajickova, A. Zemanova, M. Hill, I. Zofkova, "Is A163G polymorphism in the osteoprotegerin gene associated with heel velocity of sound in postmenopausal women?," Physiological Research , vol. 57, no. 1, pp. S153-S157, 2008.

[43] G. Assmann, J. Koenig, M. Pfreundschuh, J. T. Epplen, J. Kekow, K. Roemer, S. Wieczorek, "Genetic variations in genes encoding RANK, RANKL, and OPG in rheumatoid arthritis: a case-control study," Journal of Rheumatology , vol. 37, no. 5, pp. 900-904, 2010.

[44] Y. M. Hussien, A. Shehata, R. A. Karam, S. S. Alzahrani, H. Magdy, A. M. El-Shafey, "Polymorphism in vitamin D receptor and osteoprotegerin genes in Egyptian rheumatoid arthritis patients with and without osteoporosis," Molecular Biology Reports , vol. 40, no. 5, pp. 3675-3680, 2013.

[45] A. J. L. Brambila-Tapia, J. Duran-Gonzalez, L. Sandoval-Ramirez, J. P. Mena, M. Salazar-Paramo, J. I. Gamez-Nava, L. Gonzalez-Lopez, B. Lazalde-Medina B, N. O. Davalos, V. Peralta-Leal, M. V. Del Mercado, C. P. Beltran-Miranda, I. P. Davalos, "MTHFR C677T, MTHFR A1298C, and OPG A163G polymorphisms in Mexican patients with rheumatoid arthritis and osteoporosis," Disease Markers , vol. 32, no. 2, pp. 109-114, 2012.