Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 09 July 2021
Sec. Experimental Therapeutics
This article is part of the Research Topic Current State and Future Directions of Cranial Focused Ultrasound Therapy View all 10 articles

Peripheral Nerve Focused Ultrasound Lesioning—Visualization and Assessment Using Diffusion Weighted Imaging

\nMatthew R. Walker,Matthew R. Walker1,2Jidan ZhongJidan Zhong2Adam C. Waspe,Adam C. Waspe3,4Karolina PiorkowskaKarolina Piorkowska3Lananh N. NguyenLananh N. Nguyen5Dimitri J. Anastakis,,Dimitri J. Anastakis1,2,6James M. Drake,,James M. Drake1,3,7Mojgan Hodaie,,
Mojgan Hodaie1,2,8*
  • 1Institute of Medical Science, University of Toronto, Toronto, ON, Canada
  • 2Division of Brain, Imaging & Behaviour, Krembil Research Institute, University Health Network, Toronto, ON, Canada
  • 3Centre for Image Guided Innovation and Therapeutic Intervention, Hospital for Sick Children, Toronto, ON, Canada
  • 4Department of Medical Imaging, University of Toronto, Toronto, ON, Canada
  • 5Laboratory Medicine Program, University Health Network and University of Toronto, Toronto, ON, Canada
  • 6Department of Surgery, Toronto Western Hospital, University Health Network and University of Toronto, Toronto, ON, Canada
  • 7Department of Neurosurgery, Hospital for Sick Children, Toronto, ON, Canada
  • 8Department of Neurosurgery, Toronto Western Hospital, University Health Network, Toronto, ON, Canada

Objectives: Magnetic resonance-guided focused ultrasound (MRgFUS) is a non-invasive targeted tissue ablation technique that can be applied to the nervous system. Diffusion weighted imaging (DWI) can visualize and evaluate nervous system microstructure. Tractography algorithms can reconstruct fiber bundles which can be used for treatment navigation and diffusion tensor imaging (DTI) metrics permit the quantitative assessment of nerve microstructure in vivo. There is a need for imaging tools to aid in the visualization and quantitative assessment of treatment-related nerve changes in MRgFUS. We present a method of peripheral nerve tract reconstruction and use DTI metrics to evaluate the MRgFUS treatment effect.

Materials and Methods: MRgFUS was applied bilaterally to the sciatic nerves in 6 piglets (12 nerves total). T1-weighted and diffusion images were acquired before and after treatment. Tensor-based and constrained spherical deconvolution (CSD) tractography algorithms were used to reconstruct the nerves. DTI metrics of fractional anisotropy (FA), and mean (MD), axial (AD), and radial diffusivities (RD) were measured to assess acute (<1–2 h) treatment effects. Temperature was measured in vivo via MR thermometry. Histological data was collected for lesion assessment.

Results: The sciatic nerves were successfully reconstructed in all subjects. Tract disruption was observed after treatment using both CSD and tensor models. DTI metrics in the targeted nerve segments showed significantly decreased FA and increased MD, AD, and RD. Transducer output power was positively correlated with lesion volume and temperature and negatively correlated with MD, AD, and RD. No correlations were observed between FA and other measured parameters.

Conclusions: DWI and tractography are effective tools for visualizing peripheral nerve segments for targeting in non-invasive surgical methods and for assessing the microstructural changes that occur following MRgFUS treatment.

Introduction

Magnetic resonance-guided focused ultrasound (MRgFUS) is a technique to thermally ablate targeted tissue using MR imaging for navigation (1). It is non-invasive and does not involve ionizing radiation. Current clinically approved indications include essential tremor (2), prostate cancer (3), uterine fibroids (4), and bone metastases (5). There is great potential in extending the use of MRgFUS to the peripheral nervous system for treating conditions such as spasticity and chronic pain (6), which is supported by observations of FUS effects on nerve conduction (7). Ablation techniques, of which MRgFUS is a potential alternative, have been studied in peripheral nerve-related conditions including painful stump neuromas (8), peripheral nerve sheath tumors (9), inguinal neuralgia (10), and lumbar degenerative disease (11). Ablation has also been used for cancer pain relief (12, 13) and to treat tumors with proximity to neural structures (14). An understanding of the nerve involvement in these ablation procedures is of great concern whether the nerves are to be targeted directly or to be avoided.

Diffusion weighted imaging (DWI), an imaging technique sensitive to the movement of water in tissue, has the ability to image nerve fibers and measure their microstructural characteristics in vivo (15). Tractography algorithms permit the reconstruction of nervous system fiber connectivity based on patterns of restricted water diffusion (16). Modeling the movement of water with tensors, called diffusion tensor imaging (DTI), includes metrics for quantitative evaluation of nerve characteristics. Fractional anisotropy (FA) reflects the directionality of diffusion and has been used as a proxy for nerve integrity (16). Mean diffusivity (MD), an average of all three orthogonal tensor indices, describes the overall magnitude of diffusion and reflects the degree of water diffusion restriction within tissues regardless of fiber orientation (17). Radial diffusivity (RD), a measure of diffusion perpendicular to the primary diffusion direction, is associated with degree of myelination (18). Axial diffusivity (AD), a measure of diffusion parallel to the primary diffusion direction, is sensitive to axonal integrity (18). These DTI metrics can be used as quantitative tools to assess peripheral nerve pathology and measure microstructural changes following treatment.

Accurate imaging is imperative in non-invasive treatments such as MRgFUS as faulty targeting could exacerbate comorbidities or depress treatment outcomes. It is also important that the imaging characteristics are well-understood in order to properly assess treatment response and inform treatment parameters. Conventional MR sequences, such as T1- and T2-weighted imaging, are limited as they cannot selectively visualize peripheral nerves or quantify nerve integrity or injury. Peripheral nerves have been visualized using MR neurography (19) and selective excitation techniques (20), however these approaches do not provide quantitative assessment of nerve changes after treatment. Targeting by atlas or structural images only is also limited in the ability to account for subject variability and specificity in identifying tracts of interest (21). Tractography can remedy these limitations by enhancing navigation in treatments such as transcranial thalamotomy for essential tremor (22, 23) and deep brain stimulation (21, 24, 25), resulting in improved targeting accuracy and patient outcomes.

Early work in DTI suggests its potential to evaluate peripheral nerve injury and regeneration in vivo (26, 27). Diffusion imaging permits longitudinal assessment and, in the case of animal studies, obviates the need for large subject numbers to be sacrificed at multiple time points, providing insight into cellular changes in lieu of histological data. Imaging of peripheral nerves has associated technical challenges due to the complexity and variability of the peripheral nervous system and surrounding muscle tissue, which is also fibrous and thus carries an anisotropic diffusion signature. Previous work has shown the importance of employing appropriate diffusion processing and tractography techniques in order to achieve anatomically accurate results (24).

There is a paucity of MRgFUS lesioning studies focusing on peripheral nerves and a greater limitation on the use of DWI to guide treatment. This is primarily due to the technical limitations associated with the accurate anatomical identification of these nerves. In this study, we investigate the use of diffusion tractography for targeting the sciatic nerves of piglets in ablative MRgFUS and DTI metrics for assessing the microstructural changes following treatment. Histological lesion analysis provides insight into cellular changes after treatment and their correlations with imaging and treatment parameters.

Materials and Methods

Animal Model

These experiments were approved by the Animal Care Committee and Laboratory Animal Services at the Hospital for Sick Children in Toronto, Ontario, Canada. This study conforms to the policies of the Canadian Council on Animal Care (CCAC).

Six male Yorkshire piglets (average weight 6.7 ± 1.3 kg, age 24 ± 4 days) were used. The animals were pre-anesthetized with ketamine solution [10 mg/kg] (Ketalean, CDMV Inc., Quebec, Canada) intramuscularly before being intubated and anesthesia maintained with 2.5% isoflurane and 2 L oxygen via MR-compatible ventilator. Hair was removed from both thighs of the animal via shaving and commercial depilatory cream to facilitate skin surface coupling with the MRgFUS system, aided by degassed ultrasound gel.

Once prepared, the animals were transported to the MR facilities for pre-treatment imaging on a standard clinical diagnostic table. Heart rate, peripheral capillary oxygen saturation, and body temperature were monitored throughout the experiment. A circulating water blanket was used to help maintain the piglet's core body temperature around 37°C. Upon completion of the experiment, the animals were euthanized while under anesthesia via intravenous injection of sodium pentobarbital [120 mg/kg] (Euthanyl, CDMV Inc., Quebec, Canada).

MR Imaging

DWI and T1-weighted images were acquired before and after treatment using a clinical Philips Achieva 3T MR scanner (Philips Healthcare, Best, Netherlands) and a 32-channel receive-only cardiac coil on a standard diagnostic table. For both imaging and treatment, the animals were placed in lateral decubitus position with thighs perpendicular to the long axis of the body and the leg of interest closest to the tabletop. Separate data sets were acquired for each leg with the lateromedial coverage extending from the outer skin surface to the contralateral spinal nerve roots. This position provided stability to limit potential movement of the animal and a clear path for the ultrasound beam to target the sciatic nerve.

Anatomical images were acquired with a three-dimensional T1 magnetization prepared rapid gradient echo (MPRAGE) sequence. Acquisition parameters included: repetition time (TR) 8.1 ms; echo time (TE) 3.7 ms; flip angle 8°; matrix 224 × 224; field of view (FOV) 224 × 224 mm; slice thickness 1 mm; slice number 70; voxel resolution 1 × 1 × 1 mm; number of signal averages (NSA) 4; SENSE reduction factor 2; acquisition time 14 min 35 s.

Diffusion images were collected with a SENSE-single shot spin-echo echo-planar-imaging (SE-EPI) sequence with a b-value of 800 s/mm2 and 128 diffusion encoding directions. Two additional baseline images with b = 0 s/mm2 were acquired, one each in forward and reverse phase-encode directions, for post-processing EPI-based susceptibility distortion corrections. Other diffusion scanning parameters include: TR 5845 ms; TE 106 ms; flip angle 90°; matrix 128 × 128; FOV 205 × 205 mm; slice thickness 1.6 mm; slice number 38; voxel resolution 1.6 × 1.6 × 1.6 mm; NSA 2; SENSE reduction factor 2; diffusion gradient pulse duration/time interval 15.7/52.9 ms; acquisition time 29 min 33 s.

Image Post-processing and Tractography

Post-processing was carried out with the FSL software library (Analysis Group, FMRIB, Oxford, UK: https://fsl.fmrib.ox.ac.uk/fsl/fslwiki) (28). Corrections were performed to remedy distortions caused by EPI and susceptibility-induced off-resonance fields using the two baseline images with opposing phase-encoding as implemented in FSL “topup” (29, 30). Susceptibility-corrected images were further processed to remove distortions associated with bulk motion and eddy currents in FSL “eddy” (31). Water diffusion was modeled from fully corrected data by fitting tensors to each image voxel (32). Tensors were used to calculate scalar maps of FA, MD, AD, and RD (33). Structural and diffusion images were linearly co-registered using FSL FLIRT and manual adjustment (34). Pre- and post-treatment images within each subject were co-registered with the same method. Image processing time required was ~40 min for each subject using a moderately powerful workstation with dedicated graphical processing unit.

Fiber tracking was performed with MRtrix (Brain Research Institute, Melbourne, Australia: http://www.brain.org.au/software) (35). The response function for a single fiber population was estimated using the default threshold of FA > 0.2 (36). This response function was then used with a basis of constrained spherical deconvolution (CSD) to estimate the fiber orientation distribution (FOD) (37). CSD has been shown to be an effective tractography method in regions of complex fiber orientations and crossing fibers (38). A deterministic tractography algorithm, “SD-Stream,” was used to generate tracks (35). Tracking seeds were delineated manually at the lumbar nerve roots based on structural T1 and FOD maps. Tracts were segmented by placing inclusion regions of interest (ROI) in terminating muscle regions.

DTI metric assessment was performed by manual placement of a 2 × 2 × 1 voxel (3.2 × 3.2 × 1.6 mm3) ROI on the sciatic nerve overlapped by the post-treatment lesion as identified by tractography and T1 data. Measurements of FA, AD, RD, and MD were pulled from co-registered pre- and post-treatment images using the same ROI mask. Visually distinct lesion zones were manually identified and measured on T1, aided by voxel intensity thresholding, as a region of hyperintensity (zone I) concentrically surrounded by hypointensity (zone II).

MRgFUS Treatment

Treatment was performed with a clinical MRgFUS system (Sonalleve V1, Profound Medical, Toronto, Canada). Animals were positioned on the treatment table housing the ultrasound transducer. The same lateral decubitus position was used in pre- and post-treatment imaging. Degassed ultrasound gel was applied to the skin surface with a 20 mm gel pad (Aquaflex, Parker Laboratories, New Jersey, USA) placed between the animal and treatment table to facilitate acoustic coupling with the transducer.

Imaging during treatment was accomplished using one element of the cardiac coil (16 channels) placed on top of the animal. Skin bubble images were acquired using a three-dimensional spoiled gradient echo (FFE) sequence to confirm that no air bubbles were present which may interfere with acoustic beam propagation. T1-weighted FFE images were acquired for immediate target identification and treatment cell placement using the Sonalleve planning software.

Cells were placed on the sciatic nerve, posterior to the proximal head of the adjacent femur, with guidance from pre-treatment tractography (Figure 1A). The sciatic nerves in these piglets are ~3 mm wide therefore treatment cell diameters of 4 and 8 mm were chosen to cover the whole nerve. Sonication times are fixed based on cell diameter leading to 20 and 27 s treatments for these respective cell sizes. The distance along the beam path from skin surface to treatment cell center was measured. A single sonication was used for each treatment in order to isolate the effect of a given cell size and transducer output power. All ultrasound treatments were delivered with a frequency of 1.2 MHz. Before each therapeutic exposure, one to two test sonications at reduced power of 10 W and 20 s were performed for calibration of beam focus location. Temperature was measured simultaneously with sonication via MR thermometry (Figure 1B). Specific treatment parameters for each nerve can be found in Table 1. Immediately following treatment in both legs, the animal was repositioned on the diagnostic table for post-treatment imaging (1–2 h after sonication). Identical scanning sequences were used as in pre-treatment imaging.

FIGURE 1
www.frontiersin.org

Figure 1. Pre-treatment planning images and in-treatment MR thermometry in axial, sagittal, and coronal view. (A) T1-weighted images are shown with the MRgFUS treatment cell placed on the sciatic nerve. (B) In vivo MR thermometry displaying heat deposition during treatment.

TABLE 1
www.frontiersin.org

Table 1. MR-guided focused ultrasound treatment parameters.

Histology

Animals were euthanized immediately after post-treatment imaging (<2 h following treatment). The treated area was identified by reflecting the biceps femoris muscle to expose the sciatic nerve and lesion on surrounding muscle tissue. The treated portion of the sciatic nerve and adjacent sections of biceps femoris and semitendinosus muscles were collected. Treated muscle specimens were sampled superficially, deep or mirror to the treated area. Nerve specimens were sampled longitudinally along the axis of the beam path. Internal control specimens of both nerve and muscle were sampled several centimeters away from treated areas, confirmed to free from temperature changes via MR thermometry maps and tissue damage via gross inspection. Samples were immersion fixed with 10% neutral buffered formalin and cooled in a 4°C refrigerator for 48 h before routine histological processing. The tissue was sectioned at a 5-micron thickness at 200-micron levels and stained with hematoxylin and eosin (H&E) and luxol fast blue (LFB). Control samples were exposed to the same conditions and processing as for treated samples except for MRgFUS ablation. Thus, any differences observed between treated and control tissue would be related to MRgFUS treatment. All specimens were evaluated by a neuropathologist who was blinded to the treatment conditions of each specimen.

Statistical Analysis

DTI metric comparisons from pre- and post-treatment lesion ROIs were carried out using paired, two-tailed t-tests. Significance was taken as P < 0.05. Pearson correlations were performed between maximum temperature, temperature change from baseline, lesion volumes (zones I and II), DTI metrics, and output power (absolute and normalized by skin-to-focus distance). Statistical analysis was performed using SPSS version 23 (IBM, Inc.).

Results

Peripheral Nerve Imaging

On T1-weighted imaging, the sciatic nerve presented as a moderate hyperintensity associated with suppressed signal from nearby blood vessels. Surrounding fat and connective tissue of similar signal intensity to the nerve introduced uncertainty in assessing nerve position at all points in its caudal trajectory from the lumbar plexus. However, the sciatic nerves were successfully reconstructed via both single tensor and CSD tractography in all subjects (Figure 2). Tract models were observed extending from the vertebral roots, through the lumbar plexus, and terminating in the muscles of the leg posterior to the femur.

FIGURE 2
www.frontiersin.org

Figure 2. Tractography before and after MRgFUS treatment in sagittal view. Tracts are reconstructed using (A,B) constrained spherical deconvolution (CSD) and (C,D) diffusion tensor imaging (DTI) methods. Tracts are overlaid on T1-weighted images and colored by orientation: anterior-posterior (blue); superior-inferior (red); left-right (green). Arrowhead indicates treatment target.

MRgFUS Treatment Assessment

Tract abnormalities following MRgFUS treatment were observed in both single tensor and CSD reconstruction models (Figure 2). For the tensor model, tracts were discontinuous within the lesion in 8 of the 12 nerves imaged. DTI tracts that did extend into the lesion displayed narrowing and decreased fiber density. In the CSD model, tract continuity was maintained through the lesion in all cases with narrowing and decreased fiber density within the treatment area.

Lesions were identified on T1 images by regions of hyperintensity (zone I) surrounded by hypointensity (zone II) (Figure 3). Zone I ranged 20–531 mm3 with an average of volume of 207 mm3. Zone II ranged 81 to 1,270 mm3 with an average of 492 mm3. Total volume ranged 101–1,700 mm3 with an average of 699 mm3. Individual volume measurements and treatment parameters are listed in Table 1.

FIGURE 3
www.frontiersin.org

Figure 3. Post-treatment T1-weighted images in axial, sagittal, and coronal view. (A) The lesion on the sciatic nerve shows distinct zones of hyperintense necrotic core and surrounding hypointense edema. (B) Lesion voxels were segmented into central zone I (red) and surrounding zone II (blue).

Significant correlations were observed between output power and maximum lesion temperature (correlation coefficient r = 0.81, P = 0.001), temperature difference (r = 0.86, P = 0.0004), lesion zone I volume (r = 0.82, P = 0.001), zone II volume (r = 0.86, P = 0.0003), and total lesion volume (r = 0.88, P = 0.0001). Output energy was also correlated with maximum temperature (r = 0.64, P = 0.025), temperature difference (r = 0.64, P = 0.025).

DTI metrics measured over the sciatic nerve within the lesion immediately after treatment compared to the pre-treatment baseline revealed significantly decreased FA (P = 0.00008), and increased MD (P = 0.0008), AD (P = 0.01), and RD (P = 0.0001). Results are depicted graphically in Figure 4.

FIGURE 4
www.frontiersin.org

Figure 4. Diffusion tensor metric assessment of the treated sciatic nerve segment before and after MRgFUS. Results show significant changes in fractional anisotropy (FA) and mean (MD), axial (AD), and radial diffusivity (RD). **P < 0.01; ***P < 0.001.

Significant negative correlations were also observed for output power with diffusivity changes relative to baseline: MD (r = −0.77, P = 0.004), AD (r = −0.70, P = 0.011), and RD (r = −0.70, P = 0.011). No correlations were seen between power and FA.

Maximum temperature and temperature difference were well-correlated with total lesion volume (r = 0.81, P = 0.004 and r = 0.80, P = 0.005, respectively) and zone II volume (r = 0.88, P = 0.001 and r = 0.89, P = 0.001, respectively) and relatively weakly correlated with zone I volume (r = 0.65, P = 0.023 and r = 0.67, P = 0.017, respectively). Temperature difference correlated with MD (r = −0.68, P = 0.016) and AD (r = −0.67, P = 0.016) but not significantly with RD (r = −0.56, P = 0.058) or FA (r = 0.05, P = 0.87).

Output power normalized by skin-to-focus distance (measured as distance from skin surface to center of treatment cell) was weakly correlated with only maximum temperature (r = 0.58, P = 0.047) and temperature difference (r = 0.68, P = 0.015)—less significant than correlations with absolute output power. Selected correlations are shown graphically in Figure 5.

FIGURE 5
www.frontiersin.org

Figure 5. Relationships between single sonication MRgFUS output and treatment results. Correlations are shown between transducer output power and resulting maximum temperature, temperature change from baseline, percent change in fractional anisotropy (FA), mean (MD), axial (AD), and radial diffusivity (RD), and lesion volume (zone I, II, and total).

Histological Analysis

Gross examination of the dissected tissue revealed pallor of the surrounding muscles and red discoloration of the perineural tissue, further confirming accurate and sufficient delivery of the acoustic energy to the sciatic nerve. Control specimens appeared unremarkable (i.e., no signs tissue damage observed). Microscopic examination of the muscle specimens showed pathologic changes in all areas sampled, superficially, deep or mirror, indicating sensitivity to the acoustic energy regardless of location. All specimens showed changes in a zonal or gradient pattern. Figures 6A–D shows a representative zonal pattern where on one end, minimal endomysial edema was present while on the other end, extensive edema and myofiber dropout were present. Nerve specimens (Figures 6E–H) showed marked changes. On H&E, specimens showed extensive perineurial and endoneurial edema. LFB stain highlighted loss of myelin while H&E/LFB dual stain highlighted axonal loss. Control specimens for muscle and nerve were unremarkable histologically.

FIGURE 6
www.frontiersin.org

Figure 6. Histologic specimens taken of muscle (A–D) and sciatic nerve (E–H) following single sonication MRgFUS treatment. (A) Low magnification showing zonal changes within the muscle—severe, moderate, and mild. Within this muscle specimen, colored boxes indicate areas of high magnification shown in panels (B–D), each corresponding to these zonal changes. (B) High magnification of severe muscle involvement (blue) shows diffuse edema, myofiber vacuolation, and myofiber dropout. (C) High magnification of moderate muscle involvement (green) shows perimysial and endomysial edema, myofiber vacuolation, and myofiber necrosis. (D) High magnification of mild muscle involvement (black) shows mild myofiber atrophy and endomysial edema. (E) High magnification of involved nerve shows endoneurial edema, demyelination, and axonal fragmentation as highlighted by (F) LFB stain. (G) High magnification of control nerve H&E stain and (H) LFB stain. H&E, hematoxylin and eosin; LFB, luxol fast blue. Scale bar is 2,000 μm in (A) and 100 μm for (B–H).

Discussion

In this study, we have demonstrated the feasibility of using diffusion MR tractography to identify and visually reconstruct peripheral nerves and guide their ablative treatment with a clinical MRgFUS system. Further, we have shown the ability of DTI metrics to quantitatively assess the acute changes following treatment. MRgFUS informed by tractography is capable of producing thermal lesions focused on peripheral nerves with minimal damage to surrounding tissue. DTI metrics demonstrate significant microstructural changes to nervous tissue in the form of decreased FA and increased MD, AD, and RD, the latter three of which were found to negatively correlate with transducer output power. Histological analysis verified the damage to the nerves and sharp transition zones from lesions and adjacent untreated tissue.

Experiments of constriction injury in rabbit sciatic nerves have reported microstructural alterations consistent with decreased FA and increased MD and RD (39, 40). However, we observed increased AD while these authors cited significant (40) and non-significant (39) decreases in AD. In constriction injury, it appears that the chronic compression primarily restricts water motion parallel with the nerve and secondarily drives an inflammatory response of axonal swelling and loosening of the myelin sheath. In acute MRgFUS, the thermal ablation results in axonal fragmentation and the accumulation of cell debris which, along with the influx of inflammatory fluid and loosening myelin, contributes to a change toward isotropic diffusion within the nerve. RD, which is associated with degree of fiber myelination, experienced the largest change relative to pre-treatment baseline. Thus, we suspect the myelin disruption and widening of the periaxonal space as seen on histology to be the driving factor behind the results of increased diffusivity and decreased directionality. This is supported by a previous report (7) which identified axonotmesis, where the axons and myelin sheath are damaged but the gross structure (i.e., epineurium, perineurium, and endoneurium) remains intact, in rat sciatic nerves following FUS exposure. Our tractography results corroborate this further as fiber reconstruction with CSD within the lesion was still possible after treatment. Narrowing of reconstructed fibers and decreased fiber density within the lesion is indicative of some axonal damage and decreased diffusion directionality.

In addition to significant changes in DTI metrics at the group level, negative correlations were observed for MD, AD, and RD with transducer output power. That is, with increasing power more modest alterations to diffusivities were measured. Similar correlations were seen between temperature change from baseline and MD and AD with the RD correlation failing to reach statistical significance (P = 0.058). FA changes did not correlate with any other observable. We interpret these FA findings as a possible ceiling effect to the FA decrease in which the nervous tissue has become sufficiently disordered such that further acoustic energy will not further affect diffusion directionality within the nerve. Negative correlations between the diffusivities and both output power and temperature were unexpected, however. We hypothesize that the damage associated with higher sonication intensity and lesion temperature results in increasing damage and accumulation of cellular debris which acts to inhibit the increased motion of water within the nerve. RD exhibited greater relative change with power than both MD and AD. Thus, RD would be the most substantially affected in this regard as myelin disruption-associated RD changes are suspected to be the driving factor behind the diffusivity increases overall.

DWI has been used previously as a monitoring tool immediately after the MRgFUS treatment of bone in ex vivo lamb legs (41). Correlations were reported between applied energy, maximum temperature, and lesion volume, which was expected and is confirmed in the present study. Giles et al. also observed a positive correlation between apparent diffusion coefficient (ADC) and applied energy, maximum temperature, and lesion volume in muscle tissue adjacent to treated bone at time points ranging from <1–50 min after treatment. While this result differs from the anti-correlation currently presented, several important differences exist between studies: Giles et al. measure only ADC in muscle tissue adjacent to bone in room temperature ex vivo subjects at time points <50 min post-treatment. Conversely, we measured multiple DTI metrics in nervous tissue in vivo 1–2 h after sonication. The differences in acoustic energy attenuation between bone and soft tissue (42), thermal diffusion dynamics due to blood flow (43), and differing time points complicates direct comparison of the two studies. Further investigation is needed regarding the time evolution of MRgFUS lesions and microstructural dynamics of nerve and muscle tissue.

DTI metrics have been used longitudinally to evaluate sciatic nerve repair following peripheral nerve injury in varying levels of severity where nerve continuity is maintained (26, 27, 39) to complete disruption of axonal and surrounding connective tissue (44, 45). DTI has great potential in this regard as measurements are taken in vivo, precluding the need for large numbers of animals to be sacrificed at multiple time points. Due to limitations in data acquisition, these authors were able to reconstruct only small segments of the sciatic nerve and included spurious fibers. We present here robust tractography of the branches of the sciatic nerve from the dorsal root ganglion to their respective terminal muscle destinations, showing that accurate peripheral nerve tractography is achievable with appropriate acquisition parameters and post-processing techniques. Some of the above authors have used track length, normalized to a pre-surgical baseline, as a marker of nerve integrity. This can be problematic as the angle of nerve trajectory is a significant factor in fiber reconstruction and body position is not easily reproduced in repeated imaging sessions and between subjects. Nerve assessment via DTI metrics from a ROI, as we present here, may be a more robust method of evaluating nerve integrity.

The results of this study are consistent with previous experiments who reported similar observations of disrupted myelin and axon swelling in the sciatic nerves of pigs (19, 46) and rats (7) after treatment with FUS. These authors' observations are based on histology but are reinforced by our analysis of both histology and DTI metrics. This important validation of histological analysis by in vivo imaging obviates the need to sacrifice animals in order to assess the microstructural response to FUS treatment. DTI has been demonstrated previously as a useful adjunct to histological analysis (47, 48) and this study extends its application to MRgFUS.

We note that both Huisman et al. and Kaye et al. used large adult pigs (50–75 and 19–40 kg, respectively) with relatively large sciatic nerves (8–10 mm wide) (19, 46). Here we demonstrate the ability to image small animals (5–8 kg) and target small nerves (3–4 mm wide) with similar hardware. These smaller nerves reflect a similar size to potential human clinical MRgFUS targets such as the pudendal nerve (49) or posterior femoral cutaneous nerve (50), as noted by Huisman et al. (19).

Clinically, MRgFUS has potential as a non-invasive alternative to other ablation modalities, including radiofrequency ablation. Peripheral nerve ablation has been studied in a variety of cancerous (8, 9, 12, 13) and non-cancerous conditions (6, 10, 11). Ablation has also been used to treat tumors with close proximity to neural structures (14). For pain relief, a nerve conduction block may be indicated before ablative tract disruption is considered, as suggested by the findings of Foley et al. on FUS-related changes in nerve conduction (7). While conduction block can be successful in reducing pain, the effect can be limited in duration and nerve modulation or ablation may be further considered (51). Whether nerve conduction block or ablation is indicated, an in vivo imaging-based understanding of FUS-related nerve changes is important for characterizing treatments and predicting outcomes. This study demonstrates the utility of DWI to visualize and assess nerve changes in these treatments, thus bolstering the potential clinical utility of MRgFUS treatments.

We have demonstrated specificity on segmenting portions of the sciatic nerve based on vertebral origin and terminating muscle innervation (Supplementary Figure 1). The same seed ROIs at the nerve roots were used while disparate inclusion ROIs were placed at distal branches of the nerve. Similar techniques have been validated in the central nervous system for segmenting specific tract locations in thalamic nuclei by utilizing anatomically and functionally distinct cortical and subcortical regions in patients with essential tremor (25, 52). This demonstrates the potential to target specific portions of peripheral nerves in non-invasive treatments while avoiding branches of non-interest.

Limitations

This is an acute study and thus we are not able to determine the longitudinal effect of single-exposure MRgFUS treatment including effects at the lesion site and distal portions of the treated nerve. Previous studies have suggested that Wallerian-type degeneration may cause disruption in nearby white matter and along the length of the lesioned tract (53). Further studies incorporating multiple imaging time points are needed to understand the evolution of MRgFUS lesions and microstructural changes at both the primary injury site and distal tract segments. This study also uses ablative treatments of peripheral nerves, whereas lower intensity, non-ablative sonication may be indicated, as in the case of nerve conduction blocks (7). Future longitudinal studies may thus also incorporate lower power MRgFUS treatments to assess sub-lesional nerve changes and the inflection point between ablative and non-ablative therapies.

Technical challenges limit three-dimensional tractograms from being visualized directly on the MRgFUS targeting software. As such, specific fiber positions must be inferred in vivo relative to nearby anatomical structures and with reference to pre-treatment planning images. There exists the potential for targeting errors due to operator bias and positioning reproducibility. This problem is not limited to MRgFUS treatments but with surgical interventions incorporating the visualization of the central and peripheral nervous systems in general. Stereotactic frames in cranial applications can aid target registration but such devices are not generally available for peripheral nerve treatments. Further development in adequate forms of repeatable limb positioning and immobilization is needed for coupling tractography directly with non-invasive surgical approaches to minimize the effect of movement and maximize targeting accuracy.

Conclusions

This study demonstrates the potential of DWI and tractography for the in vivo targeting and assessment of changes following MRgFUS treatment of peripheral nerves. Significant changes in DTI metrics of FA, MD, AD, and RD were observed in the sciatic nerve following single exposure MRgFUS treatment. These DTI metric changes were correlated with sonication parameters. Robust visual reconstruction of the sciatic nerve via tractography was achieved. Diffusion imaging may thus be a valuable tool in optimizing peripheral nerve treatments with MRgFUS and evaluating the effects of treatment.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The animal study was reviewed and approved by Animal Care Committee and Laboratory Animal Services at the Hospital for Sick Children, Toronto, Ontario, Canada.

Author Contributions

MW, JZ, AW, KP, JD, and MH participated in the experimental design. Data was collected by MW, JZ, AW, and KP. The manuscript was written by MW. Manuscript critique and interpretation was contributed by JZ, AW, KP, LN, DA, JD, and MH. All authors contributed to this work and contributed to the analysis.

Funding

MH and JD are supported by a Brain Canada Multi-Investigator Research Initiative and MW by a Mitacs-Accelerate Canada Fellowship (IT05779).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank Marvin Estrada, Anson Lam, and the rest of the Laboratory Animal Services department at the Hospital for Sick Children for their animal care and handling for these experiments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.673060/full#supplementary-material

Supplementary Figure 1. CSD tractography of the sciatic nerve segmented by fiber bundle vertebral origin and distal termination (individual tracts separated by color). Full reconstructed tracts visible from (A) axial, (B) sagittal, and (C) coronal view. (D) Axial views showing only tract segments within the visible slice centered on the MRgFUS lesion. (E) Lesion zones I (red) and II (blue) are highlighted. CSD, constrained spherical deconvolution; MRgFUS, magnetic resonance-guided focused ultrasound.

References

1. Ghanouni P, Pauly KB, Elias WJ, Henderson J, Sheehan J, Monteith S, et al. Transcranial MRI-Guided Focused Ultrasound: A Review of the Technologic and Neurologic Applications. Am J Roentgenol. (2015) 205:150–9. doi: 10.2214/AJR.14.13632

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Elias WJ, Lipsman N, Ondo WG, Ghanouni P, Kim YG, Lee W, et al. A randomized trial of focused ultrasound thalamotomy for essential tremor. N Engl J Med. (2016) 375:730–9. doi: 10.1056/NEJMoa1600159

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chaussy CG, Thuroff S. High-intensity focused ultrasound for the treatment of prostate cancer: a review. J Endourol. (2017) 31:S-30-S-37. doi: 10.1089/end.2016.0548

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hesley GK, Felmlee JP, Gebhart JB, Dunagan KT, Gorny KR, Kesler JB, et al. Noninvasive treatment of uterine fibroids: early Mayo Clinic experience with magnetic resonance imaging-guided focused ultrasound. Mayo Clin Proc. (2006) 81:936–42. doi: 10.4065/81.7.936

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Liberman B, Gianfelice D, Inbar Y, Beck A, Rabin T, Shabshin N, et al. Pain palliation in patients with bone metastases using MR-guided focused ultrasound surgery: a multicenter study. Ann Surg Oncol. (2009) 16:140–6. doi: 10.1245/s10434-008-0011-2

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Foley JL, Little JW, Starr FL, Frantz C, Vaezy S. Image-guided HIFU neurolysis of peripheral nerves to treat spasticity and pain. Ultrasound Med Biol. (2004) 30:1199–207. doi: 10.1016/j.ultrasmedbio.2004.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Foley JL, Little JW, Vaezy S. Effects of high-intensity focused ultrasound on nerve conduction. Muscle Nerve. (2008) 37:241–50. doi: 10.1002/mus.20932

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Pu S, Wu J, Han Q, Zhang X, Lv Y, Xu Y, et al. Ultrasonography-guided radiofrequency ablation for painful stump neuromas to relieve postamputation pain: a pilot study. J Pain Res. (2020) 13:3437–45. doi: 10.2147/JPR.S283986

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hirbe AC, Jennings J, Saad N, Giardina JD, Tao Y, Luo J, et al. A Phase II Study of Tumor Ablation in Patients with Metastatic Sarcoma Stable on Chemotherapy. Oncologist. (2018) 23:760–e76. doi: 10.1634/theoncologist.2017-0536

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Kastler A, Aubry S, Barbier-Brion B, Jehl J, Kastler B. Radiofrequency neurolysis in the management of inguinal neuralgia: Preliminary study. Radiology. (2012) 262:701–7. doi: 10.1148/radiol.11110727

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kim HS, Wu PH, Jang IT. Lumbar degenerative disease part 1: anatomy and pathophysiology of intervertebral discogenic pain and radiofrequency ablation of basivertebral and sinuvertebral nerve treatment for chronic discogenic back pain: a prospective case series and review of lite. Int J Mol Sci. (2020) 21:1483. doi: 10.3390/ijms21041483

CrossRef Full Text | Google Scholar

12. Raslan AM, Cetas JS, Mccartney S, Burchiel KJ. Destructive procedures for control of cancer pain: the case for cordotomy: a review. J Neurosurg. (2011) 114:155–70. doi: 10.3171/2010.6.JNS10119

CrossRef Full Text | Google Scholar

13. Bentley JN, Viswanathan A, Rosenberg WS, Patil PG. Treatment of medically refractory cancer pain with a combination of intrathecal neuromodulation and neurosurgical ablation: Case series and literature review. Pain Med. (United States). (2014) 15:1488–95. doi: 10.1111/pme.12481

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kurup AN, Morris JM, Schmit GD, Atwell TD, Weisbrod AJ, Murthy NS, et al. Neuroanatomic considerations in percutaneous tumor ablation. Radiographics. (2013) 33:1195–215. doi: 10.1148/rg.334125141

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ciccarelli O, Catani M, Johansen-Berg H, Clark C, Thompson A. Diffusion-based tractography in neurological disorders: concepts, applications, and future developments. Lancet Neurol. (2008) 7:715–27. doi: 10.1016/S1474-4422(08)70163-7

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Basser PJ, Pierpaoli C. Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J Magn Reson B. (1996) 111:209–19. doi: 10.1006/jmrb.1996.0086

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hagmann P, Jonasson L, Maeder P, Thiran JP, Wedeen VJ, Meuli R. Understanding diffusion MR imaging techniques: from scalar diffusion-weighted imaging to diffusion tensor imaging and beyond. Radiographics. (2006) 26(Suppl 1):S205–23. doi: 10.1148/rg.26si065510

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Song S-K, Sun S-W, Ju W-K, Lin S-J, Cross AH, Neufeld AH. Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia. Neuroimage. (2003) 20:1714–22. doi: 10.1016/j.neuroimage.2003.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Huisman M, Staruch RM, Ladouceur-Wodzak M, van den Bosch MA, Burns DK, Chhabra A, et al. Non-Invasive Targeted Peripheral Nerve Ablation Using 3D MR Neurography and MRI-Guided High-Intensity Focused Ultrasound. (MR-HIFU): Pilot Study in a Swine Model. PLoS One. (2015) 10:e0144742. doi: 10.1371/journal.pone.0144742

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Budzik J-F, Verclytte S, Lefebvre G, Monnet A, Forzy G, Cotten A. Assessment of reduced field of view in diffusion tensor imaging of the lumbar nerve roots at 3 T. Eur Radiol. (2013) 23:1361–6. doi: 10.1007/s00330-012-2710-0

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kincses ZT, Szabó N, Valálik I, Kopniczky Z, Dézsi L, Klivényi P, et al. Target identification for stereotactic thalamotomy using diffusion tractography. PLoS One. (2012) 7:e29969. doi: 10.1371/journal.pone.0029969

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Tian Q, Wintermark M, Jeffrey Elias W, Ghanouni P, Halpern CH, Henderson JM, et al. Diffusion MRI tractography for improved transcranial MRI-guided focused ultrasound thalamotomy targeting for essential tremor. NeuroImage Clin. (2018) 19:572–80. doi: 10.1016/j.nicl.2018.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Chazen JL, Sarva H, Stieg PE, Min RJ, Ballon DJ, Pryor KO, et al. Clinical improvement associated with targeted interruption of the cerebellothalamic tract following MR-guided focused ultrasound for essential tremor. J Neurosurg. (2018) 129:1–9. doi: 10.3171/2017.4.JNS162803

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Sammartino F, Krishna V, King NKK, Lozano AM, Schwartz ML, Huang Y, et al. Tractography-Based Ventral Intermediate Nucleus Targeting: Novel Methodology and Intraoperative Validation. Mov Disord. (2016) 31:1217–25. doi: 10.1002/mds.26633

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Pouratian N, Zheng Z, Bari AA, Behnke E, Elias WJ, DeSalles AAF. Multi-institutional evaluation of deep brain stimulation targeting using probabilistic connectivity-based thalamic segmentation. J Neurosurg. (2011) 115:995–1004. doi: 10.3171/2011.7.JNS11250

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lehmann HC, Zhang J, Mori S, Sheikh KA. Diffusion tensor imaging to assess axonal regeneration in peripheral nerves. Exp Neurol. (2010) 223:238–44. doi: 10.1016/j.expneurol.2009.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Takagi T, Nakamura M, Yamada M, Hikishima K, Momoshima S, Fujiyoshi K, et al. Visualization of peripheral nerve degeneration and regeneration: Monitoring with diffusion tensor tractography. Neuroimage. (2009) 44:884–92. doi: 10.1016/j.neuroimage.2008.09.022

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. Neuroimage. (2012) 62:782–90. doi: 10.1016/j.neuroimage.2011.09.015

CrossRef Full Text | Google Scholar

29. Andersson JLR, Skare S, Ashburner J. How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. Neuroimage. (2003) 20:870–88. doi: 10.1016/S1053-8119(03)00336-7

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TEJ, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. (2004) 23 Suppl 1:S208–19. doi: 10.1016/j.neuroimage.2004.07.051

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Andersson JLR, Sotiropoulos SN. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage. (2016) 125:1063–78. doi: 10.1016/j.neuroimage.2015.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Veraart J, Sijbers J, Sunaert S, Leemans A, Jeurissen B. Weighted linear least squares estimation of diffusion MRI parameters: strengths, limitations, and pitfalls. Neuroimage. (2013) 81:335–46. doi: 10.1016/j.neuroimage.2013.05.028

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Westin C-F, Peled S, Gudbjartsson H, Kikinis R, Jolesz F. Geometrical diffusion measures for mri from tensor basis analysis. In: Proceedings of the 5th Annual Meeting of ISMRM'97. (2000). p. 1742.

Google Scholar

34. Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage. (2002) 17:825–41. doi: 10.1006/nimg.2002.1132

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Tournier J-D, Calamante F, Connelly A. MRtrix: diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol. (2012) 22:53–66. doi: 10.1002/ima.22005

CrossRef Full Text | Google Scholar

36. Tournier J-D, Calamante F, Connelly A. Determination of the appropriate b value and number of gradient directions for high-angular-resolution diffusion-weighted imaging. NMR Biomed. (2013) 26:1775–86. doi: 10.1002/nbm.3017

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Tournier J-D, Calamante F, Connelly A. Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage. (2007) 35:1459–72. doi: 10.1016/j.neuroimage.2007.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Tournier J-D, Yeh C-H, Calamante F, Cho K-H, Connelly A, Lin C-P. Resolving crossing fibres using constrained spherical deconvolution: Validation using diffusion-weighted imaging phantom data. Neuroimage. (2008) 42:617–25. doi: 10.1016/j.neuroimage.2008.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Yamasaki T, Fujiwara H, Oda R, Mikami Y, Ikeda T, Nagae M, et al. In vivo evaluation of rabbit sciatic nerve regeneration with diffusion tensor imaging. (DTI): correlations with histology and behavior. Magn Reson Imaging. (2015) 33:95–101. doi: 10.1016/j.mri.2014.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Wu W, Niu Y, Kong X, Liu D, Long X, Shu S, et al. Application of diffusion tensor imaging in quantitatively monitoring chronic constriction injury of rabbit sciatic nerves: correlation with histological and functional changes. Br J Radiol. (2017) 91:1083. doi: 10.1259/bjr.20170414

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Giles SL, Winfield JM, Collins DJ, Rivens I, Civale J, Ter Haar GR, et al. Value of diffusion-weighted imaging for monitoring tissue change during magnetic resonance-guided high-intensity focused ultrasound therapy in bone applications: an ex-vivo study. Eur Radiol Exp. (2018) 2:10. doi: 10.1186/s41747-018-0041-x

PubMed Abstract | CrossRef Full Text | Google Scholar

42. van Rhoon GC, Samaras T, Yarmolenko PS, Dewhirst MW, Neufeld E, Kuster N. CEM43C thermal dose thresholds: a potential guide for magnetic resonance radiofrequency exposure levels? Eur Radiol. (2013) 23:2215–27. doi: 10.1007/s00330-013-2825-y

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kolios MC, Sherar MD, Hunt JW. Blood flow cooling and ultrasonic lesion formation. Med Phys. (1996) 23:1287–98. doi: 10.1118/1.597694

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Li X, Chen J, Hong G, Sun C, Wu X, Peng MJ, et al. In vivo DTI longitudinal measurements of acute sciatic nerve traction injury and the association with pathological and functional changes. Eur J Radiol. (2013) 82:e707–14. doi: 10.1016/j.ejrad.2013.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Boyer RB, Kelm ND, Riley DC, Sexton KW, Pollins AC, Shack RB, et al. 4.7-T diffusion tensor imaging of acute traumatic peripheral nerve injury. Neurosurg Focus. (2015) 39:E9. doi: 10.3171/2015.6.FOCUS1590

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Kaye EA, Gutta NB, Monette S, Gulati A, Loh J, Srimathveeravalli G, et al. Feasibility Study on MR-Guided High-Intensity Focused Ultrasound Ablation of Sciatic Nerve in a Swine Model: preliminary results. Cardiovasc Intervent Radiol. (2015) 38:985–92. doi: 10.1007/s00270-015-1141-0

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Schilling KG, Janve V, Gao Y, Stepniewska I, Landman BA, Anderson AW. Histological validation of diffusion MRI fiber orientation distributions and dispersion. Neuroimage. (2018) 165:200–21. doi: 10.1016/j.neuroimage.2017.10.046

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Seehaus A, Roebroeck A, Bastiani M, Fonseca L, Bratzke H, Lori N, et al. Histological validation of high-resolution DTI in human post mortem tissue. Front Neuroanat. (2015) 9:98. doi: 10.3389/fnana.2015.00098

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Rofaeel A, Peng P, Louis I, Chan V. Feasibility of Real-Time Ultrasound for Pudendal Nerve Block in Patients with Chronic Perineal Pain. Reg Anesth Pain Med. (2008) 33:139–45. doi: 10.1016/j.rapm.2007.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Fritz J, Bizzell C, Kathuria S, Flammang AJ, Williams EH, Belzberg AJ, et al. High-resolution magnetic resonance-guided posterior femoral cutaneous nerve blocks. Skeletal Radiol. (2013) 42:579–86. doi: 10.1007/s00256-012-1553-8

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Choi EJ, Choi YM, Jang EJ, Kim JY, Kim TK, Kim KH. Neural ablation and regeneration in pain practice. Korean J Pain. (2016) 29:3–11. doi: 10.3344/kjp.2016.29.1.3

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Behrens TEJ, Johansen-Berg H, Woolrich MW, Smith SM, Wheeler-Kingshott CAM, Boulby PA, et al. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat Neurosci. (2003) 6:750–7. doi: 10.1038/nn1075

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Reginold W, Sam K, Poublanc J, Fisher J, Crawley A, Mikulis DJ. Impact of white matter hyperintensities on surrounding white matter tracts. Neuroradiology. (2018) 60:933–44. doi: 10.1007/s00234-018-2053-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: magnetic resonance-guided focused ultrasound, high intensity focused ultrasound, diffusion weighted imaging, diffusion tensor imaging, tractography, peripheral nerves, neuromodulation

Citation: Walker MR, Zhong J, Waspe AC, Piorkowska K, Nguyen LN, Anastakis DJ, Drake JM and Hodaie M (2021) Peripheral Nerve Focused Ultrasound Lesioning—Visualization and Assessment Using Diffusion Weighted Imaging. Front. Neurol. 12:673060. doi: 10.3389/fneur.2021.673060

Received: 26 February 2021; Accepted: 18 June 2021;
Published: 09 July 2021.

Edited by:

Vibhor Krishna, The Ohio State University, United States

Reviewed by:

J. Levi Chazen, Cornell University, United States
Rao P. Gullapalli, University of Maryland, United States

Copyright © 2021 Walker, Zhong, Waspe, Piorkowska, Nguyen, Anastakis, Drake and Hodaie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mojgan Hodaie, mojgan.hodaie@uhn.ca

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.