Chromatin Accessibility of Human Mitral Valves and Functional Assessment of MVP Risk Loci.

RATIONALE: Mitral valve prolapse (MVP) is a common valvopathy that leads to mitral insufficiency, heart failure, and sudden death. Functional genomic studies in mitral valves are needed to better characterize MVP-associated variants and target genes. OBJECTIVE: To establish the chromatin accessibility profiles and assess functionality of variants and narrow down target genes at MVP loci. METHODS AND RESULTS: We mapped the open chromatin regions in nuclei from 11 human pathogenic and 7 nonpathogenic mitral valves by an assay for transposase-accessible chromatin with high-throughput sequencing. Open chromatin peaks were globally similar between pathogenic and nonpathogenic valves. Compared with the heart tissue and cardiac fibroblasts, we found that MV-specific assay for transposase-accessible chromatin with high-throughput sequencing peaks are enriched near genes involved in extracellular matrix organization, chondrocyte differentiation, and connective tissue development. One of the most enriched motifs in MV-specific open chromatin peaks was for the nuclear factor of activated T cells family of TFs (transcription factors) involved in valve endocardial and interstitial cell formation. We also found that MVP-associated variants were significantly enriched (P<0.05) in mitral valve open chromatin peaks. Integration of the assay for transposase-accessible chromatin with high-throughput sequencing data with risk loci, extensive functional annotation, and gene reporter assay suggest plausible causal variants for rs2641440 at the SMG6/SRR locus and rs6723013 at the IGFBP2/IGFBP5/TNS1 locus. CRISPR-Cas9 deletion of the sequence including rs6723013 in human fibroblasts correlated with increased expression only for TNS1. Circular chromatin conformation capture followed by high-throughput sequencing experiments provided evidence for several target genes, including SRR, HIC1, and DPH1 at the SMG6/SRR locus and further supported TNS1 as the most likely target gene on chromosome 2. CONCLUSIONS: Here, we describe unprecedented genome-wide open chromatin profiles from human pathogenic and nonpathogenic MVs and report specific gene regulation profiles, compared with the heart. We also report in vitro functional evidence for potential causal variants and target genes at MVP risk loci involving established and new biological mechanisms. Graphic Abstract: A graphic abstract is available for this article.

Photocurable bioresorbable adhesives as functional interfaces between flexible bioelectronic devices and soft biological tissues.

Flexible electronic/optoelectronic systems that can intimately integrate onto the surfaces of vital organ systems have the potential to offer revolutionary diagnostic and therapeutic capabilities relevant to a wide spectrum of diseases and disorders. The critical interfaces between such technologies and living tissues must provide soft mechanical coupling and efficient optical/electrical/chemical exchange. Here, we introduce a functional adhesive bioelectronic-tissue interface material, in the forms of mechanically compliant, electrically conductive, and optically transparent encapsulating coatings, interfacial layers or supporting matrices. These materials strongly bond both to the surfaces of the devices and to those of different internal organs, with stable adhesion for several days to months, in chemistries that can be tailored to bioresorb at controlled rates. Experimental demonstrations in live animal models include device applications that range from battery-free optoelectronic systems for deep-brain optogenetics and subdermal phototherapy to wireless millimetre-scale pacemakers and flexible multielectrode epicardial arrays. These advances have immediate applicability across nearly all types of bioelectronic/optoelectronic system currently used in animal model studies, and they also have the potential for future treatment of life-threatening diseases and disorders in humans.

Spatial organization of acute myocardial ischemia.

INTRODUCTION: Myocardial ischemia is a pathological condition initiated by supply and demand imbalance of the blood to the heart. Previous studies suggest that ischemia originates in the subendocardium, i.e., that nontransmural ischemia is limited to the subendocardium. By contrast, we hypothesized that acute myocardial ischemia is not limited to the subendocardium and sought to document its spatial distribution in an animal preparation. The goal of these experiments was to investigate the spatial organization of ischemia and its relationship to the resulting shifts in ST segment potentials during short episodes of acute ischemia. METHODS: We conducted acute ischemia studies in open-chest canines (N=19) and swines (N=10), which entailed creating carefully controlled ischemia using demand, supply or complete occlusion ischemia protocols and recording intramyocardial and epicardial potentials. Elevation of the potentials at 40% of the ST segment between the J-point and the peak of the T-wave (ST40%) provided the metric for local ischemia. The threshold for ischemic ST segment elevations was defined as two standard deviations away from the baseline values. RESULTS: The relative frequency of occurrence of acute ischemia was higher in the subendocardium (78% for canines and 94% for swines) and the mid-wall (87% for canines and 97% for swines) in comparison with the subepicardium (30% for canines and 22% for swines). In addition, acute ischemia was seen arising throughout the myocardium (distributed pattern) in 87% of the canine and 94% of the swine episodes. Alternately, acute ischemia was seen originating only in the subendocardium (subendocardial pattern) in 13% of the canine episodes and 6% of the swine episodes (p<0.05). CONCLUSIONS: Our findings suggest that the spatial distribution of acute ischemia is a complex phenomenon arising throughout the myocardial wall and is not limited to the subendocardium.

Experimental Data and Geometric Analysis Repository-EDGAR.

INTRODUCTION: The "Experimental Data and Geometric Analysis Repository", or EDGAR is an Internet-based archive of curated data that are freely distributed to the international research community for the application and validation of electrocardiographic imaging (ECGI) techniques. The EDGAR project is a collaborative effort by the Consortium for ECG Imaging (CEI, ecg-imaging.org), and focused on two specific aims. One aim is to host an online repository that provides access to a wide spectrum of data, and the second aim is to provide a standard information format for the exchange of these diverse datasets. METHODS: The EDGAR system is composed of two interrelated components: 1) a metadata model, which includes a set of descriptive parameters and information, time signals from both the cardiac source and body-surface, and extensive geometric information, including images, geometric models, and measure locations used during the data acquisition/generation; and 2) a web interface. This web interface provides efficient, search, browsing, and retrieval of data from the repository. RESULTS: An aggregation of experimental, clinical and simulation data from various centers is being made available through the EDGAR project including experimental data from animal studies provided by the University of Utah (USA), clinical data from multiple human subjects provided by the Charles University Hospital (Czech Republic), and computer simulation data provided by the Karlsruhe Institute of Technology (Germany). CONCLUSIONS: It is our hope that EDGAR will serve as a communal forum for sharing and distribution of cardiac electrophysiology data and geometric models for use in ECGI research.

Sensitivity of epicardial electrical markers to acute ischemia detection.

INTRODUCTION: We hypothesize that electrocardiographic measurements from the intramyocardial space contain more sensitive markers of ischemia than those detectable on the epicardium. The goal of this study was to evaluate different electrical markers for their potential to detect the earliest phases of acute myocardial ischemia. METHODS: We conducted acute ischemia studies in open chest animal, by creating finely controlled demand or supply ischemic episodes and recording intramyocardial and epicardial potentials. RESULTS: Under the conditions of mild perfusion deficit, acute ischemia induced changes in the T wave that were larger and could be detected earlier on the epicardial surface than ST-segment changes. CONCLUSIONS: Our findings indicate that in the setting of very acute ischemia, epicardial T waves have higher sensitivity to mild degrees of acute ischemia than epicardial ST potentials. These results suggest that changes in the T wave shape may augment shifts in ST segments to improve ECG based localization of ischemia.

Effect of Spatial Resolution on Accurate Detection and Localization of Arrhythmia Rotors in Human Right Ventricular Tachycardia.

The goal of this study was to identify the spatial resolution requirements for accurate rotor detection and localization in human right ventricular tachyarrhythmias. Poor spatial resolution is often cited as a reason for the inaccuracy of cardiac mapping catheters in detecting and localizing arrhythmia rotors. High-resolution (0.7 mm) arrhythmia data from optical recordings obtained from human donor hearts (n = 12) were uniformly downsampled to lower resolutions (1.4-7 mm) to approximate the spatial resolution (4 mm) of clinical mapping catheters. Rotors were tracked at various subresolutions and compared to the rotors in the original data by computing F1-scores to create accuracy profiles for both rotor detection and localization. Further comparisons were made according to arrhythmia type, donor sex, anatomical region, and mapped surface: endocardium or epicardium. For a spatial resolution of 4.2 mm, the accuracies of rotor detection and localization were 57% ± 4% and 61% ± 7%, respectively. Arrhythmia type affected the accuracy of rotor detection (monomorphic ventricular tachycardia, 58% ± 4%; ventricular fibrillation, 56% ± 8%) and localization (monomorphic ventricular tachycardia, 70% ± 4%; ventricular fibrillation, 54% ± 13%). However, donor sex, anatomical region (right ventricular outflow tract, mid, and apical), and mapped surface (epicardium and endocardium) did not significantly affect rotor detection or localization accuracy. To achieve rotor detection accuracy of 80%, a spatial resolution of 1.4 mm or better is needed. The accuracy profiles provided here serve as a guideline for future mapping device development.

Hardware-Mappable Cellular Neural Networks for Distributed Wavefront Detection in Next-Generation Cardiac Implants.

Artificial intelligence algorithms are being adopted to analyze medical data, promising faster interpretation to support doctors' diagnostics. The next frontier is to bring these powerful algorithms to implantable medical devices. Herein, a closed-loop solution is proposed, where a cellular neural network is used to detect abnormal wavefronts and wavebrakes in cardiac signals recorded in human tissue is trained to achieve >96% accuracy, >92% precision, >99% specificity, and >93% sensitivity, when floating point precision weights are assumed. Unfortunately, the current hardware technologies for floating point precision are too bulky or energy intensive for compact standalone applications in medical implants. Emerging device technologies, such as memristors, can provide the compact and energy-efficient hardware fabric to support these efforts and can be reliably embedded with existing sensor and actuator platforms in implantable devices. A distributed design that considers the hardware limitations in terms of overhead and limited bit precision is also discussed. The proposed distributed solution can be easily adapted to other medical technologies that require compact and efficient computing, like wearable devices and lab-on-chip platforms.

Secretome of atrial epicardial adipose tissue facilitates reentrant arrhythmias by myocardial remodeling.

BACKGROUND: Epicardial adipose tissue (EAT) accumulation is associated with cardiac arrhythmias. The effect of EAT secretome (EATs) on cardiac electrophysiology remains largely unknown. OBJECTIVE: The purpose of this study was to investigate the arrhythmogenicity of EATs and its underlying molecular and electrophysiological mechanisms. METHODS: We collected atrial EAT and subcutaneous adipose tissue (SAT) from 30 patients with atrial fibrillation (AF), and EAT from 3 donors without AF. The secretome was collected after a 24-hour incubation of the adipose tissue explants. We cultured neonatal rat ventricular myocytes (NRVMs) with EATs, subcutaneous adipose tissue secretome (SATs), and cardiomyocytes conditioned medium (CCM) for 72 hours. We implemented the electrophysiological changes observed after EATs incubation into a model of human left atrium and tested arrhythmia inducibility. RESULTS: Incubation of NRVMs with EATs decreased expression of the potassium channel subunit Kcnj2 by 26% and correspondingly reduced the inward rectifier K+ current IK1 by 35% compared to incubation with CCM, resulting in a depolarized resting membrane of cardiomyocytes. EATs decreased expression of connexin43 (29% mRNA, 46% protein) in comparison to CCM. Cells incubated with SATs showed no significant differences in Kcnj2 or Gja1 expression in comparison to CCM, and their resting potential was not depolarized. Cardiomyocytes incubated with EATs showed reduced conduction velocity and increased conduction heterogeneity compared to SATs and CCM. Computer modeling of human left atrium revealed that the electrophysiological changes induced by EATs promote sustained reentrant arrhythmias if EAT partially covers the myocardium. CONCLUSION: EAT slows conduction, depolarizes the resting potential, alters electrical cell-cell coupling, and facilitates reentrant arrhythmias.

Electrophysiology and Arrhythmogenesis in the Human Right Ventricular Outflow Tract.

BACKGROUND: Right ventricular outflow tract (RVOT) is a common source of ventricular tachycardia, which often requires ablation. However, the mechanisms underlying the RVOT's unique arrhythmia susceptibility remain poorly understood due to lack of detailed electrophysiological and molecular studies of the human RVOT. METHODS: We conducted optical mapping studies in 16 nondiseased donor human RVOT preparations subjected to pharmacologically induced adrenergic and cholinergic stimulation to evaluate susceptibility to arrhythmias and characterize arrhythmia dynamics. RESULTS: We found that under control conditions, RVOT has shorter action potential duration at 80% repolarization relative to the right ventricular apical region. Treatment with isoproterenol (100 nM) shortened action potential duration at 80% repolarization and increased incidence of premature ventricular contractions (P=0.003), whereas acetylcholine (100 µM) stimulation alone had no effect on action potential duration at 80% repolarization or premature ventricular contractions. However, acetylcholine treatment after isoproterenol stimulation reduced the incidence of premature ventricular contractions (P=0.034) and partially reversed action potential duration at 80% repolarization shortening (P=0.029). Immunolabeling of RVOT (n=4) confirmed the presence of cholinergic marker VAChT (vesicular acetylcholine transporter) in the region. Rapid pacing revealed RVOT susceptibility to both concordant and discordant alternans. Investigation into transmural arrhythmia dynamics showed that arrhythmia wave fronts and phase singularities (rotors) were relatively more organized in the endocardium than in the epicardium (P=0.006). Moreover, there was a weak but positive spatiotemporal autocorrelation between epicardial and endocardial arrhythmic wave fronts and rotors. Transcriptome analysis (n=10 hearts) suggests a trend that MAPK (mitogen-activated protein kinase) signaling, calcium signaling, and cGMP-PKG (protein kinase G) signaling are among the pathways that may be enriched in the male RVOT, whereas pathways of neurodegeneration may be enriched in the female RVOT. CONCLUSIONS: Human RVOT electrophysiology is characterized by shorter action potential duration relative to the right ventricular apical region. Cholinergic right ventricular stimulation attenuates the arrhythmogenic effects of adrenergic stimulation, including increase in frequency of premature ventricular contractions and shortening of wavelength. Right ventricular arrhythmia is characterized by positive spatial-temporal autocorrelation between epicardial-endocardial arrhythmic wave fronts and rotors that are relatively more organized in the endocardium.

A transient, closed-loop network of wireless, body-integrated devices for autonomous electrotherapy.

Temporary postoperative cardiac pacing requires devices with percutaneous leads and external wired power and control systems. This hardware introduces risks for infection, limitations on patient mobility, and requirements for surgical extraction procedures. Bioresorbable pacemakers mitigate some of these disadvantages, but they demand pairing with external, wired systems and secondary mechanisms for control. We present a transient closed-loop system that combines a time-synchronized, wireless network of skin-integrated devices with an advanced bioresorbable pacemaker to control cardiac rhythms, track cardiopulmonary status, provide multihaptic feedback, and enable transient operation with minimal patient burden. The result provides a range of autonomous, rate-adaptive cardiac pacing capabilities, as demonstrated in rat, canine, and human heart studies. This work establishes an engineering framework for closed-loop temporary electrotherapy using wirelessly linked, body-integrated bioelectronic devices.

Validating defibrillation simulation in a human-shaped phantom.

BACKGROUND: We previously developed a computational model to aid clinicians in positioning implantable cardioverter-defibrillators (ICDs), especially in the case of abnormal anatomies that commonly arise in pediatric cases. We have validated the model clinically on the body surface; however, validation within the volume of the heart is required to establish complete confidence in the model and improve its use in clinical settings. OBJECTIVE: The goal of this study was to use an animal model and thoracic phantom to record the ICD potential field within the heart and on the torso to validate our defibrillation simulation system. METHODS: We recorded defibrillator shock potentials from an ICD suspended together with an animal heart in a human-shaped torso tank and compared them with simulated values. We also compared the scaled distribution threshold, an analog to the defibrillation threshold, calculated from the measured and simulated electric fields within the myocardium. RESULTS: ICD potentials recorded on the tank and cardiac surface and within the myocardium agreed well with those predicted by the simulation. A quantitative comparison of the recorded and simulated potentials yielded a mean correlation of 0.94 and a relative error of 19.1%. The simulation can also predict scaled distribution thresholds similar to those calculated from the measured potential fields. CONCLUSION: We found that our simulation could predict potential fields with high correlation with the measured values within the heart and on the torso surface. These results support the use of this model for the optimization of ICD placements.

Correction: Genetic algorithm-based personalized models of human cardiac action potential.

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

Genetic algorithm-based personalized models of human cardiac action potential.

We present a novel modification of genetic algorithm (GA) which determines personalized parameters of cardiomyocyte electrophysiology model based on set of experimental human action potential (AP) recorded at different heart rates. In order to find the steady state solution, the optimized algorithm performs simultaneous search in the parametric and slow variables spaces. We demonstrate that several GA modifications are required for effective convergence. Firstly, we used Cauchy mutation along a random direction in the parametric space. Secondly, relatively large number of elite organisms (6-10% of the population passed on to new generation) was required for effective convergence. Test runs with synthetic AP as input data indicate that algorithm error is low for high amplitude ionic currents (1.6±1.6% for IKr, 3.2±3.5% for IK1, 3.9±3.5% for INa, 8.2±6.3% for ICaL). Experimental signal-to-noise ratio above 28 dB was required for high quality GA performance. GA was validated against optical mapping recordings of human ventricular AP and mRNA expression profile of donor hearts. In particular, GA output parameters were rescaled proportionally to mRNA levels ratio between patients. We have demonstrated that mRNA-based models predict the AP waveform dependence on heart rate with high precision. The latter also provides a novel technique of model personalization that makes it possible to map gene expression profile to cardiac function.

Granger Causality-Based Analysis for Classification of Fibrillation Mechanisms and Localization of Rotational Drivers.

BACKGROUND: The mechanisms sustaining myocardial fibrillation remain disputed, partly due to a lack of mapping tools that can accurately identify the mechanism with low spatial resolution clinical recordings. Granger causality (GC) analysis, an econometric tool for quantifying causal relationships between complex time-series, was developed as a novel fibrillation mapping tool and adapted to low spatial resolution sequentially acquired data. METHODS: Ventricular fibrillation (VF) optical mapping was performed in Langendorff-perfused Sprague-Dawley rat hearts (n=18), where novel algorithms were developed using GC-based analysis to (1) quantify causal dependence of neighboring signals and plot GC vectors, (2) quantify global organization with the causality pairing index, a measure of neighboring causal signal pairs, and (3) localize rotational drivers (RDs) by quantifying the circular interdependence of neighboring signals with the circular interdependence value. GC-based mapping tools were optimized for low spatial resolution from downsampled optical mapping data, validated against high-resolution phase analysis and further tested in previous VF optical mapping recordings of coronary perfused donor heart left ventricular wedge preparations (n=12), and adapted for sequentially acquired intracardiac electrograms during human persistent atrial fibrillation mapping (n=16). RESULTS: Global VF organization quantified by causality pairing index showed a negative correlation at progressively lower resolutions (50% resolution: P=0.006, R2=0.38, 12.5% resolution, P=0.004, R2=0.41) with a phase analysis derived measure of disorganization, locations occupied by phase singularities. In organized VF with high causality pairing index values, GC vector mapping characterized dominant propagating patterns and localized stable RDs, with the circular interdependence value showing a significant difference in driver versus nondriver regions (0.91±0.05 versus 0.35±0.06, P=0.0002). These findings were further confirmed in human VF. In persistent atrial fibrillation, a positive correlation was found between the causality pairing index and presence of stable RDs (P=0.0005,R2=0.56). Fifty percent of patients had RDs, with a low incidence of 0.9±0.3 RDs per patient. CONCLUSIONS: GC-based fibrillation analysis can measure global fibrillation organization, characterize dominant propagating patterns, and map RDs using low spatial resolution sequentially acquired data.

Catheter-integrated soft multilayer electronic arrays for multiplexed sensing and actuation during cardiac surgery.

The rigidity and relatively primitive modes of operation of catheters equipped with sensing or actuation elements impede their conformal contact with soft-tissue surfaces, limit the scope of their uses, lengthen surgical times and increase the need for advanced surgical skills. Here, we report materials, device designs and fabrication approaches for integrating advanced electronic functionality with catheters for minimally invasive forms of cardiac surgery. By using multiphysics modelling, plastic heart models and Langendorff animal and human hearts, we show that soft electronic arrays in multilayer configurations on endocardial balloon catheters can establish conformal contact with curved tissue surfaces, support high-density spatiotemporal mapping of temperature, pressure and electrophysiological parameters and allow for programmable electrical stimulation, radiofrequency ablation and irreversible electroporation. Integrating multimodal and multiplexing capabilities into minimally invasive surgical instruments may improve surgical performance and patient outcomes.

Critical Volume of Human Myocardium Necessary to Maintain Ventricular Fibrillation.

BACKGROUND: Abnormal QT intervals, long QT or short QT, have been epidemiologically linked with sudden cardiac death because of ventricular fibrillation (VF). Consequently, Food and Drug Administration recommends testing all pharmacological agents for QT toxicity as a risk factor for cardiac toxicity. Such tests assess QT/QTc interval, which represents ventricular depolarization and repolarization. However, the current QT toxicity analysis does not account for the well-known anisotropy in cardiac tissue conductivity. Mines demonstrated in 1913 that cardiac wavelength (λ) determines inducibility of reentrant arrhythmia, where both repolarization time or action potential duration and conduction velocity determine λ=action potential duration×conduction velocity. We aimed to determine the role of anisotropic wavelength in inducibility of VF in explanted human left ventricular preparations. We tested the hypothesis that 3-dimensional cardiac wavelength, which takes into account anisotropic cardiac tissue conductivity, can accurately predict VF sustainability. METHODS: We conducted panoramic optical mapping of coronary perfused human left ventricular wedge preparations subjected to pharmacologically induced shortening and prolongation of action potential duration, by IK,ATP agonist pinacidil and antagonist glybenclamide, respectively. This measured action potential duration, conduction velocity, and thus determined pacing cycle length-dependent wavelengths in longitudinal (λL), transverse (λTV), and transmural (λTM) directions using S1S1 pacing protocol, from which wavelength volume (Vλ) was determined, as Vλ=λL×λTV×λTM, and compared with tissue volume. We tested a hypothesis that tissue volume/Vλ ratio can predict VF sustainability. RESULTS: At baseline, at pacing rate of 240 beats per minute, the wavelengths were λL=9.6±0.6 cm, λTV=4.2±0.3 cm, and λTM=5.8±0.2 cm, respectively (n=7), and thus Vλ=246.4±42.1 cm3. Administration of pinacidil at escalating concentrations progressively decreased Vλ, and VF became sustained, when tissue volume/Vλ was above safety factor κ=4.4±0.6 (n=9) during rapid pacing. Treatment with glybenclamide decreased VT/Vλ below κ at any pacing rate and prevented VF sustainability. CONCLUSIONS: Sustained VF was only sustained in ventricular volume exceeding critical Vλ=λL×λTV×λTM.

A Framework for Image-Based Modeling of Acute Myocardial Ischemia Using Intramurally Recorded Extracellular Potentials.

The biophysical basis for electrocardiographic evaluation of myocardial ischemia stems from the notion that ischemic tissues develop, with relative uniformity, along the endocardial aspects of the heart. These injured regions of subendocardial tissue give rise to intramural currents that lead to ST segment deflections within electrocardiogram (ECG) recordings. The concept of subendocardial ischemic regions is often used in clinical practice, providing a simple and intuitive description of ischemic injury; however, such a model grossly oversimplifies the presentation of ischemic disease-inadvertently leading to errors in ECG-based diagnoses. Furthermore, recent experimental studies have brought into question the subendocardial ischemia paradigm suggesting instead a more distributed pattern of tissue injury. These findings come from experiments and so have both the impact and the limitations of measurements from living organisms. Computer models have often been employed to overcome the constraints of experimental approaches and have a robust history in cardiac simulation. To this end, we have developed a computational simulation framework aimed at elucidating the effects of ischemia on measurable cardiac potentials. To validate our framework, we simulated, visualized, and analyzed 226 experimentally derived acute myocardial ischemic events. Simulation outcomes agreed both qualitatively (feature comparison) and quantitatively (correlation, average error, and significance) with experimentally obtained epicardial measurements, particularly under conditions of elevated ischemic stress. Our simulation framework introduces a novel approach to incorporating subject-specific, geometric models and experimental results that are highly resolved in space and time into computational models. We propose this framework as a means to advance the understanding of the underlying mechanisms of ischemic disease while simultaneously putting in place the computational infrastructure necessary to study and improve ischemia models aimed at reducing diagnostic errors in the clinic.

RHYTHM: An Open Source Imaging Toolkit for Cardiac Panoramic Optical Mapping.

Fluorescence optical imaging techniques have revolutionized the field of cardiac electrophysiology and advanced our understanding of complex electrical activities such as arrhythmias. However, traditional monocular optical mapping systems, despite having high spatial resolution, are restricted to a two-dimensional (2D) field of view. Consequently, tracking complex three-dimensional (3D) electrical waves such as during ventricular fibrillation is challenging as the waves rapidly move in and out of the field of view. This problem has been solved by panoramic imaging which uses multiple cameras to measure the electrical activity from the entire epicardial surface. However, the diverse engineering skill set and substantial resource cost required to design and implement this solution have made it largely inaccessible to the biomedical research community at large. To address this barrier to entry, we present an open source toolkit for building panoramic optical mapping systems which includes the 3D printing of perfusion and imaging hardware, as well as software for data processing and analysis. In this paper, we describe the toolkit and demonstrate it on different mammalian hearts: mouse, rat, and rabbit.