Search

Abstract For 4D Flow MRI of mean and turbulent flow a compromise between spatiotemporal undersampling and velocity encodings needs to be found. Assuming a fixed scan time budget, the impact of trading off spatiotemporal undersampling versus velocity encodings on quantification of velocity and turbulence for aortic 4D Flow MRI was investigated. For this purpose, patient-specific mean and turbulent aortic flow data were generated using computational fluid dynamics which were embedded into the patient-specific background image data to generate synthetic MRI data with corresponding ground truth flow. Cardiac and respiratory motion were included. Using the synthetic MRI data as input, 4D Flow MRI was subsequently simulated with undersampling along pseudo-spiral Golden angle Cartesian trajectories for various velocity encoding schemes. Data were reconstructed using a locally low rank approach to obtain mean and turbulent flow fields to be compared to ground truth. Results show that, for a 15-min scan, velocity magnitudes can be reconstructed with good accuracy relatively independent of the velocity encoding scheme ( $$SSI{M}_{U}=0.938\pm 0.003)$$ S S I M U = 0.938 ± 0.003 ) , good accuracy ( $$SSI{M}_{U}\ge 0.933$$ S S I M U ≥ 0.933 ) and with peak velocity errors limited to 10%. Turbulence maps on the other hand suffer from both lower reconstruction quality ( $$SSI{M}_{TKE}\ge 0.323$$ S S I M TKE ≥ 0.323 ) and larger sensitivity to undersampling, motion and velocity encoding strengths ( $$SSI{M}_{TKE}=0.570\pm 0.110)$$ S S I M TKE = 0.570 ± 0.110 ) when compared to velocity maps. The best compromise to measure unwrapped velocity maps and turbulent kinetic energy given a fixed 15-min scan budget was found to be a 7-point multi- $${V}_{enc}$$ V enc acquisition with a low $${V}_{enc}$$ V enc tuned for best sensitivity to the range of expected intra-voxel standard deviations and a high $${V}_{enc}$$ V enc larger than the expected peak velocity.

Abstract Background Standardised performance assessment of image acquisition, reconstruction and processing methods is limited by the absence of images paired with ground truth reference values. To this end, we propose MRXCAT2.0 to generate synthetic data, covering healthy and pathological function, using a biophysical model. We exemplify the approach by generating cardiovascular magnetic resonance (CMR) images of healthy, infarcted, dilated and hypertrophic left-ventricular (LV) function. Method In MRXCAT2.0, the XCAT torso phantom is coupled with a statistical shape model, describing population (patho)physiological variability, and a biophysical model, providing known and detailed functional ground truth of LV morphology and function. CMR balanced steady-state free precession images are generated using MRXCAT2.0 while realistic image appearance is ensured by assigning texturized tissue properties to the phantom labels. Finding Paired CMR image and ground truth data of LV function were generated with a range of LV masses (85–140 g), ejection fractions (34–51%) and peak radial and circumferential strains (0.45 to 0.95 and − 0.18 to − 0.13, respectively). These ranges cover healthy and pathological cases, including infarction, dilated and hypertrophic cardiomyopathy. The generation of the anatomy takes a few seconds and it improves on current state-of-the-art models where the pathological representation is not explicitly addressed. For the full simulation framework, the biophysical models require approximately two hours, while image generation requires a few minutes per slice. Conclusion MRXCAT2.0 offers synthesis of realistic images embedding population-based anatomical and functional variability and associated ground truth parameters to facilitate a standardized assessment of CMR acquisition, reconstruction and processing methods.

Abstract We propose to synthesize patient-specific 4D flow MRI datasets of turbulent flow paired with ground truth flow data to support training of inference methods. Turbulent blood flow is computed based on the Navier–Stokes equations with moving domains using realistic boundary conditions for aortic shapes, wall displacements and inlet velocities obtained from patient data. From the simulated flow, synthetic multipoint 4D flow MRI data is generated with user-defined spatiotemporal resolutions and reconstructed with a Bayesian approach to compute time-varying velocity and turbulence maps. For MRI data synthesis, a fixed hypothetical scan time budget is assumed and accordingly, changes to spatial resolution and time averaging result in corresponding scaling of signal-to-noise ratios (SNR). In this work, we focused on aortic stenotic flow and quantification of turbulent kinetic energy (TKE). Our results show that for spatial resolutions of 1.5 and 2.5 mm and time averaging of 5 ms as encountered in 4D flow MRI in practice, peak total turbulent kinetic energy downstream of a 50, 75 and 90% stenosis is overestimated by as much as 23, 15 and 14% (1.5 mm) and 38, 24 and 23% (2.5 mm), demonstrating the importance of paired ground truth and 4D flow MRI data for assessing accuracy and precision of turbulent flow inference using 4D flow MRI exams.

To improve risk stratification in extracranial internal carotid artery disease (CAD), patients who would benefit maximally from revascularization must be identified. In cardiology, the fractional flow reserve (FFR) has become a reference standard for evaluating the functional severity of coronary artery stenosis, and noninvasive surrogates thereof relying on computational fluid dynamics (CFD) have been developed. Here, we present a CFD-based workflow using digital twins of patients’ carotid bifurcations derived from computed tomography angiography for the noninvasive functional assessment of CAD. We reconstructed patient-specific digital twins of 37 carotid bifurcations. We implemented a CFD model using common carotid artery peak systolic velocity (PSV) acquired with Doppler ultrasound (DUS) as inlet boundary condition and a two-element Windkessel model as oulet boundary condition. The agreement between CFD and DUS on the PSV in the internal carotid artery (ICA) was then compared. The relative error for the agreement between DUS and CFD was 9% ± 20% and the intraclass correlation coefficient was 0.88. Furthermore, hyperemic simulations in a physiological range were feasible and unmasked markedly different pressure drops along two ICA stenoses with similar degree of narrowing under comparable ICA blood flow. Hereby, we lay the foundation for prospective studies on noninvasive CFD-based derivation of metrics similar to the FFR for the assessment of CAD.

Objectives: Paravalvular leakage (PVL) and eccentric aortic regurgitation remain a major clinical concern in patients receiving transcatheter aortic valve replacement (TAVR), and regurgitant volume remains the main readout parameter in clinical assessment. In this work we investigate the effect of jet origin and trajectory of mild aortic regurgitation on left ventricular hemodynamics in a porcine model.Methods: A pig model of mild aortic regurgitation/PVL was established by transcatheter piercing and dilating the non-coronary (NCC) or right coronary cusp (RCC) of the aortic valve close to the valve annulus. The interaction between regurgitant blood and LV hemodynamics was assessed by 4D flow cardiovascular MRI.Results: Six RCC, six NCC, and two control animals were included in the study and with one dropout in the NCC group, the success rate of model creation was 93%. Regurgitant jets originating from NCC were directed along the ventricular side of the anterior mitral leaflet and integrated well into the diastolic vortex forming in the left ventricular outflow tract. However, jets from the RCC were orientated along the septum colliding with flow within the vortex, and progressing down to the apex. As a consequence, the presence as well as the area of the vortex was reduced at the site of impact compared to the NCC group. Impairment of vortex formation was localized to the area of impact and not the entire vortex ring. Blood from the NCC jet was largely ejected during the following systole, whereas ejection of large portion of RCC blood was protracted.Conclusions: Even for mild regurgitation, origin and trajectory of the regurgitant jet does cause a different effect on LV hemodynamics. Septaly oriented jets originating from RCC collide with the diastolic vortex, reduce its size, and reach the apical region of the left ventricle where blood resides extendedly. Hence, RCC jets display hemodynamic features which may have a potential negative impact on the long-term burden to the heart.

Phenylketonuria (PKU) is an inherited metabolic disease characterized by abnormally high concentrations of the essential amino acid L-phenylalanine (Phe) in blood plasma caused by reduced activity of phenylalanine hydroxylase (PAH). While numerous studies have shown association between high plasma Phe concentration and intellectual impairment, it is not clear whether increased Phe fluctuations also observed in PKU affect the brain as well. To investigate this, time-resolved in vivo data on Phe and competing large neutral amino acid (LNAA) concentrations in neurons are needed, but cannot be acquired readily with current methods. We have used in silico modeling as an alternative approach to characterize the interactive dynamics of Phe and competing LNAAs (CL) in the neurovascular unit (NVU). Our results suggest that plasma Phe fluctuations can propagate into the NVU cells and change there the concentration of LNAAs, with the highest magnitude of this effect observed at low frequency and high amplitude-to-mean ratio of the plasma Phe concentration fluctuations. Our model further elucidates the effect of therapeutic LNAA supplementation in PKU, showing how abnormal concentrations of Phe and CL in the NVU move thereby toward normal physiologic levels.

The homeostatic regulation of large neutral amino acid (LNAA) concentration in the brain interstitial fluid (ISF) is essential for proper brain function. LNAA passage into the brain is primarily mediated by the complex and dynamic interactions between various solute carrier (SLC) transporters expressed in the neurovascular unit (NVU), among which SLC7A5/LAT1 is considered to be the major contributor in microvascular brain endothelial cells (MBEC). The LAT1-mediated trans-endothelial transport of LNAAs, however, could not be characterized precisely by available in vitro and in vivo standard methods so far. To circumvent these limitations, we have incorporated published in vivo data of rat brain into a robust computational model of NVU-LNAA homeostasis, allowing us to evaluate hypotheses concerning LAT1-mediated trans-endothelial transport of LNAAs across the blood brain barrier (BBB). We show that accounting for functional polarity of MBECs with either asymmetric LAT1 distribution between membranes and/or intrinsic LAT1 asymmetry with low intraendothelial binding affinity is required to reproduce the experimentally measured brain ISF response to intraperitoneal (IP) L-tyrosine and L-phenylalanine injection. On the basis of these findings, we have also investigated the effect of IP administrated L-tyrosine and L-phenylalanine on the dynamics of LNAAs in MBECs, astrocytes and neurons. Finally, the computational model was shown to explain the trans-stimulation of LNAA uptake across the BBB observed upon ISF perfusion with a competitive LAT1 inhibitor.

Purpose: To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction. Methods: CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided. Results: The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain. Conclusion: CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research. [ABSTRACT FROM AUTHOR]