ORIGINAL RESEARCH article
Characterizing aging in the human brainstem using quantitative multimodal MRI analysis
- 1Clinical Neuroscience, St George’s University of London, London, UK
- 2Wellcome Trust Centre for Neuroimaging, UCL Institute of Neurology, University College London, London, UK
- 3Center for Neural Science, New York University, New York, NY, USA
- 4Department of Experimental Psychology, University of Oxford, Oxford, UK
- 5LREN, Département des Neurosciences Cliniques, CHUV, Université de Lausanne, Lausanne, Switzerland
Aging is ubiquitous to the human condition. The MRI correlates of healthy aging have been extensively investigated using a range of modalities, including volumetric MRI, quantitative MRI (qMRI), and diffusion tensor imaging. Despite this, the reported brainstem related changes remain sparse. This is, in part, due to the technical and methodological limitations in quantitatively assessing and statistically analyzing this region. By utilizing a new method of brainstem segmentation, a large cohort of 100 healthy adults were assessed in this study for the effects of aging within the human brainstem in vivo. Using qMRI, tensor-based morphometry (TBM), and voxel-based quantification (VBQ), the volumetric and quantitative changes across healthy adults between 19 and 75 years were characterized. In addition to the increased R2* in substantia nigra corresponding to increasing iron deposition with age, several novel findings were reported in the current study. These include selective volumetric loss of the brachium conjunctivum, with a corresponding decrease in magnetization transfer and increase in proton density (PD), accounting for the previously described “midbrain shrinkage.” Additionally, we found increases in R1 and PD in several pontine and medullary structures. We consider these changes in the context of well-characterized, functional age-related changes, and propose potential biophysical mechanisms. This study provides detailed quantitative analysis of the internal architecture of the brainstem and provides a baseline for further studies of neurodegenerative diseases that are characterized by early, pre-clinical involvement of the brainstem, such as Parkinson’s and Alzheimer’s diseases.
Aging is ubiquitous to the human condition. Previously, the cortical and subcortical correlates of normal aging have been extensively investigated using several modalities including volumetric MRI (Jernigan et al., 2001; Walhovd et al., 2005, 2011), quantitative MRI (qMRI) (Armstrong et al., 2004; Draganski et al., 2011; Bilgic et al., 2012), cortical thickness measures (Thambisetty et al., 2010), and diffusion tensor imaging (DTI) (Pfefferbaum et al., 2010; Vaillancourt et al., 2012). These studies have established that the normal aging process is associated with distinctive morphological changes including volumetric loss of the cortex, predominantly in the prefrontal, parietal, and temporal regions in addition to the amygdala, hippocampus, striatum, and cerebellum (Woodruff-Pak et al., 2010; Draganski et al., 2011; Walhovd et al., 2011). Corresponding with this, there is a reduction in cortical white-matter myelination, reflected by falling magnetization transfer (MT) values, and increasing iron deposition particularly in basal ganglia structures, notably the substantia nigra, as reflected by increasing R2* values (Pfefferbaum et al., 2009, 2010; Haacke et al., 2010; Draganski et al., 2011; Bilgic et al., 2012). Despite these widespread effects, reported brainstem related changes remain sparse due to technical limitations of imaging, segmenting, and statistically analyzing data from this region (Luft et al., 1999; Raz et al., 2001; Lee et al., 2009). Aging is known to negatively impact on several brainstem-mediated functions, for example the sleep-wake cycle (Hut and Van der Zee, 2011), sympathetic outflow (Samuels and Szabadi, 2008), vestibular-ocular reflexes (Baloh et al., 1993), and cardiovascular reflexes (Vita et al., 1986). Post-mortem reports indicate that changes within the brainstem sub-nuclei do take place during aging (Alvarez et al., 2000; Samuels and Szabadi, 2008), some of which may be early precursors for subclinical neurodegenerative disease (Tsopelas et al., 2011).
Quantitative MRI produces quantitative MR parameters that can be used as biomarkers of tissue microstructure (Tofts, 2003; Draganski et al., 2011). Examples of these parameters include MT, Proton Density (PD), R2*, and R1. MT emerges from hydrogen in motionally restricted macromolecules and more directly relates to macromolecular content. Biologically, it is a reflection of the quantity of myelin within a voxel (Helms et al., 2008b). It is important to note that in the literature MT contrast is often reported as “MT ratio” (Dousset et al., 1992), a parameter that shows residual T1 dependence. However, by using semi-quantitative parameter, the MT saturation, MT and T1 effect can be separated (Helms et al., 2008b). PD refers to the concentration of MRI-visible water (Tofts, 2003). R2* (=1/T2*) is the relaxation rate of the transverse magnetization, and is linearly correlated with tissue iron concentration (Yao et al., 2009). Finally R1 (=1/T1) is the longitudinal relaxation rate, and arises from a mix of water content, iron, and tissue macromolecule fraction (Rooney et al., 2007). Increases in iron (Rooney et al., 2007), decreases in PD (Gelman et al., 2001), and increases in lipid content (Stanisz et al., 2005) all cause an increase in the measured R1 signal.
The aim of this cross-sectional study was to characterize volumetric and tissue parameter changes associated with aging within the human brainstem in vivo using qMRI, tensor-based morphometry (TBM), and voxel-based quantification (VBQ). This was motivated by the observation that certain neurodegenerative diseases, such as Parkinson’s and Alzheimer’s diseases, are characterized by early pre-clinical involvement of the brainstem (Simic et al., 2009; Hawkes et al., 2010). In order to further study these effects, changes associated with the normal aging process must first be better characterized.
Materials and Methods
Multiparametric maps were acquired from previous studies (FitzGerald et al., 2012; Lambert et al., 2012; Chowdhury et al., 2013). In total, imaging data for 100 healthy adults (47 males, aged 40.3 ± 21.2 years; 53 females, aged 48.2 ± 22.7 years; age ranged 19–75 years) who had MRI scanning at the Wellcome Trust Centre for Neuroimaging was used. Due to the selection criteria, these predominately fell into a bimodal distribution (as is shown Figure 5), with a younger cohort less than 60 years [n = 58, mean age 25.8 years (SD 7.6 years)] and an older cohort above 60 years [n = 42, mean age 69.1 years (SD 3.5 years)]. Subjects above 60 years had a normal neurological examination performed by a physician (Rumana Chowdhury), MMSE score>28 and a normal performance (within 1.5 SD of age-related norm) on a range of neuropsychological tests. Involvement of human volunteers was approved by the local ethics committee, and each subject provided written informed consent prior to MRI examination.
MR imaging was performed on a 3T whole-body MRI system (Magnetom TIM Trio, Siemens Healthcare, Erlangen) operated with a whole-body transmit radio-frequency (RF) coil and a 32-channel RF receive coil. MR data of 21 subjects was acquired on a second identical MRI system located within the same department. Each participant underwent a multiparametric mapping (MPM) scanning protocol for quantitative mapping of multiple MR parameters. MT, R1 (=1/T1), and PD weighted images were acquired using 3D multi-echo FLASH (fast low-angle shot) acquisitions (Helms et al., 2008a,b). Full imaging parameters are summarized in Table 1. The image resolution was 1 mm isotropic. The total acquisition time was 19 min. For each subject quantitative MT, R1, PD, and R2* (=1/T2*) maps were extracted from the acquired images using in-house MATLAB program. An additional dataset was acquired on each subject for mapping of the RF transmit field B1+ over the brain (4 mm isotropic resolution, acquisition time 3 min) using the 3D EPI SE/STE method described in (Lutti et al., 2010, 2012). A B0-field map was also acquired to correct for the distortions of the EPI images acquired for the B1+ mapping acquisition (acquisition time 2 min). The resulting B1+ maps were used to correct the MPM maps for RF transmit field inhomogeneity effects (Helms et al., 2008a). The PD maps were corrected for the spatially varying sensitivity profile of the receive coil using the UNICORT algorithm (Weiskopf et al., 2011). The resulting “flattened” signal amplitude maps were converted into PD maps by scaling the voxels by the expected average PD of white-matter [69% (Tofts, 2003)].
The MT maps were initially segmented into gray, white, and CSF tissue classes whilst maintaining the native resolution (1 mm) by using the unified segmentation within SPM8 (http://www.fil.ion.ucl.ac.uk/spm/) (Ashburner and Friston, 2005). For each individual, the total intracranial volume (TIV) was calculated by summing the gray matter, white-matter, and CSF segmentations at a threshold of 0.1. The segmentations were then registered diffeomorphically to a common group-average 1 mm isotropic template (using the “Shoot toolbox” of SPM, Ashburner and Friston, 2011) to produce deformation fields that were subsequently used in the segmentation step to warp the brainstem tissue priors to individual subject space.
A method for generating brainstem specific tissue probability maps (TPMs) and subsequent segmentation had been previously developed (Lambert et al., 2013). In brief, the brainstem TPMs were defined using a modified multivariate Mixture of Gaussians (mmMoG) to generate spatial TPMs from 0.8 mm isotropic MT and PD maps, which were masked to include the brainstem from the origin of the cerebral aqueduct to the level of the foramen magnum. Four tissue classes were identified, three gray matter and one white-matter. These were labeled for descriptive purposes according to the predominant tissue type present: tissue class one included the substantia nigra, locus coeruleus and raphe nuclei and hence was designated “monoaminergic gray matter,” though regions that contained dorsal cranial nerve nuclei were also included. Tissue class two consisted mainly of the nucleus reticularis throughout its length and the pontine nuclei, and hence was designated “reticulated gray matter.” Tissue class three was specific for the periaqueductal gray (PAG) and labeled as such. Tissue class four was the brainstem white-matter. Examples of these tissue classes projected onto the corresponding sections from Duvernoy’s 9.4T MRI brainstem atlas (Naidich and Duvernoy, 2009) have been provided in Figure 1 (reprinted from Lambert et al., 2013 with permission from Elsevier).
Figure 1. Comparison of brainstem tissue classes against three corresponding ex vivo brainstem sections from “MR microscopy at 9.4T” (taken from Duvernoy’s Atlas of the Human Brain Stem and Cerebellum with permission). Figure reproduced from Lambert et al. (2013) with permission from Elsevier.
The previously calculated tissue priors for four brainstem and one non-brainstem tissue classes were aligned with the calculated group-average template. This was achieved by realigning a whole brain template in the same space as the TPMs with the corresponding new group-average template, re-slicing to achieve 1 mm isotropic voxel size with a 4th degree spline interpolation and ensuring each individual voxel probability value summed to one across all the maps.
The actual brainstem segmentation step is summarized in Figure 2. The realigned probability maps were warped to individual subject space and used in SPM8 “New Segment” on individual MT and PD maps that had been cropped with a set-bounding box to include only the brainstem region. Specifically, five tissue classes were used; four within brainstem and one for everything else. Two Gaussians were used to model each tissue class except for the PAG matter (one Gaussian) and non-brainstem (eight Gaussians). Individual level segmentations were generated from the MT and PD images. All the images from each of the four tissue classes were then cropped using a common bounding box to increase computational speed and visually checked to ensure good quality. All of the brainstem tissue classes were then re-warped back to a group-average using a diffeomorphic warping algorithm (geodesic shooting) to allow re-estimation of the Jacobian determinants. The deformation fields were then used to warp each parametric map (MT, PD, R1, R2*) to the common template. The resulting images were masked using the non-brainstem tissue class as an exclusion mask.
Figure 2. Summary of segmentation pipeline. This includes examples of the Multiparameter Maps (MPMs) as shown at the top of the image.
This work utilized two techniques to quantitatively analyze qMRI parameter maps. The first was TBM, which is a method to characterize volumetric change in vivo. It is similar to the already well-described voxel-based morphometry (Ashburner and Friston, 2000), however instead of statistically analyzing warped modulated tissue segmentations, the Jacobian determinant images are used (Ashburner and Friston, 2001). The second was VBQ. This is a pipeline to allow unbiased mass univariate statistical analysis of quantitative parameter maps whilst controlling for multiple comparisons with family-wise error (Draganski et al., 2011). The original VBQ method (Draganski et al., 2011) was modified slightly for this study. Because the brainstem lacks gyrification, and due to the highly accurate warping algorithm used, smoothing of the warped quantitative data was avoided, hence negating the need to produce “warped-weighted average” images that were required for the published approach. This had several advantages. First, only one statistical analysis was required for each quantitative map (i.e., the warped quantitative map) rather than one for each tissue class (i.e., four), which would have been necessary with the warped-weighted average approach. Additionally, because analysis was carried out directly on the warped quantitative maps without smoothing, a higher degree of spatial accuracy could be achieved.
All statistical analysis was carried out using SPM8 in MATLAB 2010b. A two-sample t-test was initially used to check for scanner-associated differences in the acquisitions controlling for age, sex, and TIV, at FWE < 0.05 correction for multiple comparisons. No voxels survived correction. A design matrix for multiple linear regression model was then constructed including age, sex, and TIV as covariates. The TIV was centered around the mean, and the remaining covariates remained un-centered. An intercept was included in the model but no normalization was used. Using this design matrix, TBM was performed by analyzing the Jacobian determinant images (Ashburner and Friston, 2001), and VBQ by analyzing each warped quantitative map. For each analysis, the non-brainstem tissue class was used as an exclusion mask, ensuring only voxels that contained brainstem tissue were included. Each image type was assessed in a single design matrix to negate potential problems associated with uneven variances across the different quantitative maps. Each map was assessed for significant positive and negative correlates with age. For each contrast, FWE < 0.05 was reported. Finally, for each VBQ analysis where results were significant at FWE < 0.05, the T-maps were binarized for visualization at p < 0.001 uncorrected. These were used to create three-dimensional renderings, and also to assess the overlap of significant results between different modalities to better understand the interaction between the different measures. Finally, the same binarized images were used to extract the mean quantitative value for each individual. These values were plotted against the corresponding age, and the best linear fit indicated.
For anatomical reference, a high resolution (0.5 mm isotropic) ex vivo combined MT T2* MRI with anatomical annotations has been provided in Figure 3 (figure adapted from Lambert et al., 2013 with permission from Elsevier).
Figure 3. High resolution ex vivo combined MT T2* MRI with anatomical annotations for reference. Figure adapted from Lambert et al. (2013) with permission from Elsevier. Abbreviations: A8, dopaminergic center (approximate location), CP, cerebral peduncle (anterior to posterior: consisting of frontopontine, corticonuclear, corticospinal, and parietotemporal pontine tracts); CST, corticospinal tract; CTT, central tegmental tract; ICP, inferior cerebellar peduncle; ML, medial lemniscus; MLF, medial longitudinal fasciculus; PAG, periaqueductal gray; SCT, spinocerebellar tract; SNpc, substantia nigra pars compacta; SNpr, substantia nigra pars reticulata; TST, tectospinal tract; VTA, ventral tegmental area. *Artifact due to fixation.
Highly localized volumetric decreases in tissue volume were found symmetrically within the brachium conjunctivum (superior cerebellar peduncle) bilaterally. These volumetric decreases were significant at FWE < 0.05 and are summarized in Figures 4 and 5. There were no significant positive TBM correlates with age at FWE < 0.05 and p < 0.001.
Figure 4. Areas of significant regional decreases in tissue volume and magnetization transfer, binarized at p < 0.001 uncorrected.
Figure 5. Scatter plots of mean individual values extracted from binarized T-maps at p < 0.001 uncorrected. Linear fit shown in red. Increases: (A) Proton Density; (B) R1; (C) Magnetization Transfer; (D) R2*. Decreases: (E). Magnetization Transfer; (F). TBM. All Pearson’s correlation coefficients significant at p < 1 × 10−8.
Voxel-Based Quantification Analysis
Negative correlates with age
Positive correlates with age
As shown in Figures 5–7, there were more widespread significant increases in qMRI values with age. Figure 6 shows the individual parameter maps increases at six axial slices through the brainstem. Figure 7 summarizes the regional MPM increases and examines the overlap between the different parameter maps. It also demonstrates the relative proportions of each tissue type that overlap with one another.
Figure 6. Axial, coronal, and sagittal slice-wise images showing positive correlates with age (R2s, MT, A, R1) projected onto group-average MT map.
Figure 7. Significant regional increases binarized at p < 0.001 uncorrected. The legend classifies the regions of significant increase across the MPM, allowing overlapping significant increases to be simultaneously visualized. The histogram takes the total number of all significant voxels shown in the rendering, and then classifies them according to the MPMs in which those voxels are found in as per the legend.
R1 increases with age. Increases in R1 intensity values with age were found throughout the brainstem, but were confined to gray matter structures. These were found in bilateral inferior olivary nuclei, pontine nuclei, dorsal raphe nucleus, and bilateral substantia nigra. To better understand this finding, a sub analysis on these significant regions was performed, looking at the remaining parameters. Whilst there may be some selection bias with this approach (Vul et al., 2009), the correlation analysis was strictly performed to compare the behavior of the different MRI modalities within the regions of R1 change, to better understand what may be contributing to this R1 increase. A significant positive correlation with R2* was found within these regions, even when the iron rich midbrain structures were excluded (Pearson’s Correlation Coefficient within pons and medulla = 0.51, p < 1 × 10−8), implying that it is increasing iron deposition that significantly contributes to the observed R1 signal increases within brainstem gray matter structures.
R2* increases with age. Increases in R2* intensity values were localized to midbrain structures. Specifically, the substantia nigra, ventral tegmental area, and red nuclei. The pontine and medullary R1 regions described above did not show increases in R2s when analyzing the entire brainstem. This discrepancy would be due to correction for multiple comparisons at the whole brainstem level compared to a small volume region of interest.
Proton density increases with age. Significant increases in PD were observed in the brainstem white-matter. This included (from rostral to caudal): Fasciculus cerebellothalamicus, brachium conjunctivum, corticospinal tract, superior cerebellar peduncle, and medial longitudinal fasciculus. Additionally, the inferior portion of the substantia nigra also exhibited increased PD.
MT increases with age. Only one small region demonstrated MT increases: the corticobulbar portion of the cerebral peduncle bilaterally.
In this study, we have demonstrated the age effect of MT, R1, R2*, and PD in human brainstem using an automated, non-biased approach that is able to resolve the internal structure of the brainstem in unprecedented detail at 3T.
Midbrain Atrophy in Healthy Aging
Previous studies examining age-related volumetric decline in the brainstem have found no overall volume loss (Raz et al., 2001; Lee et al., 2009) but significant midbrain atrophy (Luft et al., 1999). This had previously been attributed to shrinkage of the substantia nigra (Raz, 1996), however there is sparse histopathological evidence to support this hypothesis (McCormack et al., 2004; Collier et al., 2007). Our study agrees with previous work in that brainstem aging-associated atrophy seems to be confined to the midbrain. However, our results indicate that it is predominantly volume loss of the superior cerebellar fiber bundles (brachium conjunctivum, fasciculus cerebellothalamicus (Haroian et al., 1981)] that are responsible for this finding. Additionally the decreasing myelin content, reflected by decreasing MT and increasing PD, indicates this is a regional effect due to axonal loss rather than an artifact. What is striking is the regional specificity of these findings. As with previous studies (Draganski et al., 2011), we found a marked increase in iron concentration with age in the substantia nigra and red nucleus. Though this increase is largely sequestered in neuromelanin in vivo, this substance will be released by dying neurons and hence can contribute to the regional damage observed (Chiueh, 2001; Papanikolaou and Pantopoulos, 2005). Not only is iron directly toxic to axons by causing rapid lipid peroxidation, it can also induce neurotoxic microglial factors to be released locally that will potentiate the regional insult (Zecca et al., 2004). These findings correlate with previously described cerebellar atrophy (Andersen et al., 2003; Draganski et al., 2011; Walhovd et al., 2011) and the resultant cortical and subcortical disconnection (Taniwaki et al., 2007; Alalade et al., 2011) hence would better account for these observations.
Brainstem Gray Matter Changes
Our results agree with previous studies in the observation that the pons and medulla do not show volumetric loss in aging (Sullivan et al., 2004; Lee et al., 2009). Despite this, we identified widespread and specific changes in the qMRI maps within the gray matter structures of these regions. Specifically, increasing R1 signal was observed within the pontine nuclei and inferior olive, in addition to the substantia nigra and dorsal raphe nuclei within the mesencephalon. This increase in R1 signal was significantly correlated with R2*, suggesting that increasing iron content within these structures can at least partly account for the observed gray matter changes. Within the pons, as with elsewhere in the brain, iron is predominately found in oligodendrocytes and to a lesser extent in the microglia and astrocytes (Ozawa et al., 1994). However the amount of iron in the former remains constant throughout aging and instead increases are found within the microglia and astrocytes (Zecca et al., 2004), which may relate to changes in vascular permeability. It has been speculated that the accumulation of immunoreactive iron in the microglia may cause or predispose to a neuroinflammatory response, as seen in Parkinson’s and Alzheimer’s disease (Zecca et al., 2004).
Chemical exchange brings a significant contribution to the transfer of magnetization between free water and macromolecules (Henkelman et al., 1993; Kucharczyk et al., 1994). Evidence has been presented which also suggests a significant impact of chemical exchange on the resonance frequency of water protons (Shmueli et al., 2011; Wharton and Bowtell, 2012). Physiological factors that impact chemical exchange such as tissue oxygenation (O2 and CO2 exchange), temperature, and pH (Kucharczyk et al., 1994; Liepinsh and Otting, 1996; van Zijl et al., 2003) may alter the MT and R2* estimates presented here. However only small variations of MT with pH have previously been reported within biologically plausible values (Kucharczyk et al., 1994). T1 relaxation is primarily driven by dipolar interactions between water and macromolecular protons (Koenig, 1995). Additionally the small variations in MT and T1 values of brain tissue with temperature (Lewa and Majewska, 1980; Graham et al., 1999) are unlikely to have a significant impact on our results within a range of biologically reasonable temperatures. Decreases in cerebral blood flow (CBF) that vary from region to region have been reported in normal aging (Aanerud et al., 2012). Modulated by the physiological parameters listed above, CBF may have an impact on the MPM measurements. However, the precise interaction between these remains poorly understood (Zauner and Muizelaar, 1997) and their impact on aging sparsely characterized. Several findings in this study suggest that the impact of CBF on healthy aging is minimal in the brainstem. First, biologically plausible concordance was accomplished between different parameter maps, for example the decreasing MT with increasing PD in areas of axonal and volume loss. Second, the results demonstrate respect for known anatomical boundaries such as the brachium conjunctivum with increased PD. Finally, our results are in agreement with previous histological observations, such as the increases in iron content in the substantia nigra (Zecca et al., 2001). Future work is expected to truly clarify and disambiguate the effects of CBF on MPM measurements. Not only will this improve the biophysical interpretation of these sequences in healthy tissue, but also allow better understanding of cerebral pathology where the vascular permeability will also change (Mooradian, 1988).
There are several limitations with our study. First, to accrue 100 normal control MPMs, we utilized scans that were acquired through previous studies (FitzGerald et al., 2012; Lambert et al., 2012; Chowdhury et al., 2013). There were few subjects between the ages of 35–65. Whilst this does not invalidate the findings, it makes it impossible to better characterize the temporal characteristics of the changes i.e., are they linear or non-linear. Additionally, as the upper age limit is 75, it is unclear how these changes extrapolate to those over that age. These features also bias the study toward those who are ambulatory, independent, and self motivating. Whilst it could also be argued that the latter criticism ensures that the experimental group consists only of healthy normal controls, it must also be acknowledged that declining mobility is a recognized feature of normal aging that may be unrelated to cortical changes, and this would certainly represent a selection bias in this work. Finally, it is currently unclear how these changes map to individual function, so further work is required to better understand this aspect.
In this work, the impact of physiological noise on image quality and the segmentation results was not explored. However, recent work using these techniques in the cortex highlights its robustness and also sensitivity to tissue microarchitecture (Dick et al., 2012; Sereno et al., 2013). Physiological noise has mostly been addressed in the context of fMRI, where image stability is paramount (Glover et al., 2000; Hutton et al., 2011). These methods cannot be directly implemented in anatomical imaging due to the different type of image acquisition. Potentially beneficial techniques include phase-navigator correction methods (Hu and Kim, 1994; Barry et al., 2008) although they may reduce the efficiency of the FLASH acquisitions. Alternatively real-time shimming methods for correction of respiratory-induced effects (Van Gelderen et al., 2007) or optical systems for fast prospective correction of subject motion (Zaitsev et al., 2006) may yield a significant reduction of physiological effects on anatomical scans.
Regarding the scanning parameters, 1 mm isotopic volumes are still reasonably large for certain brainstem structures, so our ability to fully characterize the changes are likewise limited. However, many structures in the brainstem are well above the 1 mm3 threshold such as the facial nerve nucleus (mean volume = 12.95 mm3) and hypoglossal nerve nucleus (mean volume = 14.39 mm3) (Sherwood et al., 2005), and hence this current study is of sufficient resolution for these nuclei. Additionally, we have demonstrated that widespread and structure specific changes can be found that also correspond to known features of brainstem aging.
This work provided a baseline of qMRI changes with aging in the brainstem. Further work is required to characterize the exact temporal dynamics of these changes, and how they extrapolate to those beyond the age of 75. Combining these techniques with diffusion-weighted imaging would allow further characterization of the changes within the white-matter and connectivity properties, and could be used to further refine the segmentations. Future work will examine these parameters in the context of the Parkinsonian disorders, such as idiopathic Parkinson’s disease, progressive supranuclear palsy, and multisystem atrophy. These conditions are good models of brainstem disease with which to develop the imaging techniques, as they all show post-mortem histological changes within the brainstem that, to date, are only evident on MRI at advanced stages of disease. Additionally, further longitudinal work is required to understand how individuals whose qMRI parameters lay well outside their expected normal age-matched distributions develop over time, to identify whether these parameters can serve as early biomarkers of neurodegenerative disease, an important prerequisite for any future neuroprotective therapy.
In conclusion, we have characterized changes in the brainstem due to normal healthy aging using qMRI and volumetric analysis. We replicate previous findings of midbrain shrinkage, and put forward a new hypothesis as to the underlying mechanism based on statistically significant regional changes in the qMRI maps. Specifically, axonal loss in the ascending cerebellar fiber bundles is responsible for the decreased brainstem volume loss rather than atrophy of the substantia nigra, as previously speculated. Additionally, it may be that the regional increases in nigral and rubral iron content, as reflected by R2* signal, underpin these observations. Finally we demonstrate widespread brainstem changes during aging evidenced by increasing PD in the white-matter, and increasing R1 in the gray matter, the latter being significantly driven by increasing iron deposition. This work provides a baseline from which brainstem pathology can be better explored in vivo using 3T MRI, and non-invasive biomarkers of different neurodegenerative conditions.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by Wellcome Trust Grant 075696/Z/04/Z (R.S.J.F., Sarah Tabrizi, John Ashburner). The Wellcome Trust Centre for Neuroimaging is supported by core funding from the Wellcome Trust 091593/Z/10/Z. Christian Lambert is supported by the Academy of Medical Sciences Starter Grants for Clinical Lecturers. Bogdan Draganski is supported by the Swiss National Science Foundation (320030_135679 and NCCR Synapsy), Foundation Parkinson Switzerland, Deutsche Forschungsgemeinschaft (KFO 247/0), Novartis Foundation for medical-biological research and the Synapsis Foundation. We thank all participants in our study and the radiographers at the Functional Imaging Laboratory for their assistance acquiring data.
Aanerud, J., Borghammer, P., Chakravarty, M. M., Vang, K., Rodell, A. B., Jónsdottir, K. Y., et al. (2012). Brain energy metabolism and blood flow differences in healthy aging. J. Cereb. Blood Flow Metab. 32, 1177–1187. doi:10.1038/jcbfm.2012.18
Alvarez, J. C., Díaz, C., Suárez, C., Fernández, J. A., González del Rey, C., Navarro, A., et al. (2000). Aging and the human vestibular nuclei: morphometric analysis. Mech. Ageing Dev. 114, 149–172. doi:10.1016/S0047-6374(00)00098-1
Armstrong, C. L., Traipe, E., Hunter, J. V., Haselgrove, J. C., Ledakis, G. E., Tallent, E. M., et al. (2004). Age-related, regional, hemispheric, and medial-lateral differences in myelin integrity in vivo in the normal adult brain. AJNR Am. J. Neuroradiol. 25, 977–984.
Barry, R. L., Martyn Klassen, L., Williams, J. M., and Menon, R. S. (2008). Hybrid two-dimensional navigator correction: a new technique to suppress respiratory-induced physiological noise in multi-shot echo-planar functional MRI. Neuroimage 39, 1142–1150. doi:10.1016/j.neuroimage.2007.09.060
Bilgic, B., Pfefferbaum, A., Rohlfing, T., Sullivan, E. V., and Adalsteinsson, E. (2012). MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. Neuroimage 59, 2625–2635. doi:10.1016/j.neuroimage.2011.08.077
Chowdhury, R., Guitart-Masip, M., Lambert, C., Dolan, R. J., and Düzel, E. (2013). Structural integrity of the substantia nigra and subthalamic nucleus predicts flexibility of instrumental learning in older-age individuals. Neurobiol. Aging 34, 2261–2270. doi:10.1016/j.neurobiolaging.2013.03.030
Collier, T. J., Lipton, J., Daley, B. F., Palfi, S., Chu, Y., Sortwell, C., et al. (2007). Aging-related changes in the nigrostriatal dopamine system and the response to MPTP in nonhuman primates: diminished compensatory mechanisms as a prelude to parkinsonism. Neurobiol. Dis. 26, 56–65. doi:10.1016/j.nbd.2006.11.013
Dick, F., Tierney, A. T., Lutti, A., Josephs, O., Sereno, M. I., and Weiskopf, N. (2012). In vivo functional and myeloarchitectonic mapping of human primary auditory areas. J. Neurosci. 32, 16095–16105. doi:10.1523/JNEUROSCI.1712-12.2012
Dousset, V., Grossman, R. I., Ramer, K. N., Schnall, M. D., Young, L. H., Gonzalez-Scarano, F., et al. (1992). Experimental allergic encephalomyelitis and multiple sclerosis: lesion characterization with magnetization transfer imaging. Radiology 182, 483–491.
Draganski, B., Ashburner, J., Hutton, C., Kherif, F., Frackowiak, R. S., Helms, G., et al. (2011). Regional specificity of MRI contrast parameter changes in normal ageing revealed by voxel-based quantification (VBQ). Neuroimage 55, 1423–1434. doi:10.1016/j.neuroimage.2011.01.052
FitzGerald, T. H., Friston, K. J., and Dolan, R. J. (2012). Action-specific value signals in reward-related regions of the human brain. J. Neurosci. 32, 16417–16423. doi:10.1523/JNEUROSCI.3254-12.2012
Gelman, N., Ewing, J. R., Gorell, J. M., Spickler, E. M., and Solomon, E. G. (2001). Interregional variation of longitudinal relaxation rates in human brain at 3.0 T: relation to estimated iron and water contents. Magn. Reson. Med. 45, 71–79. doi:10.1002/1522-2594(200101)45:1<71::AID-MRM1011>3.0.CO;2-2
Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. 44, 162–167. doi:10.1002/1522-2594(200007)44:1<162::AID-MRM23>3.3.CO;2-5
Graham, S. J., Stanisz, G. J., Kecojevic, A., Bronskill, M. J., and Henkelman, R. M. (1999). Analysis of changes in MR properties of tissues after heat treatment. Magn. Reson. Med. 42, 1061–1071. doi:10.1002/(SICI)1522-2594(199912)42:6<1061::AID-MRM10>3.0.CO;2-T
Griswold, M. A., Jakob, P. M., Heidemann, R. M., Nittka, M., Jellus, V., Wang, J., et al. (2002). Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson. Med. 47, 1202–1210. doi:10.1002/mrm.10171
Haacke, E. M., Miao, Y., Liu, M., Habib, C. A., Katkuri, Y., Liu, T., et al. (2010). Correlation of putative iron content as represented by changes in R2* and phase with age in deep gray matter of healthy adults. J. Magn. Reson. Imaging 32, 561–576. doi:10.1002/jmri.22293
Haroian, A. J., Massopust, L. C., and Young, P. A. (1981). Cerebellothalamic projections in the rat: an autoradiographic and degeneration study. J. Comp. Neurol. 197, 217–236. doi:10.1002/cne.901970205
Helms, G., Dathe, H., Kallenberg, K., and Dechent, P. (2008b). High-resolution maps of magnetization transfer with inherent correction for RF inhomogeneity and T1 relaxation obtained from 3D FLASH MRI. Magn. Reson. Med. 60, 1396–1407. doi:10.1002/mrm.21732
Henkelman, R. M., Huang, X., Xiang, Q. S., Stanisz, G. J., Swanson, S. D., and Bronskill, M. J. (1993). Quantitative interpretation of magnetization transfer. Magn. Reson. Med. 29, 759–766. doi:10.1002/mrm.1910290607
Hutton, C., Josephs, O., Stadler, J., Featherstone, E., Reid, A., Speck, O., et al. (2011). The impact of physiological noise correction on fMRI at 7T. Neuroimage 57, 101–112. doi:10.1016/j.neuroimage.2011.04.018
Jernigan, T. L., Archibald, S. L., Fennema-Notestine, C., Gamst, A. C., Stout, J. C., Bonner, J., et al. (2001). Effects of age on tissues and regions of the cerebrum and cerebellum. Neurobiol. Aging 22, 581–594. doi:10.1016/S0197-4580(01)00217-2
Kucharczyk, W., Macdonald, P. M., Stanisz, G. J., and Henkelman, R. M. (1994). Relaxivity and magnetization transfer of white matter lipids at MR imaging: importance of cerebrosides and pH. Radiology 192, 521–529.
Lambert, C., Lutti, A., Helms, G., Frackowiak, R., and Ashburner, J. (2013). Multiparametric brainstem segmentation using a modified multivariate mixture of Gaussians. Neuroimage (Amst.) 2, 684–694. doi:10.1016/j.nicl.2013.04.017
Lambert, C., Zrinzo, L., Nagy, Z., Lutti, A., Hariz, M., Foltynie, T., et al. (2012). Confirmation of functional zones within the human subthalamic nucleus: patterns of connectivity and sub-parcellation using diffusion weighted imaging. Neuroimage 60, 83–94. doi:10.1016/j.neuroimage.2011.11.082
Lee, N. J., Park, I. S., Koh, I., Jung, T. W., and Rhyu, I. J. (2009). No volume difference of medulla oblongata between young and old Korean people. Brain Res. 1276, 77–82. doi:10.1016/j.brainres.2009.04.027
Luft, A. R., Skalej, M., Schulz, J. B., Welte, D., Kolb, R., Bürk, K., et al. (1999). Patterns of age-related shrinkage in cerebellum and brainstem observed in vivo using three-dimensional MRI volumetry. Cereb. Cortex 9, 712–721. doi:10.1093/cercor/9.7.712
Lutti, A., Hutton, C., Finsterbusch, J., Helms, G., and Weiskopf, N. (2010). Optimization and validation of methods for mapping of the radiofrequency transmit field at 3T. Magn. Reson. Med. 64, 229–238. doi:10.1002/mrm.22421
Lutti, A., Stadler, J., Josephs, O., Windischberger, C., Speck, O., Bernarding, J., et al. (2012). Robust and fast whole brain mapping of the RF transmit field B1 at 7T. PLoS ONE 7:e32379. doi:10.1371/journal.pone.0032379
McCormack, A. L., Di Monte, D. A., Delfani, K., Irwin, I., DeLanney, L. E., Langston, W. J., et al. (2004). Aging of the nigrostriatal system in the squirrel monkey. J. Comp. Neurol. 471, 387–395. doi:10.1002/cne.20036
Naidich, T. P., and Duvernoy, H. M. (2009). Duvernoy’s Atlas of the Human Brain Stem and Cerebellum: High-Field MRI: Surface Anatomy, Internal Structure, Vascularization and 3D Sectional Anatomy, Illustrated Ed. Wien, New York: Springer.
Ozawa, H., Nishida, A., Mito, T., and Takashima, S. (1994). Development of ferritin-containing cells in the pons and cerebellum of the human brain. Brain Dev. 16, 92–95. doi:10.1016/0387-7604(94)90041-8
Pfefferbaum, A., Adalsteinsson, E., Rohlfing, T., and Sullivan, E. V. (2009). MRI estimates of brain iron concentration in normal aging: comparison of field-dependent (FDRI) and phase (SWI) methods. Neuroimage 47, 493–500. doi:10.1016/j
Pfefferbaum, A., Adalsteinsson, E., Rohlfing, T., and Sullivan, E. V. (2010). Diffusion tensor imaging of deep gray matter brain structures: effects of age and iron concentration. Neurobiol. Aging 31, 482–493. doi:10.1016/j.neurobiolaging.2008.04.013
Raz, N., Gunning-Dixon, F., Head, D., Williamson, A., and Acker, J. D. (2001). Age and sex differences in the cerebellum and the ventral pons: a prospective MR study of healthy adults. AJNR Am. J. Neuroradiol. 22, 1161–1167.
Rooney, W. D., Johnson, G., Li, X., Cohen, E. R., Kim, S. G., Ugurbil, K., et al. (2007). Magnetic field and tissue dependencies of human brain longitudinal 1H2O relaxation in vivo. Magn. Reson. Med. 57, 308–318. doi:10.1002/mrm.21122
Samuels, E. R., and Szabadi, E. (2008). Functional neuroanatomy of the noradrenergic locus coeruleus: its roles in the regulation of arousal and autonomic function part II: physiological and pharmacological manipulations and pathological alterations of locus coeruleus activity in humans. Curr. Neuropharmacol. 6, 254. doi:10.2174/157015908785777193
Sherwood, C. C., Hof, P. R., Holloway, R. L., Semendeferi, K., Gannon, P. J., Frahm, H. D., et al. (2005). Evolution of the brainstem orofacial motor system in primates: a comparative study of trigeminal, facial, and hypoglossal nuclei. J. Hum. Evol. 48, 45–84. doi:10.1016/j.jhevol.2004.10.003
Simic, G., Stanic, G., Mladinov, M., Jovanov-Milosevic, N., Kostovic, I., and Hof, P. R. (2009). Does Alzheimer’s disease begin in the brainstem? Neuropathol. Appl. Neurobiol. 35, 532–554. doi:10.1111/j.1365-2990.2009.01038.x
Stanisz, G. J., Odrobina, E. E., Pun, J., Escaravage, M., Graham, S. J., Bronskill, M. J., et al. (2005). T1, T2 relaxation and magnetization transfer in tissue at 3T. Magn. Reson. Med. 54, 507–512. doi:10.1002/mrm.20605
Sullivan, E. V., Rosenbloom, M., Serventi, K. L., and Pfefferbaum, A. (2004). Effects of age and sex on volumes of the thalamus, pons, and cortex. Neurobiol. Aging 25, 185–192. doi:10.1016/S0197-4580(03)00044-7
Taniwaki, T., Okayama, A., Yoshiura, T., Togao, O., Nakamura, Y., Yamasaki, T., et al. (2007). Age-related alterations of the functional interactions within the basal ganglia and cerebellar motor loops in vivo. Neuroimage 36, 1263–1276. doi:10.1016/j.neuroimage.2007.04.027
Thambisetty, M., Wan, J., Carass, A., An, Y., Prince, J. L., and Resnick, S. M. (2010). Longitudinal changes in cortical thickness associated with normal aging. Neuroimage 52, 1215–1223. doi:10.1016/j.neuroimage.2010.04.258
Tsopelas, C., Stewart, R., Savva, G. M., Brayne, C., Ince, P., Thomas, A., et al. (2011). Neuropathological correlates of late-life depression in older people. Br. J. Psychiatry 198, 109–114. doi:10.1192/bjp.bp.110.078816
Vaillancourt, D. E., Spraker, M. B., Prodoehl, J., Zhou, X. J., and Little, D. M. (2012). Effects of aging on the ventral and dorsal substantia nigra using diffusion tensor imaging. Neurobiol. Aging 33, 35–42. doi:10.1016/j.neurobiolaging.2010.02.006
Van Gelderen, P., De Zwart, J. A., Starewicz, P., Hinks, R. S., and Duyn, J. H. (2007). Real-time shimming to compensate for respiration-induced B0 fluctuations. Magn. Reson. Med. 57, 362–368. doi:10.1002/mrm.21136
van Zijl, P. C., Zhou, J., Mori, N., Payen, J. F., Wilson, D., and Mori, S. (2003). Mechanism of magnetization transfer during on-resonance water saturation. A new approach to detect mobile proteins, peptides, and lipids. Magn. Reson. Med. 49, 440–449. doi:10.1002/mrm.10398
Vita, G., Princi, P., Calabro, R., Toscano, A., Manna, L., and Messina, C. (1986). Cardiovascular reflex tests: assessment of age-adjusted normal range. J. Neurol. Sci. 75, 263–274. doi:10.1016/0022-510X(86)90074-2
Vul, E., Harris, C., Winkielman, P., and Pashler, H. (2009). Puzzlingly high correlations in fMRI studies of emotion, personality, and social cognition. Perspect. Psychol. Sci. 4, 274–290. doi:10.1111/j.1745-6924.2009.01132.x
Walhovd, K. B., Fjell, A. M., Reinvang, I., Lundervold, A., Dale, A. M., Eilertsen, D. E., et al. (2005). Effects of age on volumes of cortex, white matter and subcortical structures. Neurobiol. Aging 26, 1261–1270. doi:10.1016/j.neurobiolaging.2005.05.020
Walhovd, K. B., Westlye, L. T., Amlien, I., Espeseth, T., Reinvang, I., Raz, N., et al. (2011). Consistent neuroanatomical age-related volume differences across multiple samples. Neurobiol. Aging 32, 916–932. doi:10.1016/j.neurobiolaging.2009.05.013
Weiskopf, N., Lutti, A., Helms, G., Novak, M., Ashburner, J., and Hutton, C. (2011). Unified segmentation based correction of R1 brain maps for RF transmit field inhomogeneities (UNICORT). Neuroimage 54, 2116–2124. doi:10.1016/j.neuroimage.2010.10.023
Woodruff-Pak, D. S., Foy, M. R., Akopian, G. G., Lee, K. H., Zach, J., Nguyen, K. P., et al. (2010). Differential effects and rates of normal aging in cerebellum and hippocampus. Proc. Natl. Acad. Sci. U.S.A. 107, 1624–1629. doi:10.1073/pnas.0914207107
Yao, B., Li, T. Q., Gelderen, P. V., Shmueli, K., de Zwart, J. A., and Duyn, J. H. (2009). Susceptibility contrast in high field MRI of human brain as a function of tissue iron content. Neuroimage 44, 1259–1266. doi:10.1016/j.neuroimage.2008.10.029
Zaitsev, M., Dold, C., Sakas, G., Hennig, J., and Speck, O. (2006). Magnetic resonance imaging of freely moving objects: prospective real-time motion correction using an external optical motion tracking system. Neuroimage 31, 1038–1050. doi:10.1016/j.neuroimage.2006.01.039
Zauner, A., and Muizelaar, J. P. (1997). “Chapter 5: brain metabolism and cerebral blood flow,” in Head Injury: Pathophysiology and Management of Severe Closed Injury, eds P. Reilly and R. Bullock (London: Chapman and Hall), 89–99.
Zecca, L., Gallorini, M., Schünemann, V., Trautwein, A. X., Gerlach, M., Riederer, P., et al. (2001). Iron, neuromelanin and ferritin content in the substantia nigra of normal subjects at different ages: consequences for iron storage and neurodegenerative processes. J. Neurochem. 76, 1766–1773. doi:10.1046/j.1471-4159.2001.00186.x
Keywords: brainstem, quantitative MRI, tensor-based morphometry, voxel-based quantification, aging
Citation: Lambert C, Chowdhury R, FitzGerald THB, Fleming SM, Lutti A, Hutton C, Draganski B, Frackowiak R and Ashburner J (2013) Characterizing aging in the human brainstem using quantitative multimodal MRI analysis. Front. Hum. Neurosci. 7:462. doi: 10.3389/fnhum.2013.00462
Received: 02 May 2013; Accepted: 25 July 2013;
Published online: 20 August 2013.
Edited by:Simon Baudrexel, Goethe University Frankfurt, Germany
Reviewed by:Pierre-Louis Bazin, Max Planck Institute for Human Cognitive and Brain Sciences, Germany
Berkin Bilgic, Martinos Center for Biomedical Imaging, USA
Lijuan Zhang, SIAT; Chinese Academy of Sciences, China
Copyright: © 2013 Lambert, Chowdhury, FitzGerald, Fleming, Lutti, Hutton, Draganski, Frackowiak and Ashburner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christian Lambert, St George’s University of London, Cranmer Terrace, SW17 0RE, London, UK e-mail: firstname.lastname@example.org