- Review
- Open access
- Published:
Temporal variations in the pattern of breathing: techniques, sources, and applications to translational sciences
The Journal of Physiological Sciences volume 72, Article number: 22 (2022)
Abstract
The breathing process possesses a complex variability caused in part by the respiratory central pattern generator in the brainstem; however, it also arises from chemical and mechanical feedback control loops, network reorganization and network sharing with nonrespiratory motor acts, as well as inputs from cortical and subcortical systems. The notion that respiratory fluctuations contain hidden information has prompted scientists to decipher respiratory signals to better understand the fundamental mechanisms of respiratory pattern generation, interactions with emotion, influences on the cortical neuronal networks associated with cognition, and changes in variability in healthy and disease-carrying individuals. Respiration can be used to express and control emotion. Furthermore, respiration appears to organize brain-wide network oscillations via cross-frequency coupling, optimizing cognitive performance. With the aid of information theory-based techniques and machine learning, the hidden information can be translated into a form usable in clinical practice for diagnosis, emotion recognition, and mental conditioning.
Introduction
Humans breathe throughout their lives; from birth to death, our breaths are always changing, and no single breath is identical to another. Where do these variations originate from? What information can we extract from the variations in breathing? What are the physiological and pathophysiological implications of these variations? Breath-to-breath variations in breathing patterns can occur as uncorrelated random variations (white noise), correlated random changes, periodic variations, or nonrandom, nonperiodic fluctuations [1]. This review focuses on these temporal variations in the pattern of breathing, including the complex variability of breathing and characteristic breathing patterns.
Definition of terms The term “variability” is used in different ways. For example, heart rate variability refers to small fluctuations in the time interval between heartbeats. However, the term is also used to describe more diverse variations in the pattern of breathing. In this review, the term “respiratory variability” is used in a broader sense, as equivalent to “breath-to-breath variations in breathing patterns” described by Bruce [1]. The term “complexity”, synonymous with “complex variability”, has a meaning that is qualitatively and quantitatively distinguishable from traditional concepts and metrics of “variability” [2]. Although no formal definition of this term exists, it may be said that “complexity” is a measure of the amount of information and unpredictability, which is quantified with the information theory-based techniques mentioned in "Techniques for analyzing respiratory variability" section. As revealed in "Respiratory variability in health and disease" and "Future directions in translational sciences" sections, the concept of complex variability can be usefully applied to translational sciences. The term “nonlinearity” refers to the relationship between system inputs and outputs where the latter is not directly proportional to the former; thus, small input perturbations can cause large effects on the outputs [2]. A linear system, but not a nonlinear system, can be described by an autoregressive model or frequency response characteristics.
Recognizing the presence of 'hidden information' in physiological time series necessitates the use of fluctuation analysis techniques in statistical physics [3,4,5,6,7], which are unfamiliar to physiologists. Therefore, this review begins by summarizing techniques for analyzing complex variability and then discusses the sources of respiratory variability and translations of the information hidden therein to health-related sciences.
Techniques for analyzing respiratory variability
A common measure of gross variability is the coefficient of variation (CV), which is defined as the ratio of the standard deviation to the mean. When datasets with varied units or considerably different means are compared, the CV should be used instead of the standard deviation. The root mean square successive difference (RMSSD) measures the extent of variability between successive time points; for example, the RMSSD of the interheartbeat interval can be used as an index for monitoring changes in parasympathetic activity [8].
Respiratory variability is produced through an integrated process that involves multifunctional control mechanisms in the brain; therefore, the characteristics of present breaths are correlated with the characteristics of past breaths. The autocorrelation (AR), which is defined as the correlation between a signal and a delayed copy of itself as a function of the delay, is a common correlation metric for discrete time-series data. For example, the white noise exhibits zero correlation with any nonzero time lag. On the other hand, if a signal is not random, one or more of the autocorrelations remains significantly nonzero. Previous studies on respiratory variability have typically used the AR at one breath lag [9].
In well-controlled experimental settings, respiratory control systems can be regarded as stationary. In such cases, the mean, standard deviation, CV, and AR remain invariant throughout the observation period. However, this is not always the case; in particular, during long-term observations, e.g., overnight monitoring, respiratory control systems become nonstationary. Detrended fluctuation analysis is a scaling analysis approach that was originally designed to quantify long-range power-law correlations in signals; however, it can also be used to investigate both long-term (LTCs) and short-term correlations (STCs) in nonstationary systems [10, 11]. Detrended fluctuation analysis uses the integrated fluctuation of a signal; the integrated time series is divided into N epochs of length n, and each epoch is detrended with a least squares fit, yielding a locally detrended time-series segment \({x}_{k}^{n}(t), k=1, \dots ,\) N. The average fluctuation for a given epoch is calculated as
The LTCs or STCs can then be extracted by the scaling exponent, α, which is the slope of \(\frac{\mathrm{log}(F\left(n\right))}{\mathrm{log}(n)}\) for a specific range of n [12].
Deterministic dynamical systems that follow a unique path or evolution can exhibit complex behaviors. One example is the Rössler attractor (Fig. 1), which is a system composed of three nonlinear ordinary differential equations. In such a system, a small difference in initial conditions can result in significantly different behaviors, and predicting future behavior becomes progressively more difficult with time, a traditional indicator of ‘chaos’. The largest Lyapunov exponent (LLE) measures the predictability of a system’s behavior; a positive LLE indicates that the attractor diverges, i.e., chaotic behavior [13].
Visualization of the Rössler attractor, a continuous-time dynamical system that exhibits chaotic dynamics with fractal properties. The system is composed of three nonlinear differential equations and thus has three variables. a Temporal fluctuations of the x variable. Although each variable oscillates within a fixed range of values, its amplitude is highly variable. b Three-dimensional trajectory of the Rössler attractor visualized using the values of all three variables. c Three-dimensional plot of the Rössler attractor reconstructed using delay-time embedding (\(\tau =30\,\mathrm{ a}.\mathrm{u}.\)) of a single variable (x). Note that the complexity of the system is well preserved compared to b, in which all three variables are used. d Recurrence plot of the reconstructed attractor using delay-time embedding (\(m=3, \tau =30\,\mathrm{ a}.\mathrm{u}.\), r = 1.0) consisting of diagonal lines with different lengths, indicating that the trajectories pass neighborhoods of past trajectories with limited durations. This means that the reconstructed attractor exhibits a fractal structure or self-similarity
The complexity of the variability of a system can be visualized through state space reconstructions. These visualizations can be obtained from a single observable variable using delay-time embedding (Fig. 1c). If we have time-series data \(u\left(i\right),i=1, 2,\dots ,N\), the state space can be reconstructed as an m-dimensional state, \({X}_{i}^{m}=\{u\left(i\right), u\left(i+\tau \right), \dots , u\left(i+(m-1)\tau \right)\}\), where m is the embedding dimension and \(\tau\) is the time delay. If we choose m and \(\tau\) properly, then the complexity of the variability is well characterized by the reconstructed attractor [14, 15].
As seen in three-dimensional plots of the Rössler attractor (Fig. 1b), the trajectories never return to past trajectories but do pass near them. When we plot the neighborhoods \({X}_{j}\) of \({X}_{i}\) whose norm is within \(r\), we can visualize these state recurrences [16]. This is called a recurrence plot, and it is mapped as
where \(\Vert \cdot \Vert\) is a norm, \(r\) is a threshold distance, and \(\Theta\) is the Heaviside function, which has a value of 1 if the expression in the inner parenthesis has a value greater than zero and a value of zero otherwise. Diagonal lines of different lengths that appear in the recurrence plots (Fig. 1d) indicate that a state visits the same region of the attractor at different times [17]. Eckmann et al. [16] showed that the length of these short upward diagonal lines is inversely proportional to the LLE. Based on the probability of diagonal lines with specific lengths appearing, the amount of information or uncertainty can be calculated using an entropy measure, the Shannon entropy (ShEn) [17], which is formulated as
where \(P\left({L}_{i}\right) is\) is the probability that diagonal lines of length \({L}_{i}\) appear.
Attractors that originate from complex systems often exhibit fractal structure or self-similarity, where similar patterns appear at increasingly small scales. A widely used index for characterizing fractal structures is the correlation dimension, which was proposed by Grassberger and Procaccia [18]. For any positive number r and embedding dimension m, the correlation sum \(C\left(r\right)\), which is the discrete version of the correlation integral, is defined as the fraction of pairs in the time delay embedding vector \({X}_{i}^{m}\) that have distances smaller than r,
where \({N}^{^{\prime}}=N-(m-1)\tau\), \(\Vert \cdot \Vert\) is the Euclidian norm, and \(\Theta\) is the Heaviside function [19]. If an attractor reconstructed by delay-time embedding exhibits a fractal structure, there exists a region in \(r\) where \(\frac{\mathrm{log}(C\left(r\right))}{\mathrm{log}(r)}\) is linear. The correlation dimension is defined as the slope of this scaling region. Figure 2 shows the attractor reconstruction, recurrence plot, and correlation dimension estimation from an actual respiratory signal.
Visualization and quantification of the complexity of actual experimental data. The data, obtained from an open dataset [225], contain a human respiratory signal while the participant watched a scary video. The respiratory signal is the temporal change in chest dimension associated with the expansion and contraction of the chest cavity, which was measured using a Hall effect sensor placed high on the torso. a The temporal change in chest dimension, sampled at 10 Hz, was used as the x variable for attractor reconstruction (b) and to generate the recurrence plot (c). b Three-dimensional plots of an attractor reconstructed using delay-time embedding (\(\tau =0.6\,\mathrm{ s})\). Variables y and z represent respiratory signals 0.6 s and 1.2 s advanced relative to the x variable, respectively. c Recurrence plot of human respiration while the participant watched a scary video (\(m=3, \tau =0.6\,\mathrm{ s}\), r = 0.8). d The log–log plot of a distance r versus the correlation sum C(r) has a linear scaling region, indicative of a fractal property. e The slope of the log–log plot (\(\frac{\mathrm{log}(C\left(r\right))}{\mathrm{log}(r)}\)) converges to 1.5 as the embedding dimension grows
The amount of information and the unpredictability of fluctuations in time-series data can be calculated directly using another entropy measure, the approximate entropy (ApEn) [9, 20]. The algorithm first estimates the appearance frequency of a similar pattern of sequences by evaluating whether a sequence of data points with length m is similar to other sequences in the data, with an allowed distance of \(r\) between the points. For each sequence \({X}_{i}^{m}, i=1, \dots , N-m+1\), the appearance frequency is defined as
where m is the length of the sequence, N is the number of data points, \(\Vert \cdot \Vert\) is the norm, and \(\Theta\) is the Heaviside function. Then, ApEn can be measured using the average logarithmic appearance frequency:
The sample entropy (SampEn) is a modification of the ApEn that is used to assess the complexity of physiological time-series signals and to diagnose disease states [21]. SampEn is defined as
Correlated variability arises from either a stochastic process (e.g., Brownian motion) or a deterministic nonlinear process. Due to various measurement restrictions, respiratory data have a limited length. As a result, any nonlinear dynamics that originate from deterministic processes cannot be distinguished from correlated variations caused by stochastic events [22]. To detect nonlinearity in the observed data, the surrogate test has often been used. In this test, surrogate data are generated based on a model or a combination of Fourier and inverse Fourier transformations. The surrogate data share the same statistical properties, such as the AR and power spectrum, as the original data, but do not retain their nonlinear properties [23]. Then, a null hypothesis against the presence of nonlinearity is tested using Monte Carlo methods [23, 24]. Namely, a discriminating statistic of nonlinearity is calculated for the original and all surrogate data, and the null hypothesis is rejected if the value for the original dataset is significantly different from that of the surrogate dataset.
A major obstacle for detecting the nonlinearity of respiratory variability is that nonlinear characterization methods, such as Lyapunov exponents and correlation analyses, are sensitive to both uncorrelated and correlated noise [25]. In addition, although the surrogate test infers the presence of nonlinearity, the results are not sufficient for determining the presence of chaos [26]. An alternative method that can circumvent these limitations is noise titration. In this method, noise is gradually added to the observed data by increasing the standard deviation of the data. The noise limit (NL) is defined as the standard deviation at which nonlinearity cannot be detected by the Volterra–Wiener algorithm, with NL > 0 indicating the presence of chaos [25]. There are some circumstances in which noise titration fails to distinguish colored noise from chaos; however, a ‘remedy’ has been provided to address these circumstances [27].
Sources of respiratory variability
The respiratory control system is complex; while the primary goal of the respiratory system is to maintain arterial blood gas homeostasis, if the work of breathing is too costly, the homeostatic maintenance function can be compromised, and the magnitude and pattern of the respiratory motor output are optimized in terms of a cost function [28, 29]. Swallowing, coughing, sighing, and other nonrespiratory motor acts reset the respiratory rhythm produced by the respiratory central pattern generator (rCPG) in the brainstem. Respiratory neurons in the Bötzinger complex and ventral respiratory group are involved in generating the spatiotemporally organized activities associated with coughing and swallowing, and some respiratory neuronal networks are shared by nonrespiratory networks [30, 31]. Spatial network reorganization, i.e., the expansion and contraction of the active network during the inspiratory phase of breathing, occurs within the rCPG via a balance between excitation and inhibition [32]. This network sharing and reorganization contribute to the flexibility and variability of breathing [33, 34]. Furthermore, volitional and emotional controls of breathing can take control over the pattern of breathing, either consciously or unconsciously, via direct projections to respiratory motoneurons and projections to diverse respiratory control areas in the midbrain, pons, and medulla oblongata [35,36,37,38]. This multifaceted control produces variations in breathing that are not random but have some structure inherited from past breaths. In this section, sources of respiratory variability at each level—from automatic control by respiratory rhythmic cores in the brainstem to respiratory control by cortical and subcortical systems—are reviewed.
Complex variability intrinsic to the rCPG
The characteristics of a given breath were found to be dependent on the characteristics of the immediately preceding breath in paralyzed, artificially ventilated, and vagotomized cats whose spinal cords were cut at the T1 level [39]. This suggests that the nonrandom variability of breathing originates at least in part in the rCPG. The pre-Bötzinger complex (preBötC), which is located in the ventrolateral medulla, is the respiratory rhythmic core that generates the bursting activity that triggers inspiration [40,41,42]. Another respiratory oscillator, the parafacial respiratory group (pFRG, also known as the lateral parafacial), which is located ventrolateral to the facial nucleus, is an expiratory rhythm generator that causes active expiration during states of elevated respiratory drive [40, 43, 44]. Neurons in the preBötC autonomously exhibit periodic bursting activity due to their channel properties even if they are isolated and produce periodic inspiration-related synchronized activity in slice preparations [45, 46].
Del Negro et al. [47] analyzed respiratory variability at the neuron population level in a highly reduced, but rhythmically active preparation by measuring the integrated preBötC activity. A Poincaré map, which is an intersection of the state space, was obtained by plotting \(\int {PBC}_{n+1}\) against \(\int {PBC}_{n}\), where \(\int {PBC}_{n}\) represents the n-th integrated preBötC activity. They observed discrete transitions in the Poincaré maps from periodic oscillations to mixed-mode periodicity, quasiperiodicity, and finally disorganized aperiodic activity, with progressive increases in neuronal excitability, suggesting that the preBötC produces neural activity characteristic of a nonlinear dynamic system. Koshiya et al. [48] applied a voltage-sensitive dye imaging technique to rhythmically active slices and recorded spatiotemporal preBötC activity. They observed that the center of activity, which is calculated according to the magnitude of fluorescence intensity, moved within the preBötC during neuronal population bursts. The state space reconstructed from the moving speed of the center of activity was quasiperiodic, and a correlation dimension analysis and surrogate test suggested the presence of nonlinear dynamics.
Carroll and Ramirez [49] investigated cycle-to-cycle variability during preBötC neuronal recruitment using a multielectrode recording technique and found that respiratory neurons were stochastically activated with each burst. Furthermore, they found that the burst onset variability could not be reproduced in fully interconnected computational models but could be reproduced in sparsely connected network models with as little as 1% connectivity. However, this estimate was based on a randomly connected network that included only excitatory neurons; the preBötC neuronal network includes both excitatory and inhibitory neurons [50, 51]. A higher burst onset variability may be achieved with a higher fraction of all-to-all synaptic connections in a more realistic network with both excitatory and inhibitory neurons since the interburst interval becomes variable when inhibitory connections are included [52].
The pFRG was first identified as a presumptive rhythm generator that triggers the inspiratory pattern generator [53] in the brainstem–spinal cord preparations of neonatal rodents developed by Suzue [54]. The pFRG partially overlaps with the retrotrapezoid nucleus (RTN), which is located ventromedial to the facial nucleus, contains chemosensitive cells and distributes a CO2-dependent excitatory drive to the respiratory network [45, 55,56,57]. In more mature, intact preparations, neurons in the pFRG appear to be quiescent; however, they generate late-expiratory bursts of action potentials when they are disinhibited or activated [58]. Therefore, the pFRG is postulated to be a conditional oscillator for active expiration [40, 43]. Neurons in the preBötC and pFRG are bidirectionally connected; therefore, these inspiratory and expiratory oscillators are coupled. When the excitability of the preBötC network decreases, the inspiratory bursts skip their expected timings in an unpredictable manner; thus, quantal slowing of the respiratory rhythm, a phenomenon in which the respiratory rhythm jumps nondeterministically to integer multiples of the control period, occurs [59]. This quantal slowing could be caused by transmission failure from the pFRG to preBötC networks due to suppressed or stochastic excitatory synaptic transmission [59, 60]. Alternatively, quantal slowing could result from a breakdown of synchronized bursting in the preBötC [61]. In either case, the coupled oscillators in the preBötC and pFRG produce respiratory variability similar to atrioventricular blocks in the heart.
Normal breathing consists of three phases: inspiration, postinspiration, and late expiration, all of which are believed to originate in the rCPG in the brainstem [41, 62]. Sequential transection experiments have shown that the three-phase rhythm requires the integrity of the pontine–medullary respiratory network [63]. The pons plays two major roles in the rCPG, with both mediated by neuronal circuitry within the Kölliker-Fuse (KF) area [64]. First, the pons provides an inspiratory off-switch that causes an inspiratory-to-expiratory phase transition in conjunction with sensory feedback from slowly adapting pulmonary stretch receptors. Second, the pons regulates postinspiration, adjusting upper airway resistance during the respiratory cycle. Lesioning the pons results in a longer and more irregular inspiratory phase [65,66,67]. Stimulus after-effects (the prolongation of the inspiratory duration) have been found to be augmented following lesioning, suggesting that the pons plays an important role in the stability of the rCPG [67].
Yu et al. [68] investigated the effects of changing the input to the pons using a conductance-based model of the four different types of cells in the rCPG. The model shows that reduced pontine input causes longer inspiratory phases, reduced respiratory rate (RR), and increased breath-to-breath variability, consistent with the experimental findings. Furthermore, they investigated how channel noise affects neural dynamics at the circuit level. The model predicted that the expiratory phase is more variable than the inspiratory phase when the channel number is small, and vice versa when the channel number is large. Among the four different types of cells, the pacemaker cell exhibited the highest sensitivity to channel noise.
Astrocytic contributions to respiratory variability
Astrocytes respond to changes in neuronal network activity across brain states and behaviors [69] and modulate central pattern-generating motor circuits [70]. Astrocytes can even modulate brain-wide oscillations by transmitting the oxygenation status to higher cortical areas [69]. In the preBötC, the astrocyte glutamate–glutamine cycle and the supply of glutamine to neuronal glutamatergic terminals are essential for rhythm generation [71]. Furthermore, blocking the vesicular release of preBötC astrocytes reduces the resting breathing rate, lowers the frequency of periodic sighs, decreases rhythm variability, impairs respiratory responses to hypoxia and hypercapnia, and reduces exercise capacity [72].
Rhythmic, inwardly directed currents attributed to neuronal population bursts have been found in 10% of preBötC astrocytes [73]. Okada et al. [74] recorded the spatiotemporal activities of neurons and astrocytes in the preBötC using a calcium imaging technique. They found that a subset of astrocytes exhibited preinspiratory increases in the intracellular calcium concentration that were irregularly coupled with inspiratory neuronal bursts. In addition, they found that optogenetic stimulation of the astrocytes triggered action potentials in inspiratory neurons in the preBötC. Similar irregular coupling between the calcium activities of neurons and astrocytes has been reported in organotypic cultures of preBötC slices and pFRG slices [75, 76]. Network structure analyses with the cross-correlation technique and graph theory [77] revealed three separate but interconnected subnetworks: the glial, neuronal, and glial-neuronal networks [76]. These networks are organized into a small-world network structure commonly observed in biological networks [78]. On the other hand, half of the preBötC astrocytes showed synchronized low-frequency (0.023 Hz) oscillations; thus, a subset of the astrocytes forms a slow oscillator (Fig. 3) [73, 74, 79]. Therefore, the neurons and astrocytes in the preBötC and pFRG may form coupled slow and fast oscillators and can mutually interact, thus producing complex behaviors (Fig. 4a) [79, 80].
Representative calcium transients of EGFP-positive and EGFP-negative cells [79]. a Astrocytes are identified with GFAP promoter-controlled expression of EGFP. Scale bar: 40 µm. b Cross-correlation image between each pixel and the low-frequency calcium oscillation of cell #4 indicates that low-frequency astrocytic calcium oscillations are highly synchronized. c The top trace shows the respiratory rhythmic calcium fluctuation, and the other traces show bandpass filtered calcium signals of representative cells. For cells #1–8, the location of each cell is indicated in the left panels (a and b). For EGFP-positive cells (#1–5), the calcium signals of the cell (green) and its vicinity (light blue) are shown
Schematic drawings of coupled slow and fast oscillators that exhibit cross-frequency coupling (CFC). a Slow glial (red circles) and fast neuronal (black circles) oscillators are organized within the pFRG and preBötC and modulate their intrinsic rhythms [75]. The pFRG and preBötC conditionally couple, resulting in complex variability, such as quantal slowing. b The slow respiratory oscillator (black circle) orchestrates fast neuronal oscillations (blue circles) in the anterior cingulate, premotor, insular, and hippocampal cortices via CFC [151, 152]. pFRG parafacial respiratory group, preBötC pre-Bötzinger complex, ACC anterior cingulate cortex, PMC premotor cortex, INS insular cortex, HPC hippocampal cortex
Coupled oscillators desynchronize for sufficiently small couplings and then bifurcate to partially synchronized states when the coupling increases above a critical value [81]. Cross-frequency coupling (CFC), which is the interaction between oscillations in different frequency bands, is a widely observed phenomenon in the brain that may play a functional role in neuronal computation, communication, and learning [82]. There are various types of CFC, including phase-phase, phase-frequency, phase-amplitude, amplitude-amplitude, frequency-frequency, and amplitude-frequency (Fig. 5) [83]. In particular, phase-amplitude coupling changes quickly in response to sensory, motor, and cognitive events and correlates with performance in learning tasks (see "Control of respiration during cognition" section) [82]. However, the coupling between glial and neuronal oscillators may be more complex, since the coupling is not compatible with any known type of CFC. Therefore, this interaction may be revealed only through nonlinear analyses, such as cross-recurrence plot analysis [84, 85].
Representative modes of cross frequency interactions. Consider a slow oscillator (a) that interacts with a fast oscillator (b). In c, the power of the fast oscillation is modulated by the phase of the slow oscillation. This mode of interaction can be observed between the phase of respiration (slow oscillation) and the power of the gamma-band (40–150 Hz) oscillations during attention to breathing [151], and between the phase of respiration and the power of alpha-band (8–13 Hz) oscillations during a near-threshold spatial detection task [147]. In d, the frequency of the fast oscillation is modulated by the phase of the slow oscillation
Correlated variability and oscillations originating from the chemical control system
The integrated phrenic nerve activity (\(\int \mathrm{Phr}\)) responses to the electrical stimulation of a carotid sinus nerve (CSN) in anesthetized, paralyzed, vagotomized, and glomectomized cats have two distinct components: a rapid increase that accounts for approximately half of the full response and a gradual increase that eventually reaches a steady-state plateau [86]. After the stimulation stops, short-term potentiation (STP) of respiration is observed, in which the level of \(\int \mathrm{Phr}\) decreases rapidly but remains higher than the prestimulation level [86]. STP has also been observed after acute exposure to hypoxia, which is known as post-hypoxic persistent respiratory augmentation (PHRA) [87, 88]. Arundic acid, which is an astrocyte inhibitor, has been shown to suppress PHRA; however, astrocyte-specific Trpa1 knockout did not abolish PHRA, indicating that astrocytes mediate PHRA through mechanisms other than the putative ventilatory hypoxia sensor TRPA1 channels [89].
In addition to STP, several different types of memory effects have been shown to be associated with chemoreceptor afferent activation [90, 91]. CSN stimulation or hypoxia exposure increases \(\int \mathrm{Phr}\); however, the RR is below that of the steady-state baseline level. This phenomenon is known as the short-term depression (STD) of respiration. Injections of muscimol, a GABAA agonist, into the ventrolateral pons, where low current pulses evoke short-latency inhibition of phrenic nerve activity, have been shown to abolish STD [92]. This result suggests that the integrity of the ventrolateral pons is required for STD [92]. While STD involves a decrease in RR that occurs when hypoxic exposure is sustained for tens of seconds to a few minutes, hypoxic ventilatory depression (HVD) involves a decrease in tidal volume (VT) that occurs when moderate hypoxemia is sustained for 5–30 min [90]. Although the mechanism of HVD is unknown, recent findings have shown that astrocytes in the preBötC play an important role in counteracting HVD by releasing ATP to stimulate ventilation by activating P2Y1 receptors [93]. Repeated CSN stimulation or hypoxia exposure has been shown to result in long-term potentiation (LTP) of respiration in rats, where \(\int \mathrm{Phr}\) remained above that of the controls for at least 30 min. Vermalectomy eliminated LTP, which suggests that the cerebellar vermis plays a role in LTP [91].
The above experiments on STP, STD, and LTP were all conducted under open loop conditions, i.e., without chemical feedback control; therefore, these memory effects must be caused by neural mechanisms rather than chemical feedback mechanisms. However, correlated respiratory activity can also arise from chemical feedback [1, 94]. The VT of each breath in anesthetized, vagotomized, and spontaneously breathing rats was correlated with that immediately preceding breathing; however, such a correlation was not observed in the \(\int \mathrm{Phr}\) of anesthetized, vagotomized, paralyzed, and artificially ventilated rats [94]. Therefore, the autocorrelation originates from chemical feedback mechanisms.
Hypoxia, in combination with sleep and hypocapnia, can induce periodic breathing characterized by repeated clusters of two to five breaths interspaced with regularly spaced expiratory pauses [95]. During wakefulness, isocapnic hypoxia increases gross respiratory variability but decreases the autocorrelation coefficient at a lag of one breath for minute ventilation. The increase in respiratory variability can be decomposed into a random component and an oscillatory component, indicating that hypoxia induces hidden oscillations even in the absence of hypocapnia in healthy awake subjects [96]. The PaCO2 level also affects respiratory variability. Hypercapnia reduces respiratory variability and increases LLE, whereas hypocapnia increases respiratory variability and decreases LLE [97].
Phase resetting of the respiratory rhythm
The response dynamics of chemical feedback and mechanoreceptor feedback systems differ, with the latter depending on when the stimulus is presented [98]. For example, lung inflation during inspiration shortens the inspiratory time, whereas lung inflation during expiration prolongs the expiratory time. This suggests that the predicted timing of a certain phase, e.g., the onset of inspiration, can shift depending on the timing of a stimulus. This phenomenon is referred to as the phase resetting characteristic [99]. The phase resetting characteristics of the rCPG were first investigated by electrically stimulating the midbrain [100], which facilitates inspiration; however, these characteristics were most extensively studied by electrically stimulating the vagus nerve and superior laryngeal nerve (SLN) [1], which suppress inspiratory activity during most phases of the respiratory cycle. A plot of the time between the onset of the preceding inspiration and a stimulus (termed the old phase) versus the time between the stimulus and the onset of the following inspiration (termed the cophase) determines the topographical type of phase resetting [100]. If the stimulus is weak, a net change in the cophase as the old phase moves through one respiratory cycle becomes one respiratory cycle, and the topographical type is classified as type-1. If the stimulus is strong, the net change becomes 0, and the phase resetting is classified as type-0.
A brief SLN stimulation produces after-effects that can last for several cycles [101]. The stimulus after-effects depend on the timing of a given stimulus. Brief vagal nerve stimulation delivered during each respiratory cycle at mid-inspiration, mid-expiration, and late expiration near the expiratory-to-inspiratory phase transition can result in complex breath-to-breath variability [102]. Dhingra et al. [103] investigated vagal-dependent respiratory variability using information theory-based techniques and surrogate data testing. They found that the vagal afferent and dorsolateral pons both contribute to nonlinear variability in the pattern of breathing and are mutually dependent.
Repeated stimulation entrains respiration to the stimulation; however, the trajectories that return to the baseline state are not the same, resulting in respiratory variability. The stochastic noise associated with ion channel gating and synaptic neurotransmission affects the entrainment of respiratory rhythms to external periodic inputs [104]. Periodic vagal nerve stimulation (VNS) entrains the respiratory rhythm similar to the Hering-Breuer reflex; however, VNS increases respiratory variability because noise interacts with the input, leading to phase slips. Suppressing the activity of the KF region enhances entrainment and reduces rhythm variability during VNS, which suggests that the KF region regulates respiratory variability by controlling the gain of the Hering-Breuer reflex [104].
Swallowing
Swallowing is a physiological perturbation of the respiratory rhythm [105, 106]. In humans, swallowing during inspiration terminates the inspiratory phase, and respiration resumes with expiration, with a shorter duration than that of the control [107]. On the other hand, swallowing during expiration interrupts the expiratory phase; however, in this case, respiration normally resumes with expiration, increasing the total duration of the expiratory phase [107]. In general, swallowing acts as a strong perturbation in the rCPG, resulting in type-0 phase resetting [106, 107]. All rCPG neurons appear to be affected by swallowing, irrespective of the type of neuron. Inspiratory-augmenting, inspiratory-decrementing, and expiratory-augmenting neurons are all inhibited [31, 108, 109], while a majority of expiratory-decrementing neurons are activated in rats [108], and inspiratory-to-expiratory phase-spanning neurons are activated in guinea pigs [109]. Lesioning the KF area with ibotenic acid eliminates the respiratory phase resetting caused by swallowing, which suggests that the KF region plays an important role in coordinating breathing and swallowing [110, 111]. Furthermore, KF inhibition attenuates tonic postinspiratory vagal nerve activity and lowers the threshold for evoking swallowing. Therefore, the KF region plays a role in the airway-defensive laryngeal adductor reflex and gates the initiation of swallowing [112]. Since swallowing initiation is also inhibited by vagal feedback, dual peripheral and central gating mechanisms are involved in the coordination between breathing and swallowing [113].
Sighing
Sighs originate from a small ensemble of preBötC neurons [33, 114]. Physiological sighing requires peptidergic inputs from RTN/pFRG neurons, which express the bombesin-like neuropeptide neuromedin B or gastrin-releasing peptides [40, 114]. Physiological sighing is believed to be important in preventing the alveoli from collapsing (atelectasis), improving gas exchange [115, 116], and reducing hypoxia and hypercapnia [117]. However, a sigh is not only an augmented breath that maximally inflates the lung but also signals brain state changes, controls arousal, and regulates homeostasis of respiratory variability [118]. Vlemincx et al. [119] found that sighing increases autocorrelated respiratory variability and relieves mental stress. They hypothesized that sighs serve as psychophysiological resetters, restoring respiratory regulation by resetting the nonrandom respiratory variability when it becomes too low or too random [120,121,122]. Based on the theory of stochastic resonance [123], they postulated that an inappropriate level of respiratory variability compromises flexible and adaptive responsiveness or jeopardizes stability and hypothesized that a sigh acts as noise to restore healthy balanced respiratory variability [119]. Meanwhile, because autocorrelated respiratory variability arises from chemical feedback control [1, 94], the level of autocorrelated variability may reflect the relative contribution of the chemical feedback control to the total respiratory variability. Thus, sighing may shift respiratory control from the cortical and subcortical systems to the brainstem autonomic control system.
Breathing controlled by the limbic system
Emotion induces various physiological responses, including changes in heart rate, blood pressure, body temperature, and respiratory patterns, by activating the autonomic nervous system [124]. Among the changes in respiratory patterns, changes in the RR have been investigated extensively. Negative emotions, such as anxiety [125], fear [126], and sadness [127] increase the RR. Positive emotions, such as happiness, increase the RR [126, 127], while relief decreases the RR [125]. In addition to changes in the RR, each emotion appears to accompany a characteristic pattern of breathing that, to some extent, may overlap with the pattern from other emotions [128]. For example, fast, deep breaths are associated with excitement, while rapid, shallow breaths are associated with concentration, fear, and panic. Discriminant analyses indicate that four emotions (anger, fear, happiness, and sadness) can be adequately classified using heart rate variability, respiratory sinus arrhythmia, the mean RR, and respiratory variability, suggesting that distinct patterns of peripheral physiological activity are associated with different emotions [126].
The response to an emotion in the respiratory pattern is affected by the personality of the subject [129]. In subjects with high levels of anxiety, increases in the RR are more dominant responses to mental stress than changes in the VT [130], and changes in expiratory duration are dependent on anxiety scores [131]. Personality even affects respiratory parameters of subjects at rest [132]. The VT is smaller and the RR is higher in subjects with more anxiety and higher states on the State-Trait Anxiety Inventory [132].
The RR increases during the interval between alert presentation and actual stimulation, irrespective of changes in O2 consumption [133], which implies that breathing is controlled by the cortical or subcortical areas associated with emotion during the anticipatory anxiety period. Masaoka et al. [134] analyzed electroencephalogram (EEG) data during the anticipatory anxiety period and observed positive waves in cycle-triggered averaged EEG signals approximately 350 ms after the onset of inspiration, which are known as respiratory-related anxiety potentials (RAPs). A dipole tracing analysis based on a scalp-skull-brain head model identified the source location of the RAPs as the right temporal pole, while in the most anxious subject, it was the temporal pole and the amygdala [134]. A blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) study [135] demonstrated that the insula is essential for dyspnea perception. In addition, activation of the anterior cingulate cortex was correlated with the Breathlessness Catastrophizing Scale during dyspnea anticipation [136].
The respiration–emotion relationship is bidirectional. Deep and slow breathing (DSB) reduces anxiety and skin conductance levels in alcohol-dependent young adults [137]. Older adults also benefit from DSB in terms of vagal tone and anxiety [138]. The ameliorating effects of DSB on anxiety are believed to be mediated by reinforcement of the vagal tone, which balances sympathetic and parasympathetic activity [139]. Philippot and Blairy [128] tested whether respiratory changes affect emotions. When subjects mimicked a breathing pattern characteristic of joy, anger, fear, or sadness, the emotional state characterized by that breathing pattern was evoked, suggesting that alterations in breathing patterns can induce emotion. Masaoka et al. [140] showed that odors associated with autobiographical memories can trigger DSB and pleasant emotional experiences.
Control of respiration during cognition
Grassmann et al. [141] conducted a systematic review on respiration and cognitive loads. They found that in general, the cognitive load increases the RR and minute ventilation while not considerably impacting the VT. The end-tidal CO2 level decreased, which suggests that subjects were hyperventilated; however, oxygen consumption and CO2 release were also elevated. Changes in respiratory variability depend on the type of cognitive load [142]. Total variability in the RR decreases during sustained attention tasks, while during an arithmetic load, the autocorrelated variability decreases while the random variability increases. In addition, the frequency of sighing increased during sustained attention tasks but after arithmetic tasks, suggesting that the need for the respiratory control system to reset differs depending on the type of load [142]. Honma et al. [143] found that, compared with reading on paper, reading on a smartphone elicited fewer sighs and promoted brain overactivity in the prefrontal cortex, resulting in reduced comprehension.
There is growing evidence that human subjects can adjust their respiratory cycle to the onset of cognitive tasks, even if the tasks are not olfactory in nature [144,145,146,147,148]. Johannknecht et al. [146] found that subjects tend to align their respiratory cycle to the experimental paradigm, inhaling when the stimulus is presented and exhaling when submitting their responses. Respiratory timing affects cognitive task performance [145, 146, 148,149,150]. Zelano et al. [150] recorded intracranial EEG (iEEG) signals in patients with epilepsy and found that natural breathing synchronized electrical activity in the pyriform cortex, amygdala, and hippocampus. Fear discrimination and memory retrieval were enhanced during the inspiratory phase when the oscillatory power peaked. Cognitive performance was modulated during nose breathing but not during mouth breathing. Furthermore, Herrero et al. [151] demonstrated that coherence between the iEEG signal and breathing increased in the frontotemporal-insular network during volitionally paced breathing, whereas attention to breathing increased coherence in the anterior cingulate, premotor, insular, and hippocampal cortices. They proposed that breathing can organize neuronal oscillations throughout the brain [151]. In addition, Kluger et al. [152] applied phase-amplitude analysis to magnetoencephalography (MEG) data from quiet breathing humans and demonstrated the presence of respiration-mediated CFCs, termed respiration-modulated brain oscillations, across all major frequency bands in a widespread network of cortical and subcortical areas (Fig. 4b). Furthermore, they showed that occipital alpha power was coupled with respiration during near-threshold spatial detection tasks and that this respiration-alpha coupling was maximized with a respiratory phase lag of − 30°, indicating that the coupling occurs before behavioral consequences [147]. Time–frequency analyses revealed that compared with the alpha power prior to the presentation of undetected targets, the alpha power prior to the presentation of detected targets was significantly suppressed. Based on these results and the ‘active sensing’ concept for the functional role of olfaction [153], Kluger et al. [147] suggested that respiration actively adjusts the timing of sensory information sampling with transient oscillatory cycles of heightened cortical excitability to optimize performance. These respiratory acts are conceivably regulated by higher brain networks, albeit unconsciously; however, it has been shown that a subset of preBötC neurons regulates the balance between calm and arousal behaviors in a bottom-up fashion [154].
Respiratory variability in health and disease
Goldberger [4] presented an innovative concept: the output of a ‘healthy’ control system is not constant but instead fluctuates in a complex manner. Techniques that measure the complexity of the output have indicated that physiological control systems operate far from equilibrium and that maintaining constancy is not the goal. For example, fluctuations in the heart rate of healthy humans are chaotic [155] and multifractal, as they can be decomposed into multiple scaling regions [156], whereas fluctuations in the heart rates of patients with chronic heart failure (CHF) are less chaotic [155] and monofractal [156]. Aging also decreases the complexity of heart rate fluctuations, as quantified by SampEn [157]. Based on these observations, one might predict that the complex variability in respiratory signals would be greater in healthy individuals than in individuals with a disease. While this prediction is true in some cases, this is not always the case.
In human neonates with mild respiratory distress syndrome, the RR and VT exhibited increased complexity with increasing weight and gestational age; however, this complexity was observed only in terms of pattern matching-based entropies and not in the ApEn and SampEn, which are based on conditional probabilities [158]. This implies that respiratory fluctuations become increasingly complex with maturation. On the other hand, the LTC in the interbreath interval (IBI) time series decreases with age [159], similar to the changes in the complex variability of the heart rate [157]. Furthermore, sex affects changes in respiratory complex variability over the course of aging. The scaling exponents of the IBI time series are significantly lower (indicating decreased correlations) in healthy older males than in young males, young females, and older females [159]. The correlation dimensions of respiratory movement are lowest during slow wave sleep (stage IV) and highest during rapid eye movement (REM) sleep, with both correlated with the correlation dimension of the EEG signals [160]. STCs in the VT and minute ventilation, which may indicate the chemical control of breathing, have been observed during both non-REM and REM sleep, while LTCs have been observed only during REM sleep [12]. Exercise has opposite effects on the complex variability of respiration and heart rate, inducing a decrease in STCs in the IBI and an increase in LTCs of heart rate variability [161].
Various diseases affect respiratory variability. Breathing variability is remarkably augmented in patients with anxiety disorders, such as panic disorder [162,163,164,165]. Ventilatory complexity is also increased in patients with hyperventilation disorder; however, their respiratory control stability, which is assessed based on the loop gain, is not impaired [166]. In patients with breathing pattern disorders, a prevalent cause of exertional dyspnea, the ApEn of the VT and minute ventilation during the cardiopulmonary exercise test was significantly greater than that of controls [167].
The ventilatory flow of healthy, quietly breathing subjects exhibits nonlinear dynamics that are indicative of chaos [26]. In critically ill patients, switching from assist-controlled mechanical ventilation to inspiratory pressure support reduced the CV and eliminated nonlinear dynamics that are detectable using the noise titration technique [168]. Therefore, the chaotic feature of respiratory variability is neurogenic and is either intrinsic to the rCPG, a result of respiratory control processes driven by perturbations, or both, with little contribution from lung mechanics, if any [168]. However, changes in lung mechanics and gas exchange affect the gross variability [169, 170] and complexity of respiratory fluctuations [171,172,173]. The VT, RR, and minute ventilation are greater and the CV of the inspiratory time and minute ventilation are lower in patients with chronic obstructive pulmonary disease (COPD) than in healthy controls [170]. The random fraction of the breath variability is reduced, and the nonrandom, correlated fraction is greater in patients with restrictive lung disease than in healthy controls [169]. Interestingly, small variations in the average resting VT led to marked increases in dyspnea in these patients [169]. Compared with healthy subjects, the ApEn was significantly reduced in patients with asthma, which was correlated with the spirometric indices of airway obstruction [173]. Furthermore, the SampEn was significantly reduced in patients with COPD, which was also correlated with the spirometric indices of airway obstruction [171]. In addition, acute bronchodilation increased ventilatory complexity, as quantified by the noise titration technique during resting breathing in patients with stable COPD [172]. These changes in respiratory variability in patients with lung diseases may be a direct consequence of feedback from chemoreceptors and mechanoreceptors in the lung and airway. Alternatively, these changes could be due to alterations in autonomic function. Sympathetic nerve activity is increased in patients with COPD [174, 175], which is associated with morbidity and mortality [176]. Sympathetic neural overactivity may be a consequence of chronic hypoxia exposure [177]; however, slow breathing reduces elevated sympathetic activity in patients with COPD [178].
The most remarkable form of respiratory variability caused by instability in the chemical feedback loop is a type of periodic breathing known as Cheyne-Stokes respiration (CSR). CSR is a specific form of central sleep apnea characterized by the waxing and waning of the VT in 50 ~ 90 s intervals [179]. CSR has been observed in patients with CHF, particularly during stages 1 and 2 of non-REM sleep [180]. CSR has also been observed during the day and is more closely correlated with the severity of CHF [181]. The mechanism underlying CSR is thought to be a failure in chemical feedback control [179, 182, 183]. In control theory, the stability of a feedback system is defined by the controller gain, the plant gain, and the loop gain (LG). The controller gain is a measure of how much the controller responds to a given change in blood gas tension (\(\Delta {\dot{V}}_{E}\)/\(\Delta \mathrm{Pa{CO}}_{2}\)), and the plant gain is a measure of how much the blood gas tension changes for a given change in ventilation (\(\Delta {\mathrm{PaCO}}_{2}\)/\(\Delta {\dot{V}}_{E}\)). The LG, which is the product of the controller gain and the plant gain, represents the ratio of the ventilatory response to the ventilatory disturbance. An LG of less than 1 indicates stable breathing, whereas an LG of greater than 1 in combination with a prolonged circulatory delay results in periodic breathing. In CHF patients, a decrease in cardiac output can lead to prolonged circulatory delays and mild hypoxemia, increasing controller gain. These factors, in combination with an increase in plant gain due to sleep hypoventilation and the subsequent elevation of PaCO2, destabilize chemical feedback control, resulting in CSR [179]. Moreover, CHF patients tend to hyperventilate and become hypocapnic during wakefulness. Subsequently, the withdrawal of the wakefulness stimulus upon sleep leads to apnea [184, 185]. Nasal continuous positive airway pressure [186, 187] and inhalation of 3% CO2 [188] have been shown to ameliorate CSR by reducing plant gain.
Obstructive sleep apnea (OSA) is a common breathing disorder that involves periodic breathing with repetitive narrowing and closing of the upper airway during sleep [189]. The primary cause of OSA is an anatomically collapsible upper airway; however, additional nonanatomical factors, such as inadequate responsiveness of the upper airway dilator muscles during sleep, waking prematurely due to airway narrowing, and a high LG, characterize different phenotypes of OSA [190]. During apnea, both the plant gain and the controller gain increase; thus, the increased LG at the end of the obstruction is not the cause but the result of the obstructive event [191]. The chemical LG measured while the upper airway is stable is moderately elevated in some OSA patients; however, the increase is insufficient for causing instability in the absence of a collapsible upper airway [192, 193].
The breathing pattern of CHF patients is typically characterized by unstable respiration, such as rapid, irregular, and nonperiodic respiration with transient sighing or apnea, rather than CSR [194]. Respiratory instability is unlikely to be related to the negative feedback system of chemical respiratory control; rather, it might be caused by the stimulation of afferent vagal nerve endings due to lung edema [195, 196]. Asanoi et al. [195] developed a quantitative measure of respiratory instability (RSI) based on the frequency distribution of respiratory spectral components and the very low-frequency components. They found that patients who died from cardiac causes had a lower RSI and suggested that an RSI < 20 predicts a higher probability of subsequent all-cause and cardiovascular death. Okamoto et al. [197] analyzed stable airflow data before the onset of sleep to quantify breathing irregularities using the ShEn in patients with relatively mild CHF, ischemic disease, or atrial fibrillation. They found that the ShEn of the airflow signals in these patients was significantly greater than the ShEn of patients without heart disease.
Future directions in translational sciences
The idea of extracting hidden information from respiratory signals and utilizing these data in clinical practice and daily life is attractive. This process could be easily carried out in intensive care units, where continuous monitoring of breathing is the standard protocol. In intensive care units, respiratory variability may have predictive value for successful weaning from mechanical ventilation [198]. For example, Wysocki et al. [199] showed that the reduced CVs of the TV/inspiratory time and inspiratory time/respiratory period can be used to predict successful weaning cases. Additionally, El-Khatib et al. [200] showed that spontaneous breathing patterns during minimal mechanical ventilatory support are more chaotic in patients who failed extubation trials than in patients who passed them. Similarly, Engoren et al. [201] showed that the RR and ApEn of the VT increase upon spontaneous ventilation in weaning trials for patients who require mechanical ventilation. Nonlinear dynamics analyses can also be used to diagnose specific diseases. Miyata et al. [202] showed that the correlation dimension of chest movement with a brief period during wakefulness may be a useful index for identifying patients with OSA. Raoufy et al. [203] showed that nonlinear analyses (LLE, LTC, and SampEn) of breathing patterns have diagnostic value in asthma and can be used to differentiate uncontrolled and controlled asthma as well as nonatopic and atopic asthma using receiver operating characteristic (ROC) curve analysis.
The gold standard technique for staging sleep is polysomnography; however, this method requires expensive equipment with constrained sensors for recording and human resources for analysis. Therefore, there is a need for an automated sleep staging system that ideally uses an inexpensive, wearable or noncontact sensor. Breathing patterns in infants are considerably different between active sleep (equivalent to adult REM sleep) and quiet sleep (equivalent to adult non-REM sleep) [204]. Haddad et al. [205] reported that the CV of the IBI can be used to adequately distinguish active and quiet sleep stages in newborn infants. Harper et al. [206] applied machine learning techniques to identify sleep stages in newborn infants according to cardiorespiratory variables. Terrill et al. attempted to use nonlinear analyses of respiratory variability for sleep staging [207]. They showed that features extracted from recurrence plots of the IBI using recurrence quantification analysis [17, 208] can be used to classify sleep stages in infants. Recently, machine learning techniques have been applied not only to identify sleep stages but also to detect respiratory events (apnea, hypopnea, and CSR) during sleep [209,210,211].
Another promising research direction is an application toward emotion recognition. Emotion is tightly coupled with physiological changes that are specific to each emotion (see "Breathing controlled by the limbic system" section) [126, 212]. Advances in wearable sensors for measuring physiological signals and machine learning techniques have allowed e-health research to focus on emotion recognition [213,214,215,216]. Emotion recognition technology is expected to be applied in various fields, such as mental health conditioning, man–machine interfaces, marketing, and education [217,218,219,220,221]. Studies on determining human emotions in the engineering field generally use a two-dimensional model known as Russell’s circumplex model of affect for emotion classification [222, 223] since it can easily be used with classification algorithms. The circumplex model assumes that all affective states arise from two fundamental neurophysiological systems: one related to valence (a pleasure–displeasure continuum) and another related to arousal, or alertness. Although only a few studies have used respiration for emotion recognition to date [214], a study based on deep learning algorithms applied to the dataset DEAP [224] showed valence and arousal accuracies of 73% and 81%, respectively [221].
Conclusions
Respiratory variability contains a veritable treasure trove of hidden information. The elucidation of the mechanisms underlying this variability is undoubtedly important; however, deep learning techniques and information theory-based quantification of complex variability have allowed us to use this variability for inference and decision making without knowing its precise sources and mechanisms. On the other hand, since the structure of the model is not expected to mimic the actual system in conventional deep learning techniques, the techniques cannot be applied to elucidate the sources and mechanisms of respiratory variability. Rather, these techniques in combination with smart sensors and devices should be used to improve the health and quality of life of everyone.
Availability of data and materials
Not applicable.
Abbreviations
- CV:
-
Coefficient of variation
- RMSSD:
-
Root mean square successive difference
- AR:
-
Autocorrelation
- LTC:
-
Long-term correlation
- STC:
-
Short-term correlation
- LLE:
-
Largest Lyapunov exponent
- ShEn:
-
Shannon entropy
- ApEn:
-
Approximate entropy
- SampEn:
-
Sample entropy
- NL:
-
Noise limit
- rCPG:
-
Respiratory central pattern generator
- preBötC:
-
Pre-Bötzinger complex
- pFRG:
-
Parafacial respiratory group
- RTN:
-
Retrotrapezoid nucleus
- KF:
-
Kölliker-Fuse
- CFC:
-
Cross-frequency coupling
- CSN:
-
Carotid sinus nerve
- STP:
-
Short-term potentiation
- STD:
-
Short-term depression
- LTD:
-
Long-term depression
- PHRA:
-
Post-hypoxic persistent respiratory augmentation
- HVD:
-
Hypoxic ventilatory depression
- SLN:
-
Superior laryngeal nerve
- VNS:
-
Vagal nerve stimulation
- RR:
-
Respiratory rate
- VT:
-
Tidal volume
- EEG:
-
Electroencephalogram
- RAP:
-
Respiratory-related anxiety potential
- BOLD:
-
Blood oxygen level dependent
- fMRI:
-
Functional magnetic resonance imaging
- DSB:
-
Deep and slow breathing
- iEEG:
-
Intracranial electroencephalogram
- MEG:
-
Magnetoencephalography
- CHF:
-
Chronic heart failure
- IBI:
-
Interbreath interval
- REM:
-
Rapid eye movement
- COPD:
-
Chronic obstructive pulmonary disease
- CSR:
-
Cheyne-Stokes respiration
- LG:
-
Loop gain
- OSA:
-
Obstructive sleep apnea
- RSI:
-
Quantitative measure of respiratory instability
- ROC:
-
Receiver operating characteristic
References
Bruce EN (1996) Temporal variations in the pattern of breathing. J Appl Physiol 80:1079–1087
Goldberger AL, Moody GB, Costa MD. Variability vs. complexity. https://archive.physionet.org/tutorials/cv/. Accessed 8 June 2022
Frey U, Maksym G, Suki B (2011) Temporal complexity in clinical manifestations of lung disease. J Appl Physiol 110:1723–1731
Goldberger AL, Amaral LA, Hausdorff JM, Ivanov P, Peng CK, Stanley HE (2002) Fractal dynamics in physiology: alterations with disease and aging. Proc Natl Acad Sci 99(Suppl 1):2466–2472
Macklem PT (2008) Emergent phenomena and the secrets of life. J Appl Physiol 104:1844–1846
Suki B, Bates JH, Frey U (2011) Complexity and emergent phenomena. Compr Physiol 1:995–1029
Thamrin C, Frey U, Kaminsky DA, Reddel HK, Seely AJ, Suki B, Sterk PJ (2016) Systems biology and clinical practice in respiratory medicine. The Twain shall meet. Am J Respir Crit Care Med 194:1053–1061
Minarini G (2020) Root mean square of the successive differences as marker of the parasympathetic system and difference in the outcome after ANS stimulation. In: Aslanidis T (ed) Autonomic nervous system monitoring - heart rate variability. IntechOpen, London
van den Bosch OFC, Alvarez-Jimenez R, de Grooth HJ, Girbes ARJ, Loer SA (2021) Breathing variability-implications for anaesthesiology and intensive care. Crit Care 25:280
Chen Z, Ivanov PC, Hu K, Stanley HE (2002) Effect of nonstationarities on detrended fluctuation analysis. Phys Rev E 65:041107
Morariu VV, Buimaga-Iarinca L, Vamos C, Soltuz SM (2007) Detrended fluctuation analysis of autoregressive processes. Fluct Noise Lett 7:L249–L255
Rostig S, Kantelhardt JW, Penzel T, Cassel W, Peter JH, Vogelmeier C, Becker HF, Jerrentrup A (2005) Nonrandom variability of respiration during sleep in healthy humans. Sleep 28:411–417
Wolf A (1986) Quantifying chaos with Lyapunov exponents. In: Holden AV (ed) Chaos. Princeton University Press, New Jersey, pp 273–290
Deyle ER, Sugihara G (2011) Generalized theorems for nonlinear state space reconstruction. PLoS ONE 6:e18295
Takens F. (1981) Detecting strange attractors in turbulence. In: Proc. Warwick Symp. 1980 Dynamical Systems and Turbulence, eds. Rand DA and Young BS (Springer, Berlin). Lect Notes Math 898:366–81.
Eckmann JP, Kamphorst SO, Ruelle D (1987) Recurrence plots of dynamical systems. Europhys Lett 5:973–977
Webber CL, Zbilut JP (1994) Dynamical assessment of physiological systems and states using recurrence plot strategies. J Appl Physiol 76:965–973
Grassberger P, Procaccia A (1983) Measuring the strangeness of strange attractors. Physica D 9:189–208
Pritchard WS, Duke DW (1995) Measuring “chaos” in the brain: a tutorial review of EEG dimension estimation. Brain Cogn 27:353–397
Pincus SM (1991) Approximate entropy as a measure of system complexity. Proc Natl Acad Sci USA 88:2297–2301
Richman JS, Moorman JR (2000) Physiological time-series analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol 278:H2039–H2049
Stone L (1992) Coloured noise or low-dimensional chaos? Proc Biol Sci 250:77–81
Lancaster G, Iatsenko D, Pidde A, Ticcinelli V, Stefanovska A (2018) Surrogate data for hypothesis testing of physical systems. Phys Rep 748:1–60
Theiler J, Eubank S, Longtin A, Galdrikian B, Farmer JD (1992) Testing for nonlinearity in time series: the method of surrogate data. Physica D 58:77–94
Poon CS, Barahona M (2001) Titration of chaos with added noise. Proc Natl Acad Sci USA 98:7107–7112
Wysocki M, Fiamma MN, Straus C, Poon CS, Similowski T (2006) Chaotic dynamics of resting ventilatory flow in humans assessed through noise titration. Respir Physiol Neurobiol 153:54–65
Gao JB, Hu J, Mao X, Tung WW (2012) Detecting low-dimensional chaos by the “noise titration” technique: possible problems and remedies. Chaos Soliton Fract 45:213–223
Cherniack NS (1987) Potential role of optimization in alveolar hypoventilation and respiratory instability. In: von Euler C, Lagercrantz H (eds) Neurobiology of the control of breathing. Raven Press, New York, pp 45–50
Poon CS (1987) Ventilatory control in hypercapnia and exercise—optimization hypothesis. J Appl Physiol 62:2447–2459
Jakus J, Tomori Z, Stransky A (1985) Activity of bulbar respiratory neurones during cough and other respiratory tract reflexes in cats. Physiol Bohemoslov 34:127–136
Oku Y, Tanaka I, Ezure K (1994) Activity of bulbar respiratory neurons during fictive coughing and swallowing in the decerebrate cat. J Physiol-London 480:309–324
Baertsch NA, Severs LJ, Anderson TM, Ramirez JM (2019) A spatially dynamic network underlies the generation of inspiratory behaviors. Proc Natl Acad Sci USA 116:7493–7502
Lieske SP, Thoby-Brisson M, Telgkamp P, Ramirez JM (2000) Reconfiguration of the neural network controlling multiple breathing patterns: eupnea, sighs and gasps [see comment]. Nat Neurosci 3:600–607
Lindsey BG, Rybak IA, Smith JC (2012) Computational models and emergent properties of respiratory neural networks. Compr Physiol 2:1619–1670
Butler JE (2007) Drive to the human respiratory muscles. Respir Physiol Neurobiol 159:115–126
Rikard-Bell GC, Bystrzycka EK, Nail BS (1985) Cells of origin of corticospinal projections to phrenic and thoracic respiratory motoneurones in the cat as shown by retrograde transport of HRP. Brain Res Bull 14:39–47
Trevizan-Bau P, Dhingra RR, Furuya WI, Stanic D, Mazzone SB, Dutschmann M (2021) Forebrain projection neurons target functionally diverse respiratory control areas in the midbrain, pons, and medulla oblongata. J Comp Neurol 529:2243–2264
Yang CF, Kim EJ, Callaway EM, Feldman JL (2020) Monosynaptic projections to excitatory and inhibitory preBötzinger complex neurons. Front Neuroanat 14:58
Benchetrit G, Bertrand F (1975) A short-term memory in the respiratory centres: statistical analysis. Respir Physiol 23:147–158
Del Negro CA, Funk GD, Feldman JL (2018) Breathing matters. Nat Rev Neurosci 19:351–367
Richter DW, Smith JC (2014) Respiratory rhythm generation in vivo. Physiology 29:58–71
Smith JC, Ellenberger HH, Ballanyi K, Richter DW, Feldman JL (1991) Pre-Bötzinger Complex—a brain-stem region that may generate respiratory rhythm in mammals. Science 254:726–729
Feldman JL, Del Negro CA, Gray PA (2013) Understanding the rhythm of breathing: so near, yet so far. Annu Rev Physiol 75:423–452
Takakura AC, Malheiros-Lima MR, Moreira TS (2021) Excitatory and inhibitory modulation of parafacial respiratory neurons in the control of active expiration. Respir Physiol Neurobiol 289:103657
Ikeda K, Kawakami K, Onimaru H, Okada Y, Yokota S, Koshiya N, Oku Y, Iizuka M, Koizumi H (2017) The respiratory control mechanisms in the brainstem and spinal cord: integrative views of the neuroanatomy and neurophysiology. J Physiol Sci 67:45–62
Koshiya N, Smith JC (1999) Neuronal pacemaker for breathing visualized in vitro. Nature 400:360–363
Del Negro CA, Wilson CG, Butera RJ, Rigatto H, Smith JC (2002) Periodicity, mixed-mode oscillations, and quasiperiodicity in a rhythm-generating neural network. Biophys J 82:206–214
Koshiya N, Oku Y, Yokota S, Oyamada Y, Yasui Y, Okada Y (2014) Anatomical and functional pathways of rhythmogenic inspiratory premotor information flow originating in the pre-Bötzinger complex in the rat medulla. Neuroscience 268:194–211
Carroll MS, Ramirez JM (2013) Cycle-by-cycle assembly of respiratory network activity is dynamic and stochastic. J Neurophysiol 109:296–305
Kuwana S, Tsunekawa N, Yanagawa Y, Okada Y, Kuribayashi J, Obata K (2006) Electrophysiological and morphological characteristics of GABAergic respiratory neurons in the mouse pre-Bötzinger complex. Eur J Neurosci 23:667–674
Winter SM, Fresemann J, Schnell C, Oku Y, Hirrlinger J, Hulsmann S (2009) Glycinergic interneurons are functionally integrated into the inspiratory network of mouse medullary slices. Pflug Arch Eur J Phy 458:459–469
Lal A, Oku Y, Someya H, Miwakeichi F, Tamura Y (2016) Emergent network topology within the respiratory rhythm-generating kernel evolved in silico. PLoS ONE 11:e0154049
Onimaru H, Homma I (2003) A novel functional neuron group for respiratory rhythm generation in the ventral medulla. J Neurosci 23:1478–1486
Suzue T (1984) Respiratory rhythm generation in the in vitro brain stem-spinal cord preparation of the neonatal rat. J Physiol 354:173–183
Guyenet PG, Bayliss DA (2015) Neural control of breathing and CO2 homeostasis. Neuron 87:946–961
Guyenet PG, Stornetta RL, Souza G, Abbott SBG, Shi Y, Bayliss DA (2019) The retrotrapezoid nucleus: central chemoreceptor and regulator of breathing automaticity. Trends Neurosci 42:807–824
Zoccal DB, Silva JN, Barnett WH, Lemes EV, Falquetto B, Colombari E, Molkov YI, Moreira TS, Takakura AC (2018) Interaction between the retrotrapezoid nucleus and the parafacial respiratory group to regulate active expiration and sympathetic activity in rats. Am J Physiol Lung Cell Mol Physiol 315:L891–L909
Janczewski WA, Feldman JL (2006) Distinct rhythm generators for inspiration and expiration in the juvenile rat. J Physiol-London 570:407–420
Mellen NM, Janczewski WA, Bocchiaro CM, Feldman JL (2003) Opioid-induced quantal slowing reveals dual networks for respiratory rhythm generation. Neuron 37:821–826
Wittmeier S, Song G, Duffin J, Poon CS (2008) Pacemakers handshake synchronization mechanism of mammalian respiratory rhythmogenesis. Proc Natl Acad Sci USA 105:18000–18005
Lal A, Oku Y, Hulsmann S, Okada Y, Miwakeichi F, Kawai S, Tamura Y, Ishiguro M (2011) Dual oscillator model of the respiratory neuronal network generating quantal slowing of respiratory rhythm. J Comput Neurosci 30:225–240
Richter DW, Spyer KM (2001) Studying rhythmogenesis of breathing: comparison of in vivo and in vitro models. Trends Neurosci 24:464–472
Smith JC, Abdala APL, Koizumi H, Rybak IA, Paton JFR (2007) Spatial and functional architecture of the mammalian brain stem respiratory network: a hierarchy of three oscillatory mechanisms. J Neurophysiol 98:3370–3387
Dutschmann M, Dick TE (2012) Pontine Mechanisms of Respiratory Control. Compr Physiol 2:2443–2469
Fung ML, Wang W, St John WM (1994) Involvement of pontile NMDA receptors in inspiratory termination in rat. Respir Physiol 96:177–188
Ling L, Karius DR, Speck DF (1985) (1994) Role of N-methyl-d-aspartate receptors in the pontine pneumotaxic mechanism in the cat. J Appl Physiol 76:1138–1143
Oku Y, Dick TE (1992) Phase resetting of the respiratory cycle before and after unilateral pontine lesion in cat. J Appl Physiol 72:721–730
Yu H, Dhingra RR, Dick TE, Galan RF (2017) Effects of ion channel noise on neural circuits: an application to the respiratory pattern generator to investigate breathing variability. J Neurophysiol 117:230–242
Van Horn MR, Benfey NJ, Shikany C, Severs LJ, Deemyad T (2021) Neuron-astrocyte networking: astrocytes orchestrate and respond to changes in neuronal network activity across brain states and behaviors. J Neurophysiol 126:627–636
Turk AZ, Bishop M, Adeck A, SheikhBahaei S (2022) Astrocytic modulation of central pattern generating motor circuits. Glia 70(8):1506–1519
Hulsmann S, Oku Y, Zhang W, Richter DW (2000) Metabolic coupling between glia and neurons is necessary for maintaining respiratory activity in transverse medullary slices of neonatal mouse. Eur J Neurosci 12:856–862
Sheikhbahaei S, Turovsky EA, Hosford PS, Hadjihambi A, Theparambil SM, Liu B, Marina N, Teschemacher AG, Kasparov S, Smith JC, Gourine AV (2018) Astrocytes modulate brainstem respiratory rhythm-generating circuits and determine exercise capacity. Nat Commun 9:370
Schnell C, Fresemann J, Hulsmann S (2011) Determinants of functional coupling between astrocytes and respiratory neurons in the pre-Bötzinger complex. PLoS ONE 6:e26309
Okada Y, Sasaki T, Oku Y, Takahashi N, Seki M, Ujita S, Tanaka KF, Matsuki N, Ikegaya Y (2012) Preinspiratory calcium rise in putative pre-Bötzinger complex astrocytes. J Physiol-London 590:4933–4944
Forsberg D, Herlenius E (2019) Astrocyte networks modulate respiration—sniffing glue. Respir Physiol Neurobiol 265:3–8
Forsberg D, Ringstedt T, Herlenius E (2017) Astrocytes release prostaglandin E2 to modify respiratory network activity. Elife 6:e29566
Smedler E, Malmersjo S, Uhlen P (2014) Network analysis of time-lapse microscopy recordings. Front Neural Circuit 8:111
Watts DJ, Strogatz SH (1998) Collective dynamics of “small-world” networks. Nature 393:440–442
Oku Y, Fresemann J, Miwakeichi F, Hulsmann S (2016) Respiratory calcium fluctuations in low-frequency oscillating astrocytes in the pre-Bötzinger complex. Respir Physiol Neurobiol 226:11–17
Sakaguchi H, Okita T (2016) Cooperative dynamics in coupled systems of fast and slow phase oscillators. Phys Rev E 93:022212
Strogatz SH (2000) From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D 143:1–20
Canolty RT, Knight RT (2010) The functional role of cross-frequency coupling. Trends Cogn Sci 14:506–515
Jensen O, Colgin LL (2007) Cross-frequency coupling between neuronal oscillations. Trends Cogn Sci 11:267–269
Marwan N, Kurths J (2002) Nonlinear analysis of bivariate data with cross recurrence plots. Phys Lett A 302:299–307
Marwan N, Thiel M, Nowaczyk NR (2002) Cross recurrence plot based synchronization of time series. Nonlinear Proc Geoph 9:325–331
Wagner PG, Eldridge FL (1991) Development of short-term potentiation of respiration. Respir Physiol 83:129–139
Mateika JH, Syed Z (2013) Intermittent hypoxia, respiratory plasticity and sleep apnea in humans: present knowledge and future investigations. Respir Physiol Neurobiol 188:289–300
Mitchell GS, Baker TL, Nanda SA, Fuller DD, Zabka AG, Hodgeman BA, Bavis RW, Mack KJ, Olson EB Jr (2001) Invited review: intermittent hypoxia and respiratory plasticity. J Appl Physiol 90:2466–2475
Fukushi I, Takeda K, Pokorski M, Kono Y, Yoshizawa M, Hasebe Y, Nakao A, Mori Y, Onimaru H, Okada Y (2021) Activation of astrocytes in the persistence of post-hypoxic respiratory augmentation. Front Physiol 12:757731
Powell FL, Milsom WK, Mitchell GS (1998) Time domains of the hypoxic ventilatory response. Respir Physiol 112:123–134
Hayashi F, Coles SK, Bach KB, Mitchell GS, McCrimmon DR (1993) Time-dependent phrenic nerve responses to carotid afferent activation: intact vs. decerebellate rats. Am J Physiol 265:R811–R819
Coles SK, Dick TE (1996) Neurones in the ventrolateral pons are required for post-hypoxic frequency decline in rats. J Physiol 497(Pt 1):79–94
Rajani V, Zhang Y, Jalubula V, Rancic V, SheikhBahaei S, Zwicker JD, Pagliardini S, Dickson CT, Ballanyi K, Kasparov S, Gourine AV, Funk GD (2018) Release of ATP by pre-Bötzinger complex astrocytes contributes to the hypoxic ventilatory response via a Ca2+-dependent P2Y(1) receptor mechanism. J Physiol-London 596:3245–3269
Khatib MF, Oku Y, Bruce EN (1991) Contribution of chemical feedback loops to breath-to-breath variability of tidal volume. Respir Physiol 83:115–127
Berssenbrugge A, Dempsey J, Iber C, Skatrud J, Wilson P (1983) Mechanisms of hypoxia-induced periodic breathing during sleep in humans. J Physiol 343:507–524
Jubran A, Tobin MJ (2000) Effect of isocapnic hypoxia on variational activity of breathing. Am J Respir Crit Care Med 162:1202–1209
Fiamma MN, Straus C, Thibault S, Wysocki M, Baconnier P, Similowski T (2007) Effects of hypercapnia and hypocapnia on ventilatory variability and the chaotic dynamics of ventilatory flow in humans. Am J Physiol Regul Integr Comp Physiol 292:R1985–R1993
Song G, Poon CS (2004) Functional and structural models of pontine modulation of mechanoreceptor and chemoreceptor reflexes. Respir Physiol Neurobiol 143:281–292
Glass L, Mackey M (1988) From clocks to chaos. Princeton University Press, New Jersey
Paydarfar D, Eldridge FL (1987) Phase resetting and dysrhythmic responses of the respiratory oscillator. Am J Physiol 252:R55-62
Lewis J, Bachoo M, Polosa C, Glass L (1990) The effects of superior laryngeal nerve stimulation on the respiratory rhythm: phase-resetting and aftereffects. Brain Res 517:44–50
Sammon M, Romaniuk JR, Bruce EN (1993) Bifurcations of the respiratory pattern produced with phasic vagal-stimulation in the rat. J Appl Physiol 75:912–926
Dhingra RR, Jacono FJ, Fishman M, Loparo KA, Rybak IA, Dick TE (2011) Vagal-dependent nonlinear variability in the respiratory pattern of anesthetized, spontaneously breathing rats. J Appl Physiol 111:272–284
Dhingra RR, Dutschmann M, Galan RF, Dick TE (2017) Kölliker-Fuse nuclei regulate respiratory rhythm variability via a gain-control mechanism. Am J Physiol Regul Integr Comp Physiol 312:R172–R188
Oku Y (2020) Coordination of swallowing and breathing: how is the respiratory control system connected to the swallowing system? In: Yamaguchi K (ed) Structure-function relationships in various respiratory systems. Springer Nature, Singapore, pp 37–52
Paydarfar D, Gilbert RJ, Poppel CS, Nassab PF (1995) Respiratory phase resetting and air-flow changes induced by swallowing in humans. J Physiol-London 483:273–288
Yagi N, Oku Y, Nagami S, Yamagata Y, Kayashita J, Ishikawa A, Domen K, Takahashi R (2017) Inappropriate timing of swallow in the respiratory cycle causes breathing-swallowing discoordination. Front Physiol 8:676
Saito Y, Ezure K, Tanaka I, Osawa M (2003) Activity of neurons in ventrolateral respiratory groups during swallowing in decerebrate rats. Brain Dev 25:338–345
Sugiyama Y, Shiba K, Mukudai S, Umezaki T, Hisa Y (2014) Activity of respiratory neurons in the rostral medulla during vocalization, swallowing, and coughing in guinea pigs. Neurosci Res 80:17–31
Bonis JM, Neumueller SE, Krause KL, Pan LG, Hodges MR, Forster HV (2013) Contributions of the Kölliker-Fuse nucleus to coordination of breathing and swallowing. Respir Physiol Neurobiol 189:10–21
Bonis JM, Neumueller SE, Marshall BD, Krause KL, Qian B, Pan LG, Hodges MR, Forster HV (2011) The effects of lesions in the dorsolateral pons on the coordination of swallowing and breathing in awake goats. Respir Physiol Neurobiol 175:272–282
Bautista TG, Dutschmann M (2014) Ponto-medullary nuclei involved in the generation of sequential pharyngeal swallowing and concomitant protective laryngeal adduction in situ. J Physiol 592:2605–2623
Horton KK, Segers LS, Nuding SC, O’Connor R, Alencar PA, Davenport PW, Bolser DC, Pitts T, Lindsey BG, Morris KF, Gestreau C (2018) Central respiration and mechanical ventilation in the gating of swallow with breathing. Front Physiol 9:785
Li P, Janczewski WA, Yackle K, Kam K, Pagliardini S, Krasnow MA, Feldman JL (2016) The peptidergic control circuit for sighing. Nature 530:293
Davis GM, Moscato J (1994) Changes in lung mechanics following sighs in premature newborns without lung disease. Pediatr Pulmonol 17:26–30
Ferris BG Jr, Pollard DS (1960) Effect of deep and quiet breathing on pulmonary compliance in man. J Clin Invest 39:143–149
Cherniack NS, von Euler C, Glogowska M, Homma I (1981) Characteristics and rate of occurrence of spontaneous and provoked augmented breaths. Acta Physiol Scand 111:349–360
Ramirez JM (2014) The Integrative Role of the Sigh in Psychology, Physiology, Pathology, and Neurobiology. Prog Brain Res 209:91–129
Vlemincx E, Taelman J, Van Diest I, Van den Bergh O (2010) Take a deep breath: The relief effect of spontaneous and instructed sighs. Physiol Behav 101:67–73
Vlemincx E, Abelson JL, Lehrer PM, Davenport PW, Van Diest I, Van den Bergh O (2013) Respiratory variability and sighing: a psychophysiological reset model. Biol Psychol 93:24–32
Vlemincx E, Van Diest I, Lehrer PM, Aubert AE, Van den Bergh O (2010) Respiratory variability preceding and following sighs: a resetter hypothesis. Biol Psychol 84:82–87
Vlemincx E, Van Diest I, Van den Bergh O (2016) A sigh of relief or a sigh to relieve: the psychological and physiological relief effect of deep breaths. Physiol Behav 165:127–135
Moss F, Ward LM, Sannita WG (2004) Stochastic resonance and sensory information processing: a tutorial and review of application. Clin Neurophysiol 115:267–281
Kreibig SD (2010) Autonomic nervous system activity in emotion: a review. Biol Psychol 84:394–421
Blechert J, Lajtman M, Michael T, Margraf J, Wilhelm FH (2006) Identifying anxiety states using broad sampling and advanced processing of peripheral physiological information. Biomed Sci Instrum 42:136–141
Rainville P, Bechara A, Naqvi N, Damasio AR (2006) Basic emotions are associated with distinct patterns of cardiorespiratory activity. Int J Psychophysiol 61:5–18
Rottenberg J, Salomon K, Gross JJ, Gotlib IH (2005) Vagal withdrawal to a sad film predicts subsequent recovery from depression. Psychophysiology 42:277–281
Philippot P, Chapelle G, Blairy S (2002) Respiratory feedback in the generation of emotion. Cogn Emot 16:605–627
Homma I, Masaoka Y (2008) Breathing rhythms and emotions. Exp Physiol 93:1011–1021
Masaoka Y, Homma I (1997) Anxiety and respiratory patterns: their relationship during mental stress and physical load. Int J Psychophysiol 27:153–159
Masaoka Y, Homma I (1999) Expiratory time determined by individual anxiety levels in humans. J Appl Physiol 86:1329–1336
Kato A, Takahashi K, Homma I (2018) Relationships between trait and respiratory parameters during quiet breathing in normal subjects. J Physiol Sci 68:369–376
Masaoka Y, Homma I (2001) The effect of anticipatory anxiety on breathing and metabolism in humans. Respir Physiol 128:171–177
Masaoka Y, Homma I (2000) The source generator of respiratory-related anxiety potential in the human brain. Neurosci Lett 283:21–24
Evans KC, Banzett RB, Adams L, McKay L, Frackowiak RSJ, Corfield DR (2002) BOLD fMRI identifies limbic, paralimbic, and cerebellar activation during air hunger. J Neurophysiol 88:1500–1511
Stoeckel MC, Esser RW, Gamer M, Buchel C, von Leupoldt A (2018) Dyspnea catastrophizing and neural activations during the anticipation and perception of dyspnea. Psychophysiology 55:e13004
Clark ME, Hirschman R (1990) Effects of paced respiration on anxiety reduction in a clinical population. Biofeedback Self Regul 15:273–284
Magnon V, Dutheil F, Vallet GT (2021) Benefits from one session of deep and slow breathing on vagal tone and anxiety in young and older adults. Sci Rep 11:19267
Porges SW (2007) The polyvagal perspective. Biol Psychol 74:116–143
Masaoka Y, Sugiyama H, Katayama A, Kashiwagi M, Homma I (2012) Slow breathing and emotions associated with odor-induced autobiographical memories. Chem Senses 37:379–388
Grassmann M, Vlemincx E, von Leupoldt A, Mittelstadt JM, Van den Bergh O (2016) Respiratory changes in response to cognitive load: a systematic review. Neural Plast 2016:8146809
Vlemincx E, Taelman J, De Peuter S, Van Diest I, Van den Bergh O (2011) Sigh rate and respiratory variability during mental load and sustained attention. Psychophysiology 48:117–120
Honma M, Masaoka Y, Iizuka N, Wada S, Kamimura S, Yoshikawa A, Moriya R, Kamijo S, Izumizaki M (2022) Reading on a smartphone affects sigh generation, brain activity, and comprehension. Sci Rep-Uk 12:1589
Boyadzhieva A, Kayhan E (2021) Keeping the breath in mind: respiration, neural oscillations, and the free energy principle. Front Neurosci 15:647579
Grund M, Al E, Pabst M, Dabbagh A, Stephani T, Nierhaus T, Gaebler M, Villringer A (2022) Respiration, heartbeat, and conscious tactile perception. J Neurosci 42:643–656
Johannknecht M, Kayser C (2022) The influence of the respiratory cycle on reaction times in sensory-cognitive paradigms. Sci Rep 12:2586
Kluger DS, Balestrieri E, Busch NA, Gross J (2021) Respiration aligns perception with neural excitability. Elife 10:e70907
Perl O, Ravia A, Rubinson M, Eisen A, Soroka T, Mor N, Secundo L, Sobel N (2019) Human non-olfactory cognition phase-locked with inhalation. Nat Hum Behav 3:501–512
Nakamura NH, Fukunaga M, Oku Y (2018) Respiratory modulation of cognitive performance during the retrieval process. PLoS ONE 13:e0204021
Zelano C, Jiang HD, Zhou GY, Arora N, Schuele S, Rosenow J, Gottfried JA (2016) Nasal respiration entrains human limbic oscillations and modulates cognitive function. J Neurosci 36:12448–12467
Herrero JL, Khuvis S, Yeagle E, Cerf M, Mehta AD (2018) Breathing above the brain stem: volitional control and attentional modulation in humans. J Neurophysiol 119:145–159
Kluger D, Gross J (2021) Respiration modulates oscillatory neural network activity at rest. PLoS Biol 19:e3001457
Wachowiak M (2011) All in a sniff: olfaction as a model for active sensing. Neuron 71:962–973
Yackle K, Schwarz LA, Kam K, Sorokin JM, Huguenard JR, Feldman JL, Luo L, Krasnow MA (2017) Breathing control center neurons that promote arousal in mice. Science 355:1411–1415
Poon CS, Merrill CK (1997) Decrease of cardiac chaos in congestive heart failure. Nature 389:492–495
Ivanov PC, Amaral LA, Goldberger AL, Havlin S, Rosenblum MG, Struzik ZR, Stanley HE (1999) Multifractality in human heartbeat dynamics. Nature 399:461–465
Kaplan DT, Furman MI, Pincus SM, Ryan SM, Lipsitz LA, Goldberger AL (1991) Aging and the complexity of cardiovascular dynamics. Biophys J 59:945–949
Engoren M, Courtney SE, Habib RH (1985) (2009) Effect of weight and age on respiratory complexity in premature neonates. J Appl Physiol 106:766–773
Peng CK, Mietus JE, Liu Y, Lee C, Hausdorff JM, Stanley HE, Goldberger AL, Lipsitz LA (2002) Quantifying fractal dynamics of human respiration: age and gender effects. Ann Biomed Eng 30:683–692
Burioka N, Cornelissen G, Halberg F, Kaplan DT (2001) Relationship between correlation dimension and indices of linear analysis in both respiratory movement and electroencephalogram. Clin Neurophysiol 112:1147–1153
Busha BF (2010) Exercise modulation of cardiorespiratory variability in humans. Respir Physiol Neurobiol 172:72–80
Grassi M, Caldirola D, Vanni G, Guerriero G, Piccinni M, Valchera A, Perna G (2013) Baseline respiratory parameters in panic disorder: a meta-analysis. J Affect Disord 146:158–173
Martinez JM, Kent JM, Coplan JD, Browne ST, Papp LA, Sullivan GM, Kleber M, Perepletchikova F, Fyer AJ, Klein DF, Gorman JM (2001) Respiratory variability in panic disorder. Depress Anxiety 14:232–237
Pine DS, Coplan JD, Papp LA, Klein RG, Martinez JM, Kovalenko P, Tancer N, Moreau D, Dummit ES 3rd, Shaffer D, Klein DF, Gorman JM (1998) Ventilatory physiology of children and adolescents with anxiety disorders. Arch Gen Psychiatry 55:123–129
Yeragani VK, Radhakrishna RK, Tancer M, Uhde T (2002) Nonlinear measures of respiration: respiratory irregularity and increased chaos of respiration in patients with panic disorder. Neuropsychobiology 46:111–120
Bokov P, Fiamma MN, Chevalier-Bidaud B, Chenivesse C, Straus C, Similowski T, Delclaux C (1985) (2016) Increased ventilatory variability and complexity in patients with hyperventilation disorder. J Appl Physiol 120:1165–1172
Bansal T, Haji GS, Rossiter HB, Polkey MI, Hull JH (2018) Exercise ventilatory irregularity can be quantified by approximate entropy to detect breathing pattern disorder. Respir Physiol Neurobiol 255:1–6
Mangin L, Fiamma MN, Straus C, Derenne JP, Zelter M, Clerici C, Similowski T (2008) Source of human ventilatory chaos: lessons from switching controlled mechanical ventilation to inspiratory pressure support in critically ill patients. Respir Physiol Neurobiol 161:189–196
Brack T, Jubran A, Tobin MJ (2002) Dyspnea and decreased variability of breathing in patients with restrictive lung disease. Am J Respir Crit Care Med 165:1260–1264
Loveridge B, West P, Anthonisen NR, Kryger MH (1984) Breathing patterns in patients with chronic obstructive pulmonary disease. Am Rev Respir Dis 130:730–733
Dames KK, Lopes AJ, de Melo PL (2014) Airflow pattern complexity during resting breathing in patients with COPD: effect of airway obstruction. Respir Physiol Neurobiol 192:39–47
Teulier M, Fiamma MN, Straus C, Similowski T (2013) Acute bronchodilation increases ventilatory complexity during resting breathing in stable COPD: toward mathematical biomarkers of ventilatory function? Respir Physiol Neurobiol 185:477–480
Veiga J, Lopes AJ, Jansen JM, Melo PL (1985) (2011) Airflow pattern complexity and airway obstruction in asthma. J Appl Physiol 111:412–419
Chhabra SK, Gupta M, Ramaswamy S, Dash DJ, Bansal V, Deepak KK (2015) Cardiac sympathetic dominance and systemic inflammation in COPD. COPD 12:552–559
Heindl S, Lehnert M, Criee CP, Hasenfuss G, Andreas S (2001) Marked sympathetic activation in patients with chronic respiratory failure. Am J Respir Crit Care Med 164:597–601
Andreas S, Haarmann H, Klarner S, Hasenfuss G, Raupach T (2014) Increased sympathetic nerve activity in COPD is associated with morbidity and mortality. Lung 192:235–241
Hansen J, Sander M (2003) Sympathetic neural overactivity in healthy humans after prolonged exposure to hypobaric hypoxia. J Physiol 546:921–929
Raupach T, Bahr F, Herrmann P, Luethje L, Heusser K, Hasenfuss G, Bernardi L, Andreas S (2008) Slow breathing reduces sympathoexcitation in COPD. Eur Respir J 32:387–392
Khoo MC, Kronauer RE, Strohl KP, Slutsky AS (1982) Factors inducing periodic breathing in humans: a general model. J Appl Physiol Respir Environ Exerc Physiol 53:644–659
Naughton MT (1998) Pathophysiology and treatment of Cheyne-Stokes respiration. Thorax 53:514–518
Brack T, Thuer I, Clarenbach CF, Senn O, Noll G, Russi EW, Bloch KE (2007) Daytime Cheyne-Stokes respiration in ambulatory patients with severe congestive heart failure is associated with increased mortality. Chest 132:1463–1471
Cherniack NS, Longobardo G, Evangelista CJ (2005) Causes of Cheyne-Stokes respiration. Neurocrit Care 3:271–279
Longobardo GS, Cherniack NS, Fishman AP (1966) Cheyne-Stokes breathing produced by a model of the human respiratory system. J Appl Physiol 21:1839–1846
Hanly P, Zuberi N, Gray R (1993) Pathogenesis of Cheyne-Stokes respiration in patients with congestive heart failure. Relationship to arterial PCO2. Chest 104:1079–1084
Naughton M, Benard D, Tam A, Rutherford R, Bradley TD (1993) Role of hyperventilation in the pathogenesis of central sleep apneas in patients with congestive heart failure. Am Rev Respir Dis 148:330–338
Smith LA, Vennelle M, Gardner RS, McDonagh TA, Denvir MA, Douglas NJ, Newby DE (2007) Auto-titrating continuous positive airway pressure therapy in patients with chronic heart failure and obstructive sleep apnoea: a randomized placebo-controlled trial. Eur Heart J 28:1221–1227
Takasaki Y, Orr D, Popkin J, Rutherford R, Liu P, Bradley TD (1989) Effect of nasal continuous positive airway pressure on sleep apnea in congestive heart failure. Am Rev Respir Dis 140:1578–1584
Steens RD, Millar TW, Su X, Biberdorf D, Buckle P, Ahmed M, Kryger MH (1994) Effect of inhaled 3% CO2 on Cheyne-Stokes respiration in congestive heart failure. Sleep 17:61–68
Young T, Palta M, Dempsey J, Skatrud J, Weber S, Badr S (1993) The occurrence of sleep-disordered breathing among middle-aged adults. N Engl J Med 328:1230–1235
Eckert DJ, White DP, Jordan AS, Malhotra A, Wellman A (2013) Defining phenotypic causes of obstructive sleep apnea. Identification of novel therapeutic targets. Am J Respir Crit Care Med 188:996–1004
Younes M (2014) CrossTalk proposal: elevated loop gain is a consequence of obstructive sleep apnoea. J Physiol 592:2899–2901
Wellman A, Jordan AS, Malhotra A, Fogel RB, Katz ES, Schory K, Edwards JK, White DP (2004) Ventilatory control and airway anatomy in obstructive sleep apnea. Am J Respir Crit Care Med 170:1225–1232
Younes M, Ostrowski M, Thompson W, Leslie C, Shewchuk W (2001) Chemical control stability in patients with obstructive sleep apnea. Am J Respir Crit Care Med 163:1181–1190
Bruce EN (1995) Mechanisms and analysis of ventilatory stability. In: Dempsey JA, Pack AI (eds) Regulation of breathing. Marcel Dekker Inc., New York, pp 285–313
Asanoi H, Harada D, Oda Y, Ueno H, Takagawa J, Ishise H, Goso Y, Joho S, Inoue H (2017) Independent prognostic importance of respiratory instability and sympathetic nerve activity in patients with chronic heart failure. J Cardiol 70:476–483
Roberts AM, Bhattacharya J, Schultz HD, Coleridge HM, Coleridge JCG (1986) Stimulation of pulmonary vagal afferent C-fibers by lung edema in dogs. Circ Res 58:512–522
Okamoto S, Ishii M, Hibi S, Akishita M, Yamaguchi Y (2021) Breathing irregularities before sleep onset on polysomnography in patients with heart diseases. Sleep Breath. 26(2):605–612
Seely AJ, Bravi A, Herry C, Green G, Longtin A, Ramsay T, Fergusson D, McIntyre L, Kubelik D, Maziak DE, Ferguson N, Brown SM, Mehta S, Martin C, Rubenfeld G, Jacono FJ, Clifford G, Fazekas A, Marshall J, Canadian Critical Care Trials G (2014) Do heart and respiratory rate variability improve prediction of extubation outcomes in critically ill patients? Crit Care 18:R65
Wysocki M, Cracco C, Teixeira A, Mercat A, Diehl JL, Lefort Y, Derenne JP, Similowski T (2006) Reduced breathing variability as a predictor of unsuccessful patient separation from mechanical ventilation. Crit Care Med 34:2076–2083
El-Khatib M, Jamaleddine G, Soubra R, Muallem M (2001) Pattern of spontaneous breathing: potential marker for weaning outcome. Spontaneous breathing pattern and weaning from mechanical ventilation. Intensive Care Med 27:52–58
Engoren M (1998) Approximate entropy of respiratory rate and tidal volume during weaning from mechanical ventilation. Crit Care Med 26:1817–1823
Miyata M, Burioka N, Sako T, Suyama H, Fukuoka Y, Tomita K, Higami S, Shimizu E (2004) A short daytime test using correlation dimension for respiratory movement in OSAHS. Eur Respir J 23:885–890
Raoufy MR, Ghafari T, Darooei R, Nazari M, Mahdaviani SA, Eslaminejad AR, Almasnia M, Gharibzadeh S, Mani AR, Hajizadeh S (2016) Classification of asthma based on nonlinear analysis of breathing pattern. PLoS ONE 11:e0147976
Finer NN, Abroms IF, Taeusch HW Jr (1976) Ventilation and sleep states in newborn infants. J Pediatr 89:100–108
Haddad GG, Jeng HJ, Lai TL, Mellins RB (1987) Determination of sleep state in infants using respiratory variability. Pediatr Res 21:556–562
Harper RM, Schechtman VL, Kluge KA (1987) Machine classification of infant sleep state using cardiorespiratory measures. Electroencephalogr Clin Neurophysiol 67:379–387
Terrill PI, Wilson SJ, Suresh S, Cooper DM, Dakin C (2010) Attractor structure discriminates sleep states: recurrence plot analysis applied to infant breathing patterns. IEEE Trans Biomed Eng 57:1108–1116
Marwan N, Romano MC, Thiel M, Kurths J (2007) Recurrence plots for the analysis of complex systems. Phys Rep 438:237–329
Gaiduk M, Penzel T, Ortega JA, Seepold R (2018) Automatic sleep stages classification using respiratory, heart rate and movement signals. Physiol Meas 39:124008
Nikkonen S, Korkalainen H, Leino A, Myllymaa S, Duce B, Leppanen T, Toyras J (2021) Automatic respiratory event scoring in obstructive sleep apnea using a long short-term memory neural network. IEEE J Biomed Health Inform 25:2917–2927
Yeo M, Byun H, Lee J, Byun J, Rhee HY, Shin W, Yoon H (2022) Respiratory event detection during sleep using electrocardiogram and respiratory related signals: using polysomnogram and patch-type wearable device data. IEEE J Biomed Health Inform 26:550–560
Ekman P, Levenson RW, Friesen WV (1983) Autonomic nervous system activity distinguishes among emotions. Science 221:1208–1210
Ayata D, Yaslan Y, Kamasak ME (2020) Emotion recognition from multimodal physiological signals for emotion aware healthcare systems. J Med Biol Eng 40:149–157
Egger M, Ley M, Hanke S (2018) Emotion recognition from physiological signal analysis: a review. Electr Notes Theor Comput Sci 343:35–55
Kyamakya K, Al-Machot F, Mosa AH, Bouchachia H, Chedjou JC, Bagula A (2021) Emotion and stress recognition related sensors and machine learning technologies. Sensors-Basel 21:2273
Raheel A, Majid M, Alnowami M, Anwar SM (2020) Physiological sensors based emotion recognition while experiencing tactile enhanced multimedia. Sensors 20:4037
Cheng XF, Yue W, Shicheng D, Pengjun Z, Qifa L (2019) Heart sound signals can be used for emotion recognition. Sci Rep 9:6486
Dzedzickis A, Kaklauskas A, Bucinskas V (2020) Human emotion recognition: review of sensors and methods. Sensors-Basel 20:592
Jerath R, Beveridge C (2020) Respiratory rhythm, autonomic modulation, and the spectrum of emotions: the future of emotion recognition and modulation. Front Psychol 11:1980
Suzuki K, Laohakangvalvit T, Matsubara R, Sugaya M (2021) Constructing an emotion estimation model based on EEG/HRV indexes using feature extraction and feature selection algorithms. Sensors-Basel 21:2910
Zhang Q, Chen XX, Zhan QY, Yang T, Xia SH (2017) Respiration-based emotion recognition with deep learning. Comput Ind 92–93:84–90
Posner J, Russell JA, Peterson BS (2005) The circumplex model of affect: an integrative approach to affective neuroscience, cognitive development, and psychopathology. Dev Psychopathol 17:715–734
Russell JA (1980) A circumplex model of affect. J Pers Soc Psychol 39:1161–1178
Koelstra S, Muhl C, Soleymani M, Lee JS, Yazdani A, Ebrahimi T, Pun T, Nijholt A, Patras I (2012) DEAP: a database for emotion analysis using physiological signals. IEEE T Affect Comput 3:18–31
Sharma K, Castellini C, van den Broek EL, Albu-Schaeffer A, Schwenker F (2019) A dataset of continuous affect annotations and physiological signals for emotion analysis. Sci Data 6:196
Acknowledgements
Not applicable.
Funding
This work was supported by JSPS KAKENHI Grant Number 20K06917.
Author information
Authors and Affiliations
Contributions
YO wrote the whole manuscript. The author read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The author declares no conflicts of interest in association with this study.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Oku, Y. Temporal variations in the pattern of breathing: techniques, sources, and applications to translational sciences. J Physiol Sci 72, 22 (2022). https://doi.org/10.1186/s12576-022-00847-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12576-022-00847-z