Unsupervised deep learning with convolutional neural networks for static parallel transmit design: A retrospective study

PURPOSE

To mitigateB 1 + $$ {B}_1^{+} $$ inhomogeneity at 7T for multi-channel transmit arrays using unsupervised deep learning with convolutional neural networks (CNNs).

METHODS

Deep learning parallel transmit (pTx) pulse design has received attention, but such methods have relied on supervised training and did not use CNNs for multi-channelB 1 + $$ {B}_1^{+} $$ maps. In this work, we introduce an alternative approach that facilitates the use of CNNs with multi-channelB 1 + $$ {B}_1^{+} $$ maps while performing unsupervised training. The multi-channelB 1 + $$ {B}_1^{+} $$ maps are concatenated along the spatial dimension to enable shift-equivariant processing amenable to CNNs. Training is performed in an unsupervised manner using a physics-driven loss function that minimizes the discrepancy of the Bloch simulation with the target magnetization, which eliminates the calculation of reference transmit RF weights. The training database comprises 3824 2D sagittal, multi-channelB 1 + $$ {B}_1^{+} $$ maps of the healthy human brain from 143 subjects.B 1 + $$ {B}_1^{+} $$ data were acquired at 7T using an 8Tx/32Rx head coil. The proposed method is compared to the unregularized magnitude least-squares (MLS) solution for the target magnetization in static pTx design.

RESULTS

The proposed method outperformed the unregularized MLS solution for RMS error and coefficient-of-variation and had comparable energy consumption. Additionally, the proposed method did not show local phase singularities leading to distinct holes in the resulting magnetization unlike the unregularized MLS solution.

CONCLUSION

Proposed unsupervised deep learning with CNNs performs better than unregularized MLS in static pTx for speed and robustness.

Reproducibility of 3D pH-weighted chemical exchange saturation transfer contrast in the healthy cervical spinal cord

Spinal cord ischemia and hypoxia can be caused by compression, injury, and vascular alterations. Measuring ischemia and hypoxia directly in the spinal cord noninvasively remains challenging. Ischemia and hypoxia alter tissue pH, providing a physiologic parameter that may be more directly related to tissue viability. Chemical exchange saturation transfer (CEST) is an MRI contrast mechanism that can be made sensitive to pH. More specifically, amine/amide concentration independent detection (AACID) is a recently developed endogenous CEST contrast that has demonstrated sensitivity to intracellular pH at 9.4 T. The goal of this study was to evaluate the reproducibility of AACID CEST measurements at different levels of the healthy cervical spinal cord at 3.0 T incorporating B1 correction. Using a 3.0 T MRI scanner, two 3D CEST scans (saturation pulse train followed by a 3D snapshot gradient-echo readout) were performed on 12 healthy subjects approximately 10 days apart, with the CEST volume centered at the C4 level for all subjects. Scan-rescan reproducibility was evaluated by examining between and within-subject coefficients of variation (CVs) and absolute AACID value differences. The C4 level of the spinal cord demonstrated the lowest within-subject CVs (4.1%-4.3%), between-subject CVs (5.6%-6.3%), and absolute AACID percent difference (5.8-6.1%). The B1 correction scheme significantly improved reproducibility (adjusted p-value = 0.002) compared with the noncorrected data, suggesting that implementing B1 corrections in the spinal cord is beneficial. It was concluded that pH-weighted AACID measurements, incorporating B1-inhomogeneity correction, were reproducible within subjects along the healthy cervical spinal cord and that optimal image quality was achieved at the center of the 3D CEST volume.

Fast online spectral-spatial pulse design for subject-specific fat saturation in cervical spine and foot imaging at 1.5 T

OBJECTIVE

To compensate subject-specific field inhomogeneities and enhance fat pre-saturation with a fast online individual spectral-spatial (SPSP) single-channel pulse design.

METHODS

The RF shape is calculated online using subject-specific field maps and a predefined excitation k-space trajectory. Calculation acceleration options are explored to increase clinical viability. Four optimization configurations are compared to a standard Gaussian spectral selective pre-saturation pulse and to a Dixon acquisition using phantom and volunteer (N = 5) data at 1.5 T with a turbo spin echo (TSE) sequence. Measurements and simulations are conducted across various body parts and image orientations.

RESULTS

Phantom measurements demonstrate up to a 3.5-fold reduction in residual fat signal compared to Gaussian fat saturation. In vivo evaluations show improvements up to sixfold for dorsal subcutaneous fat in sagittal cervical spine acquisitions. The versatility of the tailored trajectory is confirmed through sagittal foot/ankle, coronal, and transversal cervical spine experiments. Additional measurements indicate that excitation field (B1) information can be disregarded at 1.5 T. Acceleration methods reduce computation time to a few seconds.

DISCUSSION

An individual pulse design that primarily compensates for main field (B0) inhomogeneities in fat pre-saturation is successfully implemented within an online "push-button" workflow. Both fat saturation homogeneity and the level of suppression are improved.

Fluid suppression in amide proton transfer-weighted (APTw) CEST imaging: New theoretical insights and clinical benefits

PURPOSE

Amide proton transfer-weighted (APTw) MRI at 3T provides a unique contrast for brain tumor imaging. However, APTw imaging suffers from hyperintensities in liquid compartments such as cystic or necrotic structures and provides a distorted APTw signal intensity. Recently, it has been shown that heuristically motivated fluid suppression can remove such artifacts and significantly improve the readability of APTw imaging.

THEORY AND METHODS

In this work, we show that the fluid suppression can actually be understood by the known concept of spillover dilution, which itself can be derived from the Bloch-McConnell equations in comparison to the heuristic approach. Therefore, we derive a novel post-processing formula that efficiently removes fluid artifact, and explains previous approaches. We demonstrate the utility of this APTw assessment in silico, in vitro, and in vivo in brain tumor patients acquired at MR scanners from different vendors.

RESULTS

Our results show a reduction of the CEST signals from fluid environments while keeping the APTw-CEST signal intensity almost unchanged for semi-solid tissue structures such as the contralateral normal appearing white matter. This further allows us to use the same color bar settings as for conventional APTw imaging.

CONCLUSION

Fluid suppression has considerable value in improving the readability of APTw maps in the neuro-oncological field. In this work, we derive a novel post-processing formula from the underlying Bloch-McConnell equations that efficiently removes fluid artifact, and explains previous approaches which justify the derivation of this metric from a theoretical point of view, to reassure the scientific and medical field about its use.

Dynamic parallel transmit diffusion MRI at 7T

Diffusion MRI (dMRI) is inherently limited by SNR. Scanning at 7 T increases intrinsic SNR but 7 T MRI scans suffer from regions of signal dropout, especially in the temporal lobes and cerebellum. We applied dynamic parallel transmit (pTx) to allow whole-brain 7 T dMRI and compared with circularly polarized (CP) pulses in 6 subjects. Subject-specific 2-spoke dynamic pTx pulses were designed offline for 8 slabs covering the brain. We used vendor-provided B0 and B1+ mapping. Spokes positions were set using the Fourier difference approach, and RF coefficients optimized with a Jacobi-matrix high-flip-angle optimizer. Diffusion data were analyzed with FSL. Comparing whole-brain averages for pTx against CP scans: mean flip angle error improved by 15% for excitation (2-spoke-VERSE 15.7° vs CP 18.4°, P = 0.012) and improved by 14% for refocusing (2-spoke-VERSE 39.7° vs CP 46.2°, P = 0.008). Computed spin-echo signal standard deviation improved by 14% (2-spoke-VERSE 0.185 vs 0.214 CP, P = 0.025). Temporal SNR increased by 5.4% (2-spoke-VERSE 8.47 vs CP 8.04, P = 0.004) especially in the inferior temporal lobes. Diffusion fitting uncertainty decreased by 6.2% for first fibers (2-spoke VERSE 0.0655 vs CP 0.0703, P < 0.001) and 1.3% for second fibers (2-spoke VERSE 0.139 vs CP 0.141, P = 0.01). In conclusion, dynamic parallel transmit improves the uniformity of 7 T diffusion-weighted imaging. In future, less restrictive SAR limits for parallel transmit scans are expected to allow further improvements.

Improved characterization of lenticulostriate arteries using compressed sensing time-of-flight at 7T

OBJECTIVES

To evaluate the feasibility of 0.2-mm isotropic lenticulostriate arteries (LSAs) imaging using compressed sensing time-of-flight (CS TOF) at around 10 min on 7T, and compare the delineation and characterization of LSAs using conventional TOF and CS TOF.

METHODS

Thirty healthy volunteers were examined with CS TOF and conventional TOF at 7T for around 10 min each. CS TOF was optimized to achieve 0.2-mm isotropic LSA imaging. The numbers of LSA stems and branches were counted and compared on a vascular skeleton. The length and distance were measured and compared on the most prominent branch in each hemisphere. Another patient with intracranial artery stenosis was studied to compare LSA delineation in CS TOF and digital subtraction angiography (DSA).

RESULTS

The number of stems visualized with CS TOF was significantly higher than with conventional TOF in both left (p = 0.002, ICC = 0.884) and right (p < 0.001, ICC = 0.938) hemispheres. The number of branches visualized by conventional TOF was significantly lower than that by CS TOF in both left (p < 0.001, ICC = 0.893) and right (p < 0.001, ICC = 0.896) hemispheres. The lengths were statistically higher in CS TOF than in conventional TOF (left: p < 0.001, ICC = 0.868; right: p < 0.001, ICC = 0.876).

CONCLUSIONS

The high-resolution CS TOF improves the delineation and characterization of LSAs over conventional TOF. High-resolution LSA imaging using CS TOF can be a promising tool for clinical research and applications in patients with neurologic diseases.

KEY POINTS

• 0.2-mm isotropic LSA imaging for around 10 min using CS TOF at 7T is feasible. • More stems and branches of LSAs with longer lengths can be delineated with CS TOF than with conventional TOF at the same scan time. • High-resolution CS TOF can be a promising tool for research and applications on LSA.

Snapshot CEST++: Advancing rapid whole-brain APTw-CEST MRI at 3 T

APTw CEST MRI suffers from long preparation times and consequently long acquisition times (~5 min). Recently, a consensus on the preparation module for clinical APTw CEST at 3 T was found in the community, and we present a fast whole-brain APTw CEST MRI sequence following this consensus preparation of pulsed RF irradiation of 2 s duration at 90% RF duty-cycle and a B1,rms of 2 μT. After optimization of the snapshot CEST approach for APTw imaging regarding flip angle, voxel size and frequency offset sampling, we extend it by undersampled GRE acquisition and compressed sensing reconstruction. This allows 2 mm isotropic whole-brain APTw imaging for clinical research at 3 T below 2 min. With this sequence, a fast snapshot APTw imaging method is now available for larger clinical studies of brain tumors.

CEST 2022 - Differences in APT-weighted signal in T1 weighted isointense lesions, black holes and normal-appearing white matter in people with relapsing-remitting multiple sclerosis

PURPOSE

To evaluate amide proton transfer weighted (APTw) signal differences between multiple sclerosis (MS) lesions and contralateral normal-appearing white matter (cNAWM). Cellular changes during the demyelination process were also assessed by comparing APTw signal intensity in T1weighted isointense (ISO) and hypointense (black hole -BH) MS lesions in relation to cNAWM.

METHODS

Twenty-four people with relapsing-remitting MS (pw-RRMS) on stable therapy were recruited. MRI/APTw acquisitions were undertaken on a 3 T MRI scanner. The pre and post-processing, analysis, co-registration with structural MRI maps, and identification of regions of interest (ROIs) were all performed with Olea Sphere 3.0 software. Generalized linear model (GLM) univariate ANOVA was undertaken to test the hypotheses that differences in mean APTw were entered as dependent variables. ROIs were entered as random effect variables, which allowed all data to be included. Regions (lesions and cNAWM) and/or structure (ISO and BH) were the main factor variables. The models also included age, sex, disease duration, EDSS, and ROI volumes as covariates. Receiver operating characteristic (ROC) curve analyses were performed to evaluate the diagnostic performance of these comparisons.

RESULTS

A total of 502 MS lesions manually identified on T2-FLAIR from twenty-four pw-RRMS were subcategorized as 359 ISO and 143 BH with reference to the T1-MPRAGE cerebral cortex signal. Also, 490 ROIs of cNAWM were manually delineated to match the MS lesion positions. A two-tailed t-test showed that mean APTw values were higher in females than in males (t = 3.52, p < 0.001). Additionally, the mean APTw values of MS lesions were higher than those of cNAWM after accounting for covariates (mean lesion = 0.44, mean cNAWM = 0.13, F = 44.12, p < 0.001).The mean APTw values of ISO lesions were higher than those of cNAWM after accounting for covariates (mean ISO lesions = 0.42, mean cNAWM = 0.21, F = 12.12, p < 0.001). The mean APTw values of BH were also higher than those of cNAWM (mean BH lesions = 0.47, mean cNAWM = 0.033, F = 40.3, p < 0.001). The effect size (i.e., difference between lesion and cNAWM) for BH was found to be higher than for ISO (14 vs. 2). Diagnostic performance showed that APT was able to discriminate between all lesions and cNAWM with an accuracy of >75% (AUC = 0.79, SE = 0.014). Discrimination between ISO lesions and cNAWM was accomplished with an accuracy of >69% (AUC = 0.74, SE = 0.018), while discrimination between BH lesions and cNAWM was achieved at an accuracy of >80% (AUC = 0.87, SE = 0.021).

CONCLUSIONS

Our results highlight the potential of APTw imaging for use as a non-invasive technique that is able to provide essential molecular information to clinicians and researchers so that the stages of inflammation and degeneration in MS lesions can be better characterized.

CEST 2022 - Fast WASABI post-processing: Access to rapid B0 and B1 correction in clinical routine for CEST MRI

CEST MRI methods, such as APT and NOE imaging reveal biomarkers with significant diagnostic potential due to their ability to access molecular tissue information. Regardless of the technique used, CEST MRI data are affected by static magnetic B₀ and radiofrequency B₁ field inhomogeneities that degrade their contrast. For this reason, the correction of B₀ field-induced artefacts is essential, whereas accounting for B₁ field inhomogeneities have shown significant improvements in image readability. In a previous work, an MRI protocol called WASABI was presented, which can map simultaneously B₀ and B₁ field inhomogeneities, while maintaining the same sequence and readout types as used for CEST MRI. Despite the highly satisfactory quality of B₀ and B₁ maps computed from the WASABI data, the post-processing method is based on an exhaustive search of a four-parameter space and an additional four-parameter non-linear model fitting step. This leads to long post-processing times that are prohibitive in clinical practice. This work provides a new method for fast post-processing of WASABI data with outstanding acceleration of the parameter estimation procedure and without compromising its stability. The resulting computational acceleration makes the WASABI technique suitable for clinical use. The stability of the method is demonstrated on phantom data and clinical 3 Tesla in vivo data.

Linear projection-based chemical exchange saturation transfer parameter estimation

Isolated evaluation of multiparametric in vivo chemical exchange saturation transfer (CEST) MRI often requires complex computational processing for both correction of B0 and B1 inhomogeneity and contrast generation. For that, sufficiently densely sampled Z-spectra need to be acquired. The list of acquired frequency offsets largely determines the total CEST acquisition time, while potentially representing redundant information. In this work, a linear projection-based multiparametric CEST evaluation method is introduced that offers fast B0 and B1 inhomogeneity correction, contrast generation and feature selection for CEST data, enabling reduction of the overall measurement time. To that end, CEST data acquired at 7 T in six healthy subjects and in one brain tumor patient were conventionally evaluated by interpolation-based inhomogeneity correction and Lorentzian curve fitting. Linear regression was used to obtain coefficient vectors that directly map uncorrected data to corrected Lorentzian target parameters. L1-regularization was applied to find subsets of the originally acquired CEST measurements that still allow for such a linear projection mapping. The linear projection method allows fast and interpretable mapping from acquired raw data to contrast parameters of interest, generalizing from healthy subject training data to unseen healthy test data and to the tumor patient dataset. The L1-regularization method shows that a fraction of the acquired CEST measurements is sufficient to preserve tissue contrasts, offering up to a 2.8-fold reduction of scan time. Similar observations as for the 7-T data can be made for data from a clinical 3-T scanner. Being a fast and interpretable computation step, the proposed method is complementary to neural networks that have recently been employed for similar purposes. The scan time acceleration offered by the L1-regularization ("CEST-LASSO") constitutes a step towards better applicability of multiparametric CEST protocols in a clinical context.

The effects of RF coils and SAR supervision strategies for clinically applicable nonselective parallel-transmit pulses at 7 T

PURPOSE

To investigate the effects of using different parallel-transmit (pTx) head coils and specific absorption rate (SAR) supervision strategies on pTx pulse design for ultrahigh-field MRI using a 3D-MPRAGE sequence.

METHODS

The PTx universal pulses (UPs) and fast online-customized (FOCUS) pulses were designed with pre-acquired data sets (B₀ , B₁ ⁺ maps, specific absorption rate [SAR] supervision data) from two different 8 transmit/32 receive head coils on two 7T whole-body MR systems. For one coil, the SAR supervision model consisted of per-channel RF power limits. In the other coil, SAR estimations were done with both per-channel RF power limits as well as virtual observation points (VOPs) derived from electromagnetic field (EMF) simulations using three virtual human body models at three different positions. All pulses were made for nonselective excitation and inversion and evaluated on 132 B₀ , B₁ ⁺ , and SAR supervision datasets obtained with one coil and 12 from the other. At both sites, 3 subjects were examined using MPRAGE sequences that used UP/FOCUS pulses generated for both coils.

RESULTS

For some subjects, the UPs underperformed when simulated on a different coil from which they were derived, whereas FOCUS pulses still showed acceptable performance in that case. FOCUS inversion pulses outperformed adiabatic pulses when scaled to the same local SAR level. For the self-built coil, the use of VOPs showed reliable overestimation compared with the ground-truth EMF simulations, predicting about 52% lower local SAR for inversion pulses compared with per-channel power limits.

CONCLUSION

FOCUS inversion pulses offer a low-SAR alternative to adiabatic pulses and benefit from using EMF-based VOPs for SAR estimation.

A numerical human brain phantom for dynamic glucose-enhanced (DGE) MRI: On the influence of head motion at 3T

PURPOSE

Dynamic glucose-enhanced (DGE) MRI relates to a group of exchange-based MRI techniques where the uptake of glucose analogues is studied dynamically. However, motion artifacts can be mistaken for true DGE effects, while motion correction may alter true signal effects. The aim was to design a numerical human brain phantom to simulate a realistic DGE MRI protocol at 3T that can be used to assess the influence of head movement on the signal before and after retrospective motion correction.

METHODS

MPRAGE data from a tumor patient were used to simulate dynamic Z-spectra under the influence of motion. The DGE responses for different tissue types were simulated, creating a ground truth. Rigid head movement patterns were applied as well as physiological dilatation and pulsation of the lateral ventricles and head-motion-induced B0 -changes in presence of first-order shimming. The effect of retrospective motion correction was evaluated.

RESULTS

Motion artifacts similar to those previously reported for in vivo DGE data could be reproduced. Head movement of 1 mm translation and 1.5 degrees rotation led to a pseudo-DGE effect on the order of 1% signal change. B0 effects due to head motion altered DGE changes due to a shift in the water saturation spectrum. Pseudo DGE effects were partly reduced or enhanced by rigid motion correction depending on tissue location.

CONCLUSION

DGE MRI studies can be corrupted by motion artifacts. Designing post-processing methods using retrospective motion correction including B0 correction will be crucial for clinical implementation. The proposed phantom should be useful for evaluation and optimization of such techniques.

DeepCEST 7 T: Fast and homogeneous mapping of 7 T CEST MRI parameters and their uncertainty quantification

PURPOSE

In this work, we investigated the ability of neural networks to rapidly and robustly predict Lorentzian parameters of multi-pool CEST MRI spectra at 7 T with corresponding uncertainty maps to make them quickly and easily available for routine clinical use.

METHODS

We developed a deepCEST 7 T approach that generates CEST contrasts from just 1 scan with robustness against B₁ inhomogeneities. The input data for a neural feed-forward network consisted of 7 T in vivo uncorrected Z-spectra of a single B₁ level, and a B₁ map. The 7 T raw data were acquired using a 3D snapshot gradient echo multiple interleaved mode saturation CEST sequence. These inputs were mapped voxel-wise to target data consisting of Lorentzian amplitudes generated conventionally by 5-pool Lorentzian fitting of normalized, denoised, B₀ - and B₁ -corrected Z-spectra. The deepCEST network was trained with Gaussian negative log-likelihood loss, providing an uncertainty quantification in addition to the Lorentzian amplitudes.

RESULTS

The deepCEST 7 T network provides fast and accurate prediction of all Lorentzian parameters also when only a single B₁ level is used. The prediction was highly accurate with respect to the Lorentzian fit amplitudes, and both healthy tissues and hyperintensities in tumor areas are predicted with a low uncertainty. In corrupted cases, high uncertainty indicated wrong predictions reliably.

CONCLUSION

The proposed deepCEST 7 T approach reduces scan time by 50% to now 6:42 min, but still delivers both B₀ - and B₁ -corrected homogeneous CEST contrasts along with an uncertainty map, which can increase diagnostic confidence. Multiple accurate 7 T CEST contrasts are delivered within seconds.

Parallel transmit (pTx) with online pulse design for task-based fMRI at 7 T

PURPOSE

Parallel transmission (pTx) is an approach to improve image uniformity for ultra-high field imaging. In this study, we modified an echo planar imaging (EPI) sequence to design subject-specific pTx pulses online. We compared its performance against EPI with conventional circularly polarised (CP) pulses.

METHODS

We compared the pTx-EPI and CP-EPI sequences in a short EPI acquisition protocol and for two different functional paradigms in six healthy volunteers (2 female, aged 23-36 years, mean age 29.2 years). We chose two paradigms that are typically affected by signal dropout at 7 T: a visual objects localiser to determine face/scene selective brain regions and a semantic-processing task.

RESULTS

Across all subjects, pTx-EPI improved whole-brain mean temporal signal-to-noise ratio (tSNR) by 11.0% compared to CP-EPI. We also compared the ability of pTx-EPI and CP-EPI to detect functional activation for three contrasts over the two paradigms: face > object and scene > object for the visual objects localiser and semantic association > pattern matching for the semantic-processing paradigm. Across all three contrasts, pTx-EPI showed higher median z-scores and detected more active voxels in relevant areas, as determined from previous 3 T studies.

CONCLUSION

We have demonstrated a workflow for EPI acquisitions with online per-subject pulse calculations. We saw improved performance in both tSNR and functional acquisitions from pTx-EPI. Thus, we believe that online calculation pTx-EPI is robust enough for future fMRI studies, especially where activation is expected in brain areas liable to significant signal dropout.

Interleaved binomial kT-Points for water-selective imaging at 7T

PURPOSE

We present a time-efficient water-selective, parallel transmit RF excitation pulse design for ultra-high field applications.

METHODS

The proposed pulse design method achieves flip angle homogenization at ultra-high fields by employing spatially nonselective k T $$ {\mathrm{k}}_T $$ -points pulses. In order to introduce water-selection, the concept of binomial pulses is applied. Due to the composite nature of k T $$ {\mathrm{k}}_T $$ -points, the pulse can be split into multiple binomial subpulse blocks shorter than half the precession period of fat, that are played out successively. Additional fat precession turns, that would otherwise impair the spectral response, can thus be avoided. Bloch simulations of the proposed interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses were carried out and compared in terms of duration, homogeneity, fat suppression and pulse energy. For validation, in vivo MP-RAGE and 3D-EPI data were acquired.

RESULTS

Simulation results show that interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses achieve shorter total pulse durations, improved flip angle homogeneity and more robust fat suppression compared to available methods. Interleaved binomial k T $$ {\mathrm{k}}_T $$ -points can be customized by changing the number of k T $$ {\mathrm{k}}_T $$ -points, the subpulse duration and the order of the binomial pulse. Using shorter subpulses, the number of k T $$ {\mathrm{k}}_T $$ -points can be increased and hence better homogeneity is achieved, while still maintaining short total pulse durations. Flip angle homogenization and fat suppression of interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses is demonstrated in vivo at 7T, confirming Bloch simulation results.

CONCLUSION

In this work, we present a time efficient and robust parallel transmission technique for nonselective water excitation with simultaneous flip angle homogenization at ultra-high field.

Parallel transmit pulse design for saturation homogeneity (PUSH) for magnetization transfer imaging at 7T

PURPOSE

This work proposes a novel RF pulse design for parallel transmit (pTx) systems to obtain uniform saturation of semisolid magnetization for magnetization transfer (MT) contrast in the presence of transmit field B1+ inhomogeneities. The semisolid magnetization is usually modeled as being purely longitudinal, with the applied B1+ field saturating but not rotating its magnetization; thus, standard pTx pulse design methods do not apply.

THEORY AND METHODS

Pulse design for saturation homogeneity (PUSH) optimizes pTx RF pulses by considering uniformity of root-mean squared B1+ , B1rms , which relates to the rate of semisolid saturation. Here we considered designs consisting of a small number of spatially non-selective sub-pulses optimized over either a single 2D plane or 3D. Simulations and in vivo experiments on a 7T Terra system with an 8-TX Nova head coil in five subjects were carried out to study the homogenization of B1rms and of the MT contrast by acquiring MT ratio maps.

RESULTS

Simulations and in vivo experiments showed up to six and two times more uniform B1rms compared to circular polarized (CP) mode for 2D and 3D optimizations, respectively. This translated into 4 and 1.25 times more uniform MT contrast, consistently for all subjects, where two sub-pulses were enough for the implementation and coil used.

CONCLUSION

The proposed PUSH method obtains more uniform and higher MT contrast than CP mode within the same specific absorption rate (SAR) budget.

The Value of APTw CEST MRI in Routine Clinical Assessment of Human Brain Tumor Patients at 3T

Background. With fast-growing evidence in literature for clinical applications of chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI), this prospective study aimed at applying amide proton transfer-weighted (APTw) CEST imaging in a clinical setting to assess its diagnostic potential in differentiation of intracranial tumors at 3 tesla (T). Methods. Using the asymmetry magnetization transfer ratio (MTRasym) analysis, CEST signals were quantitatively investigated in the tumor areas and in a similar sized region of the normal-appearing white matter (NAWM) on the contralateral hemisphere of 27 patients with intracranial tumors. Area under curve (AUC) analyses were used and results were compared to perfusion-weighted imaging (PWI). Results. Using APTw CEST, contrast-enhancing tumor areas showed significantly higher APTw CEST metrics than contralateral NAWM (AUC = 0.82; p < 0.01). In subgroup analyses of each tumor entity vs. NAWM, statistically significant effects were yielded for glioblastomas (AUC = 0.96; p < 0.01) and for meningiomas (AUC = 1.0; p < 0.01) but not for lymphomas as well as metastases (p > 0.05). PWI showed results comparable to APTw CEST in glioblastoma (p < 0.01). Conclusions. This prospective study confirmed the high diagnostic potential of APTw CEST imaging in a routine clinical setting to differentiate brain tumors.

Measuring Arterial Pulsatility With Dynamic Inflow Magnitude Contrast

The pulsatility of blood flow through cerebral arteries is clinically important, as it is intrinsically associated with cerebrovascular health. In this study we outline a new MRI approach to measuring the real-time pulsatile flow in cerebral arteries, which is based on the inflow phenomenon associated with fast gradient-recalled-echo acquisitions. Unlike traditional phase-contrast techniques, this new method, which we dub dynamic inflow magnitude contrast (DIMAC), does not require velocity-encoding gradients as sensitivity to flow velocity is derived purely from the inflow effect. We achieved this using a highly accelerated single slice EPI acquisition with a very short TR (15 ms) and a 90° flip angle, thus maximizing inflow contrast. We simulate the spoiled GRE signal in the presence of large arteries and perform a sensitivity analysis. The sensitivity analysis demonstrates that in the regime of high inflow contrast, DIMAC shows much greater sensitivity to flow velocity over blood volume changes. We support this theoretical prediction with in-vivo data collected in two separate experiments designed to demonstrate the utility of the DIMAC signal contrast. We perform a hypercapnia challenge experiment in order to experimentally modulate arterial tone within subjects, and thus modulate the arterial pulsatile flow waveform. We also perform a thigh-cuff release challenge, designed to induce a transient drop in blood pressure, and demonstrate that the continuous DIMAC signal captures the complex transient change in the pulsatile and non-pulsatile components of flow. In summary, this study proposes a new role for a well-established source of MR image contrast and demonstrates its potential for measuring both steady-state and dynamic changes in arterial tone.

CEST MRI provides amide/amine surrogate biomarkers for treatment-naïve glioma sub-typing

PURPOSE

Accurate glioma classification affects patient management and is challenging on non- or low-enhancing gliomas. This study investigated the clinical value of different chemical exchange saturation transfer (CEST) metrics for glioma classification and assessed the diagnostic effect of the presence of abundant fluid in glioma subpopulations.

METHODS

Forty-five treatment-naïve glioma patients with known isocitrate dehydrogenase (IDH) mutation and 1p/19q codeletion status received CEST MRI (B_1rms = 2μT, T_sat = 3.5 s) at 3 T. Magnetization transfer ratio asymmetry and CEST metrics (amides: offset range 3-4 ppm, amines: 1.5-2.5 ppm, amide/amine ratio) were calculated with two models: 'asymmetry-based' (AB) and 'fluid-suppressed' (FS). The presence of T2/FLAIR mismatch was noted.

RESULTS

IDH-wild type had higher amide/amine ratio than IDH-mutant_1p/19q^codel (p < 0.022). Amide/amine ratio and amine levels differentiated IDH-wild type from IDH-mutant (p < 0.0045) and from IDH-mutant_1p/19q^ret (p < 0.021). IDH-mutant_1p/19q^ret had higher amides and amines than IDH-mutant_1p/19q^codel (p < 0.035). IDH-mutant_1p/19q^ret with AB/FS mismatch had higher amines than IDH-mutant_1p/19q^ret without AB/FS mismatch ( < 0.016). In IDH-mutant_1p/19q^ret, the presence of AB/FS mismatch was closely related to the presence of T2/FLAIR mismatch (p = 0.014).

CONCLUSIONS

CEST-derived biomarkers for amides, amines, and their ratio can help with histomolecular staging in gliomas without intense contrast enhancement. T2/FLAIR mismatch is reflected in the presence of AB/FS CEST mismatch. The AB/FS CEST mismatch identifies glioma subgroups that may have prognostic and clinical relevance.

Rigid motion-resolved B 1 + prediction using deep learning for real-time parallel-transmission pulse design

Purpose

Tailored parallel-transmit (pTx) pulses produce uniform excitation profiles at 7 T, but are sensitive to head motion. A potential solution is real-time pulse redesign. A deep learning framework is proposed to estimate pTx B₁⁺ distributions following within-slice motion, which can then be used for tailored pTx pulse redesign.

Methods

Using simulated data, conditional generative adversarial networks were trained to predict B₁⁺ distributions in the head following a displacement. Predictions were made for two virtual body models that were not included in training. Predicted maps were compared with groundtruth (simulated, following motion) B₁⁺ maps. Tailored pTx pulses were designed using B₁⁺ maps at the original position (simulated, no motion) and evaluated using simulated B₁⁺ maps at displaced position (ground-truth maps) to quantify motion-related excitation error. A second pulse was designed using predicted maps (also evaluated on ground-truth maps) to investigate improvement offered by the proposed method.

Results

Predicted B₁⁺ maps corresponded well with ground-truth maps. Error in predicted maps was lower than motion-related error in 99% and 67% of magnitude and phase evaluations, respectively. Worst-case flip-angle normalized RMS error due to motion (76% of target flip angle) was reduced by 59% when pulses were redesigned using predicted maps.

Conclusion

We propose a framework for predicting B₁⁺ maps online with deep neural networks. Predicted maps can then be used for real-time tailored pulse redesign, helping to overcome head motion–related error in pTx.

Whole-brain quantitative CEST MRI at 7T using parallel transmission methods and B 1 + correction

PURPOSE

To enable whole-brain quantitative CEST MRI at ultra-high magnetic field strengths (B₀ ≥ 7T) within short acquisition times.

METHODS

Multiple interleaved mode saturation (MIMOSA) was combined with fast online-customized (FOCUS) parallel transmission (pTx) excitation pulses and B 1 + correction to achieve homogenous whole-brain coverage. Examinations of 13 volunteers were performed on a 7T MRI system with 3 different types of pulse sequences: (1) saturation in circular polarized (CP) mode and CP mode readout, (2) MIMOSA and CP readout, and (3) MIMOSA and FOCUS readout. For comparison, the inverse magnetic transfer ratio metric for relayed nuclear Overhauser effect and amide proton transfer were calculated. To investigate the number of required acquisitions for a good B 1 + correction, 4 volunteers were measured with 6 different B₁ amplitudes. Finally, time point repeatability was investigated for 6 volunteers.

RESULTS

MIMOSA FOCUS sequence using B 1 + correction, with both single and multiple points, reduced inhomogeneity of the CEST contrasts around the occipital lobe and cerebellum. Results indicate that the most stable inter-subject coefficient of variation was achieved using the MIMOSA FOCUS sequence. Time point repeatability of MIMOSA FOCUS with single-point B 1 + correction showed a maximum coefficient of variation below 8% for 3 measurements in a single volunteer.

CONCLUSION

A combination of MIMOSA FOCUS with a single-point B 1 + correction can be used to achieve quantitative CEST measurements at ultra-high magnetic field strengths. Compared to previous B 1 + correction methods, acquisition time can be reduced as additional scans required for B 1 + correction can be omitted.

Fast online-customized (FOCUS) parallel transmission pulses: A combination of universal pulses and individual optimization

PURPOSE

To mitigate spatial flip angle (FA) variations under strict specific absorption rate (SAR) constraints for ultra-high field MRI using a combination of universal parallel transmit (pTx) pulses and fast subject-specific optimization.

METHODS

Data sets consisting of B₀ , B 1 + maps, and virtual observation point (VOP) data were acquired from 72 subjects (study groups of 48/12 healthy Europeans/Asians and 12 Europeans with pathological or incidental findings) using an 8Tx/32Rx head coil on a 7T whole-body MR system. Combined optimization values (COV) were defined as combination of spiral-nonselective (SPINS) trajectory parameters and an energy regularization weight. A set of COV was optimized universally by simulating the individual RF pulse optimizations of 12 training data sets (healthy Europeans). Subsequently, corresponding universal pulses (UPs) were calculated. Using COV and UPs, individually optimized pulses (IOPs) were calculated during the sequence preparation phase (maximum 15 s). Two different UPs and IOPs were evaluated by calculating their normalized root-mean-square error (NRMSE) of the FA and SAR in simulations of all data sets. Seven additional subjects were examined using an MPRAGE sequence that uses the designed pTx excitation pulses and a conventional adiabatic inversion.

RESULTS

All pTx pulses resulted in decreased mean NRMSE compared to a circularly polarized (CP) pulse (CP = ~28%, UPs = ~17%, and IOPs = ~12%). UPs and IOPs improved homogeneity for all subjects. Differences in NRMSE between study groups were much lower than differences between different pulse types.

CONCLUSION

UPs can be used to generate fast online-customized (FOCUS) pulses gaining lower NRMSE and/or lower SAR values.

High resolution time-of-flight MR-angiography at 7 T exploiting VERSE saturation, compressed sensing and segmentation

BACKGROUND

3D Time-of-Flight (TOF) MR-angiography (MRA) substantially benefits from ultra-high magnetic field strengths (≥7 T) due to increased Signal-to-Noise ratio and improved contrast. However, high-resolution TOF-MRA usually requires long acquisition times. In addition, specific absorption rate constraints limit the choice of optimal pulse sequence parameters, especially if venous saturation is employed.

PURPOSE

To implement and evaluate an arterial TOF-MRA for accelerated high-resolution angiography at ultra-high magnetic field strength.

FIELD STRENGTHS/SEQUENCE

7 T modified gradient-echo TOF sequence including venous saturation using Variable-Rate Selective Excitation (VERSE), Compressed Sensing (CS) and sparse application of saturation pulses, called segmentation, were included for acceleration.

ASSESSMENT

To analyze the acceleration techniques all volunteers were examined with the same protocols. CS with different sampling patterns and regularization factors as well as segmentation were applied for acceleration. For comparison, conventional acceleration techniques were applied (GRAPPA PAT 3 and Partial Fourier (6/8 in slice/phase encoding)). Images were co-registered and 40 mm thick transversal maximum intensity projections were created to calculate the relative number of vessels. To analyze the visibility of small vessels, the lenticulostriate arteries (LSA) were examined. This was done via multiscale vessel enhancement filtering in a VOI and quantification via Fiji ImageJ as well as qualitatively evaluation by two radiologists. Additionally, the venous/arterial vessel-to-background ratios (vVBR/aVBR) were calculated for chosen protocols.

RESULTS

For the acceleration of a high resolution TOF-MRA (0.31 mm isotropic), under-sampling of 9.6 showed aliasing artifacts, whereas 7.2 showed no aliasing. The regularization factor R had a strong impact on the image quality according to smoothing (R = 0.01 to R = 0.005) and noise (R = 0.0005 to R = 0.00005). With the alternating sampling patterns it was shown that the k-space center should not be under-sampled too much. Additionally segmentation could be verified to be feasible for stronger acceleration with sufficient venous suppression.

CONCLUSION

The combination of several independent techniques (VERSE, CS with acceleration factor 7.2, R = 0.001, Poisson disc radius of 80%, 3 segments) enables the application of high-resolution (0.31 mm isotropic) TOF-MRA with venous saturation at 7 T in clinical time settings (TA ≈ 5 min) and within the SAR limits.