MR venography of the human brain using susceptibility weighted imaging at very high field strength

Visualizing Cerebral Small Vessel Degeneration During Aging and Diseases Using Magnetic Resonance Imaging

Cerebral small vessel disease is a major contributor to brain disorders in older adults. It is associated with a much higher risk of stroke and dementia. Due to a lack of clinical and fluid biomarkers, diagnosing and grading small vessel disease are highly dependent on magnetic resonance imaging. In the past, researchers mostly used brain parenchymal imaging markers to represent small vessel damage, but the relationships between these surrogate markers and small vessel pathologies are complex. Recent progress in high-resolution magnetic resonance imaging methods, including time-of-flight MR angiography, phase-contrast MR angiography, black blood vessel wall imaging, susceptibility-weighted imaging, and contrast-enhanced methods, allow for direct visualization of cerebral small vessel structures. They could be powerful tools for understanding aging-related small vessel degeneration and improving disease diagnosis and treatment. This article will review progress in these imaging techniques and their application in aging and disease studies. Some challenges and future directions are also discussed. EVIDENCE LEVEL: 4. TECHNICAL EFFICACY: 3.

U-Net based vessel segmentation for murine brains with small micro-magnetic resonance imaging reference datasets

Identification and quantitative segmentation of individual blood vessels in mice visualized with preclinical imaging techniques is a tedious, manual or semiautomated task that can require weeks of reviewing hundreds of levels of individual data sets. Preclinical imaging, such as micro-magnetic resonance imaging (μMRI) can produce tomographic datasets of murine vasculature across length scales and organs, which is of outmost importance to study tumor progression, angiogenesis, or vascular risk factors for diseases such as Alzheimer's. Training a neural network capable of accurate segmentation results requires a sufficiently large amount of labelled data, which takes a long time to compile. Recently, several reasonably automated approaches have emerged in the preclinical context but still require significant manual input and are less accurate than the deep learning approach presented in this paper-quantified by the Dice score. In this work, the implementation of a shallow, three-dimensional U-Net architecture for the segmentation of vessels in murine brains is presented, which is (1) open-source, (2) can be achieved with a small dataset (in this work only 8 μMRI imaging stacks of mouse brains were available), and (3) requires only a small subset of labelled training data. The presented model is evaluated together with two post-processing methodologies using a cross-validation, which results in an average Dice score of 61.34% in its best setup. The results show, that the methodology is able to detect blood vessels faster and more reliably compared to state-of-the-art vesselness filters with an average Dice score of 43.88% for the used dataset.

Laminar functional magnetic resonance imaging in vision research

Magnetic resonance imaging (MRI) scanners at ultra-high magnetic fields have become available to use in humans, thus enabling researchers to investigate the human brain in detail. By increasing the spatial resolution, ultra-high field MR allows both structural and functional characterization of cortical layers. Techniques that can differentiate cortical layers, such as histological studies and electrode-based measurements have made critical contributions to the understanding of brain function, but these techniques are invasive and thus mainly available in animal models. There are likely to be differences in the organization of circuits between humans and even our closest evolutionary neighbors. Thus research on the human brain is essential. Ultra-high field MRI can observe differences between cortical layers, but is non-invasive and can be used in humans. Extensive previous literature has shown that neuronal connections between brain areas that transmit feedback and feedforward information terminate in different layers of the cortex. Layer-specific functional MRI (fMRI) allows the identification of layer-specific hemodynamic responses, distinguishing feedback and feedforward pathways. This capability has been particularly important for understanding visual processing, as it has allowed researchers to test hypotheses concerning feedback and feedforward information in visual cortical areas. In this review, we provide a general overview of successful ultra-high field MRI applications in vision research as examples of future research.

Using high spatial resolution fMRI to understand representation in the auditory network

Following rapid methodological advances, ultra-high field (UHF) functional and anatomical magnetic resonance imaging (MRI) has been repeatedly and successfully used for the investigation of the human auditory system in recent years. Here, we review this work and argue that UHF MRI is uniquely suited to shed light on how sounds are represented throughout the network of auditory brain regions. That is, the provided gain in spatial resolution at UHF can be used to study the functional role of the small subcortical auditory processing stages and details of cortical processing. Further, by combining high spatial resolution with the versatility of MRI contrasts, UHF MRI has the potential to localize the primary auditory cortex in individual hemispheres. This is a prerequisite to study how sound representation in higher-level auditory cortex evolves from that in early (primary) auditory cortex. Finally, the access to independent signals across auditory cortical depths, as afforded by UHF, may reveal the computations that underlie the emergence of an abstract, categorical sound representation based on low-level acoustic feature processing. Efforts on these research topics are underway. Here we discuss promises as well as challenges that come with studying these research questions using UHF MRI, and provide a future outlook.

Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI

Purpose

Susceptibility Weighted Imaging (SWI) has become established in the clinical investigation of stroke, microbleeds, tumor vascularization, calcification and iron deposition, but suffers from a number of shortcomings and artefacts. The goal of this study was to reduce the sensitivity of SWI to strong B₁ and B₀ inhomogeneities at ultra-high field to generate homogeneous images with increased contrast and free of common artefacts. All steps in SWI processing have been addressed −coil combination, phase unwrapping, image combination over echoes, phase filtering and homogeneity correction −and applied to an efficient bipolar multi-echo acquisition to substantially improve the quality of SWI.

Principal results

Our findings regarding the optimal individual processing steps lead us to propose a Contrast-weighted, Laplace-unwrapped, bipolar multi-Echo, ASPIRE-combined, homogeneous, improved Resolution SWI, or CLEAR-SWI. CLEAR-SWI was compared to two other multi-echo SWI methods and standard, single-echo SWI with the same acquisition time at 7 T in 10 healthy volunteers and with single-echo SWI in 13 patients with brain tumors. CLEAR-SWI had improved contrast-to-noise and homogeneity, reduced signal dropout and was not compromised by the artefacts which affected standard SWI in 10 out of 13 cases close to tumors (as assessed by expert raters), as well as generating T2* maps and phase images which can be used for Quantitative Susceptibility Mapping. In a comparison with other multi-echo SWI methods, CLEAR-SWI had the fewest artefacts, highest SNR and generally higher contrast-to-noise.

Major conclusions

CLEAR-SWI eliminates the artefacts common in standard, single-echo SWI, reduces signal dropouts and improves image homogeneity and contrast-to-noise. Applied clinically, in a study of brain tumor patients, CLEAR-SWI was free of the artefacts which affected standard, single-echo SWI.

Penumbra quantification from MR SWI-DWI mismatch and its comparison with MR ASL PWI-DWI mismatch in patients with acute ischemic stroke

In acute-ischemic-stroke patients, penumbra assessment plays a significant role in treatment outcome. MR perfusion-weighted imaging (PWI) and diffusion-weighted imaging (DWI) mismatch ratio can provide penumbra assessment. Recently reported studies have shown the potential of susceptibility-weighted imaging (SWI) in the qualitative assessment of penumbra. We hypothesize that quantitative penumbra assessment using SWI-DWI can provide an alternative to the PWI-DWI approach and this can also reduce the overall scan-time. The purpose of the current study was to develop a framework for accurate quantitative assessment of penumbra using SWI-DWI and its validation with PWI-DWI-based quantification. In the current study, the arterial-spin-labelling (ASL) technique has been used for PWI. This retrospective study included 25 acute-ischemic-stroke patients presenting within 24 hours of the last noted baseline condition of stroke onset. Eleven patients also had follow-up MRI within 48 hours. MRI acquisition comprised DWI, SWI, pseudo-continuous-ASL (pCASL), FLAIR and non-contrast-angiography sequences. A framework was developed for the enhancement of prominent hypo-intense vein signs followed by automatic segmentation of the SWI penumbra ROI. Apparent-diffusion-coefficient (ADC) maps and cerebral-blood-flow (CBF) maps were computed. The infarct core ROI from the ADC map and the ASL penumbra ROI from CBF maps were segmented semiautomatically. The infarct core volume, SWI penumbra volume (SPV) and pCASL penumbra volume were computed and used to calculate mismatch ratios MRSWIADC and MRCBFADC . The Dice coefficient between the SWI penumbra ROI and ASL penumbra ROI was 0.96 ± 0.07. MRSWIADC correlated well (r = 0.90, p < 0.05) with MRCBFADC , which validates the hypothesis of accurate penumbra assessment using the SWI-DWI mismatch ratio. Moreover, a significant association between high SPV and the presence of vessel occlusion in the MR angiogram was observed. Follow-up data showed salvation of penumbra tissue (location and volumes predicted by proposed framework) by treatments. Additionally, functional-outcome analysis revealed 93.3% of patients with MRSWIADC > 1 benefitted from revascularization therapy. Overall, the proposed automated quantitative assessment of penumbra using the SWI-DWI mismatch ratio performs equivalently to the ASL PWI-DWI mismatch ratio. This approach provides an alternative to the perfusion sequence required for penumbra assessment, which can reduce scan time by 17% for the protocol without a perfusion sequence.

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.

Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

10.5 T MRI static field effects on human cognitive, vestibular, and physiological function

PURPOSE

To perform a pilot study to quantitatively assess cognitive, vestibular, and physiological function during and after exposure to a magnetic resonance imaging (MRI) system with a static field strength of 10.5 Tesla at multiple time scales.

METHODS

A total of 29 subjects were exposed to a 10.5 T MRI field and underwent vestibular, cognitive, and physiological testing before, during, and after exposure; for 26 subjects, testing and exposure were repeated within 2-4 weeks of the first visit. Subjects also reported sensory perceptions after each exposure. Comparisons were made between short and long term time points in the study with respect to the parameters measured in the study; short term comparison included pre-vs-isocenter and pre-vs-post (1-24 h), while long term compared pre-exposures 2-4 weeks apart.

RESULTS

Of the 79 comparisons, 73 parameters were unchanged or had small improvements after magnet exposure. The exceptions to this included lower scores on short term (i.e. same day) executive function testing, greater isocenter spontaneous eye movement during visit 1 (relative to pre-exposure), increased number of abnormalities on videonystagmography visit 2 versus visit 1 and a mix of small increases (short term visit 2) and decreases (short term visit 1) in blood pressure. In addition, more subjects reported metallic taste at 10.5 T in comparison to similar data obtained in previous studies at 7 T and 9.4 T.

CONCLUSION

Initial results of 10.5 T static field exposure indicate that 1) cognitive performance is not compromised at isocenter, 2) subjects experience increased eye movement at isocenter, and 3) subjects experience small changes in vital signs but no field-induced increase in blood pressure. While small but significant differences were found in some comparisons, none were identified as compromising subject safety. A modified testing protocol informed by these results was devised with the goal of permitting increased enrollment while providing continued monitoring to evaluate field effects.

Enhanced POCS reconstruction for partial Fourier imaging in multi-echo and time-series acquisitions

PURPOSE

To improve partial Fourier (PF) imaging reconstruction in time-series or multi-echo acquisitions.

METHODS

Many PF methods use a phase estimate to restore Hermitian symmetry before filling missing k-space entries with measured data from the opposite half. This estimate is obtained from the symmetrically sampled, central part of k-space and its low-resolution results in artifacts near high-frequency phase effects (eg, tissue boundaries, vessels), limiting PF undersampling. Enhanced projection onto convex sets (POCS) uses full-resolution phase estimates and relies on alternating the half of k-space that is acquired in time series or multi-echo acquisitions. This enables full-resolution phase estimates to be calculated for each volume/echo, which are fed into the POCS framework. We apply enhanced POCS to high-resolution multi-echo FLASH and 3D-EPI functional MRI time-series data.

RESULTS

Reconstruction errors and their bias dramatically reduce compared with existing methods, without leading to temporal blurring in time-series acquisitions. This allows for higher PF acceleration factors at virtually no cost.

CONCLUSION

Enhanced POCS results in superior PF reconstructions. Furthermore, as the resolution of the phase estimate used for symmetry correction no longer depends on the PF factor, enhanced POCS is more robust against larger PF omission.

First in-vivo human imaging at 10.5T: Imaging the body at 447 MHz

PURPOSE

To investigate the feasibility of imaging the human torso and to evaluate the performance of several radiofrequency (RF) management strategies at 10.5T.

METHODS

Healthy volunteers were imaged on a 10.5T whole-body scanner in multiple target anatomies, including the prostate, hip, kidney, liver, and heart. Phase-only shimming and spoke pulses were used to demonstrate their performance in managing the B 1 + inhomogeneity present at 447 MHz. Imaging protocols included both qualitative and quantitative acquisitions to show the feasibility of imaging with different contrasts.

RESULTS

High-quality images were acquired and demonstrated excellent overall contrast and signal-to-noise ratio. The experimental results matched well with predictions and suggested good translational capabilities of the RF management strategies previously developed at 7T. Phase-only shimming provided increased efficiency, but showed pronounced limitations in homogeneity, demonstrating the need for the increased degrees of freedom made possible through single- and multispoke RF pulse design.

CONCLUSION

The first in-vivo human imaging was successfully performed at 10.5T using previously developed RF management strategies. Further improvement in RF coils, transmit chain, and full integration of parallel transmit functionality are needed to fully realize the benefits of 10.5T.

Structure tensor informed fibre tractography at 3T

Structure tensor informed fibre tractography (STIFT) based on informing tractography for diffusion-weighted images at 3T and by utilising the structure tensor obtained from gradient-recalled echo (GRE) images at 7T is able to delineate fibres when seed voxels are placed close to the fibre boundaries. However, incorporating data from two different field strengths limits the applicability of STIFT. In this study, STIFT was implemented with both diffusion-weighted images and GRE images acquired at 3T. Instead of using the magnitude GRE data directly for STIFT as in the previous work, the utility of T₂ * maps and quantitative susceptibility maps derived from complex-valued GRE data to improve fibre delineation was explored. Single-seed tractography was performed and the results show that the optic radiation reconstructed with STIFT is more distinguishable from the inferior longitudinal fasciculus/inferior fronto-occipital fasciculus complex when compared to standard diffusion-weighted imaging tractography. We further investigated the quantitative effects of STIFT in a group of five healthy volunteers and evaluated its impact on measures of structural connectivity. The framework was extended to evaluate implementations of STIFT based on T₂ *-weighted and quantitative susceptibility-weighted images in a whole-brain connectivity study. In terms of connectivity, no systematic differences were found between STIFT and diffusion-weighted imaging tractography, suggesting that local improvements in tractography are not translated to the atlas-based structural connectivity analysis. Nevertheless, the reduction in the number of statistically significant connections in the STIFT connectivity matrix suggests that STIFT can potentially reduce the false-positive connections in fibre tractography.

Studying brain microstructure with magnetic susceptibility contrast at high-field

A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.

Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field

The SNR and CNR benefits of ultra-high field (UHF) have helped push the envelope of achievable spatial resolution in MRI. For applications based on susceptibility contrast where there is a large CNR gain, high quality sub-millimeter resolution imaging is now being routinely performed, particularly in fMRI and phase imaging/QSM. This has enabled the study of structure and function of very fine-scale structures in the brain. UHF has also helped push the spatial resolution of many other MRI applications as will be outlined in this review. However, this push in resolution comes at a cost of a large encoding burden leading to very lengthy scans. Developments in parallel imaging with controlled aliasing and the move away from 2D slice-by-slice imaging to much more SNR-efficient simultaneous multi-slice (SMS) and 3D acquisitions have helped address this issue. In particular, these developments have revolutionized the efficiency of UHF MRI to enable high spatiotemporal resolution imaging at an order of magnitude faster acquisition. In addition to describing the main approaches to these techniques, this review will also outline important key practical considerations in using these methods in practice. Furthermore, new RF pulse design to tackle the B1+ and SAR issues of UHF and the increased SAR and power requirement of SMS RF pulses will also be touched upon. Finally, an outlook into new developments of smart encoding in more dimensions, particularly through using better temporal/across-contrast encoding and reconstruction will be described. Just as controlled aliasing fully exploits spatial encoding in parallel imaging to provide large multiplicative gains in accelerations, the complimentary use of these new approaches in temporal and across-contrast encoding are expected to provide exciting opportunities for further large gains in efficiency to further push the spatiotemporal resolution of MRI.

Automated vessel exclusion technique for quantitative assessment of hepatic iron overload by R2*-MRI

BACKGROUND

Extraction of liver parenchyma is an important step in the evaluation of R2*-based hepatic iron content (HIC). Traditionally, this is performed by radiologists via whole-liver contouring and T2*-thresholding to exclude hepatic vessels. However, the vessel exclusion process is iterative, time-consuming, and susceptible to interreviewer variability.

PURPOSE

To implement and evaluate an automatic hepatic vessel exclusion and parenchyma extraction technique for accurate assessment of R2*-based HIC.

STUDY TYPE

Retrospective analysis of clinical data.

SUBJECTS

Data from 511 MRI exams performed on 257 patients were analyzed.

FIELD STRENGTH/SEQUENCE

All patients were scanned on a 1.5T scanner using a multiecho gradient echo sequence for clinical monitoring of HIC.

ASSESSMENT

An automated method based on a multiscale vessel enhancement filter was investigated for three input data types-contrast-optimized composite image, T2* map, and R2* map-to segment blood vessels and extract liver tissue for R2*-based HIC assessment. Segmentation and R2* results obtained using this automated technique were compared with those from a reference T2*-thresholding technique performed by a radiologist.

STATISTICAL TESTS

The Dice similarity coefficient was used to compare the segmentation results between the extracted parenchymas, and linear regression and Bland-Altman analyses were performed to compare the R2* results, obtained with the automated and reference techniques.

RESULTS

Mean liver R2* values estimated from all three filter-based methods showed excellent agreement with the reference method (slopes 1.04-1.05, R² > 0.99, P < 0.001). Parenchyma areas extracted using the reference and automated methods had an average overlap area of 87-88%. The T2*-thresholding technique included small vessels and pixels at the vessel/tissue boundaries as parenchymal area, potentially causing a small bias (<5%) in R2* values compared to the automated method.

DATA CONCLUSION

The excellent agreement between reference and automated hepatic vessel segmentation methods confirms the accuracy and robustness of the proposed method. This automated approach might improve the radiologist's workflow by reducing the interpretation time and operator dependence for assessing HIC, an important clinical parameter that guides iron overload management.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;47:1542-1551.

Removal of Arterial Vessel Contributions in Susceptibility-Weighted Images for Quantification of Normalized Visible Venous Volume in Children with Sickle Cell Disease

Purpose

To evaluate a new postprocessing framework that eliminates arterial vessel signal contributions in the quantification of normalized visible venous volume (NVVV, a ratio between venous and brain volume) in susceptibility-weighted imaging (SWI) exams in patients with sickle cell disease (SCD).

Materials and Methods

We conducted a retrospective study and qualitatively reviewed for hypointense arterial vessel contamination in SWI exams from 21 children with SCD. We developed a postprocessing framework using magnetic resonance angiography in combination with SWI to provide a more accurate quantification of NVVV. NVVV was calculated before and after removing arterial vessel contributions to determine the error from hypointense arterial vessels in quantifying NVVV.

Results

Hypointense arterial vessel contamination was observed in 86% SWI exams and was successfully corrected by the proposed method. The contributions of hypointense arterial vessels in the original SWI were significant and accounted for approximately 33% of the NVVV [uncorrected NVVV = 0.012 ± 0.005 versus corrected NVVV = 0.008 ± 0.003 (mean ± SD), P < 0.01].

Conclusion

Hypointense arterial vessel contamination occurred in the majority of SWI exams and led to a sizeable overestimation of the visible venous volume. A prospective longitudinal study is needed to evaluate if quantitation of NVVV was improved and to assess the role of NVVV as a biomarker of SCD severity or stroke risk.

Ultra-high field upper extremity peripheral nerve and non-contrast enhanced vascular imaging

OBJECTIVE

The purpose of this study was to explore the efficacy of Ultra-high field [UHF] 7 Tesla [T] MRI as compared to 3T MRI in non-contrast enhanced [nCE] imaging of structural anatomy in the elbow, forearm, and hand [upper extremity].

MATERIALS AND METHOD

A wide range of sequences including T1 weighted [T1] volumetric interpolate breath-hold exam [VIBE], T2 weighted [T2] double-echo steady state [DESS], susceptibility weighted imaging [SWI], time-of-flight [TOF], diffusion tensor imaging [DTI], and diffusion spectrum imaging [DSI] were optimized and incorporated with a radiofrequency [RF] coil system composed of a transverse electromagnetic [TEM] transmit coil combined with an 8-channel receive-only array for 7T upper extremity [UE] imaging. In addition, Siemens optimized protocol/sequences were used on a 3T scanner and the resulting images from T1 VIBE and T2 DESS were compared to that obtained at 7T qualitatively and quantitatively [SWI was only qualitatively compared]. DSI studio was utilized to identify nerves based on analysis of diffusion weighted derived fractional anisotropy images. Images of forearm vasculature were extracted using a paint grow manual segmentation method based on MIPAV [Medical Image Processing, Analysis, and Visualization].

RESULTS

High resolution and high quality signal-to-noise ratio [SNR] and contrast-to-noise ratio [CNR]-images of the hand, forearm, and elbow were acquired with nearly homogeneous 7T excitation. Measured [performed on the T1 VIBE and T2 DESS sequences] SNR and CNR values were almost doubled at 7T vs. 3T. Cartilage, synovial fluid and tendon structures could be seen with higher clarity in the 7T T1 and T2 weighted images. SWI allowed high resolution and better quality imaging of large and medium sized arteries and veins, capillary networks and arteriovenous anastomoses at 7T when compared to 3T. 7T diffusion weighted sequence [not performed at 3T] demonstrates that the forearm nerves are clearly delineated by fiber tractography. The proper digital palmar arteries and superficial palmar arch could also be clearly visualized using TOF nCE 7T MRI.

CONCLUSION

Ultra-high resolution neurovascular imaging in upper extremities is possible at 7T without use of renal toxic intravenous contrast. 7T MRI can provide superior peripheral nerve [based on fiber anisotropy and diffusion coefficient parameters derived from diffusion tensor/spectrum imaging] and vascular [nCE MRA and vessel segmentation] imaging.

Post Mortem Validation of MRI-Identified Veins on the Surface of the Cerebral Cortex as Potential Landmarks for Neurosurgery

Background and Objective: Image-guided neurosurgery uses information from a wide spectrum of methods to inform the neurosurgeon's judgement about which tissue to resect and which to spare. Imaging data are registered to the patient's skull so that they correspond to the intraoperative macro- and microscopic view. The correspondence between imaging and optical systems breaks down during surgery, however, as a result of cerebro-spinal fluid drain age, tissue resection, and gravity-based brain shift. In this work we investigate whether a map of surface veins, automatically segmented from MRI, could serve as additional reference system.

Methods: Gradient-echo based T2*-weighted imaging was performed on two human cadavers heads using a 7 Tesla MRI scanner. Automatic vessel segmentation was performed using the Frangi vesselness filter, and surface renderings of vessels compared with photographs of the surface of the brain following craniotomy.

Results: A high level of correspondence was established between vessel maps and the post autopsy photographs. Corresponding veins, including the prominent superior anastomotic veins, could be identified in all brain lobes.

Conclusion: Automatic surface vessel segmentation is feasible and the high correspondence to post autopsy photographs indicates that they could be used as an additional reference system for image-guided neurosurgery in order to maintain the correspondence between imaging and optical systems.This has the advantage over a skull-based reference system that veins are clearly visible to the surgeon and move and deform with the underlying tissue, potentially making this surface net of landmarks robust to brain shift.

MAVEN: An Algorithm for Multi-Parametric Automated Segmentation of Brain Veins From Gradient Echo Acquisitions

Cerebral vein analysis provides a chance to study, from an unusual viewpoint, an entire class of brain diseases, including neurodegenerative disorders and traumatic brain injuries. Manual segmentation approaches can be used to assess vascular anatomy, but they are observer-dependent and time-consuming; therefore, automated approaches are desirable, as they also improve reproducibility. In this paper, a new, fully automated algorithm, based on structural, morphological, and relaxometric information, is proposed to segment the entire cerebral venous system from MR images. The algorithm for multi-parametric automated segmentation of brain VEiNs (MAVEN) is based on a combined investigation of multi-parametric information that allows for rejection of false positives and detection of thin vessels. The method is tested on gradient echo brain data sets acquired at 1.5, 3, and 7 T. It is compared to previous methods against manual segmentation, and its inter-scan reproducibility is assessed. The achieved accuracy and reproducibility are good, meaning that MAVEN outperforms previous methods on both quantitative and qualitative analyses. It is usable at all the field strengths explored, showing comparable accuracy scores, with no need for algorithm parameter adjustments, and thus, it is a promising candidate for the characterization of the venous tree topology.

An illustrated comparison of processing methods for MR phase imaging and QSM: combining array coil signals and phase unwrapping

Phase imaging benefits from strong susceptibility effects at very high field and the high signal-to-noise ratio (SNR) afforded by multi-channel coils. Combining the information from coils is not trivial, however, as the phase that originates in local field effects (the source of interesting contrast) is modified by the inhomogeneous sensitivity of each coil. This has historically been addressed by referencing individual coil sensitivities to that of a volume coil, but alternative approaches are required for ultra-high field systems in which no such coil is available. An additional challenge in phase imaging is that the phase that develops up to the echo time is "wrapped" into a range of 2π radians. Phase wraps need to be removed in order to reveal the underlying phase distribution of interest. Beginning with a coil combination using a homogeneous reference volume coil - the Roemer approach - which can be applied at 3 T and lower field strengths, we review alternative methods for combining single-echo and multi-echo phase images where no such reference coil is available. These are applied to high-resolution data acquired at 7 T and their effectiveness assessed via an index of agreement between phase values over channels and the contrast-to-noise ratio in combined images. The virtual receiver coil and COMPOSER approaches were both found to be computationally efficient and effective. The main features of spatial and temporal phase unwrapping methods are reviewed, placing particular emphasis on recent developments in temporal phase unwrapping and Laplacian approaches. The features and performance of these are illustrated in application to simulated and high-resolution in vivo data. Temporal unwrapping was the fastest of the methods tested and the Laplacian the most robust in images with low SNR. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Accelerated mapping of magnetic susceptibility using 3D planes-on-a-paddlewheel (POP) EPI at ultra-high field strength

With the advent of ultra-high field MRI scanners in clinical research, susceptibility based MRI has recently gained increasing interest because of its potential to assess subtle tissue changes underlying neurological pathologies/disorders. Conventional, but rather slow, three-dimensional (3D) spoiled gradient-echo (GRE) sequences are typically employed to assess the susceptibility of tissue. 3D echo-planar imaging (EPI) represents a fast alternative but generally comes with echo-time restrictions, geometrical distortions and signal dropouts that can become severe at ultra-high fields. In this work we assess quantitative susceptibility mapping (QSM) at 7 T using non-Cartesian 3D EPI with a planes-on-a-paddlewheel (POP) trajectory, which is created by rotating a standard EPI readout train around its own phase encoding axis. We show that the threefold accelerated non-Cartesian 3D POP EPI sequence enables very fast, whole brain susceptibility mapping at an isotropic resolution of 1 mm and that the high image quality has sufficient signal-to-noise ratio in the phase data for reliable QSM processing. The susceptibility maps obtained were comparable with regard to QSM values and geometric distortions to those calculated from a conventional 4 min 3D GRE scan using the same QSM processing pipeline. Copyright © 2016 John Wiley & Sons, Ltd.

Overview of quantitative susceptibility mapping

Magnetic susceptibility describes the magnetizability of a material to an applied magnetic field and represents an important parameter in the field of MRI. With the recently introduced method of quantitative susceptibility mapping (QSM) and its conceptual extension to susceptibility tensor imaging (STI), the non-invasive assessment of this important physical quantity has become possible with MRI. Both methods solve the ill-posed inverse problem to determine the magnetic susceptibility from local magnetic fields. Whilst QSM allows the extraction of the spatial distribution of the bulk magnetic susceptibility from a single measurement, STI enables the quantification of magnetic susceptibility anisotropy, but requires multiple measurements with different orientations of the object relative to the main static magnetic field. In this review, we briefly recapitulate the fundamental theoretical foundation of QSM and STI, as well as computational strategies for the characterization of magnetic susceptibility with MRI phase data. In the second part, we provide an overview of current methodological and clinical applications of QSM with a focus on brain imaging. Copyright © 2016 John Wiley & Sons, Ltd.

Ultra-high-field RF coil development for evaluating upper extremity imaging applications

The purpose of this study is to develop and evaluate a custom-designed 7  T MRI coil and explore its use for upper extremity applications. An RF system composed of a transverse electromagnetic transmit coil and an eight-channel receive-only array was developed for 7  T upper extremity applications. The RF system was characterized and evaluated using scattering parameters and B₁⁺ mapping. Finite difference time domain simulations were performed to evaluate the B₁⁺ field distribution and specific absorption rate for the forearm region of the upper extremity. High-resolution 7  T images were acquired and compared with those at 3 T. The simulation and experimental results show very good B₁⁺ field homogeneity across the forearm. High-resolution images of musculotendinous, osseocartilaginous, and neurovascular structures in the upper extremity are presented with T₁ volumetric interpolated breath-hold examination, T₂ double-echo steady state, T₂ * susceptibility weighted imaging (SWI), diffusion tensor imaging, and time-of-flight sequences. Comparison between 3  T and 7  T is shown. Intricate contextual anatomy can be delineated in synovial, fibrocartilaginous, interosseous, and intraosseous trabecular structures of the forearm, as well as palmar and digital vascular anatomy (including microvascular detail in SWI). Ultra-high-field 7  T imaging holds great potential in improving the sensitivity and specificity of upper extremity imaging, especially in wrist and hand pathology secondary to bone, ligament, nerve, vascular, and other soft or hard tissue etiology.

MRI venous architecture of the thalamus

PURPOSE

Our purpose is to describe the thalamic veins using a novel approach named venous gliography in cases with primary or secondary gliomas of the thalamus. Venous gliography is defined by authors as a method to visualize veins on MRI Brain T1-weighted post contrast scans containing gliomas which have induced regional venous congestion.

METHODS

Routine clinical MR Imaging studies were reviewed to assess the presence of thalamic veins in 29 glioma cases. In addition, confocal reconstruction techniques (Anatom-e and Osirix) were used in cases that had thin sections (1.0-1.5mm) post contrast T1 weighted sequences. Multiplanar MIP and confocal volume rendered images were generated to evaluate the thalamic veins in those cases.

RESULTS

Using venous gliography and confocal reconstruction techniques, two patterns in the venous architecture of the thalamus were documented. First, the branching pattern created by the tributaries of the internal cerebral vein, namely the superior thalamic vein and the anterior thalamic vein, which together formed the superior group of thalamic veins. Second, the pattern created by the un-branched vertically oriented veins, namely the inferior thalamic veins and the posterior thalamic veins, which joined the basal vein of Rosenthal and constituted the inferior group of thalamic veins.

CONCLUSIONS

Venous gliography combined with the use of confocal reconstruction techniques provided a novel approach to display the thalamic veins that are usually not seen. The understanding of the venous architecture is mandated by the recent research where veins have taken on an important role in the perivenular spread of gliomas.

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.

Rapid brain MRI acquisition techniques at ultra-high fields

Ultra-high-field MRI provides large increases in signal-to-noise ratio (SNR) as well as enhancement of several contrast mechanisms in both structural and functional imaging. Combined, these gains result in a substantial boost in contrast-to-noise ratio that can be exploited for higher-spatial-resolution imaging to extract finer-scale information about the brain. With increased spatial resolution, however, there is a concurrent increased image-encoding burden that can cause unacceptably long scan times for structural imaging and slow temporal sampling of the hemodynamic response in functional MRI - particularly when whole-brain imaging is desired. To address this issue, new directions of imaging technology development - such as the move from conventional 2D slice-by-slice imaging to more efficient simultaneous multislice (SMS) or multiband imaging (which can be viewed as "pseudo-3D" encoding) as well as full 3D imaging - have provided dramatic improvements in acquisition speed. Such imaging paradigms provide higher SNR efficiency as well as improved encoding efficiency. Moreover, SMS and 3D imaging can make better use of coil sensitivity information in multichannel receiver arrays used for parallel imaging acquisitions through controlled aliasing in multiple spatial directions. This has enabled unprecedented acceleration factors of an order of magnitude or higher in these imaging acquisition schemes, with low image artifact levels and high SNR. Here we review the latest developments of SMS and 3D imaging methods and related technologies at ultra-high field for rapid high-resolution functional and structural imaging of the brain. Copyright © 2016 John Wiley & Sons, Ltd.

Susceptibility-weighted imaging using inter-echo-variance channel combination for improved contrast at 7 tesla

PURPOSE

To implement and optimize a new approach for susceptibility-weighted image (SWI) generation from multi-echo multi-channel image data and compare its performance against optimized traditional SWI pipelines.

MATERIALS AND METHODS

Five healthy volunteers were imaged at 7 Tesla. The inter-echo-variance (IEV) channel combination, which uses the variance of the local frequency shift at multiple echo times as a weighting factor during channel combination, was used to calculate multi-echo local phase shift maps. Linear phase masks were combined with the magnitude to generate IEV-SWI. The performance of the IEV-SWI pipeline was compared with that of two accepted SWI pipelines-channel combination followed by (i) Homodyne filtering (HPH-SWI) and (ii) unwrapping and high-pass filtering (SVD-SWI). The filtering steps of each pipeline were optimized. Contrast-to-noise ratio was used as the comparison metric. Qualitative assessment of artifact and vessel conspicuity was performed and processing time of pipelines was evaluated.

RESULTS

The optimized IEV-SWI pipeline (σ = 7 mm) resulted in continuous vessel visibility throughout the brain. IEV-SWI had significantly higher contrast compared with HPH-SWI and SVD-SWI (P < 0.001, Friedman nonparametric test). Residual background fields and phase wraps in HPH-SWI and SVD-SWI corrupted the vessel signal and/or generated vessel-mimicking artifact. Optimized implementation of the IEV-SWI pipeline processed a six-echo 16-channel dataset in under 10 min.

CONCLUSION

IEV-SWI benefits from channel-by-channel processing of phase data and results in high contrast images with an optimal balance between contrast and background noise removal, thereby presenting evidence of importance of the order in which postprocessing techniques are applied for multi-channel SWI generation.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:1113-1124.

A Comparison Study of Single-Echo Susceptibility Weighted Imaging and Combined Multi-Echo Susceptibility Weighted Imaging in Visualizing Asymmetric Medullary Veins in Stroke Patients

Background

Asymmetric medullary veins (AMV) are frequently observed in stroke patients and single-echo susceptibility weighted imaging (SWI_s) is the main technique in detecting AMV. Our study aimed to investigate which echo time (TE) on single-echo susceptibility is the optimal echo for visualizing AMV and to compare the ability in detecting AMV in stroke patients between SWI_s and multi-echo susceptibility weighted imaging (SWI_c).

Materials and Methods

Twenty patients with middle cerebral artery stroke were included. SWI was acquired by using a multi-echo gradient-echo sequence with six echoes ranging from 5 ms to 35.240 ms. Three different echoes of SWI_s including SWI_s1 (TE = 23.144 ms), SWI_s2 (TE = 29.192 ms) and SWI_s3 (TE = 35.240 ms) were reconstructed. SWI_c was averaged using the three echoes of SWI_s. Image quality and venous contrast of medullary veins were compared between SWI_s and SWI_c using peak signal-to-noise ratio (PSNR), mean opinion score (MOS), contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR). The presence of AMV was evaluated in each SWI_{s (1–3)} and SWI_c.

Results

SWI_s2 had the highest PSNR, MOS and CNR and SWI_s1 had the highest SNR among three different echoes of SWI_s. No significant difference was found in SNR between SWI_s1 and SWI_s2. PSNR, MOS and CNR in SWI_c were significantly increased by 27.9%, 28.2% and 17.2% compared with SWI_s2 and SNR in SWI_c was significantly increased by 32.4% compared with SWI_s1. 55% of patients with AMV were detected in SWI_s2, SWI_s3 and SWI_c, while 50% AMV were found in SWI_s1.

Conclusions

SWI_s using TE around 29ms was optimal in visualizing AMV. SWI_c could improve image quality and venous contrast, but was equal to SWI_s using a relative long TE in evaluating AMV. These results provide the technique basis for further research of AMV in stroke.

Automatic 3D Segmentation and Quantification of Lenticulostriate Arteries from High-Resolution 7 Tesla MRA Images

We propose a novel hybrid approach for automatic 3D segmentation and quantification of high-resolution 7 Tesla magnetic resonance angiography (MRA) images of the human cerebral vasculature. Our approach consists of two main steps. First, a 3D model-based approach is used to segment and quantify thick vessels and most parts of thin vessels. Second, remaining vessel gaps of the first step in low-contrast and noisy regions are completed using a 3D minimal path approach, which exploits directional information. We present two novel minimal path approaches. The first is an explicit approach based on energy minimization using probabilistic sampling, and the second is an implicit approach based on fast marching with anisotropic directional prior. We conducted an extensive evaluation with over 2300 3D synthetic images and 40 real 3D 7 Tesla MRA images. Quantitative and qualitative evaluation shows that our approach achieves superior results compared with a previous minimal path approach. Furthermore, our approach was successfully used in two clinical studies on stroke and vascular dementia.

Combining phase images from array coils using a short echo time reference scan (COMPOSER)

PURPOSE

To develop a simple method for combining phase images from multichannel coils that does not require a reference coil and does not entail phase unwrapping, fitting or iterative procedures.

THEORY AND METHODS

At very short echo time, the phase measured with each coil of an array approximates to the phase offset to which the image from that coil is subject. Subtracting this information from the phase of the scan of interest matches the phases from the coils, allowing them to be combined. The effectiveness of this approach is quantified in the brain, calf and breast with coils of diverse designs.

RESULTS

The quality of phase matching between coil elements was close to 100% with all coils assessed even in regions of low signal. This method of phase combination was similar in effectiveness to the Roemer method (which needs a reference coil) and was superior to the rival reference-coil-free approaches tested.

CONCLUSION

The proposed approach-COMbining Phase data using a Short Echo-time Reference scan (COMPOSER)-is a simple and effective approach to reconstructing phase images from multichannel coils. It requires little additional scan time, is compatible with parallel imaging and is applicable to all coils, independent of configuration. Magn Reson Med 77:318-327, 2017. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a novel MRI based technique that relies on estimates of the magnetic field distribution in the tissue under examination. Several sophisticated data processing steps are required to extract the magnetic field distribution from raw MRI phase measurements. The objective of this review article is to provide a general overview and to discuss several underlying assumptions and limitations of the pre-processing steps that need to be applied to MRI phase data before the final field-to-source inversion can be performed. Beginning with the fundamental relation between MRI signal and tissue magnetic susceptibility this review covers the reconstruction of magnetic field maps from multi-channel phase images, background field correction, and provides an overview of state of the art QSM solution strategies.

Automatic SWI Venography Segmentation Using Conditional Random Fields

Susceptibility-weighted imaging (SWI) venography can produce detailed venous contrast and complement arterial dominated MR angiography (MRA) techniques. However, these dense reversed-contrast SWI venograms pose new segmentation challenges. We present an automatic method for whole-brain venous blood segmentation in SWI using Conditional Random Fields (CRF). The CRF model combines different first and second order potentials. First-order association potentials are modeled as the composite of an appearance potential, a Hessian-based shape potential and a non-linear location potential. Second-order interaction potentials are modeled using an auto-logistic (smoothing) potential and a data-dependent (edge) potential. Minimal post-processing is used for excluding voxels outside the brain parenchyma and visualizing the surface vessels. The CRF model is trained and validated using 30 SWI venograms acquired within a population of deep brain stimulation (DBS) patients (age range [Formula: see text] years). Results demonstrate robust and consistent segmentation in deep and sub-cortical regions (median kappa = 0.84 and 0.82), as well as in challenging mid-sagittal and surface regions (median kappa = 0.81 and 0.83) regions. Overall, this CRF model produces high-quality segmentation of SWI venous vasculature that finds applications in DBS for minimizing hemorrhagic risks and other surgical and non-surgical applications.

Quantitative susceptibility mapping: current status and future directions

Quantitative susceptibility mapping (QSM) is a new technique for quantifying magnetic susceptibility. It has already found various applications in quantifying in vivo iron content, calcifications and changes in venous oxygen saturation. The accuracy of susceptibility mapping is dependent on several factors. In this review, we evaluate the entire process of QSM from data acquisition to individual data processing steps. We also show preliminary results of several new concepts introduced in this review in an attempt to improve the quality and accuracy for certain steps. The uncertainties in estimating susceptibility differences using susceptibility maps, phase images, and T2* maps are analyzed and compared. Finally, example clinical applications are presented. We conclude that QSM holds great promise in quantifying iron and becoming a standard clinical tool.

Evaluation of SWI in children with sickle cell disease

BACKGROUND AND PURPOSE

SWI is a powerful tool for imaging of the cerebral venous system. The SWI venous contrast is affected by blood flow, which may be altered in sickle cell disease. In this study, we characterized SWI venous contrast in patients with sickle cell disease and healthy control participants and examined the relationships among SWI venous contrast, and hematologic variables in the group with sickle cell disease.

MATERIALS AND METHODS

A retrospective review of MR imaging and hematologic variables from 21 patients with sickle cell disease and age- and sex-matched healthy control participants was performed. A Frangi vesselness filter was used to quantify the attenuation of visible veins from the SWI. The normalized visible venous volume was calculated for quantitative analysis of venous vessel conspicuity.

RESULTS

The normalized visible venous volume was significantly lower in the group with sickle cell disease vs the control group (P < .001). Normalized visible venous volume was not associated with hemoglobin, percent hemoglobin F, percent hemoglobin S, absolute reticulocyte count, or white blood cell count. A hypointense arterial signal on SWI was observed in 18 of the 21 patients with sickle cell disease and none of the 21 healthy control participants.

CONCLUSIONS

This study demonstrates the variable and significantly lower normalized visible venous volume in patients with sickle cell disease compared with healthy control participants. Decreased venous contrast in sickle cell disease may reflect abnormal cerebral blood flow, volume, velocity, or oxygenation. Quantitative analysis of SWI contrast may be useful for investigation of cerebrovascular pathology in patients with sickle cell disease, and as a tool to monitor therapies. However, future studies are needed to elucidate physiologic mechanisms of decreased venous conspicuity in sickle cell disease.

Scanning fast and slow: current limitations of 3 Tesla functional MRI and future potential

Functional MRI at 3T has become a workhorse for the neurosciences, e.g., neurology, psychology, and psychiatry, enabling non-invasive investigation of brain function and connectivity. However, BOLD-based fMRI is a rather indirect measure of brain function, confounded by physiology related signals, e.g., head or brain motion, brain pulsation, blood flow, intermixed with susceptibility differences close or distant to the region of neuronal activity. Even though a plethora of preprocessing strategies have been published to address these confounds, their efficiency is still under discussion. In particular, physiological signal fluctuations closely related to brain supply may mask BOLD signal changes related to "true" neuronal activation. Here we explore recent technical and methodological advancements aimed at disentangling the various components, employing fast multiband vs. standard EPI, in combination with fast temporal ICA. Our preliminary results indicate that fast (TR <0.5 s) scanning may help to identify and eliminate physiologic components, increasing tSNR and functional contrast. In addition, biological variability can be studied and task performance better correlated to other measures. This should increase specificity and reliability in fMRI studies. Furthermore, physiological signal changes during scanning may then be recognized as a source of information rather than a nuisance. As we are currently still undersampling the complexity of the brain, even at a rather coarse macroscopic level, we should be very cautious in the interpretation of neuroscientific findings, in particular when comparing different groups (e.g., age, sex, medication, pathology, etc.). From a technical point of view our goal should be to sample brain activity at layer specific resolution with low TR, covering as much of the brain as possible without violating SAR limits. We hope to stimulate discussion toward a better understanding and a more quantitative use of fMRI.

High-resolution anatomy of the human brain stem using 7-T MRI: improved detection of inner structures and nerves?

INTRODUCTION

The purpose of this paper is to assess the value of 7 Tesla (7 T) MRI for the depiction of brain stem and cranial nerve (CN) anatomy.

METHODS

Six volunteers were examined at 7 T using high-resolution SWI, MPRAGE, MP2RAGE, 3D SPACE T2, T2, and PD images to establish scanning parameters targeted at optimizing spatial resolution. Direct comparisons between 3 and 7 T were performed in two additional subjects using the finalized sequences (3 T: T2, PD, MPRAGE, SWAN; 7 T: 3D T2, MPRAGE, SWI, MP2RAGE). Artifacts and the depiction of structures were evaluated by two neuroradiologists using a standardized score sheet.

RESULTS

Sequences could be established for high-resolution 7 T imaging even in caudal cranial areas. High in-plane resolution T2, PD, and SWI images provided depiction of inner brain stem structures such as pons fibers, raphe, reticular formation, nerve roots, and periaqueductal gray. MPRAGE and MP2RAGE provided clear depiction of the CNs. 3D T2 images improved depiction of inner brain structure in comparison to T2 images at 3 T. Although the 7-T SWI sequence provided improved contrast to some inner structures, extended areas were influenced by artifacts due to image disturbances from susceptibility differences.

CONCLUSIONS

Seven-tesla imaging of basal brain areas is feasible and might have significant impact on detection and diagnosis in patients with specific diseases, e.g., trigeminal pain related to affection of the nerve root. Some inner brain stem structures can be depicted at 3 T, but certain sequences at 7 T, in particular 3D SPACE T2, are superior in producing anatomical in vivo images of deep brain stem structures.

Phase reconstruction from multiple coil data using a virtual reference coil

PURPOSE

This study develops a method to obtain optimal estimates of absolute magnetization phase from multiple-coil MRI data.

THEORY AND METHODS

The element-specific phases of a multi-element receiver coil array are accounted for by using the phase of a real or virtual reference coil that is sensitive over the entire imaged volume. The virtual-reference coil is generated as a weighted combination of measurements from all receiver coils. The phase-corrected multiple coil complex images are combined using the inverse covariance matrix. These methods are tested on images of an agar phantom, an in vivo breast, and an anesthetized rabbit obtained using combinations of four, nine, and three receiver channels, respectively.

RESULTS

The four- and three-channel acquisitions require formation of a virtual-reference receiver coil while one channel of the nine-channel receive array has a sensitivity profile covering the entire imaged volume. Referencing to a real or virtual coil gives receiver phases that are essentially identical except for the individual receiver channel noise. The resulting combined images, which account for receiver channel noise covariance, show the expected reduction in phase variance.

CONCLUSION

The proposed virtual reference coil method determines a phase distribution for each coil from which an optimal phase map can be obtained.

A new susceptibility-weighted image reconstruction method for the reduction of background phase artifacts

PURPOSE

To significantly reduce the background phase effects, especially at the air-tissue interface, and to enhance the desirable local structures of veins in susceptibility-weighted imaging.

METHODS

In the proposed reconstruction method called Magnitude of Complex Filtering, a complex-valued magnetic resonance image is acquired using a flow-compensated high-resolution 3D gradient-echo sequence and the magnitude of the complex-valued image is set to 1 so that the phase information, which contains details of the local susceptibility, is emphasized. Then, the nonlinear filter of the Magnitude of Complex Filtering method is applied to the complex-valued image with a constant magnitude. This filter utilizes the magnitude of the low-pass and high-pass filtered complex data to selectively reduce the background phase effects while enhancing the local structures. The filter output is then processed to generate a susceptibility-weighted image.

RESULTS

Compared with the conventional susceptibility-weighted images generated by a homodyne high-pass filter, the susceptibility-weighted images from the proposed Magnitude of Complex Filtering method show significant improvement; the undesirable artifacts at the air-tissue interface regions and the brain boundaries are significantly reduced, while the contrast of the local structures of veins is enhanced.

CONCLUSION

The Magnitude of Complex Filtering method successfully reduced most background phase effects without requiring additional processing or scan time.

Computer-aided detection of radiation-induced cerebral microbleeds on susceptibility-weighted MR images

Recent interest in exploring the clinical relevance of cerebral microbleeds (CMBs) has motivated the search for a fast and accurate method to detect them. Visual inspection of CMBs on MR images is a lengthy, arduous task that is highly prone to human error because of their small size and wide distribution throughout the brain. Several computer-aided CMB detection algorithms have recently been proposed in the literature, but their diagnostic accuracy, computation time, and robustness are still in need of improvement. In this study, we developed and tested a semi-automated method for identifying CMBs on minimum intensity projected susceptibility-weighted MR images that are routinely used in clinical practice to visually identify CMBs. The algorithm utilized the 2D fast radial symmetry transform to initially detect putative CMBs. Falsely identified CMBs were then eliminated by examining geometric features measured after performing 3D region growing on the potential CMB candidates. This algorithm was evaluated in 15 patients with brain tumors who exhibited CMBs on susceptibility-weighted images due to prior external beam radiation therapy. Our method achieved heightened sensitivity and acceptable amount of false positives compared to prior methods without compromising computation speed. Its superior performance and simple, accelerated processing make it easily adaptable for detecting CMBs in the clinic and expandable to a wide array of neurological disorders.