Differential sensitivity to longitudinal and transverse stretch mediates transcriptional responses in mouse neonatal ventricular myocytes.

To identify how cardiomyocyte mechanosensitive signaling pathways are regulated by anisotropic stretch, micropatterned mouse neonatal cardiomyocytes were stretched primarily longitudinally or transversely to the myofiber axis. Four hours of static, longitudinal stretch induced differential expression of 557 genes, compared with 30 induced by transverse stretch, measured using RNA-seq. A logic-based ordinary differential equation model of the cardiac myocyte mechanosignaling network, extended to include the transcriptional regulation and expression of 784 genes, correctly predicted measured expression changes due to anisotropic stretch with 69% accuracy. The model also predicted published transcriptional responses to mechanical load in vitro or in vivo with 63-91% accuracy. The observed differences between transverse and longitudinal stretch responses were not explained by differential activation of specific pathways but rather by an approximately twofold greater sensitivity to longitudinal stretch than transverse stretch. In vitro experiments confirmed model predictions that stretch-induced gene expression is more sensitive to angiotensin II and endothelin-1, via RhoA and MAP kinases, than to the three membrane ion channels upstream of calcium signaling in the network. Quantitative cardiomyocyte gene expression differs substantially with the axis of maximum principal stretch relative to the myofilament axis, but this difference is due primarily to differences in stretch sensitivity rather than to selective activation of mechanosignaling pathways.NEW & NOTEWORTHY Anisotropic stretch applied to micropatterned neonatal mouse ventricular myocytes induced markedly greater acute transcriptional responses when the major axis of stretch was parallel to the myofilament axis than when it was transverse. Analysis with a novel quantitative network model of mechanoregulated cardiomyocyte gene expression suggests that this difference is explained by higher cell sensitivity to longitudinal loading than transverse loading than by the activation of differential signaling pathways.

Signaling network model of cardiomyocyte morphological changes in familial cardiomyopathy.

Familial cardiomyopathy is a precursor of heart failure and sudden cardiac death. Over the past several decades, researchers have discovered numerous gene mutations primarily in sarcomeric and cytoskeletal proteins causing two different disease phenotypes: hypertrophic (HCM) and dilated (DCM) cardiomyopathies. However, molecular mechanisms linking genotype to phenotype remain unclear. Here, we employ a systems approach by integrating experimental findings from preclinical studies (e.g., murine data) into a cohesive signaling network to scrutinize genotype to phenotype mechanisms. We developed an HCM/DCM signaling network model utilizing a logic-based differential equations approach and evaluated model performance in predicting experimental data from four contexts (HCM, DCM, pressure overload, and volume overload). The model has an overall prediction accuracy of 83.8%, with higher accuracy in the HCM context (90%) than DCM (75%). Global sensitivity analysis identifies key signaling reactions, with calcium-mediated myofilament force development and calcium-calmodulin kinase signaling ranking the highest. A structural revision analysis indicates potential missing interactions that primarily control calcium regulatory proteins, increasing model prediction accuracy. Combination pharmacotherapy analysis suggests that downregulation of signaling components such as calcium, titin and its associated proteins, growth factor receptors, ERK1/2, and PI3K-AKT could inhibit myocyte growth in HCM. In experiments with patient-specific iPSC-derived cardiomyocytes (MLP-W4R;MYH7-R723C iPSC-CMs), combined inhibition of ERK1/2 and PI3K-AKT rescued the HCM phenotype, as predicted by the model. In DCM, PI3K-AKT-NFAT downregulation combined with upregulation of Ras/ERK1/2 or titin or Gq protein could ameliorate cardiomyocyte morphology. The model results suggest that HCM mutations that increase active force through elevated calcium sensitivity could increase ERK activity and decrease eccentricity through parallel growth factors, Gq-mediated, and titin pathways. Moreover, the model simulated the influence of existing medications on cardiac growth in HCM and DCM contexts. This HCM/DCM signaling model demonstrates utility in investigating genotype to phenotype mechanisms in familial cardiomyopathy.

A deep learning approach for fully automated cardiac shape modeling in tetralogy of Fallot.

BACKGROUND: Cardiac shape modeling is a useful computational tool that has provided quantitative insights into the mechanisms underlying dysfunction in heart disease. The manual input and time required to make cardiac shape models, however, limits their clinical utility. Here we present an end-to-end pipeline that uses deep learning for automated view classification, slice selection, phase selection, anatomical landmark localization, and myocardial image segmentation for the automated generation of three-dimensional, biventricular shape models. With this approach, we aim to make cardiac shape modeling a more robust and broadly applicable tool that has processing times consistent with clinical workflows. METHODS: Cardiovascular magnetic resonance (CMR) images from a cohort of 123 patients with repaired tetralogy of Fallot (rTOF) from two internal sites were used to train and validate each step in the automated pipeline. The complete automated pipeline was tested using CMR images from a cohort of 12 rTOF patients from an internal site and 18 rTOF patients from an external site. Manually and automatically generated shape models from the test set were compared using Euclidean projection distances, global ventricular measurements, and atlas-based shape mode scores. RESULTS: The mean absolute error (MAE) between manually and automatically generated shape models in the test set was similar to the voxel resolution of the original CMR images for end-diastolic models (MAE = 1.9 ± 0.5 mm) and end-systolic models (MAE = 2.1 ± 0.7 mm). Global ventricular measurements computed from automated models were in good agreement with those computed from manual models. The average mean absolute difference in shape mode Z-score between manually and automatically generated models was 0.5 standard deviations for the first 20 modes of a reference statistical shape atlas. CONCLUSIONS: Using deep learning, accurate three-dimensional, biventricular shape models can be reliably created. This fully automated end-to-end approach dramatically reduces the manual input required to create shape models, thereby enabling the rapid analysis of large-scale datasets and the potential to deploy statistical atlas-based analyses in point-of-care clinical settings. Training data and networks are available from cardiacatlas.org.

Le Cœur en Sabot: shape associations with adverse events in repaired tetralogy of Fallot.

BACKGROUND: Maladaptive remodelling mechanisms occur in patients with repaired tetralogy of Fallot (rToF) resulting in a cycle of metabolic and structural changes. Biventricular shape analysis may indicate mechanisms associated with adverse events independent of pulmonary regurgitant volume index (PRVI). We aimed to determine novel remodelling patterns associated with adverse events in patients with rToF using shape and function analysis. METHODS: Biventricular shape and function were studied in 192 patients with rToF (median time from TOF repair to baseline evaluation 13.5 years). Linear discriminant analysis (LDA) and principal component analysis (PCA) were used to identify shape differences between patients with and without adverse events. Adverse events included death, arrhythmias, and cardiac arrest with median follow-up of 10 years. RESULTS: LDA and PCA showed that shape characteristics pertaining to adverse events included a more circular left ventricle (LV) (decreased eccentricity), dilated (increased sphericity) LV base, increased right ventricular (RV) apical sphericity, and decreased RV basal sphericity. Multivariate LDA showed that the optimal discriminative model included only RV apical ejection fraction and one PCA mode associated with a more circular and dilated LV base (AUC = 0.77). PRVI did not add value, and shape changes associated with increased PRVI were not predictive of adverse outcomes. CONCLUSION: Pathological remodelling patterns in patients with rToF are significantly associated with adverse events, independent of PRVI. Mechanisms related to incident events include LV basal dilation with a reduced RV apical ejection fraction.

Right-left ventricular shape variations in tetralogy of Fallot: associations with pulmonary regurgitation.

BACKGROUND: Relationships between right ventricular (RV) and left ventricular (LV) shape and function may be useful in determining optimal timing for pulmonary valve replacement in patients with repaired tetralogy of Fallot (rTOF). However, these are multivariate and difficult to quantify. We aimed to quantify variations in biventricular shape associated with pulmonary regurgitant volume (PRV) in rTOF using a biventricular atlas. METHODS: In this cross-sectional retrospective study, a biventricular shape model was customized to cardiovascular magnetic resonance (CMR) images from 88 rTOF patients (median age 16, inter-quartile range 11.8-24.3 years). Morphometric scores quantifying biventricular shape at end-diastole and end-systole were computed using principal component analysis. Multivariate linear regression was used to quantify biventricular shape associations with PRV, corrected for age, sex, height, and weight. Regional associations were confirmed by univariate correlations with distances and angles computed from the models, as well as global systolic strains computed from changes in arc length from end-diastole to end-systole. RESULTS: PRV was significantly associated with 5 biventricular morphometric scores, independent of covariates, and accounted for 12.3% of total shape variation (p < 0.05). Increasing PRV was associated with RV dilation and basal bulging, in conjunction with decreased LV septal-lateral dimension (LV flattening) and systolic septal motion towards the RV (all p < 0.05). Increased global RV radial, longitudinal, circumferential and LV radial systolic strains were significantly associated with increased PRV (all p < 0.05). CONCLUSION: A biventricular atlas of rTOF patients quantified multivariate relationships between left-right ventricular morphometry and wall motion with pulmonary regurgitation. Regional RV dilation, LV reduction, LV septal-lateral flattening and increased RV strain were all associated with increased pulmonary regurgitant volume. Morphometric scores provide simple metrics linking mechanisms for structural and functional alteration with important clinical indices.

Atlas-based methods for efficient characterization of patient-specific ventricular activation patterns.

AIMS: Ventricular activation patterns can aid clinical decision-making directly by providing spatial information on cardiac electrical activation or indirectly through derived clinical indices. The aim of this work was to derive an atlas of the major modes of variation of ventricular activation from model-predicted 3D bi-ventricular activation time distributions and to relate these modes to corresponding vectorcardiograms (VCGs). We investigated how the resulting dimensionality reduction can improve and accelerate the estimation of activation patterns from surface electrogram measurements. METHODS AND RESULTS: Atlases of activation time (AT) and VCGs were derived using principal component analysis on a dataset of simulated electrophysiology simulations computed on eight patient-specific bi-ventricular geometries. The atlases provided significant dimensionality reduction, and the modes of variation in the two atlases described similar features. Utility of the atlases was assessed by resolving clinical waveforms against them and the VCG atlas was able to accurately reconstruct the patient VCGs with fewer than 10 modes. A sensitivity analysis between the two atlases was performed by calculating a compact Jacobian. Finally, VCGs generated by varying AT atlas modes were compared with clinical VCGs to estimate patient-specific activation maps, and the resulting errors between the clinical and atlas-based VCGs were less than those from more computationally expensive method. CONCLUSION: Atlases of activation and VCGs represent a new method of identifying and relating the features of these high-dimensional signals that capture the major sources of variation between patients and may aid in identifying novel clinical indices of arrhythmia risk or therapeutic outcome.

Improved compressed sensing and super-resolution of cardiac diffusion MRI with structure-guided total variation.

PURPOSE: Structure-guided total variation is a recently introduced prior that allows reconstruction of images using knowledge of the location and orientation of edges in a reference image. In this work, we demonstrate the advantages of a variant of structure-guided total variation known as directional total variation (DTV), over traditional total variation (TV), in the context of compressed-sensing reconstruction and super-resolution. METHODS: We compared TV and DTV in retrospectively undersampled ex vivo diffusion tensor imaging and diffusion spectrum imaging data from healthy, sham, and hypertrophic rat hearts. RESULTS: In compressed sensing at an undersampling factor of 8, the RMS error of mean diffusivity and fractional anisotropy relative to the fully sampled ground truth were 44% and 20% lower in DTV compared with TV. In super-resolution, these values were 29% and 14%, respectively. Similarly, we observed improvements in helix angle, transverse angle, sheetlet elevation, and sheetlet azimuth. The RMS error of the diffusion kurtosis in the undersampled data relative to the ground truth was uniformly lower (22% on average) with DTV compared to TV. CONCLUSION: Acquiring one fully sampled non-diffusion-weighted image and 10 diffusion-weighted images at 8× undersampling would result in an 80% net reduction in data needed. We demonstrate efficacy of the DTV algorithm over TV in reducing data sampling requirements, which can be translated into higher apparent resolution and potentially shorter scan times. This method would be equally applicable in diffusion MRI applications outside the heart.

Regional variations in ex-vivo diffusion tensor anisotropy are associated with cardiomyocyte remodeling in rats after left ventricular pressure overload.

BACKGROUND: Pressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner. METHODS: We investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts. RESULTS: We observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p < 0.001, R2 = 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios. CONCLUSIONS: These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.

Quantification of model and data uncertainty in a network analysis of cardiac myocyte mechanosignalling.

Cardiac myocytes transduce changes in mechanical loading into cellular responses via interacting cell signalling pathways. We previously reported a logic-based ordinary differential equation model of the myocyte mechanosignalling network that correctly predicts 78% of independent experimental results not used to formulate the original model. Here, we use Monte Carlo and polynomial chaos expansion simulations to examine the effects of uncertainty in parameter values, model logic and experimental validation data on the assessed accuracy of that model. The prediction accuracy of the model was robust to parameter changes over a wide range being least sensitive to uncertainty in time constants and most affected by uncertainty in reaction weights. Quantifying epistemic uncertainty in the reaction logic of the model showed that while replacing 'OR' with 'AND' reactions greatly reduced model accuracy, replacing 'AND' with 'OR' reactions was more likely to maintain or even improve accuracy. Finally, data uncertainty had a modest effect on assessment of model accuracy. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Unmasking of oestrogen-dependent changes in left ventricular structure and function in aged female rats: a potential model for pre-heart failure with preserved ejection fraction.

KEY POINTS: Heart failure with preserved ejection fraction (HFpEF) is seen more frequently in older women; risk factors include age, hypertension and excess weight. No female animal models of early stage remodelling (pre-HFpEF) have examined the effects that the convergence of such factors have on cardiac structure and function. In this study, we demonstrate that ageing can lead to the development of mild chamber remodelling, diffuse fibrosis and loss of diastolic function. The loss of oestrogens further aggravates such changes by leading to a notable drop in cardiac output (while preserving normal ejection fraction) in the presence of diffuse fibrosis that is more predominant in endocardium and is accompanied by papillary fibrosis. Excess weight did not markedly aggravate such findings. This animal model recapitulates many of the features recognized in older, female HFpEF patients and thus, may serve to examine the effects of candidate therapeutic agents. ABSTRACT: Two-thirds of patients with heart failure with preserved ejection fraction (HFpEF) are older women, and risk factors include hypertension and excess weight/obesity. Pathophysiological factors that drive early disease development (before heart failure ensues) remain obscure and female animal models are lacking. The study evaluated the intersecting roles of ageing, oestrogen depletion and excess weight on altering cardiac structure/function. Female, 18-month-old, Fischer F344 rats were divided into an aged group, aged + ovariectomy (OVX) and aged + ovariectomy + 10% fructose (OVF) in drinking water (n = 8-16/group) to induce weight gain. Left ventricular (LV) structure/function was monitored by echocardiography. At 22 months of age, animals were anaesthetized and catheter-based haemodynamics evaluated, followed by histological measures of chamber morphometry and collagen density. All aged animals developed hypertension. OVF animals increased body weight. Echocardiography only detected mild chamber remodelling with ageing while intraventricular pressure-volume loop analysis showed significant (P < 0.05) decreases vs. ageing in stroke volume (13% OVX and 15% for OVF), stroke work (34% and 52%) and cardiac output (29% and 27%), and increases in relaxation time (10% OVX) with preserved ejection fraction. Histology indicated papillary and interstitial fibrosis with ageing, which was higher in the endocardium of OVX and OVF groups. With ageing, ovariectomy leads to the loss of diastolic and global LV function while preserving ejection fraction. This model recapitulates many cardiovascular features present in HFpEF patients and may help understand the roles that ageing and oestrogen depletion play in early (pre-HFpEF) disease development.

Predictive model identifies key network regulators of cardiomyocyte mechano-signaling.

Mechanical strain is a potent stimulus for growth and remodeling in cells. Although many pathways have been implicated in stretch-induced remodeling, the control structures by which signals from distinct mechano-sensors are integrated to modulate hypertrophy and gene expression in cardiomyocytes remain unclear. Here, we constructed and validated a predictive computational model of the cardiac mechano-signaling network in order to elucidate the mechanisms underlying signal integration. The model identifies calcium, actin, Ras, Raf1, PI3K, and JAK as key regulators of cardiac mechano-signaling and characterizes crosstalk logic imparting differential control of transcription by AT1R, integrins, and calcium channels. We find that while these regulators maintain mostly independent control over distinct groups of transcription factors, synergy between multiple pathways is necessary to activate all the transcription factors necessary for gene transcription and hypertrophy. We also identify a PKG-dependent mechanism by which valsartan/sacubitril, a combination drug recently approved for treating heart failure, inhibits stretch-induced hypertrophy, and predict further efficacious pairs of drug targets in the network through a network-wide combinatorial search.

Evaluation of non-Gaussian diffusion in cardiac MRI.

PURPOSE: The diffusion tensor model assumes Gaussian diffusion and is widely applied in cardiac diffusion MRI. However, diffusion in biological tissue deviates from a Gaussian profile as a result of hindrance and restriction from cell and tissue microstructure, and may be quantified better by non-Gaussian modeling. The aim of this study was to investigate non-Gaussian diffusion in healthy and hypertrophic hearts. METHODS: Thirteen rat hearts (five healthy, four sham, four hypertrophic) were imaged ex vivo. Diffusion-weighted images were acquired at b-values up to 10,000 s/mm2 . Models of diffusion were fit to the data and ranked based on the Akaike information criterion. RESULTS: The diffusion tensor was ranked best at b-values up to 2000 s/mm2 but reflected the signal poorly in the high b-value regime, in which the best model was a non-Gaussian "beta distribution" model. Although there was considerable overlap in apparent diffusivities between the healthy, sham, and hypertrophic hearts, diffusion kurtosis and skewness in the hypertrophic hearts were more than 20% higher in the sheetlet and sheetlet-normal directions. CONCLUSION: Non-Gaussian diffusion models have a higher sensitivity for the detection of hypertrophy compared with the Gaussian model. In particular, diffusion kurtosis may serve as a useful biomarker for characterization of disease and remodeling in the heart. Magn Reson Med 78:1174-1186, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Mechanotransduction in cardiac hypertrophy and failure.

Cardiac muscle cells have an intrinsic ability to sense and respond to mechanical load through a process known as mechanotransduction. In the heart, this process involves the conversion of mechanical stimuli into biochemical events that induce changes in myocardial structure and function. Mechanotransduction and its downstream effects function initially as adaptive responses that serve as compensatory mechanisms during adaptation to the initial load. However, under prolonged and abnormal loading conditions, the remodeling processes can become maladaptive, leading to altered physiological function and the development of pathological cardiac hypertrophy and heart failure. Although the mechanisms underlying mechanotransduction are far from being fully elucidated, human and mouse genetic studies have highlighted various cytoskeletal and sarcolemmal structures in cardiac myocytes as the likely candidates for load transducers, based on their link to signaling molecules and architectural components important in disease pathogenesis. In this review, we summarize recent developments that have uncovered specific protein complexes linked to mechanotransduction and mechanotransmission within the sarcomere, the intercalated disc, and at the sarcolemma. The protein structures acting as mechanotransducers are the first step in the process that drives physiological and pathological cardiac hypertrophy and remodeling, as well as the transition to heart failure, and may provide better insights into mechanisms driving mechanotransduction-based diseases.

Connexin defects underlie arrhythmogenic right ventricular cardiomyopathy in a novel mouse model.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) termed a 'disease of the desmosome' is an inherited cardiomyopathy that recently underwent reclassification owing to the identification of left-dominant and biventricular disease forms. Homozygous loss-of-function mutations in the desmosomal component, desmoplakin, are found in patients exhibiting a biventricular form of ARVC; however, no models recapitulate the postnatal hallmarks of the disease as seen in these patients. To gain insights into the homozygous loss-of-function effects of desmoplakin in the heart, we generated cardiomyocyte-specific desmoplakin-deficient mice (DSP-cKO) using ventricular myosin light chain-2-Cre mice. Homozygous DSP-cKO mice are viable but display early ultrastructural defects in desmosomal integrity leading to a cardiomyopathy reminiscent of a biventricular form of ARVC, which includes cell death and fibro-fatty replacement within the ventricle leading to biventricular dysfunction, failure and premature death. DSP-cKO mice also exhibited ventricular arrhythmias that are exacerbated with exercise and catecholamine stimulation. Furthermore, DSP-cKO hearts exhibited right ventricular conduction defects associated with loss of connexin 40 expression and electrical wavefront propagation defects associated with loss of connexin 43 expression. Dose-dependent assessment of the effects of loss of desmoplakin in neonatal ventricular cardiomyocytes revealed primary loss of connexin 43 levels, phosphorylation and function independent of the molecular dissociation of the mechanical junction complex and fibro-fatty manifestation associated with ARVC, suggesting a role for desmoplakin as a primary stabilizer of connexin integrity. In summary, we provide evidence for a novel mouse model, which is reminiscent of the postnatal onset of ARVC while highlighting mechanisms underlying a biventricular form of human ARVC.

Non-invasive, model-based measures of ventricular electrical dyssynchrony for predicting CRT outcomes.

AIMS: Left ventricular activation delay due to left bundle branch block (LBBB) is an important determinant of the severity of dyssynchronous heart failure (DHF). We investigated whether patient-specific computational models constructed from non-invasive measurements can provide measures of baseline dyssynchrony and its reduction after CRT that may explain the degree of long-term reverse ventricular remodelling. METHODS AND RESULTS: LV end-systolic volume reduction (ΔESVLV) measured by 2D trans-thoracic echocardiography in eight patients following 6 months of CRT was significantly (P < 0.05) greater in responders (26 ± 20%, n = 4) than non-responders (11 ± 16%, n = 4). LV reverse remodelling did not correlate with baseline QRS duration or its change after biventricular pacing, but did correlate with baseline LV endocardial activation measured by electroanatomic mapping (R2 = 0.71, P < 0.01). Patient-specific models of LBBB ventricular activation with parameters obtained by matching model-computed vectorcardiograms (VCG) to those derived from standard patient ECGs yielded LV endocardial activation times that correlated well with those measured from endocardial maps (R2 = 0.90). Model-computed 3D LV activation times correlated strongly with the reduction in LVESV (R2 = 0.93, P < 0.001). Computed decreases due to simulated CRT in the time delay between LV septal and lateral activation correlated strongly with ΔESVLV (R2 = 0.92, P < 0.001). Models also suggested that optimizing VV delays may improve resynchronization by this measure of activation delay. CONCLUSIONS: Patient-specific computational models constructed from non-invasive measurements can compute estimates of LV dyssynchrony and their changes after CRT that may be as good as or better than electroanatomic mapping for predicting long-term reverse remodelling.

Atlas-Based Ventricular Shape Analysis for Understanding Congenital Heart Disease.

Congenital heart disease is associated with abnormal ventricular shape that can affect wall mechanics and may be predictive of long-term adverse outcomes. Atlas-based parametric shape analysis was used to analyze ventricular geometries of eight adolescent or adult single-ventricle CHD patients with tricuspid atresia and Fontans. These patients were compared with an "atlas" of non-congenital asymptomatic volunteers, resulting in a set of z-scores which quantify deviations from the control population distribution on a patient-by-patient basis. We examined the potential of these scores to: (1) quantify abnormalities of ventricular geometry in single ventricle physiologies relative to the normal population; (2) comprehensively quantify wall motion in CHD patients; and (3) identify possible relationships between ventricular shape and wall motion that may reflect underlying functional defects or remodeling in CHD patients. CHD ventricular geometries at end-diastole and end-systole were individually compared with statistical shape properties of an asymptomatic population from the Cardiac Atlas Project. Shape analysis-derived model properties, and myocardial wall motions between end-diastole and end-systole, were compared with physician observations of clinical functional parameters. Relationships between altered shape and altered function were evaluated via correlations between atlas-based shape and wall motion scores. Atlas-based shape analysis identified a diverse set of specific quantifiable abnormalities in ventricular geometry or myocardial wall motion in all subjects. Moreover, this initial cohort displayed significant relationships between specific shape abnormalities such as increased ventricular sphericity and functional defects in myocardial deformation, such as decreased long-axis wall motion. These findings suggest that atlas-based ventricular shape analysis may be a useful new tool in the management of patients with CHD who are at risk of impaired ventricular wall mechanics and chamber remodeling.

Modifying Ventricular Fibrillation by Targeted Rotor Substrate Ablation: Proof-of-Concept from Experimental Studies to Clinical VF.

INTRODUCTION: Recent work has suggested a role for organized sources in sustaining ventricular fibrillation (VF). We assessed whether ablation of rotor substrate could modulate VF inducibility in canines, and used this proof-of-concept as a foundation to suppress antiarrhythmic drug-refractory clinical VF in a patient with structural heart disease. METHODS AND RESULTS: In 9 dogs, we introduced 64-electrode basket catheters into one or both ventricles, used rapid pacing at a recorded induction threshold to initiate VF, and then defibrillated after 18±8 seconds. Endocardial rotor sites were identified from basket recordings using phase mapping, and ablation was performed at nonrotor (sham) locations (7 ± 2 minutes) and then at rotor sites (8 ± 2 minutes, P = 0.10 vs. sham); the induction threshold was remeasured after each. Sham ablation did not alter canine VF induction threshold (preablation 150 ± 16 milliseconds, postablation 144 ± 16 milliseconds, P = 0.54). However, rotor site ablation rendered VF noninducible in 6/9 animals (P = 0.041), and increased VF induction threshold in the remaining 3. Clinical proof-of-concept was performed in a patient with repetitive ICD shocks due to VF refractory to antiarrhythmic drugs. Following biventricular basket insertion, VF was induced and then defibrillated. Mapping identified 4 rotors localized at borderzone tissue, and rotor site ablation (6.3 ± 1.5 minutes/site) rendered VF noninducible. The VF burden fell from 7 ICD shocks in 8 months preablation to zero ICD therapies at 1 year, without antiarrhythmic medications. CONCLUSIONS: Targeted rotor substrate ablation suppressed VF in an experimental model and a patient with refractory VF. Further studies are warranted on the efficacy of VF source modulation.

Timing and magnitude of systolic stretch affect myofilament activation and mechanical work.

Dyssynchronous activation of the heart leads to abnormal regional systolic stretch. In vivo studies have suggested that the timing of systolic stretch can affect regional tension and external work development. In the present study, we measured the direct effects of systolic stretch timing on the magnitude of tension and external work development in isolated murine right ventricular papillary muscles. A servomotor was used to impose precisely timed stretches relative to electrical activation while a force transducer measured force output and strain was monitored using a charge-couple device camera and topical markers. Stretches taking place during peak intracellular Ca(2+) statistically increased peak tension up to 270%, whereas external work due to stretches in this interval reached values of 500 J/m. An experimental analysis showed that time-varying elastance overestimated peak tension by 100% for stretches occurring after peak isometric tension. The addition of the force-velocity relation explained some effects of stretches occurring before the peak of the Ca(2+) transient but had no effect in later stretches. An estimate of transient deactivation was measured by performing quick stretches to dissociate cross-bridges. The timing of transient deactivation explained the remaining differences between the model and experiment. These results suggest that stretch near the start of cardiac tension development substantially increases twitch tension and mechanical work production, whereas late stretches decrease external work. While the increased work can mostly be explained by the time-varying elastance of cardiac muscle, the decreased work in muscles stretched after the peak of the Ca(2+) transient is largely due to myofilament deactivation.

Biomechanics of cardiac electromechanical coupling and mechanoelectric feedback.

Cardiac mechanical contraction is triggered by electrical activation via an intracellular calcium-dependent process known as excitation-contraction coupling. Dysregulation of cardiac myocyte intracellular calcium handling is a common feature of heart failure. At the organ scale, electrical dyssynchrony leads to mechanical alterations and exacerbates pump dysfunction in heart failure. A reverse coupling between cardiac mechanics and electrophysiology is also well established. It is commonly referred as cardiac mechanoelectric feedback and thought to be an important contributor to the increased risk of arrhythmia during pathological conditions that alter regional cardiac wall mechanics, including heart failure. At the cellular scale, most investigations of myocyte mechanoelectric feedback have focused on the roles of stretch-activated ion channels, though mechanisms that are independent of ionic currents have also been described. Here we review excitation-contraction coupling and mechanoelectric feedback at the cellular and organ scales, and we identify the need for new multicellular tissue-scale model systems and experiments that can help us to obtain a better understanding of how interactions between electrophysiological and mechanical processes at the cell scale affect ventricular electromechanical interactions at the organ scale in the normal and diseased heart.

Successful Cardiac Resynchronization Therapy Reduces Negative Septal Work in Patient-Specific Models of Dyssynchronous Heart Failure.

In patients with dyssynchronous heart failure (DHF), cardiac conduction abnormalities cause the regional distribution of myocardial work to be non-homogeneous. Cardiac resynchronization therapy (CRT) using an implantable, programmed biventricular pacemaker/defibrillator, can improve the synchrony of contraction between the right and left ventricles in DHF, resulting in reduced morbidity and mortality and increased quality of life. Since regional work depends on wall stress, which cannot be measured in patients, we used computational methods to investigate regional work distributions and their changes after CRT. We used three-dimensional multi-scale patient-specific computational models parameterized by anatomic, functional, hemodynamic, and electrophysiological measurements in eight patients with heart failure and left bundle branch block (LBBB) who received CRT. To increase clinical translatability, we also explored whether streamlined computational methods provide accurate estimates of regional myocardial work. We found that CRT increased global myocardial work efficiency with significant improvements in non-responders. Reverse ventricular remodeling after CRT was greatest in patients with the highest heterogeneity of regional work at baseline, however the efficacy of CRT was not related to the decrease in overall work heterogeneity or to the reduction in late-activated regions of high myocardial work. Rather, decreases in early-activated regions of myocardium performing negative myocardial work following CRT best explained patient variations in reverse remodeling. These findings were also observed when regional myocardial work was estimated using ventricular pressure as a surrogate for myocardial stress and changes in endocardial surface area as a surrogate for strain. These new findings suggest that CRT promotes reverse ventricular remodeling in human dyssynchronous heart failure by increasing regional myocardial work in early-activated regions of the ventricles, where dyssynchrony is specifically associated with hypoperfusion, late systolic stretch, and altered metabolic activity and that measurement of these changes can be performed using streamlined approaches.