Search

Later puberty associates with lower areal bone mineral density (aBMD), and both are risk factors for osteoporosis. However, the association between puberty timing-associated genetic variants and aBMD during development, and the causal relationship between puberty timing and aBMD, remain uncharacterized. We constructed sex-specific polygenic risk scores (GRS) consisting of 333 genetic variants associated with later puberty in European-descent children in the Bone Mineral Density in Childhood Study (BMDCS), consisting of a longitudinal cohort with up to seven assessments (n = 933) and a cross-sectional cohort (n = 486). These GRS were tested for associations with age- and sex-specific aBMD Z-scores at the lumbar spine (LS), femoral neck (FN), total hip, and distal radius, accounting for clinical covariates using sex-stratified linear mixed models. The causal relationship between puberty timing and aBMD was tested in the BMDCS and in publicly available adult data (GEFOS consortium) using two-sample Mendelian randomization (MR). The puberty-delaying GRS was associated with later puberty and lower LS-aBMD in the BMDCS in both sexes (combined beta ± SE = -0.078 ± 0.024; p = 0.0010). In the MR framework, the puberty-delaying genetic instrument also supported a causal association with lower LS-aBMD and FN-aBMD in adults of both sexes. Our results suggest that pubertal timing is causal for diminished aBMD in a skeletal site- and sex-specific manner that tracks throughout life, potentially impacting later risk for osteoporosis, which should be tested in future studies. © 2017 American Society for Bone and Mineral Research., (© 2017 American Society for Bone and Mineral Research.)

Failure to achieve optimal bone mineral accretion during childhood and adolescence results in subsequent suboptimal peak bone mass, contributing to osteoporosis risk later in life. To identify novel genetic factors that influence pediatric bone mass at discrete skeletal sites, we performed a sex-stratified genomewide association study of areal bone mineral density (BMD) measured by dual-energy X-ray absorptiometry at the 1/3 distal radius, spine, total hip, and femoral neck in a cohort of 933 healthy European American children. We took forward signals with p < 5 × 10 -5 and minor allele frequency (MAF) >5% into an independent cohort of 486 European American children in search of replication. In doing so, we identified five loci that achieved genome wide significance in the combined cohorts (nearest genes: CPED1, IZUMO3, RBFOX1, SPBT, and TBPL2), of which the last four were novel and two were sex-specific (SPTB in females and IZUMO3 in males), with all of them yielding associations that were particularly strong at a specific skeletal site. Annotation of potential regulatory function, expression quantitative trait loci (eQTL) effects and pathway analyses identified several potential target genes at these associated loci. This study highlights the importance of sex-stratified analyses at discrete skeletal sites during the critical period of bone accrual, and identifies novel loci for further functional follow-up to pinpoint key genes and better understand the regulation of bone development in children. © 2017 American Society for Bone and Mineral Research., (© 2017 American Society for Bone and Mineral Research.)

A recent whole-genome sequencing study in search of variation associated with adult areal bone mineral density (aBMD) identified rare variants near EN1, with markedly large effect sizes, and a common variant near SOX6. To understand the developmental effects of these loci, we sought to determine if they were associated with pediatric dual-energy X-ray absorptiometry-derived aBMD and bone mineral content (BMC) and if the associations were modified by sex. Our sample comprised 733 females and 685 males of European ancestry enrolled in the longitudinal Bone Mineral Density in Childhood Study (up to 7 annual study visits). Sex- and age-specific Z-scores, adjusted for height, were calculated for the total hip, femoral neck, spine, and distal radius. Total body less head (TBLH) BMC Z-scores were also calculated. The previously reported single nucleotide polymorphisms (SNPs) near EN1 and SOX6 were derived from our imputed data set. Linear mixed-effects models were used to test associations between each SNP and bone Z-scores, plus interactions with sex were explored. The rare T allele of lead EN1 SNP rs11692564 was associated with higher aBMD Z-score for total hip (beta = 0.62, p = 9.0 × 10(-4) ) and femoral neck (beta = 0.53, p = 0.010). In sex-stratified analyses, this variant was associated with higher bone Z-scores in females only, with the associations being strongest for total hip (sex interaction p = 1.9 × 10(-4) ; beta females = 0.86, p = 6.6 × 10(-6) ) and femoral neck (sex interaction p = 0.016; beta females = 0.73, p = 0.001). The common G allele of SOX6 SNP rs11024028 was associated with higher aBMD Z-score for total hip (beta = 0.12, p = 0.009), femoral neck (beta = 0.13, p = 0.003), and TBLH-BMC (beta = 0.09, p = 0.007); furthermore, this association strengthened in males in the sex-stratified analyses. Our findings reveal that rare genetic variation near EN1 and common variation near SOX6 operates in childhood and has implications for the lifelong risk of osteoporosis and fracture. The sex differences observed need to be independently replicated. © 2016 American Society for Bone and Mineral Research., (© 2016 American Society for Bone and Mineral Research.)

Both genetics and physical activity (PA) contribute to bone mineral density (BMD), but it is unknown if the benefits of physical activity on childhood bone accretion depend on genetic risk. We, therefore, aimed to determine if PA influenced the effect of bone fragility genetic variants on BMD in childhood. Our sample comprised US children of European ancestry enrolled in the Bone Mineral Density in Childhood Study (N = 918, aged 5 to 19 years, and 52.4% female). We used a questionnaire to estimate hours per day spent in total, high-, and low-impact PA. We calculated a BMD genetic score (% BMD lowering alleles) using adult genome-wide association study (GWAS)-implicated BMD variants. We used dual-energy X-ray absorptiometry to estimate femoral neck, total hip, and spine areal-BMD and total body less head (TBLH) bone mineral content (BMC) Z-scores. The BMD genetic score was negatively associated with each bone Z-score (eg, TBLH-BMC: estimate = -0.03, p = 1.3 × 10(-6) ). Total PA was positively associated with bone Z-scores; these associations were driven by time spent in high-impact PA (eg, TBLH-BMC: estimate = 0.05, p = 4.0 × 10(-10) ) and were observed even for children with lower than average bone Z-scores. We found no evidence of PA-adult genetic score interactions (p interaction > 0.05) at any skeletal site, and there was no evidence of PA-genetic score-Tanner stage interactions at any skeletal site (p interaction > 0.05). However, exploratory analyses at the individual variant level revealed that PA statistically interacted with rs2887571 (ERC1/WNT5B) to influence TBLH-BMC in males (p interaction = 7.1 × 10(-5) ), where PA was associated with higher TBLH-BMC Z-score among the BMD-lowering allele carriers (rs2887571 AA homozygotes: estimate = 0.08 [95% CI 0.06, 0.11], p = 2.7 × 10(-9) ). In conclusion, the beneficial effect of PA on bone, especially high-impact PA, applies to the average child and those genetically predisposed to lower adult BMD (based on GWAS-implicated BMD variants). Independent replication of our exploratory individual variant findings is warranted. © 2016 American Society for Bone and Mineral Research., (© 2016 American Society for Bone and Mineral Research.)

Using adult identified bone mineral density (BMD) loci, we calculated genetic risk scores (GRS) to determine if they were associated with changes in BMD during childhood. Longitudinal data from the Bone Mineral Density in Childhood Study were analyzed (N = 798, 54% female, all European ancestry). Participants had up to 6 annual dual energy X-ray scans, from which areal BMD (aBMD) Z-scores for the spine, total hip, and femoral neck were estimated, as well as total body less head bone mineral content (TBLH-BMC) Z-scores. Sixty-three single-nucleotide polymorphisms (SNPs) were genotyped, and the percentage of BMD-lowering alleles carried was calculated (overall adult GRS). Subtype GRS that include SNPs associated with fracture risk, pediatric BMD, WNT signaling, RANK-RANKL-OPG, and mesenchymal stem cell differentiation were also calculated. Linear mixed effects models were used to test associations between each GRS and bone Z-scores, and if any association differed by sex and/or chronological age. The overall adult, fracture, and WNT signaling GRS were associated with lower Z-scores (eg, spine aBMD Z-score: βadult  = -0.04, p = 3.4 × 10(-7) ; βfracture = -0.02, p = 8.9 × 10(-6) ; βWNT  = -0.01, p = 3.9 × 10(-4) ). The overall adult GRS was more strongly associated with lower Z-scores in females (p-interaction ≤ 0.05 for all sites). The fracture GRS was more strongly associated with lower Z-scores with increasing age (p-interaction ≤ 0.05 for all sites). The WNT GRS associations remained consistent for both sexes and all ages (p-interaction > 0.05 for all sites). The RANK-RANKL-OPG GRS was more strongly associated in females with increasing age (p-interaction < 0.05 for all sites). The mesenchymal stem cell GRS was associated with lower total hip and femoral neck Z-scores, in both boys and girls, across all ages. No associations were observed between the pediatric GRS and bone Z-scores. In conclusion, adult identified BMD loci associated with BMD and BMC in the pediatric setting, especially in females and in loci involved in fracture risk and WNT signaling., (© 2015 American Society for Bone and Mineral Research.)