Time-resolved MRI Research Articles

Free-breathing time-resolved 4D MRI with improved T1-weighting contrast.

This work proposes MP-Grasp4D (magnetization-prepared golden-angle radial sparse parallel 4D) MRI, a free-breathing, inversion recovery (IR)-prepared, time-resolved 4D MRI technique with improved T1-weighted contrast. MP-Grasp4D MRI acquisition incorporates IR preparation into a radial gradient echo sequence. MP-Grasp4D employs a golden-angle navi-stack-of-stars sampling scheme, where imaging data of rotating radial stacks and navigator stacks (acquired at a consistent rotation angle) are alternately acquired. The navigator stacks are used to estimate a temporal basis for low-rank subspace-constrained reconstruction. This allows for the simultaneous capture of both IR-induced contrast changes and respiratory motion. One temporal frame of the imaging volume in MP-Grasp4D MRI is reconstructed from a single stack and an adjacent navigator stack on average, resulting in a nominal temporal resolution of 0.16 seconds per volume. Images corresponding to the optimal inversion time (TI) can be retrospectively selected for providing the best image contrast. Reader studies were conducted to assess the performance of MP-Grasp4D MRI in liver imaging across 30 subjects in comparison with standard Grasp4D MRI without IR preparation. MP-Grasp4D MRI received significantly higher scores (P < 0.05) than Grasp4D in all assessment categories. There was a moderate to almost perfect agreement (kappa coefficient from 0.42 to 0.9) between the two readers for image quality assessment. When the scan time is reduced, MP-Grasp4D MRI preserves image contrast and quality, demonstrating additional acceleration capability. MP-Grasp4D MRI improves T1-weighted contrast for free-breathing time-resolved 4D MRI and eliminates the need for explicit motion compensation. This method is expected to be valuable in different MRI applications such as MR-guided radiotherapy.

Ascending aortic impedance in young endurance athletes: a time-resolved phase-contrast MRI study.

Ventricular-vascular coupling in endurance athletes remains incompletely understood. The purpose of this study was to determine the ascending aortic impedance in endurance athletes and explore its associations with traditional cardiovascular measurements. In 15 young male endurance runners and 19 young healthy men, time-resolved (CINE) two-dimensional (2-D) phase-contrast MRI quantified the ascending aortic flow while the pressure waveform was simultaneously collected via a generalized transfer function. The aortic impedance modulus and phase were calculated in the frequency domain while characteristic impedance (ZcF) was calculated by averaging moduli between the 4th and 8th heart rate (HR) harmonics. Stroke volume (SV), left ventricular (LV) morphometry, double product, aortic compliance, and total peripheral resistance (TPR) were also measured. Endurance athletes had higher SV, slower HR, greater LV end-diastolic volume and mass, and lower double product than sedentary participants (all P < 0.05). ZcF was significantly lower in athletes than in sedentary participants (73.3 ± 19.2 vs. 93.4 ± 19.0 dyn·s/cm5, P = 0.005). Furthermore, ZcF was negatively correlated with SV (r = -0.691) and aortic compliance (r = -0.601) but was positively correlated with double product (r = 0.445) and TPR (r = 0.458; all P < 0.05). Multivariate analysis revealed that ZcF was the strongest predictor of SV followed by TPR and HR (adjusted R2 = 0.788, P < 0.001). Therefore, our findings collectively suggest that LV afterload quantified by aortic ZcF is significantly lower in endurance athletes than in sedentary adults. The lower pulsatile LV afterload may contribute to greater SV in endurance athletes.NEW & NOTEWORTHY This is the first study to investigate aortic impedance with the noninvasive, simultaneous recordings of aortic pressure using SphygmoCor XCEL and flow using phase-contrast MRI. We found that the characteristic impedance (Zc) is significantly lower in endurance athletes than sedentary adults, is the strongest predictor of stroke volume (SV), and is inversely associated with aortic compliance. These findings suggest that aortic impedance is a key determinant of the ventricular-vascular coupling adapted to long-term training in endurance athletes.

Quantitative MRI of Gd-DOTA Accumulation in the Mouse Brain After Intraperitoneal Administration: Validation by Mass Spectrometry.

In mice, intraperitoneal (ip) contrast agent (CA) administration is convenient for mapping microvascular parameters over a long-time window. However, continuous quantitative MRI of CA accumulation in brain over hours is still missing. To validate a quantitative time-resolved MRI technique for mapping the CA kinetics in brain upon ip administration. Prospective, animal model. 25 C57Bl/6JRj mice underwent MRI. 7-T, gradient echo sequence. Gd-DOTA concentration was monitored by MRI (25 s/repetition) over 135 minutes with (N = 15) and without (N = 10) ip mannitol challenge (5 g/kg). After the final repetition, the brains were sampled to quantify gadolinium by mass spectrometry (MS). Upon manual brain segmentation, the average gadolinium concentration was compared with the MS quantification in transcardially perfused (N = 20) and unperfused (N = 5) mice. Precontrast T1-maps were acquired in 8 of 25 mice. One-tailed Spearman and Pearson correlation between gadolinium quantification by MRI and by MS, D'Agostino-Pearson test for normal distribution, Bland-Altman analysis to evaluate the agreement between MRI and MS. Significance was set at P-value <0.05. MRI showed that ip administered CA reached the blood compartment (>5 mM) within 10 minutes and accumulated continuously for 2 hours in cerebrospinal fluid (>1 mM) and in brain tissue. The MRI-derived concentration maps showed interindividual differences in CA accumulation (from 0.47 to 0.81 mM at 2 hours) with a consistent distribution resembling the pathways of the glymphatic system. The average in-vivo brain concentration 2 hours post-CA administration correlated significantly (r = 0.8206) with the brain gadolinium quantification by MS for N = 21 paired observations available. The presented experimental and imaging protocol may be convenient for monitoring the spatiotemporal pattern of CA uptake and clearance in the mouse brain over 2 hours. The quantification of the CA from the MRI signal in brain is corroborated by MS. N/A TECHNICAL EFFICACY: Stage 1.

Latent signal models: Learning compact representations of signal evolution for improved time‐resolved, multi‐contrast MRI

To improve time-resolved reconstructions by training auto-encoders to learn compact representations of Bloch-simulated signal evolution and inserting the decoder into the forward model. Building on model-based nonlinear and linear subspace techniques, we train auto-encoders on dictionaries of simulated signal evolution to learn compact, nonlinear, latent representations. The proposed latent signal model framework inserts the decoder portion of the auto-encoder into the forward model and directly reconstructs the latent representation. Latent signal models essentially serve as a proxy for fast and feasible differentiation through the Bloch equations used to simulate signal. This work performs experiments in the context of T2 -shuffling, gradient echo EPTI, and MPRAGE-shuffling. We compare how efficiently auto-encoders represent signal evolution in comparison to linear subspaces. Simulation and in vivo experiments then evaluate if reducing degrees of freedom by incorporating our proxy for the Bloch equations, the decoder portion of the auto-encoder, into the forward model improves reconstructions in comparison to subspace constraints. An auto-encoder with 1 real latent variable represents single-tissue fast spin echo, EPTI, and MPRAGE signal evolution to within 0.15% normalized RMS error, enabling reconstruction problems with 3 degrees of freedom per voxel (real latent variable + complex scaling) in comparison to linear models with 4-8 degrees of freedom per voxel. In simulated/in vivo T2 -shuffling and in vivo EPTI experiments, the proposed framework achieves consistent quantitative normalized RMS error improvement over linear approaches. From qualitative evaluation, the proposed approach yields images with reduced blurring and noise amplification in MPRAGE-shuffling experiments. Directly solving for nonlinear latent representations of signal evolution improves time-resolved MRI reconstructions.

Liver-ultrasound based motion modelling to estimate 4D dose distributions for lung tumours in scanned proton therapy

Motion mitigation strategies are crucial for scanned particle therapy of mobile tumours in order to prevent geometrical target miss and interplay effects. We developed a patient-specific respiratory motion model based on simultaneously acquired time-resolved volumetric MRI and 2D abdominal ultrasound images. We present its effects on 4D pencil beam scanned treatment planning and simulated dose distributions. Given an ultrasound image of the liver and the diaphragm, principal component analysis and Gaussian process regression were applied to infer dense motion information of the lungs. 4D dose calculations for scanned proton therapy were performed using the estimated and the corresponding ground truth respiratory motion; the differences were compared by dose difference volume metrics. We performed this simulation study on 10 combined CT and 4DMRI data sets where the motion characteristics were extracted from 5 healthy volunteers and fused with the anatomical CT data of two lung cancer patients. Median geometrical estimation errors below 2 mm for all data sets and maximum dose differences of = 43.2% and = 16.3% were found. Moreover, it was shown that abdominal ultrasound imaging allows to monitor organ drift. This study demonstrated the feasibility of the proposed ultrasound-based motion modelling approach for its application in scanned proton therapy of lung tumours.

Time-Resolved MRI Assessment of Convection-Enhanced Delivery by Targeted and Nontargeted Nanoparticles in a Human Glioblastoma Mouse Model.

Convection-enhanced delivery (CED) provides direct access of infusates to brain tumors; however, clinical translation of this technology has not been realized because of the inability to accurately visualize infusates in real-time and lack of targeting modalities against diffuse cancer cells. In this study, we use time-resolved MRI to reveal the kinetics of CED processes in a glioblastoma (GBM) model using iron oxide nanoparticles (NP) modified with a glioma-targeting ligand, chlorotoxin (CTX). Mice bearing orthotopic human GBM tumors were administered a single dose of targeted CTX-conjugated NP (NPCP-CTX) or nontargeted NP (NPCP) via CED. High-resolution T2-weighted, T2*-weighted, and quantitative T2 MRI were utilized to image NP delivery in real time and determined the volume of distribution (VD) of NPs at multiple time points over the first 48 hours post-CED. GBM-specific targeting was evaluated by flow cytometry and intracellular NP localization by histologic assessment. NPCP-CTX produced a VD of 121 ± 39 mm3 at 24 hours, a significant increase compared with NPCP, while exhibiting GBM specificity and localization to cell nuclei. Notably, CED of NPCP-CTX resulted in a sustained expansion of VD well after infusion, suggesting a possible active transport mechanism, which was further supported by the presence of NPs in endothelial and red blood cells. In summary, we show that time-resolved MRI is a suitable modality to study CED kinetics, and CTX-mediated CED facilitates extensive distribution of infusate and specific targeting of tumor cells. SIGNIFICANCE: MRI is used to monitor convection-enhanced delivery in real time using a nanoparticle-based contrast agent, and glioma-specific targeting significantly improves the volume of distribution in tumors.

Time-resolved volumetric MRI in MRI-guided radiotherapy: an in silico comparative analysis

MRI-treatment units enable 2D cine-MRI centred in the tumour for motion detection in radiotherapy, but they lack 3D information due to spatio-temporal limits. To derive time-resolved 3D information, different approaches have been proposed in the literature, but a rigorous comparison among these strategies has not yet been performed. The goal of this study is to quantitatively investigate five published strategies that derive time-resolved volumetric MRI in MRI-guided radiotherapy: Propagation, out-of-plane motion compensation, Fayad model, ROI-based model and Stemkens model. Comparisons were performed using an MRI digital phantom generated with six different patient-derived motion signals and tumour-shapes. An average 4D cycle was generated as well as 2D cine-MRI data with corresponding 3D in-room ground truth. Quantitative analysis was performed by comparing the estimated 3D volume to the ground truth available for each 2D cine-MRI sample. A grouped patient statistical analysis was performed to evaluate the performance of the selected methods, in case of tumour tracking or motion estimation of the whole anatomy. Analyses were also performed based on patient characteristics. Quantitative ranking of the investigated methods highlighted that Propagation and ROI-based model strategies achieved an overall median tumour centre of mass 3D distance from the ground truth of 1.1 mm and 1.3 mm, respectively, and a diaphragm distance below 1.6 mm. Higher errors and variabilities were instead obtained for other methods, which lack the ability to compensate for in-room variations and to account for regional changes. These results were especially evident when further analysing patient characteristics, where errors above 2 mm/5 mm in tumour/diaphragm were found for more irregular breathing patterns in case of out-of-plane motion compensation, Fayad and Stemkens models. These findings suggest the potential of the proposed in silico framework to develop and compare strategies to estimate time-resolved 3DMRI in MRI-guided radiotherapy.

Dynamic cerebellar herniation in Chiari patients during the cardiac cycle evaluated by dynamic magnetic resonance imaging.

Cerebellar herniation in Chiari patients can be dynamic, following the cerebrospinal fluid pulsatility during the cardiac cycle. We present a voxel intensity distribution method (VIDM) to automatically extract the pulsatility-dependent herniation in time-resolved MRI (CINE MRI) and compare it to the simple linear measurements. The degree of herniation is furthermore compared on CINE and static sequences, and the cerebellar movement is correlated to the presence of hydrocephalus and syringomyelia. The cerebellar movement in 27 Chiari patients is analyzed with VIDM and the results were compared to linear measurements on an image viewer (visual inspection, VI) using a paired t test. Second, an ANOVA test is applied to compare the degree of herniation on static 3D MRI and CINE. Finally, the Pearson's correlation coefficient is calculated for the correlation between cerebellar movement and the presence of hydrocephalus and syringomyelia. VIDM showed significant movement in 85% of our patients. Assuming that movement < 1mm cannot be detected reliably on an image viewer, VI identified movement in 29.6% of the patients (p = 0.002). The herniation was greater on static sequences than on CINE in most cases, but this was not statistically significant. The cerebellar movement was not correlated with hydrocephalus or syringomyelia (Pearson's coefficient < 0.3). VIDM is a sensitive method to detect tissue movement on CINE MRI and could be used for Chiari patients, but also for the evaluation of cyst membranes, ventriculostomies, etc. The cerebellar movement appears not to correlate with hydrocephalus and syringomyelia in Chiari patients.

OC-0409 Comparison of different strategies to derive time-resolved volumetric MRI in MRI-guided radiotherapy

OC-0409 Comparison of different strategies to derive time-resolved volumetric MRI in MRI-guided radiotherapy

Elevated diastolic wall shear stress in regurgitant semilunar valvular lesions.

Alterations in wall shear stress (WSS) assessed using 4D flow MRI have been shown to play a role in various vascular pathologies, such as bicuspid aortic valve aortopathy. Most studies have focused on systolic WSS, whereas altered diastolic hemodynamics in regurgitant semilunar valvular lesions have not so far been well characterized. To investigate diastolic WSS in aortic and pulmonary regurgitation. Retrospective data analysis. Thirty tetralogy of Fallot patients, 19 bicuspid aortic valve patients, 11 healthy volunteers. 5 T, 3D time-resolved phase-contrast MRI with 3D velocity encoding. Estimation of WSS and its axial and circumferential vector components along cardiac cycle timeframes in the proximal main pulmonary artery in pulmonary regurgitation (PR) and in the proximal ascending aorta in aortic regurgitation (AR) as well as in healthy volunteers. Wilcoxon matched pairs test was used for intra-group comparisons and Mann-Whitney test for intergroup comparisons. Correlations were assessed using Spearman correlation. WSS along the entire cardiac cycle was higher in PR and AR in comparison with controls (mean WSS 0.381 ± 0.070 vs. 0.220 ± 0.018, P < 0.0001; 0.361 ± 0.099 vs. 0.212 ± 0.030, P < 0.0001; respectively). Peak diastolic WSS was significantly higher than the mean WSS in AR and PR (P < 0.0001-0.005). The severity of PR correlated with the peak diastolic axial WSS (Spearman's r s = 0.454, P = 0.018), whereas the severity of AR correlated with both peak systolic and diastolic tangential WSS (Spearman's r s = 0.458, P = 0.049; r s = 0.539, P = 0.017, respectively). Elevated diastolic WSS is a component of the altered flow hemodynamics in AR and PR. This may give more insight into the pathophysiologic role of WSS in vascular remodeling in AR and PR. 4 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2019;50:763-770.

On the relationship between instantaneous phase synchrony and correlation-based sliding windows for time-resolved fMRI connectivity analysis

Correlation-based sliding window analysis (CSWA) is the most commonly used method to estimate time-resolved functional MRI (fMRI) connectivity. However, instantaneous phase synchrony analysis (IPSA) is gaining popularity mainly because it offers single time-point resolution of time-resolved fMRI connectivity. We aim to provide a systematic comparison between these two approaches, on temporal, topological and anatomical levels.For this purpose, we used resting-state fMRI data from two separate cohorts with different temporal resolutions (45 healthy subjects from Human Connectome Project fMRI data with repetition time of 0.72 s and 25 healthy subjects from a separate validation fMRI dataset with a repetition time of 3 s). For time-resolved functional connectivity analysis, we calculated tapered CSWA over a wide range of different window lengths that were compared to IPSA.We found a strong association in connectivity dynamics between IPSA and CSWA when considering the absolute values of CSWA. The association between CSWA and IPSA was stronger for a window length of ∼20 s (shorter than filtered fMRI wavelength) than ∼100 s (longer than filtered fMRI wavelength), irrespective of the sampling rate of the underlying fMRI data. Narrow-band filtering of fMRI data (0.03–0.07 Hz) yielded a stronger relationship between IPSA and CSWA than wider-band (0.01–0.1 Hz). On a topological level, time-averaged IPSA and CSWA nodes were non-linearly correlated for both short (∼20 s) and long (∼100 s) windows, mainly because nodes with strong negative correlations (CSWA) displayed high phase synchrony (IPSA). IPSA and CSWA were anatomically similar in the default mode network, sensory cortex, insula and cerebellum.Our results suggest that IPSA and CSWA provide comparable characterizations of time-resolved fMRI connectivity for appropriately chosen window lengths. Although IPSA requires narrow-band fMRI filtering, it does not mandate a (semi-)arbitrary choice of window length and window overlap. A code for calculating IPSA is provided.

Time-resolved contrast-enhanced MR angiography with single-echo Dixon fat suppression.

Dixon-based fat suppression has recently gained interest for dynamic contrast-enhanced MRI, but multi-echo techniques require longer scan times and reduce temporal resolution compared to single-echo alternatives without fat suppression. The purpose of this work is to demonstrate accelerated single-echo Dixon imaging with high spatial and temporal resolution. Real-valued water and fat images can be obtained from a single measurement if the shared initial phase and that due to ΔB0 are assumed known a priori. An expression for simultaneous sensitivity encoding (SENSE) unfolding and fat-water separation is derived for the general undersampling case, and simplified under the special case of uniform Cartesian undersampling. In vivo experiments were performed in extremities and brain with SENSE acceleration factors of up to R = 8. Single-echo Dixon reconstruction of highly undersampled data was successfully demonstrated. Dynamic contrast-enhanced water and fat images provided high spatial and temporal resolution dynamic images with image update times shorter than previous single-echo Dixon work. Time-resolved contrast-enhanced MRI with single-echo Dixon fat suppression shows high image quality, improved vessel delineation, and reduced sensitivity to motion when compared to time-subtraction methods.

Optimization of the Surgical Approach in AVMs Using MRI and 4D DSA Fusion Technique : Atechnical note.

Previously published data demonstrated the possibility of displaying the angioarchitecture of intracranial vascular malformations using time-resolved 3D imaging (4D digital subtraction angiography [DSA]). The purpose of our study was to prove the technical feasibility of creating fused images of time-resolved 3D reconstructions and MPRAGE MRI data sets and to check the reliability of the correct anatomical display of the angioma nidus and the venous drainage in the fused images of patients with intracranial arteriovenous malformations (AVM). In this study 20patients with intracranial AVM underwent pretherapeutic DSA and time-resolved 3D DSA in addition to MRI including MPRAGE sequences. The images were post-processed with the fusion software tool on adedicated research workstation. The fusion of both imaging modalities was done semi-automatically with automatic co-registration software followed by amanual co-registration. Co-registered DSA/MRI data sets of 20 untreated AVMs were evaluated independently by two reviewers. Image fusion was successful in all 20cases with an acceptable additional set-up time. The fused images were highly scored by the two raters in respect to their congruency of the dedicated regions. Precise anatomical localization of the nidus, the feeding arteries and the draining veins were possible with the merged images. Creating fused images of time-resolved 3D DSA and contrast-enhanced T1-weighted MPRAGE MR images might be beneficial for the preoperative and intrasurgical workflow in patients with AVMs. This new software tool fulfils the required quality and accuracy of the merged images. The clinical validation has to be proven in further studies.

Matrix pencil decomposition of time-resolved proton MRI for robust and improved assessment of pulmonary ventilation and perfusion.

To present an improved and robust method of pulmonary function assessment from time-resolved proton MRI using a matrix pencil (MP) method in combination with a linear least squares analysis. Simulations of the signal time course in lung parenchyma were performed to compare the accuracy of Fourier decomposition (FD) and MP methods for the estimation of respiratory and cardiac amplitudes. Series of two-dimensional time-resolved lung images were acquired in healthy volunteers at 1.5 T using ultra-fast steady-state free precession. Qualitative lung ventilation- and perfusion-weighted images as well as a quantitative map of fractional ventilation, perfusion, and blood arrival time were calculated using the proposed MP method and compared with the contemporary FD technique. A region-of-interest analysis was performed on the quantitative data. The signal analysis performed using MP decomposition resulted in reduced variability of the estimated respiratory and cardiac amplitudes in comparison with FD for both simulated and in vivo data. MP decomposition provides an automatic, robust, and more accurate estimation of amplitudes of respiratory and cardiac signal modulations in the lung parenchyma than the contemporary FD technique. Magn Reson Med 77:336-342, 2017. © 2016 Wiley Periodicals, Inc.

Assessment of intracardiac flow and vorticity in the right heart of patients after repair of tetralogy of Fallot by flow-sensitive 4D MRI.

To comprehensively and quantitatively analyse flow and vorticity in the right heart of patients after repair of tetralogy of Fallot (rTOF) compared with healthy volunteers. Time-resolved flow-sensitive 4D MRI was acquired in 24 rTOF patients and 12 volunteers. Qualitative flow evaluation was based on consensus reading of two observers. Quantitative analysis included segmentation of the right atrium (RA) and ventricle (RV) in a four-chamber view to extract volumes and regional haemodynamic information for computation of regional mean and peak vorticity. Right heart intra-atrial, intraventricular and outflow tract flow patterns differed considerably between rTOF patients and volunteers. Peak RA and mean RV vorticity was significantly higher in patients (p = 0.02/0.05). Significant negative correlations were found between patients' maximum and mean RV and RA vorticity and ventricular volumes (p < 0.05). The main pulmonary artery (MPA) regurgitant flow was associated with higher RA and RV vorticity, which was significant for RA maximum and RV mean vorticity (p = 0.01/0.03). The calculation of vorticity based on 4D flow data is an alternative approach to assess intracardiac flow changes in rTOF patients compared with qualitative flow visualization. Alterations in intracardiac vorticity could be relevant with regard to the development of RV dilation and impaired function. • 4D flow MRI with vorticity calculation enables a novel approach to assess intracardiac flow. • Significantly higher intracardiac vorticity occurred in patients after repair of tetralogy of Fallot. • Regurgitant flow in the main pulmonary artery is associated with higher right heart vorticity.

Improved assessment of aortic 3D blood flow with combined k-t accelerated 3D CINE bSSFP &amp; 4D flow MRI

Background To determine hemodynamic parameters from 4D flow MRI, vessel boundaries are typically depicted by calculating a 3D phase contrast MR angiogram (PC-MRA) from 4D flow magnitude and flow images. However, this approach is limited by 1) low blood-tissue contrast of the magnitude images and 2) velocity weighting of PC-MRA may not fully depict areas of slow or swirling (vortex) flow. As a result, 3D flow analysis at the aortic root has been challenging due to the presence of vortex flow distal to the aortic valve. Fast 3D balanced steady state free precession (bSSFP) time-resolved MRI sequences offer a promising alternative, providing anatomical images with improved blood-tissue contrast. The goal of this feasibility study was to integrate results from 3D CINE bSSFP with 4D flow MRI to improve the segmentation of vessel anatomy and thus 3D blood flow assessment compared to standard segmentation using PC-MRA data alone.

Velocity Measurement in Carotid Artery: Quantitative Comparison of Time-Resolved 3D Phase-Contrast MRI and Image-based Computational Fluid Dynamics.

Background:Understanding hemodynamic environment in vessels is important for realizing the mechanisms leading to vascular pathologies.Objectives:Three-dimensional velocity vector field in carotid bifurcation is visualized using TR 3D phase-contrast magnetic resonance imaging (TR 3D PC MRI) and computational fluid dynamics (CFD). This study aimed to present a qualitative and quantitative comparison of the velocity vector field obtained by each technique.Subjects and Methods:MR imaging was performed on a 30-year old male normal subject. TR 3D PC MRI was performed on a 3 T scanner to measure velocity in carotid bifurcation. 3D anatomical model for CFD was created using images obtained from time-of-flight MR angiography. Velocity vector field in carotid bifurcation was predicted using CFD and PC MRI techniques. A statistical analysis was performed to assess the agreement between the two methods.Results:Although the main flow patterns were the same for the both techniques, CFD showed a greater resolution in mapping the secondary and circulating flows. Overall root mean square (RMS) errors for all the corresponding data points in PC MRI and CFD were 14.27% in peak systole and 12.91% in end diastole relative to maximum velocity measured at each cardiac phase. Bland-Altman plots showed a very good agreement between the two techniques. However, this study was not aimed to validate any of methods, instead, the consistency was assessed to accentuate the similarities and differences between Time-resolved PC MRI and CFD.Conclusion:Both techniques provided quantitatively consistent results of in vivo velocity vector fields in right internal carotid artery (RCA). PC MRI represented a good estimation of main flow patterns inside the vasculature, which seems to be acceptable for clinical use. However, limitations of each technique should be considered while interpreting results.

Velocity Measurement in Carotid Artery: Quantitative Comparison of Time-Resolved 3D Phase-Contrast MRI and Image-based Computational Fluid Dynamics

Background: Understanding hemodynamic environment in vessels is important for realizing the mechanisms leading to vascular pathologies. Objectives: Three-dimensional velocity vector field in carotid bifurcation is visualized using TR 3D phase-contrast magnetic resonance imaging (TR 3D PC MRI) and computational fluid dynamics (CFD). This study aimed to present a qualitative and quantitative comparison of the velocity vector field obtained by each technique. Subjects and Methods: MR imaging was performed on a 30-year old male normal subject. TR 3D PC MRI was performed on a 3 T scanner to measure velocity in carotid bifurcation. 3D anatomical model for CFD was created using images obtained from time-of-flight MR angiography. Velocity vector field in carotid bifurcation was predicted using CFD and PC MRI techniques. A statistical analysis was performed to assess the agreement between the two methods. Results: Although the main flow patterns were the same for the both techniques, CFD showed a greater resolution in mapping the secondary and circulating flows. Overall root mean square (RMS) errors for all the corresponding data points in PC MRI and CFD were 14.27% in peak systole and 12.91% in end diastole relative to maximum velocity measured at each cardiac phase. Bland-Altman plots showed a very good agreement between the two techniques. However, this study was not aimed to validate any of methods, instead, the consistency was assessed to accentuate the similarities and differences between Time-resolved PC MRI and CFD. Conclusion: Both techniques provided quantitatively consistent results of in vivo velocity vector fields in right internal carotid artery (RCA). PC MRI represented a good estimation of main flow patterns inside the vasculature, which seems to be acceptable for clinical use. However, limitations of each technique should be considered while interpreting results.

A Small Surrogate for the Golden Angle in Time-Resolved Radial MRI Based on Generalized Fibonacci Sequences.

In golden angle radial magnetic resonance imaging a constant azimuthal radial profile spacing of 111.246...(°) guarantees a nearly uniform azimuthal profile distribution in k-space for an arbitrary number of radial profiles. Even though this profile order is advantageous for various real-time imaging methods, in combination with balanced steady-state free precession (SSFP) sequences the large azimuthal angle increment may lead to strong image artifacts, due to the varying eddy currents introduced by the rapidly switching gradient scheme. Based on a generalized Fibonacci sequence, a new sequence of smaller irrational angles is introduced ( 49.750...(°), 32.039...(°), 27.198...(°), 23.628...(°), ... ). The subsequent profile orders guarantee the same sampling efficiency as the golden angle if at least a minimum number of radial profiles is used for reconstruction. The suggested angular increments are applied for dynamic imaging of the heart and the temporomandibular joint. It is shown that for balanced SSFP sequences, trajectories using the smaller golden angle surrogates strongly reduce the image artifacts, while the free retrospective choice of the reconstruction window width is maintained.

Diffusion of Drug Delivery Nanoparticles into Biogels Using Time-Resolved MicroMRI.

Nanoparticle-based therapeutic agents can in some cases provide selective delivery to tumors, yet this field would greatly benefit from more detailed understanding of particle transport into and within tumor tissue. To provide fundamental information for optimizing interstitial transport of polymeric nanoparticles, we have developed a quantitative approach employing real-time analysis of nanoparticle diffusion into bulk biological hydrogels using microMRI. We use two distinct imaging approaches to probe the migration of two novel "theranostic" polymeric agents (combining drug delivery and contrast agent functions) into bulk hydrogels. Theranostic agent diffusion measured using time-resolved MRI agrees well with diffusion measured for simple probe particles using fluorescence spectroscopies. Furthermore, compared with established fluorescence techniques, which are restricted by sample thickness, our approach provides a three-dimensional diffusion rate and concentration distribution of nanoparticles over macroscopic distances in biological media. These results carry implications for in vivo tracking of theranostic nanoparticles into tumor interstitium.