Introduction

Huntington’s disease (HD) is an autosomal, dominantly inherited, slowly progressive, neurodegenerative disease localised on chromosome four. The pathophysiological mechanism is complex and although impaired energy metabolism and neuronal excitotoxicity are named as major contributors [9, 32], it is still not fully understood, but the main result is neuronal dysfunction and loss in the brain. The disease is characterised by clinical symptoms in the motor, cognitive, and psychiatric domains.

Loss of brain cells and their connections, referred to as atrophy, can be visualised with magnetic resonance imaging (MRI) techniques and has been the focus of many studies in search of objective biomarkers for tracking disease progression. The striatum has been the primary region of interest, due to previous findings in pathological studies that demonstrated the most severe loss of neurons in this region [28]. In manifest HD, the evidence for widespread atrophy is extensive, whereby several MRI studies showed striatal [4, 6, 8, 11, 13, 16, 20, 26, 34] and cortical atrophy [29, 31]. The striatal atrophy occurs in the premanifest stages of HD, with atrophy of the caudate nucleus and putamen present up to a decade or more before clinical manifestation [5, 7, 14, 17, 25, 40].

Other subcortical grey matter structures have been examined to a lesser degree. Thalamic [6, 8, 13, 15] and pallidum [2, 8, 30, 34] atrophy is apparent in the manifest stages of HD; however, findings are controversial in the premanifest stage, as not all studies are in agreement [7, 10, 14, 17, 30, 40]. Atrophy of the hippocampus, amygdala, accumbens nucleus, and brainstem has only been reported in the manifest stage in a very limited number of studies [8, 30]. Whole brain atrophy has been demonstrated in HD [12, 39]. However, the rate of volume reduction of the subcortical structures as compared to whole brain atrophy has not been reported in any of the above mentioned studies and is of importance as this regional atrophy may very well only be a reflection of whole brain atrophy. Correlation studies have been performed for volume reductions and global clinical measures in manifest HD [8, 20, 30]. However, reports on the relationship of atrophy with clinical measures of all these structures in different disease stages in one large cohort has not been previously reported.

With the possibilities of new therapies for HD, there has been increasing interest in the development of robust, reproducible biomarkers to assess the efficacy of potential treatments. The TRACK-HD study is an international multi-centre study investigating candidate biomarkers for HD and encompasses a large, well defined, sample of HD gene carriers [39]. Application of a fully automated MRI-analysis tool to this study sample provides the opportunity to study brain atrophy in detail within a very large HD cohort. The aim of this study is to provide a more precise overview of the specific atrophy of the subcortical grey matter structures in different disease stages, whilst taking into account whole brain atrophy and provide better targets for biomarker research. Secondly, we aim to examine the relationship of all subcortical grey matter structure volumes to clinical measures for HD.

Methods

Subjects

The TRACK-HD study recruited 366 participants from four centres (The National Hospital for Neurology and Neurosurgery in London, The Department of Medical Genetics at University of British Columbia in Vancouver, The Department of Genetics and Cytogenetics at the Hôpital de la Salpetrière-Université Pierre and Marie Curie in Paris and the Department of Neurology at Leiden University Medical Center in Leiden), totalling 123 early manifest HD, 120 premanifest HD and 123 control subjects. Inclusion criteria for the early manifest HD group consisted of genetic confirmation of an expanded CAG repeat of ≥40, presence of motor disturbances on the Unified Huntington’s Disease Rating Scale (UHDRS), defined as a diagnostic motor confidence score of 4 and a total motor score (TMS) of >5, and a total functional capacity (TFC) of ≥7. Inclusion criteria for the premanifest group consisted of genetic confirmation of an expanded CAG repeat of ≥40, absence of motor disturbances on the UHDRS, defined as a TMS of ≤5. Furthermore, premanifest gene carriers required a burden of pathology score >250 [27], based on CAG length and age, to ensure a premanifest HD sample close to onset. Spouses or partners of premanifest and affected subjects or gene-negative siblings were recruited as healthy control subjects. The control group was age- and gender-matched to the combined premanifest and manifest HD group. Further subdivision of the premanifest group was performed on the basis of predicted years to diagnosis into preHD-A (10.8 or more years to disease onset) and preHD-B (closer than 10.8 years to disease onset) based on Langbehn’s et al. survival analysis formula [18, 19]. The early manifest HD subjects were also divided into two groups based on disease stage by means of TFC score, the HD1 group (TFC score 11–13) and the HD2 group (TFC score 7–10) [35]. The study was approved by the local ethical committees and written informed consent was obtained from each subject. For full details on study design see [39].

Clinical measures

From the total TRACK-HD assessment battery, the UHDRS motor scale and TFC were used for this study. The UHDRS TMS is a representation of the severity of motor disturbances, scores potentially range from 0 to 124, with higher scores indicating more severe motor abnormalities. The TFC is a scale used for assessment of global impairments within daily activities, ranging from 0 to 13, with lower scores indicating impaired global functioning.

MRI acquisition

All participants underwent 3 T MRI scanning. T1-weighted image volumes were acquired using a 3D MPRAGE acquisition sequence on a 3 T Siemens or a 3 T Phillips whole body scanner with the following imaging parameters: TR = 2200 ms (Siemens)/7.7 ms (Philips), TE = 2.2 ms (Siemens)/3.5 ms (Philips), FA = 10° (Siemens)/8° (Philips), FOV = 28 cm (Siemens)/24 cm (Philips), matrix size 256 × 256 (Siemens)/224 × 224 (Philips), 208 (Siemens)/164 (Philips), sagittal slices to cover the entire brain with a slice thickness of 1.0 mm with no gap between slices.

Post-processing

All T1-weighed scans were analysed using software provided by FMRIB’s software library (FSL) [37]. Eight subcortical regions were assessed using FMRIB’s integrated registration and segmentation tool (FIRST) [2123], namely the accumbens nucleus, amygdala, brainstem, caudate nucleus, hippocampus, pallidum, putamen, and the thalamus. Total brain tissue volume was estimated with SIENAX [36, 38]. For details on MRI post-processing, please see the Supplementary data. The whole-brain volume was used to calculate ratios of atrophy of the subcortical structures to whole-brain atrophy. The obtained ratio was divided by the control group value in order to obtain an easily interpretable ratio whereby 1.00 is the control (normal) value ratio. The following formula was used:

$$ [{\text{volume structure}}/\left( {{\text{whole-brain volume}} - {\text{volume structure}}} \right)]/{\text{control group value}} $$

A value below 1.00 demonstrates that the structure displays more atrophy than the overall whole-brain atrophy. A value greater than 1.00 demonstrates that the overall whole-brain atrophy is greater than the atrophy of the structure at hand.

Statistics

Analysis of variance for group comparison was performed for volumes and ratios of the eight different subcortical structures, whilst correcting for intracranial volume, study site, age and gender. (Note that each site used a different MRI scanner model, thus, potential scanner effects can not be differentiated from site effects.) Comparisons of each structure were made for each disease group as compared to the control group, and as compared to the previous disease stage group, e.g. the preHD-B group was, besides comparison to the control group, also compared to the preHD-A group. We calculated partial correlations between the various subcortical structures and total motor score and total functional capacity from the UHDRS battery. These correlations used volumes corrected for intracranial volume and were further controlled for age, site, and gender by “partialling” these covariates out during the correlation calculations. The UHDRS TMS and TFC analyses were performed on the basis of a two-way group division, namely the premanifest gene carriers and the manifest HD group in order to keep sufficient power to perform this analysis. The significance level was set at p value <0.05.

Results

Characteristics

The clinical characteristics of the five study groups are shown (Table 1). The progressive nature of HD inherently leads to a higher age of the manifest HD groups. Except for the HD2 group, which displayed lower levels of education, there were no significant differences in education level between the disease groups and the controls. For full details on group characteristics we refer to the recent baseline TRACK-HD paper [39].

Table 1 Group characteristics

Volumes

Examples of a typical output from the segmentation by the FIRST and SIENAX tools are depicted (Fig. 1), for the FIRST examples only the left hemisphere segmentations are shown; however, both sides were segmented for the analysis.

Fig. 1
figure 1

a Example segmentations by FSL FIRST. Dark blue Brainstem, Blue Pallidum, Light blue Thalamus, Green Caudate nucleus, Dark red Accumbens nucleus, Red Putamen, Copper Hippocampus, Pink Amygdala. b Example output from the FSL SIENAX tool. Red color depicts the final sienax segmentation results

Volume comparisons for eight subcortical structures in five groups are shown (Table 2). Furthermore, Fig. 2 shows the amount of volume reduction in terms of percentage of the control value. The putamen and the caudate nucleus show a similar pattern of a smaller volume in the preHD-A group and in every subsequent disease stage group with a continuous decline in volume. The accumbens nucleus and pallidum show a decreased volume in the preHD-B group and from every following disease stage group as compared to the controls as well as to the previous disease stage group except for the accumbens nucleus between HD1 and HD2.

Table 2 Volumes
Fig. 2
figure 2

Percentage of volume reduction. Percentage of volume for eight structures as compared to the control value

The brainstem and thalamus show signs of reduced volume in preHD-B, HD1 and HD2 as compared to the control group.

The amygdala showed no reduced volume as compared to the control group. The hippocampus showed volume reduction in stages preHD-B, HD1 and HD2 as compared to controls. In terms of percentages of volume reduction as compared to the control group value it is apparent that the pallidum and accumbens nucleus show marked atrophy with only 62 and 60% of the normal volume remaining in the HD2 group, respectively. These structures are closely followed by the putamen and the caudate nucleus with 67 and 70% remaining volume, respectively. Percentages for all structures in the four disease stages are also available in the supplementary data.

Ratios

Comparisons for volume ratio of eight structures as compared to whole-brain atrophy are shown (Table 3). A smaller ratio indicates that the volume of this specific structure declines faster than the overall brain atrophy (Table 4).

Table 3 Ratios
Table 4 Volume correlations with motor score and functional capacity

The putamen atrophy ratio showed a continuous smaller value across subsequent disease stage groups. The caudate nucleus atrophy ratio showed a similar pattern, except within the manifest stage where this ratio did not show a significant difference across HD1 and HD2. The accumbens nucleus and the pallidum atrophy ratio showed a decrease over disease stages from the preHD-B group onwards.

The ratio of whole-brain atrophy to amygdala atrophy is comparable or slightly higher to the control group ratio. The atrophy ratio for the brainstem and thalamus shows no or slightly higher ratios as compared to the control group, stating that atrophy of these structures is comparable or slower to whole brain atrophy. Finally, the hippocampus atrophy ratio shows a drop in value in HD2 as compared to control and HD1.

Clinical measures

The UHDRS TMS in the premanifest stage is most strongly correlated to the atrophy of the caudate nucleus (r = 0.272, p = 0.003) and putamen (r = −0.240, p = 0.009), but also to the atrophy of accumbens nucleus, pallidum and thalamus. The TFC demonstrates a strong ceiling effect in the premanifest stage and no significant correlations are to be found. In the manifest stage of the disease, the UHDRS TMS is most strongly related to the putamen (r = −0.530, p < 0.0001), pallidum (r = −0.390, p < 0.0001) and accumbens nucleus (r < −0.380, p < 0.0001), but also to the hippocampus, thalamus, and brainstem. The TFC is most strongly correlated to atrophy of the putamen (r = 0.270, p = 0.004), hippocampus (r = 0.240, p = 0.01) and pallidum (r = 0.230, p = 0.01) and a weak correlation to accumbens nucleus and brainstem.

Discussion

The aim of this study was to provide a comprehensive overview of the atrophy of subcortical structures, with respect to disease stage while taking into account whole-brain atrophy. The main findings show that all subcortical grey matter structures at some point have a reduced absolute volume. However, the accumbens nucleus, the caudate nucleus, the pallidum and the putamen display a true progressive decline in volume reduction which is disproportionate to overall whole-brain atrophy and this is already visible in the premanifest stages of HD. Conversely, the amygdala, brainstem and thalamus do show reduced volumes, but this merely reflects the overall whole brain atrophy rate. The hippocampus is exceptional in the fact that only in the manifest stage is a true rapid decline in this structure observed. Furthermore, the volume reductions of the accumbens nucleus, putamen, pallidum and hippocampus have a severe impact on clinical measures in manifest HD.

The atrophy of the caudate nucleus and putamen in the premanifest stages confirms earlier studies. The fact that this happens a decade or more before estimated time to onset is in accordance with the findings of Aylward et al. [3], the findings of the PREDICT-HD study [25] and with the results reported in TRACK-HD baseline paper [39].

Atrophy of the pallidum and thalamus in the premanifest stages of HD is controversial. Campodonico et al. and Harris et al. report no significant reduction in volume of the pallidum [7, 10], while others do report this finding [14, 17, 40]. With regard to the thalamus in the premanifest stage, reports of reduced volume are available in favour [10, 40] and in dispute [14, 24] of thalamic involvement. Our findings, using 3 T MRI and in the largest cohort reported thus far, confirms that atrophy in the pallidum is apparent in the premanifest stage closer to predicted disease onset, but not in the further from predicted onset stage. Although thalamic atrophy is also apparent in the premanifest close to onset stage, this merely reflects the whole brain atrophy.

The brainstem, thalamus, and amygdala follow the normal atrophy rate of the whole brain in the manifest stages of the disease. In fact, in HD stage 1 the brainstem and the amygdala show signs of being relatively spared. The hippocampus shows a different pattern of atrophy as it starts to show signs of atrophy in the preHD-B stage, but really becomes significant in terms of volume loss in HD stage 2. To some degree the above described findings have been reported by Rosas et al. and Douaud et al. [8, 30]. The added value of our study lies not only in the large sample size facilitating the further specification of the premanifest and manifest groups, but particularly in the fact that the atrophy described takes into account whole-brain atrophy. It should be noted that some caution is appropriate when interpreting the ratio to whole-brain atrophy, as this method does not necessarily imply that the volume reduction of a certain structure has to be greater than the whole brain atrophy to be of significance in terms of clinical symptoms.

The clinical impact of atrophy of subcortical grey matter is difficult to assess as many brain structures are involved in complex functionality. We have demonstrated that of the subcortical grey matter structures, specifically volume reductions of the accumbens nucleus, putamen and pallidum have the greatest influence on predicting the motor disturbances in manifest HD, and perhaps surprisingly, the impact of the caudate nucleus is minimal. The role of the pallidum in motor function is underestimated and surpasses the caudate nucleus and the thalamus. However, the putamen is the most important structure for this clinical symptom. Our study shows that decreased global functioning of manifest HD was mainly related to volume loss of the hippocampus, putamen and pallidum. Cognitive and motor deficiencies have been shown in a number of separate studies, all examining only part of the deep grey matter structures, to relate to volume reductions of the basal ganglia, thalamus and frontal lobe [1, 4, 6, 7, 10, 14, 26, 33]. Within one small cohort examination of clinical measures with the TFC was carried out by Rosas et al. who provided evidence for an association between each separate structure and the TFC in the manifest stage [30], which, apart from some differences between studies, is largely confirmed by our study. All these studies show that there is an association between volumes and clinical measures; however, we hope that by examining all structures within one large, well defined cohort, we can provide evidence for which structure has the most contribution (or at least the strongest correlation) to global functioning and motor disturbances.

TRACK-HD is a longitudinal observational biomarker study aimed at providing essential methodology for the assessment of therapeutic interventions. With respect to this goal, we can add that the accumbens nucleus and pallidum show similar biomarker potential in addition to the well recognised structures putamen and caudate nucleus. The hippocampus is also of importance as it has a high correlation with clinical symptoms. The important role of the caudate nucleus in the premanifest stage of the disease seems to diminish after manifestation, as was also suggested by the TRACK-HD baseline paper [39]. It is possibly fair to say that at different disease stages different structures play a role. A biomarker would ideally show sensitivity across all disease stages and would reflect the underlying pathologic processes. From this cross-sectional analysis it is not possible to draw any definitive conclusions; however, the putamen possibly shows the greatest potential. Longitudinal evaluation has to be performed to confirm the putamen’s true potential. Furthermore we would like to stress how taking whole brain atrophy into account aids our interpretation of regional atrophy rates and in this way their potential as biomarkers can be assessed more effectively.

A limitation of our study could be the use of a relatively new software package from FSL used for segmentation, which has not specifically been validated in HD. Visual inspection, however, did not reveal any significant mismatches. The method applied gives structure segmentation on an individual basis and can therefore be used to compare groups. In contrast to this limitation several clear reasons exist in favor of using FIRST—first of all, compared to manual segmentation the automated segmentation uses voxel intensity in contrast to the sometimes difficult visual contrast differences, reducing a rater dependent bias. Secondly, this automated technique is suitable for implantation on large datasets, whereas manual segmentation is labour intensive and prone to human error.

In conclusion, we have demonstrated that atrophy of the pallidum and the accumbens nucleus exists in the premanifest stages of the disease and confirmed the well known atrophy of the caudate nucleus and putamen. Furthermore, we have shown the important correlation of the accumbens nucleus, putamen, pallidum and hippocampus with clinical symptoms. The importance of the remaining subcortical structures should be regarded when taking into account that they show similar amounts of atrophy as the brain as a whole. These findings have implications for biomarker research, as several subcortical structures now show great potential for use as a disease progression measure.