You are here

Article Text

Article menu
Basic and translational research
Extended report
Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region
  1. Marjolein J Peters1,2,
  2. Linda Broer3,
  3. Hanneke L D M Willemen4,
  4. Gudny Eiriksdottir5,
  5. Lynne J Hocking6,
  6. Kate L Holliday7,
  7. Michael A Horan8,
  8. Ingrid Meulenbelt9,
  9. Tuhina Neogi10,
  10. Maria Popham11,
  11. Carsten O Schmidt12,
  12. Anushka Soni13,
  13. Ana M Valdes11,
  14. Najaf Amin3,
  15. Elaine M Dennison14,15,
  16. Niels Eijkelkamp16,
  17. Tamara B Harris17,
  18. Deborah J Hart11,
  19. Albert Hofman3,
  20. Frank J P M Huygen18,
  21. Karen A Jameson14,
  22. Gareth T Jones19,
  23. Lenore J Launer17,
  24. Hanneke J M Kerkhof1,2,
  25. Marjolein de Kruijf1,2,18,
  26. John McBeth7,
  27. Margreet Kloppenburg20,21,
  28. William E Ollier22,
  29. Ben Oostra23,
  30. Antony Payton22,
  31. Fernando Rivadeneira1–,3,
  32. Blair H Smith24,
  33. Albert V Smith5,25,
  34. Lisette Stolk1,2,
  35. Alexander Teumer26,
  36. Wendy Thomson7,
  37. André G Uitterlinden1–,3,
  38. Ke Wang10,
  39. Sophie H van Wingerden3,
  40. Nigel K Arden14,27,
  41. Cyrus Cooper14,27,
  42. David Felson10,
  43. Vilmundur Gudnason5,25,
  44. Gary J Macfarlane19,
  45. Neil Pendleton8,
  46. P Eline Slagboom2,9,
  47. Tim D Spector11,
  48. Henry Völzke12,
  49. Annemieke Kavelaars4,
  50. Cornelia M van Duijn3,
  51. Frances M K Williams11,
  52. Joyce B J van Meurs1,2
  1. 1Department of Internal Medicine, Erasmus Medical Center Rotterdam, Rotterdam, The Netherlands
  2. 2The Netherlands Genomics Initiative-sponsored Netherlands Consortium for Healthy Aging (NGI-NCHA), Leiden/Rotterdam, The Netherlands
  3. 3Department of Epidemiology, Erasmus Medical Center Rotterdam, Rotterdam, The Netherlands
  4. 4Laboratory of Neuroimmunology and Developmental Origins of Disease, University Medical Center Utrecht, The Netherlands
  5. 5Icelandic Heart Association Research Institute, Kopavogur, Iceland
  6. 6Aberdeen Pain Research Collaboration (Musculoskeletal Research), University of Aberdeen, Aberdeen, UK
  7. 7Arthritis Research UK Epidemiology Unit, Manchester Academic Health Science Centre, University of Manchester, Manchester, UK
  8. 8Mental Health and Neurodegeneration Group, School Community Based Medicine, University of Manchester, Manchester, UK
  9. 9Department of Medical Statistics and Bioinformatics, Section of Molecular Epidemiology, Leiden University Medical Centre, Leiden, The Netherlands
  10. 10Clinical Epidemiology Unit, Boston University School of Medicine, Boston, Massachusetts, USA
  11. 11Department of Twin Research and Genetic Epidemiology, King's College London, London, UK
  12. 12Institute for Community Medicine, University of Greifswald, Greifswald, Germany
  13. 13NIHR Musculoskeletal Biomedical Research Unit, University of Oxford, Oxford, UK
  14. 14MRC Lifecourse Epidemiology Unit, University of Southampton, Southampton General Hospital, Southampton, UK
  15. 15School of Biological Sciences, Victoria University of Wellington, Wellington, New Zealand
  16. 16Molecular Nociception Group, University College London, London, UK
  17. 17Intramural Research Program, Laboratory of Epidemiology, Demography, and Biometry, National Institute on Aging, Bethesda, Maryland, USA
  18. 18Department of Anaesthesiology, Erasmus Medical Center Rotterdam, Rotterdam, The Netherlands
  19. 19Aberdeen Pain Research Collaboration (Epidemiology Group), University of Aberdeen, Aberdeen, UK
  20. 20Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands
  21. 21Department of Clinical Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
  22. 22Centre for Integrated Genomic Medical Research, University of Manchester, Manchester, UK
  23. 23Department of Clinical Genetics, Erasmus Medical Center Rotterdam, Rotterdam, The Netherlands
  24. 24Medical Research Institute, University of Dundee, Dundee, UK
  25. 25Department of Medicine, University of Iceland, Reykjavik, Iceland
  26. 26Institute of Functional Genomics, Ernst Moritz Arndt University Greifswald, University of Greifswald, Greifswald, Germany
  27. 27NIHR Biomedical Research Unit, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Oxford, UK
  1. Correspondence to Dr Joyce B J van Meurs, Genetic Laboratory Department of Internal Medicine, Erasmus Medical Center Rotterdam, Room Ee579b, Erasmus MC, PO Box 2040, Rotterdam 3000 CA, The Netherlands; j.vanmeurs{at}erasmusmc.nl

Abstract

Background and objectives Chronic widespread pain (CWP) is a common disorder affecting ∼10% of the general population and has an estimated heritability of 48–52%. In the first large-scale genome-wide association study (GWAS) meta-analysis, we aimed to identify common genetic variants associated with CWP.

Methods We conducted a GWAS meta-analysis in 1308 female CWP cases and 5791 controls of European descent, and replicated the effects of the genetic variants with suggestive evidence for association in 1480 CWP cases and 7989 controls. Subsequently, we studied gene expression levels of the nearest genes in two chronic inflammatory pain mouse models, and examined 92 genetic variants previously described associated with pain.

Results The minor C-allele of rs13361160 on chromosome 5p15.2, located upstream of chaperonin-containing-TCP1-complex-5 gene (CCT5) and downstream of FAM173B, was found to be associated with a 30% higher risk of CWP (minor allele frequency=43%; OR=1.30, 95% CI 1.19 to 1.42, p=1.2×10−8). Combined with the replication, we observed a slightly attenuated OR of 1.17 (95% CI 1.10 to 1.24, p=4.7×10−7) with moderate heterogeneity (I2=28.4%). However, in a sensitivity analysis that only allowed studies with joint-specific pain, the combined association was genome-wide significant (OR=1.23, 95% CI 1.14 to 1.32, p=3.4×10−8, I2=0%). Expression levels of Cct5 and Fam173b in mice with inflammatory pain were higher in the lumbar spinal cord, not in the lumbar dorsal root ganglions, compared to mice without pain. None of the 92 genetic variants previously described were significantly associated with pain (p>7.7×10−4).

Conclusions We identified a common genetic variant on chromosome 5p15.2 associated with joint-specific CWP in humans. This work suggests that CCT5 and FAM173B are promising targets in the regulation of pain.

  • Gene Polymorphism
  • Fibromyalgis/Pain Syndromes
  • Epidemiology

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/

https://doi.org/10.1136/annrheumdis-2012-201742

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Introduction

    Chronic widespread pain (CWP) is a common disorder, affecting about 10% of the general population.1 The prevalence of CWP increases with age for both men and women, but is more common in women at any age.1 CWP represents a major underestimated health problem and is associated with substantial impairment and a reduced quality of life. It has been related to a number of physical and affective symptoms such as fatigue, psychological distress and somatic symptoms.1 ,2 Chronic musculoskeletal pain is one of the most common conditions seen in rheumatology clinics and accounts for 6.2% of the total healthcare costs in The Netherlands every year.3 Further research is needed to be able to understand the causal mechanisms and optimal treatment for CWP patients.

    CWP causally relates to an initial local pain stimulus, such as an acute injury or athletic injuries, or another pain state such as low back pain or local pain due to osteoarthritis (OA) or rheumatic arthritis (RA).4–6 However, most injured subjects do not develop CWP, and only a proportion of patients with OA or RA develop CWP. We therefore hypothesise that several discrete stimuli may initiate CWP via a common final pathway that involves the generation of a central pain state through the sensitisation of second order spinal neurons.

    CWP is a complex trait since both environmental and genetic factors play a role in the aetiology. Heritability estimates of twin studies suggest that 48–52% of the variance in CWP occurrence is due to genetic factors, implying a strong genetic component.7 A number of studies have examined genetic variants for CWP. These candidate gene studies examined polymorphisms in genes involved in both the peripheral and the central nervous system.8 In particular, genes involved in neurotransmission (pathway of dopamine and serotonin9–19), and genes important for the hypothalamic–pituitary–adrenal axis have been considered.20 A number of genetic variants in these candidate genes were found to be associated with CWP, individual pain sites or experimental pain. However, no consistent significant associations have been demonstrated.

    The most studied gene in relation to pain is catechol-O-methyltransferase (COMT), an enzyme that degrades neurotransmitters including dopamine. The variant allele of rs4680 (or V158M) results in reduced enzymatic activity due to its effect on thermostability,21 and has been associated with reduced opioid activity in response to painful stimuli resulting in increased pain sensitivity.22 But also for COMT, no consistent results have been observed in genetic association studies.13 ,23–29

    Overall, the results have been conflicting, which is likely due to the modest sample sizes used and paucity of replication. In general, candidate studies are biased by previous knowledge of the aetiology of the disease under study. Since knowledge about the pathophysiology of CWP is poor, the chances of success using this approach are low. Therefore our objective was to identify genetic variants involved in CWP by means of a large-scale hypothesis-free genome-wide association study (GWAS) meta-analysis including 2788 cases and 13 780 controls. To our knowledge, this is the first study presenting a large-scale GWAS meta-analysis of chronic pain. The prevalence of CWP is approximately two times higher in women than in men and there is strong evidence that women tolerate less thermal and pressure pain than men.30 Therefore only women were included in this study to reduce heterogeneity and thereby increase power.

    Materials and methods

    We performed a meta-analysis (stage 1) of GWAS data of 1308 female Caucasian CWP cases and 5791 female Caucasian controls, derived from five studies, and focused our follow-up efforts on the single-nucleotide polymorphisms (SNPs) with suggestive evidence of association (p<1×10−5) with CWP (stage 2). The study outline is summarised in figure 1.

    Figure 1
    Figure 1

    Study outline. CWP, chronic widespread pain; GWAS, genome-wide association study.

    Phenotype

    CWP was defined as subjects having pain in the left side of the body, in the right side of the body, above the waist, below the waist, and in the axial skeleton (following the Fibromyalgia Criteria of the American College of Rheumatology2). Controls were defined as subjects not having CWP. Subjects using analgesics (ATC code: N0231) were excluded from the control group. Detailed descriptions of the study specific inclusion criteria are presented in supplementary table S1.

    Study design summary

    We combined the summary statistics of GWAS in a meta-analysis comprising 1308 CWP female Caucasian cases and 5791 female Caucasian controls (stage 1). We focused our follow-up efforts on the SNPs with suggestive evidence of association (p<1×10−5) with CWP in 1480 CWP cases and 7989 controls available for replication (stage 2).

    Subjects

    A full detailed description of all study cohorts is presented in table 1 and in the supplementary methods section. For the stage 1 analysis, we included studies from The Netherlands (the Erasmus Rucphen Family study (ERF study),32 Rotterdam Study I, II and III (RS-I, RS-II and RS-III)33), and the UK (TwinsUK34 ,35). All studies were approved by their institutional ethics review committees and all participants provided written informed consent. For our stage 2 analysis, we sought follow-up samples with pre-existing GWAS in silico data (stage 2a) as well as de novo genotyping (stage 2b). The studies are from the UK (the British 1958 Birth Cohort (1958BC),23 ,36–38 the Chingford Study (CHINGFORD),39 ,40 the Dyne Steel DNA Bank for Aging and Cognition (DSDBAC),41 the EPIdemiological study of FUNctional Disorders (EPIFUND),20 and the Hertfordshire Cohort Study (HCS)42); from Iceland (the Age, Gene/Environment Susceptibility Study (AGES)43); from the USA (the Framingham Osteoarthritis Study (FOA)44); from The Netherlands (the Genetics osteoARthritis and Progression Study (GARP)45); and from Germany (the Study of Health In Pomerania (SHIP)46 ,47). All studies were approved by the local ethics committees and all participants provided written informed consent.

    View this table:
    Table 1

    Overview of all participating studies

    Genotyping, quality control and imputation

    Genotyping of the stage 1 cohorts was done by Illumina Infinium HumanHap550 Beadchip (RS-I and RS-II), the Illumina Infinium HumanHap610 (RS-II, RS-III, and TwinsUK), or the Illumina Infinium HumanHap300 (ERF and TwinsUK). More details about the genotyping, quality control (QC), and imputation are shown in the supplementary methods section. Complete information on genotyping protocols and QC measures for all stage 1 cohorts is described in the supplementary material (see supplementary table S2). Detailed descriptions of the QC and imputation procedures are provided in the supplementary material (see supplementary table S3).

    Genotypes of the stage 2a studies (1958BC, AGES, DSDBAC, FOA, GARP and SHIP) were obtained from SNP arrays and imputed data. Where unavailable, proxy SNPs were selected based on high linkage disequilibrium (LD). The stage 2b studies (CHINGFORD, EPIFUND and HCS) performed de novo genotyping, using both Sequenom iPLEX and TaqMan-based assays (supplementary methods). Genotyping platforms, calling algorithms, quality control before imputation, imputation methods and analysis software used were all study-specific (see supplementary tables S4 and S5). The explicit number of follow-up SNPs genotyped in the different studies and whether the original or a proxy SNP was used is summarised in supplementary table S6.

    GWAS analysis in the stage 1 studies

    CWP was analysed as a binary trait (cases vs controls) using logistic regression under an additive model with adjustment for age and body mass index (see supplementary table S7). To adjust for population substructure, we included the four most important PCs as covariates in the regression analysis of RS-I, RS-II and RS-III. These PCs were derived from a multidimensional scaling analysis of identity-by-state distances, using PLINK software.48 Detailed descriptions of the GWAS methods are provided in supplementary table S8).

    Stage 1: GWAS meta-analysis

    p Values for association were combined using the Meta-Analysis Tool for genome-wide association scans (METAL).49 The genomic control method50 as implemented in METAL was used to correct for any residual population stratification or relatedness not accounted for by the four most important PCs. A p value <5×10−8 was considered genome-wide significant while a p value <1×10−5 was considered suggestive.51 Power calculations were performed using CaTS software (http://www.sph.umich.edu/csg/abecasis/CaTS/). Using Bonferroni correction (p<5×10−8), power calculations showed that we had approximately 80% power to detect an OR of 1.30 for SNPs with a minor allele frequency (MAF) of 0.43, given a disease prevalence of 10% for 1308 cases and 5791 controls in the discovery group. Using a p value <1×10−5, we had 80% power to detect an OR of 1.25.

    SNP selection for replication

    We aimed to select SNPs for replication (stage 2) that were enriched for signals of association with CWP. All SNPs with suggestive evidence for association in the stage 1 analyses were selected and separated into independent loci by taking the most significantly associated SNP and eliminating all SNPs that have a HapMap CEU pairwise correlation coefficient r2>0.8 with that SNP using the PLINK software.

    Meta-analysis of stage 1 and stage 2 results

    We combined the stage 1 and stage 2 association results to derive a combined meta-analysis for the suggestively associated loci. METAL was used to conduct a fixed-effects meta-analysis as in stage 1. Estimated heterogeneity variance and forest plots were generated using comprehensive meta-analysis (http://www.meta-analysis.com).

    Functional analysis of associated SNPs

    To determine whether the associated SNPs have any regulatory effect on gene expression levels, we checked their effect (and the effect of the linked SNPs) on the expression levels of their neighbouring genes. We used the 1000 genomes data in the SNAP software52 ,53 to identify those SNPs having LD thresholds of r2>0.1. We searched two publicly available eQTL databases: the NCBI GTEx (Genotype-Tissue Expression) eQTL browser (http://www.ncbi.nlm.nih.gov/gtex/GTEX2/gtex.cgi) and the expression Quantitative Trait Loci database (http://eqtl.uchicago.edu/cgi-bin/gbrowse/eqtl/). We used SIFT54 to predict whether the coding non-synonymous variant causing an amino acid substitution affects protein function.

    RNA expression analyses in mice

    For functional follow-up, two independent mouse models of inflammatory pain were studied. The first model was based on carrageenan injections; female C57Bl/6 mice received an intraplantar injection of 20 μl λ-carrageenan (2% (w/v), Sigma Aldrich, Zwijndrecht, the Netherlands) in saline in both hind paws.55 The second model was based on Complete Freund's Adjuvant (CFA) injections; male C57Bl/6 mice (Harlan Laboratories) received an intraplantar injection of 20 μl CFA (Sigma-Aldrich) in saline in both hind paws.56 Controls were injected with saline only. At day 3 (after CFA injection) or day 6 (after carrageenan injection), thermal sensitivity (heat withdrawal latency time) was measured using the Hargreaves (IITC Life Science, Woodland Hills, California, USA) test as described.57 Intensity of the light beam was chosen to induce heat withdrawal latency time of approximately 8 s at baseline.

    After measurement the mice were sacrificed and the lumbar (L2–L5) spinal cord and the dorsal root ganglions (DRG) (L2–L5) were isolated. These areas of spinal cord and DRG were selected because pain transmission from the hind paws is mediated via primary sensory neurons that have their cell bodies in the lumbar DRG, and transmit the signal to the lumbar spinal cord through sensory fibres in the dorsal roots. Total RNA was isolated and mRNA levels of Cct5 and Fam173b were measured in the spinal cord and the DRG. For more details, see the supplementary methods section.

    All experiments were performed in accordance with international guidelines and approved by the experimental animal committee of the University Medical Center Utrecht (carrageenan experiment) or the UK Home Office Animals (Scientific Procedures) Act 1986 (CFA experiment). Mice used for the carrageenan experiment were bred and maintained in the animal facility of the University of Utrecht (The Netherlands).

    Systemic review of genetic variants previously described

    We systematically searched for associations earlier reported with pain in the HugeNavigator PhenoPedia database.58 We used the search term ‘pain’ and checked all publications for genes and SNPs associated with pain at least twice. Genes and SNPs associated with drug therapy, facial pain, migraine and postoperative pain were excluded. For all reported SNPs, we examined their association with CWP in our stage 1 meta-analysis. The significance threshold was set at p<8×10−4 using Bonferroni correction for 65 independent genetic loci. Again, power calculations were performed using CaTS software (http://www.sph.umich.edu/csg/abecasis/CaTS/). With an α level of 8×10−4, power calculations showed that we had approximately 80% power to detect an OR of 1.22 for SNPs with a minor allele frequency of 20% or higher.

    Results

    GWAS meta-analysis for CWP

    The Manhattan plot and quantile–quantile plot of the initial stage 1 meta-analysis are presented in figure 2. In total, 2 224 068 SNPs (directly genotyped or imputed) were tested for association. The overall genomic control lambda (λGC) was 1.007, indicating no significant population stratification. We identified two SNPs which were genome-wide significant (p<5×10−8), and another 39 SNPs with suggestive evidence for association (p<1×10−5) located in 10 independent genomic regions. The most significant association was observed for two imputed highly correlated SNPs (r2=0.97) located upstream of the chaperonin-containing-TCP1-complex-5 gene (CCT5) and downstream of the FAMily with sequence similarity 173, member B gene (FAM173B) (rs13361160, p=1.2×10−8 and rs2386592, p=2.6×10−8). For both SNPs, the minor allele (MAF=43%) was associated with a 30% higher risk for CWP (OR=1.30, 95% CI 1.19 to 1.42).

    Figure 2
    Figure 2

    Genome-wide association results for chronic widespread pain (CWP) (stage 1). (A) Manhattan plot showing the p value of association tests for about 2 million SNPs with CWP in the stage 1 meta-analysis. SNPs are plotted on the x-axis according to their position on each chromosome. On the y-axis, the association p values with CWP are shown (as −log 10 p values). The grey solid horizontal line represents the p value threshold of 5×10−8 (genome-wide significance). The grey dashed horizontal line represents the p value threshold of 1×10−5 (the level for suggestive evidence): SNPs in loci reaching 1×10−5 were tested for replication. (B) Quantile–quantile (QQ) plot of SNPs. The blue area represents the 95% CI around the test statistics. A QQ plot compares the additive model statistics to those expected under the null distribution using fixed effects for all analysed HapMAP CEU imputed SNPs passing quality control criteria. This figure is only reproduced in colour in the online version.

    Meta-analysis of GWAS replication

    For the 10 independent SNPs with suggestive evidence, we pursued in silico replication data in six studies (stage 2a: 1203 CWP cases and 5032 controls) and performed de novo genotyping in subjects from three additional studies (stage 2b: 277 CWP cases and 2957 controls) (a detailed description of the studies is presented in table 1 and supplementary methods). The summary results of the stage 1 and 2 meta-analysis are presented in table 2. After combining the results of stage 1 and stage 2, the top SNP was rs13361160 (OR=1.17, 95% CI 1.10 to 1.24, p=4.7×10−7, I2=28.4%). Figure 3 shows a forest plot of the association of rs13361160 with CWP across the stage 1 and stage 2 studies. The overall effect in the replication studies (stage 2 studies) was in a consistent direction but not significant (OR=1.06, 95% CI 0.98 to 1.16, p=0.16). In the combined analysis, moderate heterogeneity was observed (I2=28.4%). supplementary table S1 shows the different pain assessment methods used in the different studies to define CWP. Since four out of five stage 1 studies included joint-specific pain only (ERF, RS-I, RS-II and RS-III), we performed a sensitivity analysis in which stage 2 cohorts using non-joint pain were excluded (1958BC, DSDBAC, EPIFUND, HCS and SHIP). This resulted in a combined OR of 1.23 (95% CI 1.14 to 1.32, p=3.4×10−8, I2 = 0%). An overview of the results of the combined meta-analysis and the separate stage 1 and stage 2 analyses is presented in table 3.

    View this table:
    Table 2

    Association results of the 10 top hits

    View this table:
    Table 3

    Top hit association results

    Figure 3
    Figure 3

    Forest plot of the association of rs13361160 SNP with chronic widespread pain (CWP). Study specific estimates and summary association between rs13361160 and CWP are shown. This figure is only reproduced in colour in the online version.

    Functional analysis of rs13361160 and rs2386592

    The SNPs rs13361160 and rs2386592 (r2=0.97) are annotated to the 5p15.2-region and located 81 kb upstream of CCT5 and 57 kb downstream of FAM173B (figure 4). We tested whether rs13361160 and rs2386592 and their linked SNPs (r2>0.1) affected gene expression levels of CCT5 or FAM173B. In total, we identified 130 SNPs in LD with our top SNPs, of which two SNPs were located in the coding region: one synonymous SNP rs1042392 in the CCT5 gene (r2=0.16, D′=0.85) and one non-synonymous SNP rs2438652 in the FAM173B gene (r2 = 0.17, D′=1.0) (see supplementary table S9). The minor allele of rs2438652 causes a threonine-to-methionine substitution (T75M) which is thought to be functionally neutral. SNPs rs13361160 and rs2386592 were not recorded as influencing the expression levels of CCT5 and FAM173B, however the linked intronic SNP rs2445871 (r2=0.14 for both) had a direct eQTL effect on FAM173B expression levels in liver tissue.59

    Figure 4
    Figure 4

    Regional plot of locus 5p15.2. On the x-axis, SNPs are plotted according to their position in a 400-kb window around rs13361160. On the y-axis, the association p values with chronic widespread pain are shown (as −log 10 p values). The purple diamond highlights the most significant SNP rs13361160. Blue peaks indicate recombination sites, and the SNPs surrounding the most significant SNP are colour coded to identify their strength of linkage disequilibrium with the most significant SNP (pairwise r2 values of the HapMap CEU samples). Genes and the direction of transcription are shown at the bottom of the plot. This figure is only reproduced in colour in the online version.

    RNA expression analysis in mice

    We studied gene expression levels of the two nearest genes, Cct5 and Fam173b, in the lumbar spinal cord and the DRG in two independent mouse models of chronic inflammatory pain. In both the carrageenan treated group and the CFA treated group, mice had shorter heat withdrawal latency times than mice injected with saline only, confirming enhanced pain sensitivity (p<0.001) (see supplementary figure S1).

    The results from the multivariate analysis using the two genes (Cct5 and Fam173b examined as dependent variables), the different treatments (saline, carrageenan and CFA) and the different tissues (DRG and spinal cord) confirmed that there is a significant treatment effect for Cct5 (F(2,25)=3.399, p=0.0049), as well as for Fam173b (F(2,25)=4.911, p=0.016). Moreover, both genes showed a significant tissue effect (Cct5: F(1,25)=13.595, p=0.001, and Fam173b: F(1,25)=13.522, p=0.001), as well as a significant interaction between tissue and treatment (Cct5: F(2,25)=6.424, p=0.006, and Fam173b: F(1,25)=4.196, p=0.027) (figure 5). These findings indicate that in spinal cord but not in DRG, both Fam173b and Cct5 expression levels were up-regulated in response to two different inducers of inflammatory pain. DRG Fam173b and Cct5 expression levels in CFA/carrageenan-treated mice were indistinguishable from saline-treated mice.

    Figure 5
    Figure 5

    Quantitative PCR analysis of gene expression levels in the lumbar (L2–L5) spinal cord (A) and the dorsal root ganglions (DRG) (B) of mice after intraplantar saline (n=3), carrageenan (n=4), and Complete Freund's Adjuvant (CFA) (n=4) injection. Spinal cord and DRG were collected and analysed for RNA levels of Cct5 and Fam173b. Data were normalised for Gapdh and β-actin (housekeeping genes) expression. Data are expressed as mean±SEM, *=p<0.05.

    Candidate SNPs previously associated with chronic pain

    We examined whether genetic variants previously described for association with pain were associated with CWP in our large stage 1 meta-analysis. We identified a total of 44 genes, of which 136 SNPs had been reported at least twice with any pain phenotype (excluding facial pain, migraine, postoperative pain and response to drug therapy), and we examined the association of these 136 SNPs with CWP in the GWAS stage 1 meta-analysis. Out of 136 candidate SNPs, we were able to check 92 common SNPs (MAF>5%) in 65 independent genetic loci (see supplementary table S10). Five SNPs had a too low MAF (<=5%) and 39 SNPs were not genotyped or imputed in our meta-analysis. None of the earlier reported SNPs passed the significance threshold (p<8×10−4). Interestingly, the strongest associated SNPs are located in three genes that have been reported to be associated with pain phenotypes most frequently: COMT, GCH1 (GTP cyclo-hydrolase 1) and OPRM1 (mu opioid receptor). The effects of the SNPs in GCH1 are in the same direction as reported earlier60–62: individuals having the minor allele for rs10483639, rs4411417 or rs752688 have 15% less pain than those exhibiting the common alleles. The effect of the SNP rs599548 in OPRM1 is also in the same direction as reported earlier63: those having the minor allele for rs599548 have 19% more pain than those exhibiting the major allele. The two COMT SNPs are in weak LD with the well-known amino acid changing variant rs4860, but previously have not been reported to be significantly associated with pain.23 ,64 We have found a protective effect for the minor allele of rs2020917 (those having a minor allele have 15% less pain) and an adverse effect for the minor allele of rs5993883: those having the minor allele of rs5993883 have 14% more pain.

    Discussion

    In this study, we identified a genetic variant near CCT5 and FAM173B to be associated with CWP. Chronic pain coincided with higher RNA expression of Cct5 and Fam173b in the lumbar spinal cord of mouse models of inflammatory pain. This finding indicates that both genes in the 5p15.2 region are regulated in the context of inflammatory pain.

    Interestingly, Bouhouche et al65 reported a human pedigree in which a CCT5 mutation caused hereditary sensory neuropathy (Online Mendelian Inheritance in Man (OMIM) ID=610150), a syndrome characterised by a sensory deficit in the distal portion of the lower extremities, chronic perforating ulcerations of the feet and progressive destruction of underlying bones. Symptoms can include pain and numbness, tingling in the hands, legs or feet, and extreme sensitivity to touch. CCT5 is a subunit of the chaperonin containing t-complex polypeptide 1 (TCP-1) which assists in protein folding and assembly in the brain.66 CCT5 interacts with the serine/threonine-protein phosphatase 4 catalytic subunit PP4C.67–69 Zhang et al70 confirmed that protein phosphates like PPP4C may have a regulatory effect on the central sensitisation of nociceptive transmission in the spinal cord. Interestingly, sensitisation is thought to contribute to chronic inflammatory pain.71 Since the function of the FAM173B gene is not yet known, it is difficult to postulate the mechanism by which this gene could influence CWP. Further research into the genes in this locus is needed to ascertain whether either or both CCT5 and FAM173B are driving the observed association.

    By combining the effects across the different stage 2 studies, moderate heterogeneity was observed in the meta-analysis. This heterogeneity might be caused by different pain assessment methods used by the stage 2 cohorts. In particular, four cohorts asked the participants about joint pain specifically, while the other five also included non-joint pain. When the non-joint pain phenotype were excluded, the heterogeneity across the cohorts reduced to 0% and the overall p value for rs13361160 now reached genome-wide significance by combining the stage 1 and stage 2 effects. This might suggest that indeed phenotype heterogeneity was introduced by including non-joint pain. In general, it is anticipated that pain is a very complex trait, with different aetiological pathways introducing phenotypic heterogeneity.

    A limitation of our study is that we were not able to examine possible phenotype subgroups, such as individuals with RA, a chronic systemic inflammatory disorder that principally affects the synovial joints. Stratifying these groups of individuals might serve to increase power to find genetic loci. We here decided to analyse all CWP cases together, based on the hypothesis that several discrete stimuli need to initiate CWP via a common final pathway that involves the generation of a central pain state through the sensitisation of second order spinal neurons. In addition, the prevalence of RA is very low (about 0.5–1%),72 and the earlier defined GWAS hits for RA (ie, the HLA locus)73 were not in our top list. So, we assume the results were not dominated by this small number of individuals with RA.

    It would be helpful to dissect the phenotype of pain into quantitative sub-phenotypes, for example by measuring pain sensitivity and pain thresholds for temperature or pressure,74 or by examining functional MRIs.75 The use of quantitative and possibly more objective pain measurements in response to painful stimuli (rather than reported pain) will be of pivotal importance for future pain research. Because we have focused on the clinical pain definition using questionnaires and pain homunculus, we accept that we may have missed true pain susceptibility alleles. However, this study represents the largest genome-wide meta-analysis looking into the genetics of human CWP to date. The experiments in two independent mouse models of chronic inflammatory pain showed that the expression of Cct5 and Fam173b was higher in the lumbar spinal cord of mice with chronic inflammatory pain but not in DRG. In the spinal cord, the expression profiles of both genes were up-regulated in response to two different inducers of inflammatory pain. These findings indicate that both genes in the 5p15.2 region are co-regulated in the spinal cord during inflammation-induced pain in both independent pain models, thereby possibly contributing to the neurobiology of pain. In the lumbar DRG, containing the cell bodies of the primary sensory neurons that detect pain signals from the hind paws, Cct5 and Fam173b gene expression levels did not change by inflammation. Because of these complementary results from the two independent tissues (spinal cord and DRG), we hypothesise that the 5p15.2 region is likely to play a role in spinal central pain processing and not in regulating primary sensory neuron responses.

    In the study of candidate genes previously reported to be associated with a pain phenotype, we showed that none of the 92 studied variants were significantly associated with CWP in our GWAS meta-analysis. This can be explained by the fact that many of the previous reported loci were studied in relative modest sample sizes and in a large variety of pain phenotypes.76 Power calculations show that we had approximately 80% power to detect an OR as low as 1.22 for SNPs with an allele frequency of 20% or higher. So, even in this large meta-analysis, power was still modest to detect small ORs and we therefore cannot exclude smaller effect sizes of the tested variants, resulting in lack of reproducibility.77 This lack of reproducibility of SNPs in candidate genes in large GWAS meta-analyses has been shown before for other phenotypes such as bone mineral density (BMD).78 It is interesting to note that among the candidate SNPs, the strongest associated ones were located in the three most studied pain genes, COMT, GCH1 and OPRM1. The directions of the effects of these SNPs were the same as reported earlier, which would support true associations.

    In conclusion, our study reports a GWAS meta-analysis on CWP. We identified the genetic variant rs13361160 at the 5p15.2 locus, located 81 kb upstream of the CCT5 gene and 57 kb downstream of the FAM173B gene, to be associated with CWP. We showed an increase in expression levels of Cct5 and Fam173b in the spinal cord of inflammatory pain models of mice, and since these genes both seem to influence the central mechanism of sensitisation, they may represent a novel pathway involved in pain sensation.

    References

    1. ↵
      1. Croft P,
      2. Schollum J,
      3. Silman A
      . Population study of tender point counts and pain as evidence of fibromyalgia. BMJ 1994;309:696–9.
      OpenUrlAbstract/FREE Full Text
    2. ↵
      1. Wolfe F,
      2. Smythe HA,
      3. Yunus MB,
      4. et al
      . The American College of Rheumatology 1990 criteria for the classification of fibromyalgia. Report of the Multicenter Criteria Committee. Arthritis Rheum 1990;33:160–72.
      OpenUrlCrossRefPubMedWeb of Science
    3. ↵
      1. Meerding WJ,
      2. Bonneux L,
      3. Polder JJ,
      4. et al
      . Demographic and epidemiological determinants of healthcare costs in Netherlands: cost of illness study. BMJ 1998;317:111–15.
      OpenUrlAbstract/FREE Full Text
    4. ↵
      1. Buskila D,
      2. Mader R
      . Trauma and work-related pain syndromes: risk factors, clinical picture, insurance and law interventions. Best Pract Res Cl Rh 2011;25:199–207.
      OpenUrlCrossRef
    5. ↵
      1. Leffler AS,
      2. Kosek E,
      3. Lerndal T,
      4. et al
      . Somatosensory perception and function of diffuse noxious inhibitory controls (DNIC) in patients suffering from rheumatoid arthritis. Eur J Pain-London 2002;6:161–76.
      OpenUrlCrossRef
    6. ↵
      1. Imamura M,
      2. Imamura ST,
      3. Kaziyama HHS,
      4. et al
      . Impact of nervous system hyperalgesia on pain, disability, and quality of life in patients with knee osteoarthritis: a controlled analysis. Arthrit Rheum-Arthr 2008;59:1424–31.
      OpenUrlCrossRef
    7. ↵
      1. Kato K,
      2. Sullivan PF,
      3. Evengard B,
      4. et al
      . Importance of genetic influences on chronic widespread pain. Arthritis Rheum 2006;54:1682–6.
      OpenUrlCrossRefPubMedWeb of Science
    8. ↵
      1. Limer KL,
      2. Nicholl BI,
      3. Thomson W,
      4. et al
      . Exploring the genetic susceptibility of chronic widespread pain: the tender points in genetic association studies. Rheumatology (Oxford) 2008;47:572–7.
      OpenUrlAbstract/FREE Full Text
    9. ↵
      1. Nicholl BI,
      2. Holliday KL,
      3. Macfarlane GJ,
      4. et al
      . Association of HTR2A polymorphisms with chronic widespread pain and the extent of musculoskeletal pain: results from two population-based cohorts. Arthritis Rheum 2011;63:810–18.
      OpenUrlCrossRefPubMedWeb of Science
    10. ↵
      1. Frank B,
      2. Niesler B,
      3. Bondy B,
      4. et al
      . Mutational analysis of serotonin receptor genes: HTR3A and HTR3B in fibromyalgia patients. Clin Rheumatol 2004;23:338–44.
      OpenUrlPubMedWeb of Science
    11. ↵
      1. Bondy B,
      2. Spaeth M,
      3. Offenbaecher M,
      4. et al
      . The T102C polymorphism of the 5-HT2A-receptor gene in fibromyalgia. Neurobiol Dis 1999;6:433–9.
      OpenUrlCrossRefPubMedWeb of Science
    12. ↵
      1. Gursoy S,
      2. Erdal E,
      3. Herken H,
      4. et al
      . Association of T102C polymorphism of the 5-HT2A receptor gene with psychiatric status in fibromyalgia syndrome. Rheumatol Int 2001;21:58–61.
      OpenUrlCrossRefPubMedWeb of Science
    13. ↵
      1. Tander B,
      2. Gunes S,
      3. Boke O,
      4. et al
      . Polymorphisms of the serotonin-2A receptor and catechol-O-methyltransferase genes: a study on fibromyalgia susceptibility. Rheumatol Int 2008;28:685–91.
      OpenUrlCrossRefPubMedWeb of Science
    14. ↵
      1. Su SY,
      2. Chen JJ,
      3. Lai CC,
      4. et al
      . The association between fibromyalgia and polymorphism of monoamine oxidase A and interleukin-4. Clin Rheumatol 2007;26:12–16.
      OpenUrlCrossRefPubMedWeb of Science
    15. ↵
      1. Offenbaecher M,
      2. Bondy B,
      3. de Jonge S,
      4. et al
      . Possible association of fibromyalgia with a polymorphism in the serotonin transporter gene regulatory region. Arthritis Rheum 1999;42:2482–8.
      OpenUrlCrossRefPubMedWeb of Science
    16. ↵
      1. Cohen H,
      2. Buskila D,
      3. Neumann L,
      4. et al
      . Confirmation of an association between fibromyalgia and serotonin transporter promoter region (5- HTTLPR) polymorphism, and relationship to anxiety-related personality traits. Arthritis Rheum 2002;46:845–7.
      OpenUrlPubMedWeb of Science
    17. ↵
      1. Gursoy S
      . Absence of association of the serotonin transporter gene polymorphism with the mentally healthy subset of fibromyalgia patients. Clin Rheumatol 2002;21:194–7.
      OpenUrlCrossRefPubMedWeb of Science
    18. ↵
      1. Buskila D,
      2. Cohen H,
      3. Neumann L,
      4. et al
      . An association between fibromyalgia and the dopamine D4 receptor exon III repeat polymorphism and relationship to novelty seeking personality traits. Mol Psychiatry 2004;9:730–1.
      OpenUrlCrossRefPubMedWeb of Science
    19. ↵
      1. Potvin S,
      2. Larouche A,
      3. Normand E,
      4. et al
      . DRD3 Ser9Gly polymorphism is related to thermal pain perception and modulation in chronic widespread pain patients and healthy controls. J Pain 2009;10:969–75.
      OpenUrlCrossRefPubMedWeb of Science
    20. ↵
      1. Holliday KL,
      2. Nicholl BI,
      3. Macfarlane GJ,
      4. et al
      . Genetic variation in the hypothalamic-pituitary-adrenal stress axis influences susceptibility to musculoskeletal pain: results from the EPIFUND study. Ann Rheum Dis 2010;69:556–60.
      OpenUrl
    21. ↵
      1. Syvanen AC,
      2. Tilgmann C,
      3. Rinne J,
      4. et al
      . Genetic polymorphism of catechol-O-methyltransferase (COMT): correlation of genotype with individual variation of S-COMT activity and comparison of the allele frequencies in the normal population and parkinsonian patients in Finland. Pharmacogenetics 1997;7:65–71.
      OpenUrlCrossRefPubMedWeb of Science
    22. ↵
      1. Zubieta JK,
      2. Heitzeg MM,
      3. Smith YR,
      4. et al
      . COMT val158met genotype affects mu-opioid neurotransmitter responses to a pain stressor. Science 2003;299:1240–3.
      OpenUrlAbstract/FREE Full Text
    23. ↵
      1. Hocking LJ,
      2. Smith BH,
      3. Jones GT,
      4. et al
      . 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;149:143–51.
      OpenUrlCrossRefPubMedWeb of Science
    24. ↵
      1. Gursoy S,
      2. Erdal E,
      3. Herken H,
      4. et al
      . Significance of catechol-O-methyltransferase gene polymorphism in fibromyalgia syndrome. Rheumatol Int 2003;23:104–7.
      OpenUrlPubMedWeb of Science
    25. ↵
      1. Cohen H,
      2. Neumann L,
      3. Glazer Y,
      4. et al
      . The relationship between a common catechol-O-methyltransferase (COMT) polymorphism val(158) met and fibromyalgia. Clin Exp Rheumatol 2009;27(5 Suppl 56):S51–6.
      OpenUrlPubMed
    26. ↵
      1. Vargas-Alarcon G,
      2. Fragoso JM,
      3. Cruz-Robles D,
      4. et al
      . Catechol-O-methyltransferase gene haplotypes in Mexican and Spanish patients with fibromyalgia. Arthritis Res Ther 2007;9:R110.
      OpenUrlCrossRefPubMed
    27. ↵
      1. van Meurs JB,
      2. Uitterlinden AG,
      3. Stolk L,
      4. et al
      . A functional polymorphism in the catechol-O-methyltransferase gene is associated with osteoarthritis-related pain. Arthritis Rheum 2009;60:628–9.
      OpenUrlCrossRefPubMedWeb of Science
    28. ↵
      1. Hagen K,
      2. Pettersen E,
      3. Stovner LJ,
      4. et al
      . No association between chronic musculoskeletal complaints and Val158Met polymorphism in the Catechol-O-methyltransferase gene. The HUNT study. BMC Musculoskelet Disord 2006;7:40.
      OpenUrlCrossRefPubMed
    29. ↵
      1. Nicholl BI,
      2. Holliday KL,
      3. Macfarlane GJ,
      4. et al
      . No evidence for a role of the catechol-O-methyltransferase pain sensitivity haplotypes in chronic widespread pain. Ann Rheum Dis 2010;69:2009–12.
      OpenUrl
    30. ↵
      1. Racine M,
      2. Tousignant-Laflamme Y,
      3. Kloda LA,
      4. et al
      . A systematic literature review of 10years of research on sex/gender and experimental pain perception - Part 1: Are there really differences between women and men? Pain 2011;153:61935.
      OpenUrlWeb of Science
    31. ↵
      1. Skrbo A,
      2. Begovic B,
      3. Skrbo S
      . [Classification of drugs using the ATC system (Anatomic, Therapeutic, Chemical Classification) and the latest changes] Klasificiranje lijekova po ATC sistemu i najnovije izmjene. Med Arh 2004;58(1 Suppl 2):138–41.
      OpenUrlPubMed
    32. ↵
      1. Aulchenko YS,
      2. Heutink P,
      3. Mackay I,
      4. et al
      . Linkage disequilibrium in young genetically isolated Dutch population. Eur J Hum Genet 2004;12:527–34.
      OpenUrlCrossRefPubMedWeb of Science
    33. ↵
      1. Hofman A,
      2. van Duijn CM,
      3. Franco OH,
      4. et al
      . The Rotterdam Study: 2012 objectives and design update. Eur J Epidemiol 2011;26:657–86.
      OpenUrlCrossRefPubMed
    34. ↵
      1. Spector TD,
      2. MacGregor AJ
      . The St. Thomas’ UK Adult Twin Registry. Twin Res 2002;5:440–3.
      OpenUrlCrossRefPubMedWeb of Science
    35. ↵
      1. Spector TD,
      2. Williams FM
      . The UK Adult Twin Registry (TwinsUK). Twin Res Hum Genet 2006;9:899–906.
      OpenUrlCrossRefPubMedWeb of Science
    36. ↵
      1. Power C,
      2. Atherton K,
      3. Manor O
      . Co-occurrence of risk factors for cardiovascular disease by social class: 1958 British birth cohort. J Epidemiol Community Health 2008;62:1030–5.
      OpenUrlAbstract/FREE Full Text
    37. ↵
      Consortium UKPsD, Wellcome Trust Case Control C, Spencer CCet al. Dissection of the genetics of Parkinson's disease identifies an additional association 5′ of SNCA and multiple associated haplotypes at 17q21. Hum Mol Genet 2011;20:345–53.
      OpenUrlAbstract/FREE Full Text
    38. ↵
      Wellcome Trust Case Control C. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 2007;447:661–78.
      OpenUrlCrossRefPubMedWeb of Science
    39. ↵
      1. Hart DJ,
      2. Spector TD
      . Cigarette smoking and risk of osteoarthritis in women in the general population: the Chingford study. Ann Rheum Dis 1993;52:93–6.
      OpenUrlAbstract/FREE Full Text
    40. ↵
      1. Hart DJ,
      2. Spector TD
      . The relationship of obesity, fat distribution and osteoarthritis in women in the general population: the Chingford Study. J Rheumatol 1993;20:331–5.
      OpenUrlPubMedWeb of Science
    41. ↵
      1. Rabbitt P,
      2. McInnes L,
      3. Diggle P,
      4. et al
      . The University of Manchester longitudinal study of cognition in normal healthy old age, 1983 through 2003. Aging, Neuropsychology and Cognition 2004;11:P245–79.
      OpenUrlCrossRef
    42. ↵
      1. Syddall HE,
      2. Aihie Sayer A,
      3. Dennison EM,
      4. et al
      . Cohort profile: the Hertfordshire cohort study. Int J Epidemiol 2005;34:1234–42.
      OpenUrlFREE Full Text
    43. ↵
      1. Harris TB,
      2. Launer LJ,
      3. Eiriksdottir G,
      4. et al
      . Age, Gene/Environment Susceptibility-Reykjavik Study: multidisciplinary applied phenomics. Am J Epidemiol 2007;165:1076–87.
      OpenUrlAbstract/FREE Full Text
    44. ↵
      1. Felson DT,
      2. Zhang Y,
      3. Hannan MT,
      4. et al
      . The incidence and natural history of knee osteoarthritis in the elderly. The Framingham Osteoarthritis Study. Arthritis Rheum 1995;38:1500–5.
      OpenUrlCrossRefPubMedWeb of Science
    45. ↵
      1. Riyazi N,
      2. Meulenbelt I,
      3. Kroon HM,
      4. et al
      . Evidence for familial aggregation of hand, hip, and spine but not knee osteoarthritis in siblings with multiple joint involvement: the GARP study. Ann Rheum Dis 2005;64:438–43.
      OpenUrl
    46. ↵
      1. John U,
      2. Greiner B,
      3. Hensel E,
      4. et al
      . Study of Health In Pomerania (SHIP): a health examination survey in an east German region: objectives and design. Soz Praventivmed 2001;46:186–94.
      OpenUrlCrossRefPubMedWeb of Science
    47. ↵
      1. Volzke H,
      2. Alte D,
      3. Schmidt CO,
      4. et al
      . Cohort profile: the study of health in Pomerania. Int J Epidemiol 2011;40:294–307.
      OpenUrlFREE Full Text
    48. ↵
      1. Purcell S,
      2. Neale B,
      3. Todd-Brown K,
      4. et al
      . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–75.
      OpenUrlCrossRefPubMed
    49. ↵
      1. Willer CJ,
      2. Li Y,
      3. Abecasis GR
      . METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010;26:2190–1.
      OpenUrlAbstract/FREE Full Text
    50. ↵
      1. Devlin B,
      2. Roeder K
      . Genomic control for association studies. Biometrics 1999;55:997–1004.
      OpenUrlCrossRefPubMedWeb of Science
    51. ↵
      1. Panagiotou OA,
      2. Ioannidis JP
      , for the Genome-Wide Significance P. What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations. Int J Epidemiol 2012;41:27386.
      OpenUrlAbstract/FREE Full Text
    52. ↵
      1. Johnson AD,
      2. Handsaker RE,
      3. Pulit SL,
      4. et al
      . SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 2008;24:2938–9.
      OpenUrlAbstract/FREE Full Text
    53. ↵
      Genomes Project C. A map of human genome variation from population-scale sequencing. Nature 2010;467:1061–73.
      OpenUrlCrossRefPubMedWeb of Science
    54. ↵
      1. Kumar P,
      2. Henikoff S,
      3. Ng PC
      . Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc 2009;4:1073–81.
      OpenUrlCrossRefPubMedWeb of Science
    55. ↵
      1. Willemen HL,
      2. Eijkelkamp N,
      3. Wang H,
      4. et al
      . Microglial/macrophage GRK2 determines duration of peripheral IL-1beta-induced hyperalgesia: contribution of spinal cord CX3CR1, p38 and IL-1 signaling. Pain 2010;150:550–60.
      OpenUrlCrossRefPubMed
    56. ↵
      1. Zhao J,
      2. Yuan G,
      3. Cendan CM,
      4. et al
      . Nociceptor-expressed ephrin-B2 regulates inflammatory and neuropathic pain. Mol Pain 2010;6:77.
      OpenUrlCrossRefPubMed
    57. ↵
      1. Hargreaves K,
      2. Dubner R,
      3. Brown F,
      4. et al
      . A new and sensitive method for measuring thermal nociception in cutaneous hyperalgesia. Pain 1988;32:77–88.
      OpenUrlCrossRefPubMedWeb of Science
    58. ↵
      1. Yu W,
      2. Clyne M,
      3. Khoury MJ,
      4. et al
      . Phenopedia and Genopedia: disease-centered and gene-centered views of the evolving knowledge of human genetic associations. Bioinformatics 2010;26:145–6.
      OpenUrlAbstract/FREE Full Text
    59. ↵
      1. Schadt EE,
      2. Molony C,
      3. Chudin E,
      4. et al
      . Mapping the genetic architecture of gene expression in human liver. PLoS Biol 2008;6:e107.
      OpenUrlCrossRefPubMed
    60. ↵
      1. Tegeder I,
      2. Costigan M,
      3. Griffin RS,
      4. et al
      . GTP cyclohydrolase and tetrahydrobiopterin regulate pain sensitivity and persistence. Nat Med 2006;12:1269–77.
      OpenUrlCrossRefPubMedWeb of Science
    61. ↵
      1. Campbell CM,
      2. Edwards RR,
      3. Carmona C,
      4. et al
      . Polymorphisms in the GTP cyclohydrolase gene (GCH1) are associated with ratings of capsaicin pain. Pain 2009;141:114–18.
      OpenUrlCrossRefPubMedWeb of Science
    62. ↵
      1. Holliday KL,
      2. Nicholl BI,
      3. Macfarlane GJ,
      4. et al
      . Do genetic predictors of pain sensitivity associate with persistent widespread pain? Mol Pain 2009;5:56.
      OpenUrlCrossRefPubMed
    63. ↵
      1. Hayashida M,
      2. Nagashima M,
      3. Satoh Y,
      4. et al
      . Analgesic requirements after major abdominal surgery are associated with OPRM1 gene polymorphism genotype and haplotype. Pharmacogenomics 2008;9:1605–16.
      OpenUrlCrossRefPubMedWeb of Science
    64. ↵
      1. Kim H,
      2. Lee H,
      3. Rowan J,
      4. et al
      . Genetic polymorphisms in monoamine neurotransmitter systems show only weak association with acute post-surgical pain in humans. Mol Pain 2006;2:24.
      OpenUrlCrossRefPubMed
    65. ↵
      1. Bouhouche A,
      2. Benomar A,
      3. Bouslam N,
      4. et al
      . 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;43:441–3.
      OpenUrlAbstract/FREE Full Text
    66. ↵
      1. Kubota H,
      2. Hynes G,
      3. Willison K
      . The chaperonin containing t-complex polypeptide 1 (TCP-1). Multisubunit machinery assisting in protein folding and assembly in the eukaryotic cytosol. Eur J Biochem 1995;230:3–16.
      OpenUrlPubMedWeb of Science
    67. ↵
      1. Gingras AC,
      2. Caballero M,
      3. Zarske M,
      4. et al
      . A novel, evolutionarily conserved protein phosphatase complex involved in cisplatin sensitivity. Mol Cell Proteomics 2005;4:1725–40.
      OpenUrlAbstract/FREE Full Text
    68. ↵
      1. Chen GI,
      2. Tisayakorn S,
      3. Jorgensen C,
      4. et al
      . PP4R4/KIAA1622 forms a novel stable cytosolic complex with phosphoprotein phosphatase 4. J Biol Chem 2008;283:29273–84.
      OpenUrlAbstract/FREE Full Text
    69. ↵
      1. Glatter T,
      2. Wepf A,
      3. Aebersold R,
      4. et al
      . An integrated workflow for charting the human interaction proteome: insights into the PP2A system. Mol Syst Biol 2009;5:237.
      OpenUrlAbstract/FREE Full Text
    70. ↵
      1. Zhang X,
      2. Ozawa Y,
      3. Lee H,
      4. et al
      . Histone deacetylase 3 (HDAC3) activity is regulated by interaction with protein serine/threonine phosphatase 4. Genes Dev 2005;19:827–39.
      OpenUrlAbstract/FREE Full Text
    71. ↵
      1. Latremoliere A,
      2. Woolf CJ
      . Central sensitization: a generator of pain hypersensitivity by central neural plasticity. J Pain 2009;10:895–926.
      OpenUrlCrossRefPubMedWeb of Science
    72. ↵
      1. Silman AJ,
      2. Pearson JE
      . Epidemiology and genetics of rheumatoid arthritis. Arthritis Res 2002;4(Suppl 3):S265–72.
      OpenUrlCrossRefPubMed
    73. ↵
      1. Bax M,
      2. van Heemst J,
      3. Huizinga TW,
      4. et al
      . Genetics of rheumatoid arthritis: what have we learned? Immunogenetics 2011;63:459–66.
      OpenUrlCrossRefPubMedWeb of Science
    74. ↵
      1. Rolke R,
      2. Baron R,
      3. Maier C,
      4. et al
      . Quantitative sensory testing in the German Research Network on Neuropathic Pain (DFNS): standardized protocol and reference values. Pain 2006;123:231–43.
      OpenUrlCrossRefPubMedWeb of Science
    75. ↵
      1. deCharms RC,
      2. Maeda F,
      3. Glover GH,
      4. et al
      . Control over brain activation and pain learned by using real-time functional MRI. Proc Natl Acad Sci USA 2005;102:18626–31.
      OpenUrlAbstract/FREE Full Text
    76. ↵
      1. Moller AT,
      2. Jensen TS
      . Pain and genes: Genetic contribution to pain variability, chronic pain and analgesic responses. Eur J Pain Supplements 2010;4:197–201.
      OpenUrlCrossRef
    77. ↵
      1. Ioannidis JP
      . Why most published research findings are false. PLoS Med 2005;2:e124.
      OpenUrlCrossRefPubMed
    78. ↵
      1. Richards JB,
      2. Kavvoura FK,
      3. Rivadeneira F,
      4. et al
      . Collaborative meta-analysis: associations of 150 candidate genes with osteoporosis and osteoporotic fracture. Ann Intern Med 2009;151:528–37.
      OpenUrlCrossRefPubMedWeb of Science

    Supplementary materials

    • Supplementary Data

      This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

      Files in this Data Supplement:

    Footnotes

    • MJP, LB, HLDMW are shared first authors. AK, CMD, FMKW, JBJM are shared last authors.

    • Correction notice This article has been corrected since it was published Online First. All occurrences of ‘DSDBAG’ have been corrected to ‘DSDBAC’. The license of this article has also changed since publication to CC BY 4.0.

    • Acknowledgements For acknowledgments, please see the supplements.

    • Contributors Performed the experiments: AH, AMV, AP, AS, AT, AVS, BHS, BO, COS, DJH, EMD, FJPMH, FR, GE, GTJ, HJMK, HLDMW, IM, JM, KAJ, KLH, KW, LB, LJH, LJL, LS, MAH, MJP, MK, MK2, MP, NA, NE, SWW, TBH, TN, WO, WT. Analysed the data: AGU, AH, AK, AMV, AP, AS, AT, AVS, BHS, BO, CC, CMD, COS, DF, DJH, EMD, FJPMH, FMK, FR, GE, GJM, GTJ, HJMK, HLDMW, HV, IM, JBJM, JM, KAJ, KLH, KW, LB, LJH, LJL, LS, MAH, MJP, MK, MK2, MP, NA, NE, NKA, NP, PES, SWW, TBH, TDS, TN, VG, WO, WT. Contributed reagents/materials/analysis tools: AGU, AK, CC, CMD, DF, FMK, GJM, HV, JBJM, NE, NKA, NP, PES, TDS, VG. Wrote the manuscript: AK, AMV, AS, CC, CMD, COS, DF, FMK, GE, GJM, HJMK, HLDMW, HV, IM, JBJM, KLH, LB, LJH, LS, MAH, MJP, MK, NKA, NP, PES, TDS, TN, VG. Final approval of the manuscript: AGU, AH, AK, AMV, AP, AS, AT, AVS, BHS, BO, CC, CMD, COS, DF, DJH, EMD, FJPMH, FMK, FR, GE, GJM, GTJ, HJMK, HLDMW, HV, IM, JBJM, JM, KAJ, KLH, KW, LB, LJH, LJL, LS, MAH, MJP, MK, MK2, MP, NA, NE, NKA, NP, PES, SWW, TBH, TDS, TN, VG, WO, WT.

    • Competing interests HK is currently affiliated at Pfizer. The work conducted for this study was done during her employment at the Erasmus Medical Center. There are no other competing interests.

    • Funding All funders are specified in the acknowledgements section of the supplements.

    • Ethics approval Medical ethics committees of the cohort studies (see supplements).

    • Provenance and peer review Not commissioned; externally peer reviewed.