A kinematic model-based analysis framework for 3D Cine-DENSE-validation with an axially compressed gel phantom and application in sheep before and after antero-apical myocardial infarction.

PURPOSE: Myocardial strain is increasingly used to assess left ventricular (LV) function. Incorporation of LV deformation into finite element (FE) modeling environment with subsequent strain calculation will allow analysis to reach its full potential. We describe a new kinematic model-based analysis framework (KMAF) to calculate strain from 3D cine-DENSE (displacement encoding with stimulated echoes) MRI. METHODS: Cine-DENSE allows measurement of 3D myocardial displacement with high spatial accuracy. The KMAF framework uses cine cardiovascular magnetic resonance (CMR) to facilitate cine-DENSE segmentation, interpolates cine-DENSE displacement, and kinematically deforms an FE model to calculate strain. This framework was validated in an axially compressed gel phantom and applied in 10 healthy sheep and 5 sheep after myocardial infarction (MI). RESULTS: Excellent Bland-Altman agreement of peak circumferential (Ecc ) and longitudinal (Ell ) strain (mean difference = 0.021 ± 0.04 and -0.006 ± 0.03, respectively), was found between KMAF estimates and idealized FE simulation. Err had a mean difference of -0.014 but larger variation (±0.12). Cine-DENSE estimated end-systolic (ES) Ecc , Ell and Err exhibited significant spatial variation for healthy sheep. Displacement magnitude was reduced on average by 27%, 42%, and 56% after MI in the remote, adjacent and MI regions, respectively. CONCLUSIONS: The KMAF framework allows accurate calculation of 3D LV Ecc and Ell from cine-DENSE.

Identification of intra-individual variation in intracranial arterial flow by MRI and the effect on computed hemodynamic descriptors.

OBJECTIVES: To determine the intra-individual flow variation in serially acquired studies, and the influence of this variation on subsequent hemodynamic simulations using the inlet flow as a boundary condition. Author: Kindly check and confirm whether the corresponding authors are correctly identified.Confirmed. MATERIALS AND METHODS: This prospective study included 51 patients (37 females and 14 males) with unruptured intracranial aneurysms who have received more than three times follow-up of 2D phase-contrast MR. The flow and velocity parameters were extracted to calculate the reproducibility and variation. Patient-specific computational fluid dynamics simulations were performed using the measured flows. RESULTS: Intraclass correlation coefficients for mean and maximum velocity and flow parameters ranged from 0.77 to 0.90. A 10% CV of mean flow was identified. Variations of 10% in inlet flow resulted in hemodynamic changes including 41.41% of peak systolic wall shear stress; 39.13% of end-diastolic wall shear stress; 2.79% of low shear area at peak systole; 2.12% of low shear area at end diastole: 47.57% of time-averaged wall shear stress; and 0.17% of oscillatory shear index. CONCLUSION: This study identified 10% of intra-individual mean flow variation on phase-contrast MR. Intra-individual flow variation resulted in a non-negligible variation in wall shear stress, but relatively small variation in low shear area in hemodynamic calculations.

Computer-aided quantification of non-contrast 3D black blood MRI as an efficient alternative to reference standard manual CT angiography measurements of abdominal aortic aneurysms.

BACKGROUND: Non-contrast 3D black blood MRI is a promising tool for abdominal aortic aneurysm (AAA) surveillance, permitting accurate aneurysm diameter measurements needed for patient management. PURPOSE: To evaluate whether automated AAA volume and diameter measurements obtained from computer-aided segmentation of non-contrast 3D black blood MRI are accurate, and whether they can supplant reference standard manual measurements from contrast-enhanced CT angiography (CTA). MATERIALS AND METHODS: Thirty AAA patients (mean age, 71.9 ± 7.9 years) were recruited between 2014 and 2017. Participants underwent both non-contrast black blood MRI and CTA within 3 months of each other. Semi-automatic (computer-aided) MRI and CTA segmentations utilizing deformable registration methods were compared against manual segmentations of the same modality using the Dice similarity coefficient (DSC). AAA lumen and total aneurysm volumes and AAA maximum diameter, quantified automatically from these segmentations, were compared against manual measurements using Pearson correlation and Bland-Altman analyses. Finally, automated measurements from non-contrast 3D black blood MRI were evaluated against manual CTA measurements using the Wilcoxon test, Pearson correlation and Bland-Altman analyses. RESULTS: Semi-automatic segmentations had excellent agreement with manual segmentations (lumen DSC: 0.91 ± 0.03 and 0.94 ± 0.03; total aneurysm DSC: 0.92 ± 0.02 and 0.94 ± 0.03, for black blood MRI and CTA, respectively). Automated volume and maximum diameter measurements also had excellent correlation to their manual counterparts for both black blood MRI (volume: r = 0.99, P < 0.001; diameter: r = 0.97, P < 0.001) and CTA (volume: r = 0.99, P < 0.001; diameter: r = 0.97, P < 0.001). Compared to manual CTA measurements, bias and limits of agreement (LOA) for automated MRI measurements (lumen volume: 1.49, [-4.19 7.17] cm3; outer wall volume: -2.46, [-14.05 9.13] cm3; maximal diameter: 0.08, [-6.51 6.67] mm) were largely equivalent to those of manual MRI measurements, particularly for maximum AAA diameter (lumen volume: 0.73, [-6.47 7.93] cm3; outer wall volume: 0.98, [-10.54 12.5] cm3; maximal diameter: 0.08, [-3.67 3.83] mm). CONCLUSION: Semi-automatic segmentation of non-contrast 3D black blood MRI efficiently provides reproducible morphologic AAA assessment yielding accurate AAA diameters and volumes with no clinically relevant differences compared to either automatic or manual measurements based on CTA.

Clinical validation of three cardiovascular magnetic resonance techniques to measure strain and torsion in patients with suspected coronary artery disease.

BACKGROUND: Several cardiovascular magnetic resonance (CMR) techniques can measure myocardial strain and torsion with high accuracy. The purpose of this study was to compare displacement encoding with stimulated echoes (DENSE), tagging and feature tracking (FT) for measuring circumferential and radial myocardial strain and myocardial torsion in order to assess myocardial function and infarct scar burden both at a global and at a segmental level. METHOD: 116 patients with a high likelihood of coronary artery disease (European SCORE > 15%) underwent CMR examination including cine images, tagging, DENSE and late gadolinium enhancement (LGE) in the short axis direction. In total, 97 patients had signs of myocardial disease and 19 had no abnormalities in terms of left ventricular (LV) wall mass index, LV ejection fraction, wall motion, LGE or a history of myocardial infarction. Thirty-four patients had myocardial infarct scar with a transmural LGE extent (transmurality) that exceeded 50% of the wall thickness in at least one segment. Global circumferential strain (GCS) and global radial strain (GRS) was analyzed using FT of cine loops, deformation of tag lines or DENSE displacement. RESULTS: DENSE and tagging both showed high sensitivity (82% and 71%) at a specificity of 80% for the detection of segments with > 50% LGE transmurality, and receiver operating characteristics (ROC) analysis showed significantly higher area under the curve-values (AUC) for DENSE (0.87) than for tagging (0.83, p < 0.001) and FT (0.66, p = 0.003). GCS correlated with global LGE when determined with DENSE (r = 0.41), tagging (r = 0.37) and FT (r = 0.15). GRS had a low but significant negative correlation with LGE; DENSE r = - 0.10, FT r = - 0.07 and tagging r = - 0.16. Torsion from DENSE and tagging had a weak correlation (- 0.20 and - 0.22 respectively) with global LGE. CONCLUSION: Circumferential strain from DENSE detected segments with > 50% scar with a higher AUC than strain determined from tagging and FT at a segmental level. GCS and torsion computed from DENSE and tagging showed similar correlation with global scar size, while when computed from FT, the correlation was lower.

Left ventricular geometry during unloading and the end-systolic pressure volume relationship: Measurement with a modified real-time MRI-based method in normal sheep.

The left ventricular (LV) end-systolic (ES) pressure volume relationship (ESPVR) is the cornerstone of systolic LV function analysis. We describe a 2D real-time (RT) MRI-based method (RTPVR) with separate software tools for 1) semi-automatic level set-based shape prior method (LSSPM) of the LV, 2) generation of synchronized pressure area loops and 3) calculation of the ESPVR. We used the RTPVR method to measure ventricular geometry, ES pressure area relationship (ESPAR) and ESPVR during vena cava occlusion (VCO) in normal sheep. 14 adult sheep were anesthetized and underwent measurement of LV systolic function. Ten of the 14 sheep underwent RTMRI and eight of the 14 underwent measurement with conductance catheter; 4 had both RTMRI and conductance measurements. 2D cross sectional RTMRI were performed at apex, mid-ventricle and base levels during separate VCOs. The Dice similarity coefficient was used to compare LSSPM and manual image segmentation and thus determine LSSPM accuracy. LV cross-sectional area, major and minor axis length, axis ratio, major axis orientation angle and ESPAR were measured at each LV level. ESPVR was calculated with a trapezoidal rule. The Dice similarity coefficient between LSSPM and manual segmentation by two readers was 87.31±2.51% and 88.13±3.43%. All cross sections became more elliptical during VCO. The major axis orientation shifted during VCO but remained in the septo-lateral direction. LV chamber obliteration at the apical level occurred during VCO in 7 of 10 sheep that underwent RTMRI. ESPAR was non-linear at all levels. Finally, ESPVR was non-linear because of apical collapse. ESPVR measured by conductance catheter (EES,Index = 2.23±0.66 mmHg/ml/m2) and RT (EES,Index = 2.31±0.31 mmHg/ml/m2) was not significantly different. LSSPM segmentation of 2D RT MRI images is accurate and allows calculation of LV geometry, ESPAR and ESPVR during VCO. In the future, RTPVR will facilitate determination of regional systolic material parameters underlying ESPVR.

Myocardial injection of a thermoresponsive hydrogel with reactive oxygen species scavenger properties improves border zone contractility.

The decrease in contractility in myocardium adjacent (border zone; BZ) to a myocardial infarction (MI) is correlated with an increase in reactive oxygen species (ROS). We hypothesized that injection of a thermoresponsive hydrogel, with ROS scavenging properties, into the MI would decrease ROS and improve BZ function. Fourteen sheep underwent antero-apical MI. Seven sheep had a comb-like copolymer synthesized from N-isopropyl acrylamide (NIPAAm) and 1500 MW methoxy poly(ethylene glycol) methacrylate, (NIPAAm-PEG1500), injected (20 × 0.5 mL) into the MI zone 40 min after MI (MI + NIPAAm-PEG1500) and seven sheep were MI controls. Cardiac MRI was performed 2 weeks before and 6 weeks after MI + NIPAAm-PEG1500. BZ wall thickness at end systole was significantly higher for MI + NIPAAm-PEG1500 (12.32 ± 0.51 mm/m2 MI + NIPAAm-PEG1500 vs. 9.88 ± 0.30 MI; p = .023). Demembranated muscle force development for BZ myocardium 6 weeks after MI was significantly higher for MI + NIPAAm-PEG1500 (67.67 ± 2.61 mN/m2 MI + NIPAAm-PEG1500 vs. 40.53 ± 1.04 MI; p < .0001) but not significantly different from remote myocardium or BZ or non-operated controls. Levels of ROS in BZ tissue were significantly lower in the MI + NIPAAm-PEG1500 treatment group (hydroxyl p = .0031; superoxide p = .0182). We conclude that infarct injection of the NIPAAm-PEG1500 hydrogel with ROS scavenging properties decreased ROS and improved contractile protein function in the border zone 6 weeks after MI.

Quantitative measurement of atheroma burden: reproducibility in serial studies of atherosclerotic femoral arteries.

OBJECTIVES: This study aims to evaluate the reproducibility of measures of plaque morphology in serially acquired black-blood MRI of untreated atherosclerotic femoral arteries. MATERIALS AND METHODS: MR studies was obtained from 42 timepoints, on 12 patients with known femoral artery atherosclerosis. Images with a 3D isotropic FLASH with DANTE-prepared black blood contrast (DASH) at a 3-T scanner were acquired at baseline, within 1 week, and at 1 month. Six of the patients were scanned additionally at 6 months. Inter-scan and inter-observer variations of arterial area/volume measurements were evaluated. RESULTS: Measurement of vessel area, lumen area, wall area and wall volume showed inter-scan intraclass correlation coefficients (ICC) ranging from 0.92 to 0.97 for 3 scans, 0.91-0.97 for 4 scans, and inter-observer ICCs of 0.89-0.96. Among 3 scans, the coefficients of variance (CV) for the vessel area, lumen area, wall area and wall volume were 4.1%, 6.5%, 7.5%, and 4.4%. CVs among 4 scans ranged from 4.4% to 7.9%, and interobserver CVs ranged from 6.1% to 11.8% for the different area/volume measurements. CONCLUSION: DASH MRI is useful for quantifying atherosclerotic vessel area and volume of femoral arteries with low variability among serial repeated scans and between observers.

Sound Measurement in Patient-Specific 3D Printed Bench Models of Venous Pulsatile Tinnitus.

HYPOTHESIS: We hypothesize patient-specific flow models to be an adequate in vitro surrogate to allow for characterization of pulsatile tinnitus (PT) that affects three to five million Americans. BACKGROUND: PT, rhythmic sounds without an extracorporeal source that patients appreciate, can be caused by aberrant blood flow in large cerebral veins near the cochlea. To investigate the sound production mechanism, we created 3D printed flow models based on patient-specific cerebral venous anatomies. METHODS: Magnetic resonance angiography datasets from two patients with PT were used to generate patient-specific 3D printed flow models. A flow circuit connecting the patient-specific models to a pulsatile, continuous flow pump simulating cardiac cycle was created. Sound recordings were made along the surface of the models using an electronic stethoscope. Peak-to-rms amplitude, and area under the power spectral density (PSD) curve values were computed to evaluate the sound measurements. Wilcoxon rank sum test was used to statistically determine the differences in measurements between the patient-specific models. RESULTS: In patient-1, the recordings (peak-to-rms) from the internal jugular vein stenosis of baseline model (4.29 ±â€Š1.26 for 146 samples) were significantly louder (p < 0.001) than that of the altered model (3.29 ±â€Š0.96 for 143 samples). In patient-2, the sound measured at the transverse sinus stenosis in the pre-lumbar puncture model (4.84 ±â€Š1.11 for 148 samples) was significantly louder (p < 0.0001) than that of the post-lumbar puncture model (3.14 ±â€Š0.87 for 135 samples). CONCLUSIONS: The models are able to generate sounds very similar to those appreciated by patients and examiners in the cases of objective PT.

Shape-appearance constrained segmentation and separation of vein and artery in pulsatile tinnitus patients based on MR angiography and flow MRI.

This study reports on the development and evaluation of a novel segmentation method for extracting the internal jugular vein and the adjacent carotid artery from magnetic resonance (MR) images of patients with pulsatile tinnitus. A narrow band level set method with combined shape and appearance constraints was developed and applied to high-resolution MR images from 17 pulsatile tinnitus patients (age 52 ±â€¯23 years, 10 females). The proposed method was validated by comparing with the manual segmentation as well as by identifying the jugular vein and carotid artery based on 4D flow MRI in which the two types of vessels have opposing flow. Our study showed that the vein and artery are in contact with each other on 30.2% of all the slices. Dice value, Peak signal-to-noise ratio (PSNR), Hausdorff distance and mean sum of square distance (MSSD) between automatic and manual segmentation were 89.13 ±â€¯2.84%, 27.36 ±â€¯2.39%, 17.2 ±â€¯6.9 mm, 7.4 ±â€¯5.5 mm, demonstrating good segmentation accuracy. The average Dice similarity coefficient and the coefficient of variation compared with 4D flow MRI was 91.42 ±â€¯1.63% and 89.28 ±â€¯4.54% for the internal jugular vein and the carotid artery. The present pipeline for automatic internal jugular vein quantification holds promise for efficient image interpretation in large-scale cohort studies.

Visualizing wall enhancement over time in unruptured intracranial aneurysms using 3D vessel wall imaging.

BACKGROUND: Few studies directed at assessing the visualization of the walls of unruptured aneurysms have used higher-resolution 3D MRI vessel wall imaging. Prospective longitudinal studies are also needed to screen vessel wall changes in unruptured aneurysms. PURPOSE: To compare the aneurysm wall visualization on pre- and post-3D isotropic T1 -weighted Sampling Perfection with Application-optimized Contrasts by using different flip angle Evolutions (SPACE) images and to explore whether there is a change in wall enhancement at follow up. STUDY TYPE: Prospective. POPULATION: Twenty-nine patients with a total of 35 unruptured intracranial aneurysms. SEQUENCE: 3D T1 -weighted pre- and postcontrast SPACE (0.5 mm isotropic) at 3.0T. ASSESSMENT: The aneurysm wall visibility (0-5 scale) between pre- and postcontrast images as well as the wall enhancement (0-5 scale) between follow-up and baseline studies (6-30 months, average 12.7 months) were compared. Differences in wall visibility and enhancement were also investigated as a function of aneurysm diameter and location. STATISTICAL TEST: The Wilcoxon signed rank paired test was used to compare the wall visibility score between pre- and postcontrast SPACE images, as well as wall enhancement between follow-up and baseline. The Mann-Whitney and Kruskal-Wallis tests were used to investigate the enhancement difference between different diameters and locations. RESULTS: Postcontrast images had significantly higher wall visibility (P = 0.01). A wall enhancement score ≥2 was found in 71% of the aneurysms. Changes in levels of wall enhancement were found in 17% of the aneurysms at follow-up studies, but those changes were small. Wall visibility and enhancement scores of large aneurysms were significantly higher than small ones (P < 0.001). DATA CONCLUSION: 3D T1 -weighted higher resolution SPACE can be used to assess changes in enhancement at follow-up studies. Contrast SPACE image provides better aneurysm wall visibility and improves visualization of the aneurysm wall. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019;50:193-200.

Wall stress on ascending thoracic aortic aneurysms with bicuspid compared with tricuspid aortic valve.

OBJECTIVE: Guidelines for repair of bicuspid aortic valve-associated ascending thoracic aortic aneurysms have been changing, most recently to the same criteria as tricuspid aortic valve-ascending thoracic aortic aneurysms. Rupture/dissection occurs when wall stress exceeds wall strength. Recent studies suggest similar strength of bicuspid aortic valve versus tricuspid aortic valve-ascending thoracic aortic aneurysms; thus, comparative wall stress may better predict dissection in bicuspid aortic valve versus tricuspid aortic valve-ascending thoracic aortic aneurysms. Our aim was to determine whether bicuspid aortic valve-ascending thoracic aortic aneurysms had higher wall stresses than their tricuspid aortic valve counterparts. METHODS: Patients with bicuspid aortic valve- and tricuspid aortic valve-ascending thoracic aortic aneurysms (bicuspid aortic valve = 17, tricuspid aortic valve = 19) greater than 4.5 cm underwent electrocardiogram-gated computed tomography angiography. Patient-specific 3-dimensional geometry was reconstructed and loaded to systemic pressure after accounting for prestress geometry. Finite element analyses were performed using the LS-DYNA solver (LSTC Inc, Livermore, Calif) with user-defined fiber-embedded material model to determine ascending thoracic aortic aneurysm wall stress. RESULTS: Bicuspid aortic valve-ascending thoracic aortic aneurysms 99th-percentile longitudinal stresses were 280 kPa versus 242 kPa (P = .028) for tricuspid aortic valve-ascending thoracic aortic aneurysms in systole. These stresses did not correlate to diameter for bicuspid aortic valve-ascending thoracic aortic aneurysms (r = -0.004) but had better correlation to tricuspid aortic valve-ascending thoracic aortic aneurysms diameter (r = 0.677). Longitudinal stresses on sinotubular junction were significantly higher in bicuspid aortic valve-ascending thoracic aortic aneurysms than in tricuspid aortic valve-ascending thoracic aortic aneurysms (405 vs 329 kPa, P = .023). Bicuspid aortic valve-ascending thoracic aortic aneurysm 99th-percentile circumferential stresses were 548 kPa versus 462 kPa (P = .033) for tricuspid aortic valve-ascending thoracic aortic aneurysms, which also did not correlate to bicuspid aortic valve-ascending thoracic aortic aneurysm diameter (r = 0.007). CONCLUSIONS: Circumferential and longitudinal stresses were greater in bicuspid aortic valve- than tricuspid aortic valve-ascending thoracic aortic aneurysms and were more pronounced in the sinotubular junction. Peak wall stress did not correlate with bicuspid aortic valve-ascending thoracic aortic aneurysm diameter, suggesting diameter alone in this population may be a poor predictor of dissection risk. Our results highlight the need for patient-specific aneurysm wall stress analysis for accurate dissection risk prediction.

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms.

OBJECTIVES: To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging. MATERIALS AND METHODS: Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared. RESULTS: Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: -3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, -1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6). CONCLUSION: The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

Assessment of Reynolds stress components and turbulent pressure loss using 4D flow MRI with extended motion encoding.

PURPOSE: To measure the Reynolds stress tensor using 4D flow MRI, and to evaluate its contribution to computed pressure maps. METHODS: A method to assess both velocity and Reynolds stress using 4D flow MRI is presented and evaluated. The Reynolds stress is compared by cross-sectional integrals of the Reynolds stress invariants. Pressure maps are computed using the pressure Poisson equation-both including and neglecting the Reynolds stress. RESULT: Good agreement is seen for Reynolds stress between computational fluid dynamics, simulated MRI, and MRI experiment. The Reynolds stress can significantly influence the computed pressure loss for simulated (eg, -0.52% vs -15.34% error; P < 0.001) and experimental (eg, 306 ± 11 vs 203 ± 6 Pa; P < 0.001) data. A 54% greater pressure loss is seen at the highest experimental flow rate when accounting for Reynolds stress (P < 0.001). CONCLUSION: 4D flow MRI with extended motion-encoding enables quantification of both the velocity and the Reynolds stress tensor. The additional information provided by this method improves the assessment of pressure gradients across a stenosis in the presence of turbulence. Unlike conventional methods, which are only valid if the flow is laminar, the proposed method is valid for both laminar and disturbed flow, a common presentation in diseased vessels. Magn Reson Med 79:1962-1971, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Gated thoracic magnetic resonance angiography at 3T: noncontrast versus blood pool contrast.

Both noncontrast and contrast-enhanced approaches to gated thoracic magnetic resonance angiography (MRA) for aortic root evaluation have been reported at 3T. We compare qualitative and quantitative image quality measures for the two approaches, and assess the reproducibility of standard aortic measurements. Respiratory and cardiac gated MRA of the chest was performed at 3T in 45 patients: 23 after administration of iron-based blood pool contrast, and 22 without contrast. Image quality was assessed with a 5-point Likert scale, vessel lumen-to-muscle contrast ratios, and vessel wall sharpness. Two reviewers measured the ascending aorta diameter and valve annulus area. Interrater agreement was assessed using Bland-Altman plots and coefficient of variation (CV). Qualitative image quality was better with blood pool contrast in all principal vessels of the chest (mean Likert of 4.20 ± 0.79 vs. 2.60 ± 0.77, p < 0.001). Quantitative assessment was also improved with higher contrast ratios in all vessels (5.26 ± 3.3 vs. 1.90 ± 0.53, p < 0.001), and greater sharpness of the aortic annulus and ascending aorta (0.70 ± 0.16 vs. 0.56 ± 0.14 mm-1, p < 0.001, and 0.87 ± 0.16 vs. 0.62 ± 0.16 mm-1, p = 0.008, respectively). Reproducibility of measurement was marginally better for the ascending aorta diameter (CV of 2.80 vs. 3.23%), but substantially increased for the aortic valve annulus area with blood pool contrast (CV of 4.93 vs. 7.32%). The use of a blood pool contrast agent for gated thoracic MRA improves image quality compared to a noncontrast technique, and provides more reproducible measurements of the aortic valve annulus area.

Segmentation of lumen and outer wall of abdominal aortic aneurysms from 3D black-blood MRI with a registration based geodesic active contour model.

Segmentation of the geometric morphology of abdominal aortic aneurysm is important for interventional planning. However, the segmentation of both the lumen and the outer wall of aneurysm in magnetic resonance (MR) image remains challenging. This study proposes a registration based segmentation methodology for efficiently segmenting MR images of abdominal aortic aneurysms. The proposed methodology first registers the contrast enhanced MR angiography (CE-MRA) and black-blood MR images, and then uses the Hough transform and geometric active contours to extract the vessel lumen by delineating the inner vessel wall directly from the CE-MRA. The proposed registration based geometric active contour is applied to black-blood MR images to generate the outer wall contour. The inner and outer vessel wall are then fused presenting the complete vessel lumen and wall segmentation. The results obtained from 19 cases showed that the proposed registration based geometric active contour model was efficient and comparable to manual segmentation and provided a high segmentation accuracy with an average Dice value reaching 89.79%.

Suture Forces for Closure of Transapical Transcatheter Aortic Valve Replacement: A Mathematical Model.

BACKGROUND: Transcatheter aortic valve replacement (TAVR) has revolutionized the treatment of severe aortic stenosis in intermediate, high-risk, and inoperable patients. TAVR has multiple access routes, including transfemoral (TF), transapical (TA), direct aortic (DA), axillary, transcarotid, and transcaval. The most commonly applied algorithm is a TF-first approach, where only when patients are unsuitable for TF are alternatives such as TA considered. An infrequent - but dreaded - risk is left ventricular (LV) apical bleeding from tearing or rupture with the TA approach. With burgeoning transcatheter mitral technology that requires a TA approach, the study aim was to develop a mathematical model to determine suture forces for TA closure. METHODS: Preoperative cine-cardiac magnetic resonance imaging (MRI) was used to acquire three-dimensional (3D) LV geometry at end-systole and end-diastole. Endocardial and epicardial boundaries were manually contoured using MeVisLab, a surface reconstruction software. 3D surfaces of endocardium and epicardium were reconstructed, and surfaces at end-systole were used to create a 3D LV finite element (FE) mesh. TA access was mimicked by developing a 10-mm defect within the LV FE model. The LV apex was closed using a virtual suture technique in FE analysis with the application of two virtual sutures. After virtual closure, a FE analysis was performed of LV model diastolic filling and systolic contraction. RESULTS: Proof of concept was achieved to develop an LV transapical access site and perform FE analysis to achieve closure. The FE method of virtual suture technique successfully approximated the LV apical defect. The peak axial forces on virtual sutures at end-diastole and end-systole were 0.445N and 0.736N, respectively. CONCLUSIONS: A LV TA access model was mathematically developed that could be used to evaluate the suture tension of the TA closure process. Further development of this approach may be useful to risk-stratify patients in the future for LV apical tearing. Video 1: Cine cardiac magnetic resonance imaging of the left ventricle. Video 2: Slow motion animation of left ventricular baseline simulation. Video 3: Animation of the virtual suturing process.

Assessment of turbulent viscous stress using ICOSA 4D Flow MRI for prediction of hemodynamic blood damage.

Flow-induced blood damage plays an important role in determining the hemodynamic impact of abnormal blood flow, but quantifying of these effects, which are dominated by shear stresses in highly fluctuating turbulent flow, has not been feasible. This study evaluated the novel application of turbulence tensor measurements using simulated 4D Flow MRI data with six-directional velocity encoding for assessing hemodynamic stresses and corresponding blood damage index (BDI) in stenotic turbulent blood flow. The results showed that 4D Flow MRI underestimates the maximum principal shear stress of laminar viscous stress (PLVS), and overestimates the maximum principal shear stress of Reynolds stress (PRSS) with increasing voxel size. PLVS and PRSS were also overestimated by about 1.2 and 4.6 times at medium signal to noise ratio (SNR) = 20. In contrast, the square sum of the turbulent viscous shear stress (TVSS), which is used for blood damage index (BDI) estimation, was not severely affected by SNR and voxel size. The square sum of TVSS and the BDI at SNR >20 were underestimated by less than 1% and 10%, respectively. In conclusion, this study demonstrated the feasibility of 4D Flow MRI based quantification of TVSS and BDI which are closely linked to blood damage.

High resolution imaging of the intracranial vessel wall at 3 and 7 T using 3D fast spin echo MRI.

OBJECTIVES: High resolution MRI of the intracranial vessel wall provides important insights in the assessment of intracranial vascular disease. This study aims to refine high resolution 3D MRI techniques for intracranial vessel wall imaging at both 3 and 7 T using customized flip angle train design, and to explore their comparative abilities. MATERIALS AND METHODS: 11 patients with intracranial artery disease (four atherosclerotic plaques, six aneurysms and one reversible cerebral vasoconstriction syndrome) were imaged at 3 and 7 T with a 3D T 1-weighted fast-spin-echo sequence (SPACE) both pre and post Gd contrast injection. Wall to lumen contrast ratio (CRwall-lumen), contrast enhancement ratio (ER) and the sharpness of the vessel wall were quantified. Two experienced radiologists evaluated the image quality on a 0-5 scale. RESULTS: Both 3 and 7 T achieved good image quality with high resolution (nominal 0.5 mm isotropic) and whole brain coverage. The CRwall-lumen and the ER measurements were comparable (p > 0.05). The 7 T images were significantly sharper (sharpness: 2.69 ± 0.50 vs. 1.88 ± 0.53 mm(-1), p < 0.001) with higher image quality (reader 1 score: 3.5 ± 1.1 vs. 2.4 ± 1.1, p = 0.002) compared to 3 T. CONCLUSIONS: 3D T 1-weighted SPACE can be used for intracranial vessel wall evaluation at both 3 and 7 T. 7 T provides significantly better image quality and improves the confidence of diagnosis.

Isotropic 3D black blood MRI of abdominal aortic aneurysm wall and intraluminal thrombus.

INTRODUCTION: The aortic wall and intraluminal thrombus (ILT) have been increasingly studied as potential markers of progressive disease with abdominal aortic aneurysms (AAAs). Our goal was to develop a high resolution, 3D black blood MR technique for AAA wall and ILT imaging within a clinically acceptable scan time. METHODS: Twenty two patients with AAAs (maximal diameter 4.3±1.0cm), along with five healthy volunteers, were imaged at 3T with a 3D T1-weighted fast-spin-echo sequence using variable flip angle trains (SPACE) with a preparation pulse (DANTE) for suppressing blood signal. Volunteers and ten patients were also scanned with SPACE alone for comparison purposes. The signal to noise ratio (SNR) and the aortic wall/ILT to lumen contrast to noise ratio (CNR) were measured. Qualitative image scores (1-4 scale) assessing the inner lumen and outer wall boundaries of AAA were performed by two blinded reviewers. In patients with ILT, the ratio of ILT signal intensity (ILTSI) over psoas muscle SI (MuscleSI) was calculated, and the signal heterogeneity of ILT was quantified as standard deviation (SD) over the mean. RESULTS: All subjects were imaged successfully with an average scan time of 7.8±0.7minutes. The DANTE preparation pulse for blood suppression substantially reduced flow artifacts in SPACE with lower lumen SNR (8.8 vs. 21.4, p<0.001) and improved the wall/ILT to lumen CNR (9.9 vs. 6.3, p<0.001) in patients. Qualitative assessment showed improved visualization of lumen boundaries (73% higher scores on average, p=0.01) and comparable visualization of outer wall boundary (p>0.05). ILT was present in ten patients, with relatively high signal and a wide SD (average ILTSI/MuscleSI 1.42±0.48 (range 0.75-2.11)) and with SD/mean of 27.7%±6.6% (range 19.6%-39.4%). CONCLUSION: High resolution, 3D black blood MRI of AAAs can be achieved in a clinical accepted scan time with reduction of flow artifacts using the DANTE preparation pulse. Signal characteristics of ILT can be quantified and may be used for improved patient-specific risk stratification.

Ascending thoracic aortic aneurysm wall stress analysis using patient-specific finite element modeling of in vivo magnetic resonance imaging.

OBJECTIVES: Rupture/dissection of ascending thoracic aortic aneurysms (aTAAs) carries high mortality and occurs in many patients who did not meet size criteria for elective surgery. Elevated wall stress may better predict adverse events, but cannot be directly measured in vivo, rather determined from finite element (FE) simulations. Current computational models make assumptions that limit accuracy, most commonly using in vivo imaging geometry to represent zero-pressure state. Accurate patient-specific wall stress requires models with zero-pressure three-dimensional geometry, material properties, wall thickness and residual stress. We hypothesized that wall stress calculated from in vivo imaging geometry at systemic pressure underestimates that using zero-pressure geometry. We developed a novel method to derive zero-pressure geometry from in vivo imaging at systemic pressure. The purpose of this study was to develop the first patient-specific aTAA models using magnetic resonance imaging (MRI) to assess material properties and zero-pressure geometry. Wall stress results from FE models using systemic pressure were compared with those from models using zero-pressure correction. METHODS: Patients with aTAAs <5 cm underwent ECG-gated computed tomography angiography (CTA) and displacement encoding with stimulated echo (DENSE)-MRI. CTA lumen geometry was used to create surface contour meshes of aTAA geometry. DENSE-MRI measured cyclic aortic wall strain from which wall material property was derived. Zero- and systemic pressure geometries were created. Simulations were loaded to systemic pressure using the ABAQUS FE software. Wall stress analyses were compared between zero-pressure-corrected and systemic pressure geometry FE models. RESULTS: Peak first principal wall stress (primarily aligned in the circumferential direction) at systolic pressure for the zero-pressure correction models was 430.62 ± 69.69 kPa, whereas that without zero-pressure correction was 312.55 ± 39.65 kPa (P = 0.004). Peak second principal wall stress (primarily aligned in the longitudinal direction) at systolic pressure for the zero-pressure correction models was 200.77 ± 43.13 kPa, whereas that without zero-stress correction was 156.25 ± 25.55 kPa (P = 0.02). CONCLUSIONS: Previous FE aTAA models from in vivo CT and MRI have not accounted for zero-pressure geometry or patient-specific material property. We demonstrated that zero-pressure correction significantly impacts wall stress results. Future computational models that use wall stress to predict aTAA adverse events must take into account zero-pressure geometry and patient material property for accurate wall stress determination.