Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Profile of Cerebrospinal microRNAs in Fibromyalgia

  • Jan L. Bjersing ,

    jan.bjersing@rheuma.gu.se

    Affiliation Department of Rheumatology and Inflammation Research, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Göteborg, Sweden

  • Christopher Lundborg,

    Affiliation Department of Anaesthesiology and Intensive Care Medicine, Institute of Clinical Sciences, The Sahlgrenska Academy, University of Gothenburg, Göteborg, Sweden

  • Maria I. Bokarewa,

    Affiliation Department of Rheumatology and Inflammation Research, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Göteborg, Sweden

  • Kaisa Mannerkorpi

    Affiliations Department of Rheumatology and Inflammation Research, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, Göteborg, Sweden, University of Gothenburg Centre for Person-centred Care (GPCC), Sahlgrenska Academy, University of Gothenburg, Göteborg, Sweden

Abstract

Introduction

Fibromyalgia (FM) is characterized by chronic pain and reduced pain threshold. The pathophysiology involves disturbed neuroendocrine function, including impaired function of the growth hormone/insulin-like growth factor-1 axis. Recently, microRNAs have been shown to be important regulatory factors in a number of diseases.

The aim of this study was to try to identify cerebrospinal microRNAs with expression specific for FM and to determine their correlation to pain and fatigue.

Methods

The genome-wide profile of microRNAs in cerebrospinal fluid was assessed in ten women with FM and eight healthy controls using real-time quantitative PCR. Pain thresholds were examined by algometry. Levels of pain (FIQ pain) were rated on a 0-100 mm scale (fibromyalgia impact questionnaire, FIQ). Levels of fatigue (FIQ fatigue) were rated on a 0-100 mm scale using FIQ and by multidimensional fatigue inventory (MFI-20) general fatigue (MFIGF).

Results

Expression levels of nine microRNAs were significantly lower in patients with FM patients compared to healthy controls. The microRNAs identified were miR-21-5p, miR-145-5p, miR-29a-3p, miR-99b-5p, miR-125b-5p, miR-23a-3p, 23b-3p, miR-195-5p, miR-223-3p.

The identified microRNAs with significantly lower expression in FM were assessed with regard to pain and fatigue. miR-145-5p correlated positively with FIQ pain (r=0.709, p=0.022, n=10) and with FIQ fatigue (r=0.687, p=0.028, n=10).

Conclusion

To our knowledge, this is the first study to show a disease-specific pattern of cerebrospinal microRNAs in FM.

We have identified nine microRNAs in cerebrospinal fluid that differed between FM patients and healthy controls. One of the identified microRNAs, miR-145 was associated with the cardinal symptoms of FM, pain and fatigue.

Introduction

Fibromyalgia (FM) is a disorder characterized by chronic pain and reduced pain threshold[1]. The pathophysiology of FM involves disturbed neuroendocrine function, including impaired function of the growth hormone/insulin-like growth factor-1 axis[2-4]. FM is also associated with pronounced fatigue[5-7]. FM is regarded as a multifactorial disorder, possibly involving both environmental and genetic factors.

Clustering of FM in families is common[8,9] and a number of candidate FM genes have been proposed. So far, most of the genetic studies have been fairly small and findings have been partly divergent [10]. However, a large scale candidate study was recently performed that identified three novel candidate genes[11] associated with the disorder.

Since FM is a chronic condition, it is likely that the ongoing disease process and environmental influences leads to long-term changes in gene expression. Evidence is emerging for the involvement of such epigenetic processes in chronic pain[12] but has not been studied specifically in patients with FM.

MicroRNAs have in recent years been identified as important modulators of gene expression in disease processes and physiological pathways. The microRNAs are highly evolutionary conserved, short non-coding RNA molecules, approximately 20-22 nucleotides in length. They inhibit gene expression post-transcriptionally either by inhibition of translation or by degradation of the target messenger RNA. MicroRNAs are believed to regulate at least 30% of human genes[13] and individual microRNAs can repress hundreds of genes[14]. Conversely, one messenger RNA can be the target of several microRNAs. These regulatory networks can control complex programs of gene expression[15] and allows precise adjustments of protein output[14]. Changes in a single or small group of microRNAs may therefore sensitively reflect changes involving a large number of different messenger RNAs.

MicroRNAs are important in the regulation of many processes in the central nervous system (CNS). They have been implicated in several disease states in the CNS[16] and may play important roles in the regulation of neuronal plasticity under stress[17].

To our knowledge expression of cerebrospinal microRNAs in FM has not been previously studied.

The aim of this study was to try to identify cerebrospinal microRNAs with expression specific for FM and also to determine their relation to the cardinal FM symptoms of fatigue and pain.

Materials and Methods

Study design and subjects

Ten patients with FM were compared to eight age-matched healthy controls. Cerebrospinal fluid (CSF) was collected at rest by lumbar puncture through the L3/L4 interspace. Serum was also collected at the same time from eight of the ten FM patients. Collected CSF and serum samples were centrifuged at 800 g for 3 minutes, aliquoted, and stored frozen at -70°C until use. Samples of FM patients were taken during participation in a regular outdoors exercise program, of low to moderate levels, as reported earlier [18]. Patients did not participate in exercise on the day the sample was collected. For characteristics of patients and healthy controls, see Table 1.

Healthy controls (n=8)FM patients (n=10)
Age, years58.5 (38.0 to 68.3)48.5 (37.8 to 53.0)
Symptom duration, years-9 (6-16)
Tender points, n-15 (13-16)
Pain threshold-198.3 (168.0 to 220.6)
FIQ pain-70.3 (53.4 to 83.9)
FIQ fatigue-89 (66.5 to 92.0)
MFIGF-18.5 (13.6 to 20.0)
Antidepressants, yes0 (0%)6 (60%)
Analgesics, yes2(25%)9 (90%)
BMI, kg/m226.5 (22.5 to 28.7)28.4 (25.8 to 29.7)

Table 1. Characteristics and pharmacologic treatment of FM patients and healthy controls.

Median values and interquartile range are indicated
CSV
Download CSV

Healthy controls who had no current pain or history of long-term pain were recruited among patients receiving spinal anaesthesia for gynaecological or urological interventions (n=8). Patients not speaking or reading Swedish, or having any severe somatic or psychiatric disease were excluded.

FM patients were recruited by advertisements in local newspapers. Criteria for inclusion of FM patients: Women with FM, aged 20 to 60 years, and with interest in exercising outdoors twice a week for 15 weeks, willing to participate at blood tests[18]. FM was defined by the ACR 1990 criteria[1]: a history of long-lasting generalized pain and with pain in at least 11 of 18 tender points examined by manual palpation. Criteria for exclusion of FM patients: Patients not speaking or reading Swedish. Having any severe somatic or psychiatric disease. The participation of these FM patients in an outdoors exercise program indicates reasonably good somatic health.

The median age of FM patients was 48.5 (37.8 to 53.0) years and of healthy controls 58.5 (38.0 to 68.3, n.s.) years. The symptom duration of FM patients was 9 (6-16) years and the average number of tender points was 15 (13-16). Ninety percent of FM patients were taking analgesics during the study and 60% were taking antidepressants or sedatives. The median body mass index (BMI) tended to be higher in FM patients 28.4 kg/m2 (25.8 to 29.7) compared to healthy controls 26.5 (22.5 to 28.7 n.s).

Clinical measurements

The pain threshold was examined by using an algometer (Somedic Production AB, Sollentuna, Sweden), measured in kilopascals (kPa)[19]. The pain threshold was measured in two tender point locations in the upper and lower extremities, respectively. The mean value was applied, and a higher value indicates better health. Levels of pain (FIQ pain) were rated on a visual analogue scale (0-100 mm, Fibromyalgia Impact Questionnaire, FIQ)[20]. Levels of fatigue (FIQ fatigue) were also rated on a visual analogue scale (0-100 mm) using FIQ which gives an estimation of global fatigue, as well as with the Multidimensional Fatigue Inventory (MFI-20)[21] subscale of General Fatigue (MFIGF, range 4-20), which estimates fatigue by questions related to feeling “fit”, “tired” and “rested” . Both instruments reflect fatigue during the last week, and a higher score indicates more severe fatigue.

Sample preparation of circulating microRNA

Isolation of total RNA from CSF was conducted at Exiqon Services, Denmark. Total RNA was extracted from serum using the Qiagen miRNeasy Mini Kit. CSF was thawed on ice and centrifuged at 3000 x g for 5 min in a 4°C microcentrifuge. An aliquot of 200 μL of CSF per sample was transferred to a new microcentrifuge tube and 750 μl of a Qiazol mixture containing 1.25 μg/mL of MS2 bacteriophage RNA was added to the CSF. The tube was mixed and incubated for 5 min followed by the addition of 200 μL chloroform. The tube was mixed, incubated for 2 min and centrifuged at 12,000 x g for 15 min in a 4°C microcentrifuge. The upper aqueous phase was transferred to a new microcentrifuge tube and 1.5 volume of 100% ethanol was added. The contents were mixed thoroughly and 750 μL of the sample was transferred to a Qiagen RNeasy Mini spin column in a collection tube followed by centrifugation at 15,000 x g for 30 sec at room temperature. The process was repeated until all remaining sample had been loaded. The Qiagen RNeasy Mini spin column was rinsed with 700 μL Qiagen RWT buffer and centrifuged at 15,000 x g for 1 min at room temperature followed by another rinse with 500 μL Qiagen RPE buffer and centrifuged at 15,000 x g for 1 min at room temperature. A rinse step (500 μL Qiagen RPE buffer) was repeated twice. The Qiagen RNeasy Mini spin column was transferred to a new collection tube and centrifuged at 15,000 x g for 2 min at room temperature. The Qiagen RNeasy Mini spin column was transferred to a new microcentrifuge tube and the lid was left uncapped for 1 min to allow the column to dry. Total RNA was eluted by adding 50 μL of RNase-free water to the membrane of the Qiagen RNeasy mini spin column and incubating for 1 min before centrifugation at 15,000 x g for 1 min at room temperature. The RNA was stored in a -80°C freezer.

microRNA real-time qPCR

Amplification and detection of microRNAs was conducted at Exiqon Services, Denmark. RNA was reverse transcribed in 40μl reactions using the miRCURY LNA Universal RT microRNA PCR, Polyadenylation and cDNA synthesis kit (Exiqon). An artificial RNA (RNA spike-in) was added to the reverse transcription step. This control was used to confirm that the reverse transcription and amplification occurred with equal efficiency in all samples. cDNA was diluted and assayed in PCR reactions according to the protocol for miRCURY LNA Universal RT microRNA PCR; each microRNA was assayed by qPCR on the microRNA Ready-to-Use PCR microchip (Exiqon, Denmark, Cat No 203608) for detection of human microRNAs. Negative controls excluding template from the reverse transcription reaction were included and profiled like the samples. The amplification was performed in a LightCycler 480 Real-Time PCR System (Roche) in 384 well plates. The amplification curves were analyzed using the Roche LC software, both for determination of crossing point (Cp) value, by the 2nd derivative method, and for melting curve analysis.

Data analysis

A total of 742 human microRNAs and three snRNA reference genes were assayed in each sample. The three small nuclear RNA (snRNA) reference genes were not included in the subsequent analysis. MicroRNAs with a Cp<40 in the negative control were excluded from analysis. A detection threshold of Cp<38 was defined to achieve optimal information quality[22]. MicroRNAs with no detectable expression in any samples were omitted from further analysis.

Sixty microRNAs were present in 50% or more of subjects reflecting low levels of microRNAs detectable in the CSF and indicating a CSF-specific pattern of microRNA expression. Serum samples were available for eight of ten FM patients. Therefore, it was possible to determine the presence or absence in serum of the nine microRNAs that differed in CSF between FM and healthy controls. All these nine microRNAs were also present in the serum of these eight FM patients.

Normalisation

Normalization of the microRNA RT-qPCR data was performed with rank normalisation by calculation of fractional rank [23]. First, the Cp value of each individual microRNA species was ranked within each sample. The highest amount of microRNA in one sample was assigned the highest rank in that sample. In the next step, the rank of one microRNA was divided by the total number of detectable microRNAs in the sample. This was the fractional rank.

Thus, for each microRNA species in a sample the fractional rank was calculated using the following formula: Fractional rank = individual rank (rank of one individual microRNA species in the sample)/(number of all expressed microRNAs in sample +1)

Statistical analysis

Fractional ranks of FM patients and healthy controls were compared using with Mann-Whitney U Test. Group comparisons between healthy controls and FM patients were performed for microRNAS expressed with a Cp<38 in 50% or more of subjects. Thus, sixty microRNAs were compared. Bonferroni was used to define the critical p-value for multiple comparisons. The critical p-value is 0.0008333 for 60 comparisons and nine microRNAs were significantly different between FM patients and healthy controls (table 2). Subsequently, relations between these nine identified microRNAs and typical FM symptoms were examined with the Spearman correlation coefficient. Statistical analysis were performed using SPSS (version 19.0.0 for Mac).

Healthy controlsFMComparison of groups
Fractional rankFractional rankp-valuea
miR-21-5p0.974 (0.964 to 0.977)0.715 (0.306 to 0.804)4.6x10-5
miR-145-5p0.956 (0.880 to 0.966)0.080 (0.026 to 0.262)4.6x10-5
miR-29a-3p0.240 (0.133 to 0.563)0.021 (0.017 to 0.037)9.1x10-5
miR-99b-5p0.636 (0.476 to 0.725)0.164 (0.015 to 0.359)1.8x10-4
miR-125b-5p0.981 (0.971 to 0.988)0.912 (0.900 to 0.947)3.2x10-4
miR-23a-3p0.937 (0.883 to 0.956)0.362 (0.215 to 0.654)3.2x10-4
miR-23b-3p0.693 (0.581 to 0.772)0.223 (0.017 to 0.329)5.5x10-4
miR-195-5p0.703 (0.344 to 0.773)0.025 (0.019 to 0.128)5.5x10-4
miR-223-3p0.715 (0.560 to 0.885)0.054 (0.020 to 0.303)5.5x10-4

Table 2. List of microRNAs with significantly different expression between FM patients and healthy controls.

MicroRNAs are listed in order of significance, starting with the lowest p-values. Normalized levels of microRNAs are expressed as the fractional rank. Higher fractional rank represents higher levels of microRNA. Median fractional rank values and interquartile range are indicated. a Mann-Whitney U test.
CSV
Download CSV

Ethics

The study was approved by the ethics committee of Sahlgrenska University Hospital. Written and verbal information was given to all patients, and written consent was obtained from all patients.

Trial registration: ClinicalTrials.gov NCT00643006

Results

Differences in cerebrospinal microRNA expression between FM patients and healthy controls

MicroRNA expression was compared between FM patients and healthy controls. Nine microRNAs in CSF were significantly different in expression between patients with FM and healthy controls after applying Bonferroni correction for multiple testing (Table 2 and Table S1).

Because the controls were older, a separate comparison was made excluding the two oldest patients in the control group. The median age of this subset of healthy controls (n=6) was 47.5 (37.8 to 53) years compared to 48.5 (37.8 to 53.0, n=10, n.s.) years in FM patients. The median BMI of this subset of healthy controls was 27.5 (23.2 to 31.6, n=6) kg/m2 compared to 28.4 kg/m2 (25.8 to 29.7, n=10, n.s.) in patients with FM.

The same nine microRNAs that initially were identified as differently expressed (Table 2) were also the nine highest ranking microRNAs (ranked according to p-value) after matching groups with regard to age, by using the reduced subset of six healthy controls. Four of these microRNAs, miR-21-5p, miR-145-5p, miR-223-3p and miR-29a-3p remained statistically significant after applying Bonferroni correction for multiple testing. The remaining five microRNAS, miR-99b-5p (p<0.001), miR-125b-5p (p<0.002), miR-23a-3p (p<0.002), miR-23b-3p (p<0.003) and miR-195-5p (p<0.003) did not reach statistical significance.

Disease associated microRNAs and symptom severity

The relation of the nine identified microRNAs to pain and fatigue was subsequently tested in the FM patients. The correlations are shown in Table 3. The identified FM-associated microRNAs in cerebrospinal fluid were assessed with regard to pain and fatigue. miR-145-5p was correlated with FIQ pain (r=0.709, p=0.022, n=10) and with FIQ fatigue (r=0.687, p=0.028, n=10).

Pain thresholdFIQ painFIQ fatigueMFIGF
miR-21-5prs=-0.382rs=0.018rs=-0.103rs=0.138
p=0.276p=0.960p=0.776p=0.704
n=10n=10n=10n=10
miR-145-5prs=-0.115rs=0.709rs=0.687rs=0.276
p=0.751p=0.022p=0.028p=0.440
n=10n=10n=10n=10
miR-29a-3prs=-0.298rs=-0.055rs=-0.192rs=0.277
p=0.403p=0.881p=0.595p=0.439
n=10n=10n=10n=10
miR-99b-5prs=-0.553rs=-0.237rs=0.052rs=0.145
p=0.097p=0.510p=0.887p=0.690
n=10n=10n=10n=10
miR-125b-5prs=0.103rs=0.201rs=0.174rs=0.258
p=0.776p=0.578p=0.631p=0.472
n=10n=10n=10n=10
miR-23a-3prs=0.018rs=-0.103rs=0.036rs=0.088
p=0.960p=0.777p=0.920p=0.809
n=10n=10n=10n=10
miR-23b-3prs=-0.164rs=-0.285rs=0.055rs=0.32
p=0.651p=0.425p=0.881p=0.367
n=10n=10n=10n=10
miR-195-5prs=-0.370rs=-0.018rs=0.395rs=0.464
p=0.293p=0.960p=0.258p=0.177
n=10n=10n=10n=10
miR-223-3prs=-0.286rs=-0.310rs=0.137rs=0.447
p=0.424p=0.383p=0.705p=0.196
n=10n=10n=10n=10

Table 3. Relation of disease-associated CSF microRNAs with cardinal symptoms in FM.

Spearmans rho (rs) of miRNA relative expression levels vs. symptom severity levels , two-tailed significance (p , and number of subjects (n) are presented.
CSV
Download CSV

Discussion

The pathophysiology of chronic pain in FM involves disturbed neuroendocrine function with contribution of the growth hormone axis. The pathogenesis of fatigue in FM is less well understood. MicroRNAs are increasingly recognised as important regulatory factors with diagnostic and therapeutic potential in many human diseases.

In this study we identified several microRNAs that were differently expressed in the CSF of FM patients and the healthy controls. These difference may reflect pathological changes in the central nervous system. Changes in relative microRNA levels may also represent adaptive responses to minimise the effects of the disease and its symptoms.

miR-21 is involved in glial[24,25] and neuronal[26,27] responses to injury in the central nervous system[28]. Identifying alterations in cerebrospinal microRNA networks in FM may give new leads to the pathogenesis of this condition.

miR-145-5p correlated with pain and fatigue in the FM patients. miR-145 inhibits growth[29] and is a potential modulator of the insulin-like growth factor pathway[30]. miR-145 was downregulated in endometriosis patients[31], in ulcerative colitis patients[32] and in experimental studies of prenatal stress[33]. However, in this small sample of patients with FM, post-hoc analysis did not show any significant correlation between miR-145-5p and serum-free or total insulin-like growth factor 1 (data not shown).

miR-29 is upregulated during aging in mice[34] and its inhibition reduced cellular aging[35]. miR-29 is decreased in Alzheimers disease[36], Parkinsons disease[37] and myotonic dystrophy type 1[38]. Together with miR-125, miR-99, miR-223, and miR-145, miR-29 was upregulated during endometriosis[39]. In an exercise study in animals miR-29 was increased leading to improved cardiac function after exercise training[40].

miR-23b is a regulator of µ-opioid receptor expression and responds to long-term morphine treatment[41]. MiR-23b is repressed in several autoimmune conditions including multiple sclerosis, systemic lupus erythematosus and experimental autoimmune encephalitis[42]. Conversely, miR-23b is anti-inflammatory and its overexpression suppresses autoimmune pathogenesis under experimental conditions[42].

miR-195-5p was reduced in FM patients compared to healthy controls. miR195 is involved in energy metabolism[43] and growth[44,45] and appears to suppress BDNF protein levels in the prefrontal cortex[46]. MiR-195 protects against hypoperfusion induced dementia by inhibition of beta-amyloid production[47,48].

miR-223 also showed neuroprotective activity by the targeting of glutamate receptors[49] and is changed in the mouse prefrontal cortex after inflammatory pain[50].

All the FM-associated microRNAs were significantly reduced compared to healthy controls. Among FM patients, miR-145-5p was correlated with pain and fatigue.

A limitation of this study is the small number of subjects. Therefore, these results are indicative and validation in other groups of patients and with larger study groups would be valuable.

The relative amount of miR-145-5p was strongly downregulated in FM patients. Nevertheless, among FM patients higher levels of pain and fatigue corresponded to higher amounts of miR-145-5p. An explanation is that higher levels of miR-145-5p in FM patients represent an adaptive but ineffective response to pain and fatigue. This hypothesis should be investigated in longitudinal studies with interventions that change pain and fatigue in FM.

Conclusion

This study shows for the first time a disease-specific pattern of cerebrospinal microRNAs in FM patients. One of the identified microRNAs, miR-145 was associated with the cardinal symptoms of FM, pain and fatigue.

Supporting Information

Table S1.

Complete list of microRNAs included in the comparative analysis between FM patients and healthy controls. All microRNAs included in the group comparison are listed in order of significance, starting with the lowest p-values. Normalized levels of microRNAs are expressed as the fractional rank. Higher fractional rank represents higher levels of microRNA. Median fractional rank values and interquartile range are indicated. aMann-Whitney U test.

https://doi.org/10.1371/journal.pone.0078762.s001

(DOC)

Acknowledgments

We thank Lena Nordeman, Åsa Cider, Gunilla Jonsson, Annie Palstam and members of the GAU-study group for recruiting, examining or supervising the exercise groups.

Author Contributions

Conceived and designed the experiments: JB KM. Performed the experiments: JB. Analyzed the data: JB KM MB CL. Contributed reagents/materials/analysis tools: KM MB CL. Wrote the manuscript: JB KM MB CL.

References

  1. 1. Wolfe F, Smythe HA, Yunus MB, Bennett RM, Bombardier C et al. (1990) The American College of Rheumatology 1990 Criteria for the Classification of Fibromyalgia. Report of the Multicenter Criteria Committee. Arthritis Rheum 33: 160-172. doi:https://doi.org/10.1002/art.1780330203. PubMed: 2306288.
  2. 2. Bennett RM (2002) Adult growth hormone deficiency in patients with fibromyalgia. Curr Rheumatol Rep 4: 306-312. doi:https://doi.org/10.1007/s11926-002-0039-4. PubMed: 12126582.
  3. 3. Bennett RM, Cook DM, Clark SR, Burckhardt CS, Campbell SM (1997) Hypothalamic-pituitary-insulin-like growth factor-I axis dysfunction in patients with fibromyalgia. J Rheumatol 24: 1384-1389. PubMed: 9228141.
  4. 4. Bjersing JL, Dehlin M, Erlandsson M, Bokarewa MI, Mannerkorpi K (2012) Changes in pain and insulin-like growth factor 1 in fibromyalgia during exercise: the involvement of cerebrospinal inflammatory factors and neuropeptides. Arthritis Res Ther 14: R162. doi:https://doi.org/10.1186/ar3902. PubMed: 22776095.
  5. 5. Aaron LA, Burke MM, Buchwald D (2000) Overlapping conditions among patients with chronic fatigue syndrome, fibromyalgia, and temporomandibular disorder. Arch Intern Med 160: 221-227. doi:https://doi.org/10.1001/archinte.160.2.221. PubMed: 10647761.
  6. 6. Mease PJ, Clauw DJ, Christensen R, Crofford LJ, Gendreau RM et al. (2011) Toward development of a fibromyalgia responder index and disease activity score: OMERACT module update. J Rheumatol 38: 1487-1495. doi:https://doi.org/10.3899/jrheum.110277. PubMed: 21724721.
  7. 7. Staud R (2012) Peripheral and central mechanisms of fatigue in inflammatory and noninflammatory rheumatic diseases. Curr Rheumatol Rep 14: 539-548. doi:https://doi.org/10.1007/s11926-012-0277-z. PubMed: 22802155.
  8. 8. Arnold LM, Hudson JI, Hess EV, Ware AE, Fritz DA et al. (2004) Family study of fibromyalgia. Arthritis Rheum 50: 944-952. doi:https://doi.org/10.1002/art.20042. PubMed: 15022338.
  9. 9. Markkula R, Järvinen P, Leino-Arjas P, Koskenvuo M, Kalso E et al. (2009) Clustering of symptoms associated with fibromyalgia in a Finnish Twin Cohort. Eur J Pain 13(7): 744-750. doi:https://doi.org/10.1016/j.ejpain.2008.09.007. PubMed: 18938094.
  10. 10. Lee YH, Choi SJ, Ji JD, Song GG (2012) Candidate gene studies of fibromyalgia: a systematic review and meta-analysis. Rheumatol Int 32(2): 417-426. doi:https://doi.org/10.1007/s00296-010-1678-9. PubMed: 21120487.
  11. 11. Smith SB, Maixner DW, Fillingim RB, Slade G, Gracely RH et al. (2012) Large candidate gene association study reveals genetic risk factors and therapeutic targets for fibromyalgia. Arthritis Rheum 64: 584-593. doi:https://doi.org/10.1002/art.33338. PubMed: 21905019.
  12. 12. Denk F, McMahon SB (2012) Chronic pain: emerging evidence for the involvement of epigenetics. Neuron 73: 435-444. doi:https://doi.org/10.1016/j.neuron.2012.01.012. PubMed: 22325197.
  13. 13. Lewis BP, Burge CB, Bartel DP (2005) Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120: 15-20. doi:https://doi.org/10.1016/j.cell.2004.12.035. PubMed: 15652477.
  14. 14. Baek D, Villén J, Shin C, Camargo FD, Gygi SP et al. (2008) The impact of microRNAs on protein output. Nature 455: 64-71. doi:https://doi.org/10.1038/nature07242. PubMed: 18668037.
  15. 15. Zhu H, Fan GC (2011) Extracellular/circulating microRNAs and their potential role in cardiovascular disease. Am J Cardiovascular Dis 1: 138-149. PubMed: 22059153.
  16. 16. Meza-Sosa KF, Valle-García D, Pedraza-Alva G, Pérez-Martínez L (2012) Role of microRNAs in central nervous system development and pathology. J Neurosci Res 90: 1-12. doi:https://doi.org/10.1002/jnr.22701. PubMed: 21922512.
  17. 17. Schouten M, Aschrafi A, Bielefeld P, Doxakis E, Fitzsimons CP (2013) microRNAs and the regulation of neuronal plasticity under stress conditions. Neuroscience 241C: 188-205. PubMed: 23500144.
  18. 18. Mannerkorpi K, Nordeman L, Cider A, Jonsson G (2010) Does moderate-to-high intensity Nordic walking improve functional capacity and pain in fibromyalgia? A prospective randomized controlled trial. Arthritis Res Ther 12: R189. doi:https://doi.org/10.1186/ar3159. PubMed: 20942911.
  19. 19. Kosek E, Ekholm J, Nordemar R (1993) A comparison of pressure pain thresholds in different tissues and body regions. Long-term reliability of pressure algometry in healthy volunteers. Scand J Rehabil Med 25: 117-124. PubMed: 8248762.
  20. 20. Burckhardt CS, Clark SR, Bennett RM (1991) The fibromyalgia impact questionnaire: development and validation. J Rheumatol 18: 728-733. PubMed: 1865419.
  21. 21. Ericsson A, Mannerkorpi K (2007) Assessment of fatigue in patients with fibromyalgia and chronic widespread pain. Reliability and validity of the Swedish version of the MFI-20. Disabil Rehabil 29: 1665-1670. doi:https://doi.org/10.1080/09638280601055782. PubMed: 17852297.
  22. 22. Jensen SG, Lamy P, Rasmussen MH, Ostenfeld MS, Dyrskjøt L et al. (2011) Evaluation of two commercial global miRNA expression profiling platforms for detection of less abundant miRNAs. BMC Genomics 12: 435. doi:https://doi.org/10.1186/1471-2164-12-435. PubMed: 21867561.
  23. 23. Qiu X, Wu H, Hu R (2013) The impact of quantile and rank normalization procedures on the testing power of gene differential expression analysis. BMC Bioinformatics 14: 124. doi:https://doi.org/10.1186/1471-2105-14-124. PubMed: 23578321.
  24. 24. Bhalala OG, Pan L, Sahni V, McGuire TL, Gruner K et al. (2012) microRNA-21 regulates astrocytic response following spinal cord injury. J Neurosci Off J Soc Neurosci 32(50): 17935-17947. doi:https://doi.org/10.1523/JNEUROSCI.3860-12.2012. PubMed: 23238710.
  25. 25. Zhang L, Dong LY, Li YJ, Hong Z, Wei WS (2012) miR-21 represses FasL in microglia and protects against microglia-mediated neuronal cell death following hypoxia/ischemia. Glia 60: 1888-1895. doi:https://doi.org/10.1002/glia.22404. PubMed: 22907769.
  26. 26. Strickland IT, Richards L, Holmes FE, Wynick D, Uney JB et al. (2011) Axotomy-induced miR-21 promotes axon growth in adult dorsal root ganglion neurons. PLOS ONE 6: e23423. doi:https://doi.org/10.1371/journal.pone.0023423. PubMed: 21853131.
  27. 27. Buller B, Liu X, Wang X, Zhang RL, Zhang L et al. (2010) MicroRNA-21 protects neurons from ischemic death. FEBS J 277: 4299-4307. doi:https://doi.org/10.1111/j.1742-4658.2010.07818.x. PubMed: 20840605.
  28. 28. Shi Y, Zhang X, Tang X, Wang P, Wang H et al. (2012) MiR-21 is continually elevated long-term in the brain after exposure to ionizing radiation. Radiat Res 177: 124-128. doi:https://doi.org/10.1667/RR2764.1. PubMed: 22034847.
  29. 29. Sachdeva M, Mo YY (2010) miR-145-mediated suppression of cell growth, invasion and metastasis. Am J Transl Res 2: 170-180. PubMed: 20407606.
  30. 30. Law PT, Ching AK, Chan AW, Wong QW, Wong CK et al. (2012) MiR-145 modulates multiple components of the insulin-like growth factor pathway in hepatocellular carcinoma. Carcinogenesis 33: 1134-1141. doi:https://doi.org/10.1093/carcin/bgs130. PubMed: 22431718.
  31. 31. Wang WT, Zhao YN, Han BW, Hong SJ, Chen YQ (2013) Circulating microRNAs identified in a genome-wide serum microRNA expression analysis as noninvasive biomarkers for endometriosis. J Clin Endocrinol Metab 98: 281-289. doi:https://doi.org/10.1210/jc.2012-2415. PubMed: 23118427.
  32. 32. Pekow JR, Dougherty U, Mustafi R, Zhu H, Kocherginsky M et al. (2012) miR-143 and miR-145 are downregulated in ulcerative colitis: putative regulators of inflammation and protooncogenes. Inflamm Bowel Dis 18: 94-100. doi:https://doi.org/10.1097/00054725-201212001-00232. PubMed: 21557394.
  33. 33. Zucchi FC, Yao Y, Ward ID, Ilnytskyy Y, Olson DM et al. (2013) Maternal stress induces epigenetic signatures of psychiatric and neurological diseases in the offspring. PLOS ONE 8: e56967. doi:https://doi.org/10.1371/journal.pone.0056967. PubMed: 23451123.
  34. 34. Ugalde AP, Ramsay AJ, de la Rosa J, Varela I, Mariño G et al. (2011) Aging and chronic DNA damage response activate a regulatory pathway involving miR-29 and p53. EMBO J 30: 2219-2232. doi:https://doi.org/10.1038/emboj.2011.124. PubMed: 21522133.
  35. 35. Martinez I, Cazalla D, Almstead LL, Steitz JA, DiMaio D (2011) miR-29 and miR-30 regulate B-Myb expression during cellular senescence. Proc Natl Acad Sci U S A 108: 522-527. doi:https://doi.org/10.1073/pnas.1017346108. PubMed: 21187425.
  36. 36. Shioya M, Obayashi S, Tabunoki H, Arima K, Saito Y et al. (2010) Aberrant microRNA expression in the brains of neurodegenerative diseases: miR-29a decreased in Alzheimer disease brains targets neurone navigator 3. Neuropathol Appl Neurobiol 36: 320-330. doi:https://doi.org/10.1111/j.1365-2990.2010.01076.x. PubMed: 20202123.
  37. 37. Margis R, Margis R, Rieder CR (2011) Identification of blood microRNAs associated to Parkinsonis disease. J Biotechnol 152: 96-101. doi:https://doi.org/10.1016/j.jbiotec.2011.01.023. PubMed: 21295623.
  38. 38. Perbellini R, Greco S, Sarra-Ferraris G, Cardani R, Capogrossi MC et al. (2011) Dysregulation and cellular mislocalization of specific miRNAs in myotonic dystrophy type 1. Neuromuscul Disord NMD 21: 81-88. doi:https://doi.org/10.1016/j.nmd.2010.11.012. PubMed: 21169019.
  39. 39. Ohlsson Teague EM, Van der Hoek KH, Van der Hoek MB, Perry N, Wagaarachchi P et al. (2009) MicroRNA-regulated pathways associated with endometriosis. Mol Endocrinol 23: 265-275. PubMed: 19074548.
  40. 40. Soci UP, Fernandes T, Hashimoto NY, Mota GF, Amadeu MA et al. (2011) MicroRNAs 29 are involved in the improvement of ventricular compliance promoted by aerobic exercise training in rats. Physiol Genomics 43: 665-673. doi:https://doi.org/10.1152/physiolgenomics.00145.2010. PubMed: 21447748.
  41. 41. Rodriguez RE (2012) Morphine and microRNA Activity: Is There a Relation with Addiction? Front Genet 3: 223. PubMed: 23162566.
  42. 42. Zhu S, Pan W, Song X, Liu Y, Shao X et al. (2012) The microRNA miR-23b suppresses IL-17-associated autoimmune inflammation by targeting TAB2, TAB3 and IKK-alpha. Nature medicine 18:1077-1086.
  43. 43. Fei X, Qi M, Wu B, Song Y, Wang Y et al. (2012) MicroRNA-195-5p suppresses glucose uptake and proliferation of human bladder cancer T24 cells by regulating GLUT3 expression. FEBS Lett 586: 392-397. doi:https://doi.org/10.1016/j.febslet.2012.01.006. PubMed: 22265971.
  44. 44. van Rooij E, Sutherland LB, Liu N, Williams AH, McAnally J et al. (2006) A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc Natl Acad Sci U S A 103: 18255-18260. doi:https://doi.org/10.1073/pnas.0608791103. PubMed: 17108080.
  45. 45. Hui W, Yuntao L, Lun L, WenSheng L, ChaoFeng L et al. (2013) MicroRNA-195 inhibits the proliferation of human glioma cells by directly targeting cyclin D1 and cyclin E1. PLOS ONE 8: e54932. doi:https://doi.org/10.1371/journal.pone.0054932. PubMed: 23383003.
  46. 46. Mellios N, Huang HS, Grigorenko A, Rogaev E, Akbarian S (2008) A set of differentially expressed miRNAs, including miR-30a-5p, act as post-transcriptional inhibitors of BDNF in prefrontal cortex. Hum Mol Genet, 17(19): 3030-3042. doi:https://doi.org/10.1093/hmg/ddn201. PubMed: 18632683.
  47. 47. Ai J, Sun LH, Che H, Zhang R, Zhang TZ et al. (2013) MicroRNA-195 protects against dementia induced by chronic brain hypoperfusion via its anti-amyloidogenic effect in rats. J Neurosci Off J Soc Neurosci 33: 3989-4001. doi:https://doi.org/10.1523/JNEUROSCI.1997-12.2013. PubMed: 23447608.
  48. 48. Zhu HC, Wang LM, Wang M, Song B, Tan S et al. (2012) MicroRNA-195 downregulates Alzheimer's disease amyloid-beta production by targeting BACE1. Brain research bulletin 88:596-601.
  49. 49. Harraz MM, Eacker SM, Wang X, Dawson TM, Dawson VL (2012) MicroRNA-223 is neuroprotective by targeting glutamate receptors. Proc Natl Acad Sci U S A 109: 18962-18967. doi:https://doi.org/10.1073/pnas.1121288109. PubMed: 23112146.
  50. 50. Poh KW, Yeo JF, Ong WY (2011) MicroRNA changes in the mouse prefrontal cortex after inflammatory pain. Eur J Pain 15:801: e801-e812. PubMed: 21397537.