Vitamin D and fibromyalgia: a meta-analysis.

Vitamin D is a cofactor responsible for autoimmune disorders. There is no agreement in the studies investigating the association between vitamin D and fibromyalgia. This study aims to combine the conflicting results of the primary studies which compared these patients with control groups regarding the serum concentration of vitamin D. This meta-analysis has been designed based on PRISMA guidelines. Relevant keywords were searched in PubMed, Science direct, Scopus, Cochrane, and Google scholar and primary studies were selected. After screening the eligible studies according to inclusion/exclusion criteria, we investigated the risk of bias in the selected studies and also the heterogeneity between the primary results using Cochrane (Q) and I-squared (I2) indices. The primary results were combined using inverse variance method and Cohen statistics as well as a random effects model. Publication bias was assessed using Egger test. Sensitivity analysis was applied to investigate the influence of each primary study on the final result of the meta-analysis. Suspected factors in the heterogeneity were assessed using meta-regression models. We entered 12 eligible studies in the meta-analysis including 851 cases compared with 862 controls. The standardized mean difference of Vitamin D between the two groups was -0.56 (95% confidence interval: -1.05, -0.08). Our meta-analysis showed that vitamin D serum levels of patients with fibromyalgia was significantly lower than that of control group.

INTRODUCTION

Fibromyalgia (FM) is a chronic non-articular rheumatologic disorder specified with chronic systemic musculoskeletal pains with unknown etiology, morning stiffness, fatigue, sleep disorder, multiple tender points, low pain threshold [12], symptoms of depression and anxiety [345] and intestinal dysmotility [456].

FMS affects 4–5 percent of the 30–60-year-old population 85–95 percent of whom are female [13]. The prevalence of FM in women is nine folds higher than that in men [6].

The pathogenesis of FM is still unknown [5]. However, several studies reported the effects of abnormal levels of neurotransmitters and disorders of the endocrine system [678910]. In other words, the interaction between neuroendocrine, metabolic, and immunologic factors in developing FM cannot be ignored. In addition, other factors such as genetics, serotonin, melatonin, substance P, endorphin and vitamin D have been reported to be responsible for the pathophysiology of the disease [35].

Vitamin D is an important regulator of the immune system which acts as an active cofactor in the incidence of autoimmune disorders such as multiple sclerosis, rheumatoid arteritis, insulin dependent diabetes mellitus and irritable bowel syndrome. Vitamin D deficiency is defined as the serum level of 25 hydroxy vitamin D being lower than 20 ng/ml (50 nmol/ml) [11121314]. It seems that more than one billion people suffer from vitamin D deficiency [3]. This problem can lead to proximal muscle weakness, skeletal mineralization insufficiency, risk of falling in the elderly, and widespread body pain [3].

There is no consensus regarding the association between vitamin D and FM. Several studies have mentioned the correlation between a low concentration of vitamin D and non-specific musculoskeletal pains [31213], while, some other studies reported the probable association between vitamin D deficiency and clinical manifestations of FM [134514]. However, other studies, did not find any significant association between FM and vitamin D [1516].

According to the above-mentioned discrepancies between the results of primary studies, we aimed to compare the serum levels of vitamin D between patients with FM and control groups using meta-analysis.

MATERIALS AND METHODS

This study was designed according to the PRISMA guideline [17] as well as a pre-defined protocol, however, this protocol has not been registered.

1. Inclusion criteria

All case-control studies published up to 30 August 2016 in English, carried out among women with fibromyalgia as cases compared with control groups, were entered into the meta-analysis. These studies had to report the mean and standard deviation of blood vitamin D levels in both groups.

2. Search strategy and data banks

To find required primary studies, PubMed, Science direct, Scopus, Cochrane, Ovid and Google scholar were searched from 22 September 2016 to 7 October 2016 by two independent researchers. Any disagreements were settled by a third researcher. The search was conducted without any geographical restriction using the following keywords: Fibromyalgia, FM, Chronic pain syndrome, 25-hydroxyvitamin D, 25OHD, 25-OHD, Vitamin D. Case-control.

3. Study selection

Two experts independently reviewed the studies regarding the inclusion criteria. First, all titles were reviewed and duplicates were excluded. Then, abstracts and full texts were investigated and irrelevant papers were omitted.

Kappa agreement statistics were used during the full text investigation based on the guideline suggested by Landis and Koch [18]. These statistics were categorized into slight agreement (0–0.20), fair agreement (0.21–0.40), moderate agreement (0.41–0.60), substantial agreement (0.61–0.80) and perfect agreement (more than 0.80).

4. Quality assessment

The quality of the primary studies was assessed using the Cochrane risk of bias tool [1920] by two independent researchers and an additional investigator in the case of any conflicts. The risk of bias for each study was evaluated according to selection bias, performance bias, detection bias, attrition bias, reporting bias, and other sources of bias (e.g. bias of study design, extreme baseline imbalance).

Each of these biases were classified as high risk (score 0), low risk (score 2) and unclear risk of bias (score 1). The total risk of bias was calculated by a summation of all categories.

5. Data extraction

Two researchers extracted the required information from the selected studies including first author's name, date of publication, study design, country of the study conduction, sample size in each group, rationale of the sample size estimation, sampling design, diagnostic criteria, mean (SD) age of the participants in each group, mean (SD) of vitamin D in each group and P value for the difference of vitamin D between cases and controls. Discrepancies were investigated by a third researcher and final decisions were made.

6. Data analysis

All analyses were performed using Stata version 14 package (StataCorp, College station, TX, USA) as well as Cochrane Collaboration Review Manager Software (RevMan Version 5.3). Combining the results of primary studies was conducted based on inverse variance and Cohen statistics. The heterogeneity between the results was assessed using Cochrane (Q) and I-square indices considering P < 0.1 as the significance level.

As suggested by Higgins and Thompson [21], I-squared results were categorized into low heterogeneity (less than 26%), moderate heterogeneity (25%–75%) and high heterogeneity (more than 75%). According to the degree of heterogeneity, fixed or random effects models were applied for meta-analysis.

Point and pooled standardized mean differences (95% confidence intervals) of vitamin D levels between cases and controls were illustrated by forest plots. Any presence of publication bias was assessed using the Egger test (P value < 0.1 as significance level).

To assess the effect of each primary study on the total mean difference, sensitivity analysis was conducted. In addition, suspected factors in heterogeneity were investigated using meta-regression models.

RESULTS

Primary search showed 1579 studies (1507 records were identified through database searching and 72 records were identified through non-electronic search of the references) 921 of which were duplicate evidence and were removed. Reviewing the titles and abstracts revealed 644 irrelevant papers. Investigating the full texts of the remaining articles [13514162223242526272829], two irrelevant papers [2227] were also excluded. Others were quality assessed and entered into the final meta-analysis (Fig. 1). The Kappa statistics agreement between the researchers for full texts evaluation based on the inclusion criteria was 0.88 indicating perfect agreement.

Fig. 1

Flowchart of study selection.

All selected studies were case controls carried out among women and published between 2009 and 2016. Of these, 10 studies reported the matching variables (Table 1). The FMS diagnostic criteria were based on the American College of Rheumatology criteria. Mean age of participants in both groups was similar in nine studies. In the other three studies cases were significantly older than the controls.

Table 1

Charactristics of Primary Studies Included to Meta-analysis

NS: Non-significant.

Fig. 2 and 3 represented the risk of bias in the primary studies based on the Cochrane risk of bias tool. There was high risk of selection bias in two studies, unclear risk of performance bias in two other studies, unclear risk of attrition bias in 12 studies, unclear detection bias in one study and unclear reporting bias in one study. Other sources of bias (design and extreme baseline imbalance) were observed in two studies (one study was high risk and in one study the risk was unclear).

Fig. 2

Risk of bias summary: review authors' judgements about each risk of bias item for each included study. ‘+’ denotes low risk of bias, ‘−’ represents high risk of bias and “?” indicated unclear bias.

Fig. 3

Risk of bias graph: review authors' judgements about each risk of bias item presented as percentages across all included studies.

All studies were carried out among 851 cases and 862 controls. In eight studies, the mean level of vitamin D was lower in cases than controls, six of which reported significant results. Of four studies reporting higher levels of vitamin D among cases, only one of the differences was statistically significant. Combining the results of these primary studies using a random effects model (Q = 244.9, P < 0.001; I-squared: 95.5%), the mean difference of vitamin D between patients and controls was estimated as −0.56 (95% confidence interval: −1.05, −0.08).

To assess the effect of the mean age of the participants as a factor for heterogeneity, meta-regression models showed non-significant results (β = 9.9, P = 0.10). We conducted a subgroup analysis based on the difference between the mean age of cases and controls and found different results between studies in which their case/control groups had the same age and those that did not (Fig. 4). In addition, another subgroup analysis was conducted based on the menopausal status of the participants which showed no difference in the pooled estimates between the premenopausal population and total population (Fig. 5).

Fig. 4

Point and pooled standardized mean difference of blood vitamin D level between subjects with and without fibromyalgia according to the age-matching of the cases/controls. In studies whose cases and controls had the same mean age, the pooled estimation of standardized mean difference of vitamin D differed between cases and controls, while in the studies which cases had different mean ages than controls, no considerable difference was observed between the vitamin D concentration of these two groups.

Fig. 5

Point and pooled standardized mean difference of blood levels of vitamin D between cases and controls based on menopausal status of cases/controls. There is no significant difference in the standardized mean difference of vitamin D between studies conducted among populations in premenopausal, postmenopausal and unknown menopausal status.

Sensitivity analysis indicated that exclusion of none of the primary studies influenced the pooled estimates (Fig. 6). Finally, the Egger test showed no evidence of publication bias (β = −17.3, P = 0.23).

Fig. 6

Sensitivity analysis of the effect of each primary study on the pooled results.

DISCUSSION

Review of the primary studies selected for meta-analysis showed that most of them found lower levels of vitamin D among FM patients than the control groups. But, four studies indicated that these patients had higher levels of vitamin D than healthy individuals. However, pooled estimates showed that FM patients in comparison with the control groups had, on average, 0.56 units lower concentration of vitamin D. Most of these studies compared age-matched cases/controls and half of these studies investigated only women in the premenopausal stage.

Vitamin D deficiency is common worldwide especially among children, women and the elderly population [3]. Half of the eight studies indicating lower concentrations of vitamin D among patients with FM had been carried out in Turkey. It seems that in such a country patients with FMS had low sunlight exposure due to their limited functional capacity.

Olama et al. [1] investigated the serum levels of vitamin D in the individuals with and without FM and found lower blood vitamin D levels in patients. They also showed that vitamin D level is correlated with pain (according to the VAS score), Beck Depression score, and lumbar bone mineral density.

Okyay et al. [3] showed severe vitamin D deficiency among patients with FMS, which was similar to our results as well as those reported by Olama et al. [1], but was in contrast to the results observed by Okumus et al. [16], which was conducted among premenopausal women.

Olama et al. [1], Okyay et al. [3], Labeeb et al. [4], Baygutalp et al. [5], Altindag et al. [14], and Okumus et al. [16] reported significant correlations between vitamin D blood level and pain intensity. In addition, significant impacts of vitamin D on the quality of life in patients with FM based on the FIQ (Fibromyalgia Impact Questionnaire) were observed by Okyay et al. [3], Labeeb et al. [4], and Altindag et al. [14], which was in contrast to the results of Baygutlap et al. [5].

The largest sample size (205 cases and controls) was used in the study conducted by Mateos et al. [25] while the Baygutalp study [5] with 19 cases/controls included the smallest sample size. Both of these studies reported lower blood vitamin D levels in FM patients than in controls.

One of the limitations of the current meta-analysis was the considerable heterogeneity between the results of the primary studies and we could not detect the source of this heterogeneity by meta-regression models. However, there were differences in the inclusion/exclusion criteria between the primary studies which might be responsible for a part of this heterogeneity. In addition, some of the studies had considered the season of the study conduction but most of them had not taken it into consideration. This might be another explanation for the observed heterogeneity.

In conclusion, our meta-analysis revealed credible evidence that vitamin D can be a determinant factor for fibromyalgia. It seems that the policymakers should consider vitamin D supplementation among women as a preventive strategy. Moreover, further studies are recommended to investigate the role of other factors such as season, obesity, menopausal status, history of hysterectomy, diabetes mellitus and other metabolic disorders, malignancy, pregnancy, sexual hormone levels, and history of OCP consumption in the heterogeneity of the primary studies.