Natural biological variation of white matter
microstructure is accentuated in Huntington's disease

Sarah Gregory, Helen Crawford, Kiran Seunarine, Blair Leavitt, Alexandra Durr, Raymund Ac Roos, Rachael I. Scahill, Sarah J Tabrizi, Geraint Rees, Douglas Langbehn, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Sarah Gregory, Helen Crawford, Kiran Seunarine, Blair Leavitt, Alexandra Durr, et al. Natural biological variation of white matter microstructure is accentuated in Huntington's disease. Human Brain Mapping, 2018, 39 (9), pp.3516-3527. 10.1002/hbm.24191. hal-01867947

HAL Id: hal-01867947

https://hal.sorbonne-universite.fr/hal-01867947

Submitted on 4 Sep 2018

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# Natural biological variation of white matter microstructure is accentuated in Huntington's disease

Sarah Gregory<fnref n="1" /> ID

Helen Crawford<fnref n="1" />

Kiran Seunarine<fnref n="2" />

Blair Leavitt<fnref n="3" />

Alexandra Durr<fnref n="4" />

Raymund A. C. Roos<fnref n="5" />

Rachael I. Scahill<fnref n="1" /> ID

Sarah J. Tabrizi<fnref n="1" />

Geraint Rees<fnref n="6" />

Douglas Langbehn<fnref n="7" />

Michael Orth<fnref n="8" /> ID

**Correspondence**

Sarah Gregory, Huntington's Disease Research Centre, UCL Institute of Neurology, 2nd Floor, Russell Square House, 10-12 Russell Square, London WC1B 5EH.

Email: s.gregory@ucl.ac.uk

**Funding information**

This work was funded by the CHDI Foundation, the Wellcome Trust (GR), and the Medical Research Council (ST). Some of this work was also undertaken at UCLH/UCL who acknowledge support from the Department of Health's NIHR Biomedical Research Centre. SG, RS, GR, and ST receive support from a Wellcome Trust Collaborative Award (200181/Z/15/Z)

**Abstract**

Huntington's disease (HD) is a monogenic neurodegenerative disorder caused by a CAG-repeat expansion in the Huntingtin gene. Presence of this expansion signifies certainty of disease onset, but only partly explains age at which onset occurs. Genome-wide association studies have shown that naturally occurring genetic variability influences HD pathogenesis and disease onset. Investigating the influence of biological traits in the normal population, such as variability in white matter properties, on HD pathogenesis could provide a complementary approach to understanding disease modification. We have previously shown that while white matter diffusivity patterns in the left sensorimotor network were similar in controls and HD gene-carriers, they were more extreme in the HD group. We hypothesized that the influence of natural variation in diffusivity on effects of HD pathogenesis on white matter is not limited to the sensorimotor network but extends to cognitive, limbic, and visual networks. Using tractography, we investigated 32 bilateral pathways within HD-related networks, including motor, cognitive, and limbic, and examined diffusivity metrics using principal components analysis. We identified three independent patterns of diffusivity common to controls and HD gene-carriers that predicted HD status. The first pattern involved almost all tracts, the second was limited to sensorimotor tracts, and the third encompassed cognitive network tracts. Each diffusivity pattern was associated with network specific performance. The consistency in diffusivity patterns across both groups coupled with their association with disease status and task performance indicates that naturally-occurring patterns of diffusivity can become accentuated in the presence of the HD gene mutation to influence clinical brain function.

<footnote n="1">¹Huntington's Disease Research Centre, UCL Institute of Neurology, London, United Kingdom</footnote>

<footnote n="2">²Developmental Imaging and Biophysics Section, UCL Institute of Child Health, London, WC1N 1EH, United Kingdom</footnote>

<footnote n="3">³Centre for Molecular Medicine and Therapeutics, Department of Medical Genetics, University of British Columbia, Vancouver, BC V5Z 4H4, Canada</footnote>

<footnote n="4">⁴APHP Department of Genetics, Groupe Hospitalier Pitié-Salpêtrière, and Institut du Cerveau et de la Moelle, INSERM U1127, CNRS UMR7225, Sorbonne Universités – UPMC Université Paris VI UMR_S1127, Paris, France</footnote>

<footnote n="5">⁵Department of Neurology, Leiden University Medical Centre, Leiden, 2300RC, The Netherlands</footnote>

<footnote n="6">⁶Wellcome Trust Centre for Neuroimaging, University College London, London, United Kingdom</footnote>

<footnote n="7">⁷Departments of Psychiatry and Biostatistics, University of Iowa, Iowa City, Iowa</footnote>

<footnote n="8">⁸Department of Neurology, Ulm University Hospital, Ulm, Germany</footnote>

<footnote>*TrackOn-HD is funded by the CHDI foundation, a not for profit organization dedicated to finding treatments for Huntington's disease. There are no conflicts of interest regarding the research in the manuscript.</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2018 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.</footnote>
### KEYWORDS

Diffusion Tensor Imaging, Huntington's disease, natural variability

## 1 | INTRODUCTION

Huntington's disease (HD) is caused by a CAG-repeat-expansion in the Huntingtin gene (HTT). The expansion length is variable and explains much of the heterogeneity in the age at which mutation-carriers receive a clinical diagnosis of manifest disease (Duyao et al., 1993). A large portion of variability in age at onset, however, remains unexplained by the gene-mutation and is likely the result of other factors, both genetic and environmental (Wexler et al., 2004). Identifying these factors and their associated biological processes could provide additional targets for therapeutic interventions that may modify disease progression.

Approximately, half of the unexplained variability for age at onset may be heritable. It is therefore unsurprising that much focus has centered on the role of potential genetic modifiers that do not cause disease but nonetheless influence the effects of a disease-causing mutated gene (GeM-HD Consortium, 2015; Lee et al., 2012). Genome-wide association analysis identified possible independent modifiers that may influence the age at clinical onset in HD (GeM-HD Consortium, 2015). However, further work is required to understand the mechanisms of modifier genes in a comprehensive biological context. As a complementary approach, it would be valuable to first identify biological traits in the normal population that may influence HD pathogenesis and then investigate the underlying e.g., genetic basis. Similar to genomic variability the expression of biological traits, e.g., height or eye color, varies without causing disease. The expression of that biological trait could, however, influence the vulnerability to the effects of a mutated gene. TRACK-HD and PREDICT-HD, for example, reported changes in brain structure that were independent of *HTT* CAG repeat length and age (Paulsen et al., 2014; Tabrizi et al., 2013). The identification of biological traits could help reveal an interaction between the underlying biology and pathogenesis, or the timing of clinical manifestations of HD. Such traits, and their genetic bases, could propose a route to disease modification.

Macrostructural brain changes are a consistent phenotype in HD. In addition, considerable evidence points towards changes in white matter microstructure in HD cohorts. Diffusivity is an indirect marker of white matter tract organization. Increases of diffusivity across both the whole brain and in select white matter pathways, such as within the sensorimotor network, are suggestive of white matter degeneration in premanifest (preHD) and early manifest HD (Della Nave et al., 2010; Douaud et al., 2009; Dumas et al., 2012; Klöppel et al., 2008; Matsui et al., 2014; Novak et al., 2014; Odish et al., 2015; Poudel et al., 2014, 2015). In a recent study of the sensorimotor network, we observed a broad structural HD phenotype that encompassed gray and white matter volume, cortical thickness, and altered diffusivity in left white matter tracts in the sensorimotor network, which was associated with CAG repeat length (Orth et al., 2016). In addition, we identified an inverse relationship between levels of axial diffusivity (AD; water movement in

the direction of the main tract) and radial diffusivity (RD; water movement perpendicular to the main tract) that predicted both motor performance and disease stage in HD and was independent of CAG repeat length. Interestingly, this same diffusivity pattern was also evident in controls, suggesting that it may reflect natural biological variation in white matter microstructure, which in the presence of the HD gene-mutation independently contributes to HD pathogenesis.

Given the extant evidence of widespread white matter changes, we have investigated the potential role of natural biological variation on white matter microstructural alterations beyond those that we had described within the left sensorimotor network in HD. We tested to what extent the effects of natural variability on HD pathogenesis are specific only to sensorimotor white matter tracts most likely affected by HD pathology or whether there is a more widespread effect across the whole brain. We performed tractography analyses to investigate diffusivity in 32 tracts within networks likely to be affected as part of the HD phenotype: sensorimotor, cognitive, limbic (neuropsychiatric), and visual and subsequently employed principal components analysis (PCA) on diffusivity metrics extracted from these tracts to identify and parsimoniously describe diffusivity patterns in both controls and HD gene-mutation carriers. First, we examined whether patterns differed between controls, preHD and earlyHD or if there was some consistency within the control and HD populations. Second, we aimed to distinguish diffusivity patterns linked to key contributors to HD pathology, e.g., CAG repeat length, from those that are independent but still influence HD manifestation. Third, we investigated the relative contribution of the three diffusivity measures i.e., RD, AD, and FA to each pattern of diffusivity or principal component. Finally, we examined the relationship between each PC and both cognitive and motor performance as indexed by the global cognitive composite and grip force, respectively. We predicted that tracts associated with networks more directly affected in HD, in particular those connecting striatal regions would likely show more disruption than cortico-cortical tracts and that the main patterns of diffusivity would discriminate robustly between controls, preHD and earlyHD gene-mutation carriers. We also predicted that patterns of natural variation would be associated with HD-related biological variables, including volume changes and CAG-repeat length extending to structural connections in other brain networks and therefore representing a wider structural HD phenotype.

## 2 | METHODS AND MATERIALS

### 2.1 | Participants

We analyzed brain images from the TrackOn-HD study. Participants were recruited into the TrackOn-HD study at four study sites (London, Paris, Leiden, Vancouver) as previously described (Klöppel et al., 2015; Orth et al., 2016). For the present analyses, we used data of those participants who had complete Diffusion Tensor Imaging (DTI) data. Sixty-
one individuals (F = 44%; mean age ± SD: 43.1 ± 9.1) carried the *HTT* gene with a CAG trinucleotide repeat expansion mutation of ≥ 39 but did not have a clinical diagnosis of HD. Thirteen participants (F = 62%, mean age ± SD: 43.1 ± 3.3) had early motor manifest HD (UHDRS diagnostic confidence level of 3 or 4), and 79 were controls (F = 61%, mean age ± SD: 49.1 ± 9.7). A global cognitive composite score was derived from nine cognitive tasks that were completed in testing sessions separate to the MRI procedures: Stroop Word Reading test, Symbol Digit Modality Test, Paced Tapping, Circle Tracing (two conditions), Map Search test, Cancelation task, the Spot the Change visual working memory task, Mental Rotation task (Jones et al., 2014). The United Huntington's Disease Rating Scale (UHDRS) motor examination was administered to all participants to derive the UHDRS total motor score. Grip force as a sensorimotor test was selected as a further marker of motor performance; for details see (Kloppel et al., 2015; Orth et al., 2016). All participants were right-handed. Exclusion criteria included age below 18 or above 65 (unless previously in Track-HD study), major psychiatric, neurological or medical disorder or a history of severe head injury (Kloppel et al., 2015). The study was approved by the local ethics committees, and all participants gave written informed consent according to the Declaration of Helsinki.

## 2.2 | MRI data acquisition and analysis

Standardisation of data acquisition across sites was performed based on previous suggestions (Glover et al., 2012; Kloppel et al., 2015). 3T-MRI data were acquired on two different scanner systems (Philips Achieva at Leiden and Vancouver and Siemens TIM Trio at London and Paris). Diffusion-weighted images were collected with 42 unique gradient directions ($b = 1,000 \text{ s/mm}^2$) with eight images with no diffusion weighting ($b = 0 \text{ s/mm}^2$) (Siemens) and one image with no diffusion weighting ($b = 0 \text{ s/mm}^2$) (Philips). Acquisition parameters were TE = 88 ms, TR = 13 s and voxel size $2 \times 2 \times 2$ mm (Siemens); TE = 56ms and TR = 11s and voxel size $1.96 \times 1.96 \times 2.75$ mm (Philips). The diffusion data were preprocessed using standard FSL pipelines (Smith et al., 2004).

The T1 scan was segmented into gray and white matter using the VBM8 toolbox (http://dbm.neuro.uni-jena.de/vbm/) and combined to create an improved anatomical scan for DTI data registration. Prior to analysis, we screened each DTI dataset for artifacts, signal drop-out and motion. Data were then corrected for distortions caused by eddy currents and motion using eddyCorrect in FSL, and vector gradient information updated accordingly. The no-gradient (B0) image (an averaged B0 image was used for the Siemens data) was then skull-stripped using the Brain Extraction Tool (BET) and manually corrected. For registration purposes, we also applied BET to the structural T1 image. To improve the quality of the brain mask, we combined and dilated a thresholded segmented image with an eroded brain-extracted T1 mask, which was then applied to the original brain-extracted T1 image. We then linearly registered the resultant T1 image to the B0 image using FLIRT (Jenkinson & Smith, 2001) with standard parameters. Diffusion tensors were fit to the corrected data using dtifit and fractional anisotropy (FA) axial diffusivity (AD) and radial diffusivity (RD) values derived.

<table><caption>TABLE 1 Tracts included in principal component analysis</caption><thead><tr><th>Network</th><th>Tract</th></tr></thead><tbody><tr><td rowspan="5">Sensorimotor network</td><td>Primary Motor Cortex (M1) – Putamen</td></tr><tr><td>Primary Motor Cortex (M1) – Motor Thalamus</td></tr><tr><td>Premotor Cortex – Putamen</td></tr><tr><td>Premotor Cortex – Motor Thalamus</td></tr><tr><td>Primary Somatosensory Cortex (S1) – Putamen</td></tr><tr><td rowspan="2">Visual network</td><td>Primary Somatosensory Cortex (S1) – Somatosensory thalamus</td></tr><tr><td>Inferior Fronto-Occipital Fasciculus</td></tr><tr><td rowspan="3">Limbic network</td><td>Primary Visual Cortex (V1) – Visual Thalamus</td></tr><tr><td>Uncinate Fasciculus</td></tr><tr><td>Mid Cingulum</td></tr><tr><td rowspan="3">Cognitive network</td><td>Posterior Cingulum</td></tr><tr><td>Dorsolateral Prefrontal Cortex – Caudate</td></tr><tr><td>Dorsolateral Prefrontal Cortex – Prefrontal Thalamus</td></tr><tr><td rowspan="2">"Non-HD"</td><td>Posterior Parietal Cortex – Parietal Thalamus</td></tr><tr><td>Inferior Lateral Fasciculus</td></tr><tr><td rowspan="2">Interhemispheric tracts</td><td>Mid corpus callosum</td></tr><tr><td>Posterior corpus callosum</td></tr></tbody></table>

Tractography was performed for all listed tracts. Diffusivity measures were extracted from each tract and included in our PCA analysis.

We modeled within-voxel crossing fibers using a Bayesian probabilistic method implemented in Bedpostx (Behrens et al., 2003).

Probabilistic tractography was performed for a series of tracts using probtrackx (Behrens, Berg, Jbabdi, Rushworth, & Woolrich, 2007); more information can be found at Table 1. For the sensorimotor network, this included tracts in the right hemisphere connecting the primary motor cortex (M1) and the motor thalamus; the premotor cortex (PMC) and the motor thalamus; and the primary somatosensory cortex (S1) and the somatosensory thalamus (for completeness we also included the equivalent left-sided tracts from our previous study in the final analyses). This was then repeated replacing the thalamic region with the putamen in both hemispheres. Tracking was then performed on bilateral tracts within the cognitive network connecting the dorsolateral prefrontal cortex (DLPFC) and the DLPFC thalamic parcellation and caudate, respectively; and the posterior parietal cortex (PPC) and the parietal thalamus. We also performed tractography on tracts within “limbic circuits,” the middle and posterior cingulum and the uncinate fasciculus (UF) bilaterally. Finally, we performed fiber-tracking between the primary visual cortex (V1) and the visual thalamus, the inferior fronto-occipital fasciculus (IFOF), the inferior longitudinal fasciculus (ILF) and the middle and posterior corpus callosum (see Supporting Information Table S1 for further details). Masks were created using the Anatomy Toolbox, the WFU PickAtlas (Maldjian, Laurienti, Kraft, & Burdette, 2003) or FSL MNI template and the JHU White Matter Labels
<figure><img src="image_3.png" /><figcaption>FIGURE 1 Study workflow. After processing data from the DTI sequence, tractography was employed to discern 15 tracts bilaterally, the mid corpus callosum and posterior corpus callosum so that, altogether, we derived 32 tracts. From each tract we extracted information about axial (AD) and radial diffusivity (RD) as well as fractional anisotropy (FA). Hence, we included 96 diffusivity measures in our PCA. The first 3 principle components, each of which contains information on all 3 diffusivity measures (AD, RD, FA) explained 56% of the variability in the data. The heatmap illustrates the relative contribution of each tract's diffusivity measures to the respective PC to show the diffusivity patterns for each PC across all tracts. For the three main PCs, we then calculated a principle component score for each participant followed by group comparisons of PCs [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

atlas (Wakana et al., 2007). All masks were defined in standard MNI space and warped into native space for each participant. Exclusion masks were used to exclude any streamlines that may track via the contralateral hemisphere or outside of the anatomically-defined tract. A white matter termination mask was also used to ensure tracts stopped at the gray/white matter interface and did not extend beyond the white matter into the gray matter, CSF or dura. The resulting tracts were then warped into diffusion space using FLIRT. All tracts were visually inspected following probtrackx analyses and warping into

diffusion space. Due to poor tracking in a high number of participants, the following tracts were excluded from the final analyses: the PPC and the parietal thalamus, the lateral orbitofrontal cortex to the caudate, the anterior cingulum and the anterior corpus callosum. Mean fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) values were extracted from the mask of each tract for each participant weighted by the contribution of each voxel by the number of streamlines that pass through it (as opposed to assuming an equal contribution from each voxel).
<figure><img src="image_4.png" /><figcaption>FIGURE 2 Principal component analysis. Scree plot showing variance explained by the first 10 principal components</figcaption></figure>

## 2.3 | Statistical analysis

We performed PCA using centered, standardized (standard deviation = 1) versions of the tract-specific DTI measures. PCA is a method to identify patterns of correlation among multiple measures and is based on the correlations among those measures. To the extent that the first few principal components account for much of the combined variance in a set of variables, those components serve as a useful approximation of the overall correlation patterns. Principal components are described in terms of their correlations with their constituent variables, and the pattern represented by a component can be understood in terms of these correlations. Furthermore, participants' original data can be transformed to scores on each principal component. We can then analyze a few component scores rather than a large number of original measures, for example, to examine group differences. This reduces multiple comparisons and ideally contributes to insight regarding systematic patterns that may be related to underlying biology and pathology.

Initially, we observed that separate PCA of the controls and of the combined preHD and earlyHD group yielded very similar results (see Supporting Information Table S2 for further details). Subsequently, the final PCA was performed on all three groups combined. The data were first screened for univariate outliers with the only cases identified from tracts subsequently removed from the analyses due to general QC issues as described above. As a check of stability and potential multivariate outlier influence, PCA was also performed using Hubert's robust PCA (Hubert, Rousseeuw, & Vanden Branden, 2005) as implemented in the R package *rrcov* (Todorov & Filzmoser, 2009). These results (not shown) were in excellent agreement with conventional PCA and influence plots suggested no substantial outliers. Differences in PC scores among groups were assessed via ANOVA. Group differences and their statistical significance were calculated from the ANOVA-based means and standard errors. Statistical significance was defined as *p* < .05, after

Tukey-Kramer adjustment for multiple group comparisons. Predictive models for PC scores were all estimated by linear regression. We used logistic regression for multivariate prediction of earlyHD versus preHD clinical status. We used Pearson correlation (equivalent to linear regression), with and without control for age and CAG, to assess relationships between DTI scores and cognitive or motor measures. With the exception noted above, all analyses were performed with SAS/STAT 14.1 (SAS Institute Inc., Cary, NC).

## 3 | RESULTS

We investigated diffusivity in 32 tracts within networks likely to be affected as part of the HD phenotype: sensorimotor, cognitive, limbic (neuropsychiatric), and visual (Table 1; Figure 1). We employed PCA in all participants combined to identify diffusivity patterns and addressed three main questions. First, we examined the diffusivity patterns across the 32 white matter tracts and asked whether such patterns differed between controls, preHD and earlyHD. Second, we aimed to distinguish diffusivity patterns linked to key contributors to HD pathology, e.g., CAG repeat length, from those that are independent but still influence HD manifestation. Third, we investigated the relative contribution of the three diffusivity measures i.e., RD, AD, and FA to each pattern of diffusivity or principal component. Finally, we examined the relationship between each PC and both cognitive and motor performance as indexed by the global cognitive composite and grip force, respectively.

### 3.1 | Diffusivity patterns (PCA) across both controls and HD gene-carriers

We first examined diffusivity patterns, i.e., PC (principal components) across all 32 tracts of the combined control and HD gene-carrier groups. PCA revealed that the first three principal components
<figure><img src="image_5.png" /><figcaption>FIGURE 3 Principal component analysis of diffusivity. Heat map of correlation coefficients for each measure of diffusivity with dimensions derived from PCA performed in controls and HD participants. For each tract, there are two rows within the heat map; the top represents the left hemisphere, the bottom the right hemisphere. Abbreviations: DLPFC, dorsolateral prefrontal cortex; IFOF, inferior fronto-occipital fasciculus; M1, primary motor cortex; PMC, premotor cortex; PPC, posterior parietal cortex; S1, primary sensorimotor cortex; V1, primary visual cortex [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

explained 56.8% of the overall variance in diffusivity (Figure 2). Thus, the three corresponding patterns of diffusion metric correlations describe a substantial portion of the overall diffusion metric variability within the underlying white-matter regions. As noted in the Methods, the first three PCs of the joint data closely resembled the corresponding PCs when the controls and HD groups were analyzed separately. The similarity with the PCA of controls alone suggests that the overall variability in the data forms patterns similar to the patterns of correlated natural variability within the healthy white matter of the controls. Nonetheless, within the joint analysis there were notable differences in the mean PC scores. On the first and third PCs, the early HD group had significantly different mean scores compared to the preHD and control groups (see below).

For PC1, there was a pattern across all tracts which demonstrated an inverse relationship between FA and the two diffusivity measures, AD and RD, i.e., the lower the FA, the higher the AD and RD and vice-versa (Figure 3; Table 2). This was particularly the case in connections from the PPC to the thalamus, limbic tracts (the mid and posterior cingulum and the uncinate fasciculus), visual tracts (IFOF, V1 to thalamus), the corpus callosum and the ILF. Additionally, AD and RD were correlated with PC1 in all sensorimotor tracts (bilateral M1, PMC, and S1 to putamen and thalamus) and cognitive tracts (DLPFC to caudate and

thalamus, respectively). However, FA reductions in the sensorimotor and cognitive tracts were not as high as those in the limbic and visual tracts, the corpus callosum or the ILF (Figure 3). PC1 scores distinguished clinical groups. Scores were higher in early HD compared to controls (0.61 versus -0.24, *p* = .0002) and to the preHD group (0.61 versus 0.05, *p* = .0395) while the scores of controls and preHD were similar (-0.24 versus 0.04, *p* = .238). Thus, the pattern of diffusivity in PC1 was significantly more pronounced in earlyHD (Figure 3).

PC2 revealed a diffusivity pattern of increased FA and AD focused mainly on sensorimotor tracts (Figure 3; Table 2), while RD showed little association. PC2 scores were significantly higher in earlyHD than controls (0.48 versus -0.19, *p* = .006) and also higher than those in preHD, although the difference was not statistically significant (0.48 versus -0.04, *p* = .148). Again, there was no significant difference in scores between controls and preHD (-0.19 versus 0.04, *p* = .395). Thus, the pattern of diffusivity in PC2 was significantly more pronounced in earlyHD.

Finally, PC3 revealed a diffusivity pattern of increased FA and reduced RD in bilateral connections between the DLPFC and caudate and thalamus, respectively (Table 2), while increased AD correlated with PC3 in the right mid and bilateral posterior cingulum. PC3 scores were lower in earlyHD than both controls (-0.06 versus 0.03 *p* = .0016) and preHD (-0.06 versus 0.34, *p* = .0001) and there was no significant difference between controls and preHD (0.03 versus 0.34, *p* = .08). Thus the pattern of diffusivity in PC3 was lesser in earlyHD.

## 3.2 | Associations between diffusivity patterns and HD pathology

As age and CAG-repeat length interactions influence HD pathogenesis, we used linear regression to test whether age and CAG-repeat length were associated with our diffusivity patterns. In the combined preHD and early HD group, the diffusivity pattern evident in PC1 was strongly predicted by both age (0.083 (0.02), *t* = 4.16, *p* = .0001) and CAG-repeat length (0.291 (0.08), *t* = 3.62, *p* = .0005) but not an age-by-CAG interaction. There was a suggestion of a relationship between PC2 and CAG-repeat length (0.094 (0.06), *t* = 1.71, *p* = .092), but no relationship with age; PC3 was predicted by age (0.039 (0.01), *t* = 2.76, *p* = .007) but not CAG-repeat length. For controls, PC3 was significantly predicted by age (-0.027 (0.01), *t* = 2.52, *p* = .014), and there was a suggestion of age related difference in PC1 scores (0.02 (0.01), *t* = 1.86, *p* = .07). PC2 scores were not predicted by age.

Brain volume loss is a hallmark of HD progression and so we assessed whether diffusivity patterns distinguished preHD from earlyHD independently of whole brain, gray matter, white matter, putamen, caudate, and ventricular volumes. A preliminary model selection demonstrated that only the caudate and lateral ventricle volumes acted as joint predictors of HD status. Adding PC scores to that model, the diffusivity pattern of PC3 was also a significant predictor of disease status (log odds ratio = 1.078 per SD (0.410), *X*² = 6.90, *p* = .009), but PC1 (log odds ratio = 0.166 per SD (0.38) *X*² = 0.192, *p* = .66) and PC2 (log odds ratio = 0.448 per SD (0.4) *X*² = 1.13, *p* = .29) were not.
<table><caption>TABLE 2 PCA results [Color table can be viewed at wileyonlinelibrary.com]</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="3">PC1</th><th colspan="3">PC2</th><th colspan="3">PC3</th></tr><tr><th>FA</th><th>AD</th><th>RD</th><th>FA</th><th>AD</th><th>RD</th><th>FA</th><th>AD</th><th>RD</th></tr></thead><tbody><tr><td rowspan="2">M1-Motor Thalamus</td><td>L</td><td></td><td>0.4944</td><td>0.71558</td><td>0.67955</td><td>0.70197</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.43883</td><td>0.67178</td><td>0.64759</td><td>0.71228</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">M1-Putamen</td><td>L</td><td></td><td>0.59931</td><td>0.77047</td><td>0.78749</td><td>0.61072</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.5695</td><td>0.66923</td><td>0.77874</td><td>0.66393</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">PMC-Motor Thalamus</td><td>L</td><td></td><td>0.43661</td><td>0.61912</td><td>0.73059</td><td>0.71341</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.45098</td><td>0.64774</td><td>0.68891</td><td>0.68221</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">PMC-Putamen</td><td>L</td><td></td><td>0.63268</td><td>0.74575</td><td>0.79596</td><td>0.58385</td><td>-0.43091</td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.59244</td><td>0.67459</td><td>0.77046</td><td>0.55995</td><td>-0.46135</td><td></td><td></td><td></td></tr><tr><td rowspan="2">S1-Somatosensory Thalamus</td><td>L</td><td></td><td>0.51495</td><td>0.72899</td><td>0.70752</td><td>0.68283</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.46951</td><td>0.66197</td><td>0.61133</td><td>0.70453</td><td>-0.43991</td><td></td><td></td><td>0.46866</td></tr><tr><td rowspan="2">S1-Putamen</td><td>L</td><td></td><td>0.59981</td><td>0.77641</td><td>0.77865</td><td>0.61116</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.62317</td><td>0.76785</td><td>0.75927</td><td>0.60292</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">DLPFC-Prefrontal Thalamus</td><td>L</td><td></td><td>0.69285</td><td></td><td></td><td></td><td></td><td>0.65991</td><td></td><td>-0.67064</td></tr><tr><td>R</td><td></td><td>0.6591</td><td></td><td></td><td></td><td></td><td>0.74847</td><td></td><td>-0.74165</td></tr><tr><td rowspan="2">DLPFC-Caudate</td><td>L</td><td></td><td>0.68251</td><td>0.68763</td><td>0.48431</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.66408</td><td>0.69064</td><td>0.47344</td><td></td><td></td><td>0.59726</td><td></td><td>-0.42632</td></tr><tr><td rowspan="2">PPC - Parietal Thalamus</td><td>L</td><td>-0.57012</td><td>0.6499</td><td>0.83471</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.54872</td><td>0.55881</td><td>0.75936</td><td></td><td>0.41155</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Uncinate Fasciculus</td><td>L</td><td>-0.4484</td><td>0.43608</td><td>0.61483</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td></td><td>0.50387</td><td>0.50397</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Mid Cingulum</td><td>L</td><td>-0.54713</td><td>0.47198</td><td>0.71255</td><td></td><td>0.4385</td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.50473</td><td>0.56258</td><td>0.77297</td><td></td><td>0.40841</td><td></td><td></td><td>0.4084</td><td></td></tr><tr><td rowspan="2">Posterior Cingulum</td><td>L</td><td>-0.5396</td><td></td><td>0.76733</td><td></td><td>0.42102</td><td></td><td></td><td>0.44304</td><td></td></tr><tr><td>R</td><td>-0.53816</td><td></td><td>0.76758</td><td></td><td>0.4282</td><td></td><td></td><td>0.41278</td><td></td></tr><tr><td rowspan="2">V1 - Visual Thalamus</td><td>L</td><td>-0.52886</td><td>0.44216</td><td>0.74869</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.53247</td><td></td><td>0.65093</td><td></td><td>0.44095</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">IFOF</td><td>L</td><td>-0.65731</td><td>0.5482</td><td>0.84464</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.60571</td><td>0.49602</td><td>0.81886</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">ILF</td><td>L</td><td>-0.63733</td><td>0.50695</td><td>0.77461</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.61165</td><td>0.47702</td><td>0.79844</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mid Corpus Callosum</td><td></td><td>-0.67649</td><td>0.57568</td><td>0.84759</td><td></td><td>0.41806</td><td></td><td></td><td></td><td></td></tr><tr><td>Posterior Corpus Callosum</td><td></td><td>-0.6755</td><td>0.61222</td><td>0.82682</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Correlation co-efficients for each PC showing the contribution of each measure of diffusivity to each PC for all tracts; for visualization purposes, these have only been shown for co-efficients that are smaller than or greater than 0.4. Tracts are color-coded according to the network to which they putatively belong. Blue: Sensorimotor Network tracts; Green: Cognitive network tracts; Pink: Limbic network tracts; Orange: Visual network tracts; Purple: Control (non-HD phenotype); Brown: Interhemispheric tracts. AD: axial diffusivity; DLPFC: dorsolateral prefrontal cortex; FA: fractional anisotropy; IFOF: inferior fronto-occipital Fasciculus; ILF: inferior lateral Fasciculus; M1: primary motor cortex; PMC: premotor cortex; PPC: posterior parietal cortex; RD: radial diffusivity; S1: primary somatosensory cortex; V1: primary visual cortex.

## 3.3 | Contribution of diffusivity measures to principal component patterns

PC coefficients represent correlations between diffusivity measures and their contributions to each pattern of diffusivity. We investigated

these associations in a subset of representative tracts for each PC, plotting diffusivity metric values against PC scores with individual data points color-coded for group membership. As discussed earlier, there was considerable overlap in distribution across the groups (Figure 4).
<figure><img src="image_6.png" /><figcaption>FIGURE 4 Relationships between PC scores. Scatterplots showing the relationships between individual PC scores for manifest HD (red dots) and preHD (blue crosses) (A-C) and manifest HD and controls (green crosses) (D-F). Manifest HD participants tend to have higher PC1 and PC2 scores, and lower PC3 scores, than preHD or control participants [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

For PC1, all three groups displayed correlations between PC1 scores and RD and AD in sensorimotor (bilateral PMC and both thalamus and putamen, respectively) and visual network (bilateral IFOF) tracts plus bilateral ILF (Figure 3; Table 2), suggesting that both AD and RD contributed strongly to PC1 regardless of the type of tract (Figure 3). Dif- fusivity values and mean scores were highest in earlyHD (Figure 4). For PC2 all three groups displayed a strong positive correlation between PC2 scores and both FA and AD in bilateral PMC and both thalamic and putaminal tracts with a weak negative correlation with RD. The dif- ferentiation between early HD and the other groups is more difficult to appreciate. Finally, for PC3, all three groups displayed a moderate neg- ative correlation between RD and PC3 with the earlyHD group show- ing higher levels of RD and lower PC3 scores. AD also correlated moderately positively with PC3 scores in the bilateral mid and posterior cingulum and here, the earlyHD group had lower levels of FA and AD and lower PC scores compared to the control and preHD groups (Figure 4).

3.4 | Relationship between diffusivity patterns and behavior

Global cognitive composite scores and grip force index (as a marker of motor performance) were correlated with overall PC scores for both controls and a combined HD group (Table 3). The diffusivity pattern of PC1 was negatively correlated with global cognition and positively cor- related with motor performance for the combined HD group and to a lesser extent the control group, suggesting that widespread increased diffusivity was associated with diminished cognitive and motor per- formance. Associations for cognition were only slightly attenuated

when controlling for the age-CAG interaction and thus cannot be explained by the common causal influence of the HD gene. The sensorimotor-based diffusivity pattern for PC2 was positively corre- lated with motor performance only in both groups with little effect of age or CAG. Finally, the cognition-based diffusivity pattern for PC3, was not associated with either cognition or motor performance for either controls or HD.

Looking at specific tracts, there were correlations between increased AD in PMC tracts, and reduced motor performance bilaterally for HD and in the right hemisphere for controls, that were little dimin- ished when controlling for age and CAG repeat length. RD in the DLPFC—caudate tract was positively correlated with the global cogni- tion in both groups. This was attenuated when controlling for age and CAG and was large enough in controls to suggest that changes may be due to ageing. However, many of the attenuated correlations remained significant in HD. There was also a similar correlation for the DLPFC– thalamic tract for controls but not HD.

4 | DISCUSSION

In this study, we examined the relationships between AD, RD, and FA within white matter tracts in HD gene expansion mutation-carriers and healthy controls and identified three independent patterns of diffusiv- ity relationships common to both. Absolute diffusivity levels varied between individuals, while the relationship between AD, RD, and FA was maintained potentially representing natural biological variation of white matter microstructure. All three patterns were more pronounced in HD compared to controls; with two closely related to biological fac- tors associated with the HD gene-mutation. The first described a
<table><caption>TABLE 3 Correlations between PC scores and performance</caption><thead><tr><th></th><th colspan="2">HD</th><th colspan="2">Controls</th></tr><tr><th></th><th></th><th>Adjusted</th><th></th><th>Adjusted</th></tr></thead><tbody><tr><td><strong>Global Cognitive Score</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>PC1 Score (74)</td><td>-0.449</td><td>-0.396</td><td>-0.199</td><td>-0.121</td></tr><tr><td><em>p val</em></td><td>&lt;.0001</td><td>0.0006</td><td>0.0787</td><td>0.293</td></tr><tr><td>PC2 Score (74)</td><td>-0.133</td><td>-0.130</td><td>0.106</td><td>0.076</td></tr><tr><td><em>p val</em></td><td>0.2593</td><td>0.2785</td><td>0.354</td><td>0.5063</td></tr><tr><td>PC3 Score (74)</td><td>0.093</td><td>0.057</td><td>0.227</td><td>0.119</td></tr><tr><td><em>p val</em></td><td>0.4315</td><td>0.6396</td><td>0.0444</td><td>0.3008</td></tr><tr><td>RD: Left DLPFC-Caud(90)</td><td>-0.304</td><td>-0.230</td><td>-0.311</td><td>-0.141</td></tr><tr><td><em>p val</em></td><td>0.0035</td><td>0.0321</td><td>0.0036</td><td>0.2012</td></tr><tr><td>RD: Right DLPFC-Caud(90)</td><td>-0.284</td><td>-0.198</td><td>-0.271</td><td>-0.111</td></tr><tr><td><em>p val</em></td><td>0.0066</td><td>0.066</td><td>0.0115</td><td>0.3135</td></tr><tr><td>RD: Left DLPFC-Thal(90)</td><td>-0.165</td><td>-0.103</td><td>-0.282</td><td>-0.210</td></tr><tr><td><em>p val</em></td><td>0.1204</td><td>0.3422</td><td>0.009</td><td>0.0547</td></tr><tr><td>RD: Right DLPFC-Thal(90)</td><td>-0.101</td><td>-0.039</td><td>-0.218</td><td>-0.126</td></tr><tr><td><em>p val</em></td><td>0.3412</td><td>0.7196</td><td>0.0439</td><td>0.254</td></tr><tr><td><strong>Grip Force</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>PC1 Score (74)</td><td>0.361</td><td>0.349</td><td>0.196</td><td>0.203</td></tr><tr><td><em>p val</em></td><td>0.0016</td><td>0.0028</td><td>0.0839</td><td>0.0742</td></tr><tr><td>PC2 Score (74)</td><td>0.253</td><td>0.232</td><td>0.219</td><td>0.218</td></tr><tr><td><em>p val</em></td><td>0.0297</td><td>0.0517</td><td>0.053</td><td>0.0551</td></tr><tr><td>PC3 Score (74)</td><td>-0.015</td><td>0.006</td><td>0.117</td><td>0.117</td></tr><tr><td><em>p val</em></td><td>0.8993</td><td>0.9578</td><td>0.306</td><td>0.3077</td></tr><tr><td>AD: Left PMC-Thal(74)</td><td>0.420</td><td>0.397</td><td>0.168</td><td>0.171</td></tr><tr><td><em>p val</em></td><td>0.0002</td><td>0.0006</td><td>0.1377</td><td>0.1334</td></tr><tr><td>AD: Right PMC-Thal(90)</td><td>0.318</td><td>0.265</td><td>0.267</td><td>0.269</td></tr><tr><td><em>p val</em></td><td>0.0022</td><td>0.0257</td><td>0.0131</td><td>0.0127</td></tr><tr><td>AD: Left PMC-Put(90)</td><td>0.407</td><td>0.464</td><td>0.141</td><td>0.146</td></tr><tr><td><em>p val</em></td><td>&lt;.0001</td><td>&lt;.0001</td><td>0.1946</td><td>0.1825</td></tr><tr><td>AD: Right PMC-Put(90)</td><td>0.297</td><td>0.237</td><td>0.204</td><td>0.208</td></tr><tr><td><em>p val</em></td><td>0.0044</td><td>0.0464</td><td>0.0594</td><td>0.0557</td></tr><tr><td>RD: Left PMC-Thal(74)</td><td>0.138</td><td>0.112</td><td>0.064</td><td>0.064</td></tr><tr><td><em>p val</em></td><td>0.2415</td><td>0.3507</td><td>0.5757</td><td>0.5791</td></tr><tr><td>RD: Right PMC-Thal(90)</td><td>0.060</td><td>0.105</td><td>0.025</td><td>0.031</td></tr><tr><td><em>p val</em></td><td>0.5734</td><td>0.3857</td><td>0.8164</td><td>0.7758</td></tr><tr><td>RD: Left PMC-Put(90)</td><td>0.126</td><td>0.147</td><td>0.108</td><td>0.128</td></tr><tr><td><em>p val</em></td><td>0.2355</td><td>0.2216</td><td>0.3223</td><td>0.2448</td></tr><tr><td>RD: Right PMC-Put(90)</td><td>0.094</td><td>0.122</td><td>0.009</td><td>0.023</td></tr><tr><td><em>p val</em></td><td>0.378</td><td>0.3096</td><td>0.9324</td><td>0.8334</td></tr></tbody></table>

Pearson correlations between behavioral measures and PC Scores and individual diffusivity measures in specified tracts for HD and Control groups. Correlations are presented unadjusted and adjusted for each group; for HD the correlations are adjusted for Age, CAG, Age*CAG, for Controls the correlations are adjusted for Age. AD, axial diffusivity; DLPFC, Dorsolateral Prefrontal Cortex; PMC, premotor cortex; Put, Putamen; RD, radial diffusivity.
widespread pattern of increased AD, RD and reduced FA, while the second was restricted to increased AD and FA in sensorimotor network tracts. In the case of the HD group, these patterns were to some extent driven by CAG repeat length, the primary predictor of HD onset, indicating that HD pathogenesis may accentuate pre-existing diffusivity patterns that are independent of the HD gene mutation and also present in controls. Further, both patterns were linked to widespread HD-related brain volume loss reflecting key elements of the structural phenotype of HD. The third pattern, consisting mainly of RD changes in tracts associated with cognitive function, however, while associated with age was independent of both HD-mediated volume loss and CAG repeat length, but nonetheless, a predictor of HD diagnostic status independently contributing to HD clinical manifestations.

In this study, we examined diffusivity across 32 white matter tracts from a number of functional networks and identified three independent patterns of AD, RD and FA relationships common to controls, preHD and earlyHD gene-carriers. Having previously shown that in left-sided sensorimotor tracts diffusivity naturally varies for both controls and HD gene-carriers independently i.e., absolute diffusivity levels differs between individuals but the relationship between AD and RD is maintained (Orth et al., 2016), the current findings extend this concept of variability within naturally-occurring patterns of diffusivity in HD to include pathways within cognitive, limbic and visual networks. This subsequently supports the idea that HD pathogenesis does not necessarily abolish pre-existing patterns of diffusivity, or even generate new ones, but instead modifies existing patterns throughout the brain.

The patterns of diffusivity that we previously identified in the sensorimotor network, not only reflect natural variability but are also influenced by HD pathology and could be clinically relevant since they were associated with manifestations of unequivocal motor signs of HD (Orth et al., 2016). In the current study, therefore, we explored whether the exacerbation of such diffusivity patterns could potentially be considered a factor in HD pathology. Despite the commonalities between diffusivity patterns across the three groups, it is important to note that the distribution of the data along the PC axes, which describe the relationship between data, i.e., the pattern, differed between the earlyHD and control groups. Thus, differences in PC scores between groups indicate their placement at opposite ends of the axis. Indeed all three diffusivity patterns were accentuated in earlyHD, while remaining broadly similar for controls and preHD. For example, for PC1, where we saw a widespread pattern of increased AD and RD and to a lesser extent decreased FA across all 32 white matter tracts, the earlyHD group displayed higher scores mainly due to higher levels of RD compared to both controls and preHD. Similarly, the earlyHD group showed higher levels of both RD in cognitive tracts (PC3) and AD in bilateral sensorimotor tracts (PC2). In all cases, differences were only found in earlyHD, and there were no detectable differences in the preHD group. The evidence of both widespread and localized changes in diffusivity over and above those patterns seen in controls and preHD suggests that in the manifest stages of HD there is ongoing disorganization in white matter pathways which underscores a series of functions including motor and cognitive – both key elements of the phenotype. It is interesting to note that while there were no

differences between preHD and controls, the direction of diffusivity changes in these naturally-occurring patterns in the preHD was comparable with that seen in earlyHD. It is likely that the white matter changes in the earlyHD group begin many years prior to diagnosis but are not necessarily detectable in comparisons with non-CAG expanded individuals. There is robust evidence to suggest that macrostructural white matter changes in the cortex occur prior to those in the gray matter with similar white matter microstructural changes (Paulsen et al., 2010; Tabrizi et al., 2012; Wu et al., 2017).

We also examined the extent to which the diffusivity patterns predict HD disease status independent of genetic and structural markers of HD progression. As such, despite significantly discriminating between preHD and earlyHD, both patterns of diffusivity identified in PC1 and PC2 were closely associated with CAG repeat length and whole brain, gray matter, white matter, putamen, and caudate volume loss and therefore, not independently associated with HD clinical group classifications. Here, it would appear that widespread increased radial diffusivity (as seen in PC1) and increased axial diffusivity in the sensorimotor network (as seen in PC2) were both closely related to the overall effects of HD pathogenesis, ultimately contributing to the structural phenotype of HD.

The third DLPFC-based diffusivity pattern associated with cognitive function did, however, remain associated with disease status even when controlling for volumetric markers of HD pathology and CAG repeat length. Changes in this pattern may be due to factors that do not cause HD and the diffusivity properties represented may instead independently influence the effects of the HD mutation. As this particular pattern of diffusivity was associated with age for both controls and HD gene-carriers, changes here may be related to normal aging. In other words, factors associated with aging may influence the diffusivity properties of frontal lobe white matter tracts involved in networks subserving cognition. In healthy individuals, frontal lobe gray and white matter structure and connectivity have been associated with age-related cognitive performance loss (Kievit et al., 2014; Serbruyns et al., 2016; Zhao et al., 2015) that may be associated with reduced myelination of the frontal lobe. Given the known benefits of cognitive training on myelination in healthy individuals (Caeyenberghs et al., 2016; Takeuchi et al., 2010), similar therapeutic interventions in HD gene-carriers may counteract the effects of age and thus positively influence HD pathogenesis before disease onset, akin to the disease-modifying effects of environmental enrichment in HD animal models (Dersi et al., 2016; Nithianantharajah and Hannan, 2006; Xu et al., 2016). Although this pattern of diffusivity correlated with global cognition in controls only, increased RD in the DLPFC-caudate tract specifically correlated with cognitive performance in both controls and HD. Furthermore, this effect was reduced once correcting for age and CAG suggesting that while influenced by natural variability, the presence of the HD gene mutation has an independent effect that requires further investigation.

HD may not, however, affect all aspects of white matter microstructure in equal measure. The first pattern of widespread change was driven mainly by RD suggesting that the HD gene-mutation may contribute to widespread demyelination (Fox et al., 2011), which given the correlation of this pattern with age in both controls and gene-carriers
could be an amplification of a pattern seen in normal ageing. This over-
all effect is reflected in the correlation between the PC score and both
cognitive and motor performance. The difference in the strength of
these associations echoes the accentuated effect of naturally-occurring
age-related increased diffusivity in those with the HD gene. Increased
AD in sensorimotor tracts, coupled with an apparent inverse relation-
ship between RD and AD was similar to that identified previously (Orth
et al., 2016). Although data were from the same cohort and therefore,
some correspondence in findings is expected, this pattern is neverthe-
less consistent across both hemispheres of the sensorimotor network
with earlyHD showing the highest levels of AD. Again, this was
reflected in the strong correlations between increased AD, not RD, in
the sensorimotor tracts and motor performance. Given that motor dys-
function is one of the key clinical presentations of HD, it would appear
that AD increases are associated with motor deterioration seen in HD
(Müller et al., 2016) and may reflect degeneration in the main fiber
direction within motor pathways and subsequent disorganization lead-
ing to deficits in motor performance. In both cases however, we must
proceed with caution given the considerable uncertainty associated
with the biological interpretation of diffusivity changes (Jones et al.,
2013). Taken together, our data suggest that HD pathogenesis
increases RD across widespread networks, and AD in sensorimotor
connections as part of a wider structural phenotype of earlyHD. Inde-
pendent of key contributors to HD pathogenesis, RD increases in cog-
nitive network-based connections distinguished earlyHD from preHD
or controls. This suggests that the HD mutation exerts effects on top
of naturally-occurring white matter tract microstructural variability. In
addition, age-related effects on cognitive tract microstructure inde-
pendently modify HD pathogenesis, the identification of which may be
valuable in understanding disease modification.

ACKNOWLEDGMENTS

The authors thank the Track-On HD study participants, the CHDI/
High Q Foundation, a not-for-profit organization dedicated to find-
ing treatments for HD.

AUTHOR CONTRIBUTIONS

Track-On HD Investigators: A Coleman, J Decolongon, M Fan, T
Koren (University of British Columbia, Vancouver); C Jauffret, D
Justo, S Lehericy, K Nigaud, R Valabrègue (ICM and APHP, Pitié-
Salpêtrière University Hospital, Paris); S Klöppel, E Scheller, L Min-
kova (University of Freiberg, Freiberg), A Schoonderbeek, E P 't Hart
(Leiden University Medical Centre, Leiden); C Berna, M Desikan, R
Ghosh, D Hensman Moss, E Johnson, P McColgan, G Owen, M
Papoutsi, A Razi, J Read, (University College London, London); D
Craufurd (Manchester University, Manchester); R Reilmann, N
Weber (George Huntington Institute, Munster); H Johnson, J D
Long, J Mills (University of Iowa, Iowa); J Stout, I Labuschagne
(Monash University, Melbourne); G B Landwehrmeyer, I Mayer (Ulm
University, Ulm).

ORCID

Sarah Gregory http://orcid.org/0000-0003-4620-6963

Rachael I. Scahill http://orcid.org/0000-0002-1832-106X

Michael Orth http://orcid.org/0000-0002-6913-3536

REFERENCES

Behrens, T. E., Berg, H. J., Jbabdi, S., Rushworth, M. F., & Woolrich, M. W. (2007). Probabilistic diffusion tractography with multiple fibre orientations: What can we gain?. *NeurolImage*, 34(1), 144–155.

Behrens, T. E., Woolrich, M. W., Jenkinson, M., Johansen-Berg, H., Nunes, R. G., Clare, S., . . . Smith, S. M. (2003). Characterization and propagation of uncertainty in diffusion-weighted MR imaging. *Magnetoic Resonance in Medicine*, 50(5), 1077–1088.

Caeyenberghs, K., Metzler-Baddeley, C., Foley, S., & Jones, D. K. (2016). Dynamics of the human structural connectome underlying working memory training. *Journal of Neuroscience*, 36, 4056–4066.

Della Nave, R., Ginestroni, A., Tessa, C., Giannelli, M., Piacentini, S., Filipppi, M., & Mascalchi, M. (2010). Regional distribution and clinical correlates of white matter structural damage in Huntington disease: a tract-based spatial statistics study. *American Journal of Neuroradiology*, 31(9), 1675–1681.

Dezsi, G., Ozturk, E., Salzberg, M. R., Morris, M., O'Brien, T. J., & Jones, N. C. (2016). Environmental enrichment imparts disease-modifying and transgenerational effects on genetically-determined epilepsy and anxiety. *Neurobiology of Disease*, 93, 129–136.

Douaud, G., Behrens, T. E., Poupon, C., Cointepas, Y., Jbabdi, S., Gaura, V., . . . Remy, P. (2009). In vivo evidence for the selective subcortical degeneration in Huntington's disease. *NeurolImage*, 46(4), 958–966.

Dumas, E. M., van den Bogaard, S. J., Ruber, M. E., Reilman, R. R., Stout, J. C., Craufurd, D., . . . Roos, R. A. (2012). Early changes in white matter pathways of the sensorimotor cortex in premanifest Huntington's disease. *Human Brain Mapping*, 33(1), 203–212.

Duyao, M., Ambrose, C., Myers, R., Novelletto, A., Persichetti, F., Frontali, M., . . . Abbott, M. (1993). Trinucleotide repeat length instability and age of onset in Huntington's disease. *Nature Genetics*, 4(4), 387–392.

Fox, R. J., Cronin, T., Lin, J., Wang, X., Sakaie, K., Ontaneda, D., . . . Philips, M. D. (2011). Measuring myelin repair and axonal loss with diffusion tensor imaging. *American Journal of Neuroradiology*, 32, 85–91.

GeM-HD Consortium. (2015). Identification of genetic factors that modify clinical onset of Huntington's disease. *Cell*, **162**, 516–526.

Glover, G. H., Mueller, B. A., Turner, J. A., van Erp, T. G., Liu, T. T., Greve, D. N., . . . Potkin, S. G. (2012). Function biomedical informatics research network recommendations for prospective multicenter functional MRI studies. *Journal of Magnetic Resonance Imaging*, 36(1), 39–54.

Hubert, M., Rousseeuw, P. J., & Vanden Branden, K. (2005). ROBPCA: A new approach to robust principal component analysis. *Technometrics*, 47(1), 64–79.

Jenkinson, M., & Smith, S. (2001). A global optimisation method for robust affine registration of brain images. *Medical Image Analysis*, 5 (2), 143–156.

Jones, D. K., Knosche, T. R., & Turner, R. (2013). White matter integrity, fiber count, and other fallacies: the do's and don'ts of diffusion MRI. *NeuroImage*, 73, 239–254.

Jones, R., Stout, J. C., Labuschagne, I., Say, M., Justo, D., Coleman, A., . . . Frost, C. (2014). The potential of composite cognitive scores for tracking progression in Huntington's disease. *Journal of Huntington's Disease*, 3, 197–207.

Kievit, R. A., Davis, S. W., Mitchell, D. J., Taylor, J. R., Duncan, J., Henson, R. N., . . . Cam-CAN Research Team (2014). Distinct aspects of
frontal lobe structure mediate age-related differences in fluid intelligence and multitasking. *Nature Communications*, **5**, 5658.

Klöppel, S., Draganski, B., Golding, C. V., Chu, C., Nagy, Z., Cook, P. A., ... Frackowiak, R. S. J. (2008). White matter connections reflect changes in voluntary-guided saccades in pre-symptomatic Huntington's disease. *Brain*, **131**(1), 196-204.

Kloppel, S., Gregory, S., Scheller, E., Minkova, L., Razi, A., Durr, A., ... Tabrizi, S. J. (2015). Compensation in preclinical Huntington's disease: Evidence from the Track-On HD study. *EBioMedicine*, **2**(10), 1420-1429.

Lee, J. M., Ramos, E. M., Lee, J. H., Gillis, T., Mysore, J. S., Hayden, M. R., ... Gusella, J. F. (2012). CAG repeat expansion in Huntington disease determines age at onset in a fully dominant fashion. *Neurology*, **78**(10), 690-695.

Maldjian, J. A., Laurienti, P. J., Kraft, R. A., & Burdette, J. H. (2003). An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. *NeurolImage*, **19**(3), 1233-1239.

Matsui, J. T., Vaidya, J. G., Johnson, H. J., Magnotta, V. A., Long, J. D., Mills, J. A., ... Paulsen, J. S. (2014). Diffusion weighted imaging of prefrontal cortex in prodromal huntington's disease. *Human Brain Mapping*, **35**(4), 1562-1573.

Müller, H.-P., Gorges, M., Grön, G., Kassubek, J., Landwehrmeyer, G. B., Süßmuth, S. D., ... Orth, M. (2016). Motor network structure and function are associated with motor performance in Huntington's disease. *Journal of Neurology*, **263**(3), 539-549.

Nithianantharajah, J., & Hannan, A. J. (2006). Enriched environments, experience-dependent plasticity and disorders of the nervous system. *Nature Reviews. Neuroscience*, **7**, 697-709.

Novak, M. J. U., Seunarine, K. K., Gibbard, C. R., Hobbs, N. Z., Scahill, R. I., Clark, C. A., & Tabrizi, S. J. (2014). White matter integrity in pre-manifest and early Huntington's disease is related to caudate loss and disease progression. *Cortex*, **52**, 98-112.

Odish, O. F., Leemans, A., Reijntjes, R. H., van den Bogaard, S. J., Dumas, E. M., Wolterbeek, R., ... Roos, R. A. (2015). Microstructural brain abnormalities in Huntington's disease: A two-year follow-up. *Human Brain Mapping*, **36**(6), 2061-2074.

Orth, M., Gregory, S., Scahill, R. I., Mayer, I. S., Minkova, L., Kloppel, S., ... TRACK-On Investigators. (2016). Natural variation in sensory-motor white matter organization influences manifestations of Huntington's disease. *Human Brain Mapping*, **37**, 4615-4628.

Paulsen, J. S., Long, J. D., Ross, C. A., Harrington, D. L., Erwin, C. J., Williams, J. K., ... Barker, R. A. Investigators, P.-H., Coordinators of the Huntington Study, G. (2014). Prediction of manifest Huntington's disease with clinical and imaging measures: a prospective observational study. *Lancet Neurology*, **13**(12), 1193-1201.

Paulsen, J. S., Nopoulos, P. C., Aylward, E., Ross, C. A., Johnson, H., Magnotta, V. A., ... Investigators, P.-H. Coordinators of the Huntington's Study, G. (2010). Striatal and white matter predictors of estimated diagnosis for Huntington disease. *Brain Research Bulletin*, **82**(3-4), 201-207.

Poudel, G. R., Stout, J. C., Domínguez D, J. F., Salmon, L., Churchyard, A., Chua, P., ... Egan, G. F. (2014). White matter connectivity reflects clinical and cognitive status in Huntington's disease. *Neurobiology of Disease*, **65**, 180-187.

Poudel, G. R., Stout, J. C., Domínguez, D. J., Churchyard, A., Chua, P., Egan, G. F., & Georgiou-Karistianis, N. (2015). Longitudinal change in white matter microstructure in Huntington's disease: The IMAGE-HD study. *Neurobiology of Disease*, **74**, 406-412.

Serbruyns, L., Leunissen, I., van Ruitenbeek, P., Pauwels, L., Caeyenberghs, K., Solesio-Jofre, E., ... Swinnen, S. P. (2016). Alterations in brain white matter contributing to age-related slowing of task switching performance: The role of radial diffusivity and magnetization transfer ratio. *Human Brain Mapping*, **37**, 4084-4098.

Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., ... Matthews, P. M. (2004). Advances in functional and structural MR image analysis and implementation as FSL. *NeurolImage*, **23**, S208-S219.

Tabrizi, S. J., Reilmann, R., Roos, R. A., Durr, A., Leavitt, B., Owen, G., ... Langbehn, D. R. (2012). Potential endpoints for clinical trials in pre-manifest and early Huntington's disease in the TRACK-HD study: analysis of 24 month observational data. *Lancet Neurology*, **11**, 42-53.

Tabrizi, S. J., Scahill, R. I., Owen, G., Durr, A., Leavitt, B. R., Roos, R. A., ... Langbehn, D. R. (2013). 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. *The Lancet Neurology*, **12**, 637-649.

Takeuchi, H., Sekiguchi, A., Taki, Y., Yokoyama, S., Yomogida, Y., Komuro, N., ... Kawashima, R. (2010). Training of working memory impacts structural connectivity. *Journal of Neuroscience*, **30**, 3297-3303.

Todorov, V., & Filzmoser, P. (2009). An object-oriented framework for robust multivariate analysis. *Journal of Statistical Software*, **32**(3), 1-47.

Wakana, S., Caprihan, A., Panzenboeck, M. M., Fallon, J. H., Perry, M., Gollub, R. L., ... Mori, S. (2007). Reproducibility of quantitative trac-tography methods applied to cerebral white matter. *NeurolImage*, **36**(3), 630-644.

Wexler, N. S., Lorimer, J., Porter, J., Gomez, F., Moskowitz, C., Shackell, E., ... Landwehrmeyer, B. Project, U.S.-V.C.R. (2004). Venezuelan kind-dreds reveal that genetic and environmental factors modulate Huntington's disease age of onset. *Proceedings of the National Academy of Sciences of the United States of America*, **101**, 3498-3503.

Wu, D., Faria, A. V., Younes, L., Mori, S., Brown, T., Johnson, H., ... Miller, M. I. Investigators, P.-H., Coordinators of the Huntington Study, G. (2017). Mapping the order and pattern of brain structural MRI changes using change-point analysis in premanifest Huntington's disease. *Human Brain Mapping*, **38**(10), 5035-5050.

Xu, H., Gelyana, E., Rajsombath, M., Yang, T., Li, S., & Selkoe, D. (2016). Environmental enrichment potently prevents microglia-mediated neuroinflammation by human amyloid beta-protein oligomers. *Journal of Neuroscience*, **36**, 9041-9056.

Zhao, T. D., Cao, M., Zuo, X. N., Dong, Q., He, Y., & Shu, N. (2015). Age-related changes in the topological organization of white matter structural networks across the human lifespan. *International Journal of Developmental Neuroscience*, **47**, 26-27.

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the supporting information tab for this article.

**How to cite this article:** Gregory S, Crawford H, Seunarine K, et al. Natural biological variation of white matter microstructure is accentuated in Huntington's disease. *Hum Brain Mapp.* 2018;39:3516-3527. https://doi.org/10.1002/hbm.24191