Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Bone health in children with inflammatory bowel disease : adjusting for bone age

OBJECTIVES: Clinical results of bone mineral density for children with inflammatory bowel disease are commonly reported using reference data for chronological age. It is known that these children, particularly those with Crohn disease, experience delayed growth and maturation. Therefore, it is more appropriate to compare clinical results with bone age rather than chronological age. MATERIALS AND METHODS: Areal bone mineral density (aBMD) was measured using dual energy x-ray absorptiometry, and bone age was assessed using the Tanner-Whitehouse 3 method from a standard hand/wrist radiograph. Results were available for 44 children ages 7.99 to 16.89 years. Areal bone mineral density measurements were converted to z scores using both chronological and bone ages for each subject. RESULTS: Areal bone mineral density z scores calculated using bone age, as opposed to chronological age, were significantly improved for both the total body and lumbar spine regions of interest. When subjects were grouped according to diagnosis, bone age generated z scores remained significantly improved for those with Crohn disease but not for those diagnosed with ulcerative colitis. Grouping of children with Crohn disease into younger and older ages produced significantly higher z scores using bone age compared with chronological for the older age group, but not the younger age group. CONCLUSIONS: Our findings, in accordance with those presented in the literature, suggest that aBMD results in children with Crohn disease should include the consideration of bone age, rather than merely chronological age. Bone size, although not as easily available, would also be an important consideration for interpreting results in paediatric populations. © 2009 by European Society for Pediatric Gastroenterology, Hepatology, and Nutrition and North American Society for Pediatric Gastroenterology, Hepatology, and Nutrition.

'Sink or swim': An evaluation of the clinical characteristics of individuals with high bone mass

Summary High bone mineral density on routine dual energy X-ray absorptiometry (DXA) may indicate an underlying skeletal dysplasia. Two hundred fifty-eight individuals with unexplained high bone mass (HBM), 236 relatives (41% with HBM) and 58 spouses were studied. Cases could not float, had mandible enlargement, extra bone, broad frames, larger shoe sizes and increased body mass index (BMI). HBM cases may harbour an underlying genetic disorder. Introduction High bone mineral density is a sporadic incidental finding on routine DXA scanning of apparently asymptomatic individuals. Such individuals may have an underlying skeletal dysplasia, as seen in LRP5 mutations. We aimed to characterize unexplained HBM and determine the potential for an underlying skeletal dysplasia. Methods Two hundred fifty-eight individuals with unexplained HBM (defined as L1 Z-score ≥ +3.2 plus total hip Z-score ≥ +1.2, or total hip Z-score ≥ +3.2) were recruited from 15 UK centres, by screening 335,115 DXA scans. Unexplained HBM affected 0.181% of DXA scans. Next 236 relatives were recruited of whom 94 (41%) had HBM (defined as L1 Z-score + total hip Z-score ≥ +3.2). Fifty-eight spouses were also recruited together with the unaffected relatives as controls. Phenotypes of cases and controls, obtained from clinical assessment, were compared using random-effects linear and logistic regression models, clustered by family, adjusted for confounders, including age and sex. Results Individuals with unexplained HBM had an excess of sinking when swimming (7.11 [3.65, 13.84], p < 0.001; adjusted odds ratio with 95% confidence interval shown), mandible enlargement (4.16 [2.34, 7.39], p < 0.001), extra bone at tendon/ligament insertions (2.07 [1.13, 3.78], p = 0.018) and broad frame (3.55 [2.12, 5.95], p < 0.001). HBM cases also had a larger shoe size (mean difference 0.4 [0.1, 0.7] UK sizes, p = 0.009) and increased BMI (mean difference 2.2 [1.3, 3.1] kg/m 2, p < 0.001). Conclusion Individuals with unexplained HBM have an excess of clinical characteristics associated with skeletal dysplasia and their relatives are commonly affected, suggesting many may harbour an underlying genetic disorder affecting bone mass.

Estrogen receptor α regulates area-adjusted bone mineral content in late pubertal girls

Context: Whether the action of estrogen in skeletal development depends on estrogen receptor α as encoded by the ESR1 gene is unknown. Objectives: The aim of this study was to establish whether the gain in area-adjusted bone mineral content (ABMC) in girls occurs in late puberty and to examine whether the magnitude of this gain is related to ESR1 polymorphisms. Design: We conducted a cross-sectional analysis. Setting: The study involved the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based prospective study. Participants: Participants included 3097 11-yr-olds with DNA samples, dual x-ray absorptiometry measurements, and pubertal stage information. Outcomes: Outcome measures included separate prespecified analyses in boys and girls of the relationship between ABMC derived from total body dual x-ray absorptiometry scans and Tanner stage and of the interaction between ABMC, Tanner stage, and ESR1 polymorphisms. Results: Total body less head and spinal ABMC were higher in girls in Tanner stages 4 and 5, compared with those in Tanner stages 1, 2, and 3. In contrast, height increased throughout puberty. No differences were observed in ABMC according to Tanner stage in boys. For rs2234693 (PvuII) and rs9340799 (XbaI) polymorphisms, differences in spinal ABMC in late puberty were 2-fold greater in girls who were homozygous for the C and G alleles, respectively (P = 0.001). For rs7757956, the difference in total body less head ABMC in late puberty was 50% less in individuals homozygous or heterozygous for the A allele (P = 0.006). Conclusions: Gains in ABMC in late pubertal girls are strongly associated with ESR1 polymorphisms, suggesting that estrogen contributes to this process via an estrogen receptor α-dependent pathway.

Analysis of body composition in individuals with high bone mass reveals a marked increase in fat mass in women but not men

Context: High bone mass (HBM), detected in 0.2% of dual-energy x-ray absorptiometry (DXA) scans, is characterized by raised body mass index, the basis for which is unclear. Objective: To investigate why body mass index is elevated in individuals with HBM, we characterized body composition and examined whether differences could be explained by bone phenotypes, eg, bone mass and/or bone turnover. Design, Setting, and Participants: We conducted a case-control study of 153 cases with unexplained HBM recruited from 4 UK centers by screening 219 088 DXA scans. Atotal of 138 first-degree relatives (of whom 51 had HBM) and 39 spouses were also recruited. Unaffected individuals served as controls. Main Outcome Measures: We measured fat mass, by DXA, and bone turnover markers. Results: Amongwomen, fat mass was inversely related to age in controls (P<.01), but not in HBM cases (P<.96) in whom mean fat mass was 8.9 [95% CI 4.7, 13.0] kg higher compared with controls (fully adjusted mean difference, P<.001). Increased fat mass in male HBM cases was less marked (gender interaction P = .03). Compared with controls, lean mass was also increased in female HBM cases (by 3.3 [1.2, 5.4] kg; P<.002); however, lean mass increases wereless marked than fat mass increases, resulting in 4.5% lower percentage lean mass in HBM cases (P<.001). Osteocalcin was also lower in female HBM cases compared with controls (by 2.8 [0.1, 5.5]μg/L; P = .04). Differences in fat mass were fully attenuated after hip bone mineral density (BMD) adjustment (P = .52) but unchanged after adjustment for bone turnover (P < .001), whereas the greater hip BMD in female HBM cases was minimally attenuated by fat mass adjustment (P<.001). Conclusions: HBM is characterized by a marked increase in fat mass in females, statistically explained by their greater BMD, but not by markers of bone turnover. Copyright © 2013 by The Endocrine Society.

Elevated circulating sclerostin concentrations in individuals with high bone mass, with and without LRP5 mutations

CONTEXT: The role and importance of circulating sclerostin is poorly understood. High bone mass (HBM) caused by activating LRP5 mutations has been reported to be associated with increased plasma sclerostin concentrations; whether the same applies to HBM due to other causes is unknown. OBJECTIVE: Our objective was to determine circulating sclerostin concentrations in HBM. DESIGN AND PARTICIPANTS: In this case-control study, 406 HBM index cases were identified by screening dual-energy x-ray absorptiometry (DXA) databases from 4 United Kingdom centers (n = 219 088), excluding significant osteoarthritis/artifact. Controls comprised unaffected relatives and spouses. MAIN MEASURES: Plasma sclerostin; lumbar spine L1, total hip, and total body DXA; and radial and tibial peripheral quantitative computed tomography (subgroup only) were evaluated. RESULTS: Sclerostin concentrations were significantly higher in both LRP5 HBM and non-LRP5 HBM cases compared with controls: mean (SD) 130.1 (61.7) and 88.0 (39.3) vs 66.4 (32.3) pmol/L (both P < .001, which persisted after adjustment for a priori confounders). In combined adjusted analyses of cases and controls, sclerostin concentrations were positively related to all bone parameters found to be increased in HBM cases (ie, L1, total hip, and total body DXA bone mineral density and radial/tibial cortical area, cortical bone mineral density, and trabecular density). Although these relationships were broadly equivalent in HBM cases and controls, there was some evidence that associations between sclerostin and trabecular phenotypes were stronger in HBM cases, particularly for radial trabecular density (interaction P < .01). CONCLUSIONS: Circulating plasma sclerostin concentrations are increased in both LRP5 and non-LRP5 HBM compared with controls. In addition to the general positive relationship between sclerostin and DXA/peripheral quantitative computed tomography parameters, genetic factors predisposing to HBM may contribute to increased sclerostin levels.

An experimental and finite element investigation of the biomechanics of vertebral compression fractures

Osteoporosis is a disease characterized by low bone mass and micro-architectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fracture. Osteoporosis affects over 200 million people worldwide, with an estimated 1.5 million fractures annually in the United States alone, and with attendant costs exceeding $10 billion dollars per annum. Osteoporosis reduces bone density through a series of structural changes to the honeycomb-like trabecular bone structure (micro-structure). The reduced bone density, coupled with the microstructural changes, results in significant loss of bone strength and increased fracture risk. Vertebral compression fractures are the most common type of osteoporotic fracture and are associated with pain, increased thoracic curvature, reduced mobility, and difficulty with self care. Surgical interventions, such as kyphoplasty or vertebroplasty, are used to treat osteoporotic vertebral fractures by restoring vertebral stability and alleviating pain. These minimally invasive procedures involve injecting bone cement into the fractured vertebrae. The techniques are still relatively new and while initial results are promising, with the procedures relieving pain in 70-95% of cases, medium-term investigations are now indicating an increased risk of adjacent level fracture following the procedure. With the aging population, understanding and treatment of osteoporosis is an increasingly important public health issue in developed Western countries. The aim of this study was to investigate the biomechanics of spinal osteoporosis and osteoporotic vertebral compression fractures by developing multi-scale computational, Finite Element (FE) models of both healthy and osteoporotic vertebral bodies. The multi-scale approach included the overall vertebral body anatomy, as well as a detailed representation of the internal trabecular microstructure. This novel, multi-scale approach overcame limitations of previous investigations by allowing simultaneous investigation of the mechanics of the trabecular micro-structure as well as overall vertebral body mechanics. The models were used to simulate the progression of osteoporosis, the effect of different loading conditions on vertebral strength and stiffness, and the effects of vertebroplasty on vertebral and trabecular mechanics. The model development process began with the development of an individual trabecular strut model using 3D beam elements, which was used as the building block for lattice-type, structural trabecular bone models, which were in turn incorporated into the vertebral body models. At each stage of model development, model predictions were compared to analytical solutions and in-vitro data from existing literature. The incremental process provided confidence in the predictions of each model before incorporation into the overall vertebral body model. The trabecular bone model, vertebral body model and vertebroplasty models were validated against in-vitro data from a series of compression tests performed using human cadaveric vertebral bodies. Firstly, trabecular bone samples were acquired and morphological parameters for each sample were measured using high resolution micro-computed tomography (CT). Apparent mechanical properties for each sample were then determined using uni-axial compression tests. Bone tissue properties were inversely determined using voxel-based FE models based on the micro-CT data. Specimen specific trabecular bone models were developed and the predicted apparent stiffness and strength were compared to the experimentally measured apparent stiffness and strength of the corresponding specimen. Following the trabecular specimen tests, a series of 12 whole cadaveric vertebrae were then divided into treated and non-treated groups and vertebroplasty performed on the specimens of the treated group. The vertebrae in both groups underwent clinical-CT scanning and destructive uniaxial compression testing. Specimen specific FE vertebral body models were developed and the predicted mechanical response compared to the experimentally measured responses. The validation process demonstrated that the multi-scale FE models comprising a lattice network of beam elements were able to accurately capture the failure mechanics of trabecular bone; and a trabecular core represented with beam elements enclosed in a layer of shell elements to represent the cortical shell was able to adequately represent the failure mechanics of intact vertebral bodies with varying degrees of osteoporosis. Following model development and validation, the models were used to investigate the effects of progressive osteoporosis on vertebral body mechanics and trabecular bone mechanics. These simulations showed that overall failure of the osteoporotic vertebral body is initiated by failure of the trabecular core, and the failure mechanism of the trabeculae varies with the progression of osteoporosis; from tissue yield in healthy trabecular bone, to failure due to instability (buckling) in osteoporotic bone with its thinner trabecular struts. The mechanical response of the vertebral body under load is highly dependent on the ability of the endplates to deform to transmit the load to the underlying trabecular bone. The ability of the endplate to evenly transfer the load through the core diminishes with osteoporosis. Investigation into the effect of different loading conditions on the vertebral body found that, because the trabecular bone structural changes which occur in osteoporosis result in a structure that is highly aligned with the loading direction, the vertebral body is consequently less able to withstand non-uniform loading states such as occurs in forward flexion. Changes in vertebral body loading due to disc degeneration were simulated, but proved to have little effect on osteoporotic vertebra mechanics. Conversely, differences in vertebral body loading between simulated invivo (uniform endplate pressure) and in-vitro conditions (where the vertebral endplates are rigidly cemented) had a dramatic effect on the predicted vertebral mechanics. This investigation suggested that in-vitro loading using bone cement potting of both endplates has major limitations in its ability to represent vertebral body mechanics in-vivo. And lastly, FE investigation into the biomechanical effect of vertebroplasty was performed. The results of this investigation demonstrated that the effect of vertebroplasty on overall vertebra mechanics is strongly governed by the cement distribution achieved within the trabecular core. In agreement with a recent study, the models predicted that vertebroplasty cement distributions which do not form one continuous mass which contacts both endplates have little effect on vertebral body stiffness or strength. In summary, this work presents the development of a novel, multi-scale Finite Element model of the osteoporotic vertebral body, which provides a powerful new tool for investigating the mechanics of osteoporotic vertebral compression fractures at the trabecular bone micro-structural level, and at the vertebral body level.

Mechanical considerations in fracture fixation

Bone’s capacity to repair following trauma is both unique and astounding. However, fractures sometimes fail to heal. Hence, the goal of fracture treatment is the restoration of bone’s structure, composition and function. Fracture fixation devices should provide a favourable mechanical and biological environment for healing to occur. The use of internal fixation is increasing as these devices may be applied with less invasive techniques. Recent studies suggest however that, internal fixation devices may be overly stiff and suppresses callus formation. The degree of mechanical stability influences the healing outcome. This is determined by the stiffness of the fixation device and the degree of limb loading. This project aims to characterise the fixation stability of an internal plate fixation device and the influence of modifications to its configuration on implant stability. As there are no standardised methods for the determination of fixation stiffness, the first part of this project aims to compares different methodologies and determines the most appropriate method to characterise the stiffness of internal plate fixators. The stiffness of a fixation device also influences the physiological loads experienced by the healing bone. Since bone adapts to this applied load by undergoing changes through a remodelling process, undesirable changes could occur during the period of treatment with an implant. The second part of this project aims to develop a methodology to quantify remodelling changes. This quantification is expected to aid our understanding of the changes in pattern due to implant related remodelling and on the factors driving the remodelling process. Knowledge gained in this project is useful to understand how the configuration of internal fixation devices can promote timely healing and prevent undesirable bone loss.

Osteopenia in a teenage boy presenting with previously undiagnosed guanidinoacetate methyltransferase (GAMT) deficiency and response to creatine supplementation

A 16 y.o. fully ambulant boy born to consanguineous Indian parents, presented for assessment of a fragility femoral neck fracture sustained against a background of autism and moderately severe intellectual disability. He had a past history of infantile eczema, and epilepsy treated with anticonvulsants from 2 to 10 years of age, with no further seizures following cessation of anticonvulsants. He had a thin body habitus (see Table 1) with long fingers and a high arched palate. He had no speech and negligible social interaction, but physical examination was otherwise unremarkable. Positive investigations revealed an undetectable serum creatinine and a urinary metabolic screen which showed an elevated GUA:Phe of 160 (< 36) and a decreased creatinine of 0.3 mmol/l (1.2–29.5) consistent with the diagnosis of guanidinoacetate methyltransferase(GAMT) deficiency. He was commenced on oral creatine 5 g three times daily. Despite improvement in physical activity, height and bone density, there was no discernable improvement in his intellectual functioning. Proton and phosphorous brain and leg magnetic resonance spectroscopy(MRS) was performed at baseline and showed an increased inorganic phosphorus peak and decreased phosphocreatine synthesis in brain and decreased creatine concentration in muscle. Following creatine treatment total brain creatine(1H-MRS) and phosphocreatine/ATP ratio (31P-MRS) content increased to 30% and 60% of control values, respectively. Brain GUA returned to normal levels.

Five years of GWAS discovery

The past five years have seen many scientific and biological discoveries made through the experimental design of genome-wide association studies (GWASs). These studies were aimed at detecting variants at genomic loci that are associated with complex traits in the population and, in particular, at detecting associations between common single-nucleotide polymorphisms (SNPs) and common diseases such as heart disease, diabetes, auto-immune diseases, and psychiatric disorders. We start by giving a number of quotes from scientists and journalists about perceived problems with GWASs. We will then briefly give the history of GWASs and focus on the discoveries made through this experimental design, what those discoveries tell us and do not tell us about the genetics and biology of complex traits, and what immediate utility has come out of these studies. Rather than giving an exhaustive review of all reported findings for all diseases and other complex traits, we focus on the results for auto-immune diseases and metabolic diseases. We return to the perceived failure or disappointment about GWASs in the concluding section. © 2012 The American Society of Human Genetics.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Genetic studies of osteoporosis: A rethink required

It has been 10 years since the seminal paper by Morrison and colleagues reporting the association of alleles of the vitamin D receptor and bone density [1], a paper which arguably kick-started the study of osteoporosis genetics. Since that report there have been literally thousands of osteoporosis genetic studies published, and large numbers of genes have been reported to be associated with the condition [2]. Although some of these reported associations are undoubtedly true, this snow-storm of papers and abstracts has clouded the field to such a great extent that it is very difficult to be certain of the veracity of most genetic associations reported hereto. The field needs to take stock and reconsider the best way forward, taking into account the biology of skeletal development and technological and statistical advances in human genetics, before more effort and money is wasted on continuing a process in which the primary achievement could be said to be a massive paper mountain. I propose in this review that the primary reasons for the paucity of success in osteoporosis genetics has been: •the absence of a major gene effect on bone mineral density (BMD), the most commonly studied bone phenotype; •failure to consider issues such as genetic heterogeneity, gene–environment interaction, and gene–gene interaction; •small sample sizes and over-optimistic data interpretation; and •incomplete assessment of the genetic variation in candidate genes studied.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.