This paper describes a statistical model for volumetric MRI data analysis, which identifies the "change-point" when brain atrophy begins in premanifest Huntington's disease. Whole-brain mapping of the change-points is achieved based on brain volumes obtained using an atlas-based segmentation pipeline of T1-weighted images.
Recent advances in MRI offer a variety of useful markers to identify neurodegenerative diseases. In Huntington's disease (HD), regional brain atrophy begins many years prior to the motor onset (during the "premanifest" period), but the spatiotemporal pattern of regional atrophy across the brain has not been fully characterized. Here we demonstrate an online cloud-computing platform, "MRICloud", which provides atlas-based whole-brain segmentation of T1-weighted images at multiple granularity levels, and thereby, enables us to access the regional features of brain anatomy. We then describe a regression model that detects statistically significant inflection points, at which regional brain atrophy starts to be noticeable, i.e. the "change-point", with respect to a disease progression index. We used the CAG-age product (CAP) score to index the disease progression in HD patients. Change-point analysis of the volumetric measurements from the segmentation pipeline, therefore, provides important information of the order and pattern of structural atrophy across the brain. The paper illustrates the use of these techniques on T1-weighted MRI data of premanifest HD subjects from a large multicenter PREDICT-HD study. This design potentially has wide applications in a range of neurodegenerative diseases to investigate the dynamic changes of brain anatomy.
Magnetic resonance imaging (MRI) has substantially enhanced our ability to examine the brain anatomy and functions in neurodegenerative diseases1,2,3. T1-weighted structural MRI is one of most widely adopted imaging tools in routine clinical practice to assess the brain anatomy and related pathology. Quantitative analysis of the high-resolution T1-weighted images provides useful markers to measure anatomical changes during brain degeneration. In particular, segmentation based quantification approaches effectively reduces the image dimensionality from voxel level (on the order of (106)) to anatomical structural level ((102)) for high-throughput neuroinformatics4,5. Automated brain segmentation can be achieved using atlas-based methods6,7,8,9 that map the pre-defined anatomical labels from an atlas onto the patient images. Among the atlas-based methods, multi-atlas algorithms10,11,12,13,14 have yielded superior segmentation accuracy and robustness. Our group has developed a fully automated T1 multi-atlas segmentation pipeline, with advanced diffeomorphic image registration algorithms15, multi-atlas fusion methods16,17, and rich multi-atlas libraries18. The pipeline has been distributed on a cloud-computing platform, MRICloud19, since 2015, and it has been used to study neurodegenerative diseases, such as Alzheimer's disease (AD)20,21, Primary Progressive Aphasia22, and Huntington's Disease23.
Once the high-resolution images are segmented into brain structures, regional features, such as volumes, can be used to establish mathematical models to characterize the neuroanatomical changes. A change-point analysis method was recently established by our group to analyze the temporal order, in which statistically significant brain morphometric changes occur, based on longitudinal and/or cross-sectional MRI data. This statistical model was first developed to quantify shape-based diffeomorphometry over age in AD patients21,24; and it was later adapted to investigate brain structural changes in Huntington's disease (HD), as well as to describe brain developmental changes in neonatal brains25. In HD patients, the change-point was defined in with respect to the CAG-age product (CAP) score, as an indicator of the extent of exposure to the CAG expansion in HTT 26. It is well-known that striatal atrophy is one of the earliest markers in HD, followed by the globus pallidus27. Yet, the changes in striatum in relation to other gray and white matter structures across the brain remains unclear. Such relation is crucial for us to understand the disease progression. Change-point analysis of volumetric changes in all brain structures will likely provide systematic information of brain atrophy in premanifest phase of HD.
Here we demonstrate the procedures to perform whole-brain segmentation using MRICloud (www.mricloud.org), and steps to perform change-point analysis of volumetric data in premanifest HD subjects. The MRI data were collected from a large population multicenter PREDICT-HD study28,29 with approximately 400 controls and premanifest HD subjects. The combination of atlas-based segmentation and change-point analysis brings unique information about the spatiotemporal order of the brain structural changes and the disease progression pattern across the brain. The techniques are potentially applicable to a range of neurodegenerative diseases with various biomarkers to map the brain degeneration.
1. Atlas-based Whole Brain Segmentation
2. Change-point Analysis of Volumetric Data
NOTE: Please refer the theory and mathematical details of the change-point model to our previous publications21,24,23. Briefly, for HD data, a regression model is established to define a critical change-point (in terms of CAP score), starting from which, the linear regression between CAP and volume becomes statistical significant, while correcting for the effects of patient's age, gender, and intracranial volume. In-house Matlab scripts, along with a GUI (Figure 2A, lower panel), was developed to compute the change-points for individual structures and perform statistical analysis. The scripts are available to users upon request.
Using the procedures described in 1.1-1.3, whole brain segmentation maps can be obtained from MRICloud. In the current version of atlas (V9B), 283 parcels are segmented at the finest granularity (level 5), which can be grouped to different levels of granularity, e.g., from hemisphere to lobules and parcels, according to specific ontology definitions. Figure 3 shows two types of multi-level segmentations at five levels, in axial and coronal views. For example, at the coarsest level, Type-I segmentation defines classical brain regions of the telencephalon, diencephalon, mesencephalon, metencephalon, and myelencephalon (Figure 3A), while Type-II defines clinically and radiologically used conventions of hemispheres, cerebellum, and brainstem (Figure 3B). Here, we use Type-II definitions for the following analysis, and we average the volumes from the left and right hemispheres, as there is no known laterality in HD pathology.
Change-point analysis of the brain volumes is performed based on the segmentation results, according to section 2.1-2.2. In Figure 4, we take the basal ganglia as an example to examine the related change-points at different granularity levels. The 3D hierarchical structural definitions of basal ganglia are illustrated in the top row. Scatter plots of change-point analysis in individual structures are shown in the bottom row, where the blue dots denote the z-scores of the structural volumes from the premanifest HD patients (normalized to controls), after correcting for age, sex, and intracranial volumes; the black curves are the fitted volume z-scores, regressed to the change-point component in the model23; and the red lines indicates the positions of the detected change-points. At level 1, only hemispherical separation is available, and the corresponding scatter plot shows the whole hemisphere has a gradual atrophy as CAP score increases, with a change-point at CAP of 360 and 95% confidence interval of [352.8, 367.2] (p = 0.011 after FDR correction). At level 2, the basal ganglia is a part of the cerebral nuclei, and the cerebral nuclei shows a change-point at CAP of 232 [227.7, 236.3] (p <0.01 after FDR). At level 3, basal ganglia is an independent structure and it shows a change-point at CAP of 233 [228.6, 237.4] (p <0.01 after FDR). At level 4, basal ganglia is divided into the striatum and the globus pallidus, which have change-points at CAP of 230 [225.6, 234.4] and 243 [238.6, 247.4], respectively (both p <0.01 after FDR). At level 5, the striatum is further divided into caudate and putamen, which show change-points at CAP of 240 [234.9, 245.1] and 211 [206.8, 215.2], respectively (both p <0.01 after FDR).
Once the individual change-points are calculated for all brain structures, whole-brain spatial maps of the change-points can be obtained according to 2.3. Figure 5 shows the change-point maps at different granularity levels with Type-II ontological definitions. Note that only the structures with significant change-points (p <0.05 after FDR correction) are mapped. At level 1, the entire brain parenchyma (except cerebellum) shows a significant but relatively late change-point at CAP of 360 [352.8, 367.2] (p = 0.011 after FDR); whereas the cerebrospinal fluid (CSF, including lateral ventricle, III and IV ventricle, and sulci space) shows a slightly earlier change point at CAP of 319 [313.0, 325.0] (p=0.022 after FDR). The regions with significant change-points become more localized as the granularity level goes higher, and the spatial variance begins to unveil. For example, at level 3, the inferior white matter exhibits an earlier change-point (CAP = 305 [298.8, 311.2], p = 0.038 after FDR), compared to the anterior white matter (CAP = 371 [363.9, 378.1], p = 0.042 after FDR). At level 5, the differences among putamen, caudate, and globus pallidus become noticeable, with the putamen showing the earliest change-point, followed by the caudate and globus pallidus.
The change-rates are calculated as the percentage of volumetric changes (normalized to healthy controls) per CAP score, after the change-point. The change-rates of the corresponding structures are mapped in Figure 6, where the warm colors indicate volume increases in ventricles and sulci CSF and the cold colors indicate volume decreases in the brain parenchyma. At level 3, the deep gray matter structures show highest change rate, followed by the anterior white matter, and then inferior white matter. At level 5, the putamen and globus pallidus exhibit the fastest atrophy (0.1% volumetric loss per CAP), followed by the caudate (0.07% per CAP). As the granularity goes higher, the change-rates become higher in more localized regions.
Figure 1: Illustration of T1-weighted image segmentation on MRICloud. (A) Use of Dcm2Analyzed.exe to convert DICOM files to Analyzed format and perform de-identification. (B) Webpage for "T1-MultiAtlas Batch" processing. (C) Webpage for "My job status". (D) Webpage for visualization of segmentation results. Please click here to view a larger version of this figure.
Figure 2: Graphic User Interface (GUI) to perform volume analysis and change-point analysis. (A) Matlab GUI for batch extraction of volumetric data from MRICloud segmentation outputs (upper panel) and change-point analysis (lower panel). (B) MRICro interface to visualize the change-point maps. Please click here to view a larger version of this figure.
Figure 3: Atlas-based whole brain segmentation at multiple granularity levels with two types of ontological relationships. Axial and coronal views of multi-level segmentation maps are overlaid on T1-weighted anatomical images, according to Type-I (A) and Type-II (B) ontological definitions. Please click here to view a larger version of this figure.
Figure 4: Change-point analysis of the basal ganglia at multiple granularity levels. The hierarchical anatomical relations between hemisphere (level 1), cerebral nuclei (level 2), basal ganglia (level 3), striatum and globus pallidus (level 4), and putamen and caudate (level 5) are illustrated in the top row with 3D rendering. The scatter plots demonstrate change-point analysis of these structures, where the blue dots denote the z-scores of volumetric data (normalized to healthy controls) after correcting for age, sex, and intracranial volumes; the black curves are the fitted z-scores, regressed to the change-point dependent CAP component; and the red lines indicate the positions of the change-points. Please click here to view a larger version of this figure.
Figure 5: Whole brain change-point maps at multiple granularity levels. The regions that show significant change points (at 5% FDR) are mapped onto a T1-weighted image, and the colors indicate the change-point values in unit of CAP score. Abbreviations: WM:white matter; LV: lateral ventricle; AntWM: anterior white matter; InfWM: inferior white matter; BG: basal ganglia; ST: striatum; GP: globus pallidus; ALIC: anterior limb of the internal capsule; PLIC: posterior limb of the internal capsule; Caud: caudate; Put: putamen. Please click here to view a larger version of this figure.
Figure 6: Whole brain change-rate maps at multiple granularity levels. Change-rates are calculated as the percentage of volume changes per CAP, after the change-point (normalized to healthy controls), in the regions that show significant change-points (at 5% FDR), mapped onto a T1-weighted image. Abbreviations: WM: white matter; LV: lateral ventricle; AntWM: anterior white matter; InfWM: inferior white matter; BG: basal ganglia; ST: striatum; GP: globus pallidus; ALIC: anterior limb of the internal capsule; PLIC: posterior limb of the internal capsule; Caud: caudate; Put: putamen. Please click here to view a larger version of this figure.
As demonstrated in this paper, whole-brain segmentation of brain MRI can be conveniently achieved using our online platform MRICloud. T1-weighted MRI based volumetric marker has shown to be robust and sensitive to a range of neurodegenerative diseases1,2,3. The volumetric measures are used for various downstream analysis, such as mathematical modeling, and feature-selection and classification analysis to assist clinical diagnosis and prognosis. Change-point analysis of brain volumes allows quantitative characterization of brain atrophy during disease progression. This statistical analysis employs a regression model with a nonlinear component to define a change-point in the course of brain atrophy over a disease index, e.g., CAP score in HD. Compared to conventional group analysis that is widely adopted in most existing studies, the change-point model seeks for the exact onset points of brain atrophy, and therefore, offers more specific information of brain degeneration. Combining with an automated whole-brain segmentation pipeline, spatial maps of the change-points can be obtained, which reveal the spatiotemporal pattern of disease progression in HD. This is especially important in the premanifest phase of HD when therapeutic interventions can be most effective.
In the multi-atlas based image segmentation, brain segmentations at multiple granularities are provided to make the volumetric analysis flexible and scalable to various studies. For example, for Alzheimer's disease, several small brain structures are of particular interest, such as Hippocampus, Amygdala, or subdivisions of the temporal cortex; whereas for frontotemporal dementia, brains are examined at relatively coarse scales, such as the entire frontal and temporal lobes. According to our previous study33, the segmentation precision and reproducibility may be compromised at high granularity. In addition, the statistical challenge for multiple comparison would increase at higher granularity, as shown in Figure 5. In addition, an appropriate choice of atlas libraries, depending on age range and disease types of the study, is important for the segmentation accuracy18,32. The change-point analysis requires a relatively large number of data to perform the regression analysis, e.g., n >50 is usually needed based on empirical experience. Interpretation of change-point results heavily relies on the statistical test, because the model will produce a change-point for any dynamic data series, which sometimes is insignificant. We used a permutation test to evaluate whether the addition of change-point significantly reduces the model's residual error, and we used the bootstrap operation to mitigate the influence of outliers.
There can be variations to the change-point model depending on specific applications. For example, the disease index, which is CAP score for HD patients, can be replaced with other clinical variables, or simply age. Other covariates can also be included in the model, such as imaging protocols23 and other factors to be accounted for. Furthermore, the model is applicable to various types of biomarkers beyond volumetric marker, including other imaging markers (from different contrast mechanisms such as diffusion, perfusion, or functional MRI) and non-imaging markers. In addition, the mathematical formulation of the model can be adapted according to different assumptions and hypotheses. For examples, the model can involve a baseline linear regression before the change-point and an additional linear regression that occurs after the change-point, such as proposed in reference24.
The proposed change-point analysis in combination atlas-based whole-brain segmentation can have wide applications in both clinical and basic science studies of neurodegenerative diseases. Yet, there are limitations with the proposed techniques. The dynamic changes of biomarkers may not be linear, and thereby, linear component used in the current change-point model may not be optimal. Generalized linear model can potentially be integrated to handle nonlinear situations. Secondly, the change-point analysis requires large population longitudinal and/or cross-sectional data, such as the multi-center PREDICT-HD data used in this study. Different image acquisition protocols, from multi-center studies or longitudinal studies, may impose a bias in the volume measurements. Our multi-atlas image segmentation pipeline has shown to be relatively robust to protocol differences, compared to biological effects, such as age and pathology34. Still, it may be necessary to remove the protocol difference beforehand, or to include the protocol effect as one of the covariates in the change-point model as described in23. Currently, downstream volumetric analysis is not included in MRICloud, and users will need to extract volumes from the segmentation maps or use our Matlab toolbox to perform the multi-granularity volume analysis. We can potentially integrate this offline processing step into the cloud pipeline in the future, if this is identified as a common request by users.
The proposed change-point analysis in combination atlas-based whole-brain segmentation can have wide applications in both clinical and basic science studies of neurodegenerative diseases.
The authors have nothing to disclose.
We thank the PREDICT-HD investigators, particularly, Dr. Hans Johnson and Dr. Jane S. Pauslen from University of Iowa, for their generosity in sharing the MRI data and constructive discussion on the data analysis and results.
This work is supported by NIH grants R21 NS098018, P50 NS16375, NS40068, R01 NS086888, R01 NS084957, P41 EB015909, P41 EB015909, R01 EB000975, R01 EB008171, and U01 NS082085.
MATLAB | Mathworks | N/A | Version 2015b and above |
Dell Workstation | Dell | Dell Precision T5500 (Intel Xeon CPU) |