- 1Institute of Anatomy, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
- 2Jožef Stefan Institute, Ljubljana, Slovenia
- 3Department of Radiology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
- 4Clinical Institute of Radiology, University Medical Centre Ljubljana, Ljubljana, Slovenia
Introduction: Diffusion tensor imaging (DTI) is a magnetic resonance imaging (MRI) technique that measures the anisotropy of water diffusion. Clinical magnetic resonance imaging scanners enable visualization of the structural integrity of larger axonal bundles in the central nervous system and smaller structures like peripheral nerves; however, their resolution for the depiction of nerve fascicular morphology is limited. Accordingly, high-field strength MRI and strong magnetic field gradients are needed to depict the fascicular pattern. The study aimed to quantify diffusion tensor indices with high-field strength MRI within different anatomical compartments of the median nerve and determine if they correlate with nerve structure at the fascicular level.
Methods: Three-dimensional pulsed gradient spin-echo (PGSE) imaging sequence in 19 different gradient directions and b value 1,150 s/mm2 was performed on a 9.4T wide-bore vertical superconducting magnet. Nine-millimeter-long segments of five median nerve samples were obtained from fresh cadavers and acquired in sixteen 0.625 mm thick slices. Each nerve sample had the fascicles, perineurium, and interfascicular epineurium segmented. The diffusion tensor was calculated from the region-average diffusion-weighted signals for all diffusion gradient directions. Subsequently, correlations between diffusion tensor indices of segmentations and nerve structure at the fascicular level (number of fascicles, fascicular ratio, and cross-sectional area of fascicles or nerve) were assessed. The acquired diffusion tensor imaging data was employed for display with trajectories and diffusion ellipsoids.
Results: The nerve fascicles proved to be the most anisotropic nerve compartment with fractional anisotropy 0.44 ± 0.05. In the interfascicular epineurium, the diffusion was more prominent in orthogonal directions with fractional anisotropy 0.13 ± 0.02. Diffusion tensor indices within the fascicles and perineurium differed significantly between the subjects (p < 0.0001); however, there were no differences within the interfascicular epineurium (p ≥ 0.37). There were no correlations between diffusion tensor indices and nerve structure at the fascicular level (p ≥ 0.29).
Conclusion: High-field strength MRI enabled the depiction of the anisotropic diffusion within the fascicles and perineurium. Diffusion tensor indices of the peripheral nerve did not correlate with nerve structure at the fascicular level. Future studies should investigate the relationship between diffusion tensor indices at the fascicular level and axon- and myelin-related parameters.
1 Introduction
Certain medical conditions and penetrating injuries might cause individual nerve fascicles to be selectively more affected; therefore, accurate recognition of fascicle topography has the uttermost clinical importance (Härtig et al., 2018). As the pattern of fascicular involvement aid in the diagnostic workup of peripheral neuropathies, there is a great emphasis on its recognition. Clinical methods such as nerve conduction studies and electromyography are invasive, unpleasant, and give limited information. Accordingly, exploring available non-invasive radiologic modalities to extract information regarding nerve fascicular anatomy is imperative (Bilgen et al., 2005; Delgado-Martínez et al., 2016; Möller et al., 2018).
The application of clinical magnetic resonance imaging (MRI) for the depiction of peripheral nerves is increasing as it enables the assessment of peripheral neuropathies, nerve injuries or entrapments, and even tumors of the peripheral nerves (Chhabra et al., 2013; Bäumer et al., 2014). However, in trauma-related peripheral neuropathy, clinical MRI has limited accuracy in detecting pathologies, except in cases of severe nerve stretch or where the entire cross-section is affected (Eppenberger et al., 2014). Furthermore, when evaluating neoplasms of peripheral nerve with clinical MRI, differentiation of malignant from benign lesions can sometimes be difficult to achieve, even when there are characteristic signs of the malignancy (Chhabra et al., 2013). To surmount this obstacle, different strategies have been proposed. For example, the employment of advanced hardware, higher magnetic field magnets, and stronger gradients enabled the depiction of smaller structures as nerve fascicles; however, such studies are generally limited to ex vivo (Bilgen et al., 2005; Huang et al., 2015). Recently, specific MRI sequences such as fat-suppressed 3D fast low-angle shots have been proposed to improve the delineation of the nerve fascicles in healthy volunteers (Wang et al., 2023). Additionally, advanced techniques such as diffusion tensor imaging (DTI) have been exploited to further expand the options for depicting peripheral nerve pathologies (Khalil et al., 2008).
DTI enables measuring the effect of membranes on the apparent diffusion of water molecules. The specific arrangement of peripheral nerves results in diffusivity being predominantly directed along the axis of the nerve than in a perpendicular direction, thus being anisotropic. It is known that intact membranes are the primary determinant for anisotropic diffusion, with myelination having a modulating effect (Beaulieu, 2002). DTI measurements have been further utilized for deriving diffusion tensor (DT) indices that quantify the anisotropy (Pridmore et al., 2021). In rodent models, DTI has shown promise in distinguishing healthy, transected, and regenerating nerves (Lehmann et al., 2010). Moreover, even in macroscopic absent nerve discontinuity, DTI has been proven to detect minor nerve injuries (Boyer et al., 2015). DT indices also tend to correlate with behavioral changes and axonal density of rodents during the regeneration phases; therefore, they might serve as a promising tool in the future for recognizing unsuccessful nerve repair that requires further surgical intervention (Lehmann et al., 2010; Morisaki et al., 2011; Manzanera Esteve et al., 2019). In more recent research, DTI has also yielded a convenient tool for determining the severity of nerve injury in rats with the ability to distinguish different degrees of partial nerve transections (Manzanera Esteve et al., 2021).
As DT indices of peripheral nerve reflect the structural integrity of the nerve, they have attracted substantial attention for clinical application (Morisaki et al., 2011; Rangavajla et al., 2014). In clinical settings, it has been shown that demyelinating disease reflects in DT indices with the reduction of fractional anisotropy (FA), providing additional data regarding axonal degeneration in patients with peripheral neuropathies (Takagi et al., 2009; Kakuda et al., 2011). These data could complement clinical examination, electrophysiological evaluation, and conventional MRI for early recognition of patients with neuropathies who are eligible for neuroprotective therapies (Mathys et al., 2013). It has also been suggested that DTI could be an additional tool for assessing nerve compression syndromes, notably carpal tunnel syndrome (Khalil et al., 2008; Stein et al., 2009). In addition, DTI allows tractographic reconstructions; therefore, it can provide information on nerve integrity, predict the optimal approach of tumor resection, and predilect possible compromises in nerve impairment (Bruno et al., 2019).
Intra- and extra-fascicular structures of peripheral nerve likely possess different diffusion properties. It is unclear how this might affect clinical MRI scans, which provide an averaged diffusion-weighted signal. Notably, in vivo studies generally require image post-processing, which can affect the calculation of DT indices (Hiltunen et al., 2005; Stein et al., 2009). Accordingly, basic research is needed to further understand nerve DT indices obtained on clinical MRI. Depiction of the diffusion process at the fascicular level within different nerve compartments (fascicles, perineurium, and interfascicular epineurium) could enhance the understanding of nerve DT indices. High-field MRI is required to further depict such small structures. The present study aimed to obtain knowledge on the diffusion characteristics of different nerve compartments on high-field MRI and determine their relationship with nerve fascicular morphometric parameters.
2 Materials and methods
2.1 Sample preparation
A segment of the median nerve was obtained from the proximal upper arm of five fresh cadavers, less than 24 h postmortem. The cadavers were donated for research and educational purposes to the Institute of Anatomy, Faculty of Medicine, University of Ljubljana, through a willed cadaver donation program. Each nerve was cut into a 9 mm long segment, had carefully removed the surrounding connective tissue, and inserted in a 10-mm-diameter glass tube. To prevent sample dehydration, the glass tube was filled with perfluorinated liquid, Galden SV90 from Solvay (Brussels, Belgium), which does not produce any detectable MRI signal (Awais et al., 2022). The study was approved by the National Medical Ethics Committee of the Republic of Slovenia (Permit No: 0120–239/2020/3).
2.1.1 Nerve donors
The donors were four females and one male, with a mean age of 75 years (range 70–80). The interval from death to nerve sampling ranged from 5 to 23 h. There was limited data regarding premortem medical conditions, but all subjects had an atherosclerotic disease and/or arterial hypertension with varying degrees of severity. None of the subjects had known peripheral nerve disease.
2.2 Magnetic Resonance Microscopy image acquisition
Magnetic Resonance Microscopy (MRM) was performed on a 9.4T (400 MHz proton frequency) wide-bore vertical superconducting magnet (Jastec Superconductor Technology, Tokyo, Japan) connected to an NMR/MRI spectrometer (Tecmag, Houston TX, United States). Before the imaging, the tube with the sample was inserted in a Micro 2.5 gradient system with a 10 mm RF probe (Bruker, Ettlingen, Germany) of the magnet.
DTI of the nerves was performed using a three-dimensional (3D) pulsed gradient spin-echo (PGSE) imaging sequence with diffusion gradients in 19 different directions; however, all with the same b value of 1,150 s/mm2. The selected b value was chosen given the preliminary results, whereas we have tested various b values up to 1,800 s/mm2. The selected value provided optimal conditions for measuring the leading eigenvalue within the nerve fascicles. The theory also supports the selected b value for the two-point experiment with b1 = 0 and b2 = b > 0, where the optimal b value is equal to b = 1.1/D (Xing et al., 1997). Acquisition of an additional reference T2-weighted image with no diffusion weighting (b = 0) was needed for DTI calculation. The images were acquired with the following parameters: TE/TR = 36/880 ms; δ = 3 ms; ∆ = 27 ms; G0 = 0.26 T/m; field of view 9 × 4.5 × 10 mm3; matrix size, 256 × 128 × 16; and 4 signal averages. The image resolution along the in-plane directions was 35 μm. Scanning was performed at room temperature of 21°C with a total acquisition time of 1 day 16 h.
2.3 Image analysis and nerve morphometry
The nerve segments of peripheral nerves, acquired in 16 continuous slices of 0.625 mm thickness, were identified on reference T2-weighted images (Figure 1A). Quality assessment of slices was performed, and slices with artifacts and partial volume effect were excluded from further analysis. In each included slice, fascicles, interfascicular epineurium, perineurium, and nerve cross-sectional area (CSA) were segmented. The fascicles were defined as intraneural hypointense oval- or round-shaped tissue circumferentially surrounded by a markedly hyperintense line representing the perineurium. The latter served as a reliable segmentation border (Figure 1B). The perineurium was segmented with a single measurement by two parallel lines, as shown in Figure 1C. The hyperintense tissue between the fascicles was defined as interfascicular epineurium (Figure 1D). The nerve was segmented to include the entire nerve but a minimal proportion of the background (Figure 1E). Segmentations were performed manually with the image processing software ImageJ (National Institutes of Health, Bethesda, Maryland, United States). The area was recorded for each segmentation and expressed as CSA for the nerve and fascicles. The fascicular ratio (FR) was calculated as a net fascicular CSA/nerve CSA, the ratio of perineurium as a net perineurium/nerve CSA, and the ratio of interfascicular epineurium as net interfascicular epineurium/nerve CSA (Tagliafico and Tagliafico, 2014).
FIGURE 1. T2-weighed images, using b-value 0 s/mm2, displaying 0.625 mm thick representative slices of five analyzed median nerves. Note that each figure represents one out of sixteen slices. In figure (A) fascicles are sharply demarcated with a hyperintense line representing the perineurium. Further images depict the segmentation of (B) eight nerve fascicles, (C) thin layer of perineurium with two parallel lines, (D) interfascicular epineurium, and (E) nerve cross-sectional area.
The diffusion tensor was calculated from the acquired three-dimensional data as described previously (Basser et al., 1994; Awais et al., 2022). For each image voxel, the calculated diffusion tensor was diagonalized, which yielded maps of the tensor eigenvalues D1, D2, and D3 and of the corresponding eigenvectors (
2.4 Intra-observer agreement
A subset of 10 nerve fascicles, interfascicular epineurium areas, perineurial areas, and nerve CSA were randomly selected and segmented again by the same observer 30 days after the primary segmentation to assess intra-observer agreement. Intraclass correlation coefficient (ICC) was calculated from the FA (Koo and Li, 2016). FA was chosen for ICC calculation as this index is most commonly used DTI readout parameter in clinical environment reflecting the degree of cellular structure alignment (Kronlage et al., 2018).
2.5 Trajectory, diffusion ellipsoids, and color-coded ellipsoid/fiber orientation display
The DT data of the acquired nerve segments were displayed with trajectories and diffusion ellipsoids rendered with POV-Ray software (Persistence of Vision Pty. Ltd., version 3.7, Williamstown, Victoria, Australia) (Awais et al., 2022). The software generated a tractography display of the entire peripheral nerve length using the components of the first eigenvector. The subsequent slices were displayed with ellipsoids whose size and orientation correspond to the eigenvalues (size) and the eigenvectors (orientation).
2.6 Statistical analysis
Statistical analysis was performed using GraphPad Prism 9 (GraphPad Software Inc., San Diego, United States). The Shapiro-Wilk test was used to evaluate the groups for normality. Because normality and equal variance assumptions were met, the fascicular eigenvalues (D1, D2, and D3), as well as their derived indices (MD, FA, and D||/D⊥), were compared by two-way analysis of variance (ANOVA) followed by Tukey’s posthoc test when appropriate. When comparing DT indices of the perineurium, interfascicular epineurium, and nerve CSA one-way ANOVA followed by Tukey’s posthoc test was employed. The fascicular variability of FA was assessed using the coefficient of variation and then compared between fascicles and within fascicles with two-way ANOVA followed by Tukey’s posthoc test when appropriate. To determine the correlations between the DT indices and parameters at the fascicular level, linear regression was performed for each nerve sample then coefficients were compared using a one-sample t-test. The change of nerve FA in subsequent slices was calculated with linear regression. For the assessment of an intra-observer agreement, one-way ICC was used (Koo and Li, 2016). Differences were deemed statistically significant at p < 0.05. Data are given as means ± standard deviations, ranges, or percentages when appropriate.
3 Results
3.1 Nerve morphometric characteristics
After quality assessment, 56 image slices were included in the study (range of slices per nerve, 8–13). There were 7.94 ± 4.33 fascicles per slice with a mean fascicle CSA of 0.57 ± 0.66 mm2. The CSA of the nerve was 12.34 ± 3.53 mm2, and the FR was 0.46 ± 0.07. The ratio between perineurium/nerve CSA was 0.08 ± 0.03, and the ratio between interfascicular epineurium/nerve CSA was 0.46 ± 0.09.
3.2 DTI characteristics of nerve compartments
In the nerve fascicles, the eigenvalue D1, with an average of 0.81 ± 0.09·10−9 m2/s, was the highest and approximately 2-times higher than eigenvalues D2 or D3. The mean fascicular eigenvalue D1 was 27-times higher than in the interfascicular epineurium but 1.32-times lower than the mean eigenvalue D1 of the perineurium (Tables 1, 2, 3). There were significant differences between nerve samples regarding the fascicular and perineurium eigenvalues (p < 0.0001 and p < 0.0001, respectively), while there were no significant differences between nerve samples regarding the interfascicular epineurium eigenvalues.
The mean MD followed the same pattern as eigenvalues, with the highest values calculated in the perineurial compartment. Compared to the fascicles, the perineurium had a 50.60% ± 20.72% higher MD, and the interfascicular epineurium had a 94.22% ± 0.80% lower MD (Tables 1, 2, 3). There were significant differences between nerve samples regarding the fascicular and perineural MD (p < 0.0001 and p < 0.0001, respectively), while there were no significant differences between nerve samples regarding the MD of the interfascicular epineurium.
The fascicle was the most anisotropic peripheral nerve compartment, with a mean FA of 0.44 ± 0.05. The coefficients of variation of fascicular FA throughout the same fascicle on sequential slices and between fascicles on the same slice were 0.17 ± 0.06 and 0.16 ± 0.03, respectively, and showed no statistically significant difference. Compared to the fascicles, the mean FA was lower in the perineurium (−23.95% ± 4.05%) and even lower in the interfascicular epineurium (−70.38% ± 3.91%). D||/D⊥ had the highest and most anisotropic values calculated in the fascicular compartment, while the interfascicular epineurium was the most isotropic compartment with a mean D||/D⊥ of 1.15 ± 0.04. Nerve samples differed significantly in FA and D||/D⊥ within the fascicles (p < 0.0001 and p < 0.0001, respectively) and perineurium (p = 0.0001 and p < 0.0001 respectively), while there were no significant differences in FA and D||/D⊥ of interfascicular epineurium (Tables 1, 2, 3).
DTI maps of five nerve segments providing eigenvalues (D1, D2, and D3), MD, and FA are included in the Supplementary Materials (Supplementary Figures 1–5).
3.3 DTI characteristics of nerve cross-section
The mean eigenvalues D1, D2, and D3 of the nerve were 0.46 ± 0.23·10−9 m2/s, 0.30 ± 0.15·10−9 m2/s, and 0.28 ± 0.14·10−9 m2/s, respectively, and differed significantly between the samples (p < 0.0001). Compared to the fascicular eigenvalues, mean nerve eigenvalues D1, D2, and D3 were 43.33% ± 23.07%, 27.87% ± 17.33%, and 24.45% ± 19.24% lower, respectively.
The MD of the nerve was 0.34 ± 0.17·10−9 m2/s. This was approximately 11-times higher than the MD of the interfascicular epineurium but 1.51 and 2.29-times lower than the nerve fascicles and perineurium, respectively. The MD differed significantly between the nerve samples (p < 0.0001) (Figure 2D).
FIGURE 2. Diffusion tensor indices that are calculated from the cross-sectional area of different nerve samples. Figures compare (A) eigenvalue D1, (B) eigenvalue D2, (C) eigenvalue D3, (D) mean diffusivity (MD), (E) fractional anisotropy (FA), and (F) D||/D⊥. Data are presented as means and standard deviations. ****p < 0.0001 compared to nerves 1, 2, 3, and 5; ####p < 0.0001 compared to nerves 1, 2, 3, and 4.
The mean FA of the nerve was 0.28 ± 0.04. The highest FA was noted in nerve sample 5 (Figure 2E), which differed significantly from others (p < 0.0001). The mean coefficient of variation calculated from the FA of nerve samples was 0.12 ± 0.07. The nerve FA was approximately one-third and one-sixth lower than the FA of fascicles and perineurium, respectively, and 2-times higher than the FA of the interfascicular epineurium. The FA of the fascicles and perineurium measured together was 22.85% ± 7.75% higher than the nerve FA. Nerve samples 1 and 2 had no change of FA detected throughout the nerve segments; however, statistically significant correlations were observed in nerves 3–5 (Figure 3). D||/D⊥ showed the same pattern as FA with the highest values in nerve 5 (Figure 2F), which differed significantly from other samples (p < 0.0001).
FIGURE 3. Change of fractional anisotropy throughout the nerve segments. The figure depicts the fractional anisotropy of the nerve cross-sectional area in consecutive slices. (A) Nerve sample 1 and (B) nerve sample 2 had no change of FA detected throughout the segment; however, statistically significant correlations were observed in (C) nerve 3, (D) nerve 4, and (E) nerve 5. Note that the orientation of the nerve (i.e., proximal/distal part) was not tracked during the sample preparation process.
3.4 Correlations between DT indices and nerve structure at the fascicular level
When evaluating the correlation between the nerve FA and the average FA of all fascicles on the same slice, we noted a moderately strong correlation (r = 0.74, p = 0.001). Additionally, we found a correlation between nerve FA and FA of the largest fascicle within all five nerves (r ≥ 0.47, p < 0.05).
However, we found no correlations when evaluating correlations between DT indices and structures at the fascicular level. No correlation was found between the DT indices of fascicles and CSA of fascicles (p ≥ 0.29). There was also no correlation between the mean FA of all fascicles on the slice and the number of fascicles (p = 0.59). We observed no correlation between the DT indices of nerve and the number of fascicles (p = 0.88). Additionally, there was no correlation between the nerve FA and FR (p = 0.34), and no correlation between the nerve FA and CSA of all fascicles on the same slice (p = 0.68).
3.5 Intra-observer agreement
ICC calculated from FA of the fascicles, and the nerve CSA showed excellent intra-observer agreement, 0.98 and 0.99, respectively. Good intra-observer agreement was observed for the perineurium, 0.89, and interfascicular epineurium, 0.86.
3.6 Anisotropic diffusion in nerves presented with tractography and diffusion ellipsoids
DT data was demonstrated with graphic displays. Tractographic displays were generated using the fastest diffusion direction, which was oriented longitudinally along the course of the nerve. The three-dimensional representation of anisotropic nerve fibers of one nerve sample is shown in Figure 4A. As seen, nerve fascicles can be tracked along the entire nerve segment, their fiber density remains constant, and the fibers intermingle within the fascicles.
FIGURE 4. DT tractography and diffusion ellipsoids of one nerve sample. (A) The figure displays a three-dimensional representation of nerve fibers in an approximately 9 mm long segment of the median nerve. By convention, red color represents the x-direction, green color the y-direction, and blue color the z-direction; in our case, the nerve was oriented along the z-axis. (B) Diffusion ellipsoid display of the subsequent slices. The dimensions and orientations for the ellipsoids correspond to the eigenvalues and eigenvectors, while their color scheme is orientation-dependent and determined by the same principle as the colors in the tractography image.
The presentation with diffusion ellipsoids shows ellipsoids that are oriented in the direction of the fastest diffusion (along the first eigenvector). The main axis of the ellipsoid provides information about the main diffusion direction in the voxel, while the shape of the ellipsoid gives the information about the degree of anisotropy. The difference in the eigenvalues within the fascicles and perineurium enables both compartments to be depicted and separated (Figure 4B). Note different eigenvectors of perineurium and fascicles that differ in values and orientations. In contrast, the interfascicular epineurium with considerably slower and isotropic diffusion cannot be adequately visualized. This anatomical compartment can be partially visualized in the background, where it can be seen as small dots of different colors and shapes.
4 Discussion
In our study, DTI in the magnetic field strength of 9.4T was employed on fresh ex vivo human median nerves, and DT indices of all major nerve anatomical structures were quantified. The anisotropic diffusion was shown throughout the long axis of the nerve, and nerve fascicles proved to be the most anisotropic nerve compartment. The DT indices of fascicles and perineurium differed significantly among subjects, while the interfascicular epineurium with the slowest and practically isotropic diffusion had no inter-subject differences. The DT indices of the peripheral nerve and anatomical compartments of the nerve did not correlate with nerve structure at the fascicular level.
DTI of peripheral nerves in the upper and lower extremities has been previously performed (Eppenberger et al., 2014; Kronlage et al., 2018; Godel et al., 2019). In the upper extremity, most investigations have been performed on the median nerve as it is the most frequently affected nerve in upper extremity entrapment neuropathies (Snoj et al., 2022). Several studies have depicted tractographic images of healthy median nerves and showed that Wallerian degeneration reduces the attenuation of trackable nerve fibers (Hiltunen et al., 2005; Khalil et al., 2008; Takagi et al., 2009; Boyer et al., 2015). More recent studies have predominately focused on the calculation of DT indices of nerve, mainly the FA (Kronlage et al., 2018; Awais et al., 2022). The nerve FA might vary along the longer nerve segment. Yao and Gai reported no change in FA along the length of the median nerve in the carpal tunnel (Yao and Gai, 2009), while the few other researchers demonstrated a decreasing trend of FA from proximal to distal locations near the carpal tunnel (Hiltunen et al., 2005; Guggenberger et al., 2013). Stein et al. (2009) showed that FA could differ significantly in a few centimeters long nerve segment. In our study, the FA did not diminish or augment throughout the nerve segment in two subjects, while the other three subjects had a trend of changing FA in the subsequent slices.
The mean FA of healthy median nerve reported in more extensive meta-analysis was 0.58 (Rojoa et al., 2021). In the wrist of healthy individuals, the mean FA was found to be in a broader interval between 0.48–0.71 (Kabakci et al., 2007; Stein et al., 2009; Barcelo et al., 2013; Guggenberger et al., 2013). Only the scarcity of studies have measured the FA of nerves in the upper arm and reported values in a similar interval range (Kronlage et al., 2018; Godel et al., 2019). It has been previously shown that the age of the subject is an important determinant for nerve FA (Tanitame et al., 2012; Kronlage et al., 2018). When accounting for this factor, the FA of the median nerve in our study was approximately two times lower compared to the FA of healthy peripheral nerves of similarly aged subjects in a study by Kronlage et al. (2018). Several factors should be accounted for when interpreting lower FA in our study. The most important factors are likely environmental. In our study, the MRM acquisition was performed at room temperature. This caused the diffusion process to be approximately 40% lower than the diffusion process at body temperature and may partially explain the differences in the values of DT indices compared to indices measured in previous in vivo studies (Holz et al., 2000; Lehmann et al., 2010). In a few ex vivo studies, unfrozen or fixated nerves were used (Boyer et al., 2015; Awais et al., 2022). Sample preparation can have an important impact on its integrity. It has been shown that tissue fixation in formaldehyde can significantly decrease the FA of a heart muscle (Lohr et al., 2020). Conversely, Haga et al. (2019) did not observe a change in FA between fixed and non-fixed marmoset brains; however, the use of formaldehyde did significantly decrease eigenvalues and MD of the fixed brain. As there is a lack of data on how formaldehyde might affect DT indices of peripheral nerve, we have used a fresh nerve to exclude the effects of the fixative procedure or possible rupture of the nerve cells during the freezing/thawing cycle.
The outstanding FA in nerve 5 can be partly attributed to the low body mass index of this donor; nevertheless, such inter-individual differences can still be found between healthy individuals (Kronlage et al., 2018). As reported by Godel et al. (2019), increased radial diffusivity reflects damage to myelin integrity, whereas changes in axonal diffusivity might be more specific for axonal degeneration. In nerve 4, eigenvalues D1, D2, and D3 were equally increased which probably supports the hypothesis of a major role of environmental factors.
Differences in MR hardware and imaging protocols can also lead to discrepancies between studies. Most existing studies were performed on conventional 3T MRI (Stein et al., 2009; Barcelo et al., 2013; Guggenberger et al., 2013; Kronlage et al., 2018); however, some researchers have also applied 7T whole-body MRI systems (Riegler et al., 2016; Yoon et al., 2018). Although a high-field-strength system was used in our study, it has been previously shown on the brain that field strength has little effect on the FA (Zhan et al., 2013). Importantly, studies use scanners from different manufacturers, and it has been shown that FA of the median nerve in healthy individuals significantly differs within the wrist between different MR scanners (Guggenberger et al., 2013). Post-processing within in vivo studies involves applying threshold values to distinguish nerve from muscle fibers or ligaments (Hiltunen et al., 2005). In these studies, FA is calculated from the fiber tractography images, which normally have a high threshold value (Hiltunen et al., 2005; Kabakci et al., 2007; Stein et al., 2009; Barcelo et al., 2013). Thus, the impact of isotropic nerve compartments, such as interfascicular epineurium, that decreases FA of nerve CSA is excluded from the calculation to a certain degree. In our study, no threshold was applied to obtain the most reliable data. Consequently, the FA of fascicles and perineurium combined was approximately a quarter higher than the FA of the nerve, which included interfascicular epineurium.
The estimated signal-to-noise ratio (SNR) in the fascicular region was approximately 14 (b = 0 s/mm2). In comparison to the study by Yoon et al. (Yoon et al., 2018), this SNR was comparable to their in vivo 7T MRI system. However, our in-plane resolution was considerably better (35 μm vs. 120 μm), and slices were also thinner (0.625 mm vs. 2 mm). The resolution advantage in our experiment is due to the more sensitive receiver coil, which had only 10 mm in diameter; thus, the filling factor (nerve to coil diameter ratio) was high, and the signal reception was, therefore, considerably better than in the in vivo experiment at 7T where a surface coil was used. Some advantage was also due to a somewhat stronger magnetic field (9.4T vs. 7T). This comparison demonstrates that in vivo MR imagining at 7T is promising but still needs several improvements, especially in signal detection, to match MR microscopy results at higher fields and dedicated hardware (special gradient and signal receive coils).
The long scanning window of approximately 40 h posed a challenge concerning tissue desiccation. Immersion in formaldehyde would be optimal for preventing the autolytic process; however, it contains hydrogen atoms that produce a signal on MRI. Hence, each nerve was placed instantly after the excision into the fluorinated carbon liquid, which substantially reduced problems with stability and autolysis. In our previous pilot study (Awais et al., 2022), minor nerve shrinkage of one pixel was noted during the scanning of the nerve sample in this liquid. This issue was adequately addressed with an innovative scanning strategy by reordering the scanning loops, whereas all twenty images were acquired simultaneously, not sequentially. Thus, the influence of the sample volume change was evenly distributed among all images and did not importantly affect the calculations of DT indices.
In our study, fascicles proved to have the highest FA in the analysis of the nerve compartments. This result was expected due to the specific arrangement of nerve fibers in the peripheral nerve (Sunderland, 1978). It is well known that the fascicular pattern can change in submillimeter sections (Sunderland, 1945; Sunderland, 1978). When evaluating the FA within the same fascicle, we found no statistical difference between the fascicular coefficients of variation in sequential slices and the fascicles of the same slice. We hypothesize that the fibers crossing between the fascicles contributed to high fascicular coefficients of variation within the same fascicle and explain why no difference was observed (Figley et al., 2022).
The perineurium had only moderately lower FA than the fascicles. The latter provides a diffusion barrier made of concentrically flat perineural cells, which contributes to poorer membrane permeability for water molecules and could lead to higher FA (Hill and Williams, 2002). As the perineurium is thinner than the resolution of our MRM system (Reina et al., 2015), it is essential to consider the partial volume effect. The segmentation of perineurium might partially include highly anisotropic fascicles and isotropic interfascicular epineurium. The interfascicular epineurium, as a collagenous compartment of the extracellular matrix, does not form any non-permeable barrier. Hence, it had the lowest FA mean approaching isotropic diffusion when considering the influence of the SNR (Peltonen et al., 2013; Awais et al., 2022). Fascicle FA has an important impact on nerve FA; however, measured fascicular parameters did not affect the nerve FA. Previous studies have shown that higher FA correlates with increased axonal density, axonal diameter, and myelin density (Takagi et al., 2009). Thus, demyelinating disorders could contribute to differences in fascicular FA between subjects included in this study (Kakuda et al., 2011). As seen, the FA was not correlating with structures at the fascicular level; therefore, axon- and myelin-related parameters tend to have a more important role in diffusion (Takagi et al., 2009).
We acknowledge that this study had some limitations. First, we had limited clinical data on the subjects from whom we obtained nerve samples. Consequently, the differences in DT indices between the subjects might have been even greater than they would have been between healthy subjects. It would be imperative in future research to expand the sample size and provide more clinical data about pathologies that could affect nerve integrity. Another limitation of this study is the acquisition at room temperature. As the environmental temperature could not be strictly controlled, minor fluctuations in room temperature could partially impact the DT indices of individual samples and cause important differences between the samples. A third limitation is that the determination of individual compartments could occasionally be ambiguous; hence, individual segmentations of the perineurium or smaller fascicles were more challenging to implement. This could lead to partially overlapping regions. Another limitation is a relatively long scanning time resulting in the subjection of nerve samples to the autolytic process. Accordingly, this has to be considered when compared to in vivo studies or ex vivo studies with shorter acquisition time. The last limitation is the limited ability to translate our results directly into the partial trauma series; non-etheless, we believe basic knowledge is essential for future understanding of how partial nerve transection reflects in a change of DT indices within different nerve compartments. Moreover, the understanding of diffusion within the nerve compartments can be translated into the nerve entrapment syndromes (where oedema is present) or diabetic neuropathy and thus help evaluate how individual compartments contribute to the anisotropy change.
5 Conclusion
High-resolution DTI depicted highly anisotropic diffusion within the fascicles and perineurium. The interfascicular epineurium had more isotropic diffusion. As median nerve DT indices did not correlate with nerve structure at the fascicular level, future studies should investigate their relationship with axon- and myelin-related parameters.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by National Medical Ethics Committee of the Republic of Slovenia (Permit No: 0120–239/2020/3). The participants gave their written consent to be donated for research and educational purposes.
Author contributions
Conceptualization: LP, IS, NU, EC, and ŽS; methodology: IS and ŽS; software: IS; formal analysis: LP; investigation: LP; resources: IS, NU, EC, and ŽS; writing—original draft preparation: LP; supervision: IS, NU, EC, and ŽS. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the Slovenian Research Agency grants J3-9288, J3-4507, N3-0256, and P3-0043; and by the Czech Science Foundation (GACR) grant 22-02756K.
Acknowledgments
We are thankful to Majda Črnak Maasarani for technical support, and to Dr. Chiedozie Kenneth Ugwoke for manuscript proofreading.
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.
Publisher’s note
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1070227/full#supplementary-material
References
Awais, K., Snoj, Ž., Cvetko, E., and Serša, I. (2022). Diffusion tensor imaging of a median nerve by magnetic resonance: A pilot study. Life 12, 748. doi:10.3390/life12050748
Barcelo, C., Faruch, M., Lapègue, F., Bayol, M. A., and Sans, N. (2013). 3-T MRI with diffusion tensor imaging and tractography of the median nerve. Eur. Radiol. 23, 3124–3130. doi:10.1007/s00330-013-2955-2
Basser, P. J., Mattiello, J., and LeBihan, D. (1994). MR diffusion tensor spectroscopy and imaging. Biophys. J. 66, 259–267. doi:10.1016/S0006-3495(94)80775-1
Bäumer, P., Pham, M., Ruetters, M., Heiland, S., Heckel, A., Radbruch, A., et al. (2014). Peripheral neuropathy: Detection with diffusion-tensor imaging. Radiology 273, 185–193. doi:10.1148/radiol.14132837
Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. 15, 435–455. doi:10.1002/nbm.782
Bilgen, M., Heddings, A., Al-Hafez, B., Hasan, W., McIff, T., Toby, B., et al. (2005). Microneurography of human median nerve. J. Magnetic Reson. Imaging 21, 826–830. doi:10.1002/jmri.20345
Boyer, R. B., Kelm, N. D., Riley, D. C., Sexton, K. W., Pollins, A. C., Shack, R. B., et al. (2015). 4.7-T diffusion tensor imaging of acute traumatic peripheral nerve injury. Neurosurg. Focus 39, E9–E9. doi:10.3171/2015.6.FOCUS1590
Bruno, F., Arrigoni, F., Mariani, S., Patriarca, L., Palumbo, P., Natella, R., et al. (2019). Application of diffusion tensor imaging (DTI) and MR-tractography in the evaluation of peripheral nerve tumours: State of the art and review of the literature. Acta Biomed. 90, 68–76. doi:10.23750/abm.v90i5-S.8326
Chhabra, A., Thakkar, R. S., Andreisek, G., Chalian, M., Belzberg, A. J., Blakeley, J., et al. (2013). Anatomic MR imaging and functional diffusion tensor imaging of peripheral nerve tumors and tumorlike conditions. Am. J. Neuroradiol. 34, 802–807. doi:10.3174/ajnr.A3316
Delgado-Martínez, I., Badia, J., Pascual-Font, A., Rodríguez-Baeza, A., and Navarro, X. (2016). Fascicular topography of the human median nerve for neuroprosthetic surgery. Front. Neurosci. 10, 286. doi:10.3389/fnins.2016.00286
Eppenberger, P., Andreisek, G., and Chhabra, A. (2014). Magnetic resonance neurography: Diffusion tensor imaging and future directions. Neuroimaging Clin. N. Am. 24, 245–256. doi:10.1016/j.nic.2013.03.031
Figley, C. R., Uddin, M. N., Wong, K., Kornelsen, J., Puig, J., and Figley, T. D. (2022). Potential pitfalls of using fractional anisotropy, axial diffusivity, and radial diffusivity as biomarkers of cerebral white matter microstructure. Front. Neurosci. 15, 799576–799577. doi:10.3389/fnins.2021.799576
Godel, T., Pham, M., Kele, H., Kronlage, M., Schwarz, D., Brunée, M., et al. (2019). Diffusion tensor imaging in anterior interosseous nerve syndrome – functional MR Neurography on a fascicular level. Neuroimage Clin. 21, 101659. doi:10.1016/j.nicl.2019.101659
Guggenberger, R., Nanz, D., Bussmann, L., Chhabra, A., Fischer, M. A., Hodler, J., et al. (2013). Diffusion tensor imaging of the median nerve at 3.0T using different MR scanners: Agreement of FA and ADC measurements. Eur. J. Radiol. 82, e590–e596. doi:10.1016/j.ejrad.2013.05.011
Haga, Y., Hata, J., Uematsu, A., Seki, F., Komaki, Y., Mizumura, M., et al. (2019). Mr imaging properties of ex vivo common marmoset brain after formaldehyde fixation. Magnetic Reson. Med. Sci. 18, 253–259. doi:10.2463/mrms.mp.2018-0086
Härtig, F., Ross, M., Dammeier, N. M., Fedtke, N., Heiling, B., Axer, H., et al. (2018). Nerve ultrasound predicts treatment response in chronic inflammatory demyelinating polyradiculoneuropathy—A prospective follow-up. Neurotherapeutics 15, 439–451. doi:10.1007/s13311-018-0609-4
Hill, R. E., and Williams, P. E. (2002). A quantitative analysis of perineurial cell basement membrane collagen IV, laminin and fibronectin in diabetic and non-diabetic human sural nerve. J. Anat. 201, 185–192. doi:10.1046/j.1469-7580.2002.00083.x
Hiltunen, J., Suortti, T., Arvela, S., Seppä, M., Joensuu, R., and Hari, R. (2005). Diffusion tensor imaging and tractography of distal peripheral nerves at 3 T. Clin. Neurophysiol. 116, 2315–2323. doi:10.1016/j.clinph.2005.05.014
Holz, M., Heil, S. R., and Sacco, A. (2000). Temperature-dependent self-diffusion coefficients of water and six selected molecular liquids for calibration in accurate 1H NMR PFG measurements. Phys. Chem. Chem. Phys. 2, 4740–4742. doi:10.1039/b005319h
Huang, S. Y., Nummenmaa, A., Witzel, T., Duval, T., Cohen-Adad, J., Wald, L. L., et al. (2015). The impact of gradient strength on in vivo diffusion MRI estimates of axon diameter. Neuroimage 106, 464–472. doi:10.1016/j.neuroimage.2014.12.008
Kabakci, N., Gürses, B., Firat, Z., Bayram, A., Uluǧ, A. M., Kovanlikaya, A., et al. (2007). Diffusion tensor imaging and tractography of median nerve: Normative diffusion values. Am. J. Roentgenol. 189, 923–927. doi:10.2214/AJR.07.2423
Kakuda, T., Fukuda, H., Tanitame, K., Takasu, M., Date, S., Ochi, K., et al. (2011). Diffusion tensor imaging of peripheral nerve in patients with chronic inflammatory demyelinating polyradiculoneuropathy: A feasibility study. Neuroradiology 53, 955–960. doi:10.1007/s00234-010-0833-z
Khalil, C., Hancart, C., le Thuc, V., Chantelot, C., Chechin, D., and Cotten, A. (2008). Diffusion tensor imaging and tractography of the median nerve in carpal tunnel syndrome: Preliminary results. Eur. Radiol. 18, 2283–2291. doi:10.1007/s00330-008-0971-4
Koo, T. K., and Li, M. Y. (2016). A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med. 15, 155–163. doi:10.1016/j.jcm.2016.02.012
Kronlage, M., Schwehr, V., Schwarz, D., Godel, T., Uhlmann, L., Heiland, S., et al. (2018). Peripheral nerve diffusion tensor imaging (DTI): Normal values and demographic determinants in a cohort of 60 healthy individuals. Eur. Radiol. 28, 1801–1808. doi:10.1007/s00330-017-5134-z
Lehmann, H. C., Zhang, J., Mori, S., and Sheikh, K. A. (2010). Diffusion tensor imaging to assess axonal regeneration in peripheral nerves. Exp. Neurol. 223, 238–244. doi:10.1016/j.expneurol.2009.10.012
Lohr, D., Terekhov, M., Veit, F., and Schreiber, L. M. (2020). Longitudinal assessment of tissue properties and cardiac diffusion metrics of the ex vivo porcine heart at 7 T: Impact of continuous tissue fixation using formalin. NMR Biomed. 33, e4298. doi:10.1002/nbm.4298
Manzanera Esteve, I. v., Farinas, A. F., Pollins, A. C., Nussenbaum, M. E., Cardwell, N. L., Kahn, H., et al. (2021). Noninvasive diffusion MRI to determine the severity of peripheral nerve injury. Magn. Reson Imaging 83, 96–106. doi:10.1016/j.mri.2021.08.006
Manzanera Esteve, I. v., Farinas, A. F., Pollins, A. C., Nussenbaum, M. E., Cardwell, N. L., Kang, H., et al. (2019). Probabilistic assessment of nerve regeneration with diffusion MRI in rat models of peripheral nerve trauma. Sci. Rep. 9, 19686. doi:10.1038/s41598-019-56215-2
Mathys, C., Aissa, J., Zu Hörste, G. M., Reichelt, D. C., Antoch, G., Turowski, B., et al. (2013). Peripheral neuropathy: Assessment of proximal nerve integrity by diffusion tensor imaging. Muscle Nerve 48, 889–896. doi:10.1002/mus.23855
Möller, I., Miguel, M., Bong, A. D., Zaottini, F., and Martinoli, C. (2018). The peripheral nerves: Update on ultrasound and magnetic resonance imaging. Clin. Exp. Rheumatol. 114, 145–158.
Morisaki, S., Kawai, Y., Umeda, M., Nishi, M., Oda, R., Fujiwara, H., et al. (2011). In vivo assessment of peripheral nerve regeneration by diffusion tensor imaging. J. Magnetic Reson. Imaging 33, 535–542. doi:10.1002/jmri.22442
Peltonen, S., Alanne, M., and Peltonen, J. (2013). Barriers of the peripheral nerve. Tissue Barriers 1, e24956. doi:10.4161/tisb.24956
Pridmore, M. D., Glassman, G. E., Pollins, A. C., Manzanera Esteve, I. v., Drolet, B. C., Weikert, D. R., et al. (2021). Initial findings in traumatic peripheral nerve injury and repair with diffusion tensor imaging. Ann. Clin. Transl. Neurol. 8, 332–347. doi:10.1002/acn3.51270
Rangavajla, G., Mokarram, N., Masoodzadehgan, N., Pai, S. B., and Bellamkonda, R. v. (2014). Non-invasive imaging of peripheral nerves. Cells Tissues Organs 200, 69–77. doi:10.1159/000369451
Reina, M. A., Sala-Blanch, X., Arriazu, R., and Machés, F. (2015). “Microscopic morphology and ultrastructure of human peripheral nerves,” in Nerves and nerve injuries (Amsterdam, Netherlands: Elsevier), 91–106. doi:10.1016/B978-0-12-410390-0.00007-X
Riegler, G., Drlicek, G., Kronnerwetter, C., Heule, R., Bieri, O., Bodner, G., et al. (2016). High-resolution axonal bundle (fascicle) assessment and triple-echo steady-state T2 mapping of the median nerve at 7 T: preliminary experience. Invest. Radiol. 51, 529–535. doi:10.1097/RLI.0000000000000265
Rojoa, D., Raheman, F., Rassam, J., and Wade, R. G. (2021). Meta-analysis of the normal diffusion tensor imaging values of the median nerve and how they change in carpal tunnel syndrome. Sci. Rep. 11, 20935. doi:10.1038/s41598-021-00353-z
Snoj, Ž., Serša, I., Matičič, U., Plut, D., Cvetko, E., and Omejec, G. (2022). Median and ulnar nerve fascicle imaging using MR microscopy and high-resolution ultrasound. J. Neuroimaging 32, 420–429. doi:10.1111/jon.12982
Stein, D., Neufeld, A., Pasternak, O., Graif, M., Patish, H., Schwimmer, E., et al. (2009). Diffusion tensor imaging of the median nerve in healthy and carpal tunnel syndrome subjects. J. Magnetic Reson. Imaging 29, 657–662. doi:10.1002/jmri.21553
Sunderland, S. (1945). The intraneural topography of the radial, median and ulnar nerves. Brain 68, 243–299. doi:10.1093/brain/68.4.243
Sunderland, S. (1978). “The median nerve. Anatomical features,” in Nerves and nerve injuries (Sunderland: Kenhub), 737–779.
Tagliafico, A. S., and Tagliafico, G. (2014). Fascicular ratio: A new parameter to evaluate peripheral nerve pathology on magnetic resonance imaging: A feasibility study on a 3T MRI system. Medicine 93, e68. doi:10.1097/MD.0000000000000068
Takagi, T., Nakamura, M., Yamada, M., Hikishima, K., Momoshima, S., Fujiyoshi, K., et al. (2009). Visualization of peripheral nerve degeneration and regeneration: Monitoring with diffusion tensor tractography. Neuroimage 44, 884–892. doi:10.1016/j.neuroimage.2008.09.022
Tanitame, K., Iwakado, Y., Akiyama, Y., Ueno, H., Ochi, K., Otani, K., et al. (2012). Effect of age on the fractional anisotropy (FA) value of peripheral nerves and clinical significance of the age-corrected FA value for evaluating polyneuropathies. Neuroradiology 54, 815–821. doi:10.1007/s00234-011-0981-9
Wang, L., Yan, Z., Xiang, P., Yan, L., and Zhang, Z. (2023). MR microneurography of human peripheral fascicles using a clinical 3T MR scanner. J. Neuroradiol. [Epub ahead of print]. doi:10.1016/j.neurad.2022.11.004
Xing, D., Papadakis, N. G., Huang, C. L. H., Lee, V. M., Adrian Carpenter, T., and Hall, L. D. (1997). Optimised diffusion-weighting for measurement of apparent diffusion coefficient (ADC) in human brain. Magn. Reson Imaging 15, 771–784. doi:10.1016/S0730-725X(97)00037-4
Yao, L., and Gai, N. (2009). Median nerve cross-sectional area and MRI diffusion characteristics: Normative values at the carpal tunnel. Skelet. Radiol. 38, 355–361. doi:10.1007/s00256-008-0626-1
Yoon, D., Biswal, S., Rutt, B., Lutz, A., and Hargreaves, B. (2018). Feasibility of 7T MRI for imaging fascicular structures of peripheral nerves. Muscle Nerve 57, 494–498. doi:10.1002/mus.26035
Keywords: magnetic resonance microscopy, peripheral nerve anatomy, fractional anisotropy, mean diffusivity, eigenvectors, diffusion tensor imaging
Citation: Pušnik L, Serša I, Umek N, Cvetko E and Snoj Ž (2023) Correlation between diffusion tensor indices and fascicular morphometric parameters of peripheral nerve. Front. Physiol. 14:1070227. doi: 10.3389/fphys.2023.1070227
Received: 03 November 2022; Accepted: 03 February 2023;
Published: 23 February 2023.
Edited by:
Dan Benjamini, National Institute on Aging (NIH), United StatesReviewed by:
Richard Dortch, Barrow Neurological Institute (BNI), United StatesIvo Dumic-Cule, Clinical Hospital Dubrava, Croatia
Patryk Filipiak, New York University, United States
Copyright © 2023 Pušnik, Serša, Umek, Cvetko and Snoj. 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: Luka Pušnik, luka.pusnik7@gmail.com