Abstract
Background:
Diffusion magnetic resonance imaging (dMRI) has revealed microstructural changes in white matter (WM) in Huntington’s disease (HD).
Objective:
To compare the validities of different dMRI, i.e., diffusion kurtosis imaging (DKI) and diffusion tensor imaging (DTI) in HD.
Methods:
22 mutant
Results:
Compared to controls, m
Conclusions:
Microstructural changes in WM occur at early stage of HD and deteriorate over the disease progression. Integrating DKI and DTI would provide the best accuracies for differentiating early HD from control and identifying advanced HD.
Keywords
INTRODUCTION
Huntington’s disease (HD) is an autosomal dominant inherited neurodegenerative disorder caused by mutation in huntingtin gene (
To estimate the progression of HD pathology as a function of CAG repeats and time of exposure to the effects of the repeats, Penney et al. in 1997 proposed CAG-Age Product (CAP) scores [3]. CAP has been examined in various studies, including TRACK-HD and PREDICT-HD [4], as it serves as a representative marker for the disease burden. The CAP formula, used in this paper, was proposed by Zhang et al. and categorizes individuals into three groups: Low, Medium, and High, with the group thresholds set at CAP < 290 (Low), 290≤CAP≤368 (Medium), and CAP > 368 (High) [4]. HD can be divided into Pre-manifest HD (PreHD) and manifest HD (mHD) depending on whether the patient has characteristic motor symptoms. Symptoms can manifest in diverse ways at early stage of the disease, with signs potentially escaping detection during a standard neurological examination, thus requiring an experienced clinician’s expertise. These challenges in pinpointing the exact onset of HD’s motor symptoms, coupled with the fact that pathological changes and other symptoms (e.g., cognitive impairment and behavioral changes) may commence years before clinical diagnosis [5, 6], underscore the need for developing biomarkers with capability to accurately portray the pathological changes, predict the disease stages, as well as measure the disease progression in individuals carrying the abnormal gene. Neuropathological changes in HD are notably selective, characterized by significant cell loss and brain atrophy, particularly in the caudate and putamen [7, 8]. Multimodal magnetic resonance imaging (MRI) offers a non-invasive means to observe the pathological alterations in HD [9, 10]. While neuroimaging research has predominantly focused on grey matter volume loss, particularly in grey matter and subcortical structures, it is now evident that both white matter volume loss [11–13] and modified white matter diffusion metrics [14–19] are affected in both PreHD and mHD. Diffusion magnetic resonance imaging (dMRI) is a technique used to characterize the diffusion of water molecules within tissue, including diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), bi-tensor DTI, and neurite orientation dispersion and density imaging (NODDI) [20]. Fractional anisotropy (FA) is a parameter of DTI, which reflects the degree of anisotropy of the dispersion of water molecules. A recent meta-analysis study showed that both PreHD and mHD patients demonstrated significant decreased FA values in the corpus callosum enriched with white matter when compared with that in controls [21].
DKI stands out as a highly promising diffusion imaging technique [22, 23], extending its reach to quantify non-Gaussian water diffusion. Recent research has showcased the superior sensitivity of DKI measures in detecting developmental and pathological microstructural changes in neuronal tissues compared to conventional DTI [24, 25]. The mean kurtosis (MK) value is the major parameter in DKI, reflecting the average dispersion kurtosis of water molecules in all directions, with lower MK values indicating reduced microstructural complexity within the region of interest. Studies applying DKI technique in other neurodegenerative diseases, such as mild cognitive impairment (MCI), Alzheimer’s disease (AD) [26–29], and Parkinson’s disease (PD) [30, 31] have showed that microstructural changes in cerebral white matter and suggested neuronal loss in correspondent cortical regions. However, there is no study using DKI technique in patients with HD and only two animal-based studies showing significant microstructural changes in HD models [32, 33]. Both studies indicate that DKI is a sensitive method for detecting HD-associated white matter abnormalities [32, 33]. Compared to NODDI, DKI can capture information about non-Gaussian diffusion of tissues, which give rise to better sensitivity to detect more complex tissue microstructural changes in neurodegenerative diseases [20].
In this study, we aim to use DKI technology to study the microstructural changes of white matter fiber bundles in HD patients, compare DKI and DTI in detecting the microstructural changes of white matter in HD, and further explore the correlations between DKI imaging parameters and various clinical manifestations.
MATERIALS AND METHODS
Subjects
Twenty-six subjects with more than 35 CAG repeats were enrolled from the Huntington Disease Center at Beijing Tiantan Hospital, Capital Medical University, Between May 8, 2019 and February 8, 2023. They underwent multimodal MRI scans. In addition, we recruited 16 age- and sex-matched healthy controls during the same period. Ethical approval for this study was obtained from the ethics committee of the hospital. The inclusion and exclusion criteria are outlined below:
Inclusion criteria for mutation
Inclusion criteria for control group: 1) Age≥18 years; 2) without genetic or clinical evidence for the presence of a neurological disorder; 3) signed the informed consent; 4) sex- and age-matched with m
Exclusion criteria: 1) History of severe stroke, brain trauma, epilepsy, or mental illness; 2) History of brain tumor and space occupation; 3) History of other central nervous system degenerative diseases; 4) Patients with other diseases who can’t complete the magnetic resonance examination; 5) Refusal to sign informed consent.
All study procedures were conducted according to the declaration of Helsinki. All participants provided written informed consent approved by the Ethics Board of the Beijing Tiantan Hospital, Capital Medical University of China.
Clinical assessments
Detail medical history data from all participants were recorded, including basic information, height, weight, body mass index (BMI), educational background, onset time, family history, duration of motor, cognitive, and psychiatric symptoms, and the number of CAG repeats in the
We also conducted the Unified Huntington’s Disease Rating Scale (UHDRS) [36] assessments for all enrolled subjects, including Total Motor Score (TMS), Total Functional Capacity (TFC) [36], Functional Checklist Score (FCS), and independence scale. Depression was assessed using Beck Depression Inventory II (BDI-II) [37]. Cognitive changes were assessed by the Stroop Interference Test (STIT), Symbol Digit Modality Test (SDMT), Category Fluency Test of animals (CFT), and Mini-Mental State Examination (MMSE) [38]. Our clinician team had the above clinical rating scales training and certificated by Enroll HD program (https://enroll-hd.org/). Clinical assessments of participants in this study were undertaken independently by two clinicians, and consolidated by a third party (YH) in case of disparities. Demographic information and clinical assessments were also collected for all control participants in this study.
Imaging data acquisition
We conducted MRI scans on a Siemens Prisma 3T MRI scanner equipped with a 64-channel head coil. During the scanning process, participants remained awake with their eyes closed and heads immobilized, ensuring minimal motion. According to previous study [20], three b values of 0, 1000 and 2000 s/mm2 were usually selected for collection in DKI. In this study, b values of 0, 1000 and 2000 s/mm2 were adopted. Scanning parameters for each sequence were as follows:
High resolution T1 scanning parameters: Voxel size: 0.9×0.9×1.0 mm; Field of View (FoV) read: 240 mm; FoV phase: 100%; Repetition time (TR): 1560 ms; Echo time (TE): 1.69 ms; Slice resolution: 100%; Flip angle: 8 deg; Sagittal slices: 176; Acquisition Time (TA): 3 min 37 s.
DKI was acquired with 129 gradient directions and b values of 0 s/mm2, 1000, s/mm2, and 2000 s/mm2. Voxel size: 2.5 * 2.5 mm; FoV read: 220 mm; FoV phase: 100%; TR: 2900 ms; TE: 64 ms, TA: 6 min 40 s. Then, 65 gradient directions of b values of 0 and 1000 s/mm2 are extracted for DTI analysis.
Image data processing and quality control
Every MRI scan underwent visual review to exclude subjects with intracranial masses or other abnormalities. Scans with missing sequences or motion artifacts during DKI and DTI were excluded from analysis. We employed FMRIB’s Software Library (FSL) [39] for DKI and DTI measurements, involving brain extraction (BET) and correction for motion and eddy current-induced distortions through affine registration to a reference
Statistics and analysis
Data following a normal distribution were presented as mean±standard deviation, while data not adhering to a normal distribution were expressed as median±interquartile range. The Shapiro-Wilk test was used to test the normality of the data. For normally distributed data, we employed T-tests and one-way ANOVA to compare differences between groups. ANCOVA was used to eliminate the differences in age and gender among different groups. Non-parametric tests were used for data that did not exhibit a normal distribution. Bonferroni post hoc tests were conducted following ANOVA, as it could be applied to almost any situation with multiple comparisons. Corrected
RESULTS
Participant characteristics and clinical assessments
Following image processing and quality control, finally this research included 22 participants in the m
Demographic and clinical assessment details for the m
Comparison of clinical assessments and demographic information between three groups
SD, standard deviation; BMI, body mass index; STIT, Stroop Interference Test; SDMT, Symbol Dicit Modality Test; CFT, Category Fluency Test of animals; BDI-II, Beck Depression Inventory II; MMSE, Mini-Mental State Examination; CAP, CAG age product; TFC, total functional capacity; N/A, not applicable; *Cases with missing values were excluded.
Cognitive function evaluations, encompassing STIT, SDMT, and CFT, revealed that the cognitive function of the hCAP group was notably inferior to those of HC and m&lCAP groups (Table 1). MMSE was significantly different only between hCAP and HC groups. However, BDI-II scores did not show a significant difference between any two groups, only showed an upward trend among the three groups.
DKI and DTI analyses
According to the JCH ICBM-DTI-81 white matter atlas, which provides a comprehensive division of human white matter into 48 distinct fiber tracts, we extracted MK values from all 48 regions for all 36 subjects and examined differences in each region among the three groups. To eliminate potential confounding factors such as age and sex on white matter structure, we conducted covariance analysis.
Following ANCOVA analysis, we identified 42 regions with significant differences in MK values among the three groups (Fig. 1, Table 2, and Supplementary Table 1). Specifically, after

The number of brain regions with different MK or FA values in multiple comparisons of three groups
No., number; HC, healthy control; m& lCAP, medium and low CAP group; hCAP, high CAP group.
Correlation analysis of key white matter regions with clinical presentations
We identified key brain regions for correlation analysis based on significant differences in diffusion parameters observed among the three groups in the dMRI analysis. These parameters specifically encompassed MK values from the genu of the corpus callosum, the body of the corpus callosum, the right superior corona radiate, the right posterior corona radiate, and the left superior longitudinal fasciculus. We conducted correlation analyses with CAP, TMS, TFC, STIT, SDMT, and BDI-II (Fig. 2). Higher CAP and TMS scores indicate more severe disease burden and motor symptoms. Lower TFC, STIT, SDMT, and CFT indicate more severe functional impairment and cognitive impairment. Higher BDI-II scores reflect more pronounced depression.

Notably, the MK values of the genu of the corpus callosum exhibited negative correlations with CAP and TMS, while displaying positive correlations with TFC, STIT, and SDMT, all of which held true even after applying the FDR correction (Fig. 2). Similarly, MK values for the body of the corpus callosum demonstrated negative correlations with CAP and TMS, and positive correlations with TFC, STIT, and SDMT after FDR correction (Fig. 2). Moving on to the right superior corona radiate, MK values displayed negative correlations with CAP and TMS, as well as positive correlations with TFC, STIT, and SDMT, again even after FDR correction (Fig. 2). Similarly, MK values of the right posterior corona radiate exhibited negative correlations with CAP and TMS, and positive correlations with TFC, STIT, and SDMT after FDR correction (Fig. 2). Finally, the MK values for the left superior longitudinal fasciculus revealed negative correlations with CAP and TMS, along with positive correlations with TFC, STIT, and SDMT, all of which remained significant even after FDR correction (Fig. 2). No significant correlation was found between MK values in the above five regions and BDI-II scores
ROC and Youden Index analyses
To test whether microstructural changes in white matter can be used as an indicator to distinguish m& lCAP from HC patients and distinguish different stages of m
Brain regions with significant different FA or MK values between HC and m& lCAP were selected for ROC analysis to differentiate m& lCAP groups from controls. When using the MK values of splenium of corpus callosum, the maximum AUC of 0.848 was reached to distinguish participants with hCAP from m& lCAP, and the Youden index was 0.661 (Table 3, Supplementary Figure 2). When using the FA values of right posterior thalamic radiation, the maximum AUC of 0.920 was achieved to distinguish m& lCAP from HC, and the Youden index was 0.714 (Table 2, Supplementary Figure 2).
ROC and Youden Index analyses using MK and FA values
HC, healthy controls; m& lCAP, medium and low CAP; hCAP, high CAP; AUC, area under curve.
Brain regions with significant different FA or MK values between m& lCAP and hCAP were selected for ROC analysis to differentiate different stages of HD. When using the MK values of left superior longitudinal fasciculus, the maximum AUC of 0.911 was reached to distinguish participants with hCAP from m& lCAP, and the Youden index was 0.732 (Table 3, Supplementary Figure 2). When using the FA values of the left tapetum, the maximum AUC of 0.830 was reached to distinguish participants with hCAP from m& lCAP, and the Youden index was 0.625 (Table 3, Supplementary Figure 2).
Furthermore, we combined the FA and MK values with the best identification power to obtain an optimal differentiation effect. The MK of Splenium of corpus callosum combining with the FA of Posterior thalamic radiation R was able to distinguish m& lCAP group from HC group with AUC 0.911 and Youden index 0.786 (Table 2, Supplementary Figure 2). While combining the MK values of the left superior longitudinal fasciculus and the FA values of left tapetum would be able to distinguish hCAP group from m& lCAP group with AUC 0.902 and Youden index 0.804 (Table 2, Supplementary Figure 2).
DISCUSSION
This is the first study using DKI technique to detect microstructural changes of white matter in patients with HD. Furthermore, we explored differences in diffusion parameters among individuals with varying HD disease burdens in combination of DKI and DTI techniques. Utilizing ROI-based analysis, we identified notable white matter fibers alterations in five brain regions with the most significant changes using DKI. MK values in these regions demonstrated negative and gradient correlations with the disease severity and disease burden of HD. Importantly, apart from underscoring the superior sensitivity of DKI in detecting white matter microstructural changes compared to DTI, our findings also suggested combining FA with MK values would generate the best accuracy of identifying m
MK values in 42 regions were found to have significant differences among the three groups, and FA values in 20 regions were found to have significant differences among the three groups. Furthermore, MK values in five regions exhibited significant differences between any two groups. We did not identify any brain region that showed significant differences in FA values between any two of these three groups. Although FA values could reflect white matter microstructural changes in patients with advanced HD, it shows limited capacity to detect microstructural changes in white matter at low to moderate disease stages. These findings strongly suggest that MK values exhibit superior sensitivity in detecting white matter microstructural changes among individuals with varying degrees of HD disease burden.
The imaging characteristics of DKI can quantify non-Gaussian diffusion of water molecules in the tissue [40]. Due to the presence of various membrane structures in tissues, the distribution of water molecules in the brain tissues does not always follow a Gaussian distribution, so it has higher sensitivity for the microstructure changes of fiber bundles with more complex directions [41]. DKI is more sensitive to detect early changes in white matter microstructure in HD, similar to other neurodegenerative diseases [28, 42]. MK represents the mean value of diffusion kurtosis in all gradient directions under multiple b values, and is an index to evaluate the microstructure complexity of tissue in the region of interest. A decrease in MK values implies disruption of white matter microstructural integrity [25].
White matter damage in HD was thought secondary to the loss of grey matter in the form of Wallerian degeneration [43]. However, other study suggested white matte aberrations are a feature of HD independent of neuronal loss [44]. In an animal study, BACHD mice exhibited thinner myelin and decreased myelin compaction as early as 1 month old, suggesting that myelin abnormalities in HD are an early pathological event [45]. The above studies indicated that white matter changes in HD may attribute to impaired myelination due to mHTT deposition in oligodendrocytes and axonal damage secondary to neuronal loss.
Previous DTI analysis has unveiled abnormalities in neuronal fiber orientation and integrity within white matter and subcortical gray matter structures in both PreHD [17, 46] and mHD [12, 15]. DTI studies have consistently revealed alterations in FA values in patients with PreHD and early mHD, with the most prominent finding being decreased FA in the corpus callosum [12, 48]. In our study, we also found the MK and FA of corpus callosum decreased in PreHD and mHD compared with HC. In particular, the MK values of Genu of corpus callosum and Body of corpus callosum were significantly different when comparing any two groups among the three groups. Our results align with previous studies showing a widespread reduction in axonal density observed in the white matter tracts, including the corpus callosum and the basal ganglia surroundings in HD using NOODI [49], and also in agreement with Roland et al. study showing the importance of the corpus callosum for the pathophysiological processes of HD [50]. Previous studies have selected several key regions of white matter for diffusion imaging studies in HD patients or m
To explore the correlations between dMRI parameters and clinical characteristics, we analyzed the clinical features of our HD cohort (Table 1). Cognitive and psychological assessments showed significant differences between individuals in the hCAP group and other two groups. There was a noticeable trend toward cognitive decline between m& lCAP group and HC, as evidenced by the STIT, CFT and MMSE results (Table 1), consistent with previous studies showing cognitive decline and psychological abnormalities often appear during the early stage [51, 52]. In contrast to previous reports [53–55], we did not find significant differences in BMI among the three groups. These disparities could be attributed to the smaller number of subjects in the final analysis, with most subjects who completed MRI acquisitions being at an early stage of HD. Furthermore, we confirmed that the relationship of white matter microstructure changes associated with the overall disease burden, motor and cognitive functions, but not with the severity of depression, which was likely caused by biochemical or hormonal imbalance, rather than deficits in brain microstructural connectivity. The MK values of the above five regions were more significantly correlated with the disease burden (
CAP can serve as an indicator of disease burden and mHTT exposure in HD patients. TMS scores gauge motor symptoms and TFC can reflect the disease stage [57]. SDMT, a recognized vital marker for determining the disease stage of m
The ROC analysis using both MK and FA values could distinguish the m& lCAP group from HC group, as well as the hCAP group from the m& lCAP group. The results revealed that the MK value of the left superior longitudinal fasciculus showing the maximum AUC for differentiating the hCAP group from m& lCAP group, while the FA value of the right posterior thalamic radiation showing the maximum AUC for distinguishing m& lCAP group from HC. The ROC analysis demonstrated higher sensitivity when MK values were employed compared to FA values. These findings confirmed the high sensitivity of the DKI technique in terms of white matter microstructural changes in HD patients, although the specificities are secondary to FA evaluation by DTI. Our results demonstrated that integrative application of DKI and DTI technology is of great value for accurate staging differentiation of HD. Previous studies only compared the AUC area obtained by DKI parameter and DTI parameter in ROC analysis of disease diagnose [31, 61], and most of them concluded that DKI parameter has a larger AUC area, but did not further clarify the Youden index with its optimal cutoff value.
Our study first integrates the sensitive MK values with the specific FA values to achieve the largest Youden index for actual clinical practice, which means greater accuracy than either FA or MK when distinguishing advanced or prodromal HD using cutoff values. Early diagnosis of neurodegenerative diseases is of great significance for early disease intervention. At present, there is still no consensus on the staging and classification of HD [62]. HD-ISS classification recently proposed that HD should be regarded as a continuous state of the disease, rather than simply divided into manifest stage and pre-manifest stage [58]. The use of various non-invasive methods to accurately determine the severity and stages of HD patients would be the key for early diagnosis of HD. Our study provided an ideal diagnostic method in the clinical practice, which is a convenient, non-invasive and efficient approach through combination of various advanced dMRI imaging.
Our cross-sectional exploratory study may have its limitations. Apart from small sample size, our utilization of the ROI-based analysis method, while valuable, may not surpass the precision of other techniques such as tract-based spatial statistics (TBSS) [63] or fixel-based analysis (FBA) [64] in calculating diffusion parameters to detect subtle variations within the white matter skeleton. More diffusion parameters such as mean diffusion rate, radial diffusion rate, axial diffusion rate, radial kurtosis, axial kurtosis, and kurtosis anisotropy could be applied in a large longitudinal HD cohort.
In summary, we identified more extensive white matter structural alterations in HD compared to HC using DKI than DTI analysis. In addition, the white matter alterations detected by DKI technique exacerbated along the disease course of HD. Our data indicated that MK values could be used as biomarkers for monitoring the disease progression and therapeutic outcomes, and ameliorating white matter microstructural damaging might slow down the progression of HD. Furthermore, we proposed multimodal integration of diffusion MRI for better staging determination of subjects carrying m
Footnotes
ACKNOWLEDGMENTS
We would like to thank Dr. Ling-Xiao Cao, Dr. Wen-Zheng Hu, Ms. Qing Yang, Dr. Gang Du for their efforts during patient recruitment. Particularly, we thank the HD patients and healthy controls for participating in this study.
FUNDING
This project is funded by Beijing E-town Cooperation & Development Foundation (YCXJ-JZ-2022-007, YH) and National Natural Science Foundation of China (NSFC, 82071417, YH).
CONFLICT OF INTEREST
The authors have no conflict of interest to report.
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.
