Quantitative Susceptibility Mapping of Brain Iron Deposition in Patients With Recurrent Depression

Xinxiu Duan; Yuhang Xie; Xiufang Zhu; Lei Chen; Feng Li; Guoquan Feng; Lei Li

doi:10.30773/pi.2022.0110

Psychiatry Investig > Volume 19(8); 2022 > Article

Duan, Xie, Zhu, Chen, Li, Feng, and Li: Quantitative Susceptibility Mapping of Brain Iron Deposition in Patients With Recurrent Depression

Original Article

Psychiatry Investigation 2022;19(8):668-675.

Published online: August 24, 2022

DOI: https://doi.org/10.30773/pi.2022.0110

Quantitative Susceptibility Mapping of Brain Iron Deposition in Patients With Recurrent Depression

Xinxiu Duan^1,^*

, Yuhang Xie^2,^*

, Xiufang Zhu¹

, Lei Chen¹

, Feng Li²

, Guoquan Feng²

, Lei Li¹

¹Department of Radiology, The First People’s Hospital of Lianyungang, Lianyungang, China

²Department of Radiology, Affiliated Hospital of Jiangsu University, Zhenjiang, China

Correspondence: Lei Li, MD Department of Radiology, The First People’s Hospital of Lianyungang, No.6, Huadong Road, Lianyungang 222061, China Tel: +86-518-85605135, Fax: +86-518-85605135, E-mail: lileijs6969@126.com

^* These authors contributed equally to this work.

Received April 17, 2022 Revised May 23, 2022 Accepted June 8, 2022

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Objective

Recurrence is the most significant feature of depression and the relationship between iron and recurrent depression is still lack of direct evidence in vivo.

Methods

Twenty-one patients with depression and twenty control subjects were included. Gradient-recalled echo, T1 and T2 images were acquired using a 3.0T MRI system. After quantitative susceptibility mapping were reconstructed and standardized, a whole-brain and the regions of interest were respectively analyzed.

Results

Significant increases in susceptibility were found in multiple recurrent depression patients, which involved several brain regions (frontal lobes, temporal lobe structures, occipital lobes hippocampal regions, putamen, thalamus, cingulum, and cerebellum). Interestingly, no susceptibility changes after treatment compared to pre-treatment (all p>0.05) and no significant correlation between susceptibility and Hamilton Depression Rating Scale were found. Besides, it was close to significance that those with a higher relapse frequency or a longer mean duration of single episode had a higher susceptibility in the putamen, thalamus, and hippocampus. Further studies showed susceptibility across the putamen (ρ²=0.27, p<0.001), thalamus (ρ²=0.21, p<0.001), and hippocampus (ρ²=0.19, p<0.001) were strongly correlated with total course of disease onset.

Conclusion

Brain iron deposition is related to the total course of disease onset, but not the severity of depression, which suggest that brain iron deposition may be a sign of brain damage in multiple recurrent depression.

Keywords: Iron, Biomarkers, Depression, Brain, Magnetic resonance imaging

INTRODUCTION

Depression is a heterogeneous and etiologically complex psychiatric syndrome [1-3], and its pathophysiology is still unclear.

Iron, an essential microelement in the body, plays important roles in physiological functions and the development of human brain, including oxygen transport, DNA synthesis and repair, electron transport, and neurotransmitter metabolism [4,5]. Decrease of iron not only result in the accumulation of monoamine oxidase that significantly reducing catechol transmitters, but also preventing the degradation of serotonin and decreasing the density and activity of dopamine [6]. On the other hand, excessive iron deposition can prevent dopamine synthesis by inducing toxic free radicals, which is involved in emotional and cognitive processing [7,8]. Although these reports have suggested the abnormal iron might influence the metabolism of depression related neurotransmitters from the side, the relationship between iron and depression is still lack of direct evidence in vivo.

Magnetic resonance imaging (MRI) has been regarded as one of the most powerful methods of neuroscience in vivo for non-invasive detection of brain tissue structure, function and metabolism [9,10]. Using the sequences of T2 or T2*, brain iron have been detected in several neuro-diseases; however, these results were limited by significant technical deficiencies— non-quantitative measurement, in which phase signal was easily affected by the formation of the magnetic field distribution because of the nonlocal nature of the magnetic field distribution [11]. Quantitative susceptibility mapping (QSM) is an independent of field strength and object shape imaging that has been considered to be the optimum quantitative detection for brain iron [12,13]. The successful application of this imaging technology in diseases (Huntington’s disease, Parkinson’s disease [PD], etc) makes it possible to study the brain iron of patients with depression [14,15].

Here, using QSM measure and voxel-based analyses, the objective of this study was to systematically investigate the changes of brain iron deposition in the whole brain and regions of interest (ROIs) of depressed patients, and to assess the relationships between susceptibility and depression-related conditions.

METHODS

Subjects

Basing on the results from the Chinese version of the Modified Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders (4th ed), twenty-one recurrent depression patients were recruited from the department of neurology and psychiatry. Inclusion criteria for enrollment included no previous history of neurological diseases, tumors, cerebral microbleeds, major cardiovascular risk factors, chronic diseases, or head injuries. All subjects were scanned by 3T MR at baseline before treatment and follow-up after treatment (depressive symptoms disappeared, Hamilton Depression Rating Scale [HDRS]≤8), and received extended diagnostic examinations, including a review on clinical history, and a 24-item HDRS assessment. Additionally, according to the inclusion criteria, twenty healthy control subjects (matched for age, sex, education, and handedness) were also recruited. Group demographic details are summarized in Table 1. All the examinations were performed with the understanding and written consent of each subject, with approval from the ethics committee of The First People’s Hospital of Lianyungang (LYG-20180438), and in compliance with national legislation and Declaration of Helsinki guidelines.

Image acquisition

All images were acquired using a Siemens Magnetom Trio 3.0T scanner with a 16-channel head coil (Siemens Corp, Erlangen, Germany). T1-weighted 3D-MPRAGE anatomical images were also collected to independently resolve the underlying brain anatomy. The imaging parameters were as follows: repeat time (TR)=14 ms with echo time (TE)=4.92 ms, slice thickness=1 mm, flip angle=7°, bandwidth/pixel=140 Hz/pixel, voxel resolution=1×1×1 mm³, field of view (FOV)= 25 cm, matrix size=256×256×192, scan time=3 min 20 seconds, and parallel imaging factors=2. T2-weighted images were also acquired with turbo spin-echo for visual inspection to exclude brain abnormalities. The imaging parameters were as follows: TR=5,400 ms with TE=110 ms, slice thickness=2 mm, flip angle=150°, bandwidth/pixel= 40 Hz/pixel, turbo factor= 18, voxel resolution=1×1×2 mm³, FOV=25 cm, matrix size= 256×256, scan time=2 min 40 seconds, and parallel imaging factors=2. Gradient-recalled echo (GRE) data were collected as complex MRI signals from each receiver channel. The imaging parameters were as follows: TR=28 ms with TE=16 ms, slice thickness=2 mm, flip angle=7°, bandwidth/pixel=120 Hz/pixel, voxel resolution=1×1×2 mm³, FOV=25 cm, matrix size=256×256×80, and scan time=5 min and 10 seconds. To minimize motion and increase inter-subject reproducibility in positioning, a thin pillow was placed on the base of the coil surrounding the sides and back of the head.

Data preprocessing

The magnitude images were projected with the BET algorithm (threshold setting 0.2) in the FMRIB Software Library (FSL) v5.0.9 (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/) to generate brain masks [16]. Phase images were unwrapped with the Laplacian approach [17]. Applying the method of sophisticated harmonic artifact reduction on phase data filtering, the background field was removed using a kernel size with a maximum radius of 6 mm and a singular value decomposition threshold of 0.05 [11]. Finally, an iterative algorithm was used to generate QSM images [18]. As previous study suggested that reference normalization is only of a small adjustment [19], susceptibility was not adjusted with a reference region in this study. The original susceptibility maps were derived from the following equation (for a left-handed system) [20]:

△χ=FT−1{FT(−Pcor/[γB0TE]/(13−kZ2kX2+kY2+kZ2)}

wherein, P_cor is the phase distribution of the unwrapped and background-field-corrected phase map; γ is the gyromagnetic ratio for hydrogen protons; B0 is the main magnetic field strength; TE is the echo time; and kx, ky, kz, are coordinates in k-space.

Spatial standardization

To align pixels to the same anatomical position in the brains of different subjects, the following post-processing steps were performed using a Statistical Parametric Mapping Version 8 (SPM8) program (http://www.fil.ion.ucl.ac.uk/spm), which was run on matrix laboratory (MATLAB) 2013a (https://www.mathworks.com/products/matlab.html): The individual high-resolution MPRAGE-T1WIs were initially normalized to Montreal Neurological Institute (MNI), and co-registered to the QSM image after motion correction using a linear transformation. Then, the co-registered QSM images were normalized to MNI space using the normalized parameters of the MPRAGE-T1WIs. Applying the smoothing-compensation strategy proposed by a previously study [21], 3D Gaussian kernel (3 mm) was performed for correcting co-registration errors and other imperfections.

Regions of interest extraction

ROIs were extracted bilaterally in FSL according to the BN_Atlas_274_noCb_uint16 template (https://scalablebrainatlas.incf.org/human/BNA) [22]. In order to reduce the marginal effects caused by standardized defects and registration errors, ROIs were eroded in 3D by convolution with a 1-mmradius spherical kernel. Mask overlays were checked to ensure that all masks excluded spurious voxels. The iron concentration in postmortem samples as reported by Hallgren and Sourander23 was used to validate whether QSM data provided a reliable quantitative measure of iron.

Statistical methods

Whole-brain (depression versus control) permutation analysis was conducted with the Statistical nonParametric Mapping toolbox (SnPM) in SPM8. Significant clusters were determined with 20,000 data permutations. A significance level of p=0.005 was applied with correction for multiple comparisons using the family-wise error (FWE) (p<0.05) method and clusters with at least 10 contiguous voxels. As a large number of previous studies have confirmed the positive effects of demographics (age, sex, medication, and so on) on iron content in the brain, these covariates were adjusted in group comparison [24]. After analysis, the study results were warped to MNI space and viewed in XJVIEW software (http://www.alivelearn.net/xjview).

Wilcoxon rank-sum tests were used to show interhemispheric differences. Then, averaged susceptibility of ROIs was used for comparison on tissue iron data, and averaged median susceptibility of ROIs across hemispheres were computed to improve measurement stability. Moreover, depressionversus-control Wilcoxon rank-sum statistics were computed (p_Bonferroni<0.05). The change of susceptibility before and after treatment in the depression was detected by paired t test. In addition, Pearson correlation test (median ROI susceptibility with symptom severity-HDRS, relapse frequency, mean duration of single episode, and total course of disease onset) was also computed, and linear least-square fits were calculated for regions where the Pearson correlation test returned a significant result. All tests were performed as two-tailed test.

RESULTS

Individual comparison

Definite iron deposition in the brain were found in the basal ganglia and other deep-brain nuclei (the globus pallidus had a highest susceptibility), which were not so remarkable in the cortex. Moreover, the venous system and some small scattered signal inhomogeneities were found in maps. In addition, comparing with the controls, patients with recurrent depression showed significantly increased susceptibilities in the bilateral putamen and thalamus (Figure 1).

Whole-brain quantitative susceptibility mapping study

The whole-brain group results revealed that iron deposition under depression is spatially selective (Figure 2). Besides, the whole-brain QSM results revealed widespread absolute susceptibility increases in patients with depression, which involved both the cerebellar and cerebral structures (FWE corrected p<0.05). Largely bilateral abnormalities were found in the frontal lobes, temporal lobe structures, occipital lobes hippocampal regions, putamen, thalamus, cingulum, and cerebellum of the patient group.

Regions of interest extraction and validation

The average susceptibility of ROIs from the control group in this study was compared with the iron concentration in postmortem samples as reported by a previous study [23]. Although there were significant differences between the two groups that might be attributable to some objective reasons (the former was dominated by ferritin macromolecules, myelin, blood, calcium and trace elements; the latter was dominated by non-haemin iron), a linear regression was then applied to find the correlation between susceptibility values of controls in the present study and iron content in postmortem samples from Hallgren and Sourander [23], and a significant correlation (r=0.783, p<0.001) was found in the frontal lobe, temporal lobe, occipital lobe, cerebellum, putamen, and thalamus (Supplementary Figure 1 in the online-only Data Supplement).

Regional quantitative susceptibility mapping study

Regional QSM study showed that left/right measurements were comparable across controls and patients with depression in both the deep gray matter (DGM) structures and cortex (all p>0.140). Compared with the control group, the depression group showed significantly increased regional susceptibility in the frontal lobe (Z=3.6, p<0.001), temporal lobe (Z=4.2, p<0.001), occipital lobe (Z= 4.3, p<0.001), hippocampus (Z=3.0, p=0.007), putamen (Z=3.4, p<0.001), thalamus (Z=3.3, p<0.001), and cingulum (Z=2.3, p=0.01), as well as in the cerebellum (Z=2.8, p=0.004) (Supplementary Figure 2 in the online-only Data Supplement).

Although HDRS showed significant differences before and after treatment (t=16.295, p<0.001), there was no susceptibility change was found in the involved brain regions (Supplementary Table 1 in the online-only Data Supplement), and no significant correlation between HDRS and susceptibility was found (all p>0.05). Subsequent results revealed a nearly statistically significant association between the relapse frequency and the susceptibility in the putamen, thalamus, and hippocampus (Figure 3), and similar results were also found in the analysis of mean duration of single episode (Figure 4). Further studies showed that susceptibility across the putamen (ρ2=0.27, p<0.001), thalamus (ρ2=0.21, p<0.001), and hippocampus (ρ2=0.19, p<0.001) were strongly correlated with the total course of disease onset (Figure 5).

DISCUSSION

This study was the first to report the relationships between brain iron deposition and multiple factors associated with depression (treatment, symptom severity-HDRS, relapse frequency, mean duration of single episode, and total course of disease onset). Our study discovered widespread increased susceptibility across the basal ganglia and cortex in depression, which could not be reversed by treatment. More important was the finding that susceptibility across the putamen, hippocampus, and thalamus were strongly correlated with the total course of disease onset, but not the severity of depressive status-HDRS. It is suggested that brain iron deposition may be a cumulative marker of brain damage in depression, but not an immediate indicator of symptoms.

Before the QSM method was successful developed [25,26], there were two types of methods to detect iron in the brain, i.e., qualitative MRI methods (T2 [27], T2* [28], R2* [29], susceptibility weighted imaging [SWI]) [30], and semi-quantitative MRI methods (magnetic field correlation imaging [31], phase imaging) [32]. However, all the above-mentioned methods are essentially based on the characteristics of the magnetic field, which may cause blooming artifacts [33,34], so the traditional MRI quantification methods do not provide an absolute quantification and are under the effect of non-local influences [35]. Compared to the traditional techniques, QSM is believed to be able to give a more accurate and specific measurement on tissue magnetic susceptibility [36,37]. This susceptibility can be calculated using a map of the resonance frequency in each voxel, which traditionally utilizes the MR phase signal from GRE imaging and has been shown to correlate well with tissue iron concentration in most gray matter regions in the brain [38].

As iron deposition in the brain is spatially selective (significant brain iron deposition was found in the basal ganglia and other deep-brain nuclei), these ROIs are usually subjectively chosen, artificially outlined, extracted, and analyzed directly in most previous MRI researches. In fact, there are two major deficiencies with this strategy. First, it usually ignores the changes in the most important structure of the brain (the cortex). Abnormal deposition of iron in the cerebral cortex has been proved in vitro by other established methods, such as histochemical analysis, staining or transcranial sonography. Second, the extraction of artificially outlined ROIs is subject to two definite limitations: i) results are generally prone to inaccuracies in ROI definition due to the intricate anatomical substructures and weak contrast of anatomical boundaries on QSM; and ii) the poor consistency of repeated measurements depending on subjective outlining might cause a groupdependent study bias. In need of special is that, although ROIs extracted according to the MNI-template in FSL could overcome the above deficiencies, there was still limitation of the method adopted in this study that the registration might be not perfect because of lost or mixed image information.

In our study, consistent results from both a whole-brain and regional QSM study showed that brain regions with abnormal susceptibility in patients with depression included the frontal lobes, temporal lobes, occipital lobes, hippocampal regions, putamen, thalamus, cingulum, and cerebellum. We found it interesting that increased brain iron deposition occurred not only in two iron-rich areas (not all iron-rich brain areas were involved), but also in areas with minimal iron content (cortex, hippocampus, and cerebellum). This phenomenon was also reported in previous QSM whole-brain studies of iron deposition in PD and multiple sclerosis patients [19,24]. Although its intrinsic mechanism remains unclear, this distribution pattern of aberrant brain iron deposition is enough to draw our attention and strongly suggest that studies on iron in the brain should not be focused merely on DGM. So far, direct evidence of the relations between iron deposition patterns in the brain and depression has been exhibited in this study. In fact, regions of aberrant brain iron deposition we reported here were consistent with those reported by structural and functional MRI studies of depression [39,40]. That is, in a state of stress or depression, brain iron deposition may be one of physiological and pathological changes, with synergistic effects on the structure and function of the sensitive brain regions.

In subsequent analysis, we found the susceptibility in the involved brain regions could not be reversed by treatment, and was not correlated with the severity of depressive status-HDRS. The results were slightly different from the previous reports [41], which might be explained by two reasons: 1) different sample populations and 2) different methods of brain iron deposition analysis. Our further results found that not only the susceptibility in the putamen, thalamus and hippocampus presented a nearly statistically significant association with the relapse frequency and mean duration of single episode, but also a strongly correlated with the total course of disease onset. In fact, the relapse frequency and mean duration of single episode were the two determinations of the total course of disease onset. We believed that these three factors should be all closely related to the brain iron deposition in the above brain regions, and “a nearly statistically significant” might be attributed to the limited sample size in this study. Besides, it was not surprising to find that these regions of brain iron deposition were associated with the depressionrelated conditions, as previous studies had revealed that the putamen (part of the striatum) is related to motor and cognitive functions [42], and its dysfunction is known to be associated with loss of dopaminergic neurons within the corticalstriatum-thalamocortical circuit [4]. The thalamus is thought to be involved in the pathophysiology of depressive disorders and is currently drawing sustaining attention [43], and a recent research has reported that abnormal thalamocortical connectivity was found in depression [44]. Another brain region involved was the well-studied hippocampus, which has always been regarded as a targeted area in depression research because of its important role in emotion, cognition and memory. In general, brain iron deposition in the specific brain regions may be a cumulative marker of brain damage in depression, but not an immediate indicator of symptom.

We sincerely apologize for the shortcomings of the paper. 1) A limited sample size in this study could lead to statistical errors of type I and II. 2) We had not adopted a finer scale for the selection of cortical interest areas. 3) This study did not have sufficient follow-up to understand the further evolution of brain iron deposition in depression.

In conclusions, brain iron deposition is related to the total course of disease onset, but not the severity of depression, which suggest that brain iron deposition may be a sign of brain damage in multiple recurrent depression.

Supplementary Materials

The online-only Data Supplement is available with this article at https://doi.org/10.30773/pi.2022.0110.

Supplementary Table 1.

Comparison before and after treatment

pi-2022-0110-suppl1.pdf

Supplementary Figure 1.

The differences and relationship between average QSM value and iron concentration. A: Differences between averaged susceptibility of brain regions from controls in this study and the iron concentration in postmortem samples reported in Hallgren and Sourander [23]. B: The relationship between averaged susceptibility of brain regions in controls in this study and the iron concentration in postmortem samples reported in Hallgren and Sourander [23]. **p<0.01. 1, frontal lobe; 2, temporal lobe; 3, occipital lobe; 4, cerebellum; 5, putamen; 6, thalamus; QSM, quantitative susceptibility mapping.

pi-2022-0110-suppl2.pdf

Supplementary Table 1.

Comparison between groups (mean±standard error of the mean). *signifies a group difference of p_Bonferroni<0.05; **signifies a group difference of p_Bonferroni<0.01; ***signifies a group difference of p_Bonferroni<0.001.

pi-2022-0110-suppl3.pdf

Notes

Availability of Data and Material

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

Conflicts of Interest

The authors have no potential conflicts of interest to disclose.

Author Contributions

Conceptualization: Xinxiu Duan, Yuhang Xie, Lei Li. Data curation: Xiefang Zhu, Lei Chen. Formal analysis: Xiufang Zhu, Lei Chen, Feng Li. Methodology: Yuhang Xie, Xiufang Zhu, Lei Chen. Project administration: Lei Chen, Feng Li, Guoquan Feng. Resources: Lei Chen, Feng Li, Guoquan Feng, Lei Li. Supervision: Lei Li. Validation: Guoquan Feng. Visualization: Lei Chen, Feng Li, Guoquan Feng. Writing—original draft: Xinxiu Duan. Writing—review & editing: all authors.

Funding Statement

This work was supported by Project of Jiangsu Health Committee (2017090), project of Zhenjiang social development (SH2019083, SH2019063, SH2018030).

ACKNOWLEDGEMENTS

We would like to thank all respondents and participants involved in data collection and processing for this study.

Figure 1.

Patient with depression versus the control. Bright regions reflect the strong paramagnetism returned by metallic species such as tissue iron. Increased susceptibility in the bilateral putamen and thalamus could be directly found.

Figure 2.

Results of the voxel-based analysis of corrected susceptibility between the control groups and patients with depression. The results were overlaid onto the MINI274 brain template and viewed in XJVIEW (http://www.alivelearn.net/xjview). Red/yellow clusters represent T values at pFWE<0.05. QSM, quantitative susceptibility mapping.

Figure 3.

Median susceptibility for patients with depression plotted in relation to the relapse frequency. Dashed lines represent 95% confidence intervals.

Figure 4.

Median susceptibility for patients with depression plotted in relation to the mean duration of single episode. Dashed lines represent 95% confidence intervals.

Figure 5.

Median susceptibility for patients with depression plotted in relation to the total course of disease onset. Dashed lines represent 95% confidence intervals.

Table 1.

Study demographic details

	Control (N=20)	Depression (N=21)	p
Sex (male/female)	9/11	9/12	NS
Age (yr)	37.7±7.3	39.4±4.2	NS
Handedness (right/left)	21/0	21/0	NS
Education (yr)	9.4±3.5	9.7±2.1	NS
Family history	NA	NA
Outpatient/inpatient	NA	0/21
Unipolar/bipolar	NA	21/0
First-episode/relapse	NA	0/21
HDRS	NA	27.2±5.3
Relapse frequency (times)	NA	3.3±0.8
Mean duration of single episode (wk)	NA	8.4±1.1
Total course of disease onset (wk)	NA	24.4±2.6

Values are given mean±standard deviation. Total course of disease onset=cumulative single episode. HDRS, Hamilton Depression Rating Scale; NA, not applicable; NS, not significant (p>0.05)

REFERENCES

1. Moylan S, Maes M, Wray NR, Berk M. The neuroprogressive nature of major depressive disorder: pathways to disease evolution and resistance, and therapeutic implications. Mol Psychiatry 2013;18:595-606.

2. Li Y, Yan J, Zhu X, Zhu Y, Yao S, Xu Y, et al. Dilated Virchow-Robin spaces in the hippocampus impact behaviors and effects of anti-depressant treatment in model of depressed rats. J Affect Disord 2017;219:17-24.

3. Zhang N, Qin J, Yan J, Zhu Y, Xu Y, Zhu X, et al. Increased ASL-CBF in the right amygdala predicts the first onset of depression in healthy young first-degree relatives of patients with major depression. J Cereb Blood Flow Metab 2020;40:54-66.

4. Liu C, Liang MC, Soong TW. Nitric oxide, iron and neurodegeneration. Front Neurosci 2019;13:114

5. Muhoberac BB, Vidal R. Iron, ferritin, hereditary ferritinopathy, and neurodegeneration. Front Neurosci 2019;13:1195

6. Burhans MS, Dailey C, Beard Z, Wiesinger J, Murray-Kolb L, Jones BC, et al. Iron deficiency: differential effects on monoamine transporters. Nutr Neurosci 2005;8:31-38.

7. Bagnato F, Hametner S, Yao B, van Gelderen P, Merkle H, Cantor FK, et al. Tracking iron in multiple sclerosis: a combined imaging and histopathological study at 7 Tesla. Brain 2011;134(Pt 12):3602-3615.

8. Dixon SJ, Stockwell BR. The role of iron and reactive oxygen species in cell death. Nat Chem Biol 2014;10:9-17.

9. Haacke EM, Cheng NY, House MJ, Liu Q, Neelavalli J, Ogg RJ, et al. Imaging iron stores in the brain using magnetic resonance imaging. Magn Reson Imaging 2005;23:1-25.

10. Tambasco N, Paolini Paoletti F, Chiappiniello A, Lisetti V, Nigro P, Eusebi P, et al. T2*-weighted MRI values correlate with motor and cognitive dysfunction in Parkinson’s disease. Neurobiol Aging 2019;80:91-98.

11. Schweser F, Deistung A, Lehr BW, Reichenbach JR. Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: an approach to in vivo brain iron metabolism? Neuroimage 2011;54:2789-2807.

12. Wang Y, Spincemaille P, Liu Z, Dimov A, Deh K, Li J, et al. Clinical quantitative susceptibility mapping (QSM): biometal imaging and its emerging roles in patient care. J Magn Reson Imaging 2017;46:951-971.

13. de Rochefort L, Liu T, Kressler B, Liu J, Spincemaille P, Lebon V, et al. Quantitative susceptibility map reconstruction from MR phase data using bayesian regularization: validation and application to brain imaging. Magn Reson Med 2010;63:194-206.

14. Gong NJ, Dibb R, Bulk M, van der Weerd L, Liu C. Imaging beta amyloid aggregation and iron accumulation in Alzheimer’s disease using quantitative susceptibility mapping MRI. Neuroimage 2019;191:176-185.

15. Xiao B, He N, Wang Q, Cheng Z, Jiao Y, Haacke EM, et al. Quantitative susceptibility mapping based hybrid feature extraction for diagnosis of Parkinson’s disease. Neuroimage Clin 2019;24:102070

16. Smith SM. Fast robust automated brain extraction. Hum Brain Mapp 2002;17:143-155.

17. Schofield MA, Zhu Y. Fast phase unwrapping algorithm for interferometric applications. Opt Lett 2003;28:1194-1196.

18. Tang J, Liu S, Neelavalli J, Cheng YC, Buch S, Haacke EM. Improving susceptibility mapping using a threshold-based K-space/image domain iterative reconstruction approach. Magn Reson Med 2013;69:1396-1407.

19. Acosta-Cabronero J, Cardenas-Blanco A, Betts MJ, Butryn M, Valdes-Herrera JP, Galazky I, et al. The whole-brain pattern of magnetic susceptibility perturbations in Parkinson’s disease. Brain 2017;140:118-131.

20. Haacke EM, Tang J, Neelavalli J, Cheng YC. Susceptibility mapping as a means to visualize veins and quantify oxygen saturation. J Magn Reson Imaging 2010;32:663-676.

21. Betts MJ, Acosta-Cabronero J, Cardenas-Blanco A, Nestor PJ, Düzel E. High-resolution characterisation of the aging brain using simultaneous quantitative susceptibility mapping (QSM) and R2* measurements at 7T. Neuroimage 2016;138:43-63.

22. Fan L, Li H, Zhuo J, Zhang Y, Wang J, Chen L, et al. The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb Cortex 2016;26:3508-3526.

23. Hallgren B, Sourander P. The effect of age on the non-haemin iron in the human brain. J Neurochem 1958;3:41-51.

24. Acosta-Cabronero J, Betts MJ, Cardenas-Blanco A, Yang S, Nestor PJ. In vivo MRI mapping of brain iron deposition across the adult lifespan. J Neurosci 2016;36:364-374.

25. Liu W, Ge T, Leng Y, Pan Z, Fan J, Yang W, et al. The role of neural plasticity in depression: from hippocampus to prefrontal cortex. Neural Plast 2017;2017:6871089

26. Stüber C, Pitt D, Wang Y. Iron in multiple sclerosis and its noninvasive imaging with quantitative susceptibility mapping. Int J Mol Sci 2016;17:100

27. Vymazal J, Brooks RA, Zak O, McRill C, Shen C, Di Chiro G. T1 and T2 of ferritin at different field strengths: effect on MRI. Magn Reson Med 1992;27:368-374.

28. Chavhan GB, Babyn PS, Thomas B, Shroff MM, Haacke EM. Principles, techniques, and applications of T2*-based MR imaging and its special applications. Radiographics 2009;29:1433-1449.

29. Deistung A, Schäfer A, Schweser F, Biedermann U, Turner R, Reichenbach JR. Toward in vivo histology: a comparison of quantitative susceptibility mapping (QSM) with magnitude-, phase-, and R2*-imaging at ultra-high magnetic field strength. Neuroimage 2013;65:299-314.

30. Haacke EM, Makki M, Ge Y, Maheshwari M, Sehgal V, Hu J, et al. Characterizing iron deposition in multiple sclerosis lesions using susceptibility weighted imaging. J Magn Reson Imaging 2009;29:537-544.

31. Ropele S, de Graaf W, Khalil M, Wattjes MP, Langkammer C, Rocca MA, et al. MRI assessment of iron deposition in multiple sclerosis. J Magn Reson Imaging 2011;34:13-21.

32. Ogg RJ, Langston JW, Haacke EM, Steen RG, Taylor JS. The correlation between phase shifts in gradient-echo MR images and regional brain iron concentration. Magn Reson Imaging 1999;17:1141-1148.

33. Jensen JH, Chandra R, Ramani A, Lu H, Johnson G, Lee SP, et al. Magnetic field correlation imaging. Magn Reson Med 2006;55:1350-1361.

34. Li J, Chang S, Liu T, Wang Q, Cui D, Chen X, et al. Reducing the object orientation dependence of susceptibility effects in gradient echo MRI through quantitative susceptibility mapping. Magn Reson Med 2012;68:1563-1569.

35. Yablonskiy DA, Haacke EM. Theory of NMR signal behavior in magnetically inhomogeneous tissues: the static dephasing regime. Magn Reson Med 1994;32:749-763.

36. Lim IA, Faria AV, Li X, Hsu JT, Airan RD, Mori S, et al. Human brain atlas for automated region of interest selection in quantitative susceptibility mapping: application to determine iron content in deep gray matter structures. Neuroimage 2013;82:449-469.

37. Schweser F, Sommer K, Deistung A, Reichenbach JR. Quantitative susceptibility mapping for investigating subtle susceptibility variations in the human brain. Neuroimage 2012;62:2083-2100.

38. Schenck JF. The role of magnetic susceptibility in magnetic resonance imaging: MRI magnetic compatibility of the first and second kinds. Med Phys 1996;23:815-850.

39. Zhao K, Liu H, Yan R, Hua L, Chen Y, Shi J, et al. Altered patterns of association between cortical thickness and subcortical volume in patients with first episode major depressive disorder: a structural MRI study. Psychiatry Res Neuroimaging 2017;260:16-22.

40. Sacchet MD, Camacho MC, Livermore EE, Thomas EAC, Gotlib IH. Accelerated aging of the putamen in patients with major depressive disorder. J Psychiatry Neurosci 2017;42:164-171.

41. Yao S, Zhong Y, Xu Y, Qin J, Zhang N, Zhu X, et al. Quantitative susceptibility mapping reveals an association between brain iron load and depression severity. Front Hum Neurosci 2017;11:442

42. Herrero MT, Barcia C, Navarro JM. Functional anatomy of thalamus and basal ganglia. Childs Nerv Syst 2002;18:386-404.

43. Kong L, Chen K, Womer F, Jiang W, Luo X, Driesen N, et al. Sex differences of gray matter morphology in cortico-limbic-striatal neural system in major depressive disorder. J Psychiatr Res 2013;47:733-739.

44. Brown EC, Clark DL, Hassel S, MacQueen G, Ramasubbu R. Thalamocortical connectivity in major depressive disorder. J Affect Disord 2017;217:125-131.