Whole-brain three-dimensional T2-weighted BOLD functional magnetic resonance imaging at 7 Tesla

Whole-cerebrum guanidino and amide CEST mapping at 3 T by a 3D stack-of-spirals gradient echo acquisition

PURPOSE

To develop a 3D, high-sensitivity CEST mapping technique based on the 3D stack-of-spirals (SOS) gradient echo readout, the proposed approach was compared with conventional acquisition techniques and evaluated for its efficacy in concurrently mapping of guanidino (Guan) and amide CEST in human brain at 3 T, leveraging the polynomial Lorentzian line-shape fitting (PLOF) method.

METHODS

Saturation time and recovery delay were optimized to achieve maximum CEST time efficiency. The 3DSOS method was compared with segmented 3D EPI (3DEPI), turbo spin echo, and gradient- and spin-echo techniques. Image quality, temporal SNR (tSNR), and test-retest reliability were assessed. Maps of Guan and amide CEST derived from 3DSOS were demonstrated on a low-grade glioma patient.

RESULTS

The optimized recovery delay/saturation time was determined to be 1.4/2 s for Guan and amide CEST. In addition to nearly doubling the slice number, the gradient echo techniques also outperformed spin echo sequences in tSNR: 3DEPI (193.8 ± 6.6), 3DSOS (173.9 ± 5.6), and GRASE (141.0 ± 2.7). 3DSOS, compared with 3DEPI, demonstrated comparable GuanCEST signal in gray matter (GM) (3DSOS: [2.14%-2.59%] vs. 3DEPI: [2.15%-2.61%]), and white matter (WM) (3DSOS: [1.49%-2.11%] vs. 3DEPI: [1.64%-2.09%]). 3DSOS also achieves significantly higher amideCEST in both GM (3DSOS: [2.29%-3.00%] vs. 3DEPI: [2.06%-2.92%]) and WM (3DSOS: [2.23%-2.66%] vs. 3DEPI: [1.95%-2.57%]). 3DSOS outperforms 3DEPI in terms of scan-rescan reliability (correlation coefficient: 3DSOS: 0.58-0.96 vs. 3DEPI: -0.02 to 0.75) and robustness to motion as well.

CONCLUSION

The 3DSOS CEST technique shows promise for whole-cerebrum CEST imaging, offering uniform contrast and robustness against motion artifacts.

Submillimeter fMRI Acquisition Techniques for Detection of Laminar and Columnar Level Brain Activation

Since the first demonstration in the early 1990s, functional MRI (fMRI) has emerged as one of the most powerful, noninvasive neuroimaging tools to probe brain functions. Subsequently, fMRI techniques have advanced remarkably, enabling the acquisition of functional signals with a submillimeter voxel size. This innovation has opened the possibility of investigating subcortical neural activities with respect to the cortical depths or cortical columns. For this purpose, numerous previous works have endeavored to design suitable functional contrast mechanisms and dedicated imaging techniques. Depending on the choice of the functional contrast, functional signals can be detected with high sensitivity or with improved spatial specificity to the actual activation site, and the pertaining issues have been discussed in a number of earlier works. This review paper primarily aims to provide an overview of the subcortical fMRI techniques that allow the acquisition of functional signals with a submillimeter resolution. Here, the advantages and disadvantages of the imaging techniques will be described and compared. We also summarize supplementary imaging techniques that assist in the analysis of the subcortical brain activation for more accurate mapping with reduced geometric deformation. This review suggests that there is no single universally accepted method as the gold standard for subcortical fMRI. Instead, the functional contrast and the corresponding readout imaging technique should be carefully determined depending on the purpose of the study. Due to the technical limitations of current fMRI techniques, most subcortical fMRI studies have only targeted partial brain regions. As a future prospect, the spatiotemporal resolution of fMRI will be pushed to satisfy the community's need for a deeper understanding of whole-brain functions and the underlying connectivity in order to achieve the ultimate goal of a time-resolved and layer-specific spatial scale. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

Enabling brain-wide mapping of directed functional connectivity at 3T via layer-dependent fMRI with draining-vein suppression

Layer-dependent functional magnetic resonance imaging (fMRI) offers a compelling avenue for investigating directed functional connectivity (FC). To construct a comprehensive map of brain-wide directed FC, several technical criteria must be met, including sub-mm spatial resolution, adequate temporal resolution, functional sensitivity, global brain coverage, and high spatial specificity. Although gradient echo (GE)-based echo planar imaging (EPI) is commonly used for rapid fMRI acquisition, it faces significant challenges due to the draining-vein effect, particularly when utilizing blood oxygen level-dependent (BOLD) contrast. In this study, we mitigated this effect by incorporating velocity-nulling (VN) gradients into a GE-BOLD fMRI sequence, opting for a 3T magnetic field strength over 7T. We also integrated several advanced techniques, such as simultaneous multi-slice (SMS) acceleration and NORDIC denoising, to enhance temporal resolution, spatial coverage, and signal sensitivity. Collectively, the VN fMRI method exhibited notable spatial specificity, as evidenced by the identification of double-peak activation patterns within the primary motor cortex (M1) during a finger-tapping task. Additionally, the technique demonstrated BOLD sensitivity in the lateral geniculate nucleus (LGN). Furthermore, our VN fMRI technique displayed superior robustness when compared to conventional fMRI approaches across participants. Our findings of directed FC elucidate several layer-specific functional relationships between different brain regions and align closely with existing literature. Given the widespread availability of 3T scanners, this technical advancement has the potential for significant impact across multiple domains of neuroscience research.

Differential Laminar Activation Dissociates Encoding and Retrieval in the Human Medial and Lateral Entorhinal Cortex

The hierarchically organized structures of the medial temporal lobe are critically important for episodic memory function. Accumulating evidence suggests dissociable information processing pathways are maintained throughout these structures including in the medial and lateral entorhinal cortex. Cortical layers provide an additional dimension of dissociation as the primary input to the hippocampus derives from layer 2 neurons in the entorhinal cortex, whereas the deeper layers primarily receive output from the hippocampus. Here, novel high-resolution T2-prepared functional MRI methods were successfully used to mitigate susceptibility artifacts typically affecting MRI signals in this region providing uniform sensitivity across the medial and lateral entorhinal cortex. During the performance of a memory task, healthy human subjects (age 25-33 years, mean age 28.2 ± 3.3 years, 4 female) showed differential functional activation in the superficial and deep layers of the entorhinal cortex associated with task-related encoding and retrieval conditions, respectively. The methods provided here offer an approach to probe layer-specific activation in normal cognition and conditions contributing to memory impairment.SIGNIFICANCE STATEMENT This study provides new evidence for differential neuronal activation in the superficial versus deep layers of the entorhinal cortex associated with encoding and retrieval memory processes, respectively, in cognitively normal adults. The study further shows that this dissociation can be observed in both the medial and the lateral entorhinal cortex. The study was achieved by using a novel functional MRI method allowing us to measure robust functional MRI signals in both the medial and lateral entorhinal cortex that was not possible in previous studies. The methodology established here in healthy human subjects lays a solid foundation for subsequent studies investigating layer-specific and region-specific changes in the entorhinal cortex associated with memory impairment in various conditions such as Alzheimer's disease.

Evaluation of T2-prepared blood oxygenation level dependent functional magnetic resonance imaging with an event-related task: Hemodynamic response function and reproducibility

T2-prepared (T2prep) blood oxygenation level dependent (BOLD) functional MRI (fMRI) is an alternative fMRI approach developed to mitigate the susceptibility artifacts that are typically observed in brain regions near air-filled cavities, bleeding and calcification, and metallic objects in echo-planar-imaging (EPI) based fMRI images. Here, T2prep BOLD fMRI was evaluated in an event-related paradigm for the first time. Functional experiments were performed using gradient-echo (GRE) EPI, spin-echo (SE) EPI, and T2prep BOLD fMRI during an event-related visual task in 10 healthy human subjects. Each fMRI method was performed with a low (3.4 × 3.4 × 4 mm3) and a high (1.5 mm isotropic) spatial resolution on 3T and a high resolution (1.5 mm isotropic) on 7T. Robust activation were detected in the visual cortex with all three fMRI methods. In each group of fMRI scans (3T low resolution, 3T high resolution, and 7T high resolution), GRE EPI showed the highest signal change (ΔS/S), largest full-width-at-half-maximum (FWHM) and longest time-to-peak (TTP) extracted from the hemodynamic response functions (HRF), indicating substantial signal contribution from large draining veins which have longer response times than microvessels. In contrast, T2prep BOLD showed the lowest ΔS/S, smallest FWHM, and shortest TTP, suggesting that T2prep BOLD may have a purer T2-weighted BOLD contrast that is more sensitive to microvessels compared to GRE/SE EPI BOLD. This trend was more obvious in fMRI scans performed with a lower spatial resolution on a lower field (3T with a 3.4 × 3.4 × 4 mm3 voxel). Scan-rescan reproducibility in the same subjects was comparable among the three fMRI methods. The results from the current study are expected to be useful to establish T2prep BOLD as a useful alternative fMRI approach for event-related fMRI in brain regions with large susceptibility artifacts.

A revisit of the k-space filtering effects of magnetization-prepared 3D FLASH and balanced SSFP acquisitions: Analytical characterization of the point spread functions

PURPOSE

3D FLASH and balanced SSFP (bSSFP) are increasingly used in quantitative MRI after contrast preparation. The acquired k-space data are modulated by T₁ relaxation (or additional T₂ for bSSFP). Three separate sequence parameters including the number of phase-encoding steps per shot (N), flip angle (FA), and TR have made the transient state of rapid gradient echo (GRE) imaging difficult for analysis and optimization. Here we aim to analytically characterize the k-space filtering effect of magnetization-prepared FLASH and bSSFP with the point spread functions (PSF).

METHODS

The amplitude effect is characterized with the peak magnitude of the PSF, i.e. PSF(0), which, due to their approaching from transient state to steady-state for the GRE acquisitions, obeys a linear (with a slope and an intercept, not proportional) relationship with the prepared longitudinal magnetization (M_prep). The blurring effect is characterized by the FWHM of the PSF. The magnetization-prepared acquisition-dependent image contrast efficiency is characterized with the relative contrast-to-noise ratio (CNR) per unit time (ruCNR).

RESULTS

The slope of PSF(0) characterizes the relative contrast between different M_prep levels. The intercept of PSF(0) could lead to quantification bias for magnetization-prepared imaging. FLASH and bSSFP experience very little blurring effect, which is to the contrary of conventional fast spin echo (FSE). Analytical selections of N, FA, and TR are provided to optimize ruCNR for different scenarios.

CONCLUSIONS

PSFs of the FLASH and bSSFP acquisitions are analytically derived and numerically validated, and compared with the FSE acquisition, thus providing a useful tool for optimizing magnetization-prepared GRE acquisitions.

Three-dimensional reduced field-of-view imaging (3D-rFOVI)

PURPOSE

This study aimed at developing a 3D reduced field-of-view imaging (3D-rFOVI) technique using a 2D radiofrequency (RF) pulse, and demonstrating its ability to achieve isotropic high spatial resolution and reduced image distortion in echo planar imaging (EPI).

METHODS

The proposed 3D-rFOVI technique takes advantage of a 2D RF pulse to excite a slab along the conventional slice-selection direction (i.e., z-direction) while limiting the spatial extent along the phase-encoded direction (i.e., y-direction) within the slab. The slab is phase-encoded in both through-slab and in-slab phase-encoded directions. The 3D-rFOVI technique was implemented at 3T in gradient-echo and spin-echo EPI pulse sequences for functional MRI (fMRI) and diffusion-weighted imaging (DWI), respectively. 3D-rFOVI experiments were performed on a phantom and human brain to illustrate image distortion reduction, as well as isotropic high spatial resolution, in comparison with 3D full-FOV imaging.

RESULTS

In both the phantom and the human brain, image voxel dislocation was substantially reduced by 3D-rFOVI when compared with full-FOV imaging. In the fMRI experiment with visual stimulation, 3D isotropic spatial resolution of (2 × 2 × 2 mm³ ) was achieved with an adequate signal-to-noise ratio (81.5) and blood oxygen level-dependent (BOLD) contrast (2.5%). In the DWI experiment, diffusion-weighted brain images with an isotropic resolution of (1 × 1 × 1 mm³ ) was obtained without appreciable image distortion.

CONCLUSION

This study indicates that 3D-rFOVI is a viable approach to 3D neuroimaging over a zoomed region.

Quantification of cerebral blood volume changes caused by visual stimulation at 3 T using DANTE-prepared dual-echo EPI

PURPOSE

We investigate the influence of moving blood-attenuation effects when using "delay alternating with nutation for tailored excitation" (DANTE) pulses in conjunction with blood oxygen level dependent (BOLD) of functional MRI (fMRI) at 3 T. Based on the effects of including DANTE pulses, we propose quantification of cerebral blood volume (CBV) changes following functional stimulation.

METHODS

Eighteen volunteers in total underwent fMRI scans at 3 T. Seven volunteers were scanned to investigate the effects of DANTE pulses on the fMRI signal. CBV changes in response to visual stimulation were quantified in 11 volunteers using a DANTE-prepared dual-echo EPI sequence.

RESULTS

The inflow effects from flowing blood in arteries and draining vein effects from flowing blood in large veins can be suppressed by use of a DANTE preparation module. Using DANTE-prepared dual-echo EPI, we quantitatively measured intravascular-weighted microvascular CBV changes of 25.4%, 29.8%, and 32.6% evoked by 1, 5, and 10 Hz visual stimulation, respectively. The extravascular fraction (∆S/S)_extra at TE = 30 ms in total BOLD signal was determined to be 64.8 ± 3.4%, which is in line with previous extravascular component estimation at 3 T. Results show that the microvascular CBV changes are linearly dependent on total BOLD changes at TE = 30 ms with a slope of 0.113, and this relation is independent of stimulation frequency and subject.

CONCLUSION

The DANTE preparation pulses can be incorporated into a standard EPI fMRI sequence for the purpose of minimizing inflow effects and reducing draining veins effects in large vessels. Additionally, the DANTE-prepared dual-echo EPI sequence is a promising fast imaging tool for quantification of intravascular-weighted CBV change in the microvascular space at 3 T.

Functional Activities Detected in the Olfactory Bulb and Associated Olfactory Regions in the Human Brain Using T2-Prepared BOLD Functional MRI at 7T

Olfaction is a fundamental sense that plays a vital role in daily life in humans, and can be altered in neuropsychiatric and neurodegenerative diseases. Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) using conventional echo-planar-imaging (EPI) based sequences can be challenging in brain regions important for olfactory processing, such as the olfactory bulb (OB) and orbitofrontal cortex, mainly due to the signal dropout and distortion artifacts caused by large susceptibility effects from the sinonasal cavity and temporal bone. To date, few studies have demonstrated successful fMRI in the OB in humans. T2-prepared (T2prep) BOLD fMRI is an alternative approach developed especially for performing fMRI in regions affected by large susceptibility artifacts. The purpose of this technical study is to evaluate T2prep BOLD fMRI for olfactory functional experiments in humans. Olfactory fMRI scans were performed on 7T in 14 healthy participants. T2prep BOLD showed greater sensitivity than GRE EPI BOLD in the OB, orbitofrontal cortex and the temporal pole. Functional activation was detected using T2prep BOLD in the OB and associated olfactory regions. Habituation effects and a bi-phasic pattern of fMRI signal changes during olfactory stimulation were observed in all regions. Both positively and negatively activated regions were observed during olfactory stimulation. These signal characteristics are generally consistent with literature and showed a good intra-subject reproducibility comparable to previous human BOLD fMRI studies. In conclusion, the methodology demonstrated in this study holds promise for future olfactory fMRI studies in the OB and other brain regions that suffer from large susceptibility artifacts.

Laminar fMRI using T2-prepared multi-echo FLASH

Functional magnetic resonance imaging (fMRI) using blood oxygenation level dependent (BOLD) contrast at a sub-millimeter scale is a promising technique to probe neural activity at the level of cortical layers. While gradient echo (GRE) BOLD sequences exhibit the highest sensitivity, their signal is confounded by unspecific extravascular (EV) and intravascular (IV) effects of large intracortical ascending veins and pial veins leading to a downstream blurring effect of local signal changes. In contrast, spin echo (SE) fMRI promises higher specificity towards signal changes near the microvascular compartment. However, the T₂-weighted signal is typically sampled with a gradient echo readout imposing additional T₂'-weighting. In this work, we used a T₂-prepared (T₂-prep) sequence with short GRE readouts to investigate its capability to acquire laminar fMRI data during a visual task in humans at 7 T. By varying the T₂-prep echo time (TE_prep) and acquiring multiple gradient echoes (TE_GRE) per excitation, we studied the specificity of the sequence and the influence of possible confounding contributions to the shape of laminar fMRI profiles. By fitting and extrapolating the multi-echo GRE data to a TE_GRE = 0 ms condition, we show for the first time laminar profiles free of T₂'-pollution, confined to gray matter. This finding is independent of TE_prep, except for the shortest one (31 ms) where hints of a remaining intravascular component can be seen. For TE_GRE > 0 ms a prominent peak at the pial surface is observed that increases with longer TE_GRE and dominates the shape of the profiles independent of the amount of T₂-weighting. Simulations show that the peak at the pial surface is a result of static EV dephasing around pial vessels in CSF visible in GM due to partial voluming. Additionally, another, weaker, static dephasing effect is observed throughout all layers of the cortex, which is particularly obvious in the data with shortest T₂-prep echo time. Our simulations show that this cannot be explained by intravascular dephasing but that it is likely caused by extravascular effects of the intracortical and pial veins. We conclude that even for TE_GRE as short as 2.3 ms, the T₂'-weighting added to the T₂-weighting is enough to dramatically affect the laminar specificity of the BOLD signal change. However, the bulk of this corruption stems from CSF partial volume effects which can in principle be addressed by increasing the spatial resolution of the acquisition.

Impaired response of cerebral oxygen metabolism to visual stimulation in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disease caused by a CAG triplet repeat expansion in the Huntingtin gene. Metabolic and microvascular abnormalities in the brain may contribute to early physiological changes that subserve the functional impairments in HD. This study is intended to investigate potential abnormality in dynamic changes in cerebral blood volume (CBV) and cerebral blood flow (CBF), and cerebral metabolic rate of oxygen (CMRO2) in the brain in response to functional stimulation in premanifest and early manifest HD patients. A recently developed 3-D-TRiple-acquisition-after-Inversion-Preparation magnetic resonance imaging (MRI) approach was used to measure dynamic responses in CBV, CBF, and CMRO2 during visual stimulation in one single MRI scan. Experiments were conducted in 23 HD patients and 16 healthy controls. Decreased occipital cortex CMRO2 responses were observed in premanifest and early manifest HD patients compared to controls (P < 0.001), correlating with the CAG-Age Product scores in these patients (R2 = 0.4, P = 0.001). The results suggest the potential value of this reduced CMRO2 response during visual stimulation as a biomarker for HD and may illuminate the role of metabolic alterations in the pathophysiology of HD.

Research on BOLD-fMRI Data Denoising Based on Bayesian Estimation and Adaptive Wavelet Threshold

The acquisition of functional magnetic resonance imaging (fMRI) images of blood oxygen level-dependent (BOLD) effect and the signals to be analyzed is based on weak changes in the magnetic field caused by small changes in blood oxygen physiological levels, which are weak signals and complex in noise. In order to model and analyze the pathological and hemodynamic parameters of BOLD-fMRI images effectively, it is urgent to use effective signal analysis techniques to reduce the interference of noise and artifacts. In this paper, the noise characteristics of functional magnetic resonance imaging and the traditional signal denoising methods are analyzed. The Bayesian decision criterion takes into account the probability of the total occurrence of all kinds of references and the loss caused by misjudgment and has strong discriminability. So, an improved adaptive wavelet threshold denoising method based on Bayesian estimation is proposed. By using the correlation characteristics of multiscale wavelet coefficients, the corresponding wavelet components of useful signals and noises are processed differently; while retaining useful frequency information, the noise is weakened to the greatest extent. The new adaptive threshold wavelet denoising method based on Bayesian estimation is applied to the actual experiment, and the results of OEF (oxygen extraction fraction) are optimized. A series of simulation experiments are carried out to verify the effectiveness of the proposed method.

A novel spectrally selective fat saturation pulse design with robustness to B0 and B1 inhomogeneities: A demonstration on 3D T1-weighted breast MRI at 3 T

PURPOSE

Spectrally selective fat saturation (FatSat) sequence is commonly used to suppress signal from adipose tissue. Conventional SINC-shaped pulses are sensitive to B₀ off-resonance and B₁⁺ offset. Uniform fat saturation with large spatial coverage is especially challenging for the body and breast MRI. The aim of this study is to develop spectrally selective FatSat pulses that offer more immunity to B₀/B₁⁺ field inhomogeneities than SINC pulses and evaluate them in bilateral breast imaging at 3 T.

MATERIALS AND METHODS

Optimized composite pulses (OCP) were designed based on the optimal control theory with robustness to a targeted B₀/ B₁⁺ conditions. OCP pulses also allows flexible flip angles to meet different requirements. Comparisons with the vendor-provided SINC pulses were conducted by numerical simulation and in vivo scans using a 3D T₁-weighted (T₁w) gradient-echo (GRE) sequence with coverage of the whole-breast.

RESULTS

Simulation revealed that OCP pulses yielded almost half of the transition band and much less sensitivity to B₁⁺ inhomogeneity compared to SINC pulses with B₀ off-resonance within ±200 Hz and B₁⁺ scale error within ±0.3 (P < 0.001). Across five normal subjects, OCP FatSat pulses produced 25-41% lower residual fat signals (P < 0.05) with 27-36% less spatial variation (P < 0.05) than SINC.

CONCLUSION

In contrast to conventional SINC-shaped pulses, the newly designed OCP FatSat pulses mitigated challenges of wide range of B₀/ B₁⁺ field inhomogeneities and achieved more uniform fat suppression in bilateral breast T₁w imaging at 3 T.

Whole-Brain Functional and Diffusion Tensor MRI in Human Participants with Metallic Orthodontic Braces

Background MRI performed with echo-planar imaging (EPI) sequences is sensitive to susceptibility artifacts in the presence of metallic objects, which presents a substantial barrier for performing functional MRI and diffusion tensor imaging (DTI) in patients with metallic orthodontic material and other head implants. Purpose To evaluate the ability to reduce susceptibility artifacts in healthy human participants wearing metallic orthodontic braces for two alternative approaches: T2-prepared functional MRI and diffusion-prepared DTI with three-dimensional fast gradient-echo readout. Materials and Methods In this prospective study conducted from February to September 2018, T2-prepared functional MRI and diffusion-prepared DTI were performed in healthy human participants. Removable dental braces with bonding trays were used so that MRI could be performed with braces and without braces in the same participants. Results were evaluated in regions with strong (EPI dropout regions for functional MRI and the inferior fronto-occipital fasciculus for DTI) and minimal (motor cortex for functional MRI and the posterior limb of internal capsule for DTI) susceptibility artifacts. Signal-to-noise ratio (SNR), contrast-to-noise ratio for functional MRI, apparent diffusion coefficient and fractional anisotropy for DTI, and degree of distortion (quantified with the Jaccard index, which measures the similarity of geometric shapes) were compared in regions with strong or minimal susceptibility effects between the current standard EPI sequences and the proposed alternatives by using paired t test. Results Six participants were evaluated (mean age ± standard deviation, 40 years ± 6; three women). In brain regions with strong susceptibility effects from the metallic braces, T2-prepared functional MRI showed significantly higher SNR (37.8 ± 2.4 vs 15.5 ± 5.3; P < .001) and contrast-to-noise ratio (0.83 ± 0.16 vs 0.29 ± 0.10; P < .001), whereas diffusion-prepared DTI showed higher SNR (5.8 ± 1.5 vs 3.8 ± 0.7; P = .03) than did conventional EPI methods. Apparent diffusion coefficient and fractional anisotropy were consistent with the literature. Geometric distortion was substantially reduced throughout the brain with the proposed methods (significantly higher Jaccard index, 0.95 ± 0.12 vs 0.81 ± 0.61; P < .001). Conclusion T2-prepared functional MRI and diffusion-prepared diffusion tensor imaging can acquire functional and diffusion MRI, respectively, in healthy human participants wearing metallic dental braces with less susceptibility artifacts and geometric distortion than with conventional echo-planar imaging. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Dietrich in this issue.

The Additional Diagnostic Value of the Three-dimensional Volume Rendering Imaging in Routine Radiology Practice

Three-dimensional volume rendering (3DVR) is useful in a wide variety of medical-imaging applications. The increasingly advanced capabilities of CT and MRI to acquire volumetric data sets with isotropic voxels have resulted in the increased use of the 3DVR techniques for clinical applications. The two most commonly used techniques are the maximum intensity projection (MIP) and, more recently, 3DVR. Several kinds of medical imaging data could be reconstructed for 3D display, including CT, MRI, and ultrasonography (US). In particular, the 3D CT imaging has been developed, improved, and widely used of late. Understanding the mechanisms of 3DVR is essential for the accurate evaluation of the resulting images. Although further research is required to detect the efficiency of 3DVR in radiological applications, with wider availability and improved diagnostic performance, 3DVR is likely to enjoy widespread acceptance in the radiology practice going forward.

In vivo magnetic resonance imaging and spectroscopy. Technological advances and opportunities for applications continue to abound

Over the past decades, the field of in vivo magnetic resonance (MR) has built up an impressive repertoire of data acquisition and analysis technologies for anatomical, functional, physiological, and molecular imaging, the description of which requires many book volumes. As such it is impossible for a few authors to have an authoritative overview of the field and for a brief article to be inclusive. We will therefore focus mainly on data acquisition and attempt to give some insight into the principles underlying current advanced methods in the field and the potential for further innovation. In our view, the foreseeable future is expected to show continued rapid progress, for instance in imaging of microscopic tissue properties in vivo, assessment of functional and anatomical connectivity, higher resolution physiologic and metabolic imaging, and even imaging of receptor binding. In addition, acquisition speed and information content will continue to increase due to the continuous development of approaches for parallel imaging (including simultaneous multi-slice imaging), compressed sensing, and MRI fingerprinting. Finally, artificial intelligence approaches are becoming more realistic and will have a tremendous effect on both acquisition and analysis strategies. Together, these developments will continue to provide opportunity for scientific discovery and, in combination with large data sets from other fields such as genomics, allow the ultimate realization of precision medicine in the clinic.

Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications

Quantitative cerebral blood volume (CBV) fMRI has the potential to overcome several specific limitations of BOLD fMRI. It provides direct physiological interpretability and promises superior localization specificity in applications of sub-millimeter resolution fMRI applications at ultra-high magnetic fields (7T and higher). Non-invasive CBV fMRI using VASO (vascular space occupancy), however, is inherently limited with respect to its data acquisition efficiency, restricting its imaging coverage and achievable spatial and temporal resolution. This limitation may be reduced with recent advanced acceleration and reconstruction strategies that allow two-dimensional acceleration, such as in simultaneous multi-slice (SMS) 2D-EPI or 3D-EPI in combination with CAIPIRINHA field-of-view shifting. In this study, we sought to determine the functional sensitivity and specificity of these readout strategies with VASO over a broad range of spatial resolutions; spanning from low spatial resolution (3mm) whole-cortex to sub-millimeter (0.75mm) slab-of-cortex (for cortical layer-dependent applications). In the thermal-noise-dominated regime of sub-millimeter resolutions, 3D-EPI-VASO provides higher temporal stability and sensitivity to detect changes in CBV compared to 2D-EPI-VASO. In this regime, 3D-EPI-VASO unveils task activation located in the cortical laminae with little contamination from surface veins, in contrast to the cortical surface weighting of GE-BOLD fMRI. In the physiological-noise-dominated regime of lower resolutions, however, 2D-SMS-VASO shows superior performance compared to 3D-EPI-VASO. Due to its superior sensitivity at a layer-dependent level, 3D-EPI VASO promises to play an important role in future neuroscientific applications of layer-dependent fMRI.

High-Resolution CBV-fMRI Allows Mapping of Laminar Activity and Connectivity of Cortical Input and Output in Human M1

Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases.

Shuffled magnetization-prepared multicontrast rapid gradient-echo imaging

PURPOSE

To develop a novel acquisition and reconstruction method for magnetization-prepared 3-dimensional multicontrast rapid gradient-echo imaging, using Hankel matrix completion in combination with compressed sensing and parallel imaging.

METHODS

A random k-space shuffling strategy was implemented in simulation and in vivo human experiments at 7 T for 3-dimensional inversion recovery, T₂ /diffusion preparation, and magnetization transfer imaging. We combined compressed sensing, based on total variation and spatial-temporal low-rank regularizations, and parallel imaging with pixel-wise Hankel matrix completion, allowing the reconstruction of tens of multicontrast 3-dimensional images from 3- or 6-min scans.

RESULTS

The simulation result showed that the proposed method can reconstruct signal-recovery curves in each voxel and was robust for typical in vivo signal-to-noise ratio with 16-times acceleration. In vivo studies achieved 4 to 24 times accelerations for inversion recovery, T₂ /diffusion preparation, and magnetization transfer imaging. Furthermore, the contrast was improved by resolving pixel-wise signal-recovery curves after magnetization preparation.

CONCLUSIONS

The proposed method can improve acquisition efficiencies for magnetization-prepared MRI and tens of multicontrast 3-dimensional images could be recovered from a single scan. Furthermore, it was robust against noise, applicable for recovering multi-exponential signals, and did not require any previous knowledge of model parameters. Magn Reson Med 79:62-70, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Language Mapping Using T2-Prepared BOLD Functional MRI in the Presence of Large Susceptibility Artifacts-Initial Results in Patients With Brain Tumor and Epilepsy

At present, presurgical functional mapping is the most prevalent clinical application of functional magnetic resonance imaging (fMRI). Signal dropouts and distortions caused by susceptibility effects in the current standard echo planar imaging (EPI)-based fMRI images are well-known problems and pose a major hurdle for the application of fMRI in several brain regions, many of which are related to language mapping in presurgical planning. Such artifacts are particularly problematic in patients with previous surgical resection cavities, craniotomy hardware, hemorrhage, and vascular malformation. A recently developed T2-prepared (T2prep) fMRI approach showed negligible distortion and dropouts in the entire brain even in the presence of large susceptibility effects. Here, we present initial results comparing T2prep- and multiband EPI-fMRI scans for presurgical language mapping using a sentence completion task in patients with brain tumor and epilepsy. In all patients scanned, T2prep-fMRI showed minimal image artifacts (distortion and dropout) and greater functional sensitivity than EPI-fMRI around the lesions containing blood products and in air-filled cavities. This enhanced sensitivity in T2prep-fMRI was also evidenced by the fact that functional activation during the sentence completion task was detected with T2prep-fMRI but not with EPI-fMRI in the affected areas with the same statistical threshold, whereas cerebrovascular reactivity during a breath-hold task was preserved in these same regions, implying intact neurovascular coupling in these patients. Although further investigations are required to validate these findings with invasive methods such as direct cortical stimulation mapping as the gold standard, this approach provides an alternative method for performing fMRI in brain regions with large susceptibility effects.

The impact of ultra-high field MRI on cognitive and computational neuroimaging

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

Memory performance-related dynamic brain connectivity indicates pathological burden and genetic risk for Alzheimer's disease

BACKGROUND

The incidence of Alzheimer's disease (AD) strongly relates to advanced age and progressive deposition of cerebral amyloid-beta (Aβ), hyperphosphorylated tau, and iron. The purpose of this study was to investigate the relationship between cerebral dynamic functional connectivity and variability of long-term cognitive performance in healthy, elderly subjects, allowing for local pathology and genetic risk.

METHODS

Thirty seven participants (mean (SD) age 74 (6.0) years, Mini-Mental State Examination 29.0 (1.2)) were dichotomized based on repeated neuropsychological test performance within 2 years. Cerebral Aβ was measured by ¹¹C Pittsburgh Compound-B positron emission tomography, and iron by quantitative susceptibility mapping magnetic resonance imaging (MRI) at an ultra-high field strength of 7 Tesla (7T). Dynamic functional connectivity patterns were investigated by resting-state functional MRI at 7T and tested for interactive effects with genetic AD risk (apolipoprotein E (ApoE)-ε4 carrier status).

RESULTS

A relationship between low episodic memory and a lower expression of anterior-posterior connectivity was seen (F(9,27) = 3.23, p < 0.008), moderated by ApoE-ε4 (F(9,27) = 2.22, p < 0.005). Inherent node-strength was related to local iron (F(5,30) = 13.2; p < 0.022).

CONCLUSION

Our data indicate that altered dynamic anterior-posterior brain connectivity is a characteristic of low memory performance in the subclinical range and genetic risk for AD in the elderly. As the observed altered brain network properties are associated with increased local iron, our findings may reflect secondary neuronal changes due to pathologic processes including oxidative stress.

7-T MRI in Cerebrovascular Diseases: Challenges to Overcome and Initial Results

Magnetic resonance imaging (MRI) plays a key role in the investigation of cerebrovascular diseases. Compared with computed tomography (CT) and digital subtraction angiography (DSA), its advantages in diagnosing cerebrovascular pathology include its superior tissue contrast, its ability to visualize blood vessels without the use of a contrast agent, and its use of magnetic fields and radiofrequency pulses instead of ionizing radiation. In recent years, ultrahigh field MRI at 7 tesla (7 T) has shown promise in the diagnosis of many cerebrovascular diseases. The increased signal-to-noise ratio (SNR; 2.3x and 4.7x increase compared with 3 and 1.5 T, respectively) and contrast-to-noise ratio (CNR) at this higher field strength can be exploited to obtain a higher spatial resolution and higher lesion conspicuousness, enabling assessment of smaller brain structures and lesions. Cerebrovascular diseases can be assessed at different tissue levels; for instance, changes of the arteries feeding the brain can be visualized to determine the cause of ischemic stroke, regional changes in brain perfusion can be mapped to predict outcome after revascularization, and tissue damage, including old and recent ischemic infarcts, can be evaluated as a marker of ischemic burden. For the purpose of this review, we will discriminate 3 levels of assessment of cerebrovascular diseases using MRI: Pipes, Perfusion, and Parenchyma (3 Ps). The term Pipes refers to the brain-feeding arteries from the heart and aortic arch, upwards to the carotid arteries, vertebral arteries, circle of Willis, and smaller intracranial arterial branches. Perfusion is the amount of blood arriving at the brain tissue level, and includes the vascular reserve and perfusion territories. Parenchyma refers to the acute and chronic burden of brain tissue damage, which includes larger infarcts, smaller microinfarcts, and small vessel disease manifestations such as white matter lesions, lacunar infarcts, and microbleeds. In this review, we will describe the key developments in the last decade of 7-T MRI of cerebrovascular diseases, subdivided for these 3 levels of assessment.

Colocalization of cerebral iron with Amyloid beta in Mild Cognitive Impairment

Quantitative Susceptibility Mapping (QSM) MRI at 7 Tesla and 11-Carbon Pittsburgh-Compound-B PET were used for investigating the relationship between brain iron and Amyloid beta (Aβ) plaque-load in a context of increased risk for Alzheimer's disease (AD), as reflected by the Apolipoprotein E ε4 (APOE-e4) allele and mild cognitive impairment (MCI) in elderly subjects. Carriers of APOE-e4 with normal cognition had higher cortical Aβ-plaque-load than non-carriers. In MCI an association between APOE-e4 and higher Aβ-plaque-load was observable both for cortical and subcortical brain-regions. APOE-e4 and MCI was also associated with higher cortical iron. Moreover, cerebral iron significantly affected functional coupling, and was furthermore associated with increased Aβ-plaque-load (R²-adjusted = 0.80, p < 0.001) and APOE-e4 carrier status (p < 0.001) in MCI. This study confirms earlier reports on an association between increased brain iron-burden and risk for neurocognitive dysfunction due to AD, and indicates that disease-progression is conferred by spatial colocalization of brain iron deposits with Aβ-plaques.

Variable flip angle 3D-GRASE for high resolution fMRI at 7 tesla

PURPOSE

To evaluate the use of variable flip angle refocusing pulse trains in single-shot three-dimensional gradient and spin-echo (3D-GRASE) to reduce blurring and increase the spatial coverage for high spatial resolution T2 -weighted functional MRI at 7 Tesla.

METHODS

Variable flip angle refocusing schemes in 3D-GRASE were calculated based on extended phase graph theory. The blurring along the slice (partition) direction was evaluated in simulations, as well as phantom and in vivo experiments. Furthermore, temporal stability and functional sensitivity at 0.8 mm isotropic resolution were assessed.

RESULTS

Variable flip angle refocusing schemes yielded significantly reduced blurring compared with conventional refocusing schemes, with the full width at half maximum being approximately 30-40% narrower. Simultaneously, spatial coverage could be increased by 80%. The temporal signal-to-noise ratio was slightly reduced, but functional sensitivity was largely maintained due to increased functional contrast in the variable flip angle acquisitions. Signal-to-noise ratio and functional sensitivity were reduced more strongly in areas with insufficient radiofrequency transmission indicating higher sensitivity to experimental imperfections.

CONCLUSION

Variable flip angle refocusing schemes increase usability of 3D-GRASE for high-resolution functional MRI by reducing blurring and allowing increased spatial coverage. Magn Reson Med 76:897-904, 2016. © 2015 Wiley Periodicals, Inc.

Recent applications of UHF-MRI in the study of human brain function and structure: a review

The increased availability of ultra-high-field (UHF) MRI has led to its application in a wide range of neuroimaging studies, which are showing promise in transforming fundamental approaches to human neuroscience. This review presents recent work on structural and functional brain imaging, at 7 T and higher field strengths. After a short outline of the effects of high field strength on MR images, the rapidly expanding literature on UHF applications of blood-oxygenation-level-dependent-based functional MRI is reviewed. Structural imaging is then discussed, divided into sections on imaging weighted by relaxation time, including quantitative relaxation time mapping, phase imaging and quantitative susceptibility mapping, angiography, diffusion-weighted imaging, and finally magnetization-transfer imaging. The final section discusses studies using the high spatial resolution available at UHF to identify explicit links between structure and function. Copyright © 2015 John Wiley & Sons, Ltd.

Joint design of large-tip-angle parallel RF pulses and blipped gradient trajectories

PURPOSE

To design multichannel large-tip-angle kT-points and spokes radiofrequency (RF) pulses and gradient waveforms for transmit field inhomogeneity compensation in high field magnetic resonance imaging.

THEORY AND METHODS

An algorithm to design RF subpulse weights and gradient blip areas is proposed to minimize a magnitude least-squares cost function that measures the difference between realized and desired state parameters in the spin domain, and penalizes integrated RF power. The minimization problem is solved iteratively with interleaved target phase updates, RF subpulse weights updates using the conjugate gradient method with optimal control-based derivatives, and gradient blip area updates using the conjugate gradient method. Two-channel parallel transmit simulations and experiments were conducted in phantoms and human subjects at 7 T to demonstrate the method and compare it to small-tip-angle-designed pulses and circularly polarized excitations.

RESULTS

The proposed algorithm designed more homogeneous and accurate 180° inversion and refocusing pulses than other methods. It also designed large-tip-angle pulses on multiple frequency bands with independent and joint phase relaxation. Pulses designed by the method improved specificity and contrast-to-noise ratio in a finger-tapping spin echo blood oxygen level dependent functional magnetic resonance imaging study, compared with circularly polarized mode refocusing.

CONCLUSION

A joint RF and gradient waveform design algorithm was proposed and validated to improve large-tip-angle inversion and refocusing at ultrahigh field.

Slice-selective adiabatic magnetization T2 -preparation (SAMPA) for efficient T2 -weighted imaging at ultrahigh field strengths

PURPOSE

At high field, T₂ -weighted (T₂ w) imaging is limited by transmit field inhomogeneity and increased radiofrequency power deposition. In this work, we introduce SAMPA (Slice-selective Adiabatic Magnetization T₂ PrepAration) and demonstrate its use for efficient brain T₂ w imaging at 7 Tesla (T).

METHODS

SAMPA was designed by subsampling an optimized B₁ insensitive rotation (BIR4) waveform with small tip angle linear subpulses. To perform T₂ w imaging, SAMPA was inserted before a fast gradient echo acquisition. The off-resonance behavior, B₁ robustness, and slice selectivity of the novel T₂ preparation module were analyzed using Bloch simulations. The performance of SAMPA for T₂ w imaging was demonstrated in phantom experiments as well as in the brains of healthy volunteers at 7T.

RESULTS

Based on simulations, the proposed design operates at peak B₁ of 15 μT and higher, within a 400 Hz bandwidth. T₂ values were in strong agreement with spin echo-based T₂ mapping in phantom experiments. Whole brain, interleaved multislab three-dimensional imaging could be acquired with 0.8 mm³ isotropic resolution in 5:36 min per T₂ weighting.

CONCLUSION

Compared with previous adiabatic T₂ preparation techniques, SAMPA allows for slice-selectivity, which can lead to efficient and robust acquisitions for T₂ w imaging at high field. Magn Reson Med 76:1741-1749, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

Simultaneous multi-slice Turbo-FLASH imaging with CAIPIRINHA for whole brain distortion-free pseudo-continuous arterial spin labeling at 3 and 7 T

Simultaneous multi-slice (SMS) or multiband (MB) imaging has recently been attempted for arterial spin labeled (ASL) perfusion MRI in conjunction with echo-planar imaging (EPI) readout. It was found that SMS-EPI can reduce the T1 relaxation effect of the label and improve image coverage and resolution with little penalty in signal-to-noise ratio (SNR). However, EPI still suffers from geometric distortion and signal dropout from field inhomogeneity effects especially at high and ultrahigh magnetic fields. Here we present a novel scheme for achieving high fidelity distortion-free quantitative perfusion imaging by combining pseudo-continuous ASL (pCASL) with SMS Turbo-FLASH (TFL) readout at both 3 and 7 T. Bloch equation simulation was performed to characterize and optimize the TFL-based pCASL perfusion signal. Two MB factors (3 and 5) were implemented in SMS-TFL pCASL and compared with standard 2D TFL and EPI pCASL sequences. The temporal SNR of SMS-TFL pCASL relative to that of standard TFL pCASL was 0.76 ± 0.10 and 0.74 ± 0.11 at 7 T and 0.70 ± 0.05 and 0.65 ± 0.05 at 3T for MB factor of 3 and 5, respectively. By implementing background suppression in conjunction with SMS-TFL at 3T, the relative temporal SNR improved to 0.84 ± 0.09 and 0.79 ± 0.10 for MB factor of 3 and 5, respectively. Compared to EPI pCASL, significantly increased temporal SNR (p<0.001) and improved visualization of orbitofrontal cortex were achieved using SMS-TFL pCASL. By combining SMS acceleration with TFL pCASL, we demonstrated the feasibility for whole brain distortion-free quantitative mapping of cerebral blood flow at high and ultrahigh magnetic fields.

Three-dimensional acquisition of cerebral blood volume and flow responses during functional stimulation in a single scan

In addition to the BOLD scan, quantitative functional MRI studies require measurement of both cerebral blood volume (CBV) and flow (CBF) dynamics. The ability to detect CBV and CBF responses in a single additional scan would shorten the total scan time and reduce temporal variations. Several approaches for simultaneous CBV and CBF measurement during functional MRI experiments have been proposed in two-dimensional (2D) mode covering one to three slices in one repetition time (TR). Here, we extended the principles from previous work and present a three-dimensional (3D) whole-brain MRI approach that combines the vascular-space-occupancy (VASO) and flow-sensitive alternating inversion recovery (FAIR) arterial spin labeling (ASL) techniques, allowing the measurement of CBV and CBF dynamics, respectively, in a single scan. 3D acquisitions are complicated for such a scan combination as the time to null blood signal during a steady state needs to be known. We estimated this using Bloch simulations and demonstrate that the resulting 3D acquisition can detect activation patterns and relative signal changes of quality comparable to that of the original separate scans. The same was found for temporal signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). This approach provides improved acquisition efficiency when both CBV and CBF responses need to be monitored during a functional task.