Meyer, Nolan

Abstract Background Diffusion MRI (dMRI) measures are increasingly being used as biomarkers for aging and dementia in multi‐center studies. Measurements in longitudinal studies are often acquired on different systems, which can result in systematic differences of machine origin (e.g. vendor, hardware, algorithms). This may warrant harmonization prior to pooling multi‐site data. However, harmonization algorithms may introduce a risk of removing valuable information content and reducing sensitivity to physiological factors (e.g. sex, age, disease). Here, we tested the hypothesis that protocol standardization may minimize machine‐origin differences, obviating a need for post‐acquisition harmonization in a two‐site (Mississippi‐MS and Minnesota‐MN) cohort of the longitudinal Mayo Clinic Study of Aging (MCSA). Methods Two data sets were analyzed: first, a crossover dataset with 80 participants (scanned within 2 days) bridging two different protocols on GE and Siemens scanners; second, participants scanned using standardized protocols from the ongoing MCSA ( N MN‐MCSA =204; N MS‐MCSA =52). Parameters are listed in Figure 1. For crossover data, analyses were repeated with GE data resampled to 2‐mm isotropic resolution, and discarding Siemens‐only b‐shells to artificially reduce protocol disparities. Analysis encompassed visualization of overlain subjectwise histogram‐derived densities obtained from whole‐brain masks for dMRI measures; and median dMRI measures from within the genu of the corpus callosum (GCC) with nonparametric statistical testing. When available, both DTI and NODDI measures were compared. Results Pronounced visual distributional differences seen between MD and FA in crossover data were reduced with resampling of GE data and including only b‐shells shared across acquisitions. For the resampled/reshelled data, significant differences remained for all GCC‐FA comparisons, but not all GCC‐MD (Figure 2). For the ongoing MCSA data from two sites, MD and FA whole‐brain densities were more closely aligned, with GCC‐FA and GCC‐MD showing no significant differences (Figure 3). NODDI measures showed more discordant densities than DTI measures, including shifts in most probable values and differences in distributional arrangement. Conclusions Our work provides further evidence supporting protocol and processing standardization for minimizing machine‐origin differences prior to considering harmonization. Differences in dMRI metrics seen when comparing across unstandardized protocols were reduced when protocols were standardized. NODDI metrics had more apparent sensitivity to site than DTI measures.

Abstract Background Multi‐site aging and dementia studies are challenged by data heterogeneity caused by scanner (manufacturer and protocol) related variability. Quantitative analysis of white matter hyperintensity (WMH) from pooled samples, a traditional marker of cerebrovascular disease, is sensitive to machine variability. Our goal was to compare Deming (symmetric) linear regression and ComBat, a common method used for harmonizing neuroimaging measurements, for harmonizing WMH measurements on participants scanned on two different scanners (GE and Siemens) and protocols (2D FLAIR on GE and 3D FLAIR on Siemens). Methods Cognitively unimpaired (CU) ( n =53, age (mean (SD))=67 (15) years) and cognitively impaired (CI; MCI/dementia) ( n =34, age (mean (SD))=82 (9) years) participants from the Mayo Clinic Study of Aging and Mayo ADRC were scanned on GE and Siemens scanners within a week of each other (typically done on the same day). WMH as a percentage of total intracranial volume was measured using multispectral segmentation of FLAIR and T1‐weighted images (Gunter et. al. AAIC 2023). We used Deming regression as it includes errors in both input dimensions. ComBat was applied using age, sex, and diagnosis as covariates, with scanner vendor and protocol as the batch variable. We assessed the inter‐scanner similarity using Maximum Mean Discrepancy (MMD) analysis, where a value of zero indicates identical distributions, p >0.05 indicates no significant difference under H 0 . (Gretton et al., JMLR 2012). Results Our results showed that the distribution of WMH measurements, which were different in the CU group, were more similar after Deming than ComBat harmonization (Figure 1). This was confirmed by the lower MMD in the CU group from Deming regression (‐0.01 ( p =0.76) from 0.03 ( p =0.05)) compared to ComBat (0.04 ( p =0.02))(Figure 1). However, ComBat and Deming were unable to create distributions as compatible as original measurements for the CI group. Conclusion Using a unique dataset of the same participants scanned close in time on different scanner and protocols, we found that (1) MMD measurements were useful to test the need/improvement due to harmonization; (2) ComBat was not substantially better at creating compatible distributions between scanner vendors than Deming regression, and both were sometimes inadequate. Future work will require alternative models to linear scale and shift transformations.

Abstract Background Diffusion MRI (dMRI) measures are increasingly being used as biomarkers for aging and dementia in multi‐center studies. Measurements in longitudinal studies are often acquired on different systems, which can result in systematic differences of machine origin (e.g. vendor, hardware, algorithms). This may warrant harmonization prior to pooling multi‐site data. However, harmonization algorithms may introduce a risk of removing valuable information content and reducing sensitivity to physiological factors (e.g. sex, age, disease). Here, we tested the hypothesis that protocol standardization may minimize machine‐origin differences, obviating a need for post‐acquisition harmonization in a two‐site (Mississippi‐MS and Minnesota‐MN) cohort of the longitudinal Mayo Clinic Study of Aging (MCSA). Methods Two data sets were analyzed: first, a crossover dataset with 80 participants (scanned within 2 days) bridging two different protocols on GE and Siemens scanners; second, participants scanned using standardized protocols from the ongoing MCSA ( N MN‐MCSA =204; N MS‐MCSA =52). Parameters are listed in Figure 1. For crossover data, analyses were repeated with GE data resampled to 2‐mm isotropic resolution, and discarding Siemens‐only b‐shells to artificially reduce protocol disparities. Analysis encompassed visualization of overlain subjectwise histogram‐derived densities obtained from whole‐brain masks for dMRI measures; and median dMRI measures from within the genu of the corpus callosum (GCC) with nonparametric statistical testing. When available, both DTI and NODDI measures were compared. Results Pronounced visual distributional differences seen between MD and FA in crossover data were reduced with resampling of GE data and including only b‐shells shared across acquisitions. For the resampled/reshelled data, significant differences remained for all GCC‐FA comparisons, but not all GCC‐MD (Figure 2). For the ongoing MCSA data from two sites, MD and FA whole‐brain densities were more closely aligned, with GCC‐FA and GCC‐MD showing no significant differences (Figure 3). NODDI measures showed more discordant densities than DTI measures, including shifts in most probable values and differences in distributional arrangement. Conclusions Our work provides further evidence supporting protocol and processing standardization for minimizing machine‐origin differences prior to considering harmonization. Differences in dMRI metrics seen when comparing across unstandardized protocols were reduced when protocols were standardized. NODDI metrics had more apparent sensitivity to site than DTI measures.

ABSTRACT Purpose Magnetic resonance elastography (MRE) of the heart has predominantly utilized breath‐held acquisitions with limited anatomic coverage. This work investigates the feasibility of 3D, free‐breathing cardiac MRE. Methods A 3D hybrid radial and EPI acquisition is utilized and combined with retrospective binning of k ‐space via physiologic monitoring, iterative reconstruction, MRE processing and inversions, and post‐processing specialized for cardiac imaging. Results Feasibility of free‐breathing stiffness mapping of the left ventricular myocardium was demonstrated in 11 participants. Myocardial stiffness estimates were obtained throughout the left ventricle for seven temporal states of the cardiac cycle. Stiffness estimates were presented in participant‐specific images, 17‐segment bullseye plots, and stiffness plots throughout the cardiac cycle. Results indicate increased stiffness in systolic states, with localized stiffness heterogeneity observed between participants, illustrating a key benefit of 3D cardiac MRE acquisition. Mean stiffness estimates obtained through 11 subjects were 4.44 ± 0.68 kPa (begin‐systole), 5.74 ± 0.94, 5.42 ± 0.87 (end‐systole), 3.84 ± 0.71 (begin‐diastole), 3.51 ± 0.74, 3.54 ± 0.52, and 3.55 ± 0.59 kPa (end‐diastole). Statistical testing indicated significant differences in stiffness across cardiac states with one‐way, repeated‐measures ANOVA ( p < 0.001) and Friedman's test ( p < 0.001). Experiments in scanning three participants with and without applied motion indicated that stiffness changes were more robustly detected in systole than in diastole, although patterns were subject‐specific. Comparison of free‐breathing stiffness measurements against reference breath‐held measurements yielded similar systolic stiffness estimates. Conclusions This study establishes the feasibility of free‐breathing, 3D cardiac MRE in healthy volunteers. With the current framework, stiffness estimation appears more robust in systole, and the reliability of diastolic estimates is limited and subject‐specific.

Abstract: Objectives: Locally low-rank (LLR) denoising of functional magnetic resonance imaging (fMRI) time series image data is extended to multi-echo (ME) data. The proposed method extends the capabilities of non-physiologic noise suppression beyond single-echo applications with a dedicated ME algorithm. Materials and Methods: Following an institutional review board (IRB) approved protocol, resting-state fMRI data were acquired in 7 healthy subjects. A compact 3T scanner enabled whole-brain acquisition of multiband ME fMRI data at high spatial resolution (1.4 × 1.4 × 2.8 mm3) with a 1810 ms repetition time (TR). Image data were denoised with ME-LLR preceding functional processing. The results of connectivity maps generated from denoised data were compared with maps generated with equivalent processing of non-denoised images. To assess ME-LLR as a method to reduce scan time, comparisons were made between maps computed from image data with full and retrospectively truncated durations. Assessments were completed with seed-based connectivity analyses using echo-combined image data. In a feasibility assessment, nondenoised and denoised full-duration echo-combined data were equivalently processed with independent component analysis (ICA) and compared. Results: ME-LLR denoising yielded strengthened resting-state network connectivity maps after nuisance regression and seed-based connectivity analysis. In assessing ME-LLR as a scan reduction mechanism, maps generated from denoised data at half scan time showed comparable quality with maps generated from full-duration, non-denoised data, at both single subject and group levels. ME-LLR substantially increased temporal signal-to-noise ratio (tSNR) for image data respective to each individual echo and for image data after nuisance regression. Among echo-specific image volumes, increases in tSNR yielded by ME-LLR were most pronounced for image data with the longest echo time and thereby lowest SNR. ICA showed resting-state networks consistently identified between non-denoised and denoised data, with clearer demarcation of networks for ME-LLR. Conclusions: ME-LLR is demonstrated to suppress non-physiologic noise, enhance functional connectivity map quality, and could potentially facilitate scan time reduction in ME-fMRI.

In blood-oxygen-level-dependent (BOLD)-based resting-state functional (RS-fMRI) studies, usage of multi-echo echo-planar-imaging (ME-EPI) is limited due to unacceptable late echo times when high spatial resolution is used. Equipped with high-performance gradients, the compact 3T MRI system (C3T) enables a three-echo whole-brain ME-EPI protocol with smaller than 2.5 mm isotropic voxel and shorter than 1 s repetition time, as required in landmark fMRI studies. The performance of the ME-EPI was comprehensively evaluated with signal variance reduction and region-of-interest-, seed- and independent-component-analysis-based functional connectivity analyses and compared with a counterpart of single-echo EPI with the shortest TR possible. Through the multi-echo combination, the thermal noise level is reduced. Functional connectivity, as well as signal intensity, are recovered in the medial orbital sulcus and anterior transverse collateral sulcus in ME-EPI. It is demonstrated that ME-EPI provides superior sensitivity and accuracy for detecting functional connectivity and/or brain networks in comparison with single-echo EPI. In conclusion, the high-performance gradient enabled high-spatial-temporal resolution ME-EPI would be the method of choice for RS-fMRI study on the C3T.

Abstract Objective . Interleaved reverse-gradient fMRI (RG-fMRI) with a point-spread-function (PSF) mapping-based distortion correction scheme has the potential to minimize signal loss in echo-planar-imaging (EPI). In this work, the RG-fMRI is further improved by imaging protocol optimization and application of reverse Fourier acquisition. Approach . Multi-band imaging was adapted for RG-fMRI to improve the temporal and spatial resolution. To better understand signal dropouts in forward and reverse EPIs, a simple theoretical relationship between echo shift and geometric distortion was derived and validated by the reliable measurements using PSF mapping method. After examining practical imaging protocols for RG-fMRI in three subjects on both a conventional whole-body and a high-performance compact 3 T, the results were compared and the feasibility to further improve the RG-fMRI scheme were explored. High-resolution breath-holding RG-fMRI was conducted with nine subjects on the compact 3 T and the fMRI reliability improvement in high susceptibility brain regions was demonstrated. Finally, reverse Fourier acquisition was applied to RG-fMRI, and its benefit was assessed by a simulation study based on the breath-holding RG-fMRI data. Main results . The temporal and spatial resolution of the multi-band RG-fMRI became feasible for whole-brain fMRI. Echo shift measurements from PSF mapping well estimated signal dropout effects in the EPI pair and were useful to further improve the RG-fMRI scheme. Breath-holding RG-fMRI demonstrated improved fMRI reliability in high susceptibility brain regions. Reverse partial Fourier acquisition omitting the late echoes could further improve the temporal or spatial resolution for RG-fMRI without noticeable signal degradation and spatial resolution loss. Significance . With the improved imaging scheme, RG-fMRI could reliably investigate the functional mechanisms of the human brain in the temporal and frontal areas suffering from susceptibility-induced functional sensitivity loss.

Objective This study investigates a locally low-rank (LLR) denoising algorithm applied to source images from a clinical task-based functional MRI (fMRI) exam before post-processing for improving statistical confidence of task-based activation maps. Methods Task-based motor and language fMRI was obtained in eleven healthy volunteers under an IRB approved protocol. LLR denoising was then applied to raw complex-valued image data before fMRI processing. Activation maps generated from conventional non-denoised (control) data were compared with maps derived from LLR-denoised image data. Four board-certified neuroradiologists completed consensus assessment of activation maps; region-specific and aggregate motor and language consensus thresholds were then compared with nonparametric statistical tests. Additional evaluation included retrospective truncation of exam data without and with LLR denoising; a ROI-based analysis tracked t-statistics and temporal SNR (tSNR) as scan durations decreased. A test-retest assessment was performed; retest data were matched with initial test data and compared for one subject. Results fMRI activation maps generated from LLR-denoised data predominantly exhibited statistically significant ( p = 4.88×10 –4 to p = 0.042; one p = 0.062) increases in consensus t-statistic thresholds for motor and language activation maps. Following data truncation, LLR data showed task-specific increases in t-statistics and tSNR respectively exceeding 20 and 50% compared to control. LLR denoising enabled truncation of exam durations while preserving cluster volumes at fixed thresholds. Test-retest showed variable activation with LLR data thresholded higher in matching initial test data. Conclusion LLR denoising affords robust increases in t-statistics on fMRI activation maps compared to routine processing, and offers potential for reduced scan duration while preserving map quality.