Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

Haplotypes of P2RX7 gene polymorphisms are associated with both cold pain sensitivity and analgesic effect of fentanyl

  • 1968 Accesses

  • 11 Citations



The P2X7 receptor is a member of the P2X family of adenosine 5′-triphosphate-gated cation channels. Several recent studies have demonstrated that this receptor is involved in mechanisms related to pain and inflammation. However, unknown is whether polymorphisms of the P2RX7 gene that encodes the human P2X7 receptor influence pain sensitivity and analgesic effects of opioids. The P2RX7 gene is known to be highly polymorphic. Thus, the present study examined associations between fentanyl sensitivity and polymorphisms in the P2RX7 gene in 355 Japanese patients who underwent painful orofacial cosmetic surgery.


We first conducted linkage disequilibrium (LD) analyses for 55 reported single-nucleotide polymorphisms (SNPs) in the region within and around the P2RX7 gene using genomic samples from 100 patients. In our samples, 42 SNPs were polymorphic, and a total of five LD blocks with six Tag SNPs (rs2708092, rs1180012, rs1718125, rs208293, rs1718136, and rs7132846) were observed. Thus, we further analyzed associations between genotypes/haplotypes of these Tag SNPs and clinical data using a total of 355 samples. In the genotype-based association study, only the rs1718125 G > A SNP tended to be associated with higher pain scores on a visual analog scale 24 h after surgery (VAS24). The haplotype-based association study showed that subjects with homozygous haplotype No.3 (GTAAAC; estimated frequency: 15.0%) exhibited significantly higher cold pain sensitivity and lower analgesic effects of fentanyl for acute cold pain in the cold pressor test. Conversely, subjects who carried haplotype No.1 (ACGGAC; estimated frequency: 24.5%) tended to exhibit lower cold pain sensitivity and higher analgesic effects of fentanyl. Furthermore, subjects with homozygous haplotype No.2 (GCGGAC; estimated frequency: 22.9%) exhibited significantly lower VAS24 scores.


Cold pain sensitivity and analgesic effects of fentanyl were related to the SNP and haplotypes of the P2RX7 gene. The patients with the rs1718125 G>A SNP tended to show higher VAS24 scores. Moreover, the combination of polymorphisms from the 5′-flanking region to exon 5 recessively affected cold pain sensitivity and analgesic effects of opioids for acute cold pain. The present findings shed light on the involvement of P2RX7 gene polymorphisms in naive cold pain sensitivity and analgesic effects of fentanyl.


Extracellular adenosine 5′-triphosphate (ATP) has been recognized as a neurotransmitter and/or neuromodulater in the nervous system that specifically acts on P2 purinergic receptors on the cell surface. P2 purinergic receptors are divided into two classes. P2X receptors are ATP-gated cation channels and subdivided into seven subtypes (P2X1–7). P2Y receptors are heptahelical G-protein-coupled receptors and subdivided into eight subtypes (P2Y1, P2Y2, P2Y4, P2Y6, P2Y11, P2Y12, P2Y13, and P2Y14). Both P2X and P2Y receptors are widely expressed in the sensory nerve system and exert various effects on neuronal and glial cells [1]. Recent studies revealed that ATP and its receptors are involved in peripheral and central nociceptive transmission, including mechanisms involved in neuropathic pain [2, 3].

P2X7 receptors exhibit unique pharmacological characteristics compared with other P2X receptor subtypes. A high concentration of ATP (i.e., >100 μM) is required for P2X7 receptor activation [4]. In addition to acting as ATP-gated Ca2+-permeable cation channels, P2X7 receptors induce the formation of large nonselective pores with a 900 Da cut-off [4, 5]. Many studies have shown the involvement of P2X7 receptors in pain. P2X7 receptor knockout mice have been shown to exhibit a reduction of thermal and mechanical hypersensitivity in a partial sciatic nerve ligation model [6]. Recent developments in selective inhibitors of P2X7 receptors also showed that P2X7 receptor blockade reduced nociceptive behavior in several animal models of neuropathic and inflammatory pain [710]. Although these studies revealed an important role for P2X7 receptors in neuropathic and inflammatory pain development in animal models, the involvement of P2X7 receptors in the modulation of naive pain sensitivity and efficacy of analgesics in humans is still unclear.

Opioid analgesics, such as fentanyl and morphine, are widely used for the treatment of moderate to severe pain. However, the analgesic efficacy of opioids is well known to vary widely among individuals [11]. Individual differences may be related to various genetic and nongenetic factors, including gender, age, ethnic origin, hepatic or renal function, and mental status [12]. Several studies that used mice that lack the μ-opioid receptor (MOP) [1315] have shown that analgesia produced by opioids crucially depends on the level of MOP expression. Furthermore, several single-nucleotide polymorphisms (SNPs) in the OPRM1 gene, which encodes the human MOP protein, have been reported to lead to differences in the analgesic efficacy of opioids [16]. Several gene-association studies have also reported that the analgesic efficacy of opioids could be affected by other molecules [1721].

Many gene polymorphisms, most of which are SNPs, reportedly exist in the genes that encode P2X and P2Y receptors. The gene that encodes the human P2X7 receptor (P2RX7) is known to be highly polymorphic. Some SNPs in the P2RX7 gene have been shown to cause changes in receptor function [2224]. Only a few studies have tested associations with human pain sensitivity [24], and whether genetic polymorphisms in the P2RX7 gene exhibit associations with pain sensitivity or opioid analgesia is still unclear. In contrast to animal studies that use standardized pain tests, the analgesic effects of opioids in humans are usually evaluated in patients with actual pain, particularly cancer pain or acute postoperative pain [16]. Patients with acute postoperative pain following standardized surgical procedures may be more optimal subjects for investigating gene-opioid effect relationships [11, 17, 25]. Furthermore, because subjects prior to cosmetic orthognathic surgery have no spontaneous pain, the analgesic effects of opioids in humans can be evaluated under more optimal conditions. Therefore, the present study examined whether SNPs and haplotypes in the P2RX7 gene affect cold pain sensitivity and the analgesic effects of fentanyl, one of the most commonly used opioid analgesics, evaluated by a standardized pain test and fentanyl requirements in healthy Japanese subjects who underwent uniform surgical procedures.

Materials and methods

Ethics statement

The study protocol was approved by the Institutional Review Board, Tokyo Dental College, Chiba, Japan, and the Institutional Review Board, Tokyo Metropolitan Institute of Medical Science, Tokyo, Japan. Written informed consent was obtained from all of the patients and also from parents if required.


Enrolled in the study were 355 healthy patients (American Society of Anesthesiologists Physical Status I, age 15–52 years, 125 males and 230 females [the same patients who served as subjects in our previous report] [17]) who were scheduled to undergo cosmetic orthognathic surgery (mandibular sagittal split ramus osteotomy) for mandibular prognathism at Tokyo Dental College Suidoubashi Hospital. Patients with chronic pain, those taking pain medication, and those who had experienced Raynaud’s phenomenon were excluded.

Preoperative cold pressor-induced pain test

The patients were premedicated with oral diazepam, 5 mg, and oral famotidine, 150 mg, 90 min before the induction of anesthesia. The patients had an intravenous (i.v.) line on the forearm on their nondominant side. The temperature in the operating room was maintained at 26°C. The cold pressor-induced pain test was then performed before and 3 min after an i.v. bolus injection of fentanyl, 2 μg/kg, as previously described [25, 26]. Briefly, crushed ice cubes and cold water were blended 15 min before the test in a 1 L isolated tank, and the mixture was stirred immediately before each test to ensure uniform temperature distribution (0°C) within the tank. The dominant hand was immersed up to the wrist. Patients were instructed to keep the hand calm in the ice-cold water and withdraw it as soon as they perceived any pain. All of the patients were administered the test by the same investigator. The baseline latency to pain perception, defined as the time of immersion of the hand in the ice water, before an i.v. injection of fentanyl (PPLpre) was recorded. A cut-off point of 150 s was set to avoid tissue damage. The hand was warmed with a hair dryer as soon as it was withdrawn from the ice water until the sensation of cold was completely abolished. The patients then received i.v. fentanyl, 2 μg/kg. Three minutes after the injection, the pain perception latency of the dominant hand (PPLpost) was measured again. The analgesic effect of fentanyl in the preoperative cold pressor-induced pain test was evaluated simply as the difference between PPLpost and PPLpre (PPLpost - PPLpre).

Anesthesia, surgery, and postoperative pain management

Anesthesia, surgery, and postoperative pain management were performed as previously described [17]. Briefly, after the cold pressor-induced pain test ended, general anesthesia was induced with a target-controlled infusion (TCI) of propofol. After the induction of anesthesia, 10 ml of venous blood was sampled for the preparation of DNA specimens. Local anesthesia was performed on the right side of the surgical field with 8 ml of 2% lidocaine that contained epinephrine, 12.5 μg/ml, and right mandibular ramus osteotomy was performed. Local anesthesia was then performed on the left side, and left mandibular ramus osteotomy was performed. The bilateral mandibular bone segments were fixed in appropriate positions. Whenever systolic blood pressure or heart rate exceeded +20% of the preinduction value during surgery, i.v. fentanyl, 1 μg/kg, was administered. At the end of the surgery, rectal diclofenac sodium, 50 mg, and i.v. dexamethasone, 8 mg, were administered at the request of surgeons to prevent postoperative orofacial edema/swelling. After emergence from anesthesia and tracheal extubation, i.v. patient-controlled analgesia (PCA) with a fentanyl-droperidol combination (2 mg fentanyl and 5 mg droperidol diluted in normal saline in a total volume of 50 ml) commenced using a CADD-Legacy PCA pump (Smiths Medical Japan, Tokyo, Japan). A bolus dose of fentanyl, 20 μg, on demand and a lockout time of 10 min were set. Patient-controlled analgesia continued for 24 h postoperatively. The intensity of spontaneous pain was assessed 3 and 24 h postoperatively using a 100-mm visual analog scale (VAS), with 0 mm indicating no pain and 100 mm indicating the worst pain imaginable. Intraoperative fentanyl use and postoperative PCA fentanyl use during the first 24 h postoperative period were recorded. Doses of fentanyl administered intraoperatively and postoperatively were normalized to body weight. Additionally, perioperative fentanyl use was calculated as the sum of intraoperative fentanyl use and postoperative fentanyl use.

Genotyping procedures and linkage disequilibrium analysis

Genomic DNA was extracted from whole-blood samples using standard procedures. The extracted DNA was dissolved in TE buffer (10 mM Tris–HCl, 1 mM ethylenediaminetetraacetic acid, pH 8.0). The DNA concentration was adjusted to 5–50 ng/μl for genotyping individual SNPs or 100 ng/μl for whole-genome genotyping using a NanoDrop ND-1000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).

For the analysis of SNPs within and around the P2RX7 gene region, genotype data from whole-genome genotyping were used. Briefly, whole-genome genotyping was performed using Infinium assay II and an iScan system (Illumina, San Diego, CA) according to the manufacturer’s instructions. Five kinds of BeadChips were used to genotype 40, 67, 6, 119, and 123 samples, respectively: HumanHap300 (total markers: 317,503), HumanHap300-Duo (total markers: 318,237), Human610-Quad v1 (total markers: 620,901), Human1M v1.0 (total markers: 1,072,820), and Human 1 M-Duo v3 (total markers: 1,199,187). Some BeadChips included a number of probes specific to copy number variation markers, but most were for SNP markers on the human autosome or sex chromosome. Approximately 300,000 SNP markers were commonly included in all of the BeadChips. After the whole-genome genotyping, the data for genotyped samples were analyzed using BeadStudio or GenomeStudio with the Genotyping module v3.3.7 (Illumina) to evaluate the quality of the results, and the genotype data for all of the SNPs with P2RX7 gene annotation were extracted. In the data-cleaning process, markers that had “Cluster sep” values (i.e., an index of genotype cluster separation) <0.4 and were separated from any of the three genotype clusters were excluded from the subsequent association study.

Single-nucleotide polymorphisms for the association studies were selected based on recently advanced tagging strategies [2729]. To identify relationships between the SNPs used in the study, linkage disequilibrium (LD) analysis was performed in 55 SNPs that were in the approximately 108 kbp region that contained the P2RX7 gene among 1,072,820 markers in the Human 1 M v1.0 BeadChip for 100 samples using Haploview v.4.2 [30]. For the estimation of LD strength between the SNPs, the commonly used D’ and r2 values were pairwise calculated using the genotype dataset of each SNP. Linkage disequilibrium blocks were defined among the SNPs that showed “strong LD,” based on the default algorithm of Gabriel et al. [31], in which the upper and lower 95% confidence limits on D’ for strong LD were set at 0.98 and 0.7, respectively. Tag SNPs in the LD block were consequently determined using Tagger software with default settings, which is incorporated in Haploview and has been detailed in a previous report [29].

Statistical analysis

Parametric and nonparametric data are expressed as mean ± SD and median [interquartile range], respectively. The statistical analysis was performed using IBM SPSS statistics v.20.0.0 (IBM, Tokyo, Japan). In the present study, none of the clinically measured endpoints that were related to pain sensitivity (i.e., PPLpre) or fentanyl analgesia (i.e., analgesia measured with the preoperative cold pressor test, perioperative fentanyl use, and VAS scores at 3 and 24 h postoperatively) were normally distributed. Therefore, nonparametric analyses, including the Mann–Whitney U-test or Kruskal-Wallis test (with Steel-Dwass multiple comparison tests), were used to detect possible associations between any of the clinical or genomic parameters (e.g., sex and genotypes of the Tag SNP) and clinical endpoints related to pain sensitivity or the analgesic effects of fentanyl. The sample size of the present nonparametric data was higher than the estimated size that possesses statistical power (1 minus type II error probability) of 99% for the Cohen’s conventional “medium” effect size of 0.25, when power analysis was performed for analysis of variance with three genotype groups using G*Power v.3.1.3 [32]. Haplotype analyses were performed using HPlus v.3.2 software with default settings (Fred Hutchinson Cancer Research Center, Seattle, WA, USA) that employs expectation-maximization with a modified progressive ligation computational algorithm to infer haplotypes [33].


To identify the LD blocks in the approximately 108 kbp region that contains the P2RX7 gene, 55 SNPs among 1,072,820 markers that were included in the whole-genome genotyping (Human 1 M v1.0 BeadChip) were tested using genomic samples of 100 patients who were randomly selected from all 355 Japanese patients. In these Japanese samples, 42 SNPs were found to be polymorphic. A total of five LD blocks (LD1-5) were observed within and around the P2RX7 gene region (the exon, intron, and approximately 45 kbp 5′-flanking region and 10 kbp 3′-flanking region of the P2RX7 gene [approximately 53 kbp]), and six representative Tag SNPs (rs2708092, rs1180012, rs1718125, rs208293, rs1718136, and rs7132846) were selected in this region (Table 1, Additional file 1). Thus, we further examined associations between clinical endpoints and these six representative Tag SNPs using genomic samples from all 355 Japanese patients.

Table 1 Allelic frequencies of SNPs in five LD blocks around the P2RX7 gene of Japanese subjects

Of the 355 Japanese patients who enrolled in the study, 353 completed the study. The distributions of six representative Tag SNP genotypes are shown in Table 2. No observed genotype frequencies were significantly different from Hardy-Weinberg equilibrium.

Table 2 Distribution of six Tag SNP genotypes examined in the P2RX7 gene

The Mann–Whitney U-test revealed that the VAS score at 24 h was significantly lower (p = 0.019; Table 3) in subjects who did not carry the minor A allele of the rs1718125 SNP than subjects who carried this allele, whereas the rs1718125 SNP had no significant association with PPLpre, PPLpost-PPLpre, 24-h postoperative fentanyl use, perioperative fentanyl use, total perioperative analgesic use, and the VAS score at 3 h. The unpaired t-test, Fisher’s exact test, and Mann–Whitney U-test revealed no significant differences in age, sex, duration of surgery, and duration of anesthesia between subjects who carried the minor A allele and subjects who did not carry this allele (p =0.079 [age: Table 3], p = 0.740 [sex: Table 3], p =0.430 [duration of surgery: AA + AG, 106 min (93, 128); GG, 105 min (92, 122)], and p = 0.384 [duration of anesthesia: AA + AG, 173 min (160, 195); GG, 173 min (157, 193)], respectively). Therefore, we did not conduct additional multivariate covariate analyses. The other five SNPs did not show any associations with any of the clinical endpoints (Table 3). When multiple-testing corrections (i.e., the standard Bonferroni correction for the number of SNPs) were applied, no significant association was found between the genotype of the rs1718125 SNP and VAS score at 24 h.

Table 3 Patients’ demographic and clinical data

We further analyzed the haplotype-based associations of the six P2RX7 gene Tag SNPs with clinical endpoints. Of the 26 estimated haplotypes, 10 haplotypes (estimated frequency >1%) are listed in Table 4. Haplotype No.1 (ACGGAC) was significantly associated with both higher PPLpre score and higher PPLpost-PPLpre score in the linear regression analysis that used the recessive model (Table 5) and log-additive model (PPLpre: coefficient = 4.532 [confidence interval (CI): 0.048, 9.017], z-score = 1.981, p =0.048; PPLpost-PPLpre: coefficient = 8.735 [CI: 1.829, 15.640], z-score = 2.479, p = 0.013) and associated with higher PPLpost-PPLpre score using the dominant model (coefficient =8.802 [CI: 0.601, 17.004], z-score = 2.104, p = 0.035), although the adjusted p values of these associations with haplotype No.1 after multiple-testing corrections (i.e., the standard Bonferroni correction for the number of clinical endpoints) were not significant. The linear regression analyses that used the recessive model revealed that haplotype No.3 (GTAAAC) was significantly associated with both lower PPLpre score and lower PPLpost-PPLpre score, even after multiple-testing corrections (Table 5). These significant associations were not evident when analyzed using either a dominant or log-additive model. Furthermore, linear regression analyses that used the recessive model also revealed that haplotype No.2 (GCGGAC) was significantly associated with lower VAS scores at 24 h, even after multiple-testing corrections (Table 5). This significant association was also not evident when analyzed using either the dominant or log-additive model.

Table 4 Estimated P2RX7 gene Tag SNP haplotypes and their frequencies
Table 5 Associations between P2RX7 gene Tag SNP haplotypes and clinical data


We studied patients who underwent mandibular sagittal split ramus osteotomy. Subjects who undergo this cosmetic surgery are usually young and healthy. The operation causes considerable perioperative pain that arises from the dissected mandibular bone, and the surgical technique is highly standardized at our institute. We conducted a standardized pain test before the induction of general anesthesia in opioid-naive subjects without pain. Using these ideal subjects and methods, we found that PPLpre and the analgesic effects of fentanyl (PPLpost-PPLpre) evaluated in the cold pressor test were significantly lower in subjects who carried homozygous haplotype No.3 (GTAAAC) in the P2RX7 gene compared with the other subjects. Subjects who carried haplotype No.1 (ACGGAC), which has a different allele in the Tag SNP on the LD1-3 regions from haplotype No.3, tended to have higher PPLpre scores and higher analgesic effects of fentanyl. These results suggest that the combination of polymorphisms from the 5′-flanking region to exon 5 (from LD1 to LD3 regions) in the P2RX7 gene could recessively affect cold pain sensitivity and the analgesic effects of opioids for acute cold pain. Furthermore, VAS scores at 24 h tended to be lower in subjects who carried homozygously the major G allele of the rs1718125 SNP (3rd Tag SNP) than subjects who carried the minor A allele. We also found that VAS scores at 24 h were significantly lower in subjects who carried homozygous haplotype No.2 (GCGGAC) in the P2RX7 gene compared with the other subjects. Thus, haplotype No.2 could recessively affect the postoperative analgesic effects of opioids. Although we analyzed five types of BeadChips that were merged together in the present study because the sample sizes were small for each of the five datasets and they all presumably had Japanese ancestry, performing a meta-analysis is usually better than merging. Thus, further studies that have a greater number of samples might be required to reveal the influences of these haplotypic effects in the P2RX7 gene. Twenty-four hour postoperative fentanyl use, perioperative fentanyl use, and total perioperative analgesic use were not associated with haplotypes in the P2RX7 gene. Although further validation is needed, the analgesic/opioid requirements for postoperative pain management might not be associated with genetic polymorphisms in the P2RX7 gene.

In the present study, subjects who carried homozygous haplotype No.3 (GTAAAC) in the P2RX7 gene exhibited higher cold pain sensitivity and lower analgesic effects of fentanyl for acute cold pain in the cold pressor test, but still unknown is whether this combination of SNPs alters gene function or expression, which may be an important limitation of the present study. Many SNPs have been identified in the P2RX7 gene, and some of the SNPs in the P2RX7 gene have been shown to cause changes in receptor function [2224]. Although our present study did not focus on these reported functional SNPs, some of them are located on the present LD blocks. The rs208294 T/C SNP (Tyr155 to His), which was described as the C489T SNP in previous reports with non-Japanese subjects, is located in the LD3 region in the present study. This SNP showed high LD (D’ =1.0) with the present 4th Tag SNP, rs208293 G/A SNP. In the previous report, functional analyses in recombinant cells that expressed P2X7 receptors with the mutation of the rs208294 SNP using [Ca2+]i influx and ethidium uptake experiments revealed that the T > C mutation in this position caused a loss of function of P2X7 receptors [34]. Although the rs208293 G/A SNP alone did not show any association with cold pain sensitivity and the analgesic effects of fentanyl, the rs208294 T/C SNP might be one of the mutations in LD1-3 that caused the present results. Interestingly, P2X7 receptors have been reported to contain several alternative splicing variants (P2X7a-k), some of which presented altered function [3538]. Because the splicing difference in more than half of these isoforms occurs from the LD1 to LD3 regions (from exon 1 to exon 5) in the P2RX7 gene, polymorphisms in these regions might affect this splicing mechanism and induce a functional change in neuronal transmission via P2X7 receptors. Further studies that focus on gene polymorphisms from LD1 to LD3 in the P2RX7 gene may reveal the functional mechanisms that affect pain sensitivity and clinical efficacy of opioids. A recent report showed a genetic association between the hypofunctional rs7958311 G/A SNP (Arg270 to His) and lower pain intensity in two cohort studies of human patients with chronic pain [24]. The rs7958311 G/A SNP is located in the gap between the present LD4 and LD5 regions. Including this SNP, some other nonsynonymous SNPs (e.g., rs3751143 A/C SNP [Glu496 to Ala] [22]) are located in the gap among the present LD regions and were not evaluated in the present study. Further association analyses of these SNPs in the P2RX7 gene may be necessary to reveal the role and mechanisms of P2X7 receptors in naive pain sensitivity and the efficacy of analgesics.

The potential of P2X7 receptors as a therapeutic target in the management of neuropathic pain and inflammation has been the subject of intensive recent investigations [2, 3, 39], but the involvement of P2X7 receptors in the modulation of naive pain sensitivity and efficacy of analgesics is still unclear. P2X7 receptors are expressed predominantly on immune cells, and are speculated to contribute to the hyperexcitability of nociceptive neurons through the release of both interleukin-1β (IL-1β) and tumor necrosis factor α (TNF-α) [4042]. These cytokines are well known to play important roles in the generation and maintenance of pain, and especially in chronic pain states [39, 42]. P2X7 receptors are also known to act not only as cation channels following brief activation by extracellular ATP but also as large nonselective pores following prolonged or repeated activation. Sorge et al. reported that mice that carried a mutation in the P2RX7 gene that causes impaired pore formation showed less hypersensitivity in neuropathic and inflammatory pain states [24]. These previous data might indicate that the pore formation of P2X7 receptors closely affects chronic pain states [43]. In the present study, haplotype No.2 (GCGGAC) and the rs1718125 SNP in the P2RX7 gene showed associations with VAS scores at 24 h but not with VAS scores at 3 h or the analgesic effects of fentanyl in the cold pressor test. Immediate pain transmission in the cold pressor test and the acute phase of postoperative pain (reflected by VAS score at 3 h) might have different pain-modulated mechanisms from subacute postoperative pain (VAS score at 24 h), and the latter might be affected by the pore formation of P2X7 receptors. Further validation of the relationship between these P2RX7 gene polymorphisms and time-dependent changes in postoperative pain states is needed to clarify this issue.

The present study revealed that PPLpre and PPLpost-PPLpre in the cold pressor test were associated with haplotypes of the SNPs (haplotype No.1 and No.3) in the P2RX7 gene, suggesting that P2X7 receptors are involved in cold pain transmission and the analgesic effects of opioids for cold pain stimuli under naive conditions. However, these haplotypes showed no significant association with clinical pain conditions (VAS scores) or fentanyl requirements. The pain stimulus in the cold pressor test may be mediated by the activation of transient receptor potential ankyrin 1 (TRPA1) channels in primary sensory neurons [44]. Postoperative pain may be mediated mainly by the activation of transient receptor potential vanilloid 1 (TRPV1) channels subsequent to incision and inflammation and not by the activation of TRPA1 channels [45]. P2X7 receptors are well known to act as ATP-gated Ca2+-permeable cation channels. Furthermore, the function of TRP channels might be regulated by external Ca2+ [46]. Thus, together with our present results, TRPA1 channels may be affected by the channel function of P2X7 receptors, and TRPV1 channels may be affected by the pore formation of P2X7 receptors. Although further investigation and validation are needed to confirm this hypothesis, P2X7 receptors are suggested to play differential roles in pain transmission, depending on the type of pain stimulus. Furthermore, the mechanisms of P2X7 receptor-induced alterations in the analgesic effects of opioids are still unclear. Only a few reports have shown that the blockade of P2X7 receptors with specific antagonists or targeting small interfering RNA (siRNA) enhanced the analgesic effects of morphine in chronic morphine-treated (i.e., morphine-tolerant) rats, although these effects were not evident in naive rats [47, 48]. Thus, neuronal modulation via P2X7 receptors may affect the analgesic effects of opioids, but these modulatory effects may depend on the type of pain stimulus or pain expression (i.e., acute or chronic).


In Japanese subjects, naive cold pain sensitivity and the analgesic effect of fentanyl were related to an SNP and haplotypes in the P2RX7 gene. The rs1718125 G > A SNP tended to be associated with higher VAS scores at 24 h after mandibular sagittal split ramus osteotomy. Furthermore, subjects with homozygous haplotype No.3 (GTAAAC) exhibited higher cold pain sensitivity and lower analgesic effects of fentanyl for acute cold pain in the cold pressor test. The combination of polymorphisms from the 5′-flanking region to exon 5 in the P2RX7 gene could recessively affect cold pain sensitivity and the analgesic effects of opioids for acute cold pain. Subjects with homozygous haplotype No.2 (GCGGAC) had lower VAS scores at 24 h after surgery. Although further validation is needed, our data may provide valuable information about the role of P2X7 receptors in pain pathways and pain treatment.


  1. 1.

    Burnstock G, Knight GE: Cellular distribution and functions of P2 receptor subtypes in different systems. Int Rev Cytol 2004, 240: 31–304.

  2. 2.

    Burnstock G: Purinergic mechanisms and pain–an update. Eur J Pharmacol 2013, 716: 24–40. 10.1016/j.ejphar.2013.01.078

  3. 3.

    Tsuda M, Tozaki-Saitoh H, Inoue K: Pain and purinergic signaling. Brain Res Rev 2010, 63: 222–232. 10.1016/j.brainresrev.2009.11.003

  4. 4.

    North RA: Molecular physiology of P2X receptors. Physiol Rev 2002, 82: 1013–1067.

  5. 5.

    Locovei S, Scemes E, Qiu F, Spray DC, Dahl G: Pannexin1 is part of the pore forming unit of the P2X(7) receptor death complex. FEBS Lett 2007, 581: 483–488. 10.1016/j.febslet.2006.12.056

  6. 6.

    Chessell IP, Hatcher JP, Bountra C, Michel AD, Hughes JP, Green P, Egerton J, Murfin M, Richardson J, Peck WL, Grahames CB, Casula MA, Yiangou Y, Birch R, Anand P, Buell GN: Disruption of the P2X7 purinoceptor gene abolishes chronic inflammatory and neuropathic pain. Pain 2005, 114: 386–396. 10.1016/j.pain.2005.01.002

  7. 7.

    Honore P, Donnelly-Roberts D, Namovic MT, Hsieh G, Zhu CZ, Mikusa JP, Hernandez G, Zhong C, Gauvin DM, Chandran P, Harris R, Medrano AP, Carroll W, Marsh K, Sullivan JP, Faltynek CR, Jarvis MF: A-740003 [N-(1-{[(cyanoimino) (5-quinolinylamino) methyl]amino}-2,2-dimethylpropyl)-2-(3,4-dimethoxyphenyl)acetamide], a novel and selective P2X7 receptor antagonist, dose-dependently reduces neuropathic pain in the rat. J Pharmacol Exp Ther 2006, 319: 1376–1385. 10.1124/jpet.106.111559

  8. 8.

    McGaraughty S, Chu KL, Namovic MT, Donnelly-Roberts DL, Harris RR, Zhang XF, Shieh CC, Wismer CT, Zhu CZ, Gauvin DM, Fabiyi AC, Honore P, Gregg RJ, Kort ME, Nelson DW, Carroll WA, Marsh K, Faltynek CR, Jarvis MF: P2X7-related modulation of pathological nociception in rats. Neuroscience 2007, 146: 1817–1828. 10.1016/j.neuroscience.2007.03.035

  9. 9.

    Hughes JP, Hatcher JP, Chessell IP: The role of P2X(7) in pain and inflammation. Purinergic Signal 2007, 3: 163–169. 10.1007/s11302-006-9031-1

  10. 10.

    Carroll WA, Donnelly-Roberts D, Jarvis MF: Selective P2X(7) receptor antagonists for chronic inflammation and pain. Purinergic Signal 2009, 5: 63–73. 10.1007/s11302-008-9110-6

  11. 11.

    Ikeda K, Ide S, Han W, Hayashida M, Uhl GR, Sora I: How individual sensitivity to opiates can be predicted by gene analyses. Trends Pharmacol Sci 2005, 26: 311–317. 10.1016/

  12. 12.

    Coulbault L, Beaussier M, Verstuyft C, Weickmans H, Dubert L, Tregouet D, Descot C, Parc Y, Lienhart A, Jaillon P, Becquemont L: Environmental and genetic factors associated with morphine response in the postoperative period. Clin Pharmacol Ther 2006, 79: 316–324. 10.1016/j.clpt.2006.01.007

  13. 13.

    Loh HH, Liu HC, Cavalli A, Yang W, Chen YF, Wei LN: μ opioid receptor knockout in mice: effects on ligand-induced analgesia and morphine lethality. Brain Res Mol Brain Res 1998, 54: 321–326. 10.1016/S0169-328X(97)00353-7

  14. 14.

    Sora I, Elmer G, Funada M, Pieper J, Li XF, Hall FS, Uhl GR: μ opiate receptor gene dose effects on different morphine actions: evidence for differential in vivo μ receptor reserve. Neuropsychopharmacology 2001, 25: 41–54. 10.1016/S0893-133X(00)00252-9

  15. 15.

    Sora I, Takahashi N, Funada M, Ujike H, Revay RS, Donovan DM, Miner LL, Uhl GR: Opiate receptor knockout mice define μ receptor roles in endogenous nociceptive responses and morphine-induced analgesia. Proc Natl Acad Sci U S A 1997, 94: 1544–1549. 10.1073/pnas.94.4.1544

  16. 16.

    Kasai S, Ikeda K: Pharmacogenomics of the human μ-opioid receptor. Pharmacogenomics 2011, 12: 1305–1320. 10.2217/pgs.11.68

  17. 17.

    Ide S, Nishizawa D, Fukuda K, Kasai S, Hasegawa J, Hayashida M, Minami M, Ikeda K: Association between genetic polymorphisms in Ca(v)2.3 (R-type) Ca2+ channels and fentanyl sensitivity in patients undergoing painful cosmetic surgery. PLoS One 2013, 8: e70694. 10.1371/journal.pone.0070694

  18. 18.

    Aoki Y, Nishizawa D, Kasai S, Fukuda K, Ichinohe T, Yamashita S, Ikeda K: Association between the variable number of tandem repeat polymorphism in the third exon of the dopamine D4 receptor gene and sensitivity to analgesics and pain in patients undergoing painful cosmetic surgery. Neurosci Lett 2013, 542: 1–4.

  19. 19.

    Moriyama A, Nishizawa D, Kasai S, Hasegawa J, Fukuda K, Nagashima M, Katoh R, Ikeda K: Association between genetic polymorphisms of the beta1-adrenergic receptor and sensitivity to pain and fentanyl in patients undergoing painful cosmetic surgery. J Pharmacol Sci 2013, 121: 48–57.

  20. 20.

    Nishizawa D, Fukuda K, Kasai S, Hasegawa J, Aoki Y, Nishi A, Saita N, Koukita Y, Nagashima M, Katoh R, Satoh Y, Tagami M, Higuchi S, Ujike H, Ozaki N, Inada T, Iwata N, Sora I, Iyo M, Kondo N, Won MJ, Naruse N, Uehara-Aoyama K, Itokawa M, Koga M, Arinami T, Kaneko Y, Hayashida M, Ikeda K: Genome-wide association study identifies a potent locus associated with human opioid sensitivity. Mol Psychiatry 2014, 19: 55–62. 10.1038/mp.2012.164

  21. 21.

    Zwisler ST, Enggaard TP, Noehr-Jensen L, Mikkelsen S, Verstuyft C, Becquemont L, Sindrup SH, Brosen K: The antinociceptive effect and adverse drug reactions of oxycodone in human experimental pain in relation to genetic variations in the OPRM1 and ABCB1 genes. Fundam Clin Pharmacol 2010, 24: 517–524.

  22. 22.

    Fuller SJ, Stokes L, Skarratt KK, Gu BJ, Wiley JS: Genetics of the P2X7 receptor and human disease. Purinergic Signal 2009, 5: 257–262. 10.1007/s11302-009-9136-4

  23. 23.

    Bradley HJ, Baldwin JM, Goli GR, Johnson B, Zou J, Sivaprasadarao A, Baldwin SA, Jiang LH: Residues 155 and 348 contribute to the determination of P2X7 receptor function via distinct mechanisms revealed by single-nucleotide polymorphisms. J Biol Chem 2011, 286: 8176–8187. 10.1074/jbc.M110.211284

  24. 24.

    Sorge RE, Trang T, Dorfman R, Smith SB, Beggs S, Ritchie J, Austin JS, Zaykin DV, Vander Meulen H, Costigan M, Herbert TA, Yarkoni-Abitbul M, Tichauer D, Livneh J, Gershon E, Zheng M, Tan K, John SL, Slade GD, Jordan J, Woolf CJ, Peltz G, Maixner W, Diatchenko L, Seltzer Z, Salter MW, Mogil JS: Genetically determined P2X7 receptor pore formation regulates variability in chronic pain sensitivity. Nat Med 2012, 18: 595–599. 10.1038/nm.2710

  25. 25.

    Fukuda K, Hayashida M, Ide S, Saita N, Kokita Y, Kasai S, Nishizawa D, Ogai Y, Hasegawa J, Nagashima M, Tagami M, Komatsu H, Sora I, Koga H, Kaneko Y, Ikeda K: Association between OPRM1 gene polymorphisms and fentanyl sensitivity in patients undergoing painful cosmetic surgery. Pain 2009, 147: 194–201. 10.1016/j.pain.2009.09.004

  26. 26.

    Bisgaard T, Klarskov B, Rosenberg J, Kehlet H: Characteristics and prediction of early pain after laparoscopic cholecystectomy. Pain 2001, 90: 261–269. 10.1016/S0304-3959(00)00406-1

  27. 27.

    Carlson CS, Eberle MA, Rieder MJ, Smith JD, Kruglyak L, Nickerson DA: Additional SNPs and linkage-disequilibrium analyses are necessary for whole-genome association studies in humans. Nat Genet 2003, 33: 518–521. 10.1038/ng1128

  28. 28.

    Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 2004, 74: 106–120. 10.1086/381000

  29. 29.

    de Bakker PI, Yelensky R, Pe’er I, Gabriel SB, Daly MJ, Altshuler D: Efficiency and power in genetic association studies. Nat Genet 2005, 37: 1217–1223. 10.1038/ng1669

  30. 30.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21: 263–265. 10.1093/bioinformatics/bth457

  31. 31.

    Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science 2002, 296: 2225–2229. 10.1126/science.1069424

  32. 32.

    Faul F, Erdfelder E, Lang AG, Buchner A: G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods 2007, 39: 175–191. 10.3758/BF03193146

  33. 33.

    Li SS, Khalid N, Carlson C, Zhao LP: Estimating haplotype frequencies and standard errors for multiple single nucleotide polymorphisms. Biostatistics 2003, 4: 513–522. 10.1093/biostatistics/4.4.513

  34. 34.

    Cabrini G, Falzoni S, Forchap SL, Pellegatti P, Balboni A, Agostini P, Cuneo A, Castoldi G, Baricordi OR, Di Virgilio F: A His-155 to Tyr polymorphism confers gain-of-function to the human P2X7 receptor of human leukemic lymphocytes. J Immunol 2005, 175: 82–89. 10.4049/jimmunol.175.1.82

  35. 35.

    Feng YH, Li X, Wang L, Zhou L, Gorodeski GI: A truncated P2X7 receptor variant (P2X7-j) endogenously expressed in cervical cancer cells antagonizes the full-length P2X7 receptor through hetero-oligomerization. J Biol Chem 2006, 281: 17228–17237. 10.1074/jbc.M602999200

  36. 36.

    Cheewatrakoolpong B, Gilchrest H, Anthes JC, Greenfeder S: Identification and characterization of splice variants of the human P2X7 ATP channel. Biochem Biophys Res Commun 2005, 332: 17–27. 10.1016/j.bbrc.2005.04.087

  37. 37.

    Georgiou JG, Skarratt KK, Fuller SJ, Martin CJ, Christopherson RI, Wiley JS, Sluyter R: Human epidermal and monocyte-derived langerhans cells express functional P2X receptors. J Invest Dermatol 2005, 125: 482–490. 10.1111/j.0022-202X.2005.23835.x

  38. 38.

    Xu XJ, Boumechache M, Robinson LE, Marschall V, Gorecki DC, Masin M, Murrell-Lagnado RD: Splice variants of the P2X7 receptor reveal differential agonist dependence and functional coupling with pannexin-1. J Cell Sci 2012, 125: 3776–3789. 10.1242/jcs.099374

  39. 39.

    Alves LA, Bezerra RJ, Faria RX, Ferreira LG, da Silva FV: Physiological roles and potential therapeutic applications of the P2X7 receptor in inflammation and pain. Molecules 2013, 18: 10953–10972. 10.3390/molecules180910953

  40. 40.

    Clark AK, Staniland AA, Marchand F, Kaan TK, McMahon SB, Malcangio M: P2X7-dependent release of interleukin-1β and nociception in the spinal cord following lipopolysaccharide. J Neurosci 2010, 30: 573–582. 10.1523/JNEUROSCI.3295-09.2010

  41. 41.

    Suzuki T, Hide I, Ido K, Kohsaka S, Inoue K, Nakata Y: Production and release of neuroprotective tumor necrosis factor by P2X7 receptor-activated microglia. J Neurosci 2004, 24: 1–7. 10.1523/JNEUROSCI.3792-03.2004

  42. 42.

    Inoue K: P2 receptors and chronic pain. Purinergic Signal 2007, 3: 135–144. 10.1007/s11302-006-9045-8

  43. 43.

    Falk S, Uldall M, Heegaard AM: The role of purinergic receptors in cancer-induced bone pain. J Osteoporos 2012, 2012: 758181.

  44. 44.

    Sawada Y, Hosokawa H, Hori A, Matsumura K, Kobayashi S: Cold sensitivity of recombinant TRPA1 channels. Brain Res 2007, 1160: 39–46.

  45. 45.

    Barabas ME, Stucky CL: TRPV1, but not TRPA1, in primary sensory neurons contributes to cutaneous incision-mediated hypersensitivity. Mol Pain 2013, 9: 9. 10.1186/1744-8069-9-9

  46. 46.

    Banke TG, Chaplan SR, Wickenden AD: Dynamic changes in the TRPA1 selectivity filter lead to progressive but reversible pore dilation. Am J Physiol Cell Physiol 2010, 298: C1457-C1468. 10.1152/ajpcell.00489.2009

  47. 47.

    Zhou D, Chen ML, Zhang YQ, Zhao ZQ: Involvement of spinal microglial P2X7 receptor in generation of tolerance to morphine analgesia in rats. J Neurosci 2010, 30: 8042–8047. 10.1523/JNEUROSCI.5377-09.2010

  48. 48.

    Chen ML, Cao H, Chu YX, Cheng LZ, Liang LL, Zhang YQ, Zhao ZQ: Role of P2X7 receptor-mediated IL-18/IL-18R signaling in morphine tolerance: multiple glial-neuronal dialogues in the rat spinal cord. J Pain 2012, 13: 945–958. 10.1016/j.jpain.2012.06.007

Download references


We acknowledge Mr. Michael Arends for his assistance with editing the manuscript. We are grateful to the volunteers for their participation in this study and the anesthesiologists and surgeons at Tokyo Dental College Suidoubashi Hospital for collecting the clinical data. This work was supported by grants from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan (Tokyo, Japan; no. 22790518, 23390377, 24659549, 24790544, 25116532, 26293347, and 26860360), the Ministry of Health, Labour and Welfare (MHLW) of Japan (Tokyo, Japan; no. H21-3jigan-ippan-011, H22-Iyaku-015, and H25-Iyaku-020), Smoking Research Foundation (Tokyo, Japan), and Astellas Foundation for Research on Metabolic Disorders (Tokyo, Japan).

Author information

Correspondence to Kazutaka Ikeda.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The study was conceived and the experiments were designed by SI, DN, and KI. SI and DN performed the experiments, performed the statistical analyses, and wrote the manuscript. Clinical surgery and the experiments were performed by KF and MH. Genotyping procedures were performed by SK and JH. MM and KI supervised the experiments and finalized the manuscript. All of the authors contributed to writing the manuscript, and all of the authors read and approved the final manuscript.

Soichiro Ide, Daisuke Nishizawa contributed equally to this work.

Electronic supplementary material

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ide, S., Nishizawa, D., Fukuda, K. et al. Haplotypes of P2RX7 gene polymorphisms are associated with both cold pain sensitivity and analgesic effect of fentanyl. Mol Pain 10, 75 (2014).

Download citation


  • P2X7 receptor
  • ATP
  • Purinergic receptor
  • Single-nucleotide polymorphism
  • Pain
  • Fentanyl
  • Cold pain
  • Haplotype analysis
  • Perioperative analgesia