Description: Purpose To evaluate the feasibility of an improved motion and flow robust methodology for imaging the pulmonary vasculature using non-contrast-enhanced, free-breathing, golden-angle radial MRI. Methods Healthy volunteers ( n  = 10, age 46 ± 11 years, 50% female) and patients ( n  = 2, ages 27 and 84, both female) were imaged at 1.5 T using a Cartesian and golden-angle radial 2D balanced SSFP pulse sequence. The acquisitions were made under free breathing without contrast agent enhancement. The radial acquisitions were reconstructed at 3 temporal footprints. All series were scored from 1 to 5 for perceived diagnostic quality, artifact level, and vessel sharpness in multiple anatomical locations. In addition, vessel sharpness and blood-to-blood clot contrast were measured. Results Quantitative measurements showed higher vessel sharpness for golden-angle radial ( n  = 76, 0.79 ± 0.11 versus 0.71 ± 0.16, p  〈 .05). Blood-to-blood clot contrast was found to be 23% higher in golden-angle radial in the 2 patients. At comparable temporal footprints, golden-angle radial was scored higher for diagnostic quality (mean ± SD, 2.3 ± 0.7 versus 2.2 ± 0.6, p  〈 .01) and vessel sharpness (2.2 ± 0.8 versus 2.1 ± 0.5, p  〈 .01), whereas the artifact level did not differ (3.0 ± 0.9 versus 3.0 ± 1.0, p  = .80). The ability to retrospectively choose a temporal resolution and perform sliding-window reconstructions was demonstrated in patients. Conclusion In pulmonary artery imaging, the motion and flow robustness of a radial trajectory does both improve image quality over Cartesian trajectory in healthy volunteers, and allows for flexible selection of temporal footprints and the ability to perform real-time sliding window reconstructions, which could potentially provide further diagnostic insight.

Description: Purpose To build and evaluate a small-footprint, lightweight, high-performance 3T MRI scanner for advanced brain imaging with image quality that is equal to or better than conventional whole-body clinical 3T MRI scanners, while achieving substantial reductions in installation costs. Methods A conduction-cooled magnet was developed that uses less than 12 liters of liquid helium in a gas-charged sealed system, and standard NbTi wire, and weighs approximately 2000 kg. A 42-cm inner-diameter gradient coil with asymmetric transverse axes was developed to provide patient access for head and extremity exams, while minimizing magnet-gradient interactions that adversely affect image quality. The gradient coil was designed to achieve simultaneous operation of 80-mT/m peak gradient amplitude at a slew rate of 700 T/m/s on each gradient axis using readily available 1-MVA gradient drivers. Results In a comparison of anatomical imaging in 16 patients using T 2 -weighted 3D fluid-attenuated inversion recovery (FLAIR) between the compact 3T and whole-body 3T, image quality was assessed as equivalent to or better across several metrics. The ability to fully use a high slew rate of 700 T/m/s simultaneously with 80-mT/m maximum gradient amplitude resulted in improvements in image quality across EPI, DWI, and anatomical imaging of the brain. Conclusions The compact 3T MRI system has been in continuous operation at the Mayo Clinic since March 2016. To date, over 200 patient studies have been completed, including 96 comparison studies with a clinical 3T whole-body MRI. The increased gradient performance has reliably resulted in consistently improved image quality.

Description: Purpose The ultimate intrinsic signal-to-noise ratio (UISNR) represents an upper bound for the achievable SNR of any receive coil. To reach this threshold a complete basis set of equivalent surface currents is required. This study systematically investigated to what extent either loop- or dipole-like current patterns are able to reach the UISNR threshold in a realistic human head model between 1.5 T and 11.7 T. Based on this analysis, we derived guidelines for coil designers to choose the best array element at a given field strength. Moreover, we present ideal current patterns yielding the UISNR in a realistic body model. Methods We distributed generic current patterns on a cylindrical and helmet-shaped surface around a realistic human head model. We excited electromagnetic fields in the human head by using eigenfunctions of the spherical and cylindrical Helmholtz operator. The electromagnetic field problem was solved by a fast volume integral equation solver. Results At 7 T and above, adding curl-free current patterns to divergence-free current patterns substantially increased the SNR in the human head (locally 〉20%). This was true for the helmet-shaped and the cylindrical surface. On the cylindrical surface, dipole-like current patterns had high SNR performance in central regions at ultra-high field strength. The UISNR increased superlinearly with B0 in most parts of the cerebrum but only sublinearly in the periphery of the human head. Conclusion The combination of loop and dipole elements could enhance the SNR performance in the human head at ultra-high field strength.

Description: Purpose To characterize and suppress stripe artifact associated with velocity-selective (VS) magnetization for unenhanced MRA. Methods Extended phase graph formalism was used to show that the stripe artifact contains multiples of the fundamental frequency that is determined by the area of unipolar VS gradient. Four VS preparation pulses whose excitation profiles are spatially shifted by quarter the fundamental period of the stripes, were applied alternately. For further suppression of the artifact, k-space data at k z  = 0 were averaged over the 4 VS preparations. The proposed schemes were tested in a chicken breast phantom and healthy human subjects. Results When the standard VS preparation scheme was used, stripe artifact was shown in all the reconstructed images and appeared as artifactual peaks in k-space that corresponded to the first and second order harmonics of the fundamental frequency. Alternate application of the 4 phase-shifted VS preparation pulses suppressed the stripes, but not completely, as evidenced by residual erroneous peaks in k-space. After the k-space averaging, the stripe artifact was nearly eliminated. Conclusion Stripe artifact in VS-MRA consists of multiples of the fundamental frequency and can be effectively suppressed through alternate application of phase-shifted VS preparations along with k-space averaging.

Description: Purpose Deranged metabolism is now recognized as a key causal factor in a variety of heart diseases, and is being studied extensively. However, invasive methods may alter metabolism, and conventional imaging techniques measure tracer uptake but not downstream metabolism. These challenges may be overcome by hyperpolarized MR, a noninvasive technique currently crossing the threshold into human trials. The aim of this study was to image metabolic changes in the heart in response to endogastric glucose bolus and to acute hypertension. Methods Five postprandial pigs were scanned with hyperpolarized [1- 13 C]pyruvate cardiac MR at baseline, after oral glucose bolus, and after infusion of angiotensin-II. Results No effect of glucose bolus was seen using hyperpolarized [1- 13 C]pyruvate MR despite changes in circulating substrates. During angiotensin-II infusion, blood pressure increased 179% ( P  = 0.008) and ejection fraction decreased from 54 ± 2% to 47 ± 6% ( P  = 0.03) The hemodynamic changes were accompanied by increases in the hyperpolarized [1- 13 C]pyruvate MR derived ratios of lactate/alanine (from 0.58 ± 0.13 to 0.78 ± 0.06, P  = 0.03) and bicarbonate/alanine (from 0.55 ± 0.12 to 0.91 ± 0.14, P  = 0.007). Conclusion Glucose loading did not alter cardiac metabolism, but during acute hypertensive stress, cardiac aerobic, carbohydrate metabolism, and pyruvate-lactate exchange was altered. Hyperpolarized MR allows noninvasive evaluation of acute changes in cardiac metabolism. However, hemodynamics must be taken into account when interpreting the results.

Description: Purpose Proton MRSI is a noninvasive modality capable of generating volumetric maps of in vivo tissue metabolism without the need for ionizing radiation or injected contrast agent. Magnetic resonance spectroscopic imaging has been shown to be a viable imaging modality for studying several neuropathologies. However, a key hurdle in the routine clinical adoption of MRSI is the presence of spectral artifacts that can arise from a number of sources, possibly leading to false information. Methods A deep learning model was developed that was capable of identifying and filtering out poor quality spectra. The core of the model used a tiled convolutional neural network that analyzed frequency-domain spectra to detect artifacts. Results When compared with a panel of MRS experts, our convolutional neural network achieved high sensitivity and specificity with an area under the curve of 0.95. A visualization scheme was implemented to better understand how the convolutional neural network made its judgement on single-voxel or multivoxel MRSI, and the convolutional neural network was embedded into a pipeline capable of producing whole-brain spectroscopic MRI volumes in real time. Conclusion The fully automated method for assessment of spectral quality provides a valuable tool to support clinical MRSI or spectroscopic MRI studies for use in fields such as adaptive radiation therapy planning.

Description: Purpose This study demonstrates a DCE-MRI estimate of tumor interstitial fluid pressure (TIFP) and hydraulic conductivity in a rat model of glioblastoma, with validation against an invasive wick-in-needle (WIN) technique. An elevated TIFP is considered a mark of aggressiveness, and a decreased TIFP a predictor of response to therapy. Methods The DCE-MRI studies were conducted in 36 athymic rats (controls and posttreatment animals) with implanted U251 cerebral tumors, and with TIFP measured using a WIN method. Using a model selection paradigm and a novel application of Patlak and Logan plots to DCE-MRI data, the MRI parameters required for estimating TIFP noninvasively were estimated. Two models, a fluid-mechanical model and a multivariate empirical model, were used for estimating TIFP, as verified against WIN-TIFP. Results Using DCE-MRI, the mean estimated hydraulic conductivity (MRI-K) in U251 tumors was (2.3 ± 3.1) × 10 −5 (mm 2 /mmHg-s) in control studies. Significant positive correlations were found between WIN-TIFP and MRI-TIFP in both mechanical and empirical models. For instance, in the control group of the fluid-mechanical model, MRI-TIFP was a strong predictor of WIN-TIFP (R 2  = 0.76, p  〈 .0001). A similar result was found in the bevacizumab-treated group of the empirical model (R 2  = 0.93, p  = .014). Conclusion This research suggests that MRI dynamic studies contain enough information to noninvasively estimate TIFP in this, and possibly other, tumor models, and thus might be used to assess tumor aggressiveness and response to therapy.

Description: Purpose Avoid formation of staircase artifacts in nonlinear diffusion-based MR image reconstruction without compromising computational speed. Methods Whereas second-order diffusion encourages the evolution of pixel neighborhood with uniform intensities, fourth-order diffusion considers smooth region to be not necessarily a uniform intensity region but also a planar region. Therefore, a controlled application of fourth-order diffusivity function is used to encourage second-order diffusion to reconstruct the smooth regions of the image as a plane rather than a group of blocks, while not being strong enough to introduce the undesirable speckle effect. Results Proposed method is compared with second- and fourth-order nonlinear diffusion reconstruction, total variation (TV), total generalized variation, and higher degree TV using in vivo data sets for different undersampling levels with application to dictionary learning-based reconstruction. It is observed that the proposed technique preserves sharp boundaries in the image while preventing the formation of staircase artifacts in the regions of smoothly varying pixel intensities. It also shows reduced error measures compared with second-order nonlinear diffusion reconstruction or TV and converges faster than TV-based methods. Conclusion Because nonlinear diffusion is known to be an effective alternative to TV for edge-preserving reconstruction, the crucial aspect of staircase artifact removal is addressed. Reconstruction is found to be stable for the experimentally determined range of fourth-order regularization parameter, and therefore not does not introduce a parameter search. Hence, the computational simplicity of second-order diffusion is retained.

Description: Purpose To investigate the feasibility of measuring the subtle disruption of blood-brain barrier (BBB) using DCE-MRI with a scan duration shorter than 10 min. Methods The extended Patlak-model (EPM) was introduced to include the effect of plasma flow ( F p ) in the estimation of vascular permeability–surface area product ( PS ). Numerical simulation studies were carried out to investigate how the reduction in scan time affects the accuracy in estimating contrast kinetic parameters. DCE-MRI studies of the rat brain were conducted with Fisher rats to confirm the results from the simulation. Intracranial F98 glioblastoma models were used to assess areas with different levels of permeability. In the normal brain tissues, the Patlak model (PM) and EPM were compared, whereas the 2-compartment-exchange-model (TCM) and EPM were assessed in the peri-tumor and the tumor regions. Results The simulation study results demonstrated that scan time reduction could lead to larger bias in PS estimated by PM (〉2000%) than by EPM (〈47%), especially when F p is low. When F p was high as in the gray matter, the bias in PM- PS (〉900%) were larger than that in EPM- PS (〈42%). The animal study also showed similar results, where the PM parameters were more sensitive to the scan duration than the EPM parameters. It was also demonstrated that, in the peri-tumor region, the EPM parameters showed less change by scan duration than the TCM parameters. Conclusion The results of this study suggest that EPM can be used to measure PS with a scan duration of 10 min or less.

Description: Purpose 3D time-resolved (4D) phase contrast MRI can be used to study muscle contraction. However, 3D coverage with sufficient spatiotemporal resolution can only be achieved by interleaved acquisitions during many repetitions of the motion task, resulting in long scan times. The aim of this study was to develop a compressed sensing accelerated 4D phase contrast MRI technique for quantification of velocities and strain rate of the muscles in the lower leg during active plantarflexion/dorsiflexion. Methods Nine healthy volunteers were scanned during active dorsiflexion/plantarflexion task. For each volunteer, we acquired a reference scan, as well as 4 different accelerated scans (k-space undersampling factors: 3.14X, 4.09X, 4.89X, and 6.41X) obtained using Cartesian Poisson disk undersampling schemes. The data was reconstructed using a compressed sensing pipeline. For each scan, velocity and strain rate values were quantified in the gastrocnemius lateralis, gastrocnemius medialis, tibialis anterior, and soleus. Results No significant differences in velocity values were observed as a function acceleration factor in the investigated muscles. The strain rate calculation resulted in one positive (s + ) and one negative (s − ) eigenvalue, whereas the third eigenvalue (s 3 ) was consistently 0 for all the acquisitions. No significant differences were observed for the strain rate eigenvalues as a function of acceleration factor. Conclusions Data undersampling combined with compressed sensing reconstruction allowed obtainment of time-resolved phase contrast acquisitions with 3D coverage and quantitative information comparable to the reference scan. The 3D sensitivity of the method can help in understanding the connection between muscle architecture and muscle function in future studies.