Citations

AbstractBackgroundQuantitative susceptibility mapping (QSM) can detect iron distribution in the brain by estimating the local tissue magnetic susceptibility of every voxel. Dysregulation of neuronal iron has been correlated with typical Alzheimer’s disease (tAD), but little is known about iron distribution in atypical clinical presentations of AD such as logopenic progressive aphasia (LPA) and posterior cortical atrophy (PCA).MethodEleven PCA patients and eight LPA patients were recruited by the Neurodegenerative Research Group at Mayo Clinic, Rochester, MN, and underwent MRI that included a five‐echo gradient echo sequence for calculation of QSM. Mean QSM signal was extracted from gray and white matter for regions‐of‐interest across the brain using the Mayo Clinic Adult Lifespan Template. Hierarchical models were fit per‐region and per‐hemisphere to compare PCA and LPA to each other, to 20 healthy controls, and to 20 patients with tAD recruited by the Alzheimer’s Disease Research Center.ResultStrong evidence (posterior probability>0.99) was observed for greater susceptibility in the left middle occipital gyrus and amygdala in both LPA and PCA, and in the right precuneus and inferior temporal gyrus in PCA compared to controls. Moderate evidence for greater susceptibility (posterior probability>0.90) was also observed in the left insula, inferior occipital gyrus, entorhinal cortex, and right inferior parietal gyrus and amygdala in both LPA and PCA, along with right angular gyrus and middle occipital gyrus in PCA and right insula and caudate in LPA, when compared to controls. Comparing between phenotypes, LPA showed moderate evidence for greater susceptibility in the caudate, hippocampus, substantia nigra, left supramarginal gyrus and posterior cingulate compared to PCA, while PCA showed greater susceptibility in the right precuneus compared to LPA. Both LPA and PCA showed moderate‐strong evidence for greater susceptibility than tAD, particularly in medial and lateral parietal regions.ConclusionUnique iron profiles were identified in LPA and PCA within cortical and subcortical structures. These changes are in line with disease‐related neurodegenerative changes in the clinical phenotypes, suggesting that QSM could be an informative marker to study iron deposition in these patients.

AbstractBackgroundFluorodeoxyglucose positron emission tomography (FDG‐PET) is an established modality for differential diagnosis of dementia. Deriving phenotypic signatures from FDG‐PET via machine learning is challenging due to the high dimensional nature of brain images relative to the generally small number of labeled examples available for training, the class imbalance among those labels, and the cooccurrence of multiple pathologies. In this study, we developed a multi‐class, multi‐label framework to address these challenges.MethodA database of clinically acquired PET/CT images from 3,000 unique patients was used to develop a latent space model using matrix decomposition. This model was then applied to images from a separate cohort of Mayo Clinic Alzheimer’s Disease Research Center participants (n=1,745) labeled as cognitively unimpaired (CU) (n=1,436) or with the following potentially co‐occurring phenotypes: Alzheimer’s disease (AD) (n=165), Lewy body dementia (DLB) (n=92), behavioral variant frontotemporal dementia (bvFTD) (n=43), semantic (svPPA) (n=10) and logopenic (lvPPA) (n=13) variant PPA, and posterior cortical atrophy (PCA) (n=17). A k‐nearest neighbors classifier that is robust to these imbalanced and overlapping labels was then trained on these examples. The resulting classifier was evaluated by area under receiver‐operator characteristic curve (ROC‐AUC) via leave one out cross validation, using clinical diagnosis as the gold standard.ResultROC curves and AUC scores for each phenotype are illustrated in Fig. 1a. Because the classifier is based on a k‐nearest neighbors connectivity matrix, it has a convenient graphical representation, where images are nodes and edges are drawn between an image and its set of nearest neighbors in latent space. A self‐organizing force directed graph constructed in this way is illustrated in Fig 1b, highlighting the strong separation of CU and degenerative images, as well as the segregation of each phenotype within the neurodegenerative region of the graph.ConclusionIn this study, we developed a machine learning framework for classification of neurodegenerative disease based on k‐nearest neighbor analysis in a low dimensional latent space projection of FDG‐PET images. By leveraging low‐dimensional representations and k‐nearest neighbors analysis, this framework is robust in multi‐class, multi‐label tasks with strong class imbalance and provides a highly interpretable graphical representation.

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.

AbstractBackgroundIn vivo tau‐positron emission tomography (PET) is an attractive biomarker for Alzheimer’s disease (AD) diagnosis and treatment. However, tau‐PET is less widely available than other modalities. In this study, we tested cross‐modality synthesis of tau‐PET brain images from fluorodeoxyglucose F‐18 (FDG)‐PET using a deep convolutional neural network (CNN).MethodParticipants (n=1,192) who had brain FDG‐PET with 18F‐FDG and tau‐PET with Flortaucipir (F‐18‐AV‐1451) were included for training and testing. This cohort spanned normal aging (ages 26‐98), pre‐clinical, and clinical AD and related disorders including the FTD and DLB spectrum. External validation was done using ADNI (n=288). The PET scans were co‐registered to the corresponding MRI and subsequently warped to Mayo Clinic Adult Lifespan Template (MCALT) space. Tau‐PET images were SUVR‐normalized to the cerebellar crus, and FDG to the pons. A 3D dense‐U‐net model was utilized as an architecture. Cross‐validation experiments were conducted using 5‐fold validations (60% training set, 20% validation set, and 20% test set) with mean squared error as the loss function.ResultOur dense‐U‐net model successfully synthesized tau‐PET from metabolic images with good correlation and low prediction error for regional SUVRs (Figure 1A‐C). The model showed a robust prediction ability, performing accurately in an independent, external ADNI cohort (Figure 1D‐F). The model‐imputed tau‐PET significantly improved performance in classifying tau positivity (mean AUROC(±SD)=0.78±0.04 and 0.85±0.03 for FDG‐PET and synthesized tau‐PET, respectively) and diagnostic groups (cognitively unimpaired with abnormal amyloid‐PET vs. cognitively impaired with abnormal amyloid‐PET) compared to the original input FDG data (mean AUROC(±SD)=0.89±0.04, 0.85±0.05 0.91±0.04 for actual tau‐PET, FDG‐PET and synthesized tau‐PET, respectively) (Figure 2), suggesting enhanced clinical utility for metabolic images. The ADNI cohort also showed similar results (for tau positivity: AUROC=0.66 and 0.78 for FDG‐PET and synthesized tau‐PET, respectively; for CU A+ vs. CI A+: AUROC=0.86, 0.62, and 0.73 for actual tau‐PET, FDG‐PET, and synthesized tau‐PET; Figure 3).ConclusionWe showed that using a CNN model to predict tau‐PET from FDG‐PET is feasible. The synthesized tau‐PET can augment the value of FDG‐PET, facilitating the multi‐modal diagnosis of AD.

AbstractBackgroundDisruption of the default mode network (DMN) is thought to be a common feature of aging and across Alzheimer’s disease (AD) clinical phenotypes, while syndromic variability could be reflected in selective degeneration of macro‐scale networks optimally supporting specific mental abilities (Fig 1). As part of a collaborative effort between Mayo Clinic and the Human Connectome Project‐Aging (HCP‐A), we aimed to delineate a biomarker of global network failure across aging and the dysexecutive (dAD) and amnestic phenotypes of AD, and biomarkers of focal degeneration specific to each AD phenotype.MethodWe used a regression framework to assess the relationships between a marker of DMN dysfunction, the network failure quotient (NFQ), and factors relevant to AD such as age, cortical thickness (global cortical thickness, inferior parietal lobe), hippocampal volume, and cognition (episodic memory, cognitive flexibility) in the HCP‐A cohort, which consists of 725 non‐demented older adults. We then fit centile curves to the HCP‐A cohort to assess whether the NFQ and focal degeneration markers could distinguish dAD (n=10) and amnestic AD (n=8) patients from the HCP‐A cohort at the group‐ and individual‐level. Participants were scanned using the HCP‐A protocol, allowing for the acquisition of high‐resolution multi‐echo structural images and longer resting‐state connectivity acquisition time compared to traditional protocols.ResultThe NFQ related to age, cortical thickness, hippocampal volume, and cognition in the HCP‐A cohort. It also strongly discriminated dAD and amnestic AD from the HCP‐A cohort but did not discriminate between AD phenotypes. In contrast, markers of focal degeneration were phenotype‐specific; parieto‐frontal areas yielded largest effect sizes for dAD patients, and hippocampal and temporal areas best discriminated patients with amnestic AD (Fig 2).ConclusionThe NFQ may reflect shared system‐level pathophysiological mechanisms across aging and the dysexecutive and amnestic phenotypes of AD whereas focal degeneration markers highlight distinct pathognomonic processes across AD phenotypes. These findings suggest that inter‐individual variability in cognitive impairment in AD would be closely tied to modular, phenotype‐specific network degeneration combined with a common global functional disruption of the DMN. Such information has crucial implications for clinical biomarkers and to provide a deeper understanding of heterogeneity in AD.

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.

AbstractBackgroundWe previously demonstrated that Alzheimer’s disease (AD) patients ≥70 y/o tend to have smaller medial temporal volumes in the presence of relatively less medial temporal 18F‐flortaucipir‐PET uptake compared to young AD patients. TAR DNA‐binding protein 43 (TDP‐43) is frequently found in the brains of AD patients and is associated with old age and medial temporal atrophy. In this study, we aimed to determine whether TDP‐43 could be contributing to the mismatch between medial temporal volumes and flortaucipir‐PET uptake in older AD patients.MethodWe identified 77 participants (median age at death 79 years [range: 51‐99]) who had been enrolled into the Mayo Clinic Alzheimer’s Disease Research Center or Neurodegenerative Research Group (NRG), had undergone flortaucipir‐PET and volumetric MRI scans close to death, and autopsy confirmed AD neuropathological changes, as well as TDP‐43 assessments. Grey matter volumes and flortaucipir standard uptake value ratios (SUVRs) were calculated for hippocampus, amygdala, entorhinal, inferior temporal, and midfrontal cortices. We used bivariate‐response linear regression modelling, adjusting for Braak neurofibrillary tangle stage, to estimate the effect of TDP‐43 and age on gray matter volume and flortaucipir uptake, and tested whether TDP‐43 effect sizes varied according to age.ResultAfter accounting for Braak stage, older individuals tended to have significantly lower flortaucipir uptake across regions and lower hippocampal and entorhinal cortex (ERC) volumes. TDP‐43 was associated with reduced hippocampus and ERC volumes but was unrelated to flortaucipir SUVRs (Figure 1). Notably, the TDP‐43 effect size for the hippocampus and ERC was largely consistent across age, i.e., at any age, TDP‐43‐positive participants had significantly smaller hippocampal (∼12% decrease) and ERC (13% decrease) volumes than TDP‐43‐negative participants. However, at older ages, the mean volumes in the whole cohort and, more evidently, in those with dementia moved closer towards the mean of those with TDP‐43 (Figure 1), likely attributable to the increasing frequency of TDP‐43 with age.ConclusionTDP‐43 is a major contributor to the mismatch between flortaucipir uptake and medial temporal lobe volumes in older AD patients where lower medial temporal regional volumes are observed in the presence of less flortaucipir uptake.

AbstractBackgroundFluorodeoxyglucose positron emission tomography (FDG‐PET) is an established modality for differential diagnosis of dementia. Deriving phenotypic signatures from FDG‐PET via machine learning is challenging due to the high dimensional nature of brain images relative to the generally small number of labeled examples available for training, the class imbalance among those labels, and the cooccurrence of multiple pathologies. In this study, we developed a multi‐class, multi‐label framework to address these challenges.MethodA database of clinically acquired PET/CT images from 3,000 unique patients was used to develop a latent space model using matrix decomposition. This model was then applied to images from a separate cohort of Mayo Clinic Alzheimer’s Disease Research Center participants (n=1,745) labeled as cognitively unimpaired (CU) (n=1,436) or with the following potentially co‐occurring phenotypes: Alzheimer’s disease (AD) (n=165), Lewy body dementia (DLB) (n=92), behavioral variant frontotemporal dementia (bvFTD) (n=43), semantic (svPPA) (n=10) and logopenic (lvPPA) (n=13) variant PPA, and posterior cortical atrophy (PCA) (n=17). A k‐nearest neighbors classifier that is robust to these imbalanced and overlapping labels was then trained on these examples. The resulting classifier was evaluated by area under receiver‐operator characteristic curve (ROC‐AUC) via leave one out cross validation, using clinical diagnosis as the gold standard.ResultROC curves and AUC scores for each phenotype are illustrated in Fig. 1a. Because the classifier is based on a k‐nearest neighbors connectivity matrix, it has a convenient graphical representation, where images are nodes and edges are drawn between an image and its set of nearest neighbors in latent space. A self‐organizing force directed graph constructed in this way is illustrated in Fig 1b, highlighting the strong separation of CU and degenerative images, as well as the segregation of each phenotype within the neurodegenerative region of the graph.ConclusionIn this study, we developed a machine learning framework for classification of neurodegenerative disease based on k‐nearest neighbor analysis in a low dimensional latent space projection of FDG‐PET images. By leveraging low‐dimensional representations and k‐nearest neighbors analysis, this framework is robust in multi‐class, multi‐label tasks with strong class imbalance and provides a highly interpretable graphical representation.

AbstractBackgroundWhite matter hyperintensities (WMH) are markers of cerebral small vessel disease and have been associated with higher risk of typical amnestic Alzheimer’s disease (AD), however, little is known about the role of WMHs in atypical variants of AD. This study aimed to determine the frequency and distribution of WMH in two atypical variants of AD: logopenic progressive aphasia (LPA) and posterior cortical atrophy (PCA), and determine whether WMH are associated with demographic and clinical outcomes in these cohorts.MethodSeventy‐five LPA and 39 PCA patients were recruited by the Neurodegenerative Research Group and underwent two‐ dimensional fluid‐attenuated inversion recovery (FLAIR) MRI and Pittsburgh Compound B (PiB) PET. All patients were PiB‐positive. All FLAIR images were segmented using an in‐house software, then were manually edited by a trained imaging analyst. WMH volume was scaled by total intracranial volume to correct for differences in head size. Adjusted WMH burdens were compared between LPA and PCA, and associations were assessed between WMH burden and age, global PiB standard uptake value ratios (SUVRs), and clinical metrics, including tests of general cognitive impairment, executive function, naming, repetition, visuospatial and visuoperceptual function, and parkinsonism. Associations between WMHs and clinical scores were adjusted for age and brain volume to correct for atrophy.ResultTotal WMH burden was similar in PCA and LPA, however, regional differences were identified with PCA having greater subcortical WMH burden in right occipital, parietal, and temporal lobes compared to LPA, and LPA having greater WMH burden in left deep grey and white matter and greater subcortical burden in left parietal lobe compared to PCA. Total WMH burden was greater with older age in both groups but was not related to global PiB SUVR either globally or locally. Greater total WMH was associated with visuoperceptual performance in both cohorts after correcting for atrophy. When assessed regionally, visuoperceptual performance was associated with right subcortical occipital WMH across the cohort.ConclusionRegional distribution of WMHs differ across clinical variants of atypical AD and WMHs may impact clinical outcomes, although further work is needed to determine whether they reflect cerebrovascular disease or regionally specific neurodegenerative changes.

AbstractBackgroundStudies of aging and dementia often acquire brain MRI and PET images, then de‐identify and share these widely with other researchers. It is known that research participants can potentially be re‐identified from brain MRI and CT using face recognition, but this has not been tested with PET imaging.MethodWe generated face reconstruction images of 182 volunteers using their existing amyloid (PiB), tau (FTP), and FDG PET scans from clinical GE PET/CT scanners, and we measured how often commercial automatic face recognition (Microsoft Azure) correctly matched these reconstructions with photographs of their faces. We additionally compared these match rates with face reconstructions from CT (from PET/CT) and MRI. When testing face recognition from PET scans, CT was used only during PET reconstruction, not for face recognition. We also validated our updated mri_reface software, originally developed to de‐identify (“de‐face”) faces in MRI, with new support for PET images. We measured face recognition rates before and after applying mri_reface to PET, CT, and MRI. In addition, we measured effects of de‐facing on biomarker measurements (global and regional SUVR) using 244 ADNI Florbetapir (FBP) amyloid PET scans with two automated pipelines: FreeSurfer 6.0 (PETSurfer) and an in‐house SPM12‐based pipeline.ResultRates of successful face recognition with PET ranged from 32‐41% (Table 1) across tracers. These were lower than MRI (98%) and CT (78%) but still concerning for participant privacy. Using mri_reface reduced face recognition to 0‐3.5% for PET, 5% for CT, and 8% for MRI. mri_reface had little effect on global amyloid SUVR measurements from ADNI PET images (Figure 1). ICC values between original and de‐faced images were 1.00 with biases and median absolute differences < 0.5%, across both pipelines. Effects on individual regions all had ICC values > 0.98 with biases and median absolute differences < 2%.ConclusionFace reconstructions from PET images have sufficient quality for potential re‐identification using publicly available face recognition technology. Research studies should consider using face de‐identification (de‐facing) software on PET images in addition to CT and structural MRI. Our free, automatic mri_reface software greatly reduced the potential for face recognition without substantially affecting SUVR measurements.