Skip to main content


We’d like to understand how you use our websites in order to improve them. Register your interest.

Low-frequency BOLD fluctuations demonstrate altered thalamocortical connectivity in diabetic neuropathic pain



In this paper we explored thalamocortical functional connectivity in a group of eight patients suffering from peripheral neuropathic pain (diabetic pain), and compared it with that of a group of healthy subjects. We hypothesized that functional interconnections between the thalamus and cortex can be altered after years of ongoing chronic neuropathic pain.


Functional connectivity was studied through a resting state functional magnetic resonance imaging (fMRI) paradigm: temporal correlations between predefined regions of interest (primary somatosensory cortex, ventral posterior lateral thalamic nucleus, medial dorsal thalamic nucleus) and the rest of the brain were systematically investigated. The patient group showed decreased resting state functional connectivity between the thalamus and the cortex.


This supports the idea that chronic pain can alter thalamocortical connections causing a disruption of thalamic feedback, and the view of chronic pain as a thalamocortical dysrhythmia.


It is already known [1] that certain brain disorders, such as schizophrenia, Alzheimer's disease, autism, epilepsy and Parkinson's are associated with abnormal neural synchronization. This suggests a correlation between abnormalities in neuronal synchronization and cognitive dysfunctions, highlighting the importance of temporal coordination within the brain.

Some recent findings posed the possibility that pain can be related to the emergence of large-scale synchrony, occurring at two different frequency bands; a synchrony that can be enhanced or destroyed by attention and distraction [2]. These results seem to suggest that what matters most, in the case of cortical dynamics during pain, are not the oscillations themselves but their synchrony [3]. Pain seems to imply a "slowing down" of the cortical dynamics [4], and this dynamical aspect of the oscillations is consistent with Llinas's view of chronic pain as a thalamocortical dysrhythmia [5], where the slow rhythm is produced by a disruption of thalamocortical feedback. Indeed, experimental animal models [6] illustrated that, after induced chronic pain, the thalamocortical somatosensory loop shows important synaptic re-organization. A growing number of studies have examined the effects of attention on pain and its neural correlates (for a review see: [710]); in particular, experimental pain can decrease cognitive performance and task-related activity [1113]. Although not completely understood, the cognitive processes underlying pain perception have been shown to be impaired in patients suffering from chronic pain [1319]. Moreover several studies have reported anatomic, functional and biochemical changes in the thalamic region in patients suffering from chronic neuropathic pain suggesting that regions of the thalamus which have lost their normal somatosensory input or that show alterations in the metabolic balance can contain neurons which exhibit abnormal spontaneous and evoked activity and that electrical stimulation of these regions can produce the sensation of burning dysesthesia [7, 2024]. More specifically Iadarola and colleagues [25] using positron emission tomography found a unilateral decrease in thalamic activity in patients with chronic neuropathic pain suggesting that this condition can cause functional alterations in thalamic pain processing circuits.

The recent discovery of low-frequency band fluctuations (LFBF) in the blood oxygenation level dependent (BOLD) signal (for a review see [26]), identified during passive, resting state functional magnetic resonance imaging (fMRI) paradigms [2730], gave researchers a new tool to explore the brain's network dynamics and pathology-related changes. The concept of functional connectivity, introduced by Karl Friston [31] as "The temporal correlations across cortical regions" could represent an index of brain function [31, 32]. Intrinsic functional brain connectivity, revealed by low-frequency spontaneous signal fluctuations in fMRI signal time courses, has received increasing attention in the neuroscience community. Functionally organized systems were first described by Biswal and colleagues [33]; such resting state networks have been subsequently identified in various systems within the cerebral cortex, and related to specific types of sensory, motor and cognitive functions (for a review see [34]), suggesting the existence of several resting state networks in the human brain [27, 35]. Recently, it has been demonstrated that resting state functional connectivity (rsFC) patterns are not artefactually produced by aliasing of cardiac and respiratory cycles, but are instead localized in the gray matter and are likely related to ongoing neuronal activity [36]. Indeed, the BOLD signal within the default mode network (DMN) cannot be explained by cardiac variation rate effects alone and is possibly related to some form of cognitive neuronal processing [27, 3739]. Moreover, it has been demonstrated [27] that these functionally organized systems display BOLD signal changes comparable to task-related BOLD changes (up to 3%) and that these networks are consistent across individuals and stable across repeated sessions.

Unfortunately, despite its central role in the functioning of the cerebral cortex [40, 41], little is known about the participation of the thalamus in the brain's intrinsic activity. Occasionally, thalamocortical patterns of connectivity have been observed in resting state studies [3436, 4246], but few attempts have been made to systematically characterize thalamic nuclei on the basis of spontaneous resting state activity or to explore the effect of pathologies on thalamic rsFC. To the best of our knowledge, only one paper has focused its attention on thalamocortical resting connectivity in normal subjects: Zhang and colleagues studied the intrinsic functional relations between human cerebral cortex and the thalamus [47] finding strong relationships between specific thalamic areas and the cortex. They also defined thalamic boundaries that may serve as a functional atlas for localizing functional areas of the thalamus. In the field of brain pathologies, Welsh et al. [48] used rsFC to show that thalamocortical functional connectivity is altered in schizophrenia; Laureys et al. [49] found that there is a restoration of thalamocortical connectivity after recovery from persistent vegetative state; White et al. [50] evidenced a disruption of functional interactions within neural networks involving the thalamus and cerebral cortex after general-anesthetic-induced unconsciousness; and Mizuno et al. [51] discovered a partially enhanced thalamocortical functional connectivity in autism. On the basis of these results and those of our previous work showing a general disruption of the rsFC in neuropathic pain [52], we hypothesized that thalamocortical reverberation can also be altered after years of ongoing chronic neuropathic pain. Thus, in the present work we investigated the temporal correlations between predefined regions of interest (ROI) in the thalamic nuclei and in the somatosensory cortex (S1) and the rest of the brain, using the method of functional connectivity. This network analysis approach [53] considers the role of a set of brain regions working in unison, assuming that the time courses of interconnected areas show higher correlations [54]. We compared a group of patients suffering from diabetic neuropathic pain with age- and gender-matched healthy controls.


Eight right-handed consecutive patients suffering from peripheral NP (diabetic pain) (four women and four men; age range = 51-78, mean age = 61 years) were enrolled from our multidisciplinary pain unit (tab. s2). All patients underwent a complete neurological and psychological examination as well as standard MR brain scanning by an experienced neuroradiologist (SD) to exclude structural/white-matter abnormalities on anatomical MR-images. Patients suffering from significant psychiatric disorders were excluded. All patients were assessed using standardized pain scales (visual analog scale - VAS, numerical rating scale - NRS, McGill Pain Questionnaire MPQ - Italian version). The spontaneous component of their pain syndrome was evaluated on the MPQ checklist. VAS readings were obtained from their clinical records both the day before and on the day of the study to check the stability of the level of pain (in tab. s2 we reported the value for the day of the study). In all cases pain was restricted to the bilateral lower limbs. Duration of pain was >2 years in all cases. Patients were washed out of their medications at least one month before imaging (opioids or cannabinoids were never administered). At the time of scanning, pain intensity had reached pre-treatment levels. Maximum care was taken to avoid situations that could actually trigger evoked pain during the imaging sessions. Eight age- and gender-matched right-handed healthy volunteers (four women and four men; age range = 47-79, mean age = 59 years) acted as controls. None suffered from any neurological or psychiatric disorders, including chronic pain of any kind, nor had a history of drug or alcohol abuse. None were on medications known to alter brain activity. All females participating in the study were menopausal. All subjects gave their informed written consent, in line with the Declaration of Helsinki, and the study was approved by the local ethics committee.

Task and image acquisition

All subjects were instructed simply to keep their eyes closed, think of nothing in particular, and not to fall asleep. After the scanning session, participants were asked if they had fallen asleep during the scan; subjects who provided a positive or dubious answer were excluded from the study.

Data acquisition was performed on a 1.5 Tesla INTERA™ scanner (Philips Medical Systems) with a SENSE high-field, high resolution (MRIDC) head coil optimized for functional imaging. Resting state functional -weighted images were acquired using echoplanar (EPI) sequences, with a repetition time (TR) of 2000 ms, an echo time (TE) of 50 ms and a 90° flip angle. The acquisition matrix was 64 × 64, the field of view (FoV) 200 mm. A total of 200 volumes were acquired; each volume consisted of 19 axial slices, parallel to the anterior-posterior (AC-PC) commissure; slice thickness was 4.5 mm with a 0.5 mm gap. Two scans were added at the beginning of functional scanning and the data discarded to reach a steady-state magnetization before acquiring the experimental data. In the same session, a set of three-dimensional high-resolution T1-weighted structural images was acquired for each participant. This data-set was acquired using a Fast Field Echo (FFE) sequence, with a TR of 25 ms, ultra-short TE and a 30° flip angle. The acquisition matrix was 256 × 256, the field of view (FoV) 256 mm. The set consisted of 160 contiguous sagittal images covering the whole brain. In-plane resolution was 1 × 1 mm2 and slice thickness 1 mm (voxel of 1 × 1 × 1 mm3).

Data analysis

BOLD imaging data were analyzed using the BrainVoyager QX software (Brain Innovation, Maastricht, Holland). Functional images were preprocessed to reduce artifacts [55] following pre-processing steps: 1) slice scan time correction, using a sinc interpolation algorithm; 2) 3D motion correction: all volumes were aligned spatially to the first volume by rigid body transformations, using a trilinear interpolation algorithm; 3) spatial smoothing using a Gaussian kernel of 6 mm FWHM for the primary sensory (S1) connectivity and 8 mm FWHM for the ventral posterior thalamic nucleus (VPL) connectivity. We used a minor spatial smoothing in the S1 connectivity data to improve the spatial resolution in the thalamic area; 4) temporal filters (i.e. linear trend removal and band pass filter of 0.01-0.1) to reduce low frequency fluctuation in the blood oxygen level dependent (BOLD) signal before the functional connectivity analysis [28, 33].

After pre-processing, a series of steps were followed in order to allow for the precise anatomical location of brain activity to facilitate inter-subject averaging. First, each subject's slice-based functional scan was co-registered on his/her 3D high-resolution structural scan. Second, the 3D structural data-set of each subject was transformed into Talairach space [56]: the cerebrum was translated and rotated into the anterior-posterior commissure plane and then the borders of the cerebrum were identified. Third, using the anatomo-functional coregistration matrix and the determined Talairach reference points, the functional time course of each subject was transformed into Talairach space and the volume time course created.

ROI selection

To obtain the S1 ROI we extracted the location of the peak voxel from a localizer study conducted on eight healthy subjects (random effect analysis, thresholded at q<0.05 corrected with false discovery rate [FDR] [57]) and compared it with anatomic atlases and anatomic landmarks; over the bilateral peak voxel locations we designed a cubic ROI of 5 × 5 × 5 mm3 (coordinates: -43, -26, 49; -46, -22, 49). The localizer study (following Moore et al. [58]) was a blocked design with interleaved epochs of tactile stimulation (brushing) and rest. For determination of thalamocortical connectivity, a seed region in the ventral posterior lateral nucleus (VPL) was determined using an anatomic atlas: we calculated the center of gravity of all voxels comprised in the VPL as mapped in the Talairach Daemon database [59]. Left and right 5 × 5 × 5 mm3 seed ROIs were drawn over these two points (coordinates: 20, -20, 5; -18, -20, 5). In addition, since in the S1 rsFC a relevant area within the medial dorsal thalamic nucleus (MDN) was found to be less connected, we also computed a seed region in the bilateral MDN as determined using an anatomic atlas: we calculated the center of gravity of all voxels comprised in the MDN as mapped in the Talairach Daemon database [59]. Left and right 5 × 5 × 5 mm3 seed ROIs were drawn over these two points (coordinates: +6, -22, 8; -6, -22, 8). All ROIs are listed in tab. s3.

Additional analyses were conducted using alternative ROIs, i.e. ROIs which were moved from the original location in the rostral, lateral, ventral, medial and caudal directions (3 mm), reduced (3 × 3 × 3 voxels) and increased (8 × 8 × 8 voxels) in dimensions. We found high similarity between the resulting maps and those obtained using the original ROIs: probabilistic maps showed high overlapping between the original and the alternative connectivity maps (see Additional file 1: figure S19). A probabilistic map describes the relative frequency (expressed in percent in our cases) in a sample of different maps with which a specific voxel or area is significantly activated.

Functional connectivity analysis

The first step in all functional connectivity (FC) analyses was to extract BOLD time courses from each ROI by averaging over voxels within each region. Several nuisance covariates were included in the analyses to control for the effects of physiological processes (such as fluctuations related to cardiac and respiratory cycles) [60, 61] and motion. Specifically, we included nine additional covariates that modeled nuisance signals from white matter (WM) and cerebrospinal fluid (CSF), global signal (GS) [62], as well as six motion parameters (three rotations and three translations). Additional file 1: Figure s21 shows the overlapping (probabilistic map) of voxels correlated with GS. All seed-based predictors were z-normalized. Temporal autocorrelation correction [63] was used to eliminate the autocorrelation induced by the temporal lowpass filter. Seed ROI-driven FC maps were computed on a voxelwise basis for each previously selected region. Individual participant multiple regression analysis was carried out using the general linear model (GLM) [64] and resulted in a t-based map (SPMt) corrected for multiple comparisons with FDR [57] (q<0.05, cluster threshold >5 voxels in the native resolution).

Group statistical map

Random effect group-level analyses and group comparisons controlling for age and gender effects were conducted using the ANCOVA analysis tool implemented in BrainVoyager QX. Corrections for multiple comparisons were carried out using the FDR [57] (q<0.05, cluster threshold >5 voxels in the native resolution).

Given the variable age in our group, we conducted random effects analyses of covariance (ANCOVAs) to assess the potential effect of this variable on the spatially specific activations revealed in the preceding analyses. To test for significant group-related differences, direct voxelwise group comparisons were performed using group-level contrasts. These contrasts computed the voxelwise statistical significance of mean group differences in FC, and produced FDR thresholded maps of those voxels that showed significant differences between the two groups. For MDN and VPL connectivity the resulting maps were then projected on a 3D representation of the brain using the BrainVoyager QX 3D cortical tool.

To evaluate the spatial consistency of functional connectivity patterns across subjects we computed probabilistic maps, that were calculated separately for each ROI-generated network.

Pairwise correlations between the four a priori ROIs

For each subject and each ROI, we extracted the time course at each voxel. Using a home-made SPSS script, we computed the mean time course for each ROI. For each subject we performed pairwise correlations between the ROIs' mean time courses, in order to assess rsFC. In particular, within each hemisphere, we assessed the rsFC between the thalamus (VPL) and S1; we also assessed the rsFC between each area and its homologous area in the opposite hemisphere (left and right thalamus; left and right S1). Parametric analyses (t-tests for independent samples) were used at the group level to determine the statistical significance of the difference between patients and controls on each paired correlation.


In connectivity comparison between the pain group and the control group we calculated the lateralization (percent of left or right voxels over the total) of the portion of functionally connected Brodmann areas (see additional methods).

Voxel distance calculations

We calculated the Euclidean distance between the center of gravity of the three pairs of ROI (MDN left and right, VPL left and right, S1 left and right) and every other voxel that reached significance in the thresholded Z-score map of positive FC (cluster significance: p < 0.05). The Euclidian distance between two voxels P = (p x , p y , p z ) and Q = (q x , q y , q z ) was computed using the formula . We computed the number of significant voxels at specific distances (from 0 to 140 mm in 4-mm bins) for each individual, then we rebinned these values in three ranges thus comprised: short- (<40 mm), medium- (40-80 mm) and long- (>80 mm) distance connections. Finally we calculated group means and compared them using a simple t-test.


The aim of the study was to compare the resting state functional connectivity (rsFC) in the thalamus and in the primary somatosensory cortex between normal and pain subjects.

Although we used age as a covariate we further inspected possible age effects on the rsFC maps performing a correlational analysis between each subject-specific ROI-generated map and the relative age using the ANCOVA analysis tool implemented in BrainVoyager QX. We found no significant correlation between them (q<0.05 FDR corrected). Movement was assessed by summing the deviations (three translations plus three rotations at a radius of 50 mm) used to compensate for head motion within fMRI. The quantity reported here is head movement root mean square averaged over subjects. This quantity was 0.30 ± 0.09 mm (mean ± standard deviation) for patients and 0.28 ± 0.10 mm (mean ± standard deviation) for age-matched controls. Patients vs. controls comparisons were not significantly different (p = 0.84). Hence, it is unlikely that movement can account for the observed differences in connectivity between the two samples.

Functional connectivity of S1

Comparing the pain and healthy groups (see fig. 1 and also Additional file 1: s10, s22), patients showed a decrease in rsFC between S1 and the left thalamus, including several different nuclei, such as the MDN, laterodorsal (LDN), and pulvinar. The MDN was the most represented nucleus, with 76 voxels found to be less connected with respect to the control group (see Additional file 1: Figure s9). As regards the connectivity between S1 and the rest of the cortex, in the patient group we found several areas of reduced connectivity (Additional file 1: Figure. S12; Table. s4, s6) localized in the bilateral middle temporal gyrus, superior temporal gyrus, pre/post central gyrus, paracentral lobule, middle frontal gyrus, middle occipital gyrus, and cingulated cortex; right putamen, cuneus, fusiform gyrus, lingual gyrus and declive; left superior frontal gyrus and inferior occipital gyrus. The wider areas were in the postcentral gyrus (BA 2, 3; Left Lateralization L = 64%), middle temporal gyrus (BA 22, 21; L = 54%) and superior temporal gyrus (BA 22, 21; L = 72%).

Figure 1

S1 resting state connectivity analysis. A: Probabilistic maps of the S1 rsFC thalamic projection relative to pain subjects. B: Probabilistic maps of the S1 rsFC thalamic projection relative to healthy subjects. C: Group comparison of the S1 rsFC thalamic projection (Two sample t-test q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution). D: Example of the time course of a healthy subject. E: Group comparison of the S1 rsFC brain projection (Two sample t-test q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution). Right middle panel represents the Seed ROI. Colors from green to white indicate an increasing spatial overlapping probability (%) (single subject correlation maps before probabilistic map creation thresholded at q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution). Colors from red to yellow indicate an increased connectivity in the pain group; colors from blue to green indicate a reduced connectivity in the pain group. All rsFC maps were projected on 3D representations of the brain using BrainVoyager QX.

An increase of rsFC in the patient group (Additional file 1: fig. s13; tab. s5, s6) was found bilaterally in the inferior and superior parietal lobule, precuneus, postcentral gyrus (this increased connectivity area is more posterior and less extended compared to the decreased connectivity area in the postcentral gyrus described above) and in the right precentral gyrus. The wider areas were in the inferior parietal lobule (BA 40, L = 69%) and the precuneus (BA 7, Right Lateralization R = 70%).

To assess the spatial consistency of seed-generated S1 correlation patterns across subjects we computed spatial probabilistic maps (see fig. 1). Overall, both groups showed a high overlapping of subjects' maps; more specifically, 62% of the healthy subjects showed S1 connections in the left thalamus and 50% in the right thalamus, while only 12% of the patients showed areas of connectivity in the left thalamus and 0% in the right thalamus.

The healthy group showed no significant differences, compared to the pain group, in mean, medium and long-distance connectivity, but only in short-distance connectivity (Additional file 1: Figure. s18). In all cases patients had fewer connected voxels.

Functional connectivity of VPL

Patients showed a generalized decrease in rsFC (see fig. 2, Additional file 1: s11, s14, s23; tab. s7, s9) between VPL and the bilateral pre/postcentral gyrus, inferior parietal lobule, superior parietal lobule, supramarginal gyrus, precuneus, insula, middle frontal gyrus, inferior frontal gyrus, superior frontal gyrus, medial frontal gyrus and cingulated cortex; right caudate and middle temporal gyrus; left putamen, pulvinar and MDN. The wider areas were in the middle frontal gyrus (BA 10, 9, 46, 6; L = 45%), inferior frontal gyrus (BA 45, 46, 47, 13; L = 50%) and superior frontal gyrus (BA 10, 9, 6, 8; L = 68%).

Figure 2

VPL resting state connectivity analysis. Group comparison of the VPL rsFC brain projection (Two sample t-test q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution). Middle panel represents the Seed ROI. Colors from red to yellow indicate an increased connectivity in the pain group; colors from blue to green indicate a reduced connectivity in the pain group. All rsFC maps were projected on 3D representations of the brain using BrainVoyager QX.

An increase of rsFC in the patient group (see fig. 2, Additional file 1: s15; tab. s8, s9) was found in the bilateral cuneus, bilateral posterior cingulated cortex and bilateral lingual gyrus. The wider area was in the lingual gyrus (BA 18, 19; L = 91%).

To assess the spatial consistency of seed-generated VPL correlation patterns across subjects we computed spatial probabilistic maps. Overall, both groups showed a high overlapping of subjects' maps; more specifically, 75% of the healthy subjects showed VPL connections in the bilateral pre/postcentral gyrus and 66% in the cingulated cortex, while only 62% of the patients showed areas of connectivity in the bilateral pre/postcentral gyrus and 37% in the cingulated cortex.

The healthy group showed no significant differences, compared to the pain group, in mean, medium and long-distance connectivity but only in short-distance connectivity (Additional file 1: fig. s18) and in all cases patients had fewer connected voxels.

Functional connectivity of MDN

Patients showed a generalized decrease of rsFC (see fig. 3, Additional file 1: s16, s23; tab s10, s12) in the bilateral frontal gyrus (inferior, middle, superior), pre/postcentral gyrus, paracentral lobule, precuneus, cingulated cortex, insula, superior temporal gyrus, inferior parietal lobule, pulvinar, parahippocampal gyrus, amygdale, cerebellar tonsil; right pyramis, uvula, nodule; left supramarginal gyrus, VPL, VLN, globus pallidus, caudate. The wider areas were in the left superior temporal gyrus (BA 41, 13, 22, 42; L = 100%) and in the cuneus (BA 30, 18, 23; L = 65%).

Figure 3

MDN resting state connectivity analysis. A: Probabilistic maps of the MDN rsFC Amygdale projection relative to healthy subjects. B: Probabilistic maps of the MDN rsFC Amygdale projection relative to pain subjects. C: Group comparison of the MDN rsFC Amygdale projection. D: Group comparison of the MDN rsFC Subgenual ACC projection. E: Group comparison of the MDN rsFC S1 projection. Colors from green to white indicate an increasing spatial overlapping probability (%) (single subject correlation maps before probability map creation thresholded at q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution). Group comparison with two sample t-test, q<0.05 FDR-corrected, minimum cluster dimension K>5 voxels in the native resolution; colors from red to yellow indicate an increased connectivity in the pain group; colors from blue to green indicate a reduced connectivity in the pain group. All rsFC maps were projected on 3D representation of the brain using BrainVoyager QX.

An increase of rsFC in the patient group (see fig. 3, Additional file 1: s17; tab s11, s12) was found in the bilateral cuneus, posterior cingulated, lingual gyrus; left insula and transverse temporal gyrus. The wider areas were in the inferior frontal gyrus (BA 47, 45, 46, 44; L = 75%), middle frontal gyrus (BA 10, 6, 11, 47; L = 79%) and in the precuneus (BA 7, 31; L = 68%).

To assess the spatial consistency of seed-generated MDN correlation patterns across subjects we computed spatial probabilistic maps. Overall, as with S1 and VPL connectivity, both groups showed a high overlapping of the subject specific maps (see fig. 3), 75% of the healthy subjects showed an area connected to MDN in the left amygdala and 74% in the right amygdala while just 20% and 22% of patients showed areas of connectivity (left and right amygdala respectively).

The healthy group showed no significant difference, compared to the pain group, in mean, short, medium and long-distance connectivity (Additional file 1: fig. s18). Patients always had fewer connected voxels but they were not statistically significant.

Pairwise correlations

Results of pairwise correlations are illustrated in fig. 4. Significant differences emerged between patients and controls in rsFC between left and right S1 (t = 2.670; p = 0.09), and, within the left hemisphere, between S1 and VPL (t = 2.039; p = 0.03). No significant differences emerged between the two thalami (VPL) and between right S1 and right VPL.

Figure 4

Pairwise correlations. Mean correlations and standard deviations for each pair of investigated ROIs in the patient and control groups. Red lines indicate statistically relevant differences in rsFC between the two groups, blue dotted lines indicate no statistically relevant differences.


The methodology we used is particularly suited to analyze brain networks as a whole so we will focus the discussion on the general increase or reduction in thalamic network connectivity rather than on changes in specific cortical areas. The main finding of the present study is decreased thalamocortical connectivity in patients with chronic painful diabetic polyneuropathy. More specifically we found a reduced synchrony between S1 and some of the thalamic nuclei involved in pain processing. Exploring specifically the VPL and MDN connectivity we demonstrated that VPL has a reduced synchrony with sensorymotor, frontal, cingulated and parietal areas while MDN more with sensory and emotional processing areas.

Chronic neuropathic pain, as the result of a predominantly small fiber sensory polyneuropathy, is a frequent complication of diabetes, and patients with impaired glucose tolerance may also suffer from pain [65]. Thus, pain may be a complication of diabetic neuropathy even in the early stages. Differently from previous studies on chronic back pain [66], we dealt with the effect of both deafferentation and neuropathic pain on cerebral functions. There is growing evidence, stemming from both neurophysiological and functional imaging studies, that lesions of the peripheral nerve may sustain plastic changes of the central nervous system, at both cortical and subcortical level. For instance, tremor and dystonia were sometimes observed after a peripheral nerve trauma as well as complex regional pain syndrome [67]. Ischemic nerve block is associated with decreased intracortical inhibition evaluated by transcranial magnetic stimulation [68]. Indeed, modification of cerebral and spinal activity may also be related to pain [69]. For instance, withdrawal reflexes as a result of modification of spinal interneuronal activity may be modulated by heterotopic painful stimuli [70], and painful electrical stimulation of the median nerve may modulate cortical excitability [71]. However, long-term modifications induced by chronic pain on brain function are still a matter of debate [72], and this uncertainty may depend on the different types of neuropsychological and functional imaging protocols employed so far. Pain is characterized by both a spontaneous and an induced dimension (i.e. allodynia). While functional imaging studies based on the analysis of brain activation after painful stimuli seem to be more suitable for evaluating the induced pain mechanism, analysis of resting state or default mode of brain function is an appealing method of evaluating modification of brain activity induced by ongoing spontaneous pain [66]. In everyday neurology patients with chronic neuropathic pain are commonly found to experience the spontaneous painful sensation more when awake and at rest than during daily activities. In a previous study on the same cohort of patients [52] we observed that during rest there is a reduced default mode network [73] connectivity of bilateral primary sensorimotor areas and cingulated cortex, as well as left temporal and occipital cortices; on the other hand, there was increased connectivity in the left dorsolateral prefrontal cortex and posterior parietal cortex, (corresponding to area 7) lateral thalamus and bilateral insula. These results point toward an opposite modification of connectivity in the cingulated cortex and the insulae which are well-recognized parts of the so-called pain matrix, a network of brain areas functionally linked and activated during processing of painful thermal stimuli in healthy subjects [2]. The increased connectivity in the left dorsolateral prefrontal cortex may be coupled with the decreased connectivity in the anterior cingulated cortex, since the former area exerts an inhibitory control on the latter as well as on the medial thalamus [74]. Thus, the reduced connectivity of the medial thalamus observed in the present study is also in keeping with the findings of our previous work. The inhibitory control of the dorsolateral prefrontal cortex on the medial-thalamus-anterior cingulated cortex loop accounts for the reduction of perceived unpleasantness [74] associated with allodynia in patients with sensitized skin. Thus, in patients who suffer frequently from allodynia, such as those affected by diabetic polyneuropathy, a modulation of the left dorsolateral prefrontal cortex-anterior cingulated cortex-medial thalamus loop could take place constitutively. This is further suggested by the findings that show a decreased regional blood flow at rest in the anterior cingulated cortex in patients with chronic pain.

Another possible interpretation is that the reduced connectivity of the anterior cingulated cortex in our patients as well as others with chronic pain is related to mood modulation. This is suggested by the results of stimulation of the motor cortex in patients with chronic intractable pain [75]: PET studies on patients with decreased pain showed an increase in cerebral blood flow in the anterior cingulate cortex and orbitofrontal cortex. The anterior cingulated cortex and medial thalamus are parts of both the system analyzing the affective dimension of pain as well as the reward system encompassing the ventral striatum. The anterior cingulate cortex has connection with both the medial prefrontal cortex and the amygdala. Indeed the present study also found the amygdala to have reduced connectivity to the medial thalamus. The amygdala is implicated in the modulation of perception of emotional stimuli, particularly the enhanced perception of salient stimuli and memory of emotional experience. It also has a role in behavior control related to emotional stimuli, indeed patients with lesions of this structure have blunting of emotional and fear perceptions [76].

Taken together these results may suggest that in our patients there is an increased connectivity of cerebral areas related to inhibitory control of pain (left dorsolateral prefrontal cortex), and a reduced connectivity of areas related to modulation of emotional response to stimuli (amygdale, anterior cingulated cortex, caudate). Although this study did not include any formal evaluation of mood, patients with diabetic painful polyneuropathy frequently suffer from depression as well as sleep disturbances [77], which interfere with quality of life.

Another interesting observation on connectivity of cortical areas pertains to the increased connectivity of the primary sensory cortex with the posterior as well as the inferior parietal cortex. The latter result may be interpreted on the basis of the known role of the posterior parietal cortex in sequential motor performance, the so-called "internal controller" function matching efferent copy of motor command with sensory feedback [78]. In the presence of an altered sensory feedback, an increased connectivity between the primary sensory cortex and the "internal controller" may represent an adaptive plastic change to overcome difficulties in complex movements. These observations may also be related to the known role of the precuneus in episodic memory as well as self-awareness[79]. It is notable that in our patients precuneus connectivity with S1 was increased, but its connectivity with the medial thalamus, which is the major subcortical target of the precuneus, was decreased. Although the present study did not include a formal evaluation of memory and self-awareness, it is tempting to speculate that in patients with continuous pain there may be a rearrangement of the structures that are involved in episodic memory and this may represent a protective adaptation. This is also in keeping with the reduced connectivity of the amygdale, which is involved in the memory of emotional information [76].

Focusing more closely on the reduced connectivity of the VPL and MDN nuclei of the thalamus observed in the present study, this is consistent with the observation reported above and with a recent study in which subjects with diabetic neuropathic pain had significant reduction in brain metabolites (NAA) in the thalamus compared with the no-pain group [80]. A PET study on complex regional pain syndrome showed that the thalamus seemed to be hypoactive at least in the chronic stage of the disease [81]. This correlates with the increased regional blood flow in the thalamus after motor cortex stimulation for intractable pain (Garda Larrea e Peyron 2007). It has been suggested that activation of the thalamus is a crucial event leading to an eventual, albeit insufficient pain relieving effect.

Modulation of thalamic activity in pain patients should also be interpreted in the light of electrophysiological studies.

Pain, in healthy controls, has the unique capability of disrupting the idling electroencephalogram (EEG) rhythm on a large-scale spatial dimension (i.e. both alpha on the occipital areas and mu-rhythm on the sensorimotor areas [82]). On the contrary, it induces gamma oscillations in the sensory cortex [83]. Since thalamus activity is increased during wakefulness with closed eyes, that is the typical state in which alpha rhythm predominates, it is conceivable that a continuous or long-standing stimulus that disrupts a typical idling rhythm such as pain, may result in the deactivation of the thalamus. This correlates with the reduced large-scale connectivity of all the thalamic regions considered in our study. Activation of the direct and indirect pathways of the basal ganglia has the effect of inducing both excitation of some cortical neuronal populations and inhibition of surrounding neurons, in order to focus the neuronal output on the desired behavior or sensation [84]. This type of neuronal central excitation/surrounding inhibition is mirrored by the already known event-related de-synchronization surrounded by hyper-synchronization [85] and is particularly applicable to the frontal-subcortical circuit encompassing the medial dorsal thalamus, which also involves the striatum and the pallidum [86]. A reduced connectivity of the thalamus may underlie a disruption of these mechanisms and thus may constitute a true dysrhythmic syndrome, as hypothesized by Llinas [87]. Neuropathic pain, as well as other neurological disorders such as Parkinson's disease, is actually considered part of the dysrhythmic syndrome. Thus, combining the results of the present and our previous study, it is conceivable to hypothesize that a dysrhythmic syndrome may be present in diabetic patients with persistent neuropathic pain.

The relative role of pain vs. deafferentation in the induction of the alteration of thalamus and cortex connectivity in our patients remains unresolved. However, some evidence points toward a crucial role of neuropathic pain. Firstly, the selective reduction of activity in the medial dorsal thalamus and ventrolateral nucleus is consistent with a specific pain-related phenomenon since both nuclei are part of the respectively medial and lateral spinothalamic tract. Secondly, even if these patients were affected by polyneuropathy with electrophysiological evidence of involvement of A beta fibers (data not shown), the influence of these alterations on brain function may be hypothesized to be more evident during tasks and may be reflected by a disruption of the so-called optimal feedback control circuit which is made by the primary sensorimotor, posterior parietal and premotor cortices [[88], Wolpert cited above]. Thirdly, Baliki et al. [66] recently demonstrated that a disruption of default mode of brain function occurs in chronic back pain. Thus, our results seem to be related to a common mechanism of disruption of default mode and dysrhythmia of the brain which is shared by various conditions associated with neuropathic pain. However, a possible compromise between the importance of deafferentation vs. that of pain in the modulation of cerebral cortex may be suggested by the theory of the introceptive system [89]. This system comprises the A delta and C fibers, regardless of their fine specificity (thermoceptive, nociceptive) or site of innervation (somatic vs. visceral), the lamin spinal I neurons, reticular nuclei of the brainstem such as the A1 and A5-7 monoaminargic group and parabrachial nucleus, as well as the medial dorsal and ventrobasal thalamus and insular and anterior cingulated cortex. This system is implicated in the evaluation of the homeostatic control of the entire body and is not only involved in sensory perception, but also modulates sympathetic responses and it is also hypothesized to modulate more complex motor behavior. From this perspective, a profound modulation of thalamic connectivity at rest as reported herein is conceivable.

Taking this view the left lateralization of the S1-thalamic connectivity reduction seems difficult to explain, however as shown by the pairwise correlation (see fig. 4) the right S1-Thalamic connectivity, while not significant, is similar to that of the left and just below the threshold: we hypothesize that with a bigger group of patients we would also find the right connectivity to be significantly reduced.

From this perspective the distinction between brain function modulation induced by pain and by deafferentation may be less meaningful.


In summary, exploring thalamocortical functional connectivity in a group of eight patients suffering from peripheral neuropathic pain (diabetic pain), and compared it with that of a group of healthy subjects we found a reduced synchrony between S1 and some of the thalamic nuclei involved in pain processing. Studying specifically the VPL and MDN connectivity we demonstrated that VPL has a reduced synchrony with sensorymotor, frontal, cingulated and parietal areas while MDN more with sensory and emotional processing areas. This supports the idea that chronic pain can alter thalamocortical connections causing a disruption of thalamic feedback, and the view of chronic pain as a thalamocortical dysrhythmia.

Methodological Consideration

This study enrolled only a limited number of patients so the inference on a general population may be partly impaired. However, as mentioned earlier, our results are consistent with previous results.

The distinction between modulation induced by pain or by deafferentation, although less meaningful in the light of the theory of the introceptive system, should be analyzed in further studies comparing patients with diabetic neuropathy with and without pain.

A major point remains unresolved: is the reduced connectivity of the medial dorsal thalamus an adaptive or misadaptive plastic change? Although it continues to be a matter of debate, the medial dorsal thalamus together with other limbic structures (fornix) seems to play an important role both in retaining and in re-learning encoded information [90, 91]. Therefore, further studies should address the issue of learning and memory as well as depression in diabetic patients with painful polyneuropathy and their relationship with alteration of default mode of brain function.



blood oxygenation level dependent


cerebrospinal fluid


global signal


low-frequency band fluctuations


medial dorsal nucleus


neuropathic pain


region of interest


functional connectivity


resting state functional connectivity


primary somatosensory area


ventral posterior lateral nucleus


white matter.


  1. 1.

    Uhlhaas PJ, Singer W: Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology. Neuron. 2006, 52 (1): 155-168. 10.1016/j.neuron.2006.09.020.

  2. 2.

    Ohara S, Crone NE, Weiss N, Lenz FA: Analysis of synchrony demonstrates 'pain networks' defined by rapidly switching, task-specific, functional connectivity between pain-related cortical structures. Pain. 2006, 123 (3): 244-253. 10.1016/j.pain.2006.02.012.

  3. 3.

    Apkarian AV, Chialvo DR: The shadows of pain. Pain. 2006, 123 (3): 221-222. 10.1016/j.pain.2006.04.028.

  4. 4.

    Chen AC, Herrmann CS: Perception of pain coincides with the spatial expansion of electroencephalographic dynamics in human subjects. Neurosci Lett. 2001, 297 (3): 183-186. 10.1016/S0304-3940(00)01696-7.

  5. 5.

    Llinas RR, Ribary U, Jeanmonod D, Kronberg E, Mitra PP: Thalamocortical dysrhythmia: A neurological and neuropsychiatric syndrome characterized by magnetoencephalography. Proc Natl Acad Sci USA. 1999, 96 (26): 15222-15227. 10.1073/pnas.96.26.15222.

  6. 6.

    Brambilla M, Manuguerra M, Valente M, Caramenti GC, Sotgiu ML, Biella GE: Chronic pain as expression of neural substrates. Issues from the neuronal dynamics and mutual relations. Arch Ital Biol. 2004, 142 (3): 275-283.

  7. 7.

    Lenz FA, Treede RD: Attention, novelty, and pain. Pain. 2002, 99 (1-2): 1-3. 10.1016/S0304-3959(02)00050-7.

  8. 8.

    Eccleston C, Crombez G: Pain demands attention: a cognitive-affective model of the interruptive function of pain. Psychol Bull. 1999, 125 (3): 356-366. 10.1037/0033-2909.125.3.356.

  9. 9.

    McCabe C, Lewis J, Shenker N, Hall J, Cohen H, Blake D: Don't look now! Pain and attention. Clin Med. 2005, 5 (5): 482-486.

  10. 10.

    Villemure C, Bushnell MC: Cognitive modulation of pain: how do attention and emotion influence pain processing?. Pain. 2002, 95 (3): 195-199. 10.1016/S0304-3959(02)00007-6.

  11. 11.

    Crombez G, Eccleston C, Baeyens F, Eelen P: Habituation and the interference of pain with task performance. Pain. 1997, 70 (2-3): 149-154. 10.1016/S0304-3959(96)03304-0.

  12. 12.

    Houlihan ME, McGrath PJ, Connolly JF, Stroink G, Allen Finley G, Dick B, Phi TT: Assessing the effect of pain on demands for attentional resources using ERPs. Int J Psychophysiol. 2004, 51 (2): 181-187. 10.1016/j.ijpsycho.2003.08.001.

  13. 13.

    Lorenz J, Beck H, Bromm B: Cognitive performance, mood and experimental pain before and during morphine-induced analgesia in patients with chronic non-malignant pain. Pain. 1997, 73 (3): 369-375. 10.1016/S0304-3959(97)00123-1.

  14. 14.

    Dick B, Eccleston C, Crombez G: Attentional functioning in fibromyalgia, rheumatoid arthritis, and musculoskeletal pain patients. Arthritis Rheum. 2002, 47 (6): 639-644. 10.1002/art.10800.

  15. 15.

    Eccleston C: Chronic pain and distraction: an experimental investigation into the role of sustained and shifting attention in the processing of chronic persistent pain. Behav Res Ther. 1995, 33 (4): 391-405. 10.1016/0005-7967(94)00057-Q.

  16. 16.

    Harman K, Ruyak P: Working through the pain: a controlled study of the impact of persistent pain on performing a computer task. Clin J Pain. 2005, 21 (3): 216-222. 10.1097/00002508-200505000-00004.

  17. 17.

    Kewman DG, Vaishampayan N, Zald D, Han B: Cognitive impairment in musculoskeletal pain patients. Int J Psychiatry Med. 1991, 21 (3): 253-262.

  18. 18.

    Park H, Sprince NL, Whitten PS, Burmeister LF, Zwerling C: Risk factors for back pain among male farmers: analysis of Iowa Farm Family Health and Hazard Surveillance Study. Am J Ind Med. 2001, 40 (6): 646-654. 10.1002/ajim.10025.

  19. 19.

    Veldhuijzen DS, Kenemans JL, van Wijck AJ, Olivier B, Kalkman CJ, Volkerts ER: Processing capacity in chronic pain patients: a visual event-related potentials study. Pain. 2006, 121 (1-2): 60-68. 10.1016/j.pain.2005.12.004.

  20. 20.

    Bruggemann J, Galhardo V, Apkarian AV: Immediate reorganization of the rat somatosensory thalamus after partial ligation of sciatic nerve. J Pain. 2001, 2 (4): 220-228. 10.1054/jpai.2001.22792.

  21. 21.

    Fukui S, Matsuno M, Inubushi T, Nosaka S: N-Acetylaspartate concentrations in the thalami of neuropathic pain patients and healthy comparison subjects measured with (1)H-MRS. Magn Reson Imaging. 2006, 24 (1): 75-79. 10.1016/j.mri.2005.10.021.

  22. 22.

    Nakabeppu Y, Nakajo M, Gushiken T, Tsuchimochi S, Tani A, Kanmura Y: Decreased perfusion of the bilateral thalami in patients with chronic pain detected by Tc-99 m-ECD SPECT with statistical parametric mapping. Ann Nucl Med. 2001, 15 (5): 459-463. 10.1007/BF02988354.

  23. 23.

    Pattany PM, Yezierski RP, Widerstrom-Noga EG, Bowen BC, Martinez-Arizala A, Garcia BR, Quencer RM: Proton magnetic resonance spectroscopy of the thalamus in patients with chronic neuropathic pain after spinal cord injury. AJNR Am J Neuroradiol. 2002, 23 (6): 901-905.

  24. 24.

    Weng HR, Lenz FA, Vierck C, Dougherty PM: Physiological changes in primate somatosensory thalamus induced by deafferentation are dependent on the spinal funiculi that are sectioned and time following injury. Neuroscience. 2003, 116 (4): 1149-1160. 10.1016/S0306-4522(02)00796-0.

  25. 25.

    Iadarola MJ, Max MB, Berman KF, Byas-Smith MG, Coghill RC, Gracely RH, Bennett GJ: Unilateral decrease in thalamic activity observed with positron emission tomography in patients with chronic neuropathic pain. Pain. 1995, 63 (1): 55-64. 10.1016/0304-3959(95)00015-K.

  26. 26.

    Raichle ME, Snyder AZ: A default mode of brain function: a brief history of an evolving idea. Neuroimage. 2007, 37 (4): 1083-1090. 10.1016/j.neuroimage.2007.02.041.

  27. 27.

    Damoiseaux JS, Rombouts SA, Barkhof F, Scheltens P, Stam CJ, Smith SM, Beckmann CF: Consistent resting-state networks across healthy subjects. Proc Natl Acad Sci USA. 2006, 103 (37): 13848-13853. 10.1073/pnas.0601417103.

  28. 28.

    Greicius MD, Krasnow B, Reiss AL, Menon V: Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci USA. 2003, 100 (1): 253-258. 10.1073/pnas.0135058100.

  29. 29.

    Hampson M, Peterson BS, Skudlarski P, Gatenby JC, Gore JC: Detection of functional connectivity using temporal correlations in MR images. Hum Brain Mapp. 2002, 15 (4): 247-262. 10.1002/hbm.10022.

  30. 30.

    Vincent JL, Patel GH, Fox MD, Snyder AZ, Baker JT, Van Essen DC, Zempel JM, Snyder LH, Corbetta M, Raichle ME: Intrinsic functional architecture in the anaesthetized monkey brain. Nature. 2007, 447 (7140): 83-86. 10.1038/nature05758.

  31. 31.

    Friston KJ, Frith CD, Liddle PF, Frackowiak RS: Functional connectivity: the principal-component analysis of large (PET) data sets. J Cereb Blood Flow Metab. 1993, 13 (1): 5-14.

  32. 32.

    Horwitz B: The elusive concept of brain connectivity. Neuroimage. 2003, 19 (2 Pt 1): 466-470. 10.1016/S1053-8119(03)00112-5.

  33. 33.

    Biswal B, Yetkin FZ, Haughton VM, Hyde JS: Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995, 34 (4): 537-541. 10.1002/mrm.1910340409.

  34. 34.

    Fox MD, Raichle ME: Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci. 2007, 8 (9): 700-711. 10.1038/nrn2201.

  35. 35.

    Beckmann CF, DeLuca M, Devlin JT, Smith SM: Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1457): 1001-1013. 10.1098/rstb.2005.1634.

  36. 36.

    De Luca M, Beckmann CF, De Stefano N, Matthews PM, Smith SM: fMRI resting state networks define distinct modes of long-distance interactions in the human brain. Neuroimage. 2006, 29 (4): 1359-1367. 10.1016/j.neuroimage.2005.08.035.

  37. 37.

    Franco AR, Pritchard A, Calhoun VD, Mayer AR: Interrater and intermethod reliability of default mode network selection. Hum Brain Mapp. 2009, 30 (7): 2293-303. 10.1002/hbm.20668.

  38. 38.

    Shehzad Z, Kelly AM, Reiss PT, Gee DG, Gotimer K, Uddin LQ, Lee SH, Margulies DS, Roy AK, Biswal BB, et al.: The Resting Brain: Unconstrained yet Reliable. Cereb Cortex. 2009, 19 (10): 2209-29. 10.1093/cercor/bhn256.

  39. 39.

    van Buuren M, Gladwin TE, Zandbelt BB, Heuvel van den M, Ramsey NF, Kahn RS, Vink M: Cardiorespiratory effects on default-mode network activity as measured with fMRI. Hum Brain Mapp. 2009, 30 (9): 3031-42. 10.1002/hbm.20729.

  40. 40.

    Jones EG: The thalamus. 2007, Cambridge: Cambridge University Press, 2

  41. 41.

    Sherman SM, Guillery RW, Sherman SMEtt: Exploring the thalamus and its role in cortical function. 2006, Cambridge, Mass; London: The MIT Press, 2

  42. 42.

    Anand A, Li Y, Wang Y, Wu J, Gao S, Bukhari L, Mathews VP, Kalnin A, Lowe MJ: Activity and connectivity of brain mood regulating circuit in depression: a functional magnetic resonance study. Biol Psychiatry. 2005, 57 (10): 1079-1088. 10.1016/j.biopsych.2005.02.021.

  43. 43.

    Dosenbach NU, Fair DA, Miezin FM, Cohen AL, Wenger KK, Dosenbach RA, Fox MD, Snyder AZ, Vincent JL, Raichle ME, et al.: Distinct brain networks for adaptive and stable task control in humans. Proc Natl Acad Sci USA. 2007, 104 (26): 11073-11078. 10.1073/pnas.0704320104.

  44. 44.

    Greicius MD, Flores BH, Menon V, Glover GH, Solvason HB, Kenna H, Reiss AL, Schatzberg AF: Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biol Psychiatry. 2007, 62 (5): 429-437. 10.1016/j.biopsych.2006.09.020.

  45. 45.

    Seeley WW, Menon V, Schatzberg AF, Keller J, Glover GH, Kenna H, Reiss AL, Greicius MD: Dissociable intrinsic connectivity networks for salience processing and executive control. J Neurosci. 2007, 27 (9): 2349-2356. 10.1523/JNEUROSCI.5587-06.2007.

  46. 46.

    Stein T, Moritz C, Quigley M, Cordes D, Haughton V, Meyerand E: Functional connectivity in the thalamus and hippocampus studied with functional MR imaging. AJNR Am J Neuroradiol. 2000, 21 (8): 1397-1401.

  47. 47.

    Zhang D, Snyder AZ, Fox MD, Sansbury MW, Shimony JS, Raichle ME: Intrinsic functional relations between human cerebral cortex and thalamus. J Neurophysiol. 2008, 100 (4): 1740-1748. 10.1152/jn.90463.2008.

  48. 48.

    Welsh RC, Chen AC, Taylor SF: Low-Frequency BOLD Fluctuations Demonstrate Altered Thalamocortical Connectivity in Schizophrenia. Schizophr Bull. 2008.

  49. 49.

    Laureys S, Faymonville ME, Luxen A, Lamy M, Franck G, Maquet P: Restoration of thalamocortical connectivity after recovery from persistent vegetative state. Lancet. 2000, 355 (9217): 1790-1791. 10.1016/S0140-6736(00)02271-6.

  50. 50.

    White NS, Alkire MT: Impaired thalamocortical connectivity in humans during general-anesthetic-induced unconsciousness. Neuroimage. 2003, 19 (2 Pt 1): 402-411. 10.1016/S1053-8119(03)00103-4.

  51. 51.

    Mizuno A, Villalobos ME, Davies MM, Dahl BC, Muller RA: Partially enhanced thalamocortical functional connectivity in autism. Brain Res. 2006, 1104 (1): 160-174. 10.1016/j.brainres.2006.05.064.

  52. 52.

    Cauda F, Sacco K, Duca S, Cocito D, D'Agata F, Geminiani GC, Canavero S: Altered resting state in diabetic neuropathic pain. PLoS ONE. 2009, 4 (2): e4542-10.1371/journal.pone.0004542.

  53. 53.

    Friston KJ, Tononi G, Reeke GN, Sporns O, Edelman GM: Value-dependent selection in the brain: simulation in a synthetic neural model. Neuroscience. 1994, 59 (2): 229-243. 10.1016/0306-4522(94)90592-4.

  54. 54.

    Friston KJ, Price CJ: Generative models, brain function and neuroimaging. Scand J Psychol. 2001, 42 (3): 167-177. 10.1111/1467-9450.00228.

  55. 55.

    Miezin FM, Maccotta L, Ollinger JM, Petersen SE, Buckner RL: Characterizing the hemodynamic response: effects of presentation rate, sampling procedure, and the possibility of ordering brain activity based on relative timing. Neuroimage. 2000, 11 (6 Pt 1): 735-759. 10.1006/nimg.2000.0568.

  56. 56.

    Talairach J, Tournoux P: Co-planar stereotaxic atlas of the human brain: 3-dimensional proportional system: an approach to cerebral imaging. 1988, Stuttgart: Thieme

  57. 57.

    Genovese CR, Lazar NA, Nichols T: Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage. 2002, 15 (4): 870-878. 10.1006/nimg.2001.1037.

  58. 58.

    Moore CI, Stern CE, Corkin S, Fischl B, Gray AC, Rosen BR, Dale AM: Segregation of somatosensory activation in the human rolandic cortex using fMRI. J Neurophysiol. 2000, 84 (1): 558-569.

  59. 59.

    Lancaster JL, Woldorff MG, Parsons LM, Liotti M, Freitas CS, Rainey L, Kochunov PV, Nickerson D, Mikiten SA, Fox PT: Automated Talairach atlas labels for functional brain mapping. Hum Brain Mapp. 2000, 10 (3): 120-131. 10.1002/1097-0193(200007)10:3<120::AID-HBM30>3.0.CO;2-8.

  60. 60.

    Bandettini PA, Bullmore E: Endogenous oscillations and networks in functional magnetic resonance imaging. Hum Brain Mapp. 2008, 29 (7): 737-739. 10.1002/hbm.20607.

  61. 61.

    Napadow V, Dhond R, Conti G, Makris N, Brown EN, Barbieri R: Brain correlates of autonomic modulation: combining heart rate variability with fMRI. Neuroimage. 2008, 42 (1): 169-177. 10.1016/j.neuroimage.2008.04.238.

  62. 62.

    Fox MD, Zhang D, Snyder AZ, Raichle ME: The Global Signal and Observed Anticorrelated Resting State Brain Networks. J Neurophysiol. 2009, 3270-83. 10.1152/jn.90777.2008. 6

  63. 63.

    Woolrich MW, Ripley BD, Brady M, Smith SM: Temporal autocorrelation in univariate linear modeling of FMRI data. Neuroimage. 2001, 14 (6): 1370-1386. 10.1006/nimg.2001.0931.

  64. 64.

    Friston KJ: Statistical parametric mapping: the analysis of functional brain images. 2007, London: Academic, vii: 647-p., [632] p. of plates: ill. (some col.); 629 cm

  65. 65.

    Llewelyn JG, Tomlinson DR, Thomas PK, Dyck PJ, Thomas PK: Diabetic Neuropathies. Peripheral neuropathy. Edited by: Dyck PJ, Thomas PK. 2005, Philadelphia, Pa; [London]: Elsevier Saunders, 1951-1991. 4

  66. 66.

    Baliki MN, Geha PY, Apkarian AV, Chialvo DR: Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics. J Neurosci. 2008, 28 (6): 1398-1403. 10.1523/JNEUROSCI.4123-07.2008.

  67. 67.

    Jankovic J, Linden Van der C: Dystonia and tremor induced by peripheral trauma: predisposing factors. J Neurol Neurosurg Psychiatry. 1988, 51 (12): 1512-1519. 10.1136/jnnp.51.12.1512.

  68. 68.

    Ziemann U, Hallett M, Cohen LG: Mechanisms of deafferentation-induced plasticity in human motor cortex. J Neurosci. 1998, 18 (17): 7000-7007.

  69. 69.

    Woolf CJ, Salter MW: Neuronal plasticity: increasing the gain in pain. Science. 2000, 288 (5472): 1765-1769. 10.1126/science.288.5472.1765.

  70. 70.

    Rossi P, Pierelli F, Parisi L, Perrotta A, Bartolo M, Amabile G, Serrao M: Effect of painful heterotopic stimulation on the cutaneous silent period in the upper limbs. Clin Neurophysiol. 2003, 114 (1): 1-6. 10.1016/S1388-2457(02)00321-8.

  71. 71.

    Trompetto C, Buccolieri A, Abbruzzese G: Intracortical inhibitory circuits and sensory input: a study with transcranial magnetic stimulation in humans. Neurosci Lett. 2001, 297 (1): 17-20. 10.1016/S0304-3940(00)01648-7.

  72. 72.

    Apkarian AV, Bushnell MC, Treede RD, Zubieta JK: Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain. 2005, 9 (4): 463-484. 10.1016/j.ejpain.2004.11.001.

  73. 73.

    Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL: A default mode of brain function. Proc Natl Acad Sci USA. 2001, 98 (2): 676-682. 10.1073/pnas.98.2.676.

  74. 74.

    Lorenz J, Minoshima S, Casey KL: Keeping pain out of mind: the role of the dorsolateral prefrontal cortex in pain modulation. Brain. 2003, 126 (Pt 5): 1079-1091. 10.1093/brain/awg102.

  75. 75.

    Garcia-Larrea L, Peyron R: Motor cortex stimulation for neuropathic pain: From phenomenology to mechanisms. Neuroimage. 2007, 37 (Suppl 1): S71-79. 10.1016/j.neuroimage.2007.05.062.

  76. 76.

    Phillips ML, Drevets WC, Rauch SL, Lane R: Neurobiology of emotion perception I: The neural basis of normal emotion perception. Biol Psychiatry. 2003, 54 (5): 504-514. 10.1016/S0006-3223(03)00168-9.

  77. 77.

    Argoff CE, Cole BE, Fishbain DA, Irving GA: Diabetic peripheral neuropathic pain: clinical and quality-of-life issues. Mayo Clin Proc. 2006, 81 (4 Suppl): S3-11.

  78. 78.

    Wolpert DM, Goodbody SJ, Husain M: Maintaining internal representations: the role of the human superior parietal lobe. Nat Neurosci. 1998, 1 (6): 529-533. 10.1038/2245.

  79. 79.

    Cavanna AE, Trimble MR: The precuneus: a review of its functional anatomy and behavioural correlates. Brain. 2006, 129 (Pt 3): 564-583. 10.1093/brain/awl004.

  80. 80.

    Sorensen L, Siddall PJ, Trenell MI, Yue DK: Differences in metabolites in pain-processing brain regions in patients with diabetes and painful neuropathy. Diabetes Care. 2008, 31 (5): 980-981. 10.2337/dc07-2088.

  81. 81.

    Garcia-Larrea L, Maarrawi J, Peyron R, Costes N, Mertens P, Magnin M, Laurent B: On the relation between sensory deafferentation, pain and thalamic activity in Wallenberg's syndrome: a PET-scan study before and after motor cortex stimulation. Eur J Pain. 2006, 10 (8): 677-688. 10.1016/j.ejpain.2005.10.008.

  82. 82.

    Ploner M, Gross J, Timmermann L, Pollok B, Schnitzler A: Pain suppresses spontaneous brain rhythms. Cereb Cortex. 2006, 16 (4): 537-540. 10.1093/cercor/bhj001.

  83. 83.

    Gross J, Schnitzler A, Timmermann L, Ploner M: Gamma oscillations in human primary somatosensory cortex reflect pain perception. PLoS Biol. 2007, 5 (5): e133-10.1371/journal.pbio.0050133.

  84. 84.

    Takakusaki K, Saitoh K, Harada H, Kashiwayanagi M: Role of basal ganglia-brainstem pathways in the control of motor behaviors. Neurosci Res. 2004, 50 (2): 137-151. 10.1016/j.neures.2004.06.015.

  85. 85.

    Pfurtscheller G, Lopes da Silva FH: Event-related desynchronization (ERD) and event-related synchronization (ERS). Electroencephalography: basic principles, clinical applications, and related fields. Edited by: Niedermeyer E, Lopes da Silva FH. 2005, Philadelphia; London: Lippincott Williams & Wilkins, 1003-1016. 5

  86. 86.

    Cummings JL: Frontal-subcortical circuits and human behavior. Arch Neurol. 1993, 50 (8): 873-880.

  87. 87.

    Llinas RR, Steriade M: Bursting of thalamic neurons and states of vigilance. J Neurophysiol. 2006, 95 (6): 3297-3308. 10.1152/jn.00166.2006.

  88. 88.

    Scott SH: Optimal feedback control and the neural basis of volitional motor control. Nat Rev Neurosci. 2004, 5 (7): 532-546. 10.1038/nrn1427.

  89. 89.

    Craig AD: Interoception: the sense of the physiological condition of the body. Curr Opin Neurobiol. 2003, 13 (4): 500-505. 10.1016/S0959-4388(03)00090-4.

  90. 90.

    Miller LA, Caine D, Watson JD: A role for the thalamus in memory for unique entities. Neurocase. 2003, 9 (6): 504-514. 10.1076/neur.9.6.504.29375.

  91. 91.

    Mitchell AS, Browning PG, Wilson CR, Baxter MG, Gaffan D: Dissociable roles for cortical and subcortical structures in memory retrieval and acquisition. J Neurosci. 2008, 28 (34): 8387-8396. 10.1523/JNEUROSCI.1924-08.2008.

Download references


We wish to thank all the subjects who participated in this study. This research was supported by the Italian Ministry of Education, University and Research (MIUR), PRIN project 2008, reg. No. 2007P2CWAZ.

Author information



Corresponding author

Correspondence to Franco Cauda.

Additional information

Authors' contributions

Experiments conceived and designed by: CF SK DS CD. Experiments performed by: CF SK DS CD. Data analyzed by: CF SK DF FM. Reagents/materials/analysis tools contributed by: CF SK DS CD FM DF GG. Paper written by: CF GI SK DF. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Cauda, F., Sacco, K., D'Agata, F. et al. Low-frequency BOLD fluctuations demonstrate altered thalamocortical connectivity in diabetic neuropathic pain. BMC Neurosci 10, 138 (2009).

Download citation


  • Neuropathic Pain
  • Functional Connectivity
  • Anterior Cingulated Cortex
  • Default Mode Network
  • Blood Oxygenation Level Dependent