Alterations in the resting-state functional connectivity and hyperperfusion of pain processing areas of the brain have been demonstrated in painful diabetic peripheral neuropathy (DPN). However, the mechanisms underlying these abnormalities are poorly understood; thus there is good rationale to explore whether there is higher energy consumption in the pain processing areas of the brain. We performed a 31P magnetic resonance spectroscopy study to explore cellular energy usage (bioenergetics) in the primary somatosensory (S1) cortex in a well-characterized cohort of participants with painful and painless DPN. S1 phosphocreatine (PCr):ATP, a measure of energy consumption, was significantly reduced in painful compared with painless DPN. This is indicative of greater S1 cortical energy consumption in painful DPN. Furthermore, S1 PCr:ATP correlated with pain intensity during the MRI. S1 PCr:ATP was also significantly lower in painful-DPN individuals with moderate/severe pain compared with those with low pain. To our knowledge, this is the first study to demonstrate higher S1 cortical energy metabolism in painful compared with painless DPN. Moreover, the relationship between PCr:ATP and neuropathic pain measures shows that S1 bioenergetics is related to the severity of neuropathic pain. S1 cortical energetics may represent a biomarker of painful DPN and could have the potential to serve as a target for therapeutic interventions.

Article Highlights

  • Energy consumption within the primary somatosensory cortex appears to be greater in painful compared with painless diabetic peripheral neuropathy.

  • The measure of energy metabolism, PCr:ATP, within the somatosensory cortex correlated with pain intensity and was lower in those with moderate/severe compared with low pain.

  • To our knowledge. this is the first study to indicate higher cortical energy metabolism in painful compared with painless diabetic peripheral neuropathy, and thus has the potential to act as a biomarker for clinical pain trials.

Diabetic peripheral neuropathy (DPN) is the most common chronic complication of diabetes and causes peripheral sensory loss, which is associated with a high risk of foot ulceration and subsequent amputation (1). Up to half of patients with DPN suffer with distressing painful neuropathic symptoms (painful DPN) that result in sleep interference, mood disorders, and a reduced quality of life (2). Unfortunately, current treatments provide suboptimal pain relief and are complicated by adverse effects (3). A greater insight into the pathophysiology of painful DPN is therefore urgently required.

There is emerging evidence for the involvement of the central nervous system, with alterations in both structure and function of the brain in DPN (46). A number of studies have reported abnormalities in patients with painful DPN, including dysfunction of the descending modulatory system (4), lower connectivity of the thalamus and primary somatosensory (S1) cortex (7), and functional reorganization (6) of the pain processing areas in patients with the painful insensate phenotype. Moreover, there is a reduction in thalamic neuronal function in patients with painless DPN, which is preserved in painful DPN and equivalent to levels in healthy volunteers (8). Recent studies also show hyperperfusion of the thalamus (9) and anterior cingulate cortex (emotional pain processing region) (10) in painful DPN, which may be normalized following the treatment of pain (10). These studies are indicative that there may be neuronal hyperactivity in patients with painful DPN in the pain processing areas of the brain. There is therefore a clear rationale to study whether there is higher energy consumption of neurons in the pain processing areas of the brain.

Under normal physiological conditions, adenosine triphosphate (ATP) is generated through mitochondrial phosphorylation whereby adenosine diphosphate (ADP) is combined with inorganic phosphate (Pi). However, during periods of increased energy expenditure or reduced ATP generation (e.g., hypoxia), phosphocreatine (PCr) donates Pi to ADP in order to maintain adequate cellular levels of ATP. PCr is considered an energy reservoir that is important for maintaining a constant level of ATP during physiological or pathological perturbations (11). The 31P magnetic resonance spectroscopy (31P-MRS) is an imaging modality that allows the detection of phosphorus-containing metabolites including ATP, PCr, and Pi, within a particular region of interest known as a voxel (11). The 31P-MRS is a robust technique (11) that provides an objective and reproducible measure of these metabolites and has been used in migraine (12), where higher energy consumption (lower PCr) was shown during a migraine attack. However, this technique has not been used to investigate DPN or neuropathic pain. Therefore, we conducted a 31P-MRS study in well-characterized cohorts of patients with painful and painless DPN. We hypothesized that there will be greater S1 cortical energy usage (i.e., lower PCr:ATP) in painful compared with painless DPN.

Participants

Thirty-two right-handed adults with type 2 diabetes and DPN, aged ≥18 years, were recruited consecutively from Sheffield Teaching Hospitals diabetes outpatient clinics to the study between 2018 and 2020. Exclusion criteria included current or past history of excessive alcohol usage (>14 units per week); nondiabetic neuropathies; other neurological disorders which may confound radiological/clinical assessments, or other factors which would preclude MRI (e.g., pacemaker/claustrophobia); pregnancy; major lower limb amputation; estimated glomerular filtration rate <45 mL/min/1.73 m2; antidepressant medication use; recurrent severe hypoglycemia; and moderate-to-severe pain from causes other than DPN. All participants gave written, informed consent before participation in the study, which had prior ethics approval by the Nottingham Research Ethics Committee, Nottingham (U.K.) (reference no. 17/EM/0430).

Peripheral Neurological and Clinical Assessments

Participants underwent detailed clinical and neurophysiological assessments; full details are available in Supplementary Methods. In brief, this included structured assessment of neuropathy symptoms and signs (Toronto Clinical Neuropathy Score [TCNS]) (13), assessment of neuropathic pain symptoms and signs (Doleur Neuropathique en 4 [DN4]) (14) and different neuropathic pain symptoms (Neuropathic Pain Symptom Inventory [NPSI]) (15), lower limb comprehensive quantitative sensory testing (16), and nerve conduction studies.

On the basis of these assessments, participants were divided into two groups (17): 1) confirmed painless DPN, consisting of participants with painless DPN, with abnormal clinical findings (TCNS >5) and at least two abnormalities on nerve conduction studies one of which had to be an attribute of the sural nerve; and 2) confirmed painful DPN, with abnormal clinical findings (TCNS >5) and at least two abnormalities on nerve conduction studies—one of which had to be an attribute of the sural nerve—DN4 score ≥4, and chronic neuropathic pain for at least 6 months (18).

The Hospital Anxiety and Depression Scale (HADS) was used to determine the presence of depression (HADS-D) and anxiety (HADS-A), using a cut-off score of >10 (19).

31P-MRS

The 31P-MRS was performed with the patients having discontinued pain-relieving medications for over 48 h. Participants had blood glucose measured before the scan to exclude hypoglycemia. An 11-point numeric rating scale (NRS, 0 = no pain to 10 = most severe pain ever) was used to record pain severity during the scan.

Participants underwent cerebral 31P-MRS [echo time 0.26 ms; repetition time 4,000 ms; acquisition technique image-selected in vivo spectroscopy (20)] at 3 tesla (Ingenia 3.0T; Phillips Healthcare, Best, the Netherlands) using a dual-tuned 1H/31P head coil (RAPID Biomedical GmbH Rimpaer, Germany). A two-dimensional survey scan was performed for voxel localization. The voxels (size 25 × 25 × 40 mm3) were placed over the foot region of the S1 cortex (Fig. 1A). The 31P-MRS spectra (Fig. 1B) were quantified using the non-least-squares algorithm function (21) within the MR user interface package (jMRUI, version 5.2) (22).

Figure 1

A: Example of the placement of spectroscopy matrix in the sagittal, coronal, and transverse planes. The white box indicates the placement of the matrix over the bilateral S1 cortices. B: MRS processing in Java-based version of the user interface package (version 5.2) with original spectra after preprocessing (in red) and spectral line fitting for metabolite quantification (in purple). The quantifiable metabolites using 31P-MRS include, from left to right, PME (phosphocholine plus phosphoethanolamine; 6.78 ppm and 6.52 ppm, respectively), Pi (4.85 ppm), PDE (glycerophosphocholine plus glycerophosphoethanolamine; 2.92 ppm and 3.42 ppm, respectively), PCr (0.05 ppm), γ-ATP (−3.19 ppm and −2.62 ppm), α-ATP (−7.75 ppm and −7.4 ppm), and β-ATP (−16.86 ppm, −16.58 ppm, and 16.24 ppm). The area under the peaks is proportionate to the quantity of the metabolite within the voxel. PDE, phosphodiesters; PME, phosphomonoesters; ppm, parts per million.

Figure 1

A: Example of the placement of spectroscopy matrix in the sagittal, coronal, and transverse planes. The white box indicates the placement of the matrix over the bilateral S1 cortices. B: MRS processing in Java-based version of the user interface package (version 5.2) with original spectra after preprocessing (in red) and spectral line fitting for metabolite quantification (in purple). The quantifiable metabolites using 31P-MRS include, from left to right, PME (phosphocholine plus phosphoethanolamine; 6.78 ppm and 6.52 ppm, respectively), Pi (4.85 ppm), PDE (glycerophosphocholine plus glycerophosphoethanolamine; 2.92 ppm and 3.42 ppm, respectively), PCr (0.05 ppm), γ-ATP (−3.19 ppm and −2.62 ppm), α-ATP (−7.75 ppm and −7.4 ppm), and β-ATP (−16.86 ppm, −16.58 ppm, and 16.24 ppm). The area under the peaks is proportionate to the quantity of the metabolite within the voxel. PDE, phosphodiesters; PME, phosphomonoesters; ppm, parts per million.

Close modal

The main study outcomes were ATP, PCr, Pi, and pH. By convention, metabolite levels were expressed as ratios: PCr:ATP as a measure of cellular energetics, with lower values indicative of greater energy usage; and Pi:ATP, which measures the status of oxidative phosphorylation within the cellular mitochondrion, with lower values suggestive of dysfunction. Metabolite ratios and tissue pH were calculated for the left, right, and mean (total) S1 cortical value.

Statistical Analysis

Details are provided in the Supplementary Material and legends.

As 31P-MRS has not previously been performed in chronic pain conditions, we used a migraine study to determine our sample size. A sample size of at least 12 provides 0.90 power (α < 0.05) to detect a 0.12 difference in Pi:ATP (12).

Data and Resource Availability

The data sets generated during the current study are available from the corresponding author upon reasonable request.

There were no significant differences in demographic, biochemical, and neurophysiological variables (Table 1). The TCNS was greater in the painful compared with the painless DPN group (t test, P = 0.017, 95% CI −7.2, −0.76). As expected, the DN4 and total NPSI scores were significantly higher in the painful DPN group. Objective measures of neuropathy, including nerve conduction parameters, and quantitative sensory testing z-scores (Supplementary Table 1) were not significantly different between groups.

Table 1

Demographic details, clinical, biochemical, and neurophysiological assessments of study participants

Painless DPN (n = 12)Painful DPN (n =20)P value
Age, years 64.1 ± 6.3 61.6 ± 8.5 0.39 
Sex, n female (%) 7 (58%) 8 (40%) 0.31* 
Duration of diabetes, years 15.8 ± 5.7 12.5 ± 8.3 0.23 
Pack years smoking 12.3 ± 21.1 15.2 ± 18.6 0.69 
BMI (kg/m232.2 ± 3.4 31.3 ± 5.2 0.60 
HbA1c (mmol/mol) 64.4 ± 17.0 65.0 ± 21.9 0.94 
HbA1c (%) 8.05 ± 1.6 8.10 ± 2.0 0.94 
Creatinine (µmol/L) 77.9 ± 12.6 71.2 ± 16.0 0.22 
Total cholesterol (mmol/L) 3.9 ± 1.0 4.3 ± 1.0 0.27 
Blood glucose during MR (mmol/L) 9.7 ± 4.7 9.7 ± 5.3 0.99 
TCNS 9.9 ± 4.7 13.9 ± 4.1 0.02 
CPN CMAP (mV) 1.9 ± 1.8 2.7 ± 2.2 0.36 
CPN MNCV (m/s) 33.9 ± 13.2 38.2 ± 5.3 0.25 
CPN MNDL (ms) 10.5 ± 10.2 6.2 ± 2.9 0.11 
Tibial MNDL (ms) 6.0 ± 1.8 6.6 ± 2.1 0.43 
Sural SNAP (mV) 1.9 ± 1.7 4.8 ± 6.3 0.26 
NPSI (total score) 0.7 ± 1.7 22.3 ± 9.9 <0.01 
DN4 0.8 ± 1.2 6.9 ± 1.5 <0.01 
Depression, n present (%) 0 (0%) 3 (15%) 0.16* 
Anxiety, n present (%) 1 (8%) 4 (20%) 0.38* 
Analgesic medications, n (%) N/A  N/A 
 α2-δ ligands  8 (40%)  
 Duloxetine  4 (20%)  
 Amitriptyline  6 (30%)  
 Opiates  3 (15%)  
 Topical agents  1 (5%)  
 Other  2 (10%)  
31P-MRS results 
 Right S1 PCr:ATP 2.63 ± 0.60 2.15 ± 0.48 0.02 
 Right S1 Pi:ATP 1.15 ± 0.62 1.02 ± 0.57 0.63 
 Right S1 pH 7.032 ± 0.05 7.043 ± 0.08 0.66 
 Left S1 PCr:ATP 2.43 ± 0.51 1.95 ± 0.38 <0.01 
 Left S1 Pi:ATP 1.06 ± 0.41 1.13 ± 0.50 0.67 
 Left S1 pH 7.036 ± 0.04 7.056 ± 0.11 0.58 
 Total S1 PCr:ATP 2.53 ± 0.52 2.05 ± 0.40 <0.01 
 Total S1 Pi:ATP 1.11 ± 0.39 1.08 ± 0.48 0.19 
 Total S1 pH 7.036 ± 0.04 7.047 ± 0.08 0.57 
Painless DPN (n = 12)Painful DPN (n =20)P value
Age, years 64.1 ± 6.3 61.6 ± 8.5 0.39 
Sex, n female (%) 7 (58%) 8 (40%) 0.31* 
Duration of diabetes, years 15.8 ± 5.7 12.5 ± 8.3 0.23 
Pack years smoking 12.3 ± 21.1 15.2 ± 18.6 0.69 
BMI (kg/m232.2 ± 3.4 31.3 ± 5.2 0.60 
HbA1c (mmol/mol) 64.4 ± 17.0 65.0 ± 21.9 0.94 
HbA1c (%) 8.05 ± 1.6 8.10 ± 2.0 0.94 
Creatinine (µmol/L) 77.9 ± 12.6 71.2 ± 16.0 0.22 
Total cholesterol (mmol/L) 3.9 ± 1.0 4.3 ± 1.0 0.27 
Blood glucose during MR (mmol/L) 9.7 ± 4.7 9.7 ± 5.3 0.99 
TCNS 9.9 ± 4.7 13.9 ± 4.1 0.02 
CPN CMAP (mV) 1.9 ± 1.8 2.7 ± 2.2 0.36 
CPN MNCV (m/s) 33.9 ± 13.2 38.2 ± 5.3 0.25 
CPN MNDL (ms) 10.5 ± 10.2 6.2 ± 2.9 0.11 
Tibial MNDL (ms) 6.0 ± 1.8 6.6 ± 2.1 0.43 
Sural SNAP (mV) 1.9 ± 1.7 4.8 ± 6.3 0.26 
NPSI (total score) 0.7 ± 1.7 22.3 ± 9.9 <0.01 
DN4 0.8 ± 1.2 6.9 ± 1.5 <0.01 
Depression, n present (%) 0 (0%) 3 (15%) 0.16* 
Anxiety, n present (%) 1 (8%) 4 (20%) 0.38* 
Analgesic medications, n (%) N/A  N/A 
 α2-δ ligands  8 (40%)  
 Duloxetine  4 (20%)  
 Amitriptyline  6 (30%)  
 Opiates  3 (15%)  
 Topical agents  1 (5%)  
 Other  2 (10%)  
31P-MRS results 
 Right S1 PCr:ATP 2.63 ± 0.60 2.15 ± 0.48 0.02 
 Right S1 Pi:ATP 1.15 ± 0.62 1.02 ± 0.57 0.63 
 Right S1 pH 7.032 ± 0.05 7.043 ± 0.08 0.66 
 Left S1 PCr:ATP 2.43 ± 0.51 1.95 ± 0.38 <0.01 
 Left S1 Pi:ATP 1.06 ± 0.41 1.13 ± 0.50 0.67 
 Left S1 pH 7.036 ± 0.04 7.056 ± 0.11 0.58 
 Total S1 PCr:ATP 2.53 ± 0.52 2.05 ± 0.40 <0.01 
 Total S1 Pi:ATP 1.11 ± 0.39 1.08 ± 0.48 0.19 
 Total S1 pH 7.036 ± 0.04 7.047 ± 0.08 0.57 

The 31P-MRS results at the primary somatosensory cortex. PCr:ATP is a measure of bioenergetic status, and Pi:ATP is a measure of oxidative phosphorylation. A lower PCr:ATP is indicative of greater energy usage. One pack year smoking is defined as 20 cigarettes × number of years of smoking. The presence of depression and anxiety was defined as a score of >10 on the HADS-D and HADS-A questionnaires, respectively. Normally distributed characteristics are presented as means ± SD. Categorical and dichotomous variables are presented as number of cases and group percentages. All tests were t test unless stated otherwise. Boldface type is used to denote group comparisons with a P value of <0.05. CPN, common peroneal nerve; CMAP, compound muscle action potential; DN4, DN4 questions; MNCV, motor nerve conduction velocity; MNDL, motor nerve distal latency; MR, magnetic resonance scan; SNAP, sensory nerve action potential; N/A, not applicable.

*

χ2 test.

31P-MRS Parameters

Table 1 shows the left, right, and total (left and right combined mean) S1 cortex 31P-MRS results for the two groups. The left, right, and total PCr:ATP were significantly lower in the painful compared with the painless DPN group (left: P = 0.005, 95% CI 0.15, 0.80; right: P = 0.019, 95% CI 0.08, 0.87; total: P = 0.006, 95% CI 0.14, 0.81) (Fig. 2A). There were no significant group differences in Pi:ATP and pH at the right, left, and total S1 cortices.

Figure 2

Box and whisker plots displaying mean PCr:ATP in bilateral S1 cortex (A) and the PCr:ATP in painful-DPN participants (B) divided into low (mild) pain (NRS <3, mean ± SD, 0.4 ± 0.9) and moderate/severe pain (NRS ≥3, mean ± SD, 5.6 ± 1.8). NRS scale from 0 (no pain) to 10 (most severe pain). Statistical tests: in A, t test; in B, Mann-Whitney U test.

Figure 2

Box and whisker plots displaying mean PCr:ATP in bilateral S1 cortex (A) and the PCr:ATP in painful-DPN participants (B) divided into low (mild) pain (NRS <3, mean ± SD, 0.4 ± 0.9) and moderate/severe pain (NRS ≥3, mean ± SD, 5.6 ± 1.8). NRS scale from 0 (no pain) to 10 (most severe pain). Statistical tests: in A, t test; in B, Mann-Whitney U test.

Close modal

Correlation Analysis

The mean S1 cortex PCr:ATP showed significant negative correlations with the DN4 score (Pearson’s r −0.465, P = 0.007), pain intensity (NRS) during the MRI scan (r −0.428, P = 0.016), and total NPSI (r −0.543, P = 0.001). Subscores of the NPSI also correlated with the S1 PCr:ATP (burning spontaneous pain, r −0.424, P = 0.016; pressing spontaneous pain, r −0.375, P = 0.034; paroxysmal pain, r −0.594, P < 0.001; evoked pain, r −0.396, P = 0.025; paresthesia/dysesthesia, r −0.486, P = 0.005). There were no correlations seen between PCr:ATP to other clinical or neurophysiological measures.

Cerebral Bioenergetics (PCr:ATP) and Pain Intensity

There were no significant differences in clinical, demographic, biochemical, or neurophysiological parameters, other than pain severity between high-pain (NRS 5.6 ± 1.8) and low-pain (NRS 0.4 ± 0.9) participants with painful DPN (see Supplementary Table 2). The PCr:ATP was significantly lower in the high-pain group [median (interquartile range) 1.87 (0.35)] compared with the low-pain group [2.35 (0.67), Mann-Whitney U test P = 0.033] (Fig. 2B).

The main finding of this first-ever 31P-MRS study in clinical DPN was that the PCr:ATP, a measure of cellular bioenergetics, was significantly lower in painful compared with painless DPN, indicating higher energy consumption. The PCr:ATP was also significantly lower in participants with moderate/high pain compared with low pain within the painful-DPN group. During periods of high energy demand, PCr acts as a buffer, donating Pi to ADP to maintain cellular ATP levels. Hence the lower levels of PCr are due to increased S1 cortical cellular activity and greater energy usage in painful DPN, specifically in those with high self-reported real-time pain scores. However, the pH and Pi:ATP were similar between groups, indicating that cellular and mitochondrial functions were not different.

The S1 cortex is one of the main components of the cerebral network for processing pain, receiving ascending impulses from important pain pathways. Increased brain activity in this region during experimental pain has been shown to be coupled with increased blood flow (23) and relates to the severity of stimulus intensity (24). These studies support our findings of depleted S1 PCr in painful DPN, likely due to the continuous barrage of afferent impulses. Moreover, our results are consistent with other studies that show evidence of increased cerebral neuronal activity in painful DPN in other brain regions, including increased cerebral blood flow in the thalamus (9) and anterior cingulate cortex (10).

This study also found a relationship between S1 cortical PCr:ATP and measures of neuropathic pain. This further strengthens the link between cerebral energy consumption and neuropathic pain. In particular, the NPSI paroxysmal pain subscore (stabbing and electric shock–type pains) was most significantly correlated with PCr:ATP, whereas there were weaker correlations with other NPSI subscores such as burning and pressing spontaneous pains. A post hoc cluster analysis of the COmbination vs Monotherapy of pregaBalin and dulOxetine in Diabetic Neuropathy (COMBO-DN) painful DPN clinical found that those patients with an overall higher pain intensity were also characterized by higher paradoxical pain scores (25). Therefore, paradoxical pain may be associated with greater S1 cortical energy usage and higher overall pain intensity.

PCr:ATP levels in patients with painless DPN were significantly higher compared with those with painful DPN, despite a similar severity of neuropathy. This signifies the relatively lower resting-state neuronal activity in the S1 cortex compared with patients with painful DPN. Other studies demonstrate that painless DPN may be associated with reduced cerebral neuronal activity, including a reduction in neuronal function and blood flow (9) in the thalamus (8). Moreover, previous studies demonstrate a reduction in S1 cortex peripheral gray matter volume in participants with DPN, but the potential pathophysiological mechanisms have been unexplained (5). To our knowledge, this is the first study to show reduced bioenergetics in patients with painless DPN, suggesting reduced neuronal activity, which may be associated with cortical atrophy or a dying-back axonopathy secondary to reduced sensory input.

The strengths of this study include the use of 31P-MRS, a unique MR modality that does not require the use of an ionizing radiation or administration of contrast/radioactive tracers (11). Pain is an individual and subjective experience, and this technique offers an objective measure of tissue metabolism to gain insight into the cerebral mechanisms of painful DPN. The sample size was adequate for a preliminary, proof-of-concept study and is similar to other cerebral imaging studies in DPN (9). The groups were well characterized using a battery of validated neurological clinical scoring systems and neurophysiological tests. Now, larger cohort studies are necessary to analyze different pain phenotypes (e.g., irritable and nonirritable nociceptors) with the inclusion of deeper assessments of neuronal function such as conditioned pain modulation and spinal inhibitory function (17,26). The limitations of this study include the cross-sectional nature of the study, which does not allow a causal association to be determined.

In conclusion, this preliminary study of 31P-MRS in DPN demonstrates increased cortical energy metabolism in the S1 cortex in participants with painful DPN—specifically in those with greater neuropathic pain severity. Our findings may have clinical relevance, as 31P metabolites have the potential to act as an objective biomarker of painful DPN or neuropathic pain in future clinical trials. Further larger studies are required to confirm these findings, including examination of other brain regions involved with nociceptive processing.

This article contains supplementary material online at https://doi.org/10.2337/figshare.22592887.

Acknowledgments. The authors acknowledge the hard work, skills, and contributions of the radiographers at the University of Sheffield Magnetic Resonance Imaging Department. The authors also greatly appreciate the study volunteers who spent considerable time participating in this study. I.W. is deceased.

Funding. S.T. thanks Sheffield Teaching Hospitals Diabetes Charitable Trust for funding the study.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. G.S. recruited participants, undertook clinical and neurophysiological assessments, researched and analyzed clinical and MR data, and wrote the manuscript. A.A., I.W., S.C., and D.S. made substantial contributions to the design, acquisition, analysis, and interpretation of data, and drafting of the article, and gave final approval of the version published. S.T. made substantial contributions to conception and designs, reviewed and revised the manuscript, and gave final approval of the version published. S.T. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

1.
Sloan
G
,
Selvarajah
D
,
Tesfaye
S
.
Pathogenesis, diagnosis and clinical management of diabetic sensorimotor peripheral neuropathy
.
Nat Rev Endocrinol
2021
;
17
:
400
420
2.
Sadosky
A
,
Mardekian
J
,
Parsons
B
,
Hopps
M
,
Bienen
EJ
,
Markman
J
.
Healthcare utilization and costs in diabetes relative to the clinical spectrum of painful diabetic peripheral neuropathy
.
J Diabetes Complications
2015
;
29
:
212
217
3.
Tesfaye
S
,
Sloan
G
,
Petrie
J
, et al.
OPTION-DM trial group
.
Comparison of amitriptyline supplemented with pregabalin, pregabalin supplemented with amitriptyline, and duloxetine supplemented with pregabalin for the treatment of diabetic peripheral neuropathic pain (OPTION-DM): a multicentre, double-blind, randomised crossover trial
.
Lancet
2022
;
400
:
680
690
4.
Segerdahl
AR
,
Themistocleous
AC
,
Fido
D
,
Bennett
DL
,
Tracey
I
.
A brain-based pain facilitation mechanism contributes to painful diabetic polyneuropathy
.
Brain
2018
;
141
:
357
364
5.
Selvarajah
D
,
Wilkinson
ID
,
Maxwell
M
, et al
.
Magnetic resonance neuroimaging study of brain structural differences in diabetic peripheral neuropathy
.
Diabetes Care
2014
;
37
:
1681
1688
6.
Selvarajah
D
,
Wilkinson
ID
,
Fang
F
, et al
.
Structural and functional abnormalities of the primary somatosensory cortex in diabetic peripheral neuropathy: a multimodal MRI study
.
Diabetes
2019
;
68
:
796
806
7.
Croosu
SS
,
Røikjer
J
,
Mørch
CD
,
Ejskjaer
N
,
Frøkjær
JB
,
Hansen
TM
.
Alterations in functional connectivity of thalamus and primary somatosensory cortex in painful and painless diabetic peripheral neuropathy
.
Diabetes Care
2023
;
46
:
173
182
8.
Gandhi
RA
,
Marques
JL
,
Selvarajah
D
,
Emery
CJ
,
Tesfaye
S
.
Painful diabetic neuropathy is associated with greater autonomic dysfunction than painless diabetic neuropathy
.
Diabetes Care
2010
;
33
:
1585
1590
9.
Selvarajah
D
,
Wilkinson
ID
,
Gandhi
R
,
Griffiths
PD
,
Tesfaye
S
.
Microvascular perfusion abnormalities of the thalamus in painful but not painless diabetic polyneuropathy: a clue to the pathogenesis of pain in type 1 diabetes
.
Diabetes Care
2011
;
34
:
718
720
10.
Watanabe
K
,
Hirano
S
,
Kojima
K
, et al
.
Altered cerebral blood flow in the anterior cingulate cortex is associated with neuropathic pain
.
J Neurol Neurosurg Psychiatry
2018
;
89
:
1082
1087
11.
Liu
Y
,
Gu
Y
,
Yu
X
.
Assessing tissue metabolism by phosphorous-31 magnetic resonance spectroscopy and imaging: a methodology review
.
Quant Imaging Med Surg
2017
;
7
:
707
726
12.
Schulz
UG
,
Blamire
AM
,
Corkill
RG
,
Davies
P
,
Styles
P
,
Rothwell
PM
.
Association between cortical metabolite levels and clinical manifestations of migrainous aura: an MR-spectroscopy study
.
Brain
2007
;
130
:
3102
3110
13.
Bril
V
,
Perkins
BA
.
Validation of the Toronto Clinical Scoring System for diabetic polyneuropathy
.
Diabetes Care
2002
;
25
:
2048
2052
14.
Bouhassira
D
,
Attal
N
,
Alchaar
H
, et al
.
Comparison of pain syndromes associated with nervous or somatic lesions and development of a new neuropathic pain diagnostic questionnaire (DN4)
.
Pain
2005
;
114
:
29
36
15.
Bouhassira
D
,
Attal
N
,
Fermanian
J
, et al
.
Development and validation of the Neuropathic Pain Symptom Inventory
.
Pain
2004
;
108
:
248
257
16.
Rolke
R
,
Baron
R
,
Maier
C
, et al
.
Quantitative sensory testing in the German Research Network on Neuropathic Pain (DFNS): standardized protocol and reference values
.
Pain
2006
;
123
:
231
243
17.
Tesfaye
S
,
Boulton
AJ
,
Dyck
PJ
, et al.
Toronto Diabetic Neuropathy Expert Group
.
Diabetic neuropathies: update on definitions, diagnostic criteria, estimation of severity, and treatments
.
Diabetes Care
2010
;
33
:
2285
2293
18.
Jensen
TS
,
Baron
R
,
Haanpää
M
, et al
.
A new definition of neuropathic pain
.
Pain
2011
;
152
:
2204
2205
19.
Zigmond
AS
,
Snaith
RP
.
The hospital anxiety and depression scale
.
Acta Psychiatr Scand
1983
;
67
:
361
370
20.
Ordidge
RJ
,
Bowley
RM
,
McHale
G
.
A general approach to selection of multiple cubic volume elements using the ISIS technique
.
Magn Reson Med
1988
;
8
:
323
331
21.
Vanhamme
L
,
van den Boogaart
A
,
Van Huffel
S
.
Improved method for accurate and efficient quantification of MRS data with use of prior knowledge
.
J Magn Reson
1997
;
129
:
35
43
22.
D
S
.
Stefan
D
,
Di Cesare
F
,
Andrasescu
A
, et al
.
Quantitation of magnetic resonance spectroscopy signals: the jMRUI software package
.
Meas Sci Technol
2009
;
20
:
104035
23.
Liu
J
,
Hao
Y
,
Du
M
, et al
.
Quantitative cerebral blood flow mapping and functional connectivity of postherpetic neuralgia pain: a perfusion fMRI study
.
Pain
2013
;
154
:
110
118
24.
Bornhövd
K
,
Quante
M
,
Glauche
V
,
Bromm
B
,
Weiller
C
,
Büchel
C
.
Painful stimuli evoke different stimulus-response functions in the amygdala, prefrontal, insula and somatosensory cortex: a single-trial fMRI study
.
Brain
2002
;
125
:
1326
1336
25.
Bouhassira
D
,
Wilhelm
S
,
Schacht
A
, et al
.
Neuropathic pain phenotyping as a predictor of treatment response in painful diabetic neuropathy: data from the randomized, double-blind, COMBO-DN study
.
Pain
2014
;
155
:
2171
2179
26.
Teh
K
,
Wilkinson
ID
,
Heiberg-Gibbons
F
, et al
.
Somatosensory network functional connectivity differentiates clinical pain phenotypes in diabetic neuropathy
.
Diabetologia
2021
;
64
:
1412
1421
Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at https://www.diabetesjournals.org/journals/pages/license.