- 1Department of Neurology, Tongren Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 2Shanghai Key Laboratory of Magnetic Resonance, School of Physics and Electronic Science, East China Normal University, Shanghai, China
Background: Neuropathological studies have revealed copper and iron accumulation in the deep gray matter (DGM) nuclei of patients with Wilson’s disease (WD). However, the association between metal accumulation and neurodegeneration in WD has not been well studied in vivo. The study was aimed to investigate whether metal accumulation in the DGM was associated with the structural and functional changes of DGM in neurological WD patients.
Methods: Seventeen neurological WD patients and 20 healthy controls were recruited for the study. Mean bulk susceptibility values and volumes of DGM were obtained from quantitative susceptibility mapping (QSM). Regions of interest including the head of the caudate nucleus, globus pallidus, putamen, thalamus, substantia nigra, red nucleus, and dentate nucleus were manually segmented. The susceptibility values and volumes of DGM in different groups were compared using a linear regression model. Correlations between susceptibility values and volumes of DGM and Unified Wilson’s Disease Rating Scale (UWDRS) neurological subscores were investigated.
Results: The susceptibility values of all examined DGM in WD patients were higher than those in healthy controls (P < 0.05). Volume reductions were observed in the head of the caudate nucleus, globus pallidus, putamen, thalamus, and substantia nigra of WD patients (P < 0.001). Susceptibility values were negatively correlated with the volumes of the head of the caudate nucleus (rp = −0.657, P = 0.037), putamen (rp = −0.667, P = 0.037), and thalamus (rp = −0.613, P = 0.046) in WD patients. UWDRS neurological subscores were positively correlated with the susceptibility values of all examined DGM. The susceptibility values of putamen, head of the caudate nucleus, and dentate nucleus could well predict UWDRS neurological subscores.
Conclusion: Our study provided in vivo evidence that paramagnetic metal accumulation in the DGM was associated with DGM atrophy and neurological impairment. The susceptibility of DGM could be used as a biomarker to assess the severity of neurodegeneration in WD.
Introduction
Wilson’s disease (WD), also known as hepatolenticular degeneration, is a copper-overload inherited disease caused by ATP7B mutations on chromosome 13. ATP7B encodes a copper-transporting P-type ATPase ATP7B, which plays a critical role in maintaining copper homeostasis (Czlonkowska et al., 2018). The dysfunction of ATP7B causes excessive copper accumulation in the liver, and then large amount of non-ceruloplasmin-bound copper is released into the circulation, leading to secondary copper accumulation in the brain, kidney, and other organs (Bandmann et al., 2015). Liver disease and neurological symptoms are the most common clinical features of WD. Copper chelators can effectively improve the hepatic symptoms of WD patients, but they are less effective in improving neurological symptoms and even cause neurological deterioration in initial treatment (Li et al., 2016; Czlonkowska and Litwin, 2017). The pathogenesis of neurological impairment in WD is still not clear. In addition to the toxic effects of copper on the central nervous system, iron accumulation may also aggravate neural damage (Dusek et al., 2015, 2017).
Magnetic resonance imaging (MRI) is the most sensitive imaging tool in the diagnosis of WD. Abnormal signals were observed in more than 90% of the neurological WD patients (Litwin et al., 2013a). The most common brain MRI findings are symmetrical or asymmetric hyperintensities in T2-weighted images in the deep gray matter (DGM) nuclei and brain stem. Hyperintensities in T2-weighted images in the brain may reflect neuropathological changes such as edema, demyelination, and gliosis, which are largely reversible after effective anti-copper treatment (King et al., 1996; Czlonkowska et al., 2018). Brain atrophy is another common imaging feature of neurological WD patients, and both DGM and white matter are affected. Recent studies suggested that total brain volume and the volumes of white matter and gray matter were associated with the severity of neurological impairment in WD patients (Smolinski et al., 2019; Dusek et al., 2020). However, it is not clear whether this correlation still exists in single DGM nuclei. Although WD is a copper-overload disease, the correlation between copper accumulation and neuropathological changes is still controversial (Poujois et al., 2017). An increasing number of studies indicate that iron accumulation in the DGM is involved in the neuropathological changes of WD. Dusek et al. (2017) reported that increased iron content was detected in the basal ganglia of WD patients, and severe putamen (Put) pathological changes were associated with an increased number of iron-containing macrophages. In addition, hypointensities in T2/T2∗-weighted images and susceptibility-weighted images (SWI) in the DGM are frequently present in WD patients, also reflecting paramagnetic metal (mainly iron) accumulation in these regions (Skowronska et al., 2013; Dusek et al., 2017). A recent MRI- histopathological study demonstrated that R2∗ values correlated with iron concentration in the DGM of WD patients (Dusek et al., 2017). However, R2∗ mapping cannot distinguish between iron and diamagnetic substances. R2∗ mapping also depends on factors such as field strength, homogeneity of iron distribution, and suffers from blooming artifacts (Bilgic et al., 2012).
Quantitative susceptibility mapping (QSM) is a recently developed MRI post-processing technique, which can quantitatively measure tissue susceptibility and reflect the concentrations of paramagnetic metal in vivo (Wang et al., 2017; Li et al., 2018). Using QSM, several recent studies found significantly increased susceptibility in the DGM of neurological WD patients, suggesting paramagnetic metal accumulation in the DGM is a character of neuropathology (Fritzsch et al., 2014; Doganay et al., 2018; Dezortova et al., 2019). In addition, QSM performed better than R2∗ in distinguishing WD from healthy individuals and showed excellent diagnostic accuracy in the diagnosis of WD (Li et al., 2020).
Intracellular ferritin and hemosiderin are considered as the main source of the DGM susceptibility, which has been confirmed in recent MRI-histopathology studies (Langkammer et al., 2012; Hametner et al., 2018). In vitro evidences indicate that the Cu2+ compounds are slightly paramagnetic and the Cu+ compounds are diamagnetic, but the magnetic features of copper in vivo are still unclear (Valdes Hernandez Mdel et al., 2012). Considering that a large amount of copper accumulates in the brains of WD patients, increased susceptibility in the DGM may be explained by copper and iron accumulation.
Several case reports reported that metal accumulation in the DGM changed with the progression or improvement of neurological symptoms (Dusek et al., 2018; Zaino et al., 2019). It is speculated that metal accumulation in the DGM is a surrogate indicator of oxidative stress and neuronal loss (Yan et al., 2018). However, the relationship between metal accumulation and neurodegeneration in WD has not been quantitatively studied in vivo. Moreover, a quantitative biomarker for the evaluation of neurodegeneration in WD is still lacking. Validated neuroimaging biomarker will be very helpful in treatment monitoring and outcome prediction.
In this study, we hypothesized that metal accumulation in the DGM was associated with the structural and functional changes of DGM in neurological WD patients. Using the recently developed MRI technique QSM, we quantitatively measured the susceptibility values of DGM, and investigated the association between metal accumulation and DGM volume changes. Furthermore, the association between metal accumulation and the severity of neurological impairment in WD was also assessed.
Materials and Methods
Subjects
Twenty neurological WD patients were recruited from the Department of Neurology, Tongren Hospital, Shanghai Jiao Tong University School of Medicine, between December 2018 and December 2019. WD was confirmed according to the Leipzig diagnostic criteria, including low serum ceruloplasmin, presence of Kayser-Fleischer rings on slit-lamp examination, high 24-hour urine copper excretion, and genetic confirmation (Ferenci et al., 2003). All WD patients were on anti-copper treatment with D-penicillamine or dimercaptosuccinic acid (DMSA). Before MRI scans, the severity of neurological impairment in WD patients was assessed by two neurologists with consensus using Unified Wilson’s Disease Rating Scale (UWDRS) neurological subscale (27 items, 208 points) (Leinweber et al., 2008; Volpert et al., 2017). Neurological WD was defined as the presence of neurological symptoms at the time of onset or anytime during the course of the disease. Neurological WD patients aged 18 years or older were included in this study. Exclusion criteria included severe neurological impairment or severe liver decompensation, history of other neuropsychiatric disorders or liver disease, and hepatic encephalopathy. T1 hyperintensities in the lentiform nucleus of WD patients are associated with portosystemic shunting and hepatic encephalopathy (Kozic et al., 2003), which may affect regional susceptibility duo to manganese accumulation and the evaluation of neurological symptoms (Lee et al., 2018). Therefore, WD patients with T1 hyperintensities in the lentiform nucleus were excluded from the study. Twenty age and gender matched healthy controls (HC) were recruited for the study. This study was approved by the Medical Ethics Committee of Tongren Hospital, Shanghai Jiao Tong University School of Medicine. All subjects were informed about the study and signed written consent before participation.
Imaging Data Acquisition and Procession
Magnetic resonance imaging scans were carried out on a clinical 3T MRI system (Magnetom Prisma Fit, Siemens Healthcare, Erlangen, Germany) equipped with a 20-channel head coil. The susceptibility maps were acquired from a 3D spoiled multi-echo gradient-echo (GRE) sequence with the following imaging parameters: repetition time (TR) = 31 ms, the first echo time (TE1) = 4.07 ms, echo spacing (ΔTE) = 4.35 ms, numbers of echoes = 6, flip angle = 12°, field of view (FOV) = 240 mm× 200 mm, in-plane resolution = 0.83 mm× 0.83 mm, slice thickness = 0.8 mm, number of slices = 192, acquisition time = 7.22 min. A generalized auto-calibrating partially parallel acquisition (GRAPPA) with an acceleration factor of two in the right-left direction and elliptical sampling were used to reduce the acquisition time. T1-weighted, T2-weighted, and T2-weighted fluid-attenuated inversion recovery (FLAIR) images were also obtained from all subjects. T1-weighted images were acquired using a magnetization prepared rapid gradient echo (MPRAGE) sequence with TR = 2530 ms, inversion time (TI) = 1100 ms, TE = 2.98 ms, FOV = 256 mm× 256 mm, matrix = 256 × 256, slice thickness = 1 mm, slice number = 192, acquisition time = 6.03 min. T2-weighted images were acquired using a fast spin echo sequence with TR = 5400 ms, TE = 105 ms, FOV = 220 mm× 220 mm, matrix = 448 × 448, slice thickness = 3 mm, slice number = 35, acquisition time = 1.44 min. T2-weighted FLAIR were obtained with TR = 5500 ms, TI = 1920 ms, TE = 83 ms, FOV = 220 mm× 220 mm, matrix = 320 × 320, slice thickness = 3 mm, slice number = 35, acquisition time = 2.47 min.
T1-weighted, T2-weighted, and FLAIR images were visually assessed by two neurologists (X-ZY and X-PW). Abnormalities in the corpus striatum, thalamus (Th), midbrain, pons, and cerebellum were recorded. Abnormal signal was defined as deviation from the conventionally accepted signal intensity of a given structure. Cortical and cerebellar atrophy was assessed according to the global cortical atrophy (GCA) score (Pasquier et al., 1996). Central atrophy was assessed by measuring the width of the third ventricle. WD patients were considered to have brain atrophy if the GCA score >1 and/or the width of the third ventricle >6 mm (Dusek et al., 2020). When disagreements arose, consensus was reached through discussion.
The susceptibility maps were reconstructed using the Morphology Enabled Dipole Inversion with automatic uniform cerebrospinal fluid (CSF) zero reference (MEDI + 0) algorithm (Liu et al., 2018). The reconstruction procedure involved the following steps: The field map was first estimated by performing a one-dimensional temporal unwrapping of the phase on each voxel followed by a weighted least squares fit of the temporally unwrapped phases in each voxel over TE (Liu et al., 2013). Then, a magnitude image guided spatial unwrapping algorithm was used to account for frequency aliasing on the field map (Cusack and Papadakis, 2002). Thirdly, the background field was removed using the projection onto dipole fields (PDF) method (Liu et al., 2011). Finally, the remaining tissue field was inverted to generate a susceptibility map using MEDI (Liu et al., 2012), and the susceptibility of ventricular CSF was used as zero reference (Straub et al., 2017). The ventricular CSF mask MCSF was determined from the brain ROI mask M and the R2∗ map by taking advantage of the low R2∗ values of CSF and by imposing connectivity (Supplementary Figure S1). The parameters for MEDI + 0 used the default settings in the MEDI toolbox, λ1 and λ2 were 1000 and 100, respectively. The threshold of R2∗ was 5 s−1 (Liu et al., 2018).
The 3D regions of interest (ROIs) including the bilateral head of the caudate nucleus (CN), globus pallidus (GP), Put, Th, substantia nigra (SN), red nucleus (RN), and dentate nucleus (DN), were manually segmented on the susceptibility maps by a rater who was blinded to subject information using ITK-SNAP1. ROIs were drawn according to their anatomical boundaries on all sections where the DGM was visible (Figure 1). The most inferior or most superior slice of ROIs were excluded to minimize partial volume effects, and voxels at the boundaries were also excluded. ROIs were confirmed by a senior neurologist (X-PW) with more than 30 years of neuroimaging experience. To account for the effect of head size variability across individuals, the volume of each ROI was normalized with the skull-scaling factor, which was obtained from T1-weighted images using FSL’s SIENAX tool (Smith et al., 2002). Susceptibility values of each ROI in the left and right hemispheres were averaged for further analysis.
Figure 1. Example images showing segmentation of regions of interest from a healthy control (28 years old, male). The QSM image (A) and example segmentations (D) include the head of the caudate nucleus (CN), globus pallidus (GP), putamen (Put) and thalamus (Th). The QSM image (B) and example segmentations (E) include the red nucleus (RN) and substantia nigra (SN). The QSM image (C) and example segmentation (F) include the dentate nucleus (DN).
Statistical Analyses
Mann-Whitney U test was used to compare the age distribution between WD patients and HC. Susceptibility values and volumes of ROIs in each group were compared using a linear regression model, and age and gender were used as covariates to correct for their potential influence on susceptibility. Correlation analysis was performed between the susceptibility values and volumes of ROIs and UWDRS neurological subscores. Since the data could not completely pass the normality test (Kolmogorov-Smirnov test), Spearman rank correlation coefficients (r) were calculated to quantify correlation strength. To control for the effects of age, gender, and disease duration, partial Spearman correlation coefficients (rp) were calculated. The P-values of multiple comparisons and correlation analysis were corrected by controlling the false discovery rate (FDR) at a level of 0.05 (Benjamini and Hochberg, 1995). To test whether susceptibility values of DGM could predict the severity of neurological impairment, a multiple regression model with UWDRS neurological subscores as the dependent variable and susceptibility values of DGM as predictor variables was performed. The multiple regression model was performed using a stepwise method, which continuously increased the number of predictor variables until the best predictive model was obtained. Collinearity diagnosis was also performed by checking tolerance and variance inflation factor (VIF). Statistical analyses were conducted using SPSS Statistics 22 (SPSS Inc., Chicago, IL, United States).
Results
Subjects Characteristics
Three WD patients with T1 hyperintensities in the lentiform nucleus were excluded. A total of 17 neurological WD patients (10 males and 7 females) and 20 age-matched HC (12 males and 8 females) were included for further analysis. The mean age of WD patients and HC was 32.23 ± 7.55 (range 20–48) years and 30.4 ± 7.59 (range 18–43) years, respectively. There were no statistically differences in age distribution between the two groups (P = 0.478). The mean disease duration (from symptom onset to the time of MRI scans) was 10.83 ± 6.74 (range 0.5–21) years, and the mean treatment duration was 8.15 ± 6.07 (range 0.5–20) years. Four WD patients were treated with D-penicillamine and 13 WD patients were treated with DMSA.
T2/FLAIR hyperintensities in lentiform nucleus, Th, midbrain, and tegmentum of pons were the most common findings and were observed in 15 of 17 (88%) WD patients (Supplementary Figure S2). Brain atrophy was observed in 11 of 17 (65%) WD patients. Demographic and clinical characteristics of all participants are summarized in Table 1.
Susceptibility and Volume Analysis
In visual observation, high signals (increased susceptibility) could be observed in the DGM of WD patients. GP, Put, and SN had the most significant susceptibility changes in WD patients compared to HC. In addition, DGM atrophy was a common imaging feature of WD patients (Figure 2). The mean susceptibility values of all examined DGM in WD patients were higher than those in HC. Significantly increased susceptibility values were found in CN, GP, Put, Th, SN, and RN in WD patients (P < 0.001). The susceptibility value of DN was slightly higher in WD patients than in HC (P = 0.031) (Figure 3). Volume reductions were found in CN, GP, Put, Th, and SN in WD patients (P < 0.001), but the volumes of RN (P = 0.114) and DN (P = 0.614) were not significantly changed (Figure 4). More information about the susceptibility and volumes of DGM can be found in Supplementary Table S1.
Figure 2. Comparison of QSM images from a WD patient (36 years old, male, A–C) and a healthy control (38 years old, male, D–F). The QSM images (A,D) include head of the caudate nucleus, globus pallidus, putamen and thalamus. The QSM images (B,E) include the red nucleus and substantia nigra. The QSM images (C,F) include the dentate nucleus. Significantly increased susceptibility was observed in the deep gray matter nuclei of the WD patient.
Figure 3. Comparison of mean susceptibility values of DGM between WD patients and HC. Susceptibility values of all examined DGM were higher in WD patients than in HC. *P < 0.05, and **P < 0.01. Error bars indicate standard deviation. CN, head of the caudate nucleus; DN, dentate nucleus; GP, globus pallidus; Put, putamen; RN, red nucleus; SN, substantia nigra; Th, thalamus.
Figure 4. Comparison of mean volumes of DGM between WD patients and HC. Volume reductions were observed in CN, GP, Put, Th, and SN in WD patients. **P < 0.01. Error bars indicate standard deviation. CN, head of the caudate nucleus; DN, dentate nucleus; GP, globus pallidus; Put, putamen; RN, red nucleus; SN, substantia nigra; Th, thalamus.
Correlation Between Susceptibility Values and Volumes of DGM
In WD patients, susceptibility values and volumes of DGM were negatively correlated in CN (r = −0.757, P = 0.003), GP (r = −0.534, P = 0.038), Put (r = −0.706, P = 0.005), Th (r = −0.613, P = 0.021), and RN (r = −0.539, P = 0.038). No such correlation was found in SN (r = 0.093, P = 0.722) and DN (r = 0.235, P = 0.424) (Figure 5). After controlling for age, gender, and disease duration, susceptibility values and volumes of DGM were still negatively correlated in CN (rp = −0.657, P = 0.037), Put (rp = −0.667, P = 0.037), and Th (rp = −0.613, P = 0.046). No correlation was found between susceptibility values and volumes of DGM in GP (rp = −0.502, P = 0.118), SN (rp = 0.063, P = 0.829), RN (rp = −0.457, P = 0.141), and DN (rp = 0.239, P = 0.478).
Figure 5. Correlation between susceptibility values and volumes of DGM in WD patients. Susceptibility values and volumes of DGM were negatively correlated in CN (A), Put (B), GP (C), Th (D), and RN (F). No such correlation was found in SN (E) and DN (G). CN, head of the caudate nucleus; DN, dentate nucleus; GP, globus pallidus; Put, putamen; RN, red nucleus; SN, substantia nigra; Th, thalamus.
In the group of HC, there was no correlation between susceptibility values and volumes of DGM in CN (r = −0.302, P = 0.308), GP (r = −0.287, P = 0.308), Put (r = −0.577, P = 0.056), Th (r = −0.005, P = 0.985), SN (r = −0.477, P = 0.119), RN (r = −0.006, P = 0.985), and DN (r = 0.414, P = 0.163). Similar results were obtained after controlling for age and gender: CN (rp = −0.313, P = 0.288), GP (rp = −0.275, P = 0.315), Put (rp = −0.313, P = 0.288), Th (rp = −0.018, P = 0.945), SN (rp = −0.381, P = 0.278), RN (rp = −0.464, P = 0.256), and DN (rp = 0.433, P = 0.256).
Correlation Between Susceptibility Values of DGM and Neurological Impairment
Unified Wilson’s Disease Rating Scale neurological subscale was used to assess the severity of neurological impairment in WD patients, with a higher score indicating a more severe condition. The susceptibility values of DGM were positively correlated with UWDRS neurological subscores in CN (r = 0.701, P = 0.006), GP (r = 0.632, P = 0.011), Put (r = 0.727, P = 0.006), Th (r = 0.656, P = 0.01), and RN (r = 0.608, P = 0.014) in WD patients. No such correlation was found in SN (r = 0.432, P = 0.097) and DN (r = 0.374, P = 0.139) (Figure 6). After controlling for age, gender, and disease duration, susceptibility values of DGM were positively correlated with UWDRS neurological subscores in all examined DGM: CN (rp = 0.669, P = 0.016), GP (rp = 0.616, P = 0.026), Put (rp = 0.76, P = 0.006), Th (rp = 0.72, P = 0.009), SN (rp = 0.588, P = 0.03), RN (rp = 0.766, P = 0.006), and DN (rp = 0.58, P = 0.03).
Figure 6. Correlation between susceptibility values of DGM and UWDRS neurological subscores. Susceptibility values were positively correlated with UWDRS neurological subscores in CN (A), GP (B), Put (C), Th (D), and RN (F). No such correlation was found in SN (E) and DN (G). CN, head of the caudate nucleus; DN, dentate nucleus; GP, globus pallidus; Put, putamen; RN, red nucleus; SN, substantia nigra; Th, thalamus.
The multiple regression model revealed that the susceptibility values of Put, CN, and DN could well predict UWDRS neurological subscores. The susceptibility value of Put was superior to that of other DGM nuclei in predicting the UWDRS neurological subscores. In the model with only the susceptibility value of Put as a predictor variable, UWDRS neurological subscores increased by 1 point for every 6 ppb increase in Put (P = 0.005; R2 = 0.415). The model using the susceptibility values of CN and DN as predictor variables had the best ability in predicting the UWDRS neurological subscores. UWDRS neurological subscores increased by 1 point for every 6 ppb increase in CN (P < 0.001) or every 8 ppb increase in DN (P = 0.001; R2 = 0.745). Details about the multiple regression model can be found in Supplementary Table S2.
Correlation Between Volumes of DGM and Neurological Impairment
There was no correlation between volumes of DGM and UWDRS neurological subscores in CN (r = −0.432, P = 0.146), GP (r = −0.625, P = 0.051), Put (r = −0.311, P = 0.262), Th (r = −0.513, P = 0.082), SN (r = −0.549, P = 0.079), RN (r = −0.367, P = 0.206), and DN (r = 0.205, P = 0.43). After controlling for age, gender, and disease duration, a negative correlation was found between the volume of Th and UWDRS neurological subscores (rp = −0.689, P = 0.045). No such correlation was found in CN (rp = −0.331, P = 0.289), GP (rp = −0.575, P = 0.11), Put (rp = −0.335, P = 0.289), SN (rp = −0.502, P = 0.118), RN (rp = −0.503, P = 0.118), and DN (rp = 0.234, P = 0.42).
Discussion
In this cross-sectional study, we found that the group average susceptibility values of all examined DGM in neurological WD patients were significantly higher than those in HC. Volume reductions were observed in CN, GP, Put, Th, and SN of WD patients, but RN in the midbrain and DN in the cerebellum were not affected. Susceptibility values and volumes of DGM were negatively correlated in CN, Put, and Th. In addition, increased DGM susceptibility was associated with high UWDRS neurological subscores in all examined DGM. There was only a weak correlation between volumes of DGM and UWDRS neurological subscores in Th. The results of our study were consistent with the theory that iron and copper accumulation in DGM was a character of neuropathology in neurological WD patients. To our knowledge, our study provided the first in vivo evidence that metal accumulation in the DGM was associated with DGM atrophy and neurological impairment in neurological WD patients.
Studies using traditional GRE MRI sequences such as T2∗-weighted images and SWI had identified iron accumulation in the brains of WD patients (Yang et al., 2015; Dusek et al., 2017). However, traditional GRE MRI sequences always suffer from blooming artifacts and the non-local effect (Li et al., 2012). Most importantly, traditional GRE MRI sequences cannot quantify susceptibility source without deconvolution, and phase changes cannot reflect the actual tissue susceptibility (Wang et al., 2017). QSM is a recently developed GRE MRI post-processing technique, which can quantitatively measure tissue susceptibility and is less affected by the external magnetic field (de Rochefort et al., 2010; Wang et al., 2017). Given iron content is about 3-fold higher than copper content in the DGM of neurological WD patients and the strong paramagnetism of iron species (Dusek et al., 2017; Dezortova et al., 2019), we attribute the increased susceptibility of DGM primarily to iron accumulation. However, we cannot exclude the influence of copper accumulation on susceptibility changes. In our study, significantly increased susceptibility was detected in all examined DGM nuclei, indicating metal accumulation in the DGM is common neuropathological feature of neurological WD patients. The highest susceptibilities values were observed in GP and SN, which was consistent with the distribution of iron species in neurological WD patients (Dusek et al., 2017). Whether iron species accumulate in the dentate nucleus in WD patients is controversial. Several postmortem studies analyzed the iron content in DN in WD patients, but the results were quite different (Litwin et al., 2013b; Dusek et al., 2017). We found that the mean susceptibility value of DN in WD patients was slightly higher than that in HC. The results of our study suggested that iron species tended to accumulate in the basal ganglia and brain stem of WD patients, while DN in the cerebellum was less affected.
Interestingly, increased susceptibility was also observed in Th, which contains large amount of white matter fiber tracts. Myelin in white matter is diamagnetic, and the increased susceptibility in Th may be caused by demyelination besides metal accumulation (Hametner et al., 2018; Dezortova et al., 2019). Demyelination is a common pathological feature of white matter in WD patients. Recent studies also revealed microstructure changes in Th and white matter of WD patients using diffusion tensor imaging (DTI) (Li et al., 2014; Dong et al., 2019). However, we cannot determine whether metal accumulation or demyelination is the main factor leading to increased susceptibility in Th. In addition, an increase in the mean bulk susceptibility of DGM may be induced by DGM atrophy. Metal concentration increases with tissue atrophy, but it is unclear whether the total metal content increases (Hernandez-Torres et al., 2019). Increased susceptibility was also observed in the DGM without atrophy (RN and DN) in WD patients. Therefore, we speculate that paramagnetic metal accumulation occurs in the DGM of WD patients. Histopathological studies are needed to confirm whether there is paramagnetic metal accumulation in the DGM of WD patients.
Brain atrophy is a common MRI feature of neurological WD patients and is associated with the severity of the neurological impairment (Smolinski et al., 2019; Dusek et al., 2020). Brain atrophy seems to be an irreversible process, but the neurological symptoms of WD patients are partly reversible after anti-copper treatment. Therefore, brain volume may be more appropriate as an assessment for chronic treated WD patients than as a marker for treatment monitoring. In WD patients, basal ganglia and brainstem are the most vulnerable regions, especially Put, which is a softening, brownish, atrophic nuclei with cavitations at macroscopic examination (Poujois et al., 2017). In this study, we found significant DGM atrophy in CN, GP, Put, Th, and SN, but RN and DN were not affected. Furthermore, we found that susceptibility values and volumes of DGM were negatively correlated in CN, Put, and Th, suggesting that metal accumulation in the DGM was associated with neuronal loss. The results of our study were consistent with the study of Dusek et al. (2017) which found Put atrophy was accompanied with increased iron content and iron-containing macrophages. However, it is not clear whether iron accumulation was merely a biomarker of neurodegeneration or a cause or contributor to the neurodegeneration itself (Yan et al., 2018). The advantage of QSM in measuring DGM volumes is that it can show the boundaries of DGM more clearly than traditional MRI sequences, which can ensure the accuracy of segmentation (Liu et al., 2013). To our knowledge, this study is the first one to investigate the DGM volume changes in WD using QSM.
Brain atrophy and iron accumulation is a common feature of neurodegenerative diseases such as Alzheimer’s disease and Parkinson’s disease (Lane et al., 2018; Jiang et al., 2019). Increased iron accumulation in the DGM may be explained by the influx of iron-containing phagocytic cells (macrophages), mitochondria damage and energy production failure, or decreased iron efflux due to ceruloplasmin dysfunction (Dusek et al., 2017, 2018). Ceruloplasmin is a ferroxidase that plays a critical role in maintaining iron homeostasis and preventing the formation of free radicals (Wang and Wang, 2019). Decreased serum ceruloplasmin may aggravate iron accumulation in WD, which is similar to the pathogenesis of aceruloplasminemia (Piperno and Alessio, 2018). As mentioned above, iron accumulation may be caused by neurodegeneration in WD, and iron accumulation also exacerbates neurodegeneration due to its redox properties.
We assessed the severity of neurological impairment of WD patients using UWDRS neurological subscale, and found that the susceptibility values of DGM were positively correlated with the severity of neurological impairment in all examined DGM. The correlation between iron accumulation and neuropathological changes of Put has been confirmed by a histopathological study (Dusek et al., 2017). Brain atrophy and brain MRI abnormalities reflecting iron accumulation, such as T2∗ and SWI hypointensities, were proposed as a part of imaging scale to assess the severity of MRI changes in WD. These chronic brain abnormalities were considered irreversible overtime, and showed positive correlation with the severity of neurological impairment (Dusek et al., 2020). Our study provided in vivo evidence that metal accumulation in the DGM was associated with DGM atrophy and neurological impairment in neurological WD patients. In addition, we found the susceptibility values of Put, CN, and DN could well predict the severity of neurological impairment. The susceptibility value of Put was superior to that of other DGM nuclei in predicting the UWDRS neurological subscores. The model using the susceptibility values of CN and DN as predictor variables had the best ability in predicting the UWDRS neurological subscores.
Few studies focused on the correlations between DGM atrophy and severity of neurological impairment in WD. Zou et al. (2019) reported that the volumes of bilateral Put and left GP were associated with the severity of neurological impairment in WD, which was assessed using modified Young scale. However, in this study, we only found a weak correlation between the volumes of DGM and the UWDRS neurological subscores in Th. The method of DGM volume assessment (manual or automatic segmentation), the selection of neurological scales, the small number of WD patients, and the high heterogeneity of WD patients may be the reasons for different results of the two studies. More studies are needed to determine the relationship between DGM atrophy and the severity of neurological impairment in WD. Compared with the volumes of DGM, the susceptibility values of DGM showed better correlations with UWDRS neurological subscores. The results of our study suggested that the susceptibility values of DGM could be used as a biomarker to assess the severity of neurodegeneration, which would be helpful in the diagnosis of WD, treatment monitoring, and as a surrogate outcome for clinical trials.
Several limitations in this study have to be mentioned. In some WD patients, DGM atrophy and local high susceptibility make the boundaries of DGM unclear, which may affect the accuracy of DGM segmentation. WD patients included in this study were treated patients with mild or moderate neurological symptoms. Patients with severe neurological impairment were excluded because they could not cooperate or remain still during MRI scans. It is unclear whether the susceptibility of DGM is still associated with neurological impairment in patients with severe neurological symptoms or in untreated patients. Due to limited sample size, the results of our study need to be confirmed in a larger cohort. In addition, a longitudinal study is required to determine whether the susceptibility of DGM changes with disease progression or symptom improvement.
In conclusion, our study found paramagnetic metal accumulation in the DGM of neurological WD patients, and provided in vivo evidence that metal accumulation in the DGM was associated with DGM atrophy and neurological impairment. The susceptibility of DGM could be used as a biomarker to assess the severity of neurodegeneration in WD.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Medical Ethics Committee of Tongren Hospital, Shanghai Jiao Tong University School of Medicine. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
X-ZY and J-LC developed the conception and design of the study, collected clinical and imaging data, and analyzed and interpreted the data. G-YL analyzed and interpreted the data, and revised the manuscript. J-QL and X-PW developed the conception and design of the study, and revised the manuscript. X-ZY wrote the first draft, and all of the authors critically reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by grants from the National Natural Science Foundation of China (81671103) and Shanghai Jiao Tong University Medical School Innovative Research Plan 2015 (TM201520).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank all participants involved in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2020.573633/full#supplementary-material
Footnotes
References
Bandmann, O., Weiss, K. H., and Kaler, S. G. (2015). Wilson’s disease and other neurological copper disorders. Lancet Neurol. 14, 103–113. doi: 10.1016/S1474-4422(14)70190-5
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300.
Bilgic, B., Pfefferbaum, A., Rohlfing, T., Sullivan, E. V., and Adalsteinsson, E. (2012). MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. Neuroimage 59, 2625–2635. doi: 10.1016/j.neuroimage.2011.08.077
Cusack, R., and Papadakis, N. (2002). New robust 3-D phase unwrapping algorithms: application to magnetic field mapping and undistorting echoplanar images. Neuroimage 16(3 Pt 1), 754–764. doi: 10.1006/nimg.2002.1092
Czlonkowska, A., and Litwin, T. (2017). Wilson disease - currently used anticopper therapy. Handb. Clin. Neurol. 142, 181–191. doi: 10.1016/B978-0-444-63625-6.00015-X
Czlonkowska, A., Litwin, T., Dusek, P., Ferenci, P., Lutsenko, S., Medici, V., et al. (2018). Wilson disease. Nat. Rev. Dis. Primers 4:21. doi: 10.1038/s41572-018-0018-3
de Rochefort, L., Liu, T., Kressler, B., Liu, J., Spincemaille, P., Lebon, V., et al. (2010). Quantitative susceptibility map reconstruction from MR phase data using bayesian regularization: validation and application to brain imaging. Magn. Reson. Med. 63, 194–206. doi: 10.1002/mrm.22187
Dezortova, M., Lescinskij, A., Dusek, P., Herynek, V., Acosta-Cabronero, J., Bruha, R., et al. (2019). Multiparametric quantitative brain MRI in neurological and hepatic forms of Wilson’s disease. J. Magn. Reson. Imaging 51, 1829–1835. doi: 10.1002/jmri.26984
Doganay, S., Gumus, K., Koc, G., Bayram, A. K., Dogan, M. S., Arslan, D., et al. (2018). Magnetic susceptibility changes in the basal ganglia and brain stem of patients with wilson’s disease: evaluation with quantitative susceptibility mapping. Magn. Reson. Med. Sci. 17, 73–79. doi: 10.2463/mrms.mp.2016-0145
Dong, T., Yang, W. M., Wu, M. C., Zhang, J., Huang, P., Xu, C. S., et al. (2019). Microstructure changes in whiter matter relate to cognitive impairment in Wilson’s disease. Biosci. Rep. 39:BSR20181651. doi: 10.1042/BSR20181651
Dusek, P., Bahn, E., Litwin, T., Jablonka-Salach, K., Luciuk, A., Huelnhagen, T., et al. (2017). Brain iron accumulation in Wilson disease: a post mortem 7 Tesla MRI - histopathological study. Neuropathol. Appl. Neurobiol. 43, 514–532. doi: 10.1111/nan.12341
Dusek, P., Roos, P. M., Litwin, T., Schneider, S. A., Flaten, T. P., and Aaseth, J. (2015). The neurotoxicity of iron, copper and manganese in Parkinson’s and Wilson’s diseases. J. Trace Elem. Med. Biol. 31, 193–203. doi: 10.1016/j.jtemb.2014.05.007
Dusek, P., Skoloudik, D., Maskova, J., Huelnhagen, T., Bruha, R., Zahorakova, D., et al. (2018). Brain iron accumulation in Wilson’s disease: a longitudinal imaging case study during anticopper treatment using 7.0T MRI and transcranial sonography. J. Magn. Reson. Imaging 47, 282–285. doi: 10.1002/jmri.25702
Dusek, P., Smolinski, L., Redzia-Ogrodnik, B., Golebiowski, M., Skowronska, M., Poujois, A., et al. (2020). Semiquantitative scale for assessing brain MRI abnormalities in wilson disease: a validation study. Mov. Disord. 35, 994–1001. doi: 10.1002/mds.28018
Ferenci, P., Caca, K., Loudianos, G., Mieli-Vergani, G., Tanner, S., Sternlieb, I., et al. (2003). Diagnosis and phenotypic classification of Wilson disease. Liver Int. 23, 139–142. doi: 10.1034/j.1600-0676.2003.00824.x
Fritzsch, D., Reiss-Zimmermann, M., Trampel, R., Turner, R., Hoffmann, K. T., and Schafer, A. (2014). Seven-tesla magnetic resonance imaging in Wilson disease using quantitative susceptibility mapping for measurement of copper accumulation. Invest. Radiol. 49, 299–306. doi: 10.1097/RLI.0000000000000010
Hametner, S., Endmayr, V., Deistung, A., Palmrich, P., Prihoda, M., Haimburger, E., et al. (2018). The influence of brain iron and myelin on magnetic susceptibility and effective transverse relaxation - A biochemical and histological validation study. Neuroimage 179, 117–133. doi: 10.1016/j.neuroimage.2018.06.007
Hernandez-Torres, E., Wiggermann, V., Machan, L., Sadovnick, A. D., Li, D. K. B., Traboulsee, A., et al. (2019). Increased mean R2∗ in the deep gray matter of multiple sclerosis patients: have we been measuring atrophy? J. Magn. Reson. Imaging 50, 201–208. doi: 10.1002/jmri.26561
Jiang, H., Song, N., Jiao, Q., Shi, L., and Du, X. (2019). Iron pathophysiology in Parkinson diseases. Adv. Exp. Med. Biol. 1173, 45–66. doi: 10.1007/978-981-13-9589-5_4
King, A. D., Walshe, J. M., Kendall, B. E., Chinn, R. J., Paley, M. N., Wilkinson, I. D., et al. (1996). Cranial MR imaging in Wilson’s disease. AJR Am. J. Roentgenol. 167, 1579–1584. doi: 10.2214/ajr.167.6.8956601
Kozic, D., Svetel, M., Petrovic, B., Dragasevic, N., Semnic, R., and Kostic, V. S. (2003). MR imaging of the brain in patients with hepatic form of Wilson’s disease. Eur. J. Neurol. 10, 587–592. doi: 10.1046/j.1468-1331.2003.00661.x
Lane, D. J. R., Ayton, S., and Bush, A. I. (2018). Iron and Alzheimer’s disease: an update on emerging mechanisms. J. Alzheimers Dis. 64, S379–S395. doi: 10.3233/JAD-179944
Langkammer, C., Schweser, F., Krebs, N., Deistung, A., Goessler, W., Scheurer, E., et al. (2012). Quantitative susceptibility mapping (QSM) as a means to measure brain iron? A post mortem validation study. Neuroimage 62, 1593–1599. doi: 10.1016/j.neuroimage.2012.05.049
Lee, S., Nam, Y., Jang, J., Na, G. H., Kim, D. G., Shin, N. Y., et al. (2018). Deep gray matter iron measurement in patients with liver cirrhosis using quantitative susceptibility mapping: relationship with pallidal T1 hyperintensity. J. Magn. Reson. Imaging 47, 1342–1349. doi: 10.1002/jmri.25841
Leinweber, B., Moller, J. C., Scherag, A., Reuner, U., Gunther, P., Lang, C. J., et al. (2008). Evaluation of the unified wilson’s disease rating scale (UWDRS) in German patients with treated Wilson’s disease. Mov. Disord. 23, 54–62. doi: 10.1002/mds.21761
Li, G., Wu, R., Tong, R., Bo, B., Zhao, Y., Gillen, K. M., et al. (2020). Quantitative measurement of metal accumulation in brain of patients with wilson’s disease. Mov. Disord. doi: 10.1002/mds.28141 [Epub ahead of print].
Li, G., Zhou, X., Xu, P., Pan, X., and Chen, Y. (2014). Microstructure assessment of the thalamus in Wilson’s disease using diffusion tensor imaging. Clin. Radiol. 69, 294–298. doi: 10.1016/j.crad.2013.10.016
Li, J., Chang, S., Liu, T., Wang, Q., Cui, D., Chen, X., et al. (2012). Reducing the object orientation dependence of susceptibility effects in gradient echo MRI through quantitative susceptibility mapping. Magn. Reson. Med. 68, 1563–1569. doi: 10.1002/mrm.24135
Li, J., Lin, H., Liu, T., Zhang, Z., Prince, M. R., Gillen, K., et al. (2018). Quantitative susceptibility mapping (QSM) minimizes interference from cellular pathology in R2∗ estimation of liver iron concentration. J. Magn. Reson. Imaging 48, 1069–1079. doi: 10.1002/jmri.26019
Li, W. J., Chen, C., You, Z. F., Yang, R. M., and Wang, X. P. (2016). Current drug managements of Wilson’s disease: from west to east. Curr. Neuropharmacol. 14, 322–325. doi: 10.2174/1570159x14666151130222427
Litwin, T., Gromadzka, G., Czlonkowska, A., Golebiowski, M., and Poniatowska, R. (2013a). The effect of gender on brain MRI pathology in Wilson’s disease. Metab. Brain Dis. 28, 69–75. doi: 10.1007/s11011-013-9378-2
Litwin, T., Gromadzka, G., Szpak, G. M., Jablonka-Salach, K., Bulska, E., and Czlonkowska, A. (2013b). Brain metal accumulation in Wilson’s disease. J. Neurol. Sci. 329, 55–58. doi: 10.1016/j.jns.2013.03.021
Liu, T., Eskreis-Winkler, S., Schweitzer, A. D., Chen, W., Kaplitt, M. G., Tsiouris, A. J., et al. (2013). Improved subthalamic nucleus depiction with quantitative susceptibility mapping. Radiology 269, 216–223. doi: 10.1148/radiol.13121991
Liu, T., Khalidov, I., de Rochefort, L., Spincemaille, P., Liu, J., Tsiouris, A. J., et al. (2011). A novel background field removal method for MRI using projection onto dipole fields (PDF). NMR Biomed. 24, 1129–1136. doi: 10.1002/nbm.1670
Liu, T., Xu, W., Spincemaille, P., Avestimehr, A. S., and Wang, Y. (2012). Accuracy of the morphology enabled dipole inversion (MEDI) algorithm for quantitative susceptibility mapping in MRI. IEEE Trans. Med. Imaging 31, 816–824. doi: 10.1109/TMI.2011.2182523
Liu, Z., Spincemaille, P., Yao, Y., Zhang, Y., and Wang, Y. (2018). MEDI+0: morphology enabled dipole inversion with automatic uniform cerebrospinal fluid zero reference for quantitative susceptibility mapping. Magn. Reson. Med. 79, 2795–2803. doi: 10.1002/mrm.26946
Pasquier, F., Leys, D., Weerts, J. G., Mounier-Vehier, F., Barkhof, F., and Scheltens, P. (1996). Inter- and intraobserver reproducibility of cerebral atrophy assessment on MRI scans with hemispheric infarcts. Eur. Neurol. 36, 268–272. doi: 10.1159/000117270
Piperno, A., and Alessio, M. (2018). Aceruloplasminemia: waiting for an efficient therapy. Front. Neurosci. 12:903. doi: 10.3389/fnins.2018.00903
Poujois, A., Mikol, J., and Woimant, F. (2017). Wilson disease: brain pathology. Handb. Clin. Neurol. 142, 77–89. doi: 10.1016/B978-0-444-63625-6.00008-2
Skowronska, M., Litwin, T., Dziezyc, K., Wierzchowska, A., and Czlonkowska, A. (2013). Does brain degeneration in Wilson disease involve not only copper but also iron accumulation? Neurol. Neurochir. Pol. 47, 542–546. doi: 10.5114/ninp.2013.39071
Smith, S. M., Zhang, Y., Jenkinson, M., Chen, J., Matthews, P. M., Federico, A., et al. (2002). Accurate, robust, and automated longitudinal and cross-sectional brain change analysis. Neuroimage 17, 479–489. doi: 10.1006/nimg.2002.1040
Smolinski, L., Litwin, T., Redzia-Ogrodnik, B., Dziezyc, K., Kurkowska-Jastrzebska, I., and Czlonkowska, A. (2019). Brain volume is related to neurological impairment and to copper overload in Wilson’s disease. Neurol. Sci. 40, 2089–2095. doi: 10.1007/s10072-019-03942-z
Straub, S., Schneider, T. M., Emmerich, J., Freitag, M. T., Ziener, C. H., Schlemmer, H. P., et al. (2017). Suitable reference tissues for quantitative susceptibility mapping of the brain. Magn. Reson. Med. 78, 204–214. doi: 10.1002/mrm.26369
Valdes Hernandez Mdel, C., Maconick, L. C., Tan, E. M., and Wardlaw, J. M. (2012). Identification of mineral deposits in the brain on radiological images: a systematic review. Eur. Radiol. 22, 2371–2381. doi: 10.1007/s00330-012-2494-2
Volpert, H. M., Pfeiffenberger, J., Groner, J. B., Stremmel, W., Gotthardt, D. N., Schafer, M., et al. (2017). Comparative assessment of clinical rating scales in Wilson’s disease. BMC Neurol. 17:140. doi: 10.1186/s12883-017-0921-3
Wang, B., and Wang, X. P. (2019). Does ceruloplasmin defend against neurodegenerative diseases? Curr. Neuropharmacol. 17, 539–549. doi: 10.2174/1570159X16666180508113025
Wang, Y., Spincemaille, P., Liu, Z., Dimov, A., Deh, K., Li, J., et al. (2017). Clinical quantitative susceptibility mapping (QSM): biometal imaging and its emerging roles in patient care. J. Magn. Reson. Imaging 46, 951–971. doi: 10.1002/jmri.25693
Yan, F., He, N., Lin, H., and Li, R. (2018). Iron deposition quantification: applications in the brain and liver. J. Magn. Reson. Imaging 48, 301–317. doi: 10.1002/jmri.26161
Yang, J., Li, X., Yang, R., Yu, X., Yu, C., Qian, Y., et al. (2015). Susceptibility-weighted imaging manifestations in the brain of Wilson’s disease patients. PLoS One 10:e0125100. doi: 10.1371/journal.pone.0125100
Zaino, D., Chiarotti, I., Battisti, C., Salvatore, S., Federico, A., and Cerase, A. (2019). Six-year clinical and MRI quantitative susceptibility mapping (QSM) follow-up in neurological Wilson’s disease under zinc therapy: a case report. Neurol. Sci. 40, 199–201. doi: 10.1007/s10072-018-3557-1
Keywords: Wilson’s disease, quantitative susceptibility mapping, iron, copper, neurodegeneration
Citation: Yuan X-Z, Li G-Y, Chen J-L, Li J-Q and Wang X-P (2020) Paramagnetic Metal Accumulation in the Deep Gray Matter Nuclei Is Associated With Neurodegeneration in Wilson’s Disease. Front. Neurosci. 14:573633. doi: 10.3389/fnins.2020.573633
Received: 17 June 2020; Accepted: 27 August 2020;
Published: 16 September 2020.
Edited by:
Mark Haacke, Wayne State University, United StatesReviewed by:
Tomasz Litwin, Institute of Psychiatry and Neurology, PolandFerdinand Schweser, University at Buffalo, United States
Petr Dušek, Charles University, Czechia
Xu Li, Kennedy Krieger Institute, United States
Copyright © 2020 Yuan, Li, Chen, Li and Wang. 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: Jian-Qi Li, anFsaUBwaHkuZWNudS5lZHUuY24=; Xiao-Ping Wang, d2FuZ3hwQHVzdGMuZWR1