Quantifying the effect of uncertainty in input parameters in a simplified bidomain model of partial thickness ischaemia

Reduced blood flow in the coronary arteries can lead to damaged heart tissue (myocardial ischaemia). Although one method for detecting myocardial ischaemia involves changes in the ST segment of the electrocardiogram, the relationship between these changes and subendocardial ischaemia is not fully understood. In this study, we modelled ST-segment epicardial potentials in a slab model of cardiac ventricular tissue, with a central ischaemic region, using the bidomain model, which considers conduction longitudinal, transverse and normal to the cardiac fibres. We systematically quantified the effect of uncertainty on the input parameters, fibre rotation angle, ischaemic depth, blood conductivity and six bidomain conductivities, on outputs that characterise the epicardial potential distribution. We found that three typical types of epicardial potential distributions (one minimum over the central ischaemic region, a tripole of minima, and two minima flanking a central maximum) could all occur for a wide range of ischaemic depths. In addition, the positions of the minima were affected by both the fibre rotation angle and the ischaemic depth, but not by changes in the conductivity values. We also showed that the magnitude of ST depression is affected only by changes in the longitudinal and normal conductivities, but not by the transverse conductivities

Comparison of diurnal variations, gestational age and gender related differences in fetal heart rate (FHR) parameters between appropriate-for-gestational-age (AGA) and small-for-gestational-age (SGA) fetuses in the home environment

Objective To assess the influence of gender, time of the day and gestational age on fetal heart rate (FHR) parameters between appropriate-for-gestational-age (AGA) and small-for-gestational age (SGA) fetuses using a portable fetal ECG monitor employed in the home setting. Methods We analysed and compared the antenatal FHR data collected in the home setting on 61 healthy pregnant women with singleton pregnancies from 24 weeks gestation. Of the 61 women, 31 had SGA fetuses (estimated fetal weight below the tenth gestational centile) and 30 were pregnant with AGA fetuses. FHR recordings were collected for up to 20 h. Two 90 min intervals were deliberately chosen retrospectively with respect to signal recording quality, one during day-time and one at night-time for comparison. Results Overall, success rate of the fetal abdominal ECG in the AGA fetuses was 75.7% compared to 48.6% in the SGA group. Based on randomly selected episodes of heart rate traces where recording quality exceeded 80% we were able to show a marginal difference between day and night-time recordings in AGA vs. SGA fetuses beyond 32 weeks of gestation. A selection bias in terms of covering different representation periods of fetal behavioural states cannot be excluded. In contrast to previous studies, we neither controlled maternal diet and activity nor measured maternal blood hormone and heart rate as all mothers were monitored in the home environment. Conclusion Based on clinically unremarkable, but statistically significant differences in the FHR parameters between the AGA and SGA group we suggest that further studies with large sample size are required to assess the clinical value of antenatal fetal ECG monitoring

3D finite element electrical model of larval zebrafish ECG signals

Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace’s equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions

Influence of ischemic core muscle fibers on surface depolarization potentials in superfused cardiac tissue preparations: a simulation study

Thin-walled cardiac tissue samples superfused with oxygenated solutions are widely used in experimental studies. However, due to decreased oxygen supply and insufficient wash out of waste products in the inner layers of such preparations, electrophysiological functions could be compromised. Although the cascade of events triggered by cutting off perfusion is well known, it remains unclear as to which degree electrophysiological function in viable surface layers is affected by pathological processes occurring in adjacent tissue. Using a 3D numerical bidomain model, we aim to quantify the impact of superfusion-induced heterogeneities occurring in the depth of the tissue on impulse propagation in superficial layers. Simulations demonstrated that both the pattern of activation as well as the distribution of extracellular potentials close to the surface remain essentially unchanged. This was true also for the electrophysiological properties of cells in the surface layer, where most relevant depolarization parameters varied by less than 5.5 %. The main observed effect on the surface was related to action potential duration that shortened noticeably by 53 % as hypoxia deteriorated. Despite the known limitations of such experimental methods, we conclude that superfusion is adequate for studying impulse propagation and depolarization whereas repolarization studies should consider the influence of pathological processes taking place at the core of tissue sample

Ketamine-Induced Oscillations in the Motor Circuit of the Rat Basal Ganglia

Oscillatory activity can be widely recorded in the cortex and basal ganglia. This activity may play a role not only in the physiology of movement, perception and cognition, but also in the pathophysiology of psychiatric and neurological diseases like schizophrenia or Parkinson's disease. Ketamine administration has been shown to cause an increase in gamma activity in cortical and subcortical structures, and an increase in 150 Hz oscillations in the nucleus accumbens in healthy rats, together with hyperlocomotion

Methods for Initialization of Activation Based Inverse Electrocardiography Using Graphs Derived from Heart Surface Geometry

The activation-based inverse problem of electrocardiography is non-linear in the desired activation times. Current solutions rely on iterative algorithms. There is considerable interest in improved initialization approaches due to the importance of the initialization used. Recent efforts include the critical point algorithm of Huiskamp and Greensite and the fastest route algorithm of van Dam, Oostendorp, and van Oosterom. In this work we analyze the relationship between these two methods. We also suggest an alternative to the shortest path approach to represent the set of likely activation patterns that may have computational advantages. We also explore modifications to these two methods exploiting their relationship and using the new activation patterns. We use epicardially stimulated data and geometries, recorded at the Cardiovascular Research and Training Institute in Utah, and geometries and forward matrix supplied with the ECGSim software, to compare results. 1