Perfusion weighted imaging using combined gradient/spin echo EPIK: Brain tumour applications in hybrid MR-PET

Rapid simultaneous estimation of relaxation rates using multi-echo, multi-contrast MRI

PURPOSE

Multi-echo, multi-contrast methods are increasingly used in dynamic imaging studies to simultaneously quantify R2∗ and R2. To overcome the computational challenges associated with nonlinear least squares (NLSQ) fitting, we propose a generalized linear least squares (LLSQ) solution to rapidly fit R2∗ and R2.

METHODS

Spin- and gradient-echo (SAGE) data were simulated across T2∗ and T2 values at high (200) and low (20) SNR. Full (four-parameter) and reduced (three-parameter) parameter fits were implemented and compared with both LLSQ and NLSQ fitting. Fit data were compared to ground truth using concordance correlation coefficient (CCC) and coefficient of variation (CV). In vivo SAGE perfusion data were acquired in 20 subjects with relapsing-remitting multiple sclerosis. LLSQ R2∗ and R2, as well as cerebral blood volume (CBV), were compared with the standard NLSQ approach.

RESULTS

Across all fitting methods, T2∗ was well-fit at high (CCC = 1, CV = 0) and low (CCC ≥ 0.87, CV ≤ 0.08) SNR. Except for short T2∗ values (5-15 ms), T2 was well-fit at high (CCC = 1, CV = 0) and low (CCC ≥ 0.99, CV ≤ 0.03) SNR. In vivo, LLSQ R2∗ and R2 estimates were similar to NLSQ, and there were no differences in R2∗ across fitting methods at high SNR. However, there were some differences at low SNR and for R2 at high and low SNR. In vivo NLSQ and LLSQ three parameter fits performed similarly, as did NLSQ and LLSQ four-parameter fits. LLSQ CBV nearly matched the standard NLSQ method for R2∗- (0.97 ratio) and R2-CBV (0.98 ratio). Voxel-wise whole-brain fitting was faster for LLSQ (3-4 min) than NLSQ (16-18 h).

CONCLUSIONS

LLSQ reliably fit for R2∗ and R2 in simulated and in vivo data. Use of LLSQ methods reduced the computational demand, enabling rapid estimation of R2∗ and R2.

Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

Development of a novel 10-echo multi-contrast sequence based on EPIK to deliver simultaneous quantification of T2 and T2 * with application to oxygen extraction fraction

PURPOSE

The simultaneous quantification of T₂ and T₂ * maps based on fast sequences for combined GE and SE acquisition has rekindled research and clinical interest by offering a wide range of attractive applications, e.g., dynamic tracking of oxygen extraction fraction (OEF). However, previously published methods based on EPI-readouts have been hindered by resolution and the number of acquired echoes.

METHODS

This work presents a novel 10-echo GE-SE EPIK (EPI with keyhole) sequence for the rapid quantification of T₂ '. T₂ /T₂ * maps from the GE-SE EPIK sequence were validated using three phantoms and 15 volunteers at 3T. The incorporation of a sliding window approach, combined with the full sampling of the k-space center inherent to EPIK, enables a high effective temporal resolution. That is, for an eight-slice breath-hold experiment, a temporal sampling rate of eight reconstructed slices per 1.1 s.

RESULTS

In comparison with repeated single-echo SE, multi-echo GE, and spectroscopy methods, the GE-SE EPIK sequence shows good agreement in quantifying T₂ /T₂ * values, while the gray matter/white matter separation yielded the expected contrast differentiation. The OEF was calculated with a view to an initial application with clinical relevance, producing results comparable to those in the literature and with good sensitivity in breath-hold experiments.

CONCLUSIONS

GE-SE EPIK provides increased resolution and more echoes, including two SEs, than comparable sequences. Moreover, GE-SE EPIK achieves this within an acquisition time of 57 s for 20 slices (matrix size = 128×128; FOV = 24 cm) and with a reasonably short TE for the final echo (114 ms). The sequence can dynamically track OEF changes in a breath-hold experiment.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Mapping of whole-cerebrum resting-state networks using ultra-high resolution acquisition protocols

Resting‐state functional magnetic resonance imaging (fMRI) has been used in numerous studies to map networks in the brain that employ spatially disparate regions. However, attempts to map networks with high spatial resolution have been hampered by conflicting technical demands and associated problems. Results from recent fMRI studies have shown that spatial resolution remains around 0.7 × 0.7 × 0.7 mm³, with only partial brain coverage. Therefore, this work aims to present a novel fMRI technique that was developed based on echo‐planar‐imaging with keyhole (EPIK) combined with repetition‐time‐external (TR‐external) EPI phase correction. Each technique has been previously shown to be effective in enhancing the spatial resolution of fMRI, and in this work, the combination of the two techniques into TR‐external EPIK provided a nominal spatial resolution of 0.51 × 0.51 × 1.00 mm³ (0.26 mm³ voxel) with whole‐cerebrum coverage. Here, the feasibility of using half‐millimetre in‐plane TR‐external EPIK for resting‐state fMRI was validated using 13 healthy subjects and the corresponding reproducible mapping of resting‐state networks was demonstrated. Furthermore, TR‐external EPIK enabled the identification of various resting‐state networks distributed throughout the brain from a single fMRI session, with mapping fidelity onto the grey matter at 7T. The high‐resolution functional image further revealed mesoscale anatomical structures, such as small cerebral vessels and the internal granular layer of the cortex within the postcentral gyrus.

Perfusion weighted imaging using combined gradient/spin echo EPIK: Brain tumour applications in hybrid MR-PET

Advanced perfusion-weighted imaging (PWI) methods that combine gradient echo (GE) and spin echo (SE) data are important tools for the study of brain tumours. In PWI, single-shot, EPI-based methods have been widely used due to their relatively high imaging speed. However, when used with increasing spatial resolution, single-shot EPI methods often show limitations in whole-brain coverage for multi-contrast applications. To overcome this limitation, this work employs a new version of EPI with keyhole (EPIK) to provide five echoes: two with GEs, two with mixed GESE and one with SE; the sequence is termed "GESE-EPIK." The performance of GESE-EPIK is evaluated against its nearest relative, EPI, in terms of the temporal signal-to-noise ratio (tSNR). Here, data from brain tumour patients were acquired using a hybrid 3T MR-BrainPET scanner. GESE-EPIK resulted in reduced susceptibility artefacts, shorter TEs for the five echoes and increased brain coverage when compared to EPI. Moreover, compared to EPI, EPIK achieved a comparable tSNR for the first and second echoes and significantly higher tSNR for other echoes. A new method to obtain multi-echo GE and SE data with shorter TEs and increased brain coverage is demonstrated. As proposed here, the workflow can be shortened and the integration of multimodal clinical MR-PET studies can be facilitated.