Association of ADIPOR1 polymorphisms with bone mineral density in postmenopausal Korean women.

Adiponectin may affect bone through interactions with two known receptors, adiponectin receptors (ADIPOR) 1 and 2. We examined the association between polymorphisms of ADIPOR1 and ADIPOR2 and bone mineral density (BMD) in postmenopausal Korean women. Six polymorphisms in ADIPOR1 and four polymorphisms in ADIPOR2 were selected and genotyped in all study participants (n = 1,329). BMD at the lumbar spine and femur neck were measured using dual-energy X-ray absorptiometry. Lateral thoracolumbar (T4-L4) radiographs were obtained for vertebral fracture assessment and the occurrence of non-vertebral fractures examined using self-reported data. P values were adjusted for multiple testing using Bonferroni correction (P(corr)). ADIPOR1 rs16850799 and rs34010966 polymorphisms were significantly associated with femur neck BMD (P(corr) = 0.036 in the dominant model; P(corr) = 0.024 and P(corr) = 0.006 in the additive and dominant models, respectively). Subjects with the rare allele of each polymorphism had lower BMD, and association of rs34010966 with BMD showed a gene dosage effect. However, ADIPOR2 single nucleotide polymorphisms and haplotypes were not associated with BMD at any site. Our results suggest that ADIPOR1 polymorphisms present a useful genetic marker for BMD in postmenopausal Korean women.

Introduction

Osteoporosis is a systemic bone disease characterized by low bone mineral density (BMD) and subsequent bone loss, leading to increased risk of fracture (1993). Obesity is strongly correlated with increased bone mineral density (BMD) (Felson et al., 1993), and increase in body weight reduces fracture risk in both genders (De Laet et al., 2005). In particular, the effects of body weight are possibly attributed to both fat mass and lean mass, of which fat mass is more important in postmenopausal women (Reid et al., 1992). Mechanical load forces may contribute to this relationship, along with other factors, such as sex hormones, glucocorticoids and insulin (Reid et al., 1993; Gennari et al., 2004). Recent studies have suggested that the association of BMD with body weight may also be mediated by hormonal factors secreted by adipocytes, including leptin and adiponectin (Zoico et al., 2003; Jurimae and Jurimae 2007).

Adiponectin is an adipocyte-derived hormone that possibly affects bone. The receptors for adiponectin, ADIPOR1 and ADIPOR2, have been identified on both osteoblasts and osteoclasts (Berner et al., 2004; Shinoda et al., 2006). To date, reports on the effects of adiponectin on bone metabolism have been inconsistent. Adiponectin has been shown to increase osteoblast proliferation and differentiation and inhibit osteoclastogenesis in vitro (Luo et al., 2005; Oshima et al., 2005; Williams et al., 2009). In support of these in vitro data, transient overexpression of adiponectin in mice increased trabecular bone mass and inhibited osteoclast number and bone resorption (Oshima et al., 2005), although conflicting results have also been reported (Shinoda et al., 2006). In contrast to in vitro and animal data, clinical studies have shown an inverse association of adiponectin with BMD in perimenopausal women (Jurimae et al., 2005), diabetic men and women (Lenchik et al., 2003), and elderly men (Basurto et al., 2009).These findings suggest that the receptors, but not adiponectin itself, are important for bone metabolism.

Although the multiple risk factors influence BMD and the pathogenesis of osteoporosis, genetic factors are mainly implicated and account for approximately 50% to 85% of the variance in BMD based on twin and family studies (Slemenda et al., 1991; Arden and Spector, 1997). Several genetic epidemiological studies have demonstrated an association of adiponectin and adiponectin receptor polymorphisms with type 2 diabetes and its related phenotypes (Damcott et al., 2005; Stefan et al., 2005). Recently, Lee et al. (2006) demonstrated that the T45G polymorphism of the adiponectin gene is significantly linked to lower lumbar spine BMD in Korean women. However, the issue of whether ADIPOR1 and ADIPOR2 polymorphisms are associated with bone metabolism is yet to be established. In the current study, we investigated the associations of genetic variations in ADIPOR1 and ADIPOR2 with bone mineral density (BMD) and the risk of osteoporotic fracture in postmenopausal women.

Results

Clinical data and correlations between BMD and age, weight, height, and YSM are listed in Supplemental Data Table S1. The mean age of participants was 59.1 ± 7.3 yr (range 45-87 yr), and mean YSM was 9.7 ± 7.8 yr (range, 1-42 yr). As expected, age and YSM were inversely correlated with BMD at both the lumbar spine and femur neck regions. Weight and height were positively correlated with BMD at both sites.

On the basis of the direct sequencing of DNA from 24 Korean individuals, we identified 18 and 17 sequence variants of ADIPOR1 and ADIPOR2, respectively (Figure 1). Among the identified polymorphisms, six SNPs from ADIPOR1 [rs2275737, rs16850799, rs34010966, rs33942950, rs1342387, rs34559546] and four SNPs in ADIPOR2 [rs1029629, -63442G > C, rs12342, rs1044471] were selected for larger-scale genotyping based on minor allele frequency (MAF ≥ 0.1), LDs, and haplotype-tagging status. The genotype frequencies of all SNPs analyzed are shown in Supplemental Data Table S2. All genotype distributions were in Hardy-Weinberg equilibrium (P > 0.05). For ADIPOR1, two haplotype blocks were constructed using Haploview version 3.2 (Barrett et al., 2005) and Lewontin's method, while one block was constructed for ADIOPR2. Three common haplotypes (frequency > 0.1) in each block were investigated in detail.

Figure 1

Gene map of the (A) ADIPOR1 and (B) ADIPOR2 genes. Coding exons are marked with black blocks, and 5'- and 3'-UTRs with white blocks. The first base of the translation start site is designated nucleotide '+1'. Asterisks (*) indicate polymorphisms genotyped in a larger population (n = 1329).

Next, the association of ADIPOR1 and ADIPOR2 polymorphisms with BMD at the lumbar spine and femur neck was analyzed. In linear regression analysis adjusted for confounding variables, two SNPs of ADIPOR1 were significantly associated with BMD at the femur neck, even after Bonferroni correction was strictly adopted for multiple comparisons (Table 1). Specifically, ADIPOR1 rs16850799 and rs34010966 were significantly associated with femur neck BMD (P = 0.036 in the dominant model; P = 0.024 and P = 0.006 in the additive and dominant model, respectively). The effects of ADIPOR1 rs16850799 on BMD were gene dose-dependent. Specifically, subjects homozygous for the common rs16850799 allele displayed highest BMD (0.738 ± 0.126 g/cm2), while heterozygotes had intermediate BMD (0.726 ± 0.119 g/cm2) and rare allele homozygotes had the lowest BMD values (0.722 ± 0.117 g/cm2). For ADIPOR1 rs34010966, subjects with rare alleles had lower BMD, compared with those with common alleles. Moreover, ADIPOR1 rs34010966 was correlated with BMD at the total femur and Ward`s triangle (P = 0.030 in the dominant model and P = 0.018 in the dominant model, respectively) (Table 2). However, the ADIPOR2 polymorphisms and haplotypes were not associated with BMD at any site.

Table 1

Regression analysis of BMD at the lumbar spine and femoral neck in relation to ADIPOR1 and ADIPOR2 polymorphisms in postmenopausal Korean women (n = 1329)

The number of subjects and means and standard deviation of BMD are shown.

*C/C, C/R, and R/R represent homozygotes for the common allele, and heterozygotes and homozygotes for the rarer allele, respectively.

†Pa, Pb, and Pc are P values of additive, dominant and recessive models for multiple regression analysis, respectively.

‡P values after Bonferroni correction.

Table 2

Regression analysis of BMD at proximal femur sites in relation to ADIPOR1 polymorphisms in postmenopausal Korean women

The number of subjects and means and standard deviation of BMD are shown.

*C/C, C/R, and R/R represent homozygotes for the common allele, and heterozygotes and homozygotes for the rarer allele, respectively.

†Pa, Pb, and Pc are P values of additive, dominant and recessive models for multiple regression analysis, respectively.

‡P values after Bonferroni correction.

The genetic effects of ADIPOR1 and ADIPOR2 polymorphisms on risk of osteoporosis fracture were analyzed using multiple logistic regression analysis adjusted for confounding variables (Table 3). Vertebral and non-vertebral fractures were observed in 99 and 73 subjects, respectively. Despite the significant association of ADIPOR1 +348T > C and ADIPOR1 rs34010966 with femur neck BMD, no association between these SNPs and osteoporotic fracture was evident. None of the ADIPOR2 polymorphisms or haplotypes were associated with increased risk of any type of osteoporotic fracture.

Table 3

Logistic regression analysis of ADIPOR1 and ADIPOR2 polymorphisms in relation to any fracture risk in Korean postmenopausal women

Genotype distributions and P values for logistic analyses of three alternative models (additive, dominant and recessive), controlling for age, weight, height, and YSM, as covariates, are shown.

*P values after Bonferroni correction.

Discussion

In the present study, we focused on the genetic association of ADIPOR1 and ADIPOR2 SNPs and haplotypes with postmenopausal osteoporosis-related phenotypes. Our findings indicate that ADIPOR1 rs16850799 and rs34010966 are significantly associated with femur neck BMD. Interestingly, ADIPOR1 rs34010966 was associated with BMD at the total femur and Ward`s triangle (P = 0.030 in the dominant model and P = 0.018 in the dominant model, respectively). However, ADIPOR2 SNPs and haplotypes were not associated with BMD at any site or with any type of osteoporotic fracture. To our knowledge, this is the first clinical report supporting a role for ADIPOR1 and ADIPOR2 polymorphisms on BMD and fracture risk.

Two SNPs in the intron of ADIPOR1 were associated with femur neck BMD, but not lumbar spine BMD. This site-specific variation of BMD heritability is frequently reported in candidate gene association studies (Yerges et al., 2009). These results may be attributed to different genetic mechanisms of BMDs at different skeletal sites (Videman et al., 2007). At present, the mechanisms through which variants at the ADIPOR1 locus influence bone phenotype are hypothetical. ADIPOR1 +348 T > C and +2852 C > G are located within the intron and therefore do not induce amino acid changes. Therefore, it is unclear whether these SNPs are functional or in linkage disequilibrium with an unidentified polymorphism. Moreover, we cannot exclude the possibility that the non-coding intronic polymorphisms mediate genetic function through changes in the alternative splicing regimen (Ast, 2004). Notably, SNPs of ADIPOR1, but not those of ADIPOR2, were associated with proximal femur BMDs in our study. These findings suggest that ADIPOR1 is more important for bone biology that ADIPOR2, although both receptors are expressed in bone cells.

Alterations in the expression or conformational changes of the adiponectin receptor associated with SNPs may impair the direct effects of adiponectin on bone. Recently, Soccio et al. (2006) reported that specific ADIPOR1 SNPs linked to increased cardiovascular risk are associated with 30-40% lower ADIPOR1 mRNA levels in blood mononuclear cells and adipose tissue biopsies. Moreover, the increase in insulin resistance due to impaired adiponectin signaling may indirectly affect bone metabolism, since insulin is a potential regulator of bone (Hickman and McElduff, 1989). Insulin exerts an anabolic effect in bone through direct effects on osteoblast proliferation (Cornish et al., 1996) and indirect effects on the production of sex hormones and their binding globulin (Reid, 2008). A recent study on non-diabetic Mexican Americans has described a strong, positive correlation between ADIPOR1 expression levels in skeletal muscle and insulin sensitivity, as determined with the glucose clamp (Civitarese et al., 2004). Thus, the association between ADIPOR1 polymorphisms and bone phenotype may be mediated by an increase in insulin resistance secondary to decreased ADIPOR1 expression resulting from polymorphisms.

Despite the significant association of ADIPOR1 rs16850799 and rs34010966 with BMD at the femur neck, SNPs were not associated with risk of fracture. There may be several explanations for this finding. Firstly, although genetic risk factors are evidently important in the etiology of fracture and BMD (Nguyen et al., 2000), these are likely to be specifically linked to each phenotype. Secondly, some of the additional genetic variants predictive of fracture may affect fracture risk through BMD-independent mechanisms, such as effects on bone geometry, bone matrix and other features of bone quality. Thirdly, fall-related environmental factors, such as postural balance (Pajala et al., 2004), muscle function (Tiainen et al., 2005), and cognitive abilities (Wright et al., 2001), play critical roles in determining fracture risk. Genetic variants related to these factors may contribute to fracture through increased fall risk. Finally, the small number of fractured subjects was insufficient to demonstrate statistical power of association between genotypes and fracture risk. Therefore, a possible role of ADIPOR1 polymorphisms as a genetic marker for bone metabolism cannot be excluded.

There are several potential limitations of this study. First, our study population was restricted to persons who visited the university hospital. These subjects may not be representative of the general population, impossibly contributing to selection bias. Second, despite the association of the ADIPOR1 rs16850799 and rs34010966 with femur neck BMD, neither was associated with risk of non-vertebral fracture. However, the heritability of fracture itself has been estimated to lie between 25% and 35% (MacGregor et al., 2000), which is much lower than the hereditability of the BMD values (Arden and Spector, 1997). In addition, the genetic component of fracture may ultimately be a combination of polygenic effects, gene-gene and genetic environmental interactions. Further multiple analyses with other possible candidate genes are essential to resolve these complex relationships. Third, we did not assess the functional effects of ADIPOR1 SNP activity, and cannot confirm whether adiponectin receptor activity varies according to the absence or presence of polymorphisms and haplotypes. Therefore, we cannot definitely assert that the genotypes are functionally relevant.

In summary, to establish the possible involvement of genetic polymorphisms of ADIPOR1 and ADIPOR2 in osteoporosis, 10 common sites were genotyped in Korean postmenopausal women. ADIPOR1 rs16850799 and rs34010966 were significantly associated with femur neck BMD, even after Bonferroni correction. Our findings suggest that ADIPOR1 polymorphisms constitute one of the genetic determinants of BMD in postmenopausal Korean women.

Methods

Subjects

The study population comprised 1,329 postmenopausal women of Korean ethnicity who visited Asan Medical Center (AMC, Seoul, Korea). All subject visited our hospital for diagnosis or treatment of osteoporosis spontaneously. Menopause was defined as the absence of menstruation for at least one year and serum follicle-stimulating hormone (FSH) > 30 IU/l. Women with premature menopause (before 40 yr of age) were excluded, along with those taking drugs with a possible effect on bone metabolism for more than 6 months or within the previous 12 months (such as glucocorticoids, sex hormones, bisphosphonate or other treatments for osteoporosis). Subjects suffering from diseases that could affect bone metabolism were additionally excluded. Women with a history of stroke or dementia were excluded from analysis owing to concerns related to their limited physical activity, in addition to patients with osteophyte formation above the fourth grade of the Nathan classification (Nathan et al., 1994) and/or severe facet joint osteoarthritis in the lumbar spine diagnosed using conventional spine radiography. The study was approved by the AMC ethics review committee, and written informed consent was obtained from all subjects.

BMD measurements

Areal BMD (g/cm2) of the anterior-posterior lumbar spine (L2-L4) and femur neck was measured in 834 women with dual energy X-ray absorptiometry using Lunar equipment (Lunar, Expert XL with software version 1.90; Madison, WI). In the remaining 495 women, BMD was measured using Hologic equipment (Hologic, QDR 4500-A with software version 4.84, Waltham, MA). Owing to upper extremity dominance, BMD at the proximal femur was measured at non-dominant sites. Short-term in vivo measurement precision of the Lunar and Hologic machines, expressed as coefficient of variation, were 0.82% and 0.85% for the lumbar spine, and 1.12% and 1.20% for the femur neck, respectively. These values were obtained by scanning 17 volunteers who were not part of the study. Each volunteer underwent five scans on the same day, getting on and off the table between examinations. To derive cross-calibration equations between the two systems, BMD values were measured at the lumbar and femur neck with the two machines in 109 healthy Korean women (55±11 yr, range 31-75 yr), and calculated as follows (Jo et al.,1999)

L2-L4 BMD (g/cm2): Lunar = 1.1287 × Hologic - 0.0027

Femur neck BMD (g/cm2): Lunar = 1.1556 × Hologic - 0.0182

We obtained additional BMD values at other proximal femur sites taken after January 2001. The Hologic machine did not measure BMD at the femur shaft. Therefore, BMD values at the femur shaft and other proximal sites (total femur, trochanter and Ward`s triangle) were available for 571 and 893 participants, respectively.

Detection of vertebral and non-vertebral fractures

We examined prevalent morphological vertebral fracture in all study subjects by obtaining lateral thoracolumbar (T4-L4) radiographs. Vertebral fractures were assessed in accordance with the recommendations of the Working Group on Vertebral Fractures (Kiel, 1995). Radiographs were assessed at AMC by expert radiologists blinded to this study A vertebral fracture was defined quantitatively as more than a 20% reduction in any measured vertebral height (i.e., anterior, middle, or posterior (Genant et al., 1993). In addition, a history of non-vertebral fracture, including those of hip, wrist, forearm, humerus, rib and pelvis, was obtained using self-administered questionnaires. Fractures that had clearly been caused by major trauma, such as a traffic accident or fall from higher than standing height, were excluded.

Sequencing analysis of ADIPOR1 and ADIPOR2

We sequenced all exons, including exon-intron boundaries, and the promoter region (~1.5 Kb) to detect single nucleotide polymorphisms (SNPs) in 24 Korean DNA samples using the ABI PRISM 3730 DNA analyzer (Applied Biosystems, Foster City, CA). Sixteen primer sets were designed for amplification and sequencing analyses based on GenBank sequences (Reference Genome Sequence, ADIPOR1: NT_004671.15, and ADIPOR2: NT_009759.15). Sequence variants were verified with automated sequencing chromatograms. SNPs were detected by multiple sequence alignment using the Phred/Phrap/Consed package and Polyphred (Ewing et al., 1998; Gordon et al., 1998).

SNP Genotyping

For genotyping of polymorphic sites, amplification primers and MGB probes for TaqMan (Livak, 1999) were designed using Primer Express (Applied Biosystems). One allelic probe was labeled with FAM dye and the other with the fluorescent dye, VIC. PCR was performed using TaqMan Universal Master mix without UNG (Applied Biosystems) with 900 nM primer and 200 nM TaqMan MGB-probe. The reactions were performed in a 384-well format in a total reaction volume of 5 µl using 20 ng genomic DNA. The plates were then placed in a thermal cycler (PE 9700, Applied Biosystems) and heated at 50℃ for 2 min and 95℃ for 10 min, followed by 40 cycles of 95℃ for 15 s and 60℃ for 1 min. The TaqMan assay plates were then transferred to a Prism 7900HT instrument (Applied Biosystems) and the fluorescence intensity in each well of the plate were read. Fluorescence data files from each plate were analyzed using automated software (SDS 2.1, Applied Biosystems). Genotyping quality control was performed in 44 samples by duplicate checking, and the rate of concordance in duplicates was 100%.

Statistical analysis

To determine whether individual variants were in equilibrium at each locus in the population (Hardy-Weinberg equilibrium), χ2 tests were applied. We examined Lewontin's D' (|D'|) and the linkage disequilibrium (LD) coefficient (r2) between all pairs of biallelic loci. Haploview version 3.2 (Whitehead Institute for Biomedical research, Cambridge, MA) was used for the structure of LD block (Barrett et al., 2005). This program uses two-marker expectation maximization to estimate the maximum-likelihood values of the four gamete frequencies from which D' and log of odds (LOD) values are derived. Haplotypes and phase probabilities of all polymorphic sites for haplotypes were calculated for each individual with PHASE software (ver 2.0) using the algorithm developed by Stephens et al. (2001). Individuals with phase probabilities of less than 97% were excluded from the analysis. The genetic effects of inferred haplotypes were analyzed in a similar way to polymorphisms.

Multiple regression analyses were performed for BMD, controlling for age, yr since menopause (YSM), weight and height as covariates. The genotype and haplotype distributions between subjects with and without vertebral and non-vertebral fractures were additionally analyzed with a logistic regression model controlling for age, YSM, weight and height. Genotypes were assigned codes of 0, 1, and 2 for the additive model, 0, 1, and 1 for the dominant model, and 0, 0, and 1 for the recessive model. P values were adjusted for multiple testing by Bonferroni correction (P). All statistical analyses were conducted using SAS (SAS Institute, Cary, NC).