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.
Large‐scale Normative Modeling of Brain Microstructure over the Human Lifespan
AbstractBackgroundNormative models (NM) of brain metrics based on large, diverse populations offer novel strategies to detect individual brain abnormalities. To create an age‐dependent statistical model of brain microstructure over the human lifespan, we built the largest multi‐site NM of white matter (WM) diffusion tensor imaging (DTI) metrics based on 54,591 subjects. We used state‐of‐the‐art tools to adjust for site‐dependent effects.MethodWe used Hierarchical Bayesian Regression (HBR) to determine the age trajectory and lifespan centile curves (ages 3‐100 years) by merging data from different sections of the lifespan. We analyzed 19 international public datasets with the ENIGMA‐DTI protocol. Mean fractional anisotropy (FA), mean, axial, and radial diffusivity (MD, AxD, RD) were extracted for 21 bilateral ROIs from the JHU‐WM atlas and the whole WM skeleton. Regressions were run with each metric per ROI as a function of age and sex, and the DTI protocol as the batch effect (37 protocols). Z‐scores were derived for 81 patients with Alzheimer’s disease (AD; mean age: 77y±8.4, 46M/35F) and 225 MCI participants (mean age: 75.1y±8.1, 127M/98F) and extreme deviations (Z>|2|) were quantified for each condition. ROI‐wise areas under the ROC curve were calculated to determine the classification accuracy of the Z‐scores.ResultThe average FA for the overall WM skeleton reached its peak value at 27 years, with minima for MD, RD and AxD at 47, 42 and 54 years, respectively. We found extreme FA deviations (Z>2) for MCI (AUC=0.67) and AD (AUC=0.66) in the genu of the corpus callosum, extreme RD deviations (Z<2) in the tapetum (AUC=0.73) and Extreme AxD deviations (Z>2) for MCI in the superior longitudinal fasciculus (AUC=0.65). See Figures 1 & 2.ConclusionNM of brain microstructure with HBR yielded lifespan trajectories for widely‐used DTI metrics. By pooling large scale multi‐site dMRI to create a statistical reference, single subject deviations are detectable across the lifespan for AD and MCI. These NM will be publicly available to the research community.
Relationships between PET and blood plasma biomarkers in corticobasal syndrome
AbstractBackgroundCorticobasal syndrome (CBS) is a progressive neurodegenerative syndrome that can result from several underlying neuropathologies, including Alzheimer’s disease (AD). Little is known about the utility of blood plasma measures to predict PET biomarker confirmed AD in CBS.MethodEighteen CBS patients (8 PET Aß+; 10 PET Aß‐) and 8 cognitively unimpaired Aß‐ individuals (CU) were recruited by the Neurodegenerative Research group from the Department of Neurology, Mayo Clinic, Rochester, MN, underwent Aß (Pittsburgh Compound‐B) and tau (18F‐flortaucipir) PET and provided a blood sample. Blood plasma analysis was performed to assess p‐tau 181, total tau, Aß 42/40, neurofilament light chain (NfL) and glial fibrillary acidic protein (GFAP) using the Quanterix Simoa platform. Plasma concentrations were compared between Aß‐, and Aß+ CBS and controls using Kruskal‐Wallis corrected for multiple comparisons using Dunn’s test, with effect size measured with area under the receiver operating characteristic curve (AUROC) analysis. Spearman correlations were calculated to assess relationships between plasma concentrations and Aß and tau PET uptake.ResultThe CBS Aß+ group showed reduced Aß42/40 ratio, with elevated p‐tau 181, GFAP and NfL plasma concentrations compared to CU, while the CBS Aß‐ group only showed elevated NfL concentration compared to CU (Figure 1). Both p‐tau 181 (AUROC=0.95) and GFAP (AUROC=0.91) plasma concentrations were able to differentiate CBS Aß‐ from CBS Aß+. Plasma p‐tau 181 and GFAP concentrations showed a positive correlation to both Aß and tau PET uptake, and the Aß42/40 ratio showed a negative correlation to Aß PET uptake (Figure 2). No correlations were observed between NfL and Aß and tau PET.ConclusionThis study supports the use of plasma p‐tau 181 and GFAP to detect AD in CBS, while NfL shows potential as a nonspecific disease biomarker of CBS regardless of underlying pathology. Plasma concentrations of these biomarkers may potentially be incorporated into clinical trials targeting CBS patients with underlying AD pathology.
Association between temporal lobe cortical NODDI measures and memory function
AbstractBackgroundCortical structural changes in the temporal lobes including global and hippocampal volume loss and decrease in cortical thickness are well‐associated findings linked to cognitive decline in MCI and dementia patients. Neurite orientation dispersion and density imaging (NODDI) is a multishell diffusion MRI method enabling assessment of tissue microstructural integrity at the axonal and dendritic level. We wished to determine whether changes in temporal cortex orientation dispersion index (ODI) & neurite density index (NDI) correlate with decreases in memory function.MethodThe MCALT_ADIR122 cortical atlas was used for ROI designation. NODDI imaging was obtained on a 3T Siemens Prisma scanner with 2 mm isotropic resolution. Participants (ages 60‐95 years) were cognitively unimpaired (CU) or had mild cognitive impairment (MCI). CU & MCI subjects were enrolled in the Mayo Clinic Study of Aging (MCSA) and underwent neuropsychological testing including 3 memory tests from which a composite memory z score was calculated. A generalized additive regression model was used to estimate z scores as a smooth function of ODI, NDI and cortical thickness or volume with adjustment for age, gender, education, MCSA cycle number, WMH and intracranial volume.ResultThe Table displays demographic data. Figures 1 and 2 show the relationship of NDI and ODI measures to memory z scores across participants in individual ROIs both before and after adjustment for cortical thickness (“structural MRI”) as well as for confounds listed above. Both higher cortical thickness and higher NDI values in the hippocampus and entorhinal cortex were associated with significantly higher memory z‐scores (p<0.001) (Figure 1). Higher ODI values in addition to greater cortical thickness in temporal lobe regions, including within the hippocampus, parahippocampal gyrus, entorhinal cortex and inferior temporal gyrus, were associated with higher memory z‐scores (Figure 2). These significant associations (p <=0.01) were maintained after adjusting for cortical thickness in all these regions except entorhinal cortex for ODI.ConclusionNODDI‐derived measures of microstructural integrity (NDI & ODI) within the temporal cortex may provide information about memory decline, beyond what cortical thickness analysis provides, and may be potentially useful for monitoring therapeutic effects of AD immunotherapy in the future.
Measuring the Potential Risk of Re‐identification of Imaging Research Participants from Open‐Source Automated Face Recognition Software
AbstractBackgroundResearch subjects can potentially be re‐identified from de‐identified MRI, CT, and PET brain scans with up to 98% accuracy using Microsoft Azure’s cloud‐based commercial facial recognition software. This showed the need to “de‐face” publicly shared research brain scans. Subsequently, Microsoft has begun restricting its face recognition services, intending to prevent misuse. This study tests a variety of popular open‐source computer vision and facial recognition software packages to measure the re‐identification risk with only unrestricted free software.MethodThis study used brain MRI from Siemens 3T scanners and face photographs from 182 participants (ages 30‐90+) in the Mayo Clinic Study of Aging. Each participant had photos taken with an iPad, which were cropped loosely around the head and converted to black and white. A total of 81 renders (2D reconstruction images; Figure 1) were made from each 3D T2‐FLAIR MRI, simulating various camera and lighting positions. These face renders approximate black and white photos, allowing use with software packages developed for matching faces/features between photographs. Open‐source Python packages (OpenCV, DeepFace, face_recognition, scikit‐learn) and their included pre‐trained models (designed for photographs) were used to match participant photographs to the labelled library of MRI‐based face renders. A correct match was recorded when the top‐scoring render or set of renders automatically chosen by the algorithm was from the MRI of the correct participant.ResultFeature matching with ORB and SWIFT (both from OpenCV) produced no correct matches, despite a prior publication reporting moderate success. The best match rates (Figure 2) were 58‐59% accuracy using face_recognition or DeepFace packages with support vector machines (SVMs). Without SVMs, these packages performed modestly worse (51%‐55%).ConclusionWhile the highest match rates obtained by widely available open‐source packages (roughly 60%) are much lower than those of a limited‐access commercial product (Microsoft Azure, with 98%), they show that actors with only access to no‐cost readily available software could still re‐identify a research subject from a brain scan most of the time. This further highlights the need for image refacing software, such as mri_reface, to protect the privacy and identifies of research volunteers.
Comparing de‐facing software MiDeFace and mri_reface according to privacy protection and effects on downstream measurements
AbstractBackgroundSeveral programs are available for automatically “de‐facing” brain images. These are designed to prevent the potential use of face imagery to re‐identify research participants. Previous works have compared available software and mri_reface is frequently among the top methods. The “Minimally Invasive DeFacing of MRI images” (MiDeFace) program was more recently added to FreeSurfer, but minimal validation of MiDeFace has been published.MethodWe ran mri_reface 0.3.3 and MiDeFace 7.3.2 on two datasets: A) 207 participants with face photographs and 3D T1‐weighted and 3D T2‐weighted‐FLAIR MRI from the Mayo Clinic Study of Aging (all unimpaired, Siemens); and B) 300 participants from ADNI3 with 3D T1‐weighted MRI (50 unimpaired + 50 with clinical AD, from each of 3 vendors). For dataset A, we tested the rate at which Microsoft Azure face recognition matched each participant’s photographs with the correct MRI, comparing unmodified images, those from mri_reface, and those from MiDeFace. For dataset B, we compared measured regional brain volume and cortical thickness from original vs. de‐faced images with SPM12, FreeSurfer 7.4.1, and the FSL‐UK Biobank pipeline. For both datasets, images were visually rated to assess removal of the face and retention of brain tissue.ResultAutomated face recognition match rates were lower with mri_reface (17/207 8.2%) than with MiDeFace (24/207 11.6%), indicating stronger identity protection with mri_reface, but both were greatly reduced vs. unmodified images (203/207 98.1%). Example images with both programs are shown in Figure 1. From visual ratings (Table 1), rates of complete success (modification of the entire face while retaining the entire brain) were much higher with mri_reface, but the great majority of failures from MiDeFace were subtle and occurred in relatively less recognizable face regions. The effects of both programs on downstream imaging measurements (Figure 2) were small and very comparable, but MiDeFace was favored mildly for SPM12 and for FSL, while mri_reface was more favored for FreeSurfer.ConclusionAltogether, mri_reface more reliably removed the complete face, leading to smaller rates of potential re‐identification (8.2% vs 11.6%). Effects on downstream image measurements were small and comparable with both methods.
A Large Public Release of Clinical and Imaging Data from the Mayo Clinic Study of Aging
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.