Abstract

Purpose: Modic changes (MC) on magnetic resonance imaging (MRI) have been associated with the development and severity of low back pain (LBP). The etiology of MC remains elusive, but it has been suggested that altered metabolism may be a risk factor. As such, this study aimed to identify metabolomic biomarkers for MC phenotypes of the lumbar spine via a combined metabolomic-genomic approach.

Methods: A population cohort of 3,584 southern Chinese underwent lumbar spine MRI. Blood samples were genotyped with single-nucleotide polymorphisms (SNP) arrays (n = 2,482) and serum metabolomics profiling using magnetic resonance spectroscopy (n  = 757), covering 130 metabolites representing three molecular windows, were assessed. Genome-wide association studies (GWAS) were performed on each metabolite, to construct polygenic scores for predicting metabolite levels in subjects who had GWAS but not metabolomic data. Associations between predicted metabolite levels and MC phenotypes were assessed using linear/logistic regression and least absolute shrinkage and selection operator (LASSO). Two-sample Mendelian randomization analysis tested for causal relationships between metabolic biomarkers and MC.

Results: 20.4% had MC (10.6% type 1, 67.2% type 2, 22.2% mixed types). Significant MC metabolomic biomarkers were mean diameter of very-low-density lipoprotein (VLDL)/low-density lipoprotein (LDL) particles and cholesterol esters/phospholipids in large LDL. Mendelian randomization indicated that decreased VLDL mean diameter may lead to MC.

Conclusions: This large-scale study is the first to address metabolomics in subject with/without lumbar MC. Causality studies implicate VLDL related to MC, noting a metabolic etiology. Our study substantiates the field of “spino-metabolomics” and illustrates the power of integrating metabolomics-genomics-imaging phenotypes to discover biomarkers for spinal disorders, paving the way for more personalized spine care for patients.