Nondestructive optical determination of fiber organization in intact myocardial wall

Oxidation-reduction imaging of myoglobin reveals two-phase oxidation in the reperfused myocardium

Myocardial infarction (MI) is a serious acute cardiovascular syndrome that causes myocardial injury due to blood flow obstruction to a specific myocardial area. Under ischemic-reperfusion settings, a burst of reactive oxygen species is generated, leading to redox imbalance that could be attributed to several molecules, including myoglobin. Myoglobin is dynamic and exhibits various oxidation-reduction states that have been an early subject of attention in the food industry, specifically for meat consumers. However, rarely if ever have the myoglobin optical properties been used to measure the severity of MI. In the current study, we develop a novel imaging pipeline that integrates tissue clearing, confocal and light sheet fluorescence microscopy, combined with imaging analysis, and processing tools to investigate and characterize the oxidation-reduction states of myoglobin in the ischemic area of the cleared myocardium post-MI. Using spectral imaging, we have characterized the endogenous fluorescence of the myocardium and demonstrated that it is partly composed by fluorescence of myoglobin. Under ischemia-reperfusion experimental settings, we report that the infarcted myocardium spectral signature is similar to that of oxidized myoglobin signal that peaks 3 h post-reperfusion and decreases with cardioprotection. The infarct size assessed by oxidation-reduction imaging at 3 h post-reperfusion was correlated to the one estimated with late gadolinium enhancement MRI at 24 h post-reperfusion. In conclusion, this original work suggests that the redox state of myoglobin can be used as a promising imaging biomarker for characterizing and estimating the size of the MI during early phases of reperfusion.

Quantifying the microstructural and biomechanical changes in the porcine ventricles during growth and remodelling

Cardiac tissue growth and remodelling (G & R) occur in response to the changing physiological demands of the heart after birth. The early shift to pulmonary circulation produces an immediate increase in ventricular workload, causing microstructural and biomechanical changes that serve to maintain overall physiological homoeostasis. Such cardiac G & R continues throughout life. Quantifying the tissue's mechanical and microstructural changes because of G & R is of increasing interest, dovetailing with the emerging fields of personalised and precision solutions. This study aimed to determine equibiaxial, and non-equibiaxial extension, stress-relaxation, and the underlying microstructure of the passive porcine ventricles tissue at four time points spanning from neonatal to adulthood. The three-dimensional microstructure was investigated via two-photon excited fluorescence and second-harmonic generation microscopy on optically cleared tissues, describing the 3D orientation, rotation and dispersion of the cardiomyocytes and collagen fibrils. The results revealed that during biomechanical testing, myocardial ventricular tissue possessed non-linear, anisotropic, and viscoelastic behaviour. An increase in stiffness and viscoelasticity was noted for the left and right ventricular free walls from neonatal to adulthood. Microstructural analyses revealed concomitant increases in cardiomyocyte rotation and dispersion. This study provides baseline data, describing the biomechanical and microstructural changes in the left and right ventricular myocardial tissue during G & R, which should prove valuable to researchers in developing age-specific, constitutive models for more accurate computational simulations. STATEMENT OF SIGNIFICANCE: There is a dearth of experimental data describing the growth and remodelling of left and right ventricular tissue. The published literature is fragmented, with data reported via different experimental techniques using tissues harvested from a variety of animals, with different gender and ages. This prevents developing a continuum of data spanning birth to death, so limiting the potential that can be leveraged to aid computational modelling and simulations. In this study, equibiaxial, non-equibiaxial, and stress-relaxation data are presented, describing directional-dependent material responses. The biomechanical data is consolidated with equivalent microstructural data, an important element for the development of future material models. Combined, these data describe microstructural and biomechanical changes in the ventricles, spanning G &R from neonatal to adulthood.

Myocardial Fiber Mapping of Rat Hearts, Using Apparent Backscatter, with Histologic Validation

Myocardial fiber architecture is a physiologically important regulator of ejection fraction, strain and pressure development. Apparent ultrasonic backscatter has been shown to be a useful method for recreating the myocardial fiber architecture in human-sized sheep hearts because of the dependence of its amplitude on the relative orientation of a myofiber to the angle of ultrasonic insonification. Thus, the anisotropy of the backscatter signal is linked to and provides information about the fiber orientation. In this study, we sought to determine whether apparent backscatter could be used to measure myofiber orientation in rodent hearts. Fixed adult-rat hearts were imaged intact, and both a transmural cylindrical core and transmural wedge of the left ventricular free wall were imaged. Cylindrical core samples confirmed that backscatter anisotropy could be measured in rat hearts. Ultrasound and histologic analysis of transmural myocardial wedge samples confirmed that the apparent backscatter could be reproducibly mapped to fiber orientation (angle of the fiber relative to the direction of insonification). These data provided a quantitative relationship between the apparent backscatter and fiber angle that was applied to whole-heart images. Myocardial fiber architecture was successfully measured in rat hearts. Quantifying myocardial fiber architecture, using apparent backscatter, provides a number of advantages, including its scalable use from rodents to man, its rapid low-cost acquisition and minimal contraindications. The method outlined in this study provides a method for investigators to begin detailed assessments of how the myocardial fiber architecture changes in preclinical disease models, which can be immediately translated into the clinic.

A framework for biomechanics simulations using four-chamber cardiac models

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly-fewer elements are needed to build the cardiac geometry-compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.

Cardiac muscle organization revealed in 3-D by imaging whole-mount mouse hearts using two-photon fluorescence and confocal microscopy

The ability to image the entire adult mouse heart at high resolution in 3-D would provide enormous advantages in the study of heart disease. However, a technique for imaging nuclear/cellular detail as well as the overall structure of the entire heart in 3-D with minimal effort is lacking. To solve this problem, we modified the benzyl alcohol:benzyl benzoate (BABB) clearing technique by labeling mouse hearts with periodic acid Schiff (PAS) stain. We then imaged the hearts with a combination of two-photon fluorescence microscopy and automated tile-scan imaging/stitching. Utilizing the differential spectral properties of PAS, we could identify muscle and nuclear compartments in the heart. We were also able to visualize the differences between a 3-month-old normal mouse heart and a mouse heart that had undergone heart failure due to the expression of cardiac myosin binding protein-C (cMyBP-C) gene mutation (t/t). Using 2-D and 3-D morphometric analysis, we found that the t/t heart had anomalous ventricular shape, volume, and wall thickness, as well as a disrupted sarcomere pattern. We further validated our approach using decellularized hearts that had been cultured with 3T3 fibroblasts, which were tracked using a nuclear label. We were able to detect the 3T3 cells inside the decellularized intact heart tissue, achieving nuclear/cellular resolution in 3-D. The combination of labeling, clearing, and two-photon microscopy together with tiling eliminates laborious and time-consuming physical sectioning, alignment, and 3-D reconstruction.

Optical clearing based cellular-level 3D visualization of intact lymph node cortex

Lymph node (LN) is an important immune organ that controls adaptive immune responses against foreign pathogens and abnormal cells. To facilitate efficient immune function, LN has highly organized 3D cellular structures, vascular and lymphatic system. Unfortunately, conventional histological analysis relying on thin-sliced tissue has limitations in 3D cellular analysis due to structural disruption and tissue loss in the processes of fixation and tissue slicing. Optical sectioning confocal microscopy has been utilized to analyze 3D structure of intact LN tissue without physical tissue slicing. However, light scattering within biological tissues limits the imaging depth only to superficial portion of LN cortex. Recently, optical clearing techniques have shown enhancement of imaging depth in various biological tissues, but their efficacy for LN are remained to be investigated. In this work, we established optical clearing procedure for LN and achieved 3D volumetric visualization of the whole cortex of LN. More than 4 times improvement in imaging depth was confirmed by using LN obtained from H2B-GFP/actin-DsRed double reporter transgenic mouse. With adoptive transfer of GFP expressing B cells and DsRed expressing T cells and fluorescent vascular labeling by anti-CD31 and anti-LYVE-1 antibody conjugates, we successfully visualized major cellular-level structures such as T-cell zone, B-cell follicle and germinal center. Further, we visualized the GFP expressing metastatic melanoma cell colony, vasculature and lymphatic vessels in the LN cortex.

Analyzing Remodeling of Cardiac Tissue: A Comprehensive Approach Based on Confocal Microscopy and 3D Reconstructions

Microstructural characterization of cardiac tissue and its remodeling in disease is a crucial step in many basic research projects. We present a comprehensive approach for three-dimensional characterization of cardiac tissue at the submicrometer scale. We developed a compression-free mounting method as well as labeling and imaging protocols that facilitate acquisition of three-dimensional image stacks with scanning confocal microscopy. We evaluated the approach with normal and infarcted ventricular tissue. We used the acquired image stacks for segmentation, quantitative analysis and visualization of important tissue components. In contrast to conventional mounting, compression-free mounting preserved cell shapes, capillary lumens and extracellular laminas. Furthermore, the new approach and imaging protocols resulted in high signal-to-noise ratios at depths up to 60 µm. This allowed extensive analyzes revealing major differences in volume fractions and distribution of cardiomyocytes, blood vessels, fibroblasts, myofibroblasts and extracellular space in control vs. infarct border zone. Our results show that the developed approach yields comprehensive data on microstructure of cardiac tissue and its remodeling in disease. In contrast to other approaches, it allows quantitative assessment of all major tissue components. Furthermore, we suggest that the approach will provide important data for physiological models of cardiac tissue at the submicrometer scale.

Biomechanical properties and microstructure of human ventricular myocardium

In the multidisciplinary field of heart research it is of utmost importance to identify accurate myocardium material properties for the description of phenomena such as mechano-electric feedback or heart wall thickening. A rationally-based material model is required to understand the highly nonlinear mechanics of complex structures such as the passive myocardium under different loading conditions. Unfortunately, to date there are no experimental data of human heart tissues available to estimate material parameters and to develop adequate material models. This study aimed to determine biaxial extension and triaxial shear properties and the underlying microstructure of the passive human ventricular myocardium. Using new state-of-the-art equipment, planar biaxial extension tests were performed to determine the biaxial extension properties of the passive ventricular human myocardium. Shear properties of the myocardium were examined by triaxial simple shear tests performed on small cubic specimens excised from an adjacent region of the biaxial extension specimens. The three-dimensional microstructure was investigated through second-harmonic generation (SHG) microscopy on optically cleared tissues, which emphasized the 3D orientation and dispersion of the myofibers and adjacent collagen fabrics. The results suggest that the passive human LV myocardium under quasi-static and dynamic multiaxial loadings is a nonlinear, anisotropic (orthotropic), viscoelastic and history-dependent soft biological material undergoing large deformations. Material properties of the tissue components along local microstructural axes drive the nonlinear and orthotropic features of the myocardium. SHG microscopy investigation revealed detailed information about the myocardial microstructure due to its high resolution. It enabled the identification of structural parameters such as the fiber and the sheet orientations and corresponding dispersions. With this complete set of material data, a sophisticated material model and associated material parameters can be defined for a better description of the biomechanical response of the ventricular myocardium in humans. Such a model will lead to more accurate computational simulations to better understand the fundamental underlying ventricular mechanics, a step needed in the improvement of medical treatment of heart diseases.

STATEMENT OF SIGNIFICANCE

Unfortunately, to date there are no experimental data of human heart tissues available for material parameter estimation and the development of adequate material models. In this manuscript novel biaxial tensile and shear test data at different specimen orientations are presented, which allowed to adequately capture the direction-dependent material response. With these complete sets of mechanical data, combined with their underlying microstructural data (also presented herein), sophisticated material models and associated material parameters can be defined for the description of the mechanical behavior of the ventricular myocardium in humans. Such models will lead to accurate computational simulations to better understand the fundamental underlying ventricular mechanics, a step needed in the improvement of medical treatment of heart diseases.

Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

BACKGROUND

Cardiovascular magnetic resonance (CMR) can through the two methods 3D FLASH and diffusion tensor imaging (DTI) give complementary information on the local orientations of cardiomyocytes and their laminar arrays.

METHODS

Eight explanted rat hearts were perfused with Gd-DTPA contrast agent and fixative and imaged in a 9.4T magnet by two types of acquisition: 3D fast low angle shot (FLASH) imaging, voxels 50 × 50 × 50 μm, and 3D spin echo DTI with monopolar diffusion gradients of 3.6 ms duration at 11.5 ms separation, voxels 200 × 200 × 200 μm. The sensitivity of each approach to imaging parameters was explored.

RESULTS

The FLASH data showed laminar alignments of voxels with high signal, in keeping with the presumed predominance of contrast in the interstices between sheetlets. It was analysed, using structure-tensor (ST) analysis, to determine the most (v1(ST)), intermediate (v2(ST)) and least (v3(ST)) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e1(DTI)), intermediate (e2(DTI)) and least (e3(DTI)) orthogonal eigenvectors of extent of diffusion. The correspondence between the FLASH and DTI methods was measured and appraised. The most extended direction of FLASH signal (v1(ST)) agreed well with that of diffusion (e1(DTI)) throughout the left ventricle (representative discrepancy in the septum of 13.3 ± 6.7°: median ± absolute deviation) and both were in keeping with the expected local orientations of the long-axis of cardiomyocytes. However, the orientation of the least directions of FLASH signal continuity (v3(ST)) and diffusion (e3(ST)) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v3(ST)) and DTI (e3(DTI)) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v3(ST) and the directly measured FLASH laminar array normal was of 9 ± 7° for the lateral wall and 7 ± 9° for the septum (median ± inter quartile range), and for DTI the discrepancy between the calculated v3(DTI) and the directly measured FLASH laminar array normal was 22 ± 14° and 61 ± 53.4°. DTI was relatively insensitive to the number of diffusion directions and to time up to 72 hours post fixation, but was moderately affected by b-value (which was scaled by modifying diffusion gradient pulse strength with fixed gradient pulse separation). Optimal DTI parameters were b = 1000 mm/s(2) and 12 diffusion directions. FLASH acquisitions were relatively insensitive to the image processing parameters explored.

CONCLUSIONS

We show that ST analysis of FLASH is a useful and accurate tool in the measurement of cardiac microstructure. While both FLASH and the DTI approaches appear promising for mapping of the alignments of myocytes throughout myocardium, marked discrepancies between the cross myocyte anisotropies deduced from each method call for consideration of their respective limitations.

Cardiac response to low-energy field pacing challenges the standard theory of defibrillation

BACKGROUND

The electric response of myocardial tissue to periodic field stimuli has attracted significant attention as the basis for low-energy antifibrillation pacing, potentially more effective than traditional single high-energy shocks. In conventional models, an electric field produces a highly nonuniform response of the myocardial wall, with discrete excitations, or hot spots (HS), occurring at cathodal tissue surfaces or large coronary vessels. We test this prediction using novel 3-dimensional tomographic optical imaging.

METHODS AND RESULTS

Experiments were performed in isolated coronary perfused pig ventricular wall preparations stained with near-infrared voltage-sensitive fluorescent dye DI-4-ANBDQBS. The 3-dimensional coordinates of HS were determined using alternating transillumination. To relate HS formation with myocardial structures, we used ultradeep confocal imaging (interrogation depths, >4 mm). The peak HS distribution is located deep inside the heart wall, and the depth is not significantly affected by field polarity. We did not observe the strong colocalization of HS with major coronary vessels anticipated from theory. Yet, we observed considerable lateral displacement of HS with field polarity reversal. Models that de-emphasized lateral intracellular coupling and accounted for resistive heterogeneity in the extracellular space showed similar HS distributions to the experimental observations.

CONCLUSIONS

The HS distributions within the myocardial wall and the significant lateral displacements with field polarity reversal are inconsistent with standard theories of defibrillation. Extended theories based on enhanced descriptions of cellular scale electric mechanisms may be necessary. The considerable lateral displacement of HS with field polarity reversal supports the hypothesis of biphasic stimuli in low-energy antifibrillation pacing being advantageous.

Embryonic heart morphogenesis from confocal microscopy imaging and automatic segmentation

Embryonic heart morphogenesis (EHM) is a complex and dynamic process where the heart transforms from a single tube into a four-chambered pump. This process is of great biological and clinical interest but is still poorly understood for two main reasons. On the one hand, the existing imaging modalities for investigating EHM suffered from either limited penetration depth or limited spatial resolution. On the other hand, current works typically adopted manual segmentation, which was tedious, subjective, and time consuming considering the complexity of developing heart geometry and the large size of images. In this paper, we propose to utilize confocal microscopy imaging with tissue optical immersion clearing technique to image the heart at different stages of development for EHM study. The imaging method is able to produce high spatial resolution images and achieve large penetration depth at the same time. Furthermore, we propose a novel convex active contour model for automatic image segmentation. The model has the ability to deal with intensity fall-off in depth which is characterized by confocal microscopy images. We acquired the images of embryonic quail hearts from day 6 to day 14 of incubation for EHM study. The experimental results were promising and provided us with an insight view of early heart growth pattern and also paved the road for data-driven heart growth modeling.

Visualizing the complex 3D geometry of the perfusion border zone in isolated rabbit heart

Myocardial infarction, caused by a major blockage of a coronary artery, creates a border zone (BZ) between perfused and nonperfused tissue, which is believed to be the origin of fatal cardiac arrhythmias. We used a combination of optical clearing and polarization-sensitive optical coherence tomography to visualize a three-dimensional organization of the BZ in isolated rabbit hearts (n=5) at the microscopic level with a high spatial resolution. We found that the BZ has a complex three-dimensional structure with nonperfused areas penetrating into perfused tissue with finger-like projections. These "fingers" may play an important role in the initiation and maintenance of ventricular arrhythmias.

Extracting surface activation time from the optically recorded action potential in three-dimensional myocardium

Optical mapping has become an indispensible tool for studying cardiac electrical activity. However, due to the three-dimensional nature of the optical signal, the optical upstroke is significantly longer than the electrical upstroke. This raises the issue of how to accurately determine the activation time on the epicardial surface. The purpose of this study was to establish a link between the optical upstroke and exact surface activation time using computer simulations, with subsequent validation by a combination of microelectrode recordings and optical mapping experiments. To simulate wave propagation and associated optical signals, we used a hybrid electro-optical model. We found that the time of the surface electrical activation (t(E)) within the accuracy of our simulations coincided with the maximal slope of the optical upstroke (t(F)*) for a broad range of optical attenuation lengths. This was not the case when the activation time was determined at 50% amplitude (t(F50)) of the optical upstroke. The validation experiments were conducted in isolated Langendorff-perfused rat hearts and coronary-perfused pig left ventricles stained with either di-4-ANEPPS or the near-infrared dye di-4-ANBDQBS. We found that t(F)* was a more accurate measure of t(E) than was t(F50) in all experimental settings tested (P = 0.0002). Using t(F)* instead of t(F50) produced the most significant improvement in measurements of the conduction anisotropy and the transmural conduction time in pig ventricles.

Visualization and quantification of whole rat heart laminar structure using high-spatial resolution contrast-enhanced MRI

It has been shown by histology that cardiac myocytes are organized into laminae and this structure is important in function, both influencing the spread of electrical activation and enabling myocardial thickening in systole by laminar sliding. We have carried out high-spatial resolution three-dimensional MRI of the ventricular myolaminae of the entire volume of the isolated rat heart after contrast perfusion [dimeglumine gadopentate (Gd-DTPA)]. Four ex vivo rat hearts were perfused with Gd-DTPA and fixative and high-spatial resolution MRI was performed on a 9.4T MRI system. After MRI, cryosectioning followed by histology was performed. Images from MRI and histology were aligned, described, and quantitatively compared. In the three-dimensional MR images we directly show the presence of laminae and demonstrate that these are highly branching and are absent from much of the subepicardium. We visualized these MRI volumes to demonstrate laminar architecture and quantitatively demonstrated that the structural features observed are similar to those imaged in histology. We showed qualitatively and quantitatively that laminar architecture is similar in the four hearts. MRI can be used to image the laminar architecture of ex vivo hearts in three dimensions, and the images produced are qualitatively and quantitatively comparable with histology. We have demonstrated in the rat that: 1) laminar architecture is consistent between hearts; 2) myolaminae are absent from much of the subepicardium; and 3) although localized orthotropy is present throughout the myocardium, tracked myolaminae are branching structures and do not have a discrete identity.

Image-based models of cardiac structure in health and disease

Computational approaches to investigating the electromechanics of healthy and diseased hearts are becoming essential for the comprehensive understanding of cardiac function. In this article, we first present a brief review of existing image-based computational models of cardiac structure. We then provide a detailed explanation of a processing pipeline which we have recently developed for constructing realistic computational models of the heart from high resolution structural and diffusion tensor (DT) magnetic resonance (MR) images acquired ex vivo. The presentation of the pipeline incorporates a review of the methodologies that can be used to reconstruct models of cardiac structure. In this pipeline, the structural image is segmented to reconstruct the ventricles, normal myocardium, and infarct. A finite element mesh is generated from the segmented structural image, and fiber orientations are assigned to the elements based on DTMR data. The methods were applied to construct seven different models of healthy and diseased hearts. These models contain millions of elements, with spatial resolutions in the order of hundreds of microns, providing unprecedented detail in the representation of cardiac structure for simulation studies.

Models of cardiac tissue electrophysiology: progress, challenges and open questions

Models of cardiac tissue electrophysiology are an important component of the Cardiac Physiome Project, which is an international effort to build biophysically based multi-scale mathematical models of the heart. Models of tissue electrophysiology can provide a bridge between electrophysiological cell models at smaller scales, and tissue mechanics, metabolism and blood flow at larger scales. This paper is a critical review of cardiac tissue electrophysiology models, focussing on the micro-structure of cardiac tissue, generic behaviours of action potential propagation, different models of cardiac tissue electrophysiology, the choice of parameter values and tissue geometry, emergent properties in tissue models, numerical techniques and computational issues. We propose a tentative list of information that could be included in published descriptions of tissue electrophysiology models, and used to support interpretation and evaluation of simulation results. We conclude with a discussion of challenges and open questions.

Virtual sectioning of cardiac tissue relative to fiber orientation

Ventricular myocardium is composed of muscle fibers organised into a complex, branching, laminar (sheet-like) structure. The fibers run approximately parallel to the epicardial wall, but their orientation relative to the circumferential axis varies transmurally, rotating from around -70 degrees at the epicardium to +70 degrees at the endocardium. This ensures that any flat transmural imaging plane or histological section contains only a partial description of myocardial architecture, as fibers intersect with the image plane at a variety of angles depending on transmural depth. We have developed a new way of viewing microstructure that accounts for this variation. Extended-volume confocal 3-D images of normal rat left-ventricular wall have previously been acquired, with an approximate size of 4x1x1 mm3. The transmural fiber rotation is measured on planes parallel to the epicardium, and used to define a curvilinear coordinate system with a transmural axis, and a second axis defined relative to the local fiber orientation. Images extracted from the image volume on curvilinear planes derived from these axes reveal a consistent view of myocardial architecture.

Extracting intramural wavefront orientation from optical upstroke shapes in whole hearts

Information about intramural propagation of electrical excitation is crucial to understanding arrhythmia mechanisms in thick ventricular muscle. There is currently a controversy over whether it is possible to extract such information from the shape of the upstroke in optical mapping recordings. We show that even in the complex geometry of a whole guinea pig heart, optical upstroke morphology reveals the 3D wavefront orientation near the surface. To characterize the upstroke morphology, we use V(F)(*), the fractional level at which voltage-sensitive fluorescence, V(F), has maximal time derivative. Low values of V(F)(*)( approximately 0.2) indicate a wavefront moving away from the surface, high values of V(F)(*) ( approximately 0.6) a wavefront moving toward the surface, and intermediate values of V(F)(*) ( approximately 0.4) a wavefront moving parallel to the surface. We further performed computer simulations using Luo-Rudy II electrophysiology and a simplified 3D geometry. The simulated V(F)(*) maps for free wall and apical stimulations as well as for sinus rhythm are in good quantitative agreement with the averaged experimental results. Furthermore, computer simulations show that the effect of the curvature of the heart on wave propagation is negligible.