Uncovering a Genetic Polymorphism Located in Huntingtin Associated Protein 1 in Modulation of Central Pain Sensitization Signaling Pathways

Fibromyalgia syndrome (FMS) is characterized by widespread pain and increased sensitivity to nociceptive stimulus or tenderness. While familial aggregation could suggest a potential hereditary component in FMS development, isolation of genetic determinants has proven difficult due to the multi-factorial nature and complexity of the syndrome. Central sensitization is thought to be one of the key mechanisms leading to FMS in a subset of patients. Enhanced central pain signaling can be measured using the Nociceptive Flexion Reflex (NFR) or RIII threshold. We performed a genome-wide association study (GWAS) using an array to genotype 258,756 human genetic polymorphisms in 225 FMS patients and 77 healthy volunteers and searched for genetic variants associated with a lowered NFR threshold. We have identified a potential association between a single nucleotide polymorphism resulting in a common non-synonymous coding mutation in the Huntingtin associated protein 1 (HAP1) gene (rs4796604, MAF = 0.5) and the NFR threshold (p = 4.78E−06). The Hap1 protein is involved in trafficking and is particularly enriched in neurons. Our results suggest a possible involvement of the neuronal trafficking protein HAP1 in modulating pain signaling pathways and thus participate in the establishment of the NFR threshold.

Introduction

According to the American College of Rheumatology (ACR), the fibromyalgia syndrome (FMS) is defined as a chronic pain condition of musculoskeletal origin but uncertain cause (ACR, 2022). It is a condition characterized by widespread pain and tenderness generally associated with fatigue, sleep disturbance, mood disorders, and cognitive symptoms (ACR, 2022). The first generally accepted diagnostic criteria of fibromyalgia were published only in 1990 (Wolfe et al., 1990), despite the high prevalence of the disease affecting 0.5–5% of the European (Branco et al., 2010) and 2–8% of the United States (Clauw, 2014) population. A wide range of serious co-morbidities, including depression, anxiety, and sleep disorders are frequently associated with fibromyalgia, increasing the utilization of health resources and costs (Arnold et al., 2006; Skaer, 2014).

Although the pathophysiological mechanisms leading to FMS are still poorly understood, both neuronal plasticity and central sensitization are important processes involved in the onset of FMS. Functional magnetic resonance imaging (fMRI) comparing brain activity in FMS patients versus controls revealed an increased neuronal activity in pain-processing brain-regions as well as a diminution in the connections between those regions and the anti-nociceptive regions of the brain (Sluka and Clauw, 2016). Those observations suggest that sensitization of the central nervous system (CNS) results from a physiologic acquired imbalance between incoming pain stimulation and pain inhibitory feedback loops originating from the CNS (Latremoliere and Woolf, 2009). In line with this current working model, the most potent FMS treatments available today are directed at restoring the feedback inhibitory loop with serotonin-norepinephrine reuptake inhibitors (SNRIs) or decreasing excitatory signals (like NMDA glutamate receptor antagonists) (Sluka and Clauw, 2016).

A clear familial aggregation has been observed between FMS patients (Pellegrino et al., 1989; Buskila et al., 1996; Arnold et al., 2004, 2013), with a recent study estimating a sibling recurrence risk ratio reaching 13.6 (Arnold et al., 2013). This aggregation is not observed in rheumatoid arthritis patients under similar conditions (Arnold et al., 2004) and could suggest a possible genetic background to the condition. Moreover, heritability of chronic widespread pain (CWP), which includes FMS, has been estimated to reach 30–55% depending on the study (Kato et al., 2006; Rahman et al., 2021). Although candidate gene approaches focusing on single nucleotide polymorphisms (SNPs) affecting relevant neuro-transmitters activity did uncover some polymorphisms with significantly altered distribution in FMS patients versus controls (Ablin and Buskila, 2015), most genome-wide approaches of FMS underscored the difficulty of isolating single genetic determinants of FMS (Docampo et al., 2014; Ablin and Buskila, 2015).

The Nociceptive Flexion Reflex (NFR) or R-III threshold is triggered by electrical stimulations of variable intensity applied directly to the sural nerve. It is considered as an objective and reliable assessment of central pain pathways since direct stimulation of the nerve bypasses the activation of peripheral nociceptors (Willer, 1990; Sandrini et al., 2005). The correlation between decreased NFR intensities and increased pain sensitivity is well established both in healthy volunteers and chronic pain patients (Willer, 1977; Skljarevski and Ramadan, 2002). Previous studies have shown a decrease in NFR threshold in FMS patients as compared to non-FMS control subjects (Desmeules et al., 2003; Banic et al., 2004; Tanwar et al., 2019).

Although all patients with FMS diagnosis suffer from a chronic pain condition, they represent a heterogenous group with different symptom presentation and one of the current challenges in the field concerns the sub-categorization of FMS patients in the most appropriate therapeutic group. In the current study, we used a genome-wide association study (GWAS) approach to identify genetic determinants of central sensitization using the NFR threshold as objective marker of increased sensitivity in both patients and healthy controls. We uncovered an association between a SNP in the N-terminus of the Huntingtin associated protein 1 (HAP1) gene and NFR threshold levels.

Materials and Methods

Study Population

All data and genetic material were collected during three subsequent clinical trials with patients suffering from FMS conducted by the Division of Clinical Pharmacology and Toxicology of the University Hospital of Geneva (HUG) between 1998 and 2009. The first study (AquaFM) is a case-control trial evaluating different pain assessment tools (including NFR) between FMS patients and a matched control population (Desmeules et al., 2003; Cedraschi et al., 2004). This study was followed by another case-control clinical trial (“PNR”) assessing the psychological and physiological differences between FMS patients and controls in function of the COMTVal158 genotype (Desmeules et al., 2014). Finally, the Milna study is a placebo-controlled trial investigating the pharmaco-dynamic activity of milnacipran in a cohort of FMS patients (Matthey et al., 2013). Each study was approved by the Geneva Ethics Committee (CCER 98-63, 05-052, 05-215), which also approved the reutilization of the collected biological material for the GWAS analysis (CCER 2016-01811). All clinical studies were conducted in accordance with the ethical principles of the Helsinki declaration, the ICH guidelines and Swiss national law governing clinical trials. All participants signed an inform consent form and were free of disease or medical disorder requiring immediate treatment at the time of their enrollment and during the time of the study. All FMS patients were diagnosed according to the 1990 ACR criteria based on tender point examination only. Additional symptoms such as fatigue, sleep disturbance, mood disorders, or cognitive symptoms were integrated in the diagnosis of FMS at a later time point and not considered in this study (Hauser et al., 2017).

Nociceptive Flexion Reflex Measurements

The nociceptive flexion R-III reflex (NFR) was measured as previously described (Desmeules et al., 2003). Briefly, the sural nerve was electrically stimulated with single impulses (0.5 ms) delivered at 6–10 s intervals by a constant current stimulator with variable intensities (1–100 mA). Electromyographic responses were recorded using a pair of surface electrodes placed over the ipsilateral biceps femoris. The NFR response was identified as a multiphasic signal appearing 90–200 ms after stimulation with an AUC greater than 0.150 mV/ms. The NFR threshold (mA) was defined as the stimulation intensity generating 50% positive responses in a series of 30–40 stimulations using Hill’s equation. An upper threshold for input current intensities was set at 100 mA. Only baseline NFR thresholds measured before any study intervention, but after treatment interruption for patients able to stop their medication were considered.

DNA Samples Preparation

DNA purification and concentration of 302 samples (225 FMS patients and 77 controls) was performed in the Laboratory for Clinical Pharmacology and Toxicology of the Geneva University Hospital (HUG) using the QIAamp DNA blood mini kit (QIAGEN). Quantification was performed with the Qubit dsDNA BR Assay (Thermo Fischer Scientific).

Genotyping

Genotyping was performed at the iGE3 genomic platform of the University of Geneva (UniGE) using an Infinium CoreExome-24 BeadChip from Illumina containing 555,356 markers including 4000 custom markers (Mouterde et al., under review¹), among which 906 SNPs selected on the basis of their potential association to FMS and/or central pain sensitization through literature-mining (list of FMS related SNPs in Supplementary Data 1).

Genotype intensities were called using Illumina GenomeStudio™. All sample call rates (CRs) from the genotyping experiment exceeded 99%. The mean genotyping CR for the 302 samples was 99.73%. All SNPs were evaluated by the Illumina cluster quality score and SNPs with GenTrain score < 0.7 were discarded.

In the final data set, SNPs with a genotyping CR below 98%, a significant deviation from Hardy–Weinberg equilibrium (HWE) or a minor allele frequency (MAF) below 5% (CR < 98%, HWE p < 10^–5, MAF < 0.05) were discarded. All genetic coordinates refer to the Ensembl GRCh37 reference assembly.

Genome-Wide Association Study Quality Control and Bioinformatics Analysis

Quality control and genome-wide association analysis were performed using the PLINK 1.07v software (Purcell et al., 2007). Principal component analyses (PCAs) were performed in R version 3.5.3 using the Adegenet package (v2.1.1) (Jombart, 2008; Jombart and Ahmed, 2011). Manhattan and QQ-plots were drawn using the qqman R package (v0.1.4) (Turner, 2014). Lambda statistics was calculated in R version 3.5.3.

Quality control considered relatedness, gender attribution, heterozygosity and inbreeding, genotype consistency versus HapMap phase III population of European descend, allele frequency consistency and ethnicity versus the HRC r1.1 (Human Reference Consortium version r1.1 2016) EUR population data, as well as sample reproducibility (Supplementary Data 2). A final neighborhood analysis confirmed that all outliers had been detected (Guo et al., 2014). From the initial population of 302 participants, 3 were excluded due to first-degree family relationships and 15 based on ethnicity.

Following quality control, 258,756 SNPs (46.59%) in 284 samples including 212 FMS patients and 72 controls were cleared for further analysis. The genotyping rate in the final population is >0.999.

Genome-wide association study imputation of chromosome 17 was performed using the Michigan imputation server (Minimac4) and the EUR population of the HRC r1.1 as reference panel (r² filter = 0.3, phasing with Eagle v2.4). Following imputation, the dataset on chromosome 17 reached 348,128 SNP, which after QC filtration (MAF > 0.05, CR > 0.98, HWE > 1E−05 and imputation R² > 0.8) left 114,238 SNPs for association analysis (versus 7100 before imputation).

Single nucleotide polymorphisms were annotated with the Annovar software using the hg19 build (Wang et al., 2010). Linkage disequilibrium (LD) analysis was performed using the Ensembl tool LD Calculator for the GRCh37 assembly and association results filtered for LD r² < 0.3 (Hunt et al., 2018). Regional association plots for chromosome 17 were prepared using the LocusZoom tool (Boughton et al., 2021). Fine mapping was performed using the FUnctional Mapping and Annotation (FUMA) toolbox (Watanabe et al., 2017).

Molecular Biology Experiments

Expression plasmid containing the wild type C-terminally myc-DDK tagged HAP1 cDNA under the control of the constitutive CMV promoter or the corresponding Hap1^K4R (rs4796604 G > A) variant were transfected in the human neuroblastoma derived cell line, SH-SYS5 (ATCC CRL-2266) cells. Cells and supernatant were harvested 48 h after transfection. For protein analysis, cell lysis supernatant were recovered after 10 min centrifugation at 10,000 × g. Sub-cellular fractionation were performed using the Subcellular Protein Fractionation Kit for Cultured Cells (ThermoFisher Scientific). All protein preparations were separated by SDS-PAGE and analyzed by Western-blotting. mRNA expression levels were analyzed by RT-PCR with normalization to ACT1 mRNA. All molecular biology experiments were done in triplicates and similar transfection efficiency controlled using a co-transfected fluorescent construct. Human HAP1 expression data where extracted from the NCBI Gene Expression Omnibus (GEO) and analyzed using the GEO2R tool for comparison (Edgar et al., 2002; Barrett et al., 2013). Detailed protocols for cell culture experiments are available in Supplementary Methods 3.

Statistical Analysis

Descriptive statistics, group comparisons, multicollinearity variance inflation factor (vif) and linear regressions for sensitivity analyses were done in R version 3.5.3 (Johnston et al., 2018; Fox and Weisberg, 2019). The Mann–Whitney–Wilcoxon (numerical variables) and Chi-square (categorical variables) statistics were used for pairwise comparison. The Fischer’s exact test was used for higher order categorical variables, while correlation between ordinal or numerical variables were tested using the Kendall rank correlation coefficient test. The results of statistical tests with p-values ≤ 0.05 were considered significant. The Bonferroni correction for multiple testing was applied where appropriate. Sensitivity analyses were performed by comparing the changes of the estimates (beta values) resulting from calculating the linear regression parameters either in a subgroup of the study population (by removing patients taking comedication, considering only FMS patients or patients above 40 years old) or removing a variable form the analysis (for no diagnosis) (Thabane et al., 2013). The sensitivity analysis was considered successful is the deviation was below 12%.

Results

Population Data and Genome-Wide Associations Study Set-up

The characteristics of the final study population (212 FMS patients and 72 controls) used for the genomic analysis are described in Table 1. The gender distribution in the patient group is expected as FMS is known to affect primarily women and the proportion of men in the control population is similar to the patient group [5 versus 3%, χ² (1, N = 284) = 0.27, p = 0.60]. All patients were asked to stop medication susceptible to affect the NFR threshold, and while withdrawal was a condition for further participation in the original clinical studies (Desmeules et al., 2003; Cedraschi et al., 2004; Matthey et al., 2013), all participants with a complete NFR recording and DNA sampling were included in the current genomic study. Drugs affecting NFR threshold in patients unable to stop their medication (USM) were identified according to Sandrini et al. (2005). As expected from previous studies (Desmeules et al., 2003; Tanwar et al., 2019), there is a significant difference (Mann–Whitney–Wilcoxon test, p < 0.01; W = 6072) in NFR threshold values between the FMS patient and the control groups (Figure 1).

TABLE 1

Table 1. Genome-wide association study population.

FIGURE 1

Figure 1. Median NFR threshold is slightly lower in FMS patients than controls. Boxplot showing the median and IQR NFR thresholds values for FMS patients and controls included in the final GWAS population. The difference between the two populations is small but statistically significant according to the Mann–Whitney–Wilcoxon test, p-value <0.01; W = 6072.

Due to deviation from normality of the NFR threshold values distribution (Supplementary Data 4), log₁₀ transformed values were used for regression analysis. The co-variates selected for multiple regression analysis were age, gender, FMS diagnosis, cohort of origin, and presence of co-medication susceptible to affect the NFR threshold. There was no relevant correlation between any of those variables and the 10 first axes of the PCA of the genetic data (Supplementary Data 5). An additive genetic model was assumed.

Genome-Wide Association Study Results

The association p-values for SNPs located on autosomal chromosomes and the NFR threshold are shown in Figure 2. The corresponding QQ-plot displayed in Figure 3 and inflation statistics with λ_GC = 0.989, shows that there is no notable evidence of population substructure. The top hits from the linear regression analysis are shown in Supplementary Data 6. No SNPs located on the X chromosome appeared in the top 20 results and the Y chromosome was not investigated due to the low proportion of men in the cohort.

FIGURE 2

Figure 2. Manhattan plot from the linear regression analysis of the GWAS versus log₁₀(NFR thresholds). Each point represents one single nucleotide polymorphism ordered according to its genomic position (x axis). The y axis reflects the probability value [expressed as −log₁₀(value)] obtained from the logistic regression analysis of each individual SNP in the study. The red line marks the genome-wide association p-value threshold level (2E−07). The blue line marks an arbitrary level set at 1E−05.

FIGURE 3

Figure 3. QQ-plot from the linear regression analysis of the GWAS versus log₁₀(NFR thresholds). Plot of observed versus expected distribution of p-values across all SNPs of the GWAS analysis. The predicted p-value is determined as a normally distributed set of probability statistics based on the number of SNPs included in the study (Turner, 2014). The graph shows no major deviation from normality (diagonal) and suggest the absence of confounder in the data. The tendency observed for the lowest p-values indicating that we identified less highly specific targets than might have been expected from a completely random distribution, could reflect the high complexity of the trait as well as the genetic redundancy underlying fundamental neuronal processes. The genetic inflation statistics (λ_GC = 0.989, λ₁₀₀₀ = 0.962) reveals no systematic population stratification bias in our analysis.

None of the candidates from the genome-wide analysis of the cohort reaches the threshold set by the Bonferroni correction (i.e., 0.05/258,756 = 2E−07). However, both the highly polygenic nature of pain signaling, and the low number of samples included in the study could account for the low association statistics. In addition, the Bonferroni correction is often considered too stringent for genome-wide analyses as SNPs are not entirely independent from each other and using the method presented in Li et al. (2012) to evaluate the effective number of tests gives a suggestive p-value threshold of 5.95E−06 for the present study. The top hit of our analysis and only SNP to reach this threshold is rs4796604 (A > G) (p = 4.78E−06) located on chromosome 17.

To increase the SNP coverage around rs4796604, we performed the imputation on chromosome 17 and obtained a dataset of 348,128 SNP. Fine mapping of the association results using the imputed data, confirmed rs4796604 as the lead SNP of our analysis as shown on the regional association plots displayed in Figure 4.

FIGURE 4

Figure 4. Regional visualization plots of GWAS results surrounding lead SNP rs4796604. (A) Association results from the original GWAS analysis. (B) Association results following imputation of chromosome 17. (C) Zoomed in window (39,860,000–39,930,000) for imputed chr17 results.

Different sensitivity analyses were used to test the robustness of the association between rs4796604 and log₁₀(NFR). We found the association to be consistent with respect to comedication (excluding patients taking comedications), FMS diagnosis (using either FMS patients group alone or excluding the diagnosis from the linear regression model) and age (excluding participants below 40 years old).

Gene mapping shows that the SNP rs4796604 corresponds to a non-synonymous variant located inside the HAP1 gene (chr17, MAF of 0.5 in Europeans). HAP1 is mainly expressed in neuronal cells (Li et al., 1995, 1998b; Gutekunst et al., 1998). It has been associated with microtubule-dependent transporters (Twelvetrees et al., 2019), neurotrophic growth factors trafficking and recycling of internalized membrane receptors including TrkA, TrkB, and GABA_AR (Li et al., 2002; Gauthier et al., 2004; Kittler et al., 2004; Rong et al., 2006; Sheng et al., 2006, 2008; Twelvetrees et al., 2010; Wu et al., 2010; Yang et al., 2011; Xiang et al., 2014; Mele et al., 2017; Lim et al., 2018) in humans, as well as catecholamine and glutamate release from mouse chromaffin cells (Mackenzie et al., 2016).

Log₁₀ NFR Thresholds Correlate With rs4796604 Genotype, Fibromyalgia Syndrome Diagnosis, and Presence of Co-medication

We used a multiple linear regression analysis including five co-variates (age, gender, FMS diagnosis, presence of co-medication, and cohort of origin) to assess the GWAS data. The parameters of the model for the association between the rs4796604 genotype and log₁₀ transformed NFR threshold levels are displayed in Table 2. The resulting regression model equation is significant (F_6,277 = 8.047, p = 5.08E−08) with an adjusted R² of 0.130 that is consistent with the complexity of the trait. All variance inflation factors (vif) of the multicollinearity analysis are <1.2, well below accepted thresholds indicating the presence of collinearity in the data.

TABLE 2

Table 2. Multiple linear regression model parameters for the association between log₁₀(NFR threshold) and rs4796604.

The individual parameters of each variable show that the presence of co-medication and FMS diagnosis are important co-variates of the model while neither age nor gender bring significant contributions. It is noteworthy that the cohort of origin does not introduce any significant bias although NFR measurements were spread from 1998 to 2009.

The parameters of the standardized regression model suggest that the three significant parameters share similar weights, with the presence of any additional G allele at position rs4796604 increasing the NFR value by an factor of 0.261. In comparison, the presence of co-medications increases the threshold with a standardized beta of 0.222, while FMS diagnosis decreases the threshold by an standard estimate of 0.208 versus the control population. The NFR threshold values and distribution for the different groups in our study according to rs4796604 genotype, FMS diagnosis and presence of co-medication are presented in Table 3.

TABLE 3

Table 3. NFR threshold distribution according to significant study parameters.

Relationship Between Nociceptive Flexion Reflex Threshold, Fibromyalgia Syndrome Diagnosis, and rs4796604 Genotype

The correlation between the rs4796604 genotype and NFR thresholds is highly significant (r_T = 0.23, p = 6.12E−07, Kendall’s Tau statistics). Moreover, according to the Mann–Whitney U test, the rs4796604 AA genotype group has significantly lower NFR threshold values than the AG and GG groups (Figure 5A).

FIGURE 5

Figure 5. The rs4796604 A allele correlates with lower NFR threshold values. (A) Median NFR threshold distribution in function of rs4796604 genotypes. Mann–Whitney statistics were used to assess the differences between genotypes (with ***p < 0.001). There is no correlation between rs4796604 genotype and FMS diagnosis (B), nor presence of co-medication (C). n.s. = not significant.

In addition, although the correlation between the rs4796604 genotype and the NFR threshold is visible in every subgroup, there are no significant differences between the subgroups (i.e., FMS versus controls or with and without co-medication) based on rs4796604 genotype except, in the most severely affected AA group, where FMS patients have significantly lower NFR thresholds than the controls (p < 0.01) (Supplementary Data 7).

Thus, our genome-wide analysis identified the SNP rs4796604 as a putative candidate for the prediction of NFR threshold levels. However, there is no correlation between this polymorphism and FMS diagnosis or the presence of co-medication using the Fisher exact test (Figures 5B,C). We found no correlation between rs4796604 genotype and any other clinical scale available from the original studies records (Desmeules et al., 2003; Cedraschi et al., 2004; Matthey et al., 2013), such as pain or FMS severity as evaluated through quantitative sensory testing (QST) (hot and cold temperature pain thresholds, cold pressor test), or questionnaires assessing the impact of FMS on the patient’s quality-of-life (Fibromyalgia Impact Questionnaire – FIQ), Psychological General Well-Being (PGWB), or functional status (Physical Component Score of the SF-36). However, as previously reported in the literature (Neziri et al., 2011), although all mentioned QST results, except cold and warm perception thresholds, showed a highly significant difference between case and controls, they are only poorly correlated with NFR threshold values. The summary of the QST results in our cohort and their correlation with NFR values is available in Supplementary Material 8.

Effect of the rs4796604 Mutation on Huntingtin Associated Protein 1 (HAP1) Transcription, Protein Expression, and Sub-cellular Distribution

The rs4796604 A > G polymorphism is located within the first exon of the HAP1 gene and results in a lysine (K) to arginine (R) amino-acid substitution in position 4 of the encoded protein. The SNP rs4796604 is located at the very N-terminus of the protein and in close proximity to a N-myristylation site. Myristoyl is a lipid tail used for membrane or protein interaction which is often involved in sub-cellular localization of the attached protein (Martin et al., 2011). Thus, we decided to investigated the difference in protein levels and the cellular distribution of the mutant versus the wild-type protein. We therefore transfected the neuroblastoma-derived SH-SYS5 cells with plasmids constitutively overexpressing a C-terminal DDK-myc tagged version of the Hap1^K4R mutant and corresponding wild-type protein. Cells were harvested 48 h post-transfection and both the mRNA and protein content were analyzed.

RT-PCR analysis of mRNA extracts showed that the expression levels of both constructs are similar in the chosen expression system. In parallel, we compared HAP1 mRNA expression levels in function of the rs4796604 genotype in a human reference population of European descent (CEU) using global expression dataset available from the GEO public repository [(Huang et al., 2007), accession number: GSE7761] and the GEO2R tool. There were no differences in HAP1 mRNA expression based on rs4796604 polymorphism among the 82 participants included (22 AA, 44 AG, and 16 GG with a MAF of 0.54, p = 0.33).

Total protein expression and subcellular distribution were assessed using the C-terminal Flag-tag of the constructs. As shown in Supplementary Data 9, and in line with the expression data, the Hap1^K4R mutation did not affect the overall protein amount when compared to actin expression levels. As expected from previous studies, the major portion of the wild-type protein is found associated with secretory vesicles and/or cytoskeletal fraction (Li et al., 1996, 1998a; Engelender et al., 1997; Gutekunst et al., 1998). The Hap1^K4R mutation has no visible effect on subcellular protein distribution under the current experimental conditions (Supplementary Data 9).

Discussion

Our GWAS identified rs4796604, a polymorphism resulting in an amino-acid substitution in the HAP1 gene, as a potential new genetic determinant of the NFR threshold, a measure of enhanced central pain signaling. The linear regression model parameter obtained for rs4796604 shows that FMS diagnosis and the presence of concomitant medication are important co-factors, while neither age nor gender have a significant effect. The influence of FMS diagnosis is in line with previous studies reporting a lowering of the threshold in FMS patients compared to a control population (Desmeules et al., 2003; Banic et al., 2004; Tanwar et al., 2019). Similarly, many pharmacological compounds affecting neurotransmission in the CNS are known to have a strong influence on NFR thresholds (Skljarevski and Ramadan, 2002; Sandrini et al., 2005). Although gender influence on NFR thresholds has been previously reported, the low proportion of men in our cohort might mask this effect.

Central sensitization results from an acquired imbalance in pain-transmitting and pain-attenuating signals in the CNS. Many different neuro-transmitters (NT) such as catecholamines, glutamate and GABA are involved in transmission and down-regulation of nociceptive stimuli in nerve terminals (reviewed by Harte et al., 2018). Moreover, neuronal growth factors as BDNF or nerve growth factor (NGF) play an important role not only in the development but also in the maintenance of different neuronal cell populations that are important for the regulation of nociceptive stimuli (Dai and Zhou, 2014). Although the precise molecular function of Hap1 in neuronal cells remains unclear, this microtubule-associated protein is involved in intracellular trafficking in neurons and participates both in neurotransmitter (NT) release (Mackenzie et al., 2016) and plasma membrane receptor recycling via endocytic compartments (Li et al., 2002; Kittler et al., 2004; Rong et al., 2006; Sheng et al., 2008; Twelvetrees et al., 2010; Xiang et al., 2014) possibly by acting as a scaffold protein (Rong et al., 2007). HAP1 expression is particularly enriched in the hypothalamus, although HAP1 mRNA is detected in all brain tissues and the spinal cord.²

HAP1 was first identified as an interactor of neurodegenerative Huntingtin (Htt), the protein responsible for Huntington disease and was subsequently shown to protect specific neuronal cells from degradation (Li et al., 2003; Mele et al., 2017; Wroblewski et al., 2018). Moreover, HAP1 has many connections to nociceptive signal transduction. First, HAP1 plays a role in the NGF-receptor TrkA recycling to the plasma membrane (Rong et al., 2006). NGF promotes neurite outgrowth in developing cells and is required for survival of sensory neurons (Marmigere et al., 2006). Interestingly, intramuscular injections of NGF are able to induce neuronal sensitization in rats (Hoheisel et al., 2007) and some mutations in TrkA result in congenital insensitivity to pain with anhidrosis (CIPA), an autosomal recessive disorder characterized by the inability to feel pain and absence of sweating (Indo et al., 1996). HAP1 also participates in regulating neuronal excitability through the recycling of GABA_A receptors to postsynaptic membranes (Kittler et al., 2004; Sheng et al., 2006; Twelvetrees et al., 2010; Mele et al., 2017). Inhibitory mechanisms regulated by GABA and GABA-receptors in the spinal dorsal horn play a central role in modulating pain signals transmitted to the brain (Prescott, 2015), and the strength of the inhibitory signal at synapses correlates with the amount of GABA_A receptors present on the membrane (Jacob et al., 2008). Pharmacological modulation of GABA_A receptor mediated signaling is considered an interesting alternative for pain management in chronic pain conditions. However, the sedative properties of currently available GABA_A receptor agonists limit their use as analgesics (Fischer, 2017). In a slightly different register, HAP1 was also shown to participate directly in exocytosis and regulate plasma membrane docking of secretory vesicle from mouse neuroendocrine cells thereby regulating release of catecholamines and glutamate from synapses (Mackenzie et al., 2016). While glutamate is the main excitatory NT in human CNS acting through NDMA-receptors to transmit peripheral nociceptive signals to the brain, catecholamines, as serotonin or norepinephrine, are the mediators of descending pain modulatory signals transmission playing an important role in pain attenuation (Pereira and Goudet, 2018; Bravo et al., 2019). Finally, Hap1 is involved both in proBDNF anterograde transport, and endocytosis of the BDNF-TrkB receptor complex (Gauthier et al., 2004; Sheng et al., 2008; Wu et al., 2010; Yang et al., 2011; Xiang et al., 2014; Lim et al., 2018). BDNF is a neurotrophic growth factor promoting neuronal growth and survival, as well as synaptic plasticity influencing memory, learning and cognition processes (Sasi et al., 2017). BDNF acts as an interesting modulator of nociception and has been involved in central sensitization, although its pleiotropic roles at different levels of the CNS are hard to interpret (Merighi et al., 2008; Sasi et al., 2017). Indeed, internalization of the BDNF-TrkB complex triggers a variety of signaling cascades modulating both excitatory (glutamine) and inhibitory (GABA) signaling pathways (Sasi et al., 2017). Thus, many nociceptive signal transduction pathways are connected to HAP1 function. In addition, in a completely different approach, the authors of a recent immunolabeling study showed an important overlap between HAP1 and different neurochemical markers in the mouse dorsal root ganglion (DRG) neurons and suggested a role for HAP1 in pain transduction (Islam et al., 2020).

Although Hap1 rs4796604 genotype correlates with levels of the NFR threshold, we found no direct link between Hap1 and FMS diagnosis according to ACR 1990 criteria. However, many of the processes regulated by Hap1 have been linked with FMS and are believed to participate in the onset or the progression of the condition, suggesting that understanding Hap1 cellular activity could help specify the nature of the central sensitization process in different FMS patients. Indeed, reflecting the broad range of functions mediated by neurotrophic growth factors like BDNF and NGF, their implication in FMS has been postulated at many different levels, from regulation of neurotransmission to synaptic plasticity affecting learning and memory processes. However, while the involvement in pain signaling of those factors is well documented, if not well understood, their implication in FMS is less clear. Indeed, differences in BDNF and NGF circulating and cerebrospinal fluid (CSF) levels between FMS patients and healthy controls are still controversial although most studies report increased concentrations (Baumeister et al., 2019; Brietzke et al., 2019). However, the expression pattern of neurotrophic factors are not uniform across different regions of the CNS (Notaras and van den Buuse, 2020), thus the physiological implications of increased peripheral or even central BDNF or NGF levels remain unclear. On the genetic level, a study associated different BDNF related polymorphisms with FMS in the Korean population (Park et al., 2018). A recent study showed increased GABA_AR levels in cortical neurons of FMS patients (Pomares et al., 2020). Regarding NTs in FMS, both serotonin and norepinephrine levels were found to be significantly decreased in CSF of FMS patients (Russell et al., 1992), while glutamate levels are increased in specific brain areas of FMS patients (Foerster et al., 2015). In addition, the most common pharmacological treatments used in the context of FMS are inhibitors of the NMDA glutamate receptor or of serotonin and/or norepinephrine reuptake (Sluka and Clauw, 2016; Pereira and Goudet, 2018; Bravo et al., 2019).

Despite the observed heritability of pain-related diseases, the search for polymorphisms associated with chronic pain conditions in general and FMS in particular is a long-standing issue. Candidate gene approaches focused their attention on SNPs affecting neurotransmitters activity known to be relevant for pain signal transduction and processing as well as on endogenous pain modulating systems (reviewed by Knezevic et al., 2018; Janssen et al., 2021). Frequently reported genes include serotonin receptors and transporter (HTRs, SLC6A4), dopamine receptor and degradation (COMT, TAAR1), glutamate receptor (GRIA4), adrenergic receptors (ADRs) as well as the mu opioid receptor OPRM1, the NGF BDNF or sodium channels (TRVPs for instance). In addition to candidate approaches, several genome-wide association studies directed toward elucidation of the genetic background of FMS, chronic widespread pain (CWP) or related chronic pain conditions uncovered additional candidate genes and polymorphisms (Docampo et al., 2014; Johnston et al., 2019; Rahman et al., 2021). None of the polymorphism previously associated with FMS or other chronic pain conditions nor any of the more recent GWAS candidates did show significant association with NFR in our study (data not shown). This is not unexpected as the NFR threshold, while being intimately linked to pain signal transmission is not a proxy for FMS nor CWP.

Pain perception is a highly polygenic trait that is very difficult to evaluate accurately, and although measurement of NFR thresholds constitutes one of the most objective assessment of the central pain signaling network reactivity available, a number of modulating factors ranging from pharmacological interactions to mental and psychological states have been identified (Sandrini et al., 2005). Thus, taking all those considerations into account, the number of participants in the current study is low and could easily account for the poor statistical power of our GWAS results, which is clearly the major limitation of the current study. Additional limitations to the current study, include the recognized inter- and even intra-individual variability of NFR measurements that are affected by a number of environmental and psychological factors in addition to the physical parameters as well as the low proportion of men in our cohort or the use of the old ACR 1990 criteria for FMS diagnosis.

The replication of our findings in independent cohorts will be critical to confirm their validity. Considering the measured effect size, a MAF of 0.5 in a European population and a significance threshold of p = 0.05, the replication of rs4796604 association with NFR would require a sample size of about 200 participants (power calculations heatmaps according to Visscher et al. (2017) are presented in Supplementary Data 10). Moreover, characterization of the functional implication of the Hap1^Lys4Arg mutation at the molecular level will be central to unravel the importance of this variant for transmission of painful stimuli and understand whether there is a physiological relevant difference in nociceptive signal transmission between individuals expressing different forms of the protein. Although the position of the mutation at the very beginning of a N-myristoylated protein could potentially interfere with Hap1’s membrane attachment, we could not observe any differences in expression nor subcellular repartition between the wild type and mutant protein under current experimental conditions. However, N-myristoylation is often associated with dynamic membrane attachment which might not be apparent in a steady-state distribution. Alternatively, the high expression levels resulting from constitutive expression of the exogenous transcript could mask or disrupt the normal cellular pattern.

In summary, using a GWAS, we identified a new potential genetic determinant of NFR threshold level. While Hap1 is not a predictor of FMS, this protein is involved in the regulation of several signaling pathways involved in the pathophysiology of the disease, and differential regulation of synaptic transmission through modulation of intracellular vesicular transport at synapses might help to better characterize the enhanced sensitization in individual patients. Further studies will be directed at replicating the current results in independent cohorts, as well as unraveling of the functional implications of the Hap1^K4R mutation in Hap1 protein function.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://doi.org/10.26037/yareta:l35qwdy5t5b5vll3ae3hmzigy4.

Ethics Statement

The studies involving human participants were reviewed and approved by the Geneva Ethics Committee – CCER (Commission Cantonale d’Ethique de la Recherche) de la République et du Canton de Genève, Geneva, Switzerland. All patients/participants provided their written informed consent to participate in this study.

Author Contributions

YG, JD, AM, and EP designed the study. YG, KS, and MM performed the experiments. YG, KS, and GE performed the bioinformatic analysis. YG wrote the manuscripts’ draft and prepared the figures and tables. JD, MM, EK, GP, EP, CC, and GE reviewed the draft. All authors read and approved the final manuscript.

Funding

This work was supported by the Swiss National Science Foundation (grant number 320030_182686). The Swiss National Science Foundation (grant number 320030_159669) supported EP and MM was supported by the iGE3 Ph.D. salary award.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank G. Mottet-Osman and F. Doffey-Lazeyras for expert technical assistance with cell culture and transfection. T. Vujic, V. Lazarevic, and M. Girard for providing laboratory equipment as well as E. Tseligka and C. Tapparel for the SH-SYS5 cell line.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.807773/full#supplementary-material

Footnotes

1. ^ Mouterde, M., Daali, Y., Rollason, V., Čížková, M., Mulugeta, A., Al-Balushi, K. K., et al. (under review) Joint analysis of phenotypic and genomic diversity sheds light on the evolution of xenobiotic metabolism in humans. Genome Biology and Evolution.
2. ^ https://gtexportal.org

References

Ablin, J. N., and Buskila, D. (2015). Update on the genetics of the fibromyalgia syndrome. Best Pract. Res. Clin. Rheumatol. 29, 20–28. doi: 10.1016/j.berh.2015.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

ACR (2022). Glossary. Available online at: https://www.rheumatology.org/Learning-Center/Glossary/articleType/CategoryView/categoryId/11 (Accessed May 17 2022).

Google Scholar

Arnold, L. M., Fan, J., Russell, I. J., Yunus, M. B., Khan, M. A., Kushner, I., et al. (2013). The fibromyalgia family study: a genome-wide linkage scan study. Arthritis Rheum. 65, 1122–1128. doi: 10.1002/art.37842

PubMed Abstract | CrossRef Full Text | Google Scholar

Arnold, L. M., Hudson, J. I., Hess, E. V., Ware, A. E., Fritz, D. A., Auchenbach, M. B., et al. (2004). Family study of fibromyalgia. Arthritis Rheum. 50, 944–952. doi: 10.1002/art.20042

PubMed Abstract | CrossRef Full Text | Google Scholar

Arnold, L. M., Hudson, J. I., Keck, P. E., Auchenbach, M. B., Javaras, K. N., and Hess, E. V. (2006). Comorbidity of fibromyalgia and psychiatric disorders. J. Clin. Psychiatry 67, 1219–1225.

Google Scholar

Banic, B., Petersen-Felix, S., Andersen, O. K., Radanov, B. P., Villiger, P. M., Arendt-Nielsen, L., et al. (2004). Evidence for spinal cord hypersensitivity in chronic pain after whiplash injury and in fibromyalgia. Pain 107, 7–15. doi: 10.1016/j.pain.2003.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 41, D991–D995. doi: 10.1093/nar/gks1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Baumeister, D., Eich, W., Saft, S., Geisel, O., Hellweg, R., Finn, A., et al. (2019). No evidence for altered plasma NGF and BDNF levels in fibromyalgia patients. Sci. Rep. 9:13667. doi: 10.1038/s41598-019-49403-49407

CrossRef Full Text | Google Scholar

Boughton, A. P., Welch, R. P., Flickinger, M., VandeHaar, P., Taliun, D., Abecasis, G. R., et al. (2021). LocusZoom.js: interactive and embeddable visualization of genetic association study results. Bioinformatics 37, 3017–3018. doi: 10.1093/bioinformatics/btab186

PubMed Abstract | CrossRef Full Text | Google Scholar

Branco, J. C., Bannwarth, B., Failde, I., Abello Carbonell, J., Blotman, F., Spaeth, M., et al. (2010). Prevalence of fibromyalgia: a survey in five European countries. Semin. Arthritis Rheum. 39, 448–453. doi: 10.1016/j.semarthrit.2008.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo, L., Llorca-Torralba, M., Berrocoso, E., and Mico, J. A. (2019). Monoamines as drug targets in chronic pain: focusing on neuropathic pain. Front. Neurosci. 13:1268. doi: 10.3389/fnins.2019.01268

PubMed Abstract | CrossRef Full Text | Google Scholar

Brietzke, A. P., Antunes, L. C., Carvalho, F., Elkifury, J., Gasparin, A., Sanches, P. R. S., et al. (2019). Potency of descending pain modulatory system is linked with peripheral sensory dysfunction in fibromyalgia: an exploratory study. Medicine (Baltimore) 98:e13477. doi: 10.1097/MD.0000000000013477

PubMed Abstract | CrossRef Full Text | Google Scholar

Buskila, D., Neumann, L., Hazanov, I., and Carmi, R. (1996). Familial aggregation in the fibromyalgia syndrome. Semin. Arthritis Rheum 26, 605–611. doi: 10.1016/s0049-0172(96)80011-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Cedraschi, C., Desmeules, J., Rapiti, E., Baumgartner, E., Cohen, P., Finckh, A., et al. (2004). Fibromyalgia: a randomised, controlled trial of a treatment programme based on self management. Ann. Rheum. Dis. 63, 290–296. doi: 10.1136/ard.2002.004945

PubMed Abstract | CrossRef Full Text | Google Scholar

Clauw, D. J. (2014). Fibromyalgia: a clinical review. JAMA 311, 1547–1555. doi: 10.1001/jama.2014.3266

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, R. P., and Zhou, X. F. (2014). Handbook of Neurotoxicity: Neurotrophins and Pain. New York, NY: Springer.

Google Scholar

Desmeules, J. A., Cedraschi, C., Rapiti, E., Baumgartner, E., Finckh, A., Cohen, P., et al. (2003). Neurophysiologic evidence for a central sensitization in patients with fibromyalgia. Arthritis Rheum. 48, 1420–1429. doi: 10.1002/art.10893

PubMed Abstract | CrossRef Full Text | Google Scholar

Desmeules, J., Chabert, J., Rebsamen, M., Rapiti, E., Piguet, V., Besson, M., et al. (2014). Central pain sensitization, COMT Val158Met polymorphism, and emotional factors in fibromyalgia. J. Pain 15, 129–135. doi: 10.1016/j.jpain.2013.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Docampo, E., Escaramis, G., Gratacos, M., Villatoro, S., Puig, A., Kogevinas, M., et al. (2014). Genome-wide analysis of single nucleotide polymorphisms and copy number variants in fibromyalgia suggest a role for the central nervous system. Pain 155, 1102–1109. doi: 10.1016/j.pain.2014.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30, 207–210. doi: 10.1093/nar/30.1.207

PubMed Abstract | CrossRef Full Text | Google Scholar

Engelender, S., Sharp, A. H., Colomer, V., Tokito, M. K., Lanahan, A., Worley, P., et al. (1997). Huntingtin-associated protein 1 (HAP1) interacts with the p150Glued subunit of dynactin. Hum. Mol. Genet. 6, 2205–2212. doi: 10.1093/hmg/6.13.2205

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischer, B. D. (2017). GABA a receptors as targets for the management of pain-related disorders: historical perspective and update. CNS Neurol. Disord. Drug Targets 16, 658–663. doi: 10.2174/1871527316666170207155149

PubMed Abstract | CrossRef Full Text | Google Scholar

Foerster, B. R., Nascimento, T. D., DeBoer, M., Bender, M. A., Rice, I. C., Truong, D. Q., et al. (2015). Excitatory and inhibitory brain metabolites as targets of motor cortex transcranial direct current stimulation therapy and predictors of its efficacy in fibromyalgia. Arthritis Rheumatol. 67, 576–581. doi: 10.1002/art.38945

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, J., and Weisberg, S. (2019). An R Companion to Applied Regression. Thousand Oaks CA: Sage.

Google Scholar

Gauthier, L. R., Charrin, B. C., Borrell-Pages, M., Dompierre, J. P., Rangone, H., Cordelieres, F. P., et al. (2004). Huntingtin controls neurotrophic support and survival of neurons by enhancing BDNF vesicular transport along microtubules. Cell 118, 127–138. doi: 10.1016/j.cell.2004.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., He, J., Zhao, S., Wu, H., Zhong, X., Sheng, Q., et al. (2014). Illumina human exome genotyping array clustering and quality control. Nat. Protoc. 9, 2643–2662. doi: 10.1038/nprot.2014.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutekunst, C. A., Li, S. H., Yi, H., Ferrante, R. J., Li, X. J., and Hersch, S. M. (1998). The cellular and subcellular localization of huntingtin-associated protein 1 (HAP1): comparison with huntingtin in rat and human. J. Neurosci. 18, 7674–7686. doi: 10.1523/JNEUROSCI.18-19-07674.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Harte, S. E., Harris, R. E., and Clauw, D. J. (2018). The neurobiology of central sensitization. J. Appl. Behav. Res. 23:e12137. doi: 10.1111/jabr.12137

CrossRef Full Text | Google Scholar

Hauser, W., Ablin, J., Perrot, S., and Fitzcharles, M. A. (2017). Management of fibromyalgia: practical guides from recent evidence-based guidelines. Pol. Arch. Intern. Med. 127, 47–56. doi: 10.20452/pamw.3877

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoheisel, U., Unger, T., and Mense, S. (2007). Sensitization of rat dorsal horn neurons by NGF-induced subthreshold potentials and low-frequency activation. a study employing intracellular recordings in vivo. Brain Res. 1169, 34–43. doi: 10.1016/j.brainres.2007.06.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, R. S., Duan, S., Shukla, S. J., Kistner, E. O., Clark, T. A., Chen, T. X., et al. (2007). Identification of genetic variants contributing to cisplatin-induced cytotoxicity by use of a genomewide approach. Am. J. Hum. Genet. 81, 427–437. doi: 10.1086/519850

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunt, S. E., McLaren, W., Gil, L., Thormann, A., Schuilenburg, H., Sheppard, D., et al. (2018). Ensembl variation resources. Database 2018:bay119. doi: 10.1093/database/bay119

PubMed Abstract | CrossRef Full Text | Google Scholar

Indo, Y., Tsuruta, M., Hayashida, Y., Karim, M. A., Ohta, K., Kawano, T., et al. (1996). Mutations in the TRKA/NGF receptor gene in patients with congenital insensitivity to pain with anhidrosis. Nat. Genet. 13, 485–488. doi: 10.1038/ng0896-485

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, M. N., Maeda, N., Miyasato, E., Jahan, M. R., Tarif, A. M. M., Ishino, T., et al. (2020). Expression of huntingtin-associated protein 1 in adult mouse dorsal root ganglia and its neurochemical characterization in reference to sensory neuron subpopulations. IBRO Rep. 9, 258–269. doi: 10.1016/j.ibror.2020.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacob, T. C., Moss, S. J., and Jurd, R. (2008). GABA(A) receptor trafficking and its role in the dynamic modulation of neuronal inhibition. Nat. Rev. Neurosci. 9, 331–343. doi: 10.1038/nrn2370

PubMed Abstract | CrossRef Full Text | Google Scholar

Janssen, L. P., Medeiros, L. F., Souza, A., and Silva, J. D. (2021). Fibromyalgia: a review of related polymorphisms and clinical relevance. An. Acad. Bras. Cienc. 93, (Suppl. 4):e20210618. doi: 10.1590/0001-3765202120210618

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnston, K. J. A., Adams, M. J., Nicholl, B. I., Ward, J., Strawbridge, R. J., Ferguson, A., et al. (2019). Genome-wide association study of multisite chronic pain in UK biobank. PLoS Genet 15:e1008164. doi: 10.1371/journal.pgen.1008164

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnston, R., Jones, K., and Manley, D. (2018). Confounding and collinearity in regression analysis: a cautionary tale and an alternative procedure, illustrated by studies of British voting behaviour. Qual. Quant. 52, 1957–1976. doi: 10.1007/s11135-017-0584-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Jombart, T. (2008). adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405. doi: 10.1093/bioinformatics/btn129

PubMed Abstract | CrossRef Full Text | Google Scholar

Jombart, T., and Ahmed, I. (2011). adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27, 3070–3071. doi: 10.1093/bioinformatics/btr521

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, K., Sullivan, P. F., Evengard, B., and Pedersen, N. L. (2006). Importance of genetic influences on chronic widespread pain. Arthritis Rheum. 54, 1682–1686. doi: 10.1002/art.21798

PubMed Abstract | CrossRef Full Text | Google Scholar

Kittler, J. T., Thomas, P., Tretter, V., Bogdanov, Y. D., Haucke, V., Smart, T. G., et al. (2004). Huntingtin-associated protein 1 regulates inhibitory synaptic transmission by modulating gamma-aminobutyric acid type a receptor membrane trafficking. Proc. Natl. Acad. Sci. U S A. 101, 12736–12741. doi: 10.1073/pnas.0401860101

PubMed Abstract | CrossRef Full Text | Google Scholar

Knezevic, N. N., Tverdohleb, T., Knezevic, I., and Candido, K. D. (2018). The role of genetic polymorphisms in chronic pain patients. Int. J. Mol. Sci. 19:1707. doi: 10.3390/ijms19061707

PubMed Abstract | CrossRef Full Text | Google Scholar

Latremoliere, A., and Woolf, C. J. (2009). Central sensitization: a generator of pain hypersensitivity by central neural plasticity. J. Pain 10, 895–926. doi: 10.1016/j.jpain.2009.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M. X., Yeung, J. M., Cherny, S. S., and Sham, P. C. (2012). Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum. Genet. 131, 747–756. doi: 10.1007/s00439-011-1118-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. H., Gutekunst, C. A., Hersch, S. M., and Li, X. J. (1998a). Interaction of huntingtin-associated protein with dynactin P150Glued. J. Neurosci. 18, 1261–1269. doi: 10.1523/JNEUROSCI.18-04-01261.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. H., Hosseini, S. H., Gutekunst, C. A., Hersch, S. M., Ferrante, R. J., and Li, X. J. (1998b). A human HAP1 homologue. cloning, expression, and interaction with huntingtin. J. Biol. Chem. 273, 19220–19227. doi: 10.1074/jbc.273.30.19220

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. H., Yu, Z. X., Li, C. L., Nguyen, H. P., Zhou, Y. X., Deng, C., et al. (2003). Lack of huntingtin-associated protein-1 causes neuronal death resembling hypothalamic degeneration in Huntington’s disease. J. Neurosci. 23, 6956–6964. doi: 10.1523/JNEUROSCI.23-17-06956.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. J., Li, S. H., Sharp, A. H., Nucifora, F. C. Jr., Schilling, G., Lanahan, A., et al. (1995). A huntingtin-associated protein enriched in brain with implications for pathology. Nature 378, 398–402. doi: 10.1038/378398a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. J., Sharp, A. H., Li, S. H., Dawson, T. M., Snyder, S. H., and Ross, C. A. (1996). Huntingtin-associated protein (HAP1): discrete neuronal localizations in the brain resemble those of neuronal nitric oxide synthase. Proc. Natl. Acad. Sci. U S A. 93, 4839–4844. doi: 10.1073/pnas.93.10.4839

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Chin, L. S., Levey, A. I., and Li, L. (2002). Huntingtin-associated protein 1 interacts with hepatocyte growth factor-regulated tyrosine kinase substrate and functions in endosomal trafficking. J. Biol. Chem. 277, 28212–28221. doi: 10.1074/jbc.M111612200

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, Y., Wu, L. L., Chen, S., Sun, Y., Vijayaraj, S. L., Yang, M., et al. (2018). HAP1 Is Required for Endocytosis and Signalling of BDNF and Its Receptors in Neurons. Mol. Neurobiol. 55, 1815–1830. doi: 10.1007/s12035-016-0379-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackenzie, K. D., Lumsden, A. L., Guo, F., Duffield, M. D., Chataway, T., Lim, Y., et al. (2016). Huntingtin-associated protein-1 is a synapsin I-binding protein regulating synaptic vesicle exocytosis and synapsin I trafficking. J. Neurochem. 138, 710–721. doi: 10.1111/jnc.13703

PubMed Abstract | CrossRef Full Text | Google Scholar

Marmigere, F., Montelius, A., Wegner, M., Groner, Y., Reichardt, L. F., and Ernfors, P. (2006). The Runx1/AML1 transcription factor selectively regulates development and survival of TrkA nociceptive sensory neurons. Nat. Neurosci. 9, 180–187. doi: 10.1038/nn1631

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, D. D., Beauchamp, E., and Berthiaume, L. G. (2011). Post-translational myristoylation: fat matters in cellular life and death. Biochimie 93, 18–31. doi: 10.1016/j.biochi.2010.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthey, A., Cedraschi, C., Piguet, V., Besson, M., Chabert, J., Daali, Y., et al. (2013). Dual reuptake inhibitor milnacipran and spinal pain pathways in fibromyalgia patients: a randomized, double-blind, placebo-controlled trial. Pain Physician 16, E553–E562.

PubMed Abstract | Google Scholar

Mele, M., Aspromonte, M. C., and Duarte, C. B. (2017). Downregulation of GABAA Receptor Recycling Mediated by HAP1 Contributes to Neuronal Death in In Vitro Brain Ischemia. Mol. Neurobiol. 54, 45–57. doi: 10.1007/s12035-015-9661-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Merighi, A., Salio, C., Ghirri, A., Lossi, L., Ferrini, F., Betelli, C., et al. (2008). BDNF as a pain modulator. Prog. Neurobiol. 85, 297–317. doi: 10.1016/j.pneurobio.2008.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Neziri, A. Y., Curatolo, M., Nuesch, E., Scaramozzino, P., Andersen, O. K., Arendt-Nielsen, L., et al. (2011). Factor analysis of responses to thermal, electrical, and mechanical painful stimuli supports the importance of multi-modal pain assessment. Pain 152, 1146–1155. doi: 10.1016/j.pain.2011.01.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Notaras, M., and van den Buuse, M. (2020). Neurobiology of BDNF in fear memory, sensitivity to stress, and stress-related disorders. Mol. Psychiatry 25, 2251–2274. doi: 10.1038/s41380-019-0639-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, D. J., Kim, S. H., Nah, S. S., Lee, J. H., Kim, S. K., Lee, Y. A., et al. (2018). Association between brain-derived neurotrophic factor gene polymorphisms and fibromyalgia in a Korean population: a multicenter study. Arthritis Res. Ther. 20:220. doi: 10.1186/s13075-018-1726-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Pellegrino, M. J., Waylonis, G. W., and Sommer, A. (1989). Familial occurrence of primary fibromyalgia. Arch. Phys. Med. Rehabil. 70, 61–63.

PubMed Abstract | Google Scholar

Pereira, V., and Goudet, C. (2018). Emerging Trends in Pain Modulation by Metabotropic Glutamate Receptors. Front. Mol. Neurosci. 11:464. doi: 10.3389/fnmol.2018.00464

PubMed Abstract | CrossRef Full Text | Google Scholar

Pomares, F. B., Roy, S., Funck, T., Feier, N. A., Thiel, A., Fitzcharles, M. A., et al. (2020). Upregulation of cortical GABAA receptor concentration in fibromyalgia. Pain 161, 74–82. doi: 10.1097/j.pain.0000000000001707

PubMed Abstract | CrossRef Full Text | Google Scholar

Prescott, S. A. (2015). Synaptic inhibition and disinhibition in the spinal dorsal horn. Prog. Mol. Biol. Transl. Sci. 131, 359–383. doi: 10.1016/bs.pmbts.2014.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, M. S., Winsvold, B. S., Chavez Chavez, S. O., Borte, S., Tsepilov, Y. A., Sharapov, S. Z., et al. (2021). Genome-wide association study identifies RNF123 locus as associated with chronic widespread musculoskeletal pain. Ann. Rheum. Dis. 80, 1227–1235. doi: 10.1136/annrheumdis-2020-219624

PubMed Abstract | CrossRef Full Text | Google Scholar

Rong, J., Li, S. H., and Li, X. J. (2007). Regulation of intracellular HAP1 trafficking. J. Neurosci. Res. 85, 3025–3029. doi: 10.1002/jnr.21326

PubMed Abstract | CrossRef Full Text | Google Scholar

Rong, J., McGuire, J. R., Fang, Z. H., Sheng, G., Shin, J. Y., Li, S. H., et al. (2006). Regulation of intracellular trafficking of huntingtin-associated protein-1 is critical for TrkA protein levels and neurite outgrowth. J. Neurosci. 26, 6019–6030. doi: 10.1523/JNEUROSCI.1251-06.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, I. J., Vaeroy, H., Javors, M., and Nyberg, F. (1992). Cerebrospinal fluid biogenic amine metabolites in fibromyalgia/fibrositis syndrome and rheumatoid arthritis. Arthritis Rheum. 35, 550–556. doi: 10.1002/art.1780350509

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandrini, G., Serrao, M., Rossi, P., Romaniello, A., Cruccu, G., and Willer, J. C. (2005). The lower limb flexion reflex in humans. Prog. Neurobiol. 77, 353–395. doi: 10.1016/j.pneurobio.2005.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasi, M., Vignoli, B., Canossa, M., and Blum, R. (2017). Neurobiology of local and intercellular BDNF signaling. Pflugers. Arch. 469, 593–610. doi: 10.1007/s00424-017-1964-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, G., Chang, G. Q., Lin, J. Y., Yu, Z. X., Fang, Z. H., Rong, J., et al. (2006). Hypothalamic huntingtin-associated protein 1 as a mediator of feeding behavior. Nat. Med. 12, 526–533. doi: 10.1038/nm1382

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, G., Xu, X., Lin, Y. F., Wang, C. E., Rong, J., Cheng, D., et al. (2008). Huntingtin-associated protein 1 interacts with Ahi1 to regulate cerebellar and brainstem development in mice. J. Clin. Invest. 118, 2785–2795. doi: 10.1172/JCI35339

PubMed Abstract | CrossRef Full Text | Google Scholar

Skaer, T. L. (2014). Fibromyalgia: disease synopsis, medication cost effectiveness and economic burden. Pharmacoeconomics 32, 457–466. doi: 10.1007/s40273-014-0137-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Skljarevski, V., and Ramadan, N. M. (2002). The nociceptive flexion reflex in humans – review article. Pain 96, 3–8. doi: 10.1016/s0304-3959(02)00018-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sluka, K. A., and Clauw, D. J. (2016). Neurobiology of fibromyalgia and chronic widespread pain. Neuroscience 338, 114–129. doi: 10.1016/j.neuroscience.2016.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanwar, S., Mattoo, B., Kumar, U., and Bhatia, R. (2019). Can aberrant spinal nociception be a marker of chronicity of pain in fibromyalgia syndrome? J. Clin. Neurosci. 65, 17–22. doi: 10.1016/j.jocn.2019.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Thabane, L., Mbuagbaw, L., Zhang, S., Samaan, Z., Marcucci, M., Ye, C., et al. (2013). A tutorial on sensitivity analyses in clinical trials: the what, why, when and how. BMC Med. Res. Methodol. 13:92. doi: 10.1186/1471-2288-13-92

PubMed Abstract | CrossRef Full Text | Google Scholar

Turner, S. D. (2014). QQman: an R package for visualizing GWAS results using Q-Q and manhattan plots. BioRxiv [preprint] doi: 10.1101/005165

CrossRef Full Text | Google Scholar

Twelvetrees, A. E., Lesept, F., Holzbaur, E. L. F., and Kittler, J. T. (2019). The adaptor proteins HAP1a and GRIP1 collaborate to activate the kinesin-1 isoform KIF5C. J. Cell Sci. 132:jcs215822. doi: 10.1242/jcs.215822

PubMed Abstract | CrossRef Full Text | Google Scholar

Twelvetrees, A. E., Yuen, E. Y., Arancibia-Carcamo, I. L., MacAskill, A. F., Rostaing, P., Lumb, M. J., et al. (2010). Delivery of GABAARs to synapses is mediated by HAP1-KIF5 and disrupted by mutant huntingtin. Neuron 65, 53–65. doi: 10.1016/j.neuron.2009.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Visscher, P. M., Wray, N. R., Zhang, Q., Sklar, P., McCarthy, M. I., Brown, M. A., et al. (2017). 10 Years of GWAS discovery: biology. function, and translation. Am. J. Hum. Genet. 101, 5–22. doi: 10.1016/j.ajhg.2017.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 38:e164. doi: 10.1093/nar/gkq603

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, K., Taskesen, E., van Bochoven, A., and Posthuma, D. (2017). Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8:1826. doi: 10.1038/s41467-017-01261-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Willer, J. C. (1977). Comparative study of perceived pain and nociceptive flexion reflex in man. Pain 3, 69–80. doi: 10.1016/0304-3959(77)90036-7

CrossRef Full Text | Google Scholar

Willer, J. C. (1990). Clinical exploration of nociception with the use of reflexologic techniques. Neurophysiol. Clin. 20, 335–356. doi: 10.1016/s0987-7053(05)80203-5

CrossRef Full Text | Google Scholar

Wolfe, F., Smythe, H. A., Yunus, M. B., Bennett, R. M., Bombardier, C., Goldenberg, D. L., et al. (1990). The American College of Rheumatology 1990 Criteria for the Classification of Fibromyalgia. Report of the Multicenter Criteria Committee. Arthritis Rheum 33, 160–172.

Google Scholar

Wroblewski, G., Islam, M. N., Yanai, A., Jahan, M. R., Masumoto, K. H., and Shinoda, K. (2018). Distribution of HAP1-immunoreactive cells in the retrosplenial-retrohippocampal area of adult rat brain and its application to a refined neuroanatomical understanding of the region. Neuroscience 394, 109–126. doi: 10.1016/j.neuroscience.2018.10.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L. L., Fan, Y., Li, S., Li, X. J., and Zhou, X. F. (2010). Huntingtin-associated protein-1 interacts with pro-brain-derived neurotrophic factor and mediates its transport and release. J. Biol. Chem. 285, 5614–5623. doi: 10.1074/jbc.M109.073197

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, J., Yang, H., Zhao, T., Sun, M., Xu, X., Zhou, X. F., et al. (2014). Huntingtin-associated protein 1 regulates postnatal neurogenesis and neurotrophin receptor sorting. J. Clin. Invest. 124, 85–98. doi: 10.1172/JCI69206

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, M., Lim, Y., Li, X., Zhong, J. H., and Zhou, X. F. (2011). Precursor of brain-derived neurotrophic factor (proBDNF) forms a complex with Huntingtin-associated protein-1 (HAP1) and sortilin that modulates proBDNF trafficking, degradation, and processing. J. Biol. Chem. 286, 16272–16284. doi: 10.1074/jbc.M110.195347

PubMed Abstract | CrossRef Full Text | Google Scholar