Identification of the molecular mechanism and diagnostic biomarkers in the thoracic ossification of the ligamentum flavum using metabolomics and transcriptomics
BMC Molecular and Cell Biology volume 21, Article number: 37 (2020)
To establish a metabolite fingerprint of ossification of the thoracic ligamentum flavum (OTLF) patients using liquid chromatography-mass spectrometry (LC-MS) in combination with transcriptomic data and explore the potential molecular mechanism of pathogenesis.
The study cohort was composed of 25 patients with OTLF and 23 healthy volunteers as a control group. Thirty-seven metabolites were identified out by UPLC-MS including uric acid and hypoxanthine. Nine metabolites, including uric acid and hypoxanthine, were found with a Variable Importance in Projection (VIP) score over 1 (p < 0.05). Pathway enrichment indicated that purine metabolism pathways and the other four metabolism pathways were enriched. Transcriptomic data revealed that purine metabolism have a substantial change in gene expression of OTLF and that xanthine dehydrogenase (XDH) is the key regulatory factor. Receiver operating characteristic (ROC) analysis indicated that 17 metabolites, including uric acid, were found with an AUC value of over 0.7.
Uric acid might be the potential biomarker for OTLF and play an important role within the detailed pathway. XDH could affect purine metabolism by suppressing the expression of hypoxanthine and xanthine leading to low serum levels of uric acid in OTLF, which could be a focal point in developing new therapeutic methods for OTLF.
Ossification of the ligamentum flavum (OLF) primarily occurs in the thoracolumbar spine, especially at T9-T11, and almost all OLF cases were exclusively reported in East Asian countries [1, 2]. The overall occurance rate in China is 3.8% of the general population in 2010 and it is more common in females . The first case of OLF was reported by Polgar in 1920 and was recognized as a cause of myeloradiculopathy . OLF is a relatively rare disease and numerous studies have focused on its progression at both histopathological and cellular levels. OLF is often accompanied by other spinal degenerative diseases , which means the recognition of its clinical manifestation, diagnosis and progression is extremely difficult. Imaging studies such as magnetic resonance imaging (MRI) and computed tomography (CT) are the most commonly employed methods for radiological diagnosis . As for the treatment of OTLF, surgical treatments are the first choice for patients due to the progressive nature of the disease and poor outcomes from conservative management . The most common surgical strategy is posterior decompression as it is the most effective , however this still has its limitations in that patients with multilevel OTLF cannot have the same outcome as those with single- or dual-level lesions .
Metabolomics also knowns as metabolic profiling is used in the detection of smaller molecules with a molecular weight of less than 1800 Da . Metabolomics can be used to investigate disease or systemic health from bio-samples such as body fluid and tissue . It is an effective approach to detect both intra- and extra-cellular processes through analyzing the metabolites which are produced by various biochemical pathways . Nontargeted and targeted analyses are two main approaches in metabolomics detection. Targeted analysis is focused on a specific compound set such as amino acids , lipids  and energy metabolism . In contrast, untargeted approaches aim to find potential biomarkers for the diagnosis of diseases and therefore it needs to measure thousands of molecules at a time . In order to analyze the metabolites in samples, separation techniques such as chromatography (both liquid and gas) and capillary electrophoresis are needed and these techniques are often combined with mass spectrometry to confirm the molecular weight, or used in nuclear magnetic resonance (NMR) to confirm specific molecular structures.
The studies of the diagnosis and treatment methods of OTLF mostly focus on the radiological features. There have been some other studies exploring the pathogenesis of OTLF such as mechanical , degenerative  and genetic factors  but detailed clinical progression of the disease continues to be poorly understood. Although there have been several studies that have analyzed biomarkers for ossification of the spinal ligament (OSL) and found calcium-phosphate metabolism markers, bone turnover markers, sclerostin and so on, there are no definitive conclusions to date . Few studies focus on the metabolism dysfunction in OTLF, one previous research indicated that lectin metabolism might affect the progress of OTLF  and thus further research on metabolism disorder may be needed to develop earlier diagnosis and new therapeutic methods.
In our current study, UPLC-MS system based untargeted metabolomic approaches were applied to find potential biomarkers in the OTLF serum samples, and transcriptomic techniques were also ultilized to identify and study the regulator gene associated with the metabolic pathways found in the metabolomic analysis. Findings from these approaches may provide useful information for the diagnosis and treatment of OTLF.
Demographic and clinical features
The study cohort was composed of 25 patients with OTLF and 23 healthy volunteers as a control group. The mean age of OTLF patients was 53.88 ± 9.46 years and the mean age of the control group was 53.71 ± 9.77 years indicating no significant difference. Specific patient demographic and clinical features are illustrated in Table 1.
A total of 37 metabolites were found in the samples of the two groups. Twenty-five metabolites were detected under the positive mode mass spectrometry (Fig. 1a) and 12 metabolites were found under the negative mode (Fig. 1c). PCA and OPLS-DA analysis were used to identify the metabolite and its phenotype. As shown in the PCA score plots (Fig. 1b&d), the serum samples within each group were closely clustered into each other in both positive and negative ion mode, while samples from different groups were clearly separated. Nine metabolites (including uric acid, triacetin and hypoxanthine) were found to have a Variable Importance in Projection (VIP) score over 1 (p < 0.05) (Fig. 2a). In order to comprehend the relation between differential metabolites and the metabolic pathways, a metabolite sets enrichment analysis was carried out, which showed out of the five metabolic pathways that underewent enrichment analysis, the pyrimidine and purine metabolism pathways showed significant differences in OTLF (Fig. 2b).
Transcriptomic data GSE69787 was acquired from GEO platform and it was used to perform that the relationship between genes and the metabolites in process of ligament ossification. Three thousand fifty-six genes were found have notable differences between two groups (student t’test p < 0.05, FDR < 0.05). Heatmap illustrated the top ten up- and down-regulated genes in two groups (Fig. 3a) and these genes were enriched to the metabolism pathway by Gene Set Enrichment Analysis (GSEA). Density joyplot and dotplot were plotted according to the results of GSEA and it showed 30 pathways had significant impact on the process of ligament ossification (Fig. 3b) and 8 metabolic pathways were suppressed in patients while 7 were activated (Fig. 3c). Combined with metabolite results, purine metabolism was also found have a substantial change in gene expression. 346 differential expression genes were selected under the condition of p < 0.0001, FDR < 0.05 to make the intersection with genes on the purine metabolism pathway. As shown in the Venn diagram in Fig. 3d common genes were found which were PDE7B, PDE1A and XDH .
Potential biomarker and related pathway analysis
Receiver operating characteristic curves (ROC) were plotted to inspect the Area Under Curve (AUC) which is an effective and combined measure of sensitivity and specificity that describes the inherent validity of diagnostic tests. A total 17 metabolites were found with an AUC value (Table 2) of over 0.7 which included uric acid (Fig. 4) suggesting that the uric acid may be the potential biomarker for OTLF. Figure 5 illustrates that the relationship between XDH, uric acid and hypoxanthine within the purine metabolism pathway. The metabolism of hypoxanthine was regulated by the gene XDH and the concentration of uric acid was also regulated by XDH.
OLF and calcification of the ligamentum flavum (CLF) are differential diagnoses in patients with posterior extradural compressive lesions. They have similar clinical features, but different pathological manifestations. OLF tends to arise from the lateral capsular portion of the ligamentum flavum and is continuous with the bony laminae . Hypertrophic and degenerative changes of elastic fiber are accompanied by ossified lesions predominantly at the surface of the capsular portion. However, CLF is a crystal deposition disease that mainly affects the central portion of the ligamentum flavum. Calcium deposits mainly occur in the central part of the ligamentum flavum, which is surrounded by degenerated elastic fibers . CT can facilitate the detection of characteristic findings and help distinguish between CLF and OLF . In our study, all patients were diagnosed with thoracic OLF rather than CLF.
OLF and ossification of the posterior longitudinal ligament (OPLL) are diseases which occur ossification in spinal ligament. Both of them result in the compression of the spinal cord and nerve roots. Thus patients may have a similar clinical history that a slowly progressive of neurological symptoms from discomfort to severe myelopathy. And Kawaguchi et al. suggested that over half of the patients with cervical OPLL had OLF . However, there are some differences between the two diseases. The most frequent level of OPLL was at C5 vertebral level while OLF was predominant at upper and lower thoracic levels . In our study, only 2 of 25 patients combined with OPLL. Thus our conclusion only applies to OTLF. Although some studies reported the potential pathogenesis of OTLF such as mechanical, degenerative and genetic factors, the causes of the disease continues to be poorly understood. As far as we know, this study is the first to discover XDH might participate in the progress of OTLF.
In our study, we identified 37 metabolites in the serum of all participants based on the UPLC-MS platform by using untargeted approaches. Both groups demonstrated clusters separated in OPLS-DA models. Metabolite enrichment analysis indicated that the purine metabolism pathway has changed greatly in OTLF patients. ROC analysis indicated that uric acid with the AUC value of over 0.7 might have a diagnosis value in OTLF. Combined with transcriptomic data, 3056 different expression genes were detected and GSEA indicated that purine metabolism might have an impact on the process of OTLF. KEGG pathway suggested that low expression of XDH might increase the level of hypoxanthine and suppress the hypoxanthine metabolism to uric acid that leads to the low level of uric acid in the OTLF patients. All these revealed that uric acid might be the potential biomarker for OTLF.
Previous studies report that osteosis has a relation with metabolism. Davis et al. have suggested that heterotopic ossification might be affected by several metabolism pathways such as tricarboxylic acid (TCA) cycle, amino acid metabolism and lipolysis . Ma et al. suggested that amino acid metabolism and lipid metabolism plays an important role in bone resorption and bone formation through the regulation by strontium salt . Both studies indicate that metabolomics techniques can be used to explore the pathogenesis of OTLF. Using metabolomics related techniques, Sohn et al. suggested that bone mineral density has a positive correlation with OTLF . Fan et al. concluded that obesity was a risk factor for OTLF after determining that leptin-stimulated cell osteogenesis is regulated by STAT3, Runx2 and steroid receptor coactivator-1 . Several studies have focused on the element metabolites related to bone remodeling, but to the best of our knowledge, there have been no previous studies exploring serum metabolites in OTLF. We have found that the downregulation of gene XDH might be the possible mechanism of OTLF, might provide a novel direction of treatment. XDH catalyzes the successive oxidation of hypoxanthine to xanthine and xanthine to uric acid . The decrease of XDH expression leads to the accumulation of hypoxanthine and the decrease of uric acid. Uric acid has been studied extensively in other fields such as Parkinson disease , cardiovascular function . It has also been studied in bone formation, by stimulating osteoblasts and osteoclasts. Yan et al. suggested that acid was negatively correlated with bone formation markers in postmenopausal females . One possible explanation could be that the low level of uric acid might promote bone formation, especially the heterotopic ossification. Future studies could explore the relationship between uric acid and heterotopic ossification.
Our study’s main limitation is the relatively small sample size, and thus further studies with a larger sample size could be considered, especially focusing on the metabolites of the ligament tissue and targeted metabolic approaches. Despite this limitation, our study of OTLF based on metabolomics is the first study of its kind and offers a new diagnostic approach in the early stage of OTLF.
To the best of our knowledge, this is the first study of the relationship of the metabolome with OTLF. We have found that the process of OTLF may be associated with the level of hypoxanthine and uric acid. Uric acid might be the potential biomarker for OTLF and play an important role within the detailed pathway. XDH may regulate the level of uric acid through the purine metabolism pathway and this new discovery may in turn provide a new direction for the diagnosis and treatment of OTLF.
Serum samples were obtained from an experimental group of 25 patients with OTLF and a control group of 23 healthy volunteers. All samples were acquired from Peking Union Medical College Hospital. Written informed consent was obtained from all patients and the study was approved by the hospital institutional review board (JS-981). The inclusion criterion of the experimental group was a confirmed radiological diagnosis of OTLF. All patients had a slowly progressive history of neurological symptoms and signs which is suspected to be caused by thoracic stenosis (TSS). The axial plain CT scan of the thoracic vertebrae and sagittal reconstruction provides analysis of the morphology and density of ossification. MRI can show signal changes and deformation of the spinal cord. The combination of these two imaging examination methods can differentially diagnose OTLF from other causes of TSS (such as OPLL). The exclusion criteria include (1) participants had any history of spinal deformity, other spinal disease or trauma were removed; (2) participants had any systemic metabolic diseases (such as skeletal fluorosis) were removed; (3) participants who were not willing to sign the written consent form were removed; (4) the cases were removed where clinical information was lacking or missing and therefore statistical analysis could not be performed. The healthy controls were recruited from among healthy subjects during an annual health check at Peking Union Medical College Hospital. Standard meals was provided 3 days before sampling.
Morning fasting blood samples were taken from a peripheral vein and collected into ethylenediaminetetraacetic acid (EDTA) tubes. Samples were stored at 4 °C for 3 h before centrifugating at 3000 r/min at room temperature for 20 min. The supernatants were separated and stored at − 80 °C until further analyses.
Before further analysis, the samples were thawed at room temperature. 100 μL serum was extracted and mixed with 400 μL methanol. The mixture was vortexed for 5 min and then centrifuged at 13000 r/min for 15 min. 400 μL of supernatants were then taken for UPLC-MS analysis.
UPLC-MS analysis was performed using a Waters ACQUITY system (Waters Corporation, Milford, USA) combined with a Thermo Scientific high-resolution mass spectrometer (MS) system (Thermo Fisher Scientific Inc., San Jose, USA). C18 column (Acquity UPLC BEH C18–2.1 × 100 mm, 1.7 μm) was used to achieve the separation purpose. 0.2% formic acid solution was used for mobile phase A, pure acetonitrile was used for mobile phase B. Linear gradient elution methods was used for separation. Each sample was injected three times with 5 μL on each occasion. The temperature of the column oven was held at 40 °C and the flow rate was controlled at 0.4 mL/min. Both positive and negative ion mode was applied for data collection of MS.
Metabolites and Transcriptomic data analysis
MarkerLynx (Waters, USA) was used for peak finding, filtering and alignment in the original spectrums. SIMCA-P + 12.0 (Umetrics, Umea, Sweden) was used to perform statistical analyses between two groups such as principal component analysis (PCA) and orthogonal partial least squares (OPLS). MedCalc® Version 18.104.22.168 software was used to plot the receiver operating characteristic curve (ROC) analysis. The differential metabolites statistic analyses were carried out by SPSS 16.0. Metabolites were mapped to the KEGG (Kyoto Encyclopedia of Genes and Genomes Database)  and HMDB (The Human Metabolome Database) . Transcriptomic data GSE69787 was acquired from GEO (Gene Expression Omnibus) platform . R project  was used to ascertain the differentially expressed metabolites and genes. Volcano plot and heatmaps were also carried out through R project.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
Ossification of the Thoracic Ligamentum Flavum
Ossification of the Ligamentum Flavum
Ossification of the Spinal Ligament
Calcification of the Ligamentum Flavum
Ossification of the Posterior Longitudinal Ligament
Liquid Chromatography-Mass Spectrometry
Ultra-high Performance Liquid Chromatography-Mass Spectrometry
Principal Component Analysis
Orthogonal Partial Least Squares
Kyoto Encyclopedia of Genes and Genomes
Variable Importance in Projection
Receiver Operating Characteristic
Area Under Curve
Magnetic Resonance Imaging
Nuclear Magnetic Resonance
Fong SY, Wong HK. Thoracic myelopathy secondary to ligamentum flavum ossification. Ann Acad Med Singap. 2004;33(3):340–6.
Gao R, Yuan W, Yang L, Shi G, Jia L. Clinical features and surgical outcomes of patients with thoracic myelopathy caused by multilevel ossification of the ligamentum flavum. Spine J. 2013;13(9):1032–8.
Guo JJ, Luk KD, Karppinen J, Yang H, Cheung KM. Prevalence, distribution, and morphology of ossification of the ligamentum flavum: a population study of one thousand seven hundred thirty-six magnetic resonance imaging scans. Spine (Phila Pa 1976). 2010;35(1):51–6.
Inamasu J, Guiot BH. A review of factors predictive of surgical outcome for ossification of the ligamentum flavum of the thoracic spine. J Neurosurg Spine. 2006;5(2):133–9.
Zhou SY, Yuan B, Chen XS, Li XB, Zhu W, Jia LS. Imaging grading system for the diagnosis of dural ossification based on 102 segments of TOLF CT bone-window data. Sci Rep. 2017;7(1):2983.
Yang Z, Xue Y, Zhang C, Dai Q, Zhou H. Surgical treatment of ossification of the ligamentum flavum associated with dural ossification in the thoracic spine. J Clin Neurosci. 2013;20(2):212–6.
Kang KC, Lee CS, Shin SK, Park SJ, Chung CH, Chung SS. Ossification of the ligamentum flavum of the thoracic spine in the Korean population. J Neurosurg Spine. 2011;14(4):513–9.
Monteiro MS, Carvalho M, Bastos ML, Guedes de Pinho P. Metabolomics analysis for biomarker discovery: advances and challenges. Curr Med Chem. 2013;20(2):257–71.
De Cecco CN, Caruso D, Schoepf UJ, Wichmann JL, Ter Louw JR, Perry JD, Picard MM, Schaefer AR, Parker LW, Hardie AD. Optimization of window settings for virtual monoenergetic imaging in dual-energy CT of the liver: a multi-reader evaluation of standard monoenergetic and advanced imaged-based monoenergetic datasets. Eur J Radiol. 2016;85(4):695–9.
Martin FP, Ezri J, Cominetti O, Da Silva L, Kussmann M, Godin JP, Nydegger A. Urinary metabolic phenotyping reveals differences in the metabolic status of healthy and inflammatory bowel disease (IBD) children in relation to growth and disease activity. Int J Mol Sci. 2016;17(8):1310.
Jing F, Hu X, Cao Y, Xu M, Wang Y, Jing Y, Hu X, Gao Y, Zhu Z. Discriminating gastric cancer and gastric ulcer using human plasma amino acid metabolic profile. IUBMB Life. 2018;70(6):553–62.
Bhattacharyya S, Pence L, Yan K, Gill P, Luo C, Letzig LG, Simpson PM, Kearns GL, Beger RD, James LP. Targeted metabolomic profiling indicates structure-based perturbations in serum phospholipids in children with acetaminophen overdose. Toxicol Rep. 2016;3:747–55.
Zhang X, Lin Q, Chen J, Wei T, Li C, Zhao L, Gao H, Zheng H. High glucose-induced cardiomyocyte death may be linked to unbalanced branched-chain amino acids and energy metabolism. Molecules. 2018;23(4):807.
Zhang X, Quinn K, Cruickshank-Quinn C, Reisdorph R, Reisdorph N. The application of ion mobility mass spectrometry to metabolomics. Curr Opin Chem Biol. 2018;42:60–6.
Cai HX, Yayama T, Uchida K, Nakajima H, Sugita D, Guerrero AR, Yoshida A, Baba H. Cyclic tensile strain facilitates the ossification of ligamentum flavum through beta-catenin signaling pathway: in vitro analysis. Spine (Phila Pa 1976). 2012;37(11):E639–46.
Lang N, Yuan HS, Wang HL, Liao J, Li M, Guo FX, Shi S, Chen ZQ. Epidemiological survey of ossification of the ligamentum flavum in thoracic spine: CT imaging observation of 993 cases. Eur Spine J. 2013;22(4):857–62.
Liu Y, Zhao Y, Chen Y, Shi G, Yuan W. RUNX2 polymorphisms associated with OPLL and OLF in the Han population. Clin Orthop Relat Res. 2010;468(12):3333–41.
Kawaguchi Y. Biomarkers of ossification of the spinal ligament. Global Spine J. 2019;9(6):650–7.
Fan D, Chen Z, Chen Y, Shang Y. Mechanistic roles of leptin in osteogenic stimulation in thoracic ligament flavum cells. J Biol Chem. 2007;282(41):29958–66.
Ono K, Yonenobu K, Miyamoto S, Okada K. Pathology of ossification of the posterior longitudinal ligament and ligamentum flavum. Clin Orthop Relat Res. 1999;359:18–26.
Takahashi T, Hanakita J, Minami M. Pathophysiology of calcification and ossification of the ligamentum flavum in the cervical spine. Neurosurg Clin N Am. 2018;29(1):47–54.
Kawaguchi Y, Nakano M, Yasuda T, Seki S, Hori T, Kimura T. Ossification of the posterior longitudinal ligament in not only the cervical spine, but also other spinal regions: analysis using multidetector computed tomography of the whole spine. Spine (Phila Pa 1976). 2013;38(23):E1477–82.
Kawaguchi Y, Nakano M, Yasuda T, Seki S, Hori T, Suzuki K, Makino H, Kimura T. Characteristics of ossification of the spinal ligament; incidence of ossification of the ligamentum flavum in patients with cervical ossification of the posterior longitudinal ligament - analysis of the whole spine using multidetector CT. J Orthop Sci. 2016;21(4):439–45.
Davis EL, Salisbury EA, Olmsted-Davis E, Davis AR. Anaplerotic accumulation of tricarboxylic acid cycle intermediates as well as changes in other key metabolites during heterotopic ossification. J Cell Biochem. 2016;117(4):1044–53.
Ma B, Li X, Zhang Q, Wu D, Wang G, Jiye A, Sun J, Li J, Liu Y, Wang Y, et al. Metabonomic profiling in studying anti-osteoporosis effects of strontium fructose 1,6-diphosphate on estrogen deficiency-induced osteoporosis in rats by GC/TOF-MS. Eur J Pharmacol. 2013;718(1–3):524–32.
Sohn S, Yoon JW, Chung CK. Increased bone mineral density in patients with ossification of the ligamentum flavum: a case-control study. J Clin Densitom. 2014;17(1):195–9.
Wang CH, Zhang C, Xing XH. Xanthine dehydrogenase: an old enzyme with new knowledge and prospects. Bioengineered. 2016;7(6):395–405.
Kim IY, O'Reilly EJ, Hughes KC, Gao X, Schwarzschild MA, Hannan MT, Betensky RA, Ascherio A. Integration of risk factors for Parkinson disease in 2 large longitudinal cohorts. Neurology. 2018;90(19):e1646–53.
Kim SC, Shah NR, Rogers JR, Bibbo CF, Di Carli MF, Solomon DH. Assessment of coronary vascular function with cardiac PET in relation to serum uric acid. PLoS One. 2018;13(2):e0192788.
Yan DD, Wang J, Hou XH, Bao YQ, Zhang ZL, Hu C, Jia WP. Association of serum uric acid levels with osteoporosis and bone turnover markers in a Chinese population. Acta Pharmacol Sin. 2018;39(4):626–32.
Kyoto Encyclopedia of Genes and Genomes Database. http://www.genome.jp/kegg/kegg1.html. 2018. Accessed 25 Apr 2018.
The Human Metabolome Database. http://www.hmdb.ca/. 2018.
Gene Expression Omnibus. https://www.ncbi.nlm.nih.gov/geo/. 2018. Accessed 25 Apr 2018.
R project. https://www.r-project.org/. 2018. Accessed 25 Apr 2018.
We thank the patients for their participation and contributions in this present study.
This work was supported by grants from the National Natural Science Foundation of China (NO.81572093 and 81741123). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Written informed consent was obtained from all patients. The study was approved by the hospital institutional review board of Peking Union Medical College Hospital, and complied with the Declaration of Helsinki (JS-981).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, J., Yu, L., Guo, S. et al. Identification of the molecular mechanism and diagnostic biomarkers in the thoracic ossification of the ligamentum flavum using metabolomics and transcriptomics. BMC Mol and Cell Biol 21, 37 (2020). https://doi.org/10.1186/s12860-020-00280-3