AbstractBackgroundThe Mayo Clinic Study of Aging (MCSA) is a longitudinal, population‐based study of residents of Olmsted County, Minnesota, USA. MCSA is releasing de‐identified clinical and imaging data on GAAIN.org to benefit the research community.MethodWe included longitudinal clinical data from all participants of the Mayo Clinic Study of Aging 30‐90 years of age (average 70.0y) through the first 15 years of the study. Clinical data includes age, sex, self‐reported race/ethnicity, APoE4 allele status, cognitive impairment status (normal, MCI, dementia, or other), height, weight, blood pressure, multi‐domain neuropsychological scores, neuroimaging summary measurements, medications, cardiovascular risk factors, social activity, sleep, mood/anxiety, functional measures, and more. For our initial release of imaging data, we selected the first imaging timepoint acquired on 3T GE MRI scanners and amyloid (PiB) PET during the same visit, from all available imaging participants 50‐90 years of age (average 70.7y). Currently, we included 3D T1‐weighted (MPRAGE) MRI (1x1x1.2mm3), T2‐weighted FLAIR MRI (0.85x0.85x3.6mm3), diffusion MRI (2.7mm resolution with 41 directions at b=1000 and 5 b=0), and numeric cortical global SUVR and Centiloid values from PiB PET scans at 40‐60 minutes post‐injection. All images are released in DICOM format and thoroughly de‐identified using CTP while ensuring preservation of necessary fields for analyses. T1‐weighted and T2‐FLAIR MRI were also de‐faced (identifiable face imagery was removed and replaced with an average) using mri_reface version 0.3.3 to further protect participants’ privacy, and all de‐faced images were visually inspected by a trained data scientist for quality assurance.ResultThe first phase of our public data release includes clinical data from 5925 unique participants from the Mayo Clinic Study of Aging (MCSA), each ranging from 1‐12 visits at roughly 15‐month intervals, and three MRI sequences (T1‐w, T2‐w‐FLAIR, and diffusion) from 1802 unique participants. Demographic summaries are given in Table 1.ConclusionThis large, richly characterized, de‐identified clinical and imaging dataset is available to the research community at https://www.gaaindata.org/partner/MCSA, and it is suitable for a wide range of research uses. Additional data will be released over time, including de‐faced amyloid PET images.
Distinct Spatial Patterns for White Matter Microstructure Changes in an Aging and Alzheimer's Disease Cohort
AbstractBackgroundWhite matter (WM) damage is seen with neurodegenerative and cerebrovascular pathologies and contributes to cognitive dysfunction. We hypothesized that there are multiple disease progression patterns in WM microstructural changes related to aging and dementia, and these can be identified using SuStaIn (data‐driven clustering algorithm for disease subtype and stage discovery) on multishell diffusion MRI (NODDI) data. NODDI provides information on cellular tissue architecture – neurite density index (NDI) measures packing density of axons and dendrites, and free water fraction (FWF) measures unrestricted water, e.g., interstitial and cerebrospinal fluid.MethodWe identified 1136 participants (mean age=70 years, 51% male, 40% Aß‐positive, 19% MCI/AD dementia) from the Mayo ADRC and Mayo Clinic Study of Aging. We analyzed NDI and FWF in 55 WM regions. We selected a subset of WM features that monotonically changed with clinical disease progression and then applied the Z‐score SuStaIn algorithm.ResultNDI in seven regions and FWF in six regions survived selection criteria (Figure 1). Noteworthy, the NDI regions were predominantly located in the temporal lobe, while FWF regions were primarily distributed in the frontal and subcortical regions. The SuStaIn algorithm discovered two stable clusters: the NDI‐First (Temporal) and FWF‐First (Frontal) groups (Figure 2). For the NDI‐First group, temporal lobe NDI changed early followed by frontal FWF. The FWF‐First group had the reverse order. There was clear phase separation based on which measure changed first. Stages of the NDI‐First group correlated better with meta‐ROI Tau‐PET SUVR and with MMSE than the FWF‐First group (p<0.01), suggesting that NDI‐First may be indicative of a primarily AD pathway.ConclusionGiven NDI changes are likely degenerative and FWF are suggested to be other processes such as neurodegeneration or inflammatory, the observed changes in the temporal and frontal regions may be due to distinct etiologies, i.e., neurodegenerative in the temporal lobes vs. vascular dysfunction often observed in the frontal lobes. The SuStaIn algorithm showed promise in discovering disease subtypes of WM microstructural changes. Further work is needed to understand the mechanistic basis of these observed changes and integrate WM changes for differential etiology identification and disease staging.
Identifying transition points in AD biomarkers using machine learning: A supplement to the reliable worsening method
AbstractBackgroundBiomarker cutpoints have received significant importance in AD research. Reliable worsening method based on serial biomarker data has been proposed (Jack et al., 2017 Alz & Dem). Advanced machine learning models (ML) and explainable AI provide an opportunity to enhance existing techniques. We aimed to estimate biomarker cutpoints on cross‐sectional data using ML by identifying important transitions in the biomarker values.MethodWe first trained an ML model with the biomarker of interest as an input including other available biomarkers, demographics, and diagnosis, and predicted global cognition z‐score. We then obtained SHAP values (Shapley Additive exPlanations), fit a curve between the SHAP values and biomarker values, and calculated the first derivatives (slope) of the curve. Then, a CUSUM (Cumulative Sum) analysis was used to identify SHAP change points and cross‐referenced to values where the slopes were steepest near the zero‐crossing location of the curve. A zero SHAP is an important indicator of a feature’s progressive importance for model prediction because a feature’s influence on prediction significantly shifts when its SHAP values go from positive to negative or vice versa. However, relying on zero‐crossing alone may be insufficient when the relationship is non‐linear. Consequently, we calculated the second derivatives to identify inflection points (peaks and troughs) along the curve. Lastly, we compiled all the potential cutpoints and performed post‐processing as shown in Fig. 1.ResultThe algorithm identified cutpoints at 1.44 SUVR for all and at 1.42 SUVR for CU participants (Fig. 2). At age 70, both cutpoints estimated the percentage of abnormals at less than 40% (Fig. 3).ConclusionTPE amyloid‐PET SUVR cutpoints for CUs equal the reported RW values (1.42 SUVR). A strength of the TPE is its application on cross‐sectional data making it versatile for various biomarkers across different cohorts.
Plasma and MRI biomarkers predict post‐mortem tau pathology
AbstractBackgroundImaging and plasma markers are used as key indicators of disease for Alzheimer’s disease (AD) but their usefulness in predicting regional tau pathology is relatively understudied. Our objective was to construct predictive models for regional tau pathology measured on postmortem brain tissue using multiple ante‐mortem AD biomarkers. We focused on hippocampal and parietal regions that were immunostained with AT8 and 2E9 that reflect early and advanced aspects of tangle maturity, respectively.MethodsUsing autopsy data from 63 participants enrolled in the Mayo Clinic Study of Aging (MCSA) where the last clinical visit is within three years of death, we analyzed three ante‐mortem plasma markers ‐ phosphorylated tau protein at the threonine 181 residue, glial fibrillary acidic protein, total ß‐amyloid Aß42 relative to Aß40 ratio (Quanterix Simoa HD‐X Analyzer), and three regional MRI markers computed using SPM12 ‐ parietal thickness, parietal volume relative to total intracranial volume ratio, and hippocampal volume adjusted for total intracranial volume. Based on eight separate cross‐sectional linear regression models weighted by the inverse of time from last visit to death, we estimated the multiple and partial R2 predicting post‐mortem pathology measures, including AT8 and 2E9 in both the hippocampal and parietal regions. All models were adjusted for age, sex, and cardiovascular and metabolic conditions score.ResultsThe MRI models (R2 of 0.26‐0.53) had greater predictive power compared to the plasma models (R2 of 0.17‐0.26). The best multiple R2 was 0.53 for prediction of hippocampal AT8 using antemortem hippocampal volume adjusted for total intracranial volume followed by R2 of 0.43 for prediction of hippocampal 2E9 using hippocampal volume adjusted for total intracranial volume. Though the model prediction of parietal AT8 and 2E9 measurements was low (R2 of 0.21‐0.26), parietal thickness exhibited better predictive performance than parietal volume relative to total intracranial volume ratio.ConclusionsThe data we presented here may be utilized to validate biomarkers using regional pathology. Our results support that regional MRI is sensitive to tau pathology. Newer MRI methods should be tested to examine the specificity of newer quantitative measurements to tau pathology.
Comparing regional and global diffusion MRI markers as biomarkers of VCID
AbstractBackgroundQuantifying white matter using diffusion MRI (dMRI) has been proposed for measuring early microstructural tissue changes due to cerebral small vessel disease and aid in quantifying vascular contributions to cognitive impairment and dementia (VCID). Our goal was to compare the usefulness of longitudinal white matter changes in the commonly available diffusion MRI measures for VCID prevention trials.MethodWe included 718 participants over 50 years of age (mean age: 71.1(9.6) years) from the Mayo Clinic Study of Aging, a population‐based sample, with at least two dMRI scans and structural imaging. We computed the commonly available single‐shell dMRI measures (fractional anisotropy, mean diffusivity, and two MarkVCID measures ‐ free water and peak‐width skeletonized mean diffusivity [PSMD]) at each time point. We tested for voxel‐wise associations between dMRI markers and vascular risk measured by cardiovascular metabolic condition (CMC) and observed a region‐specific dependance across all dMRI measures. Using both global and regional dMRI measures, we (i) examined the longitudinal association of dMRI measures with cognition and (ii) computed sample size estimates for a hypothetical clinical trial. We also included white matter hyperintensities (WMH) and our previously proposed composite vascular white matter score (combination of WMH and fractional anisotropy of the genu) as a comparison.ResultVascular risk was associated with all single shell dMRI measures in the genu of the corpus callosum, which we included as a regional dMRI marker for comparison (Figure 1). All dMRI markers correlated with cognitive performance longitudinally (Table 1) and had comparable sample size estimates required for hypothetical VCID clinical trials (Figure 2). Further, global free water and the composite vascular white matter score had the smallest sample size estimates.ConclusionAll commonly used dMRI markers had significant frontal lobe changes due to vascular risk. Both global and regional corpus callosum dMRI markers were sensitive to longitudinal cognitive decline, suggesting their utility in measuring the slowing down of VCID. The composite vascular white matter score, global free water, and WMH show promise as VCID biomarkers. Further work is needed to validate these markers on multiple populations.
A Comparison of DTI‐Based Measures as Predictors of Future Cerebrovascular Degeneration
AbstractBackgroundThere has been a recent proliferation of various quantities based on single shell Diffusion Tensor Imaging (DTI) for capturing overall cerebrovascular health, especially Small Vessel Disease (SVD) as a single number. The existing literature has gaps in the comparison of these measures, so we evaluated them as predictors of both current and future White Matter Hyperintensity (WMH) as a proxy of SVD.MethodUsing baseline 3T MRI examinations (T1, DTI, FLAIR scans) of 598 participants (>60 years, amyloid negative, to limit to those on the SVD pathway) from the Mayo Clinic Study of Aging, we calculated five summary DTI quantities with varying degree of complexity. We estimated their partial R2 from cross‐sectional linear regression models and computed semi‐partial R2 from mixed models predicting WMH volume as a percentage of Total Intracranial Volume (TIV). All models were adjusted for age, sex, and scanner model, and the WMH and based were log‐transformed before fitting. We also conducted subset analyses with high (WMH+) vs. low (WMH‐) baseline WMH load (cut at 1.62%) to compare the performance of the measures in early vs. widespread WMH.ResultAll MD‐based quantities computed on the whole brain, including FWF, performed similarly to each other, but Genu FA underperformed at explaining WMH variance (Table 1). PSMD (skeletonized or atlas based) was more sensitive in WMH+ but dramatically less so in WMH‐. In WMH‐ and overall cohort, FWF and JHU MD (which is a simple averaging of MD in the JHU atlas) performed well.ConclusionAll MD‐based measures performed similarly to each other in the whole dataset but dichotomizing by WMH load showed interesting differences. PSMD measurements are driven by the extreme end (95th percentile) of the MD distribution and therefore more sensitive in WMH+. Genu FA, a regional measure, was not useful for predicting whole brain WMH change. We noticed some variation in how the quantities are affected by scanner model (GE Signa HDxt and Discovery 750) which needs broader investigation. Simply averaging MD over WM avoids additional processing and is more sensitive to WM degeneration that is not (yet) severe.
Chronic Medical Conditions and Dementia Risk: Brain Age Models for Quantifying Impact and Understanding Mechanisms
AbstractBackgroundChronic medical conditions increase dementia risk, and the mechanism may be captured computationally by measures of accelerated brain aging. We hypothesized that quantifying the impact of chronic conditions on GM and WM aging will aid in measuring the differential impact of each condition on brain aging and shed light on the mechanisms through which they increase dementia risk. We specifically focused on microstructural brain changes with Neurite Orientation Dispersion and Density Imaging (NODDI) based on advanced diffusion MRI ‐ neurite density index (NDI), free water fraction (FWF), and orientation dispersion index (ODI).MethodWe identified 1123 participants [913 cognitively unimpaired (CU), 114 with MCI, 96 with AD] from the Mayo Clinic Study of Aging (MCSA) and Mayo ADRC with multi‐shell diffusion MRI. MCSA participants had chronic medical conditions’ diagnosis (based on US Dept of HHS) extracted from the electronic health care records. Using the BorutaSHAP algorithm on NDI, ODI, and ISOVF features from GM and WM separately, we constructed two distinct brain age models for GM and WM. First, we evaluated the usefulness of the predicted age gap (PAG), which is the difference between predicted age by model and chronological age, for both WM and GM NODDI models in CU/MCI/AD. We then developed separate GM and WM models for each condition using samples from condition‐negative participants to compare PAG of GM and WM for that condition.ResultThere was greater PAG in GM (avg = ) than WM (avg = 9.15 yrs) in MCI/AD participants (Fig. 1A). Each condition shifted the GM PAG by 0.10 to 3.56 years and WM PAG by 0.38 to 3.72 years with stroke and congestive heart failure showing the greater PAG overall for GM and WM PAD respectively (Fig. 2). COPD, stroke, chronic kidney disease, hypertension had comparatively greater GM PAG than WM PAG; while hyperlipidemia, congestive heart failure and coronary artery disease had greater WM PAG than GM PAG (Fig. 1B).ConclusionChronic conditions are associated with increased GM and WM Brain aging models can aid in quantifying the impact of chronic conditions on dementia risk and understanding mechanisms.
TDP‐43 association with Subiculum and CA1 Hippocampal Subfield Atrophy in Primary Age‐Related Tauopathy
AbstractBackgroundThe TAR DNA‐binding protein of 43 kDa (TDP‐43) is linked to hippocampal volume loss and faster rates of hippocampal atrophy in Alzheimer's disease and primary age‐related tauopathy (PART). Hence, TDP‐43 is becoming an important player in age‐related neurodegeneration. To advance our understanding of TDP‐43 effect on hippocampus, we conducted an imaging‐pathological study to determine which specific hippocampal subfields are affected by TDP‐43 in cases of PART.MethodOne‐hundred fifteen autopsied cases from Mayo Clinic Alzheimer Disease Center, Study of Aging and Neurodegenerative Research Group were analyzed. All had completed antemortem brain MRIs, underwent neuropathological examination per NIA and Alzheimer’s Association (NIA‐AA) criteria and had additional assessment for TDP‐43 deposition. Only cases meeting criteria for PART (Thal phases 0‐2 and Braak NFT stages I‐IV) were included. Hippocampal subfield segmentation was performed using FreeSurfer‐7. Brain volumes were normalized to total intracranial volume (TIV) when comparing TDP‐43 groups. Statistical analyses including logistic regression analyses were utilized to explore associations between TDP‐43 status and hippocampal subfield volumes.ResultForty‐nine of the 115 cases (43%) with PART were female. Thirty‐seven (32%) cases were TDP‐positive. Mean age at last MRI was 84 years. There were no differences in pathological characteristics between both groups except for presence of hippocampal sclerosis which was more frequent in the TDP(+) group (P <0.001). Whole hippocampal volume was significantly smaller in the TDP(+) group (176 ±31 vs. 188 ±27, p=0.049). After stratification of hippocampus into head, body, and tail, TDP‐positive cases had smaller volumes of hippocampal body (p=0.031) and tail (p=0.021) but did not reach significance for hippocampal head (p=0.129). When regions were subdivided into hippocampal subfields, TDP(+) cases showed smaller volumes of molecular layer of the hippocampal body (p=0.046), CA1 body (p=0.031), subiculum body (p=0.023), and presubiculum body (p=0.046) compared to TDP(‐) cases. Logistic regression analysis adjusting for TIV confirmed these associations.ConclusionThis study underscores the association between TDP‐43 and the hippocampus and provides further evidence that TDP‐43 targets multiple regions and subfields of the hippocampus in patients with PART. Intriguingly, regions targeted are those where TDP‐43 deposition is microscopically observed, suggesting a direct “hit” by TDP‐43 deposition.
Variational autoencoder latent space as a robust and pragmatic clinical classification tool for dementia
Abstract Background Many proposed clinical decision support systems (CDSS) require multiple disparate data elements as input, which makes implementation difficult, and furthermore have a black‐box nature leading to low interpretability. Fluorodeoxyglucose Positron Emission Tomography (FDG‐PET) is an established modality for the diagnosis of dementia, and a CDSS that uses only an FDG‐PET image to produce a reliable and understandable result would ease both of these challenges to clinical application. Method A deep variational autoencoder (VAE) was used to extract a latent representation of each image through prior training from FDG‐PET brain images (n=2000). This unsupervised VAE has a novel graph convolutional architecture that makes it applicable to masked template space images. A logistic regression model was used to classify each image. A parametric study of the latent space revealed the imaging features that were used in the logistic regression model to differentiate each class. A separate dataset of participants labeled with their clinical diagnosis (n=1239) was used to assess the logistic regression model’s diagnostic accuracy of cognitively unimpaired (CU) (n=679) and between common dementia subtypes: Alzheimer’s disease (AD) (n=310), Lewy body dementia (DLB) (n=151), and behavioral variant frontotemporal dementia (bvFTD) (n=99). The logistic regression classifier was evaluated using receiver‐operator characteristic (ROC‐AUC) curves. Result ROC‐AUC curves for CU and each dementia subtype are illustrated in Figure 1a. A k nearest neighbors model was used to develop a graphical representation of the latent space where nodes are images and edges are drawn between nearest neighbors, illustrated in Figure 1b. The graph of the latent space demonstrates a comparison between dementia subtypes and CU images. Conclusion In this study, we developed a proof‐of‐concept VAE that can achieve differential diagnosis with 84.8% balanced accuracy. The model architecture provides visual interpretability, and shows that the hippocampi, parietal lobes, and cerebellum, were useful for distinguishing AD, the frontal lobe, posterior cingulate cortex, and hippocampi were useful for distinguishing DLB, and prefrontal pole, subcortical structures, and occipital lobe were useful for distinguishing bvFTD. These regions of metabolism align with clinical pathology and metabolic patterns.
Implementation and Validation of Face De‐Identification (de‐facing) in ADNI‐4
AbstractBackgroundRecent advances in automatic face recognition have increased the risk that de‐identified research imaging data could be re‐identified from face imagery in brain scans.MethodAn ADNI committee of independent imaging experts evaluated 11 published techniques for face‐deidentification (“de‐facing”) and selected four algorithms (FSL‐UK Biobank, HCP/XNAT, mri_reface, and BIC) for formal testing using 183 longitudinal scans of 61 racially and ethnically diverse ADNI participants, evaluated by their facial feature removal on 3D rendered surfaces (confirming sufficient privacy protection) and by comparing measurements from ADNI routine image analyses on unmodified vs. de‐faced images (confirming negligible side effects on analyses).ResultEffects of the de‐facing methods on FreeSurfer at UCSF (Figure 1‐top), ASHS at University of Pennsylvania (Figure 1 bottom‐left), and longitudinal atrophy at Mayo Clinic (Figure 1 bottom‐right) were comparable with each other and sufficiently small vs. unmodified scans to recommend de‐facing in ADNI‐4. The FSL and mri_reface methods were more thorough in removing facial features, partly because these also removed ears (Figure 2‐top). mri_reface was recommended for ADNI‐4 primarily because fewer (zero) images had modified brain (Figure 2‐bottom), and ADNI leadership approved de‐facing with mri_reface. ADNI performed a second validation of its latest version on another diverse sample of 100 ADNI‐3 participants. These analyses included: boundary shift interval at University College London, where average differences in measurements were all <0.05%; cortical white matter hyperintensities at UC Davis, where they found “minimum impact” (R>0.97); and others.ConclusionADNI‐4 de‐faces all indicated image types (Table 1) before subsequent pre‐processing, analyses, and public release. Sites upload images to the central repository (LONI), where they are queued to the MRI Core at Mayo Clinic for de‐facing. Image types requiring de‐facing are de‐faced using mri_reface. Trained analysts inspect de‐faced images to confirm complete removal of the face and complete non‐modification of the brain, and any flagged images are re‐processed manually until issues are resolved. De‐faced images are converted back to DICOM and re‐uploaded for further distribution. Image types not requiring de‐facing are indicated to proceed directly to subsequent steps. In summary, ADNI‐4 has adopted mri_reface because it successfully de‐faces without substantially affecting brain measurements.