Skip Navigation


Human Molecular Genetics Advance Access originally published online on November 10, 2004
Human Molecular Genetics 2005 14(1):135-143; doi:10.1093/hmg/ddi013
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Material
Right arrow All Versions of this Article:
14/1/135    most recent
ddi013v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (118)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Diatchenko, L.
Right arrow Articles by Maixner, W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Diatchenko, L.
Right arrow Articles by Maixner, W.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Human Molecular Genetics, Vol. 14, No. 1 © Oxford University Press 2005; all rights reserved

Genetic basis for individual variations in pain perception and the development of a chronic pain condition

Luda Diatchenko1,8,*, Gary D. Slade2, Andrea G. Nackley1, Konakporn Bhalang3, Asgeir Sigurdsson1, Inna Belfer4,7, David Goldman4, Ke Xu4, Svetlana A. Shabalina5, Dmitry Shagin6, Mitchell B. Max7, Sergei S. Makarov8 and William Maixner1

1Comprehensive Center for Inflammatory Disorders, University of North Carolina at Chapel Hill, Columbia Street, CB#7455, Chapel Hill, NC 27599, USA, 2Australian Research Centre for Population Oral Health, University of Adelaide, Frome Road, SA 5005, Adelaide, Australia, 3Department of Oral Medicine, Chulalongkorn University, 254 Phyathai Road, Bangkok 10330, Thailand, 4Laboratory of Neurogenetics, NIAAA, NIH, 12420 Parklawn Drive, Park 5 Building, Rockville, MD 20852, USA, 5NCBI, NIH, 6404 Landon Lane, 8600 Rockville Pike, Bethesda, MD 20894, USA, 6Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry RAS, Ul. Miklukho-Maklaya 16/10, 117997, Moscow, Russia, 7NIDCR, NIH, Pain and Neurosensory Mechanisms Branch, Building 10, Room 3C-405, Bethesda, MD 20892-1258, USA and 8Attagene, Inc., 7030 Kit Creek Road, Research Triangle Park, NC 27560, USA

* To whom correspondence should be addressed at: Comprehensive Center for Inflammatory Disorders, University of North Carolina at Chapel Hill, Columbia Street, CB#7455, 2190 Old Dental Building, Chapel Hill, NC 27599, USA. Tel: +1 9198432549; Fax: +1 9199663683; Email: lbdiatch{at}email.unc.edu

Received September 22, 2004; Accepted November 2, 2004


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
Pain sensitivity varies substantially among humans. A significant part of the human population develops chronic pain conditions that are characterized by heightened pain sensitivity. We identified three genetic variants (haplotypes) of the gene encoding catecholamine-O-methyltransferase (COMT) that we designated as low pain sensitivity (LPS), average pain sensitivity (APS) and high pain sensitivity (HPS). We show that these haplotypes encompass 96% of the human population, and five combinations of these haplotypes are strongly associated (P=0.0004) with variation in the sensitivity to experimental pain. The presence of even a single LPS haplotype diminishes, by as much as 2.3 times, the risk of developing myogenous temporomandibular joint disorder (TMD), a common musculoskeletal pain condition. The LPS haplotype produces much higher levels of COMT enzymatic activity when compared with the APS or HPS haplotypes. Inhibition of COMT in the rat results in a profound increase in pain sensitivity. Thus, COMT activity substantially influences pain sensitivity, and the three major haplotypes determine COMT activity in humans that inversely correlates with pain sensitivity and the risk of developing TMD.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
Pain perception is a complex process that is influenced by a variety of environmental and genetic factors (1Go). Although the relative importance of genetic versus environmental factors in human pain perception remains unclear, reported heritability for nociceptive and analgesic sensitivities in mice is estimated to range from 28 to 76% (1Go). Even though animal studies have provided a list of candidate ‘pain genes’, only a few genes have been identified that are associated with the perception of pain in humans. Single nucleotide polymorphisms (SNPs) in genes that code for melanocortin-1 receptor and neuronal cytochrome P4502D6 are associated with alterations in opioid analgesia in humans (1Go,2Go). Congenital insensitivity to pain (CIP) type I has been linked to the gene encoding a subunit of serine palmitoyltransferase, and CIP type IV has been linked to the gene encoding for a nerve growth factor-specific tyrosine kinase receptor (1Go,3Go). Furthermore, a common SNP in codon 158 (val158met) of the gene that codes for catecholamine-O-methyltransferase (COMT) has been proposed to contribute to differences in the human experience of pain (4Go). This enzyme has broad biological functions including the regulation of the levels of catecholamines and enkephalins (5Go). Thus, the examination of the association of COMT polymorphism with human pain perception and persistent pain conditions is of considerable importance.

A pathological pain condition that appears to be associated with COMT activity is myogenous temporomandibular joint disorder (TMD). This condition is characterized by persistent facial pain, impaired oral function (6Go–8Go) and heightened sensitivity to pain-evoking stimuli (e.g. mechanical, thermal and ischemic) at numerous body sites (9Go,10Go). TMD impacts 5–15% of the adult population and incurs billions of dollars in health care costs (6Go). In 1976, Marbach and Levitt (11Go) reported that patients with facial pain conditions comparable to TMD show increased urinary levels of catecholamine metabolites and express diminished erythrocytic COMT activity, suggesting a role for COMT in this persistent pain condition. In the present study, we examined the relationship between COMT polymorphism, pain sensitivity and the risk of TMD development.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
DISCUSSION
 MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
Pain sensitivity varies significantly in the human population
Data for this study were collected from 202 healthy female volunteers. The subjects participated in a 3 year prospective cohort study that was designed to identify risk factors for TMD (9Go,10Go,12Go). Only females were included in this study as they exhibit a higher prevalence for the condition relative to males (13Go). During an initial screening exam, we assessed the sensitivity of subjects to experimental noxious stimuli and collected peripheral blood samples for genetic analyses. Subjects were subsequently followed for up to 3 years, by both trimonthly interviews and annual physical examinations, to identify newly developed cases of TMD.

In order to evaluate each participant's pain sensitivity, a unique approach was used to derive a unitary measure of pain sensitivity for both cutaneous and deep muscle pain, which are transmitted and modulated by different neural mechanisms (14Go,15Go). To accomplish this, each of 16 measures of pain sensitivity was normalized to a mean of 0 and standard deviation of 1, producing a unit normal deviate (z-score) for each test procedure. A sum of these 16 scores produced a normalized single score of pain sensitivity (integral z-score) for each individual (see Supplemental Material for details). As shown in Figure 1, measures of individual pain sensitivity (integral z-scores) were distributed approximately normally (skewness=0.3 and kurtosis=–0.1), ranging from –22.4 (least responsive to painful stimuli) to 28.0 (most responsive to painful stimuli); although the majority of individuals display average pain sensitivity, the individual variability in pain sensitivity between people is substantial.



View larger version (34K):
[in this window]
[in a new window]
  		Figure 1. Distribution of a summary measure of pain sensitivity. A summary measure of pain sensitivity was derived from 16 individual pain measures, each standardized to unit normal deviates (z-scores) with a mean of 0 and standard deviation of 1. The 16 pain measures were thermal pain threshold conveyed by A{delta} afferents and both threshold and tolerance conveyed by C-fiber afferents, all measured in °C at each of three anatomical sites (arm, cheek and foot); tolerance to temporal summation of C-fiber mediated pain (as reported on 0–100 visual analog scale); right arm ischemic pain onset and tolerance (s); and mechanical pain thresholds (kg) assessed over the temporalis and masseter muscles, the temporomandibular joint and the ventral surfaces of wrists. Individuals represented at the extreme left-side of the figure are resistant to pain-evoking procedures whereas individuals represented at the extreme right-side of the figure are most sensitive to pain-evoking procedures.

 
Genotyping of the COMT locus
Genomic DNA from peripheral blood samples was genotyped for SNPs within the COMT gene locus. Six SNPs were chosen that display high polymorphism frequency in the human population (>40% prevalence). Figure 2A shows the positions of the SNPs within the COMT locus that code for two major forms of COMT enzyme: membrane bound (MB-COMT) and soluble (S-COMT). The first SNP (rs2097903) is located at position –1217 in the estrogen sensitive portion of the MB-COMT promoter region (16Go,17Go), whereas the second SNP (rs6269) is located in the promoter region of S-COMT (16Go,18Go). The next three SNPs (rs4633, rs4818 and rs4680) are located within the coding region for both S- and MB-COMT (19Go) (NCBI genome database). Variations in SNPs rs4633 and rs4818 are synonymous (i.e. do not produce a change in amino acid composition). In contrast, SNP rs4680 is non-synonymous and codes for a substitution of valine (val) to methionine (met) at codon 158. The last SNP (rs165599) is located in the very end of the 3'-UTR of the gene and its G allele has been reported to be associated with schizophrenia (18Go). Importantly, the COMT SNP map has been thoroughly constructed (NCBI database, http://www.ncbi.nlm.nih.gov/SNP/snp_ref.cgi?locusId=1312). In the coding region of the gene, there are no other SNPs with frequencies greater than 0.15. We further confirmed this by comparing more than 300 COMT EST sequences from NCBI database using the CLUSTALW program for multiple sequence alignments. As it has been discussed by Risch (20Go), alleles with polymorphic frequencies <0.15% are very unlikely to have significant impact on common population-based diseases/disorders such as TMD.



View larger version (18K):
[in this window]
[in a new window]
  		Figure 2. (A) Schematic diagram of COMT genomic organization, SNP positions and percentage distribution. COMT gene locus (chromosome 22, band q11.21) spans for 27 221 bases. (B) Schematic diagram of linkage disequilibrium between six SNP markers. Four SNPs, rs6269, rs4633, rs4818 and rs4680 (val158met), which occur within the central region of the COMT gene, were found to exhibit strong LDs with the strongest associations found between SNPs rs6269 and rs4818 (D'=0.94, R2=0.88) and between SNPs rs4633 and rs4680 (D'=0.96, R2=0.91). In contrast, SNP rs2097903, located in the 5' promoter region, and SNP rs165599, located in the 3'-UTR, did not show strong LDs (Supplementary Material, Table S4). Thus, LD analysis demonstrated that the COMT locus covers three haploblocks (B). SNP rs2097903 is situated on the first haploblock, SNPs rs6269, rs4633, rs4818 and rs4680 are situated on the second haploblock and SNP rs165599 is situated on the third haploblock. (C) Estimated frequencies of the COMT haplotypes. The sequence of alleles in each haplotype for haploblock 2 reflects the order of occurrence from 5' to 3' in the COMT gene (SNPs: rs6269, rs4633, rs4818 and rs4680, respectively). Seven haplotypes out of possible 16 were detected for these four SNPs with the most frequent haplotype (48.7%) composed of the most frequent alleles for all four markers (A_T_C_A for SNPs rs6269, rs4633, rs4818 and rs4680, respectively). The second major haplotype (36.5%) was composed of the least frequent alleles for all four markers (G_C_G_G). The third haplotype (10.5%) was composed of a combination of the most frequent alleles for SNPs rs4633 and rs4680 and the least frequent alleles for SNPs rs6269 and rs4818 (A_C_C_G). These three haplotypes accounted for 95.9% of all detected haplotypes. These findings are consistent with the previously reported LD analysis of the COMT SNPs (18Go,19Go,28Go) and with the report of Gabriel et al. (21Go) that each haploblock within human genomic DNA is usually represented by three to five major haplotypes.

 
Statistically significant associations were found between the integral z-score and two SNPs (Table 1; Supplementary Material, Fig. S1). The SNP rs6269 accounted for 6% of the variation in pain sensitivity as determined by analysis of variance (ANOVA, P<0.01), whereas SNP rs4818 accounted for 7% of the variation (ANOVA, P<0.01). For both SNPs, the homozygous genotypes were associated with significant differences in mean pain sensitivity (t-test, P<0.01). Consistent with previously reported results (4Go), the val158met SNP (rs4680) showed a marginal, but not statistically significant, relationship with pain sensitivity, accounting for 2% of the variation in the summary pain measure (ANOVA, P=0.18). Individuals homozygous for met/met tended to be more pain responsive than those homozygous for val/val, but again the association was marginal (t-test, P=0.06). SNPs A- rs2097903, rs4633 and rs16559 were not significantly associated with pain sensitivity (Table 1; Supplementary Material, Fig. S1). Multivariate analysis revealed that all possible combinations of the four SNPs in the coding region accounted for 10.6% of the variation in the summary measure of pain sensitivity. In this model, we first controlled for rs4680 (accounting for 2% of variance), then for rs4818 (7% of variance); the remaining SNPs did not individually contribute significantly (P>0.10) to the variance in pain sensitivity (Supplementary Material, Table S3).


View this table:
[in this window]
[in a new window]
  		Table 1. Variation in pain sensitivity (summed z-score) among tested SNPs and diplotypes of the COMT gene
 
COMT haplotypes determine sensitivity to pain
We next determined which combinations of alleles (haplotypes) were formed by the six SNPs. It has been shown that alleles form associations (haploblocks) of variable length with the average span of 18 kb in populations of European descent and only a few common haplotypes are observed (21Go). Three haploblocks were determined in our study sample by LD analysis (Supplementary Material, Table S4). Because the association with pain sensitivity was observed only for SNPs rs6269 and rs4818, located in the central COMT locus haploblock (Table 1, Fig. 2B; Supplementary Material, Fig. S1), we focused our analyses on this haploblock.

Seven haplotypes with a frequency >0.5% were detected, three of them representing 95.9% of all haplotypes observed in this study (Fig. 2C). Five combinations of these three haplotypes were present in 92% of subjects and were associated with marked gradients in pain responsiveness (Fig. 3). Subjects homozygous for the G_C_G_G haplotype had the lowest pain responsiveness (mean summed z-score=–5.23±1.5; Fig. 3, Table 1); thus, G_C_C_G is designated as the ‘low pain sensitivity’ (LPS) haplotype. Intermediate pain responsiveness was observed in individuals homozygous for A_T_C_A, which we refer to as the APS haplotype (mean summed z-score=1.75±1.47; Fig. 3, Table 1). The greatest pain responsiveness was observed in individuals heterozygous for A_T_C_A (APS) and A_C_C_G haplotypes (mean summed z-score=8.9±2.9; Fig. 3, Table 1). We refer to the A_C_C_G haplotype as the ‘high pain sensitivity’ (HPS) haplotype. Differences among the five combinations of haplotypes were significant (Table 1, overall ANOVA, P=0.0004) and factorial analysis demonstrated that each haplotype had independent effects on pain sensitivity (Table 2, factorial ANOVA model, P≤0.01). These haplotypes accounted for 10.4% of the variation (P<0.01) in pain sensitivity, representing virtually all of the variation (10.6%) explained by combinations of the four individual SNPs in the central haploblock.



View larger version (23K):
[in this window]
[in a new window]
  		Figure 3. Pain responsiveness categorized by three major COMT haplotype combinations. LPS: haplotype G_C_G_G, APS: haplotype A_T_C_A, HPS: haplotype A_C_C_G. The greater values reflect greater pain sensitivity. Each value represents the mean z-score with associated SEM.

 

View this table:
[in this window]
[in a new window]
  		Table 2. Factorial ANOVA model of effects of three major haplotypes on pain sensitivitya
 
COMT haplotypes determine enzyme activity
Functional polymorphism in the COMT gene has been described only for SNP rs4680 (val158met) (4Go,5Go,22Go), which codes for a substitution of val to met. The met substitution produces a COMT enzyme with lower thermostability, resulting in decreased enzyme activity (22Go). It is generally accepted that genetic variability in codon 158 is the primary source of individual variation in COMT activity in humans (5Go); however, in many cases, the observed association between several human diseases and disorders and the met allele (low enzyme activity) is relatively modest and occasionally opposite, such as for obsessive-compulsive disorder (reviewed in 23Go) and schizophrenia (reviewed in 24Go).

Although it appears that the val158met amino acid substitution in COMT can explain the greater pain responsiveness observed for individuals with the APS haplotype when compared with the LPS haplotype, because the APS haplotype codes for the less stable met variant, it cannot explain the greater pain responsiveness that we found for subjects with the HPS haplotype when compared with the LPS haplotype. This is because both the HPS and LPS haplotypes possess the G allele that codes for the more stable val variant (Fig. 2C). Thus, the val158met SNP alone cannot account for the observed variations in pain perception. Furthermore, even though polymorphism in SNPs rs6269 and rs4818 is significantly associated with pain z-scores, both pain sensitive haplotypes HPS and APS contain the A allele of rs6269 and the C allele of rs4818. Consequently, variations in these SNPs cannot explain why the HPS and APS haplotypes are associated with different levels of pain sensitivity. Instead, the interaction of the val158met SNP with other SNPs determines the functional outcomes. The other SNPs are either synonymous (i.e. code for the same amino acid) or are located in the promoter region of S-COMT. A haplotype-dependent regulation of mRNA expression is also unlikely because SNP rs6269, which is located in the promoter region of S-COMT, does not independently contribute to pain sensitivity (Supplementary Material, Table S3). Therefore, we hypothesized that haplotype-specific secondary structures of mRNA affect COMT mRNA stability and/or efficiency of protein translation (25Go).

To test these possibilities, we transiently transfected HEK 293 cells with full-length COMT cDNA clones that corresponded to the three major haplotypes. The expression of COMT protein was assessed by measuring COMT enzymatic activity in the lysate of transfected cells. Figure 4A shows that the LPS haplotype provides 4.8 times higher COMT activity when compared with the APS haplotype (P<0.01), which is in agreement with a previous report of reduced thermal stability of COMT protein containing met (APS haplotype) in comparison with COMT protein containing val (LPS haplotype) (22Go). However, our finding that the HPS haplotype provides 11.4 times lower COMT activity compared with the LPS haptotype (P<0.01) can only be attributed to the lower amount of protein produced by the HPS haplotype as these two haplotypes code for COMT protein with exactly the same amino acid composition. No differences in COMT RNA abundance were detected in the transfected cells as measured by real-time PCR (Fig. 4B); therefore, the three major haplotypes affect the efficiency of protein synthesis, but not RNA stability.



View larger version (21K):
[in this window]
[in a new window]
  		Figure 4. Effect of haplotype on the COMT activity in transfected HEK 293 cells. HEK 293 cells were transiently transfected with six full-length COMT cDNA clones that corresponded to the three major haplotypes (LPS, APS and HPS). The expression of COMT protein was assessed by measurement of COMT enzymatic activity (A) in the lysate from transfected cells. COMT activity was assessed from two independent clones for each haplotype. Two independent ELISAs for each of the two independent transfections were performed for each clone. The graph shows the average values for the eight measurements for each haplotype. The COMT activity was calculated in nanograms of NMN, synthesized during 1 h at 37°C under described enzymatic condition per 105 transfected cells. The highest COMT activity was observed for LPS haplotypes ([NMN]=34.3±3.0 ng per 105 cells; A), the lower activity was observed for APS haplotype ([NMN]=7.2±1.53 ng per 105 cells; A) and HPS haplotypes produced the lowest activity ([NMN]=3.0±2.20 ng per 105 cells; A). The relative abundance of COMT RNA (B) was assessed by real-time PCR. Each value represents the mean value with associated SEM.

 
Inhibition of COMT enhances sensitivity to noxious stimuli
Our association and in vitro studies strongly suggest that reductions in COMT enzymatic activity enhances pain sensitivity. To directly test whether decreased COMT activity enhances pain sensitivity, the COMT inhibitor OR486 was administered to naive rats. OR486 decreased paw withdrawal thresholds to mechanical and thermal stimuli (P<0.0001 and P<0.0007, respectively) and increased paw withdrawal frequency to noxious punctate mechanical stimuli (P<0.0001) (Fig. 5). The degree of mechanical and thermal hyperalgesia produced by OR486 was comparable to that produced by carrageenan-induced inflammation in the hindpaw (data not shown).



View larger version (14K):
[in this window]
[in a new window]
  		Figure 5. The effect of COMT inhibition on rat pain behavior. Baseline sensitivity to mechanical stimuli was estimated by measuring both threshold (A) and frequency (B) of paw withdrawal. Sensitivity to thermal stimuli was estimated by measuring latency of paw withdrawal (C). After establishing baseline sensitivity, animals received the COMT inhibitor OR486 (30 mg/kg i.p.) or vehicle 1 h prior to testing. Data are expressed as mean±SEM. ***P<0.001, **P<0.01 and *P<0.05 different from the control conditions by ANOVA and Bonferroni post hoc tests. N=8 rats per group.

 
The high activity COMT haplotype protects from developing TMD
To determine the clinical relevance of these findings, we examined whether COMT polymorphism is related to the incidence of TMD onset among 170 subjects with the five most common haplotype combinations who completed one or more follow-up visits. Fifty-eight of the participants had only ‘low COMT activity’ haplotypes (HPS and/or APS) and the remaining 112 subjects had at least one ‘high activity’ haplotype (LPS). We first confirmed that HPS and/or APS subjects were more sensitive to experimental pain at their baseline assessment when compared with LPS subjects (P=0.02; Fig. 6A). During the 3 year observational period, 15 new cases of TMD were diagnosed (12Go) at varying time periods ranging from 9 months to 3 years after recruitment, yielding an average incidence rate of 3.5 cases per 100 person-years of follow-up. The incidence rate was more than twice as high among individuals having only HPS and/or APS haplotypes (5.6 cases per 100 person-years) when compared with individuals with at least one LPS haplotype (2.5 cases per 100 person-years, Fig. 6B). The derived incidence density ratio of 2.3 was significant (95% confidence interval=1.1–4.8), suggesting that the HPS and/or APS haplotypes represent significant risk factors for TMD onset.



View larger version (9K):
[in this window]
[in a new window]
  		Figure 6. Pain sensitivity and TMD incidence by haplotype groupings: HPS/APS versus LPS. Because of the limited number of TMD cases in the cohort, all subjects were subdivided into two groups: LPS and HPS/APS. Subjects were assigned to the LPS group if they carried at least one LPS haplotype. Subjects were assigned to the HPS/APS group if they carried only APS and HPS haplotypes. (A) Shows that subjects from the LPS group demonstrated significantly lower pain responsiveness than those from the HPS/APS group (P=0.02, t-test). (B) Shows the number of incidence case per 100 person-years as a function of haplotype group.

 

   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
Several important conclusions can be drawn. First, COMT genotype is highly associated with human pain perception. There are three major COMT haplotypes (LPS, APS and HPS) that determine COMT enzymatic activity, encompassing ~96% of the examined genotypes. The LPS haplotype is associated with low pain sensitivity, APS is associated with higher pain sensitivity and HPS with the highest pain sensitivity. Collectively, these three haplotypes account for ~11% of the variability in pain perception. Given the inevitably polygenic nature of pain perception, the magnitude of the effect of COMT haplotypes on pain perception is substantial. Indeed, quantitative trait locus (QTL) mapping studies for related traits in mice have shown that each single QTL usually accounts for 5–25% of the overall variance in nociceptive sensitivity (2Go,26Go).

Secondly, synonymous SNPs within COMT haplotypes can have effects on protein function that exceed the effects of individual SNPs. Several recently published studies support the theory that the collective grouping of SNPs in haplotypes has a stronger association with the assessed phenotype than individual SNPs (18Go,25Go,27Go–29Go). However, almost all previously published haplotype association studies demonstrate the importance of haplotype reconstruction because combinations of SNPs result in synergistic effects on protein function. Recently, Duan et al. (25Go) showed that synonymous SNPs within haplotypes can have functional consequences drastically different from those of each isolated mutation. This effect was attributed to alternations in the secondary structure of mRNA, which results in alterations in mRNA stability. In contrast, our studies show that genomic variations of the COMT gene do not alter the amount of COMT mRNA, suggesting that the differences in enzymatic activity result from differences in protein translation. The fact that expressed cDNA constructs, which differed in only three SNPs rs4633, rs4818 and rs4680 (val158met), showed >11-fold difference in expressed enzyme activity, confirms that the observed association between haplotypes and pain sensitivity is largely caused by combinations of these three SNPs and not by other SNPs in the haploblock located in the 5' or intronic region of the COMT gene that can affect RNA transcription. Although the precise molecular mechanisms underlying this phenomenon are unknown, our ongoing studies suggest that interactions between SNPs have profound effects on the secondary mRNA structure, which controls the efficacy of protein translation. The identification of new functional haplotypes in the present study reinforces the findings presented by Shifman et al. (18Go) and Bray et al. (28Go), suggesting that haplotype reconstruction, rather than polymorphism at codon 158, needs to be considered in future studies that examine the relationship between COMT polymorphism, human pain and pain-related or affective disorders.

Thirdly, COMT inhibition results in a robust increase in pain sensitivity. These data provide evidence that COMT activity regulates pain sensitivity and strongly suggests that the observed association between COMT genotype and pain perception is unlikely to be epiphenomenal or caused by other functional polymorphisms in neighboring genes that possess high LDs with the tested COMT SNPs.

Finally, our results are of considerable clinical significance and represent the first study of its kind to demonstrate an association between a genetic polymorphism that impacts pain sensitivity and the risk for a clinical chronic pain condition. The presence of even a single high COMT activity (LPS) haplotype diminishes, by as much as 2.3 times, the risk of developing TMD, a common musculoskeletal pain condition. The risk ratio of 2.3 is of a magnitude comparable to genetic risk factors for other multifactorial conditions, such as schizophrenia (18Go), and is similar to other predictors of TMD, such as a history of chronic pain at other body sites (7Go,8Go). The clinical relevance of this finding is best quantified by the measure of population attributable risk for having HPS and/or APS, which was 29% in this cohort of women, indicating that nearly one-third of new TMD cases can be attributed to this COMT genotype.

The mechanism by which diminished COMT activity influences pain perception and the development of TMD is not known. It has been suggested that reduction in COMT activity associated with the met allele at codon 158 leads to a reduction in the content of enkephalins (endogenous opioid-like peptides) in certain regions of the CNS associated with pain and mood (4Go). Another possible mechanism is that reduced COMT activity results in elevated levels of catecholamimes, such as epinephrine, which promote the production of persistent pain states via the stimulation of ß2-adrenergic receptors in the peripheral and central nervous system (30Go). Future studies are required to establish the precise mechanisms. However, the clinical, animal and molecular data presented in this study are in complete agreement with the conclusion that COMT activity substantially influences pain sensitivity and that the three major haplotypes determine COMT activity in humans in a fashion that inversely correlates with pain sensitivity and the risk of developing TMD and possibly other chronic pain conditions.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
SUPPLEMENTARY MATERIAL
 REFERENCES
 
Detailed descriptions for each procedure are provided in Supplementary Material.

Clinical assessments and procedures
Study participants consisted of 202 healthy pain-free females aged 18–34 years who provided a blood sample and consent for genotyping from among a larger cohort of 244 females. All subjects underwent pain perception assessments and head/neck examinations. Pressure pain thresholds were assessed over the right and left temporalis muscles, masseter muscles, temporomandibular joints and ventral surfaces of the wrists with a hand-held pressure algometer (31Go). Thermal pain thresholds and tolerances (in °C) were determined with a computer-controlled thermal stimulator on the skin overlying the right masseter muscle, the skin overlying the right hairy forearm and the skin overlying the dorsal surface of the right foot (32Go). The temporal summation of heat pain was determined by delivering 15 heat pulses of 53°C to the skin overlying the thenar region of the right hand (33Go). Subjects verbally rated the intensity of each thermal pulse using a numerical scale (34Go). Ischemic pain thresholds and tolerances (in min) were assessed via the submaximal effort tourniquet procedure (9Go).

Genotyping
Genomic DNA was purified from 202 subjects using QIAampTM 96 DNA Blood Kit (Qiagen, Valencia, CA, USA) and used for either 5'exonuclease TaqManTM PCR for the first five COMT SNPs (rs2097903, rs6269, rs4633, rs4818, rs4680) (35Go) or duplex-specific nuclease assays (36Go) for SNP rs165599. The NCBI SNP databases were used to assign SNP numbers. SAS proc haplotype was used for haplotype reconstruction.

Assessment of COMT activity corresponding to different haplotypes
Full-length S-COMT cDNA clones corresponding to LPS, APS or HPS haplotypes were obtained from the IMAGE clone collection (Open Biosystems, Huntsville, AL, USA). Clones BG290167 and BG818517 represented LPS haplotype, clones BI821094 and F037202 represented APS haplotype and clones BI759217 and BF035214 represented HPS haplotype. Human embryonic kidney cells (HEK 293) were transiently transfected using SuperFect Reagent (Qiagen). To control the efficiency of transfection, pSV-ßGalactosidase vector (Promega, Madison, WI, USA) was cotransfected with the COMT-containing clones. Cells lysates were collected ~24 h post-transfection. The COMT enzymatic assay was based on the method described by Masuda et al. (37Go). Normetanephrine (NMN) amount was measured with a NMN ELISA kit [Immuno-Biological Laboratories, Inc. (IBL-America), Minneapolis, MN, USA]. COMT activity was then normalized for transfection efficiency by measuring the ß-galactosidase activity from each lysate using a ß-galactosidase enzyme system (Promega). RNA abundance was measured using a Opticon-2 Real Time Fluorescence Detection System (MJ Research, Reno, NV, USA) with a DyNAmo-SYBRGreen qPCR kit (MJ Research).

Statistical analysis of clinical data
Sixteen measures of pain sensitivity were standardized to unit normal deviates (z-scores) and summed to produce one summary measure of pain responsiveness for each subject. Associations with each of the six SNPs were evaluated for 202 genotyped subjects using ANOVA and Student's t-test (Table 1; Supplementary Material, Fig. S1). Independent effects of SNPs were evaluated in a multivariable ANOVA model. Associations between pain responsiveness and haplotypes were assessed in a factorial ANOVA model for 186 subjects carried one of five combinations of the three most prevalent haplotypes. Incidence of newly diagnosed TMD was computed for 170 of those subjects who had one or more follow-up evaluations over the 3 year study period. Incidence rates for two haplotype groups were computed as the number of TMD cases divided by the person-years of follow-up and compared using Poisson regression with the log of person-years as the offset.

Animal behavior and analysis
Sixteen adult male Sprague–Dawley rats (285–325 g; Charles River Laboratories, Wilmington, MA, USA) were placed in plexiglass cages positioned over an elevated perforated stainless-steel platform and habituated to the environment for 15–25 min prior to testing. Paw withdrawal threshold to punctate mechanical stimulation was assessed using the up–down method of Chaplan et al. (38Go). Immediately following determination of the response threshold, paw withdrawal frequency (%) to the presentations of a noxious von Frey monofilament was calculated. Thermal hyperalgesia was evaluated using the Hargreaves radiant heat method (39Go). After establishing stable baseline responsiveness to mechanical and thermal stimuli, separate groups of rats received OR486 (30 mg/kg i.p.; N=8) or vehicle (N=8) 1 h prior to behavioral testing. Responsiveness to von Frey filaments was reassessed at 30 min intervals for 2 h. Paw withdrawal latencies to radiant heat were subsequently assessed at 2.5 h into the testing procedure. Behavioral data were analyzed by ANOVA for repeated measures and post hoc comparisons were performed using Bonferroni test.


   SUPPLEMENTARY MATERIAL
TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
REFERENCES
 
Supplementary Material is available at HMG Online.


   ACKNOWLEDGEMENTS
 
We would like to thank Ollie Monbureau for developing the hardware and software used for pain phenotyping. We also want to thank Julie A. Clarke for performing the genomic DNA purification. Finally, we would also like to express our thanks to Jurgen Westermann from IBL-Hamburg, Germany, for assisting us with optimizing the Normetanephrine ELISA kit for COMT measurements in cell lysates. Supported by DE07509 and DE007333 (W.M.), AR/AI-44564, 5-P60 AR-30701-14 and AR/AI-44030, NIH Intramural Grants DE00366 and AA000301, the Comprehensive Neuroscience Program Grant USUHS G192BR-C4 and by Attagene's R&D funding.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 SUPPLEMENTARY MATERIAL
 REFERENCES
 

  1. Mogil, J.S. (1999) The genetic mediation of individual differences in sensitivity to pain and its inhibition. Proc. Natl Acad. Sci. USA, 96, 7744–7751.[Abstract/Free Full Text]

  2. Mogil, J.S., Wilson, S.G., Chesler, E.J., Rankin, A.L., Nemmani, K.V., Lariviere, W.R., Groce, M.K., Wallace, M.R., Kaplan, L., Staud, R. et al. (2003) The melanocortin-1 receptor gene mediates female-specific mechanisms of analgesia in mice and humans. Proc. Natl Acad. Sci. USA, 100, 4867–4872.[Abstract/Free Full Text]

  3. Nagasako, E.M., Oaklander, A.L. and Dworkin, R.H. (2003) Congenital insensitivity to pain: an update. Pain, 101, 213–219.[CrossRef][ISI][Medline]

  4. Zubieta, J.K., Heitzeg, M.M., Smith, Y.R., Bueller, J.A., Xu, K., Xu, Y., Koeppe, R.A., Stohler, C.S. and Goldman, D. (2003) COMT val158met genotype affects mu-opioid neurotransmitter responses to a pain stressor. Science, 299, 1240–1243.[Abstract/Free Full Text]

  5. Mannisto, P.T. and Kaakkola, S. (1999) Catechol-O-methyltransferase (COMT): biochemistry, molecular biology, pharmacology and clinical efficacy of the new selective COMT inhibitors. Pharmacol. Rev., 51, 593–628.[Abstract/Free Full Text]

  6. National Ambulatory Medical Care Survey, Report of the Panel on Communicative Disorders and Stroke Council. No: 81-1914. 6-1-1979. Public Health Service, NIH, Washington, D.C.

  7. John, M.T., Miglioretti, D.L., LeResche, L., Von Korff, M. and Critchlow, C.W. (2003) Widespread pain as a risk factor for dysfunctional temporomandibular disorder pain. Pain, 102, 257–263.[CrossRef][ISI][Medline]

  8. Von Korff, M., Le Resche, L. and Dworkin, S.F. (1993) First onset of common pain symptoms: a prospective study of depression as a risk factor. Pain, 55, 251–258.[CrossRef][ISI][Medline]

  9. Maixner, W., Fillingim, R., Booker, D. and Sigurdsson, A. (1995) Sensitivity of patients with painful temporomandibular disorders to experimentally evoked pain. Pain, 63, 341–351.[CrossRef][ISI][Medline]

  10. Maixner, W., Fillingim, R., Sigurdsson, A., Kincaid, S. and Silva, S. (1998) Sensitivity of patients with painful temporomandibular disorders to experimentally evoked pain: evidence for altered temporal summation of pain. Pain, 76, 71–81.[CrossRef][ISI][Medline]

  11. Marbach, J.J. and Levitt, M. (1976) Erythrocyte catechol-O-methyltransferase activity in facial pain patients. J. Dent. Res., 55, 711.[Free Full Text]

  12. Dworkin, S.F., Fricton, J.R., Hollender, L., Huggins, K.H., LeResche, L., Lund, J., Mohl, N., Ohrbach, R., Palla, S.F., Sommers, E.E. et al. (1992) Research diagnostic criteria for temporomandibular disorders: review, criteria, examinations and specifications, critique. J. Craniomandibular Disord. Facial Pain Oral Pain, 6, 302–355.

  13. Carlsson, G.E. and Le Resche, L. (1995) In Sessle, B.J., Bryant, P.S. and Dionne, R.A. (eds), Temporomandibular Disorders and Related Pain Conditions. IASP Press, Seattle, Vol. 4, pp. 211–226.

  14. Yu, X.-M., Hua, M. and Mense, S. (1991) The effects of intracerebroventricular injection of naloxone, phentolamine, methysergide on the transmission of nociceptive signals in the rat dorsal horn neurones with convergent cutaneous-deep input. Neuroscience, 44, 715–723.[CrossRef][ISI][Medline]

  15. Mense, S. (1993) Nociception from skeletal muscle in relation to clinical muscle pain. Pain, 54, 241–289.[CrossRef][ISI][Medline]

  16. Xie, T., Ho, S.L. and Ramsden, D. (1999) Characterization and implications of estrogenic down-regulation of human catechol-O-methyltransferase gene transcription. Mol. Pharmacol., 56, 31–38.[Abstract/Free Full Text]

  17. DeMille, M.M., Kidd, J.R., Ruggeri, V., Palmatier, M.A., Goldman, D., Odunsi, A., Okonofua, F., Grigorenko, E., Schulz, L.O., Bonne-Tamir, B. et al. (2002) Population variation in linkage disequilibrium across the COMT gene considering promoter region and coding region variation. Hum. Genet., 111, 521–537.[CrossRef][ISI][Medline]

  18. Shifman, S., Bronstein, M., Sternfeld, M., Pisante-Shalom, A., Lev-Lehman, E., Weizman, A., Reznik, I., Spivak, B., Grisaru, N., Karp, L. et al. (2002) A highly significant association between a COMT haplotype and schizophrenia. Am. J. Hum. Genet., 71, 1296–1302.[CrossRef][ISI][Medline]

  19. Li, T., Ball, D., Zhao, J., Murray, R.M., Liu, X., Sham, P.C. and Collier, D.A. (2000) Family-based linkage disequilibrium mapping using SNP marker haplotypes: application to a potential locus for schizophrenia at chromosome 22q11. Mol. Psychiatry, 5, 452.[ISI][Medline]

  20. Risch, N.J. (2000) Searching for genetic determinants in the new millennium. Nature, 405, 847–856.[CrossRef][Medline]

  21. Gabriel, S.B., Schaffner, S.F., Nguyen, H., Moore, J.M., Roy, J., Blumenstiel, B., Higgins, J., DeFelice, M., Lochner, A., Faggart, M. et al. (2002) The structure of haplotype blocks in the human genome. Science, 296, 2225–2229.[Abstract/Free Full Text]

  22. Lotta, T., Vidgren, J., Tilgmann, C., Ulmanen, I., Melen, K., Julkunen, I. and Taskinen, J. (1995) Kinetics of human soluble and membrane-bound catechol O-methyltransferase: a revised mechanism and description of the thermolabile variant of the enzyme. Biochemistry, 34, 4202–4210.[CrossRef][Medline]

  23. Azzam, A. and Mathews, C.A. (2003) Meta-analysis of the association between the catecholamine-O-methyl-transferase gene and obsessive-compulsive disorder. Am. J. Med. Genet., 123B, 64–69.

  24. Glatt, S.J., Faraone, S.V. and Tsuang, M.T. (2003) Association between a functional catechol O-methyltransferase gene polymorphism and schizophrenia: meta-analysis of case–control and family-based studies. Am. J. Psychiatry, 160, 469–476.[Abstract/Free Full Text]

  25. Duan, J., Wainwright, M.S., Comeron, J.M., Saitou, N., Sanders, A.R., Gelernter, J. and Gejman, P.V. (2003) Synonymous mutations in the human dopamine receptor D2 (DRD2) affect mRNA stability and synthesis of the receptor. Hum. Mol. Genet., 12, 205–216.[Abstract/Free Full Text]

  26. Abiola, O., Angel, J.M., Avner, P., Bachmanov, A.A., Belknap, J.K., Bennett, B., Blankenhorn, E.P., Blizard, D.A., Bolivar, V., Brockmann, G.A. et al. (2003) The nature and identification of quantitative trait loci: a community's view. Nat. Rev. Genet., 4, 911–916.[ISI][Medline]

  27. Drysdale, C.M., McGraw, D.W., Stack, C.B., Stephens, J.C., Judson, R.S., Nandabalan, K., Arnold, K., Ruano, G. and Liggett, S.B. (2000) Complex promoter and coding region beta 2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc. Natl Acad. Sci. USA, 97, 10483–10488.[Abstract/Free Full Text]

  28. Bray, N.J., Buckland, P.R., Williams, N.M., Williams, H.J., Norton, N., Owen, M.J. and O'Donovan, M.C. (2003) A haplotype implicated in schizophrenia susceptibility is associated with reduced COMT expression in human brain. Am. J. Hum. Genet., 73, 152–161.[CrossRef][ISI][Medline]

  29. Davidson, S. (2000) Research suggests importance of haplotypes over SNPs. Nat. Biotechnol., 18, 1134–1135.[CrossRef][ISI][Medline]

  30. Khasar, S.G., Green, P.G., Miao, F.J. and Levine, J.D. (2003) Vagal modulation of nociception is mediated by adrenomedullary epinephrine in the rat. Eur. J. Neurosci., 17, 909–915.[CrossRef][ISI][Medline]

  31. Jaeger, B. and Reeves, J.L. (1986) Quantification of changes in myofascial trigger point sensitivity with the pressure algometer following passive stretch. Pain, 27, 203–210.[CrossRef][ISI][Medline]

  32. Fruhstorfer, H., Lindblom, U. and Schmidt, W.G. (1976) Method for quantitative estimation of thermal thresholds in patients. J. Neurol. Neurosurg. Psychiatry, 39, 1071–1075.[Abstract]

  33. Price, D.D., Hu, J.W., Dubner, R. and Gracely, R.H. (1977) Peripheral suppression of first pain and central summation of second pain evoked by noxious heat pulses. Pain, 3, 57–68.[CrossRef][ISI][Medline]

  34. Vierck, C.J., Jr, Cannon, R.L., Fry, G., Maixner, W. and Whitsel, B.L. (1997) Characteristics of temporal summation of second pain sensations elicited by brief contact of glabrous skin by a preheated thermode. J. Neurophysiol., 78, 992–1002.[Abstract/Free Full Text]

  35. Shi, M.M., Bleavins, M.R. and de la Iglesia, F.A. (1999) Technologies for detecting genetic polymorphisms in pharmacogenomics. Mol. Diagn., 4, 343–351.[CrossRef][ISI][Medline]

  36. Shagin, D.A., Rebrikov, D.V., Kozhemyako, V.B., Altshuler, I.M., Shcheglov, A.S., Zhulidov, P.A., Bogdanova, E.A., Staroverov, D.B., Rasskazov, V.A. and Lukyanov, S. (2002) A novel method for SNP detection using a new duplex-specific nuclease from crab hepatopancreas. Genome Res., 12, 1935–1942.[Abstract/Free Full Text]

  37. Masuda, M., Tsunoda, M., Yusa, Y., Yamada, S. and Imai, K. (2002) Assay of catechol-O-methyltransferase activity in human erythrocytes using norepinephrine as a natural substrate. Ann. Clin. Biochem., 39, 589–594.[CrossRef][ISI][Medline]

  38. Chaplan, S.R., Pogrel, J.W. and Yaksh, T.L. (1994) Role of voltage-dependent calcium channel subtypes in experimental tactile allodynia. J. Pharmacol. Exp. Ther., 269, 1117–1123.[Abstract/Free Full Text]

  39. Hargreaves, K., Dubner, R., Brown, F., Flores, C. and Joris, J. (1988) A new and sensitive method for measuring thermal nociception in cutaneous hyperalgesia. Pain, 32, 77–88.[CrossRef][ISI][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 Eur J OrthodHome page
L. Sonnesen and P. Svensson
Temporomandibular disorders and psychological status in adult patients with a deep bite
Eur J Orthod, August 5, 2008; (2008) cjn044v1.
[Abstract] [Full Text] [PDF]

Home page
 Behav EcolHome page
M. Ruiz, E. Davis, and E. P. Martins
Courtship attention in sagebrush lizards varies with male identity and female reproductive state
Behav. Ecol., July 8, 2008; (2008) arn072v1.
[Abstract] [Full Text] [PDF]

Home page
 Br J AnaesthHome page
W. A. Macrae
Chronic post-surgical pain: 10 years on
Br. J. Anaesth., July 1, 2008; 101(1): 77 - 86.
[Abstract] [Full Text] [PDF]

Home page
 Palliat MedHome page
F Skorpen, E. Laugsand, P Klepstad, and S Kaasa
Variable response to opioid treatment: any genetic predictors within sight?
Palliative Medicine, June 1, 2008; 22(4): 310 - 327.
[Abstract] [PDF]

Home page
 J. Neurosci.Home page
S. G. Khasar, J. Burkham, O. A. Dina, A. S. Brown, O. Bogen, N. Alessandri-Haber, P. G. Green, D. B. Reichling, and J. D. Levine
Stress Induces a Switch of Intracellular Signaling in Sensory Neurons in a Model of Generalized Pain
J. Neurosci., May 28, 2008; 28(22): 5721 - 5730.
[Abstract] [Full Text] [PDF]

Home page
 Rheumatology (Oxford)Home page
K. L. Limer, B. I. Nicholl, W. Thomson, and J. McBeth
Exploring the genetic susceptibility of chronic widespread pain: the tender points in genetic association studies
Rheumatology, May 1, 2008; 47(5): 572 - 577.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
A. Nikoshkov, K. Drakenberg, X. Wang, M. Cs. Horvath, E. Keller, and Y. L. Hurd
Opioid neuropeptide genotypes in relation to heroin abuse: Dopamine tone contributes to reversed mesolimbic proenkephalin expression
PNAS, January 15, 2008; 105(2): 786 - 791.
[Abstract] [Full Text] [PDF]

Home page
 J. Dent. Res.Home page
G.D. Slade, L. Diatchenko, K. Bhalang, A. Sigurdsson, R.B. Fillingim, I. Belfer, M.B. Max, D. Goldman, and W. Maixner
Influence of Psychological Factors on Risk of Temporomandibular Disorders
J. Dent. Res., November 1, 2007; 86(11): 1120 - 1125.
[Abstract] [Full Text] [PDF]

Home page