Anisotropic field-of-views in radial imaging

Electric potential energy optimized 3D radial sampling trajectories for MRI

A novel method for creating "golden" 3D center-out radial MRI sampling trajectories was developed and analyzed. This method, called ELECTRO (ELECTRic potential energy Optimized), uses repulsive forces to minimize electric potential energy. An objective function [Formula: see text], the electric potential energies of all subsets of consecutive readouts in a 3D radial trajectory, and its reduced form were minimized using a multi-stage optimization strategy. A metric called normalized mean nearest neighbor angular distance (NMNA) was proposed for describing distributions of points on a sphere. ELECTRO and other relevant golden trajectories were compared in silico using NMNA and point spread function analysis. Consecutive readouts from an ELECTRO trajectory were well spread out, with consistent NMNA values across sphere sizes (σNMNA = 0.005) and between regions on the sphere (NMNA ≈ 1.49). Conversely, the supergolden trajectory had poor consistency in NMNA values (σNMNA = 0.090) and clustering (NMNA = 1.28 at the pole with 40,000 readouts) that lead to artifact in the point spread function. Multi-stage optimization was faster than single-stage and obtained lower values of [Formula: see text] (e.g., 0.87 vs. 0.91, for a sphere size of 40). In conclusion, ELECTRO trajectories are more golden than other 3D center-out radial trajectories, making them a suitable candidate for dynamic 3D MR imaging.

Whole-brain T2 mapping with radial sampling and retrospective motion correction at 3T

PURPOSE

Several barriers prevent the use of whole-brain T2 mapping in routine use despite increasing interest in this parameter. One of the main barriers is the long scan time resulting in patient discomfort and motion corrupted data. To address this challenge, a method for accurate whole-brain T2 mapping with a limited acquisition time and motion correction capabilities is investigated.

METHODS

A 3D radial multi-echo spin-echo sequence was implemented with optimized sampling trajectory enabling the estimation of intra-scan motion, subsequently used to correct the raw data. Motion corrected echo images are then reconstructed with linear subspace constrained reconstruction. Experiments were carried out on phantom and volunteers at 3T to evaluate the accuracy of the T2 estimation, the sensitivity to lesions and the efficiency of the correction on motion corrupted data.

RESULTS

Whole-brain T2 mapping acquired in less than 7 min enabled the depiction of lesions in the white matter with longer T2. Data retrospectively corrupted with typical motion traces of pediatric patients highly benefited from the motion correction by reducing the error in T2 estimates within the lesions. All datasets acquired on seven volunteers, with deliberate motion, also showed that motion corrupted T2 maps could be improved with the retrospective motion correction both at the voxel level and the structure level.

CONCLUSION

A whole-brain T2 mapping sequence with retrospective intra-scan motion correction and reasonable acquisition time is proposed. The method necessitates advanced iterative reconstruction strategies but no additional navigator, external device, or increased scan time is required.

High Spatiotemporal Resolution Radial Encoding Single-Vessel fMRI

High-field preclinical functional MRI (fMRI) is enabled the high spatial resolution mapping of vessel-specific hemodynamic responses, that is single-vessel fMRI. In contrast to investigating the neuronal sources of the fMRI signal, single-vessel fMRI focuses on elucidating its vascular origin, which can be readily implemented to identify vascular changes relevant to vascular dementia or cognitive impairment. However, the limited spatial and temporal resolution of fMRI is hindered hemodynamic mapping of intracortical microvessels. Here, the radial encoding MRI scheme is implemented to measure BOLD signals of individual vessels penetrating the rat somatosensory cortex. Radial encoding MRI is employed to map cortical activation with a focal field of view (FOV), allowing vessel-specific functional mapping with 50 × 50 µm2 in-plane resolution at a 1 to 2 Hz sampling rate. Besides detecting refined hemodynamic responses of intracortical micro-venules, the radial encoding-based single-vessel fMRI enables the distinction of fMRI signals from vessel and peri-vessel voxels due to the different contribution of intravascular and extravascular effects.

Dynamic pulmonary MRI using motion-state weighted motion-compensation (MostMoCo) reconstruction with ultrashort TE: A structural and functional study

PURPOSE

To improve the quality of structural images and the quantification of ventilation in free-breathing dynamic pulmonary MRI.

METHODS

A 3D radial ultrashort TE (UTE) sequence with superior-inferior navigators was used to acquire pulmonary data during free breathing. All acquired data were binned into different motion states according to the respiratory signal extracted from superior-inferior navigators. Motion-resolved images were reconstructed using eXtra-Dimensional (XD) UTE reconstruction. The initial motion fields were generated by registering images at each motion state to other motion states in motion-resolved images. A motion-state weighted motion-compensation (MostMoCo) reconstruction algorithm was proposed to reconstruct the dynamic UTE images. This technique, termed as MostMoCo-UTE, was compared with XD-UTE and iterative motion-compensation (iMoCo) on a porcine lung and 10 subjects.

RESULTS

MostMoCo reconstruction provides higher peak SNR (37.0 vs. 35.4 and 34.2) and structural similarity (0.964 vs. 0.931 and 0.947) compared to XD-UTE and iMoCo in the porcine lung experiment. Higher apparent SNR and contrast-to-noise ratio are achieved using MostMoCo in the human experiment. MostMoCo reconstruction better preserves the temporal variations of signal intensity of parenchyma compared to iMoCo, shows reduced random noise and improved sharpness of anatomical structures compared to XD-UTE. In the porcine lung experiment, the quantification of ventilation using MostMoCo images is more accurate than that using XD-UTE and iMoCo images.

CONCLUSION

The proposed MostMoCo-UTE provides improved quality of structural images and quantification of ventilation for free-breathing pulmonary MRI. It has the potential for the detection of structural and functional disorders of the lung in clinical settings.

Motion corrected silent ZTE neuroimaging

Purpose

To develop self‐navigated motion correction for 3D silent zero echo time (ZTE) based neuroimaging and characterize its performance for different types of head motion.

Methods

The proposed method termed MERLIN (Motion Estimation & Retrospective correction Leveraging Interleaved Navigators) achieves self‐navigation by using interleaved 3D phyllotaxis k‐space sampling. Low resolution navigator images are reconstructed continuously throughout the ZTE acquisition using a sliding window and co‐registered in image space relative to a fixed reference position. Rigid body motion corrections are then applied retrospectively to the k‐space trajectory and raw data and reconstructed into a final, high‐resolution ZTE image.

Results

MERLIN demonstrated successful and consistent motion correction for magnetization prepared ZTE images for a range of different instructed motion paradigms. The acoustic noise response of the self‐navigated phyllotaxis trajectory was found to be only slightly above ambient noise levels (<4 dBA).

Conclusion

Silent ZTE imaging combined with MERLIN addresses two major challenges intrinsic to MRI (i.e., subject motion and acoustic noise) in a synergistic and integrated manner without increase in scan time and thereby forms a versatile and powerful framework for clinical and research MR neuroimaging applications.

Real-time deep artifact suppression using recurrent U-Nets for low-latency cardiac MRI

PURPOSE

Real-time low latency MRI is performed to guide various cardiac interventions. Real-time acquisitions often require iterative image reconstruction strategies, which lead to long reconstruction times. In this study, we aim to reconstruct highly undersampled radial real-time data with low latency using deep learning.

METHODS

A 2D U-Net with convolutional long short-term memory layers is proposed to exploit spatial and preceding temporal information to reconstruct highly accelerated tiny golden radial data with low latency. The network was trained using a dataset of breath-hold CINE data (including 770 time series from 7 different orientations). Synthetic paired data were created by retrospectively undersampling the magnitude images, and the network was trained to recover the target images. In the spirit of interventional imaging, the network was trained and tested for varying acceleration rates and orientations. Data were prospectively acquired and reconstructed in real time in 1 healthy subject interactively and in 3 patients who underwent catheterization. Images were visually compared to sliding window and compressed sensing reconstructions and a conventional Cartesian real-time sequence.

RESULTS

The proposed network generalized well to different acceleration rates and unseen orientations for all considered metrics in simulated data (less than 4% reduction in structural similarity index compared to similar acceleration and orientation-specific networks). The proposed reconstruction was demonstrated interactively, successfully depicting catheters in vivo with low latency (39 ms, including 19 ms for deep artifact suppression) and an image quality comparing favorably to other reconstructions.

CONCLUSION

Deep artifact suppression was successfully demonstrated in the time-critical application of non-Cartesian real-time interventional cardiac MR.

A study of 3D radial density adapted trajectories for sodium imaging

Sodium imaging typically employs ultrashort echo time radial, density adapted and cones trajectories to capture the rapidly decaying short T2 signal. The present study considers the parameter choices involved in the use of these trajectories in terms of their impact on the resolution and signal to noise ratio. Many parameters have a strong effect on these image properties, particularly the number of spokes which impacts voxel size. The present article develops an understanding of the trade-offs involved and how to choose optimal (or at least reasonable) parameter values. This has a practical role in designing clinical protocols for imaging sodium.

Optimizing trajectory ordering for fast radial ultra-short TE (UTE) acquisitions

PURPOSE

Additional spoiler gradients are required in 3D UTE sequences with random view ordering to suppress magnetization refocusing. By leveraging the encoding gradient induced spoiling effect, the spoiler gradients could potentially be reduced or removed to shorten the TR and increase encoding efficiency. An analysis framework is built that models the gradient spoiling effects and a new ordering scheme is proposed for fast 3D UTE acquisition.

THEORY AND METHODS

UTE signal evolution and spatial encoding gradient induced spoiling effect are derived from the Bloch equations. And the concept is validated in 2D radial UTE simulation. Then an optimized ordering scheme, named reordered 2D golden angle (r2DGA) scheme, for 3D UTE acquisition is proposed. The r2DGA scheme is compared to the sequential and 3D golden angle schemes in both phantom and volunteer studies.

RESULTS

The proposed r2DGA ordering scheme was applied to two applications, single breath-holding and free breathing 3D lung MRI. With r2DGA ordering scheme, breath-holding lung MRI scan increased 60% scan efficiency by removing the spoiler gradients and the free breathing scan reduced 20% scan time compared to the 3D golden angle scheme by reducing the spoiler gradients.

CONCLUSIONS

The proposed r2DGA ordering scheme UTE acquisition reduces the need of spoiler gradients and increases the encoding efficiency, and shows improvements in both breath-holding and free breathing lung MRI applications.

Simultaneous measurements of 3D wall shear stress and pulse wave velocity in the murine aortic arch

PURPOSE

Wall shear stress (WSS) and pulse wave velocity (PWV) are important parameters to characterize blood flow in the vessel wall. Their quantification with flow-sensitive phase-contrast (PC) cardiovascular magnetic resonance (CMR), however, is time-consuming. Furthermore, the measurement of WSS requires high spatial resolution, whereas high temporal resolution is necessary for PWV measurements. For these reasons, PWV and WSS are challenging to measure in one CMR session, making it difficult to directly compare these parameters. By using a retrospective approach with a flexible reconstruction framework, we here aimed to simultaneously assess both PWV and WSS in the murine aortic arch from the same 4D flow measurement.

METHODS

Flow was measured in the aortic arch of 18-week-old wildtype (n = 5) and ApoE-/- mice (n = 5) with a self-navigated radial 4D-PC-CMR sequence. Retrospective data analysis was used to reconstruct the same dataset either at low spatial and high temporal resolution (PWV analysis) or high spatial and low temporal resolution (WSS analysis). To assess WSS, the aortic lumen was labeled by semi-automatically segmenting the reconstruction with high spatial resolution. WSS was determined from the spatial velocity gradients at the lumen surface. For calculation of the PWV, segmentation data was interpolated along the temporal dimension. Subsequently, PWV was quantified from the through-plane flow data using the multiple-points transit-time method. Reconstructions with varying frame rates and spatial resolutions were performed to investigate the influence of spatiotemporal resolution on the PWV and WSS quantification.

RESULTS

4D flow measurements were conducted in an acquisition time of only 35 min. Increased peak flow and peak WSS values and lower errors in PWV estimation were observed in the reconstructions with high temporal resolution. Aortic PWV was significantly increased in ApoE-/- mice compared to the control group (1.7 ± 0.2 versus 2.6 ± 0.2 m/s, p < 0.001). Mean WSS magnitude values averaged over the aortic arch were (1.17 ± 0.07) N/m2 in wildtype mice and (1.27 ± 0.10) N/m2 in ApoE-/- mice.

CONCLUSION

The post processing algorithm using the flexible reconstruction framework developed in this study permitted quantification of global PWV and 3D-WSS in a single acquisition. The possibility to assess both parameters in only 35 min will markedly improve the analyses and information content of in vivo measurements.

Ultrashort echo time (UTE) imaging reveals a shift in bound water that is sensitive to sub-clinical tendinopathy in older adults

OBJECTIVE

Use ultrashort echo time (UTE) magnetic resonance imaging to quantify bound water components of asymptomatic older Achilles tendons and investigate the relationship between UTE findings and imaging assessment of sub-clinical tendinopathy.

MATERIALS AND METHODS

Thirteen young (age 25 ± 4.8) and thirteen older (age 67 ± 4.7) adults were tested. A UTE sequence was used to quantify the transverse relaxation times of bound ([Formula: see text]) and free ([Formula: see text]) water and the bound water fraction (F_s) in the Achilles tendon. Anatomical images were collected and graded by a musculoskeletal radiologist to identify signs of sub-clinical tendinopathy. Two-sample t tests were used to compare [Formula: see text], [Formula: see text], and F_s between age groups and between adults with and without sub-clinical tendinopathy.

RESULTS

Older tendons exhibited a 60% longer [Formula: see text] (p = 0.004), similar [Formula: see text] (p = 0.86), and 5% smaller F_s (p = 0.048) than young tendons. Seven older adult tendons exhibited tendon thickening and increased signal intensity indicative of sub-clinical tendinopathy. This subset of tendons exhibited a 7% smaller bound water fraction (p = 0.02) and significantly longer [Formula: see text] (p < 0.001) than the normal tendons from young and older adults.

CONCLUSION

Older adult tendons exhibited unique UTE signatures that are consistent with disruption of the collagen fiber network and changes in macromolecular content. UTE imaging metrics were sensitive to early indicators of tissue degeneration identified on anatomical images and hence could provide a quantitative biomarker by which to track changes in tissue health resulting from injury, disease, and treatment.

Towards accelerated quantitative sodium MRI at 7 T in the skeletal muscle: Comparison of anisotropic acquisition- and compressed sensing techniques

PURPOSE

To compare three anisotropic acquisition schemes and three compressed sensing (CS) approaches for accelerated tissue sodium concentration (TSC) quantification using ²³Na MRI at 7 T.

MATERIALS AND METHODS

Three anisotropic 3D-radial acquisition sequences were evaluated using simulations, phantom- and in vivo TSC measurements: An anisotropic density-adapted 3D-radial sequence (3DPR-C), a 3D acquisition-weighted density-adapted stack-of-stars sampling scheme (SOS) and a SOS approach with golden-ratio rotation (SOS-GR). Eight healthy volunteers were examined at a 7 Tesla MRI system. TSC measurements of the calf were conducted with a nominal spatial resolution of Δx = (3.0 × 3.0 × 15.0) mm³ and a field of view of (156.0 × 156.0 × 240.0) mm³ for multiple undersampling factors (USF). Three CS reconstructions were evaluated: Total variation CS (TV-CS), 3D dictionary-learning compressed sensing (3D-DLCS) and TV-CS with a block matching prior (TV-BL-CS). Results of the simulations and measurements were compared to a simulated ground truth (GT) or a fully sampled reference measurement (FS), respectively. The deviation of the mean TSC evaluated in multiple ROI (mE_GT/FS) and the normalized root-mean-squared error (NRMSE) for simulations were evaluated for CS and NUFFT reconstructions.

RESULTS

In simulations, the SOS-GR yielded the lowest NRMSE and mE_GT (< 4%) with NUFFT for an acquisition time (TA) of less than 2 min. CS further improved the results. In simulations and measurements, the best TSC quantification results were obtained with 3D-DLCS and SOS-GR (lowest NRMSE, mE_GT < 2.6% in simulations, mE_GT < 10.7% for phantom measurements and mE_FS < 6% in vivo) with an USF = 4.1 (TA < 2 min). TV-CS showed no or only slight improvements to NUFFT. The results of TV-BL-CS were similar to 3D-DLCS.

DISCUSSION

The TA for TSC measurements could be reduced to less than 2 min by using adapted sequences such as SOS-GR and CS reconstruction approaches such as 3D-DLCS or TV-BL-CS, while the quantitative accuracy stays comparable to a fully sampled NUFFT reconstruction (approx. 8 min TA). In future, the lower TA could improve clinical applicability of TSC measurements.

3D variable-density SPARKLING trajectories for high-resolution T2*-weighted magnetic resonance imaging

We have recently proposed a new optimization algorithm called SPARKLING (Spreading Projection Algorithm for Rapid K-space sampLING) to design efficient compressive sampling patterns for magnetic resonance imaging (MRI). This method has a few advantages over conventional non-Cartesian trajectories such as radial lines or spirals: i) it allows to sample the k-space along any arbitrary density while the other two are restricted to radial densities and ii) it optimizes the gradient waveforms for a given readout time. Here, we introduce an extension of the SPARKLING method for 3D imaging by considering both stacks-of-SPARKLING and fully 3D SPARKLING trajectories. Our method allowed to achieve an isotropic resolution of 600 μm in just 45 seconds for T2∗-weighted ex vivo brain imaging at 7 Tesla over a field-of-view of 200 × 200 × 140 mm³ . Preliminary in vivo human brain data shows that a stack-of-SPARKLING is less subject to off-resonance artifacts than a stack-of-spirals.

Variable anisotropic FOV for 3D radial imaging with spiral phyllotaxis (VASP)

PURPOSE

To develop a new 3D radial trajectory based on the natural spiral phyllotaxis (SP), with variable anisotropic FOV.

THEORY & METHODS

A 3D radial trajectory based on the SP with favorable interleaving properties for cardiac imaging has been proposed by Piccini et al (Magn Reson Med. 2011;66:1049-1056), which supports a FOV with a fixed anisotropy. However, a fixed anisotropy can be inefficient when sampling objects with different anisotropic dimensions. We extend Larson's 3D radial method to provide variable anisotropic FOV for spiral phyllotaxis (VASP). Simulations were performed to measure distance between successive projections, analyze point spread functions, and compare aliasing artifacts for both VASP and conventional SP. VASP was fully implemented on a whole-body clinical MR scanner. Phantom and in vivo cardiac images were acquired at 1.5 tesla.

RESULTS

Simulations, phantom, and in vivo experiments confirmed that VASP can achieve variable anisotropic FOV while maintaining the favorable interleaving properties of SP. For an anisotropic FOV with 100:100:35 ratio, VASP required ~65% fewer radial projections than the conventional SP to satisfy Nyquist criteria. Alternatively, when the same number of radial projections were used as in conventional SP, VASP produced fewer aliasing artifacts for anisotropic objects within the excited imaging volumes.

CONCLUSION

We have developed a new method (VASP), which enables variable anisotropic FOV for 3D radial trajectory with SP. For anisotropic objects within the excited imaging volumes, VASP can reduce scan times and/or reduce aliasing artifacts.

Longitudinal evaluation of demyelinated lesions in a multiple sclerosis model using ultrashort echo time magnetization transfer (UTE-MT) imaging

Alterations in myelin integrity are involved in many neurological disorders and demyelinating diseases, such as multiple sclerosis (MS). Although magnetic resonance imaging (MRI) is the gold standard method to diagnose and monitor MS patients, clinically available MRI protocols show limited specificity for myelin detection, notably in cerebral grey matter areas. Ultrashort echo time (UTE) MRI has shown great promise for direct imaging of lipids and myelin sheaths, and thus holds potential to improve lesion detection. In this study, we used a sequence combining magnetization transfer (MT) with UTE ("UTE-MT", TE ​= ​76 ​μs) and with short TE ("STE-MT", TE ​= ​3000 ​μs) to evaluate spatial and temporal changes in brain myelin content in the cuprizone mouse model for MS on a clinical 7 ​T scanner. During demyelination, UTE-MT ratio (UTE-MTR) and STE-MT ratio (STE-MTR) values were significantly decreased in most white matter and grey matter regions. However, only UTE-MTR detected cortical changes. After remyelination in subcortical and cortical areas, UTE-MTR values remained lower than baseline values, indicating that UTE-MT, but not STE-MT, imaging detected long-lasting changes following a demyelinating event. Next, we evaluated the potential correlations between imaging values and underlying histopathological markers. The strongest correlation was observed between UTE-MTR and percent coverage of myelin basic protein (MBP) immunostaining (r² ​= ​0.71). A significant, although lower, correlation was observed between STE-MTR and MBP (r² ​= ​0.48), and no correlation was found between UTE-MTR or STE-MTR and gliosis immunostaining. Interestingly, correlations varied across brain substructures. Altogether, our results demonstrate that UTE-MTR values significantly correlate with myelin content as measured by histopathology, not only in white matter, but also in subcortical and cortical grey matter regions in the cuprizone mouse model for MS. Readily implemented on a clinical 7 ​T system, this approach thus holds great potential for detecting demyelinating/remyelinating events in both white and grey matter areas in humans. When applied to patients with neurological disorders, including MS patient populations, UTE-MT methods may improve the non-invasive longitudinal monitoring of brain lesions, not only during disease progression but also in response to next generation remyelinating therapies.

Fast self-navigated wall shear stress measurements in the murine aortic arch using radial 4D-phase contrast cardiovascular magnetic resonance at 17.6 T

PURPOSE

4D flow cardiovascular magnetic resonance (CMR) and the assessment of wall shear stress (WSS) are non-invasive tools to study cardiovascular risks in vivo. Major limitations of conventional triggered methods are the long measurement times needed for high-resolution data sets and the necessity of stable electrocardiographic (ECG) triggering. In this work an ECG-free retrospectively synchronized method is presented that enables accelerated high-resolution measurements of 4D flow and WSS in the aortic arch of mice.

METHODS

4D flow and WSS were measured in the aortic arch of 12-week-old wildtype C57BL/6 J mice (n = 7) with a radial 4D-phase-contrast (PC)-CMR sequence, which was validated in a flow phantom. Cardiac and respiratory motion signals were extracted from the radial CMR signal and were used for the reconstruction of 4D-flow data. Rigid motion correction and a first order B0 correction was used to improve the robustness of magnitude and velocity data. The aortic lumen was segmented semi-automatically. Temporally averaged and time-resolved WSS and oscillatory shear index (OSI) were calculated from the spatial velocity gradients at the lumen surface at 14 locations along the aortic arch. Reproducibility was tested in 3 animals and the influence of subsampling was investigated.

RESULTS

Volume flow, cross-sectional areas, WSS and the OSI were determined in a measurement time of only 32 min. Longitudinal and circumferential WSS and radial stress were assessed at 14 analysis planes along the aortic arch. The average longitudinal, circumferential and radial stress values were 1.52 ± 0.29 N/m2, 0.28 ± 0.24 N/m2 and - 0.21 ± 0.19 N/m2, respectively. Good reproducibility of WSS values was observed.

CONCLUSION

This work presents a robust measurement of 4D flow and WSS in mice without the need of ECG trigger signals. The retrospective approach provides fast flow quantification within 35 min and a flexible reconstruction framework.

Short-T2 MRI: Principles and recent advances

Among current modalities of biomedical and diagnostic imaging, MRI stands out by virtue of its versatile contrast obtained without ionizing radiation. However, in various cases, e.g., water protons in tissues such as bone, tendon, and lung, MRI performance is limited by the rapid decay of resonance signals associated with short transverse relaxation times T₂ or T₂*. Efforts to address this shortcoming have led to a variety of specialized short-T₂ techniques. Recent progress in this field expands the choice of methods and prompts fresh considerations with regard to instrumentation, data acquisition, and signal processing. In this review, the current status of short-T₂ MRI is surveyed. In an attempt to structure the growing range of techniques, the presentation highlights overarching concepts and basic methodological options. The most frequently used approaches are described in detail, including acquisition strategies, image reconstruction, hardware requirements, means of introducing contrast, sources of artifacts, limitations, and applications.

Implementation of the FLORET UTE sequence for lung imaging

PURPOSE

Magnetic resonance imaging of lungs is inherently challenging, but it has become more common with the use of UTE sequences and their relative insensitivity to motion. Spiral UTE sequences have been touted recently as having greater k-space sampling efficiencies than radial UTE, but few are designed for the shorter T2 * of the lung. In this study, FLORET (Fermat looped, orthogonally encoded trajectories), a recently developed spiral 3D-UTE sequence designed for the short T2 * species, was implemented in human lungs for the first time and the images were compared with traditional radial UTE images.

METHODS

The FLORET sequence was implemented with parameters optimized for lung imaging on healthy and diseased (cystic fibrosis) subjects. On healthy subjects, radial UTE images (3D-radial and 2D-radial with phase encoding) were acquired for comparison to FLORET. Various metrics including SNR, vasculature contrast, diaphragm sharpness, and parenchymal density ratios were acquired and compared among the separate UTE sequences.

RESULTS

The FLORET sequence performed similarly to traditional radial UTE methods with a much shorter total scan time for fully sampled images (FLORET: 1 minute 55 seconds, 3D-radial: 3 minutes 25 seconds, 2D-radial with phase encoding: 7 minutes 22 seconds). Additionally, the FLORET image obtained on the cystic fibrosis subject resulted in the observation of cystic fibrosis lung pathology similar or superior to that of the other UTE-MRI techniques.

CONCLUSION

The FLORET sequence allows for faster acquisition of high diagnostic-quality lung images and its short T2 * components without sacrificing SNR, image quality, or tissue/disease quantification.

Real-time cardiovascular MR with spatio-temporal artifact suppression using deep learning-proof of concept in congenital heart disease

Purpose

Real‐time assessment of ventricular volumes requires high acceleration factors. Residual convolutional neural networks (CNN) have shown potential for removing artifacts caused by data undersampling. In this study, we investigated the ability of CNNs to reconstruct highly accelerated radial real‐time data in patients with congenital heart disease (CHD).

Methods

A 3D (2D plus time) CNN architecture was developed and trained using synthetic training data created from previously acquired breath hold cine images from 250 CHD patients. The trained CNN was then used to reconstruct actual real‐time, tiny golden angle (tGA) radial SSFP data (13 × undersampled) acquired in 10 new patients with CHD. The same real‐time data was also reconstructed with compressed sensing (CS) to compare image quality and reconstruction time. Ventricular volume measurements made using both the CNN and CS reconstructed images were compared to reference standard breath hold data.

Results

It was feasible to train a CNN to remove artifact from highly undersampled radial real‐time data. The overall reconstruction time with the CNN (including creation of aliased images) was shown to be >5 × faster than the CS reconstruction. In addition, the image quality and accuracy of biventricular volumes measured from the CNN reconstructed images were superior to the CS reconstructions.

Conclusion

This article has demonstrated the potential for the use of a CNN for reconstruction of real‐time radial data within the clinical setting. Clinical measures of ventricular volumes using real‐time data with CNN reconstruction are not statistically significantly different from gold‐standard, cardiac‐gated, breath‐hold techniques.

MRI gradient-echo phase contrast of the brain at ultra-short TE with off-resonance saturation

Larmor-frequency shift or image phase measured by gradient-echo sequences has provided a new source of MRI contrast. This contrast is being used to study both the structure and function of the brain. So far, phase images of the brain have been largely obtained at long echo times as maximum phase signal-to-noise ratio (SNR) is achieved at TE = T2* (∼40 ms at 3T). The structures of the brain, however, are compartmentalized and complex with a wide range of signal relaxation times. At such long TE, the short-T2 components are largely attenuated and contribute minimally to phase contrast. The purpose of this study was to determine whether proton gradient-echo images of the brain exhibit phase contrast at ultra-short TE (UTE). Our data showed that UTE images acquired at 7 T without off-resonance saturation do not contain significant phase contrast between gray and white matter. However, UTE images of the brain can attain strong phase contrast even at a nominal TE of 106 μs by using off-resonance RF saturation pulses, which provide direct saturation of ultra-short-T2 components and indirect saturation of longer-T2 components via magnetization transfer. In addition, phase contrast between gray and white matter acquired at UTE with off-resonance saturation is reversed compared to that of the long-T2 signals acquired at long TEs. This finding opens up a potential new way to manipulate image phase contrast of the brain. By accessing short and ultra-short-T2 species, MRI phase images may further improve the characterization of tissue microstructure in the brain.

In vivo characterization of brain ultrashort-T2 components

PURPOSE

Recent nuclear magnetic resonance and MRI studies have measured a fast-relaxing signal component with T2∗<1 ms in white matter and myelin extracts. In ex vivo studies, evidence suggests that a large fraction of this component directly arises from bound protons in the myelin phospholipid membranes. Based on these results, this ultrashort-T₂ component in nervous tissue is a new potential imaging biomarker of myelination, which plays a critical role in neuronal signal conduction across the brain and loss or degradation of myelin is a key feature of many neurological disorders. The goal of this work was to characterize the relaxation times and frequency shifts of ultrashort-T₂ components in the human brain.

METHODS

This required development of an ultrashort echo time relaxometry acquisition strategy and fitting procedure for robust measurements in the presence of ultrashort T2∗ relaxation times and large frequency shifts.

RESULTS

We measured an ultrashort-T₂ component in healthy volunteers with a median T2∗ between 0.5-0.7 ms at 3T and 0.2-0.3 ms at 7T as well as an approximately -3 ppm frequency shift from water.

CONCLUSION

To our knowledge, this is the first time a chemical shift of the ultrashort-T₂ brain component has been measured in vivo. This chemical shift, at around 1.7 ppm, is similar to the primary resonance of most lipids, indicating that much of the ultrashort-T₂ component observed in vivo arises from bound protons in the myelin phospholipid membranes. Magn Reson Med 80:726-735, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Murine pulmonary imaging at 7T: T2* and T1 with anisotropic UTE

PURPOSE

To measure the T2* and T₁ of mouse lung at 7T via anisotropic-resolution radial ultrashort echo-time imaging with ellipsoidal k-space coverage.

METHODS

Ellipsoidal field-of-view was created by expanding uniform spherical k-space coverage. The effects of T2* and ellipsoidal sampling on image resolution were investigated by using point-spread-function analysis and resolution phantoms. Finally, this ellipsoidal sampling approach was used to measure the lung T2* and T₁ of healthy C57BL/6 mice at the increasingly common preclinical field strength of 7T.

RESULTS

Lung parenchyma T2* of 17- to 23-week-old mice at 7T was 0.395 ± 0.033 ms. T₁ of lung and left- and right-heart ventricles was 1452.5 ± 87.0 ms, 1810.5 ± 54.6 ms, and 1602.6 ± 120.9 ms, respectively. Ellipsoidal k-space sampling provides enhanced resolution for a fixed scanning time or provides equivalent (although anisotropic) spatial resolution with reduced scanning times, while simultaneously avoiding fold-in artifacts.

CONCLUSION

Using these techniques, the first T2* and T₁ measures of mouse lung at 7T are reported. Ultrashort echo-time imaging with ellipsoidal k-space coverage significantly increases in-plane resolution without increasing scanning time, or equivalently, decreases scanning time while maintaining equivalent in-plane resolution. Magn Reson Med 79:2254-2264, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Fast volumetric imaging of bound and pore water in cortical bone using three-dimensional ultrashort-TE (UTE) and inversion recovery UTE sequences

We report the three-dimensional ultrashort-TE (3D UTE) and adiabatic inversion recovery UTE (IR-UTE) sequences employing a radial trajectory with conical view ordering for bi-component T2 * analysis of bound water (T2 *(BW) ) and pore water (T2 *(PW) ) in cortical bone. An interleaved dual-echo 3D UTE acquisition scheme was developed for fast bi-component analysis of bound and pore water in cortical bone. A 3D IR-UTE acquisition scheme employing multiple spokes per IR was developed for bound water imaging. Two-dimensional UTE (2D UTE) and IR-UTE sequences were employed for comparison. The sequences were applied to bovine bone samples (n = 6) and volunteers (n = 6) using a 3-T scanner. Bi-component fitting of 3D UTE images of bovine samples showed a mean T2 *(BW) of 0.26 ± 0.04 ms and T2 *(PW) of 4.16 ± 0.35 ms, with fractions of 21.5 ± 3.6% and 78.5 ± 3.6%, respectively. The 3D IR-UTE signal showed a single-component decay with a mean T2 *(BW) of 0.29 ± 0.05 ms, suggesting selective imaging of bound water. Similar results were achieved with the 2D UTE and IR-UTE sequences. Bi-component fitting of 3D UTE images of the tibial midshafts of healthy volunteers showed a mean T2 *(BW) of 0.32 ± 0.08 ms and T2 *(PW) of 5.78 ± 1.24 ms, with fractions of 34.2 ± 7.4% and 65.8 ± 7.4%, respectively. Single-component fitting of 3D IR-UTE images showed a mean T2 *(BW) of 0.35 ± 0.09 ms. The 3D UTE and 3D IR-UTE techniques allow fast volumetric mapping of bound and pore water in cortical bone. Copyright © 2016 John Wiley & Sons, Ltd.

Free-breathing volumetric fat/water separation by combining radial sampling, compressed sensing, and parallel imaging

PURPOSE

Conventional fat/water separation techniques require that patients hold breath during abdominal acquisitions, which often fails and limits the achievable spatial resolution and anatomic coverage. This work presents a novel approach for free-breathing volumetric fat/water separation.

METHODS

Multiecho data are acquired using a motion-robust radial stack-of-stars three-dimensional GRE sequence with bipolar readout. To obtain fat/water maps, a model-based reconstruction is used that accounts for the off-resonant blurring of fat and integrates both compressed sensing and parallel imaging. The approach additionally enables generation of respiration-resolved fat/water maps by detecting motion from k-space data and reconstructing different respiration states. Furthermore, an extension is described for dynamic contrast-enhanced fat-water-separated measurements.

RESULTS

Uniform and robust fat/water separation is demonstrated in several clinical applications, including free-breathing noncontrast abdominal examination of adults and a pediatric subject with both motion-averaged and motion-resolved reconstructions, as well as in a noncontrast breast exam. Furthermore, dynamic contrast-enhanced fat/water imaging with high temporal resolution is demonstrated in the abdomen and breast.

CONCLUSION

The described framework provides a viable approach for motion-robust fat/water separation and promises particular value for clinical applications that are currently limited by the breath-holding capacity or cooperation of patients. Magn Reson Med 78:565-576, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Fast 3D ultrashort echo-time spiral projection imaging using golden-angle: A flexible protocol for in vivo mouse imaging at high magnetic field

PURPOSE

To develop a fast three-dimensional (3D) k-space encoding method based on spiral projection imaging (SPI) with an interleaved golden-angle approach and to validate this novel sequence on small animal models.

METHODS

A disk-like trajectory, in which each disk contained spirals, was developed. The 3D encoding was performed by tilting the disks with a golden angle. The sharpness was first calculated at different T2* values. Then, the sharpness was measured on phantom using variable undersampling ratios. Finally, the sampling method was validated by whole brain time-of-flight angiography and ultrasmall superparamagnetic iron oxide (USPIO) enhanced free-breathing liver angiography on mouse.

RESULTS

The in vitro results demonstrated the robustness of the method for short T2* and high undersampling ratios. In vivo experiments showed the ability to properly detect small vessels in the brain with an acquisition time shorter than 1 min. Free-breathing mice liver angiography showed the insensitivity of this protocol toward motions and flow artifacts, and enabled the visualization of liver motion during breathing.

CONCLUSIONS

The method implemented here allowed fast 3D k-space sampling with a high undersampling ratio. Combining the advantages of center-out spirals with the flexibility of the golden angle approach could have major implications for real-time imaging. Magn Reson Med 77:1831-1840, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Selective in vivo bone imaging with long-T2 suppressed PETRA MRI

PURPOSE

To design and evaluate an optimized PETRA (point-wise encoding time reduction with radial acquisition) sequence with long-T₂ suppression at 3 Tesla.

METHODS

An adiabatic inversion recovery-based scheme was used to null the long-T₂ signal. To minimize scan time, the signal was sampled multiple times after each inversion with variable excitation flip angles designed to yield constant short-T₂ signal amplitude. The excitation pulses were phase-modulated, allowing for increased flip angle and higher signal-to-noise ratio (SNR). A fast, noniterative image reconstruction algorithm was designed to minimize image artifacts due to nonuniform excitation profile.

RESULTS

Phase-modulated pulse excitation, along with the noniterative reconstruction algorithm, allows the use of larger radiofrequency pulse flip angles, resulting in effective suppression of long-T₂ protons and improved image SNR without causing image artifacts. Midtibia images representative of collagen-bound water yielded SNR of 15 at 1-mm isotropic resolution in 6.5 minutes with a standard extremity coil. Further, the technology is shown to be suited for generating multi-angle projection images of bone akin to X-ray images displaying subtle anatomic detail.

CONCLUSION

Optimized long-T₂ suppressed PETRA allows imaging of bone matrix water unencumbered by long-T₂ soft tissue and pore water protons, opening up new possibilities for anatomic bone imaging at isotropic resolution and quantification in clinically practical scan times. Magn Reson Med 77:989-997, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Combined outer volume suppression and T2 preparation sequence for coronary angiography

PURPOSE

To develop a magnetization preparation sequence for simultaneous outer volume suppression (OVS) and T2 weighting in whole-heart coronary magnetic resonance angiography.

METHODS

A combined OVS and T2 preparation sequence (OVS-T2 Prep) was designed with a nonselective adiabatic 90° tipdown pulse, two adiabatic 180° refocusing pulses, and a 2D spiral -90° tipup pulse. The OVS-T2 Prep preserves the magnetization inside an elliptic cylinder with T2 weighting, while saturating the magnetization outside the cylinder. Its performance was tested on phantoms and on 13 normal subjects with coronary magnetic resonance angiography using 3D cones trajectories.

RESULTS

Phantom studies showed expected T2 -dependent signal amplitude in the spatial passband and suppressed signal in the spatial stopband. In vivo studies with full-field-of-view cones yielded a passband-to-stopband signal ratio of 3.18 ± 0.77 and blood-myocardium contrast-to-noise ratio enhancement by a factor of 1.43 ± 0.20 (P < 0.001). In vivo studies with reduced-field-of-view cones showed that OVS-T2 Prep well suppressed the aliasing artifacts, as supported by significantly reduced signal in the regions with no tissues compared to the images acquired without preparation (P < 0.0001).

CONCLUSION

OVS-T2 Prep is a compact sequence that can accelerate coronary magnetic resonance angiography by suppressing signals from tissues surrounding the heart while simultaneously enhancing the blood-myocardium contrast.

Quantifying temperature-dependent T1 changes in cortical bone using ultrashort echo-time MRI

PURPOSE

To demonstrate the feasibility of using ultrashort echo-time MRI to quantify T1 changes in cortical bone due to heating.

METHODS

Variable flip-angle T1 mapping combined with 3D ultrashort echo-time imaging was used to measure T1 in cortical bone. A calibration experiment was performed to detect T1 changes with temperature in ex vivo cortical bone samples from a bovine femur. Ultrasound heating experiments were performed using an interstitial applicator in ex vivo bovine femur specimens, and heat-induced T1 changes were quantified.

RESULTS

The calibration experiment demonstrated that T1 increases with temperature in cortical bone. We observed a linear relationship between temperature and T1 with a linear coefficient between 0.67 and 0.84 ms/°C over a range of 25-70°C. The ultrasound heating experiments showed increased T1 changes in the heated regions, and the relationship between the temperature changes and T1 changes was similar to that of the calibration.

CONCLUSION

We demonstrated a temperature dependence of T1 in ex vivo cortical bone using a variable flip-angle ultrashort echo-time T1 mapping method.

Anisotropic field-of-view support for golden angle radial imaging

PURPOSE

To provide anisotropic field-of-view (FOV) support for golden angle radial imaging.

THEORY AND METHODS

In radial imaging, uniform spoke density leads to a circular FOV, which is excessive for objects with anisotropic dimensions. Larson et al previously showed that the angular k-space spoke density can be determined by the desired anisotropic FOV. We show that conventional golden angle sampling can be deployed in an angle-normalized space and transformed back to k-space such that the desired nonuniform spoke density is preserved for arbitrary temporal window length. Elliptical FOVs were used to illustrate this generalized mapping approach. Point-spread-function and spoke density analysis was performed. Phantom and in vivo cardiac images were acquired.

RESULTS

Simulations, phantom, and in vivo experiments confirmed that the proposed method is able to achieve anisotropic FOV while still maintaining the benefits of golden angle sampling. This approach requires 50% less spokes for elliptical FOV with major-to-minor-axis ratio of 1:0.3, when compared with isotropic FOV with the same undersampling factor.

CONCLUSION

We demonstrate a simple method for applying golden angle view ordering to anisotropic FOV radial imaging. This can reduce imaging time for objects with anisotropic dimensions while still allowing arbitrary temporal window selection. Magn Reson Med 76:229-236, 2016. © 2015 Wiley Periodicals, Inc.

High-resolution variable-density 3D cones coronary MRA

PURPOSE

To improve the spatial/temporal resolution of whole-heart coronary MR angiography by developing a variable-density (VD) 3D cones acquisition suitable for image reconstruction with parallel imaging and compressed sensing techniques.

METHODS

A VD 3D cones trajectory design incorporates both radial and spiral trajectory undersampling techniques to achieve higher resolution. This design is used to generate a VD 3D cones trajectory with 0.8 mm/66 ms isotropic spatial/temporal resolution, using a similar number of readouts as our previous fully sampled cones trajectory (1.2 mm/100 ms). Scans of volunteers and patients are performed to evaluate the performance of the VD trajectory, using non-Cartesian L1 -ESPIRiT for high-resolution image reconstruction.

RESULTS

With gridding reconstruction, the high-resolution scans experience an expected drop in signal-to-noise and contrast-to-noise ratios, but with L1 -ESPIRiT, the apparent noise is substantially reduced. Compared with 1.2 mm images, in each volunteer, the L1 -ESPIRiT 0.8 mm images exhibit higher vessel sharpness values in the right and left anterior descending arteries.

CONCLUSION

Coronary MR angiography with isotropic submillimeter spatial resolution and high temporal resolution can be performed with VD 3D cones to improve the depiction of coronary arteries.

In Vivo Quantitative MR Imaging of Bound and Pore Water in Cortical Bone

PURPOSE

To translate and evaluate an in vivo magnetic resonance (MR) imaging protocol for quantitative mapping of collagen-bound and pore water concentrations in cortical bone that involves relaxation-selective ultrashort echo time (UTE) methods.

MATERIALS AND METHODS

All HIPAA-compliant studies were performed with institutional review board approval and written informed consent. UTE imaging sequences were implemented on a clinical 3.0-T MR imaging unit and were used for in vivo imaging of bound and pore water in cortical bone. Images of the lower leg and wrist were acquired in five volunteers each (lower leg: two men and three women aged 24, 24, 49, 30, and 26 years; wrist: two men and three women aged 31, 23, 25, 24, and 26 years) to generate bound and pore water concentration maps of the tibia and radius. Each volunteer was imaged three times, and the standard error of the measurements at the region-of-interest (ROI) level was computed as the standard deviation across studies, pooled across volunteers and ROIs.

RESULTS

Quantitative bound and pore water maps in the tibia and radius, acquired in 8-14 minutes, had per-voxel signal-to-noise ratios of 18 (bound water) and 14 (pore water) and inter-study standard errors of approximately 2 mol (1)H per liter of bone at the ROI level.

CONCLUSION

The results of this study demonstrate the feasibility of quantitatively mapping bound and pore water in vivo in human cortical bone with practical human MR imaging constraints.

Non-contrast-enhanced peripheral angiography using a sliding interleaved cylinder acquisition

PURPOSE

To develop a new sequence for non-contrast-enhanced peripheral angiography using a sliding interleaved cylinder (SLINCYL) acquisition.

METHODS

A venous saturation pulse was incorporated into a three-dimensional magnetization-prepared balanced steady-state free precession sequence for non-contrast-enhanced peripheral angiography to improve artery-vein contrast. The SLINCYL acquisition, which consists of a series of overlapped thin slabs for volumetric coverage similar to the original sliding interleaved ky (SLINKY) acquisition, was used to evenly distribute the venous-suppression effects over the field of view. In addition, the thin-slab-scan nature of SLINCYL and the centric-ordered sampling geometry of its readout trajectory were exploited to implement efficient fluid-suppression and parallel imaging schemes. The sequence was tested in healthy subjects and a patient.

RESULTS

Compared to a multiple overlapped thin slab acquisition, both SLINKY and SLINCYL suppressed the venetian blind artifacts and provided similar artery-vein contrast. However, SLINCYL achieved this with shorter scan times and less noticeable artifacts from k-space amplitude modulation than SLINKY. The fluid-suppression and parallel imaging schemes were also validated. A patient study using the SLINCYL-based sequence well identified stenoses at the superficial femoral arteries, which were also confirmed with digital subtraction angiography.

CONCLUSION

Non-contrast-enhanced angiography using SLINCYL can provide angiograms with improved artery-vein contrast in the lower extremities.

Depiction of achilles tendon microstructure in vivo using high-resolution 3-dimensional ultrashort echo-time magnetic resonance imaging at 7 T

OBJECTIVES

The objective of this study was to demonstrate the feasibility of depicting the internal structure of the Achilles tendon in vivo using high-resolution 3-dimensional ultrashort echo-time (UTE) magnetic resonance imaging at 7 T.

MATERIALS AND METHODS

For our UTE imaging, a minimum-phase radiofrequency pulse and an anisotropic field-of-view 3-dimensional radial acquisition were used to minimize the echo time and scan time. A fat saturation pulse was applied every 8 spoke acquisitions to reduce blurring and chemical shift artifacts from fat and to improve the dynamic range of the tendon signal. Five healthy volunteers and 1 patient were scanned with an isotropic spatial resolution of up to 0.6 mm. Fat-suppressed UTE images were qualitatively evaluated and compared with non-fat-suppressed UTE images and longer echo-time images.

RESULTS

High-resolution UTE imaging was able to visualize the microstructure of the Achilles tendon. Fat suppression substantially improved the depiction of the internal structure. The UTE images revealed a fascicular pattern in the Achilles tendon and fibrocartilage at the tendon insertion. In a patient who had tendon elongation surgery after birth, there was a clear depiction of disrupted tendon structure.

CONCLUSIONS

High-resolution fat-suppressed 3-dimensional UTE imaging at 7 T allows for the evaluation of the Achilles tendon microstructure in vivo.

Measuring venous blood oxygenation in fetal brain using susceptibility-weighted imaging

PURPOSE

To evaluate fetal cerebral venous blood oxygenation, Yv, using principles of MR susceptometry.

MATERIALS AND METHODS

A cohort of 19 pregnant subjects, with a mean gestational age of 31.6 ± 4.7 weeks were imaged using a modified susceptibility-weighted imaging (SWI) sequence. Data quality was first assessed for feasibility of oxygen saturation measurement, and data from five subjects (mean ± std gestational age of 33.7 ± 3.6 weeks) were then chosen for further quantitative analysis. SWI phase in the superior sagittal sinus was used to evaluate oxygen saturation using the principles of MR susceptometry. Systematic error in the measured Y(v) values was studied through simulations.

RESULTS

Simulations showed that the systematic error in Yv depended upon the assumed angle of the vessel, θ, relative to the main magnetic field and the error in that vessel angle δθ. For the typical vessel angle of θ = 30° encountered in the fetal data analyzed, a δθ as large as ±20° led to an absolute error, δYv, of less than 11%. The measured mean oxygen saturation across the five fetuses was 66% ± 9.4%. This average cerebral venous blood oxygenation value is in close agreement with values in the published literature.

CONCLUSION

We have reported the first in vivo measurement of human fetal cerebral venous oxygen saturation using MRI.

Three-dimensional through-time radial GRAPPA for renal MR angiography

PURPOSE

To achieve high temporal and spatial resolution for contrast-enhanced time-resolved MR angiography exams (trMRAs), fast imaging techniques such as non-Cartesian parallel imaging must be used. In this study, the three-dimensional (3D) through-time radial generalized autocalibrating partially parallel acquisition (GRAPPA) method is used to reconstruct highly accelerated stack-of-stars data for time-resolved renal MRAs.

MATERIALS AND METHODS

Through-time radial GRAPPA has been recently introduced as a method for non-Cartesian GRAPPA weight calibration, and a similar concept can also be used in 3D acquisitions. By combining different sources of calibration information, acquisition time can be reduced. Here, different GRAPPA weight calibration schemes are explored in simulation, and the results are applied to reconstruct undersampled stack-of-stars data.

RESULTS

Simulations demonstrate that an accurate and efficient approach to 3D calibration is to combine a small number of central partitions with as many temporal repetitions as exam time permits. These findings were used to reconstruct renal trMRA data with an in-plane acceleration factor as high as 12.6 with respect to the Nyquist sampling criterion, where the lowest root mean squared error value of 16.4% was achieved when using a calibration scheme with 8 partitions, 16 repetitions, and a 4 projection × 8 read point segment size.

CONCLUSION

3D through-time radial GRAPPA can be used to successfully reconstruct highly accelerated non-Cartesian data. By using in-plane radial undersampling, a trMRA can be acquired with a temporal footprint less than 4s/frame with a spatial resolution of approximately 1.5 mm × 1.5 mm × 3 mm.

Validation of a golden angle radial sequence (GOLD) for abdominal T1 mapping during free breathing: demonstrating clinical feasibility for quantifying gastric secretion and emptying

PURPOSE

To validate a magnetic resonance imaging sequence suitable for quantitative assessment of acid suppression by a proton pump inhibitor (PPI) on gastric secretion and emptying in clinical practice.

METHODS

A golden angle radial sequence (GOLD) was validated in a series of in vitro and in vivo experiments and clinical feasibility was shown in two studies. The impact of free breathing and image plane orientation on T1 values was evaluated in a controlled in vivo experiment. The free-breathing GOLD sequence was compared against a standard breath-hold gradient echo sequence for gastric half emptying time in 23 subjects during a gastric emptying study. Pilot data from five subjects assessed the sensitivity of the GOLD sequence to detect changes in acid secretion volume produced by PPI treatment.

RESULTS

The coronal free-breathing GOLD sequence and the axial breath-hold standard gradient echo sequence showed good agreement of the gastric half emptying time (6 ± 3 min, P = 0.053). The GOLD sequence demonstrated sensitivity to reduction of gastric secretion volumes induced by PPI treatment (55 ± 5 mL, P < 0.001).

CONCLUSION

The GOLD sequence allowed for free breathing, multislice, combined imaging and T1 mapping of the stomach content. GOLD presents a promising multipurpose, noninvasive imaging tool for monitoring gastric function in clinical studies.

Measurement techniques for magnetic resonance imaging of fast relaxing nuclei

In this review article, techniques for sodium ((23)Na) magnetic resonance imaging (MRI) are presented. These techniques can also be used to image other nuclei with short relaxation times (e.g., (39)K, (35)Cl, (17)O). Twisted projection imaging, density-adapted 3D projection reconstruction, and 3D cones are preferred because of uniform k-space sampling and ultra-short echo times. Sampling density weighted apodization can be applied if intrinsic filtering is desired. This approach leads to an increased signal-to-noise ratio compared to postfiltered acquisition in cases of short readout durations relative to T 2 (*) relaxation time. Different MR approaches for anisotropic resolution are presented, which are important for imaging of thin structures such as myocardium, cartilage, and skin. The third part of this review article describes different methods to put more weighting either on the intracellular or the extracellular sodium signal by means of contrast agents, relaxation-weighted imaging, or multiple-quantum filtering.

High-resolution 3D radial bSSFP with IDEAL

Radial trajectories facilitate high-resolution balanced steady state free precession (bSSFP) because the efficient gradients provide more time to extend the trajectory in k-space. A number of radial bSSFP methods that support fat-water separation have been developed; however, most of these methods require an environment with limited B0 inhomogeneity. In this work, high-resolution bSSFP with fat-water separation is achieved in more challenging B0 environments by combining a 3D radial trajectory with the IDEAL chemical species separation method. A method to maintain very high resolution within the timing constraints of bSSFP and IDEAL is described using a dual-pass pulse sequence. The sampling of a unique set of radial lines at each echo time is investigated as a means to circumvent the longer scan time that IDEAL incurs as a multiecho acquisition. The manifestation of undersampling artifacts in this trajectory and their effect on chemical species separation are investigated in comparison to the case in which each echo samples the same set of radial lines. This new bSSFP method achieves 0.63 mm isotropic resolution in a 5-min scan and is demonstrated in difficult in vivo imaging environments, including the breast and a knee with ACL reconstruction hardware at 1.5 T.

3D imaging using magnetic resonance tomosynthesis (MRT) technique

PURPOSE

To introduce an alternative approach to three-dimensional (3D) magnetic resonance (MR) imaging using a method that is similar to x-ray tomosynthesis.

METHODS

Variable angle tilted-projection images are acquired using a multiple-oblique view (MOV) pulse sequence. Reconstruction is performed using three methods similar to that of x-ray tomosynthesis, which generate a set of tomographic images with multiple 2D projection images. The reconstruction algorithm is further modified to reformat to the practical imaging situations of MR. The procedure is therefore termed magnetic resonance tomosynthesis (MRT). To analyze the characteristics of MRT, simulations are performed. Phantom and in vivo experiments were done to suggest potential applications.

RESULTS

Simulation results show anisotropic features that are structurally dependent in terms of resolution. Partial blurrings along slice direction were observed. In phantom and in vivo experiments, the reconstruction performance is particularly noticeable in the low SNR case where improved images with lower noise are obtained. Reformatted reconstruction using thinner slice thickness and∕or extended field-of-view can increase spatial resolution partially and alleviate slice profile imperfection.

CONCLUSIONS

Results demonstrate that MRT can generate adequate 3D images using the MOV images. Various reconstruction methods in tomosynthesis were readily adapted, while allowing other tomosynthesis reconstruction algorithms to be incorporated. A reformatted reconstruction process was incorporated for applications relevant to MR imaging.

Free-breathing multiphase whole-heart coronary MR angiography using image-based navigators and three-dimensional cones imaging

Noninvasive visualization of the coronary arteries in vivo is one of the most important goals in cardiovascular imaging. Compared to other paradigms for coronary MR angiography, a free-breathing three-dimensional whole-heart iso-resolution approach simplifies prescription effort, requires less patient cooperation, reduces overall exam time, and supports retrospective reformats at arbitrary planes. However, this approach requires a long continuous acquisition and must account for respiratory and cardiac motion throughout the scan. In this work, a new free-breathing coronary MR angiography technique that reduces scan time and improves robustness to motion is developed. Data acquisition is accomplished using a three-dimensional cones non-Cartesian trajectory, which can reduce the number of readouts 3-fold or more compared to conventional three-dimensional Cartesian encoding and provides greater robustness to motion/flow effects. To further enhance robustness to motion, two-dimensional navigator images are acquired to directly track respiration-induced displacement of the heart and enable retrospective compensation of all acquired data (none discarded) for image reconstruction. In addition, multiple cardiac phases are imaged to support retrospective selection of the best phase(s) for visualizing each coronary segment. Experimental results demonstrate that whole-heart coronary angiograms can be obtained rapidly and robustly with this proposed technique.

A novel method for quantitative monitoring of transplanted islets of langerhans by positive contrast magnetic resonance imaging

The Automatic Quantitative Ultrashort Echo Time imaging (AQUTE) protocol for serial MRI allows quantitative in vivo monitoring of iron labeled pancreatic islets of Langerhans transplanted into the liver, quantifying graft implantation and persistence in a rodent model. Rats (n = 14), transplanted with iron oxide loaded cells (0-4000 islet equivalents, IEQ), were imaged using a 3D radial ultrashort echo time difference technique (dUTE) on a Siemens MAGNETOM 3T clinical scanner up to 5 months postsurgery. In vivo 3D dUTE images gave positive contrast from labeled cells, suppressing liver signal and small vessels, allowing automatic quantification. Position of labeled islet clusters was consistent over time and quantification of hyperintense pixels correlated with the number of injected IEQs (R² = 0.898, p < 0.0001), and showed persistence over time (5 months posttransplantation). Automatic quantification was superior to standard imaging and manual counting methods, due to the uniform suppressed background and high contrast, resulting in significant timesavings, reproducibility and ease of quantification. Three-dimensional coverage of the whole liver in the absence of cardiac/respiratory artifact provided further improvement over conventional imaging. This imaging protocol reliably quantifies transplanted islet mass and has high translational potential to clinical studies of transplanted pancreatic islets.

Dual half-echo phase correction for implementation of 3D radial SSFP at 3.0 T

Fat/water separation methods such as fluctuating equilibrium magnetic resonance and linear combination steady-state free precession have not yet been successfully implemented at 3.0 T due to extreme limitations on the time available for spatial encoding with the increase in magnetic field strength. We present a method to utilize a three-dimensional radial sequence combined with linear combination steady-state free precession at 3.0 T to take advantage of the increased signal levels over 1.5 T and demonstrate high spatial resolution compared to Cartesian techniques. We exploit information from the two half-echoes within each pulse repetition time to correct the accumulated phase on a point-by-point basis, thereby fully aligning the phase of both half-echoes. The correction provides reduced sensitivity to static field (B(0)) inhomogeneity and robust fat/water separation. Resultant images in the knee joint demonstrate the necessity of such a correction, as well as the increased isotropic spatial resolution attainable at 3.0 T. Results of a clinical study comparing this sequence to conventional joint imaging sequences are included.

Elliptical field-of-view PROPELLER imaging

Traditionally two-dimensional scans are designed to support an isotropic field-of-view (iFOV). When imaging elongated objects, significant savings in scan time can potentially be achieved by supporting an elliptical field-of-view (eFOV). This work presents an empirical closed-form solution to adapt the PROPELLER trajectory for an eFOV. The proposed solution is built on the geometry of the PROPELLER trajectory permitting the scan prescription and data reconstruction to remain largely similar to standard PROPELLER. The achieved FOV is experimentally validated by the point spread function (PSF) of a phantom scan. The details of potential savings in scan time and the signal-to-noise ratio (SNR) performance in comparison to iFOV scans for both phantom and in-vivo images are also described.