Ian J. LeGrice

Fibroblast Network in Rabbit Sinoatrial Node. Structural and Functional Identification of Homogeneous and Heterogeneous Cell Coupling

Abstract— Cardiomyocytes form a conducting network that is assumed to be electrically isolated from nonmyocytes in vivo. In cell culture, however, cardiac fibroblasts can contribute to the spread of excitation via functional gap junctions with cardiomyocytes. To assess the ability of fibroblasts to form gap junctions in vivo, we combine in situ detection of connexins in rabbit sinoatrial node (a tissue that is particularly rich in fibroblasts) with identification of myocytes and fibroblasts using immunohistochemical labeling and confocal microscopy. We distinguish two spatially distinct fibroblast populations expressing different connexins: fibroblasts surrounded by other fibroblasts preferentially express connexin40, whereas fibroblasts that are intermingled with myocytes largely express connexin45. Functionality of homogeneous and heterogeneous cell coupling was investigated by dye transfer in sinoatrial node tissue explants. These studies reveal spread of Lucifer yellow, predominantly along extended threads of interconnected fibroblasts (probably via connexin40), and occasionally between neighboring fibroblasts and myocytes (probably via connexin45). Our findings show that cardiac fibroblasts form a coupled network of cells, which may be functionally linked to myocytes in rabbit SAN.

Cardiac Microstructure: Implications for Electrical Propagation and Defibrillation in the Heart

Abstract— Our understanding of the electrophysiological properties of the heart is incomplete. We have investigated two issues that are fundamental to advancing that understanding. First, there has been widespread debate over the mechanisms by which an externally applied shock can influence a sufficient volume of heart tissue to terminate cardiac fibrillation. Second, it has been uncertain whether cardiac tissue should be viewed as an electrically orthotropic structure, or whether its electrical properties are, in fact, isotropic in the plane orthogonal to myofiber direction. In the present study, a computer model that incorporates a detailed three-dimensional representation of cardiac muscular architecture is used to investigate these issues. We describe a bidomain model of electrical propagation solved in a discontinuous domain that accurately represents the microstructure of a transmural block of rat left ventricle. From analysis of the model results, we conclude that (1) the laminar organization of myocytes determines unique electrical properties in three microstructurally defined directions at any point in the ventricular wall of the heart, and (2) interlaminar clefts between layers of cardiomyocytes provide a substrate for bulk activation of the ventricles during defibrillation.

Three-dimensional transmural organization of perimysial collagen in the heart

There is strong support for the view that the ventricular myocardium has a laminar organization in which myocytes are grouped into branching layers separated by cleavage planes. However, understanding of the extent and functional implications of this architecture has been limited by the lack of a systematic three-dimensional description of the organization of myocytes and associated perimysial collagen. We imaged myocytes and collagen across the left ventricular wall at high resolution in seven normal rat hearts using extended volume confocal microscopy. We developed novel reconstruction and segmentation techniques necessary for the quantitative analysis of three-dimensional myocyte and perimysial collagen organization. The results confirm that perimysial collagen has an ordered arrangement and that it defines a laminar organization. Perimysial collagen is composed of three distinct forms: extensive meshwork on laminar surfaces, convoluted fibers connecting adjacent layers, and longitudinal cords. While myolaminae are the principal form of structural organization throughout most of the wall, they are not seen in the subepicardium, where perimysial collagen is present only as longitudinal cords.

Laminar arrangement of ventricular myocytes influences electrical behavior of the heart.

The response of the heart to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. A long-held view of cardiac electrical anisotropy is that electrical conductivity is greatest along the myocyte axis allowing most rapid propagation of electrical activation in this direction, and that conductivity is isotropic transverse to the myocyte axis supporting a slower uniform spread of activation in this plane. In this context, knowledge of conductivity in two directions, parallel and transverse to the myofiber axis, is sufficient to characterize the electrical action of the heart. Here we present new experimental data that challenge this view. We have used a novel combination of intramural electrical mapping, and experiment-specific computer modeling, to demonstrate that left ventricular myocardium has unique bulk conductivities associated with three microstructurally-defined axes. We show that voltage fields induced by intramural current injection are influenced by not only myofiber direction, but also the transmural arrangement of muscle layers or myolaminae. Computer models of these experiments, in which measured 3D tissue structure was reconstructed in-silico, best matched recorded voltages with conductivities in the myofiber direction, and parallel and normal to myolaminae, set in the ratio 4:2:1, respectively. These findings redefine cardiac tissue as an electrically orthotropic substrate and enhance our understanding of how external shocks may act to successfully reset the fibrillating heart into a uniform electrical state. More generally, the mechanisms governing the destabilization of coordinated electrical propagation into ventricular arrhythmia need to be evaluated in the light of this discovery.

Three Distinct Directions of Intramural Activation Reveal Nonuniform Side-to-Side Electrical Coupling of Ventricular Myocytes

Background—The anisotropy of cardiac tissue is a key determinant of 3D electric propagation and the stability of activation wave fronts in the heart. The electric properties of ventricular myocardium are widely assumed to be axially anisotropic, with activation propagating most rapidly in the myofiber direction and at uniform velocity transverse to this. We present new experimental evidence that contradicts this view. Methods and Results—For the first time, high-density intramural electric mapping (325 electrodes at ≈4×4×1-mm spacing) from pig left ventricular tissue was used to reconstruct 3D paced activation surfaces projected directly onto 3D tissue structure imaged throughout the same left ventricular volume. These data from 5 hearts demonstrate that ventricular tissue is electrically orthotropic with 3 distinct propagation directions that coincide with local microstructural axes defined by the laminar arrangement of ventricular myocytes. The maximum conduction velocity of 0.67±0.019 ms−1 was aligned with the myofiber axis. However, transverse to this, the maximum conduction velocity was 0.30±0.010 ms−1, parallel to the myocyte layers and 0.17±0.004 ms−1 normal to them. These orthotropic conduction velocities give rise to preferential activation pathways across the left ventricular free wall that are not captured by structurally detailed computer models, which incorporate axially anisotropic electric properties. Conclusions—Our findings suggest that current views on uniform side-to-side electric coupling in the heart need to be revised. In particular, nonuniform laminar myocardial architecture and associated electric orthotropy should be included in future models of initiation and maintenance of ventricular arrhythmia.

A triaxial-measurement shear-test device for soft biological tissues.

A novel shear-test device for soft biological tissue, capable of applying simple shear deformations simultaneously in two orthogonal directions while measuring the resulting forces generated in three axes, is described. We validated the device using a synthetic gel, the properties of which were ascertained from independent tensile and rotational shear tests. Material parameters for the gel were fitted using neo-Hookean analytical solutions to the independent test data, and these matched the results from the device. Preliminary results obtained with rat septal myocardium are also presented to demonstrate the feasibility of the apparatus in determining the shear characteristics of living tissue.

Intramural Multisite Recording of Transmembrane Potential in the Heart

Heart surface optical mapping of transmembrane potentials has been widely used in studies of normal and pathological heart rhythms and defibrillation. In these studies, three-dimensional spatio-temporal events can only be inferred from two-dimensional surface potential maps. We present a novel optical system that enables high fidelity transmural recording of transmembrane potentials. A probe constructed from optical fibers is used to deliver excitation light and collect fluorescence from seven positions, each 1 mm apart, through the left ventricle wall of the rabbit heart. Excitation is provided by the 488-nm line of a water-cooled argon-ion laser. The fluorescence of the voltage-sensitive dye di-4-ANEPPS from each tissue site is split at 600 nm and imaged onto separate photodiodes for later signal ratioing. The optics and electronics are easily expandable to accommodate multiple optical probes. The system is used to record the first simultaneous measurements of transmembrane potential at a number of sites through the intact heart wall.

Method and apparatus for soft tissue material parameter estimation using tissue tagged Magnetic Resonance Imaging

We describe an experimental method and apparatus for the estimation of constitutive parameters of soft tissue using Magnetic Resonance Imaging (MRI), in particular for the estimation of passive myocardial material properties. MRI tissue tagged images were acquired with simultaneous pressure recordings, while the tissue was cyclically deformed using a custom built reciprocating pump actuator A continuous three-dimensional (3D) displacement field was reconstructed from the imaged tag motion. Cavity volume changes and local tissue microstructure were determined from phase contrast velocity and diffusion tensor MR images, respectively. The Finite Element Method (FEM) was used to solve the finite elasticity problem and obtain the displacement field that satisfied the applied boundary conditions and a given set of material parameters. The material parameters which best fit the FEM predicted displacements to the displacements reconstructed from the tagged images were found by nonlinear optimization. The equipment and method were validated using inflation of a deformable silicon gel phantom in the shape of a cylindrical annulus. The silicon gel was well described by a neo-Hookian material law with a single material parameter C1=8.71+/-0.06kPa, estimated independently using a rotational shear apparatus. The MRI derived parameter was allowed to vary regionally and was estimated as C1 =8.80+/-0.86kPa across the model. Preliminary results from the passive inflation of an isolated arrested pig heart are also presented, demonstrating the feasibility of the apparatus and method for isolated heart preparations. FEM based models can therefore estimate constitutive parameters accurately and reliably from MRI tagging data.

An image-based model of atrial muscular architecture: Effects of structural anisotropy on electrical activation

Background— Computer models that capture key features of the heterogeneous myofiber architecture of right and left atria and interatrial septum provide a means of investigating the mechanisms responsible for atrial arrhythmia. The data necessary to implement such models have not previously been available. The aims of this study were to characterize surface geometry and myofiber architecture throughout the atrial chambers and to investigate the effects of this structure on atrial activation. Methods and Results— Atrial surface geometry and myofiber orientations were reconstructed in 3D at 50×50×50-&mgr;m3 resolution from serial images acquired throughout the sheep atrial chambers. Myofiber orientations were determined by Eigen-analysis of the structure tensor. These data have been incorporated into an anatomic model that provides the first quantitative representation of myofiber architecture throughout the atrial chambers. By simulating activation on this 3D structure, we have confirmed the roles of specialized myofiber tracts such as the crista terminalis, pectinate muscles, and the Bachman bundle on the spread of activation from the sinus node. We also demonstrate how the complex myocyte arrangement in the posterior left atrium contributes to activation time dispersion adjacent to the pulmonary veins and increased vulnerability to rhythm disturbance generated by ectopic stimuli originating in the pulmonary vein sleeves. Conclusions— We have developed a structurally detailed, image-based model of atrial anatomy that provides deeper understanding of the role that myocyte architecture plays in normal and abnormal atrial electric function.Background— Computer models that capture key features of the heterogeneous myofiber architecture of right and left atria and interatrial septum provide a means of investigating the mechanisms responsible for atrial arrhythmia. The data necessary to implement such models have not previously been available. The aims of this study were to characterize surface geometry and myofiber architecture throughout the atrial chambers and to investigate the effects of this structure on atrial activation. Methods and Results— Atrial surface geometry and myofiber orientations were reconstructed in 3D at 50×50×50-μm3 resolution from serial images acquired throughout the sheep atrial chambers. Myofiber orientations were determined by Eigen-analysis of the structure tensor. These data have been incorporated into an anatomic model that provides the first quantitative representation of myofiber architecture throughout the atrial chambers. By simulating activation on this 3D structure, we have confirmed the roles of specialized myofiber tracts such as the crista terminalis, pectinate muscles, and the Bachman bundle on the spread of activation from the sinus node. We also demonstrate how the complex myocyte arrangement in the posterior left atrium contributes to activation time dispersion adjacent to the pulmonary veins and increased vulnerability to rhythm disturbance generated by ectopic stimuli originating in the pulmonary vein sleeves. Conclusions— We have developed a structurally detailed, image-based model of atrial anatomy that provides deeper understanding of the role that myocyte architecture plays in normal and abnormal atrial electric function.

Cardiac electrophysiology and tissue structure: bridging the scale gap with a joint measurement and modelling paradigm

Significant tissue structures exist in cardiac ventricular tissue that are of supracellular dimension. It is hypothesized that these tissue structures contribute to the discontinuous spread of electrical activation, may contribute to arrhymogenesis and also provide a substrate for effective cardioversion. However, the influences of these mesoscale tissue structures in intact ventricular tissue are difficult to understand solely on the basis of experimental measurement. Current measurement technology is able to record at both the macroscale tissue level and the microscale cellular or subcellular level, but to date it has not been possible to obtain large volume, direct measurements at the mesoscales. To bridge this scale gap in experimental measurements, we use tissue‐specific structure and mathematical modelling. Our models have enabled us to consider key hypotheses regarding discontinuous activation. We also consider the future developments of our intact tissue experimental programme.