Surface Reconstruction of the Pediatric Larynx via Structure from Motion Photogrammetry: A Pilot Study.

Endoscopy is the gold standard for characterizing pediatric airway disorders, however, it is limited for quantitative analysis due to lack of three-dimensional (3D) vision and poor stereotactic depth perception. We utilize structure from motion (SfM) photogrammetry, to reconstruct 3D surfaces of pathologic and healthy pediatric larynges from monocular two-dimensional (2D) endoscopy. Models of pediatric subglottic stenosis were 3D printed and airway endoscopies were simulated. 3D surfaces were successfully reconstructed from endoscopic videos of all models using an SfM analysis toolkit. Average subglottic surface error between SfM reconstructed surfaces and 3D printed models was 0.65 mm as measured by Modified Hausdorff Distance. Average volumetric similarity between SfM surfaces and printed models was 0.82 as measured by Jaccard Index. SfM can be used to accurately reconstruct 3D surface renderings of the larynx from 2D endoscopy video. This technique has immense potential for use in quantitative analysis of airway geometry and virtual surgical planning.

Computational fluid dynamics modeling aiding surgical planning in a toddler with Parkes Weber syndrome.

Parkes Weber syndrome is a fast-flow and slow-flow vascular anomaly with limb overgrowth that can lead to congestive heart failure and limb ischemia. Current management strategies have focused on symptom management with focal embolization. A pediatric case with early onset heart failure is reported. We discuss the use of computational fluid dynamics (CFD) modeling to guide a surgical management strategy in a toddler with an MAP2K1 mutation.

Shear stress associated with cardiopulmonary bypass induces expression of inflammatory cytokines and necroptosis in monocytes.

Cardiopulmonary bypass (CPB) is required during most cardiac surgeries. CBP drives systemic inflammation and multiorgan dysfunction that is especially severe in neonatal patients. Limited understanding of molecular mechanisms underlying CPB-associated inflammation presents a significant barrier to improve clinical outcomes. To better understand these clinical issues, we performed mRNA sequencing on total circulating leukocytes from neonatal patients undergoing CPB. Our data identify myeloid cells, particularly monocytes, as the major cell type driving transcriptional responses to CPB. Furthermore, IL-8 and TNF-α were inflammatory cytokines robustly upregulated in leukocytes from both patients and piglets exposed to CPB. To delineate the molecular mechanism, we exposed THP-1 human monocytic cells to CPB-like conditions, including artificial surfaces, high shear stress, and cooling/rewarming. Shear stress was found to drive cytokine upregulation via calcium-dependent signaling pathways. We also observed that a subpopulation of THP-1 cells died via TNF-α-mediated necroptosis, which we hypothesize contributes to post-CPB inflammation. Our study identifies a shear stress-modulated molecular mechanism that drives systemic inflammation in pediatric CPB patients. These are also the first data to our knowledge to demonstrate that shear stress causes necroptosis. Finally, we observe that calcium and TNF-α signaling are potentially novel targets to ameliorate post-CPB inflammation.

Analysis of Factors Influencing Accuracy of Volume Flow Measurement in Dialysis Access Fistulas Based on Duplex Ultrasound Simulation.

OBJECTIVE: We developed a duplex ultrasound simulator and used it to assess accuracy of volume flow measurements in dialysis access fistula (DAF) models. METHODS: The simulator consists of a mannequin, computer, and mock transducer. Each case is built from a patient's B-mode images that are used to create a 3-dimensional surface model of the DAF. Computational fluid dynamics is used to determine blood flow velocities based on model vessel geometry. The simulator displays real-time B-mode and color-flow images, and Doppler spectral waveforms are generated according to user-defined settings. Accuracy was assessed by scanning each case and measuring volume flow in the inflow artery and outflow vein for comparison with true volume flow values. RESULTS: Four examiners made 96 volume flow measurements on four DAF models. Measured volume flow deviated from the true value by 35 ± 36%. Mean absolute deviation from true volume flow was lower for arteries than veins (22 ± 19%, N = 48 vs. 58 ± 33%, N = 48, p < 0.0001). This finding is attributed to eccentricity of outflow veins which resulted in underestimating true cross-sectional area. Regression analysis indicated that error in measuring cross-sectional area was a predictor of error in volume flow measurement (ß = 0.948, p < 0.001). Volume flow error was reduced from 35 ± 36% to 9 ± 8% (p < 0.000001) by calculating vessel area as an ellipse. CONCLUSIONS: Duplex volume flow measurements are based on a circular vessel shape. DAF inflow arteries are circular, but outflow veins can be elliptical. Simulation-based analysis showed that error in measuring volume flow is mainly due to assumption of a circular vessel.

Genetic correlates of wall shear stress in a patient-specific 3D-printed cerebral aneurysm model.

OBJECTIVES: To study the correlation between wall shear stress and endothelial cell expression in a patient-specific, three-dimensional (3D)-printed model of a cerebral aneurysm. MATERIALS AND METHODS: A 3D-printed model of a cerebral aneurysm was created from a patient's angiogram. After populating the model with human endothelial cells, it was exposed to media under flow for 24 hours. Endothelial cell morphology was characterized in five regions of the 3D-printed model using confocal microscopy. Endothelial cells were then harvested from distinct regions of the 3D-printed model for mRNA collection and gene analysis via quantitative polymerase chain reaction (qPCR.) Cell morphology and mRNA measurement were correlated with computational fluid dynamics simulations. RESULTS: The model was successfully populated with endothelial cells, which survived under flow for 24 hours. Endothelial morphology showed alignment with flow in the proximal and distal parent vessel and aneurysm neck, but disorganization in the aneurysm dome. Genetic analysis of endothelial mRNA expression in the aneurysm dome and distal parent vessel was compared with the proximal parent vessels. ADAMTS-1 and NOS3 were downregulated in the aneurysm dome, while GJA4 was upregulated in the distal parent vessel. Disorganized morphology and decreased ADAMTS-1 and NOS3 expression correlated with areas of substantially lower wall shear stress and wall shear stress gradient in computational fluid dynamics simulations. CONCLUSIONS: Creating 3D-printed models of patient-specific cerebral aneurysms populated with human endothelial cells is feasible. Analysis of these cells after exposure to flow demonstrates differences in both cell morphology and genetic expression, which correlate with areas of differential hemodynamic stress.

Small Left Ventricular Size Is an Independent Risk Factor for Ventricular Assist Device Thrombosis.

The prevalence of ventricular assist device (VAD) therapy has continued to increase due to a stagnant donor supply and growing advanced heart failure (HF) population. We hypothesize that left ventricular (LV) size strongly influences biocompatibility and risk of thrombosis. Unsteady computational fluid dynamics (CFD) was used in conjunction with patient-derived computational modeling and virtual surgery with a standard, apically implanted inflow cannula. A dual-focus approach of evaluating thrombogenicity was employed: platelet-based metrics to characterize the platelet environment and flow-based metrics to investigate hemodynamics. Left ventricular end-diastolic dimensions (LVEDds) ranging from 4.5 to 6.5 cm were studied and ranked according to relative thrombogenic potential. Over 150,000 platelets were individually tracked in each LV model over 15 cardiac cycles. As LV size decreased, platelets experienced markedly increased shear stress histories (SHs), whereas platelet residence time (RT) in the LV increased with size. The complex interplay between increased SH and longer RT has profound implications on thrombogenicity, with a significantly higher proportion of platelets in small LVs having long RT times and being subjected to high SH, contributing to thrombus formation. Our data suggest that small LV size, rather than decreased VAD speed, is the primary pathologic mechanism responsible for the increased incidence of thrombosis observed in VAD patients with small LVs.

Left Ventricular Assist Device Inflow Cannula Angle and Thrombosis Risk.

BACKGROUND: As heart failure prevalence continues to increase in the setting of a static donor supply, left ventricular assist device (LVAD) therapy for end-stage heart failure continues to grow. Anecdotal evidence suggests that malalignment of the LVAD inflow cannula may increase thrombosis risk, but this effect has not been explored mechanistically or quantified statistically. Our objective is to elucidate the impact of surgical angulation of the inflow cannula on thrombogenicity. METHODS AND RESULTS: Unsteady computational fluid dynamics is used in conjunction with computational modeling and virtual surgery to model flow through the left ventricle for 5 different inflow cannula angulations. We use a holistic approach to evaluate thrombogenicity: platelet-based (Lagrangian) metrics to evaluate the platelet mechanical environment, combined with flow-based (Eulerian) metrics to investigate intraventricular hemodynamics. The thrombogenic potential of each LVAD inflow cannula angulation is quantitatively evaluated based on platelet shear stress history and residence time. Intraventricular hemodynamics are strongly influenced by LVAD inflow cannula angulation. Platelet behavior indicates elevated thrombogenic potential for certain inflow cannula angles, potentially leading to platelet activation. Our analysis demonstrates that the optimal range of inflow angulation is within 0±7° of the left ventricular apical axis. CONCLUSIONS: Angulation of the inflow cannula >7° from the apical axis (axis connecting mitral valve and ventricular apex) leads to markedly unfavorable hemodynamics as determined by computational fluid dynamics. Computational hemodynamic simulations incorporating Lagrangian and Eulerian metrics are a powerful tool for studying optimization of LVAD implantation strategies, with the long-term potential of improving outcomes.

Computational fluid dynamics of cerebral aneurysm coiling using high-resolution and high-energy synchrotron X-ray microtomography: comparison with the homogeneous porous medium approach.

BACKGROUND: Computational modeling of intracranial aneurysms provides insights into the influence of hemodynamics on aneurysm growth, rupture, and treatment outcome. Standard modeling of coiled aneurysms simplifies the complex geometry of the coil mass into a homogeneous porous medium that fills the aneurysmal sac. We compare hemodynamics of coiled aneurysms modeled from high-resolution imaging with those from the same aneurysms modeled following the standard technique, in an effort to characterize sources of error from the simplified model. MATERIALS: Physical models of two unruptured aneurysms were created using three-dimensional printing. The models were treated with coil embolization using the same coils as those used in actual patient treatment and then scanned by synchrotron X-ray microtomography to obtain high-resolution imaging of the coil mass. Computational modeling of each aneurysm was performed using patient-specific boundary conditions. The coils were modeled using the simplified porous medium or by incorporating the X-ray imaged coil surface, and the differences in hemodynamic variables were assessed. RESULTS: X-ray microtomographic imaging of coils and incorporation into computational models were successful for both aneurysms. Porous medium calculations of coiled aneurysm hemodynamics overestimated intra-aneurysmal flow, underestimated oscillatory shear index and viscous dissipation, and over- or underpredicted wall shear stress (WSS) and WSS gradient compared with X-ray-based coiled computational fluid dynamics models. CONCLUSIONS: Computational modeling of coiled intracranial aneurysms using the porous medium approach may inaccurately estimate key hemodynamic variables compared with models incorporating high-resolution synchrotron X-ray microtomographic imaging of complex aneurysm coil geometry.

LVAD Outflow Graft Angle and Thrombosis Risk.

This study quantifies thrombogenic potential (TP) of a wide range of left ventricular assist device (LVAD) outflow graft anastomosis angles through state-of-the-art techniques: 3D imaged-based patient-specific models created via virtual surgery and unsteady computational fluid dynamics with Lagrangian particle tracking. This study aims at clarifying the influence of a single parameter (outflow graft angle) on the thrombogenesis associated with flow patterns in the aortic root after LVAD implantation. This is an important and poorly-understood aspect of LVAD therapy, because several studies have shown strong inter and intrapatient thrombogenic variability and current LVAD implantation strategies do not incorporate outflow graft angle optimization. Accurate platelet-level investigation, enabled by statistical treatment of outliers in Lagrangian particle tracking, demonstrates a strong influence of outflow graft anastomoses angle on thrombogenicity (platelet residence times and activation state characterized by shear stress accumulation) with significantly reduced TP for acutely-angled anastomosed outflow grafts. The methodology presented in this study provides a device-neutral platform for conducting comprehensive thrombogenicity evaluation of LVAD surgical configurations, empowering optimal patient-focused surgical strategies for long-term treatment and care for advanced heart failure patients.

Development of a Duplex Ultrasound Simulator and Preliminary Validation of Velocity Measurements in Carotid Artery Models.

OBJECTIVE: Duplex ultrasound scanning with B-mode imaging and both color Doppler and Doppler spectral waveforms is relied upon for diagnosis of vascular pathology and selection of patients for further evaluation and treatment. In most duplex ultrasound applications, classification of disease severity is based primarily on alterations in blood flow velocities, particularly the peak systolic velocity (PSV) obtained from Doppler spectral waveforms. We developed a duplex ultrasound simulator for training and assessment of scanning skills. METHODS: Duplex ultrasound cases were prepared from 2-dimensional (2D) images of normal and stenotic carotid arteries by reconstructing the common carotid, internal carotid, and external carotid arteries in 3 dimensions and computationally simulating blood flow velocity fields within the lumen. The simulator displays a 2D B-mode image corresponding to transducer position on a mannequin, overlaid by color coding of velocity data. A spectral waveform is generated according to examiner-defined settings (depth and size of the Doppler sample volume, beam steering, Doppler beam angle, and pulse repetition frequency or scale). The accuracy of the simulator was assessed by comparing the PSV measured from the spectral waveforms with the true PSV which was derived from the computational flow model based on the size and location of the sample volume within the artery. RESULTS: Three expert examiners made a total of 36 carotid artery PSV measurements based on the simulated cases. The PSV measured by the examiners deviated from true PSV by 8% ± 5% (N = 36). The deviation in PSV did not differ significantly between artery segments, normal and stenotic arteries, or examiners. CONCLUSION: To our knowledge, this is the first simulation of duplex ultrasound that can create and display real-time color Doppler images and Doppler spectral waveforms. The results demonstrate that an examiner can measure PSV from the spectral waveforms using the settings on the simulator with a mean absolute error in the velocity measurement of less than 10%. With the addition of cases with a range of pathologies, this duplex ultrasound simulator will be a useful tool for training health-care providers in vascular ultrasound applications and for assessing their skills in an objective and quantitative manner.

Hemodynamic conditions in a failing peripheral artery bypass graft.

OBJECTIVE: The mechanisms of restenosis in autogenous vein bypass grafts placed for peripheral artery disease are not completely understood. We investigated the role of hemodynamic stress in a case study of a revised bypass graft that failed due to restenosis. METHODS: The morphology of the lumen was reconstructed from a custom three-dimensional ultrasound system. Scans were taken at 1, 6, and 16 months after a patch angioplasty procedure. Computational hemodynamic simulations of the patient-specific model provided the blood flow features and the hemodynamic stresses on the vessel wall at the three times studied. RESULTS: The vessel was initially free of any detectable lesions, but a 60% diameter-reducing stenosis developed during the 16-month study interval. As determined from the simulations, chaotic and recirculating flow occurred downstream of the stenosis due to the sudden widening of the lumen at the patch location. Curvature and a sudden increase in the lumen cross-sectional area induced these flow features that are hypothesized to be conducive to intimal hyperplasia. Favorable agreement was found between simulation results and in vivo Doppler ultrasound velocity measurements. CONCLUSIONS: Transitional and chaotic flow occurs at the site of the revision, inducing a complex pattern of wall shear as computed with the hemodynamic simulations. This supports the hypothesis that the hemodynamic stresses in the revised segment, produced by the coupling of vessel geometry and chaotic flow, led to the intimal hyperplasia and restenosis of the graft.