Dissecting the Genetic Basis of Low Back Pain Independent of BMI Through Genomic Structural Equation Modeling
Journal of Pain Research, 2026
Huo L., Tan L., Wang F., Sun J., Niu Y., Jiang Y., Wu M., Shi J., Hao Y., Wang J., Huang S., Chen Z.
| Disease area | Application area | Sample type | Products |
|---|---|---|---|
Neurology | Pathophysiology | Plasma | Olink Explore 3072/384 |
Abstract
Background: Low back pain (LBP) is a leading cause of disability worldwide. Although body mass index (BMI) is a well-established risk factor for LBP, a substantial proportion of patients with LBP do not present with abnormal BMI, suggesting the involvement of BMI-independent mechanisms. However, the genetic architecture underlying BMI-independent LBP remains poorly understood. This study aimed to identify and characterize genetic variants associated with LBP that are independent of BMI.
Methods: This study was a secondary analysis of publicly available genome-wide association study (GWAS) summary statistics. Genetic associations were analyzed using a Genomic Structural Equation Modeling (Genomic SEM) framework. BMI summary statistics were obtained from the Genetic Investigation of ANthropometric Traits (GIANT) consortium (~700,000 individuals of European ancestry), and LBP data were derived from the FinnGen cohort, including 60,099 cases and 440,249 controls of European ancestry, with LBP defined by the International Classification of Diseases, 10th Revision (ICD-10) code M54. Genome-wide association study by subtraction (GWAS-by-subtraction) was applied to identify BMI-independent LBP associations. Statistical fine-mapping, transcriptome-wide association studies (TWAS), proteome-wide association studies (PWAS), and colocalization analyses were subsequently performed to prioritize putative causal genes.
Results: Three independent genome-wide significant loci associated with BMI-independent LBP were identified: rs6916321 at B7NZA1 (P = 2.41 × 10− 10), rs2596501 near HLA-B (P = 6.56 × 10⁻9), and a novel locus rs4148946 at SEC24C (P = 1.74 × 10⁻8). Fine-mapping highlighted rs2596501 as a likely causal variant with a posterior inclusion probability of 0.9999. Multi-omic integration consistently prioritized CHST3 as a candidate gene (P = 4.49 × 10⁻5 in PWAS). Pathway enrichment analyses implicated neuronal signaling and immune-related pathways, with cell-type enrichment observed in differentiated neurons (P = 0.0063). These findings were derived from large population-based cohorts of European ancestry.
Conclusion: This study refines the genetic architecture of BMI-independent LBP and identifies novel loci with convergent multi-omic evidence implicating CHST3 in disease susceptibility. The results highlight biological mechanisms beyond adiposity that may contribute to LBP risk and provide a foundation for future functional and translational research.