Genome-wide association study identifies RNF123 locus as associated with chronic widespread musculoskeletal pain.

BACKGROUND AND OBJECTIVES: Chronic widespread musculoskeletal pain (CWP) is a symptom of fibromyalgia and a complex trait with poorly understood pathogenesis. CWP is heritable (48%-54%), but its genetic architecture is unknown and candidate gene studies have produced inconsistent results. We conducted a genome-wide association study to get insight into the genetic background of CWP. METHODS: Northern Europeans from UK Biobank comprising 6914 cases reporting pain all over the body lasting >3 months and 242 929 controls were studied. Replication of three independent genome-wide significant single nucleotide polymorphisms was attempted in six independent European cohorts (n=43 080; cases=14 177). Genetic correlations with risk factors, tissue specificity and colocalisation were examined. RESULTS: Three genome-wide significant loci were identified (rs1491985, rs10490825, rs165599) residing within the genes Ring Finger Protein 123 (RNF123), ATPase secretory pathway Ca2+transporting 1 (ATP2C1) and catechol-O-methyltransferase (COMT). The RNF123 locus was replicated (meta-analysis p=0.0002), the ATP2C1 locus showed suggestive association (p=0.0227) and the COMT locus was not replicated. Partial genetic correlation between CWP and depressive symptoms, body mass index, age of first birth and years of schooling were identified. Tissue specificity and colocalisation analysis highlight the relevance of skeletal muscle in CWP. CONCLUSIONS: We report a novel association of RNF123 locus and a suggestive association of ATP2C1 locus with CWP. Both loci are consistent with a role of calcium regulation in CWP. The association with COMT, one of the most studied genes in chronic pain field, was not confirmed in the replication analysis.

Enrichment of genomic pathways based on differential DNA methylation profiles associated with chronic musculoskeletal pain in older adults: An exploratory study.

Our study aimed to identify differentially methylated CpGs/regions and their enriched genomic pathways associated with underlying chronic musculoskeletal pain in older individuals. We recruited cognitively healthy older adults with (n = 20) and without (n = 9) self-reported musculoskeletal pain and collected DNA from peripheral blood that was analyzed using MethylationEPIC arrays. We identified 31,739 hypermethylated CpG and 10,811 hypomethylated CpG probes (ps ≤ 0.05). All CpG probes were clustered into 5966 regions, among which 600 regions were differentially methylated at p ≤ 0.05 level, including 294 hypermethylated regions and 306 hypomethylated regions (differentially methylated regions). Ingenuity pathway enrichment analysis revealed that the pain-related differentially methylated regions were enriched across multiple pathways. The top 10 canonical pathways were linked to cellular signaling processes related to immune responses (i.e. antigen presentation, programed cell death 1 receptor/PD-1 ligand 1, interleukin-4, OX40 signaling, T cell exhaustion, and apoptosis) and gamma-aminobutyric acid receptor signaling. Further, Weighted Gene Correlation Network Analysis revealed a comethylation network module in the pain group that was not preserved in the control group, where the hub gene was the cyclic adenosine monophosphate-dependent transcription factor ATF-2. Our preliminary findings provide new epigenetic insights into the role of aberrant immune signaling in musculoskeletal pain in older adults while further supporting involvement of dysfunctional GABAergic signaling mechanisms in chronic pain. Our findings need to be urgently replicated in larger cohorts as they may serve as a basis for developing and targeting future interventions.

Genetic Underpinnings of Musculoskeletal Pain During Treatment With Aromatase Inhibitors for Breast Cancer: A Biological Pathway Analysis.

BACKGROUND: Musculoskeletal pain (MSKP) is the most reported symptom during treatment with aromatase inhibitors (AIs) for breast cancer. The mechanisms underlying MSKP are multidimensional and not well understood. The goals of this biological pathway analysis were to (1) gain an understanding of the genetic variation and biological mechanisms underlying MSKP with AI therapy and (2) identify plausible biological pathways and candidate genes for future investigation. METHOD: Genes associated with MSKP during AI therapy or genes involved in drug metabolism of and response to AIs were identified from the literature. Studies published through February 2019 were queried in PubMed®. The genes identified from the literature were entered into QIAGEN's Ingenuity® Pathway Analysis (IPA) software to generate canonical pathways, upstream regulators, and networks through a core analysis. RESULTS: The 17 genes identified were ABCB1, ABCG1, CYP17A1, CYP19A1, CYP27B1, CYP2A6, CYP3A4, CYP3A5, ESR1, OATP1B1, OPG, RANKL, SLCO3A1, TCL1A, UGT2A1, UGT2B17, and VDR. These genes are involved in encoding bone-remodeling regulators, drug-metabolizing enzymes (cytochrome P450 family, UDP-glucuronosyltransferases family), or drug transporters (ATP-binding cassette transporters, organic anion transporters). Multiple plausible biological pathways (e.g., nicotine degradation, melatonin degradation) and candidate genes (e.g., NFKB, HSP90, AKT, ERK1/2, FOXA2) are proposed for future investigation based on the IPA results. CONCLUSION: Multiple genes and molecular-level etiologies may contribute to MSKP with AI therapy in women with breast cancer. Our innovative combination of gene identification from the literature plus biological pathway analysis allowed for the emergence of novel candidate genes and biological pathways for future investigations.

Interleukin 6 decreases nociceptor expression of the potassium channel KV1.4 in a rat model of hand-arm vibration syndrome.

Chronic muscle pain is a prominent symptom of the hand-arm vibration syndrome (HAVS), an occupational disease induced by exposure to vibrating power tools, but the underlying mechanism remains unknown. We evaluated the hypothesis that vibration induces an interleukin 6 (IL-6)-mediated downregulation of the potassium voltage-gated channel subfamily A member 4 (KV1.4) in nociceptors leading to muscle pain. Adult male rats were submitted to a protocol of mechanical vibration of the right hind limb. Twenty-four hours after vibration, muscle hyperalgesia was observed, concomitant to increased levels of IL-6 in the gastrocnemius muscle and decreased expression of KV1.4 in the dorsal root ganglia. Local injection of neutralizing antibodies against IL-6 attenuated the muscle hyperalgesia induced by vibration, whereas antisense knockdown of this channel in the dorsal root ganglia mimicked the muscle hyperalgesia observed in the model of HAVS. Finally, knockdown of the IL-6 receptor signaling subunit glycoprotein 130 (gp130) attenuated both vibration-induced muscle hyperalgesia and downregulation of KV1.4. These results support the hypothesis that IL-6 plays a central role in the induction of muscle pain in HAVS. This likely occurs through intracellular signaling downstream to the IL-6 receptor subunit gp130, which decreases the expression of KV1.4 in nociceptors.

Multi-OMICS analyses of frailty and chronic widespread musculoskeletal pain suggest involvement of shared neurological pathways.

Chronic widespread musculoskeletal pain (CWP) and frailty are prevalent conditions in older people. We have shown previously that interindividual variation in frailty and CWP is genetically determined. We also reported an association of frailty and CWP caused by shared genetic and common environmental factors. The aim of this study was to use omic approaches to identify molecular genetic factors underlying the heritability of frailty and its genetic correlation with CWP. Frailty was quantified through the Rockwood Frailty Index (FI) as a proportion of deficits from 33 binary health deficit questions in 3626 female twins. Common widespread pain was assessed using a screening questionnaire. OMICS analysis included 305 metabolites and whole-genome (>2.5 × 10 SNPs) and epigenome (∼1 × 10 MeDIP-seq regions) assessments performed on fasting blood samples. Using family-based statistical analyses, including path analysis, we examined how FI scores were related to molecular genetic factors and to CWP, taking into account known risk factors such as fat mass and smoking. Frailty Index was significantly correlated with 51 metabolites after correction for multiple testing, with 20 metabolites having P-values between 2.1 × 10 and 4.0 × 10. Three metabolites (uridine, C-glycosyl tryptophan, and N-acetyl glycine) were statistically independent and thought to exert a direct effect on FI. Epiandrosterone sulphate, previously shown to be highly inversely associated with CWP, was found to exert an indirect influence on FI. Bioinformatics analysis of genome-wide association study and EWAS showed that FI and its covariation with CWP was through genomic regions involved in neurological pathways. Neurological pathway involvement accounts for the associated conditions of aging CWP and FI.

Genetic variant rs3750625 in the 3'UTR of ADRA2A affects stress-dependent acute pain severity after trauma and alters a microRNA-34a regulatory site.

α2A adrenergic receptor (α2A-AR) activation has been shown in animal models to play an important role in regulating the balance of acute pain inhibition vs facilitation after both physical and psychological stress. To our knowledge, the influence of genetic variants in the gene encoding α2A-AR, ADRA2A, on acute pain outcomes in humans experiencing traumatic stress has not been assessed. In this study, we tested whether a genetic variant in the 3'UTR of ADRA2A, rs3750625, is associated with acute musculoskeletal pain (MSP) severity following motor vehicle collision (MVC, n = 948) and sexual assault (n = 84), and whether this influence was affected by stress severity. We evaluated rs3750625 because it is located in the seed binding region of miR-34a, a microRNA (miRNA) known to regulate pain and stress responses. In both cohorts, the minor allele at rs3750625 was associated with increased musculoskeletal pain in distressed individuals (stress*rs3750625 P = 0.043 for MVC cohort and P = 0.007 for sexual assault cohort). We further found that (1) miR-34a binds the 3'UTR of ADRA2A, (2) the amount of repression is greater when the minor (risk) allele is present, (3) miR-34a in the IMR-32 adrenergic neuroblastoma cell line affects ADRA2A expression, (4) miR-34a and ADRA2A are expressed in tissues known to play a role in pain and stress, (5) following forced swim stress exposure, rat peripheral nerve tissue expression changes are consistent with miR-34a regulation of ADRA2A. Together, these results suggest that ADRA2A rs3750625 contributes to poststress musculoskeletal pain severity by modulating miR-34a regulation.

Are Epigenetic Factors Implicated in Chronic Widespread Pain?

BACKGROUND: Chronic widespread musculoskeletal pain (CWP) is the cardinal symptom of fibromyalgia and affects about 12% of the general population. Familial aggregation of CWP has been repeatedly demonstrated with estimated heritabilities of around 50%, indicating a genetic susceptibility. The objective of the study was to explore genome-wide disease-differentially methylated positions (DMPs) for chronic widespread pain (CWP) in a sample of unrelated individuals and a subsample of discordant monozygotic (MZ) twins. METHODOLOGY/PRINCIPLE FINDINGS: A total of N = 281 twin individuals from the TwinsUK registry, including N = 33 MZ twins discordant for self-reported CWP, were part of the discovery sample. The replication sample included 729 men and 756 women from a subsample of the KORA S4 survey-an independent population-based cohort from Southern Germany. Epigenome-wide analysis of DNA methylation was conducted using the Illumina Infinium HumanMethylation 450 DNA BeadChip in both the discovery and replication sample. Of our 40 main loci that were carried forward for replication, three CPGs reached significant p-values in the replication sample, including malate dehydrogenase 2 (MDH2; p-value 0.017), tetranectin (CLEC3B; p-value 0.039), and heat shock protein beta-6 (HSPB6; p-value 0.016). The associations between the collagen type I, alpha 2 chain (COL1A2) and monoamine oxidase B (MAOB) observed in the discovery sample-both of which have been previously reported to be biological candidates for pain-could not be replicated. CONCLUSION/SIGNIFICANCE: Our results may serve as a starting point to encourage further investigation in large and independent population-based cohorts of DNA methylation and other epigenetic changes as possible disease mechanisms in CWP. Ultimately, understanding the key mechanisms underlying CWP may lead to new treatments and inform clinical practice.

Genetic variant in the osteoprotegerin gene is associated with aromatase inhibitor-related musculoskeletal toxicity in breast cancer patients.

BACKGROUND: Aromatase inhibitor (AI) therapy is associated with musculoskeletal (MS) toxicity, which adversely affects quality of life and therapy adherence. Our objective was to evaluate whether genetic variants may predict endocrine therapy-related MS pain and hot flashes in a prospective observational cohort study. PATIENTS & METHODS: 254 early breast cancer patients starting AI (n = 159) or tamoxifen therapy (n = 95) were included in this genetic biomarker study. MS and vasomotor symptoms were assessed at baseline and after 3, 6 and 12 months of therapy. AI-induced MS pain was defined as an increase in arthralgia or myalgia relative to baseline. Single nucleotide polymorphisms (SNP) in candidate genes involved in oestrogen signalling or previously associated with AI-related MS pain or oestrogen levels were selected. RESULTS: Overall, 13 SNPs in CYP19, CYP17, osteoprotegerin (OPG) and oestrogen receptor 1 exhibited an allele frequency >0.05 and were included in the analysis. Patients carrying the G allele of rs2073618 in OPG experienced significantly more AI-induced MS toxicity compared to the wildtype allele, after correction for multiple testing (P = 0.046). Furthermore, this SNP was associated with severity of pain (P = 0.018). No association was found with regard to the other SNPs, both in AI and tamoxifen-treated patients. Neither could an association with vasomotor symptoms be demonstrated. CONCLUSION: The SNP rs2073618 in OPG is associated with an increased risk of MS symptoms and pain with AI therapy, which has not been reported previously. Validation of this finding in larger cohorts and further functional studies are required.

The phenotypic and genetic signatures of common musculoskeletal pain conditions.

Musculoskeletal pain conditions, such as fibromyalgia and low back pain, tend to coexist in affected individuals and are characterized by a report of pain greater than expected based on the results of a standard physical evaluation. The pathophysiology of these conditions is largely unknown, we lack biological markers for accurate diagnosis, and conventional therapeutics have limited effectiveness. Growing evidence suggests that chronic pain conditions are associated with both physical and psychological triggers, which initiate pain amplification and psychological distress; thus, susceptibility is dictated by complex interactions between genetic and environmental factors. Herein, we review phenotypic and genetic markers of common musculoskeletal pain conditions, selected based on their association with musculoskeletal pain in previous research. The phenotypic markers of greatest interest include measures of pain amplification and 'psychological' measures (such as emotional distress, somatic awareness, psychosocial stress and catastrophizing). Genetic polymorphisms reproducibly linked with musculoskeletal pain are found in genes contributing to serotonergic and adrenergic pathways. Elucidation of the biological mechanisms by which these markers contribute to the perception of pain in these patients will enable the development of novel effective drugs and methodologies that permit better diagnoses and approaches to personalized medicine.

Dynamic variations of c-Fos expression in the spinal cord exposed to intramuscular hypertonic saline-induced muscle nociception.

BACKGROUND: Numerous studies focus on intramuscular (i.m.) injection of hypertonic saline-induced muscle pain and nociception. The spatio-temporal characteristics and dynamic variation of spinal neuronal activities elicited by i.m. hypertonic saline remain unknown. METHODS: The spatio-temporal variations of c-Fos expression in the lumbar spinal cord exposed to i.m. injection of 5.8% saline were investigated in male rats. RESULTS: After a unilateral i.m. 5.8% saline injection, c-Fos expression in dorsal horn of spinal L4-6 segments was significantly enhanced bilaterally (p < 0.05 and p < 0.001). These 5.8% saline-induced bilateral spinal Fos expression occurred rapidly 0.5 h following the i.m. injection, and reached the peak levels within 1 h, which declined gradually and returned to the control levels within 8 h. Compared with intact rats without i.m. insults, no significant influence of the spinalization on spinal c-Fos expression was found. However, the 5.8% saline-induced increases in Fos expression in intact and spinalized rats differed significantly. During muscle nociception, the c-Fos expression within the superficial layer (laminae I-II) and the deep layer (laminae V-VI) in spinalized rats were significantly lower and higher than that of in intact rats (p < 0.05 and p < 0.001). Fentanyl (20 µg/kg, intraperitoneal) completely attenuated the 5.8% saline intramuscularly induced increases in c-Fos expression in laminae III-VI (p < 0.001), but not laminae I-II. CONCLUSIONS: It is suggested that spinal nociceptive neuronal activities in superficial and deep layers may differently be modulated by endogenous descending facilitation and inhibition, respectively.

Catechol-O-methyltransferase gene polymorphism and chronic human pain: a systematic review and meta-analysis.

In human studies, low COMT (catechol-O-methyltransferase) activity has been associated with increased sensitivity to acute clinical preoperative or postoperative pain. We explored the association between the COMT genotype and three chronic pain conditions: migrainous headache, fibromyalgia, or chronic widespread pain and chronic musculoskeletal pain. Furthermore, we evaluated whether COMT genotype affects the efficacy of opioids in chronic pain. After a systematic literature review, we carried out meta-analyses on the three chronic pain conditions. The efficacy of opioids was evaluated using a systematic review only. The meta-analyses showed that fibromyalgia or chronic widespread pain is the only type of chronic pain that could be associated with the COMT single nucleotide polymorphism rs4680 (Val158Met). Met158, which results in the low-activity variant of COMT, is the risk allele. In chronic clinical pain, the effect of the COMT polymorphism depends on the pain condition. Low COMT activity is not associated with migrainous headache or chronic musculoskeletal pain conditions, but it may increase the risk for fibromyalgia or chronic widespread pain. Low COMT activity increases opioid receptors and enhances opioid analgesia and adverse effects in some cancer pains. Findings from animal studies that have utilized COMT inhibitors elucidate the mechanism behind these findings. In rodent pain models, COMT inhibitors are pronociceptive, except for neuropathic pain models, where nitecapone was found to be antiallodynic. The complex interplay between enhanced adrenergic and dopaminergic activity in different parts of the nociceptive system probably explains the complicated actions of low COMT activity.

Genetic variation in the beta-2 adrenergic receptor is associated with chronic musculoskeletal complaints in adolescents.

BACKGROUND: There is significant evidence to suggest that psychological and stress-related factors are important predictors of the onset of chronic widespread pain (CWP) and fibromyalgia (FM). The hypothalamic-pituitary-adrenal axis, together with the efferent sympathetic/adrenomedullary system, influence all body organs (including muscles) during short- and long-term threatening stimuli. The aim of this study was to investigate the relationship between genetic variants in adrenergic candidate genes and chronic musculoskeletal complaints (MSCs) in adolescents. METHODS: Adolescents from the Western Australian Pregnancy (Raine) Cohort attending the 17-year cohort review completed a questionnaire containing a broad range of psychosocial factors and pain assessment (n = 1004). Blood samples were collected for DNA extraction and genotyping. Genotype data was obtained for 14 single nucleotide polymorphisms (SNPs) in two candidate genes - beta-2 adrenergic receptor (ADRB2) and catecholamine-O-methyltransferase (COMT). Haplotypes were reconstructed for all individuals with genotype data. RESULTS AND CONCLUSION: Both female gender and poor mental health were associated with (1) an increased risk for chronic, disabling comorbid neck and low back pain (CDCP); and (2) an increase in the number of areas of pain. Of the 14 SNPs evaluated, only SNP rs2053044 (ADRB2, recessive model) displayed an association with CDCP [odds ratio (OR) = 2.49; 95% confidence interval (CI) = 1.25, 4.98; p = 0.01] and pain in three to four pain areas in the last month (OR = 1.86; 95% CI = 1.13, 3.06; p = 0.02). These data suggest that genetic variants in ADRB2 may be involved in chronic MSCs.

Aromatase inhibitors, estrogens and musculoskeletal pain: estrogen-dependent T-cell leukemia 1A (TCL1A) gene-mediated regulation of cytokine expression.

INTRODUCTION: Arthralgias and myalgias are major side effects associated with aromatase inhibitor (AI) therapy of breast cancer. In a recent genome-wide association study, we identified SNPs - including one that created an estrogen response element near the 3' end of the T-cell leukemia 1A (TCL1A) gene - that were associated with musculoskeletal pain in women on adjuvant AI therapy for breast cancer. We also showed estrogen-dependent, SNP-modulated variation in TCL1A expression and, in preliminary experiments, showed that TCL1A could induce IL-17RA expression. In the present study, we set out to determine whether these SNPs might influence cytokine expression and effect more widely, and, if so, to explore the mechanism of TCL1A-related AI-induced side effects. METHODS: The functional genomic experiments performed included determinations of TCL1A, cytokine and cytokine receptor expression in response to estrogen treatment of U2OS cells and lymphoblastoid cell lines that had been stably transfected with estrogen receptor alpha. Changes in mRNA and protein expression after gene knockdown and overexpression were also determined, as was NF-κB transcriptional activity. RESULTS: Estradiol (E2) increased TCL1A expression and, in a TCL1A SNP-dependent fashion, also altered the expression of IL-17, IL-17RA, IL-12, IL-12RB2 and IL-1R2. TCL1A expression was higher in E2-treated lymphoblastoid cell lines with variant SNP genotypes, and induction of the expression of cytokine and cytokine receptor genes was mediated by TCL1A. Finally, estrogen receptor alpha blockade with ICI-182,780 in the presence of E2 resulted in greatly increased NF-κB transcriptional activity, but only in cells that carried variant SNP genotypes. These results linked variant TCL1A SNP sequences that are associated with AI-dependent musculoskeletal pain with increased E2-dependent TCL1A expression and with downstream alterations in cytokine and cytokine receptor expression as well as NF-κB transcriptional activity. CONCLUSIONS: SNPs near the 3' terminus of TCL1A were associated with AI-dependent musculoskeletal pain. E2 induced SNP-dependent TCL1A expression, which in turn altered IL-17, IL-17RA, IL-12, IL-12RB2, and IL-1R2 expression as well as NF-κB transcriptional activity. These results provide a pharmacogenomic explanation for a clinically important adverse drug reaction as well as insights into a novel estrogen-dependent mechanism for the modulation of cytokine and cytokine receptor expression.

Influence of catechol-O-methyltransferase (COMT) gene polymorphisms in pain sensibility of Brazilian fibromialgia patients.

Fibromyalgia syndrome (FS) is a rheumatic syndrome affecting to 2-3% of individuals of productive age, mainly women. Neuroendocrine and genetic factors may play a significant role in development of the disease which is characterized by diffuse chronic pain and presence of tender points. Several studies have suggested an association between FS, especially pain sensitivity, and polymorphism of the catechol-O-methyltransferase (COMT) gene. The aim of the present study was to characterize the SNPs rs4680 and rs4818 of the COMT gene and assess its influence in pain sensitivity of patients with fibromyalgia screened by the Fibromyalgia Impact Questionnaire (FIQ). DNA was extracted from peripheral blood of 112 patients with fibromyalgia and 110 healthy individuals and was used as template in PCR for amplification of a 185-bp fragment of the COMT gene. The amplified fragment was sequenced for analyses of the SNPs rs4680 and rs4818. The frequency of mutant genotype AA of SNP rs6860 was 77.67% in patients with FS and 28.18% for the control group. For the SNP rs4818, the frequency of mutant genotype CC was 73.21 and 39.09% for patients with FS and controls, respectively. Moreover, the FIQ score was higher in patients with the homozygous mutant genotype for SNPs rs4680 (87.92 points) and rs4818 (86.14 points). These results suggest that SNPs rs4680 and rs4818 of the COMT gene may be associated with fibromyalgia and pain sensitivity in FS Brazilian patients.