Persistent Glossopharyngeal Nerve Respiratory Discharge Patterns after Ponto-Medullary Transection

Shape and size of the nasopharyngeal airway is controlled by muscles innervated facial, glossopharyngeal, vagal, and hypoglossal cranial nerves. Contrary to brainstem networks that drive facial, vagal and hypoglossal nerve activities (FNA, VNA, HNA) the discharge patterns and origins of glossopharyngeal nerve activity (GPNA) remain poorly investigated. Here, an in situ perfused brainstem preparation (n=19) was used for recordings of GPNA in relation to phrenic (PNA), FNA, VNA and HNA. Brainstem transections were performed (n=10/19) to explore the role of pontomedullary synaptic interactions in generating GPNA. GPNA generally mirrors FNA and HNA discharge patterns and displays pre-inspiratory activity relative to the PNA, followed by robust inspiratory discharge in coincidence with PNA. Postinspiratory (early expiratory) discharge was, contrary to VNA, generally absent in FNA, GPNA or HNA. As described previously FNA and HNA discharge was virtually eliminated after pontomedullary transection while an apneustic inspiratory motor discharge was maintained in PNA, VNA and GPNA. After brainstem transection GPNA displayed an increased tonic activity starting during mid-expiration and thus developed prolonged pre-inspiratory activity compared to control. In conclusion respiratory GPNA reflects FNA and HNA which implies similar function in controlling upper airway patency during breathing. That GPNA preserved its pre-inspiratory/inspiratory discharge pattern in relation PNA after pontomedullary transection suggest that GPNA premotor circuits may have a different anatomical distribution compared HNA and FNA and thus may therefore hold a unique role in in preserving airway patency.

Peritoneal sepsis caused by Escherichia coli triggers brainstem inflammation and alters the function of sympatho-respiratory control circuits

BACKGROUND

Sepsis has a high mortality rate due to multiple organ failure. However, the influence of peripheral inflammation on brainstem autonomic and respiratory circuits in sepsis is poorly understood. Our working hypothesis is that peripheral inflammation affects central autonomic circuits and consequently contributes to multiorgan failure in sepsis.

METHODS

In an Escherichia coli (E. coli)-fibrin clot model of peritonitis, we first recorded ventilatory patterns using plethysmography before and 24 h after fibrin clot implantation. To assess whether peritonitis was associated with brainstem neuro-inflammation, we measured cytokine and chemokine levels in Luminex assays. To determine the effect of E. coli peritonitis on brainstem function, we assessed sympatho-respiratory nerve activities at baseline and during brief (20 s) hypoxemic ischemia challenges using in situ-perfused brainstem preparations (PBPs) from sham or infected rats. PBPs lack peripheral organs and blood, but generate vascular tone and in vivo rhythmic activities in thoracic sympathetic (tSNA), phrenic and vagal nerves.

RESULTS

Respiratory frequency was greater (p < 0.001) at 24 h post-infection with E. coli than in the sham control. However, breath-by-breath variability and total protein in the BALF did not differ. IL-1β (p < 0.05), IL-6 (p < 0.05) and IL-17 (p < 0.04) concentrations were greater in the brainstem of infected rats. In the PBP, integrated tSNA (p < 0.05) and perfusion pressure were greater (p < 0.001), indicating a neural-mediated pathophysiological high sympathetic drive. Moreover, respiratory frequency was greater (p < 0.001) in PBPs from infected rats than from sham rats. Normalized phase durations of inspiration and expiration were greater (p < 0.009, p < 0.015, respectively), but the post-inspiratory phase (p < 0.007) and the breath-by-breath variability (p < 0.001) were less compared to sham PBPs. Hypoxemic ischemia triggered a biphasic response, respiratory augmentation followed by depression. PBPs from infected rats had weaker respiratory augmentation (p < 0.001) and depression (p < 0.001) than PBPs from sham rats. In contrast, tSNA in E. coli-treated PBPs was enhanced throughout the entire response to hypoxemic ischemia (p < 0.01), consistent with sympathetic hyperactivity.

CONCLUSION

We show that peripheral sepsis caused brainstem inflammation and impaired sympatho-respiratory motor control in a single day after infection. We conclude that central sympathetic hyperactivity may impact vital organ systems in sepsis.

Dynamics of ventilatory pattern variability and Cardioventilatory Coupling during systemic inflammation in rats

Introduction: Biometrics of common physiologic signals can reflect health status. We have developed analytics to measure the predictability of ventilatory pattern variability (VPV, Nonlinear Complexity Index (NLCI) that quantifies the predictability of a continuous waveform associated with inhalation and exhalation) and the cardioventilatory coupling (CVC, the tendency of the last heartbeat in expiration to occur at preferred latency before the next inspiration). We hypothesized that measures of VPV and CVC are sensitive to the development of endotoxemia, which evoke neuroinflammation. Methods: We implanted Sprague Dawley male rats with BP transducers to monitor arterial blood pressure (BP) and recorded ventilatory waveforms and BP simultaneously using whole-body plethysmography in conjunction with BP transducer receivers. After baseline (BSLN) recordings, we injected lipopolysaccharide (LPS, n = 8) or phosphate buffered saline (PBS, n =3) intraperitoneally on 3 consecutive days. We recorded for 4-6 h after the injection, chose 3 epochs from each hour and analyzed VPV and CVC as well as heart rate variability (HRV). Results: First, the responses to sepsis varied across rats, but within rats the repeated measures of NLCI, CVC, as well as respiratory frequency (fR), HR, BP and HRV had a low coefficient of variation, (<0.2) at each time point. Second, HR, fR, and NLCI increased from BSLN on Days 1-3; whereas CVC decreased on Days 2 and 3. In contrast, changes in BP and the relative low-(LF) and high-frequency (HF) of HRV were not significant. The coefficient of variation decreased from BSLN to Day 3, except for CVC. Interestingly, NLCI increased before fR in LPS-treated rats. Finally, we histologically confirmed lung injury, systemic inflammation via ELISA and the presence of the proinflammatory cytokine, IL-1β, with immunohistochemistry in the ponto-medullary respiratory nuclei. Discussion: Our findings support that NLCI reflects changes in the rat's health induced by systemic injection of LPS and reflected in increases in HR and fR. CVC decreased over the course to the experiment. We conclude that NLCI reflected the increase in predictability of the ventilatory waveform and (together with our previous work) may reflect action of inflammatory cytokines on the network generating respiration.

Cannabinoid Receptor mRNA Expression in Central and Peripheral Tissues in a Rodent Model of Peritonitis

Introduction: Our laboratory investigates changes in the respiratory pattern during systemic inflammation in various rodent models. The endogenous cannabinoid system (ECS) regulates cytokine production and mitigates inflammation. Inflammation not only affects cannabinoid (CB) 1 and CB2 receptor gene expression (Cnr1 and Cnr2), but also increases the predictability of the ventilatory pattern. Objectives: Our primary objective was to track ventilatory pattern variability and transcription of Cnr1 and Cnr2 mRNA, and of Il1b, Il6, and tumor necrosis factor-alpha (Tnfa) mRNAs at multiple time points in central and peripheral tissues during systemic inflammation induced by peritonitis. Methods: In male Sprague Dawley rats (n=24), we caused peritonitis by implanting a fibrin clot containing either 0 or 25×106 Escherichia coli intraperitoneally. We recorded breathing with whole-animal plethysmography at baseline and 1 h before euthanasia. We euthanized the rats at 3, 6, or 12 h after inoculation and harvested the pons, medulla, lung, and heart for gene expression analysis. Results: With peritonitis, Cnr1 mRNA more than Cnr2 mRNA was correlated to Il1b, Il6, and Tnfa mRNAs in medulla, pons, and lung and changed oppositely in the pons, medulla, and lung. These changes were associated with increased predictability of ventilatory pattern. Specifically, nonlinear complexity index correlated with increased Cnr1 mRNA in the pons and medulla, and coefficient of variation for cycle duration correlated with Cnr1 and Cnr2 mRNAs in the lung. Conclusion: The mRNAs for ECS receptors varied with time during the central and peripheral inflammatory response to peritonitis. These changes occurred in the brainstem, which contains the network that generates breathing pattern and thus, may participate in ventilatory pattern changes during systemic inflammation.

Heartbeats entrain breathing via baroreceptor-mediated modulation of expiratory activity

NEW FINDINGS

  Cardio-ventilatory coupling refers to the onset of inspiration occurring at a preferential latency following the last heartbeat (HB) in expiration. According to the cardiac-trigger hypothesis, the pulse pressure initiates an inspiration via baroreceptor activation. However, the central neural substrate mediating this coupling remains undefined. Using a combination of animal data, human data and mathematical modelling, this study tests the hypothesis that the HB, by way of pulsatile baroreflex activation, controls the initiation of inspiration that occurs through a rapid neural activation loop from the carotid baroreceptors to Bötzinger complex expiratory neurons.

ABSTRACT

Cardio-ventilatory coupling refers to a heartbeat (HB) occurring at a preferred latency prior to the next breath. We hypothesized that the pressure pulse generated by a HB activates baroreceptors that modulate brainstem expiratory neuronal activity and delay the initiation of inspiration. In supine male subjects, we recorded ventilation, electrocardiogram and blood pressure during 20-min epochs of baseline, slow-deep breathing and recovery. In in situ rodent preparations, we recorded brainstem activity in response to pulses of perfusion pressure. We applied a well-established respiratory network model to interpret these data. In humans, the latency between a HB and onset of inspiration was consistent across different breathing patterns. In in situ preparations, a transient pressure pulse during expiration activated a subpopulation of expiratory neurons normally active during post-inspiration, thus delaying the next inspiration. In the model, baroreceptor input to post-inspiratory neurons accounted for the effect. These studies are consistent with baroreflex activation modulating respiration through a pauci-synaptic circuit from baroreceptors to onset of inspiration.

Traube-Hering waves are formed by interaction of respiratory sinus arrhythmia and pulse pressure modulation in healthy men

Excessive blood pressure variation is linked to the development of hypertension and other diseases. This study assesses the relative role of respiratory sinus arrhythmia (RSA) and pulse pressure (PP) on the amplitude and timing of blood pressure variability with respiration [Traube-Hering (TH) waves]. We analyzed respiratory, electrocardiogram, and blood pressure traces from healthy, supine male subjects (n = 10, mean age = 26.7 ± 1.4) during 20-min epochs of resting, slow deep breathing (SDB), and recovery. Across all epochs, blood pressure and heart rate (HR) were modulated with respiration and the magnitude of RSA; TH waves increased during SDB. The data were deconstructed using a simple mathematical model of blood pressure to dissect the relative roles of RSA and PP on TH waves. We constructed the time series of the R-wave peaks and compared the recorded TH waves with that predicted by the model. Given that cardiac output is determined by both heart rate and stroke volume, it was surprising that the magnitude of the TH waves could be captured by only HR modulation. However, RSA alone did not accurately predict the timing of TH waves relative to the respiratory cycle. Adding respiratory modulation of PP to the model corrected the phase shift showing the expected pattern of BP rising during inspiration with the peak of the TH wave during early expiration. We conclude that short-term variability of blood pressure referred to as TH waves has at least two independent mechanisms whose interaction forms their pattern: RSA and respiratory-driven changes in PP.NEW & NOTEWORTHY Variability in blood pressure has become an important metric to consider as more is learned about the link between excessive blood pressure variability and adverse health outcomes. In this study using slow deep breathing in human subjects, we found that heart rate and pulse pressure variations have comparable effects on the amplitude of blood pressure waves, and it is the common action of the two that defines the phase relationship between respiration and blood pressure oscillations.

Peripheral-to-central immune communication at the area postrema glial-barrier following bleomycin-induced sterile lung injury in adult rats

The pathways for peripheral-to-central immune communication (P → C I-comm) following sterile lung injury (SLI) are unknown. SLI evokes systemic and central inflammation, which alters central respiratory control and viscerosensory transmission in the nucleus tractus solitarii (nTS). These functional changes coincide with increased interleukin-1 beta (IL-1β) in the area postrema, a sensory circumventricular organ that connects P → C I-comm to brainstem circuits that control homeostasis. We hypothesize that IL-1β and its downstream transcriptional target, cyclooxygenase-2 (COX-2), mediate P → C I-comm in the nTS. In a rodent model of SLI induced by intratracheal bleomycin (Bleo), the sigh frequency and duration of post-sigh apnea increased in Bleo- compared to saline- treated rats one week after injury. This SLI-dependent change in respiratory control occurred concurrently with augmented IL-1β and COX-2 immunoreactivity (IR) in the funiculus separans (FS), a barrier between the AP and the brainstem. At this barrier, increases in IL-1β and COX-2 IR were confined to processes that stained for glial fibrillary acidic protein (GFAP) and that projected basolaterally to the nTS. Further, FS radial-glia did not express TNF-α or IL-6 following SLI. To test our hypothesis, we blocked central COX-1/2 activity by intracerebroventricular (ICV) infusion of Indomethacin (Ind). Continuous ICV Ind treatment prevented Bleo-dependent increases in GFAP + and IL-1β + IR, and restored characteristics of sighs that reset the rhythm. These data indicate that changes in sighs following SLI depend partially on activation of a central COX-dependent P → C I-comm via radial-glia of the FS.

Brainstem inflammation modulates the ventilatory pattern and its variability after acute lung injury in rodents

KEY POINTS

Compared with sham rats, rats a week after acute lung injury (ALI) express more pro-inflammatory cytokines in their brainstem respiratory control nuclei, exhibit a higher respiratory frequency (fR) and breathe with a more predictable pattern. These characteristics of the respiratory pattern persist in in situ preparations even after minimizing pulmonary and chemo-afferent inputs. Interleukin (IL)-1β microinjected in the nucleus tractus solitarii increases fR and the predictability of the ventilatory pattern similar to rats with ALI. Intracerebroventricular infusion of indomethacin, an anti-inflammatory drug, mitigates the effect of ALI on fR and ventilatory pattern variability. We conclude that changes in the ventilatory pattern after ALI result not only from sensory input due to pulmonary damage and dysfunction but also from neuro-inflammation.

ABSTRACT

Acute lung injury (ALI) increases respiratory rate (fR) and ventilatory pattern variability (VPV), but also evokes peripheral and central inflammation. We hypothesized that central inflammation has a role in determining the ventilatory pattern after ALI. In rat pups, we intratracheally injected either bleomycin to induce ALI or saline as a sham control. One week later, we recorded the ventilatory pattern of the rat pups using flow-through plethysmography, then formed in situ preparations from these pups and recorded their 'fictive' patterns from respiratory motor nerves. Compared with the ventilatory pattern of the sham rat pups, injured rat pups had increased fR and predictability. Surprisingly, the fictive patterns of the in situ preparations from ALI pups retained these characteristics despite removing their lungs to eliminate pulmonary sensory inputs and perfusing them with hyperoxic artificial cerebral spinal fluid to minimize peripheral chemoreceptor input. Histological processing revealed increased immunoreactivity of the pro-inflammatory cytokine Interleukin-1β (IL-1β) in the nucleus tractus solitarii (nTS) from ALI but not sham rats. In subsequent experiments, we microinjected IL-1β in the nTS bilaterally in anaesthetized naïve adult rats, which increased fR and predictability of ventilatory pattern variability (VPV) after 2 h. Finally, we infused indomethacin intracerebroventricularly during the week of survival after ALI. This did not affect sham rats, but mitigated changes in fR and VPV in ALI rats. We conclude that neuro-inflammation has an essential role in determining the ventilatory pattern of ALI rats.

Volumetric mapping of the functional neuroanatomy of the respiratory network in the perfused brainstem preparation of rats

KEY POINTS

The functional neuroanatomy of the mammalian respiratory network is far from being understood since experimental tools that measure neural activity across this brainstem-wide circuit are lacking. Here, we use silicon multi-electrode arrays to record respiratory local field potentials (rLFPs) from 196-364 electrode sites within 8-10 mm³ of brainstem tissue in single arterially perfused brainstem preparations with respect to the ongoing respiratory motor pattern of inspiration (I), post-inspiration (PI) and late-expiration (E2). rLFPs peaked specifically at the three respiratory phase transitions, E2-I, I-PI and PI-E2. We show, for the first time, that only the I-PI transition engages a brainstem-wide network, and that rLFPs during the PI-E2 transition identify a hitherto unknown role for the dorsal respiratory group. Volumetric mapping of pontomedullary rLFPs in single preparations could become a reliable tool for assessing the functional neuroanatomy of the respiratory network in health and disease.

ABSTRACT

While it is widely accepted that inspiratory rhythm generation depends on the pre-Bötzinger complex, the functional neuroanatomy of the neural circuits that generate expiration is debated. We hypothesized that the compartmental organization of the brainstem respiratory network is sufficient to generate macroscopic local field potentials (LFPs), and if so, respiratory (r) LFPs could be used to map the functional neuroanatomy of the respiratory network. We developed an approach using silicon multi-electrode arrays to record spontaneous LFPs from hundreds of electrode sites in a volume of brainstem tissue while monitoring the respiratory motor pattern on phrenic and vagal nerves in the perfused brainstem preparation. Our results revealed the expression of rLFPs across the pontomedullary brainstem. rLFPs occurred specifically at the three transitions between respiratory phases: (1) from late expiration (E2) to inspiration (I), (2) from I to post-inspiration (PI), and (3) from PI to E2. Thus, respiratory network activity was maximal at respiratory phase transitions. Spatially, the E2-I, and PI-E2 transitions were anatomically localized to the ventral and dorsal respiratory groups, respectively. In contrast, our data show, for the first time, that the generation of controlled expiration during the post-inspiratory phase engages a distributed neuronal population within ventral, dorsal and pontine network compartments. A group-wise independent component analysis demonstrated that all preparations exhibited rLFPs with a similar temporal structure and thus share a similar functional neuroanatomy. Thus, volumetric mapping of rLFPs could allow for the physiological assessment of global respiratory network organization in health and disease.

Increasing Local Excitability of Brainstem Respiratory Nuclei Reveals a Distributed Network Underlying Respiratory Motor Pattern Formation

The core circuit of the respiratory central pattern generator (rCPG) is located in the ventrolateral medulla, especially in the pre-Bötzinger complex (pre-BötC) and the neighboring Bötzinger complex (BötC). To test the hypothesis that this core circuit is embedded within an anatomically distributed pattern-generating network, we investigated whether local disinhibition of the nucleus tractus solitarius (NTS), the Kölliker-Fuse nuclei (KFn), or the midbrain periaqueductal gray area (PAG) can similarly affect the respiratory pattern compared to disinhibition of the pre-BötC/BötC core. In arterially-perfused brainstem preparations of rats, we recorded the three-phase respiratory pattern (inspiration, post-inspiration and late-expiration) from phrenic and vagal nerves before and after bilateral microinjections of the GABA(A)R antagonist bicuculline (50 nl, 10 mM). Local disinhibition of either NTS, pre-BötC/BötC, or KFn, but not PAG, triggered qualitatively similar disruptions of the respiratory pattern resulting in a highly significant increase in the variability of the respiratory cycle length, including inspiratory and expiratory phase durations. To quantitatively analyze these motor pattern perturbations, we measured the strength of phase synchronization between phrenic and vagal motor outputs. This analysis showed that local disinhibition of all brainstem target nuclei, but not the midbrain PAG, significantly decreased the strength of phase synchronization. The convergent perturbations of the respiratory pattern suggest that the rCPG expands rostrally and dorsally from the designated core but does not include higher mid-brain structures. Our data also suggest that excitation-inhibition balance of respiratory network synaptic interactions critically determines the network dynamics that underlie vital respiratory rhythm and pattern formation.

Ventilatory pattern variability as a biometric for severity of acute lung injury in rats

We hypothesize that ventilatory pattern variability (VPV) varies with the magnitude of acute lung injury (ALI). In adult male rats, we instilled a low- or high- dose of bleomycin or saline (PBS) intratracheally. While representative samples of pulmonary tissue indicated graded lung injury, coefficient of variation (CV) of TTOT did not differ among the 3 groups. Broncho-alveolar lavage fluid (BALF), respiratory rate (fR), mutual information were greater in ALI than sham rats; but did not differ between bleomycin doses. However, nonlinear complexity index (NLCI), which is the difference in sample entropy between original and surrogate data sets was greater for high- versus low- dose; but did not differ between low-dose and sham groups. Further, NLCI correlated to an injury index based on protein concentration of BALF and failure to gain weight. Finally, Receiver Operator Curves (ROCs) indicated that both mutual information and NLCI had greater sensitivity and specificity than fR and CVTTOT in identifying ALI. Thus, nonlinear analyses of VPV can distinguish ALI and out performs fR as a biometric.

Periodicity: A Characteristic of Heart Rate Variability Modified by the Type of Mechanical Ventilation After Acute Lung Injury

We present a novel approach to quantify heart rate variability (HRV) and the results of applying this approach to synthetic and original data sets. Our approach evaluates the periodicity of heart rate by calculating the transform of Relative Shannon Entropy, the maximum value of the RR interval periodogram, and the maximum, mean values, and sample entropy of the autocorrelation function. Synthetic data were generated using a Van der Pol oscillator; and the original data were electrocardiogram (ECG) recordings from anesthetized rats after acute lung injury while on biologically variable (BVV) or continuous mechanical ventilation (CMV). Analysis of the synthetic data revealed that our measures were correlated highly to the bandwidth of the oscillator and assessed periodicity. Then, applying these analytical tools to the ECGs determined that the heart rate (HR) of BVV group had less periodicity and higher variability than the HR of the CMV group. Quantifying periodicity effectively identified a readily apparent difference in HRV during BVV and CMV that was not identified by power spectral density measures during BVV and CMV. Cardiorespiratory coupling is the probable mechanism for HRV increasing during BVV and becoming periodic during CMV. Thus, the absence or presence of periodicity in ventilation determined HRV, and this mechanism is distinctly different from the cardiorespiratory uncoupling that accounts for the loss of HRV during sepsis.

Lung-injury depresses glutamatergic synaptic transmission in the nucleus tractus solitarii via discrete age-dependent mechanisms in neonatal rats

Transition periods (TPs) are brief stages in CNS development where neural circuits can exhibit heightened vulnerability to pathologic conditions such as injury or infection. This susceptibility is due in part to specialized mechanisms of synaptic plasticity, which may become activated by inflammatory mediators released under pathologic conditions. Thus, we hypothesized that the immune response to lung injury (LI) mediated synaptic changes through plasticity-like mechanisms that depended on whether LI occurred just before or after a TP. We studied the impact of LI on brainstem 2nd-order viscerosensory neurons located in the nucleus tractus solitarii (nTS) during a TP for respiratory control spanning (postnatal day (P) 11-15). We injured the lungs of Sprague-Dawley rats by intratracheal instillation of Bleomycin (or saline) just before (P9-11) or after (P17-19) the TP. A week later, we prepared horizontal slices of the medulla and recorded spontaneous and evoked excitatory postsynaptic currents (sEPSCs/eEPSCs) in vitro from neurons in the nTS that received monosynaptic glutamatergic input from the tractus solitarii (TS). In rats injured before the TP (pre-TP), neurons exhibited blunted sEPSCs and TS-eEPSCs compared to controls. The decreased TS-eEPSCs were mediated by differences in postsynaptic α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic-acid receptors (AMPAR). Specifically, compared to controls, LI rats had more Ca2+-impermeable AMPARs (CI-AMPARs) as indicated by: 1) the absence of current-rectification, 2) decreased sensitivity to polyamine, 1-Naphthyl-acetyl-spermine-trihydrochloride (NASPM) and 3) augmented immunoreactive staining for the CI-AMPAR GluA2. Thus, pre-TP-LI acts postsynaptically to blunt glutamatergic transmission. The neuroimmune response to pre-TP-LI included microglia hyper-ramification throughout the nTS. Daily intraperitoneal administration of minocycline, an inhibitor of microglial/macrophage function prevented hyper-ramification and abolished the pre-TP-LI evoked synaptic changes. In contrast, rat-pups injured after the TP (post-TP) exhibited microglia hypo-ramification in the nTS and had increased sEPSC amplitudes/frequencies, and decreased TS-eEPSC amplitudes compared to controls. These synaptic changes were not associated with changes in CI-AMPARs, and instead involved greater TS-evoked use-dependent depression (reduced paired pulse ratio), which is a hallmark of presynaptic plasticity. Thus we conclude that LI regulates the efficacy of TS → nTS synapses through discrete plasticity-like mechanisms that are immune-mediated and depend on whether the injury occurs before or after the TP for respiratory control.

A Framework for Patient State Tracking by Classifying Multiscalar Physiologic Waveform Features

OBJECTIVE

state-of-the-art algorithms that quantify nonlinear dynamics in physiologic waveforms are underutilized clinically due to their esoteric nature. We present a generalizable framework for classifying multiscalar waveform features, designed for patient-state tracking directly at the bedside.

METHODS

an artificial neural network classifier was designed to evaluate multiscale waveform features against a fingerprint database of multifractal synthetic time series. The results are mapped into a physiologic state space for near real-time patient-state tracking.

RESULTS

the framework was validated on cardiac beat-to-beat dynamics processed with the multiscale entropy algorithm, and assessed using PhysioNet databases. We then applied our algorithm to predict 28-day mortality for sepsis patients, and found it had greater prognostic accuracy than standard clinical severity scores.

CONCLUSION

we developed a novel framework to classify multiscale features of beat-to-beat dynamics, and performed an initial clinical validation to demonstrate that our approach generates a robust quantification of a patient's state, compatible with real-time bedside implementations.

SIGNIFICANCE

the framework generates meaningful and actionable patient-specific information, and could facilitate the dissemination of a new class of "always-on" diagnostic tools.

Kölliker-Fuse nuclei regulate respiratory rhythm variability via a gain-control mechanism

Respiration varies from breath to breath. On the millisecond timescale of spiking, neuronal circuits exhibit variability due to the stochastic properties of ion channels and synapses. Does this fast, microscopic source of variability contribute to the slower, macroscopic variability of the respiratory period? To address this question, we modeled a stochastic oscillator with forcing; then, we tested its predictions experimentally for the respiratory rhythm generated by the in situ perfused preparation during vagal nerve stimulation (VNS). Our simulations identified a relationship among the gain of the input, entrainment strength, and rhythm variability. Specifically, at high gain, the periodic input entrained the oscillator and reduced variability, whereas at low gain, the noise interacted with the input, causing events known as "phase slips", which increased variability on a slow timescale. Experimentally, the in situ preparation behaved like the low-gain model: VNS entrained respiration but exhibited phase slips that increased rhythm variability. Next, we used bilateral muscimol microinjections in discrete respiratory compartments to identify areas involved in VNS gain control. Suppression of activity in the nucleus tractus solitarii occluded both entrainment and amplification of rhythm variability by VNS, confirming that these effects were due to the activation of the Hering-Breuer reflex. Suppressing activity of the Kölliker-Fuse nuclei (KFn) enhanced entrainment and reduced rhythm variability during VNS, consistent with the predictions of the high-gain model. Together, the model and experiments suggest that the KFn regulates respiratory rhythm variability via a gain control mechanism.

Effects of ion channel noise on neural circuits: an application to the respiratory pattern generator to investigate breathing variability

Neural activity generally displays irregular firing patterns even in circuits with apparently regular outputs, such as motor pattern generators, in which the output frequency fluctuates randomly around a mean value. This "circuit noise" is inherited from the random firing of single neurons, which emerges from stochastic ion channel gating (channel noise), spontaneous neurotransmitter release, and its diffusion and binding to synaptic receptors. Here we demonstrate how to expand conductance-based network models that are originally deterministic to include realistic, physiological noise, focusing on stochastic ion channel gating. We illustrate this procedure with a well-established conductance-based model of the respiratory pattern generator, which allows us to investigate how channel noise affects neural dynamics at the circuit level and, in particular, to understand the relationship between the respiratory pattern and its breath-to-breath variability. We show that as the channel number increases, the duration of inspiration and expiration varies, and so does the coefficient of variation of the breath-to-breath interval, which attains a minimum when the mean duration of expiration slightly exceeds that of inspiration. For small channel numbers, the variability of the expiratory phase dominates over that of the inspiratory phase, and vice versa for large channel numbers. Among the four different cell types in the respiratory pattern generator, pacemaker cells exhibit the highest sensitivity to channel noise. The model shows that suppressing input from the pons leads to longer inspiratory phases, a reduction in breathing frequency, and larger breath-to-breath variability, whereas enhanced input from the raphe nucleus increases breathing frequency without changing its pattern.

NEW & NOTEWORTHY

A major source of noise in neuronal circuits is the "flickering" of ion currents passing through the neurons' membranes (channel noise), which cannot be suppressed experimentally. Computational simulations are therefore the best way to investigate the effects of this physiological noise by manipulating its level at will. We investigate the role of noise in the respiratory pattern generator and show that endogenous, breath-to-breath variability is tightly linked to the respiratory pattern.

C57BL/6J mouse apolipoprotein A2 gene is deterministic for apnea

RATIONALE

Brainstem apolipoprotein AII (apoa2) mRNA expression correlates with apnea in breathing present in the adult C57Bl/6J (B6) sleep apnea model.

OBJECTIVES

To test the hypothesis that the B6 apoa2 gene contributes to the trait, we performed plethysmographic testing in apoa2 knock out (KO: -/-) mice, an in situ brainstem-spinal cord preparation comparing KO to WT (+/+) mice, and B6xDBA recombinant inbred strains (RISs).

MEASUREMENTS AND MAIN RESULTS

Apoa2 WT do, but KO and heterozygote (+/-) mice do not exhibit apnea during post-hypoxic breathing, measured in vivo. In the in situ model, pauses and instability in fictive phrenic bursting are substantially reduced in KO vs. WT preparations. In 24 RISs, apnea number in vivo was higher in strains with B6 apoa2 than with DBA apoa2 alleles.

CONCLUSIONS

The B6 apoa2 polymorphism is directly involved in breath production, and its identification suggests a novel pathway influencing risk for adult sleep apnea.

Functional connectivity in raphé-pontomedullary circuits supports active suppression of breathing during hypocapnic apnea

Hyperventilation is a common feature of disordered breathing. Apnea ensues if CO2 drive is sufficiently reduced. We tested the hypothesis that medullary raphé, ventral respiratory column (VRC), and pontine neurons have functional connectivity and persistent or evoked activities appropriate for roles in the suppression of drive and rhythm during hyperventilation and apnea. Phrenic nerve activity, arterial blood pressure, end-tidal CO2, and other parameters were monitored in 10 decerebrate, vagotomized, neuromuscularly-blocked, and artificially ventilated cats. Multielectrode arrays recorded spiking activity of 649 neurons. Loss and return of rhythmic activity during passive hyperventilation to apnea were identified with the S-transform. Diverse fluctuating activity patterns were recorded in the raphé-pontomedullary respiratory network during the transition to hypocapnic apnea. The firing rates of 160 neurons increased during apnea; the rates of 241 others decreased or stopped. VRC inspiratory neurons were usually the last to cease firing or lose rhythmic activity during the transition to apnea. Mayer wave-related oscillations (0.04-0.1 Hz) in firing rate were also disrupted during apnea. Four-hundred neurons (62%) were elements of pairs with at least one hyperventilation-responsive neuron and a correlational signature of interaction identified by cross-correlation or gravitational clustering. Our results support a model with distinct groups of chemoresponsive raphé neurons contributing to hypocapnic apnea through parallel processes that incorporate disfacilitation and active inhibition of inspiratory motor drive by expiratory neurons. During apnea, carotid chemoreceptors can evoke rhythm reemergence and an inspiratory shift in the balance of reciprocal inhibition via suppression of ongoing tonic expiratory neuron activity.

Respiratory modulation of sympathetic activity is attenuated in adult rats conditioned with chronic hypobaric hypoxia

Respiratory modulation of sympathetic nerve activity (SNA) depends on numerous factors including prior experience. In our studies, exposing naïve adult, male Sprague-Dawley rats to acute intermittent hypoxia (AIH) enhanced respiratory-modulation of splanchnic SNA (sSNA); whereas conditioning them to chronic hypobaric hypoxia (CHH) attenuated modulation. Further, AIH can evoke increased SNA in the absence phrenic long-term facilitation. We hypothesized that AIH would restore respiratory modulation of SNA in CHH rats. In anesthetized, CHH-conditioned (0.5 atm, 2 wks) rats (n=16), we recorded phrenic and sSNA before during and after AIH (8% O2 for 45s every 5min for 1h). At baseline, sSNA was not modulated with respiration. The sSNA was not recruited during a single brief exposure of hypoxia nor after 10 repetitive exposures. Further, the sSNA chemoresponse was not restored 1h after completing AIH. Thus, CHH-conditioning blocked the short-term plasticity expressed in sympatho-respiratory efferent activities and this was associated with reduced respiratory modulation of sympathetic activity and with attenuation of the sympatho-respiratory chemoresponse.

Increased cardio-respiratory coupling evoked by slow deep breathing can persist in normal humans

Slow deep breathing (SDB) has a therapeutic effect on autonomic tone. Our previous studies suggested that coupling of the cardiovascular to the respiratory system mediates plasticity expressed in sympathetic nerve activity. We hypothesized that SDB evokes short-term plasticity of cardiorespiratory coupling (CRC). We analyzed respiratory frequency (fR), heart rate and its variability (HR&HRV), the power spectral density (PSD) of blood pressure (BP) and the ventilatory pattern before, during, and after a 20-min epoch of SDB. During SDB, CRC and the relative PSD of BP at fR increased; mean arterial pressure decreased; but HR varied; increasing (n = 3), or decreasing (n = 2) or remaining the same (n = 5). After SDB, short-term plasticity was not apparent for the group but for individuals differences existed between baseline and recovery periods. We conclude that a repeated practice, like pranayama, may strengthen CRC and evoke short-term plasticity effectively in a subset of individuals.

Learning to breathe: habituation of Hering-Breuer inflation reflex emerges with postnatal brainstem maturation

The Hering-Breuer (HBR) reflex is considered a major regulatory feedback for the generation and patterning of respiratory activity. While HBR is important in neonates, its significance in adults is controversial. Previous experiments that investigated the plasticity of entrainment of the respiratory rhythm by vagal input demonstrated postnatal changes in HBR plasticity. Here we analyzed postnatal changes in the plasticity of HBR by mimicking the classic lung inflation tests with repetitive tonic vagal stimulation across different postnatal stages in an in situ perfused brainstem preparation of rat. The study shows that neonates stereotypically exhibit HBR stimulus-dependent prolongation of expiration while juvenile preparations (>postnatal day 16) showed significant habituation of HBR following repetitive stimulation. Subsequent experiments employing physiological lung inflation tests in situ confirmed HBR habituation in juveniles. We conclude that postnatal emergence of HBR habituation explains the weak contribution and high activation threshold of HBR in the regulation of eupnea.

Cardiorespiratory coupling: common rhythms in cardiac, sympathetic, and respiratory activities

Cardiorespiratory coupling is an encompassing term describing more than the well-recognized influences of respiration on heart rate and blood pressure. Our data indicate that cardiorespiratory coupling reflects a reciprocal interaction between autonomic and respiratory control systems, and the cardiovascular system modulates the ventilatory pattern as well. For example, cardioventilatory coupling refers to the influence of heart beats and arterial pulse pressure on respiration and is the tendency for the next inspiration to start at a preferred latency after the last heart beat in expiration. Multiple complementary, well-described mechanisms mediate respiration's influence on cardiovascular function, whereas mechanisms mediating the cardiovascular system's influence on respiration may only be through the baroreceptors but are just being identified. Our review will describe a differential effect of conditioning rats with either chronic intermittent or sustained hypoxia on sympathetic nerve activity but also on ventilatory pattern variability. Both intermittent and sustained hypoxia increase sympathetic nerve activity after 2 weeks but affect sympatho-respiratory coupling differentially. Intermittent hypoxia enhances sympatho-respiratory coupling, which is associated with low variability in the ventilatory pattern. In contrast, after constant hypobaric hypoxia, 1-to-1 coupling between bursts of sympathetic and phrenic nerve activity is replaced by 2-to-3 coupling. This change in coupling pattern is associated with increased variability of the ventilatory pattern. After baro-denervating hypobaric hypoxic-conditioned rats, splanchnic sympathetic nerve activity becomes tonic (distinct bursts are absent) with decreases during phrenic nerve bursts and ventilatory pattern becomes regular. Thus, conditioning rats to either intermittent or sustained hypoxia accentuates the reciprocal nature of cardiorespiratory coupling. Finally, identifying a compelling physiologic purpose for cardiorespiratory coupling is the biggest barrier for recognizing its significance. Cardiorespiratory coupling has only a small effect on the efficiency of gas exchange; rather, we propose that cardiorespiratory control system may act as weakly coupled oscillator to maintain rhythms within a bounded variability.

Physiological and pathophysiological interactions between the respiratory central pattern generator and the sympathetic nervous system

Respiratory modulation seen in the sympathetic nerve activity (SNA) implies that the respiratory and sympathetic networks interact. During hypertension elicited by chronic intermittent hypoxia (CIH), the SNA displays an enhanced respiratory modulation reflecting strengthened interactions between the networks. In this chapter, we review a series of experimental and modeling studies that help elucidate possible mechanisms of sympatho-respiratory coupling. We conclude that this coupling significantly contributes to both the sympathetic baroreflex and the augmented sympathetic activity after exposure to CIH. This conclusion is based on the following findings. (1) Baroreceptor activation results in perturbation of the respiratory pattern via transient activation of postinspiratory neurons in the Bötzinger complex (BötC). The same BötC neurons are involved in the respiratory modulation of SNA, and hence provide an additional pathway for the sympathetic baroreflex. (2) Under hypercapnia, phasic activation of abdominal motor nerves (AbN) is accompanied by synchronous discharges in SNA due to the common source of this rhythmic activity in the retrotrapezoid nucleus (RTN). CIH conditioning increases the CO2 sensitivity of central chemoreceptors in the RTN which results in the emergence of AbN and SNA discharges under normocapnic conditions similar to those observed during hypercapnia in naïve animals. Thus, respiratory-sympathetic interactions play an important role in defining sympathetic output and significantly contribute to the sympathetic activity and hypertension under certain physiological or pathophysiological conditions, and the theoretical framework presented may be instrumental in understanding of malfunctioning control of sympathetic activity in a variety of disease states.

Pontine mechanisms of respiratory control

Pontine respiratory nuclei provide synaptic input to medullary rhythmogenic circuits to shape and adapt the breathing pattern. An understanding of this statement depends on appreciating breathing as a behavior, rather than a stereotypic rhythm. In this review, we focus on the pontine-mediated inspiratory off-switch (IOS) associated with postinspiratory glottal constriction. Further, IOS is examined in the context of pontine regulation of glottal resistance in response to multimodal sensory inputs and higher commands, which in turn rules timing, duration, and patterning of respiratory airflow. In addition, network plasticity in respiratory control emerges during the development of the pons. Synaptic plasticity is required for dynamic and efficient modulation of the expiratory breathing pattern to cope with rapid changes from eupneic to adaptive breathing linked to exploratory (foraging and sniffing) and expulsive (vocalizing, coughing, sneezing, and retching) behaviors, as well as conveyance of basic emotions. The speed and complexity of changes in the breathing pattern of behaving animals implies that "learning to breathe" is necessary to adjust to changing internal and external states to maintain homeostasis and survival.

Inducible re-expression of HEXIM1 causes physiological cardiac hypertrophy in the adult mouse

AIMS

The transcription factor hexamethylene-bis-acetamide-inducible protein 1 (HEXIM1) regulates myocardial vascularization and growth during cardiogenesis. Our aim was to determine whether HEXIM1 also has a beneficial role in modulating vascularization, myocardial growth, and function within the adult heart.

METHODS AND RESULTS

To achieve our objective, we created and investigated a mouse line wherein HEXIM1 was re-expressed in adult cardiomyocytes to levels found in the foetal heart. Our findings support a beneficial role for HEXIM1 through increased vascularization, myocardial growth, and increased ejection fraction within the adult heart. HEXIM1 re-expression induces angiogenesis, that is, essential for physiological hypertrophy and maintenance of cardiac function. The ability of HEXIM1 to co-ordinate processes associated with physiological hypertrophy may be attributed to HEXIM1 regulation of other transcription factors (HIF-1-α, c-Myc, GATA4, and PPAR-α) that, in turn, control many genes involved in myocardial vascularization, growth, and metabolism. Moreover, the mechanism for HEXIM1-induced physiological hypertrophy appears to be distinct from that involving the PI3K/AKT pathway.

CONCLUSION

HEXIM1 re-expression results in the induction of angiogenesis that allows for the co-ordination of tissue growth and angiogenesis during physiological hypertrophy.

Decreased Hering-Breuer input-output entrainment in a mouse model of Rett syndrome

Rett syndrome, a severe X-linked neurodevelopmental disorder caused by mutations in the gene encoding methyl-CpG-binding protein 2 (Mecp2), is associated with a highly irregular respiratory pattern including severe upper-airway dysfunction. Recent work suggests that hyperexcitability of the Hering-Breuer reflex (HBR) pathway contributes to respiratory dysrhythmia in Mecp2 mutant mice. To assess how enhanced HBR input impacts respiratory entrainment by sensory afferents in closed-loop in vivo-like conditions, we investigated the input (vagal stimulus trains) - output (phrenic bursting) entrainment via the HBR in wild-type and MeCP2-deficient mice. Using the in situ perfused brainstem preparation, which maintains an intact pontomedullary axis capable of generating an in vivo-like respiratory rhythm in the absence of the HBR, we mimicked the HBR feedback input by stimulating the vagus nerve (at threshold current, 0.5 ms pulse duration, 75 Hz pulse frequency, 100 ms train duration) at an inter-burst frequency matching that of the intrinsic oscillation of the inspiratory motor output of each preparation. Using this approach, we observed significant input-output entrainment in wild-type mice as measured by the maximum of the cross-correlation function, the peak of the instantaneous relative phase distribution, and the mutual information of the instantaneous phases. This entrainment was associated with a reduction in inspiratory duration during feedback stimulation. In contrast, the strength of input-output entrainment was significantly weaker in Mecp2 (-/+) mice. However, Mecp2 (-/+) mice also had a reduced inspiratory duration during stimulation, indicating that reflex behavior in the HBR pathway was intact. Together, these observations suggest that the respiratory network compensates for enhanced sensitivity of HBR inputs by reducing HBR input-output entrainment.

Respiratory and behavioral dysfunction following loss of the GABAA receptor α4 subunit

γ-Aminobutyric acid type A (GABAA) receptor plasticity participates in mediating adaptation to environmental change. Previous studies in rats demonstrated that extrasynaptic GABAA receptor subunits and receptors in the pons, a brainstem region involved in respiratory control, are upregulated by exposure to sustained hypobaric hypoxia. In these animals, expression of the mRNA encoding the extrasynaptic α4 subunit rose after 3 days in sustained hypoxia, while those encoding the α6 and δ subunits increased dramatically by 2 weeks. However, the participation of extrasynaptic subunits in maintaining respiration in normoxic conditions remains unknown. To examine the importance of α4 in a normal environment, respiratory function, motor and anxiety-like behaviors, and expression of other GABAA receptor subunit mRNAs were compared in wild-type (WT) and α4 subunit-deficient mice. Loss of the α4 subunit did not impact frequency, but did lead to reduced ventilatory pattern variability. In addition, mice lacking the subunit exhibited increased anxiety-like behavior. Finally, α4 subunit loss resulted in reduced expression of other extrasynaptic (α6 and δ) subunit mRNAs in the pons without altering those encoding the most prominent synaptic subunits. These findings on subunit-deficient mice maintained in normoxia, in conjunction with earlier findings on animals maintained in chronic hypoxia, suggest that the expression and regulation of extrasynaptic GABAA receptor subunits in the pons is interdependent and that their levels influence respiratory control as well as adaptation to stress.

Control of breathing by interacting pontine and pulmonary feedback loops

The medullary respiratory network generates respiratory rhythm via sequential phase switching, which in turn is controlled by multiple feedbacks including those from the pons and nucleus tractus solitarii; the latter mediates pulmonary afferent feedback to the medullary circuits. It is hypothesized that both pontine and pulmonary feedback pathways operate via activation of medullary respiratory neurons that are critically involved in phase switching. Moreover, the pontine and pulmonary control loops interact, so that pulmonary afferents control the gain of pontine influence of the respiratory pattern. We used an established computational model of the respiratory network (Smith et al., 2007) and extended it by incorporating pontine circuits and pulmonary feedback. In the extended model, the pontine neurons receive phasic excitatory activation from, and provide feedback to, medullary respiratory neurons responsible for the onset and termination of inspiration. The model was used to study the effects of: (1) "vagotomy" (removal of pulmonary feedback), (2) suppression of pontine activity attenuating pontine feedback, and (3) these perturbations applied together on the respiratory pattern and durations of inspiration (T(I)) and expiration (T(E)). In our model: (a) the simulated vagotomy resulted in increases of both T(I) and T(E), (b) the suppression of pontine-medullary interactions led to the prolongation of T(I) at relatively constant, but variable T(E), and (c) these perturbations applied together resulted in "apneusis," characterized by a significantly prolonged T(I). The results of modeling were compared with, and provided a reasonable explanation for, multiple experimental data. The characteristic changes in T(I) and T(E) demonstrated with the model may represent characteristic changes in the balance between the pontine and pulmonary feedback control mechanisms that may reflect specific cardio-respiratory disorders and diseases.

Quantifying interactions between real oscillators with information theory and phase models: application to cardiorespiratory coupling

Interactions between oscillators can be investigated with standard tools of time series analysis. However, these methods are insensitive to the directionality of the coupling, i.e., the asymmetry of the interactions. An elegant alternative was proposed by Rosenblum and collaborators [M. G. Rosenblum, L. Cimponeriu, A. Bezerianos, A. Patzak, and R. Mrowka, Phys. Rev. E 65, 041909 (2002); M. G. Rosenblum and A. S. Pikovsky, Phys. Rev. E 64, 045202 (2001)] which consists in fitting the empirical phases to a generic model of two weakly coupled phase oscillators. This allows one to obtain the interaction functions defining the coupling and its directionality. A limitation of this approach is that a solution always exists in the least-squares sense, even in the absence of coupling. To preclude spurious results, we propose a three-step protocol: (1) Determine if a statistical dependency exists in the data by evaluating the mutual information of the phases; (2) if so, compute the interaction functions of the oscillators; and (3) validate the empirical oscillator model by comparing the joint probability of the phases obtained from simulating the model with that of the empirical phases. We apply this protocol to a model of two coupled Stuart-Landau oscillators and show that it reliably detects genuine coupling. We also apply this protocol to investigate cardiorespiratory coupling in anesthetized rats. We observe reciprocal coupling between respiration and heartbeat and that the influence of respiration on the heartbeat is generally much stronger than vice versa. In addition, we find that the vagus nerve mediates coupling in both directions.

Isoflurane and ketamine anesthesia have different effects on ventilatory pattern variability in rats

We hypothesize that isoflurane and ketamine impact ventilatory pattern variability (VPV) differently. Adult Sprague-Dawley rats were recorded in a whole-body plethysmograph before, during and after deep anesthesia. VPV was quantified from 60-s epochs using a complementary set of analytic techniques that included constructing surrogate data sets that preserved the linear structure but disrupted nonlinear deterministic properties of the original data. Even though isoflurane decreased and ketamine increased respiratory rate, VPV as quantified by the coefficient of variation decreased for both anesthetics. Further, mutual information increased and sample entropy decreased and the nonlinear complexity index (NLCI) increased during anesthesia despite qualitative differences in the shape and period of the waveform. Surprisingly mutual information and sample entropy did not change in the surrogate sets constructed from isoflurane data, but in those constructed from ketamine data, mutual information increased and sample entropy decreased significantly in the surrogate segments constructed from anesthetized relative to unanesthetized epochs. These data suggest that separate mechanisms modulate linear and nonlinear variability of breathing.

Linking Inflammation, Cardiorespiratory Variability, and Neural Control in Acute Inflammation via Computational Modeling

Acute inflammation leads to organ failure by engaging catastrophic feedback loops in which stressed tissue evokes an inflammatory response and, in turn, inflammation damages tissue. Manifestations of this maladaptive inflammatory response include cardio-respiratory dysfunction that may be reflected in reduced heart rate and ventilatory pattern variabilities. We have developed signal-processing algorithms that quantify non-linear deterministic characteristics of variability in biologic signals. Now, coalescing under the aegis of the NIH Computational Biology Program and the Society for Complexity in Acute Illness, two research teams performed iterative experiments and computational modeling on inflammation and cardio-pulmonary dysfunction in sepsis as well as on neural control of respiration and ventilatory pattern variability. These teams, with additional collaborators, have recently formed a multi-institutional, interdisciplinary consortium, whose goal is to delineate the fundamental interrelationship between the inflammatory response and physiologic variability. Multi-scale mathematical modeling and complementary physiological experiments will provide insight into autonomic neural mechanisms that may modulate the inflammatory response to sepsis and simultaneously reduce heart rate and ventilatory pattern variabilities associated with sepsis. This approach integrates computational models of neural control of breathing and cardio-respiratory coupling with models that combine inflammation, cardiovascular function, and heart rate variability. The resulting integrated model will provide mechanistic explanations for the phenomena of respiratory sinus-arrhythmia and cardio-ventilatory coupling observed under normal conditions, and the loss of these properties during sepsis. This approach holds the potential of modeling cross-scale physiological interactions to improve both basic knowledge and clinical management of acute inflammatory diseases such as sepsis and trauma.

A method for analyzing temporal patterns of variability of a time series from Poincare plots

The Poincaré plot is a popular two-dimensional, time series analysis tool because of its intuitive display of dynamic system behavior. Poincaré plots have been used to visualize heart rate and respiratory pattern variabilities. However, conventional quantitative analysis relies primarily on statistical measurements of the cumulative distribution of points, making it difficult to interpret irregular or complex plots. Moreover, the plots are constructed to reflect highly correlated regions of the time series, reducing the amount of nonlinear information that is presented and thereby hiding potentially relevant features. We propose temporal Poincaré variability (TPV), a novel analysis methodology that uses standard techniques to quantify the temporal distribution of points and to detect nonlinear sources responsible for physiological variability. In addition, the analysis is applied across multiple time delays, yielding a richer insight into system dynamics than the traditional circle return plot. The method is applied to data sets of R-R intervals and to synthetic point process data extracted from the Lorenz time series. The results demonstrate that TPV complements the traditional analysis and can be applied more generally, including Poincaré plots with multiple clusters, and more consistently than the conventional measures and can address questions regarding potential structure underlying the variability of a data set.