- Open Access
Identification of candidate genes associated with fibromyalgia susceptibility in southern Spanish women: the al-Ándalus project
Journal of Translational Medicine volume 16, Article number: 43 (2018)
Candidate-gene studies on fibromyalgia susceptibility often include a small number of single nucleotide polymorphisms (SNPs), which is a limitation. Moreover, there is a paucity of evidence in Europe. Therefore, we compared genotype frequencies of candidate SNPs in a well-characterised sample of Spanish women with fibromyalgia and healthy non-fibromyalgia women.
A total of 314 women with a diagnosis of fibromyalgia (cases) and 112 non-fibromyalgia healthy (controls) women participated in this candidate-gene study. Buccal swabs were collected for DNA extraction. Using TaqMan™ OpenArray™, we analysed 61 SNPs of 33 genes related to fibromyalgia susceptibility, symptoms, or potential mechanisms.
We observed that the rs841 and rs1799971 GG genotype was more frequently observed in fibromyalgia than in controls (p = 0.04 and p = 0.02, respectively). The rs2097903 AT/TT genotypes were also more often present in the fibromyalgia participants than in their control peers (p = 0.04). There were no differences for the remaining SNPs.
We identified, for the first time, associations of the rs841 (guanosine triphosphate cyclohydrolase 1 gene) and rs2097903 (catechol-O-methyltransferase gene) SNPs with higher risk of fibromyalgia susceptibility. We also confirmed that the rs1799971 SNP (opioid receptor μ1 gene) might confer genetic risk of fibromyalgia. We did not adjust for multiple comparisons, which would be too stringent and yield to non-significant differences in the genotype frequencies between cases and controls. Our findings may be biologically meaningful and informative, and should be further investigated in other populations. Of particular interest is to replicate the present study in a larger independent sample to confirm or refute our findings. On the other hand, by including 61 SNPs of 33 candidate-genes with a strong rationale (they were previously investigated in relation to fibromyalgia susceptibility, symptoms or potential mechanisms), the present research is the most comprehensive candidate-gene study on fibromyalgia susceptibility to date.
Fibromyalgia is a disease in which people experience chronic diffuse musculoskeletal pain that is usually accompanied by other symptoms such as fatigue, unrefreshed sleep, and cognitive problems [1,2,3,4]. This disease is nine times more common in women than in men . Family aggregation suggests but not confirm genetic susceptibility to fibromyalgia . Variations in neurotransmitter-related genes increase susceptibility to fibromyalgia via hypersensitivity to peripheral painful stimulus by the central nervous system (CNS) . Catechol-O-methyltransferase (COMT), the most widely studied gene in fibromyalgia, is involved in degrading catecholamines and several other neurotransmitters and, therefore, in modulating pain perception by the CNS. The association between COMT single nucleotide polymorphisms (SNPs) and fibromyalgia susceptibility is controversial [7, 8].
Guanosine triphosphate cyclohydrolase 1 (GCH1) and opioid receptor μ1 (OPRM1) are candidate neurotransmitter-related genes that may confer fibromyalgia susceptibility [9, 10]. GCH1 gene participates in the synthesis of dopamine and serotonin. In Korean individuals, the rs841 SNP was associated with discomfort with a tender point examination . OPRM1 gene encodes μ-opioid receptor that bids opiates. In the Turkish population, a study identified that rs1799971 SNP is associated with fibromyalgia susceptibility . However, whether GCH1 and OPRM1 genes are associated with fibromyalgia susceptibility in the Caucasian population is unknown.
Candidate-gene studies on fibromyalgia susceptibility often include a small number of SNPs (e.g., [9,10,11]), which is a limitation [12, 13]. Additionally, there is a paucity of evidence in Europe. The present candidate-gene study compared genotype frequencies of candidate SNPs in a well-characterised sample of Spanish women with fibromyalgia (cases) vs. healthy non-fibromyalgia women (controls). We, therefore, analysed SNPs that had been previously investigated in relation to fibromyalgia susceptibility, symptoms, or potential mechanisms.
The participants were recruited mainly via fibromyalgia associations from Andalusia (southern Spain). The fibromyalgia patients invited to non-fibromyalgia acquaintances with similar sociodemographic characteristics to participate in the study as controls. All participants signed an informed consent form. The Ethics Committee of the Virgen de las Nieves Hospital (Granada, Spain) approved the study. We followed the ethical guidelines of the Declaration of Helsinki.
The fibromyalgia participants had been previously diagnosed with fibromyalgia by a professional rheumatologist and met the 1990 American College of Rheumatology (ACR) criteria for fibromyalgia, which was further confirmed by a tender point examination . Controls neither had a medical diagnosis of fibromyalgia nor fulfilled the 1990 ACR criteria.
The participants were genotyped for 61 SNPs (online Additional file 1: Table S1) that had been previously investigated in relation to fibromyalgia susceptibility, symptoms, or potential mechanisms (online Additional file 1: Table S2). As described elsewhere [14, 15], we collected buccal mucosa cells and we performed DNA non-organic extraction (proteinase K and salting-out) procedures. Afterwards, we conducted a spectrophotometric quantification (NanoDrop 2000c, ThermoFisher). All the samples were standardised to 50 ng/μL and they had an A260/230 ratio between 1.7 and 1.9. The samples in low concentration values were processed by Genomiphi™ V2 DNA Amplification Kit (Sigma Aldrich). Until being processed, all the samples were stored at − 20 °C. The OpenArray® genotyping was performed according to the manufacturer’s protocol using Accufill™ system liquid robot that adds 1.2 μl of genomic DNA sample and 3.8 µl of the following reagents: 2× TaqMan® Universal PCR Master Mix, No AmpErase® UNG, 20× Primer and TaqMan® Probe (FAM™ dye) mix, and sterile-filtered water. The online Additional file 1: Table S3 shows the manufacturer thermal cycling conditions.
Plates include a NTC for each SNP in the analysis, and each plate has a total of 48 samples. Online Additional file 1: Tables S1 and S4 provide further details about the TaqMan® OpenArray® custom assay designs and the 61 analysed SNPs, respectively. We performed a TaqMan™ OpenArray™ Genotyping Plate, Custom Format 64 QuantStudio™ 12 K Flex. The data were analysed using the TaqMan® Genotyper Software and downstream analysis using the AutoCaller™ Software.
The Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium (LD = r2 > 0.5) were checked for all the SNPs using SNPStats v3.0.1 . Using SPSS for Mac v.20.0 (IBM, Armonk, NY, USA), we computed Pearson’s χ2 and logistic regression to analyse the differences between fibromyalgia and controls on SNPs genotype frequency. Significance was set at p < 0.05.
Online Additional file 1: Table S5 shows the characteristics of the participants included in the present study. Those SNPs that did not fulfil genotyping quality controls were excluded. Through this strict filtering out rule we excluded the following nine SNPs because they did not meet the HWE criterion: the rs6323, rs7911, rs136078, rs806377, rs1050450, rs3746544, rs4411417, rs7124442, and rs12620053. HWE and LD were confirmed for the remaining 52 SNPs, which were, therefore, included in the present study.
The rs841and rs1799971 GG genotypes were more frequently observed in fibromyalgia than in controls [odds ratio (OR) (95% confidence interval, 95% CI) = 1.66 (1.03–2.67), p = 0.04 and 1.72 (1.11–2.67), p = 0.02, respectively; Table 1]. The rs2097903 AT/TT genotypes were also present in the fibromyalgia participants more often than in their non-fibromyalgia peers [OR (95% CI) = 1.66 (1.02–2.71), p = 0.04; Table 1). There were no differences for the remaining 49 SNPs (online Additional file 1: Table S6).
The present candidate-gene study (including 61 SNPs of 33 genes) is the largest conducted on fibromyalgia susceptibility. We identified, for the first time, that the rs841 (GCH1 gene) and rs2097903 (COMT gene) SNPs were associated with having fibromyalgia. We also confirmed that the rs1799971 SNP (OPRM1 gene) seems to confer genetic risk of fibromyalgia. The COMT gene is involved in neurotransmitters degradation while the GCH1 gene is related to their synthesis. The μ-opioid receptor, which is encoded by the OPRM1 gene, bids both endogenous and exogenous opiates. In line with previous literature , our findings suggest that an augmented processing of pain may be involved in the genetic susceptibility to fibromyalgia.
Vargas-Alarcón et al.  showed that genotype frequencies of the rs4680, rs4818, and rs6269 SNPs, but not rs2097903, were different between fibromyalgia and control participants. In contrast, we found that the rs2097903 SNP was different between the study cohorts; 78.8% of the patients, and only 69.1% of the controls were TT genotype carriers. Vargas-Alarcón et al.  excluded fibromyalgia participants with rheumatic comorbidities and control participants with chronic pain. Of note is that we corroborated the rheumatologist diagnosis of our fibromyalgia participants by fulfilling the 1990 ACR criteria and discarding it in controls. We included fibromyalgia participants with other rheumatic conditions and controls that experienced chronic pain if they did not fulfil the fibromyalgia diagnosis. Therefore, the ecological validity of our findings is high. Moreover, a later genome-wide association study  in the Spanish population also failed to replicate the findings by Vargas-Alarcón et al. .
GCH enzyme participates in nitric oxide (NO) production. It is widely recognised that NO plays a key role in health and disease by a substantial number of paths . Among other roles, increased NO concentrations often lead to dorsal horn hyperexcitability . Overall, our results are consistent with a previous study in a Korean population in which genotype frequencies of rs841, rs752688, rs4411417, and rs3783641 SNPs in the GCH1 gene were not associated with fibromyalgia susceptibility . However, we found a significant association of the rs841 SNP and fibromyalgia susceptibility with 76.8% of the cases carrying the GG genotype vs. 66.7% of the controls. Interestingly, rs841 was the only SNP correlated to discomfort with a tender point examination in Korean patients . Currently, no study in Caucasian population is available; therefore, future research testing associations of GCH1 gene SNPs, particularly the rs841, with fibromyalgia susceptibility in Caucasians is welcome.
A hypersensitive CNS also appears to have behavioural implications. An active lifestyle is associated with better health status in fibromyalgia [19,20,21]. However, impaired functioning of the anterior cingulate cortex and amygdala, among other neural areas, may yield to processing harmless movements as painful . The activity of these structures is modulated by μ-opioid receptor availability, which is reduced in fibromyalgia . In agreement with a study in Turkey , our findings suggest that the rs1799971 SNP confers fibromyalgia susceptibility, as most of the present patients’ sample were GG genotype (55.8%) and most of the controls were AA/AG (57.7%).
Most of the non-significant findings emerged in SNPs of sodium voltage-gated channel alpha subunit 9 (SCN9A) and adrenoceptor alpha 1A (ADRA1A) genes. Genotype frequencies of a SCN9A SNP (rs6754031) were, however, significantly different between fibromyalgia and controls in Mexico . ADRA1A gene activates mitogenic responses and regulates growth and proliferation of many cells. In two subsamples (from Spain and Mexico), a previous study analysed associations of 4 ADRA1A SNPs , showing that only rs1383914 was associated with an increased risk of fibromyalgia susceptibility in the Spanish sample .
We did not include a replication sample; therefore, our findings must be interpreted with caution until they are corroborated or refuted in future studies conducted in independent samples. In order to prevent flares on pain, we collected mucosa cells using buccal swabs (i.e., a both low level of invasiveness and pain method), which on the other hand, precludes conducting a functional study (e.g., gene or protein expression). The unadjusted analyses for multiple comparisons are another limitation, which would yield to non-significant differences in the genotype frequencies of the SNPs. However, to conclude that there is no difference from a statistical point of view would be too stringent. Our findings may be biologically meaningful and informative, and should be further investigated in other populations. Additionally, the unbalanced size of the study cases and controls limits our statistical power. We included, however, more people in the smaller group (cases, n = 112) than most of previous studies .
In conclusion, we identified associations of the rs841 (GCH1 gene) and rs2097903 (COMT gene) SNPs with higher risk of fibromyalgia. We also confirmed that rs1799971 SNP (OPRM1 gene) seem to confer higher susceptibility to fibromyalgia. Further studies are needed to confirm or refute the present findings.
American College of Rheumatology
adrenoceptor alpha 1A
central nervous system
guanosine triphosphate cyclohydrolase 1
opioid receptor μ1
sodium voltage-gated channel alpha subunit 9
single nucleotide polymorphisms
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. Report of the Multicenter Criteria Committee. Arthritis Rheum. 1990; 33:160–72. http://www.ncbi.nlm.nih.gov/pubmed/2306288. Accessed 18 Jan 2015.
Wolfe F, Clauw DJ, Fitzcharles M-A, Goldenberg DL, Häuser W, Katz RS, et al. Fibromyalgia criteria and severity scales for clinical and epidemiological studies: a modification of the ACR Preliminary Diagnostic Criteria for Fibromyalgia. J Rheumatol. 2011;38:1113–22.
Wolfe F, Clauw DJ, Fitzcharles M-A, Goldenberg DL, Häuser W, Katz RL, et al. 2016 Revisions to the 2010/2011 fibromyalgia diagnostic criteria. Semin Arthritis Rheum. 2016;46:319–29.
Segura-Jiménez V, Álvarez-Gallardo IC, Carbonell-Baeza A, Aparicio VA, Ortega FB, Casimiro AJ, et al. Fibromyalgia has a larger impact on physical health than on psychological health, yet both are markedly affected: the al-Ándalus project. Semin Arthritis Rheum. 2015; 44:563–70. http://linkinghub.elsevier.com/retrieve/pii/S0049017214002261. Accessed 3 Aug 2015.
Branco JC, Bannwarth B, Failde I, Abello Carbonell J, Blotman F, Spaeth M, et al. Prevalence of fibromyalgia: a survey in five European countries. Semin Arthritis Rheum. 2010;39:448–53. http://www.ncbi.nlm.nih.gov/pubmed/19250656.
Arnold LM, Hudson JI, Hess E V., Ware AE, Fritz DA, Auchenbach MB, et al. Family study of fibromyalgia. Arthritis Rheum. 2004;50:944–52. http://doi.wiley.com/10.1002/art.20042.
Docampo E, Escaramís G, Gratacòs M, Villatoro S, Puig A, Kogevinas M, et al. Genome-wide analysis of single nucleotide polymorphisms and copy number variants in fibromyalgia suggest a role for the central nervous system. Pain. 2014;155:1102–9. http://dx.doi.org/10.1016/j.pain.2014.02.016.
Vargas-Alarcón G, Fragoso J-M, Cruz-Robles D, Vargas A, Vargas A, Lao-Villadóniga J-I, et al. Catechol-O-methyltransferase gene haplotypes in Mexican and Spanish patients with fibromyalgia. Arthritis Res Ther. 2007;9:R110.
Kim SK, Kim SH, Nah SS, Hyun Lee J, Hong SJ, Kim HS, et al. Association of guanosine triphosphate cyclohydrolase 1 gene polymorphisms with fibromyalgia syndrome in a Korean population. J Rheumatol. 2013;40:316–22.
Solak Ö, Erdoǧan MÖ, Yildiz H, Ulaşli AM, Yaman F, Terzi ESA, et al. Assessment of opioid receptor μ1 gene A118G polymorphism and its association with pain intensity in patients with fibromyalgia. Rheumatol Int. 2014;34(9):1257–61.
Estévez-López F, Aparicio VA, Ruiz JR, Martínez-González LJ, Delgado-Fernández M, Álvarez-Cubero MJ. The TT genotype of the rs6860 polymorphism of the charged multivesicular body protein 1A gene is associated with susceptibility to fibromyalgia in southern Spanish women. Rheumatol Int. 2017;15:1–3. http://www.ncbi.nlm.nih.gov/pubmed/29247346.
Klepstad P. Polymorphism in the μ-opioid receptor gene OPRM1 A118G—an example of the enigma of genetic variability behind chronic pain syndromes. Scand J Pain. 2014;5:8–9. https://doi.org/10.1016/j.sjpain.2013.11.006.
Tabor HK, Risch NJ, Myers RM. Candidate-gene approaches for studying complex genetic traits: practical considerations. Nat Rev Genet. 2002;3:391–7. http://www.ncbi.nlm.nih.gov/pubmed/11988764.
Freeman B, Smith N, Curtis C, Huckett L, Mill J, Craig IW. DNA from buccal swabs recruited by mail: evaluation of storage effects on long-term stability and suitability for multiplex polymerase chain reaction genotyping. Behav Genet. 2003;33:67–72. http://www.ncbi.nlm.nih.gov/pubmed/12645823.
Gómez-Martín A, Hernández AF, Martínez-González LJ, González-Alzaga B, Rodríguez-Barranco M, López-Flores I, et al. Polymorphisms of pesticide-metabolizing genes in children living in intensive farming communities. Chemosphere. 2015;139:534–40. http://www.ncbi.nlm.nih.gov/pubmed/26318115.
Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006; 22:1928–9. http://www.ncbi.nlm.nih.gov/pubmed/16720584.
Ghimire K, Altmann HM, Straub AC, Isenberg JS. Nitric oxide: what’s new to no? Am J Physiol Cell Physiol. 2017;312:C254–62. http://ajpcell.physiology.org/lookup/doi/10.1152/ajpcell.00315.2016.
Kim HY, Lee I, Chun SW, Kim HK. Reactive oxygen species donors increase the responsiveness of dorsal horn neurons and induce mechanical hyperalgesia in rats. Neural Plast. 2015;2015:293423. http://www.ncbi.nlm.nih.gov/pubmed/26457204.
Estévez-López F, Segura-Jiménez V, Álvarez-Gallardo ICIC, Borges-Cosic M, Pulido-Martos M, Carbonell-Baeza A, et al. Adaptation profiles comprising objective and subjective measures in fibromyalgia: the al-Ándalus project. Rheumatology (Oxford). 2017;56:2015–24. http://www.ncbi.nlm.nih.gov/pubmed/28968914.
Segura-Jiménez V, Soriano-Maldonado A, Estévez-López F, Álvarez-Gallardo ICIC, Delgado-Fernández M, Ruiz JRJR, et al. Independent and joint associations of physical activity and fitness with fibromyalgia symptoms and severity: The al-Ándalus project. J Sports Sci. 2017;35:1565–74. http://www.ncbi.nlm.nih.gov/pubmed/27618648. Accessed 30 Nov 2016.
Ellingson LD, Stegner AJ, Schwabacher IJ, Koltyn KF, Cook DB. Exercise strengthens central nervous system modulation of pain in fibromyalgia. Brain Sci. 2016;6:13. http://www.ncbi.nlm.nih.gov/pubmed/26927193. Accessed 3 Mar 2016.
Coghill RC, McHaffie JG, Yen Y-F. Neural correlates of interindividual differences in the subjective experience of pain. Proc Natl Acad Sci USA. 2003;100:8538–42. http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=166264&tool=pmcentrez&rendertype=abstract. Accessed 8 Feb 2016.
Harris RE, Clauw DJ, Scott DJ, McLean SA, Gracely RH, Zubieta J-K. Decreased central μ-opioid receptor availability in fibromyalgia. J Neurosci. 2007;27:10000–6. http://www.jneurosci.org/content/27/37/10000.abstract.
Vargas-Alarcon G, Alvarez-Leon E, Fragoso J-M, Vargas A, Martinez A, Vallejo M, et al. A SCN9A gene-encoded dorsal root ganglia sodium channel polymorphism associated with severe fibromyalgia. BMC Musculoskelet Disord. 2012;13:23. http://www.biomedcentral.com/1471-2474/13/23.
Vargas-Alarcón G, Fragoso JM, Cruz-Robles D, Vargas A, Martinez A, Lao-Villadóniga JI, et al. Association of adrenergic receptor gene polymorphisms with different fibromyalgia syndrome domains. Arthritis Rheum. 2009;60:2169–73.
Lee YHY, Choi SSJ, Ji JJD, Song GGG. Candidate gene studies of fibromyalgia: a systematic review and meta-analysis. Rheumatol Int. 2012;32:417–26. http://link.springer.com/article/10.1007/s00296-010-1678-9. Accessed 12 Feb 2015.
Conception and design of the research: FE-L, VA-A, MD-F, LJM-G, JR-R, and MJA-C; collected the data: FE-L, DC-M, VA-A, VS-J, ICA-G, ASM, MBC, and PAM; analysed the samples in the laboratory: DC-M, LJM-G, and MJA-C; analysed the data: FE-L, RG, MD-F, and JRR; interpreted the results: FE-L, RG, MD-F, LJM-G, JRR, and MJA-C; drafted the manuscript: FE-L, LJM-G, JRR, and MJA-C; critically revised the manuscript: all the authors. All authors read and approved the final manuscript.
The authors gratefully acknowledge all participants for their collaboration and enthusiasm. We would like to thank the collaborators in the al-Ándalus project and all the members of the Physical Activity for Health Promotion (PA-HELP; CTS-1018) research group. We would also like to thank Ms. Gema Garcia (Genomic unit, Genyo) for her helpful support in DNA extraction and Genotyping of the samples. We are grateful to Ms. Carmen Sainz-Quinn for assistance with the English language.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
All participants signed an informed consent form. The Ethics Committee of the Virgen de las Nieves Hospital (Granada, Spain) approved the study. We followed the ethical guidelines of the Declaration of Helsinki.
This work was supported by the Spanish Ministry of Economy and Competitiveness [I+D+i DEP2010-15639, I+D+i DEP2013-40908-R to M.D.-F.; BES-2014-067612 to F.E.-L.]; the Spanish Ministry of Education [FPU2014/02518 to M.B.-C.]; the Consejería de Turismo, Comercio y Deporte, Junta de Andalucía [CTCD-201000019242-TRA to M.D.-F.]; Consejería de Salud, Junta de Andalucía [PI-0520-2016 to M.D.-F.], and the University of Granada, Plan Propio de Investigación 2016, Excellence actions: Units of Excellence; Unit of Excellence on Exercise and Health (UCEES). This work is part of a Ph.D. Thesis conducted in the Biomedicine Doctoral Studies of the University of Granada, Spain. The funders of the present study did not have any role in the study design, data collection and analyses, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
TaqMan™ OpenArray™ custom assay designs of candidate gene SNPs included in the present study. Table S2. Rationale for the inclusion of the 64 single nucleotide polymorphisms (SNPs) in the present study. Table S3. Thermal cycling conditions. Table S4. Further details of the single nucleotide polymorphisms (SNP) included in the present study. Table S5. Socio-demographic and clinical characteristics of the study samples. Table S6. Genotype frequencies of single nucleotide polymorphisms (SNP) in fibromyalgia (FM) and non-fibromyalgia (controls, HC) women.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Estévez-López, F., Camiletti-Moirón, D., Aparicio, V.A. et al. Identification of candidate genes associated with fibromyalgia susceptibility in southern Spanish women: the al-Ándalus project. J Transl Med 16, 43 (2018). https://doi.org/10.1186/s12967-018-1416-8
- Chronic pain
- COMT gene
- GCH1 gene
- OPRM1 gene