Supplementary methods.

<div><p>Background</p><p>Osteoporosis is a major global health issue, weakening bones and increasing fracture risk. Dual-energy X-ray absorptiometry (DXA) is the standard for measuring bone mineral density (BMD) and diagnosing osteoporosis, but its costliness and complexity impede widespread screening adoption. Predictive modeling using genetic and clinical data offers a cost-effective alternative for assessing osteoporosis and fracture risk. This study aims to develop BMD prediction models using data from the UK Biobank (UKBB) and test their performance across different ethnic and geographical populations.</p><p>Methods and findings</p><p>We developed BMD prediction models for the femoral neck (FNK) and lumbar spine (SPN) using both genetic variants and clinical factors (such as sex, age, height, and weight), within 17,964 British white individuals from UKBB. Models based on regression with least absolute shrinkage and selection operator (LASSO), selected based on the coefficient of determination (R<sup>2</sup>) from a model selection subset of 5,973 individuals from British white population. These models were tested on 5 UKBB test sets and 12 independent cohorts of diverse ancestries, totaling over 15,000 individuals. Furthermore, we assessed the correlation of predicted BMDs with fragility fractures risk in 10 years in a case-control set of 287,183 European white participants without DXA-BMDs in the UKBB.</p><p>With single-nucleotide polymorphism (SNP) inclusion thresholds at 5×10<sup>−6</sup> and 5×10<sup>−7</sup>, the prediction models for FNK-BMD and SPN-BMD achieved the highest R<sup>2</sup> of 27.70% with a 95% confidence interval (CI) of [27.56%, 27.84%] and 48.28% (95% CI [48.23%, 48.34%]), respectively. Adding genetic factors improved predictions slightly, explaining an additional 2.3% variation for FNK-BMD and 3% for SPN-BMD over clinical factors alone. Survival analysis revealed that the predicted FNK-BMD and SPN-BMD were significantly associated with fragility fracture risk in the European white population (<i>P</i> < 0.001). The hazard ratios (HRs) of the predicted FNK-BMD and SPN-BMD were 0.83 (95% CI [0.79, 0.88], corresponding to a 1.44% difference in 10-year absolute risk) and 0.72 (95% CI [0.68, 0.76], corresponding to a 1.64% difference in 10-year absolute risk), respectively, indicating that for every increase of one standard deviation in BMD, the fracture risk will decrease by 17% and 28%, respectively. However, the model’s performance declined in other ethnic groups and independent cohorts. The limitations of this study include differences in clinical factors distribution and the use of only SNPs as genetic factors.</p><p>Conclusions</p><p>In this study, we observed that combining genetic and clinical factors improves BMD prediction compared to clinical factors alone. Adjusting inclusion thresholds for genetic variants (e.g., 5×10<sup>−6</sup> or 5×10<sup>−7</sup>) rather than solely considering genome-wide association study (GWAS)-significant variants can enhance the model’s explanatory power. The study highlights the need for training models on diverse populations to improve predictive performance across various ethnic and geographical groups.</p></div>