Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

The impact of intense training on endogenous estrogen and progesterone concentrations and bone mineral acquisition in adolescent rowers

The effect of 18 months of training on the ovarian hormone concentrations and bone mineral density (BMD) accrual was assessed longitudinally in 14 adolescent rowers and 10 matched controls, aged 14–15 years. Ovarian hormone levels were assessed by urinary estrone glucuronide (E1G) and pregnanediol glucuronide (PdG) excretion rates, classifying the menstrual cycles as ovulatory or anovulatory. Total body (TB), total proximal femur (PF), femoral neck (FN) and lumbar spine (LS) (L2–4) bone mass were measured at baseline and 18 months using dual-energy X-ray densitometry. Results were expressed as bone mineral content (BMC), BMD and bone mineral apparent density (BMAD). Five rowers had anovulatory menstrual cycles compared with zero prevalence for the control subjects. Baseline TB BMD was significantly higher in the ovulatory rowers, with PF BMD, FN BMD and LS BMD similar for all groups. At completion, the LS bone accrual of the ovulatory rowers was significantly greater (BMC 8.1%, BMD 6.2%, BMAD 6.2%) than that of the anovulatory rowers (BMC 1.1%, BMD 3.9%, BMAD 1.6%) and ovulatory controls (BMC 0.5%, BMD 1.1%, BMAD 1.1%). No difference in TB, PF or FN bone accrual was observed among groups. This study demonstrated an osteogenic response to mechanical loading, with the rowers accruing greater bone mass than the controls at the lumbar spine. However, the exercise-induced osteogenic benefits were less when rowing training was associated with low estrogen and progesterone metabolite excretion.

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.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Subchondral bone osteoblasts can induce chondrocyte mineralization during osteoarthritis and this process is relevant to cartilage degradation

Pathological mineralization of articular cartilage is a characteristic feature of osteoarthritis (OA); however, the underlying mechanisms, and their relevance to cartilage degeneration, are not clear. The involvement of subchondral bone changes in OA have been reported previously with the characterization of abnormal subchondral bone mineral density (BMD), osteiod volume, altered bone mechanical parameters and an increase in bone turnover markers. A number of osteoarthritic animal models have demonstrated that subchondral bone changes often precede cartilage degeneration. In this study site specific localization of mineralization markers were detected in the OA cartilage. Chondrocytes and osteoblasts derived from OA cartilage and subchondral bone showed a significant increase in the mRNA expressions of mineralization markers. Interestingly, osteoblasts from OA subchondral bone could significantly decrease cartilage matrix expression; whereas, increase mineralization of chondrocytes (Figure 1). Osteogenic factors, such as CBFA1, ALP, and type X collagen (Col-X), were detected in chondrocytes under mineralization conditions (Figure 2). Furthermore, chondrocyte mineralization was followed by increased mRNA and protein levels of MMP-2, MMP-9 and MMP-13, all of which are detrimental to cartilage integrity in vivo. The data reported here suggests that the upregulation of subchondral bone-mineralization, typical of OA progression, causes cartilage mineralization, and that the mineralization of chondrocytes induce increased MMP levels with a subsequent degradation of the articular cartilage.

Near infrared (NIR) absorption spectra correlates with subchondral bone micro-CT parameters in osteoarthritic rat models

Determining the properties and integrity of subchondral bone in the developmental stages of osteoarthritis, especially in a form that can facilitate real-time characterization for diagnostic and decision-making purposes, is still a matter for research and development. This paper presents relationships between near infrared absorption spectra and properties of subchondral bone obtained from 3 models of osteoarthritic degeneration induced in laboratory rats via: (i) menisectomy (MSX); (ii) anterior cruciate ligament transaction (ACL); and (iii) intra-articular injection of mono-ido-acetate (1 mg) (MIA), in the right knee joint, with 12 rats per model group (N = 36). After 8 weeks, the animals were sacrificed and knee joints were collected. A custom-made diffuse reflectance NIR probe of diameter 5 mm was placed on the tibial surface and spectral data were acquired from each specimen in the wavenumber range 4000–12 500 cm− 1. After spectral acquisition, micro computed tomography (micro-CT) was performed on the samples and subchondral bone parameters namely: bone volume (BV) and bone mineral density (BMD) were extracted from the micro-CT data. Statistical correlation was then conducted between these parameters and regions of the near infrared spectra using multivariate techniques including principal component analysis (PCA), discriminant analysis (DA), and partial least squares (PLS) regression. Statistically significant linear correlations were found between the near infrared absorption spectra and subchondral bone BMD (R2 = 98.84%) and BV (R2 = 97.87%). In conclusion, near infrared spectroscopic probing can be used to detect, qualify and quantify changes in the composition of the subchondral bone, and could potentially assist in distinguishing healthy from OA bone as demonstrated with our laboratory rat models.

Mice lacking β-Adrenergic receptors have increased bone mass but are not protected from deleterious skeletal effects of ovariectomy

Activation of β2-adrenergic receptors inhibits osteoblastic bone formation and enhances osteoclastic bone resorption. Whether β-blockers inhibit ovariectomy-induced bone loss and decrease fracture risk remains controversial. To further explore the role of β-adrenergic signaling in skeletal acquisition and response to estrogen deficiency, we evaluated mice lacking the three known β-adrenergic receptors (β-less). Body weight, percent fat, and bone mineral density were significantly higher in male β-less than wild-type (WT) mice, more so with increasing age. Consistent with their greater fat mass, serum leptin was significantly higher in β-less than WT mice. Mid-femoral cross-sectional area and cortical thickness were significantly higher in adult β-less than WT mice, as were femoral biomechanical properties (+28 to +49%, P < 0.01). Young male β-less had higher vertebral (1.3-fold) and distal femoral (3.5-fold) trabecular bone volume than WT (P < 0.001 for both) and lower osteoclast surface. With aging, these differences lessened, with histological evidence of increased osteoclast surface and decreased bone formation rate at the distal femur in β-less vs. WT mice. Serum tartrate-resistance alkaline phosphatase-5B was elevated in β-less compared with WT mice from 8–16 wk of age (P < 0.01). Ovariectomy inhibited bone mass gain and decreased trabecular bone volume/total volume similarly in β-less and WT mice. Altogether, these data indicate that absence of β-adrenergic signaling results in obesity and increased cortical bone mass in males but does not prevent deleterious effects of estrogen deficiency on trabecular bone microarchitecture. Our findings also suggest direct positive effects of weight and/or leptin on bone turnover and cortical bone structure, independent of adrenergic signaling.

COL1A1 C-propeptide cleavage site mutation causes high bone mass, bone fragility and jaw lesions: A new cause of gnathodiaphyseal dysplasia?

Gnathodiaphyseal dysplasia (GDD) is a rare autosomal dominant condition characterized by bone fragility, irregular bone mineral density (BMD) and fibro-osseous lesions in the skull and jaw. Mutations in Anoctamin-5 (ANO5) have been identified in some cases. We aimed to identify the causative mutation in a family with features of GDD but no mutation in ANO5, using whole exome capture and massive parallel sequencing (WES). WES of two affected individuals (a mother and son) and the mother's unaffected parents identified a mutation in the C-propeptide cleavage site of COL1A1. Similar mutations have been reported in individuals with osteogenesis imperfecta (OI) and paradoxically increased BMD. C-propeptide cleavage site mutations in COL1A1 may not only cause 'high bone mass OI', but also the clinical features of GDD, specifically irregular sclerotic BMD and fibro-osseous lesions in the skull and jaw. GDD patients negative for ANO5 mutations should be assessed for mutations in type I collagen C-propeptide cleavage sites.

Life-course determinants of bone mass in young adults from a transitional rural community in India: The Andhra Pradesh Children and Parents Study (APCAPS)

Background: Undernutrition and physical inactivity are both associated with lower bone mass. Objective: This study aimed to investigate the combined effects of early-life undernutrition and urbanized lifestyles in later life on bone mass accrual in young adults from a rural community in India that is undergoing rapid socioeconomic development. Design: This was a prospective cohort study of participants of the Hyderabad Nutrition Trial (1987–1990), which offered balanced protein-calorie supplementation to pregnant women and preschool children younger than 6 y in the intervention villages. The 2009–2010 follow-up study collected data on current anthropometric measures, bone mineral density (BMD) measured by dual-energy X-ray absorptiometry, blood samples, diet, physical activity, and living standards of the trial participants (n = 1446, aged 18–23 y). Results: Participants were generally lean and had low BMD [mean hip BMD: 0.83 (women), 0.95 (men) g/cm2; lumbar spine: 0.86 (women), 0.93 (men) g/cm2]. In models adjusted for current risk factors, no strong evidence of a positive association was found between BMD and early-life supplementation. On the other hand, current lean mass and weight-bearing physical activity were positively associated with BMD. No strong evidence of an association was found between BMD and current serum 25-hydroxyvitamin D or dietary intake of calcium, protein, or calories. Conclusions: Current lean mass and weight-bearing physical activity were more important determinants of bone mass than was early-life undernutrition in this population. In transitional rural communities from low-income countries, promotion of physical activity may help to mitigate any potential adverse effects of early nutritional disadvantage.

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.

Solid volume fraction estimation of bone:marrow replica models using ultrasound transit time spectroscopy

The acceptance of broadband ultrasound attenuation (BUA) for the assessment of osteoporosis suffers from a limited understanding of both ultrasound wave propagation through cancellous bone and its exact dependence upon the material and structural properties. It has recently been proposed that ultrasound wave propagation in cancellous bone may be described by a concept of parallel sonic rays; the transit time of each ray defined by the proportion of bone and marrow propagated. A Transit Time Spectrum (TTS) describes the proportion of sonic rays having a particular transit time, effectively describing the lateral inhomogeneity of transit times over the surface aperture of the receive ultrasound transducer. The aim of this study was to test the hypothesis that the solid volume fraction (SVF) of simplified bone:marrow replica models may be reliably estimated from the corresponding ultrasound transit time spectrum. Transit time spectra were derived via digital deconvolution of the experimentally measured input and output ultrasonic signals, and compared to predicted TTS based on the parallel sonic ray concept, demonstrating agreement in both position and amplitude of spectral peaks. Solid volume fraction was calculated from the TTS; agreement between true (geometric calculation) with predicted (computer simulation) and experimentally-derived values were R2=99.9% and R2=97.3% respectively. It is therefore envisaged that ultrasound transit time spectroscopy (UTTS) offers the potential to reliably estimate bone mineral density and hence the established T-score parameter for clinical osteoporosis assessment.