A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network.

BACKGROUND: The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice. METHODS: FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription. RESULTS: FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency. CONCLUSIONS: Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Evidence of sex differences in cancer-related cardiac complications in mouse models of pancreatic and liver cancer.

Abnormal heart rate variability (HRV) is commonly observed in cancer patients who have undergone targeted therapy and/or surgery, yet the effects of cancer itself on cardiac function remain underexplored. Specifically, there is limited knowledge about sex-specific manifestations of HRV in cancer patients. Transgenic mouse models are widely used to study different types of cancer. Here, we aimed to investigate the sex-specific effects of cancer on cardiac function using transgenic mouse models of pancreatic and liver cancers. This study used male and female transgenic mice with cancer and wild-type controls. Cardiac function was assessed by recording electrocardiograms in conscious mice. RR intervals were detected to determine HRV using time and frequency domain analyses. Histological analysis with Masson's trichrome staining was performed to determine structural changes. In females, increased HRV was observed in both pancreatic and liver cancer-bearing mice. In contrast, in males, increased HRV was observed only in the liver cancer group. Male pancreatic cancer mice demonstrated autonomic balance shift showing an increase in parasympathetic to sympathetic tone. The heart rate (HR) was higher in control and liver cancer male mice groups than in females. Histological analysis did not show significant sex differences but suggested a higher degree of remodeling in liver cancer mice than in control, specifically in the right atrium and left ventricle. This study revealed sex differences in cancer's HR modulation. Specifically, female cancer mice had lower median HR and higher HRV. These findings indicate that sex must be considered when using HRV as a cancer biomarker.

Molecular and Functional Remodeling of Superior and Inferior SAN in a Rat Model of HCM.

BACKGROUND: Recently, our laboratory presented functional and molecular evidence for the presence of 2 competing sinoatrial node (SAN) pacemakers in healthy human and rat hearts. Anatomically localized near the superior vena cava and inferior vena cava, the superior and inferior SANs (sSAN and iSAN, respectively) preferentially control fast and slow normal heart rates. However, only 1 dominant pacemaker, primarily the sSAN, was functional in the failing rat heart with hypertrophic cardiomyopathy. OBJECTIVES: This study aimed to determine the transcriptional basis of functional silencing of 1 of 2 dominant pacemakers in failing rat hearts. METHODS: Ascending aortic constriction was performed on 1-week-old male Sprague-Dawley rat pups to induce left ventricular hypertrophy and heart failure. The dominant pacemaker was anatomically mapped in adult (10-12 weeks old) healthy and failing rat hearts using optical mapping in isolated right atrial tissue preparations. RNA sequencing was used to identify regional sSAN/iSAN gene expression differences between healthy and failing rat hearts. RESULTS: In all failing rat hearts optically mapped in this study (n = 4), only the sSAN pacemaker was functional, while the iSAN was silent. Compared to healthy rat hearts, a total of 3,640 genes were downregulated, and 4,518 genes were upregulated in failing rat hearts. The functional quiescence of the iSAN in these failing rat hearts may be explained by their downregulation of sodium, potassium, and calcium ion channels as well as their downregulation of specific structural genes, including ankyrin, titin, and myosin heavy chain. Moreover, the iSAN showed predominant downregulation of several key transcription factors such as Tbx5, Tbx3, Shox2, and Smad9. CONCLUSIONS: Pressure-overload-induced heart failure resulted in significant downregulation of critical transcription factors, ion channels, and structural transcripts of the iSAN, which could explain the functional silencing of the iSAN in failing rat hearts.

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.

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

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

Evidence of Superior and Inferior Sinoatrial Nodes in the Mammalian Heart.

OBJECTIVES: This study sought to investigate the shift of leading pacemaker locations in healthy and failing mammalian hearts over the entire range of physiological heart rates (HRs), and to molecularly characterize spatial regions of spontaneous activity. BACKGROUND: A normal heartbeat originates as an action potential in a group of pacemaker cells known as the sinoatrial node (SAN), located near the superior vena cava. HRs and the anatomical site of origin of pacemaker activity in the adult heart are known to dynamically change in response to various physiological inputs, yet the mechanism of this pacemaker shift is not well understood. METHODS: Optical mapping was applied to ex vivo rat and human isolated right atrial tissues, and HRs were modulated with acetylcholine and isoproterenol. RNA sequencing was performed on tissue areas that elicited spontaneous activity, and comparisons were made to neighboring myocardial tissues. RESULTS: Functional and molecular evidence identified and confirmed the presence of 2 competing right atrial pacemakers localized near the superior vena cava and the inferior vena cava-the superior SAN (sSAN) and inferior SAN (iSAN), respectively-which preferentially control the fast and slow HRs. Both of these regions were evident in non-failing rat and human hearts and maintained spontaneous activity in the rat heart when physically separated from one another. Molecular analysis of these 2 pacemaker regions revealed unique but similar transcriptional profiles, suggesting iSAN dominance when the sSAN is silent. CONCLUSIONS: The presence of 2 spatially distinct dominant pacemakers, sSAN and iSAN, in the mammalian heart clarifies previous identification of migrating pacemakers and corresponding changes in P-wave morphology in mammalian species.

Heart slice culture system reliably demonstrates clinical drug-related cardiotoxicity.

The limited availability of human heart tissue and its complex cell composition are major limiting factors for the reliable testing of drug efficacy and toxicity. Recently, we developed functional human and pig heart slice biomimetic culture systems that preserve the viability and functionality of 300 µm heart slices for up to 6 days. Here, we tested the reliability of this culture system for testing the cardiotoxicity of anti-cancer drugs. We tested three anti-cancer drugs (doxorubicin, trastuzumab, and sunitinib) with known different mechanisms of cardiotoxicity at three concentrations and assessed the effect of these drugs on heart slice viability, structure, function and gene expression. Slices incubated with any of these drugs for 48 h showed diminished in viability as well as loss of cardiomyocyte structure and function. Mechanistically, RNA sequencing of doxorubicin-treated tissues demonstrated a significant downregulation of cardiac genes and upregulation of oxidative stress responses. Trastuzumab treatment downregulated cardiac muscle contraction-related genes consistent with its clinically known effect on cardiomyocytes. Interestingly, sunitinib treatment resulted in significant downregulation of angiogenesis-related genes, in line with its mechanism of action. Similar to hiPS-derived-cardiomyocytes, heart slices recapitulated the expected toxicity of doxorubicin and trastuzumab, however, slices were superior in detecting sunitinib cardiotoxicity and mechanism in the clinically relevant concentration range of 0.1-1 µM. These results indicate that heart slice culture models have the potential to become a reliable platform for testing and elucidating mechanisms of drug cardiotoxicity.

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.

Mitochondrial Enzymes of the Urea Cycle Cluster at the Inner Mitochondrial Membrane.

Mitochondrial enzymes involved in energy transformation are organized into multiprotein complexes that channel the reaction intermediates for efficient ATP production. Three of the mammalian urea cycle enzymes: N-acetylglutamate synthase (NAGS), carbamylphosphate synthetase 1 (CPS1), and ornithine transcarbamylase (OTC) reside in the mitochondria. Urea cycle is required to convert ammonia into urea and protect the brain from ammonia toxicity. Urea cycle intermediates are tightly channeled in and out of mitochondria, indicating that efficient activity of these enzymes relies upon their coordinated interaction with each other, perhaps in a cluster. This view is supported by mutations in surface residues of the urea cycle proteins that impair ureagenesis in the patients, but do not affect protein stability or catalytic activity. We find the NAGS, CPS1, and OTC proteins in liver mitochondria can associate with the inner mitochondrial membrane (IMM) and can be co-immunoprecipitated. Our in-silico analysis of vertebrate NAGS proteins, the least abundant of the urea cycle enzymes, identified a protein-protein interaction region present only in the mammalian NAGS protein-"variable segment," which mediates the interaction of NAGS with CPS1. Use of super resolution microscopy showed that NAGS, CPS1 and OTC are organized into clusters in the hepatocyte mitochondria. These results indicate that mitochondrial urea cycle proteins cluster, instead of functioning either independently or in a rigid multienzyme complex.