Introduction

Chronic pain is a serious, prevalent, worldwide health problem [1]. Chronic nonspecific back pain is the most common chronic pain syndrome, impacting around 38% of the global population [2]. Opioids remain the frontline treatment for chronic pain conditions, contributing to what the Centers for Disease Control and Prevention have declared to be an opioid epidemic. In concert with a steady rise in prescription opioids, nonmedical use of these medications has also increased, with opioids now reportedly the most commonly abused drugs in the country [3]. In 2019, 9.7 million people (3.5%) aged 12 or older had misused prescription pain relievers in the past year [4]. Orthopedic pain (34.8%) was the primary reason for an opioid prescription, followed by dental conditions (17.3%), back pain (14.0%), and headache (12.9%) [5]. Nearly one-third of chronic pain patients endorse opioid misuse behaviors [6], but not every chronic pain patient prescribed opioids develops a problematic pattern of use. Accordingly, identifying individuals who are at risk for opioid misuse prior to beginning opioid therapy is of significant clinical value. Unfortunately, the processes which underlie enhanced misuse and addiction risk are not currently understood.

While multiple neurotransmitter systems likely contribute to the risk of misuse of pain medications, opioid systems are of particular interest. Activation of the endogenous mu-opioid receptor (MOR) system has been long known to reduce both sensory and affective responses to pain and stress [e.g. [7, 8]], and disruptions in this system are suspected to be involved in the pathogenesis of chronic pain. Human and animal studies indicate substantial interindividual variations in MOR activation and/or endogenous opioid system function. Such differences are further influenced by chronic pain [9,10,11,12,13,14,15,16], introducing varying capacity to respond to exogenous opioid treatment or perhaps even a predisposition to develop persistent pain.

The MOR neurotransmitter system also mediates the reinforcing and hedonic effects of both natural and artificial rewards [17,18,19,20,21]. MORs robustly modulate activity within the mesolimbic pathway—a critical reward circuit consisting of dopamine neuron projections from the ventral tegmental area (VTA) to the nucleus accumbens (NAc) [22]. Activation of MORs located on VTA GABAergic interneurons produces hyperpolarization [23, 24]. This results in disinhibition of DA neurons projecting to the NAc, causing enhanced NAc dopamine release and concomitant increases in reward learning and seeking [see review, [25]]. MOR activation contributes to drug reinforcement and addiction processes, where enhanced MOR activity is associated with heightened drive, increased consumption, and enhanced hedonic reactivity to reward [25].

The present study examines the role of MOR-mediated neurotransmission in risk for prescription opioid misuse in humans. Previous reports have investigated the effects of opioid medications on the MOR system in persons with opioid use disorders (OUD) [26, 27], which likely include adaptations of MOR signaling as a result of chronic opioid use and/or binge/deprivation patterns of use; however, no study to date has examined an at-risk population without current OUD. Here, we utilized [11C]-carfentanil positron emission tomography (PET) with the radiotracer [11C]carfentanil, a selective μ-opioid receptor radioligand [28] to examine relationships between baseline MOR concentrations (non-displaceable binding potential, BPND) and endogenous opioid release in response to a standardized sustained pain challenge, as a function of misuse risk. Misuse risk was calculated from the Pain Medication Questionnaire (PMQ)—a clinical screening instrument used to assess potential for opioid medication misuse in the context of chronic pain [29]—in a sample of chronic nonspecific back pain patients. We initially hypothesized that individuals at high risk for opioid misuse would exhibit higher levels of baseline MOR BPND, possibly secondary to chronic of lower endogenous opioid function in regions (i.e., NAc, amygdala) previously shown to be dysregulated in chronic pain [15] and further associated with addiction [30].

Methods and materials

Participants

We studied 28 chronic pain patients (11 males, 17 females, mean ± SD age, 38.0 ± 10.1) recruited from a pain clinic, and 15 healthy controls (HC; 6 males, 9 females, age 40.5 ± 9.5) recruited via advertisement. [11C]-carfentanil PET data from 16 chronic pain patients have been reported previously [15] and reanalyzed here. All participants were right-handed, non-smoking adults. Patients were excluded if they endorsed lifetime substance dependence, current nicotine dependence, or use of antipsychotics, stimulants, or recreational drugs. Power analyses examining group differences in MOR BPND between persons with fibromyalgia and controls indicated n = 27 would have 90% power to detect an effect at α = 0.05. All participants provided written informed consent. All procedures adopted were approved by the Investigational Review Board and Radioactive Drug Research Committee at the University of Michigan.

Intake measures

Pain-related measures

Chronic pain subjects’ pain experiences were assessed using the McGill Pain Questionnaire (MPQ, [31]), a self-report measure used to quantify the sensory and affective qualities of pain and the Brief Pain Inventory (BPI, [32]), which characterizes pain severity and impact of pain on daily life. As pain fluctuates day to day, patients were asked to complete the BPI each day for a week. These data were averaged and used to compute pain severity and interference. For patients taking prescription opioids (n = 14), morphine milligram equivalent (MME) for total daily opioid dose was calculated using CDC guidelines (Table 1).

Table 1 Chronic back pain characteristics.

Opioid misuse risk

Opioid misuse risk was assessed prior to scanning using the PMQ, a 26-item self-report measure with good internal consistency and test-retest reliability in the context of persistent pain [33]. The PMQ has demonstrated superior predictive utility out of 14 similar assessment tools [34]. Individuals with high, relative to low, PMQ scores exhibit greater risk of OUD and likelihood to prematurely request opioid prescription refills [35]. Patients with PMQ scores >21 were classified as high risk (PMQ-H, n = 13), and those with PMQ scores ≤21 were grouped as low risk (PMQ-L, n = 15) [35].

Personality

Participants completed the Zuckerman-Kuhlman Personality Questionnaire (ZKPQ) prior to scanning to assess dimensions of impulsivity (ZKPQ-Impulsivity), a tendency to act without planning, and sensation seeking (ZKPQ-Sensation Seeking), a preference for novelty and excitement, which have both been previously associated with risk taking [36]. Three PMQ-L participants did not complete the ZKPQ and were excluded from these analyses.

Neuroimaging

PET. [11C]-carfentanil data were acquired with a Siemens HR + scanner (Knoxville, TN). Radiotracer synthesis, image acquisition, and preprocessing protocols have been described in detail elsewhere [7, 37, 38]. Briefly, fifty percent of the [11C] carfentanil dose was administered as a bolus, and the remaining 50% as a continuous infusion using a computer-controlled pump to achieve steady-state tracer levels. The total activity of [11C]-carfentanil administered during each scan was 15.6 ± 0.6 mCi with a mass injection of <0.03 μg/kg per scan (Table 2). PET data were corrected for decay, attenuation, and motion [39]. A modified Logan graphical analysis [40] using the occipital cortex as a reference region was used to transform dynamic image data on a voxel-by-voxel basis into two sets of parametric maps: a tracer transport measure (K1 ratio) and a measure of receptor availability in vivo (non-displaceable binding potential, BPND). The Logan method together with the bolus-continuous infusion allows for linearity in the plot typically 5–7 min after radiotracer administration, allowing for the calculation of BPND values early during scanning. Data obtained from 45 to 90 min post-tracer administration was utilized for the analyses presented here, as previously described [7, 15, 41, 42].

Table 2 Psychophysiological measures.

For spatial normalization, T1-weighted magnetic resonance imaging (MRI) scans were obtained on a General Electric (GE) Sigma 3 T scanner or GE Discovery 3 T scanner (number of slices = 154; voxel resolution = 1 mm3; flip angle = 15°; FOV = 250 × 260mm2; TR = 9.2 ms; TE = 1.9 ms). MRI and PET images were coregistered to each other using SPM12 (www.fil.ion.ucl.ac.uk/spm/) and normalized into Montreal Neurological Institute (MNI) space using the advanced normalization tools (ANTs) software package [43]. PET data were then smoothed with a 6-mm FWHM Gaussian kernel.

Pain-stress challenge

Participants underwent two randomized and counterbalanced PET scans: a pain scan, where participants experienced moderate levels of sustained pain over 20 min, and a baseline scan. This challenge has been described in detail previously [7, 15, 42]. Briefly, pain was induced via infusion of hypertonic saline (5%) into the left masseter muscle 45-min post-radiotracer administration. Every 15 s, participants rated their pain on a VAS scale of 0 (no pain) to 100 (most intense pain imaginable) and the infusion rate was adjusted via a computer-controlled closed-loop system targeting an average VAS rating of about 40 [44, 45]. At the completion of the pain challenge, subjects completed the MPQ to describe their experimentally-induced pain. Prior to scanning and immediately after the pain challenge, participants completed the Positive and Negative Affectivity Schedule (PANAS, [46]), and the Profile of Mood States (POMS, [47]) to changes in affective state and expressed as percent change. One subject in the PMQ-L group did not complete their pain scan and two subjects in both the PMQ-H and control group experienced technical problems and their data could not be utilized, resulting in n = 38 for pain analyses.

Image analysis

We used an a priori region of interest (ROI) approach to explore differences in MOR-related measures between groups. We focused on the NAc and amygdala due to their corresponding functionality in pain, negative affect and, more broadly, their prominent roles in addiction [7, 25, 30, 48, 49]. Additionally, alterations in MOR activity within these regions have been previously described in both animal and human models of addiction [e.g. [50, 51]]. ROIs were defined using the Harvard-Oxford atlas (http://www.cma.mgh.harvard.edu/) using a probability threshold of 25% (see Fig. 1).

Fig. 1: Regions of interest.
figure 1

Location of regions of interest. Visualizations were created with MRIcroGL (http://www.cabiatl.com/mricrogl/).

Two measures of endogenous opioid system function were examined: Baseline MOR BPND and pain-induced changes in MOR BPND, defined as the percent change in MOR BPND from baseline to the pain stress condition. The latter measure, reflects processes associated with endogenous opioid release, such as competition between the radiotracer and endogenous neurotransmitter and reduction in receptor affinity after the activation of MORs by the endogenous ligand [7, 8, 15].

We also conducted whole-brain exploratory analyses to determine whether we could detect group differences in baseline BPND or in pain-induced changes in MOR BPND outside our a priori defined regions of interest. We examined the effects of group on baseline and pain-induced changes to MOR BPND by applying a general linear model on a voxel-by-voxel basis using SPM12 (Wellcome Department of Cognitive Neurology, London, United Kingdom) for Matlab (MathWorks, Natick, Massachusetts). A cluster-forming threshold of p < 0.001 and a peak level significance threshold of p < 0.05 FWE (family-wise error) with a minimum cluster size of 10 voxels was used to define significant peaks.

Statistical analysis

All statistical analyses were conducted using SPSS 26 (Armonk, NY, USA). Correlations between MOR BPND and changes in psychophysical measures were calculated using Kendall’s Tau-b coefficients due to the skewed distribution of these measures [52]. Comparisons between groups were performed using one-way analysis of variance (ANOVA) and t-tests. ROI analyses were corrected for multiple comparisons using a Bonferroni correction, with a significance threshold to p = 0.0125 to account for the 4 ROIs.

Results

There were no significant differences in age (mean ± SD: Control, 40.5 ± 9.5, PMQ-H, 33.9 ± 10.3, PMQ-L 41.5 ± 8.6) between groups (F(2,40) = 2.7, p = 0.08). Baseline clinical pain characteristics including pain severity, pain intensity and pain interference, and morphine equivalents are shown in Table 1. The number of individuals being prescribed opioids (χ2(1,28) = 1.3, p = 0.26) also did not significantly differ between PMQ-H and PMQ-L groups.

A one-way ANOVA revealed significant group differences in ZKPQ-Sensation Seeking (F(2,37) = 4.19, p = 0.02) with Tukey’s post-hoc comparisons indicating significantly higher ZKPQ-Sensation Seeking in PMQ-H compared to PMQ-L (p = 0.03), but only a trending toward significance when compared to the HC group (p = 0.08). However, there were no significant group differences in ZKPQ-Impulsivity observed (F(2,37) = 1.07, p = 0.35).

As it may be expected that opioid misuse risk scores may track ZKPQ-Sensation Seeking and ZKPQ-Impulsivity scores, correlations between ZKPQ-Sensation Seeking, ZKPQ-Impulsivity and PMQ scores were examined. There were no significant relationships detected between PMQ scores and ZKPQ-Impulsivity (τb = 0.16, p = 0.30); however, a significant positive relationship was observed for ZKPQ-Sensation Seeking (τb = 0.40, p = 0.01).

Baseline MOR BPND in chronic pain patients and controls

PMQ

A one-way ANOVA showed significant differences in MOR BPND in the right amygdala (F(2,40) = 5.58, p = 0.01), with Tukey’s post-hoc comparisons indicating the PMQ-H group exhibited higher BPND compared to the PMQ-L group (p = 0.01) and HC group (p = 0.09) (see Fig. 2). A similar trend was noted in the left amygdala (F(2,40) = 2.98, p = 0.06). Consistently, we also noted a significant positive correlation between PMQ scores and MOR BPND in the same region (right amygdala, τb = 0.31, p = 0.02; left amygdala, τb = 0.34, p = 0.01), indicating higher baseline MOR availability was associated with higher risk for opioid misuse. No group differences in NAc MOR BPND were noted (right NAc: F(2,40) = 1.41, p = 0.26; left NAc: F(2,40) = 0.28, p = 0.76).

Fig. 2: Regional differences in Mu-opioid receptor binding between individuals at high and low opioid misuse risk.
figure 2

A Significant differences in baseline MOR BPND between groups was obtained in the right amygdala, with individuals in the high-risk for opioid misuse group (PMQ-H) exhibited the highest mu-opioid receptor MOR BPND. B Significant relationships between opioid misuse risk scores and baseline µ-opioid BP were also observed.

Exploratory whole-brain analyses revealed significant regional effects of group within an area corresponding to the left extended amygdala (x, y, z coordinates, −10, −6, −8; F = 15.72, cluster size = 90 mm3; pFWE = 0.04, Fig. 3) in line with the results of our ROI analyses. No other regional differences as baseline were noted. Post-hoc analyses in SPSS revealed significantly lower MOR BPND in the PMQ-L group relative to both the PMQ-H group (p < 0.001) and the Control group (p = 0.003). No differences were observed between PMQ-H and Control groups (p = 0.163).

Fig. 3: Whole-brain group differences in Mu-opioid receptor binding potential.
figure 3

Significant differences (pFWE < 0.05) in baseline MOR BPND between groups was observed in the left extended amygdala, with individuals in the low-risk for opioid misuse group (PMQ-L) exhibiting lower mu-opioid receptor availability relative to PMQ-H and Control groups.

Affect

Both PMQ-H (p < 0.001) and PMQ-L (p = 0.01) groups showed significantly higher POMS-TMD scores than HC participants (F(2,40) = 9.76, p < 0.001). These differences are unsurprising given the historic association between mood disturbance and chronic pain. No significant group differences in pre-scan positive affect were observed (F(2,40) = 1.45, p = 0.25) though there was a trend (F(2,40) = 3.05, p = 0.06) for higher negative affect among PMQ-H participants relative to controls (Table 2). There were no significant relationships detected between any of the affective measures and binding.

Back pain

Significant negative relationships between MPQ sensory back pain and MOR BPND in the right amygdala (τb = −0.29, p = 0.03) and left amygdala (τb = −0.31, p = 0.02) were observed. No other significant relationships were observed between MOR BPND and back pain intensity, unpleasantness, or among MPQ affective ratings.

Personality

There were no significant relationships between ZKPQ-Sensation Seeking and baseline MOR availability for any group. We did, however, observe significant positive relationships between ZKPQ-Impulsivity and MOR availability within the left NAc (τb = 0.29, p = 0.02), right NAc (τb = 0.31, p = 0.01), left amygdala (τb = 0.27, p = 0.03), and a trend in the right amygdala (τb = 0.22, p = 0.06).

Pain stress-induced changes in MOR BPND

PMQ

A one-way ANOVA indicated group differences in the capacity to activate endogenous opioid neurotransmission in response to the pain stressor (F(2,35) = 4.30, p = 0.02), with Tukey’s post-hoc comparisons indicating MOR release was lower in the PMQ-L group than in the HC group in the left NAc (p = 0.03), however, this was not significant after our stringent Bonferroni correction for multiple comparisons (p < 0.0125). There was a trend towards a positive relationship between PMQ scores and pain-induced changes in MOR BPND within the left NAc (τb = 0.28, p = 0.05). No other group differences or relationships between PMQ scores and pain-induced changes in MOR BPND were noted.

We did not observe any significant effects of group within our whole-brain pain analyses. However, though not significant after correction for multiple corrections, we did observe a peak within an area corresponding to the extended amygdala (x, y, z coordinates, −10, −4, −8; F = 11.39, cluster size = 27mm3; puncorrected < 0.001).

Affect

There were no significant group differences in the pain-induced changes in POMS-TMD (F(2,35) = 0.94, p = 0.40), positive affect (F(2,35) = 0.15, p = 0.86), or negative affect (F(2,35) = 0.83, p = 0.44) scores. We detected negative relationships between ΔPOMS-TMD and pain-induced changes in MOR BPND in the left amygdala (τb = −0.27, p = 0.02), left NAc (τb = −0.24, p = 0.04). A negative relationship between ΔPANAS negative affect and MOR release within the left NAc was also observed (τb = 0.29, p = 0.01). Hence, increases in distress and negative affect during the pain stress challenge were associated with lower endogenous opioid release. No other relationships between change in affect and MOR release were observed.

Back pain

No significant relationships were observed between pain-induced changes in MOR BPND and back pain intensity, unpleasantness, or MPQ ratings.

Experimental pain

We then examined the relationships between pain-induced changes in MOR BPND and experimental pain ratings. One-way ANOVAs indicated no significant group differences in pain intensity (F(2,35) = 2.16, p = 0.13), pain unpleasantness (F(2,35) = 0.57, p = 0.57), average 15-sec VAS ratings (F(2,35) = 2.79, p = 0.08), or MPQ sensory (F(2,35) = 0.39, p = 0.68) or MPQ affective (F(2,35) = 0.11, p = 0.89) ratings following the pain challenge. There were no significant relationships between any of the experimental pain measures and pain-induced changes in MOR BPND.

Personality

Consistent with the baseline findings, there were no significant relationships between ZKPQ-Sensation Seeking or ZKPQ-Impulsivity and pain-stress-induced change in MOR availability for any group.

Discussion

The present study provides the first evidence in humans that risk for opioid misuse is associated with variations in endogenous opioid functional measures. Here we observe that individuals at high risk to misuse opioids (PMQ-H) have higher MOR BPND at rest within regions implicated in emotion and reward processing relative to individuals at lower risk (PMQ-L). Further, that following a standardized pain stressor, lower levels of endogenous opioid release were observed among individuals in the PMQ-L group, relative to controls, within the NAc, however, this difference did not reach our stringent a priori significance threshold. We also show that subjective pain experiences and emotional states are distinctly related to in vivo measures of MOR neurotransmission. Among chronic pain patients, back pain severity ratings were inversely related to baseline MOR BPND, while enhanced negative affect during the pain stressor was related to blunted endogenous opioid system activity. These data may yield insight into the means by which chronic pain and dysregulated pain processing are reciprocally related to one another, and for the first time elucidate the neurobiological validity of a subjective measure of risk for opioid misuse.

The amygdala has long been recognized to play a significant role in pain regulation. MORs in the amygdala have been shown to inhibit nociceptive signaling in preclinical models [53, 54], and activation of endogenous opioid, MOR-mediated neurotransmission in the amygdala is associated with reductions in pain ratings in healthy human subjects [7, 38]. Previous research has posited that regional deficits in MOR-mediated neurotransmission are likely contributors to insufficient endogenous opioid pain control, amplified acute pain sensitivity, and risk for development of additional chronic pain disorders [15, 55, 56]. Indeed, we find both lower MOR availability among the PMQ-L patients relative to controls, as well as significant inverse relationships between self-reported sensory back pain symptoms and amygdala MOR BPND among chronic pain patients. Interestingly, in our ROI analyses, we noted significant differences between groups within the right amygdala but only a trend within the left amygdala. While the origin of such hemispheric differences is not known, there are well described functional asymmetries in the processing of negative emotions within the amygdala, with the right side playing a more prominent role [57, 58]. In addition, hemispheric asymmetries in MOR BPND have been previously noted in large scale PET studies with the amygdala and a NAc displaying higher MOR BPND within the left hemisphere [59].

While lower MOR BPND has been associated with deficient pain control, higher MOR BPND has been associated with risk-taking phenotypes such as high impulsivity and low deliberation [42]. In part, higher MOR BPND among PMQ-H patients may be associated with risk for substance misuse as a function of the opioid system’s influence over preferences for immediate rewards. In humans, greater preference for immediate monetary rewards and greater discounting of probabilistic gains and losses [60] is noted among individuals with opioid addictions relative to non-addicted individuals [61,62,63]—a preference which is exacerbated following mild opioid deprivation [64]. Preclinical models in rodents have demonstrated dose-dependent increases in delay discounting rates [65] and impulsive responding to the five-choice serial reaction time task [66] following morphine administration, while studies in mice have revealed reductions in motor impulsivity [67] and in perseverative responses to obtain reward [68] among MOR knockouts. Here, we observed the highest MOR BPND among those at highest risk to misuse opioids and among those reporting the highest levels of impulsivity but, interestingly, we did not observe any relationships between PMQ scores and impulsivity. This implies these constructs may be capturing distinct information about addiction and misuse vulnerability. But these data, along with previous work [42], demonstrate that even when using unique measures of risk, enhanced risk for substance misuse and addiction is related to endogenous opioid system functioning.

A comparatively lower magnitude pain-induced change in MOR BPND was observed in the left NAc among patients in the low risk group. The endogenous opioid system is heavily implicated in both pain and stress [e.g. [7, 38]]. Several lines of evidence indicate lower MOR reactivity in response to a painful challenge reflects lower capacity to effectively regulate responses to stress and may partially underlie the pathophysiology of chronic pain [15]. Here, we find that endogenous opioid activation of MOR-mediated neurotransmission in response to the pain challenge was relatively, though not significantly, blunted in the PMQ-L group relative to controls; an effect not observed in the PMQ-H group. These effects were noted in the absence of significant group differences in either sensory or affective responses to the challenge, indicating that while pain-induced endogenous opioid system activity varied, the experience of pain was similar across groups. We also noted a trend towards a positive relationship between PMQ scores and change in MOR BPND within the left NAc. This suggests that the capacity to activate the endogenous opioid system in response to pain may be an overlapping but also distinctive biomarker of chronic pain and risk for opioid misuse. Finally, we noted inverse relationships between total mood disturbance, negative affect, and endogenous opioid system activation by the experimental challenge, within the left NAc, indicating greater attenuation of negative affect is associated with greater recruitment of the MOR system during the challenge. This is consistent with previous PET studies that have shown MOR system activation in this region is involved in the modulation of not only pain, but also human affective responses [7, 38].

The present study provides the first evidence in humans that risk for opioid misuse is related to interindividual variations in the function of the endogenous opioid system. While previous PET studies have reported alterations in MOR system markers among persons with addictive disorders, including individuals with nicotine dependence [69], alcohol dependence [70], and gambling disorders [71], the present findings are the first to report altered endogenous opioid activity in persons at risk for opioid misuse. A few limitations should be noted in consideration of these results. First, these data were derived from a sample of non-neuropathic pain patients and the results may not be generalizable to populations with other forms of clinical pain, such as generalized muscular (e.g., fibromyalgia) or neuropathic pain. Second, the sample size is relatively small and cross-sectional, however well-controlled through specialized pain clinics. Future studies with larger, more heterogeneous samples will allow for more nuanced examinations of causal effects and aid in unpacking the unique contributions of psychological and neurobiological variables to risk for opioid misuse in the context of persistent pain. Additionally, our chronic pain sample included individuals at risk for opioid misuse, but did not include persons with a history of opioid misuse. Prevalence of opioid misuse within chronic pain populations is estimated to be as high as 29% [72]. It will be important for future investigations to determine how MOR signaling varies according not only to misuse risk but also according to misuse status. Despite these limitations, the present study contributes novel findings to the existing body of literature on neurobiological features of risk for opioid misuse.