Abstract
Chronic pain is a widespread health condition affecting millions of individuals worldwide. Despite its high prevalence, the underlying neural mechanisms and pathophysiology of chronic pain are still not well understood. Introduced about 35 years ago, brain diffusion MRI has emerged as a powerful tool to investigate white matter microstructure and connectivity changes in chronic pain conditions. In this review, we provide an overview of the current state of research on brain diffusion MRI in chronic pain, including the methods used to acquire and analyze diffusion MRI data, as well as the major findings in the field.
We discuss the evidence supporting the role of altered white matter microstructure and connectivity in chronic pain conditions, highlighting the importance of studying multiple chronic pain syndromes to identify common neurobiological pathways. We also explore the potential clinical applications of diffusion MRI in chronic pain research, such as identifying biomarkers for diagnosis and treatment response.
In addition, we explore the shortcomings and difficulties associated with brain diffusion MRI in chronic pain studies, such as the necessity for uniformity in data acquisition and analysis methods. Finally, we discuss emerging techniques and future directions in brain diffusion MRI research that may provide new insights into the pathophysiology of chronic pain and potential new therapeutic targets. In summary, we are optimistic that brain diffusion MRI has significant potential in furthering our knowledge of chronic pain and enhancing clinical outcomes, although further research is required due to the limited body of literature currently available, and that targeted therapeutic strategies are yet to be identified.
Introduction
The understanding of chronic pain (CP) has greatly advanced over the past 25 years, with many scientists recognizing the significant role of brain processes in its development (1). A vast number of studies have now thoroughly documented the impact of various CP conditions on brain structure and function (2). When assessed using structural Magnetic Resonance Imaging (MRI) techniques, such as T1 or T2 weighted images, some brain impairments were found to be common across different pain conditions, such as reduced gray matter density (GMD) in the hippocampus and secondary somatosensory cortex. However, others were unique to specific types of CP, such as decreased GMD in the posterior insula in low back pain (LBP) compared to decreased GMD in the anterior insula in chronic regional pain syndrome (3). These studies have led to the unexpected identification of brain regions involved in chronic pain, including areas responsible for emotional learning, memory, and decision-making, rather than just the anticipated somatosensory regions. These changes may play a role in both the development and maintenance of CP, as seen in individuals transitioning to chronic low back pain (CLBP), who present smaller amygdala and hippocampi volumes (4).
Functional MRI (fMRI) studies have also indicated that regions outside of the somatosensory areas are largely involved in chronic pain. Areas known for processing cognitive, learning, and emotional information are found to have increased activity during the transition from acute to chronic pain (5). Chronic pain has also been shown to affect the organization of the brain at the network level. When studying resting state networks in various chronic pain conditions, research has found that the default mode network is the primary network undergoing reorganization due to an increase in sustained pain signaling through the medial prefrontal cortex (6). Some functional features are also capable of forecasting transition from a subacute phase to a chronic pain phase (7). Interestingly, these properties appear to be plastic and may partially reverse and reorganize upon receiving treatment or during the progression of the condition (8–10).
However, despite the tremendous progress made in MRI technology and analysis methods, the structural changes in white matter (WM) related to CP have been studied to a lesser extent. White matter is a highly complex system that makes up over one-third of the brain’s total volume and recent research has shown that various acute and chronic neurological conditions, such as Alzheimer’s, Parkinson’s, depression, multiple sclerosis, traumatic brain injury, and others, can affect the microstructure of white matter (11–16). Currently, diffusion-weighted MRI (DWI) is the most promising non-invasive tool for investigating WM structural changes. Although dMRI methods have advanced significantly in recent years (17), one of the earliest and still commonly used methods is to model the dMRI signal using diffusion tensor imaging (DTI). Accordingly, fractional anisotropy (FA), a common DTI measure, is one of the most widely reported dMRI measures to evaluate WM integrity. As most other fields, the pain research community is searching for valid biomarkers to help people suffering from this condition. For CP, we are searching for biomarkers that can predict the transition from acute to chronic pain and identify treatment response potential, and brain imaging is considered one of the most valid approaches (18–23). However, biomarkers from diffusion MRI for CP are currently under-investigated, but they hold promise for rapid clinical application (24–28). It is likely that CP biomarkers will soon be discovered using diffusion MRI due to the critical role that the brain plays in CP and its treatment.
The purpose of this review is to critically evaluate the current research on the characteristics of brain white matter as revealed through diffusion MRI methods and its impairments in various chronic pain conditions. We have used the latest definition of chronic pain provided by the International Association for the Study of Pain and its subcategories to classify the types of chronic pain (29,30). The focus was on chronic primary pain, which encompasses five categories: 1) chronic widespread pain, 2) complex regional pain syndrome, 3) chronic primary headache or orofacial pain, 4) chronic primary visceral pain, and 5) chronic primary musculoskeletal pain. After reviewing 447 publications, only 59 met the criteria, with the highest number of articles in category 3 (21 articles), and only two articles in category 2.
Methods
Information source
PubMed and Scopus databases were searched for articles dated up to 10th March 2022 Articles were also searched using the Medical Subject Headings (MeSH) term on PubMed.
Review strategy
This review follows the recommendations from the Preferred Reporting Items for Systematic reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR) Checklist. Search strategies were developed with a librarian of the Health Sciences Library of the Université de Sherbrooke. The keywords chosen for the review were: magnetic resonance imaging, diffusion, pain, and brain. The full search strategy can be viewed in the supplementary material.
Study strategy and inclusion/exclusion criterias
The study strategy is illustrated in Figure 1. First, we screened 447 articles to exclude review articles and articles that were ineligible based on the title and abstract (i.e., articles that did not contain the words “diffusion,” “brain,” or “pain,”). Second, we excluded case reports, articles not in English, and articles where pain was not the main focus. Third, we excluded studies that did not acquire dMRI in brain regions of individuals with chronic pain and studies that included participants who were under the age of one, veterans, amputees, or had undergone deep brain stimulation (DBS), traumatic brain injury (TBI), surgery, chronic pain from a stroke or neurodegenerative disease, or experimental acute pain. Methodological articles and animal studies were also excluded. Fourth, among the remaining articles, we separated them into three categories: chronic primary pain, chronic secondary pain based on the IASP classification (29,30), or articles missing information about the chronic pain condition, and only included articles on chronic primary pain. Finally, articles were excluded if dMRI analysis was done only on a nerve, only on healthy controls, or if no results were documented. One article was added after further investigation, as it was mentioned in several of the included articles and was found to have used dMRI in chronic pain patients but did not used our search terms in either the title or the abstract. The remaining articles were analyzed.
Study selection and analyzes
The first three steps described in the section study strategy and inclusion and exclusion criteria were conducted by MS and PT. Any disagreement was resolved between MS and PT. Subsequently, the classification into primary chronic pain and secondary chronic pain was made by MS and supported by PT and GLe. The remaining articles were then separated into five categories based on the classification of IASP: 1) chronic widespread pain; 2) complex regional pain syndrome; 3) chronic primary headache or orofacial pain; 4) chronic primary visceral pain and 5) chronic primary musculosketal pain. If applicable, half of the articles in each category were separated and respectively analyzed by MS and PT. The remaining articles were analyzed by MM and GLe. Finally, each reviewer (MS, PT, MM and GLe) extracted, based on the chart developed by PT, general study design information and dMRI-specific study data. An overview of the chart used by each reviewer is outlined in table 1.
Results
From the 370 unique articles initially identified, 247 articles remained after removing reviews and articles which did not fit our inclusion criteria for the title and abstract. Subsequent removal of: case-reports, articles not available in English and articles that only reported pain anecdotally yielded a total of 175 articles. Further exclusion of articles that did not acquire dMRI in brain regions or on a specified chronic pain condition population, left 107 articles. Then, the remaining articles were separated into three categories: chronic primary pain (62 articles), chronic secondary pain (41 articles) and articles missing information about the chronic pain condition (4 articles); only chronic primary pain articles were kept. Afterwards, we excluded articles where the dMRI was acquired: on a nerve, only on healthy participants and articles without results, leaving 58 articles. Additionally, one article was added at the end because it was mentioned several times by remaining articles, giving a total of 59 articles for the final analysis.
Chronic primary headache or orofacial pain
The IASP classification defines chronic primary headache or orofacial pain as headache or orofacial pain that occurs on at least 15 days per month for a period longer than 3 months. Daily pain duration is at least 2 hours (untreated) or may occur in several shorter attacks (29).
This category yielded 21 articles from which: 11 studies were on migraine (31–41), two were on cluster headache (42,43) one on chronic headache (44), three on temporomandibular disorder (TMD) (45–47) and four on burning mouth syndrome (48– 51)
For the migraine studies, four studies used TBSS to generate whole brain FA skeleton (31,37,39,40). Neeb et coll found no difference in FA, MD, AD, and RD between migraine patients and HC. The other three studies further analyzed the TBSS results using tractography. Gomez-Beldarrain et coll., using a probabilistic white matter tract atlas, found lower FA in the right anterior insula, bilateral cingulate gyri and right uncinate fasciculus. Szabó et coll., using probabilistic tractography with a seeding in regions showing group differences from TBSS, found lower FA, higher MD and RD in the right frontal WM. Kattem Husøy et coll., using deterministic tractography with a whole brain seeding (FA > 0.25), studied patients with varying types of primary headaches, but only migraine and tension-type headache results were analyzed. They found higher AD in several areas of the skeleton, such as the CC and corticospinal tract (CST), in migraine and tension-type headaches groups. No differences were observed for FA, MD, and RD
Three studies have used TBSS to generate FA maps only (32,35,38). One study has mentioned using a linear model based (Basser and Pierpaoli, 1996) to generate FA maps (34). Zhang et coll, using SPM on diffusion metric maps registered to MNI, found no difference in FA, RD, and AD. Coppola et coll. have found notably higher RD in various WM tracts such as, the bilateral genu of the CC and the bilateral posterior limb of the IC. DaSilva et coll. have found lower FA in two regions of the brain in the posterior limb of the IC and the CR. Marciszewski et coll. have evaluated FA, AD, MD, and RD during different phases of the cycle of a migraine. Results showed variation of MD during the cycle of migraine in different regions such as the spinal trigeminal nucleus. Also, FA was higher in the region of the medial lemniscus/ventral trigeminal thalamic tract in the entire migraine cycle.
Three studies have used only tractography (33,36,41). Chong et coll. used TRACULA from Free Surfer to develop a classification model to distinguish migraine patients from persistent post-traumatic headache patients such that, no diffusion metrics were compared. Planchuelo-Gómez evaluated structural connectivity using the MRtrix tools. No diffusion metrics were compared, but they found an increase in streamline count in connections between subcortical regions, the insula, and the cingulate and orbitofrontal cortex and between the insula and the temporal region. Silvestro et coll. have used Bedpostx from FSL to investigate the brain connectome. No diffusion metrics were compared, but the connectome showed higher connection probability in precuneus, cuneus, amygdala, calcarine cortex, posterior cingulate cortex, anterior cingulate cortex, postcentral gyrus, superior parietal lobule, lingual and fusiform gyri, middle frontal gyrus and inferior and superior parietal lobules.
For the cluster headache studies, both groups used TBSS to generate whole brain FA skeleton (42,43). Szabó et coll. reported diffusion metrics on the FA skeleton (FA, MD, AD, RD) and showed a decrease of FA in several WM tracts, such as the corpus callosum (CC) and the corona radiata (CR). While Teepker et coll. reported changes in FA in different structures, such as the internal capsule (IC) and thalamus, without specifying if they observed an increase or decrease of FA. Results from these two studies are therefore difficult to compare.
For the chronic headache study, the authors used constrained spherical deconvolution (CSD) and semi-automated tractography to compute whole brain tractograms, then extracted FA, RD and AD from two specific tracks, the cingulum and the uncinate fasciculus (44). FA and RD were respectively higher and lower in chronic headaches patients compared to healthy controls (HC) in the cingulum. No difference was observed in the uncinate fasciculus and for the AD metric on both bundles.
For the TMD studies, two studies used TBSS to generate whole brain FA skeleton and probabilistic tractography (45,46). Salomons et coll. evaluated the relation between FA of a corticospinal track and behavioral aspects, however, no comparison between groups regarding diffusion metrics was performed. In the other TBSS TMD study, Moayedi et coll. have showed lower FA notably in the right IC and the right external/extreme capsule adjacent to the insula. They also found, in the CC, higher connection probability to the frontal pole and lower connection probability to the dorsolateral prefrontal cortex. The authors did also dMRI on the trigeminal nerve directly, but those results were not included in this review.
Gustin et coll. have generated FA maps using classic Basser and Pierpaoli method (52) without referencing any specific software (47). FA was computed to determine whether the primary somatosensory cortex (S1) reorganization was associated with a change in regional anatomy in TMD and painful trigeminal neuralgia (PTN). They found that the S1 face region FA was lower in the PTN compared to TMD and controls, but no difference was found between TMD and controls. For the burning mouth syndrome studies, two articles have generated whole brain FA skeleton using TBSS (49,50) Khan et coll. also did probabilistic tractography.
In these two articles, no significant group differences in FA were found between patients suffering from burning mouth and controls. The remaining article used graph theoretical metrics to evaluate changes in the brain network at local and global connectivity using Brain Connectivity Toolbox (48,51). Wada et coll. showed higher connectivity at the bilateral rostral anterior cingulate cortex, right medial orbitofrontal cortex, and left pars orbitalis. Kurokawa et coll. showed higher betweenness centrality in the left insula, right amygdala, and right lateral orbito-frontal cortex and lower betweenness centrality in the right inferotemporal cortex. For both articles, no diffusion metrics were compared.
Chronic primary visceral pain
Chronic primary visceral pain is CPP localized in the head or neck, thoracic, abdominal, or pelvic region (29). The respective anatomical location is compatible with typical referral pain patterns from specific internal organs. Chronic primary visceral pain includes pain: in the head/neck viscera of the digestive system; in the thoracic region (e.g., noncardiac chest pain and reflux hypersensitivity); in the abdominal region arising from viscera of the digestive system (e.g., epigastric pain syndrome, irritable bowel syndrome, centrally mediated abdominal pain syndrome, and biliary dyskinesia) and in the pelvic region due to involvement of the viscera of the digestive, urinary, and genital systems (e.g., bladder pain syndrome, anal spasm, chronic pelvic pain, and chronic testicular pain).
For this category of pain, 19 studies were included. From these, seven studies focus only on irritable bowel syndrome (IBS) (53–59) (precisely, Nan et coll. focused on constipation-predominant IBS), three on only urological chronic pelvic pain syndrome (UCPPS) (60– 62), four only on primary dysmenorrhea [PDM] (63–66), two on chronic prostatitis/ pelvic pain syndrome [CP/CPPS] (67,68), one on provoked vestibulodynia (PVD) (69), one on interstitial cystitis/bladder pain syndrome (IC/BPS) (70) and one focused on two conditions: UCPSS and IBS (71).
Generally, all these papers have reported fractional anisotropy (FA) and mean diffusivity (MD) values. In studies with individuals suffering from IBS, the authors have, among other things, observed a higher FA value in posterior/anterior cingulate, corpus collosum, fornix and S1 regions when compared to control groups (53,54,56). However, the study of Hubbard and collaborators reported lower FA in the right dorsal cingulum bundle compared to controls and the study of Nan and collaborators reported lower FA in the corpus callosum (55,59). Although TBSS are used to extract metrics, the ROIs targeted are not the same between study. For example, a study looked at ROIs from white matter subserving areas associated with pain implicated in IBS including lateral (VPL) and medial (MD) thalamus, ACC (mid-anterior cingulum bundle), S1, insula and adjacent external capsule, and fornix (53) whereas another study focused on cingulum region essentially, from JHU atlas (55). Thus, the comparison between the different studies becomes difficult and hazardous.
For people with UCPPS, studies have shown that patients had lower FA and higher MD in brain regions commonly associated with perception and integration of pain information such as thalamus, basal ganglia and splenium of the corpus callosum (62,71). In contrast with the precedent studies, Huang and collaborators had observed that UCPPS patients exhibited higher FA in a portion of the left corticospinal track projecting through the precentral gyrus adjacent to the motor homunculus compared to control participants (61). Interestingly, a multi-site study has shown significant spatial variations in FA and the apparent diffusion coefficient (ADC, mathematically almost identical to MD) in statistical parametric maps for all healthy volunteers and UCPPS and positive control participants with adequate quality DTI when evaluated between sites (60). These observations thus highlight some of the consequences that multi-site studies can have for neuroimaging-based clinical phenotyping in UCPPS.
Moreover, compared to healthy controls, patients with PDM had increased FA along with decreased MD in the corpus callosum, superior longitudinal fasciculus, corona radiata, internal and external capsule (66). While on the other hand, some studies have shown a lower FA in superior/posterior corona radiata, external capsule, dorsal-posterior cingulum and a higher MD in these same regions in patients with PDM (63,65). Also, another study has observed a significant negative correlation between the FA values of fiber tracks of thalamus-SI pathways and pain intensity in women with PDM, suggesting that women with higher level of pain had alteration in white matter integrity (64).
Chronic prostatitis/pelvic pain syndrome is another type of visceral pain. A study from Farmer and collaborators showed that there was no difference in mean whole brain white matter FA between patients with CP/CPPS and control participants (67). Using a different methodological approach, a DTI study used Trackvis diffusion toolkit to investigate the topological organization and proprieties of white matter brain network in patients with CP/CPPS (72). Briefly, the authors observed that the architecture of white matter networks of CP/CPPS patients and healthy controls were both presenting small-world properties. However, they showed that the frontal-parietal control network was presenting altered connectivity in patients with CP/CPPS.
Regarding the other subtype of visceral pain, individuals with PVD showed significant increases in the FA of the somatosensory and basal ganglia areas as well as a higher MD in the basal ganglia (internal capsule and pallidum) when compared to healthy controls (69). For IC/BPS patients, a study observed a decreased FA in parts of the right anterior thalamic radiation, left forceps major, and right longitudinal fasciculus, and an increased FA was detected in the right superior and bilateral inferior longitudinal fasciculi (70). Besides, these brain white matter abnormalities are correlated with symptom severity associated with IC/BPS.
Chronic primary musculoskeletal pain
Chronic primary musculoskeletal pain is CPP located in the muscles, bones, joints, or tendons (29). A typical example is chronic primary low-back pain. Chronic primary musculoskeletal pain syndromes are distinguished according to location: upper (chronic primary cervical pain), middle (chronic primary thoracic pain), lower back (chronic primary low-back pain), and limbs (chronic primary limb pain).
For this category of pain, ten studies were included. From these, five were on sub-acute/chronic low back pain (SBP/CLBP) (8,73–76), three on chronic musculoskeletal pain syndrome (77–79), one on nonspecific low back pain (80) and one on chronic neck pain (81).
For the low back pain studies, two were on a dataset that investigates transition from subacute to chronic low back pain (8,73). However, they used different methods that are not comparable, one used TBSS and identified clusters with lower FA only in subjects that transitioned to develop chronic pain (one cluster in the temporal part of left superior longitudinal fasciculus, a second located in the left retro-lenticular part of the internal capsule, and a third one in the left anterior limb of the internal capsule and part of the corpus callosum including the anterior corona radiata). While the other used probabilistic tractography to construct a network that identified a higher density of connections between corticolimbic regions (nucleus accumbens, amygdala, hippocampus-parahippocampus, and the prefrontal cortex) that were present in subjects that transitioned to develop chronic pain.
For the three studies specifically on chronic low back pain (75,76,82), one investigated WM regions on a TBSS skeleton that are adjacent to regions representing lower back and finger sensory regions as predetermined by fMRI experiments and were able to observe lower FA in those regions for the CLBP patients (76). Another study used the same concept, but for fMRI seeds located in the insula (75). They did not find differences in FA in the WM ROI bordering the insula between HC and CLBP, but they showed that FA of this region increased after spine surgery or facet joint block treatment in CLBP patients. The last study (82) also used TBSS to generate whole brain FA skeleton, but then segmented the skeleton using 50 WM ROIs from the ICBM-DTI-81 white matter atlas and found reduced FA in regions such as the corpus callosum, anterior and posterior thalamic radiation, SLF and corona radiata.
For chronic musculoskeletal pain syndrome, two studies are from the same group and most likely on an overlapping dataset comprising a heterogenous chronic pain population of ∼ 50% of back pain, ∼20 % of osteoarthritis (OA) and ∼20% of post trauma patients. They either used TBSS only (78) or TBSS and a more complex multi-compartment Ball- and-Sticks modeling (77) and found regional differences in the splenium of the corpus callosum (lower FA, mode of anisotropy [MO] and Ball- and-Sticks partial volume fraction metrics) and the cingulum (lower FA and partial volume fraction metric). They also reported changes in both DTI based (increases in RD and/or AD in these cases) and Ball- and-Sticks (decreases in MO and/or partial volume fraction) metrics in internal and external capsules, cerebral peduncles and SLF. The third study on chronic musculoskeletal pain syndrome is one of the only using a recent modeling of the diffusion signal with a Neurite Orientation Dispersion and Density Imaging (NODDI) approach (79). They observed lower WM orientation dispersion index in few regions such as the anterior corona radiata, the right posterior thalamic radiation, the uncinate fasciculus and the fornix.
The participants from the nonspecific low back pain study (80) could potentially be classified as CLBP as they cite comparable studies that included CLBP participants. Nonetheless, this study is the only one using graph theoretical metrics to evaluate changes in the white matter architecture of low back pain patients and cannot be compared to any others. They found that local efficiency, a network metric that characterizes the ability of a network node to process specialized information at the local level, was lower in the low back pain patients when compared to healthy controls.
The last MSK study is about chronic neck pain (81) and compared neck pain from traumatic and idiopathic origin to healthy controls in a women-only cohort. After generating FA and diffusivities maps, they used the ICBM-DTI-81 WM label atlas to extract properties of several ROIs, they found decreased FA, and increased MD and RD in the left cingulum hippocampus and tapetum in traumatic compared to idiopathic neck pain and in the left tapetum in traumatic neck pain patients compared to controls.
Chronic widespread pain
Chronic widespread pain is diffuse musculoskeletal pain in at least 4 of 5 body regions and in at least 3 or more body quadrants (as defined by upper–lower/left–right side of the body) and axial skeleton (neck, back, chest, and abdomen) (29,83). Chronic widespread pain is characterized by the core features of CPP, such as pain persisting for at least 3 months, and associated with significant emotional distress and/or functional disability.
For this category of pain, six studies were included. All of which were made on fibromyalgia (FM) patients (84–89). Some of the earliest publications were found in this subgroup (84,89), where they were using the lowest number of gradient directions (6 and 9) and basic diffusion maps obtained from unspecified software or possibly in house pipelines. These two studies found similar differences in the thalamus (lower FA), and the Lutz study also found higher FA in several other grey matter regions such as the postcentral gyri, amygdalae, and hippocampi. One study that investigated diffusion maps extracted from the GE MRI software FuncTool did not report any differences in their diffusion measure (88) and another one, that investigated probabilistic tractography with the FSL tools, found more connection within the cerebellum in FM patients when compared to controls (86). One study compared young and old FM patients and report a lower FA in the corpus callosum only in older patients using FSL TBSS approach (85). The last study used diffusion maps extracted from ExploreDTI software to investigate the effect of hyperbaric oxygen therapy in FM patients and found an increase in FA in anterior and superior thalamic radiation, left Insula and right thalamus following the therapy (87).
Complex regional pain syndrome (CRPS)
Complex regional pain syndrome (CRPS) is a type of CPP characterized by pain in a regional distribution that usually starts distally in an extremity after trauma and that is disproportionate in magnitude or duration to the typical course (of pain) after similar tissue trauma (29,90,91). For this chronic pain condition, only two studies were found (92,93). Both studies used a whole brain TBSS approach. Hotta et coll. identified regions of the WM skeleton showing higher MD, AD and RD. They showed extensive clusters (between ∼300 and 1400 voxels) in the corpus callosum, internal capsule, corona radiata and SLF, and smaller ones (between ∼14 and 150 voxels) namely in the fornix, external capsule, posterior thalamic radiation and cingulum. As for Geha et coll., they showed lower FA in a cluster within the left callosal fiber track. In addition, they performed probabilistic tractography to trace and quantify pathways from these clusters in both groups. They observe that in the left hemisphere (ipsilateral to the seed), average number of connections at long distances (> 60 mm from seed), as well as total number of connections, were significantly lower in CRPS patients.
Discussion
The purpose of this review was to provide a critical summary of the use of brain dMRI for the study of primary chronic pain conditions. Each article was classified according to the latest IASP chronic pain definition, dMRI sequence and analysis method. The main findings of this review highlight the difficulty of delineating common white matter abnormalities for each chronic pain condition. This observation comes in part from the vast number of possibilities to analyze and further report results from dMRI data. Notably:
Twenty-three percent of studies (14/59) did not report diffusion metrics comparisons between CP and healthy controls (seven articles in the orofacial pain category, four articles in the visceral pain category, two articles in the musculoskeletal pain category and one article in the widespread pain category). Alternately, these articles used dMRI to develop prediction models, quantify graph theory properties, or correlated metrics with various behavioral and physiological parameters.
Fourteen percent of the studies (8/59) did not report any differences between the groups that were investigated. Of these, four were in the orofacial pain category (two migraine and two BMS studies), one in the musculoskeletal pain category (CLBP), two in the visceral pain category (one IBS and one CP/CPPS) and one on widespread pain (fibromyalgia).
Finally, six percent of the studies (4/59) reported differences of FA or ODI in the fornix (two in visceral pain, one in CRPS and one in MSK) even though, due to its unique location surrounded by CSF, it is most likely affected by partial volume effects (PVE) even when using state-of-the art dMRI acquisition sequences (94). Furthermore, as presented later, when using a TBSS approach, the fornix is almost absent from the TBSS FA skeleton.
Main critics of the approaches used in the 59 articles
Acquisition parameters
One of the main issues explaining the lack of dMRI reproducibility is the large variability of acquisition parameters. For example, a recent study reported “Multisite, multiscanner, and multisubject acquisitions for studying variability in diffusion weighted MRI “ and showed that with certain acquisition parameters, interscanner (same subject) variability could approach intersubject (same scanner) variability (95). As for any MRI sequences, many parameters need to be properly chosen to ensure quality diffusion sensitized images. Three main parameters which drastically influence quality and analysis outcome are discussed. 1) the b value; representing the strength, duration and timing of the diffusion encoding gradients (96,97); 2) the number of diffusions encoding gradient directions; generally representing the number of diffusion gradient directions applied over a sphere (98,99); and 3) the voxel dimensions; representing the length, width and height of the 3D image voxels, (100). The most common sequences developed for dMRI, suggested for DTI had: a b value of 1000s/mm2, 30 unique gradient directions and a 2mm isotropic resolution.
Across studies, the b-value parameter was the most stable with ∼75% of studies using a b value of 1000. Still, some of the b-values reported were either in an unusual range (700, 800, 900, 1200 and 1300) without any explanation or were not even reported. A striking observation is that only two groups performed multi-shell acquisitions, an acquisition strategy which emerged over 15 years ago (101) that allows for state-of-the-art dMRI analysis and that with modern multiband sequences is easily doable in an acceptable amount of time (between 10 and 20 minutes at most).
When investigating the number of diffusions encoding gradient directions, we found that ∼55% of the studies acquired a minimum of 30 directions (with a maximum of 99 directions). However, this parameter was also highly variable with some studies acquiring only six or nine directions, which is troublesome as the theoretical limit to properly describe a diffusion tensor is six gradient directions (102,103). In the context of tractography (and more generally fiber orientation estimation), high angular resolution diffusion imaging (HARDI) was introduced as an efficient acquisition strategy to overcome DTI limitations. Based on the acquisition of over 50 gradient directions at a single high b-value, HARDI (paired with multi-shell acquisition) is also doable in an acceptable amount of time (between 10 and 20 minutes at most).
Lastly, when investigating image resolution, the overall voxel volume ranged from 2.4mm3 to 20mm3. Forty-seven percent of studies used the conventional 2mm isotropic resolution, however, substantial variation was also noted, with some studies acquiring highly anisotropic voxels (e.g. 1.875×1.873×3 or 2×2.5×4) which can bias dMRI analysis. For example, the bigger the voxel, the more likely the voxel will contain multiple fiber populations (i.e. decrease fiber orientation homogeneity) thus impacting diffusion metrics (100) and anisotropic voxels will impact tractography algorithm, especially in branching situation (104).
Overall, the observations made in this section highlight the need for more standardized acquisition parameters and pre-processing pipelines across chronic pain dMRI studies. In addition to DTI MRI sequence recommendations and following human connectome project (HCP) guidelines (105) we recommend that new studies acquire multi-shell HARDI data at over 2mm isotropic resolution and a b=0 reverse phase-encoded image to correct susceptibility image distortions. We think that recommending this type of acquisition scheme will leave most space for protocol evolution and optimization while promoting reproducibility and replicability across datasets and analysis methods. For a recent review of dMRI preprocessing, we refer to this article (106).
Processing and analysis tools
The tools provided by the FMRIB group at Oxford were the most used, with over 70% of the papers using one of the FSL tools to process their dMRI data. Almost 50% of the papers use the Tract Based Spatial Statistic (TBSS) pipeline (107) to identify voxel-wise differences. Although TBSS is a valid approach that provides significant advantages over classic whole-brain voxel-based analysis when it comes to group comparison, this method does not use the full potential of dMRI images as it summarizes the complexity of the whole-brain white matter to a few voxels wide WM skeleton. Furthermore, this approach does not completely exclude registration errors. Therefore, almost no information coming from tracks spanning up to the cortex or tracks that are smaller or located in complex regions can be found in the TBSS maps. The other most commonly used softwares to process dMRI data were ExploreDTI (108) and MRtrix (109), both offering advanced tools specifically designed for dMRI data to generate tractography and diffusion metrics maps. To circumvent issues brought by TBSS, some studies presented here have used a clever approach to identify white matter tracts impacted by the chronic pain condition under study. They used the clusters of significant differences identified in the TBSS results as seeds to perform probabilistic tractography. Although this approach allows for the identification of actual white matter tracts that were not present in the WM skeleton, it still cannot reconcile the fact that a significant volume of white matter was not included in the original TBSS analysis. Similarly, several groups have identified GM ROI from fMRI experiments and then expanded this region to include adjacent white matter as a WM ROI to extract diffusion properties. However, this needs to be taken carefully as a WM bundle passing close to a GM region does not necessarily connect with that region. Indeed, many WM tracts travel long distances in the brain without connecting with each region they are bordering along the way. A more appropriate approach could have been to use the GM ROI to generate seeds from/to which WM fibers might connect and extract those tracks from a whole brain tractogram.
To illustrate the constraints of TBSS methods, we examined the disparities between TBSS and a track-based approach for two specific tracks, namely the fornix and the accumbofrontal (AcF) track, which are of particular significance in chronic pain research owing to the regions they connect. Notably, the fornix serves as the primary pathway for efferent signals from the hippocampus, and it plays a critical role in memory circuitry (94). Moreover, hippocampus volume has been shown to be a risk factor in the transition from acute to chronic pain (110). The AcF track connects the orbitofrontal cortex to the nucleus accumbens (111), both regions that are also implicated in the transition from acute to chronic pain (7). To demonstrate the benefits of employing targeted approaches for investigating white matter tracks involved in chronic pain, we analyzed the overlap of two streamline bundles extracted using separate techniques with a commonly used whole brain FA skeleton. For the first technique, data from a single subject from Mansour et coll. paper (73), that can be found on the OpenPain neuroimaging dataset platform (https://openpain.org) was used to generate a track segmentation of the left AcF track in native imaging space. Firstly, the anatomical (T1w structural image, 1 mm isotropic) and diffusion data (8 b0s, 60 b1000, whole brain 2 mm isotropic) were processed using tractoflow version 2.3.0 (112) outputting a particle filtered whole brain tractogram (113). Next the brainnetome atlas was registered to the subject space using the T1w image and a non-linear registration using the ANTS toolbox (114) followed by a nearest neighbor interpolation of the Brainnetome atlas (115) to the diffusion data. Streamlines for the left AcF were segmented by selecting all streamlines that joined the brainnetome ROIs (left nucleus accumbens and the left orbital frontal gyrus). The AcF track was then registered to MNI space using the previously calculated non-linear transformation. For comparison of the FA skeleton to multiple methods, the population averaged right fornix used in the RecobundlesX analysis pipeline (116) along with the aforementioned left AcF track were overlaid on the binarized whole brain FA skeleton. As a quantitative measure of overlap, voxels intersected by these tracks were extracted and the percentage of voxels overlapping the binarized FA skeleton were output for each track. The fornix and AcF tracks are displayed in Figure 2 along with the percentage of overlapping voxels showing 15% voxels overlap between TBSS and tractography approached for the fornix (figure 2 left), and only 12% voxels overlap for the AcF track (figure 2 right). This analysis was reproduced for a few other subjects from the OpanPain database and similar overlaps were obtained (not shown).
Where should the field go?
Similar to the challenge faced by fMRI for processing and analysis (103,104), dMRI is facing reproducibility and replication issues. For example, when observing results of this review, one is faced with the relatively large variation in published results, measures and methods across studies. Typically, three main parameters can be used to improve group comparisons and improve study replicability: i) Increasing group differences, ii) Increasing study sample size and, iii) decreasing intra-group variability.
Increasing group differences: To increase group differences, novel biomarkers that increase sensitivity and specificity must be found. Recent advances in MR hardware, acquisition schemes (HARDI, DSI, multi-shell) and analysis method allow the investigation of new brain markers (117). However, the CP field has yet to yield the full potential of these methods. A striking example is DTI. The DTI method has taken so much space in the field of dMRI of white matter that the two are often used interchangeably in literature. Furthermore, the latest studies about dMRI do not use DTI as it cannot accurately represent more than one fiber population per white matter voxel. Accordingly, the tensor model usually fails to capture proper microstructure in voxels where crossing fiber occurs and because this happens in about 60 to 90% of all white matter voxels in the brain (118) the DTI method is now seen as providing unanatomical sound information. There are several new methods that exists, each attempting to account for a better approximation of the underlying white matter tissue organization. Notably in the articles from this review, four studies used whole brain connectome metrics and two studies used multi-compartment local models.
Increasing study sample size: Chronic pain lacks a large dMRI imaging open dataset. Recent improvements in metadata sharing, unthresholded statistical maps sharing, standardized dataset structures, multi-site harmonization, and standardized metrics could allow the rapid expansion of open-access datasets to study chronic pain.
Decreasing intra-group variability: Both chronic pain and dMRI are subject to large intra-group variability. For chronic pain, intra-group variability can be mitigated by reducing chronic pain condition confounds (improving clinical evaluation specificity). Concerning dMRI, intra-group variability can be mitigated at each step of the study from acquisition to analysis. We would like to emphasize that, in the context of group comparisons, dMRI data presents additional challenges over other structural MR data such as region of interest positioning, registration and smoothing. Notably, due to these common dMRI pitfalls (119), application of commonly used structural MR analysis methods, such as whole brain voxel-based analysis, are discouraged. Therefore, we recommend that the acquisition parameters, pre-processing, processing and analysis be conducted with standardized pipelines specific to dMRI.
Overall, the goal of this review is not to favorize one over the other but to encourage the field to move towards more up to date modelling of dMRI signal and using methods along High Angular Resolution Diffusion Imaging (HARDI) concepts. This procedural update will allow pain researchers to capture subtle white matter abnormalities that cannot be observed with more conventional methods.
A note on acute pain
The link between white matter characteristics and acute pain perception, as often reported by pain thermal stimulation fMRI studies, is under-studied. Preliminary screening revealed that only 13 studies specifically investigated this aspect. Therefore, due to significant methodological differences compared to chronic pain, acute pain studies were excluded from this review.
Limitations
Due to evolving chronic pain defintions and dMRI nomenclature, it is likely that articles performing dMRI on primary chronic pain patients were missed in this review. For example, an article from Geha et al (93), referenced by other studies in this review, was erroneously excluded because its dMRI nomenclature did not meet our inclusion criteria. Moreover, the most recent chronic pain definition was used to select and classify the literature to be reviewed (29), therefore, some papers might have identified their participants otherwise would they have used these new definitions.
Data Availability
All data produced in the present work are contained in the manuscript
Author Approval
All authors have seen and approved the manuscript.
Conflits of interest
MD is co-founder of IMEKA inc.
All other authors declare no conflict of interest.
Funding
MS is supported by a PhD scholarship from the CIHR. GLi is supported by FRQS UNIQUE postdoctoral scholarship. PT is supported by FRQS J1 salary award and Arthritis Society star career development award. GLe is supported by FRQS J2 salary award.
Footnotes
↵# co-first authors
References
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.
- 10.↵
- 11.↵
- 12.
- 13.
- 14.
- 15.
- 16.↵
- 17.↵
- 18.↵
- 19.
- 20.
- 21.
- 22.
- 23.↵
- 24.↵
- 25.
- 26.
- 27.
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.
- 58.
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵