Back

The genetic architecture of trabecular bone score and association with fracture: a genome-wide association and meta-analysis

Lu, H.; Tibbitts, G.; LI, G. H.-Y.; Frysz, M.; Lary, C. W.; Trajanoska, K.; Das, N. R.; de Groot, L.; Cheung, C.-L.; van Schoor, N. M.; Uitterlinden, A. G.; van der Velde, N.; Tobias, J. H.; Rivadeneira, F.; Evans, D. M.; Kiel, D. P.; Kemp, J. P.; Medina-Gomez, C.

2025-10-03 genomics
10.1101/2025.10.01.678241 bioRxiv
Show abstract

BackgroundTrabecular bone score (TBS) is a texture-based measurement derived from DXA scans, which describes the distribution of mineral across the vertebral bodies. Identifying its genetic determinants is crucial for enhancing understanding of its biological basis and clarifying its relationship with fracture risk. MethodsWe conducted a large-scale genome-wide association study (GWAS) of TBS from 44,767 participants (95.08% European ancestry). Sex heterogeneity was assessed through sex-stratified GWAS. Post-GWAS analyses included functional annotation, gene-set enrichment analysis, and cis-expression quantitative trait locus (cis-eQTL) colocalization. Two-sample and multivariable Mendelian randomization (MR) were employed to investigate the causal effect of TBS on fractures. FindingsWe identified 33 independent variants associated with TBS, which explained 4.06% of TBS phenotypic variance. All the discovered signals map to previously identified BMD-associated loci. Additionally, we identified genetic variants in the RAB11FIP3 locus that exhibited a sex-heterogeneous effect, being associated with TBS exclusively in males. Functional annotation and colocalization identified functional genes related to TBS. The two-sample and multivariable MR analyses indicated that TBS potentially has an independent causal effect on fracture risk. InterpretationOur study unveiled 33 independent loci associated with TBS, all in BMD-associated loci. One locus was further identified as associated with TBS only in males. MR results suggested that genetically derived TBS may be causally associated with fracture risk at different sites, potentially beyond BMD. This study provides insights into the TBS genetic architecture and uncovers its potential clinical applications in fracture risk prediction. FundingAll funding information can be found in the Acknowledgements section. Research in contextO_ST_ABSEvidence before this studyC_ST_ABSOsteoporosis is a common disease prevalent in older adults, mainly diagnosed by low bone mineral density and disruption of bone architecture. While the genetic determinants of bone mass have been well studied, the genetic architecture of bone features beyond BMD remains largely unknown. Added value of this studyThis study presents a large genome-wide association meta-analysis of the trabecular bone score (TBS), a measurement of trabecular distribution of bone mineralization across vertebral bodies. Our study identified 33 independent variants associated with TBS and one association signal only in males. Functional annotation highlighted crucial genes related to TBS. The two-sample and multivariable Mendelian randomization (MR) analyses suggested that TBS may be causally linked with fracture risk, potentially independent of BMD. Implications of all the available evidenceOur study represents an expansion of the traditional bone density phenotype used in most previous skeletal genetic studies. Based on the findings that all of the significant loci for TBS have been previously identified in GWAS of BMD, we conclude that TBS-associated variants have pleiotropic effects affecting not only mineral density but also the distribution patterns of minerals in the cancellous bone. Continuing to perform genetic studies of new skeletal phenotypes beyond BMD could ultimately impact the management and prediction of osteoporosis patients.

Matching journals

The top 4 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.