- Research article
- Open Access
A diagnostic biomarker profile for fibromyalgia syndrome based on an NMR metabolomics study of selected patients and controls
BMC Neurology volume 17, Article number: 88 (2017)
Fibromyalgia syndrome (FMS) is a chronic pain syndrome. A plausible pathogenesis of the disease is uncertain and the pursuit of measurable biomarkers for objective identification of affected individuals is a continuing endeavour in FMS research. Our objective was to perform an explorative metabolomics study (1) to elucidate the global urinary metabolite profile of patients suffering from FMS, and (2) to explore the potential of this metabolite information to augment existing medical practice in diagnosing the disease.
We selected patients with a medical history of persistent FMS (n = 18), who described their recent state of the disease through the Fibromyalgia Impact Questionnaire (FIQR) and an in-house clinical questionnaire (IHCQ). Three control groups were used: first-generation family members of the patients (n = 11), age-related individuals without any indications of FMS or related conditions (n = 10), and healthy young (18–22 years) individuals (n = 20). All subjects were female and the biofluid under investigation was urine. Correlation analysis of the FIQR showed the FMS patients represented a well-defined disease group for this metabolomics study. Spectral analyses of urine were conducted using a 500 MHz 1H nuclear magnetic resonance (NMR) spectrometer; data processing and analyses were performed using Matlab, R, SPSS and SAS software.
Results and discussion
Unsupervised and supervised multivariate analyses distinguished all three control groups and the FMS patients, and significant increases in metabolites related to the gut microbiome (hippuric, succinic and lactic acids) were observed. We have developed an algorithm for the diagnosis of FMS consisting of three metabolites — succinic acid, taurine and creatine — that have a good level of diagnostic accuracy (Receiver Operating Characteristic (ROC) analysis — area under the curve 90%) and on the pain and fatigue symptoms for the selected FMS patient group.
Our data and comparative analyses indicated an altered metabolic profile of patients with FMS, analytically detectable within their urine. Validation studies may substantiate urinary metabolites to supplement information from medical assessment, tender-point measurements and FIQR questionnaires for an improved objective diagnosis of FMS.
Fibromyalgia syndrome (FMS) is a common chronic pain syndrome characterized by widespread musculo-skeletal pain and associated with multiple other symptoms such as cognitive impairment, disrupted sleep and chronic fatigue. The American College of Rheumtology (ACR) first published criteria for FMS in 1990  which emphasized chronic widespread musculo-skeletal pain (including pain in the axial skeleton) in the presence of pain on at least 11 of 18 specified tender point sites with digital palpation of 4 kg/cm2.
The 2010 ACR updated criteria for FMS  are applied in a 2-part, self-administered questionnaire and do not require a tender point assessment. The first part assesses the presence of pain at 19 sites on a body diagram (widespread pain index) and part 2 measures the symptom severity score (0–3) of 3 core symptoms (insomnia, fatigue and cognitive impairment) and an average score (0–3) for additional somatic symptoms. FMS is the most common cause of widespread or generalized musculo-skeletal pain and affects 2–8% of the adult population with the highest prevalence in women between 30 and 55 years. [3, 4].
FMS is currently viewed as a central sensitivity syndrome associated with abnormal pain processing. It is regarded as a “pain amplification syndrome” associated with increased sensitivity of the nervous system and decreased anti-nociception which results in the clinical phenomena of hyperalgesia and allodynia. Dysfunction in central mono-aminergic neurotransmission which involves serotonin, norepinephrine, nerve growth factor, substance P and others have been implicated in the patho-physiology of FMS. [5,6,7,8] FMS patients often have associated comorbidities such as irritable bowel syndrome, interstitial cystitis and mood disorders [9, 10].
In the absence of an objective biomarker, the diagnosis of FMS is based on a comprehensive clinical assessment. Before 2010, the diagnosis was principally based on the 1990 ACR criteria of widespread pain (including in the axial skeleton) > 3 months and at least 11 painful “tender points” with digital palpation. Although the 2010 ACR criteria do not include a “tender point” count, a musculo-skeletal clinical examination remains mandatory, to exclude other couses of widespread pain and also to identify peripheral pain generators e.g. myofascial trigger points. Selective use of laboratory testing is used to exclude other causes of widespread pain such as polymyalgia rheumatica and hypothyroidism.
The pursuit of specific and measurable biomarkers that may assist in objectively identifying susceptible individuals, confirming disease diagnosis and facilitating treatment, is a continuing endeavour in FMS research. The development of high-throughput metabolic profiling and the study of the metabolome have proven to be particularly applicable in neurological research where small molecules are key in neurochemical metabolism and in performing a role as neurotransmitters, signalling modulators and osmolytes. It is now generally anticipated that metabolomics profiling methods, linked to systems biology approaches, will emerge with well-defined metabolic phenotypes, enhancing the understanding of brain metabolism in health and disease. Recently, a few metabolomics studies have been reported on fibromyalgia, potentially disclosing novel insights into metabolic perturbations in the brain that go brain metabolic homeostasis beyond alterations of neurotransmission variations associated with neurological disorders .
In a pilot study, presented only as a poster at an Annual Meeting of the Rheumatologic Society of the UK , Richards and co-workers (2001) reported that muscle metabolites detected in the urine of fibromyalgia patients may suggest a prevailing muscle damage. Although not by definition a metabolomics study, their targeted metabolite analysis of urine by nuclear magnetic resonance (NMR) spectroscopy revealed significant levels of creatine in FMS patients and elevated (t-test p < 0.05) urinary excretion of choline, taurine, citrate and trimethylamine N-oxide (TMAO) relative to matched controls.
The first metabolomics study on FMS, reported in 2013 , used 50 μl blood samples collected on blood spot cards (Whatman 903 Protein Saver Snap Apart Card, GE Healthcare, Westborough, MA, USA) from patients diagnosed with FMS (n = 14), osteoarthritis (OA; n = 15) and rheumatoid arthritis (RA; n = 12). Samples were dried and then transported to the laboratory for mid-infrared micro-spectroscopy (IRMS) and other analyses. The RA and OA groups appeared to be metabolically similar, but different from the metabolite profile of FMS. The IRMS approach did not conclusively identify the metabolites responsible for the diagnostic spectral differentiation, although changes in tryptophan catabolism seemed to be involved.
Another metabolomics approach to FMS involved liquid chromatography/quadrupole–time-of-flight/mass spectrometry (LC/Q-TOF/MS) with multivariate statistical analysis aimed at discriminating FMS patients (n = 22) and controls (n = 21) from blood plasma analysis . Lysophosphocholine (lysoPCs), phosphocholine and ceramide lipids dominated the metabolite profile. The metabolites that discriminated the most between FMS patients and controls were identified as 1-tetradecanoyl-sn-glycero-3-phosphocholine [PC(14:0/0:0)] and 1-hexadecanoyl-sn-glycero-3-phosphocholine [PC(16:0/0:0)] — suggesting that lysoPCs may be potential biomarkers for FMS.
In addition to these metabolomics findings, a recent review on biomarkers of FMS included contributions from genetic and proteomic studies . Although genetic factors have been shown to influence predisposition to FMS, no specific genes have been confirmed as being involved in this disease . The review also listed several proteins of the immune response, cytoskeleton remodelling and the inflammatory process in FMS. Their role in FMS, however, is still controversial.
Thus the availability of biomarkers for unequivocal and objective diagnosis of FMS remains elusive in clinical practice. Yet, metabolites identified as being involved in the aetiology and pathogenesis of FMS could meanwhile contribute to insights into various presentations of FMS and provide ancillary diagnostic testing criteria to complement general diagnostic procedures. We thus present here the outcomes of a 1H NMR metabolomics study on FMS. All experimental subjects were females and provided urine samples for the study. The investigation was designed as an untargeted approach and revealed metabolite information with predictive potential to discriminate between FMS patients and healthy young controls. The outcomes thus underscore the versatility of metabolomics to provide insights into disease pathophysiology, furthering potential novel approaches to supplement existing protocols proposed for the practising clinician to assess FMS and monitor its treatment .
Experimental subjects, physical characteristics, symptoms and clinical profiles
All the patients that were included in this study were previously diagnosed with FMS by the same specialist pain clinician from his chronic pain practice in Pretoria. This practice manages the full spectrum of chronic pain disorders, with a special interest in FMS and related pain disorders. The diagnosis was based on a comprehensive clinical assessment using the 1990 criteria. All patients in the study were confirmed with FMS before 2010 and all were on a comprehensive evidence-based management programme according to international guidelines. They were only included if they continued to complain of widespread musculo-skeletal pain (including in the axial skeleton) in the presence of >11 painful tender points with musculo-skeletal assessment.
Informed consent was obtained from all the participants in this study by means of a voluntarily completed consent form; ethical approval for the study was obtained as specified under Declarations. All participants in the study were female and the sample material investigated was urine. The experimental subjects consisted of one FMS patient group (Group 1) and three control groups (Groups 2 to 4). Clinical description and serum and urine sample collection on all experimental groups commenced from 2009 to 2011. Case definition and selection for the eventual study was done by clinical and scientific group of co-workers in 2010. Following scrutinizing of the records of patients with a medical history of FMS, a group of 18 FMS patients eventually selected based on the above selection criteria as well as after excluding outliers based on statistical analysis [(S5 in the Additional file 1: Supplementary Information (SI)].
The socio-demographic, tender point and myofascial pain experience, awareness of gastro-intestinal symptoms, pain-specific medication and levels of emotional experience associated with FMS for the 18 patients was obtained through the FIQR  and the IHCQ. The questionnaires are presented in Table S1 and the response to the IHCQ are summarised in Table S2 of the SI. The IHCQ included 18 items that could be extended to a total of 30 sub-items. The questionnaire provided socio-demographic information on the patients (3 items), personal clinical experience of pain and their FMS condition (5 items) and use of medication against pain (2 questions), emotional experience (7 questions) and digestive functioning (1 item). The urine samples were provided by the patients prior to application of pressure to the tender-points (TPs). For a total of 16 of these FMS patients a complete set of data was available for the comparative analysis of the FIQR and metabolomics data, as some information on some patients had to be excluded because the data were incomplete. Some degree of comorbidity of conditions that overlap with FMS (e.g. chronic fatigue syndrome) could not be excluded, as the mean level of energy in the FIQR was rated at 7,0 and according to responses to the IHCQ, 94% of the patients experienced sleep disturbances and did not awoke refreshed. The responses to an experience of mood disturbances (58% answered “Yes”) and anxiety (52% answered “Yes”) for the FMS patients as a group were moderate. Responses on depression was inconsistent (mean FIQR-score = 5.1 with 84% “Yes” answers on the IHCQ) but 88% indicated discomforts with their gastrointestinal functions (Indicated as Irritable Bowel Syndrome (IBS) in the IHCQ). These scores were accepted as indications of the mental and physical profile of the FMS patient group and were not further clinically verified.
Three control groups were used: (1) a group of 11 subjects that were first-degree relatives of the patients, meaning that they were a mother, sister or daughter relation (Group 2: CF); (2) a group of 10 unrelated subjects, selected by physicians and defined as unrelated and age matched controls to the patients (Group 3: CO); (3) a control group of young and healthy individuals, comprising 20 randomly selected students (aged 18–22 years) of North-West University (NWU) (Group 4: CN). All individuals in the control groups showed no indications of FMS or related conditions and was not required to complete the FIQR or IHCQ.
This investigation used availability sampling on the clinically selected FMS patients and controls (CO, CF and CN). However, statistical analyses indicated that the sample sizes provided sufficient power to detect large effects at a univariate level.in the FMS and CN comparison.
Sample preparation and 1H NMR analysis
Spectral analyses were conducted according to the protocol at the NMR facility of the Translational Metabolic Laboratory at Radboud University Medical Centre in Nijmegen, the Netherlands [18, 19]. The urine samples were collected in South Africa, stored at –80 °C and transported to the Netherlands before being thawed at room temperature for analysis. A 1 ml volume of each sample was centrifuged at 3000 rpm for 10 min to remove any sediments or debris. A 70 μl volume of a deuterated solution containing 20.2 mM of trimethyl-2,2,3,3-tetradeuteropropionic acid (TSP, sodium salt; Sigma Aldrich) was added to 700 μl of the supernatant and vortexed. This internal standard (IS) solution served to lock the signal during analysis and to provide a chemical shift reference of δ = 0.00. The sample was then acidified to pH 2.5 ± 0.05, with 37% concentrated hydrochloric acid (HCl). A 650 μl aliquot of the acidified sample was then transferred to a 5 mm NMR tube (Wilmad Royal Imperial; Wilmad LabGlass, USA) and analysed on a 500 MHz Bruker Avance spectrometer (Bruker Analytische Messtechnik, Karlsruhe, Germany) (pulse angle 90o, delay time 4 s, number of scans 256, temperature 298 K). Water suppression was achieved by using gated irradiation focused on the water frequency. All samples were automatically shimmed prior to acquisition of data, using topshim from Bruker BioSpin. The resultant raw spectral data, in the form of free induction decay, were Fourier transformed. These transformed spectra were then manually corrected for phase and baseline. All the samples were normalized with reference to the creatinine CH3 peak at 3.13 ppm. We opted for two methods of spectral analysis. The first method entailed equidistant binning  using a bin width of 0.02 ppm applied to the selected region of 0.5–10 ppm, which gave a total of 461 integrated units per NMR spectrum, excluding the water region, for each individual of the four experimental groups. The second method entailed variable-sized binning. The equal-binning procedure masks subtle chemical shift differences, hides potentially significant changes of low-intensity peaks and incurs the risk of splitting peaks or spectral features between bins . To avoid these problems we also used variable bin sizes in areas of peaks above the noise level, preventing peak division between multiple bins. This approach was specifically applied for the identification and quantification of discernible and important known metabolites, generating data for univariate analysis.
Data and statistical analysis
The original normalized spectral data (presented in Additional file 1 as Table S4 in S1 – Supplementary information (SI) or Additional file 2 – Raw data matrix) were pre-processed by performing log transformation and auto-scaling. Outliers were detected through Hotelling’s T2 and PCA scores (using a 90% confidence region) analysis and resulted in the exclusion of 4 outliers from further analysis. Univariate statistical analyses, specifically the Mann–Whitney test p-values (MW) and associated effect sizes (ES), were generated for each feature. Multivariate analyses were performed using cluster analysis (Euclidean distance and Ward linkage) principal components analysis (PCA) and partial least squares discriminant analysis (PLS-DA), using a 90% confidence interval (CI). Data processing and analyses were performed using Matlab (MATLAB with Statistics and PLS Toolbox Release 2012b, The MathWorks, Inc., Natick, MA, USA); R (R version 3.2.3 downloaded from https://www.R-project.org with the corrplot package downloaded from https://cran.r-project.org/web/packages/corrplot); the SPSS software package (SPSS Inc. (2015). IBM SPSS Statistics Version 22, Release 22.0.0, © IBM Corporation and its licensors - http://www-01.ibm.com/software/analytics/spss/) and SAS (SAS Institute Inc. 2016 The SAS System for Windows Release 9.4 TS Level 1 M3, SAS Institute Inc., Cary, NC, USA). A table containing all discriminant information, i.e. the power and VIP values as generated from the PCA and PLS-DA analyses, respectively, as well as the ES and MW p-values, was constructed.
We did not test for a normal distribution of the data, given the small number of cases and used Pearson’s r and Spearman’s rho to assess correlations, analysed through SPSS version 12.0 (SPSS, Inc., Chicago, IL). All tests were one-tailed, given the positive fold changes (FC) observed for all metabolites.
Characteristics of the FMS patient group
The age profile of the patients concurs with the general agreement of FMS being uncommon in young subjects (<25–30 years), increasing with age towards the prevalence peak in middle-aged individuals, and then declines [3, 4]. According to the feedback, 88% of our patients had stable relationships with a male partner, 89% had one form or another of day-filling or employment activities, and their emotional experience was not severely affected by their disease. The pain experience and medication used resembled that generally prescribed for FMS. The mean scores and ranges of the 21 FIQR questions obtained for our patient group and those (designated as the reference group) used for the standardization of the questionnaire  are compared in Table S3.
To characterize further the relationship between questions or variables making up the FIQR questionnaire, we calculated Kendall’s tau correlation coefficients for the FMS patient group (Fig. 1). The correlation coefficients along with their associated significance levels are indicated in Table S3. The highest correlation (r = 0.817) was indicated for the relationship between pain and the symptoms for FMS. The function domain contains 9 physical functioning items related to the ability to perform relatively demanding but regular daily muscle tasks. Apart from the low score for ‘brushing hair’, all remaining items showed high correlation coefficients among each other, ranging from 0.399 to 0.778. These high values collectively substantiate the major signs and symptoms experienced by the FMS patients. The ‘overall impact’ domain contained 2 items that asked about the number of days individuals felt well (could reach their goals) and the corresponding number they were unable to work because of FMS symptoms. These again showed high correlation coefficients, ranging from 0.421 to 0.686, with the 8 items in the functional domain indicating the underlying negative impact of the FMS symptoms on the daily routine of the FMS patients. The symptoms domain contained 9 items on which patients had to rate work and physiological, psychological and environmental difficulties related to FMS. Lower correlations, ranging from 0.076 to 0.499 (mean = 0.25), were found between the 8 functional items and sleep patterns, memory, anxiety and depression, indicating little overlap within the patient group with other FMS-related conditions. Taking everything into account, we conclude the FMS patients represent a well-defined group for this explorative metabolomics study.
Data generation and case selection
Representative scaled NMR spectra from an FMS patient and from the young control group (CN) is shown in Fig. 2 to illustrate some of the discernible qualitative NMR differences observed in these selected examples. Close inspection of the spectra indicates that there were no immediately discernible qualitative differences between the two representative examples, suggesting that FMS is not associated with distinctive metabolic aberrations, as otherwise observed in monogenetic disorders such as inborn errors of metabolism. Using the equal-bins spectral data, case reduction was first applied to all four experimental groups (Additional file 1: Figure S2). Four outliers were identified using a 95% confidence region in a Hotelling’s T2 test in conjunction with the respective PCA score plots with 90% confidence regions. Cases that were identified as outliers by either method were removed. The outliers were: group 1 (FMS patients) – one outlier; group 2 (CF; family controls) – two outliers; group 3 (CO; matched controls) – no outliers; group 4 (CN; young controls) – one outlier.
Supposed changes in metabolite profiles from the FMS patients and the three control groups (excluding outliers) were established through three multivariate approaches: unsupervised Euclidian and Ward hierarchical cluster analyses presented as dendrograms, unsupervised PCA, and supervised PLS-DA models, applied to the original 461 1H NMR profiled bins for the four experimental groups.
Figure 3 shows the group separations based on the unsupervised cluster analysis, indicating the perceived closeness of spectral data encapsulated in the NMR bins. The main clusters formed between the CF family members group (Fig. 3a) and the CO age-matched group (Fig. 3b) relative to the FMS patients are heterogeneous in terms of case distribution. In contrast, two well-defined clusters were formed between the FMS patients and CN young controls (Fig. 3c), suggesting distinct differences in the spectral fingerprints between these two groups.
Next, group separations based on unsupervised PCA and supervised PLS-DA were performed. The data were log transformed and auto-scaled. The PCA between the CF family members (Fig. 3d), CO matched controls (Fig. 3e) and FMS patients complemented results from the cluster analyses. A complete separation was obtained between all three control groups and the FMS patients (Fig. 3f–i) through supervised PLS-DA. Evaluation of the PLS-DA model shown in Fig. 3i (FMS vs CN) was performed by calculating the goodness-of-fit (R2) and predictive ability (Q2) parameters. These metrics confirmed the complete separation between the FMS and CN young control groups, with good model fit (R2 = 0. 96), however this model may not generalize well (Q2 = 0.29).
From the equal binning analysis it is evident that there are bins or combinations of bins that can discriminate between our patient and control groups. However, since it is not clinically practical to measure bins, we did not investigate this data further. Instead, the metabolites potentially responsible for the separation of the FMS patients and the CN young controls were subsequently identified by analyzing variable bins from the NMR spectra and converting these measures to concentration values of the identified metabolites.
Metabolite profile of the FMS patient group
Twenty-one metabolites could be identified and quantified from the NMR spectra. From this list we selected twenty endogenous metabolites (listed in Table 1), and also included 2-hydroxyisobutyric acid of exogenous origin , with high VIP, ES and ES values, despite being present in low concentrations. The endogenous metabolites include seven amino acids (tyrosine, leucine, valine, histidine, alanine, threonine and lysine), seven metabolites directly or indirectly associated with energy metabolism (lactic acid, succinic acid, citric acid, 3-hydroxyisovaleric acid, creatine, carnitine and formic acid), three osmolytes (taurine, TMAO and dimethylglycine), a major mammalian detoxification product (hippuric acid), histamine and an N-acetyl-derivative. The N-acetyl-derivative showed a singlet at 2.03 ppm, possibly indicative of an N-acetyl group. One-dimensional spectral data suggested that aspartic acid (multiplet at 4.70 ppm) could be the moiety linked to the N-acetyl group, which, however, could not be substantiated as N-acetyl-aspartic acid by two-dimensional NMR spectral analysis Additional file 1: Figure S4). We thus designated the variable as an N-acetyl derivative (N-acetyl-X).
We subsequently performed multivariate (log and centred concentration values) and univariate (unscaled concentration values) analyses on the reduced bins (endogenous metabolites, converted to their respective metabolite concentrations) to refine our identification of the key variables that discriminate between the FMS patients and the controls. All cases were retained for this analysis as none were identified as outliers based on the concentrations. Multivariate PCA (Fig. 4a) and PLS-DA (Fig. 4b) both indicated that the 20 metabolites contained information that differentiates, but did not separate, the FMS patients from the young controls. Model performance was evaluated using the goodness of fit (R2) and goodness of prediction (Q2) parameters, which were R2(X) = 0.52 and Q2(Y) = 0.05, respectively, indicating a reasonable (R2) but not necessarily reproducible (Q2) fit between the variation in the data and the components (quantified metabolites) comprising the model for the present FMS group. It thus appears that some metabolites below the sensitivity range for quantification from the NMR spectra might be required for reproducibility (Q2) and for further differentiation between the FMS patients and young controls.
Univariate analyses using Mann–Whitney p-values and fold changes, as summarized in a volcano plot (Fig. 5a), point to important substances that cause group differentiation. The outcome of this analysis of the data set of 20 variables is presented in Fig. 5a, indicating which large-magnitude changes (fold change: |log2 FC| ˃ 1.5) are also statistically significant (Mann–Whitney test: p < 0.05). Six informative metabolites complied with these measures, with their respective VIP values shown in brackets: succinic acid (0.246), taurine (5.214), tyrosine (0.365), lactic acid (2.832), creatine (4.402) and trimethylamine N-oxide (TMAO; 2.209).
Important endogenous metabolites in the FMS patient group
A summary of the results for the univariate and multivariate statistical analyses is presented in Table 1, nine of which could be related to physiological functions that could be related to FMS.
The neurological functions of succinic acid, tyrosine and lactic acid are well known: the aerobic mitochondrial energy regeneration function, a precursor for neurotransmitters and a key metabolite in the astrocyte-neuron lactate shuttle , respectively. Taurine is an abundant β-amino acid in the mammalian brain  and has been shown to be a neurotransmitter in the substantia nigra (SN). It has been suggested from micro-dialysis experiments on Sprague-Dawley rats that osmoregulation of the nonsynaptic taurine pool of the SN could influence the nigral cell vulnerability, seen in the pathogenesis of Parkinson’s disease . Likewise, nutritional studies  suggest that TMAO may be involved in diet-induced variations in the balance of several osmolytes, including betaine, choline, creatinine and creatine, whereas creatine has also been proposed as being involved in pain experienced in FMS . Thus, we subsequently evaluated the potential diagnostic value of these six metabolites on FMS by means of a logistic regression analysis, as indicated below.
Important exogenous metabolites in the FMS patient group
The pain intensity of patients with FMS has been reported to correlate with the degree of small intestinal bacterial overgrowth [9, 10]. This clinical observation may have pathogenetic relevance for FMS, because bacterial overgrowth leads to the exposure of immune cells to luminal antigens and consequent immune modulation. An untargeted NMR metabolomics study of celiac disease, a multifactorial immune-mediated enteropathy , suggested alterations of energy metabolism - a clinical characteristic in FMS - while urine data pointed to alterations of gut microbiota. At least three metabolites observed in the urine samples of our FMS patient group suggest perturbations in their gut metabolome (Fig. 6): (1) Hippuric acid is a normal and major component of urine and appear in humans as an increased excretory product from unnatural (detoxification) and natural (gastroesophageal reflux disease in children) sources. (2) 2-Hydroxyisobutyric acid, the most discriminatory variable between our FMS group and controls (VIP = 6.2 – Table 1), is an apparent catabolic from gut microbiotica and was shown to be statistically linked to Faecalibacterium prausnitzii  an important commensal bacterium of the human gut flora proposed to be an indicator of the dynamic basis of host–microbiome symbiosis. (3) Lactic acid is a key intermediate in many biochemical processes and is a measure of critical illness in patients with poor prognosis. It may be of endogenous (L-lactate) or exogenous (D-lactate) origin and we recently proposed that the determination of its enantiomers in infectious conditions may provide a basis for substantiating the clinical significance of disease markers . The presence of these exogenous markers of gut origin provides further indications of the connectivity between disturbances in the gut microbial populations and the metabolic consequences of the altered microbial–mammalian metabolic balance influencing host disease, which will be discussed below in the context of FMS.
A putative biosignature for FMS
A combination of three selection methods (Forward, Backward and Step-wise selection) was used to identify the best metabolite predictors. Instead of using one selection method, a combined approach was chosen since each method has its advantages and disadvantages [30, 31]. Although our aim was to explore a small set of highly discriminatory endogenous metabolites, we also investigated the potential of a combination of these metabolites to function as a biosignature for the FMS patient group. We followed a forced entry approach to evaluate the combination of metabolites. Table 2 lists the methods used as well as the preferred metabolite predictors selected from the six informative metabolites. The last model (Forced entry) entered succinic acid, taurine and creatine, and produced the best model based on −2 Log Likelihood (−2LL) from the present data. Table 2 also reports other model performance measures, but -2LL was used to select the best model as it gives an indication of the variation not explained in the data, and gave the lowest -2LL value compared to the other models. The Forced entry model was also well calibrated since the Hosmer Lemeshow (HL) statistic was not significant. The model fit is reported by using the Max Rescaled R-squared value and only the Forced entry model had a satisfactory value of above 0.6.
Finally, the classification ability of each model was assessed by using a Receiver Operating Characteristic (ROC) analysis to the data mentioned. The values of the area under the ROC curve (AUC) provide a measure of how well this combination could distinguish between the two groups. A value of AUC = 1 represents a perfect test, while a cursory guide for classifying the accuracy of a diagnostic test is given by: AUC = 0.90–1 (excellent, i.e. high sensitivity and high specificity); 0.80–0.90 (good); 0.70–0.80 (fair); 0.60–0.70 (poor); 0.50–0.60 (fail). To provide some indication of how well the model would potentially generalize, the last column in Table 2 reports the classification ability when one sample is left out repeatedly — in other words, based on a leave-one-out cross-validation strategy (AUC (LOO CV)). Again the Forced entry model performed the best (AUC = 90% [0.8972]; AUC (LOOCV) = 88% [0.8750]).
Correlation between clinical and metabolic indicators
Pearson and Spearman correlation analysis was done to compare the bivariate relationships between responses to the FIQR and the three endogenous variables defining the biosignature of FMS. Specifically, correlations were assessed between the sum of all three FIQR domains as well as the sum of the functional, impact and symptoms domains and SUM-3, SUM-2, creatine, succinic acid and taurine. Finally we inspected the data for symptoms related to metabolism to be included in the bivariate correlation analysis. In this regard it should be noted: (1) The scores of the 21 questions of the FIQR corresponds to an average based on the subjective self-assessment of the FMS patients as used in the behavioural sciences (i.e., it is not empirically based). We therefore used the mean scores of fibromyalgia patients on the symptoms for experience of pain, low energy levels and tenderness to touch only as a directive to include these symptoms in the bivariate correlation analysis . Their mean values did not differed in practice from the data of a reference group of the revised FIQR (Additional file 1: Figure S1). (2) The number of FMS cases is relatively small for assessment of normality in the data distribution. We therefore included the Pearson and Spearman correlations in Table 3, but used only the Spearman’s correlations for the interpretation of the bivariate correlation analyses, with guideline values for “small”(r ≥ 0.1), “medium” (r ≥ 0.3) and “large” (r ≥ 0.5) as operational convention for the correlation coefficients .
The results shown indicate a medium and borderline significant relationship between the SUM-3 biosignature and the sum of the FIQR, with insignificant correlations for its functional and impact domains. Sum-3 and the symptoms domain showed a large and significant correlation. The relationship between SUM-2 and the sum of the FIQR and its three domains improved significantly. Taken together these results directs to a more meaningful relationship between the metabolites which comprise the biosignature and clinical symptoms related to biochemical perturbations in FMS. This impression is substantiated by the strong and significant relationship between SUM-3 and SUM-2 on the experience of pain (p = 0.004 and 0.016, respectively) and loss of energy (p = 0.006 and 0.001, respectively) in the FMS patients as a group. Notably this relationship is not shared by succinic acid (a metabolite from the Krebs cycle) and taurine (an osmolite), but a good and significant relationship was shown between creatine and the symptoms pain and energy (p = 0.024 and 0.003, respectively). The relationship between the biosignature components to tenderness to touch, the third clinical symptom evaluated, was statistically insignificant (not included in Table 2). All correlation coefficients were positive indicating that patients with high scores on the biosignature will likely also have high FIQR scores.
The results of this metabolomics study lead to three main discussion points – whether FMS presents with a unique global metabolic profile which characterizes this disease, whether metabolomics studies contributed to the advancement of an objective clinical diagnosis of FMS in patients so affected and on gut microbial–host metabolic perturbations in FMS.
As the overall health status of individuals is captured in their metabolic state, there exists a prevailing view that metabolomics results embody global biochemical changes in an individual due to a disease and neurological conditions , and supported by our results and of two other NMR metabolomics investigations. The first NMR metabolomics study evaluated the diagnostic accuracy of biomarker profiles in three neurological conditions: idiopathic intracranial hypertension, multiple sclerosis, and cerebrovascular disease relative to controls with either no or combined neurological diseases . It appeared that the metabolomics investigation identified differences in metabolite profiles in patients suffering from these three conditions. A related conclusion was drawn from the second NMR metabolomics study of FMS . Although a relatively small number of patients formed the experimental group, the metabolomics approach was successful in identifying distinct metabolic profiles for FMS patients relative to controls, supporting the concept that the Platelet Activating Factor/Platelet Activating Factor Receptor (PAF/PAFr) system plays a role in modulating pain signalling. Our results furthermore indicated the differentiation of the three control groups used (family members, an age-matched group, and young individuals) and the FMS patients (Additional file 1: Figure S3). Statistical assessment of the outcome of a supervised PLS-DA model confirmed the complete separation between the FMS and young control groups. Good model fit values substantiated some unique differences between the global metabolic profiles of the FMS patients and the healthy young controls. The metabolites principally responsible for the differentiation between our FMS patients and controls included taurine and TMAO which were also reported to be significantly increased (p < 0.05) in an FMS patient group in a preliminary targeted NMR study . In addition, we observed perturbed succinic acid suggesting altered energy metabolism in FMS. This result is linked to a study  where there was relatively elevated: glucose, the glycolytic intermediate phosphoenolpyruvate, pyruvate and nicotinamide adenine dinucleotide (NAD+) seen in dried blood spots from FMS patients. This observation was previously reported for patients with chronic widespread pain .
A common thread in the metabolomics studies on FMS discussed here is the affirmation of the ability of metabolomics to identify distinct metabolic profiles for FMS patients relative to controls. Some metabolites/biomarkers could therefore contribute to the disease phenotype by having a role in the pathogenesis of FMS. The biomarkers revealed in these metabolomics studies seemed, however, not to be metabolically closely linked, but may be due to the multi-factorial nature of FMS. Noteworthy also are the two main limitations of our own and the other two metabolomics studies: the FMS groups investigated and analytical methods used. Most metabolomics studies are limited by the number of experimental subjects available for investigation, and therefore ultimately call for follow-up validation studies with larger and better-defined experimental groups. Further, given the complexity of the human metabolome and the multi-dimensional nature of biofluids and other biological samples available for metabolomics studies, no single analytical technology can fully disclose and account for the information encapsulated in these samples. Nonetheless, metabolomics retains a promise well beyond the scope of standard clinical chemistry techniques, for affording detailed characterization of metabolic phenotypes and is believed, eventually, to lead to so-called precision medicine in which knowledge of their unique metabolic derangements explains the disease state of individual patients . A third limitation in the present study is the use of the 1990 criteria for FMS (1, 14) as the patient selection was one before publication of the revised criteria in 2011. The use of the revised criteria is now standard practice in our pain clinic.
So, can metabolomics studies contribute to the advancement of objective clinical diagnosis of FMS? The results of the present and the two other metabolomics studies on the disease imply that they can, albeit with qualifications. The analyses of blood spots from FMS patients provided information using IRMS technology that differentiated samples from FMS subjects from those with RA or OA with zero misclassifications (100% accuracy). The accuracy of the metabolomics approach was 75%, but with the advantage of disclosing a prioritized list of metabolites that may underlie the differences identified . The possible role of lysoPCs as biomarkers or as contributors to the FMS phenotype and function in the pathogenesis of this condition suggest they are potential new disease biomarkers and thereby open a new approach for the treatment in FMS . Likewise, the predictive potential of the combination of succinic acid, taurine and creatin
e proved to be excellent for discriminating between our cases of FMS and controls (AUC = 90%). The combination of creatine and succinic acid also showed a significant correlation with the characteristic symptoms of pain and fatigue in FMS. The inclusion of this predictive information on these three metabolites could in time be considered to form part of the initial evaluation of patients suspected of suffering from the disease, in anticipation of validation of FMS diagnostic markers.
Finally, the involvement of gut microbial–host metabolic perturbations in FMS may prove to contribute significantly in defining the clinical profile in FMS. In health, brain-gut interactions are crucial in the maintaining of homeostasis . It appears that neuroplasticity-related systems and neurotransmitter systems are influenced by the gut–brain axis regulation and perturbed homeostasis may contribute to risk of disease through alterations in gastrointestinal tract, central nervous, autonomic nervous and immune systems . The frequent comorbidity of fibromyalgia with stress related disorders, such as chronic fatigue and irritable bowel syndromes and some CNS related abnormalities, suggests that gut–brain axis regulation may at least be a partial common denominator for these disorders. This view may well be revealed by data from a follow-up targeted metabolomics investigation of high sensitivity, like through mass spectrometric-based technologies.
An untargeted 1H NMR metabolomics analysis of urine samples obtained from a group of clinically well-defined female FMS patients with no psychiatric co-morbidity could be fully differentiated from a group of young healthy women. The presence of metabolic indicators of perturbations in the gut microbiome (hippuric, 2-hydroxyisobutyric and lactic acids) supports the paradigm that regulation of the gut-brain axis becomes affected in stress related disorders, like FMS. Three metabolite markers (taurine, creatine and succinic acid) were important for the differentiation between FMS patients and controls and were significant indicators of the pain and fatigue symptoms in FMS. ROC analysis and odds ratios substantiated the good predictive potential of a combination of these three metabolites for FMS in the present patient group. Follow-up metabolomics research on a larger number of urine samples, including those from individuals at high risk of developing the disease, as well as longitudinal studies on FMS patients during treatment, are needed to validate the findings presented here and to potentially detect effects which would require greater statistical power. These markers may in time provide objective supplementary information together with tender-point measurements and FIQR questionnaires used to confirm FMS.
Proton nuclear magnetic resonance
American College of Rheumatology
Area under the ROC curve
Area under the ROC curve
Group of first-degree relatives of the patients
Group of healthy young subjects
Central nervous system
Group of age-matched subjects but unrelated to the patients
- CV AUROC:
Chronic widespread pain
Fibromyalgia Impact Questionnaire
In-house clinical questionnaire
liquid chromatography/quadrupole–time of flight/mass spectrometry
- LOO CV:
Nuclear Technologies in Medicine and Biosciences Initiative
Platelet Activating Factor/Platelet Activating Factor Receptor
Principal components analysis
Partial least squares discriminant analysis
- Q2 :
Predictive ability parameter
- R2 :
Receiver Operating Characteristic
Numerical sum of the concentrations of taurine, creatine and TMAO
Variable important in projection
Wolfe F, Smythe HA, Yunus MB, Bennett RM, Bombardier C, Goldenberg DL, et al. The American College of Rheumatology 1990 criteria for the classification of fibromyalgia. Arthritis Rheum. 1990;33(2):160–72.
Wolfe F, Clauw DJ, Fitzcharles MA, Goldenberg DL, Katz RS, Mease P, et al. The American College of Rheumatology preliminary diagnostic criteria for fibromyalgia and measurement of symptom severity. Arthritis Care Res. 2010;62(5):600–10.
Branco JC, Bannwarth B, Failde I, Carbonell JA, Blotman F, Spaeth M, Saraiva F, et al. Prevalence of fibromyalgia: a survey in five European countries. Semin Arthritis Rheum. 2010;39:448–53.
Lawrence RC, Felson DT, Helmick CG, Arnold LM, Choi H, Deyo RA, et al. Estimates of the prevalence of arthritis and other rheumatic conditions in the United States: part II. Arthritis Rheum. 2008;58(1):26–35.
Russell IJ, Orr MD, Littman B, Vipraio GA, Alboukrek D, Michalek JE, et al. Elevated cerebrospinal fluid levels of substance P in patients with the fibromyalgia syndrome. Arthritis Rheum. 1994;37:1593–601.
Yunus MB, Dailey JW, Aldag JC, Masi AT, Jobe PC. Plasma and urinary catecholamines in primary fibromyalgia: a controlled study. J Rheumatol. 1992;19:95–7.
Legangneux E, Mora JJ, Spreux-Varoquaux O, Thorin I, Herrou M, Alvado G, et al. Cerebrospinal fluid biogenic amine metabolites, plasma-rich platelet serotonin and [3H]imipramine reuptake in the primary fibromyalgia syndrome. Rheumatol. 2001;40:290–6.
Meyer HP. Myofascial pain syndrome and its suggested role in the pathogenesis and treatment of fibromyalgia syndrome. Curr Pain Headache Rep. 2002;6(4):274–83.
Levine JS, Burakoff R. Extraintestinal manifestations of inflammatory bowel disease. Gastroenterol Hepatol (NY). 2011;7(4):235–41.
Clauw DJ. Fibromyalgia and related conditions. Mayo Clin Proc Elsevier. 2015;90(5):680–92.
Dumas M-E, Davidovic L. Metabolic phenotyping and systems biology approaches to understanding neurological disorders. F1000Prime Rep. 2013;5:18.
Richards SCM, Bell J, Cheung Y-L, Cleare A, Scott DL. Abstract 382: muscle metabolites detected in urine in fibromyalgia and chronic fatigue syndrome may suggest ongoing muscle damage. 2001.
Hackshaw KV, Rodriguez-Saona L, Plans M, Bell LN, Buffington CAT. A bloodspot-based diagnostic test for fibromyalgia syndrome and related disorders. Analyst. 2013;138:4453–62.
Caboni P, Liori B, Kumar A, Santoru ML, Asthana S, Pieroni E, et al. Metabolomics analysis and modeling suggest a Lysophosphocholines-PAF receptor interaction in fibromyalgia. PLoS One. 2014;9(9):e107626.
Giacomelli C, Sernissi F, Rossi A, Bombardieri S, Bazzichi L. Biomarkers in fibromyalgia: a review. Curr Biomark Find. 2014;4:35–41.
Xiao Y, Russell IJ, Liu YG. A brain-derived neurotrophic factor polymorphism Val66Met identifies fibromyalgia syndrome subgroup with higher body mass index and C-reactive protein. Rheumatol Int. 2012;32:2479–85.
Bennett RM, Friend R, Jones KD, Ward R, Han BK, Ross RL. The revised fibromyalgia impact questionnaire (FIQR): validation and psychometric properties. Arthritis Res Ther. 2009;11(5):415.
Engelke UF, Liebrand-van Sambeek ML, De Jong JG, Leroy JG, Morava É, Smeitink JA, et al. N-acetylated metabolites in urine: proton nuclear magnetic resonance spectroscopic study on patients with inborn errors of metabolism. Clin Chem. 2004;50(1):58–66.
Wevers RA, Engelke U, Rotteveel JJ, Heerschap A, De Jong JG, Abeling NG, et al. 1H NMR spectroscopy of body fluids in patients with inborn errors of purine and pyrimidine metabolism. J Inherit Metab Dis. 1997;20(3):345–50.
Viant MR, Lyeth BG, Miller MG, Berman RF. An NMR metabolomic investigation of early metabolic disturbances following traumatic brain injury in a mammalian model. NMR Biomed. 2005;18:507–16.
Worley B, Powers R. Multivariate analysis in metabolomics. Curr Metabolomics. 2013;1:92–107.
Ditre CM, Griffin TD, Murphy GF, Sueki H, Telegan B, Johnson WC, et al. Effects of 2-hydroxy acids on photoaged skin: a pilot clinical, histologic, and ultrastructural study. J Am Acad Dermatol. 1996;34:187–95.
Pellerin L, Pellegri G, Bittar PG, Charnay Y, Bouras C, Martin JL, et al. Evidence supporting the existence of an activity-dependent astrocyte-neuron lactate shuttle. Dev Neurosci. 1998;20:291–9.
Huxtable RJ. Physiological actions of taurine. Physiol Rev. 1992;72:101–63.
Morales I, Dopico JG, Sabate M, Gonzalez-Hernandez T, Rodriguez M. Substantia nigra osmoregulation: taurine and ATP involvement. Am J Physiol Cell Physiol. 2007;292(5):C1934–41.
Solanky KS, Bailey NJ, Beckwith-Hall BM, Bingham S, Davis A, et al. Biofluid 1H NMR-based metabonomic techniques in nutrition research - metabolic effects of dietary isoflavones in humans. J Nutr Biochem. 2005;16(4):236–44.
Bertini I, Calabro A, De Carli V, Luchinat C, Nepi S, Porfirio B, Tenori L. The metabonomic signature of celiac disease. J Proteome Res. 2008;8(1):170–7.
Li M, Wang B, Zhang M, Rantalainen M, Wang S, Zhou H, et al. Symbiotc gut microbes modulate human metabolic phenotypes. Proc Natl Acad Sci U S A. 2008;105:2117–22.
Mason S, Reinecke CJ, Kulik W, van Cruchten A, Solomons R, van Furth MT. Cerebrospinal fluid in tuberculous meningitis exhibits only the L-enantiomer of lactic acid. BMC Infect Dis. 2016;16(1):1.
Field AP. Discovering statistics using IBM SPSS statistics: and sex and drugs and rock 'n' roll. 3rd ed. London: Sage Publications; 2009.
Thompson B. Significance, effect sizes, stepwise methods, and other issues: strong arguments move the field. J Exp Educ. 2001;70:80–93.
Cohen J. Statistical power analysis of the behavioural sciences. Revised ed. New York: Academic Press, Inc.; 1988. p. 79–81.
Quinones MP, Kaddurah-Daouk R. Metabolomics tools for identifying biomarkers for neuropsychiatric diseases. Neurobiol Dis. 2009;35(2):165–76.
Sinclair AJ, Viant MR, Ball AK, Burdon MA, Walker EA, Stewart PM, et al. NMR-based metabolomic analysis of cerebrospinal fluid and serum in neurological diseases- a diagnostic tool? NMR Biomed. 2010;23(2):123–32.
Mäntyselkä P, Miettola J, Niskanen L, Kumpusalo E. Glucose regulation and chronic pain at multiple sites. Rheumatology. 2008;47(8):1235–8.
Clish CB. Metabolomics: an emerging but powerful tool for precision medicine. Mol Case Stud. 2015;1(1):a000588.
Mayer EA, Tillisch K. The brain-gut axis in abdominal pain syndromes. Ann Rev Med. 2011;62
Cryan JF, O’Mahony SM. The microbiome-gut-brain axis: from bowel to behavior. Neurogastroenterol Motil. 2011;23(3):187–92.
Research funding for this project was provided by the Technological Innovation Agency (TIA) of the South African Department of Science and Technology (DST) and from the Nuclear Technologies in Medicine and the Biosciences Initiative (NTeMBI) of the Nuclear Energy Corporation of South Africa (NECSA). BM received a post-graduate bursary from the National Research Foundation (NRF) of South Africa and a SAVUSA-SKILL stipend from the Embassy of the Netherlands in Pretoria, enabling her to pursue the NMR analyses at the Radboud University Medical Centre in Nijmegen.
Availability of data and materials
See “Additional files” under Supplementary Information.
This investigation required a multidisciplinary approach and the inputs of all authors were essential to produce the concept and final manuscript. CR and HM defined the aim of the study, developed the experimental design and obtained ethical approval for the study; they are, respectively, the promoter and co-promoter for BM. HM performed the selection and clinical description of the FMS patients and CF and CO control groups and provided the urine samples for these individuals. CR provided it for the CN control group and was responsible for coordination and integration of inputs from the authors. BM conducted all experimental analyses and compiled the clinical information provided by HM. SM and UE gave guidance and assessment on 1H NMR data generation and spectral analyses. RW was responsible for critical evaluation of the analytical aspects of the clinical chemistry data and for their interpretation. MvR performed all the statistical analyses. All authors have given their approval of the version of the manuscript as submitted and agreed to the accountability requirements.
There are no competing interests.
Consent for publication
Ethics approval and consent to participate
Informed consent was obtained from all the participants in this study by means of a voluntarily completed consent form. Ethical approval for this study was obtained via the consortium under the Nuclear Technologies in Medicine and Biosciences Initiative (NTeMBI) (ethical approval by Pharma Ethics Pty, Ltd., reference number 11064365). Pharma Ethics confirmed the following: “The study has been accepted as complying to the Ethics Standards for Clinical Research with a new drug in participants based on FDA, ICH GCP and the Declaration of Helsinki guidelines. The Ethics Committee (IRB) granting this approval is in compliance with the Guidelines for Good Practice in the Conduct of Clinical Trials in Human Participants in South Africa (2006), ICH Harmonised Tripartite Guidelines E6: Note: for the Guidance in Good Clinical Practice (CPMP/ICH/135/95) and FDA Code of Federal Regulation Part 50, 56 and 312.”
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information (SI) (PDF providing detailed descriptions of methods and additional information – clinical, spectral and statistical – to support the manuscript). (DOCX 1270 kb)
Raw data matrix (Raw NMR spectral data (Excel format) normalized relative to the CH3 singlet of creatinine at 3.13 ppm). (XLSX 298 kb)
About this article
Cite this article
Malatji, B.G., Meyer, H., Mason, S. et al. A diagnostic biomarker profile for fibromyalgia syndrome based on an NMR metabolomics study of selected patients and controls. BMC Neurol 17, 88 (2017). https://doi.org/10.1186/s12883-017-0863-9
- Fibromyalgia syndrome
- Proton nuclear magnetic resonance (1H–NMR) spectroscopy
- Metabolite markers