SummaryPreferential locations of atherosclerotic plaque are strongly associated with the areas of low wall shear stress and disturbed haemodynamic characteristics such as flow detachment, flow recirculation and oscillatory flow. The areas of low wall shear stress are also associated with the reduced production of adenosine triphosphate in the endothelial layer, as well as the resulting reduced production of inositol trisphosphate (IP3). The subsequent variation in Ca2+ signalling and nitric oxide synthesis could lead to the impairment of the atheroprotective function played by nitric oxide. In previous studies, it has been suggested that the reduced IP3 and Ca2+ signalling can explain the correlation of atherosclerosis with induced low WSS and disturbed flow characteristics. The massively parallel implementation described in this article provides insight into the dynamics of coupled smooth muscle cells and endothelial cells mapped onto the surface of an idealised arterial bifurcation. We show that variations in coupling parameters, which model normal and pathological conditions, provide vastly different smooth muscle cell Ca2+ dynamics and wave propagation profiles. The extensibility of the coupled cells model and scalability of the implementation provide a solid framework for in silico investigations of the interaction between complex cellular chemistry and the macro‐scale processes determined by fluid dynamics. © 2016 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.

Generic solutions for misalignment correction in image stacks are based on slice-byslice registration to a reference volume. Slice misalignment in cardiac Magnetic Resonance cine series typically includes translational and rotational components. Recovery of rotation in medial down to apical slices with the basic slice-to-volume registration may provide unreliable results because of rotational symmetry of the short-axis view of the left ventricle. This paper describes a novel method for correction of stack misalignment in cardiac cine series with recovery of the rotational component. The core of the presented method is a custom spatial transform which parametrises a set of independent centred rigid transforms, one per slice in a stack. The stack alignment transform was designed to circumvent the weaknesses associated with slice-to-volume registration and improve the reliability of misalignment correction over the slice-by-slice registration approach. The image similarity metric for every iteration during the stack alignment transform optimisation is calculated on the whole image stack with all slice correction parameters contributing to the result. The presented method was evaluated on 50 clinical datasets and the results demonstrate an improvement of misalignment correction with the proposed transform. The study also provides evidence of the importance of recovery of slice rotations in this application; the presented method is general enough to be applied to similar problems in other image modalities.

No abstract

Myocardial perfusion imaging, coupled with quantitative perfusion analysis, provides an important diagnostic tool for the identification of ischaemic heart disease caused by coronary stenoses. The accurate mapping between coronary anatomy and under-perfused areas of the myocardium is important for diagnosis and treatment. However, in the absence of the actual coronary anatomy during the reporting of perfusion images, areas of ischaemia are allocated to a coronary territory based on a population-derived 17-segment (American Heart Association) AHA model of coronary blood supply. This work presents a solution for the fusion of 2D Magnetic Resonance (MR) myocardial perfusion images and 3D MR angiography data with the aim to improve the detection of ischaemic heart disease. The key contribution of this work is a novel method for the mediated spatiotemporal registration of perfusion and angiography data and a novel method for the calculation of patient-specific coronary supply territories. The registration method uses 4D cardiac MR cine series spanning the complete cardiac cycle in order to overcome the under-constrained nature of non-rigid slice-to-volume perfusion-to-angiography registration. This is achieved by separating out the deformable registration problem and solving it through phase-to-phase registration of the cine series. The use of patient-specific blood supply territories in quantitative perfusion analysis (instead of the population-based model of coronary blood supply) has the potential of increasing the accuracy of perfusion analysis. Quantitative perfusion analysis diagnostic accuracy evaluation with patient-specific territories against the AHA model demonstrates the value of the mediated spatiotemporal registration in the context of ischaemic heart disease diagnosis.

This work investigates the effect of arterial bifurcation angulation on atherosclerosis development through in-silico simulations of coupled cell dynamics. The computational model presented here combines cellular pathways, fluid dynamics, and physiologically-realistic vessel geometries as observed in the human vasculature. The coupled cells model includes endothelial cells (ECs) and smooth muscle cells (SMCs) with ion dynamics, hetero and homotypic coupling, as well as electro-diffusive coupling. Three arterial bifurcation surface models were used in the coupled cells simulations. All three simulations showed propagating waves of Ca2+ in both the SMC and EC layers, following the introduction of a luminal agonist, in this case ATP. Immediately following the introduction of ATP concentration Ca2+ waves propagate from the area of high ATP toward the areas of low ATP concentration, forming complex patterns where waves interact with eachother, collide and fade. These dynamic phenomena are repeated with a series of waves of slower velocity. The underlying motivation of this research was to examine the macro-scale phenomena, given that the characteristic length scales of atherosclerotic plaques are much larger than a single cell. The micro-scale dynamics were modeled on macro-scale arterial bifurcation surfaces containing over one million cells. The results of the simulations presented here suggest that susceptibility to atherosclerosis development depends on the bifurcation angulation. In conjunction with findings reported in the literature, the simulation results demonstrate that arterial bifurcations containing wider angles have a more prominent influence on the coupled cells pathways associated with the development of atherosclerosis, by means of disturbed flow and lower SMC Ca2+ concentrations. The discussion of the results considers the findings of this research within the context of the potential link between information transport through frequency encoding of Ca2+ wave dynamics and development of atheroprone conditions.