Kamykowski, Michael G.

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.

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.

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.

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.

Abstract Obstructive sleep apnoea (OSA) is associated with an increased risk for cognitive impairment and dementia, which likely involves Alzheimer’s disease pathology. Non-rapid eye movement slow-wave activity (SWA) has been implicated in amyloid clearance, but it has not been studied in the context of longitudinal amyloid accumulation in OSA. This longitudinal retrospective study aims to investigate the relationship between polysomnographic and electrophysiological SWA features and amyloid accumulation. From the Mayo Clinic Study of Aging cohort, we identified 71 participants ≥60 years old with OSA (mean baseline age = 72.9 ± 7.5 years, 60.6% male, 93% cognitively unimpaired) who had at least 2 consecutive Amyloid Pittsburgh Compound B (PiB)-PET scans and a polysomnographic study within 5 years of the baseline scan and before the second scan. Annualized PiB-PET accumulation [global ΔPiB(log)/year] was estimated by the difference between the second and first log-transformed global PiB-PET uptake estimations divided by the interval between scans (years). Sixty-four participants were included in SWA analysis. SWA was characterized by the mean relative spectral power density (%) in slow oscillation (SO: 0.5–0.9 Hz) and delta (1–3.9 Hz) frequency bands and by their downslopes (SO-slope and delta-slope, respectively) during the diagnostic portion of polysomnography. We fit linear regression models to test for associations among global ΔPiB(log)/year, SWA features (mean SO% and delta% or mean SO-slope and delta-slope), and OSA severity markers, after adjusting for age at baseline PiB-PET, APOE ɛ4 and baseline amyloid positivity. For 1 SD increase in SO% and SO-slope, global ΔPiB(log)/year increased by 0.0033 (95% CI: 0.0001; 0.0064, P = 0.042) and 0.0069 (95% CI: 0.0009; 0.0129, P = 0.026), which were comparable to 32% and 59% of the effect size associated with baseline amyloid positivity, respectively. Delta-slope was associated with a reduction in global ΔPiB(log)/year by −0.0082 (95% CI: −0.0143; −0.0021, P = 0.009). Sleep apnoea severity was not associated with amyloid accumulation. Regional associations were stronger in the pre-frontal region. Both slow-wave slopes had more significant and widespread regional associations. Annualized PiB-PET accumulation was positively associated with SO and SO-slope, which may reflect altered sleep homeostasis due to increased homeostatic pressure in the setting of unmet sleep needs, increased synaptic strength, and/or hyper-excitability in OSA. Delta-slope was inversely associated with PiB-PET accumulation, suggesting it may represent residual physiological activity. Further investigation of SWA dynamics in the presence of sleep disorders before and after treatment is necessary for understanding the relationship between amyloid accumulation and SWA physiology.

Abstract Disruption of the default mode network is a hallmark of Alzheimer’s disease, which has not been extensively examined in atypical phenotypes. We investigated cross-sectional and 1-year longitudinal changes in default mode network sub-systems in the visual and language variants of Alzheimer’s disease, in relation to age and tau. Sixty-one amyloid-positive Alzheimer’s disease participants diagnosed with posterior cortical atrophy (n = 33) or logopenic progressive aphasia (n = 28) underwent structural MRI, resting-state functional MRI and [18F]flortaucipir PET. One-hundred and twenty-two amyloid-negative cognitively unimpaired individuals and 60 amyloid-positive individuals diagnosed with amnestic Alzheimer’s disease were included as controls and as a comparison group, respectively, and had structural and resting-state functional MRI. Forty-one atypical Alzheimer’s disease participants, 26 amnestic Alzheimer’s disease participants and 40 cognitively unimpaired individuals had one follow-up functional MRI ∼1–2 years after the baseline scan. Default mode network connectivity was calculated using the dual regression method for posterior, ventral, anterior ventral and anterior dorsal sub-systems derived from independent component analysis. A global measure of default mode network connectivity, the network failure quotient, was also calculated. Linear mixed-effects models and voxel-based analyses were computed for each connectivity measure. Both atypical and amnestic Alzheimer’s disease participants had lower cross-sectional posterior and ventral and higher anterior dorsal connectivity and network failure quotient relative to cognitively unimpaired individuals. Age had opposite effects on connectivity in Alzheimer’s disease participants and cognitively unimpaired individuals. While connectivity declined with age in cognitively unimpaired individuals, younger Alzheimer’s disease participants had lower connectivity than the older ones, particularly in the ventral default mode network. Greater baseline tau-PET uptake was associated with lower ventral and anterior ventral default mode network connectivity in atypical Alzheimer’s disease. Connectivity in the ventral default mode network declined over time in atypical Alzheimer’s disease, particularly in older participants, with lower tau burden. Voxel-based analyses validated the findings of higher anterior dorsal default mode network connectivity, lower posterior and ventral default mode network connectivity and decline in ventral default mode network connectivity over time in atypical Alzheimer’s disease. Visuospatial symptoms were associated with default mode network connectivity disruption. In summary, default mode connectivity disruption was similar between atypical and amnestic Alzheimer’s disease variants, and discriminated Alzheimer’s disease from cognitively unimpaired individuals, with decreased posterior and increased anterior connectivity and with disruption more pronounced in younger participants. The ventral default mode network declined over time in atypical Alzheimer’s disease, suggesting a shift in default mode network connectivity likely related to tau pathology.

AbstractBackgroundPredictive models of tau‐PET accumulation have largely relied on the functional connectome to estimate its spatial patterns in Alzheimer’s disease (AD). Patterns of global functional organization represent a promising avenue as they encapsulate many biological properties relevant to neurodevelopment, cognition and degeneration. We compared approaches based on local connectivity (from task‐free fMRI) and global organization (from FDG‐PET) to predict regional tau‐PET across AD phenotypes.MethodWe included 430 cognitively unimpaired (CU) individuals (214 A+T‐, 216 A+T+) and 62 amnestic mild cognitive impairment (aMCI), 64 amnestic AD, 21 dysexecutive AD (dAD), 41 logopenic progressive aphasia (LPA) and 49 posterior cortical atrophy (PCA) patients (all A+T+) from the Alzheimer’s Disease Research Center, Mayo Clinic Study of Aging, and Neurodegenerative Research Groups (referred to as “Mayo cohort”). The connectivity approach consisted of assessing the relationship between a region’s tau level and its functional connectivity to the region of highest tau (epicenter). The global organization approach consisted of 1) performing a principal component analysis on 3000 FDG‐PET images that did not overlap with the Mayo cohort, 2) projecting FDG‐PET images of the Mayo cohort onto the low‐dimensional space generated in Step 1 and extracting corresponding weighted values across the first 100 components, 3) using these weighted values as predictors of regional tau‐PET in a linear model and determining optimal parameters based on predicted R 2 values, and 4) using the parameters from the model built in Step 3 to predict tau‐PET in the ADNI cohort (n = 86) (Fig.2).ResultThe connectivity‐based approach yielded adjusted R 2 values ranging from 0.19‐0.45 across diagnostic groups (Fig.1). The global organization approach yielded averaged adjusted R 2 and predicted R 2 values of 0.71 and 0.63, respectively, in the Mayo cohort. Out‐of‐sample median R 2 in the ADNI cohort was 0.41 (range 0.10‐0.60) (Fig.3).ConclusionAn unbiased approach based on the global functional organization of the brain and using an imaging modality widely used in clinical settings outperformed a connectome‐based one to predict regional tau‐PET across the AD phenotypic spectrum. Patterns of tau accumulation may better relate to the large‐scale physiology of the brain underlying specific mental functions rather than local connectivity factors.

AbstractBackgroundPredictive models of tau‐PET accumulation have largely relied on the functional connectome to estimate its spatial patterns in Alzheimer’s disease (AD). Patterns of global functional organization represent a promising avenue as they encapsulate many biological properties relevant to neurodevelopment, cognition and degeneration. We compared approaches based on local connectivity (from task‐free fMRI) and global organization (from FDG‐PET) to predict regional tau‐PET across AD phenotypes.MethodWe included 430 cognitively unimpaired (CU) individuals (214 A+T‐, 216 A+T+) and 62 amnestic mild cognitive impairment (aMCI), 64 amnestic AD, 21 dysexecutive AD (dAD), 41 logopenic progressive aphasia (LPA) and 49 posterior cortical atrophy (PCA) patients (all A+T+) from the Alzheimer’s Disease Research Center, Mayo Clinic Study of Aging, and Neurodegenerative Research Groups (referred to as “Mayo cohort”). The connectivity approach consisted of assessing the relationship between a region’s tau level and its functional connectivity to the region of highest tau (epicenter). The global organization approach consisted of 1) performing a principal component analysis on 3000 FDG‐PET images that did not overlap with the Mayo cohort, 2) projecting FDG‐PET images of the Mayo cohort onto the low‐dimensional space generated in Step 1 and extracting corresponding weighted values across the first 100 components, 3) using these weighted values as predictors of regional tau‐PET in a linear model and determining optimal parameters based on predicted R2 values, and 4) using the parameters from the model built in Step 3 to predict tau‐PET in the ADNI cohort (n = 86) (Fig.2).ResultThe connectivity‐based approach yielded adjusted R2 values ranging from 0.19‐0.45 across diagnostic groups (Fig.1). The global organization approach yielded averaged adjusted R2 and predicted R2 values of 0.71 and 0.63, respectively, in the Mayo cohort. Out‐of‐sample median R2 in the ADNI cohort was 0.41 (range 0.10‐0.60) (Fig.3).ConclusionAn unbiased approach based on the global functional organization of the brain and using an imaging modality widely used in clinical settings outperformed a connectome‐based one to predict regional tau‐PET across the AD phenotypic spectrum. Patterns of tau accumulation may better relate to the large‐scale physiology of the brain underlying specific mental functions rather than local connectivity factors.

AbstractFrom a complex systems perspective, clinical syndromes emerging from neurodegenerative diseases are thought to result from multiscale interactions between aggregates of misfolded proteins and the disequilibrium of large-scale networks coordinating functional operations underpinning cognitive phenomena. Across all syndromic presentations of Alzheimer’s disease, age-related disruption of the default mode network is accelerated by amyloid deposition. Conversely, syndromic variability may reflect selective neurodegeneration of modular networks supporting specific cognitive abilities. In this study, we leveraged the breadth of the Human Connectome Project-Aging cohort of non-demented individuals (N = 724) as a normative cohort to assess the robustness of a biomarker of default mode network dysfunction in Alzheimer’s disease, the network failure quotient, across the aging spectrum. We then examined the capacity of the network failure quotient and focal markers of neurodegeneration to discriminate patients with amnestic (N = 8) or dysexecutive (N = 10) Alzheimer’s disease from the normative cohort at the patient level, as well as between Alzheimer’s disease phenotypes. Importantly, all participants and patients were scanned using the Human Connectome Project-Aging protocol, allowing for the acquisition of high-resolution structural imaging and longer resting-state connectivity acquisition time. Using a regression framework, we found that the network failure quotient related to age, global and focal cortical thickness, hippocampal volume, and cognition in the normative Human Connectome Project-Aging cohort, replicating previous results from the Mayo Clinic Study of Aging that used a different scanning protocol. Then, we used quantile curves and group-wise comparisons to show that the network failure quotient commonly distinguished both dysexecutive and amnestic Alzheimer’s disease patients from the normative cohort. In contrast, focal neurodegeneration markers were more phenotype-specific, where the neurodegeneration of parieto-frontal areas associated with dysexecutive Alzheimer’s disease, while the neurodegeneration of hippocampal and temporal areas associated with amnestic Alzheimer’s disease. Capitalizing on a large normative cohort and optimized imaging acquisition protocols, we highlight a biomarker of default mode network failure reflecting shared system-level pathophysiological mechanisms across aging and dysexecutive and amnestic Alzheimer’s disease and biomarkers of focal neurodegeneration reflecting distinct pathognomonic processes across the amnestic and dysexecutive Alzheimer’s disease phenotypes. These findings provide evidence that variability in inter-individual cognitive impairment in Alzheimer’s disease may relate to both modular network degeneration and default mode network disruption. These results provide important information to advance complex systems approaches to cognitive aging and degeneration, expand the armamentarium of biomarkers available to aid diagnosis, monitor progression and inform clinical trials.

AbstractBackgroundNeurite Orientation Dispersion and Density Imaging (NODDI) is becoming increasingly common as multiband acceleration becomes available for clinical scanning, but fitting the model can be extremely slow. Since the release of the original MATLAB implementation in 2012, three alternatives have become popular: Accelerated Microstructure Imaging via Convex Optimization (AMICO), the Microstructure Diffusion Toolbox (MDT), and the CUDA Diffusion Modelling Toolbox (cuDIMOT). They accelerate fitting through reformulation and/or using GPUs, but the implementation differences raise the question of which approach is best. To our knowledge all previous comparisons lacked knowledge of the true values and used the original implementation as the nominally correct one, with in vivo data. We remedy that using simulations with known NODDI values, in addition to in vivo scans.MethodWe created a software phantom with Neurite Density Index (NDI) values linearly increasing from 0.1 to 1.0 in x, Orientation Dispersion (ODI) from 0 to 1 in y, and Isotropic water Volume Fraction (ISOVF) from 0 to 0.9 in z. We simulated five scans using the ADNI3 advanced diffusion MRI protocol. The phantom had 20 voxels in each direction, and the fiber bundle principal direction randomly varied between the simulations. The simulations were fit by the original, AMICO, MDT, and cuDIMOT implementations, and root mean square errors were calculated as shown in Table 1. A discrepancy measure was also calculated using the brain voxels of three in vivo human scans acquired using the same protocol on a Siemens Prisma, substituting the original implementation’s values for the truth since the latter is unknown for in vivo scans.ResultcuDIMOT had the least error (for ODI < 0.8) and was fastest. The original implementation was a close second in accuracy but much slower. cuDIMOT’s ODI underestimation can be corrected with a transfer function in the noise‐free case, but that approach fails with actual scans.ConclusionThe original or cuDIMOT implementations should be used, with cuDIMOT the favorite when at least one CUDA‐capable GPU is available and the number of scans are large. Correcting cuDIMOT’s ODI underestimation using a different initialization or form of Dawson’s integral should be investigated.