Abstract Computational fluid dynamic methods are currently being used clinically to simulate blood flow and pressure and predict the functional significance of atherosclerotic lesions in ...patient-specific models of the coronary arteries extracted from noninvasive coronary computed tomography angiography (cCTA) data. One such technology, FFR CT , or noninvasive fractional flow reserve derived from CT data, has demonstrated high diagnostic accuracy as compared to invasively measured fractional flow reserve (FFR) obtained with a pressure wire inserted in the coronary arteries during diagnostic cardiac catheterization. However, uncertainties in modeling as well as measurement results in differences between these predicted and measured hemodynamic indices. Uncertainty in modeling can manifest in two forms – anatomic uncertainty resulting in error of the reconstructed 3D model and physiologic uncertainty resulting in errors in boundary conditions or blood viscosity. We present a data-driven framework for modeling these uncertainties and study their impact on blood flow simulations. The incompressible Navier–Stokes equations are used to model blood flow and an adaptive stochastic collocation method is used to model uncertainty propagation in the Navier–Stokes equations. We perform uncertainty quantification in two geometries, an idealized stenosis model and a patient specific model. We show that uncertainty in minimum lumen diameter (MLD) has the largest impact on hemodynamic simulations, followed by boundary resistance, viscosity and lesion length. We show that near the diagnostic cutoff ( FFR CT = 0.8 ), the uncertainty due to the latter three variables are lower than measurement uncertainty, while the uncertainty due to MLD is only slightly higher than measurement uncertainty. We also show that uncertainties are not additive but only slightly higher than the highest single parameter uncertainty. The method presented here can be used to output interval estimates of hemodynamic indices and visualize patient-specific maps of sensitivities.
Top-down prefrontal cortex inputs to the hippocampus have been hypothesized to be important in memory consolidation, retrieval, and the pathophysiology of major psychiatric diseases; however, no such ...direct projections have been identified and functionally described. Here we report the discovery of a monosynaptic prefrontal cortex (predominantly anterior cingulate) to hippocampus (CA3 to CA1 region) projection in mice, and find that optogenetic manipulation of this projection (here termed AC-CA) is capable of eliciting contextual memory retrieval. To explore the network mechanisms of this process, we developed and applied tools to observe cellular-resolution neural activity in the hippocampus while stimulating AC-CA projections during memory retrieval in mice behaving in virtual-reality environments. Using this approach, we found that learning drives the emergence of a sparse class of neurons in CA2/CA3 that are highly correlated with the local network and that lead synchronous population activity events; these neurons are then preferentially recruited by the AC-CA projection during memory retrieval. These findings reveal a sparsely implemented memory retrieval mechanism in the hippocampus that operates via direct top-down prefrontal input, with implications for the patterning and storage of salient memory representations.
The goal of understanding living nervous systems has driven interest in high-speed and large field-of-view volumetric imaging at cellular resolution. Light sheet microscopy approaches have emerged ...for cellular-resolution functional brain imaging in small organisms such as larval zebrafish, but remain fundamentally limited in speed. Here, we have developed SPED light sheet microscopy, which combines large volumetric field-of-view via an extended depth of field with the optical sectioning of light sheet microscopy, thereby eliminating the need to physically scan detection objectives for volumetric imaging. SPED enables scanning of thousands of volumes-per-second, limited only by camera acquisition rate, through the harnessing of optical mechanisms that normally result in unwanted spherical aberrations. We demonstrate capabilities of SPED microscopy by performing fast sub-cellular resolution imaging of CLARITY mouse brains and cellular-resolution volumetric Ca2+ imaging of entire zebrafish nervous systems. Together, SPED light sheet methods enable high-speed cellular-resolution volumetric mapping of biological system structure and function.
Display omitted
•Light sheet microscopy speed is increased by extending the detection depth of field•A simple, scalable method is developed for extending the axial point spread function•Rapid, cellular-resolution nervous system mapping across the entire larval zebrafish•Fast automated identification of co-active neurons across the nervous system
By harnessing optical mechanisms that normally result in unwanted spherical aberrations, SPED light sheet microscopy allows high-speed mapping of biological structures such as the entire vertebrate nervous system and its activity at a cellular resolution.
In the cardiovascular system, blood flow rates, blood velocities and blood pressures can be modeled using the Navier–Stokes equations. Inputs to the system are typically uncertain, such as (a) the ...geometry of the arterial tree, (b) clinically measured blood pressure and viscosity, (c) boundary resistances, among others. Due to a large number of such parameters, efficient quantification of uncertainty in solution fields in this multi-parameter space is challenging. We use an adaptive stochastic collocation method to quantify the impact of uncertainty in geometry in patient-specific models. We develop a novel subdivision method to define the stochastic space of geometries. To accelerate convergence and make the problem tractable, we use a machine learning approach to approximate the simulation-based solution. Towards this, a reduced order model of the Navier–Stokes equations is developed using a segmental resistance analog boundary conditions (ratio of pressure to flow). Using an offline database of pre-computed solutions, we compute a map (rule) from the features to solution fields. We achieve significant speed-up (of a few orders of magnitude) by approximating the simulation-based solution using a machine learning predictor. A bootstrap aggregated decision tree was found to be the best predictor among many candidate regressors (correlation coefficient of training set was 0.94). We demonstrate stochastic space convergence using the adaptive stochastic collocation method, and also show robustness to the choice of geometry parameterization. The sensitivities to geometry obtained using machine learning had a correlation coefficient of 0.92 with the values obtained using finite element simulations. Segments with significant disease in the larger arteries had the highest sensitivities. Terminal segments are more sensitive to dilation and proximal healthy segments are more sensitive to erosion. Sensitivity to geometry is highest when geometric resistance is comparable to net downstream resistance.
Simulations of blood flow in both healthy and diseased vascular models can be used to compute a range of hemodynamic parameters including velocities, time varying wall shear stress, pressure drops, ...and energy losses. The confidence in the data output from cardiovascular simulations depends directly on our level of certainty in simulation input parameters. In this work, we develop a general set of tools to evaluate the sensitivity of output parameters to input uncertainties in cardiovascular simulations. Uncertainties can arise from boundary conditions, geometrical parameters, or clinical data. These uncertainties result in a range of possible outputs which are quantified using probability density functions (PDFs). The objective is to systemically model the input uncertainties and quantify the confidence in the output of hemodynamic simulations. Input uncertainties are quantified and mapped to the stochastic space using the stochastic collocation technique. We develop an adaptive collocation algorithm for Gauss-Lobatto-Chebyshev grid points that significantly reduces computational cost. This analysis is performed on two idealized problems--an abdominal aortic aneurysm and a carotid artery bifurcation, and one patient specific problem--a Fontan procedure for congenital heart defects. In each case, relevant hemodynamic features are extracted and their uncertainty is quantified. Uncertainty quantification of the hemodynamic simulations is done using (a) stochastic space representations, (b) PDFs, and (c) the confidence intervals for a specified level of confidence in each problem.
We present a computational framework for multiscale modeling and simulation of blood flow in coronary artery bypass graft (CABG) patients. Using this framework, only CT and non-invasive clinical ...measurements are required without the need to assume pressure and/or flow waveforms in the coronaries and we can capture global circulatory dynamics. We demonstrate this methodology in a case study of a patient with multiple CABGs. A patient-specific model of the blood vessels is constructed from CT image data to include the aorta, aortic branch vessels (brachiocephalic artery and carotids), the coronary arteries and multiple bypass grafts. The rest of the circulatory system is modeled using a lumped parameter network (LPN) 0 dimensional (0D) system comprised of resistances, capacitors (compliance), inductors (inertance), elastance and diodes (valves) that are tuned to match patient-specific clinical data. A finite element solver is used to compute blood flow and pressure in the 3D (3 dimensional) model, and this solver is implicitly coupled to the 0D LPN code at all inlets and outlets. By systematically parameterizing the graft geometry, we evaluate the influence of graft shape on the local hemodynamics, and global circulatory dynamics. Virtual manipulation of graft geometry is automated using Bezier splines and control points along the pathlines. Using this framework, we quantify wall shear stress, wall shear stress gradients and oscillatory shear index for different surgical geometries. We also compare pressures, flow rates and ventricular pressure–volume loops pre- and post-bypass graft surgery. We observe that PV loops do not change significantly after CABG but that both coronary perfusion and local hemodynamic parameters near the anastomosis region change substantially. Implications for future patient-specific optimization of CABG are discussed.
Vein maladaptation, leading to poor long-term patency, is a serious clinical problem in patients receiving coronary artery bypass grafts (CABGs) or undergoing related clinical procedures that subject ...veins to elevated blood flow and pressure. We propose a computational model of venous adaptation to altered pressure based on a constrained mixture theory of growth and remodeling (G&R). We identify constitutive parameters that optimally match biaxial data from a mouse vena cava, then numerically subject the vein to altered pressure conditions and quantify the extent of adaptation for a biologically reasonable set of bounds for G&R parameters. We identify conditions under which a vein graft can adapt optimally and explore physiological constraints that lead to maladaptation. Finally, we test the hypothesis that a gradual, rather than a step, change in pressure will reduce maladaptation. Optimization is used to accelerate parameter identification and numerically evaluate hypotheses of vein remodeling.
Objectives
Our goal was to evaluate the efficacy of a fully automated method for assessing the image quality (IQ) of coronary computed tomography angiography (CCTA).
Methods
The machine learning ...method was trained using 75 CCTA studies by mapping features (noise, contrast, misregistration scores, and un-interpretability index) to an IQ score based on manual ground truth data. The automated method was validated on a set of 50 CCTA studies and subsequently tested on a new set of 172 CCTA studies against visual IQ scores on a 5-point Likert scale.
Results
The area under the curve in the validation set was 0.96. In the 172 CCTA studies, our method yielded a Cohen’s kappa statistic for the agreement between automated and visual IQ assessment of 0.67 (
p
< 0.01). In the group where good to excellent (
n
= 163), fair (
n
= 6), and poor visual IQ scores (
n
= 3) were graded, 155, 5, and 2 of the patients received an automated IQ score > 50 %, respectively.
Conclusion
Fully automated assessment of the IQ of CCTA data sets by machine learning was reproducible and provided similar results compared with visual analysis within the limits of inter-operator variability.
Key points
•
The proposed method enables automated and reproducible image quality assessment.
•
Machine learning and visual assessments yielded comparable estimates of image quality.
•
Automated assessment potentially allows for more standardised image quality.
•
Image quality assessment enables standardization of clinical trial results across different datasets.
Kawasaki Disease (KD) is the leading cause of acquired pediatric heart disease. A subset of KD patients develops aneurysms in the coronary arteries, leading to increased risk of thrombosis and ...myocardial infarction. Currently, there are limited clinical data to guide the management of these patients, and the hemodynamic effects of these aneurysms are unknown. We applied patient-specific modeling to systematically quantify hemodynamics and wall shear stress in coronary arteries with aneurysms caused by KD. We modeled the hemodynamics in the aneurysms using anatomic data obtained by multi-detector computed tomography (CT) in a 10-year-old male subject who suffered KD at age 3 years. The altered hemodynamics were compared to that of a reconstructed normal coronary anatomy using our subject as the model. Computer simulations using a robust finite element framework were used to quantify time-varying shear stresses and particle trajectories in the coronary arteries. We accounted for the cardiac contractility and the microcirculation using physiologic downstream boundary conditions. The presence of aneurysms in the proximal coronary artery leads to flow recirculation, reduced wall shear stress within the aneurysm, and high wall shear stress gradients at the neck of the aneurysm. The wall shear stress in the KD subject (2.95–3.81 dynes/sq cm) was an order of magnitude lower than the normal control model (17.10–27.15 dynes/sq cm). Particle residence times were significantly higher, taking 5 cardiac cycles to fully clear from the aneurysmal regions in the KD subject compared to only 1.3 cardiac cycles from the corresponding regions of the normal model. In this novel quantitative study of hemodynamics in coronary aneurysms caused by KD, we documented markedly abnormal flow patterns that are associated with increased risk of thrombosis. This methodology has the potential to provide further insights into the effects of aneurysms in KD and to help risk stratify patients for appropriate medical and surgical interventions.
Fractional flow reserve (FFR) is commonly used to assess the functional significance of coronary artery disease but is theoretically limited in evaluating individual stenoses in serially diseased ...vessels. We sought to characterize the accuracy of assessing individual stenoses in serial disease using invasive FFR pullback and the noninvasive equivalent, fractional flow reserve by computed tomography (FFR
). We subsequently describe and test the accuracy of a novel noninvasive FFR
-derived percutaneous coronary intervention (PCI) planning tool (FFR
) in predicting the true significance of individual stenoses.
Patients with angiographic serial coronary artery disease scheduled for PCI were enrolled and underwent prospective coronary CT angiography with conventional FFR
-derived post hoc for each vessel and stenosis (FFR
). Before PCI, the invasive hyperemic pressure-wire pullback was performed to derive the apparent FFR contribution of each stenosis (FFR
). The true FFR attributable to individual lesions (FFR
) was then measured following PCI of one of the lesions. The predictive accuracy of FFR
, FFR
, and the novel technique (FFR
) was then assessed against FFR
. From the 24 patients undergoing the protocol, 19 vessels had post hoc FFR
and FFR
calculation. When assessing the distal effect of all lesions, FFR
correlated moderately well with invasive FFR ( R=0.71; P<0.001). For lesion-specific assessment, there was significant underestimation of FFR
using FFR
(mean discrepancy, 0.06±0.05; P<0.001, representing a 42% error) and conventional trans-lesional FFR
(0.05±0.06; P<0.001, 37% error). Using FFR
, stenosis underestimation was significantly reduced to a 7% error (0.01±0.05; P<0.001).
FFR pullback and conventional FFR
significantly underestimate true stenosis contribution in serial coronary artery disease. A novel noninvasive FFR
-based PCI planner tool more accurately predicts the true FFR contribution of each stenosis in serial coronary artery disease.