Spatiotemporal characterization of breathing-induced B0 field fluctuations in the cervical spinal cord at 7T

Dixon-based B0 self-navigation in radial stack-of-stars multi-echo gradient echo imaging

PURPOSE

To develop a Dixon-basedB 0 $$ {\mathrm{B}}_0 $$ self-navigation approach to estimate and correct temporalB 0 $$ {\mathrm{B}}_0 $$ variations in radial stack-of-stars gradient echo imaging for quantitative body MRI.

METHODS

The proposed method estimates temporalB 0 $$ {\mathrm{B}}_0 $$ variations using aB 0 $$ {\mathrm{B}}_0 $$ self-navigator estimated by a graph-cut-based water-fat separation algorithm on the oversampled k-space center. TheB 0 $$ {\mathrm{B}}_0 $$ self-navigator was employed to correct for phase differences between radial spokes (one-dimensional [1D] correction) and to perform a motion-resolved reconstruction to correct spatiotemporal pseudo-periodicB 0 $$ {\mathrm{B}}_0 $$ variations (three-dimensional [3D] correction). Numerical simulations, phantom experiments and in vivo neck scans were performed to evaluate the effects of temporalB 0 $$ {\mathrm{B}}_0 $$ variations on the field-map, proton density fat fraction (PDFF) andT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ map, and to validate the proposed method.

RESULTS

TemporalB 0 $$ {\mathrm{B}}_0 $$ variations were found to cause signal loss and phase shifts on the multi-echo images that lead to an underestimation ofT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ , while PDFF mapping was less affected. TheB 0 $$ {\mathrm{B}}_0 $$ self-navigator captured slowly varying temporalB 0 $$ {\mathrm{B}}_0 $$ drifts and temporal variations caused by respiratory motion. While the 1D correction effectively correctedB 0 $$ {\mathrm{B}}_0 $$ drifts in phantom studies, it was insufficient in vivo due to 3D spatially varying temporalB 0 $$ {\mathrm{B}}_0 $$ variations with amplitudes of up to 25 Hz at 3 T near the lungs. The proposed 3D correction locally improved the correction of field-map andT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ and reduced image artifacts.

CONCLUSION

TemporalB 0 $$ {\mathrm{B}}_0 $$ variations particularly affectT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ mapping in radial stack-of-stars imaging. The self-navigation approach can be applied without modifying the MR acquisition to correct forB 0 $$ {\mathrm{B}}_0 $$ drift and physiological motion-inducedB 0 $$ {\mathrm{B}}_0 $$ variations, especially in the presence of fat.

Thermal stimulus task fMRI in the cervical spinal cord at 7 Tesla

Although functional magnetic resonance imaging (fMRI) is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0 . Increasing field strength enables higher spatial resolution and improved sensitivity to blood oxygenation level-dependent (BOLD) signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D echo-planar imaging (EPI) protocols, which differ in sensitivity to spatial and temporal B0 inhomogeneity. The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. All protocols used 3 mm slice thickness. For each protocol, the BOLD response to 13 10-s noxious thermal stimuli applied to the right thumb was acquired in a 10-min fMRI run. Image quality, temporal signal to noise ratio (SNR), and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols. Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed. Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment. PRACTITIONER POINTS: First stimulus task fMRI results in the spinal cord at 7 T. Single-shot 0.75 mm 2D EPI produced the highest mean z-statistic. Multi-shot 0.60 mm 2D EPI provided the best-localized activation and least distortion.

Prospective motion correction for cervical spinal cord MRS

PURPOSE

To develop prospective motion correction for single-voxel MRS in the human cervical spinal cord.

METHODS

A motion MR navigator was implemented using reduced field-of-view 2D-selective RF excitation together with EPI readout. A short-echo semi-LASER sequence (TE  = 30 ms) was updated to incorporate this real-time image-based motion navigator, as well as real-time shim and frequency navigators. Five healthy participants were studied at 3 T with a 64-channel head-neck receive coil. Single-voxel MRS data were measured in a voxel located at the C3-5 vertebrae level. The motion navigator was used to correct for translations in the X-Y plane and was validated by assessing spectral quality with and without prospective correction in the presence of subject motion.

RESULTS

Without prospective correction, motion resulted in severe lipid contamination in the MR spectra. With prospective correction, the quality of spinal cord MR spectra in the presence of motion was similar to that obtained in the absence of motion, with comparable spectral signal-to-noise ratio and linewidth and no significant lipid contamination.

CONCLUSION

Prospective motion and B0 correction allow acquisition of good-quality MR spectra in the human cervical spinal cord in the presence of motion. This new technique should facilitate reliable acquisition of spinal cord MR spectra in both research and clinical settings.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Thermal Stimulus Task fMRI in the Cervical Spinal Cord at 7 Tesla

PURPOSE

Although functional MRI is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0. Increasing field strength enables higher spatial resolution and improved sensitivity to BOLD signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first stimulus task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D EPI protocols, as they differ in sensitivity to spatial and temporal B0 inhomogeneity.

METHODS

The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. For each protocol, the BOLD response to thirteen 10-second noxious thermal stimuli applied to the right thumb was acquired in a 10-minute fMRI run. Image quality, temporal SNR, and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols.

RESULTS

Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed.

CONCLUSION

Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Characterization of motion dependent magnetic field inhomogeneity for DWI in the kidneys

PURPOSE

Diffusion-weighted imaging (DWI) of the abdomen has increased dramatically for both research and clinical purposes. Motion and static field inhomogeneity related challenges limit image quality of abdominopelvic imaging with the most conventional echo-planar imaging (EPI) pulse sequence. While reversed phase encoded imaging is increasingly used to facilitate distortion correction, it typically assumes one motion independent magnetic field distribution. In this study, we describe a more generalized workflow for the case of kidney DWI in which the field inhomogeneity at multiple respiratory phases is mapped and used to correct all images in a multi-contrast DWI series.

METHODS

In this HIPAA-compliant and IRB-approved prospective study, 8 volunteers (6 M, ages 28-51) had abdominal imaging performed in a 3 T MRI system (MAGNETOM Prisma; Siemens Healthcare, Erlangen, Germany) with ECG gating. Coronal oblique T2-weighted HASTE images were collected for anatomical reference. Sagittal phase-contrast (PC) MRI images through the left renal artery were collected to determine systolic and diastolic phases. Cardiac triggered oblique coronal DWI were collected at 10 b-values between 0 and 800 s/mm2 and 12 directions. DWI series were distortion corrected using field maps generated by forward and reversed phase encoded b = 0 images collected over the full respiratory cycle and matched by respiratory phase. Morphologic accuracy, intraseries spatial variability, and diffusion tensor imaging (DTI) metrics mean diffusivity (MD) and fractional anisotropy (FA) were compared for results generated with no distortion correction, correction with only one respiratory bin, and correction with multiple respiratory bins across the breathing cycle.

RESULTS

Computed field maps showed significant variation in static field with kidney laterality, region, and respiratory phase. Distortion corrected images showed significantly better registration to morphologic images than uncorrected images; for the left kidney, the multiple bin correction outperformed one bin correction. Line profile analysis showed significantly reduced spatial variation with multiple bins than one bin correction. DTI metrics were mostly similar between correction methods, with some differences observed in MD between uncorrected and corrected datasets.

CONCLUSIONS

Our results indicate improved morphology of kidney DWI and derived parametric maps as well as reduced variability over the full image series using the motion-resolved distortion correction. This work highlights some morphologic and quantitative metric improvements can be obtained for kidney DWI when distortion correction is performed in a respiratory-resolved manner.

Real-time shimming with FID navigators

PURPOSE

To implement a method for real-time field control using rapid FID navigator (FIDnav) measurements and evaluate the efficacy of the proposed approach for mitigating dynamic field perturbations and improvingT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality.

METHODS

FIDnavs were embedded in a gradient echo sequence and a subject-specific linear calibration model was generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. To confirm the accuracy of FID-navigated field updates, phantom and volunteer scans were performed with online updates of the scanner B0 shim settings. To evaluate improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality with real-time shimming, 10 volunteers were scanned at 3T while performing deep-breathing and nose-touching tasks designed to modulate the B0 field. Quantitative image quality metrics were compared with and without FID-navigated field control. An additional volunteer was scanned at 7T to evaluate performance at ultra-high field.

RESULTS

Applying measured FIDnav shim updates successfully compensated for applied global and linear field offsets in phantoms and across all volunteers. FID-navigated real-time shimming led to a substantial reduction in field fluctuations and a consequent improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality in volunteers performing deep-breathing and nose-touching tasks, with 7.57% ± 6.01% and 8.21% ± 10.90% improvement in peak SNR and structural similarity, respectively.

CONCLUSION

FIDnavs facilitate rapid measurement and application of field coefficients for slice-wise B0 shimming. The proposed approach can successfully counteract spatiotemporal field perturbations and substantially improvesT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality, which is important for a variety of clinical and research applications, particularly at ultra-high field.

Impact of autocalibration method on accelerated EPI of the cervical spinal cord at 7 T

PURPOSE

The spinal cord contains sensorimotor neural circuits of scientific and clinical interest. However, spinal cord functional MRI (fMRI) is significantly more technically demanding than brain fMRI, due primarily to its proximity to the lungs. Accelerated echo-planar imaging (EPI) at 7 T is particularly vulnerable to k-space phase inconsistencies induced by motion or B0 fluctuation, during either autocalibration signal (ACS) or time-series acquisition. For 7 T brain fMRI, sensitivity to motion and B0 fluctuation can be reduced using a re-ordered segmented EPI ACS based on the fast low-angle excitation echo-planar technique (FLEET). However, respiration-induced B0 fluctuations (exceeding 100 Hz at C7) are greater, and fewer k-space lines per slice are required for cervical spinal cord fMRI at 7 T, necessitating a separate evaluation of ACS methods.

METHODS

We compared 24-line single-shot EPI with 48-line two-shot segmented EPI, two-shot FLEET, and gradient echo (GRE)-based ACS acquisition methods, performed under various physiological conditions, in terms of temporal signal-to-noise ratio and prevalence of artifacts in generalized autocalibrating partially parallel acquisition (GRAPPA)-accelerated EPI of the cervical spinal cord at 7 T.

RESULTS

Segmented EPI and FLEET ACS produce images with nearly identical patterns of severe image artifacts. GRE and single-shot EPI ACS consistently produce images free from significant artifacts, and temporal signal-to-noise ratio is significantly greater for GRE ACS, particularly in lower slices where through-slice dephasing is most severe.

CONCLUSIONS

GRE and single-shot EPI-ACS acquisition methods, which are robust to respiration-induced phase errors between k-space segments, produce images with fewer and less severe artifacts than either FLEET or conventionally segmented EPI for accelerated EPI of the cervical spinal cord at 7 T.

Optimized multi-echo gradient-echo magnetic resonance imaging for gray and white matter segmentation in the lumbosacral cord at 3 T

Atrophy in the spinal cord (SC), gray (GM) and white matter (WM) is typically measured in-vivo by image segmentation on multi-echo gradient-echo magnetic resonance images. The aim of this study was to establish an acquisition and analysis protocol for optimal SC and GM segmentation in the lumbosacral cord at 3 T. Ten healthy volunteers underwent imaging of the lumbosacral cord using a 3D spoiled multi-echo gradient-echo sequence (Siemens FLASH, with 5 echoes and 8 repetitions) on a Siemens Prisma 3 T scanner. Optimal numbers of successive echoes and signal averages were investigated comparing signal-to-noise (SNR) and contrast-to-noise ratio (CNR) values as well as qualitative ratings for segmentability by experts. The combination of 5 successive echoes yielded the highest CNR between WM and cerebrospinal fluid and the highest rating for SC segmentability. The combination of 3 and 4 successive echoes yielded the highest CNR between GM and WM and the highest rating for GM segmentability in the lumbosacral enlargement and conus medullaris, respectively. For segmenting the SC and GM in the same image, we suggest combining 3 successive echoes. For SC or GM segmentation only, we recommend combining 5 or 3 successive echoes, respectively. Six signal averages yielded good contrast for reliable SC and GM segmentation in all subjects. Clinical applications could benefit from these recommendations as they allow for accurate SC and GM segmentation in the lumbosacral cord.

Quantitative 23 Na-MRI of the intervertebral disk at 3 T

Monitoring the tissue sodium content (TSC) in the intervertebral disk geometry noninvasively by MRI is a sensitive measure to estimate changes in the proteoglycan content of the intervertebral disk, which is a biomarker of degenerative disk disease (DDD) and of lumbar back pain (LBP). However, application of quantitative sodium concentration measurements in ²³ Na-MRI is highly challenging due to the lower in vivo concentrations and smaller gyromagnetic ratio, ultimately yielding much smaller signal relative to ¹ H-MRI. Moreover, imaging the intervertebral disk geometry imposes higher demands, mainly because the necessary RF volume coils produce highly inhomogeneous transmit field patterns. For an accurate absolute quantification of TSC in the intervertebral disks, the B1 field variations have to be mitigated. In this study, we report for the first time quantitative sodium concentration in the intervertebral disks at clinical field strengths (3 T) by deploying ²³ Na-MRI in healthy human subjects. The sodium B1 maps were calculated by using the double-angle method and a double-tuned (¹ H/²³ Na) transceive chest coil, and the individual effects of the variation in the B1 field patterns in tissue sodium quantification were calculated. Phantom measurements were conducted to evaluate the quality of the Na-weighted images and B1 mapping. Depending on the disk position, the sodium concentration was calculated as 161.6 mmol/L-347 mmol/L, and the mean sodium concentration of the intervertebral disks varies between 254.6 ± 54 mmol/L and 290.1 ± 39 mmol/L. A smoothing effect of the B1 correction on the sodium concentration maps was observed, such that the standard deviation of the mean sodium concentration was significantly reduced with B1 mitigation. The results of this work provide an improved integration of quantitative ²³ Na-MRI into clinical studies in intervertebral disks such as degenerative disk disease and establish alternative scoring schemes to existing morphological scoring such as the Pfirrmann score.

Respiratory-triggered quantitative MR spectroscopy of the human cervical spinal cord at 7 T

PURPOSE

Ultra-high field 1 H MR spectroscopy (MRS) is of great interest to help characterizing human spinal cord pathologies. However, very few studies have been reported so far in this small size structure at these fields due to challenging experimental difficulties caused by static and radiofrequency field heterogeneities, as well as physiological motion. In this work, in line with the recent developments proposed to strengthen spinal cord MRS feasibility at 7 T, a respiratory-triggered acquisition approach was optimized to compensate for dynamic B 0 field heterogeneities and to provide robust cervical spinal cord MRS data.

METHODS

A semi-LASER sequence was purposely used, and a dedicated raw data processing algorithm was developed to enhance MR spectral quality by discarding corrupted scans. To legitimate the choices done during the optimization stage, additional tests were carried out to determine the impact of breathing, voluntary motion, body mass index, and fitting algorithm. An in-house quantification tool was concomitantly designed for accurate estimation of the metabolite concentration ratios for choline, N-acetyl-aspartate (NAA), myo-inositol and glutathione. The method was tested on a cohort of 14 healthy volunteers.

RESULTS

Average water linewidth and NAA signal-to-noise ratio reached 0.04 ppm and 11.01, respectively. The group-average metabolic ratios were in good agreement with previous studies and showed intersession reproducibility variations below 30%.

CONCLUSION

The developed approach allows a rise of the acquired MRS signal quality and of the quantification robustness as compared to previous studies hence offering strengthened possibilities to probe the metabolism of degenerative and traumatic spinal cord pathologies.

Multi-shot acquisitions for stimulus-evoked spinal cord BOLD fMRI

PURPOSE

To demonstrate the feasibility of 3D multi-shot magnetic resonance imaging acquisitions for stimulus-evoked blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) in the human spinal cord in vivo.

METHODS

Two fMRI studies were performed at 3T. The first study was a hypercapnic gas challenge where data were acquired from healthy volunteers using a multi-shot 3D fast field echo (FFE) sequence as well as single-shot multi-slice echo-planar imaging (EPI). In the second study, another cohort of healthy volunteers performed an upper extremity motor task while fMRI data were acquired using a 3D multi-shot acquisition.

RESULTS

Both 2D-EPI and 3D-FFE were shown to be sensitive to BOLD signal changes in the cervical spinal cord, and had comparable contrast-to-noise ratios in gray matter. FFE exhibited much less signal drop-out and weaker geometric distortions compared to EPI. In the motor paradigm study, the mean number of active voxels was highest in the ventral gray matter horns ipsilateral to the side of the task and at the spinal level associated with innervation of finger extensors.

CONCLUSIONS

Highly multi-shot acquisition sequences such as 3D-FFE are well suited for stimulus-evoked spinal cord BOLD fMRI.

Feasibility of human spinal cord perfusion mapping using dynamic susceptibility contrast imaging at 7T: Preliminary results and identified guidelines

PURPOSE

To explore the feasibility of dynamic susceptibility contrast MRI at 7 Tesla for human spinal cord perfusion mapping and fill the gap between brain and spinal cord perfusion mapping techniques.

METHODS

Acquisition protocols for high-resolution single shot EPI in the spinal cord were optimized for both spin-echo and gradient-echo preparations, including cardiac gating, acquisition times and breathing cycle recording. Breathing-induced MRI signal fluctuations were investigated in healthy volunteers. A specific image- and signal-processing pipeline was implemented to address them. Dynamic susceptibility contrast was then evaluated in 3 healthy volunteers and 5 patients. Bolus depiction on slice-wise signal within cord was investigated, and maps of relative perfusion indices were computed.

RESULTS

Signal fluctuations were increased by 1.9 and 2.3 in free-breathing compared to apnea with spin-echo and gradient-echo, respectively. The ratio between signal fluctuations and bolus peak in healthy volunteers was 5.0% for spin-echo and 3.8% for gradient-echo, allowing clear depiction of the bolus on every slice and yielding relative blood flow and volume maps exhibiting the expected higher perfusion of gray matter. However, signal fluctuations in patients were increased by 4 in average (using spin-echo), compromising the depiction of the bolus in slice-wise signal. Moreover, 3 of 18 slices had to be discarded because of fat-aliasing artifacts.

CONCLUSION

Dynamic susceptibility contrast MRI at 7 Tesla showed great potential for spinal cord perfusion mapping with a reliability never achieved thus far for single subject and single slice measurements. Signal stability needs to be improved in acquisition conditions associated with patients; guidelines to achieve that have been identified and shared.

Cortico-spinal imaging to study pain

Functional magnetic resonance imaging of the brain has helped to reveal mechanisms of pain perception in health and disease. Recently, imaging approaches have been developed that allow recording neural activity simultaneously in the brain and in the spinal cord. These approaches offer the possibility to examine pain perception in the entire central pain system and in addition, to investigate cortico-spinal interactions during pain processing. Although cortico-spinal imaging is a promising technique, it bears challenges concerning data acquisition and data analysis strategies. In this review, we discuss studies that applied simultaneous imaging of the brain and spinal cord to explore central pain processing. Furthermore, we describe different MR-related acquisition techniques, summarize advantages and disadvantages of approaches that have been implemented so far and present software that has been specifically developed for the analysis of spinal fMRI data to address challenges of spinal data analysis.

Dynamic distortion correction for functional MRI using FID navigators

PURPOSE

To develop a method for slice-wise dynamic distortion correction for EPI using rapid spatiotemporal B₀ field measurements from FID navigators (FIDnavs) and to evaluate the efficacy of this new approach relative to an established data-driven technique.

METHODS

A low-resolution reference image was used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B₀ inhomogeneity variations up to second order from measured FIDnavs. Five volunteers were scanned at 3 T using a 64-channel coil with FID-navigated EPI. The accuracy of voxel shift measurements and geometric distortion correction was assessed for experimentally induced magnetic field perturbations. The temporal SNR was evaluated in EPI time-series acquired at rest and with a continuous nose-touching action, before and after image realignment.

RESULTS

Field inhomogeneity coefficients and voxel shift maps measured using FIDnavs were in excellent agreement with multi-echo EPI measurements. The FID-navigated distortion correction accurately corrected image geometry in the presence of induced magnetic field perturbations, outperforming the data-driven approach in regions with large field offsets. In functional MRI scans with nose touching, FIDnav-based correction yielded temporal SNR gains of 30% in gray matter. Following image realignment, which accounted for global image shifts, temporal SNR gains of 3% were achieved.

CONCLUSIONS

Our proposed application of FIDnavs enables slice-wise dynamic distortion correction with high temporal efficiency. We achieved improved signal stability by leveraging the encoding information from multichannel coils. This approach can be easily adapted to other EPI-based sequences to improve temporal SNR for a variety of clinical and research applications.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

Rapid measurement and correction of spatiotemporal B0 field changes using FID navigators and a multi-channel reference image

PURPOSE

To measure spatiotemporal B₀ field changes in real time using FID navigators (FIDnavs) and to demonstrate the efficacy of retrospectively correcting high-resolution T 2 * -weighted images using a novel FIDnav framework.

METHODS

A forward model of the complex FIDnav signals was generated by simulating the effect of changes in the underlying B₀ inhomogeneity coefficients, with spatial encoding provided by a multi-channel reference image. Experiments were performed at 3T to assess the accuracy of B₀ field estimates from FIDnavs acquired from a 64-channel head coil under different shim settings and in 5 volunteers performing deep-breathing and nose-touching tasks designed to modulate the B₀ field. Second-order, in-plane spherical harmonic (SH) inhomogeneity coefficients estimated from FIDnavs were incorporated into an iterative reconstruction to retrospectively correct 2D gradient-echo images acquired in both axial and sagittal planes.

RESULTS

Spatiotemporal B₀ field changes measured from rapidly acquired FIDnavs were in good agreement with the results of second-order SH fitting to the measured field maps. FIDnav field estimates accounted for a significant proportion of the ΔB₀ variance induced by deep breathing (64 ± 21%) and nose touching (67 ± 34%) across all volunteers. Ghosting, blurring, and intensity modulation artifacts in T 2 * -weighted images, induced by spatiotemporal field changes, were visibly reduced following retrospective correction with FIDnav inhomogeneity coefficients.

CONCLUSIONS

Spatially resolved B₀ inhomogeneity changes up to second order can be characterized in real time using the proposed approach. Retrospective FIDnav correction substantially improves T 2 * -weighted image quality in the presence of strong B₀ field modulations, with potential for real-time shimming.

A method for correcting breathing-induced field fluctuations in T2*-weighted spinal cord imaging using a respiratory trace

PURPOSE

Spinal cord MRI at ultrahigh field is hampered by time-varying magnetic fields associated with the breathing cycle, giving rise to ghosting artifacts in multi-shot acquisitions. Here, we suggest a correction approach based on linking the signal from a respiratory bellows to field changes inside the spinal cord. The information is used to correct the data at the image reconstruction level.

METHODS

The correction was demonstrated in the context of multi-shot T2*-weighted imaging of the cervical spinal cord at 7T. A respiratory trace was acquired during a high-resolution multi-echo gradient-echo sequence, used for structural imaging and quantitative T2* mapping, and a multi-shot EPI time series, as would be suitable for fMRI. The coupling between the trace and the breathing-induced fields was determined by a short calibration scan in each individual. Images were reconstructed with and without trace-based correction.

RESULTS

In the multi-echo acquisition, breathing-induced fields caused severe ghosting in images with long TE, which led to a systematic underestimation of T2* in the spinal cord. The trace-based correction reduced the ghosting and increased the estimated T2* values. Breathing-related ghosting was also observed in the multi-shot EPI images. The correction largely removed the ghosting, thereby improving the temporal signal-to-noise ratio of the time series.

CONCLUSIONS

Trace-based retrospective correction of breathing-induced field variations can reduce ghosting and improve quantitative metrics in multi-shot structural and functional T2*-weighted imaging of the spinal cord. The method is straightforward to implement and does not rely on sequence modifications or additional hardware beyond a respiratory bellows.

Spinal Cord Imaging in Amyotrophic Lateral Sclerosis: Historical Concepts-Novel Techniques

Amyotrophic lateral sclerosis (ALS) is the most common adult onset motor neuron disease with no effective disease modifying therapies at present. Spinal cord degeneration is a hallmark feature of ALS, highlighted in the earliest descriptions of the disease by Lockhart Clarke and Jean-Martin Charcot. The anterior horns and corticospinal tracts are invariably affected in ALS, but up to recently it has been notoriously challenging to detect and characterize spinal pathology in vivo. With recent technological advances, spinal imaging now offers unique opportunities to appraise lower motor neuron degeneration, sensory involvement, metabolic alterations, and interneuron pathology in ALS. Quantitative spinal imaging in ALS has now been used in cross-sectional and longitudinal study designs, applied to presymptomatic mutation carriers, and utilized in machine learning applications. Despite its enormous clinical and academic potential, a number of physiological, technological, and methodological challenges limit the routine use of computational spinal imaging in ALS. In this review, we provide a comprehensive overview of emerging spinal cord imaging methods and discuss their advantages, drawbacks, and biomarker potential in clinical applications, clinical trial settings, monitoring, and prognostic roles.

Influence of respiration-induced B0 variations in real-time phase-contrast echo planar imaging of the cervical cerebrospinal fluid

PURPOSE

Respiration induces temporal variations of the main magnetic field B₀ along the spinal cord. These variations are typically not compensated for in velocity quantifications using phase-contrast MRI. The goal of this study was to analyze errors caused by respiration-induced B₀ variations in real-time phase-contrast echo planar imaging (PCEPI) of cervical cerebrospinal fluid (CSF) velocity measurements and to evaluate this effect for various sequence parameters using numerical simulations.

METHODS

Real-time B₀ measurements with double gradient echo sequence and PCEPI measurements were acquired in the cervical CSF of 10 healthy subjects. Dynamic phase offsets attributed to respiration-induced B₀ variations were analyzed by quantifying amplitudes and comparing the temporal behavior with respiratory signals. In experiments and simulations, the influence of the echo time (TE) and the delay between PCEPI images (Δt) with respect to respiration on the dynamic phase offsets were investigated.

RESULTS

A good agreement was found between phase offsets extracted from both acquisition types. Furthermore, respiratory signals qualitatively matched the temporal behavior of the measured phase offsets showing a dependency on subject-dependent local B₀ distribution and respiration physiology. Simulations revealed residual background phases in PCEPI velocity quantification varying with TE and Δt.

CONCLUSION

Respiration-induced B₀ variations result in dynamic background phases in real-time PCEPI velocity quantifications of the CSF in the cervical spine. The current work underlines that these background phases need to be corrected to avoid confounding effects.

Data consistency-driven determination of B 0 -fluctuations in gradient-echo MRI

PURPOSE

Introduce a method to estimate B 0 -fluctuations based on the analysis of raw k-space data, without sequence modifications or external hardware, and correct for their detrimental effects in gradient-echo MRI.

THEORY

Inconsistencies in multi-channel raw k-space data can be used to estimate B 0 -fluctuations by exploiting coil-sensitivity information.

METHODS

The proposed method, dubbed consistency navigation, is used to extract B 0 -fluctuations from T 2 * -weighted 3D gradient-echo data. These results are compared with the results from an MR phase navigator and respiratory bellows. The spatial variation of the B 0 -fluctuation amplitude is derived using the sensitivity maps of the coil array and compared with direct measurements based on dynamic 2D gradient-echo data.

RESULTS

B 0 -fluctuations derived from the consistency navigator and MR phase navigator are highly correlated. Images corrected for these fluctuations show marked improvements in homogeneity and tissue delineation. The spatial variation of the B 0 -fluctuation amplitude follows closely the variation directly measured from time-resolved 2D scans.

CONCLUSIONS

Based on the consistency navigator, an accurate estimation of the spatiotemporal characteristics of B 0 -fluctuations and correction of T 2 * -weighted images has been demonstrated.

Functional Neuroimaging of Nociceptive and Pain-Related Activity in the Spinal Cord and Brain: Insights From Neurovascular Coupling Studies

Spinal cord and brain processes underlie pain perception, which produces systemic cardiovascular changes. In turn, the autonomic nervous system regulates vascular function in the spinal cord and brain in order to adapt to these systemic changes, while neuronal activity induces local vascular changes. Thus, autonomic regulation and pain processes in the brain and spinal cord are tightly linked and interrelated. The objective of this topical review is to discuss work on neurovascular coupling during nociceptive processing in order to highlight supporting evidence and limitations for the use of cerebral and spinal fMRI to investigate pain mechanisms and spinal nociceptive processes. Work on functional neuroimaging of pain is presented and discussed in relation to available neurovascular coupling studies and related issues. Perspectives on future work are also discussed with an emphasis on differences between the brain and the spinal cord and on different approaches that may be useful to improve current methods, data analyses and interpretation. In summary, this review highlights the lack of data on neurovascular coupling during nociceptive stimulation and indicates that hemodynamic and BOLD responses measured with fMRI may be biased by nonspecific vascular changes. Future neuroimaging studies on nociceptive and pain-related processes would gain further understanding of neurovascular coupling in the brain and spinal cord and should take into account the effects of systemic vascular changes that may affect hemodynamic responses. Anat Rec, 301:1585-1595, 2018. © 2018 Wiley Periodicals, Inc.