A multipurpose, adolescent idiopathic scoliosis-specific, short MRI protocol: A feasibility study in volunteers.

INTRODUCTION: Visualization of scoliosis typically requires ionizing radiation (radiography and CT) to visualize bony anatomy. MRI is often additionally performed to screen for neural axis abnormalities. We propose a 14-minutes radiation-free scoliosis-specific MRI protocol, which combines MRI and MRI-based synthetic CT images to visualize soft and osseous structures in one examination. We assess the ability of the protocol to visualize landmarks needed to detect 3D patho-anatomical changes, screen for neural axis abnormalities, and perform surgical planning and navigation. METHODS: 18 adult volunteers were scanned on 1.5 T MR-scanner using 3D T2-weighted and synthetic CT sequences. A predefined checklist of relevant landmarks was used for the parameter assessment by three readers. Parameters included Cobb angles, rotation, torsion, segmental height, area and centroids of Nucleus Pulposus and Intervertebral Disc. Precision, reliability and agreement between the readers measurements were evaluated. RESULTS: 91 % of Likert-based questions scored ≥ 4, indicating moderate to high confidence. Precision of 3D dot positioning was 1.0 mm. Precision of angle measurement was 0.6° (ICC 0.98). Precision of vertebral and IVD height measurements was 0.4 mm (ICC 0.99). Precision of area measurement for NP was 8 mm2 (ICC 0.55) and for IVD 18 mm2 (ICC 0.62) for IVD. Precision of centroid measurement for NP was 1.3 mm (ICC 0.88-0.92) and for IVD 1.1 mm (ICC 0.88-91). CONCLUSIONS: The proposed MRI protocol with synthetic CT reconstructions, has high precision, reliability and agreement between the readers for multiple scoliosis-specific measurements. It can be used to study scoliosis etiopathogenesis and to assess 3D spinal morphology.

Indirect 1H-[13C] MRS of the human brain at 7 T using a 13C-birdcage coil and eight transmit-receive 1H-dipole antennas with a 32-channel 1H-receive array.

The neuronal tricarboxylic acid and glutamate/glutamine (Glu/Gln) cycles play important roles in brain function. These processes can be measured in vivo using dynamic 1H-[13C] MRS during administration of 13C-labeled glucose. Proton-observed carbon-edited (POCE) MRS enhances the signal-to-noise ratio (SNR) compared with direct 13C-MRS. Ultra-high field further boosts the SNR and increases spectral dispersion; however, even at 7 T, Glu and Gln 1H-resonances may overlap. Further gain can be obtained with selective POCE (selPOCE). Our aim was to create a setup for indirect dynamic 1H-[13C] MRS in the human brain at 7 T. A home-built non-shielded transmit-receive 13C-birdcage head coil with eight transmit-receive 1H-dipole antennas was used together with a 32-channel 1H-receive array. Electromagnetic simulations were carried out to ensure that acquisitions remained within local and global head SAR limits. POCE-MRS was performed using slice-selective excitation with semi-localization by adiabatic selective refocusing (sLASER) and stimulated echo acquisition mode (STEAM) localization, and selPOCE-MRS using STEAM. Sequences were tested in a phantom containing non-enriched Glu and Gln, and in three healthy volunteers during uniformly labeled 13C-glucose infusions. In one subject the voxel position was alternated between bi-frontal and bi-occipital placement within one session. [4-13C]Glu-H4 and [4-13C]Gln-H4 signals could be separately detected using both STEAM-POCE and STEAM-selPOCE in the phantom. In vivo, [4,5-13C]Glx could be detected using both sLASER-POCE and STEAM-POCE, with similar sensitivities, but [4,5-13C]Glu and [4,5-13C]Gln signals could not be completely resolved. STEAM-POCE was alternately performed bi-frontal and bi-occipital within a single session without repositioning of the subject, yielding similar results. With STEAM-selPOCE, [4,5-13C]Glu and [4,5-13C]Gln could be clearly separated. We have shown that with our setup indirect dynamic 1H-[13C] MRS at 7 T is feasible in different locations in the brain within one session, and by using STEAM-selPOCE it is possible to separate Glu from Gln in vivo while obtaining high quality spectra.

Automatic generation of subject-specific finite element models of the spine from magnetic resonance images.

The generation of subject-specific finite element models of the spine is generally a time-consuming process based on computed tomography (CT) images, where scanning exposes subjects to harmful radiation. In this study, a method is presented for the automatic generation of spine finite element models using images from a single magnetic resonance (MR) sequence. The thoracic and lumbar spine of eight adult volunteers was imaged using a 3D multi-echo-gradient-echo sagittal MR sequence. A deep-learning method was used to generate synthetic CT images from the MR images. A pre-trained deep-learning network was used for the automatic segmentation of vertebrae from the synthetic CT images. Another deep-learning network was trained for the automatic segmentation of intervertebral discs from the MR images. The automatic segmentations were validated against manual segmentations for two subjects, one with scoliosis, and another with a spine implant. A template mesh of the spine was registered to the segmentations in three steps using a Bayesian coherent point drift algorithm. First, rigid registration was applied on the complete spine. Second, non-rigid registration was used for the individual discs and vertebrae. Third, the complete spine was non-rigidly registered to the individually registered discs and vertebrae. Comparison of the automatic and manual segmentations led to dice-scores of 0.93-0.96 for all vertebrae and discs. The lowest dice-score was in the disc at the height of the implant where artifacts led to under-segmentation. The mean distance between the morphed meshes and the segmentations was below 1 mm. In conclusion, the presented method can be used to automatically generate accurate subject-specific spine models.

In vivo phosphorus magnetic resonance spectroscopic imaging of the whole human liver at 7 T using a phosphorus whole-body transmit coil and 16-channel receive array: Repeatability and effects of principal component analysis-based denoising.

Quantitative three-dimensional (3D) imaging of phosphorus (31 P) metabolites is potentially a promising technique with which to assess the progression of liver disease and monitor therapy response. However, 31 P magnetic resonance spectroscopy has a low sensitivity and commonly used 31 P surface coils do not provide full coverage of the liver. This study aimed to overcome these limitations by using a 31 P whole-body transmit coil in combination with a 16-channel 31 P receive array at 7 T. Using this setup, we determined the repeatability of whole-liver 31 P magnetic resonance spectroscopic imaging (31 P MRSI) in healthy subjects and assessed the effects of principal component analysis (PCA)-based denoising on the repeatability parameters. In addition, spatial variations of 31 P metabolites within the liver were analyzed. 3D 31 P MRSI data of the liver were acquired with a nominal voxel size of 20 mm isotropic in 10 healthy volunteers twice on the same day. Data were reconstructed without denoising, and with PCA-based denoising before or after channel combination. From the test-retest data, repeatability parameters for metabolite level quantification were determined for 12 31 P metabolite signals. On average, 31 P MR spectra from 100 ± 25 voxels in the liver were analyzed. Only voxels with contamination from skeletal muscle or the gall bladder were excluded and no voxels were discarded based on (low) signal-to-noise ratio (SNR). Repeatability for most quantified 31 P metabolite levels in the liver was good to excellent, with an intrasubject variability below 10%. PCA-based denoising increased the SNR ~ 3-fold, but did not improve the repeatability for mean liver 31 P metabolite quantification with the fitting constraints used. Significant spatial heterogeneity of various 31 P metabolite levels within the liver was observed, with marked differences for the phosphomonoester and phosphodiester metabolites between the left and right lobe. In conclusion, using a 31 P whole-body transmit coil in combination with a 16-channel 31 P receive array at 7 T allowed 31 P MRSI acquisitions with full liver coverage and good to excellent repeatability.

A plug-and-play, lightweight, single-axis gradient insert design for increasing spatiotemporal resolution in echo planar imaging-based brain imaging.

The goal of this study was to introduce and evaluate the performance of a lightweight, high-performance, single-axis (z-axis) gradient insert design primarily intended for high-resolution functional magnetic resonance imaging, and aimed at providing both ease of use and a boost in spatiotemporal resolution. The optimal winding positions of the coil were obtained using a genetic algorithm with a cost function that balanced gradient performance (minimum 0.30 mT/m/A) and field linearity (≥16 cm linear region). These parameters were verified using field distribution measurements by B0 -mapping. The correction of geometrical distortions was performed using theoretical field distribution of the coil. Simulations and measurements were performed to investigate the echo planar imaging echo-spacing reduction due to the improved gradient performance. The resulting coil featured a 16-cm linear region, a weight of 45 kg, an installation time of 15 min, and a maximum gradient strength and slew rate of 200 mT/m and 1300 T/m/s, respectively, when paired with a commercially available gradient amplifier (940 V/630 A). The field distribution measurements matched the theoretically expected field. By utilizing the theoretical field distribution, geometrical distortions were corrected to within 6% of the whole-body gradient reference image in the target region. Compared with a whole-body gradient set, a maximum reduction in echo-spacing of a factor of 2.3 was found, translating to a 344 µs echo-spacing, for a field of view of 192 mm, a receiver bandwidth of 920 kHz and a gradient amplitude of 112 mT/m. We present a lightweight, single-axis gradient insert design that can provide high gradient performance and an increase in spatiotemporal resolution with correctable geometrical distortions while also offering a short installation time of less than 15 min and minimal system modifications.

PCA denoising and Wiener deconvolution of 31 P 3D CSI data to enhance effective SNR and improve point spread function.

PURPOSE: This study evaluates the performance of 2 processing methods, that is, principal component analysis-based denoising and Wiener deconvolution, to enhance the quality of phosphorus 3D chemical shift imaging data. METHODS: Principal component analysis-based denoising increases the SNR while maintaining spectral information. Wiener deconvolution reduces the FWHM of the voxel point spread function, which is increased by Hamming filtering or Hamming-weighted acquisition. The proposed methods are evaluated using simulated and in vivo 3D phosphorus chemical shift imaging data by 1) visual inspection of the spatial signal distribution; 2) SNR calculation of the PCr peak; and 3) fitting of metabolite basis functions. RESULTS: With the optimal order of processing steps, we show that the effective SNR of in vivo phosphorus 3D chemical shift imaging data can be increased. In simulations, we show we can preserve phosphorus-containing metabolite peaks that had an SNR < 1 before denoising. Furthermore, using Wiener deconvolution, we were able to reduce the FWHM of the voxel point spread function with only partially reintroducing Gibb-ringing artifacts while maintaining the SNR. After data processing, fitting of the phosphorus-containing metabolite signals improved. CONCLUSION: In this study, we have shown that principal component analysis-based denoising in combination with regularized Wiener deconvolution allows increasing the effective spectral SNR of in vivo phosphorus 3D chemical shift imaging data, with reduction of the FWHM of the voxel point spread function. Processing increased the effective SNR by at least threefold compared to Hamming weighted acquired data and minimized voxel bleeding. With these methods, fitting of metabolite amplitudes became more robust with decreased fitting residuals.

Field drift correction of proton resonance frequency shift temperature mapping with multichannel fast alternating nonselective free induction decay readouts.

PURPOSE: To demonstrate that proton resonance frequency shift MR thermometry (PRFS-MRT) acquisition with nonselective free induction decay (FID), combined with coil sensitivity profiles, allows spatially resolved B0 drift-corrected thermometry. METHODS: Phantom experiments were performed at 1.5T and 3T. Acquisition of PRFS-MRT and FID were performed during MR-guided high-intensity focused ultrasound heating. The phase of the FIDs was used to estimate the change in angular frequency δωdrift per coil element. Two correction methods were investigated: (1) using the average δωdrift over all coil elements (0th-order) and (2) using coil sensitivity profiles for spatially resolved correction. Optical probes were used for independent temperature verification. In-vivo feasibility of the methods was evaluated in the leg of 1 healthy volunteer at 1.5T. RESULTS: In 30 minutes, B0 drift led to an apparent temperature change of up to -18°C and -98°C at 1.5T and 3T, respectively. In the sonicated area, both corrections had a median error of 0.19°C at 1.5T and -0.54°C at 3T. At 1.5T, the measured median error with respect to the optical probe was -1.28°C with the 0th-order correction and improved to 0.43°C with the spatially resolved correction. In vivo, without correction the spatiotemporal median of the apparent temperature was at -4.3°C and interquartile range (IQR) of 9.31°C. The 0th-order correction had a median of 0.75°C and IQR of 0.96°C. The spatially resolved method had the lowest median at 0.33°C and IQR of 0.80°C. CONCLUSION: FID phase information from individual receive coil elements allows spatially resolved B0 drift correction in PRFS-based MRT.

Shortening of apparent transverse relaxation time of inorganic phosphate as a breast cancer biomarker.

Phosphorus MRS offers a non-invasive tool for monitoring cell energy and phospholipid metabolism and can be of additional value in diagnosing cancer and monitoring cancer therapy. In this study, we determined the transverse relaxation times of a number of phosphorous metabolites in a group of breast cancer patients by adiabatic multi-echo spectroscopic imaging at 7 T. The transverse relaxation times of phosphoethanolamine, phosphocholine, inorganic phosphate (Pi ), glycerophosphocholine and glycerophosphatidylcholine were 184 ± 8 ms, 203 ± 17 ms, 87 ± 8 ms, 240 ± 56 ms and 20 ± 10 ms, respectively. The transverse relaxation time of Pi in breast cancer tissue was less than half that of healthy fibroglandular tissue. This effect is most likely caused by an up-regulation of glycolysis in breast cancer tissue that leads to interaction of Pi with the GAPDH enzyme, which forms part of the reversible pathway of exchange of Pi with gamma-adenosine tri-phosphate, thus shortening its apparent transverse relaxation time. As healthy breast tissue shows very little glycolytic activity, the apparent T2 shortening of Pi due to malignant transformation could possibly be used as a biomarker for cancer.

Improved fat suppression of the breast using discretized frequency shimming.

PURPOSE: Robust fat suppression is essential in bilateral breast MRI at 7 Tesla. The lack of good fat suppression can result in errors when calculating the enhancement curve from dynamic contrast-enhanced acquisitions. In this work we propose discretized frequency shimming to improve the quality of fat suppression by exploiting the capabilities of a parallel transmit setup. METHODS: With discretized frequency shimming, different resonance frequencies are used for the different transmitters of the radiofrequency coil. This can be particularly useful for frequency-selective radiofrequency pulses, such as water excitation. The potential of this method has been investigated by simulations on seven acquired B0 maps. Prospective experiments on phantoms and four healthy volunteers were performed to illustrate the performance of discretized frequency shimming. RESULTS: Simulations show a decrease in off-resonance excitation for all subjects when discretized frequency shimming is used. This decrease is highly subject-dependent, but can result in up to 30% better fat suppression, as observed in phantom and volunteer experiments. CONCLUSIONS: Discretized frequency shimming improves fat suppression in bilateral breast MRI at 7 Tesla. Magn Reson Med 79:593-599, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A comparison of navigators, snap-shot field monitoring, and probe-based field model training for correcting B0 -induced artifacts in T2*-weighted images at 7 T.

PURPOSE: To compare methods for estimating B0 maps used in retrospective correction of high-resolution anatomical images at ultra-high field strength. The B0 maps were obtained using three methods: (1) 1D navigators and coil sensitivities, (2) field probe (FP) data and a low-order spherical harmonics model, and (3) FP data and a training-based model. METHODS: Data from nine subjects were acquired while they performed activities inducing B0 field fluctuations. Estimated B0 fields were compared with reference data, and the reductions of artifacts were compared in corrected T2* images. RESULTS: Reduction of sum-of-squares difference relative to a reference image was evaluated, and Method 1 yielded the largest artifact reduction: 27 ± 15%, 20 ± 18% (mean ± 1 standard deviation) for deep breathing and combined deep breathing and hand motion activities. Method 3 performed almost as well (24 ± 18%, 15 ± 17%), provided that adequate training data were used, and Method 2 gave a similar result (21 ± 16%, 19 ± 17%). CONCLUSION: This study confirms that all of the investigated methods can be used in retrospective image correction. In terms of image quality, Method 1 had a small advantage, whereas the FP-based methods measured the B0 field slightly more accurately. The specific strengths and weaknesses of FPs and navigators should therefore be considered when determining which B0 -estimation method to use. Magn Reson Med 78:1373-1382, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Erratum to: Dynamic contrast-enhanced breast MRI at 7T and 3T: an intra-individual comparison study.

[This corrects the article DOI: 10.1186/s40064-015-1654-7.].

Dynamic contrast-enhanced breast MRI at 7T and 3T: an intra-individual comparison study.

The aim of this study is to compare the current state of lesion identification, the BI-RADS classification and the contrast-enhancement behavior at 7T and 3T breast MRI in the same patient group. Twenty-seven patients with thirty suspicious lesions were selected for this prospective study and underwent both 7T and 3T MRI. All examinations were rated by two radiologists (R1 and R2) independently on image quality, lesion identification and BI-RADS classification. We assessed sensitivity, specificity, NPV and PPV, observer agreement, lesion sizes, and contrast-enhancement-to-noise ratios (CENRs) of mass lesions. Fifteen of seventeen histopathological proven malignant lesions were detected at both field strengths. Image quality of the dynamic series was good at 7T, and excellent at 3T (P = 0.001 for R1 and P = 0.88 for R2). R1 found higher rates of specificity, NPV and PPV at 7T when compared to 3T, while R2 found the same results for sensitivity, specificity, NPV and PPV for both field strengths. The observers showed excellent agreement for BI-RADS categories at 7T (κ = 0.86) and 3T (κ = 0.93). CENRs were higher at 7T (P = 0.015). Lesion sizes were bigger at 7T according to R2 (P = 0.039). Our comparison study shows that 7T MRI allows BI-RADS conform analysis. Technical improvements, such as acquisition of T2w sequences and adjustment of B1+ field inhomogeneity, are still necessary to allow clinical use of 7T breast MRI.

Fat suppression techniques for obtaining high resolution dynamic contrast enhanced bilateral breast MR images at 7T.

OBJECTIVES: To compare water selective excitation (WSE) and Dixon fat suppression in the context of high-resolution dynamic contrast enhanced MRI of the breast at 7T. METHODS: Ten healthy volunteers and one patient with a malignant breast lesion were scanned at 7T. The MRI protocol contained 3D T1-weighted gradient echo images obtained with both WSE fat suppression, multi echo Dixon fat suppression, and without fat suppression. Images were acquired at a (0.8mm)(3) or (0.7mm)(3) isotropic resolution with equal field of view and optimized such to obtain a maximal SNR. Image quality was scored qualitatively on overall image quality, sharpness of anatomical details, presence of artifacts, inhomogeneous fat suppression and the presence of water-fat shift. A quantitative scoring was obtained from the signal to noise ratio and contrast to noise ratio. RESULTS: WSE scored significantly better in terms of overall image quality and the absence of artifacts. No significant difference in contrast to noise ratio was found between the two fat suppression methods. CONCLUSION: When maximizing temporal and spatial resolution of high resolution DCE MRI of the breast, water selective excitation provides better image quality than multi echo Dixon at 7T.

Multiparametric MRI With Dynamic Contrast Enhancement, Diffusion-Weighted Imaging, and 31-Phosphorus Spectroscopy at 7 T for Characterization of Breast Cancer.

OBJECTIVES: To describe and to correlate tumor characteristics on multiparametric 7 tesla (T) breast magnetic resonance imaging (MRI) with prognostic characteristics from postoperative histopathology in patients with breast cancer. MATERIALS AND METHODS: Institutional review board approval and written informed consent of 15 women (46-70 years) with 17 malignant lesions were obtained. In this prospective study (March 2013 to March 2014), women were preoperatively scanned using dynamic contrast-enhanced MRI, diffusion-weighted imaging, and 31-phosphorus spectroscopy (¹³P-MRS). The value of the protocol was assessed to quantify tumor differentiation and proliferation. Dynamic contrast-enhanced MRI was assessed according to the American College of Radiology Breast Imaging Reporting and Data System-MRI lexicon. Apparent diffusion coefficients (ADCs) were calculated from diffusion-weighted imaging. On ¹³P-MRS, at the location of the tumor, the amount of phosphorus components was obtained in a localized spectrum. In this spectrum, the height of phosphodiester (PDE) and phosphomonoester (PME) peaks was assessed to serve as a measure for metabolic activity, stratifying tumors into a PDE > PME, PDE = PME, or PDE < PME group. Tumor grade and mitotic count from resection specimen were compared with the MRI characteristics using explorative analyses. RESULTS: On dynamic contrast-enhanced MRI, the mean tumor size was 24 mm (range, 6-55 mm). An inverse trend was seen between ADC and tumor grade (P = 0.083), with mean ADC of 867 × 10⁻6 mm²/s for grade 1 (N = 4), 751 × 10⁻6 mm²/s for grade 2 (N = 6), and 659 × 10⁻6 mm²/s for grade 3 (N = 2) tumors. Between P-MR spectra and mitotic count, a relative increase of PME over PDE showed significant association with increasing mitotic counts (P = 0.02); a mean mitotic count of 6 was found in the PDE greater than PME group (N = 7), 8 in the PDE = PME group (N = 1), and 17 in the PDE < PME group (N = 3). CONCLUSIONS: Multiparametric 7 T breast MRI is feasible in clinical setting and shows association between ADC and tumor grade, and between ¹³P-MRS and mitotic count.

Radiofrequency configuration to facilitate bilateral breast (31) P MR spectroscopic imaging and high-resolution MRI at 7 Tesla.

PURPOSE: High-resolution MRI combined with phospholipid detection may improve breast cancer grading. Currently, configurations are optimized for either high-resolution imaging or (31) P spectroscopy. To be able to perform both imaging as well as spectroscopy in a single session, we integrated a (1) H receiver array into a (1) H-(31) P transceiver at 7T. To ensure negligible signal loss due to coupling between elements, we investigated the use of a floating decoupling loop to enable bilateral MRI and (31) P MRS. METHODS: Two quadrature double-tuned radiofrequency coils were designed for bilateral breast MR with active detuning at the (1) H frequency. The two coils were placed adjacent to each other and decoupled for both frequencies with a single resonant floating loop. Sensitivity of the bilateral configuration, facilitating space for a 26-element (1) H receive array, was compared with a transceiver configuration. RESULTS: The floating loop was able to decouple the elements over 20 dB for both frequencies. Enlargement of the elements, to provide space for the receivers, and the addition of detuning electronics altered the (31) P sensitivity by 0.4 dB. CONCLUSION: Dynamic contrast-enhanced scans of 0.7 mm isotropic, diffusion-weighted imaging, and (31) P MR spectroscopic imaging can be acquired at 7T in a single session as demonstrated in a patient with invasive ductal carcinoma.

MRI of the carotid artery at 7 Tesla: quantitative comparison with 3 Tesla.

PURPOSE: To evaluate the 7 Tesla (T) MRI of the carotid arteries, as quantitatively compared with 3T. MATERIALS AND METHODS: The 7T MRI of the carotid arteries was performed in six healthy subjects and in two patients with carotid stenosis. The healthy group was scanned at 3T and at 7T, using current coil setups at both field strengths. T1 and T2 values of the normal carotid vessel wall were assessed at both field strengths. B1 (+) maps and signal to noise ratio (SNR) maps were obtained, as well as T1 weighted images with a resolution as high as 0.4 × 0.4 × 1.5 mm(3) . RESULTS: The T1 of the normal carotid vessel wall was found to be 1227 ± 47 ms at 3T and 1628 ± 130 ms at 7T, while a T2 of 55 ± 11 ms at 3T and 46 ± 4 ms at 7T was found. A twofold average gain in SNR at the carotid arteries was found with 7T. T1 weighted images showed an increased SNR at 7T for all subjects. CONCLUSION: Evaluation between 3T and 7T carotid MRI with optimized setups at both field strengths showed improved SNR at 7T, an increase in vessel wall T1 and a decrease in vessel wall T2.