BMD loci contribute to ethnic and developmental differences in skeletal fragility across populations: Assessment of evolutionary selection pressures

Bone mineral density (BMD) is a highly heritable trait used both for the diagnosis of osteoporosis in adults and to assess bone health in children. Ethnic differences in BMD have been documented, with markedly higher levels in individuals of African descent, which partially explain disparity in osteoporosis risk across populations. To date, 63 independent genetic variants have been associated with BMD in adults of Northern-European ancestry. Here, we demonstrate that at least 61 of these variants are predictive of BMD early in life by studying their compound effect within two multiethnic pediatric cohorts. Furthermore, we show that within these cohorts and across populations worldwide the frequency of those alleles associated with increased BMD is systematically elevated in individuals of Sub-Saharan African ancestry. The amount of differentiation in the BMD genetic scores among Sub-Saharan and non-Sub-Saharan populations together with neutrality tests, suggest that these allelic differences are compatible with the hypothesis of selective pressures acting on the genetic determinants of BMD. These findings constitute an explorative contribution to the role of selection on ethnic BMD differences and likely a new example of polygenic adaptation acting on a human trait

Trabecular Bone Score Reference Values for Children and Adolescents According to Age, Sex, and Ancestry.

Trabecular bone score (TBS) is used for fracture prediction in adults, but its utility in children is limited by absence of appropriate reference values. We aimed to develop reference ranges for TBS by age, sex, and population ancestry for youth ages 5 to 20 years. We also investigated the association between height, body mass index (BMI), and TBS, agreement between TBS and lumbar spine areal bone mineral density (aBMD) and bone mineral apparent density (BMAD) Z-scores, tracking of TBS Z-scores over time, and precision of TBS measurements. We performed secondary analysis of spine dual-energy X-ray absorptiometry (DXA) scans from the Bone Mineral Density in Childhood Study (BMDCS), a mixed longitudinal cohort of healthy children (n = 2014) evaluated at five US centers. TBS was derived using a dedicated TBS algorithm accounting for tissue thickness rather than BMI. TBS increased only during ages corresponding to pubertal development with an earlier increase in females than males. There were no differences in TBS between African Americans and non-African Americans. We provide sex-specific TBS reference ranges and LMS values for calculation of TBS Z-scores by age and means and SD for calculation of Z-scores by pubertal stage. TBS Z-scores were positively associated with height Z-scores at some ages. TBS Z-scores explained only 27% and 17% of the variance of spine aBMD and BMAD Z-scores. Tracking of TBS Z-scores over 6 years was lower (r = 0.47) than for aBMD or BMAD Z-scores (r = 0.74 to 0.79), and precision error of TBS (2.87%) was greater than for aBMD (0.85%) and BMAD (1.22%). In sum, TBS Z-scores provide information distinct from spine aBMD and BMAD Z-scores. Our robust reference ranges for TBS in a well-characterized pediatric cohort and precision error estimates provide essential tools for clinical assessment using TBS and determination of its value in predicting bone fragility in childhood and adolescence. © 2022 American Society for Bone and Mineral Research (ASBMR)

CYP11B1 variants influence skeletal maturation via alternative splicing.

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10 <sup>-12</sup> ). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10 <sup>-40</sup> ), exon 4 inclusion (P = 4.29 × 10 <sup>-34</sup> ), and decreased exon 3 and 5 splicing (P = 7.85 × 10 <sup>-43</sup> ), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism