Impact of locus coeruleus and its projections on memory and aging
Jason Langley, Sana Hussain, Daniel E. Huddleston, Ilana J. Bennett, Xiaoping P. Hu
LLangley, et al.
LC Tractography 1
Impact of locus coeruleus and its projections on memory and aging
Jason Langley , Sana Hussain , Daniel E. Huddleston , Ilana J. Bennett , and Xiaoping P. Hu * Center for Advanced Neuroimaging, University of California Riverside, Riverside, CA 2.
Department of Bioengineering, University of California Riverside, Riverside, CA 3.
Department of Neurology, Emory University, Atlanta, GA 4.
Department of Psychology, University of California Riverside, Riverside, CA * Correspondence to: Xiaoping P. Hu, Ph.D. E-mail: xhu ’at’ engr.ucr.edu
Abstract
Introduction:
Locus coeruleus (LC) is the primary source of norepinephrine to the brain and its efferent projections innervate many brain regions, including the thalamus. LC degrades with normal aging, but not much is known regarding whether its structural connectivity evolves with age or predicts aspects of cognition.
Methods:
Here, we use high-resolution diffusion tensor imaging (DTI)-based tractography to examine structural connectivity between LC and the thalamus in younger and older adults.
Results:
We found LC projections to be bundled in a fiber tract anatomically consistent with the central tegmental tract (CTT) and branch from this tract into the thalamus. The older cohort exhibited a significant reduction in mean and radial diffusivity within CTT as compared to the young cohort. We also observed a significant correlation between CTT mean, axial, and radial diffusivities and memory performance (delayed recall) in the older adult cohort.
Discussion:
These observations suggest that although LC projections degrade with age, the degree of degradation is associated with cognitive abilities in older adults.
Keywords
Locus coeruleus, DTI, aging, central tegmental tract, memory
Significance Statement
Locus coeruleus modulates several cognitive processes including modulates arousal, attention modulation, and memory. Sustaining the integrity of locus coeruleus neurons is hypothesized to play a key role in staving off age-related cognitive decline. However, less is known about how efferent projections of locus coeruleus change with age or cognition. Here, we examine how age effects microstructure of the central tegmental tract, a fiber tract in which locus coeruleus efferent projections are bundled, and whether age-related changes in the microstructure of this tract are associated with cognitive decline. angley, et al.
LC Tractography 2
1. Introduction
Locus coeruleus (LC) is a small cylindrical structure located along the dorsal edge of the brainstem anterior to the 4th ventricle that is filled with catecholaminergic neurons. As the primary source of norepinephrine to much of the brain, it innervates many areas including the neocortex, thalamus, and hippocampus (Aston-Jones and Cohen, 2005). The release of norepinephrine by neurons in LC is hypothesized to reduce inflammation and protect neurons from oxidative stress (Feinstein, et al., 2002; Robertson, 2013; Troadec, et al., 2001), with the density of LC neurons correlating with cognitive performance in older adults (Kelly, et al., 2017; Wilson, et al., 2013). As such, LC is thought to play a pivotal role in forestalling cognitive decline, such as what occurs in normal aging. Neuromelanin-sensitive contrast, associated with incidental (Sasaki, et al., 2006) or explicit (Chen, et al., 2014; Huddleston, et al., 2018; Schwarz, et al., 2011) magnetization transfer effects, has been found to colocalize with melanized neurons in LC (Keren, et al., 2009; Keren, et al., 2015). LC regions of interest (ROIs) derived from magnetization transfer contrast have been employed to examine how LC modulates attention (Krebs, et al., 2018; Murphy, et al., 2014), memory encoding (Clewett, et al., 2014; Clewett, et al., 2018; Hammerer, et al., 2018), or examine resting state functional connectivity between LC and memory networks (Jacobs, et al., 2015; Jacobs, et al., 2018). Other studies have found that structural measures within LC, such as magnetization transfer contrast (Clewett, et al., 2016; Dahl, et al., 2019) or diffusivity (Langley, et al., 2020b), correlate with memory performance in older adults. Together, these studies suggest that MRI can be used as a proxy for assessing norepinephrine function (Cassidy, et al., 2019) and its contribution to cognition in adults across the lifespan. Histological studies have found a rostral-caudal gradient of age-related cell loss in LC, with rostral portions of LC exhibiting the greatest neuronal loss (Chan-Palay and Asan, 1989; Manaye, et al., 1995). This is consistent with imaging studies that have found reduced neuromelanin-sensitive contrast in rostral LC ROIs of older adults (Betts, et al., 2017; Dahl, et al., 2019; Liu, et al., 2018). Using retrograde tracers in macaques, noradrenergic fibers arising from the rostral portions of LC have been shown to enter the dorsomedial area of the central tegmental tract (CTT) and a group of fibers branch out to the dorsal thalamus. The remaining fibers in the noradrenergic bundle move adjacent to the ventral surface of the thalamus, merge into the medial forebrain bundle, and branch through the substantia innominata toward the amygdala (Bowden, et al., 1978; Tanaka, et al., 1982). Age-related degradation of the noradrenergic bundle may account for decreased catecholaminergic innervation observed in several brain regions by tracer studies in older primates (Goldman-Rakic and Brown, 1981). However, less is known about how LC projections in the noradrenergic bundle vary with age or relate to memory in humans. Diffusion tensor imaging (DTI) tractography can be used to reconstruct fiber tracts in vivo (Mori and van Zijl, 2002) and their integrity can be examined using measures such as fractional anisotropy (FA), which gives an estimate of restricted diffusivity in the fiber bundle. In this study, we build on prior work examining age-related differences in LC microstructure and explore how age affects microstructure of the CTT using a high resolution DTI protocol. Our focus here is on the CCT given that this first segment of the noradrenergic bundle has defined endpoints (LC, thalamus) and is well established in histology. We further assess the relationship between white matter integrity of the CTT and memory performance. Finally, we investigate the rostral-caudal selectivity of age-related gray matter microstructural differences in LC and their effect on memory.
2. Materials and Methods
Sixty-six participants (30 older and 36 younger participants) were recruited for this study from the Riverside, California community and the student population at University of California, Riverside, respectively. All participants in the study gave written informed consent in accordance with local institutional review board (IRB) regulations. Prior to enrollment in the study, potential participants were angley, et al.
LC Tractography 3 screened over the phone for normal global cognition as defined by a score of >17 on a subset of the Montreal Cognitive Assessment (MOCA) adapted for phone screening (Pendlebury, et al., 2013). They were also screened for contraindications to MRI imaging and diagnosed neuropsychiatric condition, such as depression or stroke, that could influence their performance. One younger subject was excluded from the diffusion-weighted MRI analysis due to problems with data acquisition and three participants (2 older and 1 younger) were excluded from the analyses because of significant motion artifacts.
The final sample size was 28 older and 35 younger participants. Demographic information (gender, age), MOCA, and Rey Auditory Verbal Learning Test (RAVLT) scores were collected on each subject. Group means for age, MOCA, and RAVLT scores are given in Table 1.
Table 1.
Demographic data for the final sample are presented as mean ± standard deviation. MOCA - Montreal cognitive assessment scoring; RAVLT - Rey Auditory Verbal Learning Test. Variable Younger (n=35) Older (n=28) p Value Gender (M/F) 15/20 13/15 0.19 Age (years) 20.7 ± ± -4 MOCA score 27.5 ± ± ± ± ± ± -4 RAVLT Delay 11.1 ± ± Imaging data were acquired on a 3 T MRI scanner (Prisma, Siemens Healthineers, Malvern, PA). Excitation was performed using the body coil on the scanner and signal was received using a 32-channel receive only coil at the Center for Advanced Neuroimaging at University of California, Riverside. Anatomic images were acquired with an MP-RAGE sequence (echo time (TE)/repetition time (TR)/inversion time=3.02/2600/800 ms, GRAPPA acceleration factor=2, flip angle=8°, voxel size=0.8 × × ) for registration from subject space to common space. High-resolution diffusion-weighted MRI data were collected with a diffusion-weighted single-shot spin-echo, echo planar imaging sequence with the following parameters: TE / TR = 78 / 4500 ms, FOV = 194 ×
168 mm , matrix size of 204 × × × , multiband factor = 2, and 64 slices with no gap, covering the brain from the middle of the cerebellum to the striatum. Monopolar diffusion-encoding gradients were applied in 30 directions with b values of 500 s/mm and 2000 s/mm . Two sets of b =0 images were acquired, with the two sets having opposite polarities of phase-encoding direction for the correction of susceptibility distortion (Andersson, et al., 2003). Imaging data were analyzed with FMRIB Software Library (FSL). A transformation was derived between individual subject space to Montreal Neurological Institute (MNI) 152 T -weighted space using FMRIB’s Linear Image Registration Tool (FLIRT) and FMRIB’s Nonlinear Image Registration Tool (FNIRT) in the FSL software package using the following steps (Smith, et al., 2004; Woolrich, et al., 2009). (1) The T -weighted image was skull stripped using the brain extraction tool (BET) in FSL, (2) brain extracted T -weighted images were aligned with the MNI brain extracted image using an affine transformation, and (3) a nonlinear transformation (FNIRT) was used to generate a transformation from individual T -weighted images to T -weighted MNI152 common space. Diffusion-weighted data were analyzed with FSL (Jenkinson, et al., 2002; Jenkinson and Smith, 2001; Smith, et al., 2004) and MATLAB (The Mathworks, Natick, MA). Standard preprocessing steps were applied to correct susceptibility induced distortions in the diffusion MR data. Diffusion MR data angley, et al.
LC Tractography 4 were first corrected for eddy-current distortion, motion, and for susceptibility distortion using eddy in FSL (Andersson, et al., 2003; Andersson and Sotiropoulos, 2016). Next, skull stripping of the T -weighted image and susceptibility corrected b =0 image was performed using the brain extraction tool in the FSL software package (Smith, 2002). Measures derived from the diffusion MR data, including FA, mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) were calculated using the dtifit tool in FSL. A bi-exponential model of water diffusion was used to estimate the fractional contribution of freely diffusing (i.e., unrestricted) water to the overall signal per voxel (Hoy, et al., 2014; Pasternak, et al., 2009). The free water elimination algorithm was implemented in DiPy (Garyfallidis, et al., 2014). Diffusion-weighted data from each subject was assessed for gross head motion and ghosting artifacts from pulsatile motion. Three subjects were excluded from the analysis due to significant motion artifacts and had root-mean-square (RMS) translation estimate between consecutive volumes greater than 4 voxels. For the remaining subjects used in the analysis, the maximal RMS translation estimate between consecutive volumes was 1.53 voxels (standard deviation: 0.74 voxels). Registration between diffusion-weighted and T images was derived using the b =0 image with a rigid body transformation with boundary based registration cost function. The quality of each registration was assessed for each subject and no significant misregistrations were observed across the sample. Prior to fiber tracking, principal diffusion directions and crossing fibers were estimated using dtifit and bedpostX. Next, brain extraction of the b=
0 image was performed using BET (Smith, 2002). Segmentation of grey matter (GM), white matter (WM), cerebrospinal fluid and their partial volume estimates were created from the anatomical scans using FAST (Zhang, et al., 2001). Probabilistic tractography, as implemented in FSL, was used to track the fiber tracts originating from LC to thalamus. A LC atlas was derived from magnetization transfer images in 31 healthy participants (mean age: 26.1 years) using processes similar to those described previously (Langley, et al., 2020b). This atlas was used to define a LC seed mask for tractography separately for each hemisphere. A bilateral thalamus atlas from the Harvard-Oxford atlas (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/atlases) was selected as the target mask. The procedures used to transform bilateral seed and target masks in diffusion space were as follows: Each mask was transformed to each subject’s DTI-space using the nonlinear transformation described above concatenated with a rigid-body transformation using applywarp. The boundary between gray matter GM in the thalamus was found using the following procedure. First, the thalamus and WM partial volume estimates were transformed to DTI space. Next, the partial volume estimates of WM were thresholded at a level of 0.5, dilated, and binarized. Finally, the overlap of the thalamus mask and processed WM maps were defined as the thalamus and WM boundary mask. For each subject, probabilistic tractography was carried out from all voxels in each seed mask for both hemispheres using a model accommodating two fiber orientations (Behrens, et al., 2007) in DTI space. Probability density functions were calculated by drawing 5000 samples per voxel in the LC seed mask from the uncertainty distribution of the fiber orientations between the LC seed mask and the thalamus and WM boundary mask. The number of these samples that reached the thalamus and WM target mask was counted as the waytotal. After each subject’s LC-thalamus tracts were tracked, the pathways were normalized by the waytotal number, and the normalized value in each voxel was interpreted as the tract probability given the prior knowledge that the tract exists in that voxel. Then, each subject’s tract probabilities were transformed to MNI space using a nonlinear transformation. Next, all transformed tracts in MNI space were averaged and multiplied by 100 to generate a percentile track probability for LC-thalamus tract in older and young populations. The process for creating the tract atlas is summarized in Figure 1. Finally, the young LC-thalamus tract was then transformed back to native diffusion space and mean FA, MD, AD, RD, and free water were measured in the resulting bilateral LC-thalamus population tract for each subject. angley, et al.
LC Tractography 5
Figure 1.
Schematic describing the creation of the CTT population atlas.
For each younger adult participant, normalized LC-thalamus pathways were transformed to MNI space, thresholded, and binarized. Next, the binarized tracts were summed and normalized by the total number of young participants. the resulting value in each voxel of the population atlas represents the overlap across all participants.
The LC atlas was divided into rostral and caudal ROIs with each ROI having a volume of 50% of the total volume of the LC atlas. Each ROI was transformed to diffusion space, thresholded, and binarized in an identical procedure to that of the LC seed mask described above. The mean volume of the bilateral rostral and caudal ROIs in diffusion space was each approximately 30 voxels. Mean FA, MD, RD, and AD were measured in caudal and rostral LC ROIs for each subject.
A nigrostriatal tract atlas was created to test the specificity of the relationship between CTT diffusivity and memory performance as well as diffusivity of LC subregions. A substantia nigra pars compacta atlas was derived from magnetization transfer images in 31 healthy participants (mean age: 26.1 years) using processes similar to those described previously (Langley, et al., 2020a; Langley, et al., 2017) was used to define the nigral seed mask. A bilateral putamen atlas from the Harvard-Oxford atlas (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/atlases) was selected as the target mask and tractography was performed in the younger adult group as described in section 2.6. The nigrostriatal tract atlas was transformed to native diffusion space, thresholded at a level of 0.6, and binarized. For each subject, mean FA, MD, RD, and AD were measured in the nigrostriatal atlas. Results for the nigrostriatal analysis are summarized in the Supplemental Materials section.
All statistical analyses were performed using IBM SPSS Statistics software version 24 (IBM Corporation, Somers, NY, USA) and results are reported as mean ± standard deviation. A p value of 0.05 was considered significant for all statistical tests performed in this work. Normality of diffusion data was assessed using the Shapiro-Wilk test for each group and all data was found to be normal. We hypothesized that degradation of CTT, as measured by DTI indices, is correlated with cognitive measures in the older adult group but not in the younger adult group. Thus, to assess the impact of age-related decline in this tract, we performed a Spearman rank correlation of mean CTT DTI indices with cognitive measures (RAVLT delayed recall scores) separately in younger and older adult groups. Age group differences in the magnitude of these relationships were assessed using a one-tailed Fisher’s z-test for independent correlations. We further hypothesized that the rostral LC subregion has stronger angley, et al. LC Tractography 6 correlations with cognitive measures than the caudal LC ROI in older adults. This hypothesis was tested by performing a Spearman rank correlation of RAVLT delayed recall with rostral or caudal LC subregions separately in the older adult group. LC subregion ROI differences in the magnitude of these relationships were assessed using a one-tailed Fisher’s z-test for independent correlations. The effects of age group (younger, older) and LC subregion (caudal, rostral) were tested with separate two-way analyses of variance (ANOVA) for each diffusion index (MD, AD, RD). If the interaction was significant, post hoc comparisons between the age groups were performed for each LC subregion using respective two-tailed t -tests. Age group comparisons of the LC-thalamus tract (CTT) were made using separate two-tailed t- tests for each diffusion parameter. Finally, to examine the relationship between CTT microstructure and microstructure of LC subregions, Spearman rank correlations were performed on each CTT DTI index to the corresponding LC subregion DTI index (i.e. CTT MD was correlated with rostral LC MD). LC subregion differences in the magnitude of these relationships were assessed using a one-tailed Fisher’s z-test for independent correlations.
3. Results
The existence and variability of fiber tracts between LC and thalamus across individuals were assessed by measuring the maximum intensity value at any given voxel. The intensity value at each voxel in Figure 2 represents the amount of overlap across individuals for that voxel since the tracts for each individual were binarized and then summed to develop group maps for each tract. The maximum probability of the fiber tract in each group was 1. The variability of fiber tracts between LC and thalamus was further probed by examining the total volume of the tract with intensity values greater than 0.75. Approximately 53% of the tract volume had a probability value of 0.75 or greater in the younger group. Both results indicate significant overlap and consistency between individual tractography results. Importantly, tracts leaving the LC seed mask were contained within a region anatomically consistent with the CTT and projecting to the dorsal portion of the thalamus (see Figure 2).
Figure 2.
Sagittal (A; top row) and axial (A; bottom row) of population maps showing the tract from LC to thalamus in the young group. The lines in sagittal view at X=5 mm identify the location of the axial slices shown in the bottom row. The population map was generated by taking individual tracts, thresholding at a level of 0.15 (corresponding to 15% of the total streamlines sent from LC), binarizing, and transforming to standard space. The tract is anatomically consistent with the CTT. A three dimensional rendering of the CTT is shown in B. angley, et al.
LC Tractography 7
Relative to young adults, higher FA (older: 0.31 ± ± t =-2.6, p= × -4 mm /s ± × -5 mm /s; young: 3.31 × -4 mm /s ± × -5 mm /s; t =2.4, p= × -4 mm /s ± × -5 mm /s; young: 3.31 × -4 mm /s ± × -5 mm /s; t =3.0, p= × -4 mm /s ± × -5 mm /s; young: 4.31 × -4 mm /s ± × -5 mm /s; t =1.3, p =0.197) or free water (older: 0.34 ± ± t= p =0.158) was seen between cohorts. These results are summarized in Figure 3. Figure 3 . Group comparisons between younger and older adult groups in CTT MD (A), RD (B), and AD (C). Statistically significant decreases in mean tract MD and mean tract RD were observed in the older group as compared to the younger group. Correlations between delayed recall and CTT MD (D), RD (E), and AD (F) were observed in the older group.
As shown in Figure 3, significant correlations were seen between diffusion metrics in the CTT and RAVLT delayed recall in the older group (mean tract MD: r =0.477, p =0.005; mean tract RD: r =0.384, p =0.022; mean tract AD: r =0.459, p =0.007) but not in the young group (mean tract MD: r =-0.268, p =0.104; mean tract RD: r =-0.309, p =0.06; mean tract AD: r =-0.173, p =0.282), with higher diffusivity relating to better memory performance. Fisher’s z tests revealed that these relationships were significantly weaker in the young compared to older group for all diffusion measures (mean tract MD: z =-2.79, p =0.003; mean tract RD: z =-2.71, p =0.007; mean tract AD: z =-2.51, p =0.012). No correlation was observed between CTT free water and RAVLT delayed recall in the older group ( r =-0.098; p =0.347) but a significant correlation was observed between CTT free water and RAVLT delayed recall in the younger group ( r =-0.392; p = 0.032). Controlling for age, education, and gender did not substantially change the correlations between CTT diffusion metrics and RAVLT delayed recall in the older adult group (mean tract MD: r =0.373, p =0.028; mean tract RD: r =0.296, p =0.046; mean tract AD: r =0.428, p =0.013) or in younger adult group (mean tract MD: r =-0.267, p =0.111; mean tract RD: r =-0.310, p =0.062; mean tract AD: r =-0.175, p =0.299). The effect of age and LC subregion on diffusivity was assessed using separate 2 Age Group x 2 LC subregion ANOVAs for each diffusion parameter. Results revealed significant main effects of LC subregion for all diffusion measures (MD: F =60.9, p <10 -4 ; RD: F =61.3, p <10 -4 ; AD: F =70.9, p <10 -4 ), with increased diffusivity in the rostral versus caudal LC subregion for MD (rostral: 3.89 × -4 mm /s ± × -5 mm /s; caudal: 3.50 × -4 mm /s ± × -5 mm /s), RD (rostral: 3.35 × -4 mm /s ± × -5 mm /s; caudal: 2.98 × -4 mm /s ± × -5 mm /s), and AD (rostral: 4.99 × -4 mm /s ± × -5 mm /s; caudal: 4.51 × -4 mm /s ± × -5 mm /s). Results for the comparison between rostral and caudal LC subregions are summarized in Figure 4. Main effects of age group were significant for RD ( F =4.8, angley, et al. LC Tractography 8 p <0.032), trending for MD ( F =3.5; p =0.066), and not significant for AD ( F =0.001; p =0.98), with increased diffusivity in younger versus older adults for RD (younger: 3.23 × -4 mm /s ± × -5 mm /s; older: 3.03 × -4 mm /s ± × -5 mm /s) and MD (younger: 3.76 × -4 mm /s ± × -5 mm /s; older: 3.57 × -4 mm /s ± × -5 mm /s). There were no significant interactions between age group and LC subregion for any diffusion measure (MD: F =0.001, p= F =0.102, p= F =0.119, p =0.731). Figure 4.
Group comparisons from the LC subregion analysis. Group comparisons for MD, RD, and AD in the rostral and caudal LC subregions from the older adult cohort are shown in A, B, and C, respectively. Comparisons for MD, RD, and AD in the rostral and caudal LC subregions from the younger adult cohort are shown in D, E, and F respectively. Significant decreases in MD, RD, and AD of the caudal LC subregion are seen in both cohorts.
LC integrity relates to memory performance.
For rostral LC, significant positive correlations were found between the RAVLT delayed recall score and MD ( r =0.645; p <10 -4 ), RD ( r =0.614; p= × -4 ), and AD ( r =0.626; p= × -4 ) in the older group, with higher diffusivity being associated with better memory performance. In the young group, the RAVLT delayed recall score did not significantly correlate with MD ( r =-0.141; p =0.199), RD ( r =-0.173; p =0.150), or AD ( r =-0.117; p =0.242). Fisher’s z tests revealed that these relationships were significantly smaller in the young compared to older group for all diffusion measures (mean tract MD: z =-3.42, p =0.0001; mean tract RD: z =-3.50, p =0.0005; mean tract AD: z =-3.19, p =0.0014). These correlations for rostral LC are shown in Figure 5. Figure 5.
Rostral LC diffusivity and memory performance. The top row contains correlations between RAVLT delayed recall and MD (A), RD (B), and AD (C) in the older adult group. The bottom row contains correlations between RAVLT delayed recall and MD (D), RD (E), and AD (F) in rostral LC subregion of the younger adult group. Significant correlations are seen between rostral LC diffusivity measures and RAVLT delayed recall score in the older adult group but not in the younger adult group. angley, et al.
LC Tractography 9
For caudal LC, significant positive correlations were also found between the RAVLT delayed recall score and MD ( r =0.606; p= × -4 ), RD ( r =0.594; p= × -4 ), and AD ( r =0.643; p< -4 ) in the older group. However, no significant correlations were observed between the RAVLT delayed recall score and any diffusion measure (MD: r =-0.009, p =0.478; RD: r =0.078, p =0.466; AD: r =-0.010, p =0.486) in the younger group. Fisher’s z tests revealed that these relationships were significantly smaller in the young compared to older group for all diffusion measures (mean tract MD: z =-2.67, p =0.0076; mean tract RD: z =-2.27, p =0.0232; mean tract AD: z =-2.51, p =0.0036). These correlations for caudal LC are shown in Figure 6. Figure 6.
Caudal LC diffusivity and memory performance.
The top row contains correlations between RAVLT delayed recall and MD (A), RD (B), and AD (C) in the older adult group. The bottom row contains correlations between RAVLT delayed recall and MD (D), RD (E), and AD (F) in caudal LC subregion of the younger adult group. Significant correlations are seen between caudal LC diffusivity measures and RAVLT delayed recall score in the older adult group but not in the younger adult group.
LC subregion diffusivity relates to CTT diffusivity.
For diffusivity indices in rostral LC, significant positive correlations were found between diffusivity indices in CTT (MD: r =0.656, p= × -4 ; RD: r =0.507, p= r =0.679, p= × -5 ) in the older group, with higher diffusivity in rostral LC being associated with higher diffusivity in CTT. However, no correlation was observed between caudal LC MD and CTT MD ( r =0.292, p= r =0.165, p= r =0.444, p= r =0.787, p =4.7 × -9 ; RD: r =0.669, p =6.1 × -10 ; AD: r =0.822, p =1.5 × -17 ) and caudal (MD: r =0.492, p =0.002; RD: r =0.351, p =0.004; AD: r =0.561, p =8.1 × -7 ) LC subregions. Fisher’s z tests revealed that these relationships were significantly smaller in the caudal LC compared to the rostral LC for all diffusivities in the older (MD: z=- p =0.006; RD: z=- p =0.020; AD: z= -2.0, p =0.049) and younger group (MD: z=- p =0.007; RD: z=- p =0.023; AD: z= -2.7, p =0.009). Correlations for the older cohort are shown in Figure 7.
4. Discussion
In furthering our current understanding of the role LC plays in cognition and age-related cognitive decline, it would be advantageous to map the white matter fiber tracts connecting LC with the rest of the brain. To this end, we focused on mapping LC efferent projections in the CTT, first segment of the noradrenergic bundle, since this connection has been well mapped by invasive tracer studies of non-human primates (Bowden, et al., 1978; Samuels and Szabadi, 2008; Tanaka, et al., 1982). After careful inspection of Duvernoy’s Atlas of the Human Brain Stem and Cerebellum (Naidich, et al., 2009), we found LC efferent projections to be bundled in a fiber tract anatomically consistent with the CTT and branched from this tract into the thalamus. Age-related reductions in RD and increases in FA were observed in this tract and we found significant correlations between CTT diffusivity measures and angley, et al.
LC Tractography 10
Figure 7.
Relationship between LC subregion diffusivity and CTT diffusivity in the older adult group.
The top row contains correlations between rostral LC diffusivity indices and CTT diffusivity indices. The bottom row contains correlations between caudal LC diffusivity indices and CTT diffusivity indices. Significant correlations are seen between all rostral diffusivity indices and CTT diffusivity indices while only AD in the caudal LC correlates with CTT AD.
RAVLT delayed recall in older, but not younger, adults. A complementary analysis of LC subregions found significant age-related microstructural differences in the rostral and caudal portions of LC (reduced MD and RD in older adults), with microstructure in both subregions correlating with RAVLT delayed recall in older, but not younger, adults. For both age groups, higher LC diffusivity was associated with higher CTT diffusivity in rostral, compared to caudal, subregions. Of note, none of these effects were observed in a control tract, the NST (see Supplemental Materials). Taken together, these findings suggest that axial connectivity and axonal size in CTT white matter play a role in the retention of cognitive abilities in aging. White matter fiber tracts in CTT were identified by tractography from a seed region in LC to target region in the thalamus. Individuals demonstrated a high degree of overlap between tracts connecting LC and thalamus in standard space as indicated by the maximum overlap intensity and percent of tract volume shared by at least 75% of individuals in the younger adult group. This overlap measure is more accurately seen as a probability of tract existence. In particular, this type of probability measure reflects the anatomical variability that occurs across individuals and should not be thought of as a measure of connectivity strength (Hua, et al., 2009). In addition to finding that LC efferent projections are bundled in the CTT (Bowden, et al., 1978; Tanaka, et al., 1982), we observed higher FA and lower RD in this tract in the older adult group. Efferent projections of LC are spatially organized with axons in the rostral portion projecting to the thalamus, hypothalamus, cortex, and forebrain while axons in the caudal portion project into the cerebellum and spinal cord (Bowden, et al., 1978; Tanaka, et al., 1982; Westlund, et al., 1984). Our analysis examining the relationship between CTT microstructural integrity and microstructural integrity of LC subregions found CTT microstructural integrity more closely related to microstructural integrity in the rostral LC subregion than in the caudal LC subregion in both cohorts. In particular, in the older adult cohort, which was expected to experience some age-related neuronal loss, we found significant correlations between CTT diffusivity and diffusivity in the rostral LC subregion but no relationship was observed between diffusivity in the caudal LC subregion and CTT diffusivity. These results accord with the aforementioned tracer studies (Bowden, et al., 1978; Tanaka, et al., 1982; Westlund, et al., 1984) and indicate that CTT diffusivity may be a proxy for LC integrity when scan resolution prohibits accurate LC segmentation in older adults. The shape and density of cells changes along the rostral-caudal axis of LC. Cells are less dense and circularly shaped in the rostral portion of LC while the shape becomes less circular and density of cells increases in the caudal portions of LC (Chan-Palay and Asan, 1989; Manaye, et al., 1995). High neuronal density in the caudal LC should affect LC microstructure by increasing barriers to diffusion in that region (Beaulieu, 2002), which is expected to lead to decreases in diffusivity. Interestingly, rostral angley, et al.
LC Tractography 11 and caudal portions of LC in both cohorts had different microstructural properties with the caudal LC ROI exhibiting decreased diffusivity as compared to the rostral LC ROI. This result may reflect the differences in density of cells along the rostral-caudal axis of LC. Age-related microstructural differences were observed in both rostral and caudal subregions of LC, with older adults showing reduced mean and radial diffusivity compared to younger adults. In contrast to our predictions, however, these age effects were not regionally selective. This observation is not in agreement with earlier neuromelanin-sensitive imaging studies which found an age-related reduction in neuromelanin sensitive contrast in the rostral portions of LC (Betts, et al., 2017; Dahl, et al., 2019; Liu, et al., 2018). The lack of regional selectivity with diffusion measures is not in accordance with several histology studies that found the rostral LC to exhibit extensive age-related neuronal loss (Chan-Palay and Asan, 1989; Manaye, et al., 1995), suggesting that these diffusivity measures may not be sensitive enough to detect subtle neuronal loss in this region. However, results from histology are not in complete agreement with several studies finding no age-related neuronal loss in the rostral portions of LC (Mouton et al., 1994; Ohm et al., 1997). These studies have speculated that the neuronal loss seen in other histology may be due to Alzheimer’s disease since the rostral portion of LC is selectively affected by Alzheimer’s disease (German et al., 1992; Hoogendijk et al., 1995). Significant correlations were observed between reduced AD and RD in this tract and worse performance on RAVLT delayed recall in the older adult group. Reduced RD may indicate a winnowing of axons in CTT (Barazany, et al., 2009; Wheeler-Kingshott and Cercignani, 2009). We further found that lower performance on RAVLT delayed recall scores in the older adult group was correlated with reduced diffusivity in both rostral and caudal LC subregion ROIs. Interestingly, the correlations between RAVLT delayed recall and diffusivity in rostral and caudal LC were of similar strength and may suggest that there is no selectivity in participation of LC neurons in memory. This interpretation accords with both histology (Kelly, et al., 2017; Wilson, et al., 2013) and imaging (Clewett, et al., 2016) in older adults which found a reduction in total neuron count or contrast from melanized neurons in LC correlates with decreased cognitive performance. In contrast, only a correlation between mean CTT free water and RAVLT delayed recall was seen in the younger cohort, with lower RAVLT delayed recall scores associated with higher CTT free water. Clinical studies assessing free water in white matter tracts have linked the free water measure to cellular density, atrophy, or neuroinflammation (Pasternak, et al., 2014; Pasternak, et al., 2012). Lower performance on RAVLT delayed recall is likely driven by individual differences in density of CTT as it is unlikely that there is significant atrophy or neuroinflammation in the CTT of a population of young control subjects. No correlation was observed between RAVLT delayed recall and caudal or rostral LC diffusivity measures in the younger group. This result stands in contrast to a recently published imaging study that found higher neuromelanin-sensitive contrast in medial and caudal LC ROIs to be correlated with improved performance on an emotional memory task in younger adults (Jacobs, et al., 2020). RAVLT delayed recall is sensitive measure for diagnosis of cognitive impairment in older adults (Zhao, et al., 2015). However, younger adults, on average, perform better on RAVLT delayed recall than older adults (Messinis, et al., 2016) and younger adults tend to have similar performance on RAVLT delayed recall (Badcock, et al., 2011). Our data accords with these studies and the lack of correlation between diffusivity in LC subregions with RAVLT delayed recall may be due to the inability of RAVLT delayed recall to detect subtle differences in normal cognition in younger adults. There are some caveats in the present study. First, neuromelanin-sensitive data were not acquired in the same cohort so it is impossible to examine microstructural changes in rostral and caudal LC subregion ROIs alongside contrast from melanized neurons. Second, we only considered fiber tracts verified by anatomical tracer studies in non-human primate and limited our tractography analysis to the CTT. Referencing anatomical tracer studies in non-human primates is necessary to interpret tractography angley, et al.
LC Tractography 12 results given the complexity and structural proximity of multiple pathways in the brainstem and ensure that no false tracts were identified. Finally, it is impossible to exclude the possibility that some of the healthy controls in the older adult cohort have an undiagnosed neurological condition which may affect the microstructure of LC or CTT.
5. Conclusion
In this work, a high-resolution diffusion-weighted MRI protocol was employed to examine age-related microstructural and compositional differences in LC subregions as well as in CTT. Older age was most associated with axonal thinning in rostral portions of LC and CTT. Strong correlations were observed between diffusivity and RAVLT delayed recall score in CTT and both subregions of LC. In particular, reduced diffusivity was correlated with worse memory performance in older adults and may suggest that connectivity and axonal size in white matter between LC and thalamus play a role in the retention of cognitive abilities.
Author Contributions
JL contributed to the conception and design of the study, data analysis, conducted statistical analyses, and drafted the manuscript. SH contributed in data and statistical analyses. DH contributed to the interpretation of data. IJB contributed to the conception and design of the study, conducted statistical analyses, interpretation of data, and acquired financial support for the study. XH contributed to the conception, design of the study, and interpretation of data. All authors contributed to revisions of the manuscript.
Funding
This work was supported by R00 AG047334 (Bennett) and R21 AG054804 (Bennett) from the National Institutes of Health/National Institute on Aging, 1K23NS105944-01A1 (Huddleston) from the National Institutes of Health/National Institute of Neurological Diseases and Stroke, and the American Parkinson’s Disease Foundation Center for Advanced Research at Emory University (Huddleston).
Disclosures
The authors have no conflicts of interest to disclose
References
Andersson, J.L., Skare, S., Ashburner, J. (2003) How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. NeuroImage, 20:870-88. Andersson, J.L.R., Sotiropoulos, S.N. (2016) An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage, 125:1063-1078. Aston-Jones, G., Cohen, J.D. (2005) Adaptive gain and the role of the locus coeruleus–norepinephrine system in optimal performance. J Comp Neurol, 493:99-110. Badcock, J.C., Dragovic, M., Dawson, L., Jones, R. (2011) Normative data for Rey's Auditory Verbal Learning Test in individuals with schizophrenia. Arch Clin Neuropsychol, 26:205-13. Barazany, D., Basser, P.J., Assaf, Y. (2009) In vivo measurement of axon diameter distribution in the corpus callosum of rat brain. Brain, 132:1210-20. angley, et al.
LC Tractography 13 Beaulieu, C. (2002) The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed, 15:435-55. Behrens, T.E., Berg, H.J., Jbabdi, S., Rushworth, M.F., Woolrich, M.W. (2007) Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? NeuroImage, 34:144-55. Betts, M.J., Cardenas-Blanco, A., Kanowski, M., Jessen, F., Duzel, E. (2017) In vivo MRI assessment of the human locus coeruleus along its rostrocaudal extent in young and older adults. Neuroimage, 163:150-159. Bowden, D.M., German, D.C., Poynter, W.D. (1978) An autoradiographic, semistereotaxic mapping of major projections from locus coeruleus and adjacent nuclei in Macaca mulatta. Brain Res, 145:257-76. Cassidy, C.M., Zucca, F.A., Girgis, R.R., Baker, S.C., Weinstein, J.J., Sharp, M.E., Bellei, C., Valmadre, A., Vanegas, N., Kegeles, L.S., Brucato, G., Kang, U.J., Sulzer, D., Zecca, L., Abi-Dargham, A., Horga, G. (2019) Neuromelanin-sensitive MRI as a noninvasive proxy measure of dopamine function in the human brain. Proc Natl Acad Sci U S A, 116:5108-5117. Chan-Palay, V., Asan, E. (1989) Quantitation of catecholamine neurons in the locus coeruleus in human brains of normal young and older adults and in depression. J Comp Neurol, 287:357-72. Chen, X., Huddleston, D.E., Langley, J., Ahn, S., Barnum, C.J., Factor, S.A., Levey, A.I., Hu, X. (2014) Simultaneous imaging of locus coeruleus and substantia nigra with a quantitative neuromelanin MRI approach. Magn Reson Imaging, 32:1301-6. Clewett, D., Schoeke, A., Mather, M. (2014) Locus coeruleus neuromodulation of memories encoded during negative or unexpected action outcomes. Neurobiol Learn Mem, 111:65-70. Clewett, D.V., Huang, R., Velasco, R., Lee, T.H., Mather, M. (2018) Locus Coeruleus Activity Strengthens Prioritized Memories Under Arousal. J Neurosci, 38:1558-1574. Clewett, D.V., Lee, T.H., Greening, S., Ponzio, A., Margalit, E., Mather, M. (2016) Neuromelanin marks the spot: identifying a locus coeruleus biomarker of cognitive reserve in healthy aging. Neurobiol Aging, 37:117-26. Dahl, M.J., Mather, M., Duzel, S., Bodammer, N.C., Lindenberger, U., Kuhn, S., Werkle-Bergner, M. (2019) Rostral locus coeruleus integrity is associated with better memory performance in older adults. Nat Hum Behav, 3:1203-1214. Feinstein, D.L., Heneka, M.T., Gavrilyuk, V., Dello Russo, C., Weinberg, G., Galea, E. (2002) Noradrenergic regulation of inflammatory gene expression in brain. Neurochem Int, 41:357-65. Garyfallidis, E., Brett, M., Amirbekian, B., Rokem, A., van der Walt, S., Descoteaux, M., Nimmo-Smith, I., Dipy, C. (2014) Dipy, a library for the analysis of diffusion MRI data. Front Neuroinform, 8:8. Goldman-Rakic, P.S., Brown, R.M. (1981) Regional changes of monoamines in cerebral cortex and subcortical structures of aging rhesus monkeys. Neuroscience, 6:177-87. Hammerer, D., Callaghan, M.F., Hopkins, A., Kosciessa, J., Betts, M., Cardenas-Blanco, A., Kanowski, M., Weiskopf, N., Dayan, P., Dolan, R.J., Duzel, E. (2018) Locus coeruleus integrity in old age is selectively related to memories linked with salient negative events. Proc Natl Acad Sci U S A, 115:2228-2233. Hoy, A.R., Koay, C.G., Kecskemeti, S.R., Alexander, A.L. (2014) Optimization of a free water elimination two-compartment model for diffusion tensor imaging. Neuroimage, 103:323-333. Hua, K., Oishi, K., Zhang, J., Wakana, S., Yoshioka, T., Zhang, W., Akhter, K.D., Li, X., Huang, H., Jiang, H., van Zijl, P., Mori, S. (2009) Mapping of functional areas in the human cortex based on connectivity through association fibers. Cereb Cortex, 19:1889-95. Huddleston, D.E., Langley, J., Dusek, P., He, N., Faraco, C.C., Crosson, B., Factor, S., Hu, X.P. (2018) Imaging parkinsonian pathology in substantia nigra with MRI. Curr Radiol Rep, 6:15. Jacobs, H.I., Priovoulos, N., Poser, B.A., Pagen, L.H., Ivanov, D., Verhey, F.R., Uludag, K. (2020) Dynamic behavior of the locus coeruleus during arousal-related memory processing in a multi-modal 7T fMRI paradigm. Elife, 9. angley, et al.
LC Tractography 14 Jacobs, H.I., Wiese, S., van de Ven, V., Gronenschild, E.H., Verhey, F.R., Matthews, P.M. (2015) Relevance of parahippocampal-locus coeruleus connectivity to memory in early dementia. Neurobiol Aging, 36:618-26. Jacobs, H.I.L., Muller-Ehrenberg, L., Priovoulos, N., Roebroeck, A. (2018) Curvilinear locus coeruleus functional connectivity trajectories over the adult lifespan: a 7T MRI study. Neurobiol Aging, 69:167-176. Jenkinson, M., Bannister, P., Brady, M., Smith, S. (2002) Improved optimization for the robust and accurate linear registration and motion correction of brain images. NeuroImage, 17:825-41. Jenkinson, M., Smith, S. (2001) A global optimisation method for robust affine registration of brain images. Med Image Anal, 5:143-56. Kelly, S.C., He, B., Perez, S.E., Ginsberg, S.D., Mufson, E.J., Counts, S.E. (2017) Locus coeruleus cellular and molecular pathology during the progression of Alzheimer's disease. Acta Neuropathol Commun, 5:8. Keren, N.I., Lozar, C.T., Harris, K.C., Morgan, P.S., Eckert, M.A. (2009) In vivo mapping of the human locus coeruleus. NeuroImage, 47:1261-7. Keren, N.I., Taheri, S., Vazey, E.M., Morgan, P.S., Granholm, A.C., Aston-Jones, G.S., Eckert, M.A. (2015) Histologic validation of locus coeruleus MRI contrast in post-mortem tissue. Neuroimage, 113:235-245. Krebs, R.M., Park, H.R., Bombeke, K., Boehler, C.N. (2018) Modulation of locus coeruleus activity by novel oddball stimuli. Brain Imaging Behav, 12:577-584. Langley, J., Huddleston, D.E., Crosson, B., Song, D.D., Factor, S.A., Hu, X. (2020a) Multimodal assessment of nigrosomal degeneration in Parkinson's disease. Parkinsonism & related disorders, 80:102-107. Langley, J., Huddleston, D.E., Sedlacik, J., Boelmans, K., Hu, X.P. (2017) Parkinson's disease-related increase of T2*-weighted hypointensity in substantia nigra pars compacta. Mov Disord, 32:441-449. Langley, J., Hussain, S., Flores, J.J., Bennett, I.J., Hu, X. (2020b) Characterization of age-related microstructural changes in locus coeruleus and substantia nigra pars compacta. Neurobiol Aging, 87:89-97. Liu, K.Y., Acosta-Cabronero, J., Cardenas-Blanco, A., Loane, C., Berry, A.J., Betts, M.J., Kievit, R.A., Henson, R.N., Duzel, E., Cam, C.A.N., Howard, R., Hammerer, D. (2018) In vivo visualization of age-related differences in the locus coeruleus. Neurobiol Aging, 74:101-111. Manaye, K.F., McIntire, D.D., Mann, D.M., German, D.C. (1995) Locus coeruleus cell loss in the aging human brain: a non-random process. J Comp Neurol, 358:79-87. Messinis, L., Nasios, G., Mougias, A., Politis, A., Zampakis, P., Tsiamaki, E., Malefaki, S., Gourzis, P., Papathanasopoulos, P. (2016) Age and education adjusted normative data and discriminative validity for Rey's Auditory Verbal Learning Test in the elderly Greek population. J Clin Exp Neuropsychol, 38:23-39. Mori, S., van Zijl, P.C. (2002) Fiber tracking: principles and strategies - a technical review. NMR Biomed, 15:468-80. Murphy, P.R., O'Connell, R.G., O'Sullivan, M., Robertson, I.H., Balsters, J.H. (2014) Pupil diameter covaries with BOLD activity in human locus coeruleus. Hum Brain Mapp, 35:4140-54. Naidich, T.P., Duvernoy, H.M., Delman, B.N., Sorensen, A.G., Kollias, S.S., Haacke, E.M. (2009) Duvernoy's Atlas of the Human Brain Stem and Cerebellum. Vienna, Austria. Springer-Verlag Wien. 876 p. Pasternak, O., Koerte, I.K., Bouix, S., Fredman, E., Sasaki, T., Mayinger, M., Helmer, K.G., Johnson, A.M., Holmes, J.D., Forwell, L.A., Skopelja, E.N., Shenton, M.E., Echlin, P.S. (2014) Hockey Concussion Education Project, Part 2. Microstructural white matter alterations in acutely concussed ice hockey players: a longitudinal free-water MRI study. J Neurosurg, 120:873-81. Pasternak, O., Sochen, N., Gur, Y., Intrator, N., Assaf, Y. (2009) Free water elimination and mapping from diffusion MRI. Magn Reson Med, 62:717-30. Pasternak, O., Westin, C.F., Bouix, S., Seidman, L.J., Goldstein, J.M., Woo, T.U., Petryshen, T.L., Mesholam-Gately, R.I., McCarley, R.W., Kikinis, R., Shenton, M.E., Kubicki, M. (2012) Excessive extracellular volume reveals a neurodegenerative pattern in schizophrenia onset. J Neurosci, 32:17365-72. angley, et al.
LC Tractography 15 Pendlebury, S.T., Welch, S.J., Cuthbertson, F.C., Mariz, J., Mehta, Z., Rothwell, P.M. (2013) Telephone assessment of cognition after transient ischemic attack and stroke: modified telephone interview of cognitive status and telephone Montreal Cognitive Assessment versus face-to-face Montreal Cognitive Assessment and neuropsychological battery. Stroke, 44:227-9. Robertson, I.H. (2013) A noradrenergic theory of cognitive reserve: implications for Alzheimer's disease. Neurobiol Aging, 34:298-308. Samuels, E.R., Szabadi, E. (2008) Functional neuroanatomy of the noradrenergic locus coeruleus: its roles in the regulation of arousal and autonomic function part I: principles of functional organisation. Curr Neuropharmacol, 6:235-53. Sasaki, M., Shibata, E., Tohyama, K., Takahashi, J., Otsuka, K., Tsuchiya, K., Takahashi, S., Ehara, S., Terayama, Y., Sakai, A. (2006) Neuromelanin magnetic resonance imaging of locus ceruleus and substantia nigra in Parkinson's disease. Neuroreport, 17:1215-8. Schwarz, S.T., Rittman, T., Gontu, V., Morgan, P.S., Bajaj, N., Auer, D.P. (2011) T1-Weighted MRI shows stage-dependent substantia nigra signal loss in Parkinson's disease. Mov Disord, 26:1633–38. Smith, S.M. (2002) Fast robust automated brain extraction. Hum Brain Mapp, 17:143-55. Smith, S.M., Jenkinson, M., Woolrich, M.W., Beckmann, C.F., Behrens, T.E., Johansen-Berg, H., Bannister, P.R., De Luca, M., Drobnjak, I., Flitney, D.E., Niazy, R.K., Saunders, J., Vickers, J., Zhang, Y., De Stefano, N., Brady, J.M., Matthews, P.M. (2004) Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage, 23 Suppl 1:S208-19. Tanaka, C., Ishikawa, M., Shimada, S. (1982) Histochemical mapping of catecholaminergic neurons and their ascending fiber pathways in the rhesus monkey brain. Brain Res Bull, 9:255-70. Troadec, J.D., Marien, M., Darios, F., Hartmann, A., Ruberg, M., Colpaert, F., Michel, P.P. (2001) Noradrenaline provides long-term protection to dopaminergic neurons by reducing oxidative stress. J Neurochem, 79:200-10. Westlund, K.N., Bowker, R.M., Ziegler, M.G., Coulter, J.D. (1984) Origins and terminations of descending noradrenergic projections to the spinal cord of monkey. Brain Res, 292:1-16. Wheeler-Kingshott, C.A., Cercignani, M. (2009) About "axial" and "radial" diffusivities. Magn Reson Med, 61:1255-60. Wilson, R.S., Nag, S., Boyle, P.A., Hizel, L.P., Yu, L., Buchman, A.S., Schneider, J.A., Bennett, D.A. (2013) Neural reserve, neuronal density in the locus ceruleus, and cognitive decline. Neurology, 80:1202-8. Woolrich, M.W., Jbabdi, S., Patenaude, B., Chappell, M., Makni, S., Behrens, T., Beckmann, C., Jenkinson, M., Smith, S.M. (2009) Bayesian analysis of neuroimaging data in FSL. NeuroImage, 45:S173-86. Zhang, Y.Y., Brady, M., Smith, S. (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Transactions on Medical Imaging, 20:45-57. Zhao, Q., Guo, Q., Liang, X., Chen, M., Zhou, Y., Ding, D., Hong, Z. (2015) Auditory Verbal Learning Test is Superior to Rey-Osterrieth Complex Figure Memory for Predicting Mild Cognitive Impairment to Alzheimer's Disease. Curr Alzheimer Res, 12:520-6. angley, et al.
LC Tractography 16
Supplementary Material
S1. Introduction
A nigrostriatal tract (NST) atlas was created using procedure similar to that described in section 2.7 and used as a control tract to examine the specificity of central tegmental tract (CTT) diffusivity on memory performance and locus coeruleus (LC) subregion diffusivity.
S2. Statistical Analysis
All statistical analyses were performed using IBM SPSS Statistics software version 24 (IBM Corporation, Somers, NY, USA) and results are reported as mean ± standard deviation. A p value of 0.05 was considered significant for all statistical tests performed in this work. Normality of diffusion data was assessed using the Shapiro-Wilk test for each group and all data was found to be normal. Unpaired t-tests were used to assess age group differences in each DTI index (FA, MD, AD, RD) from the NST. We performed separate Spearman rank correlations between mean NST DTI indices and memory performance (RAVLT delayed recall scores) separately for younger and older adults. Separate Spearman rank correlations were also performed between each NST DTI index and the corresponding DTI index from each LC subregion (rostral, caudal; i.e. NST MD was correlated with rostral LC MD) for each age group. S3. Results
The older adult cohort exhibited higher FA (older: 0.28 ± ± t =-2.6, p =0.007) and lower MD (older: 3.03 × -4 mm /s ± × -5 mm /s; young: 3.45 × -4 mm /s ± × -5 mm /s; t =4.2, p <0.001), lower RD (older: 2.65 × -4 mm /s ± × -5 mm /s; young: 3.01 × -4 mm /s ± × -5 mm /s; t =3.9, p <0.001), and lower AD (older: 3.82 × -4 mm /s ± × -5 mm /s; young: 4.44 × -4 mm /s ± × -5 mm /s; t =4.6, p <0.001) relative to young adults. To probe the potential contribution of motion to the significant CTT-memory correlations observed within older adults, we tested for similar correlations between NST diffusion metrics and RAVLT delayed recall . No significant correlations were observed between diffusion metrics in NST and RAVLT delayed recall in the older group (mean tract MD: r =0.209, p =0.276; mean tract RD: r =0.165, p =0.392; mean tract AD: r =0.266, p =0.162) or in the younger group (mean tract MD: r =0.279, p =0.094; mean tract RD: r =0.280, p =0.094; mean tract AD: r =0.269, p =0.107). These results suggest that correlations between CTT diffusivity and memory in older adults are not due to motion. Group comparisons and correlations between NST and memory in the older adult group are summarized in Figure S1. The specificity of the relationship between CTT diffusivity and LC subregion diffusivity was assessed by correlating NST and LC subregion diffusivity. No significant correlations were observed between rostral LC diffusivity indices and NST diffusivity (MD: r =0.209, p =0.276; RD: r =0.165, p =0.392; AD: r =0.266, p =0.162) in the older group or in the younger group (MD: r =0.046, p =0.787; RD: r =-0.001, p =0.957; AD: r =0.139, p =0.367). For the caudal LC ROI, no significant correlation was seen between caudal LC diffusivity indices and NST diffusivity (MD: r =0.195, p =0.247; RD: r =0.135, p =0.424; AD: r =0.247, p =0.140) in the younger group or in the older group (MD: r =0.111, p =0.938; RD: r =-0.089, p =0.666; AD: r =0.183, p =0.360). These results suggest that the correlation between CTT diffusivity and LC subregion diffusivity is specific to CTT. angley, et al. LC Tractography 17
Figure S1. Group comparisons between younger and older adult groups in nigrostriatal (NST) MD (A), RD (B), and AD (C). Statistically significant decreases in mean tract MD, RD and AD were observed in the older group as compared to the younger group. No correlations between delayed recall and NST MD (D), RD (E), and AD (F) were observed in the older group.