Journal of Applied Physiology Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 97: 1830-1839, 2004. First published June 18, 2004; doi:10.1152/japplphysiol.00298.2004
8750-7587/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
97/5/1830    most recent
00298.2004v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (8)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Baldwin, D. N.
Right arrow Articles by Frey, U.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Baldwin, D. N.
Right arrow Articles by Frey, U.

Effect of sighs on breathing memory and dynamics in healthy infants

D. N. Baldwin,1,3 B. Suki,2 J. J. Pillow,3 H. L. Roiha,1 S. Minocchieri,1 and U. Frey1

1Pediatric Respiratory Medicine, University Children's Hospital, Bern 3010, Switzerland; 2Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215; and 3Clinical Sciences Division, Institute for Child Health and Centre for Child Health Research, University of Western Australia, Perth 6008, Australia

Submitted 18 March 2004 ; accepted in final form 16 June 2004


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Deep inspirations (sighs) play a significant role in altering lung mechanical and airway wall function; however, their role in respiratory control remains unclear. We examined whether sighs act via a resetting mechanism to improve control of the respiratory regulatory system. Effects of sighs on system variability, short- and long-range memory, and stability were assessed in 25 healthy full-term infants at 1 mo of age [mean 36 (range 28–57) days] during quiet sleep. Variability was examined using moving-window coefficient of variation, short-range memory using autocorrelation function, and long-range memory using detrended fluctuation analysis. Stability was examined by studying the behavior of the attractor with use of phase-space plots. Variability of tidal volume (VT) and minute ventilation (E) increased during the initial 15 breaths after a sigh. Short-range memory of VT decreased during the 50 breaths preceding a sigh, becoming uncorrelated (random) during the 10-breath presigh window. Short-range memory increased after a sigh for the entire 50 breaths compared with the randomized data set and for 20 breaths compared with the presigh window. Similar, but shorter duration, changes were noted in E. No change in long-range memory was seen after a sigh. Coefficient of variation and range of points located within a defined attractor segment increased after a sigh. Thus control of breathing in healthy infants shows long-range stability and improvement in short-range memory and variability after a sigh. These results add new evidence that the role of sighs is not purely mechanical.

nonlinear dynamics; feedback; control; systems theory


INFANTS DEMONSTRATE DYNAMIC regulation of lung function over time. Lung volume is actively elevated above that determined by tissue elastic properties and a compliant chest wall (26). Recent studies suggested that dynamic regulation of airway tone is particularly important in infants (16). It is likely that this dynamic regulatory process represents a tightly controlled negative-feedback system, whereby information from peripheral and central chemoreceptors, as well as pulmonary stretch receptors, is provided to respiratory regulatory neurons in the brain stem, with continuously variable output to the respiratory system musculature (43, 44).

Biological feedback systems such as those utilized in, e.g., respiratory control, are reflective of engineering control systems and may be modeled using similar mathematical concepts (9, 35). From an engineering point of view, a control system requires an appropriate balance between system stability and sensitivity, permitting responsiveness to fluctuations while maintaining function within tightly regulated limits. Excess sensitivity of such a system may lead to dramatic "overshoot"-type responses, instability, and positive feedback (24). Reduction in sensitivity and excessive stability may result in diminished capacity for response to environmental change and poor feedback control. On the other hand, a certain amount of stability is required in order for the biological system to return to a steady state after a perturbation. Variability observed in the output parameters of a system is often the result of these two competing mechanisms. Long- and short-range memory within such a dynamical system are responsible for maintaining this delicate equilibrium by incorporating information from previous events into the current system state. An ideal error-tolerant system returns quickly to a steady state, whereby long-range memory is not altered. Long-range memory may therefore be considered a feature of stability.

From engineering systems, it is known that many control loop feedback systems incorporate a resetting mechanism. There is evidence that spontaneous deep inspirations (sighs) play an important role in resetting the mechanical properties of the lung tissue (10, 12) and airway walls (7); however, there is little information in the literature as to whether sighs also act as a "resetting" mechanism for the control loop system. In adults, sighs during sleep are associated with increased inspiratory and expiratory flows and increased inspiratory time (TI) as a fraction of respiratory cycle time (TI/TT), indicating that they may play an important role in the control process (31). Sighs are a more frequently observed phenomenon in infant than in adult breathing (19, 31) and, given the increased contribution of dynamic processes to respiratory control in this age group, may make a major contribution to homeokinesis. Fleming et al. (13) used a linear approach to investigate the developmental changes in oscillatory responses of breath-to-breath minute ventilation (E) after a sigh in ≤7-mo-old infants. Newborn infants had a stable system with delayed responses to the sigh stimulus that evolved with age to exhibit a stable, but more responsive system suggestive of improvement in feedback control. Animal studies demonstrate no difference in respiratory pattern after a sigh between rapid eye movement (REM) and non-REM (NREM) sleep; however, a smaller tidal volume (VT) and longer expiratory time (TE) are evident, resolving on most occasions by the third postsigh breath (20). The factors responsible for inducing a sigh and the effects of this event on the stability of the system remain unclear.

We hypothesized that sighs may play a crucial role not only for lung mechanics but also as a resetting mechanism for the negative-feedback control system. To examine whether full-term healthy infants demonstrate alterations in regulation of breathing after sighs during quiet natural sleep, we analyzed the effect of sighs on tidal breathing indexes in terms of respiratory variability, short- and long-range memory, and stability.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Study Design

This cross-sectional observational study examined regulation of breathing in relation to sighs in healthy infants during quiet natural sleep. Breath-by-breath VT, E, and minimum exhaled O2 and end-tidal CO2, as well as the classical parameters of respiratory drive [TI/TT, VT/TI, time to peak tidal expiratory flow (tPTEF)/TE, TI, TE, and maximal inspiratory flow (Imax)], were examined for 50 breaths before and after a spontaneous sigh. Data were analyzed in terms of variability, short- and long-range memory, and stability of VT and E. Variability was examined using moving-window coefficient of variation (CV). Short-range memory was assessed for data sets of 50 breaths using an autocorrelation function (c) (23). Detrended fluctuation analysis (28) examined long-range memory over 100 breaths. Stability of the system was assessed using phase-space plots (23).

Subjects

Data were obtained from a subgroup of 25 infants from a total cohort of 80 healthy full-term infants recruited antenatally to a separate large epidemiological study. Data from each infant were visually examined for the presence of a sigh. A sigh was defined as a tidal breath at least double the mean VT of the preceding 10 breaths (10, 40). A sigh breath was identified in 40 infants; however, 15 were excluded because of the absence of NREM sleep state, insufficient data adjacent to the sigh (<50 breaths), significant data artifact or when a change in sleep state, movement, or arousal occurred during the study. This left data from 25 infants that were suitable for analysis. Demographic data are provided in Table 1. All infants were born after unremarkable pregnancies, and delivery was normal. None of the infants had respiratory infections between birth and date of the study. The study was approved by the Ethics Committee of the University Hospital of Berne. Written informed parental consent was obtained before study, and parents were most often present at the time of measurement.


View this table:
[in this window]
[in a new window]
 
Table 1. Demographic data

 
Measurements

The technique for measurement of tidal breathing parameters in spontaneously sleeping infants has been published in detail previously (8). Briefly, infants were studied in the supine position in a cot or in the mother's arms during quiet natural sleep. Sleep state was defined clinically using the criteria of Prechtl (33), with closed eyes, absence of gross body or limb movements, and regular respiratory pattern. Heart rate and O2 saturation were monitored continuously for the duration of study (Biox 3700, Datex-Ohmeda, Helsinki, Finland). Measurements were made over periods of up to 10 min using methodology and equipment compliant with recently published international standards for measurement of lung function in infancy (5, 17, 18).

Tidal breathing was measured by application of a compliant silicone face mask to the infant's mouth and nose (size 1, Homedica, Cham, Switzerland) after the face mask was warmed to skin temperature (17). The mask was connected in series with a CO2 analyzer, combined nitric oxide and O2 sensor (Exhalyser, EcoMedics, Duernten, Switzerland), and prototype ultrasonic flowmeter (Spiroson model M30.8001, EcoMedics) containing a bacterial filter (Spirette) surrounding a dead space (VD) reduction apparatus (size 1). The measurement setup was connected to a bias flow of 14 l/min.

Combined instrumental VD without the face mask was 3 ml. Total VD of the silicone mask was 15 ml. The effective total measurement VD was 10.5 ml, with 50% of the mask space assumed to be occupied by the infant's face. This gives a mean VD of ~2.2 ml/kg body wt, which is at the upper limit of current recommendations (17) for short-range measurements. End-tidal CO2 was monitored for the entire test procedure and did not increase during the measurements. Signals were 12-bit analog-to-digital converted and sampled at 200 Hz using a custom-designed data acquisition and analysis package (WBreath version 3.7.6.0 [EC] , Firmware v3.06, EcoMedics). The sampling time delay was corrected for each signal.

Data Analysis

Data were analyzed using custom software written in Matlab (Matlab 6, version 12, Mathworks). VT was integrated from the flow signal after correction for any technical signal drift. E was calculated by multiplying VT by the respiratory frequency (f) on a breath-by-breath basis. A representative example of VT and associated sigh used for analysis is shown as a function of time in Fig. 1A. Breath-to-breath VT was plotted as a function of breath number for a 10-min time series (Fig. 1B). The sigh breath itself was excluded from the analysis.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 1. Sigh data from a representative infant. A: raw tidal flow data. B: breath-by-breath tidal volume (VT).

 
Classical tidal breathing parameters associated with respiratory drive.   We analyzed the mean, along with 95% confidence interval (CI), of TI/TT, VT/TI, tPTEF/TE, TI, TE, and Imax (31, 41) for 50 breaths before and after the sigh and compared them using paired t-tests. Moving-window CV analysis was used to determine whether these parameters were predictive for the sigh or changed after the sigh (see below).

Change in minimum exhaled O2 and end-tidal CO2 in response to a sigh.   To examine the possible influence of alterations in O2 and CO2 on sigh and postsigh control of breathing, breath-to-breath changes in minimum exhaled O2 and end-tidal CO2 were investigated in 18 of 25 and 25 of 25 infants, respectively. In seven cases, O2 was not measured, inasmuch as an analyzer was not present within our equipmental setup. Mean, along with 95% CI, minimum exhaled O2 and end-tidal CO2 values were calculated for the 10-breath period immediately before the sigh and compared with values for each breath after the sigh for a period of 20 breaths.

Variability of VT and E.   Variability was examined using an established moving-window algorithm whereby CV was obtained by dividing the SD by the mean for repeated windows of 10 breaths (23, 29). A 50% window overlap was used for each subsequent measurement until the end of the 50-breath series. This technique provided 8 windows of 10 breath periods for each of the 50-breath pre- and postsigh data series. Comparisons were made between changes of CV in the period leading to the sigh, alteration of CV in response to the sigh, and the period required for return to the baseline level of CV after the sigh. Results for all 25 infants were pooled and mean (95% CI) for each window was displayed as a function of window number before and after the sigh.

Short-range memory.   Within each analysis window containing 10 breaths as described above, the breath-to-breath c (Matlab 6, version 12, Mathworks) was also calculated with a lag of 1 breath [c(1)], which represents the correlation coefficient between 2 consecutive breaths; c(1) = 1 indicates a 100% correlation between adjacent breaths, whereas c(1) = 0 indicates that they are totally uncorrelated or random. Statistical analysis was performed for the entire group (n = 25).

To test whether correlations represented the presence of true memory, data within an individual window were shuffled 10 times (i.e., the order of data points within the window was randomized), with the autocorrelation calculation repeated on each occasion. The average c(1) for the 10 shuffles was then obtained for each window. By definition, a randomized data set such as this is uncorrelated, and c(1) should approximate 0. The existence of a significant difference in c(1) for the whole group of infants before and after shuffling indicated the presence of true short-range memory.

Long-range memory.   In a subgroup of 12 infants, 100 tidal breaths were available on either side of the sigh, such that long-range correlations could be examined using a modified form of detrended fluctuation analysis (28). Calculation of the correlation exponent ({alpha}) with this algorithm uses a modified root-mean-square analysis, a detailed description of which has been reported previously (8, 28, 30). Briefly, the breath-by-breath VT and E time series before and after the sigh were analyzed as follows. The breath-to-breath series (xi = x1, x2,..., xN) was first integrated [y(k)]

The series was then divided into equally sized nonoverlapping windows of length n, and a linear regression line was fit through the data in each window, giving the local trend of the data y(k). The fluctuation function F(n) of the time series corresponding to each window length (n) was calculated after detrending of data within each window by subtraction of the local trend [yn(k)]

If the functional form of F(n) follows a power law, F(n) ~ n{alpha}, the process is said to exhibit long-range correlation. The exponent {alpha} can be obtained as the slope of a straight-line fit to F(n) vs. n on a log-log plot (Fig. 2). A random process with uncorrelated data is represented by {alpha} = 0.5. A positively correlated signal occurs when {alpha} > 0.5.



View larger version (7K):
[in this window]
[in a new window]
 
Fig. 2. Detrended fluctuation analysis from a representative infant (presigh VT). Fluctuation function [F(n)] was calculated at increasing window sizes (n) after local detrending was performed within each window. Correlation exponent ({alpha}) is the gradient on a log-log plot and represents long-range memory of the system. Data from 100 breaths before and after a sigh were examined, and {alpha} values were compared for the entire group (n = 12) before and after a sigh using paired t-test.

 
Stability.   The dynamics of a system can be characterized by its attraction toward a steady state within a setting of continuously occurring intrinsic and extrinsic fluctuations that tend to perturb it away from this stable state or "basin" of attraction. We consider a sigh to be a local nonstationarity within the respiratory control system and observe the return of the system to a steady state in response to this transient stimulus. We used phase-space plots of raw tidal flow data to assess changes in system stability after a sigh (23, 37). To undertake such an analysis, the scalar data series must first be converted to a geometric reconstruction of the system dynamics. A graphical representation of this form is known as an embedding within phase space. The plotted variables constitute the trajectory of the system behavior as it evolves in time. Consequently, the points within the phase space represent all possible states of the system for the period of observation and comprise the attractor (4).

The attractor was reconstructed using the technique of time-delay embedding (39). For a given time series si = (s1,s2, ..., sN) of N measurements sampled at a fixed rate, a data point si within the scalar time series is plotted against a point si + {xi} that is "delayed" in time within the series by a given value {xi}, called the embedding delay. The appropriate embedding of the vector time series (Si) is given by the relation {xi} = m{tau}, where m represents the number of dimensions in which the data are to be embedded and {tau} is the time lag. The dimension is therefore represented by the number of delayed points used in the reconstruction of the attractor. Each time a point within the scalar series is related to another point located a given delayed time distant, a dimension is added to the reconstruction. Hence, a reconstruction in three dimensions relates a given reference point with two other points located at specified (and different) delays from the point of interest

where i = 1, 2, ..., N – (m – 1){xi}. In each infant, 50 breaths of raw tidal flow data (sampled every 5 ms) were examined visually within a two- or three-dimensional embedded phase space. Each pre- and postsigh plot therefore comprised 15,000–18,000 data points. Different embedding time delays were systematically examined to obtain the best unfolding of the attractor dynamics in phase space. In most instances, the embedding delay was similar between infants, and in each infant the same embedding delay was used to examine pre- and postsigh data. As would be expected from a periodic or quasi-periodic signal, the embedding delay was related to the breathing frequency, in most cases on the order of 400 ms. This delay is similar in magnitude to that reported previously for respiratory data embedding (37) and is representative of the dominant period relation observable in low-dimensional attractor systems in which a self-excited or forced periodic behavior is present. Within such systems, a useful starting point for identification of the time delay is around one-quarter of the dominant period (1). We have selected an embedding delay of ~0.25 of normal respiratory cycle period to provide an acceptable unfolding of the attractor dynamics.

To quantify the dynamics before and after a sigh, we examined the variability between cycles of the attractor within phase space plotted in two dimensions. A small region of one limb of the attractor cycles was selected to obtain at least one point from each embedded cycle. Data from the selected region were normally distributed in each infant, and variability of the cycles around the mean was assessed in terms of CV and range.

Statistical Analysis

Comparisons of indexes of variability, autocorrelation, and detrended fluctuation analysis, as well as pre- and postsigh attractor dynamics, were made for the grouped data using paired t-test. Values are means and 95% CI unless otherwise stated. Results were considered statistically significant when P < 0.05.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Classical Tidal Breathing Parameters Associated With Respiratory Drive

When examined over the 50-breath period, mean VT, TI, and TI/TT decreased after a sigh (Table 2). Other parameters of respiratory drive (TE, VT/TI, tPTEF/TE, and Imax) were unchanged. More detailed analysis using a moving-window CV algorithm showed that most of the observed alteration in TI and T/TT occurred during the first one to two windows (10–15 breaths) after a sigh.


View this table:
[in this window]
[in a new window]
 
Table 2. Classical indexes of respiratory drive

 
Change in Minimum Exhaled O2 and End-Tidal CO2 in Response to a Sigh

The minimum exhaled O2 and end-tidal CO2 concentration data demonstrated a high level of stability during the 50-breath period leading to a sigh. Consequently, the mean, along with 95% CI, value for the 10 breaths immediately before a sigh were used to compare data after a sigh (Fig. 3). A significant increase in minimum exhaled O2 was evident for the first two breaths after a sigh compared with the presigh values when examined for the entire group. Changes in end-tidal CO2 lasted slightly longer, with a significant reduction to the fifth postsigh breath.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 3. Response of minimum exhaled O2 and end-tidal CO2 to a sigh. Data are group means with 95% confidence intervals (CI). Postsigh breaths were compared with the mean of 10 breaths immediately before the sigh. Minimum exhaled O2 increased significantly from before the sigh for the first 2 postsigh breaths: *P < 0.05. The 4th postsigh breath achieved marginal significance: §P = 0.05. End-tidal CO2 returned to presigh levels within 6 breaths after the sigh.

 
Variability of VT and E

CV of VT increased in 23 of 25 infants after a sigh when examined over the entire 50-breath period. For the moving-window analysis, cross-sectional data were normally distributed, allowing the use of mean and 95% CI for presentation of data. Individual VT variability data for each analysis window were grouped for all 25 infants and are presented in Fig. 4A. Baseline VT variability was not significantly different between measurement windows during the presigh period (P = 0.22), implying that no indication of the approaching sigh was evident in the presigh variability. Compared with the 10-breath window immediately before the sigh, there was a significant increase in variability during the first two analysis windows after the sigh, representing a period of 15 breaths. Similar changes in variability were observed in E (Fig. 4B).



View larger version (10K):
[in this window]
[in a new window]
 
Fig. 4. Change in variability [coefficient of variation (CV)] of VT (A) and minute ventilation (E, B) in response to a sigh. Data are group (n = 25) means with 95% CI for each analysis window before and after a sigh. Each analysis window examined 10 breaths with a 50% overlap between windows. A: compared with presigh window, variability of VT increased for the first 2 postsigh windows (15 breaths). B: variability of E also increased for a period of 2 analysis windows (15 breaths) after the sigh before returning to baseline levels. *P < 0.05. **P < 0.001.

 
Short-Range Memory

Autocorrelation data for windows located at an increasing distance after the sigh and for the window immediately before the sigh are compared in Fig. 5A. Grouped mean CV of c(1) for VT was significantly increased during the first three postsigh windows, representing a period of 20 breaths (Fig. 5A). In contrast, c(1) of E was significantly higher for the first postsigh window only; the second window had a borderline significance of P = 0.05 (Fig. 5B).



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 5. Change in autocorrelation [c(1)] of VT (A) and E (B) with a sigh. Data are group (n = 25) means with 95% CI for each analysis window before and after a sigh. Each analysis window examined 10 breaths with a 50% overlap between windows. A: c(1) of VT increased during the first 3 postsigh windows compared with the presigh window. This represents a period of 20 breaths over which breath-to-breath short-range memory was increased after the sigh. After this period, memory properties returned to values not significantly different from before the sigh. B: c(1) of E also increased after a sigh; however, changes were shorter lasting, being significantly different for the 1st postsigh window only (10 breaths). The 2nd postsigh window achieved marginal significance (P = 0.05) compared with the presigh window. *P < 0.05. **P < 0.001.

 
Linear regression analysis through the group means for each analysis window before the sigh showed a significant decrease of correlation of VT (r2 = 0.59, P = 0.02) and E (r2 = 0.53, P = 0.04) compared with a randomized series during the period leading to the sigh (Fig. 6). The VT and E series demonstrated a weak breath-to-breath correlation that was significantly different from a randomized series but became less significant as the sigh approached. In the case of the VT series, the original and randomized data were not significantly different for the presigh window (P = 0.08), indicating a loss of deterministic correlation. After the sigh, VT and E showed a significant increase in correlation compared with their respective randomized series, which persisted for the duration of eight postsigh analysis windows (i.e., the full 50 breaths).



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 6. Linear regressions for group mean c(1) values for VT and E at each presigh analysis window. Linear regressions through the grouped mean c(1) for VT and E demonstrated a decrease in correlation toward the impending sigh. Regression for VT (solid line): c(1) = –0.0123x + 0.155, r2 = 0.59. Regression for E (dashed line): c(1) = –0.0132x + 0.192, r2 = 0.53. Data from each infant were shuffled within each analysis window to remove any correlation properties associated with short-range memory in the data set. Compared with the randomized, uncorrelated data series, regression slopes were significantly different for VT (P = 0.02) and E (P = 0.04).

 
Long-Range Memory

The exponent {alpha}, characterizing the long-range correlation of VT and E, was unchanged in response to the sigh (Table 3). The observed trend (P = 0.06) toward change in long-range correlation (over 100 breaths) of VT is consistent with a contribution from short-range correlations observed over 50 breaths.


View this table:
[in this window]
[in a new window]
 
Table 3. Long-range regulation of VT and I before and after a sigh

 
Stability

An example of a typical phase-space embedding in two and three dimensions from a single infant is shown in Fig. 7. The significant alteration in attractor dynamics after a sigh is clearly evident. Similar visual changes were observed in all infants. CV of data points within the analyzed attractor region increased significantly after the sigh in 18 of 25 infants [0.18 (95% CI = 0.09–0.27) and 0.37 (95% CI = 0.22–0.52) before and after the sigh, respectively, P < 0.001]. The range increased in 22 of 25 infants [0.03 (95% CI = 0.02–0.04) and 0.05 (95% CI = 0.04–0.06) before and after the sigh, respectively, P < 0.02]. Importantly, on each occasion, the respiratory behavior returned to the same position within the phase space after the sigh, indicating overall stability of the system.



View larger version (42K):
[in this window]
[in a new window]
 
Fig. 7. Phase-space embedded plots from a representative infant: presigh (A) and postsigh (B) embeddings in 3 dimensions. Each embedding comprises ~15,000 data points from the raw tidal flow signal. A given point within the series is plotted against points with delays of 100 points (500 ms) and 120 points (600 ms). The same delays were used before and after the sigh. Although the location of the attractor within the phase space is unchanged after the sigh, dynamics of the system are visually markedly different. Phase-space embeddings in 2 dimensions were used to quantify changes in dynamics before (C) and after (D) the sigh. In this case, the appropriate delay was 80 points (400 ms). Variability of the postsigh attractor was significantly higher than before the sigh when examined in terms of CV (P < 0.001) and range (P < 0.02).

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Sighs play an important role in altering the mechanical properties of the lungs and airway walls. In bronchoconstricted adults, sighs lead to an increase in specific airway conductance and bronchodilatation secondary to an increase in transpulmonary pressure that is not related to change in lung volume (11). In preterm infants, mechanical effects of sighs include maintenance of lung volume (32), increased dynamic lung compliance, and reduced lung resistance (10), confirming their role in altering lung mechanics.

The exact origin of spontaneous sighs remains controversial. They may result from activation of lung and chest wall receptors in response to reduced pulmonary compliance associated with alveolar atelectasis (6). Animal studies have demonstrated increasing sigh frequency in response to hypoxia (3) and hypercapnia (34), whereas interruption of vagal pathways leads to loss of these feedback responses (3). Vagal neural feedback to the brain stem respiratory center plays a significant role within the reflex pathway responsible for mediating sigh activity (25). For example, the increase in neurally mediated end-expiratory volume and respiratory frequency after sighs in cats is not observed after passive inflation maneuvers (38). The shorter time required for end-expiratory volume than for lung compliance to return to baseline in these cats is further evidence of dissociation between neural and mechanical effects of sighs. The increased sigh frequency reported in infants with an immature breathing control system (e.g., preterm infants) (2, 19) highlights the likely importance of sighs within the neuroregulatory feedback loop.

Recently, regulatory properties of the respiratory negative-feedback loop have been analyzed using mathematical techniques derived from engineering systems control theory (24, 35). The behavior of the control loop after a sigh may be reflective of the maturity and functional integrity of neurorespiratory feedback control. Using these techniques, we investigated whether sighs are important for the regulation and resetting of the neurorespiratory controller. Our data showed that a highly regulated and stable negative-feedback loop system exists in healthy infants and that sighs represent a mechanism for improving the memory associated with neurorespiratory control of breathing within this system.

Classical Tidal Breathing Parameters Associated With Respiratory Drive

This study suggests that established indexes of respiratory drive and control, such as TI/TT, VT/TI, tPTEF/TE, TI, TE, and Imax, are not predictive of an impending sigh. Most of these parameters remained unchanged during the 50-breath postsigh period, although a small, but statistically significant, decrease in TI and TI/TT occurred during the two window periods immediately after the sigh compared with the presigh window. Ueda et al. (41) showed that tPTEF was one factor that identified reduced respiratory drive in infants of smoking mothers, suggesting its potential usefulness for the examination of control. The reason for the insensitivity of most of the classical drive parameters to changes in respiratory control induced by a sigh in this study is unclear. Although TI/TT and TI are higher in adults during presigh breaths than immediately after a sigh (31), physiological differences make comparisons between infants and adults particularly difficult. Techniques permitting analysis of the dynamic characteristics of respiratory drive, such as those derived from system control theory, may reveal more information about the continually changing behavior of this control system than is provided by standard linear methods.

Change in Minimum Exhaled O2 and End-Tidal CO2 in Response to a Sigh

The correlation between arterial O2 and CO2 values and exhaled gas concentrations is known to be limited (42); however, because of the noninvasive nature of this study, we were unable to analyze blood gases. We did, however, monitor O2 and CO2 in exhaled gas continuously. Breath-to-breath minimum exhaled O2 and end-tidal CO2 revealed highly stable values over 50 breaths before and after a sigh, with exhaled concentrations of both gases returning to presigh values within six breaths (Fig. 3). These changes were of a shorter duration than those observed in variability or short-range memory of VT and E after a sigh. Changes in flow and volume during the first few breaths after a sigh may impact on the concentrations of exhaled gases. When examined over a longer period, the mean minimum O2 and end-tidal CO2 were not different from presigh values.

Variability of VT and E

A typical feature of an adaptive control system operating in a noisy environment is variability. A system locked into a steady state is likely to be insensitive to small perturbations within its external environment. Reduction of physiological variability has been associated with morbidity and mortality, e.g., in the reduced heart rate variability observed with age (22, 30) and heart disease (29) and in extreme cases preceding life-threatening arrhythmias (27, 36). In the absence of changes in external noise, changes in the variability of breathing induced by a sigh may be considered an indicator of changes in neurorespiratory feedback control. Variability in VT or E does not predict the occurrence of a sigh. In contrast, the sigh itself introduces a significant increase in variability of VT and E for a period of 15 breaths (Fig. 4), after which the system returns rapidly to a baseline level of fluctuation.

Memory in Respiratory Control

The existence of memory within a feedback system allows the recent history (or state) to play a role in the response to a transient perturbation, thus assisting with return to steady-state behavior. Many biological systems exhibit long-range memory characterized by a power law in which decay of correlations slowly fades with time. These systems may be said to demonstrate scale-invariant or fractal characteristics. Although this has been shown for heart rate (29) and breathing (30) in adults, we recently demonstrated such long-range memory in infant breathing (8). Detrended fluctuation analysis (28) performed on E during this study identified no change of long-range correlation after a sigh in healthy infants. A trend toward an increase in {alpha} of VT after sigh might be related to the observed changes in short-range memory described below. The presence of a power law memory is likely to be beneficial for the long-range stability of the control system. However, such analysis does not reveal the rapid dynamics after a perturbation of the system.

The analysis techniques used in this study allowed us to follow rapid changes in system dynamics. In the context of an engineering control system, the presence of short-range memory may act like a low-pass filter for a sudden perturbation, allowing a smoother response to the stimulus. The loss of short-range memory [decrease in c(1)] as the sigh approached and apparent random behavior of VT immediately before the sigh suggest that the complete loss of correlations in VT may predict the next sigh. After the sigh, c(1) was significantly increased compared with the randomized series for the 50 breaths after the sigh but steadily decreased probably toward the onset of the next sigh. These data are compatible with a recovery of short-range memory of the respiratory control system after the sigh, which is consistent with the hypothesis that sighs may act as a resetting mechanism for the control loop. Our findings are supported by the data of Franco et al. (15), who recently demonstrated differences in the response to a sigh between healthy infants and those who later succumbed to sudden infant death syndrome. Together, these two studies show that sighs may play an important role in resetting autonomic tone in NREM sleep.

Stability of the Control System

To characterize the clear increase in fluctuations of the attractor after a sigh, we analyzed a segment of the attractor path within a predetermined region of two-dimensional phase space (Fig. 7). Compared with the locked-in and unresponsive behavior of the presigh attractor (Fig. 7A), there was a clear increase in fluctuations of the attractor after a sigh (Fig. 7B). The comparatively low variability indexes observed before a sigh occurred alongside a gradual reduction in short-range memory of the system, and we speculate that this reduced memory and apparent reduction in system responsiveness are intimately related. The significant increase in path fluctuation suggests that the attractor was less "locked-in" and more responsive to perturbations (Fig. 7B). This marked increase in cycle fluctuation within phase space is confirmed by the significant increases in the range and variability of data points obtained from the attractor segments. Importantly, the breathing control system returned to the same basin of attraction during the period after the sigh in all infants, as indicated by the unchanged position of the attractor within the phase-space plot. Thus, although the system is more variable and, perhaps, better regulated, it remains stable. It is worth noting that stability and variability represent different phenomena in system behavior and often coexist. Furthermore, an increase in variability (e.g., CV) does not necessarily imply an alteration of correlation properties. Consequently, a system may be intrinsically stable while demonstrating significant variability in response to, e.g., extrinsic fluctuations (24).

Technical Issues and Limitations of the Study

Because the aim of this study was to assess the alterations in neuroregulatory control after a sigh, we examined the output parameters of this complex system in terms of tidal breathing indexes. The techniques employed during this study allowed us to investigate the global properties of this feedback system. A limitation of these methods was the inability to separate and define the contributions of the various input parameters to this regulatory system, e.g., those from chemoreceptors and pulmonary stretch receptors.

From a physiological point of view, the conditions under which measurements are made are important. Measurement conditions were standardized as far as possible, while the noninvasive nature of this study was maintained. Nevertheless, even the pressure and temperature of the face mask can influence control of breathing (14). The mask was warmed to skin temperature before application, and infants were given time to adapt to the new face-mask conditions before measurements were begun. Although infants were studied in a cot or in the mother's arms, in all cases, position of the body and head were standardized. A major strength of this study lies in the homogeneous population of healthy infants, all of whom were measured under conditions of quiet unsedated sleep. Such studies must be performed carefully, because sigh behavior of infants during different sleep stages may vary (19). In this study, we undertook careful prospective evaluation of behaviorally defined sleep patterns (33) to avoid the possibility of observing alterations in sleep state after a sigh.

Sighing occurs about every 50–100 breaths in healthy infants. Thus the information on long-range correlations is limited to a maximum of ~100 breaths. Detrended fluctuation analysis was designed for much longer data sets, and a limited number of data points can introduce bias. Using a correlated noise signal with a known {alpha} of 1.0, we have tested the effect of reducing the number of data points from 1,000 to 100. Using a data set size 100, we found that the difference in {alpha} was ≤15% compared with the larger series. To partly overcome this limitation of data set size, we always compared the same number of data points before and after a sigh. Phase-space plot analysis is dependent on the correct choice of raw data, the embedding dimension, and the time delay. We used raw tidal flow time series data, which represent the most basic form of information obtained from the respiratory system with use of our equipment. Consequently, the attractor may be visualized in phase space in its "original" form, before data manipulations such as integration. In addition, the accuracy of the reconstruction of the attractor is dependent on the amount of data included. Attractor reconstructions normally incorporate at least several thousand points (1, 21, 23, 37); hence, short data series of 50–100 breaths (e.g., VT) are insufficient. For our purposes, a time series of 50 breaths consisted of ~15,000–18,000 data points, providing much clearer information regarding the attractor behavior and, therefore, about the global system dynamics.

Reconstruction of the attractor in phase space is dependent on the nature of the individual signal under study. Numerous techniques have been developed to assist with accurate representation of the intrinsic dynamics of the system graphically. It is possible to obtain a suitable value of the embedding delay ({xi}) in several ways (1, 23). In this study, the embedding delay was increased systematically, and the attractor trajectories were viewed for consistency of appearance at different delays. That is, the unfolding of the attractor trajectories was examined at many delays, with an appropriate delay being one that provided the best visual representation of the system dynamics.

Summary

We found an increased variability of tidal breathing after a sigh in healthy infants during quiet sleep. Short-range breath-to-breath memory of the neurorespiratory control system decreased toward the sigh. After the sigh, this memory was recovered before steadily returning to presigh conditions. Long-range memory and the location of the system attractor within phase space were unchanged after the sigh, consistent with the presence of long-range stability. Attractor dynamics were markedly different, with a relatively fixed, locked-in behavior before the sigh changing to a more adaptive pattern with increased fluctuation after the sigh. These findings support the hypothesis that sighs in infants not only act to reset mechanics but also to change the properties of the neurorespiratory control system. Finally, an intriguing finding was the observation that loss of short-range memory appeared to be a predictor of a sigh.

The ability to sigh may be an important mechanism of control of breathing. The response to a sigh may be different in various groups of infants at risk for inadequate control of breathing, such as premature infants, those with neurological impairment, or those at risk for sudden infant death syndrome. Further studies should investigate sigh behavior during maturation and in disease to develop a new sensitive marker of control of breathing in these infants at risk.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study was supported in part by Swiss National Fund Grant 32-68025.02. Partial support was provided to D. N. Baldwin by Swiss Government Grant 2002-0041 and European Respiratory Society Long-Term Research Fellowship 2002-032.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank the parents of the infants who participated in the study and G. Wirz, M. Graf, C. Becher, and H. Gehr for help with enrollment of infants and assistance in performing measurements.


    FOOTNOTES
 

Address for reprint requests and other correspondence: D. N. Baldwin, Pediatric Respiratory Medicine, Univ. Children's Hospital of Berne, Berne CH-3010, Switzerland (E-mail: dn_baldwin{at}hotmail.com).

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 

  1. Addison P. Fractals and Chaos. Bristol, UK: Institute of Physics Publishing, 1997.
  2. Alvarez JE, Bodani J, Fajardo CA, Kwiatkowski K, Cates DB, and Rigatto H. Sighs and their relationship to apnea in the newborn infant. Biol Neonate 63: 139–146, 1993.
  3. Bartlett D Jr. Origin and regulation of spontaneous deep breaths. Respir Physiol 12: 230–238, 1971.
  4. Bassingthwaite JB, Liebovitch LS, and West BJ. Fractal Physiology. New York: Oxford University Press, 1994.
  5. Bates JH, Schmalisch G, Filbrun D, and Stocks J. Tidal breath analysis for infant pulmonary function testing. ERS/ATS Task Force on Standards for Infant Respiratory Function Testing European Respiratory Society/American Thoracic Society. Eur Respir J 16: 1180–1192, 2000.
  6. Bendixon HH, Smith GM, and Mead J. Pattern of ventilation in young adults. J Appl Physiol 19: 195–198, 1964.
  7. Brown RH and Mitzner W. Airway response to deep inspiration: role of inflation pressure. J Appl Physiol 91: 2574–2578, 2001.
  8. Cernelc M, Suki B, Reinmann B, Hall GL, and Frey U. Correlation properties of tidal volume and end-tidal O2 and CO2 concentrations in healthy infants. J Appl Physiol 92: 1817–1827, 2002.
  9. Cleave JP, Levine MR, and Fleming PJ. The control of ventilation: a theoretical analysis of the response to transient disturbances. J Theor Biol 108: 261–283, 1984.
  10. Davis GM and Moscato J. Changes in lung mechanics following sighs in premature newborns without lung disease. Pediatr Pulmonol 17: 26–30, 1994.
  11. Duggan CJ, Chan J, Whelan AJ, and Berend N. Bronchodilatation induced by deep breaths in relation to transpulmonary pressure and lung volume. Thorax 45: 930–934, 1990.
  12. Ferris BG Jr and Pollard DS. Effect of deep and quiet breathing on pulmonary compliance in man. J Clin Invest 39: 143–149, 1960.
  13. Fleming PJ, Goncalves AL, Levine MR, and Woollard S. The development of stability of respiration in human infants: changes in ventilatory responses to spontaneous sighs. J Physiol 347: 1–16, 1984.
  14. Fleming PJ, Levine MR, and Goncalves A. Changes in respiratory pattern resulting from the use of a facemask to record respiration in newborn infants. Pediatr Res 16: 1031–1034, 1982.
  15. Franco P, Verheulpen D, Valente F, Kelmanson I, de Broca A, Scaillet S, Groswasser J, and Kahn A. Autonomic responses to sighs in healthy infants and in victims of sudden infant death. Sleep Med 4: 569–577, 2003.
  16. Frey U, Jackson AC, and Silverman M. Differences in airway wall compliance as a possible mechanism for wheezing disorders in infants. Eur Respir J 12: 136–142, 1998.
  17. Frey U, Stocks J, Coates A, Sly P, and Bates J. Specifications for equipment used for infant pulmonary function testing. ERS/ATS Task Force on Standards for Infant Respiratory Function Testing European Respiratory Society/American Thoracic Society. Eur Respir J 16: 731–740, 2000.
  18. Frey U, Stocks J, Sly P, and Bates J. Specification for signal processing and data handling used for infant pulmonary function testing. ERS/ATS Task Force on Standards for Infant Respiratory Function Testing European Respiratory Society/American Thoracic Society. Eur Respir J 16: 1016–1022, 2000.
  19. Hoch B, Bernhard M, and Hinsch A. Different patterns of sighs in neonates and young infants. Biol Neonate 74: 16–21, 1998.
  20. Issa FG and Porostocky S. Effect of sleep on changes in breathing pattern accompanying sigh breaths. Respir Physiol 93: 175–187, 1993.
  21. Jeong J, Gore JC, and Peterson BS. Detecting determinism in short time series, with an application to the analysis of a stationary EEG recording. Biol Cybern 86: 335–342, 2002.
  22. Jokinen V, Syvanne M, Makikallio TH, Airaksinen KE, and Huikuri HV. Temporal age-related changes in spectral, fractal and complexity characteristics of heart rate variability. Clin Physiol 21: 273–281, 2001.
  23. Kantz H and Schreiber T. Nonlinear Time Series Analysis. Cambridge: Cambridge University Press, 2002.
  24. Khoo MC. Determinants of ventilatory instability and variability. Respir Physiol 122: 167–182, 2000.
  25. Knowlton GC and Larrabee MG. A unitary analysis of pulmonary volume receptors. Am J Physiol 147: 100–114, 1947.
  26. Kosch PC and Stark AR. Dynamic maintenance of end-expiratory lung volume in full-term infants. J Appl Physiol 57: 1126–1133, 1984.
  27. Makikallio TH, Koistinen J, Jordaens L, Tulppo MP, Wood N, Golosarsky B, Peng CK, Goldberger AL, and Huikuri HV. Heart rate dynamics before spontaneous onset of ventricular fibrillation in patients with healed myocardial infarcts. Am J Cardiol 83: 880–884, 1999.
  28. Peng CK, Havlin S, Stanley HE, and Goldberger AL. Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. Chaos 5: 82–87, 1995.
  29. Peng CK, Mietus J, Hausdorff JM, Havlin S, Stanley HE, and Goldberger AL. Long-range anticorrelations and non-Gaussian behavior of the heartbeat. Phys Rev Lett 70: 1343–1346, 1993.
  30. Peng CK, Mietus JE, Liu Y, Lee C, Hausdorff JM, Stanley HE, Goldberger AL, and Lipsitz LA. Quantifying fractal dynamics of human respiration: age and gender effects. Ann Biomed Eng 30: 683–692, 2002.
  31. Perez-Padilla R, West P, and Kryger MH. Sighs during sleep in adult humans. Sleep 6: 234–243, 1983.
  32. Poets CF, Rau GA, Neuber K, Gappa M, and Seidenberg J. Determinants of lung volume in spontaneously breathing preterm infants. Am J Respir Crit Care Med 155: 649–653, 1997.
  33. Prechtl HF. The behavioural states of the newborn infant. Brain Res 76: 185–212, 1974.
  34. Reininger EJ and Segall P. Effect of inhalation of O2 and CO2 gas mixtures on spontaneous gasps and apnea in unaesthetized dog (Abstract). Physiologist 13: 290, 1970.
  35. Revow M, England SJ, O'Beirne H, and Bryan AC. A model of the maturation of respiratory control in the newborn infant. IEEE Trans Biomed Eng 36: 414–423, 1989.
  36. Skinner JE, Pratt CM, and Vybiral T. A reduction in the correlation dimension of heartbeat intervals precedes imminent ventricular fibrillation in human subjects. Am Heart J 125: 731–743, 1993.
  37. Small M, Judd K, Lowe M, and Stick S. Is breathing in infants chaotic? Dimension estimates for respiratory patterns during quiet sleep. J Appl Physiol 86: 359–376, 1999.
  38. Szereda-Przestaszewska M, Bartlett D Jr, and Wise JC. Changes in respiratory frequency and end-expiratory volume accompanying augmented breaths in cats. Pflügers Arch 364: 29–33, 1976.
  39. Takens F. Detecting Strange Attractors in Turbulence. New York: Springer, 1981.
  40. Thach BT and Taeusch HW Jr. Sighing in newborn human infants: role of inflation-augmenting reflex. J Appl Physiol 41: 502–507, 1976.
  41. Ueda Y, Stick SM, Hall GL, and Sly PD. Control of breathing in infants born to smoking mothers. J Pediatr 135: 226–232, 1999.
  42. Wu CH, Chou HC, Hsieh WS, Chen WK, Huang PY, and Tsao PN. Good estimation of arterial carbon dioxide by end-tidal carbon dioxide monitoring in the neonatal intensive care unit. Pediatr Pulmonol 35: 292–295, 2003.
  43. Younes M and Riddle W. A model for the relation between respiratory neural and mechanical outputs. I. Theory. J Appl Physiol 51: 963–978, 1981.
  44. Younes M and Riddle W. Relation between respiratory neural output and tidal volume. J Appl Physiol 56: 1110–1119, 1984.



This article has been cited by other articles:


Home page
Eur Respir JHome page
M. Kohler, S. Kriemler, E. M. Wilhelm, H. Brunner-LaRocca, M. Zehnder, and K. E. Bloch
Children at high altitude have less nocturnal periodic breathing than adults
Eur. Respir. J., July 1, 2008; 32(1): 189 - 197.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
J. L. Fisher and S. S. Margulies
Modeling the effect of stretch and plasma membrane tension on Na+-K+-ATPase activity in alveolar epithelial cells
Am J Physiol Lung Cell Mol Physiol, January 1, 2007; 292(1): L40 - L53.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
97/5/1830    most recent
00298.2004v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (8)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Baldwin, D. N.
Right arrow Articles by Frey, U.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Article