Precision mapping of the default network reveals common and distinct (inter) activity for autobiographical memory and theory of mind.

The default network is widely implicated as a common neural substrate for self-generated thought, such as remembering one's past (autobiographical memory) and imagining the thoughts and feelings of others (theory of mind). Findings that the default network comprises subnetworks of regions, some commonly and some distinctly involved across processes, suggest that one's own experiences inform their understanding of others. With the advent of precision functional MRI (fMRI) methods, however, it is unclear if this shared substrate is observed instead due to traditional group analysis methods. We investigated this possibility using a novel combination of methodological strategies. Twenty-three participants underwent multi-echo resting-state and task fMRI. We used their resting-state scans to conduct cortical parcellation sensitive to individual variation while preserving our ability to conduct group analysis. Using multivariate analyses, we assessed the functional activation and connectivity profiles of default network regions while participants engaged in autobiographical memory, theory of mind, or a sensorimotor control condition. Across the default network, we observed stronger activity associated with both autobiographical memory and theory of mind compared to the control condition. Nonetheless, we also observed that some regions showed preferential activity to either experimental condition, in line with past work. The connectivity results similarly indicated shared and distinct functional profiles. Our results support that autobiographical memory and theory of mind, two theoretically important and widely studied domains of social cognition, evoke common and distinct aspects of the default network even when ensuring high fidelity to individual-specific characteristics.NEW & NOTEWORTHY We used cutting-edge precision functional MRI (fMRI) methods such as multi-echo fMRI acquisition and denoising, a robust experimental paradigm, and individualized cortical parcellation across 23 participants to provide evidence that remembering one's past experiences and imagining the thoughts and feelings of others share a common neural substrate. Evidence from activation and connectivity analyses indicate overlapping and distinct functional profiles of these widely studied episodic and social processes.

Dynamic T 2 * relaxometry of hyperpolarized [1-13 C]pyruvate MRI in the human brain and kidneys.

PURPOSE: This study aimed to quantify T 2 * $$ {T}_2^{\ast } $$ for hyperpolarized [1-13 C]pyruvate and metabolites in the healthy human brain and renal cell carcinoma (RCC) patients at 3 T. METHODS: Dynamic T 2 * $$ {T}_2^{\ast } $$ values were measured with a metabolite-specific multi-echo spiral sequence. The dynamic T 2 * $$ {T}_2^{\ast } $$ of [1-13 C]pyruvate, [1-13 C]lactate, and 13 C-bicarbonate was estimated in regions of interest in the whole brain, sinus vein, gray matter, and white matter in healthy volunteers, as well as in kidney tumors and the contralateral healthy kidneys in a separate group of RCC patients. T 2 * $$ {T}_2^{\ast } $$ was fit using a mono-exponential function; and metabolism was quantified using pyruvate-to-lactate conversion rate maps and lactate-to-pyruvate ratio maps, which were compared with and without an estimated T 2 * $$ {T}_2^{\ast } $$ correction. RESULTS: The T 2 * $$ {T}_2^{\ast } $$ of pyruvate was shown to vary during the acquisition, whereas the T 2 * $$ {T}_2^{\ast } $$ of lactate and bicarbonate were relatively constant through time and across the organs studied. The T 2 * $$ {T}_2^{\ast } $$ of lactate was similar in gray matter (29.75 ± 1.04 ms), white matter (32.89 ± 0.9 ms), healthy kidney (34.61 ± 4.07 ms), and kidney tumor (33.01 ± 2.31 ms); and the T 2 * $$ {T}_2^{\ast } $$ of bicarbonate was different between whole-brain (108.17 ± 14.05 ms) and healthy kidney (58.45 ± 6.63 ms). The T 2 * $$ {T}_2^{\ast } $$ of pyruvate had similar trends in both brain and RCC studies, reducing from 75.56 ± 2.23 ms to 22.24 ± 1.24 ms in the brain and reducing from 122.72 ± 9.86 ms to 57.38 ± 7.65 ms in the kidneys. CONCLUSION: Multi-echo dynamic imaging can quantify T 2 * $$ {T}_2^{\ast } $$ and metabolism in a single integrated acquisition. Clear differences were observed in the T 2 * $$ {T}_2^{\ast } $$ of metabolites and in their behavior throughout the timecourse.

Fast and High-Resolution T2 Mapping Based on Echo Merging Plus k-t Undersampling with Reduced Refocusing Flip Angles (TEMPURA) as Methods for Human Renal MRI.

PURPOSE: To develop a highly accelerated multi-echo spin-echo method, TEMPURA, for reducing the acquisition time and/or increasing spatial resolution for kidney T2 mapping. METHODS: TEMPURA merges several adjacent echoes into one k-space by either combining independent echoes or sharing one echo between k-spaces. The combined k-space is reconstructed based on compressed sensing theory. Reduced flip angles are used for the refocusing pulses, and the extended phase graph algorithm is used to correct the effects of indirect echoes. Two sequences were developed: a fast breath-hold sequence; and a high-resolution sequence. The performance was evaluated prospectively on a phantom, 16 healthy subjects, and two patients with different types of renal tumors. RESULTS: The fast TEMPURA method reduced the acquisition time from 3-5 min to one breath-hold (18 s). Phantom measurements showed that fast TEMPURA had a mean absolute percentage error (MAPE) of 8.2%, which was comparable to a standardized respiratory-triggered sequence (7.4%), but much lower than a sequence accelerated by purely k-t undersampling (21.8%). High-resolution TEMPURA reduced the in-plane voxel size from 3 × 3 to 1 × 1 mm2, resulting in improved visualization of the detailed anatomical structure. In vivo T2 measurements demonstrated good agreement (fast: MAPE = 1.3%-2.5%; high-resolution: MAPE = 2.8%-3.3%) and high correlation coefficients (fast: R = 0.85-0.98; high-resolution: 0.82-0.96) with the standardized method, outperforming k-t undersampling alone (MAPE = 3.3-4.5%, R = 0.57-0.59). CONCLUSION: TEMPURA provides fast and high-resolution renal T2 measurements. It has the potential to improve clinical throughput and delineate intratumoral heterogeneity and tissue habitats at unprecedented spatial resolution.

DeepEMC-T2 mapping: Deep learning-enabled T2 mapping based on echo modulation curve modeling.

PURPOSE: Echo modulation curve (EMC) modeling enables accurate quantification of T2 relaxation times in multi-echo spin-echo (MESE) imaging. The standard EMC-T2 mapping framework, however, requires sufficient echoes and cumbersome pixel-wise dictionary-matching steps. This work proposes a deep learning version of EMC-T2 mapping, called DeepEMC-T2 mapping, to efficiently estimate accurate T2 maps from fewer echoes. METHODS: DeepEMC-T2 mapping was developed using a modified U-Net to estimate both T2 and proton density (PD) maps directly from MESE images. The network implements several new features to improve the accuracy of T2/PD estimation. A total of 67 MESE datasets acquired in axial orientation were used for network training and evaluation. An additional 57 datasets acquired in coronal orientation with different scan parameters were used to evaluate the generalizability of the framework. The performance of DeepEMC-T2 mapping was evaluated in seven experiments. RESULTS: Compared to the reference, DeepEMC-T2 mapping achieved T2 estimation errors from 1% to 11% and PD estimation errors from 0.4% to 1.5% with ten/seven/five/three echoes, which are more accurate than standard EMC-T2 mapping. By incorporating datasets acquired with different scan parameters and orientations for joint training, DeepEMC-T2 exhibits robust generalizability across varying imaging protocols. Increasing the echo spacing and including longer echoes improve the accuracy of parameter estimation. The new features proposed in DeepEMC-T2 mapping all enabled more accurate T2 estimation. CONCLUSIONS: DeepEMC-T2 mapping enables simplified, efficient, and accurate T2 quantification directly from MESE images without dictionary matching. Accurate T2 estimation from fewer echoes allows for increased volumetric coverage and/or higher slice resolution without prolonging total scan times.

QRAGE-Simultaneous multiparametric quantitative MRI of water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength.

PURPOSE: To introduce quantitative rapid gradient-echo (QRAGE), a novel approach for the simultaneous mapping of multiple quantitative MRI parameters, including water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength. METHODS: QRAGE leverages a newly developed multi-echo MPnRAGE sequence, facilitating the acquisition of 171 distinct contrast images across a range of inversion and TE points. To maintain a short acquisition time, we introduce MIRAGE2, a novel model-based reconstruction method that exploits prior knowledge of temporal signal evolution, represented as damped complex exponentials. MIRAGE2 minimizes local Block-Hankel and Casorati matrices. Parameter maps are derived from the reconstructed contrast images through postprocessing steps. We validate QRAGE through extensive simulations, phantom studies, and in vivo experiments, demonstrating its capability for high-precision imaging. RESULTS: In vivo brain measurements show the promising performance of QRAGE, with test-retest SDs and deviations from reference methods of < 0.8% for water content, < 17 ms for T1, and < 0.7 ms for T2*. QRAGE achieves whole-brain coverage at a 1-mm isotropic resolution in just 7 min and 15 s, comparable to the acquisition time of an MP2RAGE scan. In addition, QRAGE generates a contrast image akin to the UNI image produced by MP2RAGE. CONCLUSION: QRAGE is a new, successful approach for simultaneously mapping multiple MR parameters at ultrahigh field.

JUST-Net: Jointly unrolled cross-domain optimization based spatio-temporal reconstruction network for accelerated 3D myelin water imaging.

PURPOSE: We introduced a novel reconstruction network, jointly unrolled cross-domain optimization-based spatio-temporal reconstruction network (JUST-Net), aimed at accelerating 3D multi-echo gradient-echo (mGRE) data acquisition and improving the quality of resulting myelin water imaging (MWI) maps. METHOD: An unrolled cross-domain spatio-temporal reconstruction network was designed. The main idea is to combine frequency and spatio-temporal image feature representations and to sequentially implement convolution layers in both domains. The k-space subnetwork utilizes shared information from adjacent frames, whereas the image subnetwork applies separate convolutions in both spatial and temporal dimensions. The proposed reconstruction network was evaluated for both retrospectively and prospectively accelerated acquisition. Furthermore, it was assessed in simulation studies and real-world cases with k-space corruptions to evaluate its potential for motion artifact reduction. RESULTS: The proposed JUST-Net enabled highly reproducible and accelerated 3D mGRE acquisition for whole-brain MWI, reducing the acquisition time from fully sampled 15:23 to 2:22 min within a 3-min reconstruction time. The normalized root mean squared error of the reconstructed mGRE images increased by less than 4.0%, and the correlation coefficients for MWI showed a value of over 0.68 when compared to the fully sampled reference. Additionally, the proposed method demonstrated a mitigating effect on both simulated and clinical motion-corrupted cases. CONCLUSION: The proposed JUST-Net has demonstrated the capability to achieve high acceleration factors for 3D mGRE-based MWI, which is expected to facilitate widespread clinical applications of MWI.

Multiband accelerated 2D EPI for multi-echo brain QSM at 3 T.

PURPOSE: Data for QSM are typically acquired using multi-echo 3D gradient echo (GRE), but EPI can be used to accelerate QSM and provide shorter acquisition times. So far, EPI-QSM has been limited to single-echo acquisitions, which, for 3D GRE, are known to be less accurate than multi-echo sequences. Therefore, we compared single-echo and multi-echo EPI-QSM reconstructions across a range of parallel imaging and multiband acceleration factors. METHODS: Using 2D single-shot EPI in the brain, we compared QSM from single-echo and multi-echo acquisitions across combined parallel-imaging and multiband acceleration factors ranging from 2 to 16, with volume pulse TRs from 21.7 to 3.2 s, respectively. For single-echo versus multi-echo reconstructions, we investigated the effect of acceleration factors on regional susceptibility values, temporal noise, and image quality. We introduce a novel masking method based on thresholding the magnitude of the local field gradients to improve brain masking in challenging regions. RESULTS: At 1.6-mm isotropic resolution, high-quality QSM was achieved using multi-echo 2D EPI with a combined acceleration factor of 16 and a TR of 3.2 s, which enables functional applications. With these high acceleration factors, single-echo reconstructions are inaccurate and artefacted, rendering them unusable. Multi-echo acquisitions greatly improve QSM quality, particularly at higher acceleration factors, provide more consistent regional susceptibility values across acceleration factors, and decrease temporal noise compared with single-echo QSM reconstructions. CONCLUSION: Multi-echo acquisition is more robust for EPI-QSM across parallel imaging and multiband acceleration factors than single-echo acquisition. Multi-echo EPI can be used for highly accelerated acquisition while preserving QSM accuracy and quality relative to gold-standard 3D-GRE QSM.

Accelerated 3D multi-echo spin-echo sequence with a subspace constrained reconstruction for whole mouse brain T 2 mapping.

PURPOSE: To accelerate whole-brain quantitative T 2 $$ {\mathrm{T}}_2 $$ mapping in preclinical imaging setting. METHODS: A three-dimensional (3D) multi-echo spin echo sequence was highly undersampled with a variable density Poisson distribution to reduce the acquisition time. Advanced iterative reconstruction based on linear subspace constraints was employed to recover high-quality raw images. Different subspaces, generated using exponential or extended-phase graph (EPG) simulations or from low-resolution calibration images, were compared. The subspace dimension was investigated in terms of T 2 $$ {\mathrm{T}}_2 $$ precision. The method was validated on a phantom containing a wide range of T 2 $$ {\mathrm{T}}_2 $$ and was then applied to monitor metastasis growth in the mouse brain at 4.7T. Image quality and T 2 $$ {\mathrm{T}}_2 $$ estimation were assessed for 3 acceleration factors (6/8/10). RESULTS: The EPG-based dictionary gave robust estimations of a large range of T 2 $$ {\mathrm{T}}_2 $$ . A subspace dimension of 6 was the best compromise between T 2 $$ {\mathrm{T}}_2 $$ precision and image quality. Combining the subspace constrained reconstruction with a highly undersampled dataset enabled the acquisition of whole-brain T 2 $$ {\mathrm{T}}_2 $$ maps, the detection and the monitoring of metastasis growth of less than 500 µ m 3 $$ \mu {\mathrm{m}}^3 $$ . CONCLUSION: Subspace-based reconstruction is suitable for 3D T 2 $$ {\mathrm{T}}_2 $$ mapping. This method can be used to reach an acceleration factor up to 8, corresponding to an acquisition time of 25 min for an isotropic 3D acquisition of 156 µ $$ \mu $$ m on the mouse brain, used here for monitoring metastases growth.

Quantification of Cardiac Iron Overload at 3 T MRI in a Rabbit Model Utilizing ME-GRE T2* Sequence.

BACKGROUND: Myocardial iron overload can lead to myocardial dysfunction, muscle cell injury, and end-stage heart failure. The enhanced signal-to-noise ratio and technical advancements have made 3 T magnetic resonance imaging (MRI) more accessible in clinical settings. However, 3 T assessments for early diagnosis of myocardial iron overload are scarce. PURPOSE: To evaluate the feasibility of myocardial iron quantification using 3 T MRI in a rabbit model of iron overload. STUDY TYPE: Animal model. ANIMAL MODEL: Overall, 40 male New Zealand white rabbits were categorized into control (N = 8; no treatment) and experimental (N = 32; weekly 200 mg/kg iron dextran injections) groups. SEQUENCE: 3 T MRI with multi-echo gradient echo (ME-GRE) T2* sequence. ASSESSMENT: Each week, two experimental rabbits were randomly selected for blood collection to determine serum iron (SI) levels; their tissue was harvested to assess myocardial and hepatic iron deposition. STATISTICAL TESTS: Spearman's rank correlation tests were used to evaluate the correlations among R2*, cardiac iron concentration (CIC), liver iron concentration (LIC), total amount of iron injected, and SI levels. P ≤ 0.05 was considered statistically significant. RESULTS: The myocardial T2* value in the experimental group was significantly lower than that of the control group. An excellent correlation was observed between R2* values and CIC (r = 0.854). CIC moderately correlated with LIC (r = 0.712) and the total amount of iron injected (r = 0.698). A strong correlation was observed between the total amount of iron injected and LIC (r = 0.866). SI levels poorly correlated with the total amount of iron injected (r = 0.205, P = 0.277) and LIC (r = 0.170, P = 0.370) but fairly correlated with CIC (r = 0.415, P = 0.022). DATA CONCLUSION: A 3 T MRI with an ME-GRE sequence may serve as a noninvasive method for evaluating cardiac iron content. EVIDENCE LEVEL: N/A TECHNICAL EFFICACY: Stage 1.

Comparing the efficacy of data-driven denoising methods for a multi-echo fMRI acquisition at 7T.

In functional magnetic resonance imaging (fMRI) of the brain the measured signal is corrupted by several (e.g. physiological, motion, and thermal) noise sources and depends on the image acquisition. Imaging at ultrahigh field strength is becoming increasingly popular as it offers increased spatial accuracy. The latter is of particular benefit in brainstem neuroimaging given the small cross-sectional area of most nuclei. However, physiological noise scales with field strength in fMRI acquisitions. Although this problem is in part solved by decreasing voxel size, it is clear that adequate physiological denoising is of utmost importance in brainstem-focused fMRI experiments. Multi-echo sequences have been reported to facilitate highly effective denoising through TE-dependence of Blood Oxygen Level Dependent (BOLD) signals, in a denoising method referred to as multi-echo independent component analysis (ME-ICA). It has not been explored previously how ME-ICA compares to other data-driven denoising approaches at ultrahigh field strength. In the current study, we compared the efficacy of several denoising methods, including anatomical component based correction (aCompCor), Automatic Removal of Motion Artifacts (ICA-AROMA) aggressive and non-aggressive options, ME-ICA, and a combination of ME-ICA and aCompCor. We assessed several data quality metrics, including temporal signal-to-noise ratio (tSNR), delta variation signal (DVARS), spectral density of the global signal, functional connectivity and Shannon spectral entropy. Moreover, we looked at the ability of each method to uncouple the global signal and respiration. In line with previous reports at lower field strengths, we demonstrate that after applying ME-ICA, the data is best post-processed in order to remove spatially diffuse noise with a method such as aCompCor. Our findings indicate that ME-ICA combined with aCompCor and the aggressive option of ICA-AROMA are highly effective denoising approaches for multi-echo data acquired at 7T. ME-ICA combined with aCompCor potentially preserves more signal-of-interest as compared to the aggressive option of ICA-AROMA.

LARO: Learned acquisition and reconstruction optimization to accelerate quantitative susceptibility mapping.

Quantitative susceptibility mapping (QSM) involves acquisition and reconstruction of a series of images at multi-echo time points to estimate tissue field, which prolongs scan time and requires specific reconstruction technique. In this paper, we present our new framework, called Learned Acquisition and Reconstruction Optimization (LARO), which aims to accelerate the multi-echo gradient echo (mGRE) pulse sequence for QSM. Our approach involves optimizing a Cartesian multi-echo k-space sampling pattern with a deep reconstruction network. Next, this optimized sampling pattern was implemented in an mGRE sequence using Cartesian fan-beam k-space segmenting and ordering for prospective scans. Furthermore, we propose to insert a recurrent temporal feature fusion module into the reconstruction network to capture signal redundancies along echo time. Our ablation studies show that both the optimized sampling pattern and proposed reconstruction strategy help improve the quality of the multi-echo image reconstructions. Generalization experiments show that LARO is robust on the test data with new pathologies and different sequence parameters. Our code is available at https://github.com/Jinwei1209/LARO-QSM.git.

Multi-echo dipole inversion for magnetic susceptibility mapping.

PURPOSE: Reconstructing tissue magnetic susceptibility (QSM) from MRI phase data involves solving multiple consecutive ill-posed inverse problems such as phase unwrapping, background field removal, and field-to-source inversion. Multi-echo acquisitions present an additional challenge, as the magnetization field is typically computed from the multiple phase data prior to reconstructing the susceptibility map. Processing the multiple phase data introduces errors during the field estimation, violating assumptions of the subsequent inverse problems, manifesting as streaking artifacts in the susceptibility map. To address this challenge, we propose a multi-echo field-to-source forward model that forgoes the field estimation step. Moreover, we propose a fully general underestimation correction step to recover susceptibility sources that were regularized away during the field-to-source inversion. METHODS: The multi-echo forward model and correction step were validated on the QSM Challenge 2.0 datasets and compared to the standard single field-to-source model in in vivo human brains using different types of deconvolution algorithms. RESULTS: On the QSM Challenge 2.0 datasets the multi-echo forward model and correction step attain state-of-the-art results on all metrics by a wide margin. Experiments in in vivo brains show that the multi-echo model is in agreement with the single field-to-source model and that the proposed forward model and correction step can be used with any available dipole inversion method. CONCLUSION: A multi-echo field-to-source forward model forgoes the need to fit multi-echo phase data and achieves state-of-the-art results on the QSM Challenge 2.0 data. Underestimated low-frequency susceptibility distributions can be partially recovered using a correction step.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

PURPOSE: To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction. METHODS: Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF), R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study. RESULTS: Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05, R 2 $$ {R}^2 $$ = 0.999) and Cartesian without motion ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.32 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04, R 2 $$ {R}^2 $$ = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524, R 2 $$ {R}^2 $$ = 0.977) compared to motion-averaged reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572, R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM. CONCLUSION: Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Lung parenchyma transverse relaxation rates at 0.55 T.

PURPOSE: To determine R2 and R 2 ' $$ {R}_2^{\prime } $$ transverse relaxation rates in healthy lung parenchyma at 0.55 T. This is important in that it informs the design and optimization of new imaging methods for 0.55T lung MRI. METHODS: Experiments were performed in 3 healthy adult volunteers on a prototype whole-body 0.55T MRI, using a custom free-breathing electrocardiogram-triggered, single-slice echo-shifted multi-echo spin echo (ES-MCSE) pulse sequence with respiratory navigation. Transverse relaxation rates R2 and R 2 ' $$ {R}_2^{\prime } $$ and off-resonance ∆f were jointly estimated using nonlinear least-squares estimation. These measurements were compared against R2 estimates from T2 -prepared balanced SSFP (T2 -Prep bSSFP) and R 2 * $$ {R}_2^{\ast } $$ estimates from multi-echo gradient echo, which are used widely but prone to error due to different subvoxel weighting. RESULTS: The mean R2 and R 2 ' $$ {R}_2^{\prime } $$ values of lung parenchyma obtained from ES-MCSE were 17.3 ± 0.7 Hz and 127.5 ± 16.4 Hz (T2  = 61.6 ± 1.7 ms; T 2 ' $$ {\mathrm{T}}_2^{\prime } $$  = 9.5 ms ± 1.6 ms), respectively. The off-resonance estimates ranged from -60 to 30 Hz. The R2 from T2 -Prep bSSFP was 15.7 ± 1.7 Hz (T2  = 68.6 ± 8.6 ms) and R 2 * $$ {R}_2^{\ast } $$ from multi-echo gradient echo was 131.2 ± 30.4 Hz ( T 2 * $$ {\mathrm{T}}_2^{\ast } $$  = 8.0 ± 2.5 ms). Paired t-test indicated that there is a significant difference between the proposed and reference methods (p < 0.05). The mean R2 estimate from T2 -Prep bSSFP was slightly smaller than that from ES-MCSE, whereas the mean R 2 ' $$ {R}_2^{\prime } $$ and R 2 * $$ {R}_2^{\ast } $$ estimates from ES-MCSE and multi-echo gradient echo were similar to each other across all subjects. CONCLUSIONS: Joint estimation of transverse relaxation rates and off-resonance is feasible at 0.55 T with a free-breathing electrocardiogram-gated and navigator-gated ES-MCSE sequence. At 0.55 T, the mean R2 of 17.3 Hz is similar to the reported mean R2 of 16.7 Hz at 1.5 T, but the mean R 2 ' $$ {R}_2^{\prime } $$ of 127.5 Hz is about 5-10 times smaller than that reported at 1.5 T.

Subtractionless compressed-sensing-accelerated whole-body MR angiography using two-point Dixon fat suppression with single-pass half-reduced contrast dose: feasibility study and initial experience.

PURPOSE: To investigate the feasibility and clinical utility of a compressed-sensing-accelerated subtractionless whole-body MRA (CS-WBMRA) protocol with only contrast injection for suspected arterial diseases, by comparison to conventional dual-pass subtraction-based whole-body MRA (conventional-WBMRA) and available computed tomography angiography (CTA). MATERIALS AND METHODS: This prospective study assessed 86 patients (mean age, 56 years ± 16.4 [standard deviation]; 25 women) with suspected arterial diseases from May 2021 to December 2022, who underwent CS-WBMRA (n = 48, mean age, 55.9 years ± 16.4 [standard deviation]; 25 women) and conventional-WBMRA (n = 38, mean age, 48 years ± 17.4 [standard deviation]; 20 women) on a 3.0 T MRI after random group assignment based on the chronological order of enrolment. Of all enrolled patients administered the CS-WBMRA protocol, 35% (17/48) underwent CTA as required by clinical demands. Two experienced radiologists independently scored the qualitative image quality and venous enhancement contamination. Quantitative image assessment was carried out by determining and comparing the apparent signal-to-noise ratios (SNRs) and contrast-to-noise ratios (CNRs) of four representative arterial segments. The total examination time and contrast-dose were also recorded. The independent samples t-test or the Wilcoxon rank sum test was used for statistical analysis. RESULTS: The overall scores of CS-WBMRA outperformed those of conventional-WMBRA (3.40 ± 0.60 vs 3.22 ± 0.55, P < 0.001). In total, 1776 and 1406 arterial segments in the CS-WBMRA and conventional-WBMRA group were evaluated. Qualitative image scores for 7 (of 15) vessel segments in the CS-WMBRA group had statistically significantly increased values compared to those of the conventional-WBMRA groups (P < 0.05). Scores from the other 8 segments showed similar image quality (P > 0.05) between the two protocols. In the quantitative analysis, overall apparent SNRs were significantly higher in the conventional-WBMRA group than in the CS-WBMRA group (214.98 ± 136.05 vs 164.90 ± 118.05; P < 0.001), while overall apparent CNRs were not significantly different in these two groups (CS vs conventional: 107.13 ± 72.323 vs 161.24 ± 118.64; P > 0.05). In the CS-WBMRA group, 7 of 1776 (0.4%) vessel segments were contaminated severely by venous enhancement, while in the convention-WBMRA group, 317 of 1406 (23%) were rated as severe contamination. In the CS-WBMRA group, total examination and reconstruction times were only 7 min and 10 min, respectively, vs 20 min and < 30 s for the conventional WBMRA group, respectively. The contrast agent dose used in the CS-WBMRA protocol was reduced by half compared to conventional-WBMRA protocol (18.7 ± 3.5 ml vs 37.2 ± 5.4 ml, P = 0.008). CONCLUSION: The CS-WBMRA protocol provides excellent image quality and sufficient diagnostic accuracy for whole-body arterial disease, with relatively faster workflow and half-dose reduction of contrast agent, which has greater potential in clinical practice compared with conventional-WBMRA.

Age differences in the functional architecture of the human brain.

The intrinsic functional organization of the brain changes into older adulthood. Age differences are observed at multiple spatial scales, from global reductions in modularity and segregation of distributed brain systems, to network-specific patterns of dedifferentiation. Whether dedifferentiation reflects an inevitable, global shift in brain function with age, circumscribed, experience-dependent changes, or both, is uncertain. We employed a multimethod strategy to interrogate dedifferentiation at multiple spatial scales. Multi-echo (ME) resting-state fMRI was collected in younger (n = 181) and older (n = 120) healthy adults. Cortical parcellation sensitive to individual variation was implemented for precision functional mapping of each participant while preserving group-level parcel and network labels. ME-fMRI processing and gradient mapping identified global and macroscale network differences. Multivariate functional connectivity methods tested for microscale, edge-level differences. Older adults had lower BOLD signal dimensionality, consistent with global network dedifferentiation. Gradients were largely age-invariant. Edge-level analyses revealed discrete, network-specific dedifferentiation patterns in older adults. Visual and somatosensory regions were more integrated within the functional connectome; default and frontoparietal control network regions showed greater connectivity; and the dorsal attention network was more integrated with heteromodal regions. These findings highlight the importance of multiscale, multimethod approaches to characterize the architecture of functional brain aging.

Improved Resting-State Functional MRI Using Multi-Echo Echo-Planar Imaging on a Compact 3T MRI Scanner with High-Performance Gradients.

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.

Optimization of Gradient-Echo Echo-Planar Imaging for T2* Contrast in the Brain at 0.5 T.

Gradient-recalled echo (GRE) echo-planar imaging (EPI) is an efficient MRI pulse sequence that is commonly used for several enticing applications, including functional MRI (fMRI), susceptibility-weighted imaging (SWI), and proton resonance frequency (PRF) thermometry. These applications are typically not performed in the mid-field (<1 T) as longer T2* and lower polarization present significant challenges. However, recent developments of mid-field scanners equipped with high-performance gradient sets offer the possibility to re-evaluate the feasibility of these applications. The paper introduces a metric "T2* contrast efficiency" for this evaluation, which minimizes dead time in the EPI sequence while maximizing T2* contrast so that the temporal and pseudo signal-to-noise ratios (SNRs) can be attained, which could be used to quantify experimental parameters for future fMRI experiments in the mid-field. To guide the optimization, T2* measurements of the cortical gray matter are conducted, focusing on specific regions of interest (ROIs). Temporal and pseudo SNR are calculated with the measured time-series EPI data to observe the echo times at which the maximum T2* contrast efficiency is achieved. T2* for a specific cortical ROI is reported at 0.5 T. The results suggest the optimized echo time for the EPI protocols is shorter than the effective T2* of that region. The effective reduction of dead time prior to the echo train is feasible with an optimized EPI protocol, which will increase the overall scan efficiency for several EPI-based applications at 0.5 T.

A geometric view of signal sensitivity metrics in multi-echo fMRI.

In multi-echo fMRI (ME-fMRI), two metrics have been widely used to measure the performance of various acquisition and analysis approaches. These are temporal SNR (tSNR) and differential contrast-to-noise ratio (dCNR). A key step in ME-fMRI is the weighted combination of the data from multiple echoes, and prior work has examined the dependence of tSNR and dCNR on the choice of weights. However, most studies have focused on only one of these two metrics, and the relationship between the two metrics has not been examined. In this work, we present a geometric view that offers greater insight into the relation between the two metrics and their weight dependence. We identify three major regimes: (1) a tSNR robust regime in which tSNR is robust to the weight selection with most weight variants achieving close to optimal performance, whereas dCNR shows a pronounced dependence on the weights with most variants achieving suboptimal performance; (2) a dCNR robust regime in which dCNR is robust to the weight selection with most weight variants achieving close to optimal performance, while tSNR exhibits a strong dependence on the weights with most variants achieving significantly lower than optimal performance; and (3) a within-type robust regime in which both tSNR and dCNR achieve nearly optimal performance when the form of the weights are variants of their respective optimal weights and exhibit a moderate decrease in performance for other weight variants. Insight into the behavior observed in the different regimes is gained by considering spherical representations of the weight dependence of the components used to form each metric. For multi-echo acquisitions, dCNR is shown to be more directly related than tSNR to measures of CNR and signal-to-noise ratio (SNR) for task-based and resting-state fMRI scans, respectively.

Evaluating the efficacy of multi-echo ICA denoising on model-based fMRI.

fMRI is an indispensable tool for neuroscience investigation, but this technique is limited by multiple sources of physiological and measurement noise. These noise sources are particularly problematic for analysis techniques that require high signal-to-noise ratio for stable model fitting, such as voxel-wise modeling. Multi-echo data acquisition in combination with echo-time dependent ICA denoising (ME-ICA) represents one promising strategy to mitigate physiological and hardware-related noise sources as well as motion-related artifacts. However, most studies employing ME-ICA to date are resting-state fMRI studies, and therefore we have a limited understanding of the impact of ME-ICA on complex task or model-based fMRI paradigms. Here, we addressed this knowledge gap by comparing data quality and model fitting performance of data acquired during a visual population receptive field (pRF) mapping (N = 13 participants) experiment after applying one of three preprocessing procedures: ME-ICA, optimally combined multi-echo data without ICA-denoising, and typical single echo processing. As expected, multi-echo fMRI improved temporal signal-to-noise compared to single echo fMRI, with ME-ICA amplifying the improvement compared to optimal combination alone. However, unexpectedly, this boost in temporal signal-to-noise did not directly translate to improved model fitting performance: compared to single echo acquisition, model fitting was only improved after ICA-denoising. Specifically, compared to single echo acquisition, ME-ICA resulted in improved variance explained by our pRF model throughout the visual system, including anterior regions of the temporal and parietal lobes where SNR is typically low, while optimal combination without ICA did not. ME-ICA also improved reliability of parameter estimates compared to single echo and optimally combined multi-echo data without ICA-denoising. Collectively, these results suggest that ME-ICA is effective for denoising task-based fMRI data for modeling analyzes and maintains the integrity of the original data. Therefore, ME-ICA may be beneficial for complex fMRI experiments, including voxel-wise modeling and naturalistic paradigms.