Impact of chronic fetal hypoxia and inflammation on cardiac pacemaker cell development
11 Impact of chronic fetal hypoxia and inflammation on cardiac pacemaker cell development
Martin G. Frasch , Dino A. Giussani Department of Obstetrics and Gynecology, University of Washington, Seattle, WA, USA Center on Human Development and Disability, University of Washington, Seattle, WA, USA Department of Physiology, Development & Neuroscience, University of Cambridge, Cambridge, UK
Corresponding author:
Dr. Martin G. Frasch Department of Obstetrics and Gynecology University of Washington 1959 NE Pacific St Box 356460 Seattle, WA 98195 Email: [email protected]
Abstract
Chronic fetal hypoxia and infection are examples of adverse conditions during complicated pregnancy, which impact cardiac myogenesis and increase the lifetime risk of heart disease. However, the effects that chronic hypoxic or inflammatory environments exert on cardiac pacemaker cells are poorly understood. Here, we review the current evidence and novel avenues of bench-to-bed research in this field of perinatal cardiogenesis as well as its translational significance for early detection of future risk for cardiovascular disease.
Introduction
A suboptimal prenatal environment, such as fetal exposure to diminished nutrient or oxygen supply, impacts myocardial development predisposing to cardiometabolic disorders in later life.[1,2] There is however no direct data during development so far on the timing of the emergence of the cardiac pacemaker cell synchronization (CPS) nor on the effects of adverse intrauterine conditions on cardiac pacemaker cell development and alterations of CPS; the basis for an altered fetal heartbeat. Interestingly, some of the recent insight into cardiac pacemaker cell development comes from bioengineering. This has focused on the creation of biological pacemakers [3,4], the biology of the human pluripotent stem cells differentiating into adult cardiomyocytes [5], as well as transcriptomics associated with cardiac mature pacemaker cell function at single-cell resolution.[6] The analysis of fetal heart rate variability (HRV) has been the mainstay of noninvasive fetal health monitoring for decades, yet its utility to identify endangered fetuses has remained limited.[7,8] In this review, we will discuss how the CPS provides the subcellular, cellular and multicellular substrate for the emergence of an intrinsic fetal HRV (iHRV). We hope that expanding research into the iHRV as a proxy to fetal cardiac health will not only help in the diagnosis of fetuses at risk in adverse pregnancy, but that it will also inspire future studies into this aspect of cardiac cellular development and its significance for the early detection of cardiac diseases in later life.[9,10] The heartbeat is generated by the sinoatrial node, a small region in the right atrium of the heart, comprised of approximately 10,000 cells that interact with each other to set the normal cardiac rhythm. The fluctuations in heart rate result from a combination of its own intrinsic variability and the modulation of the pacemaker activity by several neural and endocrine control mechanisms (Fig. 1).[11] Remarkably, the intrinsic fluctuations of the heartbeat are seen not only in the isolated node, devoid of all neural and hormonal inputs but also at the single pacemaker cell level, isolated from the node.[12] The electrical activity in such a single cell is generated by ions flowing through discrete channels within the cell membrane. Based on experimental data, in an in silico pacemaker model consisting of roughly 6000 channels, the intrinsic dynamics of the irregular beating seen experimentally in single nodal cells was shown to be due to the (pseudo)random opening and closing of single channels.[13] In addition to synchronization between single-channel activity, pacemaker cells rely on an interstitial tissue architecture for their isolation from the rest of the myocardium and a reservoir of brown adipose tissue for ready access to energy substrates.[14–16]
Together, these three components (Fig. 1), made up of ion channels, isolated tissue and the energy reservoir form a dynamic system responsible for the generation of the heartbeat and its vulnerability to developmental challenges, such as exposure to hypoxia or infection, which in turn may determine future susceptibility to dysfunction in later life. 1) How is cardiac pacemaker synchronization (CPS) achieved?
There is no universally accepted mechanism explaining CPS. However, several important contributing components have been identified. Each channel in a pacemaker cell has one or more gates, all of which must be open to allow current to flow through the channel. Since these gates are thought to open and close in a random manner, each gate can be modeled by a Markov process, assigning a pseudorandom number to each gate every time that it changes state from open to closed or vice versa. This number, in conjunction with the classical voltage ‐ dependent Hodgkin Huxley ‐ like rate constants that control the speed with which a gate will open or close, then determines when that gate will next change state. In the study by Guevara and Lewis [13], the model’s behavior was consistent with the hypothesis that the irregular beating seen experimentally in single nodal cells is due to the (pseudo)random opening and closing of single channels. However, since the pseudorandom number generator used in the simulations is deterministic, one cannot state that the activity in the model is random (or stochastic). If one or more electric shocks are delivered directly to the heart near the intrinsic pacemaker, the heart rhythm is modified transiently but re-establishes itself with the same frequency as before within a few seconds.[17] Such behavior represents a memory of the pattern of the intrinsic CPS, an imprint of the physiological relationship between the pacemaker cells and their dynamic properties, at the channel, cellular and multicellular levels. In addition to the complex temporal dynamics of CPS based on the channel physiology and self-organization, there is an emerging role for compartmentalized (i.e., associated with distinct, spatially-confined microdomains) organization of pacemaker signaling complexes in the sinoatrial node pacemaker cells.[6] Upon structural remodeling of the sinoatrial node, disruption in subcellular targeting of pacemaker proteins and associated signaling molecules may affect their biophysical properties, neurohormonal regulation as well as protein–protein interactions within the pacemaker signaling complex. This will disturb rhythmic generation of action potentials, thereby contributing to the pathophysiology of the sinus node function. Evidence for these relationships is clear, derived from patients and animal models with genetic defects of scaffolding proteins, which are closely associated with sinus node disease. For instance, via indirect changes of key components in the coupled-clock systems in terms of protein expression, function and membrane localization. This extends beyond the classical concept of electrical remodeling, according to which dysfunction can be explained by straightforward increases or decreases in protein expression alone, and adds a new dimension to cardiovascular disease. It thus introduces a novel framework for therapeutic approaches for the treatment of pacemaker dysfunction, targeted at preventing the degradation of cardiac cytoarchitecture or abnormalities in the compartmentalization of cardiac pacemaker microdomains, incurred as a result of the intrauterine insults.[18] The single-cell resolution in situ study of cardiac conduction system has recently yielded an unprecedented and detailed transcriptional biomarker landscape of the cardiac pacemaker cells that includes ion channels (Hcn4, Hcn1, hyperpolarization-activated cyclic nucleotide-gated channels 4 and 1), transcription factors (Isl1, Tbx3, Tbx18, Shox2), gap junction gene Gjc1, Igfbp5 (insulin-like growth factor–binding protein 5), Smoc2 (SPARC-related modular calcium-binding protein), Ntm (Neurotrimin), Cpne5 (Copine 5) and Rgs6 (regulator of G-protein signaling type 6).[6,19] These represent targets for future studies of CPS development in utero in healthy and complicated pregnancy (Fig. 1). The automaticity of action potentials within the sinoatrial node’s pacemaker cells is driven by a dynamic, bidirectionally coupled interaction of two clock systems. These are represented by the spontaneous rhythmic local Ca releases generated by a Ca clock and the activity of the electrogenic surface membrane molecules, especially the ion channels (the M clock). The coupling between the intracellular Ca and the fluctuations of the surface membrane potential is mediated by the Ca –cAMP–protein kinase A (PKA) signaling. Consequently, beta-adrenergic receptor stimulation can restore or further enhance this coupling and the CPS, as it increases the cAMP concentrations.[10] The study highlights how dependent CPS is on energy availability and complements in this sense the work by Gu and colleagues[20], discussed further below. Future developmental studies of the expression profiles of these biomarkers should yield insights into how physiological and pathophysiological stimuli during gestation impact the maturation and synchronization of the cardiac pacemaker cells.
2) How is CPS impacted by developmental hypoxia and infection?
CPS memory goes beyond re-establishing its intrinsic rhythm. Cardiac disease can also impair CPS leaving its “memory imprint” on it.[9] That is, impaired intrinsic properties of pacemaker cells become manifest in an altered heart rate and HRV in the context of heart disease.[9] Ischemia slows pacemaker activity in rabbit sinoatrial node pacemaker cells because two inward currents, I
Ca,T and a presumed I
NCX , are diminished.[21] Infection, in particular, gram-negative bacterial infection impacts CPS and hence HRV via direct effects on pacemaker channels resulting in ionic current remodeling.[22,23] This was studied in the human hyperpolarization-activated cyclic nucleotide-gated channel 2 (hHCN2), which is stably expressed in HEK293 and cardiac pacemaker cells, exposed to lipopolysaccharides (LPS). The ‘funny’ current ( I f) of pacemaker cells is the main membrane mechanism generating the diastolic depolarization; the phase of the action potential responsible for cardiac spontaneous activity. I f is also responsible for the modulation of the heart rate by the autonomic nervous system (ANS).[24] LPS released from Gram ‐ negative bacteria during severe sepsis reduces both native I f and HCN2 ‐ or HCN4 ‐ mediated currents and shifts their activation to more negative potentials.[25] The LPS effect is not mediated by any of the intracellular modulatory pathways affecting HCN channels, but is instead due to a direct interaction with the channel.[22] Klöckner et al . have shown that alterations of the HCN current require the integrity of the LPS molecule, since neither the pro ‐ inflammatory lipid A alone nor the O ‐ chain polysaccharidic region alone were effective. LPS release during sepsis can modify the I f properties by shifting them out of their physiological range, which may impair normal rate control and decrease the cardiac safety factor, making the system more prone to life ‐ threatening perturbations (Fig. 1).[25] In a recent study, the vulnerability of sinus node cells to lipotoxicity has been highlighted.[26] Future studies should investigate how this vulnerability relates to fetal exposures to malnutrition, maternal obesity, chronic hypoxia and IUGR. Energy availability plays a role in CPS via modulation of the mitochondrial dynamics.[20] While the data in induced cardiac pacemaker cells should be considered with caution, as it compares to native cardiac pacemaker cells, this approach represents one of the rare opportunities to study this rare subpopulation of human cardiac myocytes with rhythm generating properties. The isolated, native sinoatrial node cells exhibited a remarkable resilience to prolonged hypoxia followed by reoxygenation [27–29], a condition that would irreversibly damage chamber cardiomyocytes. This is due to a significantly lower global metabolic demand in the cardiac pacemaker cells compared to cardiomyocytes. It remains to be investigated whether and from what time point in gestation the fetal cardiac pacemaker cells exhibit similar ischemia tolerance and what consequences such ischemic events have on programming CPS properties and susceptibility to future insults, e.g., with increased risk for arrhythmia development. Overall, ischemia and infection modulate CPS via altering the intrinsic properties of the electrical activity of pacemaker cells. This may lead to changes in HRV and arrhythmia.[30] Three major components with varying degrees of vulnerabilities to adverse intrauterine conditions emerge, namely HCN2/HCN4 channels, brown adipose tissue and the intercellular matrix (Fig. 1). It is likely that interference by LPS, ischemia, hypoxia or fibrosis with any of these constituents of the cardiac pacemaker cell environment, would yield delays in CPS maturation and function.
3) CPS and fetal development?
What are the most relevant currents to focus on that contribute to CPS and iHRV?[26,31] The literature suggests that these are the If , mediated in particular via HCN4. How do these behave during fetal and early postnatal development? In the rabbit, the sodium current (I(Na)) contributes to the sinoatrial node automaticity in the newborn but not in the adult period, when heart rate is slower. In contrast, heart rate is high and the I(Na) is functional in the adult mouse sinoatrial node. [32] Canine sinoatrial node cells have TTX-sensitive I(Na), which decreases with postnatal age. The current does not contribute to normal automaticity in isolated adult cells but can be recruited to sustain excitability if nodal cells are hyperpolarized. This is particularly relevant in the fetal and newborn periods, when the I(Na) is large and the autonomic balance favors vagal tone. The pacemaker current, I f, is expressed in ventricular myocytes from neonatal rats and progressively disappears with advancing age. When present, it shows electrophysiological properties similar to I f re-expressed in hypertrophied adult rat ventricular myocytes. Thus, it is likely that the occurrence of I f in ventricular myocytes of hypertrophied and failing hearts is due to the re-expression of a fetal gene.[33] In the heart, both during development and in adulthood, the kinetic properties and expression levels of I f match the need of the cells to perform pacemaker activity. While at immature stages of embryonic development all cardiomyocytes are autorhythmic and express a robust I f current activated within the physiological range of voltages, shortly after birth the pacemaker activity and I f expression become restricted to cells of the conduction system only, while in the working myocardium the current is downregulated and its activation is shifted to more negative, non ‐ physiological voltages.[24] The characteristics of native I f are the results of a complex interplay between different factors. The kinetic properties and expression levels of the f ‐ channels, encoded by the HCN genes (HCN1–4), are finely regulated by auxiliary proteins (caveolin ‐
3, MiRP ‐
1, KCR ‐
1, SAP97, TRP8) and lipids (phosphatidylinositol 4,5 ‐ bisphosphate) as well as by post ‐ translational modifications, such as phosphorylation and glycosylation. Any alteration of these regulatory pathways may modify the contribution of I f to the action potential, leading to an arrhythmogenic phenotype.[24] It is not clear when and how gradually during gestation the CPS develops, so that we can speak of a “healthy” intrinsic pacemaker rhythm. These unknowns, as well as the genetic program of pacemaker cell development, are under investigation.[3,34,35] Our recent data suggest that fetal intrinsic beat-to-beat HRV is present in the healthy fetus near-term and shows a memory of previous exposure to chronic hypoxia.[36] The two-clock coupling mechanism discussed in section 1 [10] may provide at least a partial explanation for this memory to chronic hypoxic conditions which, in addition to lowering energy availability, is known to program sympathetic hyperactivity in the offspring.[36–38]
4) What is iHRV?
From studies of healthy adult subjects during exercise and studies of heart-transplant patients, the field is aware that intrinsic components of cardiac rhythm can contribute substantially to HRV.[11,39,40] The iHRV is thought to result from a three-scale level system of synchronization occurring at the ion channel, cell and multi-cell layers of organization.[41] At each of these levels, various physiological and pathophysiological stimuli shape, modulate and program the synchronization properties and the development CPS in ways we are only just beginning to unravel. These intrinsic fluctuations at multiple time scales combined with a modulatory input from the ANS explain, at least in part, the fractal and complexity properties of iHRV, respectively.[11,23,40] With regard to modulation by the ANS, the focus so far has been on sympathetic beta-adrenergic signaling. Future work should also explore the modulation by vagal components. We have recently discovered that iHRV originates in fetal life and that chronic fetal hypoxia significantly alters it.[36] Interestingly, iHRV and its memory of chronic fetal hypoxia are characterized by nonlinear properties reflecting its recurrent state and chaotic behavior; properties inherent of fractal dynamics.[42] The significant relationship between nonlinear measures of fetal iHRV and left ventricular end-diastolic pressure (LVEDP), which is elevated in fetuses from hypoxic pregnancy, suggests that such iHRV indices could reflect fetal myocardial dysfunction, particularly during cardiac diastole. Therefore, alterations in fetal iHRV may prove clinically useful as a biomarker of impaired cardiac reserve and fetal myocardial decompensation during antepartum or intrapartum monitoring of the fetal heart in a high-risk pregnancy or complicated labor. What are the mechanisms contributing to fetal iHRV in late gestation in a normal pregnancy? The above review of the mechanisms underlying the CPS and influences on its activity suggests a number of factors capable of impacting the ion channel activity directly, such as LPS or hypoxia, via energy availability or via the alterations in the compartmentalization architecture of the sinoatrial node pacemaker cells. We describe in Section 1 twelve genes involved in the programmed differentiation of the cardiac pacemaker cells throughout embryonic development. The disruption of this complex spatiotemporal genetic program is likely to impact the CPS and hence leave a pathophysiological imprint on the susceptibility of the cardiac pacemaker cells to future insults, in particular cardiac arrhythmia or arrest. The study of the two-clock coupling-based generation of CPS underscores the intricate interplay between the cardiac pacemaker cells, neural autonomic input and energy availability as another layer of the putative transfer mechanisms of in utero hypoxia upon iHRV.[10] This likely dependence upon sympathetic influences links fetal adaptation to chronic stressors with the autochthonic responses within the cardiac pacemaker cell network, emerging from the adaptive processes within the excitatory cells themselves in response to chronic hypoxia. Therefore, there are several putative transfer mechanisms by which in utero chronic hypoxia may imprint upon iHRV. A systematic exploration of such candidate pathways would be a rich future avenue of research (summarized in Fig. 1). 0
5) Can the impact of fetal hypoxia or infection on CPS be captured by HRV monitoring?
In the above sections, we discuss the contributions to CPS, iHRV and HRV of the complex, multi-layered activity of the cardiac pacemaker cells. This opens a multitude of questions regarding the clinical application, which we list below in the hope that they will inspire future investigations: a) Are there suitable HRV properties likely to capture iHRV in vivo? In [36] we report such putative HRV measures representing recurrent states, chaotic and fractal dynamics. We invite the interested reader to explore this study and test the identified measurable outcomes in their own investigations; b) Could we distinguish between the effect of chronic fetal hypoxia versus inflammation on CPS from HRV analysis? Our findings ex vivo indicate an impact of chronic hypoxia on CPS and iHRV and studies in vivo identified HRV signatures capable of tracking fetal systemic and gut- and brain-specific inflammatory responses over a period of days [43–45]. However, it is yet not clear to what extent the HRV measures identified in these studies reflect contributions from iHRV. The mechanisms discussed in Section 2 support this possibility. To date, no data exist on the specificity of HRV outcomes to selective challenges, such as cardiac inflammation. This remains the subject of future studies; c) What is the relationship between CPS and myocardial development? In a recent study, we found a significant correlation between fetal iHRV and two measures of fetal cardiac diastolic dysfunction, namely LVEDP and the minimum rate of change of ventricular pressure (dP/dt) in fetuses from pregnancy affected by chronic hypoxia.[36] Could this reflect a dual pathological impact of hypoxia on both cardiomyocyte and pacemaker cell development? Alternatively, could these relationships reflect a functional physiological adaptive response relationship between these cell populations? These questions remain open to debate; d) The HRV code has been proposed as an overarching concept incorporating the various multi-scale contributions of inter-organ communication reflected in HRV.[8,46,47] We suggest that an understanding of the impact on CPS of chronic fetal hypoxia or infection during pregnancy should be sought within the integrative framework of the HRV code; e) Could we use mathematical modeling to derive predictions for HRV properties, likely to work as iHRV biomarkers in the human clinical setting? Mathematical models exist for CPS and fetal cardiovascular responses to labor, for example.[48,49] Such models could be combined and explored in silico to derive numerical predictions of HRV properties specific to fetal iHRV and their quantitative contribution to fetal cardiac function in vivo. 1
6) Conclusions and outlook
The three major components of the cardiac pacemaker system, comprised of the cellular ion channels, isolated tissue and the energy reservoir, form a dynamic network responsible for the generation of iHRV (Fig. 1). Future studies should determine the spatiotemporal vulnerability of this system, i.e., time periods during embryonic and fetal development that are particularly sensitive and susceptible to adverse intrauterine insults, such as chronic hypoxia or infection. This developmental vulnerability, in turn, may determine future susceptibility to cardiac dysfunction or limited tolerance in later life. Specifically, the following molecular candidate targets could provide interesting studies of CPS development in utero and postnatally: ion channels (Hcn4, Hcn1), transcription factors (Isl1, Tbx3, Tbx18, Shox2), gap junction gene Gjc1, Igfbp5, Smoc2, Ntm, Cpne5 and Rgs6. Functional studies in pacemaker cells should explore the developmental profile of CPS dynamics in vitro or using the Langendorff preparation ex vivo with attention to the vulnerability of the CPS to noxious stimuli such as ischemia, hypoxia or inflammation. To gauge the impact of energy availability during intrauterine development on CPS maturation, the effects of alterations in nutrition on CPS development should be investigated. Such studies could involve maternal malnutrition or maternal obesity during pregnancy with or without chronic hypoxia and/or infection to assess partial and combined effects. Such studies would also permit the spatial and qualitative characterization of changes to the isolating tissue, such as fibrosis, and the micro-domain organization within the cardiac pacemaker cells. Acknowledgements
MGF has been funded by the Canadian Institutes of Health Research (CIHR). DAG has been funded by the British Heart Foundation (BHF). 3
References
1. McMillen, I.C.; Robinson, J.S. Developmental origins of the metabolic syndrome: prediction, plasticity, and programming.
Physiol. Rev. , , 571–633. 2. Giussani, D.A.; Davidge, S.T. Developmental programming of cardiovascular disease by prenatal hypoxia. J. Dev. Orig. Health Dis. , , 328–337. 3. Vedantham, V. New Approaches to Biological Pacemakers: Links to Sinoatrial Node Development. Trends Mol. Med. , , 749–761. 4. Protze, S.I.; Liu, J.; Nussinovitch, U.; Ohana, L.; Backx, P.H.; Gepstein, L.; Keller, G.M. Sinoatrial node cardiomyocytes derived from human pluripotent cells function as a biological pacemaker. Nat. Biotechnol. , , 56–68. 5. van den Heuvel, N.H.L.; van Veen, T.A.B.; Lim, B.; Jonsson, M.K.B. Lessons from the heart: mirroring electrophysiological characteristics during cardiac development to in vitro differentiation of stem cell derived cardiomyocytes. J. Mol. Cell. Cardiol. , , 12–25. 6. Goodyer, W.R.; Beyersdorf, B.M.; Paik, D.T.; Tian, L.; Li, G.; Buikema, J.W.; Chirikian, O.; Choi, S.; Venkatraman, S.; Adams, E.L.; et al. Transcriptomic Profiling of the Developing Cardiac Conduction System at Single-Cell Resolution. Circ. Res. , , 379–397. 7. Georgieva, A.; Abry, P.; Chudáček, V.; Djurić, P.M.; Frasch, M.G.; Kok, R.; Lear, C.A.; Lemmens, S.N.; Nunes, I.; Papageorghiou, A.T.; et al. Computer-based intrapartum fetal monitoring and beyond: a review of the 2nd Workshop on Signal Processing and Monitoring in Labor (Oct 2017, Oxford UK). Acta Obstet. Gynecol. Scand. . 8. Frasch, M.G. Saving the brain one heartbeat at a time.
J. Physiol. , . 9. Yaniv, Y.; Lyashkov, A.E.; Lakatta, E.G. Impaired signaling intrinsic to sinoatrial node pacemaker cells affects heart rate variability during cardiac disease. J Clin Trials , . 10. Tsutsui, K.; Monfredi, O.J.; Sirenko-Tagirova, S.G.; Maltseva, L.A.; Bychkov, R.; Kim, M.S.; Ziman, B.D.; Tarasov, K.V.; Tarasova, Y.S.; Zhang, J.; et al. A coupled-clock system drives the automaticity of human sinoatrial nodal pacemaker cells. Sci. Signal. , . 11. Yaniv, Y.; Lyashkov, A.E.; Lakatta, E.G. The fractal-like complexity of heart rate variability beyond neurotransmitters and autonomic receptors: signaling intrinsic to sinoatrial node pacemaker cells. Cardiovasc Pharm Open Access , . 12. Wilders, R.; Jongsma, H.J. Beating irregularity of single pacemaker cells isolated from the rabbit sinoatrial node. Biophys. J. , , 2601–2613. 13. Guevara, M.R.; Lewis, T.J. A minimal single-channel model for the regularity of beating in the sinoatrial node. Chaos , , 174–183. 14. Bressan, M.; Henley, T.; Louie, J.D.; Liu, G.; Christodoulou, D.; Bai, X.; Taylor, J.; Seidman, C.E.; Seidman, J.G.; Mikawa, T. Dynamic Cellular Integration Drives Functional Assembly of the Heart’s Pacemaker Complex. Cell Rep. , , 2283–2291. 15. Unudurthi, S.D.; Wolf, R.M.; Hund, T.J. Role of sinoatrial node architecture in maintaining a balanced source-sink relationship and synchronous cardiac pacemaking. Front. Physiol. , , 446. 16. Fahrenbach, J.P.; Mejia-Alvarez, R.; Banach, K. The relevance of non-excitable cells for cardiac pacemaker function. J. Physiol. , , 565–578. 17. Glass, L. Synchronization and rhythmic processes in physiology. Nature , , 277–284. 18. Lang, D.; Glukhov, A.V. Functional Microdomains in Heart’s Pacemaker: A Step Beyond Classical Electrophysiology and Remodeling. Frontiers in Physiology . 19. Hoogaars, W.M.H.; Engel, A.; Brons, J.F.; Verkerk, A.O.; de Lange, F.J.; Wong, L.Y.E.; Bakker, M.L.; Clout, D.E.; Wakker, V.; Barnett, P.; et al. Tbx3 controls the sinoatrial node gene program and imposes pacemaker function on the atria. Genes Dev. , , 1098–1112. 4 20. Gu, J.-M.; Grijalva, S.I.; Fernandez, N.; Kim, E.; Foster, D.B.; Cho, H.C. Induced cardiac pacemaker cells survive metabolic stress owing to their low metabolic demand. Exp. Mol. Med. , , 105. 21. Du, Y.-M.; Nathan, R.D. Ionic basis of ischemia-induced bradycardia in the rabbit sinoatrial node. J. Mol. Cell. Cardiol. , , 315–325. 22. Klockner, U.; Rueckschloss, U.; Grossmann, C.; Matzat, S.; Schumann, K.; Ebelt, H.; Muller-Werdan, U.; Loppnow, H.; Werdan, K.; Gekle, M. Inhibition of cardiac pacemaker channel hHCN2 depends on intercalation of lipopolysaccharide into channel-containing membrane microdomains. J. Physiol. , , 1199–1211. 23. Papaioannou, V.E.; Verkerk, A.O.; Amin, A.S.; de Bakker, J.M. Intracardiac origin of heart rate variability, pacemaker funny current and their possible association with critical illness. Curr. Cardiol. Rev. , , 82–96. 24. DiFrancesco, D. The role of the funny current in pacemaker activity. Circ. Res. , , 434–446. 25. Barbuti, A.; DiFrancesco, D. The “funny” side of sepsis. J. Physiol. , , 1171. 26. Linscheid, N.; Logantha, S.J.R.J.; Poulsen, P.C.; Zhang, S.; Schrölkamp, M.; Egerod, K.L.; Thompson, J.J.; Kitmitto, A.; Galli, G.; Humphries, M.J.; et al. Quantitative proteomics and single-nucleus transcriptomics of the sinus node elucidates the foundation of cardiac pacemaking. Nat. Commun. , , 2889. 27. Kohlhardt, M.; Mnich, Z.; Maier, G. Alterations of the excitation process of the sinoatrial pacemaker cell in the presence of anoxia and metabolic inhibitors. J. Mol. Cell. Cardiol. , , 477–488. 28. Senges, J.; Mizutani, T.; Pelzer, D.; Brachmann, J.; Sonnhof, U.; Kübler, W. Effect of hypoxia on the sinoatrial node, atrium, and atrioventricular node in the rabbit heart. Circ. Res. , , 856–863. 29. Nishi, K.; Yoshikawa, Y.; Sugahara, K.; Morioka, T. Changes in electrical activity and ultrastructure of sinoatrial nodal cells of the rabbit’s heart exposed to hypoxic solution. Circ. Res. , , 201–213. 30. Yaniv, Y.; Tsutsui, K.; Lakatta, E.G. Potential effects of intrinsic heart pacemaker cell mechanisms on dysrhythmic cardiac action potential firing. Front. Physiol. , , 47. 31. Lakatta, E.G.; DiFrancesco, D. What keeps us ticking: a funny current, a calcium clock, or both? J. Mol. Cell. Cardiol. , , 157–170. 32. Protas, L.; Oren, R.V.; Clancy, C.E.; Robinson, R.B. Age-dependent changes in Na current magnitude and TTX-sensitivity in the canine sinoatrial node. J. Mol. Cell. Cardiol. , , 172–180. 33. Cerbai, E.; Pino, R.; Sartiani, L.; Mugelli, A. Influence of postnatal-development on I(f) occurrence and properties in neonatal rat ventricular myocytes. Cardiovasc. Res. , , 416–423. 34. Hu, W.; Xin, Y.; Zhao, Y.; Hu, J. Shox2: The Role in Differentiation and Development of Cardiac Conduction System. Tohoku J. Exp. Med. , , 177–186. 35. Vedantham, V.; Galang, G.; Evangelista, M.; Deo, R.C.; Srivastava, D. RNA sequencing of mouse sinoatrial node reveals an upstream regulatory role for Islet-1 in cardiac pacemaker cells. Circ. Res. , , 797–803. 36. Frasch, M.G.; Herry, C.L.; Niu, Y.; Giussani, D.A. First evidence that intrinsic fetal heart rate variability exists and is affected by hypoxic pregnancy. J. Physiol. . 37. Giussani, D.A.; Camm, E.J.; Niu, Y.; Richter, H.G.; Blanco, C.E.; Gottschalk, R.; Blake, E.Z.; Horder, K.A.; Thakor, A.S.; Hansell, J.A.; et al. Developmental programming of cardiovascular dysfunction by prenatal hypoxia and oxidative stress.
PLoS One , , e31017. 38. Rouwet, E.V.; Tintu, A.N.; Schellings, M.W.M.; van Bilsen, M.; Lutgens, E.; Hofstra, L.; Slaaf, D.W.; Ramsay, G.; Le Noble, F.A.C. Hypoxia induces aortic hypertrophic growth, left 5 ventricular dysfunction, and sympathetic hyperinnervation of peripheral arteries in the chick embryo. Circulation , , 2791–2796. 39. Bernardi, L.; Salvucci, F.; Suardi, R.; Solda, P.L.; Calciati, A.; Perlini, S.; Falcone, C.; Ricciardi, L. Evidence for an intrinsic mechanism regulating heart rate variability in the transplanted and the intact heart during submaximal dynamic exercise? Cardiovasc. Res. , , 969–981. 40. Yaniv, Y.; Ahmet, I.; Liu, J.; Lyashkov, A.E.; Guiriba, T.R.; Okamoto, Y.; Ziman, B.D.; Lakatta, E.G. Synchronization of sinoatrial node pacemaker cell clocks and its autonomic modulation impart complexity to heart beating intervals. Heart Rhythm , , 1210–1219. 41. Binah, O.; Weissman, A.; Itskovitz-Eldor, J.; Rosen, M.R. Integrating beat rate variability: from single cells to hearts. Heart Rhythm , , 928–932. 42. Webber, C.L., Jr Recurrence quantification of fractal structures. Front. Physiol. , , 382. 43. Durosier, L.D.; Herry, C.L.; Cortes, M.; Cao, M.; Burns, P.; Desrochers, A.; Fecteau, G.; Seely, A.J.E.; Frasch, M.G. Does heart rate variability reflect the systemic inflammatory response in a fetal sheep model of lipopolysaccharide-induced sepsis? Physiol. Meas. , , 2089–2102. 44. Herry, C.L.; Cortes, M.; Wu, H.-T.; Durosier, L.D.; Cao, M.; Burns, P.; Desrochers, A.; Fecteau, G.; Seely, A.J.E.; Frasch, M.G. Temporal Patterns in Sheep Fetal Heart Rate Variability Correlate to Systemic Cytokine Inflammatory Response: A Methodological Exploration of Monitoring Potential Using Complex Signals Bioinformatics. PLoS One , , e0153515. 45. Liu, H.L.; Garzoni, L.; Herry, C.; Durosier, L.D.; Cao, M.; Burns, P.; Fecteau, G.; Desrochers, A.; Patey, N.; Seely, A.J.E.; et al. Can Monitoring Fetal Intestinal Inflammation Using Heart Rate Variability Analysis Signal Incipient Necrotizing Enterocolitis of the Neonate? Pediatr. Crit. Care Med. , , e165–76. 46. Herry, C.L.; Burns, P.; Desrochers, A.; Fecteau, G.; Durosier, L.D.; Cao, M.; Seely, A.J.E.; Frasch, M.G. Vagal contributions to fetal heart rate variability: an omics approach. Physiol. Meas. , , 065004. 47. Frasch, M.G. Heart rate variability code: Does it exist and can we hack it? arXiv [q-bio.TO] Am. J. Physiol. Heart Circ. Physiol. , , H397–421. 49. Wang, Q.; Gold, N.; Frasch, M.G.; Huang, H.; Thiriet, M.; Wang, X. Mathematical Model of Cardiovascular and Metabolic Responses to Umbilical Cord Occlusions in Fetal Sheep. Bull. Math. Biol. , , 2264–2293. 6 Figure legends