Abstract
Background:
Impaired myelination may contribute to Huntington’s disease (HD) pathogenesis.
Objective:
This study assessed differences in white matter (WM) microstructure between HD patients and controls, and tested whether drumming training stimulates WM remodelling in HD. Furthermore, it examined whether training-induced microstructural changes are related to improvements in motor and cognitive function.
Methods:
Participants undertook two months of drumming exercises. Working memory and executive function were assessed before and post-training. Changes in WM microstructure were investigated with diffusion tensor magnetic resonance imaging (DT-MRI)-based metrics, the restricted diffusion signal fraction (Fr) from the composite hindered and restricted model of diffusion (CHARMED) and the macromolecular proton fraction (MPF) from quantitative magnetization transfer (qMT) imaging. WM pathways linking putamen and supplementary motor areas (SMA-Putamen), and three segments of the corpus callosum (CCI, CCII, CCIII) were studied using deterministic tractography. Baseline MPF differences between patients and controls were assessed with tract-based spatial statistics.
Results:
MPF was reduced in the mid-section of the CC in HD subjects at baseline, while a significantly greater change in MPF was detected in HD patients relative to controls in the CCII, CCIII, and the right SMA-putamen post-training. Further, although patients improved their drumming and executive function performance, such improvements did not correlate with microstructural changes. Increased MPF suggests training-induced myelin changes in HD.
Conclusion:
Though only preliminary and based on a small sample size, these results suggest that tailored behavioural stimulation may lead to neural benefits in early HD, that could be exploited for delaying disease progression.
INTRODUCTION
The pathology of Huntington’s disease
Huntington’s disease (HD) is a genetic, neurodegenerative disorder caused by an expansion of the CAG repeat within the huntingtin gene, leading to debilitating cognitive, psychiatric, and motor symptoms. In addition to striatal grey matter (GM) degeneration [1], HD pathology has been linked to white matter (WM) changes [2–9]. Additionally, an increasing body of research suggests that myelin-associated biological processes at the cellular and molecular level contribute to WM abnormalities [10–15]. Myelin is a multi-layered membrane sheath wrapping axons and is produced by oligodendrocytes. Axon myelination is vital during brain development and critical for healthy brain function [16]. Oligodendrocyte/myelin dysfunction can slow down or stop otherwise fast axonal transport, which in turn can result in synaptic loss and eventually axonal degeneration [17].
Interventions and brain plasticity
As HD is caused by a single-gene, it is an ideal model to study neurodegeneration as a whole, and test for possible beneficial interventions that can slow or suppress disease onset. Despite this, no disease-modifying treatment is approved for patients with HD at present. Recent developments in gene therapy have generated much excitement. However, these have yet to be proved to lead to measurable changes in disease progression [18]. Furthermore, a number of questions linger, for example on the relative strength of different approaches and the possible side effects of each therapy [19]. Importantly, these treatments aim to modify and not cure the disease, and while symptomatic therapies for HD are present, and are used for treating chorea and some of the psychiatric symptoms, their effectiveness varies between patients and may lead to clinically significant side-effects [18]. This stresses the need to develop better symptomatic therapies to aid patients and manage HD symptoms.
Environmental stimulation and behavioural interventions may have the potential to reduce disease progression and delay disease onset [20–22]. Furthermore, previous studies have detected training-related changes in the WM of both healthy controls [23, 24] and patients, including subjects with HD [25]. For example, DT MRI studies have shown microstructural WM changes following balance training in healthy [24] and traumatic brain injury young adults [26]. Other imaging studies have shown DT MRI changes as a result of juggling [27], abacus training [28], extensive piano practice [29, 30], working memory training [31], reasoning training [32], and meditation training [33].
Converging evidence implicates myelin plasticity as one of the routes by which experience shapes brain structure and function [27, 34–38]. Plastic changes in myelination may be implicated in early adaptation and longer-term consolidation and improvement in motor tasks [39–42]. Changes in myelin-producing oligodendrocytes and in GM and WM microstructure have been reported within the first hours of skill acquisition [43, 44], implying that experience can be quickly translated into adaptive changes in the brain.
This study assessed whether two months of drumming training, involving practising drumming patterns in ascending order of difficulty, could trigger WM microstructural changes, and potentially myelin remodelling, in individuals with HD. Specifically, we hypothesised that changes in microstructural metrics would be more marked in patients than in healthy subjects, based on reports of larger training-associated changes in structural MRI metrics in patient populations than in healthy controls [34].
The drumming intervention was designed to target cognitive and motor functions known to be mediated by cortico-basal ganglia loops. More specifically the training focused on the learning of novel motor sequences and their rhythm and timing, engaging executive processes known to be impaired in HD [45]. These included focused attention (e.g., paying attention to the drumming sequence), multi-tasking attention (e.g., listening and movement execution), movement-switching (e.g., switching between dominant and non-dominant hand) and response speed [25]. At the anatomical level, attention and executive functions rely on cortico-basal ganglia loops involving the striatum, which also plays a fundamental role in motor control and motor learning [46, 47]. This shared reliance on overlapping cortico-basal ganglia networks may contribute to the beneficial effects of physical exercise on executive functioning in healthy older adults [48], and in patients with Parkinson’s disease [49]. Additionally, in a previous pilot study assessing the feasibility of the present drumming training, we observed WM changes and improvements in executive functions in HD patients following the intervention [25].
Previous WM plasticity neuroimaging studies [27, 50] have predominantly employed indices from diffusion tensor MRI (DT-MRI) [51]. However, while sensitive, such measures are not specific to changes in specific sub-compartments of WM microstructure, challenging the interpretation of any observed change in DT-MRI indices [52, 53]. To improve compartmental specificity beyond DT-MRI, the present study explored changes in the macromolecular proton fraction (MPF) from quantitative magnetization transfer (qMT) [54] imaging and the restricted diffusion signal fraction (Fr) from the composite hindered and restricted model of diffusion (CHARMED) [55]. Fractional anisotropy (FA) and radial diffusivity (RD) from DT-MRI [51] were included for comparability with previous training studies [27, 57].
The MPF has been proposed as a proxy MRI marker of myelin [58]. Accordingly, histology studies show that this measure reflects demyelination accurately in Shiverer mice [59], is sensitive to de-myelination in multiple sclerosis patients [60] and reflects WM myelin content in post-mortem studies of multiple sclerosis brains [61]. Fr, on the other hand, represents the fraction of signal-attenuation that can be attributed to restricted diffusion, which is presumed to be predominantly intra-axonal, and therefore provides a proxy measure of axonal density [62].
Training effects were investigated in WM pathways linking the putamen and the supplementary motor area (SMA-Putamen), and within three segments of the corpus callosum (CCI, CCII, CCIII). The SMA has efferent and afferent projections to the primary motor cortex and is involved in movement execution, and previous work has reported altered DT-MRI metrics in the putamen-motor tracts of symptomatic HD patients [63]. The anterior and anterior-mid sections of the corpus callosum contain fibres connecting the motor, premotor and supplementary motor areas in each hemisphere [64]. Previous work has demonstrated a thinning of the corpus callosum in post-mortem HD brains [65], altered callosal DT-MRI metrics in both pre-symptomatic and symptomatic HD patients [66, 67], and a correlation between these metrics and performance on motor function tests [68]. Given previous reports of an effect of motor learning on myelin plasticity [38], we expected changes following training to be more marked in MPF, as compared to the other non-myelin sensitive metrics assessed in this study. We also investigated the relationship between training-associated changes in MRI measures, and changes in drumming performance and cognitive/executive function. Finally, as previous evidence has shown widespread reductions in MPF in premanifest and manifest HD patients [69], we used tract-based spatial statistics (TBSS) [70] to investigate patient-control differences in MPF before training, across the whole brain; this aided the interpretation of the post-training microstructure changes we detected.
MATERIALS AND METHODS
Participants
The study was approved by the local National Health Service (NHS) Research Ethics Committee (Wales REC 1 13/WA/0326) and all participants provided written informed consent. All subjects were drumming novices and none had taken part in our previously-reported pilot study [25]. Fifteen HD patients were recruited from HD clinics in Cardiff and Bristol. Genetic testing confirmed the presence of the mutant huntingtin allele. Thirteen age, sex, and education-matched healthy controls were recruited from the School of Psychology community panel at Cardiff University and from patients’ spouses, carers or family members. The inclusion criteria were the following: no history of head injury, stroke, cerebral haemorrhages or any other neurological condition; eligible for MRI scanning; stable medication for at least four weeks prior to the study.
Of the recruited sample, two patients were not MRI compatible, four withdrew during the study and one patient’s MRI data had to be excluded due to excessive motion. Therefore, while drumming performance and cognitive data from 11 patients were assessed, only 8 patients had a complete MRI dataset. One control participant was excluded due to an incidental MRI finding, two participants dropped out of the study and a fourth participant was not eligible for MRI. In total, we assessed drumming and cognitive tests performance in 8 controls, while MRI data from nine controls were available for analyses. Table 1 summarizes patients demographic and background clinical characteristics. Most patients were at early disease stages, however two were more advanced, as shown by their Total Motor Score (TMS; 69 and 40, respectively) and Functional Assessment Score (FAS; 18 and 17, respectively). Table 2 summarizes demographic variables and performance in the Montreal Cognitive Assessment (MoCA) [71] and in the revised National Adult Reading Test (NART-R) (Nelson, 1991) for patients and controls. While the groups did not differ significantly in age, controls were on average slightly older, performed significantly better on the MoCA, and had a significantly higher NART-IQ than patients.
Demographics and background clinical information of the patients’ cohort
Based on TMS and FAS, most of the patients were at early disease stages, however two of them were more advanced. CAG, cytosine-adenine-guanine; TMS, Total Motor Score out of 124 (the higher the scores the more impaired the performance); FAS, Functional Assessment Score out of 25 (the higher the scores the better the performance); SD, Standard Deviation.
Demographics and general cognitive profile of patients and controls
Both groups were matched for age, sex and years of education but the patient group had a lower NART-IQ and performed less well than the control group in the MoCA. NART-IQ, verbal IQ estimate based on the National Adult Reading Test; MoCA, Montreal Cognitive Assessment score out of 30.
Training intervention: Drumming-based rhythm exercises
The rhythm exercise and drumming training previously described in [25] was applied. Participants were provided with twenty-two 15 min training sessions on CDs, a pair of Bongo drums and a drumming diary and could practise at home. They were asked to exercise for 15 min per day, 5 times per week, for 2 months and to record the date and time of each exercise in their diary. Each training session introduced a drumming pattern based on one of the following rhythms: Brazilian samba, Spanish rumba, West-African kuku and Cuban son. After a brief warm up, trainees were encouraged to drum along with the instructor, initially with each hand separately and then with both hands alternating, starting with the dominant hand first and then reversing the order of the hands. The first exercises were based on very simple, slow, and regular patterns but the level of complexity and speed increased over the training sessions.
Importantly, each individual progressed through the training adaptively at their own pace i.e., as long as they exercised for the specified time, they could repeat each session as often as they felt necessary to master it. To maintain engagement and motivation, the training incorporated pieces of music based on rhythms participants had learned and could drum along to. The researcher (JBT) supervised the first training sessions and then remained in regular telephone contact (at least once a week) with each participant throughout the intervention. Whenever possible, carers and/or spouses were involved in the study to support the training. Control participants started with Session 3 since the first two exercises were built on a very low level of complexity, with slow, regular patterns of movement required. Patients, on the other hand, started with simpler exercises, but could progress to the following sessions whenever they felt comfortable.
Drumming assessment
Progress in drumming ability was assessed by digitally recording participants’ drumming performance for three patterns of ascending levels of difficulty (easy, medium and hard), which were not part of the training sessions, at baseline and after the training. Each recording was judged by an independent rater, blind to group and time, according to an adopted version of the Trinity College London marking criteria for percussion (2016) (http://www.trinitycollege.com).
Cognitive assessments
Different aspects of cognition and executive function were assessed before and after the training as previously described [25]. Multi-tasking was assessed with a dual task requiring simultaneous box crossing and digit sequences repetition [73]. Attention switching was assessed with the trails test (VT) requiring the verbal generation of letter and digit sequences in alternate order relative to a baseline condition of generating letter or digit sequences only [69]. Distractor suppression was tested with the Stroop task involving the naming of incongruent ink colours of colour words. Verbal and category fluency were tested using the letter cues “F”, “A”, “S” and “M”, “C”, “R” as well as the categories of “animals” and “boys’ names” and “supermarket items” and “girls’ names” respectively [74]. In total, we assessed 7 outcome variables, and percentage change scores in performance were computed for each of these variables (Table 3).
Cognitive outcome variables assessed in this study
Tests were carried out before and after the training, and a percentage change score was computed for each variable.
MRI data acquisition
MRI data were acquired on a 3 Tesla General Electric HDx MRI system (GE Medical Systems, Milwaukee) using an eight channel receive-only head RF coil at the Cardiff University Brain Research Imaging Centre (CUBRIC). The MRI protocol comprised the following images sequences: a high-resolution fast spoiled gradient echo (FSPGR) T1-weighted (T1-w) sequence for registration; a diffusion-weighted spin-echo echo-planar sequence (SE∖EPI) with 60 uniformly distributed directions (b = 1200 s/mm2), according to an optimized gradient vector scheme [75]; a CHARMED acquisition with 45 gradient orientations distributed on 8 shells (maximum b-value = 8700 s/mm2) [55]; and a 3D MT-weighted fast spoiled gradient recalled-echo (FSPGR) sequence [76]. The acquisition parameters of all scan sequences are reported in Table 4. Diffusion data acquisition was peripherally gated to the cardiac cycle. The off-resonance irradiation frequencies (⊝) and their corresponding saturation pulse amplitude (ΔSAT) for the 11 Magnetization transfer (MT) weighted images were optimized using Cramer-Rao lower bound optimization [76].
Scan parameters. All sequences were acquired at 3T. For each of the sequences, the main acquisition parameters are provided
T1-w, T1-weighted; MT-w, MT-weighted; FSPGR, fast spoiled gradient echo; SE, spin-echo; EPI, echo-planar imaging; SPGR, spoiled gradient recalled-echo; FoV, field of view; TE, echo time; TR, repetition time.
MRI data processing
The diffusion-weighted data were corrected for distortions induced by the diffusion-weighted gradients, artefacts due to head motion and EPI-induced geometrical distortions by registering each image volume to the T1-w anatomical images [77], with appropriate reorientation of the encoding vectors [78], all done in ExploreDTI (Version 4.8.3) [79]. A two-compartment model was fitted to derive maps of FA and RD in each voxel [80]. CHARMED data were corrected for motion and distortion artefacts according to the extrapolation method of [81] The number of distinct fiber populations (1, 2, or 3) in each voxel was obtained using a model selection approach [52], and Fr was calculated per voxel with an in-house software [52] coded in MATLAB (The MathWorks, Natick, MA). MT-weighted SPGR volumes for each participant were co-registered to the MT-volume with the most contrast using an affine (12 degrees of freedom, mutual information) registration to correct for inter-scan motion using Elastix [82]. The 11 MT-weighted SPGR images and T1 map were modelled by the two pool Ramani’s pulsed MT approximation [83, 84], which included corrections for amplitude of B0 field inhomogeneities. This approximation provided MPF maps, which were nonlinearly warped to the T1-w images using the MT-volume with the most contrast as a reference using Elastix (normalized mutual information cost function) [82].
Deterministic tractography
Training-related changes in FA, RD, Fr, and MPF were quantified using a tractography approach in pathways interconnecting the putamen and the supplementary motor area bilaterally (SMA-Putamen), and within three segments of the corpus callosum (CCI, CCII, and CCIII) [64] (Fig. 1).
Pathway regions of interest. Sagittal views of the reconstructed WM pathways displayed on a T1-weighted image for one control participant. (A) CCI, CCII, and CCIII [64]: CCI is the most anterior portion of the CC and maintains prefrontal connections between both hemispheres; CCII is the portion that maintains connections between premotor and supplementary motor areas of both hemispheres. CCIII maintains connections between primary motor cortices of both hemispheres. (B) SMA-putamen pathway: this pathway has efferent and afferent projections to the primary motor cortex and is involved in movement execution.
Whole brain tractography was performed for each participant in their native space using the damped Richardson-Lucy algorithm [85], which allows the recovery of multiple fiber orientations within each voxel including those affected by partial volume. The tracking algorithm estimated peaks in the fiber orientation density function (fODF) by selecting seed points at the vertices of a 2×2×2 mm grid superimposed over the image and propagated in 0.5 mm steps along these axes re-estimating the fODF peaks at each new location [86]. Tracks were terminated if the fODF threshold fell below 0.05 or the direction of pathways changed through an angle greater than 45° between successive 0.5 mm steps. This procedure was then repeated by tracking in the opposite direction from the initial seed-points.
The WM tracts of interest were extracted from the whole-brain tractograms by applying way-point regions of interest (ROI) [87]. These were drawn manually by one operator (JBT) blind to the identity of each dataset on color-coded fiber orientation maps in native space guided by the following anatomical landmark protocols (Fig. 2).
Sagittal views of the tractography protocols. (A) CCI, CCII and CCIII; (B) SMA - putamen pathway. Booleian logic OR waypoint regions of interest gates are illustrated in blue; AND gates in green. M, Midline.
Corpus callosum
Reconstruction of the CC segments followed the protocol of Hofer and Frahm [64] as illustrated in Fig. 2A. Segment reconstructions were visually inspected and, if necessary, additional gates were placed to exclude streamlines inconsistent with the known anatomy of the CC.
SMA-putamen pathway
One axial way-point ROI was placed around the putamen and one axial ROI around the supplementary motor cortex [88] (Fig. 2B). A way-point gate to exclude fibres projecting to the brain stem was placed inferior to the putamen.
Statistical analyses
Statistical analyses were carried out in R Statistical Software (Foundation for Statistical Computing, Vienna, Austria).
Assessment of training effects on drumming performance
Improvements in drumming performance were analysed with a two-way mixed analysis of variance (ANOVA) testing for the effects of group (HD/controls), time of assessment (before/after the training) and group by time interaction effects. We also confirmed detected effects to the ones obtained by running a robust mixed ANOVA, using the bwtrim R function from the WRS2 package [89]. This implements robust methods for statistical estimation and therefore provides a good option to deal with data presenting small sample sizes, skewed distributions and outliers [90]. Significant effects were further explored with post-hoc paired and independent t-tests. The reliability of the post-hoc analyses was assessed with bootstrap analysis based on 1000 samples and the 95% confidence interval (CI) of the mean difference is provided for each significant comparison.
Assessment of group differences in the effect of training on cognitive performance
Performance measures in executive function tasks have been shown to share underlying cognitive structures [91]. Therefore, PCA was employed to reduce the complexity of the cognitive data and hence the problem of multiple comparisons as well as to increase experimental power. PCA was run on change scores for all participants across both groups. Due to the relatively small sample size, we first confirmed with the Kaiser-Meyer-Olkin (KMO) test that our data was suited for PCA. Subsequently, we followed guidelines to limit the number of extracted components [92, 93], as follows: first, we employed the Kaiser criterion of including all components with an eigenvalue greater than 1; second, we inspected the Cattell scree plot [94] to identify the minimal number of components that accounted for most variability in the data; third, we assessed each component’s interpretability. A PCA procedure with orthogonal Varimax rotation of the component matrix was used. Loadings that exceeded a value of 0.5 were considered as significant.
Next, we assessed group differences in the component scores with permutation analyses, to understand whether the training had differentially affected HD patients as compared to controls. Permutation testing relies only on minimal assumptions and can therefore be applied when the assumptions of a parametric approach are untenable such as in the case of small sample sizes. Significant group differences were tested using 5,000 permutations and the effect sizes of significant differences were assessed with Cohen’s d [95]. Multiple comparison correction was based on a 5% false discovery rate (FDR) using the Benjamini-Hochberg procedure [96].
Training effects on WM microstructure
Median measures of FA, RD, Fr and MPF were derived for each of the reconstructed tracts in ExploreDTI [79]. A percentage change score in these measures between baseline and post-training was calculated in each tract (CCI, CCII, CCIII, left and right SMA-Putamen).
Previous research has shown that variation in the microstructural properties of WM may represent a global effect, rather than being specific to individual tracts, and that WM measures are highly correlated across WM areas [56, 98]. Therefore, we inspected the inter-tract correlation for each of microstructural metric and found that MPF values were highly correlated, whereas this was not true for the other metrics (Fig. 3). Hence, percentage change scores in MPF across the different tracts were transformed with PCA to extract meaningful anatomical properties, following the procedure described above for the PCA of cognitive change scores. PC scores for each participant were then used as dependent variables in a permutation-based analysis using 5,000 permutations to assess group differences in training associated changes in MPF. Finally, as a post-hoc exploration, we looked for between-groups differences in MPF changes in the individual tracts using 5000 permutations.
Correlation matrices for the MRI metrics investigated across the different WM pathways. Colour intensity and the size of the circles are proportional to the strength of the correlation. *p < 0.05, **p < 0.01, ***p < 0.001. The absolute correlation coefficient is plotted. MPF values were highly correlated across tracts, whereas this was not true for the other metrics.
Training-associated changes in FA, Fr and RD were investigated with permutation analyses separately for each tract. Significant group differences in these measures were tested using 5,000 permutations. Multiple comparison correction was based on a 5% FDR using the Benjamini-Hochberg procedure [96]. Cohen’s d [95] was used to assess the effect size for those changes found to be significant.
Training effects on WM microstructure
TBSS [70] was carried out to investigate baseline differences in MPF between HD subjects and healthy controls, to gain a better insight into differences in training-associated changes. To produce significance maps, a voxel-wise analysis was performed on the MPF projected 4D data for all voxels with FA≥0.20 to exclude peripheral tracts where significant inter-subject variability exists. Inference based on permutations (5,000 permutations) and threshold-free-cluster-enhancement was used. The significance level was set at p < 0.05 and corrected by multiple comparisons (family-wise error, FWE).
Relationship between changes in MRI measures and changes in drumming and cognitive performance
We computed percentage change scores for the drumming performance, in the same way cognitive change scores were calculated. Scores were computed for the easy test pattern in patients and for the medium test pattern in controls, as these training patterns showed a significant improvement in the two groups, respectively. Spearman correlation coefficients were calculated between drumming and cognitive performance, and microstructural components that showed significant group differences, to assess whether microstructural changes were related to any drumming and/or cognitive benefits of the training.
Exploration of the possible confounding effects of differences in training-compliance and IQ
We examined the training diaries of each participant to assess compliance with training. Each session was marked as completed if the whole 15-min session had been carried out. Each participant was assigned a score representing the number of training sessions they performed (e.g., a score of 40 if 40 sessions had been carried out). We then assessed group differences with permutation analyses, to test whether there was a significant difference in the amount of training sessions carried out by the groups, and hence understand whether this variable had to be accounted for in the analysis. Significant group differences were tested using 5,000 permutations.
Finally, as there was a significant difference in IQ between patients and controls (Table 2), we investigated whether any training-associated change might be due to differences in premorbid intelligence. The sample size of this experiment was very small, and therefore it was not possible to perform a multiple regression or an analysis of covariance (ANCOVA), to understand the possible influence of IQ as confounding variable. Accordingly, the statistical power to establish the incremental validity of a covariate in explaining an outcome has been shown to be extremely low, and therefore to require large sample sizes [99]. As a potential solution to this issue, we instead performed separate non-parametric Spearman correlation analyses between NART-IQ scores (as this test showed the largest difference between groups), and MRI and cognitive measures showing significant training effects; this allowed us to gain some insight into whether there was a significant association between premorbid intelligence and training-associated changes.
RESULTS
Training effects on drumming performance
The mixed ANOVA of drumming performance for the easy and medium test pattern showed a significant effect of group [easy: F(1,17) = 22.3, p < 0.001; medium: F(1,17) = 13.1, p = 0.002] and time [easy: F(1,17) = 12.83, p = 0.004; medium: F(1,17) = 13.4, p = 0.002] but no interaction (easy: p = 0.8; medium: p = 0.3). For the hard test pattern there was only a significant effect of group [F (1,17) = 9.95, p = 0.006] but not of time (p = 0.1) and there was no interaction (p = 0.4). Results from the robust mixed ANOVA were largely consistent with the above. Specifically, the easy and medium test patterns showed a significant effect of group (easy: p = 0.002; medium: p = 0.02) and time (easy: p = 0.04; medium: p = 0.049) but no interaction (easy: p = 0.45; medium: p = 0.69). The hard test pattern showed a significant effect of group (p = 0.02) but not of time (p = 0.22) and no interaction (p = 0.8). Figure 4 summarises the average drumming performance per group and time point. Overall patients’ drumming performance was poorer than controls. Patients improved their drumming performance significantly for the easy pattern [t (10) = 2.7, p = 0.02; 95% CI of mean difference: 1.5–7.8] and controls for the medium pattern [t (7) = 3.8, p = 0.01; 95% CI of mean difference: 2.8–8.5].
Mean ratings for drumming performance according to the Trinity College London marking criteria for percussion (2016) as a function of group and time point. Patients improved their drumming performance significantly for the easy test pattern and controls for the medium difficult test pattern. *p < 0.05, **p < 0.01, bootstrapping based on 1000 samples.
Group differences in the effect of training on cognitive performance
Three components, accounting for 79% of the variance in performance improvement in the cognitive benchmark tests, were extracted. The first component loaded highly on performance changes in the dual task (total number of boxes identified under dual task condition), the Stroop task (Stroop interference score), and the trails making task (Trail test switching). Since these variables all measure executive functions including focused attention and distractor suppression, the first component was labelled “executive” component. The second component loaded on variables reflecting the ability to correctly recall digits sequences (i.e., number of correct digits recalled under single and dual task condition) and was therefore labelled “working memory capacity” component. Finally, the third extracted component loaded highly on verbal and category fluency and was therefore labelled “fluency” component (Table 5).
Rotated Component Loadings on Change in the Cognitive Benchmark Tests
Significant loadings (>0.5) are highlighted in bold.
We tested whether the two groups differed in terms of post-training cognition changes, by running permutation analyses on the individual scores for the three extracted components. The two groups differed in the executive component (t = –1.03, p = 0.008, FDR-corrected p = 0.024, d = 1.15). However, no significant group differences were detected in the other two components (Working Memory capacity: t = –0.22, p = 0.3296, FDR-corrected p = 0.3296; Fluency: t = –0.39, p = 0.242 FDR corrected p =0.3296).
Training effects on WM microstructure
Table 6 reports a summary of the training associated changes in FA, RD, Fr and MPF, across the different tracts.
Summary statistics for the permutation analysis of training effects on FA, RD, Fr and MPF, across the investigated tracts
Training-associated group differences in FA
Permutation analyses of FA changes across the different tracts revealed no significant differences between HD and control groups [CCI: t = 1.22, p = 0.91 (FDR-corrected); CCII: t = 2.65, p = 0.91 (FDR-corrected); CCIII: t = 0.325, p = 0.13 (FDR-corrected); right SMA-Putamen: t = –9.54, p = 0.10 (FDR-corrected); left SMA-Putamen: t = 5.16, p = 0.77 (FDR-corrected).
Training-associated group differences in RD
There were no significant differences in RD changes following training between HD patients and controls [CCI: t = –0.48, p = 0.45 (FDR-corrected); CCII: t = –1.29, p = 0.45 (FDM-corrected); CCIII: t = –1.04, p = 0.45 (FDR-corrected); right SMA-Putamen, t = 4.01, p = 0.81 (FDR-corrected); left SMA-Putamen, t = –3.68, p = 0.39 (FDR-corrected).
Training-associated group differences in Fr
Permutation analyses of Fr changes across the different tracts revealed no significant differences between HD and control groups [CCI: t = 3.39, p = 0.82 (FDR-corrected; CCII: t = –0.17, p = 0.82 FDR-corrected; CCIII: t = 3.08, p = 0.82 (FDR-corrected); right SMA-Putamen: t = –5.24, p = 0.82 (FDR-corrected); left SMA-Putamen: t = 1.05, p = 0.82 (FDR-corrected)].
Training-associated group differences in MPF
PCA of change scores in MPF revealed one single component explaining 70.2% of the variance. This component presented high loadings from all the tracts investigated. A significant group difference was found for the MPF change-score component, indicating that HD patients presented significantly greater MPF changes in response to training, as compared to controls [t (14) = –1.743, n = 17, p = 0.03, d = 1.796].
Finally, we found a significant difference in mean MPF change scores between the two groups for CCII [t (14) = –20.72, p = 0.04, d = 0.93], CCIII [t (14) = –25.87, p = 0.04, d = 1.07], and the right SMA-putamen pathway [t (14) = –25.48, p = 0.04, d = 1.15] after FDR correction, therefore indicating that there was a differential group effect of training on MPF within these tracts (Fig. 5).
MPF changes scores: PCA scree plot (A); plot summarising how each variable is accounted for in every principal component - colour intensity and the size of the circles are proportional to the loading: PC1 loads on CCI, CCII, CCIII and right SMA-Putamen, while PC2 loads mostly on the left SMA-Putamen; the absolute correlation coefficient is plotted (B); correlation circle, interpreted as follows: 1) positively correlated variables are grouped together, 2) negatively correlated variables are positioned on opposite sides of the plot origin (opposite quadrants), 3) the distance between variables and the origin measures the quality of the variable on the factor map. Variables that are away from the origin are well represented on the factor map (C); Bar graph of the percentage change in MPF across the inspected tracts; Error bars represent the standard error; training was associated with a significantly greater change in MPF in CCII, CCIII, and right SMA-Putamen; *p < 0.05, results corrected for multiple comparisons with FDR (D).
Relationship between training-associated changes in MRI measures, and changes in drumming and cognitive performance.
We did not find a significant association between the ‘MPF’ component scores and improvement in drumming performance (PC1: rho = –0.14, p > 0.05). Moreover, although no significant correlation was observed between the ‘Executive’ and the ‘MPF’ component scores, there was a positive trend (rho = 0.348, p = 0.171).
Baseline differences in MPF
Baseline MPF was reduced in the HD group compared to controls, in the midbody of the CC (t = 3.13, p = 0.05, FWE corrected). Figure 6 shows the areas with reduced MPF in HD patients, in blue.
TBSS analysis of baseline MPF values (A). Light blue areas show a significant reduction of MPF in patients with HD compared to controls (p < 0.05, FWE corrected). The midbody of the CC was mostly found to be affected, which carries connections to the premotor, supplementary motor and motor areas of the brain. Tracts showing significantly greater MPF changes in HD patients post-training as compared to controls (B). Areas showing significant MPF reductions at baseline overlap with tracts showing significant changes post-training (i.e., CCII and CCIII).
Exploration of the possible confounding effects of differences in training-compliance and IQ
The permutation analysis of training compliance revealed that there was not a significant difference in the number of training sessions between HD patients (mean = 38.1, SD = 4.2) and healthy controls (mean = 39.6, SD = 1.2), t = 1.45, p = 0.76. Therefore, the total time spent training should not have influenced any training-associated change in microstructure and/or cognitive measures. Furthermore, we did not detect a significant association of participants’ NART-IQ scores with ‘MPF’ component scores (rho = –0.4, p = 0.10), nor with ‘Executive’ component scores (rho = –0.34, p = 0.18), suggesting that premorbid intelligence should not have had an influence on training-associated effects.
DISCUSSION
Based on evidence that myelin impairment contributes to WM damage in HD [3], and the suggestion that myelin plasticity underlies the learning of new motor skills [27, 38], the present study explored whether two months of drumming training would result in changes in WM microstructure in HD patients. Specifically, we expected to detect changes in MPF, as marker of WM myelin plasticity, in HD patients relative to healthy controls.
Firstly, we demonstrated a behavioural effect of the training by showing a significant improvement in drumming performance in patients (easy test pattern) and controls (medium test pattern). We did not detect any group differences in training-associated changes in the diffusion-based indices of FA, RD and Fr. However, as hypothesised, we found a group difference in training-induced changes in the MPF PCA component. Specifically, HD patients showed significantly higher increases in MPF relative to controls. Furthermore, through exploratory post-hoc investigations, we detected significantly higher training-induced MPF changes within the CCII, CCIII and the right SMA-putamen pathway between patients and controls. Additionally, TBSS analysis of baseline differences in MPF suggested partial overlap of WM areas showing significant MPF reductions at baseline with areas showing changes post-training (i.e., CCII and CCIII).
MPF can be affected by inflammation [100] and in advanced HD it is likely that inflammation goes hand-in-hand with myelin breakdown [101]. However, a recent CSF biomarker study found no evidence of neuro-inflammation in early-manifest HD [102]. Furthermore, recent evidence shows that this measure may be inconsistent when investigated in relatively small WM areas, presumably because of the effect of spatial heterogeneity in myelin thickness [103]. Nevertheless, the within-subjects design employed in the present study should have helped to minimise noise due to the spatial inconsistency of this measure. Therefore, though preliminary and based on a small sample size, these findings suggest that two months of drumming and rhythm exercises may result in myelin remodelling in patients with early HD.
It is plausible that this group difference arose due to WM microstructural differences between patients and controls before the training. Accordingly, the HD group showed a significantly lower baseline MPF, consistent with lower myelin content [3]. Furthermore, previous studies have reported that training-associated percentage changes in MRI measures tend to be higher in patients than in healthy subjects [34]. One possibility is that in the healthy brain, neural networks may be optimally myelinated, and further increasing myelin may not improve performance [104–106]. Hence, the MPF changes in patients relative to controls might depict mainly a ‘catch-up effect’ to the better baseline status of the control group. However, disentangling the impact of prior WM microstructural differences on microstructural plasticity during learning is beyond the scope of the current work.
Notably, the behavioral effect of drumming training and cognition differed between patients and controls. Patients improved in the easy drumming test pattern, and controls improved in the medium test pattern. Furthermore, consistent with evidence from our pilot study [25], patients showed increases in the executive function components whilst control participants did not show improvements in their cognition. Therefore, inter-group differences in microstructural changes might not only be due to baseline WM microstructural differences, but also to a different behavioral effect of the task between HD subjects and controls. For instance, control participants performed close to ceiling in the easy test pattern, and as the training was tailored to patients’ needs, some of the earlier practice sessions may not have optimally challenged them. The fact that the training seemed more taxing for patients than controls may also explain why improvements in executive functions and changes in MPF were only observed for the patients. Interestingly, baseline IQ was significantly different between the two groups, and this might have had some influence on training performance and on training-associated effects. Here, we failed to find an association between IQ scores and changes in MRI and cognitive measures. However, future studies with larger sample sizes might allow to utilize more advanced statistical approaches, such as an ANCOVA, to better model the possible influence of premorbid IQ on training effects, both at the neural and cognitive level.
A critical question relevant to all training studies concerns the functional significance of any observed neural changes. If, and to what degree adaptive alterations in myelin content can facilitate behavioural change remains poorly understood [105]. In the present study, no significant relationships between changes in MRI measures and changes in drumming proficiency or performance in cognitive tests were found. This might have been due to non-specific training-related neural responses. Specifically, while the training exercise might have triggered changes in brain structure, training-induced changes may not necessarily co-vary with improvements in performance. Alternatively, it might be that the study was insufficiently powered to detect brain-function correlations. The minimum sample size required to detect a correlation was calculated to be 64 people (α= 0.05; 80% power; medium effect size; GPower 3 software). Therefore, these results need replication in larger samples. A lack of correlation between structural and functional changes after training has been reported in other studies (including well-powered studies) and may suggest that these processes follow different time courses and/or may occur in different brain regions [107].
It is important to note that our study did not include a non-intervention patient group. Within the 12-month time period of this study it was not possible to recruit a sufficiently large number of well-matched patient controls. Therefore, we cannot disentangle the effects of the training on WM microstructure from HD-associated pathological changes. However, given that HD is a progressive neurodegenerative disease associated with demyelination [3], it is unlikely that increases in MPF observed in the patient group were due to the disease itself. Finally, while the majority of training studies assess brain structural changes between baseline and post-training [34], presumably on account of cost and participant compliance, we suggest that acquiring intermittent scans during the training period could have helped to better capture and understand changes in WM microstructure observed in this study. Accordingly, future studies might be able to provide greater insights into the complex nonlinear relationships between structural changes and behaviour [108].
To conclude, we have demonstrated that two months of drumming and rhythm exercises result in a significantly greater change in a proxy MRI measure of myelin in patients with HD relative to healthy controls. Whilst the current results require replication in a larger patient group with an appropriately matched patient control group, they suggest that behavioural stimulation may result in neural benefits in HD that could be exploited for future therapeutics aiming to delay disease progression.
CONFLICT OF INTEREST
The authors have no conflict of interest to report.
Footnotes
ACKNOWLEDGMENTS
The present research was funded by a Wellcome Trust Institutional Strategic Support Fund Award (ref: 506408) to CMB, AR, and DKJ and a Wellcome Trust PhD studentship to CC (ref: 204005/Z/16/Z); DKJ is supported by a New Investigator Award from the Wellcome Trust (ref: 096646/Z/11/Z). We would like to thank Candace Ferman and Louise Gethin for collating the patients’ clinical details and Jilu Mole for assistance with data analysis. This manuscript has been released as a pre-print at bioRxiv: 
.