Abstract
BACKGROUND AND PURPOSE: In spastic paraplegia type 5, spinal cord atrophy and white matter signal abnormalities in the brain are the main MR imaging alterations. However, the specific mechanism remains unclear. We explored the microstructural changes occurring in spastic paraplegia type 5 and assessed the relation between MR imaging and clinical data.
MATERIALS AND METHODS: Seventeen patients with spastic paraplegia type 5 and 17 healthy controls were scanned with DTI and T1 mapping on a 3T MR imaging scanner. Fractional anisotropy, mean diffusivity, radial diffusivity, axial diffusivity, and T1 values were obtained using Tract-Based Spatial Statistics and the Spinal Cord Toolbox. Neurofilament light and myelin basic protein in the CSF were measured. The differences in MR imaging and biochemical data between patients with spastic paraplegia type 5 and healthy controls were compared using the Student t test.
RESULTS: A widespread reduction of fractional anisotropy values and an elevation of mean diffusivity, T1, and radial diffusivity values were found in most cervical, T4, and T5 spinal cords; corona radiata; optic radiations; and internal capsules in spastic paraplegia type 5. A variation in axial diffusivity values was shown only in C2, C6, and the corona radiata but not in the gray matter. The levels of neurofilament light and myelin basic protein were higher in those with spastic paraplegia type 5 than in healthy controls (myelin basic protein, 3507 [SD, 2291] versus 127 [SD, 219] pg/mL; neurofilament light, 617 [SD, 207] versus 265 [SD, 187] pg/mL; P < .001). No correlation was found between the clinical data and MR imaging–derived measures.
CONCLUSIONS: Multiparametric MR imaging and biochemical indicators demonstrated that demyelination (mainly) and axonal loss led to the white matter integrity loss without gray matter injury in spastic paraplegia type 5.
ABBREVIATIONS:
- AD
- axial diffusivity
- FA
- fractional anisotropy
- HC
- healthy control
- HSP
- hereditary spastic paraplegia
- MBP
- myelin basic protein
- MD
- mean diffusivity
- NFL
- neurofilament light
- RD
- radial diffusivity
- SPG5
- spastic paraplegia type 5
- SPRS
- Spastic Paraplegia Rating Scale
Hereditary spastic paraplegia (HSP) is a neurodegenerative disorder characterized by retrograde axonal degeneration of the corticospinal tracts.1 Spastic paraplegia type 5 (SPG5) is a rare subtype of HSP for which treatment with cholesterol-lowering drugs can be attempted,2 but the target of drug action remains uncertain. Recent evidence has indicated that SPG5 is caused by a recessive mutation in the oxysterol-7α-hydroxylase gene CYP7B1, which leads to the accumulation of neurotoxic oxysterols, especially 27-hydroxycholesterol (27-OHC), which was found to impair the viability of human cortical neurons in SPG5.2⇓-4 Given that axons in the white matter are known to exert essential functions in lipid transport,5 we hypothesized that patients with SPG5 might have white matter integrity loss. An improved understanding of white matter integrity loss patterns (ie, demyelination and/or dying-back axonopathy) would be of substantial clinical relevance for enabling personalized treatment approaches for patients with SPG5.6,7
DTI is a robust method to detect in vivo white matter abnormalities. Specifically, radial diffusivity (RD) values can provide a quantitative evaluation of myelin sheaths, while axial diffusivity (AD) values support the assessment of axon status.8,9 Accordingly, DTI has been used to evaluate common subtypes of HSP, such as SPG4 and SPG11.10⇓⇓⇓-14 However, to date, no study has reported using DTI to assess both the brain and spinal cord of SPG5. One challenge with the application of DTI to HSPs is that physiologic motion and susceptibility artifacts in the spinal cord can challenge DTI quality. T1 mapping, which is relatively less susceptible to artifacts, can reflect myelin conditions and could potentially complement DTI to explore the myelin injury in HSP; the working principle of T1 mapping exploits the relatively faster relaxation of myelin water compared with nonmyelin water.15 Moreover, neurofilament light (NFL) and myelin basic protein (MBP) in the CSF are promising biomarkers that can accurately reflect neuroaxonal loss and myelin damage, respectively.16,17
In this study, we performed multiparametric structural MR imaging, including DTI and T1 mapping, to investigate the brain and spinal cord microstructural alterations in 17 patients with genetically confirmed SPG5 and conducted biochemical measurements to verify observed trends in myelin and axonal damage. Subsequently, we explored the correlations of these indicators with clinical manifestations.
MATERIALS AND METHODS
Subject Selection
Participants were prospectively recruited into the registered cohort spastic paraplegia study (National Clinical Trials: 04006418) from 2019 to 2020. This case-control parallel cohort study received approval from the local ethical committees (The First Affiliated Hospital of Fujian Medical University). All the participants agreed to participate and signed an informed consent form.
Seventeen patients with SPG5 with genetic confirmation were included (Online Supplemental Data), and 17 age- and sex-matched healthy controls (HCs) were also recruited for imaging assessment. Eleven inpatients without neurologic diseases requiring a subarachnoid block in anesthesia served as controls for CSF biochemical evaluation (Online Supplemental Data).
Clinical and Biochemical Measurements
All participants completed neurologic examinations and CSF collection on the same day as the MR imaging examination. Disease severity was quantified with the Spastic Paraplegia Rating Scale (SPRS).18 CSF NFL was measured using Simoa NF-light Advantage Kits (Quanterix) on a Simoa HD-1 Analyzer instrument. MBP was measured by a Human MBP DuoSet ELISA (R&D Systems).
MR Imaging Acquisition
All participants underwent an examination on a 3T MR scanner (Magnetom Skyra; Siemens) equipped with a 20-channel head-neck coil and a 24-channel spine-array coil. The following sequences were performed: 1) axial T2WI, T2-FLAIR, and sagittal 3D T1WI covering the whole brain; 2) sagittal 3D-T2WI and 3D-T1WI covering the cervical and thoracic spinal cord; 3) DTI with a single-shot echo-planar imaging sequence for both brain and spinal cord (32 and 20 different gradient directions, respectively); and 4) a T1-mapping sequence based on B1 inhomogeneity–corrected variable flip angle methods for the spinal cord. DTI and T1 mapping consisted of 2 axial slabs in the spinal cord: cervical (C2–C7) and thoracic (T1–T5) slabs. The parameters are listed in the Online Supplemental Data.
Image Processing
MR Imaging Analysis.
DWI preprocessing was performed on the basis of FSL (http://fsl.fmrib.ox.ac.uk/fsl) and the Diffusion Toolkit (DTK; http://www.trackvis.org/dtk/). The 3D T1WIs were preprocessed using the CAT12 toolbox of Statistical Parametric Mapping, Version 12 (SPM; http://www.fil.ion.ucl.ac.uk/spm/software/spm12). Spinal cord images were analyzed on the basis of the open-source software Spinal Cord Toolbox (Version 4.01; https://spinalcordtoolbox.com/) (Online Supplemental Data).19
Brain Image Analysis.
The DTI analysis of the brain was performed with Tract-Based Spatial Statistics (TBSS; http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/TBSS). Voxel-based whole-brain white matter measures (DTI metrics) were assessed with TBSS using FSL 5.0 (http://www.fmrib.ox.ac.uk/fsl).
Voxel-based morphometry analysis was performed to analyze the alterations in the gray matter using the CAT12 toolbox.
Spinal Cord Image Analysis.
Atlas-based analyses of spinal cord DTI and T1-mapping data were performed. The analysis pipeline is shown in the Online Supplemental Data. The spinal cord segments that met the quality control requirements are summarized in the Online Supplemental Data.
Statistical Analysis.
Statistical analyses were performed with SPSS (Version 26; IBM). The sample size was determined by the number of eligible patients willing to participate without prospective calculation. Demographic, clinical, and biochemical variables were compared between SPG5 and HCs groups by means of Student t and Pearson χ2 tests. Bonferroni adjustment was used to correct for multiple comparisons. Correlations were analyzed, counting age and sex as nuisance covariates.
TBSS statistical analysis was performed using a nonparametric permutation inference tool, FSL Randomize (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Randomise/UserGuide). DTI metrics were compared between 2 groups using a 2-sample independent t test. A voxelwise 2-sample t test was conducted to detect the differences in the gray matter between patients with SPG5 and HCs. The total intracranial volume of each participant was entered as a covariate of no interest. Statistical inferences were made at P < .001 (uncorrected) and for clusters larger than κ = 20.
RESULTS
Demographic and Clinical Features
The demographic and clinical data of all participants are shown in the Table. According to the Harding criteria,20 16 patients had pure HSP and 1 patient had complicated HSP. All 17 patients presented with lower limb spasticity and weakness and dorsal column dysfunction. Sixteen patients had a severely reduced or nonexistent vibration sense in the lower limbs.
Characteristic | HC1 | HC2 | SPG5 | P Value1 | P Value2 |
---|---|---|---|---|---|
Participantsb | 17 | 11 | 17 | >.99c | <.01c |
Sexb | .99c | .98c | |||
Men | 11 | 8 | 11 | ||
Women | 6 | 3 | 6 | ||
Age (yr)e | 29 (SD, 10) (14–49) | 39 (SD, 18) (11–69) | 30 (SD, 10) (13–49) | .93d | .13d |
SPRSe | 16 (SD, 9) (2–38) | ||||
Disease duration (yr)e | 18 (SD, 10) (6–37) |
↵a 1 = Healthy controls recruited for imaging assessment, 2 = hospitalized patients requiring subarachnoid block in anesthesia without neurologic disease enrolled as controls for CSF biochemical assessment.
↵b Data are numbers of participants.
↵c Pearson χ2 test.
↵d Student t test.
↵e Data are mean (SD). Data in parentheses are ranges.
The CSF MBP levels and CSF NFL levels in patients with SPG5 (CSF MBP, 3507 [SD, 2291] versus 127 [SD, 219] pg/mL, P < .001; CSF NFL, 617 [SD, 207] versus 265 [SD, 187] pg/mL; P < .001) were both higher than those in the control group. However, the increased degree in MBP was much higher than that in NFL in patients with SPG5 (compared with controls, the MBP increased by 27.61 [3507/127] folds on average, while the NFL increased only by 2.33 [617/265] in patients with SPG5). Conventional MR imaging showed no apparent signal abnormalities in the brain or spinal cord of those with SPG5 or HCs.
Brain MR Imaging Analysis
TBSS Analysis.
Compared with HCs, an extensive and bilateral symmetric reduction of fractional anisotropy (FA) and an elevation of mean diffusivity (MD) and RD were found in the cerebral peduncles, optic radiation, internal capsules, corona radiata, corpus callosum, posterior thalamic radiation, cingulum, and superior and inferior longitudinal fasciculi of patients with SPG5. The areas of elevated RD and MD closely mirrored the areas featuring reduced FA. However, an elevation of AD reached a statistically significant level only in the corona radiata (Fig 1).
Voxel-Based Morphometry Analysis.
There were no statistically significant differences in gray matter volume between patients with SPG5 and HCs with a less conservative threshold (P < .001, uncorrected).
Spinal Cord Analysis.
Compared with HCs, patients with SPG5 had reduced FA values and elevated RD, MD, and T1 values at the C2–T5 segments in the white matter, dorsal columns, and bilateral lateral corticospinal tract. Significant differences of FA, RD, and T1 values in the white matter regions between patients with SPG5 and HCs were observed in most cervical spinal cords but not in the thoracic segments, except T4 and T5. In contrast, no evident differences of AD values were observed in most spinal cord segments, and marked reduced AD values were found only at specific segments (eg, C2 and C6) (Fig 2 and Online Supplemental Data). Notably, no significant differences in DTI metrics or T1 values were detected in spinal cord gray matter between patients with SPG5 and HCs (Online Supplemental Data).
Correlation Analysis.
No correlations were found between DTI parameters and disease severity (SPRS scores) or disease duration using TBSS. DTI or T1 values at all spinal cord segments in patients with SPG5 did not correlate with the clinical data (P > .05, after adjustment for multiple comparisons).
DISCUSSION
This prospective study investigated microstructural alterations occurring in patients with SPG5. First, the FA values of the brain and spinal cord white matter were reduced diffusely, indicative of widespread white matter injury in SPG5. Second, the RD, AD, T1 values, and biochemical indicators (NFL and MBP) demonstrated that the white matter injury pattern of SPG5 was demyelination (mainly). Third, there were no correlations between the degree of white matter injury and the clinical data in SPG5.
Although case reports showed scattered high signal intensity in the cerebral white matter,21,22 our patients with SPG5 demonstrated a normal appearance on conventional MR imaging (16 of 17 were pure HSP), like the ‘‘pure’’ phenotype HSP that has been reported without white matter abnormalities.23,24 However, reduced FA in the brain and spinal cord revealed white matter integrity loss in patients with SPG5. Meanwhile, corticospinal tract and dorsal column involvement was consistent with the clinical symptoms: lower limb spasticity, weakness, and dorsal column dysfunction. In addition, we noted that the optic radiations were also affected in SPG5, though there were no clinical manifestations. Therefore, SPG5-specific widespread white matter defects were possibly related to the accumulation of 27-OHC22 rather than a genetic cause of damage to a specific site. On the basis of MR imaging, we have demonstrated the clinical feasibility of treating patients with SPG5 by lowering the pathologically elevated levels of oxysterols.25
Reduced FA could be caused by axonal loss and/or myelin breakdown. RD and AD values derived from DTI can, respectively, indicate myelin sheath and axon status.26 Our findings showed that RD increased in both the spinal cord and brain; AD values varied in only several cervical spinal cord segments and brain regions. These findings suggest that the white matter damage in SPG5 was more related to myelin destruction than to axonal loss. Similar to our results, an MR spectroscopy study showed an elevated mIns/Cr ratio in the occipitoparietal region in 2 patients, also indicating demyelination in SPG5.27 Moreover, considering the limited stability of spinal cord DTI, we used the T1-mapping sequence to detect the possible myelin injury.28 A prolonged T1-relaxation time in the spinal cord may indicate demyelination changes in SPG5 to some extent.
Furthermore, the NFL and MBP levels in the CSF were also analyzed to explore the myelin and axonal degeneration in SPG5. NFL is a component of the cytoskeleton of myelinated axons and can reflect axonal injury,16 while the MBP holds the layers of myelin together, so quantifying MBP supports the assessment of myelin destruction.17 Significant differences of both NFL and MBP were found between HCs and those with SPG5. Additionally, the mean MBP level in the CSF (3507 pg/mL) in our study was much higher than that in neuromyelitis optica (706 pg/mL) or multiple sclerosis (106 pg/mL);29 the mean NFL level (617 pg/mL) was lower than that in Guillain-Barré syndrome (1308 pg/mL)30 or multiple sclerosis (884 pg/mL).31
Thus, our DTI evidence and biochemical measurements together support the cautious deduction that the copresence of demyelination (mainly) and axonal loss resulted in white matter degeneration in SPG5. On the other hand, no gray matter abnormalities were found in the brain or spinal cord, in good agreement with the results of a previous MR spectroscopy study: no increased or reduced Cho/Cr ratio in the gray matter of the occipital region.27 In brief, demyelination in both white matter and preserved gray matter in patients with SPG5 was consistent with the characteristics of pure HSP as reported in the literature12,14,32
We found no correlations between clinical data (disease duration or SPRS scores) and MR imaging metrics. As previously reported,11 the ceiling effect may be a good explanation. Patients with SPG5 in our study were disabled and had a long disease duration (the mean SPRS of our patients was 16, the mean disease duration was 18 years); even though neurodegeneration was still taking place, the correlation analyses failed to capture clinical deterioration.33 Meanwhile, considering the slowly progressive nature of the disease, white matter injury may be compensated by neuroaxonal redundancy or adaptative neuroplasticity, which led to a discrepancy between the functional evaluation of SPG5 patients and DTI structural assessment in SPG5.34 Additionally, cholesterol is a vital component of myelin; we assumed that abnormalities in cholesterol metabolism could influence the early development of myelin before the degeneration.
Our study had several limitations. First, the interpretations of AD, RD, and T1 values are still contested.35 Diffusional kurtosis imaging and neurite orientation dispersion and density imaging may be the optimal MR imaging techniques to strengthen the findings of this study.35 Second, lumbar puncture is an invasive examination; the HCs who underwent the MR imaging examinations were not the same ones who underwent the lumbar puncture to get the biochemical measurements. Third, the data reliability of thoracic spinal cord DTI may be a problem in this study, which is expected to be solved by more sequence improvement.
CONCLUSIONS
Multiparametric structural MR imaging and biochemical measurements revealed that the neurodegeneration progression of SPG5 was related to white matter integrity loss (demyelination mainly) rather than neuronopathy. Our identification of the widespread myelin damage can support the development of SPG5 treatment regimens in the future.
Acknowledgments
We thank all patients with SPG5 and healthy controls participating in this study.
Footnotes
Ying Liu and Zhixian Ye contributed equally to this study.
This work was supported by a Startup Fund for Scientific Research, Fujian Medical University (2017XQ1016).
Disclosure forms provided by the authors are available with the full text and PDF of this article at www.ajnr.org.
References
- Received July 16, 2021.
- Accepted after revision September 10, 2021.
- © 2022 by American Journal of Neuroradiology