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.

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.

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.

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.