Search

In atrial fibrillation (AF), the ECG P-wave, which represents atrial depolarization, is replaced with chaotic and irregular fibrillation waves (f waves). The f-wave frequency, Ff, shows significant variations over time. Cardiorespiratory interactions regulated by the autonomic nervous system have been suggested to play a role in such variations. We conducted a simulation study to test whether the spatiotemporal release pattern of the parasympathetic neurotransmitter acetylcholine (ACh) modulates the frequency of atrial reentrant circuits. Understanding parasympathetic involvement in AF may guide more effective treatment approaches and could help to design autonomic markers alternative to heart rate variability (HRV), which is not available in AF patients. 2D tissue and 3D whole-atria models of human atrial electrophysiology in persistent AF were built. Different ACh release percentages (8% and 30%) and spatial ACh release patterns, including spatially random release and release from ganglionated plexi (GPs) and associated nerves, were considered. The temporal pattern of ACh release, ACh(t), was simulated following a sinusoidal waveform of frequency 0.125 Hz to represent the respiratory frequency. Different mean concentrations (ACh̄) and peak-to-peak ranges of ACh (ΔACh) were tested. We found that temporal variations in Ff, Ff(t), followed the simulated temporal ACh(t) pattern in all cases. The temporal mean of Ff(t), F̄f, depended on the fibrillatory pattern (number and location of rotors), the percentage of ACh release nodes and ACh̄. The magnitude of Ff(t) modulation, ΔFf, depended on the percentage of ACh release nodes and ΔACh. The spatial pattern of ACh release did not have an impact on F̄f and only a mild impact on ΔFf. The f-wave frequency, being indicative of vagal activity, has the potential to drive autonomic-based therapeutic actions and could replace HRV markers not quantifiable from AF patients.

Summary: The procedures used routinely for collagen and lipofuscin evaluation are, in many cases, qualitative, observer dependent, and disregard spatial distribution. Here, we present a protocol for automatic quantification and spatial characterization of collagen and lipofuscin from label-free microscopy images of human ventricular tissues. We describe the steps for sample collection, tissue processing, image acquisition, and quantification of collagen and lipofuscin. This protocol avoids discrepancies between observers and can be adapted to other tissues and species.For complete details on the use and execution of this protocol, please refer to García-Mendívil et al. (2022).1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Objective: Chronic kidney disease patients have a decreased ability to maintain normal electrolyte concentrations in their blood, which increases the risk for ventricular arrhythmias and sudden cardiac death. Non-invasive monitoring of serum potassium and calcium concentration, [K+] and [Ca2+], can help to prevent arrhythmias in these patients. Electrocardiogram (ECG) markers that significantly correlate with [K+] and [Ca2+] have been proposed, but these relations are highly variable between patients. We hypothesized that inter-individual differences in cell type distribution across the ventricular wall can help to explain this variability.Methods: A population of human heart-torso models were built with different proportions of endocardial, midmyocardial and epicardial cells. Propagation of ventricular electrical activity was described by a reaction-diffusion model, with modified Ten Tusscher-Panfilov dynamics. [K+] and [Ca2+] were varied individually and in combination. Twelve-lead ECGs were simulated and the width, amplitude and morphological variability of T waves and QRS complexes were quantified. Results were compared to measurements from 29 end-stage renal disease (ESRD) patients undergoing hemodialysis (HD).Results: Both simulations and patients data showed that most of the analyzed T wave and QRS complex markers correlated strongly with [K+] (absolute median Pearson correlation coefficients, r, ranging from 0.68 to 0.98) and [Ca2+] (ranging from 0.70 to 0.98). The same sign and similar magnitude of median r was observed in the simulations and the patients. Different cell type distributions in the ventricular wall led to variability in ECG markers that was accentuated at high [K+] and low [Ca2+], in agreement with the larger variability between patients measured at the onset of HD. The simulated ECG variability explained part of the measured inter-patient variability.Conclusion: Changes in ECG markers were similarly related to [K+] and [Ca2+] variations in our models and in the ESRD patients. The high inter-patient ECG variability may be explained by variations in cell type distribution across the ventricular wall, with high sensitivity to variations in the proportion of epicardial cells.Significance: Differences in ventricular wall composition help to explain inter-patient variability in ECG response to [K+] and [Ca2+]. This finding can be used to improve serum electrolyte monitoring in ESRD patients.

We present the meshfree mixed collocation method (MCM) for cardiac electrophysiology simulation. Capitalizing on the meshfree property of MCM, we introduce an immersed grid approach for automated generation of meshfree node grids from medical image data. This approach allows us to avoid the time-consuming mesh generation and processing that mesh-based methods like the finite element method (FEM) require. We employ the MCM to solve the cardiac monodomain model considering electrical propagation in 2D tissue sheets, 3D tissue slabs, and a realistic biventricular anatomy. We demonstrate that the solutions obtained by the MCM are in good agreement with the FEM, particularly when immersed grid is used. These findings confirm the suitability of the MCM for cardiac electrophysiology simulation and make the MCM a promising alternative to the FEM for cardiac electrical investigations.

Connexin 43 (CX43) is one of the major components of gap junctions, the structures responsible for the intercellular communication and transmission of the electrical impulse in the left ventricle. There is limited information on the histological changes of CX43 with age and their effect on electrophysiology, especially in humans. Here, we analyzed left ventricular biopsies from living donors starting at midlife to characterize age-related CX43 remodeling. We assessed its quantity, degree of lateralization, and spatial heterogeneity together with fibrotic deposition. We observed no significant age-related remodeling of CX43. Only spatial heterogeneity increased slightly with age, and this increase was better explained by biological age than by chronological age. Importantly, we found that CX43 features varied considerably among individuals in our population with no relevant relationship to age or fibrosis content, in contrast to animal species. We used our experimental results to feed computational models of human ventricular electrophysiology and to assess the effects of interindividual differences in specific features of CX43 and fibrosis on conduction velocity, action potential duration, and arrhythmogenicity. We found that larger amounts of fibrosis were associated with the highest arrhythmic risk, with this risk being increased when fibrosis deposition was combined with a reduction in CX43 amount and/or with an increase in CX43 spatial heterogeneity. These mechanisms underlying high arrhythmic risk in some individuals were not associated with age in our study population. In conclusion, our data rule out CX43 remodeling as an age-related arrhythmic substrate in the population beyond midlife, but highlight its potential as a proarrhythmic factor at the individual level, especially when combined with increased fibrosis.

The adaptation lag of the QT interval after heart rate (HR) has been proposed as an arrhythmic risk marker. Most studies have quantified the QT adaptation lag in response to abrupt, step-like changes in HR induced by atrial pacing, in response to tilt test or during ambulatory recordings. Recent studies have introduced novel methods to quantify the QT adaptation lag to gradual, ramp-like HR changes in stress tests by evaluating the differences between the measured QT series and an estimated, memoryless QT series obtained from the instantaneous HR. These studies have observed the QT adaptation lag to progressively reduce when approaching the stress peak, with the underlying mechanisms being still unclear. This study analyzes the contribution of β-adrenergic stimulation to QT interval rate adaptation in response to gradual, ramp-like HR changes. We first quantify the QT adaptation lag in Coronary Artery Disease (CAD) patients undergoing stress test. To uncover the involved mechanisms, we use biophysically detailed computational models coupling descriptions of human ventricular electrophysiology and β-adrenergic signaling, from which we simulate ventricular action potentials and ECG signals. We characterize the adaptation of the simulated QT interval in response to the HR time series measured from each of the analyzed CAD patients. We show that, when the simulated ventricular tissue is subjected to a time-varying β-adrenergic stimulation pattern, with higher stimulation levels close to the stress peak, the simulated QT interval presents adaptation lags during exercise that are more similar to those measured from the patients than when subjected to constant β-adrenergic stimulation. During stress test recovery, constant and time-varying β-adrenergic stimulation patterns render similar adaptation lags, which are generally shorter than during exercise, in agreement with results from the patients. In conclusion, our findings support the role of time-varying β-adrenergic stimulation in contributing to QT interval adaptation to gradually increasing HR changes as those seen during the exercise phase of a stress test.

Abstract The two most common modes of death among chronic heart failure (CHF) patients are sudden cardiac death (SCD) and pump failure death (PFD). Periodic repolarization dynamics (PRD) quantifies low-frequency oscillations in the T wave vector of the electrocardiogram (ECG) and has been postulated to reflect sympathetic modulation of ventricular repolarization. This study aims to evaluate the prognostic value of PRD to predict SCD and PFD in a population of CHF patients. 20-min high-resolution (1000 Hz) ECG recordings from 569 CHF patients were analyzed. Patients were divided into two groups, $$\hbox {PRD}^+$$ PRD + and $$\hbox {PRD}^-$$ PRD - , corresponding to PRD values above and below the optimum cutoff point of PRD in the study population. Univariate Cox regression analysis showed that SCD risk in the $$\hbox {PRD}^+$$ PRD + group was double the risk in the $$\hbox {PRD}^-$$ PRD - group [hazard ratio (95% CI) 2.001 (1.127–3.554), $$\hbox {p}

Cardiovascular diseases cause a high number of deaths nowadays. To improve these statistics, new strategies to better understand the electrical and mechanical abnormalities underlying them are urgently required. This study focuses on the development of a sensor to measure tissue stretch in excised tissues, enabling improved knowledge of biomechanical properties and allowing greater control in real time. A system made of biocompatible materials is described, which is based on two cantilevered platforms that integrate an optical fiber inside them to quantify the amount of stretch the tissues are exposed to with a precision of μm. The operating principle of the sensor is based on the variation of the optical path with the movement of the platforms onto which the samples are fixed. The conducted tests highlight that this system, based on a simple topology and technology, is capable of achieving the desired purpose (a resolution of ∼1 μm), enabling the tissue to be bathed in any medium within the system.

Abstract We investigated the ability of time-warping-based ECG-derived markers of T-wave morphology changes in time ( $$d_{w}$$ d w ) and amplitude ( $$d_a$$ d a ), as well as their non-linear components ( $${d_w^{{\mathrm{NL}}}}$$ d w NL and $${d_a^{\mathrm{NL}}}$$ d a NL ), and the heart rate corrected counterpart ( $$d_{w,c}$$ d w , c ), to monitor potassium concentration ( $$[K^{+}]$$ [ K + ] ) changes ( $$\Delta [K^+]$$ Δ [ K + ] ) in end-stage renal disease (ESRD) patients undergoing hemodialysis (HD). We compared the performance of the proposed time-warping markers, together with other previously proposed $$[K^{+}]$$ [ K + ] markers, such as T-wave width ( $$T_w$$ T w ) and T-wave slope-to-amplitude ratio ( $$T_{S/A}$$ T S / A ), when computed from standard ECG leads as well as from principal component analysis (PCA)-based leads. 48-hour ECG recordings and a set of hourly-collected blood samples from 29 ESRD-HD patients were acquired. Values of $$d_w$$ d w , $$d_a$$ d a , $${d_w^{\mathrm{NL}}}$$ d w NL , $${d_a^{\mathrm{NL}}}$$ d a NL and $$d_{w,c}$$ d w , c were calculated by comparing the morphology of the mean warped T-waves (MWTWs) derived at each hour along the HD with that from a reference MWTW, measured at the end of the HD. From the same MWTWs $$T_w$$ T w and $$T_{S/A}$$ T S / A were also extracted. Similarly, $$\Delta [K^+]$$ Δ [ K + ] was calculated as the difference between the $$[K^{+}]$$ [ K + ] values at each hour and the $$[K^{+}]$$ [ K + ] reference level at the end of the HD session. We found that $$d_{w}$$ d w and $$d_{w,c}$$ d w , c showed higher correlation coefficients with $$\Delta [K^+]$$ Δ [ K + ] than $$T_{S/A}$$ T S / A —Spearman’s ( $$\rho$$ ρ ) and Pearson’s (r)—and $$T_w$$ T w —Spearman’s ( $$\rho$$ ρ )—in both SL and PCA approaches being the intra-patient median $$\rho \ge 0.82$$ ρ ≥ 0.82 and $$r \ge 0.87$$ r ≥ 0.87 in SL and $$\rho \ge 0.82$$ ρ ≥ 0.82 and $$r \ge 0.89$$ r ≥ 0.89 in PCA respectively. Our findings would point at $$d_{w}$$ d w and $$d_{w,c}$$ d w , c as the most suitable surrogate of $$\Delta [K^+]$$ Δ [ K + ] , suggesting that they could be potentially useful for non-invasive monitoring of ESRD-HD patients in hospital, as well as in ambulatory settings. Therefore, the tracking of T-wave morphology variations by means of time-warping analysis could improve continuous and remote $$[K^{+}]$$ [ K + ] monitoring of ESRD-HD patients and flagging risk of $$[K^{+}]$$ [ K + ] -related cardiovascular events.

Summary: Age-related fibrosis in the left ventricle (LV) has been mainly studied in animals by assessing collagen content. Using second-harmonic generation microscopy and image processing, we evaluated amount, aggregation and spatial distribution of LV collagen in young to old pigs, and middle-age and elder living donors. All collagen features increased when comparing adult and old pigs with young ones, but not when comparing adult with old pigs or middle-age with elder individuals. Remarkably, all collagen parameters strongly correlated with lipofuscin, a biological age marker, in humans. By building patient-specific models of human ventricular tissue electrophysiology, we confirmed that amount and organization of fibrosis modulated arrhythmia vulnerability, and that distribution should be accounted for arrhythmia risk assessment. In conclusion, we characterize the age-associated changes in LV collagen and its potential implications for ventricular arrhythmia development. Consistency between pig and human results substantiate the pig as a relevant model of age-related LV collagen dynamics.

Background: The health benefits of sports and exercise training are well known. However, an acute bout of exercise transiently increases the risk of sudden cardiac death (SCD). To minimize the cardiovascular risks of exercise, more insight into the prevention and causes of SCD is needed. Methods: The observatory for the prevention of sudden death in sports, PREMUBID, was created with the aim of fostering research to assess the benefits and risks of exercise at different volumes and intensities and to get insight into the underlying mechanisms of potentially cardiac (mal) adaptations. Results: The observatory gathers researchers from a wide range of disciplines working at institutions in Europe and North America. The guiding principles of PREMUBID are to broaden the understanding of SCD in sports, strengthening collaborative research across the globe, and to develop, implement and evaluate robust pre-participation screening and emergency care strategies to further reduce the number of fatal cardiac events in sport events. During the inaugural meeting of the observatory, members and affiliated researchers discussed possibilities to initiate collaborative research projects and to exchange staff and students to share information and practices to prevent SCD. The final goal is to translate the obtained knowledge to understandable messages for the general population and healthcare workers to ensure that the population at large benefits from it. Conclusions: The PREMUBID consortium aims to produce novel knowledge and insights in SCD prevention, in order to maximize the health benefits associated with acute and long-term exercise training.

The autonomic nervous system (ANS) plays an essential role in the generation and maintenance of cardiac arrhythmias. The cardiac ANS can be divided into its extrinsic and intrinsic components, with the latter being organized in an epicardial neural network of interconnecting axons and clusters of autonomic ganglia called ganglionated plexi (GPs). GP ablation has been associated with a decreased risk of atrial fibrillation (AF) recurrence, but the accurate location of GPs is required for ablation to be effective. Although GP stimulation triggers both sympathetic and parasympathetic ANS branches, a predominance of parasympathetic activity has been shown. This study aims was to develop a method to locate atrial parasympathetic innervation sites based on measurements from a grid of electrograms (EGMs). Electrophysiological models representative of non-AF, paroxysmal AF (PxAF), and persistent AF (PsAF) tissues were developed. Parasympathetic effects were modeled by increasing the concentration of the neurotransmitter acetylcholine (ACh) in randomly distributed circles across the tissue. Different circle sizes of ACh and fibrosis geometries were considered, accounting for both uniform diffuse and non-uniform diffuse fibrosis. Computational simulations were performed, from which unipolar EGMs were computed in a 16 × 1 6 electrode mesh. Different distances of the electrodes to the tissue (0.5, 1, and 2 mm) and noise levels with signal-to-noise ratio (SNR) values of 0, 5, 10, 15, and 20 dB were tested. The amplitude of the atrial EGM repolarization wave was found to be representative of the presence or absence of ACh release sites, with larger positive amplitudes indicating that the electrode was placed over an ACh region. Statistical analysis was performed to identify the optimal thresholds for the identification of ACh sites. In all non-AF, PxAF, and PsAF tissues, the repolarization amplitude rendered successful identification. The algorithm performed better in the absence of fibrosis or when fibrosis was uniformly diffuse, with a mean accuracy of 0.94 in contrast with a mean accuracy of 0.89 for non-uniform diffuse fibrotic cases. The algorithm was robust against noise and worked for the tested ranges of electrode-to-tissue distance. In conclusion, the results from this study support the feasibility to locate atrial parasympathetic innervation sites from the amplitude of repolarization wave.

Computational models of cardiac electrophysiology are promising tools for reducing the rates of non-response patients suitable for cardiac resynchronization therapy (CRT) by optimizing electrode placement. The majority of computational models in the literature are mesh-based, primarily using the finite element method (FEM). The generation of patient-specific cardiac meshes has traditionally been a tedious task requiring manual intervention and hindering the modeling of a large number of cases. Meshless models can be a valid alternative due to their mesh quality independence. The organization of challenges such as the CRT-EPiggy19, providing unique experimental data as open access, enables benchmarking analysis of different cardiac computational modeling solutions with quantitative metrics. We present a benchmark analysis of a meshless-based method with finite-element methods for the prediction of cardiac electrical patterns in CRT, based on a subset of the CRT-EPiggy19 dataset. A data assimilation strategy was designed to personalize the most relevant parameters of the electrophysiological simulations and identify the optimal CRT lead configuration. The simulation results obtained with the meshless model were equivalent to FEM, with the most relevant aspect for accurate CRT predictions being the parameter personalization strategy (e.g., regional conduction velocity distribution, including the Purkinje system and CRT lead distribution).

Objective: Non-invasive estimation of serum potassium, [K+], and calcium, [Ca2+], can help to prevent life-threatening ventricular arrhythmias in patients with advanced renal disease, but current methods for estimation of electrolyte levels have limitations. We aimed to develop new markers based on the morphology of the QRS complex of the electrocardiogram (ECG). Methods: ECG recordings from 29 patients undergoing hemodialysis (HD) were processed. Mean warped QRS complexes were computed in two-minute windows at the start of an HD session, at the end of each HD hour and 48 h after it. We quantified QRS width, amplitude and the proposed QRS morphology-based markers that were computed by warping techniques. Reference [K+] and [Ca2+] were determined from blood samples acquired at the time points where the markers were estimated. Linear regression models were used to estimate electrolyte levels from the QRS markers individually and in combination with T wave morphology markers. Leave-one-out cross-validation was used to assess the performance of the estimators. Results: All markers, except for QRS width, strongly correlated with [K+] (median Pearson correlation coefficients, r, ranging from 0.81 to 0.87) and with [Ca2+] (r ranging from 0.61 to 0.76). QRS morphology markers showed very low sensitivity to heart rate (HR). Actual and estimated serum electrolyte levels differed, on average, by less than 0.035 mM (relative error of 0.018) for [K+] and 0.010 mM (relative error of 0.004) for [Ca2+] when patient-specific multivariable estimators combining QRS and T wave markers were used. Conclusion: QRS morphological markers allow non-invasive estimation of [K+] and [Ca2+] with low sensitivity to HR. The estimation performance is improved when multivariable models, including T wave markers, are considered. Significance: Markers based on the QRS complex of the ECG could contribute to non-invasive monitoring of serum electrolyte levels and arrhythmia risk prediction in patients with renal disease.

Centenarians are the paradigm of human extreme longevity and healthy aging, because they have postponed, if not avoided, mayor age-related diseases. The purpose of this study was to investigate potential differences in resting heart rate variability (HRV) between young adults, octogenarians, and centenarians and assess whether HRV variables are predictors of all-cause mortality in centenarians. To this end, three groups of participants: young adults (N = 20; 20.6 ± 2.3 years), octogenarians (N = 18; 84.1 ± 2.6 years), and centenarians (N = 17; 101.9 ± 1.9 years) were monitored for 15 min at rest (seated, without moving or talking) to measure RR intervals, from which HRV was evaluated. Our results showed a clear decrease with age in the main parasympathetic HRV variables, as well as in the standard deviation (SD) of the RR series [SD of normal-to-normal interval (SDNN)] and in low frequency (LF) heart rate (HR) oscillations, although differences between octogenarians and centenarians did not reach statistical significance. In 14 centenarians followed until death, only SDNN showed significant correlation (ρ = 0.536; p = 0.048) with survival prognosis. Additionally, SDNN

Ion channels have recently attracted attention as potential mediators of skin disease. Here, we explored the consequences of genetically encoded induction of the cell volume-regulating Ca2+-activated KCa3.1 channel (Kcnn4) for murine epidermal homeostasis. Doxycycline-treated mice harboring the KCa3.1+-transgene under the control of the reverse tetracycline-sensitive transactivator (rtTA) showed 800-fold channel overexpression above basal levels in the skin and solid KCa3.1-currents in keratinocytes. This overexpression resulted in epidermal spongiosis, progressive epidermal hyperplasia and hyperkeratosis, itch and ulcers. The condition was accompanied by production of the pro-proliferative and pro-inflammatory cytokines, IL-β1 (60-fold), IL-6 (33-fold), and TNFα (26-fold) in the skin. Treatment of mice with the KCa3.1-selective blocker, Senicapoc, significantly suppressed spongiosis and hyperplasia, as well as induction of IL-β1 (-88%) and IL-6 (-90%). In conclusion, KCa3.1-induction in the epidermis caused expression of pro-proliferative cytokines leading to spongiosis, hyperplasia and hyperkeratosis. This skin condition resembles pathological features of eczematous dermatitis and identifies KCa3.1 as a regulator of epidermal homeostasis and spongiosis, and as a potential therapeutic target.

Background and Objectives: Recent studies in humans and dogs have shown that ventricular repolarization exhibits a low-frequency (LF) oscillatory pattern following enhanced sympathetic activity, which has been related to arrhythmic risk. The appearance of LF oscillations in ventricular repolarization is, however, not immediate, but it may take up to some minutes. This study seeks to characterize the time course of the action potential (AP) duration (APD) oscillatory behavior in response to sympathetic provocations, unveil its underlying mechanisms and establish a potential link to arrhythmogenesis under disease conditions.Materials and Methods: A representative set of human ventricular computational models coupling cellular electrophysiology, calcium dynamics, β-adrenergic signaling, and mechanics was built. Sympathetic provocation was modeled via phasic changes in β-adrenergic stimulation (β-AS) and mechanical stretch at Mayer wave frequencies within the 0.03–0.15 Hz band.Results: Our results show that there are large inter-individual differences in the time lapse for the development of LF oscillations in APD following sympathetic provocation, with some cells requiring just a few seconds and other cells needing more than 3 min. Whereas, the oscillatory response to phasic mechanical stretch is almost immediate, the response to β-AS is much more prolonged, in line with experimentally reported evidences, thus being this component the one driving the slow development of APD oscillations following enhanced sympathetic activity. If β-adrenoceptors are priorly stimulated, the time for APD oscillations to become apparent is remarkably reduced, with the oscillation time lapse being an exponential function of the pre-stimulation level. The major mechanism underlying the delay in APD oscillations appearance is related to the slow IKs phosphorylation kinetics, with its relevance being modulated by the IKs conductance of each individual cell. Cells presenting short oscillation time lapses are commonly associated with large APD oscillation magnitudes, which facilitate the occurrence of pro-arrhythmic events under disease conditions involving calcium overload and reduced repolarization reserve.Conclusions: The time course of LF oscillatory behavior of APD in response to increased sympathetic activity presents high inter-individual variability, which is associated with different expression and PKA phosphorylation kinetics of the IKs current. Short time lapses in the development of APD oscillations are associated with large oscillatory magnitudes and pro-arrhythmic risk under disease conditions.

BackgroundRecent clinical, experimental and modeling studies link oscillations of ventricular repolarization in the low frequency (LF) (approx. 0.1 Hz) to arrhythmogenesis. Sympathetic provocation has been shown to enhance both LF oscillations of action potential duration (APD) and beat-to-beat variability (BVR) in humans. We hypothesized that beta-adrenergic blockade would reduce LF oscillations of APD and BVR of APD in humans and that the two processes might be linked.Methods and ResultsTwelve patients with normal ventricles were studied during routine electrophysiological procedures. Activation-recovery intervals (ARI) as a conventional surrogate for APD were recorded from 10 left and 10 right ventricular endocardial sites before and after acute beta-adrenergic adrenergic blockade. Cycle length was maintained constant with right ventricular pacing. Oscillatory behavior of ARI was quantified by spectral analysis and BVR as the short-term variability. Beta-adrenergic blockade reduced LF ARI oscillations (8.6 ± 4.5 ms2 vs. 5.5 ± 3.5 ms2, p = 0.027). A significant correlation was present between the initial control values and reduction seen following beta-adrenergic blockade in LF ARI (rs = 0.62, p = 0.037) such that when initial values are high the effect is greater. A similar relationship was also seen in the beat-to beat variability of ARI (rs = 0.74, p = 0.008). There was a significant correlation between the beta-adrenergic blockade induced reduction in LF power of ARI and the witnessed reduction of beat-to-beat variability of ARI (rs = 0.74, p = 0.01). These clinical results accord with recent computational modeling studies which provide mechanistic insight into the interactions of LF oscillations and beat-to-beat variability of APD at the cellular level.ConclusionBeta-adrenergic blockade reduces LF oscillatory behavior of APD (ARI) in humans in vivo. Our results support the importance of LF oscillations in modulating the response of BVR to beta-adrenergic blockers, suggesting that LF oscillations may play role in modulating beta-adrenergic mechanisms underlying BVR.

Ventricular arrhythmias and sudden cardiac death during long-term space missions are a major concern for space agencies. Long-duration spaceflight and its ground-based analog head-down bed rest (HDBR) have been reported to markedly alter autonomic and cardiac functioning, particularly affecting ventricular repolarization of the electrocardiogram (ECG). In this study, novel methods are developed, departing from previously published methodologies, to quantify the index of Periodic Repolarization Dynamics (PRD), an arrhythmic risk marker that characterizes sympathetically-mediated low-frequency oscillations in the T-wave vector. PRD is evaluated in ECGs from 42 volunteers at rest and during an orthostatic tilt table test recorded before and after 60-day –6° HDBR. Our results indicate that tilt test, on top of enhancing sympathetic regulation of heart rate, notably increases PRD, both before and after HDBR, thus supporting previous evidence on PRD being an indicator of sympathetic modulation of ventricular repolarization. Importantly, long-term microgravity exposure is shown to lead to significant increases in PRD, both when evaluated at rest and, even more notably, in response to tilt test. The extent of microgravity-induced changes in PRD has been associated with arrhythmic risk in prior studies. An exercise-based, but not a nutrition-based, countermeasure is able to partially reverse microgravity-induced effects on PRD. In conclusion, long-term exposure to microgravity conditions leads to elevated low-frequency oscillations of ventricular repolarization, which are potentiated following sympathetic stimulation and are related to increased risk for repolarization instabilities and arrhythmias. Tested countermeasures are only partially effective in counteracting microgravity effects.

Background and Objectives: Enhanced beat-to-beat variability of ventricular repolarization (BVR) has been linked to arrhythmias and sudden cardiac death. Recent experimental studies on human left ventricular epicardial electrograms have shown that BVR closely interacts with low-frequency (LF) oscillations of activation recovery interval during sympathetic provocation. In this work human ventricular computational cell models are developed to reproduce the experimentally observed interactions between BVR and its LF oscillations, to assess underlying mechanisms and to establish a relationship with arrhythmic risk.Materials and Methods: A set of human ventricular action potential (AP) models covering a range of experimental electrophysiological characteristics was constructed. These models incorporated stochasticity in major ionic currents as well as descriptions of β-adrenergic stimulation and mechanical effects to investigate the AP response to enhanced sympathetic activity. Statistical methods based on Automatic Relevance Determination and Canonical Correlation Analysis were developed to unravel individual and common factors contributing to BVR and LF patterning of APD in response to sympathetic provocation.Results: Simulated results reproduced experimental evidences on the interactions between BVR and LF oscillations of AP duration (APD), with replication of the high inter-individual variability observed in both phenomena. ICaL, IKr and IK1 currents were identified as common ionic modulators of the inter-individual differences in BVR and LF oscillatory behavior and were shown to be crucial in determining susceptibility to arrhythmogenic events.Conclusions: The calibrated family of human ventricular cell models proposed in this study allows reproducing experimentally reported interactions between BVR and LF oscillations of APD. Ionic factors involving ICaL, IKr and IK1 currents are found to underlie correlated increments in both phenomena in response to sympathetic provocation. A link to arrhythmogenesis is established for concomitantly elevated levels of BVR and its LF oscillations.

Background: End-stage renal disease patients undergoing hemodialysis (ESRD-HD) therapy are highly susceptible to malignant ventricular arrhythmias caused by undetected potassium concentration ([K+]) variations (Δ[K+]) out of normal ranges. Therefore, a reliable method for continuous, noninvasive monitoring of [K+] is crucial. The morphology of the T-wave in the electrocardiogram (ECG) reflects Δ[K+] and two time-warping-based T-wave morphological parameters, dw and its heart-rate corrected version dw,c, have been shown to reliably track Δ[K+] from the ECG. The aim of this study is to derive polynomial models relating dw and dw,c with Δ[K+], and to test their ability to reliably sense and quantify Δ[K+] values. Methods: 48-hour Holter ECGs and [K+] values from six blood samples were collected from 29 ESRD-HD patients. For every patient, dw and dw,c were computed, and linear, quadratic, and cubic fitting models were derived from them. Then, Spearman’s (ρ) and Pearson’s (r) correlation coefficients, and the estimation error (ed) between Δ[K+] and the corresponding model-estimated values (Δ^[K+]) were calculated. Results and Discussions: Nonlinear models were the most suitable for Δ[K+] estimation, rendering higher Pearson’s correlation (median 0.77 ≤r≤ 0.92) and smaller estimation error (median 0.20 ≤ed≤ 0.43) than the linear model (median 0.76 ≤r≤ 0.86 and 0.30 ≤ed≤ 0.40), even if similar Spearman’s ρ were found across models (median 0.77 ≤ρ≤ 0.83). Conclusion: Results support the use of nonlinear T-wave-based models as Δ[K+] sensors in ESRD-HD patients.

This work aims to validate the Polar H7 heart rate (HR) sensor for heart rate variability (HRV) analysis at rest and during various exercise intensities in a cohort of male volunteers with different age, body composition and fitness level. Cluster analysis was carried out to evaluate how these phenotypic characteristics influenced HR and HRV measurements. For this purpose, sixty-seven volunteers performed a test consisting of the following consecutive segments: sitting rest, three submaximal exercise intensities in cycle-ergometer and sitting recovery. The agreement between HRV indices derived from Polar H7 and a simultaneous electrocardiogram (ECG) was assessed using concordance correlation coefficient (CCC). The percentage of subjects not reaching excellent agreement (CCC > 0.90) was higher for high-frequency power (PHF) than for low-frequency power (PLF) of HRV and increased with exercise intensity. A cluster of unfit and not young volunteers with high trunk fat percentage showed the highest error in HRV indices. This study indicates that Polar H7 and ECG were interchangeable at rest. During exercise, HR and PLF showed excellent agreement between devices. However, during the highest exercise intensity, CCC for PHF was lower than 0.90 in as many as 60% of the volunteers. During recovery, HR but not HRV measurements were accurate. As a conclusion, phenotypic differences between subjects can represent one of the causes for disagreement between HR sensors and ECG devices, which should be considered specifically when using Polar H7 and, generally, in the validation of any HR sensor for HRV analysis.

Cardiomyocytes’ geometry and connexin 43 (CX43) amount and distribution are structural features that play a pivotal role in electrical conduction. Their quantitative assessment is of high interest in the study of arrhythmias, but it is usually hampered by the lack of automatic tools. In this work, we propose a software algorithm (Myocyte Automatic Retrieval and Tissue Analyzer, MARTA) to automatically detect myocytes from fluorescent microscopy images of cardiac tissue, measure their morphological features and evaluate the expression of CX43 and its degree of lateralization. The proposed software is based on the generation of cell masks, contouring of individual cells, enclosing of cells in minimum area rectangles and splitting of these rectangles into end-to-end and middle compartments to estimate CX43 lateral-to-total ratio. Application to human ventricular tissue images shows that mean differences between automatic and manual methods in terms of cardiomyocyte length and width are below 4 μm. The percentage of lateral CX43 also agrees between automatic and manual evaluation, with the interquartile range approximately covering from 3% to 30% in both cases. MARTA is not limited by fiber orientation and has an optimized speed by using contour filtering, which makes it run hundreds of times faster than a trained expert. Developed for CX43 studies in the left ventricle, MARTA is a flexible tool applicable to morphometric and lateralization studies of other markers in any heart chamber or even skeletal muscle. This open-access software is available online.

In this study, we teased out the dynamical mechanisms underlying the generation of arrhythmogenic early afterdepolarizations (EADs) in a three-variable model of a mammalian ventricular cell. Based on recently published studies, we consider a 1-fast, 2-slow variable decomposition of the system describing the cellular action potential. We use sweeping techniques, such as the spike-counting method, and bifurcation and continuation methods to identify parametric regions with EADs. We show the existence of isolas of periodic orbits organizing the different EAD patterns and we provide a preliminary classification of our fast–slow decomposition according to the involved dynamical phenomena. This investigation represents a basis for further studies into the organization of EAD patterns in the parameter space and the involved bifurcations.

[This corrects the article DOI: 10.1371/journal.pone.0163619.].

TRPV4 channels are calcium-permeable cation channels that are activated by several physicochemical stimuli. Accordingly, TRPV4 channels have been implicated in the regulation of osmosensing, mechanotransduction, thermosensation, and epithelial/endothelial barrier functions. Whether TRPV4 is also mechanistically implicated in melanoma cell proliferation is not clear. Here, we hypothesized that TRPV4 is expressed in human melanoma and that pharmacological activation interferes with cell proliferation.TRPV4 functions were studied in melanoma cell lines (A375, SK-MEL-28, MKTBR), immortalized non-cancer keratinocytes (HaCaT), and murine 3T3 fibroblasts by patch-clamp, qRT-PCR, intracellular calcium measurements, cell proliferation, and flow cytometric assays of apoptosis and cell cycle. The selective TRPV4-activator, GSK1016790A, elicited non-selective cation currents with TRPV4-typical current-voltage-relationship in all cell lines. GSK1016790A-induced currents were blocked by the TRPV4-blocker, HC067047. TRPV4 mRNA expression was demonstrated by qRT-PCR. In A375 cells, TRPV4 activation was frequently paralleled by co-activation of calcium/calmodulin-regulated KCa3.1 channels. Light microscopy showed that TRPV4-activation produced rapid cellular disarrangement, nuclear densification, and detachment of a large fraction of all melanoma cell lines and HaCaT cells. TRPV4-activation induced apoptosis and drastically inhibited A375 and HaCaT proliferation that could be partially prevented by HC067047.Our study showed that TRPV4 channels were functionally expressed in human melanoma cell lines and in human keratinocytes. Pharmacological TRPV4 activation in human melanoma cells and keratinocytes caused severe cellular disarrangement, necrosis and apoptosis. Pharmacological targeting of TRPV4 could be an alternative or adjuvant therapeutic strategy to treat melanoma progression and other proliferative skin disorders.

Background. Nonlinear heart rate variability (HRV) indices have extended the description of autonomic nervous system (ANS) regulation of the heart. One of those indices is approximate entropy, ApEn, which has become a commonly used measure of the irregularity of a time series. To calculate ApEn, a priori definition of parameters like the threshold on similarity and the embedding dimension is required, which has been shown to be critical for interpretation of the results. Thus, searching for a parameter-free ApEn-based index could be advantageous for standardizing the use and interpretation of this widely applied entropy measurement. Methods. A novel entropy index called multidimensional approximate entropy, MApEnmax, is proposed based on summing the contribution of maximum approximate entropies over a wide range of embedding dimensions while selecting the similarity threshold leading to maximum ApEn value in each dimension. Synthetic RR interval time series with varying levels of stochasticity, generated by both MIX(P) processes and white/pink noise, were used to validate the properties of the proposed index. Aging and congestive heart failure (CHF) were characterized from RR interval time series of available databases. Results. In synthetic time series, MApEnmax values were proportional to the level of randomness; i.e., MApEnmax increased for higher values of P in generated MIX(P) processes and was larger for white than for pink noise. This result was a consequence of all maximum approximate entropy values being increased for higher levels of randomness in all considered embedding dimensions. This is in contrast to the results obtained for approximate entropies computed with a fixed similarity threshold, which presented inconsistent results for different embedding dimensions. Evaluation of the proposed index on available databases revealed that aging was associated with a notable reduction in MApEnmax values. On the other hand, MApEnmax evaluated during the night period was considerably larger in CHF patients than in healthy subjects. Conclusion. A novel parameter-free multidimensional approximate entropy index, MApEnmax, is proposed and tested over synthetic data to confirm its capacity to represent a range of randomness levels in HRV time series. MApEnmax values are reduced in elderly patients, which may correspond to the reported loss of ANS adaptability in this population segment. Increased MApEnmax values measured in CHF patients versus healthy subjects during the night period point to greater irregularity of heart rate dynamics caused by the disease.

Cardiac electrophysiological computational models are often developed from previously published models. The new models may incorporate additional features to adapt the model to a different species or may upgrade a specific ionic formulation based on newly available experimental data. A relevant challenge in the development of a new model is the estimation of certain ionic current conductances that cannot be reliably identified from experiments. A common strategy to estimate those conductances is by means of constrained non-linear least-squares optimization. In this work, a novel methodology is proposed for estimation of ionic current conductances of cardiac electrophysiological models by using a response surface approximation-based constrained optimization with trust region management. Polynomial response surfaces of a number of electrophysiological markers were built using statistical sampling methods. These markers included action potential duration (APD), triangulation, diastolic and systolic intracellular calcium concentration, and time constants of APD rate adaptation. The proposed methodology was applied to update the Carro et al. human ventricular action potential model after incorporation of intracellular potassium ([K+]i) dynamics. While the Carro et al. model was well suited for investigation of arrhythmogenesis, it did not allow simulation of [K+]i changes. With the methodology proposed in this study, the updated Carro et al. human ventricular model could be used to simulate [K+]i changes in response to varying extracellular potassium ([K+]o) levels. Additionally, it rendered values of evaluated electrophysiological markers within physiologically plausible ranges. The optimal values of ionic current conductances in the updated model were found in a notably shorter time than with previously proposed methodologies. As a conclusion, the response surface optimization-based approach proposed in this study allows estimating ionic current conductances of cardiac electrophysiological computational models while guaranteeing replication of key electrophysiological features and with an important reduction in computational cost with respect to previously published approaches. The updated Carro et al. model developed in this study is thus suitable for the investigation of arrhythmic risk-related conditions, including those involving large changes in potassium concentration.

Atrial fibrillation (AF) usually manifests as reentrant circuits propagating through the whole atria creating chaotic activation patterns. Little is yet known about how differences in electrophysiological and ionic properties between patients modulate reentrant patterns in AF. The goal of this study is to quantify how variability in action potential duration (APD) at different stages of repolarization determines AF dynamics and their modulation by ionic block using a set of virtual whole-atria human models. Six human whole-atria models are constructed based on the same anatomical structure and fiber orientation, but with different electrophysiological phenotypes. Membrane kinetics for each whole-atria model are selected with distinct APD characteristics at 20, 50, and 90% repolarization, from an experimentally calibrated population of human atrial action potential models, including AF remodeling and acetylcholine parasympathetic effects. Our simulations show that in all whole-atria models, reentrant circuits tend to organize around the pulmonary veins and the right atrial appendage, thus leading to higher dominant frequency (DF) and more organized activation in the left atrium than in the right atrium. Differences in APD in all phases of repolarization (not only APD90) yielded quantitative differences in fibrillation patterns with long APD associated with slower and more regular dynamics. Long APD50 and APD20 were associated with increased interatrial conduction block and interatrial differences in DF and organization index, creating reentry instability and self-termination in some cases. Specific inhibitions of IK1, INaK, or INa reduce DF and organization of the arrhythmia by enlarging wave meandering, reducing the number of secondary wavelets, and promoting interatrial block in all six virtual patients, especially for the phenotypes with short APD at 20, 50, and/or 90% repolarization. This suggests that therapies aiming at prolonging the early phase of repolarization might constitute effective antiarrhythmic strategies for the pharmacological management of AF. In summary, simulations report significant differences in atrial fibrillatory dynamics resulting from differences in APD at all phases of repolarization.

Background Patients with chronic heart failure are at high risk of sudden cardiac death (SCD). Increased dispersion of repolarization restitution has been associated with SCD, and we hypothesize that this should be reflected in the morphology of the T‐wave and its variations with heart rate. The aim of this study is to propose an electrocardiogram (ECG)‐based index characterizing T‐wave morphology restitution (TMR), and to assess its association with SCD risk in a population of chronic heart failure patients. Methods and Results Holter ECGs from 651 ambulatory patients with chronic heart failure from the MUSIC (MUerte Súbita en Insuficiencia Cardiaca) study were available for the analysis. TMR was quantified by measuring the morphological variation of the T‐wave per RR increment using time‐warping metrics, and its predictive power was compared to that of clinical variables such as the left ventricular ejection fraction and other ECG‐derived indices, such as T‐wave alternans and heart rate variability. TMR was significantly higher in SCD victims than in the rest of patients (median 0.046 versus 0.039, P

Sudden cardiac death (SCD) and pump failure death (PFD) are common endpoints in chronic heart failure (CHF) patients, but prevention strategies are different. Currently used tools to specifically predict these endpoints are limited. We developed risk models to specifically assess SCD and PFD risk in CHF by combining ECG markers and clinical variables.The relation of clinical and ECG markers with SCD and PFD risk was assessed in 597 patients enrolled in the MUSIC (MUerte Súbita en Insuficiencia Cardiaca) study. ECG indices included: turbulence slope (TS), reflecting autonomic dysfunction; T-wave alternans (TWA), reflecting ventricular repolarization instability; and T-peak-to-end restitution (ΔαTpe) and T-wave morphology restitution (TMR), both reflecting changes in dispersion of repolarization due to heart rate changes. Standard clinical indices were also included.The indices with the greatest SCD prognostic impact were gender, New York Heart Association (NYHA) class, left ventricular ejection fraction, TWA, ΔαTpe and TMR. For PFD, the indices were diabetes, NYHA class, ΔαTpe and TS. Using a model with only clinical variables, the hazard ratios (HRs) for SCD and PFD for patients in the high-risk group (fifth quintile of risk score) with respect to patients in the low-risk group (first and second quintiles of risk score) were both greater than 4. HRs for SCD and PFD increased to 9 and 11 when using a model including only ECG markers, and to 14 and 13, when combining clinical and ECG markers.The inclusion of ECG markers capturing complementary pro-arrhythmic and pump failure mechanisms into risk models based only on standard clinical variables substantially improves prediction of SCD and PFD in CHF patients.

The purpose of this study is to characterize and attenuate the influence of mean heart rate (HR) on nonlinear heart rate variability (HRV) indices (correlation dimension, sample, and approximate entropy) as a consequence of being the HR the intrinsic sampling rate of HRV signal. This influence can notably alter nonlinear HRV indices and lead to biased information regarding autonomic nervous system (ANS) modulation. First, a simulation study was carried out to characterize the dependence of nonlinear HRV indices on HR assuming similar ANS modulation. Second, two HR-correction approaches were proposed: one based on regression formulas and another one based on interpolating RR time series. Finally, standard and HR-corrected HRV indices were studied in a body position change database. The simulation study showed the HR-dependence of non-linear indices as a sampling rate effect, as well as the ability of the proposed HR-corrections to attenuate mean HR influence. Analysis in a body position changes database shows that correlation dimension was reduced around 21% in median values in standing with respect to supine position (p < 0.05), concomitant with a 28% increase in mean HR (p < 0.05). After HR-correction, correlation dimension decreased around 18% in standing with respect to supine position, being the decrease still significant. Sample and approximate entropy showed similar trends. HR-corrected nonlinear HRV indices could represent an improvement in their applicability as markers of ANS modulation when mean HR changes.

Drugs that prolong the heart rate corrected QT interval (QTc) on the electrocardiogram (ECG) by blocking the hERG potassium channel and also block inward currents (late sodium or L-type calcium) are not associated with torsade de pointes (e.g. ranolazine and verapamil). Thus, identifying ECG signs of late sodium current block could aid in the determination of proarrhythmic risk for new drugs. A new cardiac safety paradigm for drug development (the "CiPA" initiative) will involve the preclinical assessment of multiple human cardiac ion channels and ECG biomarkers are needed to determine if there are unexpected ion channel effects in humans.In this study we assess the ability of eight ECG morphology biomarkers to detect late sodium current block in the presence of QTc prolongation by analyzing a clinical trial where a selective hERG potassium channel blocker (dofetilide) was administered alone and then in combination with two late sodium current blockers (lidocaine and mexiletine). We demonstrate that late sodium current block has the greatest effect on the heart-rate corrected J-Tpeak interval (J-Tpeakc), followed by QTc and then T-wave flatness. Furthermore, J-Tpeakc is the only biomarker that improves detection of the presence of late sodium current block compared to using QTc alone (AUC: 0.83 vs. 0.72 respectively, p

AIMS:Human atrial electrophysiology exhibits high inter-subject variability in both sinus rhythm (SR) and chronic atrial fibrillation (cAF) patients. Variability is however rarely investigated in experimental and theoretical electrophysiological studies, thus hampering the understanding of its underlying causes but also its implications in explaining differences in the response to disease and treatment. In our study, we aim at investigating the ability of populations of human atrial cell models to capture the inter-subject variability in action potential (AP) recorded in 363 patients both under SR and cAF conditions. METHODS AND RESULTS:Human AP recordings in atrial trabeculae (n = 469) from SR and cAF patients were used to calibrate populations of computational SR and cAF atrial AP models. Three populations of over 2000 sampled models were generated, based on three different human atrial AP models. Experimental calibration selected populations of AP models yielding AP with morphology and duration in range with experimental recordings. Populations using the three original models can mimic variability in experimental AP in both SR and cAF, with median conductance values in SR for most ionic currents deviating less than 30% from their original peak values. All cAF populations show similar variations in G(K1), G(Kur) and G(to), consistent with AF-related remodeling as reported in experiments. In all SR and cAF model populations, inter-subject variability in I(K1) and I(NaK) underlies variability in APD90, variability in I(Kur), I(CaL) and I(NaK) modulates variability in APD50 and combined variability in Ito and I(Kur) determines variability in APD20. The large variability in human atrial AP triangulation is mostly determined by I(K1) and either I(NaK) or I(NaCa) depending on the model. CONCLUSION:Experimentally-calibrated human atrial AP models populations mimic AP variability in SR and cAF patient recordings, and identify potential ionic determinants of inter-subject variability in human atrial AP duration and morphology in SR versus cAF.