Abstract
Background
Indirect traumatic optic neuropathy (ITON) is a devastating cause of permanent visual loss. Axonal degeneration, the characteristic pathological change of ITON, cannot be assessed by conventional imaging. Diffusion tensor imaging (DTI) has been widely used as a sensitive non-invasive imaging technique to obtain information on axonal integrity.
Purpose
To study time-dependent changes in ITON patients with DTI and to provide imaging evidence for clinical diagnosis and therapy.
Material and Methods
We enrolled 28 subjects with unilateral ITON who underwent detailed ocular examinations and magnetic resonance imaging (MRI) examinations. The differences between injured optic nerve (ON) and contralateral ON were tested. The patients were divided into three groups based on time (from injury to examination). Groups 1, 2, and 3 corresponded to the time taken from injury to examination: <7 days, 7–30 days, and >30 days, respectively. DT magnetic resonance imaging (MRI) differences among the groups were compared, including the relationship between diffusion parameters and mean thicknesses of the peripapillary retinal nerve fiber layer (RNFL) and ganglion cell complex (GCC) in the macular area.
Results
Compared with contralateral ON, we observed reduced fractional anisotropy (FA) of injured nerves in group 2. Reduced FA and decreased axial diffusivity (λ//) and increased radial diffusivity (λ┴) and mean diffusivity (MD) of injured nerves were observed in group 3. The mean FA value of injured nerves showed a progressive decreasing trend, and mean λ┴ value exhibited a progressive increasing trend. For injured eyes, the MD and λ┴ increases strongly correlated with the decreased mean thicknesses of RNFL and GCC. Conversely, FA was significantly associated with mean RNFL thickness.
Conclusion
DT-MRI parameters could be useful biomarkers in detecting ON changes in ITON patients.
Introduction
Traumatic optic neuropathy is a devastating cause of permanent visual loss following contusive injuries to the head, particularly to the forehead (1), and indirect traumatic optic neuropathy (ITON) is the most common form (2,3). ITON is defined as traumatic loss of vision that occurs without external or initial ophthalmoscopic evidence of injury to the eye or its nerve (4). Although the degree of visual loss after indirect traumatic optic neuropathy can be quite variable, approximately 50% of patients maintain “light perception” or “no light perception” vision (5).
No treatment, megadose corticosteroids, optic canal decompression surgery, or a combination of these are the present management options for ITON, but the optimal treatment remains controversial mainly because the lack of histologic data in the clinical examination (3,6,7). Short TI inversion recovery (STIR) magnetic resonance imaging (MRI) sequences can detect the high signal of swelling ON in early ITON, but this change is somewhat subjective and can be mild at the time of trauma or within the first few days afterwards. Apparent ON atrophy and ON sheath dilation would be detected more than 2 months post trauma (8), so conventional imaging cannot sensitively assess ITON. Furthermore, conventional imaging fails to accurately identify the underlying pathology. Therefore, the development of a non-invasive imaging technique to detect early changes and to evaluate ON axonal degeneration is important. Diffusion tensor imaging (DTI) is a sensitive technology for investigating central nervous system (CNS) white matter (WM) diseases and injuries (9–11). DTI utilizes anisotropic diffusion displacement-probability ellipsoids that are characterized by three eigenvalues (λ1, λ2, and λ3) and three eigenvectors. The average of the three eigenvalues is referred to as mean diffusivity (MD), MD = (λ1 + λ2 + λ3)/3. In CNS WM, the largest eigenvalue λ1 represents water diffusivity parallel to axonal fibers. It is referred to as λ//, or axial diffusivity. The water diffusivities perpendicular to the axonal fibers, λ2 and λ3, are averaged and referred to as λ┴ = (λ2 + λ3)/2, or radial diffusivity. Fractional anis-otropy (FA) was defined by the following equation:
Given that the axonal structural integrity reflects the random displacement of water molecules (i.e. anisotropic diffusion), it is feasible to utilize DTI as a means to obtain information about axonal state (12). Mac Donald et al. (13), in a mouse model, showed that DTI was highly sensitive to experimental pericontusional white matter injury throughout the full range of assessed time points, and DTI signal changes could determine the approximate time since injury. ITON is a closed injury produced by force imparted to the skull and transmitted into the ON, and this could cause ON damage and subsequent axonal degeneration (2); therefore, we hypothesized that ITON could be investigated with DTI.
The aims of our study were to observe time-dependent ON changes in ITON patients with DT-MRI, which could reflect the axonal state and myelin degradation, and to investigate the relationships among DT-MRI parameters and the mean thicknesses of peripapillary retinal nerve fiber layer (RNFL) and ganglion cell complex (GCC) in macular areas.
Material and Methods
Subjects
Demographic and clinical characteristics of unilateral ITON patients.
HM, hand movement; LP, light perception; NLP, no light perception.
Each patient underwent high-resolution orbital computed tomography (CT) with a 64-slice multidetector CT (Philips Brilliance CT, Philips Medical Systems, Best, The Netherlands) to rule out the presence of optic canal fractures or bone fragments in the orbital apex. The patients underwent detailed ocular examinations, including visual acuity, slit-lamp biomicroscopy, Goldmann applanation tonometry, and fundus examination. Spectral domain optical coherence tomography (SD-OCT) examination was performed with the RTVue (Optovue Inc., Fremont, CA, USA). RTVue was used to generate images of ocular microstructures and employed a scanning laser diode with a wavelength of 840 ± 10 nm, which can perform 26,000 A scans/s with a depth resolution of 5 µm. In our study, the mean thicknesses of RNFL and GCC in the macular areas of all patients were performed on the same day as the MRI examination. RTVue OCT images were obtained by blinded examiners. Only good quality images, as defined by a signal strength index of ≥30 were used for the analysis.
MRI
All patients underwent MRI examination using a 1.5-Tesla magnet (GE Healthcare, Waukesha, WI, USA) scanner with an 8-channel head coil. Subjects were asked to close their eyes to minimize any effects of deliberate eye movement during the acquisition. Head motion was minimized by restraining foam pads.
Each subject was scanned using a high-resolution T2-weighted (T2W) fluid-attenuated inversion recovery sequences (FLAIR) sequence (TR, 9000 ms; TE, 120 ms; TI, 2125 ms; slices, 32; slice thickness, 4 mm; slice gap, 0.8 mm; field of view [FOV], 240 × 210 mm2; matrix size, 256 × 222) to detect any brain abnormalities, and any patient with abnormalities was excluded.
Axial T2W images of ON (TR, 3800 ms; TE, 102 ms; slices, 12; slice thickness, 3.5 mm; slice gap, 0.3 mm; FOV, 180 × 180 mm2; matrix size, 288 × 224) was obtained as the location sequence for the DTI.
Coronal ON DTI images were obtained using spin-echo echo-planar-imaging sequences with parallel acquisition. Coronal-oblique slices were set orthogonal to the nerves. The imaging coverage was from the optic papilla to the orbital apex of the ON (Fig. 1).
Position of optic nerve slices on an axial view. There are eight slices from the anterior part (adjacent to the optic papilla) to the posterior part (near the orbital apex).
We used the following acquisition parameters for ON DTI: one b0 and six non-collinear gradient directions with b = 600 s/mm2; slices, 8; thickness, 5 mm; gap, 0; FOV, 22 × 22 cm2; acquisition matrix size, 128 × 128; NEX, 16; voxel size, 0.86 × 0.86 × 5 mm3. The signal to noise ratio (SNR) was set at 35–40 in our study. The eddy current distortions and motion artifacts in the DT-MRI data were corrected using FMRIBs diffusion toolbox (FSL, Oxford, UK, http://www.fmrib.ox.ac.uk/fsl/fdt/index.html) (14,15).
The MD, FA, and eigenvalues (λ1, λ2, and λ3) of ON were calculated voxel-by-voxel using DTI Studio (Version 3.0, MRI Studio Software, Johns Hopkins University, Baltimore, MD, USA). All diffusion-weighted (DW) images were visually inspected for apparent artifacts due to head motion, and no significant image deterioration was found. Subsequently, the eddy current was eliminated by realigning all DW and non-DW images. The regions of interest (ROIs) were defined manually on the DW image consisting of two square 2 × 2 voxels. The ROIs were placed in the center of the nerve (Fig. 2) (16–18). The six elements of the diffusion tensor were calculated for each pixel using multivariate linear fitting. After tensor diagonalization, three eigenvalues and eigenvectors were obtained, and FA maps were calculated. The eigenvector associated with the largest eigenvalue was used as an indicator for fiber orientation. The color FA map images were also calculated from DW images. The intraorbital fourth layer of the nerve (about 2.0 cm behind the optic papilla) was used for the following analysis (19).
ROI selection in the fourth slice of the optic nerve. (a) Diffusion-weighted image. (b) Fractional anisotropy map. The yellow arrow indicates the optic nerve.
We certify that our research adhered to the tenets of the Declaration of Helsinki. All institutional and governmental regulations concerning the ethical use of human subjects were followed during this study, and informed consent was obtained from each subject.
Statistical analysis
Statistical analysis was performed using SPSS v13.0 (SPSS Inc., Chicago, IL, USA). Differences between injured ON and contralateral ON were explored using paired t-tests. One-way analysis of variance (ANOVA) test was used to compare the DT-MRI parameters among three subgroups at different time points. The DT-MRI parameters of two examinations of the same one patient were compared using paired t-tests. The relationships among the diffusion parameters of ON and mean thicknesses of peripapillary RNFL and GCC in the macular area were investigated using Pearson correlation analysis. Values of P < 0.05 were considered statistically significant.
Results
Diffusion parameters
Diffusion parameters of injured nerves and contralateral nerves by ITON group.
Axial, radial, and mean diffusivities (λ║, λ┴, MD) are given in µm/ms. Fractional anisotropy is without units. All values are mean ± sd. ITON=Indirect traumatic optic neuropathy.
Group 1, the time (from injury to examination) was <7 days. Group 2, the time (from injury to examination) was 7–30 days. Group 3, the time (from injury to examination) was >30 days.
When DT-MRI parameters of three subgroups with different time points were compared, mean FA showed a progressive decreasing trend (F = 8.99, P = 0.001; Fig. 3, Table 3). Compared with Group 1, the mean MD values of Groups 2 and 3 both increased (F = 5.55, P = 0.008; Fig. 3, Table 3). The mean value of λ┴ showed a progressive increasing trend (F = 8.53, P = 0.001; Fig. 3, Table 3). The injured ON showed an initial decrease in mean λ// in the first week after injury compared with the contralateral ON (t = 1.72, P = 0.120; Table 2). The mean λ// value in injured ON showed a trend towards increasing in Group 2 and decreased again in Group 3 (F = 1.00, P = 0.379; Fig. 3, Table 3).
Dynamic diffusion parameters changes in injured optic nerve. FA showed a progressive decreasing trend. MD was normal within 7 days after injury and elevated thereafter. λ┴ showed a progressive increase trend. λ// decreased the first week following injury and then elevated from 7 days to 1 month after surgery before decreasing again at more than 1 month after injury. One-way ANOVA analysis of diffusion parameters of patients with ITON by group. Axial, radial, and mean diffusivities (λ║, λ┴, MD) are given in µm/ms. Fractional anisotropy is without units. All values are mean ± sd. P < 0.05 in Dunnett’s t multiple comparisons when compared with first group. ITON, indirect traumatic optic neuropathy.
DT-MR parameters comparison between two time points (n = 13).
Axial, radial, and mean diffusivities (λ║, λ┴, MD) are given in µm/ms. Fractional anisotropy is without units. All values are mean ± sd.
Mean RNFL and GCC thicknesses
Mean RNFL and GCC thickness of injured eyes and contralateral eyes at different time stage in ITON.
GCC, ganglion cell complex; Group 1, the time (from injury to examination) was < 7 days; Group 2, the time (from injury to examination) was 7–30 days; Group 3, the time (from injury to examination) was >30 days; ITON, indirect traumatic optic neuropathy; RNFL, retinal nerve fiber layer.
Correlation between diffusion parameters and RNFL and GCC thicknesses
Correlation matrix of diffusion parameters, RNFL, and GCC.
GCC, ganglion cell complex; RNFL, retinal nerve fiber layer.
Discussion
A previous study demonstrated that force applied to the frontal bone could transfer to the intracanalicular segment of the ON (2). Because ON vasculature in the canal is pial, compression and contusion of the nerve produces a compartment syndrome whereby swelling exacerbates ischemia. The ON also suffers shearing injury under the fixed edge of the falciform dural fold just before it enters the optic canal (3). A histopathological study of 174 patients who died after closed head injury demonstrated evidence of hemorrhage in the ON sheath (83%), in the nerve interstitium (36%), and shearing lesions and ischemic necrosis of the intracanalicular and intracranial segments of the ON (20).
Diffuse axonal injury (DAI) is a consequence of sustained acceleration and deceleration forces that can shear axons and produce microscopic changes in the brain similar to ITON. Like ITON, DAI cannot be assessed by conventional imaging. However, DTI has been used to investigate diffuse DAI, and a direct comparison between DTI-detected WM integrity changes and histologic findings in an animal model of axonal injury suggested that DTI may be a valuable imaging tool for detecting diffuse axonal injury (10,21). Therefore, we used DTI to investigate ITON patients.
Challenges associated with DT-MRI of the ON are the small diameter of the nerve, artifacts caused by eye movement, and the susceptibility effects caused by nearby sinuses. In light of these, and despite using different sequences and protocols, it is remarkable that several different groups have reported similar values in healthy controls (MD 1.0–1.3 × 10–3mm2 s–1 and FA 0.4–0.6) (16–18,22). In our study, the mean FA and MD values of 28 control eyes were 0.54 ± 0.06 and 1.54 ± 0.20 × 10–3mm2 s–1, respectively. No paper was retrieved about using DTI to investigate ITON, so we consulted many papers about DAI to explain our results.
When axonal injury occurs, centrifugal disintegration of the axonal cytoskeleton leads to sudden and rapid axon fragmentation. This stage is short-lived, lasting from several hours to days. At the same time, myelin sheaths that surround the axons become less tightly wrapped and eventually break apart and form ovoids. In the weeks following the precipitating injury, the microglia are activated and digests the axonal and myelin debris (12). Song et al. in a mouse model of retinal ischemia (23), performed serial diffusion ON measurements and showed that λ// and λ┴ could differentiate axonal from myelin damage during the course of degeneration. According to this animal study, λ// showed a significant decrease in the first few days of degeneration, which corresponded to the disintegration of the axonal microstructure, whereas myelin remained intact. Five days after the initial injury, λ┴ increased, which corresponded to myelin sheath degradation.
In Group 1, we observed that the mean λ// of injured ON was reduced compared with the contralateral ON. In Group 2, the mean FA of injured ON was reduced, and λ┴ increased relative to contralateral ON. Meanwhile, λ// of the injured ON was slightly increased when compared with Group 1 (Table 2). Histological and electron microscopy results from a mouse traumatic brain axonal injury model (13) indicated that there were three distinct temporal phases of injury in pericontusional WM. According to this animal study, the predominant pathology in the first 24 h was relatively pure axonal injury without gliosis or demyelination, and the diffusion measurements changes were that λ//↓,λ┴≈, FA↓, and MD≈/↓. In the later acute stage (4 days), axonal injury was still marked, but gliosis became more pronounced. They hypothesized that DTI in the first 4 days would show a decreased axial diffusivity attributable to extensive axonal injury. This would, in turn, cause a marked reduction in relative anisotropy and a modest reduction in MD. In the subacute stage (7 days to 1 month), axonal injury became less apparent, and myelin thinning and active demyelination by macrophages was prevalent, while gliosis remained prominent. At the subacute stage, there would be an increase in radial diffusivity attributable to demyelination and reductions in axial diffusivity as a result of persistent axonal injury (λ//≈, λ┴↑, FA↓, and MD↑) (13). Our research yielded similar results. We attributed the mean λ// reduction in Group1 (the first 7 days after injury) to axonal injury without gliosis or demyelination. This difference did not reach statistical significance, which may be due to the small sample size. In Group 2, it appears that the second stage, i.e. 7 days to 1 month after injury, was the subacute stage of ITON. The main pathologic change during this stage was demyelination.
In Group 3 (injured over 1 month), we observed significantly decreased FA and λ// values and increased MD and λ┴ values compared with contralateral ON (Table 2). Similarly, a rat spinal cord injury model showed that λ// in and surrounding a chronic spinal cord injury decreased at 12 weeks after injury. These changes may be maximal at the sites of greatest axonal loss, which would explain the correlation between decreasing λ// values and decreasing numbers of labeled axons (24). Reduction in FA at this stage is due to a reduction in parallel diffusivity and an increase in perpendicular diffusivity. Increased MD was due to a marked increase in perpendicular diffusivity, while, decreased λ// implied axonal loss. Therefore, we suggested that ITON began to enter the chronic stage after 1 month.
Although the loss of axons may be expected to increase water diffusion as it does in chronic cerebral infarcts (25,26), the λ// in traumatically injured spinal cord and brain is not elevated; rather, it is less than that observed in normal tissue (24,27). Vorisek et al. (28) provided a possible explanation when they noted that there were both astrocyte proliferation and increased chondroitin sulfate proteoglycan expression in the extracellular matrix after experimental DAI, both of which may contribute to consistently decreased water diffusion.
We found that the mean macular thickness of GCC, which is the sum of the thicknesses of the retinal nerve fiber, ganglion cell, and inner plexiform layers (29) was significantly associated with mean RNFL thickness in injured eyes. A study on glaucoma suggested a correlation between the mean GCC and RNFL thicknesses in the intact hemisphere and thought that GCC thickness could be considered to represent damage of ganglion cells and their axons (30).
Our findings demonstrate that the degree of progressive axonal dieback in injured ON (mean thicknesses of GCC and RNFL) strongly correlated with the degree of abnormality in measured DTI parameters (MD, FA, and λ┴). In a rat spinal cord injury model with greater axonal loss closer to the injury site, λ┴ increased and λ// decreased, which correlated with the greatest loss of axons (31). Therefore, it appeared that DTI parameters could reflect changes in axonal numbers.
We found that the mean FA value of injured nerves showed a progressive decreasing trend, and mean λ┴ values showed a progressive increasing trend as the time since injury increased. When following the same patient over time, mean FA decreased and mean λ┴ increased. The average thicknesses of RNFL and GCC of injured eyes were significantly thinner than those in contralateral eyes in the second and third stages. These results suggest that neuronal and axonal degeneration are difficult to reverse in these stages.
There are some limitations to this study. Histologic validation was not possible, so our discussion of histological changes is speculative. In addition, the standard MRI of patients was not complete, so we could not provide more information about nerve diameters or areas. DT-MRI sequencing of the ON needs to be improved. However, the novelty of this study is the time-dependent DTI change observed in ONs of ITON patients without other traumatic findings, which, to our knowledge, has never been done before.
In conclusion, we observed time-dependent changes of DTI parameters could reflect the progression of axonal degeneration in ITON. We suggest that ITON patients with normal conventional imaging results could benefit from this objective non-invasive test because improved assessment techniques may facilitate the development of effective therapeutic strategies, aid in the clinical management of patients, and allow more accurate prognostic statements to be made at earlier times following injury.
Footnotes
Acknowledgements
We would like to thank all the patients who participated in this study.
