Abstract
The brainstem is among the first regions to accumulate Alzheimer’s disease (AD)-related hyperphosphorylated tau pathology during aging. We aimed to examine associations between brainstem volume and neocortical amyloid-β or tau pathology in 271 middle-aged clinically normal individuals of the Framingham Heart Study who underwent MRI and PET imaging. Lower volume of the medulla, pons, or midbrain was associated with greater neocortical amyloid burden. No associations were detected between brainstem volumes and tau deposition. Our results support the hypothesis that lower brainstem volumes are associated with initial AD-related processes and may signal preclinical AD pathology.
INTRODUCTION
The pathophysiological process of Alzheimer’s disease (AD) starts two to three decades prior to the clinical symptoms. Several autopsy studies suggested that neuromodulatory nuclei in the brainstem may be among the first regions in the brain to accumulate tau pathology, starting as early as young adulthood [1–3]. Accrual of tau pathology has been identified in midbrain nuclei, including the substantia nigra and raphe nuclei, and in the pontine nuclei, including the locus coeruleus and pedunculopontine nuclei. Some of these changes may occur early in life; for example, the locus coeruleus has been reported to harbor hyperphosphorylated tau aggregates in approximately 50% of 30-year-old individuals [4, 5]. Similarly, raphe nucleus neurons accumulate tau deposition in precortical stages [6]. Previous studies reported lower grey matter volume in the rostral midbrain and pons in AD patients relative to cognitively normal individuals [7]. Cognitively healthy older individuals who later progressed to AD dementia also exhibited lower grey matter volume in the midbrain and in a specific pontine cluster colocalized to the locus coeruleus as compared to those who did not progress [8, 9]. These findings suggest that brainstem volumetric changes may be occurring early, starting in clinically normal individuals, and reflect underlying progression of AD pathology. As tau pathology is one of the major contributors to neurodegeneration, we sought to examine the association between brainstem volume and neocortical amyloid-β (Aβ) or tau pathology in the third generation (Gen 3) cohort of the Framingham Heart Study (FHS), comprising middle-aged clinically normal individuals.
METHODS
Participants
Individuals from the FHS Gen 3 who underwent 3T-MRI, 18F Flortaucipir (FTP), and 11C Pittsburgh Compound B (PiB)-PET were included for the present study (n = 271, mean age = 54.4 years (SD = 8, range = 32–73)). Participants in the FHS Gen 3 cohort are grandchildren of the original FHS cohort [10, 11]. The time between the FTP-PET visit, PIB-PET visit, and MRI was less than one year (mean = 0.5, SD = 1.0). All participants underwent a clinical evaluation to exclude medical or neurological disorders that could impact their cognitive abilities. The majority of participants have had APOE genotyping (N = 261, 96%). Participants underwent cognitive testing of which a global cognitive composite was derived (termed the PC1). This composite was obtained by extracting the first principal component from a principal component analysis forcing a single score solution including the Trails Making Test part B, the Hooper Visual Organization Test, the Logical Memory test, the Visual Reproductions test, the Paired Associate Learning and the Similarities. Measures that had a skewed distribution were log-transformed, and directionality was reversed such that higher scores reflect better performance. Study protocols were approved by Boston University School of Medicine and the Partners Human Research Committee of Massachusetts General Hospital and all participants provided written informed consent.
Structural magnetic resonance imaging
Anatomical scans of the FHS participants were collected at Boston University using a Philips 3T Achieva (Philips, Best, Netherland; repetition time =6800 ms, echo time = 3.1 ms, flip angle = 9 degrees, and a voxel size = 0.98 x 0.98 x 1.2 mm) All T1-weighted images were processed using FreeSurfer (FS) version 6 using the software package’s default, automated reconstruction protocol as described previously [12]. Briefly, each T1-weighted image was subjected to a volume-based and surface-based automated segmentation process involving intensity normalization, skull stripping, segregating left and right hemispheres, and in addition for the surface-based stream, brainstem and cerebellum were removed, topology defects were corrected, borders between grey/white matter and grey/cerebrospinal fluid were defined. Together, these resulted in a parcellation of cortical and subcortical areas. Using FS’s visualization toolbox, we visually inspected and, if necessary, edited each image. These edits included control point editing for over- or under-estimation of gray/white matter boundaries and manual editing the skull stripping procedures if needed. The brainstem volumes were segmented using a Bayesian segmentation algorithm that relies on a probabilistic atlas of the brainstem and neighboring structures, providing volumetric measures of the pons, medulla, and midbrain (the superior cerebellar pedunculus was not included in our analyses) [13]. The segmentation of the brainstem structures was inspected visually, and we detected no issues. As a control region not specific to AD (outside the brainstem and not a site of atrophy or amyloid/tau accumulation in initial AD), we included the striatum volume, calculated as the sum of the caudate and putamen volumes, and fourth ventricle volume, given that this is located in the vicinity of the brainstem. In addition, we also included hippocampal volume representing a control region vulnerable to early AD pathology. Volumes were expressed in percentage relative to the total cranial volume.
Positron emission tomography
All PiB-PET and FTP-PET data were acquired at Massachusetts General Hospital as previously reported [14]: PiB-PET was acquired with a 8.5 to 15 mCi (315–555 MBq) bolus injection followed immediately by a 60-min dynamic acquisition in 69 frames (12x15 s, 57x60 s) and FTP was acquired from 80–100 min after a 9.0 to 11.0 mCi (330–405 MBq) bolus injection in 4 x 5 min frames on a Siemens/CTI ECAT EXACT HR+scanner (n = 218) or a GE Discovery DMI 5-Ring TOF PET/CT scanner [15] (n = 53). To harmonize data across these cameras, GE Discovery images were smoothed with a 6mm Gaussian smoothing filter. In addition, PET camera (HR+ or GE) was included as covariate in all statistical analyses.
PET data was reconstructed and attenuation-corrected, and each frame was evaluated to verify adequate count statistics and absence of head motion. PiB PET data were expressed as the distribution volume ratio (DVR) with cerebellar grey as reference tissue by using the Logan graphical method applied to data from 40 to 60 min after injection [16]. The FTP PET data were expressed in standardized uptake values (SUVr) with cerebellar cortex as reference tissue. To evaluate the anatomy of cortical PiB or FTP binding, each individual PET data set was rigidly coregistered to the subject’s MPRAGE data. The FreeSurfer regions-of-interest (ROIs) were transformed into the PET native space. PET data were partial volume corrected (PVC) using the Geometrical Transfer Matrix method as implemented in FreeSurfer [17]. Neocortical PiB retention was also assessed in a large cortical ROI aggregate that included frontal, lateral temporal, and retrosplenial cortices (FLR) [18]. For the FTP data we a priori chose several regions known to harbor neurofibrillary tangles in the early Braak stages: entorhinal cortex, inferior temporal cortex, fusiform cortex and amygdala and averaged values bilaterally [19].
Statistical analyses
All analyses were done in SAS version 9.4. Group characteristics are represented in mean and standard deviation for continuous variables, or proportions for dichotomous data. Neocortical PiB was log-transformed because of its skewed distribution. Associations between adjusted brainstem volumes (or volumes of the control regions in separate models) and neocortical amyloid (PiB) or regional tau (FTP) deposition were examined using linear regression analyses, including age, sex, and PET camera as covariates. We also examined whether neocortical amyloid burden would modify the association between brainstem volumes and tau burden by including the interaction between PiB and brainstem volume.
Sensitivity analyses were done on the PET data without PVC. Only associations that were significant in both the PVC and non-PVC analyses were considered robust and interpreted. All reported p-values are two-sided.
RESULTS
A description of the demographics of our sample is provided in Table 1. Lower midbrain volume was significantly associated with greater amygdala and fusiform tau (non-PVC), but these associations were at-trend or nonsignificant for the PVC tau data. We observed no other associations between any of the brainstem volumes and the regional tau measurements (Table 2). But lower volume of the medulla, pons, and midbrain were each associated with greater neocortical amyloid burden (Table 2 and Fig. 1). In the non-PVC data, these relationships were also observed for the medulla and pons. These results did not change when adjusting the models for global cognition (PC1, Supplementary Table 1) or APOE E4 status (Supplementary Table 2), or when removing potential influential PiB-values > 2SD from the mean (Supplementary Table 3). In separate models, striatal volume or fourth ventricular volume were not associated with any of the PET variables, thus supporting the specificity of these findings. Hippocampal volume was not significantly associated with amyloid burden or tau pathology. No other associations between hippocampal volume and tau pathology were detected. We observed no interaction between neocortical amyloid burden and any of the brainstem volumes on cortical tau burden measures (Table 2).
Demographics of the sample
*given the skewed distribution of this variable, the distribution is summarized here using the median and interquartile range. PET, positron emission tomography; ICV, intracranial volume; PVC, partial volume correction; Q1/Q3, quartile 1/3; PiB, Pittsburgh Compound B; FTP, Flortaucipir.
Associations between brainstem volumes and Alzheimer’s disease pathology
β±SE (p-value) presented unless otherwise specified. For interaction analyses only the p-value is presented. All models adjusted for age, sex, and camera. Analyses were also adjusted for total cranial volume. Neocortical PiB was log transformed for normality. Top section of the table shows the partial volume corrected (PVC)-results, the bottom part the non-PVC results. *Number of observations for hippocampal volume was n = 257 for associations with neocortical PiB and n = 227 for tau measures. †Number of observations for fourth ventricle was n = 231 for associations with neocortical PiB and n = 219 for tau measures. PVC, partial volume correction; PiB, Pittsburgh Compound B; FTP, Flortaucipir.

Brainstem volumes are associated with neocortical PiB Unadjusted associations between the brainstem volume (as a % of total cranial volume) and neocortical PiB deposition (log-transformed, here referred to as amyloid FLR). FLR, fronto-lateral temporal–retrosplenial aggregate (neocortical region); PiB, Pittsburgh Compound B.
DISCUSSION
Several brainstem nuclei exhibit AD-related neurodegenerative processes early in adulthood [1], hence we examined associations between subregional brainstem volumes and the initial stages of AD pathology in a large cohort of middle-aged cognitively normal individuals. We found that lower brainstem volumes, in particular the medulla and the pons, were associated with greater Aβ deposition. Contrary to our expectations, based on the Braak staging [4] and previous MRI-work [8], we did not observe robust relationships between any of the brainstem volumes and tau deposition. Effect sizes for the amygdala and fusiform gyrus—regions with strong connections with the brainstem—are in a similar range to those of neocortical PiB, providing some initial support for the hypothesis that neurodegenerative diseases including AD start in the brainstem. Interestingly, hippocampal volume was not associated with Aβ or tau deposition in our sample. There is a plethora of studies reporting correlations between hippocampal volume and amyloid or tau pathology, or disease progression, in older or impaired individuals [20–24]. However, in younger or midlife individuals, previous studies reported no relationship between APOE genotyping and hippocampal volume [25–27], indicating that loss of grey matter tissue is a downstream effect of accumulating pathology [28]. In fact, recent work by Bussy et al. (2019) indicated that amyloid may already be accumulating at younger ages, prior to any measurable volumetric changes [25]. We observed that lower brainstem volumes were associated with initial amyloid accumulation in the neocortex, the earliest sign of preclinical AD. Based on these observations, we speculate that our associations between brainstem volumes and amyloid pathology reflect brain changes that occur prior to cortical neurodegeneration and thus may be valuable in the detection of preclinical AD. Other physiological markers of brainstem function, such as changes in heart rate variability, have also been associated with a greater risk of developing dementia. However, autopsy studies have identified a long-time lag between tau deposition in the locus coeruleus and neuronal loss [29], indicating that the observed volumetric changes in the locus coeruleus may reflect loss of projection fiber density or other elements besides locus coeruleus neuronal loss. Similarly, while tau accumulation in the raphe nucleus occurs prior to cortical tau deposition, neuronal loss in the raphe nucleus was typically detected only in later Braak stages [6, 30]. Thus, the changes in brainstem volume may reflect processes other than tau deposition in brainstem nuclei. Previous volumetric MRI studies have already reported a relationship between lower brainstem volume and worse cognition in prodromal AD, but not in normal older individuals [7, 9] and no prior study has related brainstem volumes to brain amyloid and tau burden. It is important to note that these findings were specific to the brainstem, and not observed for any of the control regions, including the hippocampus. We did not yet probe the potential subtle functional implications of these AD-related changes in brainstem volume, but plan to do so when our sample size is large enough to capture cognitive variability in this cognitively healthy middle-aged population.
Unfortunately, detecting tau aggregation in vivo in these tiny nuclei of the brainstem is not yet feasible, due to the off-target binding of current tau radiotracers to neuromelanin [31, 32], which is abundantly present in the substantia nigra and locus coeruleus. Novel emerging imaging methods may enable examination of the integrity of the locus coeruleus, substantia nigra and potentially also the other neuromodulatory subcortical systems, as well as their molecular alterations during aging or disease contributing to volumetric changes [33–35]. To better understand the temporal dynamics between volumetric changes and deposition of cortical beta-amyloid and tau, future studies should explore these associations through longitudinal studies that include repeated MRI and PET imaging.
Footnotes
ACKNOWLEDGMENTS
This work was supported by National Institutes of Health grants R01 AG049607 (to S.S.), R01 AG062559 (to H.I.L.J.), HHSN268201500001I and 75N92019D00031 (to S.S.). RFB is supported by a K99/R00 award from NIH-NIA (R00AG061238), an Alzheimer’s Association Research Fellowship and philanthropic support.
