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

* Md Shafiqur Rahman
* Bendik S Winsvold
* S.O. Chavez Chavez
* Sigrid Børte
* Yakov A. Tsepilov
* Sodbo Zh. Shapov
* HUNT All-In Pain
* Yurii Aulchenko
* Knut Hagen
* Egil A. Fors
* Kristian Hveem
* John-Anker Zwart
* J.B.J. van Meurs
* Maxim B. Freidin
* Frances M.K. Williams

## Abstract

**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 6,914 cases reporting pain all over the body lasting more than 3 months and 242,929 controls were studied. Replication of three lead genome-wide significant single nucleotide polymorphisms (SNPs) was attempted in 6 independent European cohorts (N=43,080; cases=14,177). Genetic correlations with risk factors, tissue specificity, and colocalization were examined.

**Results** Three genome-wide significant loci were identified (*rs1491985, rs10490825, rs165599*) residing within the genes *RNF123, ATP2C1*, and *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 colocalization analysis highlight the relevance of skeletal muscle in CWP.

**Conclusions** We report a novel association of *RNF123* locus with CWP and suggest a role of *ATP2C1*, consistent with a role of calcium regulation in CWP. The association to *COMT*, one of the most studied genes in chronic pain field, was not confirmed in the replication analysis.

**What is already known about this subject?**

*   Chronic widespread musculoskeletal pain (CWP) is a primary diagnostic feature of fibromyalgia.

*   CWP is moderately heritable, but precise genes involved in the pathogenesis of CWP are yet to be identified.

**What does this study add?**

*   This is the largest genetic study conducted on CWP to date and identified novel genetic risk loci (*RNF123* and *ATP2C1*).

*   The genetic signal points to peripheral pain mechanisms in CWP, and shows genetic correlation with other traits, including BMI and depression.

**How might this impact on clinical practice or future developments?**

*   The findings add to etiological basis of CWP.

Keywords
*   Chronic widespread pain
*   fibromyalgia
*   genome-wide association study
*   twin
*   *RNF123*

## Introduction

Chronic widespread musculoskeletal pain (CWP) is a common complex trait influenced by genetic and environmental factors, most of which have yet to be determined[1]. CWP and fibromyalgia syndrome are sometimes used interchangeably although the latter is generally more severe and includes other features such as sleep disturbance, fatigue and depression [2]. It is thought to represent a subgroup at the more severe end of the spectrum of CWP[3]. The prevalence of CWP is 10.6% in the world population and 14.2% in the UK population[4, 5]. It is associated with high societal cost[6]. CWP is responsible for excess mortality[7], which is thought to be attributable to cardiovascular disease, respiratory disease and cancer. Females are more affected by CWP than males[4], and the prevalence rises with age[5]. In addition to age and sex, a number of exposures have been proposed as risk factors for CWP[8, 9] but only increased body mass index (BMI) has been consistently reported across studies, including longitudinal studies[10-12].

Heritability estimates for CWP range between 48-54%, indicating a substantial genetic contribution[13]. To date, the candidate gene approach has been extensively applied to identify genetic factors in CWP[14] (see below) but few agnostic studies have been published[15]. The only genome-wide association study (GWAS) meta-analysis combining 14 studies identified a locus lying on chromosome 5 intergenic to *CCT5* and *FAM173B*[15]. *CCT5* has previously been implicated in neuropathy[16] and there is increasing evidence that small fiber neuropathy underlies a subset of fibromyalgia[17].

Genetic factors are known to be shared by chronic pain conditions[18, 19]. One of the most extensively studied chronic pain-associated genes encodes catechol-O-methyltransferase (COMT), an enzyme which regulates the production of catecholamines that act as neurotransmitters in the central nervous system pain tract. A non-synonymous change of A to G encoding a valine (Val) to methionine (Met) substitution at codon 158 (Val158Met; *rs4680*) reduces the enzymatic activity of *COMT*. This single nucleotide polymorphism (SNP) has been reported to be associated with CWP in a small study of 122 participants[20], but a subsequent association study of 3,017 participants did not confirm earlier findings[21]. An inconclusive role of *COMT* was observed for temporomandibular disorders (TMD) as well[22, 23]. Further investigation is required to identify genetic variants underlying CWP which will shed light on the pathophysiological mechanisms underlying the development of chronic pain and may reveal therapeutic targets.

## Materials and Methods

An overview of study design presented in figure 1.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F1)

Figure 1. 
Overview of study design.

### Participant selection

For the discovery analysis, we performed a GWAS of CWP using UK Biobank (UKB) comprising 249,843 participants of European descent (6,914 CWP cases and 242,929 controls). Independent SNPs passing a threshold p<5.0E-08 were submitted for replication in 43,080 individuals of European ancestry (14,177 CWP cases and 28,903 controls) from six independent cohorts originating in the UK (TwinsUK and The English Longitudinal Study of Ageing (ELSA)), the Netherlands (The Rotterdam Study I, II and III (RS-I, RS-II and RS-III)), and Norway (The Nord-Trøndelag Health Survey (HUNT)). All study participants provided written informed consent. Study-specific ethical approval was provided by the local ethics committees. The UKB dataset was used under project #18219. Description of each study cohort is presented in supplementary methods.

### Phenotype

In UKB, CWP cases were defined by combining self-reported diagnosis of pain all over the body lasting for more than 3 months; simultaneous pain in the knee, shoulder, hip and back lasting 3+ months; and doctor’s diagnosis of fibromyalgia. Controls comprised those who reported no pain in the last month or reported pain all over the body in the previous month that did not last for three months or reported only ≥ 3 months of non-musculoskeletal pain (headache, facial and abdominal pain). Those reporting a diagnosis of rheumatoid arthritis, polymyalgia rheumatica, arthritis not otherwise specified, systemic lupus erythematosus, ankylosing spondylitis or myopathy were excluded from the study (supplementary figure S1). Further phenotype details for UKB and replication cohorts are provided in supplementary methods.

### Genotyping

Genotyping and imputation methods across cohorts are summarised in supplementary table S1. In brief, UKB participants were genotyped using Applied Biosystems UKB Axiom array and Applied Biosystems UKB Lung Exome Variant Evaluation Axiom array. HUNT participants were genotyped using Illumina/HumanCoreExome12 v1.0, Illumina/HumanCoreExome12 v1.1 and UM HUNT Biobank v1.0. ELSA participants were genotyped using the Illumina/HumanOmni2.5-4v1 and Illumina/HumanOmni2.5-8v1.3. TwinsUK participants were genotyped using Illumina/HumanHap300, Illumina/HumanHap610Q, Illumina/1M-Duo and 1.2MDuo 1M. RS-I participants were genotyped using Illumina/HumanHap 550K V.3 and Illumina/HumanHap 550K V.3 DUO. RS-II participants were genotyped using Illumina/HumanHap 550V.3DUO and Illumina/HumanHap610Q. RS-III genotyping was performed using Illumina/HumanHap610Q.

### Statistical analysis and in-silico follow-up

The details of statistical analysis, and in-silico follow-up are described in supplementary methods. In brief, GWAS in the discovery sample was performed using linear mixed-effects model implemented in BOLT-LMM v2.3.2[24]. An additive genetic model for SNP effect on CWP was adjusted for age, genetically determined sex, genotyping platform, and the first ten genetic principle components[25]. A sensitivity GWAS (controls: 223,606 and CWP cases: 6,914) was performed excluding participants with chronic non-musculoskeletal pain such as headache, facial and abdominal pain from the controls. Independent SNPs at GWAS significant loci were identified using Conditional & Joint[26] (COJO) analysis and submitted for replication. Independent SNPs across all replication cohorts were meta-analysed using fixed-effects model with both sample size, and inverse-variance weighting implemented in METAL[27]. LD score regression (LDSR)[28] was used to estimate inflation (λGC) in test statistics, SNP-based heritability on the observed and liability scale and genetic correlations. We also estimated partial genetic correlations[29]. We used FUMA[30] for the functional annotation of SNPs, gene mapping, tissue specificity and gene-set enrichment. Differential expression of replicated lead SNP was assessed using the GTEx version-8 tissues. Colocalization of GWAS independent SNPs in human skeletal muscle and dorsal root ganglion (DRG) tissues was assessed utilizing publicly available data[31, 32].

## Results

Details of the discovery and replication cohorts are presented in table 1. Cases were enriched for females compared to controls in all cohorts (p<0.001) and were on average older in the discovery, and in four replication cohorts (p<0.05). In all cohorts BMI was significantly higher in cases than controls (p<0.0001) except for RS-3 where a similar but non-significant trend was observed (p=0.0827).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/T1)

Table 1. 
Sample characteristics stratified by case / control status for discovery and replication cohorts.

### Discovery genome-wide association study

Three genomic loci tagged by *rs1491985, rs10490825*, and *rs165599* passed genome-wide significance threshold of p<5E-08 (figure 2). Observed inflation in test statistics (λGC = 1.146, supplementary figure S2) was due to polygenicity (LDSR intercept: 1.002±0.0085, LDSR ratio=0.0118±0.0497) rather than confounding. SNP-heritability of CWP was 0.033±0.004 on the observed scale, and 0.22±0.02 on the liability scale meaning that the observed SNPs explain 22% of variance in CWP risk. Lead SNPs were located in the gene *RNF123* (chromosome 3, *rs1491985*, intronic variant, p=1.60E-08), *ATP2C1* (chromosome 3, *rs10490825*, intronic variant, p=1.30E-08), and *COMT* (chromosome 22, rs165599, 3’-UTR variant, p=2.50E-08), respectively (figure 3A,C; supplementary table S2). Six additional loci near or within genes *HNRNPA1P46, LRRC3B, PDE6A, DPYSL2, ANXA11*, and *AL138498*.1 were identified at suggestive GWAS threshold of p<5E-07. Sensitivity GWAS excluding participants with chronic non-musculoskeletal pain provided findings similar to original results except that *COMT* locus now became suggestively significant (p=5.3E-08) (supplementary figure S3).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F2)

Figure 2. 
Manhattan plot of summary GWAS of chronic wide-spread pain (CWP) derived from UK biobank European ancestry data. Each circle in the plots represents an SNP, which was positioned following genomic build GRCh37. Colored circles indicated linkage disequilibrium (LD). The y-axis shows the corresponding –log10 P-values based on UKB discovery cohort and the x-axis show chromosome positions along with gene names and flaking size.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F3)

Figure 3. 
Regional plots for three independent CWP associated SNPs. Independent SNPs are colored in purple. Other colored circles indicated pairwise LD. Strength of LD (r2) presented in the upper right corner of each plot.

### Replication results and meta-analysis

Results are presented in supplementary table S3, with meta-analysis of the 6 replication samples in figure 4 (supplementary table S4,S5). Given the significance threshold for replication: 0.05/3=0.017, association between CWP and *rs1491985* was considered replicated (sample-size based p=0.0002; standard-error based p=0.0003). *Rs10490825* showed suggestive association with CWP (sample-size based p=0.0227; standard-error based p=0.0490) and demonstrated a consistent direction of effect in 5 of the 6 replication samples. *Rs165599* did not replicate (sample-size based p-value=0.7300; standard-error based p-value=0.5000) and the direction of effect was not consistent across cohorts: in three cohorts, allele A was protective, while in the other three it was the risk allele. None of the three SNPs displayed statistically significant heterogeneity in the replication cohorts.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F4)

Figure 4. 
Forest plot of the association for (A) *rs1491985*, (B) *rs10490825*, and (C) *rs165599* with chronic widespread pain (CWP). X-axis shows, effect size measures are presented as beta value. The red square with horizontal black line represents the cohort-specific effect with a corresponding confidence interval for the SNP of interest. Size of the square indicates the weight of the study and reflects sample size. The vertical black line indicates “line of no effect”. Overall effect is presented as a black diamond. Test-statistics for each cohort, meta-analysis and heterogeneity are available on the left-hand side. The *rs1491985*, and *rs10490825* were not present in ELSA, therefore the *rs9870858*, and *rs1732984*8 were used as proxy, respectively (online supplementary methods).

### CWP shares genetic components with BMI, depression, age at first birth, and years of schooling

Two hundred and nine traits from sources other than UKB were examined for genetic correlation with CWP. We selected traits for which the absolute value of the correlation coefficient (rg) was >0.2, and for which the Bonferroni-corrected p-value was <0.01/209=4.78E-05. Twenty-three traits fulfilled these criteria (supplementary figure S4). The highest positive genetic correlation was observed for depressive symptoms (rg=0.65) and the highest negative correlation was observed for college completion (rg= −0.61). Many of the 23 genetically correlated traits were correlated with each other raising concerns about their independency of correlations with CWP. We therefore calculated partial genetic correlations conditionally independent of each other.

Using hierarchical clustering, based on genetic correlations we identified 7 clusters (supplementary figure S5A), with 7 traits selected to represent each cluster (BMI, triglycerides, depressive symptoms, coronary artery disease, smoking, age of first birth, and years of schooling) to quantify partial genetic correlation with CWP. We found depressive symptoms (rg=0.59), BMI (rg=0.20), age of first birth (rg=-0.26), and years of schooling (rg=-0.17) independently correlated with CWP (supplementary figure S5B, supplementary table S6).

### Tissue-specific expression of CWP mapped gene-sets

The results of functional annotation of GWAS significant loci are presented in supplementary figure S6 (supplementary results). Four different gene mapping strategies were implemented in FUMA (genome-wide gene-based association analysis, positional, eQTL, and chromatin interaction mapping) linking annotated SNPs to 89 genes of which 9 *(MST1, GMPPB, APEH, RNF123, ARVCF, AMIGO3, IP6K1, TANGO2*, and *TRAIP)* were identified using all four methods (figure 5A,B)[33]. At genomic risk loci 3p21.31, 3q22.1, and 22q11.21 a total of 54, 10 and 3 genes were prioritized by either eQTL or chromatin interaction mapping, respectively (figure 5C,D). Mapped genes were investigated for tissue-specific gene expression and for gene-set enrichment. In 54 specific GTEx tissues types, differentially expressed gene sets were enriched for skeletal muscle, several brain tissues, heart, whole blood, pancreas, and transverse colon (figure 6A, supplementary table S7). In 30 general GTEx tissue types, differentially expressed gene sets enriched for skeletal muscle, pancreas, heart, blood and brain (figure 6B, supplementary table S8). In both sets of GTEx tissues, overall enrichment for differentially expressed gene sets containing *RNF123* and *ATP2C1* gene were stronger for skeletal muscle than other tissues. *RNF123* was found to be highly expressed in skeletal muscle compared to other tissue types (figure 6C). None of the hallmark gene sets available in the molecular signature database was identified in the analysis.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F5)

Figure 5. 
(A) Manhattan plot of the genome-wide gene-based association analysis, (B) Venn diagram showing overlap of genes implicated by genome-wide gene-based analysis implemented in MAGMA, positional mapping, chromatin interaction mapping (Ci Map), and eQTL mapping, (C) & (D) The circus plot displaying chromatin interactions and eQTLs on chromosomes 3 and chromosomes 22, respectively. The y-axis shows the ─log10 transformed two-tailed p-value of each gene from a linear model and the chromosomal position on the x-axis. The red dotted line indicates the Bonferroni-corrected threshold for genome-wide significance of the gene-based test. The most outer layer of the circus plot displaying Manhattan plot with –log10 P-values for CWP associated SNPs. Each independent SNP is presented with RSID. LD relationship between independent SNPs at the locus and their proxies are indicated with red (r 2 > 0.8) and orange (r 2 > 0.6). Grey SNPs indicate minimal LD with r2 ≤0.20. The outer circle represents chromosome with genomic risk loci are highlighted in blue. Either chromatin interaction or eQTLs mapped genes are displayed on the inner circle. Chromatin interactions and eQTLs mapped genes are presented in orange or green color, respectively. Genes mapped with both approaches are colored red.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/03/2020.11.30.20241000/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2020/12/03/2020.11.30.20241000/F6)

Figure 6. 
(A) Differentially expressed gene (DEG) plots for CWP in 54 tissue types from GTEX v8, (B) DEG plots for CWP in 30 general tissue types from GTEX v8 and (C) Differential expression of *RNF123* gene across tissue types from GTEX v8 (A, B) In both plots, the y-axis represents the ─log10 transformed two-tailed P-value of the hypergeometric test. Significantly enriched DEG sets (Bonferroni corrected p-value < 0.05) are highlighted in red. (C) Y-axis represents TPM (transcripts Per Million) and X-axis represents the GTEx (Version-8) tissues. The figure was adapted from GTEx Portal ([https://www.gtexportal.org/home/gene/ENSG00000164068](https://www.gtexportal.org/home/gene/ENSG00000164068)).

### Colocalization with skeletal muscle and dorsal root ganglion eQTLs

Colocalization analysis identified 93% of probability of shared eQTL variant *rs6809879* which controls *CDHR4* expression in the skeletal muscle and CWP association signal near *RNF123* locus (supplementary table S9, supplementary figure S7A). Additionally, significant colocalization was found for *rs13093525*, which controls *APEH* expression in DRG at exon-level (72% of probability of shared variant with *RNF123* locus). Both *rs6809879* and *rs13093525* were in complete LD with lead SNP *rs1491985(R2=1)* (supplementary table S10, supplementary figure S7B). No evidence of skeletal muscle or DRG eQTL colocalization was observed for *ATP2C1* and *COMT* loci.

## Discussion

CWP is a highly prevalent condition with moderate heritability and serves as a cardinal diagnostic feature of fibromyalgia. Therefore, our findings are of importance for better understanding the genetic basis of fibromyalgia. We report here the largest GWAS of CWP to date using 249,843 participants from the UKB, identifying 3 genome-wide significant loci implicating *RNF123, ATP2C1*, and *COMT* genes. The association in *RNF123* was replicated, whereas *ATP2C1* showed a suggestive association, and the *COMT* locus did not replicate in 41,296 individuals from independent cohorts.

*RNF123* (Ring Finger Protein, also known as KPC1) gene encodes E3 ubiquitin-protein ligase, which has a role in cell cycle progression, metabolism of proteins and innate immunity[34, 35]. Gene expression was higher in skeletal muscle than other tissues. It is not clear how *RNF123* contributes to CWP but it does have a role in malignancy[36, 37], osteoclast differentiation[38] and depression[39]. Consistent with this, an increased cancer incidence and poor cancer survival[40], poor bone health[41, 42], and co-existence of depression[43] have been reported in patients with fibromyalgia suggesting a shared pathway between these conditions and CWP. However, colocalization of *RNF123* locus with skeletal muscle, and DRG eQTLs reiterates the functional relevance of this locus in CWP and is in keeping with the symptomatic skeletal muscle tenderness of which patients complain.

*ATP2C1* (ATPase secretory pathway Ca2+ transporting 1) demonstrated a suggestive association in replication (p=0.0227). There was a consistent direction of effect for *ATP2C1* locus in all six replication cohorts except for ELSA, where we used a proxy SNP, which had close to zero effect size (beta= −0.0004±0.0110). Both positional and eQTL gene mapping strategies implicated *ATP2C1* with CWP. The *ATP2C1* gene encodes for the adenosine triphosphate (ATP)-powered magnesium-dependent calcium pump protein hSPCA1, which mediates Golgi uptake of cytosolic Ca(2+) and Mg(2+)[44]. A loss of function mutation in the *ATP2C1* leads to autosomal dominant Hailey-Hailey disease (HHD), a skin condition characterized by blistering and erosion of the epidermis[45]. Interestingly, published evidence shows HHD can be treated successfully with low-dose naltrexone (LDN), an opioid receptor antagonist, which has been used in the management of fibromyalgia, although opioids are now not recommended[46, 47]. Additionally, the role of calcium regulation in pain processing is well known[48-50]. Pregabalin is an α2δ calcium channel antagonist, which is currently used in the management of fibromyalgia[51], all supporting the regulation of calcium influencing CWP/fibromyalgia. This is the first study to implicate *ATP2C1* with musculoskeletal pain using an agnostic approach.

*COMT* is one of the most studied genes in human pain[52]. Almost 30 SNPs and three haploblocks of the *COMT* gene have been studied in acute clinical, experimental and chronic pain. *Rs4680* of the *COMT* gene is extensively studied in many pain phenotypes such as pain sensitivity, TMD, and fibromyalgia[53]. Across multiple ethnic populations, *rs4680* was implicated in the development of fibromyalgia[14]. However, a meta-analysis of 8 case-control studies (589 fibromyalgia cases and 527 controls) did not confirm earlier association[54]. To date, the largest study assessed the association between *COMT* haplotypes (*rs4680, rs4818, rs4633*, and *rs6269*), and fibromyalgia included 60,367 participants (2,713 ICD-9 diagnosed fibromyalgia) and found no association[55]. They have also been refuted in other European CWP samples[21, 56] and a large candidate gene study of fibromyalgia[57]. However, we identified *rs165599*, located at 3’UTR of *COMT*, associated with CWP in the discovery sample but not in the meta-analysis or any of the six replication cohorts. This variant is not in LD with previously studied *COMT* SNPs *rs4680, rs4818, rs4633* and *rs6269*, and was found not to be associated with chronic musculoskeletal pain including CWP neither when studied as single SNP nor as a part of a haploblock[58-60]. Several explanations of our non-replication of *COMT* locus are possible. First, there was lower power pertaining to overall meta-analysis, which was estimated at 48%, based on the effect size observed in the discovery sample (n=249,843), replication sample size (n=43,080), and the number of tests conducted (n=3). Our meta-analysis did have 90% power to detect a relative risk as small as 1.04 but the estimated *COMT* effect was only 1.012 (beta=0.0027±0.004; OR=1.012, 95% CI=0.97–1.05). Our replication sample size was, however, larger than many of the earlier studies (included 122 to 3811 participants) that reported the association between *COMT* and CWP[20, 61]. Second, we observed a tendency towards non-significance for the *COMT* locus in the sensitivity GWAS due to the exclusion of participants with non-musculoskeletal pain from the control group suggesting that *COMT* predisposes to chronic pain in general. Finally, genetic factors underlying chronic pain and psychiatric co-morbidity (e.g., depression and neuroticism) are known to be shared[62]. However, previous GWAS on chronic pain[29, 63, 64], depression[65] and neuroticism[66] have failed to detect an association with *COMT*. Thus, if there is a role of *COMT* in CWP, it is likely minimal.

Epidemiological studies have consistently reported higher BMI to be associated with an increased risk of CWP[10-12]. Our analysis showed significantly higher BMI in CWP cases compared to controls (p<0.0001) in all cohorts except RS-3. In line with this, we observed a positive genetic overlap between BMI and CWP independent of genetic confounders. Similarly, genetically independent pairwise genetic correlation for depressive symptoms, age of first birth, and years of schooling was seen with CWP. These findings indicate the presence of shared molecular pathways underlying these traits.

Functional analysis showed that prioritised genes are differentially expressed in skeletal muscle, several areas of the CNS, pancreas, whole blood and heart tissues. These findings suggest the involvement of nervous, musculoskeletal, and neuroendocrine systems in CWP. These physiological systems have been implicated in fibromyalgia in previous studies[67-69]. Evidence suggests that both peripheral and central pain mechanisms influence CWP[70, 71]. We observed overall stronger enrichment for differentially expressed gene sets in skeletal muscle than in other GTEx tissues such as brain. Also, skeletal muscle and DRG eQTLs colocalize with the *RNF123* locus. These findings suggest a substantial involvement of peripheral pain mechanisms in CWP. However, we established tissue specificity in the expression of multiple genes likely related to CWP. An in-depth investigation of these genes may provide additional insight into the pathogenesis of CWP.

There are several limitations to the study. The case definition of CWP depends on self-report together with exclusion of other conditions with symptoms leading to chronic pain[72]. A clinical diagnosis of CWP would have been infeasible in a sample this large. Also, we used common SNPs to estimate the heritability of CWP, so the contribution of other variants in the heritability estimated remains unknown. Other musculoskeletal pain assessed in UKB (knee, hip, back, and neck or shoulder pain) exhibited SNP-heritability on the liability scale in the range 0.08–0.12[62], and these estimates are considerably lower than our estimate for CWP, suggesting CWP is a trait of high genetic influence.

In summary, this study identified a novel association for CWP in the gene *RNF123* and suggested the role of calcium regulation, by the involvement of the *ATP2C1* locus. The association of the *COMT* locus with CWP was not replicated, suggesting a small influence, if any. We found evidence that the epidemiological association of BMI and CWP is at least in part genetically mediated. Finally, our results suggest a profound role of peripheral mechanisms in the pathogenesis of CWP.

## Supporting information

Supplementary Material [[supplements/241000_file03.docx]](pending:yes)

## Data Availability

Summary statistics from our GWAS discovery was deposited at Zenodo (https://doi.org/10.5281/zenodo.4298764). Other data relevant to the study are included in the article or uploaded as supplementary information.

## Competing interests

YSA is co-owner of Maatschap PolyOmica and PolyKnomics BV, private organizations, providing services, research and development in the field of computational and statistical, quantitative and computational (gen)omics.

## Contributors

MDSR prepared the manuscript. Discovery GWAS: MDSR, MF and FW designed the study. MDSR, YAT, MF and FW developed phenotype definition. MDSR performed data management and GWAS in the discovery cohorts. Replication: MDSR, MF and FW developed replication protocol and coordinated the study; MDSR, BW, JAW, MF and FW designed the study; BW, SB, KH, EF, KHveem, JAZ, HUNT All-in Pain authors, SOCC, JBJVM and FW collected replication cohort’s data; MDSR, BW, SB, KHveem, JAZ, SOCC and JBJVM were involved data management; MDSR, BW, SB, KH, EF, JAZ, HUNT All-in Pain authors, MF and FW developed phenotype definition; MDSR, BW, SB, SOCC and JBJVM analysed the data. MDSR performed replication meta-analysis. Post-GWAS bioinformatics: MDSR, YAT and SSh performed post-GWAS bioinformatics analysis. MF, YA and FW provided the statistical and bioinformatics consultation in the project. MDSR, KH, SOCC, YAT, SSh, YA, JBJVM, MF and FW interpreted the findings. FW was responsible for the overall supervision of the project. All authors commented on the manuscript and agreed on the final version to be published.

## Funding

MDSR received funding from the European Union’s Horizon 2020 research and innovation program IMforFUTURE, under H2020-MSCA-ITN grant agreement number 721815. The work of YAT was supported by PolyOmica and by the Russian Foundation for Basic Research (project 19-015-00151). The work of SZS was supported by the Federal Agency of Scientific Organizations via the Institute of Cytology and Genetics (project 0324-2019-0040-C-01/AAAA-A17-117092070032-4). The work of YSA was supported by PolyOmica.

## Data availability statement

Summary statistics from our GWAS discovery was deposited at Zenodo ([https://doi.org/10.5281/zenodo.4298764](https://doi.org/10.5281/zenodo.4298764)). Other data relevant to the study are included in the article or uploaded as supplementary information.

## Patient and public involvement

Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

## Patient consent for publication

Not required

## Ethics approval

Study-specific ethical approval was provided by the local ethics committees. Details provided in the supplementary material.

## Acknowledgements

We would like to thank all the participants of UK Biobank, The Nord-Trøndelag Health Survey, English Longitudinal Study of Aging, TwinsUK and Rotterdam study I, II and III. The contribution of inhabitants, general practitioners and pharmacists of the Ommoord district to the Rotterdam Study is gratefully acknowledged.

The UK Biobank study was approved by the National Health Service National Research Ethics Service (ref.11/NW/0382) and all participants provided written informed consent. Genome-wide association analysis was performed using the UK Biobank resource under project number 18219. The English Longitudinal Study of Ageing is jointly run by University College London, Institute for Fiscal Studies, University of Manchester and National Centre for Social Research. Genetic analyses have been carried out by UCL Genomics and funded by the Economic and Social Research Council and the National Institute on Aging. Data governance was provided by the METADAC data access committee, funded by ESRC, Welcome, and MRC. (2015-2018: Grant Number MR/N01104X/1 2018-2020: Grant Number ES/S008349/1). TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust in partnership with King’s College London. The Nord-Trøndelag Health Study (The HUNT Study) is a collaboration between HUNT Research Centre (Faculty of Medicine and Health Sciences, NTNU, Norwegian University of Science and Technology), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health. The genotyping was financed by the National Institute of Health (NIH), University of Michigan, The Norwegian Research Council, and Central Norway Regional Health Authority and the Faculty of Medicine and Health Sciences, Norwegian University of Science and Technology (NTNU). The genotype quality control and imputation has been conducted by the K.G. Jebsen center for genetic epidemiology, Department of public health and nursing, Faculty of medicine and health sciences, Norwegian University of Science and Technology (NTNU). The Rotterdam Study is supported by the Erasmus MC University Medical Center and Erasmus University Rotterdam; The Netherlands Organization for Scientific Research (NWO); The Netherlands Organization for Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly (RIDE); The Netherlands Genomics Initiative (NGI); the Ministry of Education, Culture and Science; the Ministry of Health, Welfare and Sports; the European Commission (DG XII); and the Municipality of Rotterdam.

## Footnotes

*   Author Yakov A. Tsepilov name was revised. 

*   Received November 30, 2020.
*   Revision received December 3, 2020.
*   Accepted December 3, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  1.Häuser W, Ablin J, Fitzcharles M-A, Littlejohn G, Luciano JV, Usui C, et al. Fibromyalgia. Nature Reviews Disease Primers. 2015 2015/08/13; 1(1):15022.
    
    

2.  2.Wolfe F, Clauw DJ, Fitzcharles MA, Goldenberg DL, Hauser W, Katz RL, et al. 2016 Revisions to the 2010/2011 fibromyalgia diagnostic criteria. Semin Arthritis Rheum. 2016 Dec; 46(3):319–329.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.semarthrit.2016.08.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

3.  3.Shipley M. Chronic widespread pain and fibromyalgia syndrome. Medicine. 2018 2018/04/01/; 46(4):252–255.
    
    

4.  4.Fayaz A, Croft P, Langford RM, Donaldson LJ, Jones GT. Prevalence of chronic pain in the UK: a systematic review and meta-analysis of population studies. BMJ Open. 2016 Jun 20; 6(6):e010364.
    
    

5.  5.Mansfield KE, Sim J, Jordan JL, Jordan KP. A systematic review and meta-analysis of the prevalence of chronic widespread pain in the general population. Pain. 2016 Jan; 157(1):55–64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/j.pain.0000000000000314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26270591&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

6.  6.Boonen A, van den Heuvel R, van Tubergen A, Goossens M, Severens JL, van der Heijde D, et al. Large differences in cost of illness and wellbeing between patients with fibromyalgia, chronic low back pain, or ankylosing spondylitis. Annals of the Rheumatic Diseases. 2005; 64(3):396.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjg6IjY0LzMvMzk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

7.  7.Macfarlane GJ, Barnish MS, Jones GT. Persons with chronic widespread pain experience excess mortality: longitudinal results from UK Biobank and meta-analysis. Ann Rheum Dis. 2017 Nov; 76(11):1815–1822.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjEwOiI3Ni8xMS8xODE1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

8.  8.Kvalheim S, Sandven I, Hagen K, Zwart JA. Smoking as a risk factor for chronic musculoskeletal complaints is influenced by age. The HUNT study. Pain. 2013 Jul; 154(7):1073–1079.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2013.03.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23623251&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000320592200016&link_type=ISI) 

9.  9.Kvalheim S, Sandvik L, Winsvold B, Hagen K, Zwart JA. Early menarche and chronic widespread musculoskeletal complaints--Results from the HUNT study. Eur J Pain. 2016 Mar; 20(3):458–464.
    
    

10. 10.Mundal I, Gråwe RW, Bjørngaard JH, Linaker OM, Fors EA. Prevalence and long-term predictors of persistent chronic widespread pain in the general population in an 11-year prospective study: the HUNT study. BMC Musculoskelet Disord. 2014; 15:213.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2474-15-213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24951013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

11. 11.Mundal I, Grawe RW, Bjorngaard JH, Linaker OM, Fors EA. Psychosocial factors and risk of chronic widespread pain: an 11-year follow-up study--the HUNT study. Pain. 2014 Aug; 155(8):1555–1561.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2014.04.033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24813831&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340309800019&link_type=ISI) 

12. 12.Mork PJ, Vasseljen O, Nilsen TI. Association between physical exercise, body mass index, and risk of fibromyalgia: longitudinal data from the Norwegian Nord-Trondelag Health Study. Arthritis Care Res (Hoboken). 2010 May; 62(5):611–617.
    
    

13. 13.Kato K, Sullivan PF, Evengard B, Pedersen NL. Importance of genetic influences on chronic widespread pain. Arthritis Rheum. 2006 May; 54(5):1682–1686.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.21798&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16646040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237533100040&link_type=ISI) 

14. 14.Park DJ, Lee SS. New insights into the genetics of fibromyalgia. Korean J Intern Med. 2017; 32(6):984–995.
    
    

15. 15.Peters MJ, Broer L, Willemen HL, Eiriksdottir G, Hocking LJ, Holliday KL, et al. Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region. Ann Rheum Dis. 2013 Mar; 72(3):427–436.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjg6IjcyLzMvNDI3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

16. 16.Bouhouche A, Benomar A, Bouslam N, Chkili T, Yahyaoui M. Mutation in the epsilon subunit of the cytosolic chaperonin-containing t-complex peptide-1 (Cct5) gene causes autosomal recessive mutilating sensory neuropathy with spastic paraplegia. J Med Genet. 2006 May; 43(5):441–443.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiam1lZGdlbmV0IjtzOjU6InJlc2lkIjtzOjg6IjQzLzUvNDQxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. 17.Lawson VH, Grewal J, Hackshaw KV, Mongiovi PC, Stino AM. Fibromyalgia syndrome and small fiber, early or mild sensory polyneuropathy. Muscle Nerve. 2018; 58(5):625–630.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mus.26131&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29572887&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

18. 18.Vehof J, Zavos HM, Lachance G, Hammond CJ, Williams FM. Shared genetic factors underlie chronic pain syndromes. Pain. 2014 Aug; 155(8):1562–1568.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2014.05.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24879916&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340309800020&link_type=ISI) 

19. 19.Tsepilov YA, Freidin MB, Shadrina AS, Sharapov SZ, Elgaeva EE, van Zundert J, et al. Analysis of genetically independent phenotypes identifies shared genetic factors associated with chronic musculoskeletal pain at different anatomic sites. bioRxiv. 2019:810283.
    
    

20. 20.Gürsoy S, Erdal E, Herken H, Madenci E, Alaşehirli B, Erdal N. Significance of catechol-O-methyltransferase gene polymorphism in fibromyalgia syndrome. Rheumatol Int. 2003 May; 23(3):104–107.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12739038&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183552400002&link_type=ISI) 

21. 21.Hagen K, Pettersen E, Stovner LJ, Skorpen F, Zwart JA. No association between chronic musculoskeletal complaints and Val158Met polymorphism in the Catechol-O-methyltransferase gene. The HUNT study. BMC Musculoskelet Disord. 2006 May 4; 7:40.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2474-7-40&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16674809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

22. 22.Diatchenko L, Nackley AG, Slade GD, Bhalang K, Belfer I, Max MB, et al. Catechol-O-methyltransferase gene polymorphisms are associated with multiple pain-evoking stimuli. Pain. 2006 Dec 5; 125(3):216–224.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2006.05.024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16837133&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242673100006&link_type=ISI) 

23. 23.Smith SB, Parisien M, Bair E, Belfer I, Chabot-Doré A-J, Gris P, et al. Genome-wide association reveals contribution of MRAS to painful temporomandibular disorder in males. Pain. 2019; 160(3):579–591.
    
    

24. 24.Loh P-R, Tucker G, Bulik-Sullivan BK, Vilhjálmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nature Genetics. 2015 2015/03/01; 47(3):284–290.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3190&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

25. 25.Price AL, Zaitlen NA, Reich D, Patterson N. New approaches to population stratification in genome-wide association studies. Nat Rev Genet. 2010 Jul; 11(7):459–463.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2813&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20548291&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278998500008&link_type=ISI) 

26. 26.Yang J, Ferreira T, Morris AP, Medland SE, Madden PAF, Heath AC, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nature Genetics. 2012 2012/04/01; 44(4):369–375.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22426310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

27. 27.Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics; 2010:2190–2191.
    
    

28. 28.Bulik-Sullivan BK, Loh P-R, Finucane HK, Ripke S, Yang J, Patterson N, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature Genetics. 2015 2015/03/01; 47(3):291–295.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

29. 29.Freidin MB, Tsepilov YA, Palmer M, Karssen LC, Suri P, Aulchenko YS, et al. Insight into the genetic architecture of back pain and its risk factors from a study of 509,000 individuals. Pain. 2019 Jun; 160(6):1361–1373.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/j.pain.0000000000001514&link_type=DOI) 

30. 30.Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nature Communications. 2017 2017/11/28; 8(1):1826.
    
    

31. 31.Parisien M, Khoury S, Chabot-Doré A-J, Sotocinal SG, Slade GD, Smith SB, et al. Effect of Human Genetic Variability on Gene Expression in Dorsal Root Ganglia and Association with Pain Phenotypes. Cell Reports. 2017; 19(9):1940–1952.
    
    

32. 32.Consortium TG. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science. 2020; 369(6509):1318.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

33. 33.Heberle H, Meirelles GV, da Silva FR, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015 May 22; 16(1):169.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12859-015-0611-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25994840&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

34. 34.Kamura T, Hara T, Matsumoto M, Ishida N, Okumura F, Hatakeyama S, et al. Cytoplasmic ubiquitin ligase KPC regulates proteolysis of p27(Kip1) at G1 phase. Nat Cell Biol. 2004 Dec; 6(12):1229–1235.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb1194&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15531880&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225453300016&link_type=ISI) 

35. 35.Wang S, Yang YK, Chen T, Zhang H, Yang WW, Song SS, et al. RNF123 has an E3 ligase-independent function in RIG-I-like receptor-mediated antiviral signaling. EMBO Rep. 2016 Aug; 17(8):1155–1168.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiZW1ib3IiO3M6NToicmVzaWQiO3M6OToiMTcvOC8xMTU1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

36. 36.Iida Y, Ciechanover A, Marzese DM, Hata K, Bustos M, Ono S, et al. Epigenetic Regulation of KPC1 Ubiquitin Ligase Affects the NF-κB Pathway in Melanoma. Clin Cancer Res. 2017 Aug 15; 23(16):4831–4842.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6MTA6IjIzLzE2LzQ4MzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMi8wMy8yMDIwLjExLjMwLjIwMjQxMDAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

37. 37.Wang X, Bustos MA, Zhang X, Ramos RI, Tan C, Iida Y, et al. Downregulation of the Ubiquitin-E3 Ligase RNF123 Promotes Upregulation of the NF-κB1 Target SerpinE1 in Aggressive Glioblastoma Tumors. Cancers (Basel). 2020; 12(5):1081.
    
    

38. 38.Liu J, Wan W, Qiao Y, Li Y, Liu Q, Zhu Q, et al. A rare variant of RNF123 in familial ankylosing spondylitis causes dysfunction of osteoclast formation. Rheumatology. 2019; 59(5):1174–1177.
    
    

39. 39.Glahn DC, Curran JE, Winkler AM, Carless MA, Kent JW, Jr.., Charlesworth JC, et al. High dimensional endophenotype ranking in the search for major depression risk genes. Biol Psychiatry. 2012 Jan 1; 71(1):6–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2011.08.022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21982424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298137400005&link_type=ISI) 

40. 40.McBeth J, Silman AJ, Macfarlane GJ. Association of widespread body pain with an increased risk of cancer and reduced cancer survival: a prospective, population-based study. Arthritis Rheum. 2003 Jun; 48(6):1686–1692.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.10973&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12794837&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183410300028&link_type=ISI) 

41. 41.Al-Allaf AW, Mole PA, Paterson CR, Pullar T. Bone health in patients with fibromyalgia. Rheumatology. 2003; 42(10):1202–1206.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/rheumatology/keg356&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12810939&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185596700012&link_type=ISI) 

42. 42.Walker-Bone K, Harvey NC, Ntani G, Tinati T, Jones GT, Smith BH, et al. Chronic widespread bodily pain is increased among individuals with history of fracture: findings from UK Biobank. Arch Osteoporos. 2016; 11(1):1–1.
    
    

43. 43.Alciati A, Sgiarovello P, Atzeni F, Sarzi-Puttini P. Psychiatric problems in fibromyalgia: clinical and neurobiological links between mood disorders and fibromyalgia. Reumatismo. 2012 Sep 28; 64(4):268–274.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23024971&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

44. 44.Micaroni M, Giacchetti G, Plebani R, Xiao GG, Federici L. ATP2C1 gene mutations in Hailey-Hailey disease and possible roles of SPCA1 isoforms in membrane trafficking. Cell Death Dis. 2016 Jun 9; 7(6):e2259.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/cddis.2016.147&link_type=DOI) 

45. 45.Sudbrak R, Brown J, Dobson-Stone C, Carter S, Ramser J, White J, et al. Hailey-Hailey disease is caused by mutations in ATP2C1 encoding a novel Ca(2+) pump. Hum Mol Genet. 2000 Apr 12; 9(7):1131–1140.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/9.7.1131&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10767338&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086696100014&link_type=ISI) 

46. 46.Kollman N, Bass J. Generalized familial benign chronic pemphigus (Hailey-Hailey disease) treated successfully with low-dose naltrexone. JAAD Case Rep. 2018 Aug; 4(7):725–727.
    
    

47. 47.Albers LN, Arbiser JL, Feldman RJ. Treatment of Hailey-Hailey Disease With Low-Dose Naltrexone. JAMA Dermatol. 2017 Oct 1; 153(10):1018–1020.
    
    

48. 48.Park J, Luo ZD. Calcium channel functions in pain processing. Channels (Austin). 2010 Nov-Dec; 4(6):510–517.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/chan.4.6.12869&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21150297&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

49. 49.Bourinet E, Altier C, Hildebrand ME, Trang T, Salter MW, Zamponi GW. Calcium-permeable ion channels in pain signaling. Physiol Rev. 2014 Jan; 94(1):81–140.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00023.2013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24382884&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

50. 50.Younger J, Noor N, McCue R, Mackey S. Low-dose naltrexone for the treatment of fibromyalgia: findings of a small, randomized, double-blind, placebo-controlled, counterbalanced, crossover trial assessing daily pain levels. Arthritis Rheum. 2013 Feb; 65(2):529–538.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.37734&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23359310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314169400030&link_type=ISI) 

51. 51.Derry S, Cording M, Wiffen PJ, Law S, Phillips T, Moore RA. Pregabalin for pain in fibromyalgia in adults. Cochrane Database Syst Rev. 2016 Sep 29; 9(9):Cd011790.
    
    

52. 52.Mogil JS. Pain genetics: past, present and future. Trends Genet. 2012 Jun; 28(6):258–266.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tig.2012.02.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22464640&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305094000002&link_type=ISI) 

53. 53.Kambur O, Männistö PT. Catechol-O-methyltransferase and pain. Int Rev Neurobiol. 2010/11/26 ed; 2010:227–279.
    
    

54. 54.Zhang L, Zhu J, Chen Y, Zhao J. Meta-analysis reveals a lack of association between a common catechol-O-methyltransferase (COMT) polymorphism val□□met and fibromyalgia. Int J Clin Exp Pathol. 2014; 7(12):8489–8497.
    
    

55. 55.Lee C, Liptan G, Kantorovich S, Sharma M, Brenton A. Association of Catechol-O-methyltransferase single nucleotide polymorphisms, ethnicity, and sex in a large cohort of fibromyalgia patients. BMC Rheumatology. 2018 2018/12/12; 2(1):38.
    
    

56. 56.Nicholl BI, Holliday KL, Macfarlane GJ, Thomson W, Davies KA, Neill TW, et al. No evidence for a role of the *catechol-O-methyltransferase* pain sensitivity haplotypes in chronic widespread pain. Annals of the Rheumatic Diseases. 2010; 69(11):2009.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjEwOiI2OS8xMS8yMDA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDMvMjAyMC4xMS4zMC4yMDI0MTAwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

57. 57.Smith SB, Maixner DW, Fillingim RB, Slade G, Gracely RH, Ambrose K, et al. Large candidate gene association study reveals genetic risk factors and therapeutic targets for fibromyalgia. Arthritis & Rheumatism. 2012; 64(2):584–593.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.33338&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21905019&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

58. 58.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 Research & Therapy. 2007 2007/10/26; 9(5):R110.
    
    

59. 59.Hocking LJ, Smith BH, Jones GT, Reid DM, Strachan DP, Macfarlane GJ. Genetic variation in the beta2-adrenergic receptor but not catecholamine-O-methyltransferase predisposes to chronic pain: results from the 1958 British Birth Cohort Study. Pain. 2010 Apr; 149(1):143–151.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2010.01.023&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20167428&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275789500022&link_type=ISI) 

60. 60.Diatchenko L, Slade GD, Nackley AG, Bhalang K, Sigurdsson A, Belfer I, et al. Genetic basis for individual variations in pain perception and the development of a chronic pain condition. Human Molecular Genetics. 2004; 14(1):135–143.
    
    

61. 61.Tammimaki A, Mannisto PT. Catechol-O-methyltransferase gene polymorphism and chronic human pain: a systematic review and meta-analysis. Pharmacogenet Genomics. 2012 Sep; 22(9):673–691.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/FPC.0b013e3283560c46&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22722321&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307652600004&link_type=ISI) 

62. 62.Meng W, Adams MJ, Reel P, Rajendrakumar A, Huang Y, Deary IJ, et al. Genetic correlations between pain phenotypes and depression and neuroticism. European Journal of Human Genetics. 2020 2020/03/01; 28(3):358–366.
    
    

63. 63.Meng W, Chan BW, Harris C, Freidin MB, Hebert HL, Adams MJ, et al. A genome-wide association study finds genetic variants associated with neck or shoulder pain in UK Biobank. Human molecular genetics. 2020; 29(8):1396–1404.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddaa058&link_type=DOI) 

64. 64.Meng W, Adams MJ, Palmer CNA, Agee M, Alipanahi B, Bell RK, et al. Genome-wide association study of knee pain identifies associations with GDF5 and COL27A1 in UK Biobank. Communications Biology. 2019 2019/08/28; 2(1):321.
    
    

65. 65.Howard DM, Adams MJ, Clarke T-K, Hafferty JD, Gibson J, Shirali M, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience. 2019 2019/03/01; 22(3):343–352.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

66. 66.Nagel M, Jansen PR, Stringer S, Watanabe K, de Leeuw CA, Bryois J, et al. Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways. Nature Genetics. 2018 2018/07/01; 50(7):920–927.
    
    

67. 67.Olsen NJ, Park JH. Skeletal muscle abnormalities in patients with fibromyalgia. Am J Med Sci. 1998 Jun; 315(6):351–358.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00000441-199806000-00003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9638891&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

68. 68.Staud R. Autonomic dysfunction in fibromyalgia syndrome: postural orthostatic tachycardia. Curr Rheumatol Rep. 2008 Dec; 10(6):463–466.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11926-008-0076-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19007537&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

69. 69.Furlan R, Colombo S, Perego F, Atzeni F, Diana A, Barbic F, et al. Abnormalities of cardiovascular neural control and reduced orthostatic tolerance in patients with primary fibromyalgia. J Rheumatol. 2005 Sep; 32(9):1787–1793.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoianJoZXVtIjtzOjU6InJlc2lkIjtzOjk6IjMyLzkvMTc4NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzAzLzIwMjAuMTEuMzAuMjAyNDEwMDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

70. 70.Sluka KA, Clauw DJ. Neurobiology of fibromyalgia and chronic widespread pain. Neuroscience. 2016 Dec 3; 338:114–129.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F03%2F2020.11.30.20241000.atom) 

71. 71.Staud R. Peripheral pain mechanisms in chronic widespread pain. Best practice & research Clinical rheumatology. 2011; 25(2):155–164.
    
    

72. 72.Häuser W, Perrot S, Sommer C, Shir Y, Fitzcharles MA. Diagnostic confounders of chronic widespread pain: not always fibromyalgia. Pain Rep. 2017; 2(3).