Feasibility of a fast method for B1-inhomogeneity correction for FSPGR sequences

Using Phase Difference Information to Detect Errors in the Flip Angle Measured with Actual Flip Angle Imaging at 7T

Flip angle (FA) measurements using the actual flip angle imaging (AFI) method may induce significant errors in ultrahigh fields. We aimed to develop a method for detecting errors in FA measurements using phase information at 7 tesla. We performed computer simulations to elucidate the relationship between the FA calculation errors and the phase difference between the two AFI source images. We then examined whether a method based on the phase difference could detect FA calculation errors and determine the prescribed nominal FA of the scanner for accurate measurements in phantoms and healthy volunteers. The simulations confirmed that the calculated FA values erroneously decreased when the longitudinal magnetization and phase in one of the source images were inverted. Tests on phantoms and human subjects demonstrated that the phase difference information between the source images with a cut-off of 90° could readily detect FA calculation errors in the AFI method.

3D T1 relaxation time measurements in an equine model of subtle post-traumatic osteoarthritis using MB-SWIFT

The aim of this study is to assess whether articular cartilage changes in an equine model of post-traumatic osteoarthritis (PTOA), induced by surgical creation of standard (blunt) grooves, and very subtle sharp grooves, could be detected with ex vivo T1 relaxation time mapping utilizing three-dimensional (3D) readout sequence with zero echo time. Grooves were made on the articular surfaces of the middle carpal and radiocarpal joints of nine mature Shetland ponies and osteochondral samples were harvested at 39 weeks after being euthanized under respective ethical permissions. T1 relaxation times of the samples (n = 8 + 8 for experimental and n = 12 for contralateral controls) were measured with a variable flip angle 3D multiband-sweep imaging with Fourier transform sequence. Equilibrium and instantaneous Young's moduli and proteoglycan (PG) content from OD of Safranin-O-stained histological sections were measured and utilized as reference parameters for the T1 relaxation times. T1 relaxation time was significantly (p < 0.05) increased in both groove areas, particularly in the blunt grooves, compared with control samples, with the largest changes observed in the superficial half of the cartilage. T1 relaxation times correlated weakly (Rs  ≈ 0.33) with equilibrium modulus and PG content (Rs  ≈ 0.21). T1 relaxation time in the superficial articular cartilage is sensitive to changes induced by the blunt grooves but not to the much subtler sharp grooves, at the 39-week timepoint post-injury. These findings support that T1 relaxation time has potential in detection of mild PTOA, albeit the most subtle changes could not be detected.

Voxel-based mapping of five MR biomarkers in the wrist bone marrow

OBJECTIVE

MRI is a reliable and accurate technique to characterize rheumatoid arthritis. The aim of this study was to provide voxel-by-voxel 3D maps of the proton density fat fraction (PDFF), the T₁ of water (T_1W), the T₁ of fat (T_1F), the T₂* of water (T₂*_W), the T₂* of fat (T₂*_F) in the wrist bone marrow.

MATERIALS AND METHODS

The experiments were conducted on 14 healthy volunteers (mean age: 24 ± 4). The data were acquired at 1.5 T using two optimized four-echo 3D 1.2 × 1.2 × 1.2 mm³-isotropic spoiled gradient sequences. A repeatability study was carried out. The measurements were done using a homemade parametric viewer software.

RESULTS

The inter-volunteer results were, on average: PDFF = 86 ± 3%, T_1W = 441 ± 113 ms, T_1F = 245 ± 19 ms, T₂*_W = 6 ± 1 ms and T₂*_F = 16 ± 3 ms. The coefficients of variation were for fat based biomarkers CV_PDFF < 5%, CV_T1F < 15% and CV_T2*F < 10% in the repeatability study.

DISCUSSION

The protocol and quantification tool proposed in this study provide high-resolution voxel-by-voxel 3D maps of five biomarkers in the wrist in less than 4 min of acquisition.

Use of Multiple-Tube Phantom: A Method to Globally Correct Native T1 Relaxation Time Inhomogeneity in Dedicated Molecular Magnetic Resonance Breast Coil

Background:

Native T1 relaxation time (T1₀) presents an important prerequisite to reliably quantify pharmacokinetic parameter like K^trans (volume transfer constant). Native T1 value can be varied because of the inhomogeneity in the breast coil, thus influencing the K^trans measurement.

Purpose:

The current study aims to design and use a phantom with multiple tubes for both breast cuffs to assess native T1 inhomogeneity across the dedicated molecular magnetic resonance (mMR) breast coil and adopt corrective method to spatially normalize T1 values to improve homogeneity.

Materials and Methods:

Two phantoms with multiple tubes (19 tubes) specially designed and filled with contrast medium with known T1 value were placed in each mMR breast coil cuff. Native T1 at various spatial locations was calculated applying dual flip angle sequence. Correction factors were derived at various spatial locations as a function of deviation of the native T1 value from phantom and applied to correct the native T1 relaxation time.

Results:

A statistically significant difference between native T1 values of the right and left anterior (P = 0.0095), middle (P = 0.0081), and posterior (P = 0.0004) parts of the breast coil. No significant difference was seen in the corrected T1 values between anterior (P = 0.402), middle (P = 0.305), and posterior (P = 0.349) aspects of both sides of the breast coil.

Conclusion:

Inhomogeneity in the native T1 value exists in dedicated mMR breast coil, and significant improvement can be achieved using specially designed external phantom with multiple tubes.

A standard system phantom for magnetic resonance imaging

Purpose

A standard MRI system phantom has been designed and fabricated to assess scanner performance, stability, comparability and assess the accuracy of quantitative relaxation time imaging. The phantom is unique in having traceability to the International System of Units, a high level of precision, and monitoring by a national metrology institute. Here, we describe the phantom design, construction, imaging protocols, and measurement of geometric distortion, resolution, slice profile, signal‐to‐noise ratio (SNR), proton‐spin relaxation times, image uniformity and proton density.

Methods

The system phantom, designed by the International Society of Magnetic Resonance in Medicine ad hoc committee on Standards for Quantitative MR, is a 200 mm spherical structure that contains a 57‐element fiducial array; two relaxation time arrays; a proton density/SNR array; resolution and slice‐profile insets. Standard imaging protocols are presented, which provide rapid assessment of geometric distortion, image uniformity, T₁ and T₂ mapping, image resolution, slice profile, and SNR.

Results

Fiducial array analysis gives assessment of intrinsic geometric distortions, which can vary considerably between scanners and correction techniques. This analysis also measures scanner/coil image uniformity, spatial calibration accuracy, and local volume distortion. An advanced resolution analysis gives both scanner and protocol contributions. SNR analysis gives both temporal and spatial contributions.

Conclusions

A standard system phantom is useful for characterization of scanner performance, monitoring a scanner over time, and to compare different scanners. This type of calibration structure is useful for quality assurance, benchmarking quantitative MRI protocols, and to transition MRI from a qualitative imaging technique to a precise metrology with documented accuracy and uncertainty.

Analysis Protocols for MRI Mapping of Renal T1

The computation of T₁ maps from MR datasets represents an important step toward the precise characterization of kidney disease models in small animals. Here the main strategies to analyze renal T₁ mapping datasets derived from small rodents are presented. Suggestions are provided with respect to essential software requirements, and advice is provided as to how dataset completeness and quality may be evaluated. The various fitting models applicable to T₁ mapping are presented and discussed. Finally, some methods are proposed for validating the obtained results.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concept and experimental procedure.

Confounding factors in multi-parametric q-MRI protocol: A study of bone marrow biomarkers at 1.5 T

OBJECT

The MRI tissue characterization of vertebral bone marrow includes the measurement of proton density fat fraction (PDFF), T₁ and T₂* relaxation times of the water and fat components (T_1W, T_1F, T₂*_W, T₂*_F), IVIM diffusion D, perfusion fraction f and pseudo-diffusion coefficient D*. However, the measurement of these vertebral bone marrow biomarkers (VBMBs) is affected with several confounding factors. In the current study, we investigated these confounding factors including the regional variation taking the example of variation between the anterior and posterior area in lumbar vertebrae, B₁ inhomogeneity and the effect of fat suppression on f.

MATERIALS AND METHODS

A fat suppressed diffusion-weighted sequence and two 3D gradient multi-echo sequences were used for the measurements of the seven VBMBs. A turbo flash B₁ map sequence was used to estimate B₁ inhomogeneities and thus, to correct flip angle for T₁ quantification. We introduced a correction to perfusion fraction f measured with fat suppression, namely f_PDFF.

RESULTS

A significant difference in the values of PDFF, f and f_PDFF, T_1F, T₂*_W and D was observed between the anterior and posterior region. Although, little variations of flip angle were observed in this anterior-posterior direction in one vertebra but larger variations were observed in head-feet direction from L1 to L5 vertebrae.

DISCUSSION

The regional difference in PDFF, f_PDFF and T₂*_W can be ascribed to differences in the trabecular bone density and vascular network within vertebrae. The regional variation of VBMBs shows that care should be taken in reproducing the same region-of-interest location along a longitudinal study. The same attention should be taken while measuring f in fatty environment, and measuring T₁. Furthermore, the MRI-protocol presented here allows for measurements of seven VBMBs in less than 6 min and is of interest for longitudinal studies of bone marrow diseases.

Iron leakage owing to blood–brain barrier disruption in small vessel disease CADASIL

Objective

To assess the relationship among iron accumulation, blood–brain barrier (BBB) damage, and cognitive function in patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL).

Methods

We enrolled 21 patients with NOTCH3 mutations and 21 age-matched healthy controls in this cross-sectional study. All participants underwent global physical and cognitive assessments and brain MRI using voxel-based quantitative susceptibility mapping (QSM; iron deposition measure) and dynamic contrast-enhanced MRI (BBB permeability measure). We compared behavioral and imaging data between the groups and analyzed the correlations in each group.

Results

Among 21 NOTCH3 mutation carriers, 10 were symptomatic and 11 asymptomatic. Montreal Cognitive Assessment scores were significantly different among the groups (symptomatic < asymptomatic < control participants). Voxel-based QSM analysis revealed that the symptomatic group had higher QSM values than did the asymptomatic group in the putamen, caudate nucleus, temporal pole, and centrum semiovale. These QSM values were positively correlated with regional BBB permeabilities (putamen: r = 0.57, p = 0.006; caudate nucleus: r = 0.51, p = 0.019; temporal pole: r = 0.48, p = 0.030; centrum semiovale: r = 0.45, p = 0.044) and negatively correlated with Montreal Cognitive Assessment scores (caudate nucleus: r = −0.53, p = 0.012; temporal pole: r = −0.56, p = 0.008).

Conclusions

This study showed that cerebral iron burden was associated with regional BBB permeability and cognitive dysfunction in patients with CADASIL, highlighting the potential of these imaging techniques as auxiliary biomarkers to monitor the course of small vessel disease.

B1 inhomogeneity correction of RARE MRI with transceive surface radiofrequency probes

PURPOSE

The use of surface radiofrequency (RF) coils is common practice to boost sensitivity in (pre)clinical MRI. The number of transceive surface RF coils is rapidly growing due to the surge in cryogenically cooled RF technology and ultrahigh-field MRI. Consequently, there is an increasing need for effective correction of the excitation field ( B 1 + ) inhomogeneity inherent in these coils. Retrospective B₁ correction permits quantitative MRI, but this usually requires a pulse sequence-specific analytical signal intensity (SI) equation. Such an equation is not available for fast spin-echo (Rapid Acquisition with Relaxation Enhancement, RARE) MRI. Here we present, test, and validate retrospective B₁ correction methods for RARE.

METHODS

We implemented the commonly used sensitivity correction and developed an empirical model-based method and a hybrid combination of both. Tests and validations were performed with a cryogenically cooled RF probe and a single-loop RF coil. Accuracy of SI quantification and T₁ contrast were evaluated after correction.

RESULTS

The three described correction methods achieved dramatic improvements in B₁ homogeneity and significantly improved SI quantification and T₁ contrast, with mean SI errors reduced from >40% to >10% following correction in all cases. Upon correction, images of phantoms and mouse heads demonstrated homogeneity comparable to that of images acquired with a volume resonator. This was quantified by SI profile, SI ratio (error < 10%), and percentage of integral uniformity (PIU > 80% in vivo and ex vivo compared to PIU > 87% with the reference RF coil).

CONCLUSION

This work demonstrates the efficacy of three B₁ correction methods tailored for transceive surface RF probes and RARE MRI. The corrected images are suitable for quantification and show comparable results between the three methods, opening the way for T₁ measurements and X-nuclei quantification using surface transceiver RF coils. This approach is applicable to other MR techniques for which no analytical SI exists.

A pharmacokinetic model including arrival time for two inputs and compensating for varying applied flip-angle in dynamic gadoxetic acid-enhanced MR imaging

PURPOSE

Pharmacokinetic models facilitate assessment of properties of the micro-vascularization based on DCE-MRI data. However, accurate pharmacokinetic modeling in the liver is challenging since it has two vascular inputs and it is subject to large deformation and displacement due to respiration.

METHODS

We propose an improved pharmacokinetic model for the liver that (1) analytically models the arrival-time of the contrast agent for both inputs separately; (2) implicitly compensates for signal fluctuations that can be modeled by varying applied flip-angle e.g. due to B1-inhomogeneity. Orton's AIF model is used to analytically represent the vascular input functions. The inputs are independently embedded into the Sourbron model. B1-inhomogeneity-driven variations of flip-angles are accounted for to justify the voxel's displacement with respect to a pre-contrast image.

RESULTS

The new model was shown to yield lower root mean square error (RMSE) after fitting the model to all but a minority of voxels compared to Sourbron's approach. Furthermore, it outperformed this existing model in the majority of voxels according to three model-selection criteria.

CONCLUSION

Our work primarily targeted to improve pharmacokinetic modeling for DCE-MRI of the liver. However, other types of pharmacokinetic models may also benefit from our approaches, since the techniques are generally applicable.

Influence of B1-Inhomogeneity on Pharmacokinetic Modeling of Dynamic Contrast-Enhanced MRI: A Simulation Study

OBJECTIVE

To simulate the B₁-inhomogeneity-induced variation of pharmacokinetic parameters on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

MATERIALS AND METHODS

B₁-inhomogeneity-induced flip angle (FA) variation was estimated in a phantom study. Monte Carlo simulation was performed to assess the FA-deviation-induced measurement error of the pre-contrast R₁, contrast-enhancement ratio, Gd-concentration, and two-compartment pharmacokinetic parameters (K^trans, v_e, and v_p).

RESULTS

B₁-inhomogeneity resulted in -23-5% fluctuations (95% confidence interval [CI] of % error) of FA. The 95% CIs of FA-dependent % errors in the gray matter and blood were as follows: -16.7-61.8% and -16.7-61.8% for the pre-contrast R₁, -1.0-0.3% and -5.2-1.3% for the contrast-enhancement ratio, and -14.2-58.1% and -14.1-57.8% for the Gd-concentration, respectively. These resulted in -43.1-48.4% error for K^trans, -32.3-48.6% error for the v_e, and -43.2-48.6% error for v_p. The pre-contrast R₁ was more vulnerable to FA error than the contrast-enhancement ratio, and was therefore a significant cause of the Gd-concentration error. For example, a -10% FA error led to a 23.6% deviation in the pre-contrast R₁, -0.4% in the contrast-enhancement ratio, and 23.6% in the Gd-concentration. In a simulated condition with a 3% FA error in a target lesion and a -10% FA error in a feeding vessel, the % errors of the pharmacokinetic parameters were -23.7% for K^trans, -23.7% for v_e, and -23.7% for v_p.

CONCLUSION

Even a small degree of B₁-inhomogeneity can cause a significant error in the measurement of pharmacokinetic parameters on DCE-MRI, while the vulnerability of the pre-contrast R₁ calculations to FA deviations is a significant cause of the miscalculation.

Spiral Perfusion Imaging With Consecutive Echoes (SPICE™) for the Simultaneous Mapping of DSC- and DCE-MRI Parameters in Brain Tumor Patients: Theory and Initial Feasibility

Dynamic contrast-enhanced (DCE) and dynamic susceptibility contrast (DSC) magnetic resonance imaging (MRI) are the perfusion imaging techniques most frequently used to probe the angiogenic character of brain neoplasms. With these methods, T₁- and T₂/T₂*-weighted imaging sequences are used to image the distribution of gadolinium (Gd)-based contrast agents. However, it is well known that Gd exhibits combined T₁, T₂, and T₂* shortening effects in tissue, and therefore, the results of both DCE- and DSC-MRI can be confounded by these opposing effects. In particular, residual susceptibility effects compete with T₁ shortening, which can confound DCE-MRI parameters, whereas dipolar T₁ and T₂ leakage and residual susceptibility effects can confound DSC-MRI parameters. We introduce here a novel perfusion imaging acquisition and postprocessing method termed Spiral Perfusion Imaging with Consecutive Echoes (SPICE) that can be used to simultaneously acquire DCE- and DSC-MRI data, which requires only a single dose of the Gd contrast agent, does not require the collection of a precontrast T₁ map for DCE-MRI processing, and eliminates the confounding contrast agent effects due to contrast extravasation. A detailed mathematical description of SPICE is provided here along with a demonstration of its utility in patients with high-grade glioma.

Difference in the intratumoral distributions of extracellular-fluid and intravascular MR contrast agents in glioblastoma growth

Contrast enhancement by an extracellular-fluid contrast agent (CA) (Gd-DOTA) depends primarily on the blood-brain-barrier permeability (bp), and transverse-relaxation change caused by intravascular T₂ CA (superparamagnetic iron oxide nanoparticles, SPIONs) is closely associated with the blood volume (BV). Pharmacokinetic (PK) vascular characterization based on single-CA-using dynamic contrast-enhanced MRI (DCE-MRI) has shown significant measurement variation according to the molecular size of the CA. Based on this recognition, this study used a dual injection of Gd-DOTA and SPIONs for tracing the changes of bp and BV in C6 glioma growth (Days 1 and 7 after the tumor volume reached 2 mL). bp was quantified according to the non-PK parameters of Gd-DOTA-using DCE-MRI (wash-in rate, maximum enhancement ratio and initial area under the enhancement curve (IAUC)). BV was estimated by SPION-induced ΔR₂ * and ΔR₂ . With validated measurement reliability of all the parameters (coefficients of variation ≤10%), dual-contrast MRI demonstrated a different region-oriented distribution between Gd-DOTA and SPIONs within a tumor as follows: (a) the BV increased stepwise from the tumor center to the periphery; (b) the tumor periphery maintained the augmented BV to support continuous tumor expansion from Day 1 to Day 7; (c) the internal tumor area underwent significant vascular shrinkage (i.e. decreased ΔR₂ and ΔR₂ ) as the tumor increased in size; (d) the tumor center showed greater bp-indicating parameters, i.e. wash-in rate, maximum enhancement ratio and IAUC, than the periphery on both Days 1 and 7 and (e) the tumor center showed a greater increase of bp than the tumor periphery in tumor growth, as suggested to support tumor viability when there is insufficient blood supply. In the MRI-histologic correlation, a prominent BV increase in the tumor periphery seen in MRI was verified with increased fluorescein isothiocyanate-dextran signals and up-regulated immunoreactivity of CD31-VEGFR. In conclusion, the spatiotemporal alterations of BV and bp in glioblastoma growth, i.e. augmented BV in the tumor periphery and increased bp in the center, can be sufficiently evaluated by MRI with dual injection of extracellular-fluid Gd chelates and intravascular SPION.