AbstractBackgroundWe have previously shown that 3D T1‐weighted and 3D T2‐FLAIR brain MRI can potentially be re‐identified using face recognition, but this has not been tested for other MRI sequences commonly used in imaging studies of aging and dementia. Here we compare match rates of automated face recognition to re‐identify 3D T1‐weighted, 3D T2‐FLAIR, 3D T2‐weighted, 2D Axial T2‐FLAIR, multi‐echo GRE (T2*), diffusion MRI, functional MRI, ASL, and T2* GRE images.MethodWe used a dataset of 182 Mayo Clinic participants each with MRI and face photographs. MRIs used Siemens and GE vendor‐product sequences similar to ADNI3 and ADNI4 protocols. We created face reconstructions (“basic”, with minimal preprocessing, and “advanced” that automatically filled missing face regions using a template) from each MRI sequence, and we trained a Microsoft Azure PersonGroup face recognition classifier to recognize these MRI‐based faces. We then queried the classifier with five face photographs for each participant, and we recorded how often the classifier chose the correct MRI as the top‐ranked match. For image types with match rates > = 10%, we then repeated the experiment after de‐facing each MRI with mri_reface, our leading software that replaces faces in each image (“re‐faces”) with an average face.ResultResults are presented in Table 1. 3D T2‐FLAIR, 3D T1‐weighted, and 3D T2‐weighted images were matched at 96‐98%, and mri_reface reduced all to 8%. Axial T2‐FLAIR and ME‐GRE were matched at 32‐45%, and mri_reface reduced both to 1%. T2*‐GRE matched at 10%, and mri_reface reduced this to 1%. ASL matched at 8%, while dMRI and fMRI matched at only 0‐2%.ConclusionAmong current‐generation vendor‐product sequences, we recommend that de‐facing be applied before research data sharing of T1‐weighted, T2‐weighted, T2‐FLAIR, T2*, and ME‐GRE images. Minimal match rates for current‐generation product fMRI, ASL, and fMRI images suggest that they present low risk of re‐identification and can be shared without de‐facing, but this conclusion should be re‐evaluated if they are acquired without fat‐suppression, with a larger field of view, or if newer developments reduce the current levels of artifacts and distortion around the face.
Modeling the temporal evolution of plasma p‐tau in relation to amyloid beta and tau PET
Abstract INTRODUCTION The timing of plasma biomarker changes is not well understood. The goal of this study was to evaluate the temporal co‐evolution of plasma and positron emission tomography (PET) Alzheimer's disease (AD) biomarkers. METHODS We included 1408 Mayo Clinic Study of Aging and Alzheimer's Disease Research Center participants. An accelerated failure time (AFT) model was fit with amyloid beta (Aβ) PET, tau PET, plasma p‐tau217, p‐tau181, and glial fibrillary acidic protein (GFAP) as endpoints. RESULTS Individual timing of plasma p‐tau progression was strongly associated with Aβ PET and GFAP progression. In the population, GFAP became abnormal first, then Aβ PET, plasma p‐tau, and tau PET temporal meta‐regions of interest when applying cut points based on young, cognitively unimpaired participants. DISCUSSION Plasma p‐tau is a stronger indicator of a temporally linked response to elevated brain Aβ than of tau pathology. While Aβ deposition and a rise in GFAP are upstream events associated with tau phosphorylation, the temporal link between p‐tau and Aβ PET was the strongest. Highlights Plasma p‐tau progression was more strongly associated with Aβ than tau PET. Progression on plasma p‐tau was associated with Aβ PET and GFAP progression. P‐tau181 and p‐tau217 become abnormal after Aβ PET and before tau PET. GFAP became abnormal first, before plasma p‐tau and Aβ PET.
Field strength difference in extent of artifacts induced by CERTAS Plus valves in patients with idiopathic normal pressure hydrocephalus
Background and purpose : Post-shunt MRI is usually performed at 1.5T under the general assumption that shunt-related susceptibility artifacts would be greater at higher field strengths. Purpose The purpose is to show that imaging post-shunt idiopathic normal pressure hydrocephalus (iNPH) patients at 3T is feasible and with reduced artifacts as compared to 1.5T. Methods We manually measured transverse dimensions of artifact at the levels of lateral ventricles, cerebral aqueduct, and cerebellar hemisphere. Areas/volumes of artifacts were calculated assuming an elliptic/ellipsoid shape. Relative extent of shunt-related artifact between field strengths was rated by 3 readers on a 5-point Likert scale. A Wilcoxon Signed Rank Test was used to compare artifact at 1.5T vs 3T for each sequence, with a significance level set at p < 0.05. Results Artifact areas were calculated in 22 iNPH patients; artifacts were on average smaller at 3T vs 1.5T on MPRAGE, DWI, and GRE sequences. On T2 FLAIR and T2 FSE, artifacts at 3T were larger than 1.5T. On the qualitative analysis, artifact effects were less at 3T vs 1.5T on DWI, greater at 3T on T2 FSE, and had mixed results on GRE. Conclusion Our results indicate feasibility of post-shunt imaging with the CERTAS Plus valve at 3T based on shunt-related artifact that is less than or equal in extent to that on 1.5T on most standard clinical imaging sequences. Our findings, corroborated by the qualitative image review, suggest that dedicated clinical imaging sequences for devices may allow for reduction in artifact extent at both 1.5T and 3T.
Objective assessment of patients with idiopathic normal pressure hydrocephalus following ventriculoperitoneal shunt placement using activity-monitoring data: pilot study
OBJECTIVE Idiopathic normal pressure hydrocephalus (iNPH) results in significant morbidity in the elderly with symptoms of dementia, gait instability, and urinary incontinence. In well-selected patients, ventriculoperitoneal shunt (VPS) placement often results in clinical improvement. Most postshunt assessments of patients rely on subjective scales. The goal of this study was to assess the utility of remote activity monitoring to provide objective evidence of gait improvement following VPS placement for iNPH. METHODS Patients with iNPH were prospectively enrolled and fitted with 5 activity monitors (on the hip and bilateral thighs and ankles) that they wore for 4 days preoperatively within 30 days of surgery and for 4 days within 30 days postoperatively. Monitors collected continuous data for number of steps, cadence, body position (upright, prone, supine, and lateral decubitus), gait entropy, and the proportion of each day spent active or static. Data were retrieved from the devices and a comparison of pre- and postoperative movement assessment was performed. The gait data were also correlated with formal clinical gait assessments before and after lumbar puncture and with motion analysis laboratory testing at baseline and 1 month and 1 year after VPS placement. RESULTS Twenty patients fulfilled the inclusion and exclusion criteria (median age 76 years). The baseline median number of daily steps was 1929, the median percentage of the day spent inactive was 70%, the median percentage of the day with a static posture was 95%, the median gait velocity was 0.49 m/sec, and the median number of steps required to turn was 8. There was objective improvement in median entropy from pre- to postoperatively, increasing from 0.6 to 0.8 (p = 0.002). There were no statistically significant differences for any of the remaining variables measured by the activity monitors when comparing the preoperative to the 1-month postoperative time point. All variables from motion analysis testing showed statistically significant differences or a trend toward significance at 1 year after VPS placement. Among the significantly correlated variables at baseline, cadence was inversely correlated with percentage of gait cycle spent in the support phase (contact with ground vs swing phase). CONCLUSIONS This pilot study suggests that activity monitoring provides an early objective measure of improvement in gait entropy after VPS placement among patients with iNPH, although a more significant improvement was noted on the detailed clinical gait assessments. Further long-term studies are needed to determine the utility of remote monitoring for assessing gait improvement following VPS placement.
Predicting amyloid PET and tau PET stages with plasma biomarkers
Abstract Staging the severity of Alzheimer’s disease pathology using biomarkers is useful for therapeutic trials and clinical prognosis. Disease staging with amyloid and tau PET has face validity; however, this would be more practical with plasma biomarkers. Our objectives were, first, to examine approaches for staging amyloid and tau PET and, second, to examine prediction of amyloid and tau PET stages using plasma biomarkers. Participants (n = 1136) were enrolled in either the Mayo Clinic Study of Aging or the Alzheimer’s Disease Research Center; had a concurrent amyloid PET, tau PET and blood draw; and met clinical criteria for cognitively unimpaired (n = 864), mild cognitive impairment (n = 148) or Alzheimer’s clinical syndrome with dementia (n = 124). The latter two groups were combined into a cognitively impaired group (n = 272). We used multinomial regression models to estimate discrimination [concordance (C) statistics] among three amyloid PET stages (low, intermediate, high), four tau PET stages (Braak 0, 1–2, 3–4, 5–6) and a combined amyloid and tau PET stage (none/low versus intermediate/high severity) using plasma biomarkers as predictors separately within unimpaired and impaired individuals. Plasma analytes, p-tau181, Aβ1–42 and Aβ1–40 (analysed as the Aβ42/Aβ40 ratio), glial fibrillary acidic protein and neurofilament light chain were measured on the HD-X Simoa Quanterix platform. Plasma p-tau217 was also measured in a subset (n = 355) of cognitively unimpaired participants using the Lilly Meso Scale Discovery assay. Models with all Quanterix plasma analytes along with risk factors (age, sex and APOE) most often provided the best discrimination among amyloid PET stages (C = 0.78–0.82). Models with p-tau181 provided similar discrimination of tau PET stages to models with all four plasma analytes (C = 0.72–0.85 versus C = 0.73–0.86). Discriminating a PET proxy of intermediate/high from none/low Alzheimer’s disease neuropathological change with all four Quanterix plasma analytes was excellent but not better than p-tau181 only (C = 0.88 versus 0.87 for unimpaired and C = 0.91 versus 0.90 for impaired). Lilly p-tau217 outperformed the Quanterix p-tau181 assay for discriminating high versus intermediate amyloid (C = 0.85 versus 0.74) but did not improve over a model with all Quanterix plasma analytes and risk factors (C = 0.85 versus 0.83). Plasma analytes along with risk factors can discriminate between amyloid and tau PET stages and between a PET surrogate for intermediate/high versus none/low neuropathological change with accuracy in the acceptable to excellent range. Combinations of plasma analytes are better than single analytes for many staging predictions with the exception that Quanterix p-tau181 alone usually performed equivalently to combinations of Quanterix analytes for tau PET discrimination.
Unique brain iron profiles associated with logopenic progressive aphasia and posterior cortical atrophy
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.
Harmonizing machine learning imaging biomarkers – CDESH and the case for harmonization in both MR image and low‐dimensional (scalar) output space
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.
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.
Temporal Association of Imaging Biomarkers of Cerebrovascular Disease and Alzheimer’s Disease
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.
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.