Rebecca A.B. Burton

Generation of histo−anatomically representative models of the individual heart: tools and application

This paper presents methods to build histo-anatomically detailed individualized cardiac models. The models are based on high-resolution three-dimensional anatomical and/or diffusion tensor magnetic resonance images, combined with serial histological sectioning data, and are used to investigate individualized cardiac function. The current state of the art is reviewed, and its limitations are discussed. We assess the challenges associated with the generation of histo-anatomically representative individualized in silico models of the heart. The entire processing pipeline including image acquisition, image processing, mesh generation, model set-up and execution of computer simulations, and the underlying methods are described. The multifaceted challenges associated with these goals are highlighted, suitable solutions are proposed, and an important application of developed high-resolution structure–function models in elucidating the effect of individual structural heterogeneity upon wavefront dynamics is demonstrated.

Development of an anatomically detailed MRI-derived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function

Recent advances in magnetic resonance (MR) imaging technology have unveiled a wealth of information regarding cardiac histoanatomical complexity. However, methods to faithfully translate this level of fine-scale structural detail into computational whole ventricular models are still in their infancy, and, thus, the relevance of this additional complexity for simulations of cardiac function has yet to be elucidated. Here, we describe the development of a highly detailed finite-element computational model (resolution: ∼125 μm) of rabbit ventricles constructed from high-resolution MR data (raw data resolution: 43 × 43 × 36 μm), including the processes of segmentation (using a combination of level-set approaches), identification of relevant anatomical features, mesh generation, and myocyte orientation representation (using a rule-based approach). Full access is provided to the completed model and MR data. Simulation results were compared with those from a simplified model built from the same images but excluding finer anatomical features (vessels/endocardial structures). Initial simulations showed that the presence of trabeculations can provide shortcut paths for excitation, causing regional differences in activation after pacing between models. Endocardial structures gave rise to small-scale virtual electrodes upon the application of external field stimulation, which appeared to protect parts of the endocardium in the complex model from strong polarizations, whereas intramural virtual electrodes caused by blood vessels and extracellular cleft spaces appeared to reduce polarization of the epicardium. Postshock, these differences resulted in the genesis of new excitation wavefronts that were not observed in more simplified models. Furthermore, global differences in the stimulus recovery rates of apex/base regions were observed, causing differences in the ensuing arrhythmogenic episodes. In conclusion, structurally simplified models are well suited for a large range of cardiac modeling applications. However, important differences are seen when behavior at microscales is relevant, particularly when examining the effects of external electrical stimulation on tissue electrophysiology and arrhythmia induction. This highlights the utility of histoanatomically detailed models for investigations of cardiac function, in particular for future patient-specific modeling.

Axial Stretch of Rat Single Ventricular Cardiomyocytes Causes an Acute and Transient Increase in Ca2+ Spark Rate

We investigate acute effects of axial stretch, applied by carbon fibers (CFs), on diastolic Ca2+ spark rate in rat isolated cardiomyocytes. CFs were attached either to both cell ends (to maximize the stretched region), or to the center and one end of the cell (to compare responses in stretched and nonstretched half-cells). Sarcomere length was increased by 8.01±0.94% in the stretched cell fraction, and time series of XY confocal images were recorded to monitor diastolic Ca2+ spark frequency and dynamics. Whole-cell stretch causes an acute increase of Ca2+ spark rate (to 130.7±6.4%) within 5 seconds, followed by a return to near background levels (to 104.4±5.1%) within 1 minute of sustained distension. Spark rate increased only in the stretched cell region, without significant differences in spark amplitude, time to peak, and decay time constants of sparks in stretched and nonstretched areas. Block of stretch-activated ion channels (2 &mgr;mol/L GsMTx-4), perfusion with Na+/Ca2+-free solution, and block of nitric oxide synthesis (1 mmol/L L-NAME) all had no effect on the stretch-induced acute increase in Ca2+ spark rate. Conversely, interference with cytoskeletal integrity (2 hours of 10 &mgr;mol/L colchicine) abolished the response. Subsequent electron microscopic tomography confirmed the close approximation of microtubules with the T-tubular–sarcoplasmic reticulum complex (to within ≈10-8m). In conclusion, axial stretch of rat cardiomyocytes acutely and transiently increases sarcoplasmic reticulum Ca2+ spark rate via a mechanism that is independent of sarcolemmal stretch-activated ion channels, nitric oxide synthesis, or availability of extracellular calcium but that requires cytoskeletal integrity. The potential of microtubule-mediated modulation of ryanodine receptor function warrants further investigation.

Three-Dimensional Models of Individual Cardiac Histoanatomy: Tools and Challenges

Abstract:  There is a need for, and utility in, the acquisition of data sets of cardiac histoanatomy, with the vision of reconstructing individual hearts on the basis of noninvasive imaging, such as MRI, enriched by reference to detailed atlases of serial histology obtained from representative samples. These data sets would be useful not only as a repository of knowledge regarding the specifics of cardiac histoanatomy, but could form the basis for generation of individualized high‐resolution cardiac structure–function models. The current article presents a step in this general direction: it illustrates how whole‐heart noninvasive imaging can be combined with whole‐heart histology in an approach to achieve automated construction of histoanatomically detailed models of cardiac 3D structure and function at hitherto unprecedented resolution and accuracy (based on 26.4 × 26.4 × 24.4 μm MRI voxel size, and enriched by histological detail). It provides an overview of the tools used in this quest and outlines challenges posed by the approach in the light of applications that may benefit from the availability of such data and tools.

Histo-anatomical structure of the living isolated rat heart in two contraction states assessed by diffusion tensor MRI

Deformation and wall-thickening of ventricular myocardium are essential for cardiac pump function. However, insight into the histo-anatomical basis for cardiac tissue re-arrangement during contraction is limited. In this report, we describe dynamic changes in regionally prevailing cardiomyocyte (fibre) and myolaminar (sheet) orientations, using Diffusion Tensor Imaging (DTI) of ventricles in the same living heart in two different mechanical states. Hearts, isolated from Sprague–Dawley rats, were Langendorff-perfused and imaged, initially in their slack state during cardioplegic arrest, then during lithium-induced contracture. Regional fibre- and sheet-orientations were derived from DTI-data on a voxel-wise basis. Contraction was accompanied with a decrease in left-handed helical fibres (handedness relative to the baso-apical direction) in basal, equatorial, and apical sub-epicardium (by 14.0%, 17.3%, 15.8% respectively; p < 0.001), and an increase in right-handed helical fibres of the sub-endocardium (by 11.0%, 12.1% and 16.1%, respectively; p < 0.001). Two predominant sheet-populations were observed, with sheet-angles of either positive (β+) or negative (β−) polarity relative to a ‘chamber-horizontal plane’ (defined as normal to the left ventricular long-axis). In contracture, mean ‘intersection’-angle (geometrically quantifiable intersection of sheet-angle projections) between β+ and β− sheet-populations increased from 86.2 ± 5.5° (slack) to 108.3 ± 5.4° (p < 0.001). Subsequent high-resolution DTI of fixed myocardium, and histological sectioning, reconfirmed the existence of alternating sheet-plane populations. Our results suggest that myocardial tissue layers in alternating sheet-populations align into a more chamber-horizontal orientation during contraction. This re-arrangement occurs via an accordion-like mechanism that, combined with inter-sheet slippage, can significantly contribute to ventricular deformation, including wall-thickening in a predominantly centripetal direction and baso-apical shortening.

Rearrangement of Atrial Bundle Architecture and Consequent Changes in Anisotropy of Conduction Constitute the 3-Dimensional Substrate for Atrial Fibrillation

Background— Anisotropy of conduction facilitates re-entry and is, therefore, a key determinant of the stability of atrial fibrillation (AF). Little is known about the effect of AF on atrial bundle architecture and consequent changes in anisotropy of conduction and maintenance of AF. Methods and Results— Direct contact mapping was performed in left atria of goats with acute AF (n=6) or persistent AF (n=5). The degree and direction of anisotropic conduction were analyzed. Mapped tissue regions were imaged by high-resolution MRI for identification of endocardial and epicardial bundle directions. Correlation between endocardial and epicardial bundle directions and between bundle directions and anisotropic conduction was quantified. In persistent AF, epicardial bundles were oriented more perpendicularly to endocardial bundles than in acute AF (% angles <20° between epicardial and endocardial bundle directions were 7.63% and 21.25%, respectively; P<0.01). In acute AF, the direction of epicardially mapped anisotropic conduction correlated with endocardial but not with epicardial bundles. In persistent AF, the direction of anisotropic conduction correlated better with epicardial than with endocardial bundles (% angles <20° between direction of anisotropic conduction and bundle direction were 28.77% and 18.45%, respectively; P<0.01). Conclusions— During AF, atrial bundle rearrangement manifests itself in more perpendicular orientation of epicardial to endocardial bundles. Propagation of fibrillation waves is dominated by endocardial bundles in acute AF and by epicardial bundles in persistent AF. Together with the loss of endo-epicardial electrical connections, rearrangement of atrial bundles underlies endo-epicardial dissociation of electrical activity and the development of a 3-dimensional AF substrate.

Progressive changes in T₁, T₂ and left-ventricular histo-architecture in the fixed and embedded rat heart.

Chemical tissue fixation, followed by embedding in either agarose or Fomblin, is common practice in time‐intensive MRI studies of ex vivo biological samples, and is required to prevent tissue autolysis and sample motion. However, the combined effect of fixation and sample embedding may alter tissue structure and MRI properties. We investigated the progressive changes in T1 and T2 relaxation times, and the arrangement of locally prevailing cardiomyocyte orientation determined using diffusion tensor imaging, in embedded ex vivo rat hearts fixed using Karnovskys solution (glutaraldehyde–formaldehyde mix). Three embedding media were investigated: (i) standard agarose (n = 3 hearts); (ii) Fomblin (n = 4 hearts); and (iii) iso‐osmotic agarose (n = 3 hearts); in the latter, the osmolarity of the fixative and embedding medium was adjusted to 300 mOsm to match more closely that of native tissue. The T1 relaxation time in the myocardium showed a pronounced decrease over a 48‐h period following embedding in Fomblin (−11.3 ± 6.2%; mean ± standard deviation), but was stable in standard agarose‐ and iso‐osmotic agarose‐embedded hearts. The mean myocardial T2 relaxation time increased in all embedded hearts: by 35.1 ± 14.7% with standard agarose embedding, 13.1 ± 5.6% with Fomblin and 13.3 ± 1.4% with iso‐osmotic agarose. Deviation in the orientation of the primary eigenvector of the diffusion tensor occurred in all hearts (mean angular changes of 6.6°, 3.2° and 1.9° per voxel after 48 h in agarose‐, Fomblin‐ and iso‐osmotic agarose‐embedded hearts, respectively), indicative of progressive structural changes in myocardial histo‐architecture, in spite of previous exposure to fast‐acting tissue fixation. Our results suggest that progressive structural changes occur in chemically fixed myocardium, and that the extent of these changes is modulated by the embedding medium, and by osmotic gradients between the fixative in the tissue and the surrounding medium. Copyright

Fast Measurement of Sarcomere Length and Cell Orientation in Langendorff-Perfused Hearts Using Remote Focusing Microscopy

Rationale: Sarcomere length (SL) is a key indicator of cardiac mechanical function, but current imaging technologies are limited in their ability to unambiguously measure and characterize SL at the cell level in intact, living tissue. Objective: We developed a method for measuring SL and regional cell orientation using remote focusing microscopy, an emerging imaging modality that can capture light from arbitrary oblique planes within a sample. Methods and Results: We present a protocol that unambiguously and quickly determines cell orientation from user-selected areas in a field of view by imaging 2 oblique planes that share a common major axis with the cell. We demonstrate the effectiveness of the technique in establishing single-cell SL in Langendorff-perfused hearts loaded with the membrane dye di-4-ANEPPS. Conclusions: Remote focusing microscopy can measure cell orientation in complex 2-photon data sets without capturing full z stacks. The technique allows rapid assessment of SL in healthy and diseased heart experimental preparations.

Human-based approaches to pharmacology and cardiology: an interdisciplinary and intersectorial workshop

Both biomedical research and clinical practice rely on complex datasets for the physiological and genetic characterization of human hearts in health and disease. Given the complexity and variety of approaches and recordings, there is now growing recognition of the need to embed computational methods in cardiovascular medicine and science for analysis, integration and prediction. This paper describes a Workshop on Computational Cardiovascular Science that created an international, interdisciplinary and inter-sectorial forum to define the next steps for a human-based approach to disease supported by computational methodologies. The main ideas highlighted were (i) a shift towards human-based methodologies, spurred by advances in new in silico, in vivo, in vitro, and ex vivo techniques and the increasing acknowledgement of the limitations of animal models. (ii) Computational approaches complement, expand, bridge, and integrate in vitro, in vivo, and ex vivo experimental and clinical data and methods, and as such they are an integral part of human-based methodologies in pharmacology and medicine. (iii) The effective implementation of multi- and interdisciplinary approaches, teams, and training combining and integrating computational methods with experimental and clinical approaches across academia, industry, and healthcare settings is a priority. (iv) The human-based cross-disciplinary approach requires experts in specific methodologies and domains, who also have the capacity to communicate and collaborate across disciplines and cross-sector environments. (v) This new translational domain for human-based cardiology and pharmacology requires new partnerships supported financially and institutionally across sectors. Institutional, organizational, and social barriers must be identified, understood and overcome in each specific setting.

Three-dimensional histology: tools and application to quantitative assessment of cell-type distribution in rabbit heart.

Aims Cardiac histo-anatomical organization is a major determinant of function. Changes in tissue structure are a relevant factor in normal and disease development, and form targets of therapeutic interventions. The purpose of this study was to test tools aimed to allow quantitative assessment of cell-type distribution from large histology and magnetic resonance imaging- (MRI) based datasets. Methods and results Rabbit heart fixation during cardioplegic arrest and MRI were followed by serial sectioning of the whole heart and light-microscopic imaging of trichrome-stained tissue. Segmentation techniques developed specifically for this project were applied to segment myocardial tissue in the MRI and histology datasets. In addition, histology slices were segmented into myocytes, connective tissue, and undefined. A bounding surface, containing the whole heart, was established for both MRI and histology. Volumes contained in the bounding surface (called ‘anatomical volume’), as well as that identified as containing any of the above tissue categories (called ‘morphological volume’), were calculated. The anatomical volume was 7.8 cm3 in MRI, and this reduced to 4.9 cm3 after histological processing, representing an ‘anatomical’ shrinkage by 37.2%. The morphological volume decreased by 48% between MRI and histology, highlighting the presence of additional tissue-level shrinkage (e.g. an increase in interstitial cleft space). The ratio of pixels classified as containing myocytes to pixels identified as non-myocytes was roughly 6:1 (61.6 vs. 9.8%; the remaining fraction of 28.6% was ‘undefined’). Conclusion Qualitative and quantitative differentiation between myocytes and connective tissue, using state-of-the-art high-resolution serial histology techniques, allows identification of cell-type distribution in whole-heart datasets. Comparison with MRI illustrates a pronounced reduction in anatomical and morphological volumes during histology processing.