Graff‐Radford, Jonathan

AbstractBackgroundDiffusion MRI (dMRI) based quantification of microstructural changes in white matter (WM) are increasingly proposed to measure cerebrovascular disease (CVD) related changes. Our objective was to compare the proposed dMRI markers using postmortem neuropathologic CVD data and then evaluate the predictors of longitudinal dMRI markers using serial scans.MethodWe identified n=51 participants (mean age: 83.8 years, 63% males) with ante‐mortem dMRI and postmortem CVD evaluation and n=718 participants (mean age: 71.1 years, 56% males) with at least two dMRI scans from the population‐based sample of Mayo Clinic Study of Aging. We computed dMRI measures proposed for measuring CVD: Free Water (FW), Fractional Anisotropy (FA) adjusted for FW (FAadj), Peak width Skeletonized Mean Diffusivity (PSMD), and FA of the genu of corpus callosum (Genu‐FA). We used the FW and PSMD (PSMD release 1.8.1) kits from the MarkVCID consortium.Using weighted linear regression models with adjustments for MRI scan time to death, we evaluated associations between the baseline dMRI and two pathology CVD scores: Strozyk (represents the presence and number of macroscopic lesions) and Kalaria (represents the summary score of vessel wall modifications) scales. We ran linear mixed effect models with all available serial dMRI measures as outcomes and vascular risk (measured by the number of cardiovascular and metabolic conditions‐CMC), baseline imaging measures (amyloid and white matter hyperintensities (WMH)) and their interaction with time as key predictors.ResultMost dMRI markers were predictors of postmortem CVD pathology (Table 1). In longitudinal dMRI models, terms associated with amyloid explained little variability in dMRI (<1%) (Table 2). WMH was a significant predictor of baseline and decline in dMRI measures with WMH explaining a considerable percentage of variability in FW (34%) and PSMD (28%). Vascular risk measures were significant predictors of three out of the four dMRI measures.ConclusionThe diffusion measures proposed as surrogate markers of CVD map reasonably well to CVD pathological scales and do not change substantially as a function of baseline amyloidosis. Our findings shed light on variability in the proposed diffusion markers using postmortem data and longitudinal imaging data. Further work is needed to evaluate their clinical utility.

AbstractBackgroundCerebral Microbleeds (CMBs) are tracked in anti‐amyloid drug safety monitoring. T2*GRE or SWI MRI images are visually assessed for the presence of CMBs. Automated methods tuned for sensitivity often have unacceptably high false positive rates requiring extensive human cleanup.MethodAssume contemporaneous T2*GRE and T1‐weighted images. CMB detection (Figures 1,2) has three steps: 1) Segmentation outputs from the T1‐weighted image are resampled into the T2*GRE image voxel raster; 2) Normalized cross‐correlation of the T2*GRE with a CMB template identifies “candidate CMBs”; 3) a neural network using voxel patches centered on each candidate CMB rejects false positives. Patches presented to the network are 4D and include five contrasts drawn from the T2*GRE, correlation, T1‐weighted, WM probability and CSF probability images. The simple network architecture consists of convolutional blocks followed by dense blocks. The system is trained using images from ADNI, ARIC, MCSA, and the Mayo Clinic ADRC that have been previously visually graded for CMBs (Table in Figure 3). Importantly, we include participants with and without CMBs in our data reflecting the fact that the majority of people (even the elderly) have no CMBs. Data is split 70:10:20 by subject for training, validation, and testing. Various network configurations were evaluated.ResultHigh candidate detection sensitivity can be achieved at the expense of presenting hundreds to thousands of false positive candidates per input image to the rejection network. Our best performing network configuration to date has AUC 0.982 (95% bootstrap CI 0.971 to 0.990 based on the 10% validation subset). Choosing an operating point with 95% sensitivity in the validation set and admits, on average, 0.6 false positive CMB candidates per input T2*GRE image. Sensitivity is hindered by a small number of obvious errors in the ground truth visual assessment. Validation loss tracks training loss with no indication of over‐fitting. Example patches are shown in Figure 4.ConclusionExcessive false positive detection requires significant human correction and effort, preventing automated CMB detection from being widely adopted. We have an automated detection system that is sensitive (95%) and provides excellent false positive rejection (<1 false positive CMB/input image) reducing needed human effort.

AbstractBackgroundAlzheimer’s disease (AD) pathologic changes (e.g. amyloid and tau aggregation), and cerebrovascular disease often coexist, and it remains unclear whether these are mechanistically linked processes (Graff‐Radford, Brain 2019; Provenzano, JAMA Neurol 2013). Most studies evaluating these associations are cross‐sectional; however, establishing temporal relationships may provide additional insight into potential mechanistic associations. One method to assess the temporal evolution of a biomarker is the Accelerated Failure Time (AFT) model, which assumes a common pattern of progression for a given biomarker that is shifted earlier or later in time for each individual (Therneau, Neuroimage 2021). In this work we investigate the relationship of disease trajectories of cerebrovascular disease (via white matter hyperintensity, WMH, volume) and AD biomarkers (via amyloid and tau PET).MethodsWe included Mayo Clinic Study of Aging (MCSA) and Mayo Alzheimer’s Disease Research Center (ADRC) participants with amyloid PET, tau PET, and MRI, age > 50 years, and a clinical diagnosis of cognitively unimpaired (CU), mild cognitive impairment (MCI) or Alzheimer’s disease dementia. The AFT model was fit with amyloid PET (PiB) global meta‐ROI SUVR, tau PET (flortaucipir) temporal meta‐ROI SUVR, and WMH volume as the tri‐variate endpoints, with APOE genotype, sex, education, total intracranial volume, and referral/population subject as covariates. A per subject/endpoint random effect represents individual adjustments for each biomarker to indicate if the subject's biomarker became abnormal earlier or later relative to their demographic peers, accounting for the covariates.ResultsThe study included 1144 participants (Table 1). There was a strong correlation between earlier onset of amyloid and earlier onset of tau (r=0.64, CI=0.60‐0.67) (Figure 1, Table 2); early/late amyloid timing explained nearly half (0.642 = 41%) of the variation in tau onset. There were only weak associations between an individual’s WMH and amyloid PET onset adjustments (r=0.10, 1.0%) and WMH and tau PET onset adjustments (r=‐0.07, 0.5%).ConclusionThis analysis supports a strong mechanistic relationship between amyloid and tau aggregation, as reported by others. However, our findings do not support a mechanistic association between WMH and amyloid or tau PET.

AbstractBackgroundDysexecutive Alzheimer's disease (dAD) is a relatively early‐onset variant of AD primarily degrading core executive functions in the absence of predominant behavioral symptoms. Clinical observations suggest substantial variability in clinical, cognitive and neuroimaging profiles within this syndrome. In this presentation, we will discuss our ongoing work disentangling this heterogeneity using unsupervised machine learning techniques.MethodWe collected clinical, neuropsychological and multimodal imaging (MRI, FDG‐PET, amyloid‐PET, tau‐PET) data on 52 dAD patients assessed in our behavioral neurology clinic at Mayo Clinic Rochester. We first performed a spectral clustering analysis based on FDG‐PET to delineate latent patterns of network degeneration (referred as “eigenbrains”) and assessed the relationships between these eigenbrains and clinical and cognitive symptomatology. We then performed a hierarchical clustering on these eigenbrains to derive data‐driven subtypes of dAD. We compared clinical and cognitive data between subtypes, and then compared the imaging profiles to 52 amyloid‐negative age‐ and sex‐matched controls.ResultSix eigenbrains explained approximatively 48% of the variance in FDG‐PET patterns and primarily reflected heteromodal to primary motor/sensory, left‐right hemispheric asymmetry, and anterior‐to‐posterior gradients of macro‐scale cortical organization (Fig 1). These eigenbrains differentially related to reported age at symptom onset, degree of clinical impairment, and performance on a wide range of cognitive domains (executive functions, episodic memory, and visuospatial). Hierarchical clustering revealed four dAD subtypes (diffuse, biparietal, left‐hemisphere, right‐hemisphere). These subtypes differed in reported age at symptom onset and cognitive profile, where the heteromodal‐diffuse subtype exhibited an overall worse clinical picture and the biparietal had a milder profile compared to other subtypes, which was not explained by disease duration. Additionally, spatial patterns of tau distribution and neurodegeneration overlapped with patterns of FDG‐PET hypometabolism in each subtype, whereas patterns of amyloid deposition were similar across subtypes (Figs 2 & 3).ConclusionAlmost half of the variance in macro‐scale patterns of hypometabolism in this dAD cohort was accounted for by six eigenbrains. These eigenbrains captured the inter‐individual variability in age at symptom onset and cognitive impairment. Four dAD subtypes derived from these eigenbrains revealed meaningful differences in clinical, cognitive, and imaging profiles. Recognizing this heterogeneity has significant clinical implications for diagnosis, prognosis, and symptom monitoring.

AbstractBackgroundData heterogeneity due to cross‐scanner variations is a major challenge in multi‐center aging and dementia studies. Brain structural measures for the same participant can vary due to scanner differences (manufacturer, scanner age/technology, signal‐to‐noise ratio, pulse sequence design, contrast, resolution) and result in biomarker variability. In this work, we applied basic image processing and deep learning methods (style transfer (ST) and generative adversarial networks (GANs)) for harmonizing T1‐weighted MRI scans across change in vendors (GE to Siemens) and resolution (1x1x1.2mm to 0.8x0.8x0.8mm) and test whether this harmonization improved compatibility of AD biomarker measures.MethodWe utilized training (n=814) and test (n=113 same participants scanned on two scanners close‐in‐time, typically on the same day, with different resolution) data from cognitively impaired and unimpaired participants scanned on GE and Siemens scanners from Mayo Clinic Study of Aging and Mayo ADRC. The harmonization was conducted to translate GE to Siemens. All scans were affine transformed to the same template space, resampled to 1.5 mm3, and bias corrected. Histogram matching was performed using cumulative density functions. A ResNet18 was trained for the ST’s content and style layers. Two types of GANs, conditional GANs (cGAN) and cycle‐consistent GANs (cycleGAN), were constructed. We then measured regional cortical thickness from each image using FreeSurfer‐7.1.1 and compared the thickness measurements from the original and translated scans using intraclass‐correlation (ICC) and Deming regression.ResultIn the test dataset, ST and cycleGAN showed the best performances when comparing the translated GE to Siemens. With ST, the ICC increased (vs. the original scans) slightly for the temporal meta‐ROI (0.956 from 0.954) and more substantially for the frontal pole ROI (0.706 from 0.594). The cycleGAN showed an increase in ICC for the frontal pole ROI (0.610 from 0.594) but not the meta‐ROI (Fig. 1).ConclusionEach method leveraged a different strength. Histogram matching captures contrast, ST extracts textural information, and GANs learn combination of contrast, edge and texture properties. We found that deep learning is a promising approach to cross‐scanner variability. Our comparisons revealed anatomical differences across methods; thus, disease specific measure should be considered for selection of the scanner harmonization method.

AbstractBackgroundSeveral hypotheses exist regarding the neurologic substrate underlying visual hallucinations (VH) in Dementia with Lewy Bodies (DLB). Studies have suggested an association with regional cerebral hypometabolism in the primary visual or visual association cortices on FDG‐PET. We compared regional cerebral glucose metabolism and posterior cingulate island sign (CIS) ratios in DLB and in Mild Cognitive Impairment with Lewy Bodies (MCI‐LB) participants with (VH+) and without (VH‐) VH.MethodParticipants with clinically probable DLB or MCI‐LB and an FDG‐PET scan were identified from the Mayo ADRC database. Responses on the Visual Hallucinations and Delusions Questionnaire and Neuropsychiatric Inventory were used to assess presence or absence of VH. Standardized uptake value ratios (SUVr) were used to calculate regional cerebral glucose hypometabolism in the primary visual cortex (PVC), lateral occipital cortex (LOC), and whole occipital cortex (WOC). The CIS ratio was defined as posterior cingulate divided by the sum of the precuneus plus cuneus FDG‐PET uptake ratio.ResultTwo cohorts were identified (27 DLB and 18 MCI‐LB, 87% male). Eighteen DLB participants were VH+ and 9 were VH‐, whereas 5 MCI‐LB were VH+ and 13 were VH‐. There were no statistically significant differences in demographic (age, sex, education) and clinical features (MoCA, CDR, presence of other core DLB features, and UPDRS) between VH+ vs VH‐ in those with DLB, and between VH+ vs VH‐ in those with MCI‐LB. Similarly, there were no statistically significant differences in DLB between VH+ vs VH‐ in PVC (1.55±0.14 vs 1.58±0.15), LOC (1.22±0.16 vs 1.28±0.12), WOC (1.29±0.15 vs 1.34±0.11), or CIS ratio (1.16±0.08 vs 1.14±0.10), nor in MCI‐LB between VH+ vs VH‐ in PVC (1.68±0.19 vs 1.58±0.14), LOC (1.40±0.21 vs 1.32±0.14), WOC (1.45±0.20 vs 1.37±0.14), or CIS ratio (1.17±0.08 vs 1.19±0.10) (P>0.05 for all comparisons).ConclusionThese findings do not support a direct correlation between regional cerebral glucose metabolism and VH in DLB or MCI‐LB. Further work in DLB and MCI‐LB with larger numbers is warranted. Other neurologic substrates and mechanisms, such as REM sleep intrusions during wakefulness, neurotransmitter dysfunction, and disruption of functional connectivity across brain regions should be explored as contributors to VH in DLB.

AbstractBackgroundSlow walking speed is associated with numerous negative health outcomes including global cognitive decline.1 More detailed gait analysis can differentiate among different dementia subtypes,10 but the predictive power of baseline gait metrics for domain‐specific cognitive decline needs further characterization.MethodWe conducted a follow up study on 534 participants of the Mayo Clinic Study of Aging who underwent evaluations every 15 months (median follow‐up = 5.4 years). We analyzed quantitative gait metrics, acquired with a pressure‐sensing walkway (GAITRite® instrument), and cognitive scores using a linear mixed effects model to identify correlations between z‐scored baseline gait metrics (cadence, swing time, stance time, velocity, double support time (DST), stride length, stride length standard deviation (SD), stride time SD, swing time SD, stance time SD, DST SD) and cognitive domain z‐scores (memory, attention, language, visuospatial, global) at baseline and over time. Models controlled for age, sex, and probable REM sleep behavior disorder at baseline.ResultThe table includes both baseline associations and difference in slopes (Δ), relative to the average person, in cognitive domains. These slopes ranged from about ‐0.05 for visuospatial to ‐0.12 for attention.ConclusionIncreased gait velocity and stride length are associated with less severe domain‐specific and global cognitive decline while increased swing time SD and DST are associated with more severe rates of decline in some domains. These findings suggest that baseline gait characteristics could be a helpful predictor of future cognitive changes.References:1. Savica R, Wennberg AMV, Hagen C, et al. Comparison of Gait Parameters for Predicting Cognitive Decline: The Mayo Clinic Study of Aging. Journal of Alzheimer's disease : JAD. 2017;55(2):559. doi:10.3233/JAD‐1606972. Mc Ardle R, del Din S, Galna B, Thomas A, Rochester L. Differentiating dementia disease subtypes with gait analysis: feasibility of wearable sensors? Gait & Posture. 2020;76:372–376. doi:10.1016/j.gaitpost.2019.12.028

AbstractBackgroundCognitive resilience, the discrepancy observed in individuals with high Alzheimer’s disease (AD) neuropathology and near‐normal cognition, is a growing area of study, wherein multiple factors are responsible. Currently, the role of the microglial activation in cognitive resilience is unknown. Herein described is a case report characterizing the amyloid responsive microglia (ARM) in an 88‐year‐old male decedent homozygous for APOE E4, clinically cognitively intact, and found to have AD neuropathologic change (ADNC) upon autopsy.MethodThe participant was enrolled in the Mayo Clinic Study of Aging (MCSA) and continuously followed for 14 years from 73‐ to 87‐years‐old, wherein serial MRI and clinical assessments were obtained. An autopsy was performed upon death, and routine neuropathologic sections were taken. Formalin‐fixed, paraffin‐embedded sections were stained with H&E, and single immunohistochemistry was performed using antibodies against tau (AT8), beta‐amyloid (6F3D), alpha‐synuclein (LB509), TDP‐43 (pS409/410), and CD163 (10D6). Digital image analysis was performed on dorsolateral prefrontal cortex sections scanned at 40x (Leica GT450). ARM (CD163% area) and beta‐amyloid (% area) measurements were obtained using QuPath threshold classifiers. Descriptive statistics were performed (Graphpad).ResultThe decedent was an 88‐year‐old male, homozygous for APOE E4, with a history of three cardiovascular and metabolic conditions, and lifelong smoking. His MMSE (average 28/30) and global cognition (average ‐0.51) were relatively stable (Figure 1), with no clinical concerns for cognitive impairment. Cortical thickness was normal (temporal meta‐ROI>3 mm) on four serial MRIs. The participant reported high cognitive, physical activities, and education‐occupation scores (upper quartile for the 70+ population). Neuropathologic evaluation was notable for severe cerebral amyloid angiopathy, high ADNC (A3B3C3), and remote microinfarcts. CD163‐positive ARM were observed in the neocortices and perivascular spaces, and the ARM:beta‐amyloid ratio was 0.26 (Figure 1). This ARM:beta‐amyloid ratio trended upwards (>90th percentile) compared to APOE E4/E4‐matched‐ and age‐matched controls provided from the University of Pennsylvania.ConclusionThis case report demonstrates a cognitively resilient individual with an elevated ARM:amyloid ratio. Conversely, lower ARM ratios were observed in a cognitively impaired cohort, raising the possibility ARM may associate with cognition. ARM may represent a protective and phagocytic population, and further studies correlating microglia, cognition, and putative mechanism are required.

AbstractBackgroundCerebral Microbleeds (CMBs) are tracked in anti‐amyloid drug safety monitoring. T2*GRE or SWI MRI images are visually assessed for the presence of CMBs. Automated methods tuned for sensitivity often have unacceptably high false positive rates requiring extensive human cleanup.MethodAssume contemporaneous T2*GRE and T1‐weighted images. CMB detection (Figures 1,2) has three steps: 1) Segmentation outputs from the T1‐weighted image are resampled into the T2*GRE image voxel raster; 2) Normalized cross‐correlation of the T2*GRE with a CMB template identifies “candidate CMBs”; 3) a neural network using voxel patches centered on each candidate CMB rejects false positives. Patches presented to the network are 4D and include five contrasts drawn from the T2*GRE, correlation, T1‐weighted, WM probability and CSF probability images. The simple network architecture consists of convolutional blocks followed by dense blocks. The system is trained using images from ADNI, ARIC, MCSA, and the Mayo Clinic ADRC that have been previously visually graded for CMBs (Table in Figure 3). Importantly, we include participants with and without CMBs in our data reflecting the fact that the majority of people (even the elderly) have no CMBs. Data is split 70:10:20 by subject for training, validation, and testing. Various network configurations were evaluated.ResultHigh candidate detection sensitivity can be achieved at the expense of presenting hundreds to thousands of false positive candidates per input image to the rejection network. Our best performing network configuration to date has AUC 0.982 (95% bootstrap CI 0.971 to 0.990 based on the 10% validation subset). Choosing an operating point with 95% sensitivity in the validation set and admits, on average, 0.6 false positive CMB candidates per input T2*GRE image. Sensitivity is hindered by a small number of obvious errors in the ground truth visual assessment. Validation loss tracks training loss with no indication of over‐fitting. Example patches are shown in Figure 4.ConclusionExcessive false positive detection requires significant human correction and effort, preventing automated CMB detection from being widely adopted. We have an automated detection system that is sensitive (95%) and provides excellent false positive rejection (<1 false positive CMB/input image) reducing needed human effort.

AbstractBackgroundIn 2018 we presented the Computational DESH (CDESH) metric (DOI:10.1016/j.nicl.2018.11.015), a machine learning measure of abnormal CSF space patterns commonly observed on MRI in Disproportionately Enlarged Subarachnoid‐space Hydrocephalus (DESH‐type) NPH. We have shown that CDESH positive individuals often have discordant CSF and PET amyloid measures (DOI:10.1016/j.jalz.2019.06.4921). Positive CDESH scores are predictive of cognitive decline (DOI:10.1016/j.neuroimage.2021.117899, DOI: 10.1212/WNL.0000000000008616). The score is calculated from regional CSF volumes in a sulcal atlas and within the ventricles using a support vector machine (SVM). The Mayo Clinic Study of Aging (MCSA) changed scanner platforms (GE to Siemens) and T1‐weighted image resolution (1.2mm x 1.05mm x 1.05mm to 0.8mm isotropic) in 2017. We refer to these as “isotropic” or “anisotropic” scans. Our objective was to assess the dependence of a machine learning method (specifically CDESH) to image input variation and, if necessary, determine what needed corrections.MethodAll data considered is from the MCSA or Mayo Clinic ADRC. The SVM was trained on anisotropic images with laboriously labelled ground truth. We resampled the isotropic data from a crossover study (N=112, typically scanned same‐day with both platforms/resolutions) to match the anisotropic resolution. Two approaches to correction were evaluated (Figure 1). In Figure 2 we compare of the score with and without correction cross‐sectionally in 998 MCSA participants. In Figure 3 we present longitudinal series from 23 participants with multiple scans using the isotropic resolution protocol and at least one prior anisotropic image.ResultsComparisons of data from the crossover study are presented in Figure 1. Matching voxel resolution to the data used to train the SVM is critical and a subsequent mild linear score adjustment may be warranted. Cross‐sectional comparisons (Figure 2) confirm improvement matching voxel resolution with a possible mild linear score adjustment. After matching voxel resolution CDESH score trajectories (Figure 3) are markedly more stable across this scanner/protocol change.ConclusionThe SVM‐based CDESH score is sensitive to voxel resolution. Corrections are needed in image processing and possibly in the scalar output score to achieve stability across protocol and scanner changes. The resulting correction supports coherent long timeseries data and longitudinal extension of previous work.