In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging

In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging Vasileios Ioakeimidis , Marco Palombo , View ORCID Profile Chiara Casella , Lucy Layland , Carolyn B. McNabb , Robin Schubert , Philip Pallmann , Monica E. Busse , Cheney J. G. Drew , Sundus Alusi , Timothy Harrower , Jane Davies , View ORCID Profile Anne E. Rosser , View ORCID Profile Claudia Metzler-Baddeley doi: https://doi.org/10.1101/2025.03.17.25324107 Vasileios Ioakeimidis 1 Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University , , , , UK 2 Danish Research Centre for Magnetic Resonance, Department for Radiology and Nuclear Medicine, Copenhagen University Hospital Amager and Hvidovre , Copenhagen, Denmark Find this author on Google Scholar Find this author on PubMed Search for this author on this site Marco Palombo 1 Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University , , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Chiara Casella 3 Department for Forensic and Neurodevelopmental Sciences, Institute of Psychiatry, Psychology and Neuroscience, King’s College London , London, UK 4 Early life imaging research department, School of Biomedical Engineering and Imaging Sciences, King’s College London , London, UK 5 London Collaborative Ultra high field System (LoCUS), Kings College London , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Chiara Casella Lucy Layland 1 Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University , , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Carolyn B. McNabb 1 Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University , , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Robin Schubert 6 George Huntington Institut (GHI) , , , Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Philip Pallmann 7 Centre for Trials Research, School of Medicine, Cardiff University , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Monica E. Busse 7 Centre for Trials Research, School of Medicine, Cardiff University , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Cheney J. G. Drew 7 Centre for Trials Research, School of Medicine, Cardiff University , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sundus Alusi 8 The Walton Centre for Neurology and Neurosurgery , , , Liverpool, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Timothy Harrower 9 Royal Devon and Exeter NHS Trust , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jane Davies 10 Cardiff and Vale University Health Board, Main University Hospital Wales Building, Cardiff University , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Anne E. Rosser 11 Cardiff University Brain Repair Group, School of Biosciences, Cardiff University , , , UK 12 Advanced Neurotherapeutics Centre (ANTC), Department of Neurology and Psychological Medicine, School of Medicine, Cardiff University , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Anne E. Rosser Claudia Metzler-Baddeley 1 Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University , , , , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Claudia Metzler-Baddeley For correspondence: metzler-baddeleyc{at}cardiff.ac.uk Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract Background Huntington’s Disease (HD) is an inherited neurodegenerative disorder characterised by progressive cognitive and motor decline resulting from atrophy within basal ganglia networks. Although no disease-modifying therapies currently exist, several novel clinical trials are ongoing. Sensitive non-invasive imaging biomarkers are therefore essential for evaluating therapeutic effects. Soma and Neurite Density Imaging (SANDI), a multi-shell diffusion-weighted imaging model, estimates intracellular signal fractions arising from sphere-shaped soma that show promise as proxies for HD-related neurodegeneration. Although HD is rare, it offers a valuable model for understanding other neurodegenerative diseases due to its clear genetic cause and shared patterns of protein abnormalities. Objective To characterise HD-related microstructural abnormalities in the basal ganglia and thalami using SANDI and examine associations between SANDI indices, volumetric measurements, and motor performance. Methods T1-weighted anatomical and multi-shell diffusion-weighted images (b-values: 200 s/mm²– 6,000 s/mm²) were acquired using a 3T Siemens Connectom scanner (300mT/m) in 56 HD individuals (Mean Age = 46.1, SD Age = 13.8, 25 females) and 57 healthy controls (Mean Age = 45.0, SD Age = 13.8, 31 females). HD participants completed Quantitative Motor (Q-Motor) tasks, including speeded and paced finger tapping, which were reduced to one principal component of motor performance. Following standard diffusion-weighted data preprocessing, SANDI and diffusion tensor models estimated apparent soma density, apparent soma size, apparent neurite density, extracellular signal fraction, fractional anisotropy, and mean diffusivity. The caudate, putamen, pallidum, and thalamus were segmented bilaterally, and micro-structural and volumetric indices were extracted and compared. Correlations between SANDI indices, Q-Motor performance, and volumetric measures were analysed. Results HD was associated with reduced apparent soma density ( r rb = 0.32, p ≤ 0.007) and increased apparent soma size ( r rb = 0.45, p < 0.001) and extracellular signal fraction ( r rb = 0.34, p ≤ 0.003) in the basal ganglia, but not the thalami. These differences were more pronounced at HD-Integrated Staging System 0-1 than 2-3. No differences were found in apparent neurite density ( r rb = 0.18, p = 0.17). HD-related increases in fractional anisotropy and mean diffusivity in the basal ganglia were replicated. Q-Motor component scores correlated negatively with apparent soma density and positively with apparent soma size and extracellular signal fraction. SANDI indices and age explained up to 63% of striatal atrophy in HD. Conclusion SANDI measures detected HD-related neurodegeneration in the striatum, accounted significantly for striatal atrophy, and correlated with motor impairments. Reduced apparent soma density and increased apparent soma size align with ex vivo evidence of medium spiny neuron loss and glial reactivity. SANDI shows promise as an in vivo biomarker and surrogate outcome measure for clinical trials of disease-modifying therapies for HD and other neurodegenerative diseases. Introduction Huntington’s disease (HD) is an autosomal dominantly inherited neurodegenerative disorder caused by a pathogenic CAG repeat expansion of the Huntingtin gene. 1 HD is characterised by a progressive loss of cognitive and motor functions as well as psychiatric disturbances. The clinical onset of HD is commonly defined by the manifestation of motor symptoms, such as chorea, reduced voluntary motor control, bradykinesia, and difficulty maintaining rhythmic and paced movements. 2 – 4 However, changes in the brain, notably striatal atrophy in the basal ganglia (BG) 5 , 6 may precede the motor onset by up to 24 years 7 – 10 and correlate with motor and cognitive decline. 2 , 5 , 6 , 11 – 15 While HD is rare (∼12 in 100,000), it can be seen as a model neurodegenerative disorder, due to its clear genetic cause, well-characterized disease progression, and shared features of protein abnormalities with more common disorders like Alzheimer’s and Parkinson’s disease. 5 There is presently no approved disease-modifying therapy for HD, but several clinical trials are underway to test the safety and efficacy of novel therapeutics. 16 The recent surge in potential disease-modifying targets has generated a demand for surrogate outcome measures that are sensitive to HD neuropathology and allow a mechanistic assessment of therapeutic effects on striatal neurodegeneration in a timely manner. Volumetric measurements from non-invasive MRI are known to be sensitive to disease progression 17 – 20 and have been adopted into the Huntington’s Disease Integrated Staging System (HD-ISS). 21 However, volumetric measurements do not provide information about the underlying neuropathological tissue changes that lead to striatal atrophy, such as the loss of medium spiny neurons (MSN) 22 , 23 and changes in glia cell density and morphology, including enlargement of reactive astrocytes and microglia. 22 , 24 – 26 Diffusion-weighted Imaging (DWI) is widely used to investigate brain tissue microstructure in vivo by exploiting apparent water displacement due to Brownian motion. 27 , 28 Most DWI studies in HD have used diffusion tensor imaging (DTI), which models extra-cellular water diffusion as a Gaussian tensor and measures diffusion properties such as mean diffusivity (MD) and the degree of diffusion anisotropy (fractional anisotropy; FA). 29 DTI studies in HD have consistently reported increases in MD and FA in striatal grey matter, 30 , 31 that likely arise from the selective neurodegeneration of MSN connections. Advances in multi-shell and ultra-strong gradient DWI, 32 have enabled increasingly sophisticated biophysical models that require data acquisition over a range of b-values to separate extra- from intracellular diffusion signals. 33 – 35 Several approaches that model intraneurite space with cylinders or sticks have been put forward (e.g. Composite and Restricted Model of Diffusion, CHARMED 33 ; Convex optimization modelling for microstructure informed tractography; COMMIT 36 ; Neurite Orientation Dispersion and Density Index; NODDI 34 ). NODDI studies in premanifest HD have revealed reduced apparent neurite density in white mater with localised reductions in fibre orientation dispersion in the corpus callosum and basal ganglia capsules. 37 However, despite putamen volume loss, no NODDI-based differences in striatal grey matter were detected in gene-positive individuals long before motor onset. 10 Soma And Neurite Density Imaging (SANDI) 38 is a diffusion MRI model that extends biophysical multi-compartment models to capture microstructural features of cell bodies and processes in vivo . SANDI requires multi-shell acquisition protocols with b-values over 3,000 s/mm 2 to capture intracellular restricted signal fractions. Biologically, its parameters describe the relative contributions of somas, neurites (axons, dendrites and glial processes), and extracellular space to the diffusion-weighted MRI signal: the soma signal fraction (f is ) reflects the density of neuronal and glial cell bodies, the neurite fraction (f in ) captures density of axons, dendrites and glial processes, and the extracellular fraction (f extra ) and extracellular diffusivity (D e ) account for the surrounding milieu. The MR apparent soma radius (r s ) further provides an estimate of (volume-weighted) average cell body size. Mathematically, SANDI models the diffusion signal attenuation as a weighted combination of restricted diffusion in spheres (somas), restricted diffusion in cylinders of infinitesimally small radius (neurites), and Gaussian diffusion in the extra-cellular compartment, constrained by the condition that the compartment fractions sum to one. These parameters thus provide biologically interpretable and mathematically well-defined indices of cellular composition, enabling non-invasive characterization of brain microstructure beyond what can be achieved with conventional diffusion models. SANDI has been shown to provide highly reproducible and repeatable parameter estimation across grey matter regions in the human brain 39 that align closely with its known cyto- and myeloarchitecture. For instance, the gradients of apparent soma density maps were found to closely match those of Brodmann areas in human cortical regions with different soma density profiles, 38 and correlated in the mouse brain with cell density distributions from the Allen atlas. 38 , 40 These findings suggest the potential of the SANDI model for quantifying neurodegenerative processes in the grey matter of the living human brain. Clinical applications of SANDI have so far been limited to multiple sclerosis (MS) and amyotrophic lateral sclerosis (ALS). In MS, reductions in apparent soma and neurite density, together with increases in extracellular signal fraction, have been reported in both grey and white matter, consistent with demyelination, axonal loss, and neurodegeneration. 41 , 42 These abnormalities correlated with disease severity, 42 , 43 cortical and subcortical atrophy, 41 , 44 and elevated serum neurofilament light chain levels, a marker of axonal damage 44 that is also sensitive to HD progression. 45 , 46 Similar patterns, namely reduced apparent soma density and region-specific alterations in apparent soma size have recently been observed in ALS, 43 aligning with motor neuron loss and glial responses. The primary objective of this study was to determine whether SANDI indices were sensitive and specific to microstructural grey matter differences in the BG, compared with the thalami, in individuals with HD relative to healthy controls. The thalami were selected as control regions based on the established trajectory of neurodegeneration in HD, which begins with early loss of MSN in the striatum before extending to neighbouring structures such as the putamen and thalamus. At most participants were at early disease stages, we assumed the thalami would remain relatively unaffected in this sample. Secondary objectives were to explore the extent to which HD-related SANDI differences accounted for BG atrophy, assessed using volumetric measurements, and for performance differences in motor tasks, including speeded and paced finger tapping, which are known to be associated with striatal atrophy in HD. 47 Finally, relationships between SANDI indices and disease burden using the HD-ISS 21 and the CAG-Age Product (CAP 100 ) score 48 were explored. Materials and methods Participants MRI data from 56 individuals with HD and 57 age- and sex-matched healthy controls (HC) were included in the analyses. Imaging data were retrospectively pooled from three projects, all of which acquired scans on the same Siemens Connectom system at the Cardiff University Brain Research Imaging Centre (CUBRIC) using identical acquisition protocols. Thirty-eight of the HD participants were drawn from a randomised controlled feasibility trial of HD-DRUM, 49 a remote rhythmic training intervention, with ethical approval from the Wales Research Ethics Committee 2 (REC Reference: 22/WA/0147). 55 Here we report their baseline MRI and behavioural data collected prior to randomisation. An additional 18 individuals with HD and 18 age-matched HC were included from a previous study of white matter microstructure in premanifest HD (REC Reference: 18/WA/0172). 50 Further comparison data were obtained from 25 age- and sex-matched HC from the Wales Advanced Neuroimaging Database (WAND) Study 51 (REC Reference: 18.08.14.5332RA3) and 14 HC from the HD-DRUM trial. All participants provided written informed consent in accordance with the Declaration of Helsinki. HD individuals were identified and screened for eligibility in five clinics in the UK (Bristol, Birmingham, Cardiff, Exeter, and Liverpool). HC volunteers were recruited from online advertisements on the Cardiff University social network, Viva Engage, or in HD clinics as support partners or family members. HC participants were also recruited through Healthwise Wales and by word of mouth. For the WAND study, data collection was reported elsewhere. 51 Individuals aged 18 years or older with a good command of English were eligible to participate. Additional inclusion criteria for individuals with HD 49 were a positive genetic test for the presence of the mutant huntingtin allele (CAG length ≥ 36 repeats) and/or clinical diagnosis of HD, along with a Unified Huntington’s Disease Rating Scale (UHDRS) Total Functional Capacity (TFC) score between 9 and 13. 52 Exclusion criteria for all participants included an inability to provide informed consent and any contraindication for MRI (e.g. pacemakers, stents). Further exclusion criteria were a history of any other neurological condition for HD participants, and for HC, a history of neurological or psychiatric disorders, and/or alcohol or drug abuse associated with grey matter volume loss. To characterise general cognitive functioning, HD participants completed the Montreal Cognitive Assessment (MOCA). 53 Verbal intellectual ability was assessed using the Test of Premorbid Functioning (TOPF). 54 Disease burden was estimated by the TFC score and the CAG-Age Product (CAP100) score, 48 calculated as: HD-ISS 21 staging and associated clinical information for stratification were obtained from the Enroll-HD 55 observational study (former Registry; REC no 04/WSE05/89). However, due to the retrospective pooling of datasets, HD-ISS information was available for only 30 HD participants (see Supplementary Table 1 for demographic and clinical characteristics). Motor outcome measures Participants in the HD-DRUM study completed a range of motor tasks from the Quantitative-Motor (Q-Motor) test-battery, 56 – 58 which provided reliable assessments of speeded finger tapping performance in clinical HD trials. 59 , 60 Tasks included left and right 1) speeded index finger and foot tapping using force transducers, 47 2) paced finger and foot tapping 47 with a metronome-paced and memory-paced phase, using a fast (0.55s inter-onset interval; IOI) or slow (1.1s IOI) pace, 3) 3D pointing to four target locations in a predefined sequence using a position-tracking stylus with the dominant hand 61 and 4) 3D target pointing and speeded finger tapping dual task performed with dominant and non-dominant hand, respectively. Outcome measures included mean IOI (seconds) and area under the curve (AUC) (Newton-seconds) for the speeded tapping, mean absolute deviation from the metronome pace (seconds) for the paced tapping, target frequency (Hz) for the target pointing, and target frequency, mean IOI and AUC for the dual-task condition. Image acquisition All MRI data included in the analyses were acquired on the same 3T Siemens Connectom scanner (Siemens Healthcare, Erlangen, Germany) with ultra-strong magnetic gradients (300mT/m) at CUBRIC using identical acquisition protocols as detailed below. No scanner changes or upgrades were performed during the data collection periods. T1-weighted (T1w) images were acquired using a magnetisation-prepared 180-degrees radio-frequency pulses and rapid gradient-echo (MPRAGE), with the following parameters: repetition time (TR) 2,300 ms, echo time (TE) 2 ms, field of view (FOV) 256 x 256 x 192 mm, matrix size 256 x 256 x 192, resolution 1 x 1 x 1 mm 3 , flip angle 9°, inversion time (TI) 857 ms, in-plane acceleration (GeneRalised Autocalibrating Partial Parallel Acquisition; GRAPPA) factor 2, phase-encoding direction anterior to posterior (AP), and acquisition time of 6 minutes. Multi-shell High Angular Resolution Diffusion Imaging (HARDI) 62 data were obtained at b-values of 200 s/mm² (20 directions), 500 s/mm² (20 directions), 1,200 s/mm² (30 directions). 2,400 s/mm² (61 directions), 4,000 s/mm² (61 directions) and 6,000 s/mm² (61 directions) using a single-shot spin-echo, echo-planar imaging sequence with TR = 3,000 ms, TE = 59 ms, FOV 220 x 200 mm in-plane; matrix size 110 x 110 x 66; 2 mm 3 resolution, gradient pulse duration - δ = 7 ms, gradient pulses separation - Δ = 24 ms in AP phase-encoding direction with an in-plane acceleration (GRAPPA) factor of 2. Fifteen non-diffusion-weighted (b-value = 0 s/mm²) images were acquired [two initial and 11 interspersed at the 33 rd volume and every 20 th volume thereafter in AP direction and 2 images in the posterior-to-anterior (PA) direction]. The HARDI acquisition time was 18 minutes. Image processing Diffusion-weighted image preprocessing Multi-shell HARDI data were pre-processed and corrected for signal drift, susceptibility-induced distortions, motion and eddy current-induced distortions, gradient non-uniformity and Gibbs ringing artifacts using a custom in-house pipeline comprising tools from the FMRIB Software Library (FSL version 6.0.3), 63 the MRtrix software package, 64 ExploreDTI 65 (version 4.8.6) and in-house MATLAB-based scripts. 39 The FSL brain extraction tool 63 was used to mask the first non-diffusion-weighted image from each phase-encoding direction to exclude non-brain data. The diffusion-weighted MRI volumes were fitted to temporally interspersed b0 volumes to correct for within-image intensity drift by using custom code in MATLAB R2017b (MathWorks Inc., Natick, Massachusetts, USA). Slicewise outlier detection (SOLID) 66 was applied with modified Z-score thresholds of 3.5 (lower) and 10 (upper), utilising a variance-based intensity metric. FSL’s top-up tool 67 , 68 was used to estimate susceptibility-induced off-resonance fields from b0 images that were acquired in opposing phase-encoding directions (AP and PA) and then FSL’s eddy tool 69 was used to correct eddy current-induced distortions and subject movements. Gradient non-uniformity distortions were corrected using in-house code in MATLAB R2017b. Finally, Gibbs ringing correction was performed in MRtrix3 using the local subvoxel-shifts method. 70 For the purpose of comparing our results with the previous literature, 30 DTI was fitted with ExploreDTI using data with b-values of 500 s/mm² and 1,200 s/mm² to produce outcome maps for FA and MD, estimated with linearly weighted least squares regression. SANDI analysis The SANDI model 38 assumes three compartments, namely intra-neurite signal modelled as diffusion inside impermeable randomly oriented sticks, intra-soma signal modelled as restricted diffusion inside spheres, and extra-cellular signal modelled as Gaussian isotropic diffusion. The direction-averaged (or spherical mean) normalized diffusion signal has thus the following expression: where f in + f is + f ec = 1; A stick and A sphere are the normalized, directionally-averaged (or spherical mean) signals for restricted diffusion within neurites and soma, respectively and A ball is the normalized, directionally-averaged (or spherical mean) signal of the extra-cellular space. The specific expressions are given in 38 . The parameters estimated from the direction-averaged (or spherical mean) data are D in , proxy of the intra-neurite effective axial diffusivity; D e , proxy of the extracellular effective mean diffusivity; r s , a proxy of apparent soma radius as well as the signal fractions subject to the constraint f in + f is + f ec = 1, proxy respectively of the relaxation-weighted neurite, soma and extracellular volume fractions. The bulk diffusivity inside the sphere D is is fixed to 3 μm 2 /ms. The parameters were fitted using a Random Forest regression algorithm (TreeBagger Matlab®) with 200 trees, trained on simulated data, using the code publicly available at https://github.com/palombom/SANDI-Matlab-Toolbox-Latest-Release . The training data consisted of simulated signals for 10 5 parameter combinations, uniformly sampled: f in and f is ∈ [0, 1], D in ∈ [0.5, 3] μm 2 /ms, D e ∈ [0.5, 3] μm 2 /ms and r s ∈ [1, 12.5] μm. Rician noise with a distribution of standard deviations randomly sampled from the voxels within the brain mask of the noise map was obtained using Marchenko-Pastur principal component analysis (MP-PCA)-based method 71 – 73 in MRtrix3 and was added to account for realistic SNR levels and rectified noise floor. The loss function of the training was the mean squared error between predicted parameters and ground truth values. These noise maps were subsequently used to fit the SANDI model 38 to the pre-processed multi-shell diffusion data with the SANDI MATLAB Toolbox ( https://github.com/palombom/SANDI-Matlab-Toolbox-Latest-Release 38 using all the default settings. The model fitting produced maps of the intra-neurite, extracellular and intra-soma signal fractions (f in , f ec , f is ), apparent soma size (r s ; measured in μm) and intra- and extra-neurite diffusivities (D in , D e ; measured in mm 2 /ms). Post-hoc sensitivity analysis of the SANDI model parameters revealed very low sensitivity to changes in D in . Consequently, it was excluded from further analysis. T1-weighted image preprocessing The default FreeSurfer 74 (v6) recon-all pipeline was utilised to segment subcortical BG ROIs of the caudate, putamen, and globus pallidum as well as of the thalamus as control ROIs. ROIs were segmented from T1w images and were identified and labelled in each hemisphere. Extraction of microstructural metrics from regions-of-interest Median values of each microstructural index from DTI (FA, MD) and SANDI (f is , f in , f ec , r s , D e ) models were extracted for each ROI using FSL’s fslmaths. ROI masks were aligned with the diffusion space using rigid transformation with FSL’s flirt 75 before eroding the boundaries of the subcortical masks by one voxel with the default 3 × 3 × 3 kernel settings to minimise partial volume effects and then aligning all microstructural maps with the masks. Volumetric measures for each ROI and intra-cranial volume (ICV) were extracted from FreeSurfer v6 ROI volumes were normalised for ICV. The addition of brain volumes allowed exploration of the extent to which any HD-related SANDI differences accounted for BG atrophy. Statistical analysis Statistical analyses were performed in JASP (v0.18.1.0) 76 , R version 4.4.1 (2024-06-14) 77 in R-studio (2024.9.0.375) 78 and SPSS (v27) (IBM Corp). 79 Data normality was assessed using the Shapiro-Wilk test, with p < 0.05 indicating non-normal distribution. Descriptive statistics for each group were reported as percentages (%), means and standard deviations (SD). Medians of each microstructural index in each ROI were compared between the groups with Mann-Whitney-U tests because of lack of normality and unequal variance between groups. Effect sizes (ES) for group comparisons were therefore reported with rank biserial correlation ( r rb ). Multiple comparisons were corrected with Benjamini-Hochberg’s method to control a false discovery rate (FDR) of 0.05 80 and applied to all statistical tests that related to the same theoretical inference. 81 Multi-collinearity of imaging metrics was explored by calculating the intercorrelations between all SANDI, DTI, and volumetric measurements averaged across all four ROIs (caudate, putamen, pallidum, and thalamus). Spearman’s rho ( ρ ) correlation matrices were calculated in the full sample, as well as the HC and HD groups separately. Hierarchical linear regression analyses were conducted to explore SANDI predictors of the variance in volumetric measurements. Regression analyses were carried out for each ROI and each group separately. HD data were modelled by firstly accounting for age and TFC scores (available for all HD participants) simultaneously. This was followed by step-wise inclusion of all SANDI indices using an iterative forward selection and backward elimination method based on each variable’s F-statistic and p-value that aimed to maximise the adjusted R 2 -value while keeping only the most significant predictors. HD and HC data were modelled in the same way using all SANDI variables and age, except for the inclusion of TFC scores that were only available for HD participants. Principal Component Analysis (PCA) was carried out to reduce the dimensionality of HD participants’ Q-Motor data and hence the number of multiple correlations with microstructural SANDI indices. PCA followed established guidelines to limit the number of extracted components in relatively small sample sizes. 82 , 83 First the Kaiser criterion of including all components with an eigenvalue greater than 1 was applied and the Cattell scree plot was inspected to identify the minimal number of components that accounted for most variability in the data. Each extracted component was then assessed for interpretability. PCA was conducted using orthogonal Varimax rotation of the component matrix with Kaiser normalization. Loadings that were greater than 0.5 were considered to be statistically significant. Spearman’s rho ( ρ ) correlations were then calculated between HD participants’ motor component scores and the CAP100 with the SANDI indices, DTI, and volumetric measures in each ROI. HD-ISS categories were obtained either directly from Enroll-HD or calculated using the online HD-ISS calculator ( https://enroll-hd.org/calc/html_basic.htm ) based on caudate and putamen volumes derived from the FreeSurfer v6 74 cross-sectional pipeline. For participants labelled “<2” in the Enroll-HD database (reflecting missing imaging data), we distinguished between stages 0 and 1 using two complementary approaches. First, adjusted normative modelling was applied, in which caudate and putamen volumes from healthy controls were modelled using linear regression with age, sex, and ICV as covariates, and individual z-scores were generated from model residuals. Second, raw group-based normalisation was performed by z-scoring ICV-normalised volumes relative to the control mean and SD. Participants labelled “<2” were classified as HD-ISS 1 if either adjusted caudate or putamen z-scores fell ≥ 2 SD below the control mean; otherwise, they were assigned to stage 0. Both approaches produced identical classifications. Explorative analyses of HD-ISS-related differences in SANDI indices were conducted for BG ROIs. Due to small sample sizes, participants at Stages 0 and 1 were combined (HD-ISS 0-1), as were those at Stage 2 and 3 (HD-ISS 2-3). SANDI metrics for each ROI were averaged across hemispheres. Pairwise comparisons between HD-ISS 0-1 versus HC, HD-ISS 0-1 versus HD-ISS 2-3, and HD-ISS 2-3 versus HC were performed using Mann–Whitney U tests without additional FDR correction. Results Demographics Sample characteristics, including age and sex distribution for the HD and HC groups are described in Table 1 . The groups were comparable with regards to age and sex ( p >0.05). View this table: View inline View popup Download powerpoint Table 1. Demographic and clinical information of participants Information for HD-ISS calculation was available for 30 individuals. Of these 4 were stratified into Stage 0, 9 into Stage 1, 5 into Stage 2, and 12 into Stage 3; no participants were classified as Stage 4. One individual could not be assigned a stage because their clinical presentation did not align with the HD-ISS assumption that motor symptoms precede functional decline. An additional, 25 individuals could not be classified due to missing clinical data or CAG repeat lengths between 36-39. Demographic and clinical characteristics for each HD-ISS stage, as well as for unclassified participants are summarised in Supplementary Table 1. Supplementary Table 2 provides descriptive statistics for the Q-Motor measures. Imaging analysis Differences between HD and HC groups in volumetric measures Table 2 summarises the volumetric group differences across the eight ROIs with their ES. Figure 1 illustrates the ES maps of volumetric differences between HD and HC, as well as raincloud plots for each ROI. Mann-Whitney tests identified significantly reduced volumes in all BG ROIs in HD compared to HC with moderate ES ( r rb = 0.46-0.55). In contrast, no reductions were present in the right ( FDR-p = .05) and left thalamus ( FDR-p = .063) (both r rb = 0.2). Download figure Open in new tab Figure 1. Volumetric differences in the basal ganglia and thalamus between HD and healthy control (HC) groups. Regions-of-interest (ROIs) were segmented using FreeSurfer v6. All ROIs, except the left thalamus, showed significantly smaller volumes in the HD cohort after FDR correction for multiple comparisons. Colours indicate the strength of rank-biserial correlations ( r rb ) from Mann-Whitney U tests: Red = strong effect ( r rb ≥ 0.5), Yellow = medium effect (0.3 ≤ r rb < 0.5), White = small effect ( r rb < 0.3). Raincloud plots show the distribution of the volumetric measures in each ROI per group with orange for HD and green for HC participants. * p < 0.05; *** p < 0.001 View this table: View inline View popup Download powerpoint Table 2. Descriptive and Mann-Whitney Statistics for Intracranial Volume-Normalised Regions of Interest Differences between HD and HC groups in microstructural measures Descriptive and statistical microstructural results are shown in Table 3 . Figure 2 provides ES mappings and raincloud plots of the different microstructural measures for each ROI comparing the two groups. Download figure Open in new tab Figure 2. Microstructural differences in ROIs between HD and healthy control (HC) groups. Median values of each microstructural measure were extracted per ROI. A) HD individuals show reduced soma density (f is ) in the basal ganglia (BG). B) Soma radius (r s ) is elevated in the caudate and putamen but reduced in the pallidum. C) Extracellular water fraction (f ec ) is increased in BG regions in the HD group. D) Extracellular diffusivity (D e ) is higher in the putamen and pallidum. E) Fractional anisotropy (FA) is elevated in the BG, and F) mean diffusivity (MD; expressed in x10 -4 mm²/s) is increased in the striatum. Colours indicate the strength of rank-biserial correlations ( r rb ) from Mann-Whitney U tests: Red = strong effect ( r rb ≥ 0.5), Yellow = medium effect (0.3 ≤ r rb < 0.5), White = small effect ( r rb < 0.3). Raincloud plots show the distribution of the volumetric measures in each ROI per group with orange for HD and green for HC participants. * p < 0.05; ** p < 0.01; *** p < 0.001 View this table: View inline View popup Table 3. Descriptive and Mann-Whitney Statistics for Microstructural Measures in Regions of Interest The SANDI measures revealed between-group differences in apparent soma density, apparent soma size, and extracellular signal fraction with moderate ES in the BG ( r rb = 0.32-0.53) ( Fig. 2 ). HD individuals exhibited reduced apparent soma density ( Fig. 2A ) across all BG ROIs. Apparent soma size was elevated in bilateral caudate and putamen but reduced in the pallidum ( Fig. 2B ). Elevated extracellular signal fraction was observed in all BG regions ( Fig. 2C ). Extracellular diffusivity was higher in the putamen and pallidum but not in the caudate ( Fig. 2D ). No differences were observed for apparent neurite density f in in any of the ROIs ( FDR-p ≥ 0.167, absolute r rb < 0.2). Increased FA and MD ( FDR-p < 0.01) with moderate ES ( r rb = 0.30-0.61) were observed in the BG ( Fig. 2E-F ) with the exception of MD in the left pallidum ( FDR-p = 0.854). No differences were found in the thalami for any of the microstructural metrics ( rrb = 0.02-0.19). Correlations between BG microstructure and motor performance in HD Motor measures PCA extracted one principal component that explained 64% of the Q-Motor data with high loadings (> 0.5 or <-0.5) from all variables (Supplementary Table 3). Spearman’s rho correlation between Q-motor component scores and the disease burden CAP 100 score revealed a positive correlation ( ρ = 0.61, p = 0.002, N = 24), i.e., higher disease burden was associated with higher scores in the Q-Motor component reflecting slower and less accurate motor performance (Supplementary Figure 1). Spearman’s rho correlations between Q-motor component scores and SANDI microstructural measures from all ROIs and scatter plots are displayed in Figure 3 . Figure 4 displays correlations and scatter plots of the Q-Motor component with DTI and volumetric indices. Overall, negative correlations were observed between principal Q-Motor component scores and apparent soma density in the BG, apparent soma size in the pallidum, and volumetric measurements in all ROIs, indicating that lower apparent soma density and volumes were associated with impaired motor performance, reflected in higher scores in the Q-Motor component. Conversely, positive correlations were present between Q-Motor component scores and apparent soma size, extracellular signal fraction, and extracellular diffusivity in the BG, i.e. larger apparent soma size and extracellular signal were associated with impaired motor performance ( Figure 3 ). Download figure Open in new tab Figure 3. A) Correlation matrix and B-G) selected scatter plots illustrating Spearman’s rho correlations between SANDI measures and the Q-Motor principal component. A) Each cell represents the Spearman’s rho correlation strength, with pink indicating negative and green positive correlations. B-G) Each plot includes a best-fit least squares linear regression line with standard error indicated by the grey shaded area, along with the Spearman’s rho ( ρ ) and the corresponding FDR-p value. Regression lines are included for visualisation only and do not reflect variance explained (R²) or imply linear model fit. Scatter dot colours represent participants’ HD-ISS stage. Unclassified refers to those participants who could not be classified due to having CAG 36-40 or incomplete clinical data. Abbreviations: D e : Extracellular diffusivity; f ec : Extracellular signal fraction; f in : Neurite density signal fraction; f is : Soma density signal fraction; PC : Principal component; r s : Soma radius; TFC : Total functional capacity. Download figure Open in new tab Figure 4. A) Correlation matrix and B-D) selected scatter plots illustrating Spearman’s rho correlations between DTI, volumetric measures and the Q-Motor principal component. A) Each cell represents the Spearman’s rho correlation strength, with pink indicating negative and green positive correlations. B-D) Each plot includes a best-fit least squares linear regression line with standard error indicated by the grey shaded area, along with the Spearman’s rho ( ρ ) and the corresponding FDR-p value. Regression lines are included for visualisation only and do not reflect variance explained (R²) or imply linear model fit. Scatter dot colours represent participants’ HD-ISS stage. Unclassified refers to those participants who could not be classified due to having CAG 36-40 or incomplete clinical data. Abbreviations: FA : Fractional anisotropy; MD : Mean diffusivity. Microstructural predictors of BG volumes Hierarchical linear regression analyses of the HC data ( Table 4 , Figure 5 ) revealed that 17% of volumetric variation in the left caudate was accounted for by age and apparent soma density while 11% of variation in the right caudate was explained by age and extracellular diffusivity (D e ). Age accounted for 25% of volumetric differences in the left putamen, and together with D e for 29% of differences in right putamen. Age alone explained 17% of left and 11% of right thalamic volume variation. No significant regression models were present for bilateral globus pallidus. Download figure Open in new tab Figure 5. Standardised beta coefficients of SANDI microstructural metrics predicting volume (normalised for intracranial volume) in regions of interest in the healthy control group. Data were modelled by firstly accounting for age, followed by the step-wise inclusion of all SANDI indices. The figure displays the predictor variables included in the final regression models for each region of interest. Abbreviations: De : Extracellular diffusivity; fis : Soma density signal fraction; rs : Soma radius; TFC: Total functional capacity. * p < 0.05; ** p < 0.01; *** p < 0.001. View this table: View inline View popup Table 4. Hierarchical Linear Regression Predicting Normalised Volumes from SANDI Microstructural Metrics, Controlling for Age in the Healthy Control Participants In contrast, regression analyses for the HD data ( Table 5 , Figure 6 ) showed that 60% of volume variation in the left caudate and 51% of variation in the right caudate were accounted for by age and apparent soma density and size. Similarly, 57% of variation in the right putamen volume were explained by age and apparent soma size, while 63% of volume differences in the left putamen were explained by age, apparent soma size as well as extracellular signal and diffusivity. Comparable to the HC results, age alone accounted for 30% of volume variation in the left and for 35% in the right thalamus while no age effects were present for bilateral globus pallidus. However, in HD individuals 42% of volume variation in the left pallidum was explained by apparent soma size and extracellular signal and 27% of volume variation in the right pallidum by extracellular signal fraction only. No significant contributions of TFC were present. Download figure Open in new tab Figure 6. Standardised beta coefficients of SANDI microstructural metrics predicting volume (normalised for intracranial volume) in A) left and B) right hemisphere regions of interest in individuals with Huntington’s disease. Data were modelled by firstly accounting for age and Total Functional Capacity scores simultaneously, followed by the step-wise inclusion of all SANDI indices. The figure displays the predictor variables included in the final regression models for each region of interest. Abbreviations: De : Extracellular diffusivity; fec : Extracellular signal fraction; fis : Soma density signal fraction; rs : Soma radius. * p < 0.05; ** p < 0.01; *** p < 0.001. View this table: View inline View popup Table 5. Hierarchical Linear Regression Models Predicting Normalised Volumes in each Region of Interest from SANDI Microstructural Metrics, Controlling for Age and TFC in HD Participants Correlation of disease burden with microstructural and volumetric measures The CAP 100 score as an index of disease progression was correlated with brain measurements to explore the relationship between disease burden and microstructural and volumetric differences ( Figure 7 ). Figure 7A summarises correlation coefficient strengths and levels of significance. Figures 7B-J display scatter plots of significant correlations between the CAP 100 and microstructural and volumetric measures. Increased CAP 100 was negatively correlated with apparent soma density, and volume size in bilateral caudate and putamen, as well as with pallidal apparent soma size. Positive correlations were observed between CAP 100 and extracellular diffusivity and signal fraction, apparent soma size in caudate, putamen, and thalamus, apparent neurite density in putamen, FA in the BG, and MD in caudate and putamen. Download figure Open in new tab Figure 7. A) Correlation matrix and B-J) selected scatter plots illustrating Spearman’s rho correlations between SANDI, DTI, and volumetric measures with CAP 100 . Each scatter plot includes a best-fit least squares linear regression line with standard error indicated by the grey shaded area, along with the Spearman’s rho ( ρ ) and the corresponding FDR-p value. Regression lines are included for visualisation only and do not reflect variance explained (R²) or imply linear model fit. Scatter dot colours represent participants’ HD-ISS stage and those who were not classified due to having CAG<40 or incomplete clinical data. Abbreviations: De : Extracellular diffusivity; FA : Fractional anisotropy; fec : Extracellular signal fraction; fin : Neurite density signal fraction; fis : Soma density signal fraction; MD : Mean diffusivity; rs : Soma radius; vol : Normalised volume. Exploratory pairwise comparisons between HD-ISS groups and HC Exploratory comparisons of apparent soma density, apparent soma size, extracellular signal fraction and extracellular diffusivity, averaged across left and right BG ROIs, were conducted between HD-ISS 0-1 (n = 13), HD-ISS 2-3 (n = 17), and HC (Supplementary Table 1). Compared with HC, HD-ISS 0-1 showed increased extracellular signal fraction in the BG, larger apparent soma size in the caudate and putamen, and lower apparent soma density alongside larger extracellular diffusivity in the pallidum ( r rb = -0.255–0.483, p ≤ 0.033). Relative to HD-ISS 0-1, those at HD-ISS 2-3 exhibited reduced apparent soma density and increased extracellular diffusivity in the caudate and putamen, as well as smaller apparent soma size in the pallidum ( r rb = -0.439–0.447, p = 0.01–0.031). Compared with HC, the HD-ISS 2-3 group showed differences in the same directions across all metrics in all BG regions. Descriptive statistics and full statistical analysis results are provided in Supplementary Table 4, and effect sizes with 95% confidence intervals are shown in Supplementary Figure 2. Cross-correlations between microstructural and volumetric measures Cross-correlation matrices between SANDI microstructural indices, DTI metrics, and normalised volumes are shown in Supplementary Figure 3 for the full sample, and separately for HC and HD groups. Across all participants, moderate to strong correlations were observed between conceptually related diffusion-derived measures (e.g. mean diffusivity with extracellular diffusivity, fractional anisotropy with neurite-related indices), indicating partially shared variance as expected. In the HC group, SANDI metrics reflecting apparent soma properties (soma radius and soma fraction) showed distinct correlation patterns: soma and neurite fraction were strongly inversely correlated ( ρ = -0.92) but neither correlated with soma radius. All SANDI measures correlated weakly with FA and volume and moderately with MD. Group-specific matrices revealed broadly similar correlation structures in HC and HD participants. However, several associations, particularly between SANDI indices and normalised volume, were stronger in the HD group, consistent with disease-related coupling between microstructural alterations and macroscopic atrophy. Discussion This study is the first to investigate HD-related microstructural differences in the BG using SANDI, a diffusion MRI technique designed to probe grey matter microstructure. The objective of the study was to explore SANDI indices as potential non-invasive in vivo markers of HD neuropathology that may offer greater specificity than volumetric measurements. Apparent soma density and apparent soma size were sensitive to HD-related differences in the BG and together with age, accounted for up to 63% of striatal atrophy. SANDI indices also correlated with motor impairments and CAP 100 disease burden. Together, these findings highlight the potential of SANDI metrics as imaging biomarkers of disease progression and as surrogate outcome measures for assessing the neural effects of emerging disease-modifying therapies in HD and other neurodegenerative diseases. Microstructural and Volumetric Alterations in HD Well-established patterns of marked volume loss accompanied by increases in FA and MD in the BG were replicated in HD compared with healthy controls. FA increases in the caudate and putamen are thought to reflect the selective degeneration of MSN connections. 20 , 84 No microstructural differences and only trends for volumetric reduction were observed in the thalami. This pattern of macro- and microstructural differences in the BG aligns with previous reports in premanifest and early manifest HD stages. 31 , 85 Importantly, SANDI revealed novel HD-related microstructural differences in the BG. Compared with healthy controls, those with HD showed reduced apparent soma density together with increased extracellular signal fraction and diffusivity across the caudate, putamen, and pallidum, but not the thalami. Apparent soma size was also increased in the caudate and putamen, whereas it was reduced in the globus pallidus. HD-related reductions in apparent soma density in the BG are consistent with the loss of striatal MSN, the histopathological hallmark of HD, 22 , 23 and with downstream degeneration of pallidal neurons by the loss of striatal projections. Changes in apparent soma size may reflect alterations in neural and glial cell proportions and/or morphology, including astrocyte and microglia swelling in response to neurodegeneration 22 , 24 – 26 as well as soma shrinkage preceding neuronal cell death. 86 , 87 Thus, increased apparent soma size in the striatum may indicate HD-related reorganisation of cell composition driven by MSN loss and reactive glial cell swelling, whereas smaller apparent soma size in the pallidum may reflect infiltration of smaller glia cells associated with secondary neuronal loss following striatal degeneration. Exploratory analyses of BG SANDI differences across HD-ISS Stages 0-3 (Supplementary Figure 2) showed that reductions in striatal apparent soma density and increases in extracellular diffusivity were evident in HD-ISS 2-3 but not in HD-ISS 0-1, whereas increases in striatal apparent soma size and extracellular signal fraction were present in both groups. In the pallidum, reductions in apparent soma density and increases in extracellular signal fraction and diffusivity were observed in both HD-ISS groups, while reduced apparent soma size was only present in HD-ISS 2-3. These findings suggest that SANDI indices are sensitive to changes in BG cellular composition both at early (HD-ISS 0-1) and across disease stages (HD-ISS 0-3). However, these preliminary observations require validation in larger, prospective cohorts before SANDI indices can be considered reliable biomarker of disease progression. With the development of disease-modifying therapies for HD, the need for non-invasive imaging markers that can sensitively capture treatment effects has become increasingly important. Striatal volume remains the gold-standard imaging marker, and DTI has been widely used to assess tissue microstructure in both cross-sectional and longitudinal studies. 2, 13 – 15 , 88 An important question therefore concerns how SANDI indices compare with these established imaging markers and the extent to which they share variance. In this study we observed moderate effect sizes for group differences in striatal SANDI indices, comparable to those for volume, FA, and MD. ES for apparent soma size overlapped with those for volume and MD, while apparent soma density showed slightly lower but broadly similar ES to FA. Multicollinearity between microstructural and volumetric imaging measures was assessed using cross-correlations between all imaging metrics (SANDI, DTI, volume) averaged across ROIs (see Supplementary Figure 2). In HCs, apparent soma and neurite fractions were strongly inversely correlated but neither correlated with soma radius. SANDI indices showed only weak correlations with FA and volume, but moderate correlations with MD. This pattern suggests that soma density and size parameters capture microstructural features that are distinct from neurite fraction, FA, and volume, but overlap to some extent with MD, an unspecific metric influenced by multiple biological and geometrical tissue properties. 89 , 90 In HD participants, correlations between SANDI metrics and volumes were stronger, and soma radius and fraction were inversely correlated, likely reflecting disease-related reorganisation of microstructural properties associated with atrophy. Microstructural predictors of basal ganglia atrophy in HD To explore the relationship between SANDI microstructural and volumetric measures in the BG, hierarchical linear regression analyses were conducted for each ROI and for HD and HC groups separately, testing age and SANDI indices as predictors of volume. For the HD group, TFC was included as a marker of disease burden. These analyses showed that SANDI indices accounted for a substantial proportion of BG atrophy but not thalamic volume. Up to 63% of HD-related striatal atrophy was explained by apparent soma density, apparent soma size, and age. In the pallidum, apparent soma size and extracellular signal fraction accounted for 27% of atrophy on the right and 42% on the left, whereas age alone predicted thalamic volume. This pattern in the thalami mirrored findings in healthy controls, where age was the strongest predictor across all ROIs (except bilateral pallidum) and the only significant predictor in the thalami; diffusivity contributed to right-lateralised BG regions and apparent soma density to the left caudate. Together, these results demonstrate that SANDI-derived measures of apparent soma density and size capture HD-related differences in striatal grey-matter microstructure in vivo . As outlined above, reductions in apparent soma density and increases in apparent soma size are consistent with the characteristic loss of medium spiny neurons and reactive gliosis in HD and explain a significant proportion of atrophy in the caudate and putamen. Associations Between SANDI microstructural indices and Motor Function The observed correlations between SANDI metrics and motor measures provide novel insights into the functional implications of microstructural alterations in HD. Consistent with the role of the BG in motor initiation and coordination, 91 the present study showed that reduced apparent soma density in the striatum and reduced apparent soma size in the pallidum were associated with poorer Q-Motor performance, which itself was linked to greater disease burden. This was reflected in increased IOI and ACU during speeded tapping tasks, and in impaired performance on paced tapping tasks. Similarly, increased apparent soma size in the striatum, together with elevated extracellular signal and diffusivity across all BG regions, as well as striatal FA and MD, were associated with motor impairments. These findings indicate that microstructural changes related to BG neurodegeneration and glial reactivity are directly linked to subtle motor deficits in HD. Limitations HD-ISS classification was possible for only 30 of the 56 HD participants (54%) due to the retrospective pooling of MRI datasets collected before HD-ISS information was incorporated into the Enroll-HD database, missing clinical data, and HD-ISS model exclusions for individuals with CAG repeat lengths between 36-39. Of the 30 classified participants, 4 were assigned to i Stage 0, 9 to Stage 1, 5 ito Stage 2, and 12 to Stage 3; none were classified as Stage 4. Because of the small sample sizes, HD-ISS comparisons were conducted on combined groups (HD-ISS 0-1 versus HD-ISS 2-3). These combined groups showed some overlap in imaging and behavioural features ( Figures 3 , 4, and 7) and therefore do not represent discrete disease stages. Consequently, these exploratory findings should be interpreted with caution and require replication in larger, prospective cohorts before SANDI metrics can be considered markers of disease progression. In addition, several methodological considerations regarding HD-ISS classification should be noted. Classification was obtained either from Enroll-HD or calculated using the online HD-ISS tool. For participants labelled “<2” in Enroll-HD, Stages 0 and 1 were assigned based on z-scores for striatal volumes derived from the FreeSurfer v6 cross-sectional pipeline. However, HD-ISS classification and online calculator were developed using volume cut-offs generated from the longitudinal FreeSurfer v6 pipeline. Consequently, applying cross-sectional data in this study may have introduced classification bias. Furthermore, hierarchical regression models were fitted separately within HD and HC groups and were not intended to support formal statistical comparisons of regression coefficients across groups. Rather, these analyses were designed to explore within-group associations and to provide a descriptive comparison of patterns of relationships linking microstructural indices to regional atrophy. Hierarchical regression models may also be susceptible to over-fitting, and therefore the reported coefficients may not reflect true out-of-sample R 2 . The assessment of potential multicollinearity between SANDI-derived metrics was based primarily on the healthy control data, as correlations observed in the HD group likely reflect disease-related shifts in microstructural relationships rather than intrinsic model dependencies. In healthy controls, soma radius did not correlate with apparent soma or neurite fractions, whereas a pronounced inverse relationship was observed between soma and neurite fractions, consistent with the biophysical constraints of the model. Importantly, correlations between SANDI metrics and conventional DTI measures (FA) and regional volume were generally weak, indicating that apparent soma estimates capture distinct aspects of grey matter microstructure rather than redundant information. Moderate correlations with mean diffusivity were observed and are expected, given that MD is a non-specific metric sensitive to multiple microstructural features and susceptible to partial volume effects associated with tissue loss. Consequently, individual regression coefficients should not be interpreted in isolation as fully independent biological effects. As with all model-based approaches applied to correlated imaging parameters, some degree of shared variance and parameter coupling is unavoidable. As we cannot rule out some degeneracy between apparent soma density f is and apparent soma size r s , 𝑟 𝑠 should not be interpreted as absolute histological soma radii, and inferences should instead focus on relative group differences and spatial patterns. Similarly, HD-related microstructural changes (e.g., changes to membrane permeability) could affect model parameter fidelity, and thus 𝑓 𝑖𝑠 and 𝑟 𝑠 should be treated as MRI-derived effective indices rather than direct quantitative measures of neuron loss or glial hypertrophy. However, the correlation structure observed in healthy controls suggests that multicollinearity is unlikely to critically undermine model stability or interpretability in the present analyses. Importantly, all comparisons were made under an identical acquisition protocol and fitting pipeline, meaning that group-level contrasts remain informative even if absolute parameter values are biased. The primary conclusions of the present study therefore do not rely on the precise magnitude of individual coefficients but on the consistent pattern and direction of associations across regions and measures. Together, these factors support the feasibility and initial clinical relevance of SANDI for probing striatal microstructural alterations in Huntington’s disease. Clinical Implications and Future Directions The present study acquired multi-shell (max b-value = 6,000 s/mm 2 ) DWI data on a non-clinical ultra-strong gradient (300mT/m) 3T MRI system. Ultra-strong gradients improve the signal-to-noise-ratio at high b-values by enabling shorter echo times (TE), thereby enhancing sensitivity to small water displacement, 39 and reducing bias from inter-compartmental exchange. 92,93 However, for clinical translation it is important to demonstrate that HD-related SANDI differences can be detected on standard clinical MRI systems without requiring ultra-strong gradients. Although equivalent data in HD are not yet available, we have shown the feasibility of estimating SANDI metrics from multi-shell DWI acquired on a 3T MRI system with standard gradient strength (67mT/m) (maximum b-value of 6,000 s/mm 2 ) in healthy adults 94 and people with MS 42 , 44 ,94 . Additionally, Zeng et al., 43 reported significant differences in SANDI metrics between individuals with ALS and healthy controls using a 3T MRI Prisma system with a maximum b-value of 3,000 s/mm 2 , and two further studies using 3T clinical scanners with standard gradient strength showed detectable microstructural SANDI alterations in MS. Collectively, these findings indicate that SANDI can be implemented on clinical scanners, particularly as commercial systems progress toward stronger gradient capabilities such as Siemens Magnetom Cima.X. Our findings highlight the potential of SANDI-derived metrics as future markers for tracking disease progression and assessing therapeutic efficacy in HD, as well as in more prevalent neurodegenerative conditions such as Alzheimer’s and Parkinson’s disease. The sensitivity of apparent soma density and size, and extracellular water signal to microstructural changes in HD provides information complementary to volumetric measures, which are currently the most widely used imaging modality in clinical trials. Furthermore, the observed associations between SANDI metrics, motor performance, and disease burden underscore their relevance for evaluating the neural effects of emerging disease-modifying treatments. Ultimately, the choice of imaging markers depends on scientific or clinical objective: DTI and volumetric measures are informative for documenting differences or tracking longitudinal change, whereas gaining mechanistic insight into tissue microstructure and therapeutic effects requires going beyond DTI and volumes. Advanced models such as SANDI offer the potential to provide more specific and biologically meaningful characterisation of grey matter pathology. Conclusion Our study demonstrates the potential of SANDI for characterising microstructural abnormalities in individuals with HD. These abnormalities align with ex vivo neuropathological findings of striatal medium spiny neuron loss and gliosis, account for a substantial proportion of striatal atrophy, and correlate with motor impairment and disease burden. Although conventional MRI lacks the resolution to directly capture histopathology, advanced biophysical models such as SANDI may help bridge this gap by providing biologically interpretable parameters that reflect tissue composition and capture histopathological changes in vivo . This approach offers a promising avenue for advancing HD research and improving clinical care. Data availability This research utilised baseline data from the HD-DRUM project that has been endorsed by the Enroll-HD Scientific Oversight Committee (SOC) (14/11/2022). At the end of the HD-DRUM project, the coded study data will be shared and made accessible to the research community via the Enroll-HD specific data request process. WAND data are publicly available ( https://git.cardiff.ac.uk/cubric/wand ). Data Availability This research utilised baseline data from the HD–DRUM project that has been endorsed by the Enroll–HD Scientific Oversight Committee (SOC) (14/11/2022). At the end of the HD–DRUM project, the coded study data will be shared and made accessible to the research community via the Enroll–HD specific data request process. WAND data are publicly available (https://git.cardiff.ac.uk/cubric/wand). Funding This work was supported by a National Institute for Health Research (NIHR) and Health and Care Research Wales (HCRW) Advanced Fellowship to CM-B (grant number: NIHR-FS(A)- 2022). The Centre for Trials Research at Cardiff University receives infrastructure funding from HCRW. MP is supported by the UKRI Future Leaders Fellowship MR/T020296/2. CC was funded by a Wellcome Trust PhD studentship (204005/Z/16/Z) and LL by a PhD studentship of the School of Psychology at Cardiff University. The WAND project was funded by a Wellcome Trust Investigator Award (096646/Z/11/Z), a Wellcome Trust Strategic Award (104943/Z/14/Z), and a Wellcome Discovery Awards (227882/Z/23/Z and 317797/Z/24/Z). Competing interests The authors have no conflicting interests to declare. Supplementary figure legends Download figure Open in new tab Supplementary Figure 1. Scatterplot showing positive relationship between the Q-Motor principal component and the disease burden measure (CAP 100 ) with the Spearman’s rho ( ρ ) test. Scatter dot colours represent participants’ HD-ISS stage and those who were not classified due to having CAG<40 or incomplete clinical data. Download figure Open in new tab Supplementary Figure 2. Bar plot showing the effect sizes and 95% confidence intervals for exploratory pairwise comparisons between HD-ISS 0-1 (n = 13), HD-ISS 2-3 (n = 17), and healthy controls (HC). Significant comparisons are marked with * (p < 0.05), ** (p < 0.01), and *** (p < 0.001). Download figure Open in new tab Supplementary Figure 3. Correlation heatmaps showing the cross-correlation of SANDI, DTI, and volumetric (normalised for intracranial volume) measures, averaged across the caudate, putamen, pallidum, and thalamus. Heatmaps are shown separately for the full sample, healthy controls (HC), and Huntington’s disease (HD) participants. Correlations are expressed as Spearman’s rho coefficients. View this table: View inline View popup Download powerpoint Supplementary Table 1. Demographic and clinical information per HD-ISS stage View this table: View inline View popup Download powerpoint Supplementary Table 2. Descriptive statistics for motor outcome measures View this table: View inline View popup Download powerpoint Supplementary Table 3. Rotated Component Loadings on the Q-motor Outcome Measures View this table: View inline View popup Download powerpoint Supplementary Table 4. Descriptive statistics and non-parametric pairwise comparisons for SANDI indices (significant for total HD sample vs HC) between HD-ISS 0-1, HD-ISS 2-3, and Healthy Controls. Acknowledgements We would like to thank Amy Dangerfield, Allison Cooper and Sonya Foley-Bozorgzad for their help with data collection as well as Derek Jones and John Evans for their advice and support with regards to the implementation of MRI data acquisition protocols. We would like to thank the following clinical and administrative staff at the participating patient identification centres for their help with identifying suitable patients for the study: Eileen Donovan, Kim Munnery, and Jane Davies from the Cardiff HD clinic; Jessica Prado Mota from the Royal Devon University Healthcare NHS Foundation Trust in Exeter; Jenni Burns from the Walton Centre NHS Foundation Trust in Liverpool; Natalie Rosewell, Anya Soonderpershad, and Dr Liz Coulthard from the Bristol Brain Centre; Claire Tilley and Dr Hugh Rickards from the Birmingham and Solihull Mental Health NHS Foundation Trust. In addition, we would like to thank all Public Involvement contributors and the members of the Enroll-HD Scientific Oversight Committee for their input into the study as well as all participants for their generous time commitment to help us conducting this research. Footnotes This version has been revised in accordance with peer review feedback as detailed on https://doi.org/10.7554/eLife.107661.1 References 1. ↵ Rüb U , Seidel K , Heinsen H , Vonsattel JP , den Dunnen WF , Korf HW . Huntington’s disease (HD): the neuropathology of a multisystem neurodegenerative disorder of the human brain . Brain Pathology . 2016 ; 26 ( 6 ): 726 – 740 . doi: 10.1111/bpa.12426 OpenUrl CrossRef PubMed 2. ↵ Tabrizi SJ , Scahill RI , Durr A , et al. Biological and clinical changes in premanifest and early stage Huntington’s disease in the TRACK-HD study: the 12-month longitudinal analysis . Lancet Neurol . 2011 ; 10 ( 1 ): 31 – 42 . doi: 10.1016/S1474-4422(10)70276-3 OpenUrl CrossRef PubMed Web of Science 3. Roos RAC . Huntington’s disease: A clinical review . Orphanet J Rare Dis . 2010 ; 5 ( 1 ). doi: 10.1186/1750-1172-5-40 OpenUrl CrossRef PubMed 4. ↵ Walker FO . Huntington’s disease . The Lancet . 2007 ; 369 ( 9557 ): 218 – 228 . doi: 10.1016/S0140-6736(07)60111-1 OpenUrl CrossRef PubMed Web of Science 5. ↵ Ross CA , Aylward EH , Wild EJ , et al. Huntington disease: Natural history, biomarkers and prospects for therapeutics . Nat Rev Neurol.Nature Publishing Group . 2014 ; 10 ( 4 ): 204 – 216 . doi: 10.1038/nrneurol.2014.24 OpenUrl CrossRef PubMed 6. ↵ Aylward EH , Liu D , Nopoulos PC , et al. Striatal volume contributes to the prediction of onset of Huntington disease in incident cases . Biol Psychiatry . 2012 ; 71 ( 9 ): 822 – 828 . doi: 10.1016/j.biopsych.2011.07.030 OpenUrl CrossRef PubMed Web of Science 7. ↵ Tan B , Shishegar R , Poudel GR , Fornito A , Georgiou-Karistianis N . Cortical morphometry and neural dysfunction in Huntington’s disease: a review . Eur J Neurol . 2021 ; 28 ( 4 ): 1406 – 1419 . doi: 10.1111/ene.14648 OpenUrl CrossRef PubMed 8. Kinnunen KM , Schwarz AJ , Turner EC , et al. Volumetric MRI-Based Biomarkers in Huntington’s Disease: An Evidentiary Review . Front Neurol . 2021 ; 12 ( September ). doi: 10.3389/fneur.2021.712555 OpenUrl CrossRef 9. Wilkes FA , Abaryan Z , Ching CRK , et al. Striatal morphology and neurocognitive dysfunction in Huntington disease: The IMAGE-HD study . Psychiatry Res Neuroimaging . 2019 ; 291 ( July ): 1 – 8 . doi: 10.1016/j.pscychresns.2019.07.003 OpenUrl CrossRef 10. ↵ Scahill RI , Zeun BMBS P , Osborne-Crowley K , et al. Biological and Clinical Characteristics of Gene Carriers Far from Predicted Onset in the Huntington’s Disease Young Adult Study (HD-YAS): A Cross-Sectional Analysis . Vol 19 .; 2020 . 11. ↵ Paulsen JS , Langbehn DR , Stout JC , et al. Detection of Huntington’s disease decades before diagnosis: The Predict-HD study . J Neurol Neurosurg Psychiatry . 2008 ; 79 ( 8 ): 874 – 880 . doi: 10.1136/jnnp.2007.128728 OpenUrl Abstract / FREE Full Text 12. Biglan KM , Zhang Y , Long JD , et al. Refining the diagnosis of huntington disease: The PREDICT-HD study . Front Aging Neurosci . 2013 ; 5 ( APR ): 1 – 8 . doi: 10.3389/fnagi.2013.00012 OpenUrl CrossRef PubMed 13. ↵ Tabrizi SJ , Langbehn DR , Leavitt BR , et al. Biological and clinical manifestations of Huntington’s disease in the longitudinal TRACK-HD study: cross-sectional analysis of baseline data . Lancet Neurol . 2009 ; 8 ( 9 ): 791 – 801 . doi: 10.1016/S1474-4422(09)70170-X OpenUrl CrossRef PubMed Web of Science 14. Tabrizi SJ , Reilmann R , Roos RA , et al. Potential endpoints for clinical trials in premanifest and early Huntington’s disease in the TRACK-HD study: analysis of 24 month observational data . Lancet Neurol . 2012 ; 11 ( 1 ): 42 – 53 . doi: 10.1016/S1474-4422(11)70263-0 OpenUrl CrossRef PubMed Web of Science 15. ↵ Tabrizi SJ , Scahill RI , Owen G , et al. Predictors of phenotypic progression and disease onset in premanifest and early-stage Huntington’s disease in the TRACK-HD study: Analysis of 36-month observational data . Lancet Neurol . 2013 ; 12 ( 7 ): 637 – 649 . doi: 10.1016/S1474-4422(13)70088-7 OpenUrl CrossRef PubMed Web of Science 16. ↵ Tabrizi SJ , Estevez-Fraga C , van Roon-Mom WMC , et al. Potential disease-modifying therapies for Huntington’s disease: lessons learned and future opportunities . Lancet Neurol . 2022 ; 21 ( 7 ): 645 – 658 . doi: 10.1016/S1474-4422(22)00121-1 OpenUrl CrossRef PubMed 17. ↵ Vonsattel PJG , Keller C , Cortes Ramirez PE . Huntington’s disease-neuropathology . In: Handbook of Clinical Neurology . Vol 100. 2011 : 83 – 100 . 18. Liu CF , Younes L , Tong XJ , et al. Longitudinal imaging highlights preferential basal ganglia circuit atrophy in Huntington’s disease . Brain Commun . 2023 ; 5 ( 5 ). doi: 10.1093/braincomms/fcad214 OpenUrl CrossRef 19. Johnson EB , Ziegler G , Penny W , et al. Dynamics of Cortical Degeneration Over a Decade in Huntington’s Disease . Biol Psychiatry . 2021 ; 89 ( 8 ): 807 – 816 . doi: 10.1016/j.biopsych.2020.11.009 OpenUrl CrossRef PubMed 20. ↵ Estevez-Fraga C , Scahill R , Rees G , Tabrizi SJ , Gregory S . Diffusion imaging in Huntington’s disease: Comprehensive review . J Neurol Neurosurg Psychiatry . 2021 ; 92 ( 1 ): 62 – 69 . doi: 10.1136/jnnp-2020-324377 OpenUrl Abstract / FREE Full Text 21. ↵ Tabrizi SJ , Schobel S , Gantman EC , et al. A biological classification of Huntington’s disease: the Integrated Staging System . Lancet Neurol . 2022 ; 21 ( 7 ): 632 – 644 . doi: 10.1016/S1474-4422(22)00120-X OpenUrl CrossRef PubMed 22. ↵ Myers RH , Vonsattel JP , Paskevich PA , et al. Decreased Neuronal and Increased Oligodendroglial Densities in Huntington’s Disease Caudate Nucleus . J Neuropathol Exp Neurol . 1991 ; 50 ( 6 ): 729 – 742 . doi: 10.1097/00005072-199111000-00005 OpenUrl CrossRef PubMed 23. ↵ Vis JC , Schipper E , de Boer-van H , et al. Expression pattern of apoptosis-related markers in Huntington’s disease . Acta Neuropathol . 2005 ; 109 ( 3 ): 321 – 328 . doi: 10.1007/s00401-004-0957-5 OpenUrl CrossRef PubMed Web of Science 24. ↵ Vis JC , Nicholson LFB , Faull RLM , Evans WH , Severs NJ , Green CR . Connexin expression in Huntington’s diseased human brain . Cell Biol Int . 1998 ; 22 ( 11-12 ): 837 – 847 . doi: 10.1006/cbir.1998.0388 OpenUrl CrossRef PubMed Web of Science 25. Simmons DA , Casale M , Alcon B , Pham N , Narayan N , Lynch G . Ferritin accumulation in dystrophic microglia is an early event in the development of Huntington’s Disease . Glia . 2007 ; 55 ( 10 ): 1074 – 1084 . doi: 10.1002/glia.20526 OpenUrl CrossRef PubMed Web of Science 26. ↵ Sapp E , Kegel KB , Aronin N , et al. Early and Progressive Accumulation of Reactive Microglia in the Huntington Disease Brain . J Neuropathol Exp Neurol . 2001 ; 60 ( 2 ): 161 – 172 . doi: 10.1093/jnen/60.2.161 OpenUrl CrossRef PubMed 27. ↵ Le Bihan D , Turner R , Moonen CTW , Pekar J . Imaging of diffusion and microcirculation with gradient sensitization: Design, strategy, and significance . Journal of Magnetic Resonance Imaging . 1991 ; 1 ( 1 ): 7 – 28 . doi: 10.1002/jmri.1880010103 OpenUrl CrossRef PubMed Web of Science 28. ↵ Basser PJ , Mattiello J , LeBihan D . MR diffusion tensor spectroscopy and imaging . Biophys J . 1994 ; 66 ( 1 ): 259 – 267 . doi: 10.1016/S0006-3495(94)80775-1 OpenUrl CrossRef PubMed Web of Science 29. ↵ Basser PJ , Pierpaoli C . Microstructural and Physiological Features of Tissues Elucidated by Quantitative-Diffusion-Tensor MRI . J Magn Reson B . 1996 ; 111 ( 3 ): 209 – 219 . doi: 10.1006/jmrb.1996.0086 OpenUrl CrossRef PubMed Web of Science 30. ↵ Georgiou-Karistianis N , Gray MA , Domínguez D JF , et al. Automated differentiation of pre-diagnosis Huntington’s disease from healthy control individuals based on quadratic discriminant analysis of the basal ganglia: The IMAGE-HD study . Neurobiol Dis . 2013 ; 51 : 82 – 92 . doi: 10.1016/j.nbd.2012.10.001 OpenUrl CrossRef PubMed 31. ↵ Hobbs NZ , Cole JH , Farmer RE , et al. Evaluation of multi-modal, multi-site neuroimaging measures in Huntington’s disease: Baseline results from the PADDINGTON study . Neuroimage Clin . 2013 ; 2 ( 1 ): 204 – 211 . doi: 10.1016/j.nicl.2012.12.001 OpenUrl CrossRef 32. ↵ Jones DK , Alexander DC , Bowtell R , et al. Microstructural imaging of the human brain with a ‘super-scanner’: 10 key advantages of ultra-strong gradients for diffusion MRI . Neuroimage . 2018 ; 182 ( May ): 8 – 38 . doi: 10.1016/j.neuroimage.2018.05.047 OpenUrl CrossRef PubMed 33. ↵ Assaf Y , Basser PJ . Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain . Neuroimage . 2005 ; 27 ( 1 ): 48 – 58 . doi: 10.1016/j.neuroimage.2005.03.042 OpenUrl CrossRef PubMed 34. ↵ Zhang H , Schneider T , Wheeler-Kingshott CA , Alexander DC . NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain . Neuroimage . 2012 ; 61 ( 4 ): 1000 – 1016 . doi: 10.1016/j.neuroimage.2012.03.072 OpenUrl CrossRef PubMed 35. ↵ Martinez-Heras E , Grussu F , Prados F , Solana E , Llufriu S . Diffusion-Weighted Imaging: Recent Advances and Applications. Seminars in Ultrasound , CT and MRI . 2021 ; 42 ( 5 ): 490 – 506 . doi: 10.1053/j.sult.2021.07.006 OpenUrl CrossRef 36. ↵ Daducci A , Dal Palu A , Lemkaddem A , Thiran JP . COMMIT: Convex Optimization Modeling for Microstructure Informed Tractography . IEEE Trans Med Imaging . 2015 ; 34 ( 1 ): 246 – 257 . doi: 10.1109/TMI.2014.2352414 OpenUrl CrossRef PubMed 37. ↵ Zhang J , Gregory S , Scahill RI , et al. In vivo characterization of white matter pathology in premanifest huntington’s disease . Ann Neurol . 2018 ; 84 ( 4 ): 497 – 504 . doi: 10.1002/ana.25309 OpenUrl CrossRef PubMed 38. ↵ Palombo M , Ianus A , Guerreri M , et al. SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI . Neuroimage . 2020 ; 215 ( August 2019 ): 116835 . doi: 10.1016/j.neuroimage.2020.116835 OpenUrl CrossRef PubMed 39. ↵ Koller K , Rudrapatna U , Chamberland M , et al. MICRA: Microstructural image compilation with repeated acquisitions . Neuroimage . 2021 ; 225 ( October 2020 ): 117406 . doi: 10.1016/j.neuroimage.2020.117406 OpenUrl CrossRef PubMed 40. ↵ Ianuş A , Carvalho J , Fernandes FF , et al. Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas . Neuroimage . 2022 ; 254 . doi: 10.1016/j.neuroimage.2022.119135 OpenUrl CrossRef PubMed 41. ↵ Krijnen EA , Russo AW , Salim Karam E , et al. Detection of grey matter microstructural substrates of neurodegeneration in multiple sclerosis . Brain Commun . 2023 ; 5 ( 3 ). doi: 10.1093/braincomms/fcad153 OpenUrl CrossRef 42. ↵ Margoni M , Pagani E , Preziosa P , et al. In vivo quantification of brain soma and neurite density abnormalities in multiple sclerosis . J Neurol . 2023 ; 270 ( 1 ): 433 – 445 . doi: 10.1007/s00415-022-11386-3 OpenUrl CrossRef PubMed 43. ↵ Zeng JY , Huang HW , Zhuang SP , et al. Soma and neurite density imaging detects brain microstructural impairments in amyotrophic lateral sclerosis . Eur J Radiol . 2025 ; 184 . doi: 10.1016/j.ejrad.2025.111981 OpenUrl CrossRef 44. ↵ Barakovic M , Weigel M , Cagol A , et al. A novel imaging marker of cortical “cellularity” in multiple sclerosis patients . Sci Rep . 2024 ; 14 ( 1 ): 9848 . doi: 10.1038/s41598-024-60497-6 OpenUrl CrossRef PubMed 45. ↵ Byrne LM , Rodrigues FB , Johnson EB , et al. Evaluation of Mutant Huntingtin and Neurofilament Proteins as Potential Markers in Huntington’s Disease . Vol 10.; 2018 . https://www.science.org 46. ↵ Sampedro F , Pérez-Pérez J , Martínez-Horta S , et al. Cortical microstructural correlates of plasma neurofilament light chain in Huntington’s disease . Parkinsonism Relat Disord . 2021 ; 85 : 91 – 94 . doi: 10.1016/j.parkreldis.2021.03.008 OpenUrl CrossRef PubMed 47. ↵ Bechtel N , Scahill RI , Rosas HD , et al. Tapping linked to function and structure in premanifest and symptomatic Huntington disease . Neurology . 2010 ; 75 ( 24 ): 2150 – 2160 . doi: 10.1212/WNL.0b013e3182020123 OpenUrl CrossRef PubMed 48. ↵ Warner JH , Long JD , Mills JA , et al. Standardizing the CAP Score in Huntington’s Disease by Predicting Age-at-Onset . J Huntingtons Dis . 2022 ; 11 ( 2 ): 153 – 171 . doi: 10.3233/JHD-210475 OpenUrl CrossRef PubMed 49. ↵ Ioakeimidis V , Busse M , Drew CJG , et al. Protocol for a randomised controlled unblinded feasibility trial of HD-DRUM: a rhythmic movement training application for cognitive and motor symptoms in people with Huntington’s disease . BMJ Open . 2024 ; 14 ( 7 ): e082161 . doi: 10.1136/bmjopen-2023-082161 OpenUrl Abstract / FREE Full Text 50. ↵ Casella C , Chamberland M , Laguna PL , et al. Mutation-related magnetization-transfer, not axon density, drives white matter differences in premanifest Huntington disease: Evidence from in vivo ultra-strong gradient MRI . Hum Brain Mapp . 2022 ;(October 2021): 3439 – 3460 . doi: 10.1002/hbm.25859 OpenUrl CrossRef PubMed 51. ↵ McNabb CB , Driver ID , Hyde V , et al. WAND: A multi-modal dataset integrating advanced MRI, MEG, and TMS for multi-scale brain analysis . Sci Data . 2025 ; 12 ( 1 ): 220 . doi: 10.1038/s41597-024-04154-7 OpenUrl CrossRef PubMed 52. ↵ Siesling S , Van Vugt JPP , Zwinderman KAH , Kieburtz K , Roos RAC . Unified Huntington’s disease rating scale: A follow up . Movement Disorders . 1998 ; 13 ( 6 ): 915 – 919 . doi: 10.1002/mds.870130609 OpenUrl CrossRef PubMed Web of Science 53. ↵ Nasreddine ZS , Phillips NA , Bédirian V , et al. The Montreal Cognitive Assessment, MoCA: A Brief Screening Tool For Mild Cognitive Impairment . J Am Geriatr Soc . 2005 ; 53 ( 4 ): 695 – 699 . doi: 10.1111/j.1532-5415.2005.53221.x OpenUrl CrossRef PubMed Web of Science 54. ↵ Wechsler D . Test of Premorbid Functioning . UK version (TOPF UK) . Preprint posted online 2011 . 55. ↵ Schubert R , Fiedler H , Barallon P , Habbel B , Reilmann R . F55 Q-motor speeded tapping and pointing assessments: feasibility of a dual task design. In: F: Clinical Studies: Case Reports , Observational Studies and Trials. BMJ Publishing Group Ltd ; 2022 : A56 . 1 -A56. doi: 10.1136/jnnp-2022-ehdn.146 OpenUrl Abstract / FREE Full Text 56. ↵ Schubert R , Rosser A , Bachoud-Lévi AC , Craufurd D , Reilmann R . F63 Preliminary results from Q-MOTOR/Q-COG analyses in the repair-hd study . In: Clinical Studies. BMJ Publishing Group Ltd ; 2018 : A62.1 – A62 . doi: 10.1136/jnnp-2018-EHDN.164 OpenUrl CrossRef 57. Reilmann R , Schubert R . Motor Outcome Measures in Huntington Disease Clinical Trials . Vol 144 . 1st ed. Elsevier B.V. ; 2017 . doi: 10.1016/B978-0-12-801893-4.00018-3 OpenUrl CrossRef PubMed 58. ↵ Reilmann R , Rosser A , Kostyk S , et al. F41 The proof-hd phase 3 study: pridopidine’s outcome on function in huntington disease (PROOF). In: F: Clinical Studies: Case Reports , Observational Studies and Trials. BMJ Publishing Group Ltd ; 2021 : A36.1 – A36 . doi: 10.1136/jnnp-2021-EHDN.84 OpenUrl CrossRef 59. ↵ Reilmann R , Anderson KE , Feigin A , et al. Safety and efficacy of laquinimod for Huntington’s disease (LEGATO-HD): a multicentre, randomised, double-blind, placebo-controlled, phase 2 study . Lancet Neurol . 2024 ; 23 ( 3 ): 243 – 255 . doi: 10.1016/S1474-4422(23)00454-4 OpenUrl CrossRef PubMed 60. ↵ Hermle D , Schubert R , Barallon P , et al. Multifeature quantitative motor assessment of upper limb ataxia including drawing and reaching . Ann Clin Transl Neurol . 2024 ; 11 ( 5 ): 1097 – 1109 . doi: 10.1002/acn3.52024 OpenUrl CrossRef PubMed 61. ↵ Tuch DS , Reese TG , Wiegell MR , Makris N , Belliveau JW , Wedeen VJ . High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity . Magn Reson Med . 2002 ; 48 ( 4 ): 577 – 582 . doi: 10.1002/mrm.10268 OpenUrl CrossRef PubMed Web of Science 62. ↵ Smith SM . Fast robust automated brain extraction . Hum Brain Mapp . 2002 ; 17 ( 3 ): 143 – 155 . doi: 10.1002/hbm.10062 OpenUrl CrossRef PubMed Web of Science 63. ↵ Tournier JD , Smith R , Raffelt D , et al. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation . Neuroimage . 2019 ; 202 : 116137 . doi: 10.1016/j.neuroimage.2019.116137 OpenUrl CrossRef PubMed 64. ↵ Leemans A , Jeurissen B , Sijbers J , Jones DK . ExploreDTI: a graphical toolbox for processing, analyzing, and visualizing diffusion MR data . In: ISMRM . 2009 . 65. ↵ Sairanen V , Leemans A , Tax CMW . Fast and accurate Slicewise OutLIer Detection (SOLID) with informed model estimation for diffusion MRI data . Neuroimage . 2018 ; 181 : 331 – 346 . doi: 10.1016/j.neuroimage.2018.07.003 OpenUrl CrossRef PubMed 66. ↵ Andersson JLR , Skare S , Ashburner J . How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging . Neuroimage . 2003 ; 20 ( 2 ): 870 – 888 . doi: 10.1016/S1053-8119(03)00336-7 OpenUrl CrossRef PubMed Web of Science 67. ↵ Smith SM , Jenkinson M , Woolrich MW , et al. Advances in functional and structural MR image analysis and implementation as FSL . Neuroimage . 2004 ; 23 : S208 – S219 . doi: 10.1016/j.neuroimage.2004.07.051 OpenUrl CrossRef PubMed Web of Science 68. ↵ Andersson JLR , Sotiropoulos SN . An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging . Neuroimage . 2016 ; 125 : 1063 – 1078 . doi: 10.1016/j.neuroimage.2015.10.019 OpenUrl CrossRef PubMed 69. ↵ Kellner E , Dhital B , Kiselev VG , Reisert M . Gibbs-ringing artifact removal based on local subvoxel-shifts . Magn Reson Med . 2016 ; 76 ( 5 ): 1574 – 1581 . doi: 10.1002/mrm.26054 OpenUrl CrossRef PubMed 70. ↵ Cordero-Grande L , Christiaens D , Hutter J , Price AN , Hajnal J V . Complex diffusion-weighted image estimation via matrix recovery under general noise models . Neuroimage . 2019 ; 200 : 391 – 404 . doi: 10.1016/j.neuroimage.2019.06.039 OpenUrl CrossRef PubMed 71. ↵ Veraart J , Fieremans E , Novikov DS . Diffusion MRI noise mapping using random matrix theory . Magn Reson Med . 2016 ; 76 ( 5 ): 1582 – 1593 . doi: 10.1002/mrm.26059 OpenUrl CrossRef PubMed 72. Veraart J , Novikov DS , Christiaens D , Ades-aron B , Sijbers J , Fieremans E . Denoising of diffusion MRI using random matrix theory . Neuroimage . 2016 ; 142 : 394 – 406 . doi: 10.1016/j.neuroimage.2016.08.016 OpenUrl CrossRef PubMed 73. ↵ Fischl B. FreeSurfer . Neuroimage . 2012 ; 62 ( 2 ): 774 – 781 . doi: 10.1016/j.neuroimage.2012.01.021 OpenUrl CrossRef PubMed Web of Science 74. ↵ Jenkinson M . Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images . Neuroimage . 2002 ; 17 ( 2 ): 825 – 841 . doi: 10.1016/S1053-8119(02)91132-8 OpenUrl CrossRef PubMed Web of Science 75. ↵ JASP Team . JASP . Preprint posted online 2023 . 76. ↵ R Core Team . R: A Language and Environment for Statistical Computing . R Foundation for Statistical Computing . Preprint posted online 2024. 77. ↵ RStudio Team . RStudio: Integrated Development for R . RStudio, PBC . Preprint posted online 2024. 78. ↵ IBM Corp . IBM SPSS Statistics for Windows . IBM Corp . Preprint posted online 2024 . 79. ↵ Benjamini Y , Hochberg Y . Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing . J R Stat Soc Series B Stat Methodol . 1995 ; 57 ( 1 ): 289 – 300 . doi: 10.1111/j.2517-6161.1995.tb02031.x OpenUrl CrossRef PubMed 80. ↵ Lakens D . Performing high-powered studies efficiently with sequential analyses . Eur J Soc Psychol . 2014 ; 44 ( 7 ): 701 – 710 . doi: 10.1002/ejsp.2023 OpenUrl CrossRef PubMed 81. ↵ Preacher KJ , MacCallum RC . Exploratory Factor Analysis in Behavior Genetics Research: Factor Recovery with Small Sample Sizes . Behav Genet . 2002 ; 32 ( 2 ): 153 – 161 . doi: 10.1023/A:1015210025234 OpenUrl CrossRef PubMed 82. ↵ de Winter JCF , Dodou D , Wieringa PA . Exploratory Factor Analysis With Small Sample Sizes . Multivariate Behav Res . 2009 ; 44 ( 2 ): 147 – 181 . doi: 10.1080/00273170902794206 OpenUrl CrossRef PubMed Web of Science 83. ↵ Matz OC , Spocter M . The Effect of Huntington’s Disease on the Basal Nuclei: A Review . Cureus. Published online April 25 , 2022 . doi: 10.7759/cureus.24473 OpenUrl CrossRef 84. ↵ Sánchez-Castañeda C , Cherubini A , Elifani F , et al. Seeking huntington disease biomarkers by multimodal, cross-sectional basal ganglia imaging . Hum Brain Mapp . 2013 ; 34 ( 7 ): 1625 – 1635 . doi: 10.1002/hbm.22019 OpenUrl CrossRef PubMed Web of Science 85. ↵ Fricker M , Tolkovsky AM , Borutaite V , Coleman M , Brown GC . Neuronal Cell Death . Physiol Rev . 2018 ; 98 ( 2 ): 813 – 880 . doi: 10.1152/physrev.00011.2017 OpenUrl CrossRef PubMed 86. ↵ Kerr JFR , Wyllie AH , Currie AR . Apoptosis: A Basic Biological Phenomenon with Wideranging Implications in Tissue Kinetics . Br J Cancer . 1972 ; 26 ( 4 ): 239 – 257 . doi: 10.1038/bjc.1972.33 OpenUrl CrossRef PubMed Web of Science 87. ↵ Klaus A , Alves da Silva J , Costa RM. What, If, and When to Move: Basal Ganglia Circuits and Self-Paced Action Initiation . Annu Rev Neurosci . 2019 ; 42 ( 1 ): 459 – 483 . doi: 10.1146/annurev-neuro-072116-031033 OpenUrl CrossRef PubMed 88. ↵ Olesen JL , Østergaard L , Shemesh N , Jespersen SN . Diffusion time dependence, power-law scaling, and exchange in gray matter . Neuroimage . 2022 ; 251 : 118976 . doi: 10.1016/j.neuroimage.2022.118976 OpenUrl CrossRef PubMed 89. ↵ Jelescu IO , de Skowronski A , Geffroy F , Palombo M , Novikov DS . Neurite Exchange Imaging ((NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange . Neuroimage . 2022 ; 256 . doi: 10.1016/j.neuroimage.2022.119277 OpenUrl CrossRef PubMed 90. ↵ Schiavi S , Palombo M , Zacà D , et al. Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics . Hum Brain Mapp . 2023 ; 44 ( 13 ): 4792 – 4811 . doi: 10.1002/hbm.26416 OpenUrl CrossRef PubMed View the discussion thread. Back to top Previous Next Posted February 04, 2026. Download PDF Data/Code Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging Vasileios Ioakeimidis , Marco Palombo , Chiara Casella , Lucy Layland , Carolyn B. McNabb , Robin Schubert , Philip Pallmann , Monica E. Busse , Cheney J. G. Drew , Sundus Alusi , Timothy Harrower , Jane Davies , Anne E. Rosser , Claudia Metzler-Baddeley medRxiv 2025.03.17.25324107; doi: https://doi.org/10.1101/2025.03.17.25324107 Share This Article: Copy Citation Tools In vivo mapping of striatal neurodegeneration in Huntington’s disease with Soma and Neurite Density Imaging Vasileios Ioakeimidis , Marco Palombo , Chiara Casella , Lucy Layland , Carolyn B. McNabb , Robin Schubert , Philip Pallmann , Monica E. Busse , Cheney J. G. Drew , Sundus Alusi , Timothy Harrower , Jane Davies , Anne E. Rosser , Claudia Metzler-Baddeley medRxiv 2025.03.17.25324107; doi: https://doi.org/10.1101/2025.03.17.25324107 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Neurology Subject Areas All Articles Addiction Medicine (568) Allergy and Immunology (863) Anesthesia (300) Cardiovascular Medicine (4434) Dentistry and Oral Medicine (444) Dermatology (382) Emergency Medicine (608) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1509) Epidemiology (15227) Forensic Medicine (30) Gastroenterology (1124) Genetic and Genomic Medicine (6595) Geriatric Medicine (668) Health Economics (997) Health Informatics (4534) Health Policy (1368) Health Systems and Quality Improvement (1613) Hematology (540) HIV/AIDS (1264) Infectious Diseases (except HIV/AIDS) (15916) Intensive Care and Critical Care Medicine (1103) Medical Education (623) Medical Ethics (146) Nephrology (667) Neurology (6599) Nursing (346) Nutrition (998) Obstetrics and Gynecology (1144) Occupational and Environmental Health (957) Oncology (3332) Ophthalmology (974) Orthopedics (369) Otolaryngology (420) Pain Medicine (436) Palliative Medicine (130) Pathology (663) Pediatrics (1692) Pharmacology and Therapeutics (691) Primary Care Research (711) Psychiatry and Clinical Psychology (5447) Public and Global Health (9230) Radiology and Imaging (2198) Rehabilitation Medicine and Physical Therapy (1370) Respiratory Medicine (1196) Rheumatology (593) Sexual and Reproductive Health (712) Sports Medicine (530) Surgery (712) Toxicology (99) Transplantation (289) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'a0007841185f4193',t:'MTc3OTUwMTY1NQ=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();