Skip to main content

Association between a variation in the phosphodiesterase 4D gene and bone mineral density



Fragility fractures caused by osteoporosis are a major cause of morbidity and mortality in aging populations. Bone mineral density (BMD) is a useful surrogate marker for risk of fracture and is a highly heritable trait. The genetic variants underlying this genetic contribution are largely unknown.


We performed a large-scale association study investigating more than 25,000 single nucleotide polymorphisms (SNPs) located within 16,000 genes. Allele frequencies were estimated in contrasting DNA pools from white females selected for low (<0.87 g/cm2, n = 319) and high (> 1.11 g/cm2, n = 321) BMD at the lumbar spine. Significant findings were verified in two additional sample collections.


Based on allele frequency differences between DNA pools and subsequent individual genotyping, one of the candidate loci indicated was the phosphodiesterase 4D (PDE4D) gene region on chromosome 5q12. We subsequently tested the marker SNP, rs1498608, in a second sample of 138 white females with low (<0.91 g/cm2) and 138 females with high (>1.04 g/cm2) lumbar spine BMD. Odds ratios were 1.5 (P = 0.035) in the original sample and 2.1 (P = 0.018) in the replication sample. Association fine mapping with 80 SNPs located within 50 kilobases of the marker SNP identified a 20 kilobase region of association containing exon 6 of PDE4D. In a second, family-based replication sample with a preponderance of females with low BMD, rs1498608 showed an opposite relationship with BMD at different sites (p = 0.00044-0.09). We also replicated the previously reported association of the Ser37Ala polymorphism in BMP2, known to interact biologically with PDE4D, with BMD.


This study indicates that variants in the gene encoding PDE4D account for some of the genetic contribution to bone mineral density variation in humans. The contrasting results from different samples indicate that the effect may be context-dependent. PDE4 inhibitors have been shown to increase bone mass in normal and osteopenic mice, but up until now there have been no reports implicating any member of the PDE4 gene family in human osteoporosis.

Peer Review reports


The postmenopausal loss of bone mass and subsequent increased risk of low-energy (fragility) fractures is an important public health problem, especially in countries with a high proportion of elderly individuals. More than 1 million fragility fractures, primarily in postmenopausal women, occur each year in the US. The annual direct medical costs exceed US$10 billion [1]. Bone mineral density (BMD) measured with dual energy X-ray absorptiometry (DEXA) has been widely used to estimate the risk of fracture in epidemiological studies and to study treatment effects of antiresorptive agents in clinical trials. There are several well documented environmental and biological factors known to influence bone mineral density and the risk of fragility fractures including female gender, age, previous fragility fracture, low body weight, reduced lifetime exposure to estrogen, low calcium intake, physical inactivity, vitamin D deficiency, smoking, and excessive alcohol consumption [25]. There is also a strong genetic component to interindividual BMD variability, with heritability estimates ranging from 0.46 to 0.84 at different body sites [68]. Numerous candidate genes have been tested for association to BMD and fragility fractures. A polymorphism in a transcription factor-binding site of the collagen 1A1 (COL1A1) gene has shown one of the most consistent associations to osteoporosis, even if the association is generally weak for BMD and varies between populations [911]. Linkage studies have also been performed with the aim of locating genetic loci influencing BMD variation [1219]. So far, the genes responsible for the resulting linkage peaks have not been identified. Recently, linkage of a compound osteoporosis phenotype was reported to chromosome 20p12. Subsequent positional cloning efforts implicated BMP2, the gene encoding for bone morphogenetic protein 2, as responsible for the linkage [20]. Nevertheless, the associations reported thus far that have been independently validated account for only a small portion of the genetic contribution to BMD and osteoporosis.

Studies that rely on direct association approaches based on linkage disequilibrium within populations are expected to have greater statistical power and be more feasible to implement than traditional linkage studies to identify common variations that influence common, complex traits such as osteoporosis [21]. Recently, there has been increasing interest in the use of whole-genome association methods to identify genes that are involved in complex trait variation. To date, however, few such large-scale studies have been reported. In an effort to identify genes and variants that influence risk of osteoporosis, we conducted a large-scale study using more than 25,000 single nucleotide polymorphisms (SNPs) located within approximately 16,000 genes in DNA pools of unrelated females at the extremes of the lumbar spine bone mineral density distribution. A number of intriguing associations identified in this study are currently being scrutinized in further detail. In this paper we report the most advanced of these, which is the association of a variation in PDE4D, the gene encoding cyclic AMP-dependent phosphodiesterase 4D, with lumbar spine BMD. PDE4D selective inhibitors have been shown to promote osteoblast differentiation in progenitor cells [22] and to increase bone mass by promoting bone formation in normal mice [23] but the gene has not until now been implicated in human bone metabolism.



Discovery sample: unrelated females from UK twin collection

The population sample from which the discovery samples were chosen consisted of 5,436 female twins collected at the Twin Unit at St Thomas Hospital, London, England. They were selected without regard to health or trait. The volunteers had been recruited through advertisements and had undergone extensive investigation at the Twin Unit at St Thomas Hospital. Investigations included several questionnaires inquiring about present and past diseases, symptoms, family history, socio-economical factors, and medication. Subjects underwent an extensive clinical assessment including DEXA measurements of bone mineral density and anthropometric measurements [17, 24]. All individuals with data on lumbar spine BMD were considered for inclusion. To exclude relatives, the individual with the most extreme BMD was kept in each twin pair. Individuals with diseases or medication known to influence BMD were excluded, as were individuals younger than 40 years because of the observed complex relationship between age and BMD. In addition, individuals with fractures were excluded from the high BMD group. BMD values were adjusted for age, weight, BMI, self reported leisure time physical activity, smoking, and menopausal status using an ordinary least squares model including second and third order polynomial terms for age and second order terms for weight and BMI. We included BMI and weight as covariates because both were independently associated with BMD in this sample. Based on the trade-off between group sizes and separation, target sizes of 350 were chosen, resulting in a separation of approximately 1.9 SD between groups. After assessment of DNA availability sufficient for such a large scale study, group sizes were reduced to 319 and 321 individuals in the high and low groups, respectively. Lumbar spine BMD T-scores were calculated with the females between 20 and 35 years of age as the reference population. Based on this, 32% of the women in the low BMD group had osteoporosis and an additional 58% had osteopenia according to WHO criteria. The characteristics of the selected individuals are reported in Table 1.

Table 1 Characteristics of subjects in the groups selected for low and high lumbar spine bone mineral density.

Replication sample: Australian twin collection

A twin sample from Royal North Shore Hospital, Sidney, Australia was collected similarly as the UK twin collection. 731 individuals including twin pairs and singletons with lumbar spine BMD assessments were available for genotyping (Table 3). Groups of unrelated subjects corresponding to the lower and upper quartiles of the age- and BMI-adjusted lumbar spine BMD distribution were defined similarly to the discovery sample. The characteristics of the selected individuals are reported in Table 1.

Replication sample: international multi-center sib-pair study

The second replication sample was a multi-center (Australia, UK, New Zealand, Belgium) study that collected sib pairs concordant and discordant for bone mineral density [17]. Probands (BMD Z-score <-1.5 at lumbar spine, femoral neck, or hip) were identified and their siblings were contacted and underwent DEXA measurements at the lumbar spine and hip. Participants had to be between 25–85 years of age. Exclusion criteria included steroid medication, hyperparathyroidism, immobility, amenorrhea, anorexia nervosa, and unstable thyroid disease. Nine hundred and eight individual samples were genotyped from 392 families. In the present analysis we included only females older than 40 years of age. In this sample the distribution of family sizes were 164 singletons, 248 families of 2, 34 families of 3, 7 families of 4, and 3 families of 5 members (Table 3). Lumbar spine BMD levels were adjusted by age and BMI as described for the Australian sample.

Human subjects protection

All studies were approved by the appropriate research ethics committees. All participants gave their informed consent to participate in genetic studies before enrollment.

Bone mineral density

Bone mineral density was estimated from the L1-L4 vertebrae, hip, and forearm using DEXA according to the user's manual for the Hologic QDR 4500W, (Hologic, Waltham, Massachusetts, United States) at all collection sites.

SNP markers and genotyping

A set of 25,494 SNP markers was selected from a collection of 125,799 experimentally validated polymorphic variations [25]. This set was limited to SNPs located within gene coding regions, minor allele frequencies greater than 0.02 (95% have frequencies greater than 0.1), and a target inter-marker spacing of 40 kb. SNP annotation is based on NCBI dbSNP database, refSNP, build 118. Genomic annotation is based on NCBI Genome Build 34. Gene annotation is based on LocusLink genes for which NCBI was providing positions on the Mapview FTP site.

For pooled DNA assays, 25 ng of case and control DNA pools was used for amplification at each site. All PCR and MassEXTEND™ reactions were conducted using standard conditions [26]. Relative allele frequency estimates were derived from area under the peak calculations of mass spectrometry measurements from four analyte aliquots as described elsewhere [26]. For individual genotyping, the same procedure was applied except only 2.5 ng DNA was used and only one mass spectrometry measurement was taken.

Primers used to genotype rs1498608 were ATAACCTCGGGGTCCAGAAA (forward PCR primer), GAATCCCTGTTCATTCCTTG (reverse PCR primer) and CCCTAAAAACTGTTCCAGGTA (extension primer). The primers used to genotype the Ser37Ala polymorphism in BMP2 were AGCTGGGCCGCAGGAAGTTCG (forward PCR primer), TCGTCAGAGGGCTGGGATGAG (reverse PCR primer) and TGAGGGGCGGCCCGACG (extension primer).

Statistical analysis

Tests of association between adjusted lumbar spine BMD group and each SNP using pooled DNA were carried out in a similar fashion as explained elsewhere [27]. Briefly, the test statistic is based on the difference in allele frequencies between the two groups divided by the known sources of variation in each allele frequency estimate, including sampling and pool-specific measurement variation. Sources of measurement variation incorporated in the test statistic are pool formation, PCR/mass extension, and chip measurement. When three or more replicate measurements of a SNP were available within a model level, the corresponding variance component was estimated from the data. Otherwise, the following historical laboratory averages were used: pool formation = 5.0 × 10-5, PCR/mass extension = 1.7 × 10-4, and chip measurement = 1.0 × 10-4.

Tests of association using individual genotypes were carried out using a chi-square test of heterogeneity to compare allele frequencies, and Fisher's exact test to compare genotype frequencies (due to low frequency contingency table cells). Confidence intervals and P-values for odds ratios were derived using Fisher's exact test. When one or more cell counts were zero, non-infinite odds ratios were estimated by adding 0.5 to each cell [28]. In the samples that included a combination of singletons, sib pairs, and occasionally additional relatives, we estimated the relationship between genotypes and phenotypes using the generalized estimating equations (GEE) approach with a Gaussian link by clustering on family using an exchangeable correlation matrix [29]. Hypothesis testing was carried out with a Wald test statistic. The geepack implementation of GEE in the R statistical software platform was used [30]. No attempt was made to correct P-values for multiple testing. Rather, P-values are provided to compare the relative strength of association. P-values less than 0.05 are referred to as statistically significant.


Initial genome scan in UK sample

We carried out a genome-wide association study using 25,494 SNPs located within 10 kb of 15,995 LocusLink annotated genes. An overview of the investigative process is shown in Figure 1. The basic design was a two-group study in subjects from the tails of the adjusted lumbar spine BMD distribution. We selected lumbar spine BMD as the phenotype to create the contrasting groups because it had a high estimate of heritability in our twin sample (h2 = 0.82) [7]. The selected low and high BMD groups consisted of 319 and 321 individuals, respectively. The adjusted BMD range was 0.56–0.87 g/cm2 in the low BMD group and 1.11–1.60 g/cm2 in the high BMD group, corresponding to the upper and lower 22nd percentiles. Other characteristics of the samples are described in Table 1. To facilitate the screening of such a large number of SNPs, we employed a high-throughput approach using DNA pools, chip-based mass spectrometry [26, 3133], and a three-phase SNP selection strategy (Figure 1). In the first phase, we performed a single PCR and primer extension reaction for each SNP on two DNA pools consisting of equimolar amounts of DNA from each individual in the low BMD group and high BMD group, respectively. Relative allele frequencies obtained from four mass spectrometry measurements of the extension products were compared between pools. In the second phase, the 1,520 SNPs (~6%) with the most statistically significant associations (nominal P-values < 0.05) were re-measured in triplicate on each DNA pool. In the third phase, we genotyped the 140 most significant SNPs (9%) from step two (nominal P-values < 0.02) on all individuals comprising the pools. Based on the genotype results, 78 SNPs were confirmed to have statistically significant allele frequency differences between cases and controls (P < 0.05). The liberal criteria for selecting SNPs from each phase represent a practical trade-off between following up false positive versus overlooking false negative associations. We chose to follow up as many SNPs from each phase as seemed reasonable.

Figure 1
figure 1

Overview of DNA pool-based large-scale association strategy. Phases 1 and 2 are conducted using DNA pools yielding allele frequencies, all subsequent steps involve genotyping of individual samples. hME, homogeneous MassEXTEND; MS, mass spectrometry. See text for more details.

One of the associations was found with rs1498608, an A/T polymorphism within intron 5 of PDE4D on chromosome 5q12. Allele frequencies for the T allele based on genotyping were 0.91 in the low lumbar spine BMD pool and 0.88 in the high BMD pool (OR = 1.5, P = 0.035). Complete genotype counts and summary statistics are reported in Table 2. Observed genotype frequencies were consistent with expected frequencies under Hardy-Weinberg equilibrium within each study collection. Menopausal status did not have a significant influence on the effect (P = 0.87).

Table 2 Allele and genotype frequencies of rs1498608 in high and low lumbar spine BMD groups.
Table 3 Characteristics of complete genotyped replication samples.

Genome-wide studies using tens of thousands of SNPs and liberal statistical selection criteria are expected to yield a high proportion of false positive associations. To distinguish the true genetic effects from among the false positives, the 78 selected SNPs were genotyped in a second twin sample from Australia.

Replication in Australian sample

The Australian replication sample, a combination of female and male twin pairs and singletons, was analyzed in two ways. First, to create a design and carry out an analysis comparable to the discovery sample, unrelated individuals were selected from the lower and upper quartiles of the sex-specific adjusted lumbar spine BMD distribution (Table 1). A similar effect was observed for the marker SNP rs1498608 in females (OR = 2.14, P = 0.018) and males (OR = 1.55, P = 0.35) as in the original UK collection (Table 2). The second method of analysis utilized generalized estimating equation (GEE) models to take into account all of the available genotype information by carrying out a regression-type analysis while accounting for familial covariance. The characteristics of this sample are reported in Table 3. The regression of marker genotypes on adjusted lumbar spine BMD, with sex included as a covariate, found the AA genotype to be associated with significantly higher levels than the AT (β = 7.8 g/cm2, P = 0.049) or TT (β = 8.0 g/cm2, P = 0.037) genotypes, thus confirming the results observed in the unrelated tails of this sample. Similar GEE analyses carried out for femoral neck and hip BMD were not statistically significant.

Replication in international multi-center family study

Being a sample of mostly affected sib pairs, this sample was unsuitable for formation of groups with contrasting BMD because of the preponderance of individuals with low BMD (Table 3). Therefore, we restricted the analysis to using a generalized estimating equation, regressing marker genotypes on BMD values. Surprisingly, the estimates in this sample were opposite to that in the Australian sample, as the AA genotype was associated with lower adjusted lumbar spine BMD values than both the AT (β = -5.3 g/cm2, P = 0.11) and the TT (β = -5.4 g/cm2, P = 0.09) genotypes. Using Z-scores at the femoral neck (P = 0.0007 and 0.0004), total hip (P = 0.003 and 0.007), and lumbar spine (P = 0.03 and 0.02) as dependent variables confirmed this pattern of association. In all cases the AT and TT genotypes had very similar point estimates.

Association fine mapping

In order to better define the extent of the region of association and possibly identify other SNPs more strongly associated with BMD, we performed DNA pool based association fine mapping in the UK sample using 80 publicly available intronic SNPs in the 100 kb region surrounding the incident SNP (Figure 2). This analysis identified a 20 kb region of association encompassing exon 6 of PDE4D.

Figure 2
figure 2

Association fine mapping of PDE4D region. Eighty public domain SNPs in a 100-kb window around the initial marker SNP (indicated as a bold symbol at the center of the plot) were compared between the UK low and high BMD pools. Ten of 80 SNPs were significant at P = 0.05 (horizontal dashed line). The x-axis corresponds to the chromosomal position of each SNP, the y-axis to the test P-values (--log10 scale). The continuous dark line presents the results of a goodness-of-fit test for an excess of significant associations (at a 5% significance level) in a 10-kb sliding window assessed at 1-kb increments. The continuous light gray line is the result of a nonlinear smoothing function showing a weighted average of the P-values across the region. The color of each point corresponds to the minor allele frequency of each SNP in the control sample (see legend below graph). Vertical dashed lines are placed every 20 kb. The LocusLink gene annotation for NCBI genome build 34 shows the location of exon six. All investigated SNPs are located within introns.

Replication of BMP2association

As described in the discussion below, PDE4D inhibition is known to influence BMP2-induced alkaline phosphatase activity in osteoblast precursor cells. Recently, variation in the gene encoding BMP2 was found to be associated with osteoporosis in a study employing whole genome linkage and subsequent positional cloning [20]. Since we were unaware of any published independent attempts to replicate this finding, we genotyped the Ser37Ala polymorphism in our UK and international samples. In the UK sample, the allele frequency of the rare allele was 2.2% in the low BMD group and 1.4% in the high BMD group, with an odds ratio of 1.6 (P = 0.28). We tested for, but were unable to detect, an interaction between the Ser37Ala polymorphism and rs1498608 on the association with lumbar spine BMD. In the international sample we performed an allele based general estimating equation to estimate the effect of the rare allele on BMD in that sample. The allele frequency of Ala37 in this sample, mainly selected for low BMD, was 1.9%. The effect of the Ala allele was estimated to decrease the adjusted lumbar spine BMD by 0.06 g/cm2 (P = 0.0029). There were no homozygous Ala individuals in this sample.


In an association study using SNPs in nearly 16,000 genes we obtained evidence that variation in the SNP rs1498608 located within PDE4D is associated with low bone mineral density at the lumbar spine in females. PDE4D encodes cyclic AMP-dependent phosphodiesterase 4D. Phosphodiesterases are a superfamily of enzymes involved in degradation of cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) [34, 35]. cAMP and cGMP are important second messengers participating in the response of various cells to hormones. In osteoblasts, cAMP produced in response to parathyroid hormone or prostaglandins regulates osteoblastic differentiation [3639], which leads to increases in cancellous bone volume as indicated by experiments in animal models [4045]. The intracellular level of cAMP is regulated by G protein-coupled adenylyl cyclase [46], and degradation is mediated by the phosphodiesterases. The phosphodiesterase superfamily consists of seven families, PDE1-7, distinguished by substrate specificity, chromatographic behaviour during purification, and affinity for biochemical activators and inhibitors. Of these, the PDE4 family is specific for cAMP and is selectively inhibited by rolipram. Four PDE4 genes, 4A, 4B, 4C, and 4D, have been cloned from rat and humans, all of which are predicted to have multiple protein products due to alternate spicing of RNAs. PDE4 inhibitors have been shown to increase bone formation in normal mice [23] and to ameliorate loss of bone mass in animal models of osteopenia [47, 48]. PDE4A and PDE4D are expressed in two common mouse osteoblastic cell lines, ST2 and MC3T3-E1, that represent different stages in the osteoblast differentiation pathway [22]. PDE4 inhibition with rolipram increased BMP2-induced alkaline phosphatase activity, a marker of early osteoblast differentiation in ST2 cells. Furthermore, rolipram increased the expression of alkalic phosphatase, osteopontin, collagen type I and osteocalcin in the same osteoblast precursor cells [22]. In spite of these experimental data, we are not aware of any published attempts to investigate the role of PDE4 genes in human osteoporosis. However, variation in PDE4D was recently reported to be associated with the risk of ischemic stroke [49]. Given the central role of PDE4 in second messenger signalling, it is quite conceivable that PDE4D variants may have effects on the risk for different common diseases. There are other examples of genes having such pleiotropic effects, the most notable being APOE in hyperlipidemia and Alzheimer's disease [50, 51]. It should also be noted that Gretarsdottir et al found that the PDE4D association with stroke was strongest for a region in the recently extended 5' end of the gene, which is close to 1,000 kb upstream of rs1498608 [49]. Assuming a contribution of PDE4D to the risk of osteoporosis as well as stroke, it is possible that different domains are involved in the different diseases.

Given the interaction between BMP2 and PDE4 for the inhibition on osteoblastic differentiation in vitro, it is interesting to note that variants in the gene encoding for BMP2 have also been found to increase risk of osteoporosis in humans [20]. In the current study, we replicated the association between the Ser37Ala variant in BMP2 and measures of osteoporosis in an international family-based sample ascertained via low BMD probands. Although not statistically significant, this finding was supported by the results in the discovery sample of unrelated high and low spine BMD subjects. The allele frequency in the low BMD group was 2.2% and in the high group 1.4%, with an odds ratio of 1.6 (P = 0.28). The rare allele was less common in our low BMD group than the low spine BMD group (3%) in the Icelandic sample. However, our allele frequencies in the low and high BMD groups and the resulting OR corresponded well with the figures in the Danish sample (1.8% vs 1.0%, RR = 1.8) reported in the same paper [20]. We found no evidence for statistical interaction between the variations in BMP2 and PDE4D in either sample. However, given the low minor allele frequencies of each SNP, there was very little power to test for interaction effects.

The starting point of the present study was a large-scale association study of more than 25,000 SNPs located in 16,000 genes. After a stepwise selection process an association between an intronic SNP in PDE4D and lumbar spine bone mineral density was detected, providing the first evidence that a variant of this gene could contribute to the risk of osteoporosis in humans. The effect was similar in size in premenopausal and postmenopausal women, indicating that the effect would be on the attainment of peak bone mass rather than the rate of decrease in BMD after menopause. The lack of a detectable interaction with female sex hormones is supported by having observed a similar genetic effect in the small sample of males in our study. The genetic contribution to peak bone mass is possibly bigger and definitely better documented than the as yet unproven genetic influence on postmenopausal bone loss [52], and it is possible that PDE4D could contribute to this effect, especially in light of the documented anabolic effect on bone by PDE4 inhibitors.

An association with an intronic SNP provides little evidence for a change in amount or function of the protein that could explain the association. None of the 80 SNPs investigated as part of the association fine mapping were non-synonymous coding changes, which is consistent with the extraordinary lack of variation that others have reported for all PDE classes [53] and PDE4D in particular [49]. This makes it unlikely (but still possible) that the observed association would be due to a non-synonymous and disruptive single-base coding change in linkage disequilibrium with our marker SNP. Therefore it is more likely that the effect is mediated by a change in RNA splicing or expression.

Given the functional similarity between different PDE4 enzymes, we went back and scrutinized our data for associations with SNPs in the other PDE4 genes that may have been overlooked during the first stage of the scan. The only SNP in PDE4B in our assay set, rs1318475, was taken through to the second stage (Figure 1) where it was estimated to have an OR of 0.78 (P = 0.041), but failed the criteria to be taken forward to the genotyping stage. Similarly, a SNP roughly 18 kb downstream of PDE4C, rs874628, was also taken forward to the second stage where it displayed an OR of 1.3 (P = 0.08). These results suggest that further investigation into possible associations between variants in all PDE4 genes and bone mineral density may be justified.

The route by which these genetic associations were arrived at and the potential for spurious association must be considered. Recent published work has brought light to the need for proper validation to verify genetic findings for complex traits [5456]. In the current study, the initial association found between the PDE4D marker and bone mineral density was one result from over 25,000 hypothesis tests. A conservative Bonferroni adjustment to yield an experiment-wide type I error rate of 0.05 would demand a test-wise p-value on the order of 10-6. Given the modest sample size, only common variations with relatively large effects (OR > 2) would reach such significance levels. Instead, we chose to be more mindful of the role of type II error rates and apply a more liberal set of criteria in the initial phases of the study and verify true genetic effects by independent replication. The analysis of 78 selected markers in the Australian replication sample resulted in multiple associations of continuing interest, with rs1498608 displaying one of the strongest associations. A one-sided test of association comparing the results in the discovery and replication samples yields a p-value of 0.0074. This would not be considered significant on an experiment-wide level after Bonferroni adjustment, which would require a p-value on the order of 0.0006 or lower. The analysis in the international replication sample produced contradictory data in that the A allele, which in the first two samples was associated with increased lumbar spine bone mineral density, was associated with decreased BMD at all tested sites. The pattern of association evident from the first two samples, with AT and TT genotypes having very similar point estimates, was preserved in this sample, even in the face of the reverse direction of association. The highly statistically significant association between rs1498608 and femoral neck and hip BMD in this third sample and the consistency in the pattern of association would be unexpected from a spurious result. A possible explanation for the contradictory results could be the fact that the international sample consists mostly of individuals with low BMD since the probands all have a BMD Z-score < -1.5, and most of the siblings also have low BMD. It is possible that within such a selected sample the relationship between rs1498608 and BMD could be altered due to interactions with other genetic or environmental factors.

The well-documented anabolic effect on bone by PDE4 inhibitors lends indirect support for the association reported here, and it would seem that the possible role of PDE4D variants in the genetic contribution to BMD in humans merits further investigation.


The result of the present large scale association study together with data from previously published animal models suggest that genetic variation in the gene encoding PDE4D contributes to the variation in lumbar spine BMD in humans.


  1. Ray NF, Chan JK, Thamer M, Melton LJ: Medical expenditures for the treatment of osteoporotic fractures in the United States in 1995: report from the National Osteoporosis Foundation. J Bone Miner Res. 1997, 12: 24-35.

    CAS  Article  PubMed  Google Scholar 

  2. Cummings SR, Nevitt MC, Browner WS, Stone K, Fox KM, Ensrud KE, Cauley J, Black D, Vogt TM: Risk factors for hip fracture in white women. Study of Osteoporotic Fractures Research Group. N Engl J Med. 1995, 332: 767-773. 10.1056/NEJM199503233321202.

    CAS  Article  PubMed  Google Scholar 

  3. Espallargues M, Sampietro-Colom L, Estrada MD, Sola M, del Rio L, Setoain J, Granados A: Identifying bone-mass-related risk factors for fracture to guide bone densitometry measurements: a systematic review of the literature. Osteoporos Int. 2001, 12: 811-822. 10.1007/s001980170031.

    CAS  Article  PubMed  Google Scholar 

  4. Kanis JA: Diagnosis of osteoporosis and assessment of fracture risk. Lancet. 2002, 359: 1929-1936. 10.1016/S0140-6736(02)08761-5.

    Article  PubMed  Google Scholar 

  5. Jordan KM, Cooper C: Epidemiology of osteoporosis. Best Pract Res Clin Rheumatol. 2002, 16: 795-806. 10.1053/berh.2002.0264.

    CAS  Article  PubMed  Google Scholar 

  6. Arden NK, Baker J, Hogg C, Baan K, Spector TD: The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: a study of postmenopausal twins. J Bone Miner Res. 1996, 11: 530-534.

    CAS  Article  PubMed  Google Scholar 

  7. Hunter DJ, de Lange M, Andrew T, Snieder H, MacGregor AJ, Spector TD: Genetic variation in bone mineral density and calcaneal ultrasound: a study of the influence of menopause using female twins. Osteoporos Int. 2001, 12: 406-411. 10.1007/s001980170110.

    CAS  Article  PubMed  Google Scholar 

  8. Pocock NA, Eisman JA, Hopper JL, Yeates MG, Sambrook PN, Eberl S: Genetic determinants of bone mass in adults. A twin study. J Clin Invest. 1987, 80: 706-710.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Efstathiadou Z, Kranas V, Ioannidis JP, Georgiou I, Tsatsoulis A: The Sp1 COLIA1 gene polymorphism, and not vitamin D receptor or estrogen receptor gene polymorphisms, determines bone mineral density in postmenopausal Greek women. Osteoporos Int. 2001, 12: 326-331. 10.1007/s001980170123.

    CAS  Article  PubMed  Google Scholar 

  10. Efstathiadou Z, Tsatsoulis A, Ioannidis JP: Association of collagen Ialpha 1 Sp1 polymorphism with the risk of prevalent fractures: a meta-analysis. J Bone Miner Res. 2001, 16: 1586-1592.

    CAS  Article  PubMed  Google Scholar 

  11. Mann V, Hobson EE, Li B, Stewart TL, Grant SF, Robins SP, Aspden RM, Ralston SH: A COL1A1 Sp1 binding site polymorphism predisposes to osteoporotic fracture by affecting bone density and quality. J Clin Invest. 2001, 107: 899-907.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Devoto M, Shimoya K, Caminis J, Ott J, Tenenhouse A, Whyte MP, Sereda L, Hall S, Considine E, Williams CJ, Tromp G, Kuivaniemi H, Ala-Kokko L, Prockop DJ, Spotila LD: First-stage autosomal genome screen in extended pedigrees suggests genes predisposing to low bone mineral density on chromosomes 1p, 2p and 4q. Eur J Hum Genet. 1998, 6: 151-157. 10.1038/sj.ejhg.5200169.

    CAS  Article  PubMed  Google Scholar 

  13. Niu T, Chen C, Cordell H, Yang J, Wang B, Wang Z, Fang Z, Schork NJ, Rosen CJ, Xu X: A genome-wide scan for loci linked to forearm bone mineral density. Hum Genet. 1999, 104: 226-233. 10.1007/s004390050940.

    CAS  Article  PubMed  Google Scholar 

  14. Koller DL, Econs MJ, Morin PA, Christian JC, Hui SL, Parry P, Curran ME, Rodriguez LA, Conneally PM, Joslyn G, Peacock M, Johnston CC, Foroud T: Genome screen for QTLs contributing to normal variation in bone mineral density and osteoporosis. J Clin Endocrinol Metab. 2000, 85: 3116-3120. 10.1210/jc.85.9.3116.

    CAS  PubMed  Google Scholar 

  15. Deng HW, Xu FH, Liu YZ, Shen H, Deng H, Huang QY, Liu YJ, Conway T, Li JL, Davies KM, Recker RR: A whole-genome linkage scan suggests several genomic regions potentially containing QTLs underlying the variation of stature. Am J Med Genet. 2002, 113: 29-39. 10.1002/ajmg.10742.

    Article  PubMed  Google Scholar 

  16. Karasik D, Myers RH, Cupples LA, Hannan MT, Gagnon DR, Herbert A, Kiel DP: Genome screen for quantitative trait loci contributing to normal variation in bone mineral density: the Framingham Study. J Bone Miner Res. 2002, 17: 1718-1727.

    CAS  Article  PubMed  Google Scholar 

  17. Wilson SG, Reed PW, Bansal A, Chiano M, Lindersson M, Langdown M, Prince RL, Thompson D, Thompson E, Bailey M, Kleyn PW, Sambrook P, Shi MM, Spector TD: Comparison of genome screens for two independent cohorts provides replication of suggestive linkage of bone mineral density to 3p21 and 1p36. Am J Hum Genet. 2003, 72: 144-155. 10.1086/345819.

    CAS  Article  PubMed  Google Scholar 

  18. Econs MJ, Koller DL, Hui SL, Fishburn T, Conneally PM, Johnston CCJ, Peacock M, Foroud TM: Confirmation of linkage to chromosome 1q for peak vertebral bone mineral density in premenopausal white women. Am J Hum Genet. 2004, 74: 223-228. 10.1086/381401.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. Peacock M, Koller DL, Hui S, Johnston CC, Foroud T, Econs MJ: Peak bone mineral density at the hip is linked to chromosomes 14q and 15q. Osteoporos Int. 2004, 15: 489-496. 10.1007/s00198-003-1560-7.

    Article  PubMed  Google Scholar 

  20. Styrkarsdottir U, Cazier JB, Kong A, Rolfsson O, Larsen H, Bjarnadottir E, Johannsdottir VD, Sigurdardottir MS, Bagger Y, Christiansen C, Reynisdottir I, Grant SF, Jonasson K, Frigge ML, Gulcher JR, Sigurdsson G, Stefansson K: Linkage of Osteoporosis to Chromosome 20p12 and Association to BMP2. PLoS Biol. 2003, 1: E69-10.1371/journal.pbio.0000069.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Risch NJ: Searching for genetic determinants in the new millennium. Nature. 2000, 405: 847-856. 10.1038/35015718.

    CAS  Article  PubMed  Google Scholar 

  22. Wakabayashi S, Tsutsumimoto T, Kawasaki S, Kinoshita T, Horiuchi H, Takaoka K: Involvement of phosphodiesterase isozymes in osteoblastic differentiation. J Bone Miner Res. 2002, 17: 249-256.

    CAS  Article  PubMed  Google Scholar 

  23. Kinoshita T, Kobayashi S, Ebara S, Yoshimura Y, Horiuchi H, Tsutsumimoto T, Wakabayashi S, Takaoka K: Phosphodiesterase inhibitors, pentoxifylline and rolipram, increase bone mass mainly by promoting bone formation in normal mice. Bone. 2000, 27: 811-817. 10.1016/S8756-3282(00)00395-1.

    CAS  Article  PubMed  Google Scholar 

  24. Andrew T, Mak YT, Reed P, MacGregor AJ, Spector TD: Linkage and association for bone mineral density and heel ultrasound measurements with a simple tandem repeat polymorphism near the osteocalcin gene in female dizygotic twins. Osteoporos Int. 2002, 13: 745-754. 10.1007/s001980200102.

    CAS  Article  PubMed  Google Scholar 

  25. Nelson MR, Marnellos G, Kammerer S, Hoyal CR, Shi MM, Cantor CR, Braun A: Large-scale validation of single nucleotide polymorphisms in gene regions. Genome Res. 2004, 14: 1664-1668. 10.1101/gr.2421604.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Bansal A, van den Boom D, Kammerer S, Honisch C, Adam G, Cantor CR, Kleyn P, Braun A: Association testing by DNA pooling: an effective initial screen. Proc Natl Acad Sci U S A. 2002, 99: 16871-16874. 10.1073/pnas.262671399.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. Barratt BJ, Payne F, Rance HE, Nutland S, Todd JA, Clayton DG: Identification of the sources of error in allele frequency estimations from pooled DNA indicates an optimal experimental design. Ann Hum Genet. 2002, 66: 393-405. 10.1046/j.1469-1809.2002.00125.x.

    CAS  Article  PubMed  Google Scholar 

  28. Agresti A: Categorical data analysis. 2002, New York, Wiley-Interscience, xv, 710-2nd

    Chapter  Google Scholar 

  29. Slager SL, Schaid DJ, Wang L, Thibodeau SN: Candidate-gene association studies with pedigree data: controlling for environmental covariates. Genet Epidemiol. 2003, 24: 273-283. 10.1002/gepi.10228.

    CAS  Article  PubMed  Google Scholar 

  30. Yan J, Fine J: Estimating equations for association structures. Stat Med. 2004, 23: 859-74; discussion 875-7,879-80. 10.1002/sim.1650.

    Article  PubMed  Google Scholar 

  31. Kammerer S, Burns-Hamuro LL, Ma Y, Hamon SC, Canaves JM, Shi MM, Nelson MR, Sing CF, Cantor CR, Taylor SS, Braun A: Amino acid variant in the kinase binding domain of dual-specific A kinase-anchoring protein 2: a disease susceptibility polymorphism. Proc Natl Acad Sci U S A. 2003, 100: 4066-4071. 10.1073/pnas.2628028100.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. Buetow KH, Edmonson M, MacDonald R, Clifford R, Yip P, Kelley J, Little DP, Strausberg R, Koester H, Cantor CR, Braun A: High-throughput development and characterization of a genomewide collection of gene-based single nucleotide polymorphism markers by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Proc Natl Acad Sci U S A. 2001, 98: 581-584. 10.1073/pnas.021506298.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. Mohlke KL, Erdos MR, Scott LJ, Fingerlin TE, Jackson AU, Silander K, Hollstein P, Boehnke M, Collins FS: High-throughput screening for evidence of association by using mass spectrometry genotyping on DNA pools. Proc Natl Acad Sci U S A. 2002, 99: 16928-16933. 10.1073/pnas.262661399.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Manganiello VC, Murata T, Taira M, Belfrage P, Degerman E: Diversity in cyclic nucleotide phosphodiesterase isoenzyme families. Arch Biochem Biophys. 1995, 322: 1-13. 10.1006/abbi.1995.1429.

    CAS  Article  PubMed  Google Scholar 

  35. Beavo JA: Cyclic nucleotide phosphodiesterases: functional implications of multiple isoforms. Physiol Rev. 1995, 75: 725-748.

    CAS  PubMed  Google Scholar 

  36. Farndale RW, Sandy JR, Atkinson SJ, Pennington SR, Meghji S, Meikle MC: Parathyroid hormone and prostaglandin E2 stimulate both inositol phosphates and cyclic AMP accumulation in mouse osteoblast cultures. Biochem J. 1988, 252: 263-268.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Kumegawa M, Ikeda E, Tanaka S, Haneji T, Yora T, Sakagishi Y, Minami N, Hiramatsu M: The effects of prostaglandin E2, parathyroid hormone, 1,25 dihydroxycholecalciferol, and cyclic nucleotide analogs on alkaline phosphatase activity in osteoblastic cells. Calcif Tissue Int. 1984, 36: 72-76.

    CAS  Article  PubMed  Google Scholar 

  38. Ishizuya T, Yokose S, Hori M, Noda T, Suda T, Yoshiki S, Yamaguchi A: Parathyroid hormone exerts disparate effects on osteoblast differentiation depending on exposure time in rat osteoblastic cells. J Clin Invest. 1997, 99: 2961-2970.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. Partridge NC, Bloch SR, Pearman AT: Signal transduction pathways mediating parathyroid hormone regulation of osteoblastic gene expression. J Cell Biochem. 1994, 55: 321-327.

    CAS  Article  PubMed  Google Scholar 

  40. Jee WS, Ueno K, Deng YP, Woodbury DM: The effects of prostaglandin E2 in growing rats: increased metaphyseal hard tissue and cortico-endosteal bone formation. Calcif Tissue Int. 1985, 37: 148-157.

    CAS  Article  PubMed  Google Scholar 

  41. Jee WS, Ueno K, Kimmel DB, Woodbury DM, Price P, Woodbury LA: The role of bone cells in increasing metaphyseal hard tissue in rapidly growing rats treated with prostaglandin E2. Bone. 1987, 8: 171-178. 10.1016/8756-3282(87)90017-2.

    CAS  Article  PubMed  Google Scholar 

  42. High WB: Effects of orally administered prostaglandin E-2 on cortical bone turnover in adult dogs: a histomorphometric study. Bone. 1987, 8: 363-373. 10.1016/8756-3282(87)90068-8.

    CAS  Article  PubMed  Google Scholar 

  43. Whitfield JF, Morley P: Small bone-building fragments of parathyroid hormone: new therapeutic agents for osteoporosis. Trends Pharmacol Sci. 1995, 16: 382-386. 10.1016/S0165-6147(00)89079-3.

    CAS  Article  PubMed  Google Scholar 

  44. Reeve J: PTH: a future role in the management of osteoporosis?. J Bone Miner Res. 1996, 11: 440-445.

    CAS  Article  PubMed  Google Scholar 

  45. Finkelstein JS, Klibanski A, Schaefer EH, Hornstein MD, Schiff I, Neer RM: Parathyroid hormone for the prevention of bone loss induced by estrogen deficiency. N Engl J Med. 1994, 331: 1618-1623. 10.1056/NEJM199412153312404.

    CAS  Article  PubMed  Google Scholar 

  46. Casperson GF, Bourne HR: Biochemical and molecular genetic analysis of hormone-sensitive adenylyl cyclase. Annu Rev Pharmacol Toxicol. 1987, 27: 371-384. 10.1146/

    CAS  Article  PubMed  Google Scholar 

  47. Miyamoto K, Waki Y, Horita T, Kasugai S, Ohya K: Reduction of bone loss by denbufylline, an inhibitor of phosphodiesterase 4. Biochem Pharmacol. 1997, 54: 613-617. 10.1016/S0006-2952(97)00211-6.

    CAS  Article  PubMed  Google Scholar 

  48. Waki Y, Horita T, Miyamoto K, Ohya K, Kasugai S: Effects of XT-44, a phosphodiesterase 4 inhibitor, in osteoblastgenesis and osteoclastgenesis in culture and its therapeutic effects in rat osteopenia models. Jpn J Pharmacol. 1999, 79: 477-483. 10.1254/jjp.79.477.

    CAS  Article  PubMed  Google Scholar 

  49. Gretarsdottir S, Thorleifsson G, Reynisdottir ST, Manolescu A, Jonsdottir S, Jonsdottir T, Gudmundsdottir T, Bjarnadottir SM, Einarsson OB, Gudjonsdottir HM, Hawkins M, Gudmundsson G, Gudmundsdottir H, Andrason H, Gudmundsdottir AS, Sigurdardottir M, Chou TT, Nahmias J, Goss S, Sveinbjornsdottir S, Valdimarsson EM, Jakobsson F, Agnarsson U, Gudnason V, Thorgeirsson G, Fingerle J, Gurney M, Gudbjartsson D, Frigge ML, Kong A, Stefansson K, Gulcher JR: The gene encoding phosphodiesterase 4D confers risk of ischemic stroke. Nat Genet. 2003, 35: 131-138. 10.1038/ng1245.

    CAS  Article  PubMed  Google Scholar 

  50. Ghiselli G, Gregg RE, Zech LA, Schaefer EJ, Brewer HBJ: Phenotype study of apolipoprotein E isoforms in hyperlipoproteinaemic patients. Lancet. 1982, 2: 405-407. 10.1016/S0140-6736(82)90439-1.

    CAS  Article  PubMed  Google Scholar 

  51. Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, Small GW, Roses AD, Haines JL, Pericak-Vance MA: Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science. 1993, 261: 921-923.

    CAS  Article  PubMed  Google Scholar 

  52. Peacock M, Turner CH, Econs MJ, Foroud T: Genetics of osteoporosis. Endocr Rev. 2002, 23: 303-326. 10.1210/er.23.3.303.

    CAS  Article  PubMed  Google Scholar 

  53. Houslay MD, Adams DR: PDE4 cAMP phosphodiesterases: modular enzymes that orchestrate signalling cross-talk, desensitization and compartmentalization. Biochem J. 2003, 370: 1-18. 10.1042/BJ20021698.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN: Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003, 33: 177-182. 10.1038/ng1071.

    CAS  Article  PubMed  Google Scholar 

  55. Ioannidis JP, Ntzani EE, Trikalinos TA, Contopoulos-Ioannidis DG: Replication validity of genetic association studies. Nat Genet. 2001, 29: 306-309. 10.1038/ng749.

    CAS  Article  PubMed  Google Scholar 

  56. Trikalinos TA, Ntzani EE, Contopoulos-Ioannidis DG, Ioannidis JP: Establishment of genetic associations for complex diseases is independent of early study findings. Eur J Hum Genet. 2004, 12: 762-769. 10.1038/sj.ejhg.5201227.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

Download references


The authors acknowledge R. Prince, N. Gilchrist, J.Y. Reginster, I. Fogelman, I. Smith, and the clinical teams who contributed patients to the study.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Andreas Braun.

Additional information

Competing interests

The authors of this manuscript affiliated with Sequenom, Inc. (RR, SM, CH, GM, SK, MN, AB) declare competing financial interests, which may include current or prior receipt of salary and/or stock ownership, as Sequenom, Inc. may be affected financially by the publication of this manuscript.

Authors' contributions

RR drafted the manuscript and participated in study design and data analysis. SM was study project leader and participated in data analysis. SK supervised the operational aspects of the study, and was responsible for the development of the SNP marker assay set. CH and GM participated in the development of the SNP marker assay set. SW was project manager for the international multicenter study. PS was principal investigator for the Australian Twin collection. TS was principal investigator for the UK twin and international collections. MN participated in study design, performed the statistical analyses and helped draft the manuscript. AB had the overall scientific responsibility for the study.

Richard H Reneland, Steven Mah contributed equally to this work.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Reneland, R.H., Mah, S., Kammerer, S. et al. Association between a variation in the phosphodiesterase 4D gene and bone mineral density. BMC Med Genet 6, 9 (2005).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Bone Mineral Density
  • Generalize Estimate Equation
  • Lumbar Spine Bone Mineral Density
  • PDE4 Inhibitor
  • Replication Sample