Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The Effects of Involuntary Respiratory Contractions on Cerebral Blood Flow during Maximal Apnoea in Trained Divers

  • Troy J. Cross ,

    t.cross@griffith.edu.au

    Affiliations Griffith Health Institute and Heart Foundation Research Centre, Griffith University, Gold Coast Campus, Queensland, Australia, Division of Cardiovascular Diseases, Mayo Clinic, Rochester, Minnesota, United States of America

  • Justin J. Kavanagh,

    Affiliation Griffith Health Institute and Heart Foundation Research Centre, Griffith University, Gold Coast Campus, Queensland, Australia

  • Toni Breskovic,

    Affiliation Department of Physiology, University of Split School of Medicine, Split, Croatia

  • Petra Zubin Maslov,

    Affiliation Department of Physiology, University of Split School of Medicine, Split, Croatia

  • Mihajlo Lojpur,

    Affiliation Department of Anaesthesiology, Clinical Hospital Center Split, Split, Croatia

  • Bruce D. Johnson,

    Affiliation Division of Cardiovascular Diseases, Mayo Clinic, Rochester, Minnesota, United States of America

  • Zeljko Dujic

    Affiliation Department of Physiology, University of Split School of Medicine, Split, Croatia

Abstract

The effects of involuntary respiratory contractions on the cerebral blood flow response to maximal apnoea is presently unclear. We hypothesised that while respiratory contractions may augment left ventricular stroke volume, cardiac output and ultimately cerebral blood flow during the struggle phase, these contractions would simultaneously cause marked ‘respiratory’ variability in blood flow to the brain. Respiratory, cardiovascular and cerebrovascular parameters were measured in ten trained, male apnoea divers during maximal ‘dry’ breath holding. Intrathoracic pressure was estimated via oesophageal pressure. Left ventricular stroke volume, cardiac output and mean arterial pressure were monitored using finger photoplethysmography, and cerebral blood flow velocity was obtained using transcranial ultrasound. The increasingly negative inspiratory intrathoracic pressure swings of the struggle phase significantly influenced the rise in left ventricular stroke volume (R2 = 0.63, P<0.05), thereby contributing to the increase in cerebral blood flow velocity throughout this phase of apnoea. However, these contractions also caused marked respiratory variability in left ventricular stroke volume, cardiac output, mean arterial pressure and cerebral blood flow velocity during the struggle phase (R2 = 0.99, P<0.05). Interestingly, the magnitude of respiratory variability in cerebral blood flow velocity was inversely correlated with struggle phase duration (R2 = 0.71, P<0.05). This study confirms the hypothesis that, on the one hand, involuntary respiratory contractions facilitate cerebral haemodynamics during the struggle phase while, on the other, these contractions produce marked respiratory variability in blood flow to the brain. In addition, our findings indicate that such variability in cerebral blood flow negatively impacts on struggle phase duration, and thus impairs breath holding performance.

Introduction

A truly ‘maximal’ apnoea is typified by two distinct, temporal phases. The first phase is a quiescent period where respiratory neuromuscular output is voluntarily suppressed, the glottis is closed, and no movement of the chest wall occurs (i.e., ‘easy-going’ phase) [1]. The beginning of the second period, the ‘struggle’ phase, is marked by the onset of involuntary respiratory contractions that increase in both magnitude and frequency until the breaking point of apnoea is reached [2], [3]. Several factors are known to influence the apnoea breakpoint, these factors include: the degree of hypoxaemia and hypercapnia; lung volume and rate of alveolar gas-exchange; and afferent information originating within the diaphragm during contraction [4][7]. In addition to these factors, the ability to maintain adequate delivery of O2 to the cerebral tissues must necessarily impact on breath holding performance [8], [9], insofar as brain hypoxia/hypercapnia may induce cognitive impairment, unconsciousness and, in some cases, may result in fatality [10], [11]. The question arises: What mechanisms exist to ensure an adequate supply of O2 to the cerebral tissues during maximal apnoea?

It is elementary that maximal breath holding should reduce the O2 content of arterial blood [12]. To maintain cerebral O2 delivery in spite of this arterial hypoxaemia, bulk blood flow to the brain progressively increases throughout maximal breath holding [8]. The rise in cerebral blood flow during the easy-going phase of apnoea is mediated by: (i) a hypercapnia-induced vasodilation of the cerebral arterioles [13], [14]; and (ii) a progressive increase in mean arterial pressure (MAP) consequent to massive peripheral vasoconstriction [8], [12], [15]. It is of note that the arterial hypertension during the struggle phase appears solely mediated by an increasing left ventricular stroke volume (LVSV) and thus cardiac output (CO) – there is no observable change in total peripheral resistance (TPR) throughout this phase of breath holding [9], [16]. The mechanisms which contribute to the rise in LVSV during the struggle phase remain uncertain.

Previous investigators have hypothesised that forceful inspiratory contractions may act to improve left ventricular output during the struggle phase [9], [16]. Certainly, sustained negative intrathoracic pressure (ITP) is known to improve right venous return and LVSV [17][21]. And given that TPR remains unchanged during the struggle phase, any increase in left ventricular output would cause a proportional rise in MAP, therefore contributing to the increase in cerebral blood flow. Although this rationale is attractive, it remains to be shown that involuntary respiratory contractions directly influence left ventricular output during the struggle phase of maximal breath holding.

It must be considered that while involuntary respiratory contractions may facilitate haemodynamics in the manner described above, these contractions are likely to cause pronounced ‘respiratory’ variability in cerebral blood flow during the struggle phase [9]. For example, the acute increases and decreases in ITP which occur throughout the respiratory-cycle are known to impair and augment left ventricular output, respectively [22], [23]. This respiratory-modulation of LVSV causes directionally similar fluctuations in CO and MAP [24], [25], ultimately producing respiratory variability in cerebral blood flow [26], [27]. It can be reasoned that these abrupt variations in cerebral blood flow induce an ongoing mismatch between cerebral tissue O2 supply and utilisation which may, in turn, affect cyclic periods of tissue hypoxia in the vicinity of central chemoreceptors. According to this idea, respiratory variability in cerebral blood flow may contribute an additional ventilatory stimulus during apnoea, compounding the individual’s ‘urge to breathe’, directly opposing further apnoeic effort. To date, no study has directly investigated the role of involuntary respiratory contractions in the cerebral blood flow response to maximal breath holding. Moreover, the impact of respiratory variability in cerebral blood flow on apnoea duration (i.e., performance) is unknown.

The first aim of this study was to determine the effects of involuntary respiratory contractions on the cerebral blood flow response to maximal breath holding. It was hypothesised that the increasingly negative ITP swings, generated by inspiratory muscle contraction, would augment LVSV and ultimately contribute to the rise in cerebral blood flow over the struggle phase. At the same time, it was expected that these contractions would produce marked respiratory variability in LVSV, CO, MAP and thus cerebral blood flow during the struggle phase. The second aim of this study was to determine whether cerebral blood flow is a determinant of breath holding performance. It was hypothesised that struggle phase duration would be directly proportional to the rise in bulk cerebral blood flow, yet inversely related to the magnitude of its respiratory variability, throughout this phase of apnoea.

Methods

Ethics Statement

The present study conformed to the principles outlined in the Declaration of Helsinki and was approved by the research ethics board at the School of Medicine, University of Split, Croatia.

Participants and Experimental Design

Ten trained, male apnoea divers (28±2 yr; 183±2 cm, 84±3 kg) volunteered to participate in the study and provided written informed consent. The subjects underwent a pre-participatory health screening to ensure they were physically active non-smokers, with no history of cardiac, pulmonary or metabolic disease. We purposely chose to recruit trained apnoea divers because these individuals are accustomed to extended periods of breath holding. These individuals display, on average, greater absolute durations for the easy-going and struggle phases than naive participants [28].

Subjects first performed a number of pulmonary function tests in the supine position, after which they were instructed to remain supine for 10 mins before performing breath-hold manoeuvres (i.e., rest period). Breath-hold manoeuvres commenced after full inflation to total lung capacity (TLC) while wearing a nose-clip. The subjects were instructed to keep their glottis closed during each breath-hold. Subjects were asked to refrain from performing preparatory hyperventilation, and were not allowed to perform glossopharyngeal insufflations prior to the manoeuvre. Subjects completed two practice trials, followed by 2–3 experimental breath holds, separated by at least 10 mins of supine rest. The subjects were encouraged to completely relax their respiratory musculature during the early period of the breath-hold (i.e., easy-going phase) and, should the urge arise, to allow respiratory contractions to develop “naturally”. The subjects were also instructed to avoid contraction of their peripheral musculature, including neck and facial muscles. The latter point was particularly important for obtaining high-quality measurements of middle cerebral artery blood flow velocity using the transcranial Doppler probe (see below). In the event that peripheral muscle activation was observed by the investigator, or reported by the subject, breath holding was immediately terminated. Subjects were excluded from the study if, after a further two attempts, a satisfactory apnoea could not be obtained.

Pulmonary Function, Arterial O2 Saturation and Respiratory Pressures

Before each experiment, the subjects performed a forced vital capacity (FVC) manoeuvre while in the supine posture (Quark PFT, Cosmed, Rome, Italy). All pulmonary function testing was performed and reported in accordance with the American Thoracic Society guidelines [29]. Arterial O2 saturation of haemoglobin (SaO2) was measured via finger pulse oximetry (Poet II, Criticare Systems, Waukesha, WI). Oesophageal pressure was measured using a latex balloon-tip catheter (Ackrad Laboratories, CooperSurgical, Trumbull, CT, USA) advanced through the nose into the oesophagus [30]. The catheter balloon was then inflated with 1 ml of air, and the “occlusion” test was performed to ensure correct placement of the catheter in lower third of the oesophagus [31]. The catheter was connected to a differential pressure transducer (PX138-005D5V, Omega Engineering Inc., Stamford, CT, USA). The transducer was calibrated before each test using a digital manometer (HHP-90, Omega Engineering inc., Stamford, CT, USA).

Arterial Blood Pressure and Cerebral Blood Flow Velocity

Beat-by-beat arterial blood pressure was measured using a pneumatic cuff placed around the middle phalanx of the non-dominant hand, and was connected to a photoplethysmograph (Finometer, Finapress Medical Systems, Arnhem, Netherlands). The hand bearing the finger cuff was positioned at the level of the heart in order to negate hydrostatic pressure artifact. The hand was kept in this position for the duration of the experimental protocol. Cerebral blood flow velocity (CBFV) was obtained via transcranial Doppler ultrasound of the proximal middle cerebral artery, measured using a 2-MHz probe fixed at a constant angle over the right posterior temporal “window” (Transcranial Doppler, Multigon, Yonkers, NY, USA).

Data Analysis and Processing

Oesophageal pressure, finger arterial pressure, and transcranial Doppler flow velocity signals were sampled continuously at 1000 Hz (Powerlab 16SP, ADInstruments Inc., Castle Hill, Australia) and stored on a personal computer for off-line analyses. The oesophageal pressure signal was taken to represent intrathoracic pressure (ITP). The beginning of the struggle phase was identified as the first swing in ITP greater than ±2 cmH2O of the mean ITP calculated from the previous 5 s of data. The magnitudes of static inspiratory and expiratory contractions were determined as the minima and maxima of the ITP tracing, respectively. The total time elapsed between two consecutive inspiratory efforts (Ttot) was used to calculate the frequency of respiratory contraction (fR = 60/Ttot).

Mean arterial pressure (MAP) and CBFV were computed as the arithmetic means of the arterial pulse wave and transcranial Doppler flow velocity signals over each beat interval, respectively [32]. Left ventricular stroke volume (LVSV) was computed from the arterial pulse wave using the improved Modelflow algorithm [33]. Cardiac output (CO) was taken as the product of heart rate (HR) and LVSV, and total peripheral resistance (TPR) as MAP/CO. The cerebrovascular resistance index (CVRi) was taken as MAP/CBFV. The recorded ITP signal was decimated to 50 Hz. The beat-by-beat time series of HR, LVSV, CO, MAP, TPR, CBFV and CVRi were resampled to 50 Hz using cubic spline interpolation. These haemodynamic variables are under constant, dynamic autonomic control; processes that are most effective within the low-frequency range (i.e., 0.04–0.15 Hz) [27], [34]. Furthermore, these variables are also affected either directly or indirectly by respiration – such modulation of haemodynamics occurs within the broadband of respiratory frequencies (i.e., >0.15 Hz). To extract the variability in each dataset that is due solely to respiratory effort, all time series were band-pass filtered using an 8th order, zero-phase Butterworth filter with the low and high cut-off frequencies set equal to 0.20 and 0.80 Hz. This specific range of frequencies encompassed the upper- and lower-most values for fR observed during the experiment in all subjects.

The above parameters were averaged into seven time-epochs representing: the resting period; and the first, second and third tertiles of the easy-going phase (EP1, EP2 and EP3) and the struggle phase durations (SP1, SP2 and SP3), respectively. The means () of each time series were computed using raw data, while standard deviations (SD) were calculated using the filtered (‘processed’) time series (Figure 1). The respiratory variability in each time-series was quantified by the coefficient of variation: CV (%) = 100(SD/). The absolute SD of the ITP time-series was used as an index of respiratory variability, because CV cannot be computed for data where may approach zero. Lastly, respiratory variability was not calculated for TPR and CVRi. We acknowledge that previous investigators have reported high-frequency variability in TPR and CVRi in humans [32]. However, within the scope of the present research design, we could not properly discern whether such variability in TPR and CVRi was of physiological origin or simply reflected computational errors induced by slight temporal delays between measurement devices (e.g., MAP measured via the photoplethysmograph and CBFV obtained via the Doppler ultrasound probe).

thumbnail
Figure 1. Example of left ventricular stroke volume time-series of a representative subject.

The left and right panels display data obtained while the subject was resting supine, and during maximal apnoea, respectively. LVSV: left ventricular stroke volume. The vertical dashed line represents the start of the struggle phase. The raw LVSV tracing was obtained by re-sampling beat-by-beat values to 50 Hz via cubic spline interpolation. The interpolated data was band-pass filtered using an 8th order, zero-phase Butterworth filter with low and high cut-off frequencies of 0.20 and 0.80 Hz (i.e., the broadband of respiratory frequencies). This ‘processed’ data was taken to represent the respiratory variability in LVSV.

https://doi.org/10.1371/journal.pone.0066950.g001

Statistical Analyses

All data were averaged across repeated trials. One-way analyses of variance were used to evaluate the main-effect of time (rest, EP1, EP2, EP3, SP1, SP2 and SP3) on respiratory, cardiovascular and cerebrovascular variables. Pair-wise comparisons were assessed using the Bonferroni post-hoc adjustment. The role of inspiratory and expiratory ITP swings in the augmentation of LVSV during the struggle phase was examined using stepwise multiple linear regression. Path analysis was used to determine the extent to which respiratory contractions affect respiratory variations in CBFV during the struggle phase (see below). Lastly, the influence of cerebral blood flow on breath holding performance was assessed using stepwise multiple linear regression, where the overall rise and magnitude of respiratory variability in CBFV during the struggle phase were input into the regression model as independent variables, while the duration of the struggle phase (in seconds) was input as the dependent variable. Results are presented as mean ± standard error of the mean (SEM) unless stated otherwise. All data were analysed using SPSS 20.0 (SPSS, Inc., Chicago, IL, USA). Statistical analyses were considered significant if P<0.05.

Path Analysis

Path analysis is used to determine whether a set of multivariate data are consistent with a conceptual structural model [35], [36]. The structural model theorises a number of causal relationships between variables, which may be either direct or indirect in nature. The influence of one variable on another is determined by performing a series of multiple and simple linear regression analyses which, in turn, yield a set of standardised ß-coefficients for each directed path [35][37]. A ß-coefficient of –0.50 indicates that for every 1 SD increase in variable X1 there occurs a 0.50 SD decrease in variable Y. Conversely, a ß-coefficient of 1.25 indicates that for every 1 SD increase in variable X1 there occurs a 1.25 SD increase in variable Y. In the example model displayed in Figure 2, ßX1.Y represents the direct effect of X1 on Y. It is also possible that the effect of X1 on Y is mediated by other intervening variables separate from the direct path X1→Y. For example, it may be that: X1→X2→X3→Y. The total indirect effect of X1 on Y is then quantified by summing the products of ß-coefficients for all mediating or intervening paths (e.g., ßX1.X2ßX2.X3ßX3.Y). The total impact of X1 on the terminal dependent variable Y is given by the sum of its direct and indirect effects. Path analysis was conducted using SPSS AMOS 20.0 (SPSS, Inc., Chicago, IL, USA). The causal model used in this analysis was constructed a priori, using the rationale outlined below.

thumbnail
Figure 2. Path analysis.

: mean variability observed in the corresponding parameter during the struggle phase; ITP: intrathoracic pressure; HR: heart rate; LVSV: left ventricular stroke volume; CO: cardiac output; MAP: mean arterial blood pressure; TPR: total peripheral resistance; CVRi: cerebrovascular resistance index; CBFV: cerebral blood flow-velocity. The top panel depicts a hypothetical Path model to illustrate the direct and indirect effects of each parameter (X1, X2 and X3) on the terminal endogenous variable (i.e., Y). The symbol ß represents the standardised path coefficient describing the influence of one variable on another. See Methods for further explanation. The hypothesised causal model used to examine the influence of respiratory contractions on the respiratory variability in cerebral blood flow velocity is displayed in the middle-panel. The ‘reduced-fit’ causal model is presented in the bottom-panel. Only those parameters which were identified as significant (P<0.05) mediators of the terminal dependent variable (i.e., ) are retained in this panel. Solid arrows denote the direction of causality between corresponding parameters. Values located adjacent to arrows represent the standardised coefficient of the corresponding path. Values displayed in bold represent the coefficient of determination (R2) of the relationship between the corresponding parameter and all preceding paths.

https://doi.org/10.1371/journal.pone.0066950.g002

The multivariate data used in the hypothesised causal model depicted in Figure 2 were the mean variability () in respiratory, cardiovascular and cerebrovascular parameters during the struggle phase. The majority of the paths outlined in this model were justified for inclusion because of the well-known physiological dependencies between parameters. The remaining paths on the diagram were justified based on the following rationale: Path to : during tidal inspiration, the fall in ITP abruptly impairs left ventricular output and decreases LVSV, while the expiratory rise in ITP augments LVSV [22], [23]. These “within-breath” changes represent the respiratory modulation of LVSV in healthy individuals; Path to : it is known that some of the variability in spontaneous HR is attributable to central respiratory rhythm (respiratory sinus arrhythmia) [24], [25] and may be evidenced by synchronous variations in ITP and HR; Path to : changes in ITP may be transmitted through the wall of the thoracic aorta, directly modulating arterial blood pressure [38].

Results

Pulmonary Function, Arterial O2 Saturation and Breath Holding Performance

The resting pulmonary function of our participants appeared within the normal range, where FVC, forced expiratory volume in 1 s (FEV1) and FEV1/FVC were 6.57±0.36 L (120±5% predicted), 5.40±0.25 L (118±4% predicted) and 0.82±0.01 (100±1% predicted), respectively. The subjects’ average breath-hold duration was 187±7 s. The coefficient of variation (CV) for total breath-hold duration was 10±2%, indicating minor variability between repeated breath holdings efforts. The mean durations of the easy-going and struggle phases were 96±4 s (51±2% of total) and 91±5 s (49±2% of total), respectively. All subjects were able to produce at least two satisfactory apnoeic efforts without discernable activation of their peripheral musculature. Moreover, our subjects did not experience loss of consciousness, or overt neurological impairment during and after maximal apnoea.

All subjects displayed significant (P<0.05) arterial desaturation during the breath-hold manoeuvre, where SaO2 values decreased from rest (99±1%) to a nadir (87±2%) which occurred ∼30 s after the end of breath holding – this delay represents the circulatory transit time from the lungs to the finger where the oximeter probe was located.

Respiratory, Cardiovascular and Cerebrovascular Parameters During Breath Holding

The time-course changes in respiratory and haemodynamic variables during maximal breath holding are displayed in Figure 3. HR did not change significantly from resting values through breath holding. LVSV and CO decreased below resting values at the beginning of the easy-going phase (P<0.05). After this initial decline, LVSV and CO progressively increased with time (P<0.05) reaching values similar to those observed at rest by the breakpoint. MAP and CBFV dropped below resting values at the start of the easy-going phase (P<0.05), rising with time thereafter until the apnoea breakpoint (P<0.05). TPR increased abruptly at the onset of breath holding, achieving its highest value by the end of the easy-going phase (P<0.05), whereafter TPR remained stable above baseline until the end of breath holding (P<0.05). CVRi was notably higher than resting values at the start of the easy-going phase (P<0.05), but progressively decreased from this point until the middle of the struggle phase (P<0.05).

thumbnail
Figure 3. Time-course changes in respiratory, cardiovascular and cerebrovascular variables at rest and during maximal apnoea.

Vertical error bars represent means ± SEM. EP1, EP2 and EP3: first, second and third tertiles of the easy-going phase; SP1, SP2 and SP3: first second and third tertiles of the struggle phase; ITP: intrathoracic pressure; HR: heart rate; LVSV: left ventricular stroke volume; CO: cardiac output; MAP: mean arterial blood pressure; CBFV: cerebral blood flow-velocity; TPR: total peripheral resistance in units of mmHg⋅L−1⋅min; CVRi: cerebrovascular resistance index in units of mmHg⋅cm−1⋅s. The closed and open circles in the top left-hand panel represent expiratory and inspiratory ITP swings, respectively. *Significantly different from previous time-point, P<0.05. Significantly different from resting value, P<0.05.

https://doi.org/10.1371/journal.pone.0066950.g003

At the onset of involuntary respiratory contractions, inspiratory ITP swings became more negative with time from the beginning until the end of the struggle phase (P<0.05). Expiratory ITP swings significantly increased from the middle of the struggle phase until the apnoea breakpoint (P<0.05). The frequency of contractions (fR) increased from the beginning until the end of the struggle phase (SP1 = 20±5 efforts•min−1, SP2 = 25±4 efforts•min−1, SP3 = 34±4 efforts•min−1; All P-values were <0.05). The values for fR were significantly higher during this phase of breath holding than those observed at rest (12±4 efforts•min−1, P<0.05).

The Influence of Intrathoracic Pressure on Left Ventricular Stroke Volume During the Struggle Phase

Stepwise multiple regression revealed that the increasing magnitude of inspiratory (not expiratory) ITP swings influenced the rise in LVSV during the struggle phase of maximal breath holding (Figure 4, R2 = 0.63, P<0.01). The linear equation representing this relationship was: ΔLVSV = –0.565•ΔITPinsp +0.580, where ΔITPinsp is the change in the magnitude of inspiratory ITP swings from the beginning until the end of the struggle phase.

thumbnail
Figure 4. Influence of inspiratory intrathoracic pressure on left ventricular stroke volume during the struggle phase.

LVSV: left ventricular stroke volume; Δ: absolute change in values from beginning until the end of the struggle phase of breath holding; ITPinsp: inspiratory swing in intrathoracic pressure during involuntary respiratory contractions. The solid line represents the line of best-fit determined from simple linear regression. The hashed lines denote the ±95% confidence interval of the regression function.

https://doi.org/10.1371/journal.pone.0066950.g004

The‘Respiratory’ Variability in Respiratory and Haemodynamic Parameters During Breath Holding

The respiratory variability in all parameters was systematically lower during the easy-going phase compared with resting values (Figure 5; P<0.05). Thereafter, the variability in each parameter increased from the start until the middle of the struggle phase (P<0.05) – was the only variable that continued to rise until the apnoea breakpoint (P<0.05). The overall coefficient of determination (R2) for the entire ‘reduced-fit’ causal model displayed in Figure 2 was 99% (P<0.05). The direct, indirect and total effects of each parameter on are reported in Table 1. It was revealed that exerted a modest and positive influence on the respiratory variability in cerebral blood flow during the struggle phase (). The indirect path through which affected was mediated by , , and finally, . exerted the least overall influence on .

thumbnail
Figure 5. Respiratory variability in intrathoracic pressure, cardiovascular and cerebrovascular parameters at rest and during maximal apnoea.

Vertical error bars represent mean coefficients of variation (CV) ± SEM. EP1, EP2 and EP3: first, second and third tertiles of the easy-going phase; SP1, SP2 and SP3: first second and third tertiles of the struggle phase; ITP: intrathoracic pressure; HR: heart rate; LVSV: left ventricular stroke volume; CO: cardiac output; MAP: mean arterial blood pressure; CBFV: cerebral blood flow velocity; TPR: total peripheral resistance; CVRi: cerebrovascular resistance index. aThe standard deviation of ITP is presented here as a measure of respiratory variability instead of CV (see Methods). *Significantly different from previous time-point, P<0.05. Significantly different from resting value, P<0.05.

https://doi.org/10.1371/journal.pone.0066950.g005

thumbnail
Table 1. The impact of‘respiratory’ variability in intrathoracic pressure and central haemodynamics on cerebral blood flow velocity during the struggle phase.

https://doi.org/10.1371/journal.pone.0066950.t001

The Influence of Cerebral Blood Flow on Struggle Phase Duration

Stepwise multiple regression revealed that the rise in bulk cerebral blood flow (ΔCBFV) and during the struggle phase significantly influenced the duration of this phase of apnoea itself (Figure 6; R2 = 0.97, P<0.01). The equation describing this multiple linear relationship was: SPD = (2.12•ΔCBFV) – (16.19•) +99.27, where SPD is the duration of the struggle phase in seconds. Semi-partial correlation coefficients were obtained from the stepwise multiple regression. These coefficients represent the strength of influence of one independent variable on the dependent variable, after the variance due to other independent variables has been partialled out [39]. The semi-partial correlation coefficients for the input parameters ΔCBFV and were 0.507 and –0.433, respectively.

thumbnail
Figure 6. The influence of cerebral blood flow velocity on the struggle phase duration.

Δ: absolute change in values from beginning until the end of the struggle phase of breath holding; CBFV: cerebral blood flow velocity; : mean variability observed in the corresponding parameter during the struggle phase. The R2 values presented in these panels are calculated from each independent variable’s zero-order correlation coefficient. The solid lines represents the line of best-fit determined from simple linear regressions. The hashed lines denote the ±95% confidence intervals of the regression functions.

https://doi.org/10.1371/journal.pone.0066950.g006

Discussion

The present study demonstrates that involuntary respiratory contractions significantly influence the rise in CBFV via facilitation of LVSV during the struggle phase of breath holding. At the same time, these contractions generate abrupt oscillations in central haemodynamics, ultimately producing marked ‘respiratory’ variability in CBFV during the struggle phase (i.e., ). Moreover, this study is the first to provide evidence that a greater magnitude of impairs breath holding performance in man. These findings provide further insight into the effects of involuntary respiratory contractions on the cerebral blood flow response to maximal apnoea, and on breath holding performance.

Cerebral Blood Flow at the Onset of Maximal Breath Holding

LVSV, CO and MAP together decreased at the onset of maximal breath holding (Figure 3), consistent with previous observations [8], [9], [16]. Two possible reasons exist which may explain this abrupt impairment of central haemodynamics at the onset of maximal apnoea. Firstly, the large, negative inspiratory ITP swing generated immediately prior to breath holding at TLC may transiently enhance right ventricular preload and, subsequently, reduce LVSV via ventricular interdependence [40]. Secondly, after breath holding commences, the positive elastic “recoil” pressures exerted by the lungs and chest wall at TLC [41] may increase left ventricular afterload, thereby continuing to impair central haemodynamics throughout the early moments of apnoea (i.e., EP1) [42]. Importantly, we observed that CVRi increased from rest to the beginning of maximal apnoea, possibly in an effort to maintain CBFV despite the marked fall in MAP (i.e., cerebral autoregulation) [27]. Nevertheless, the magnitude of rise in CVRi at this time did not compensate for the observed reduction in arterial ‘driving’ pressure and, consequently, CBFV immediately decreased below resting values at the start of apnoea.

Cerebral Blood Flow During the Easy-going Phase

The asphyxia induced by prolonged apnoea results in vasodilation of the cerebral vasculature and a mounting arterial hypertension during the easy-going phase of breath holding [8], [13][15]. These haemodynamic responses together favour a continuous rise in cerebral blood flow [8], [14]. The arterial hypertension during the easy-going phase is primarily due to chemoreflex activation of the sympathetic nervous system, resulting in a massive vasoconstriction of the peripheral vasculature [15], [43]. Our data are in agreement with the above, wherein the subjects of the present study displayed a progressive decline in CVRi and a continuous rise in TPR, MAP and CBFV over the easy-going phase. It is of interest to note that LVSV systematically increased from the beginning until the end of the easy-going phase in our subjects. This progressive rise in LVSV affected an overall increase in systemic CO, thereby contributing to the observed arterial hypertension. It has been hypothesised that the massive peripheral vasoconstriction during maximal apnoea promotes a central-shift in blood volume, which may augment right venous return and, subsequently, left ventricular output [8], [44]. Notwithstanding this potential increase in central blood volume, the contribution of CO to the mounting arterial hypertension during the easy-going phase was relatively minor compared to the direct effect on MAP exerted by peripheral vasoconstriction.

Cerebral Blood Flow During the Struggle Phase

In the present study, CBFV continued to increase over the duration of the struggle phase, as noted previously by others [8], [16]. It has long been known that arterial CO2 is a potent regulator of bulk cerebral blood flow, wherein arterial hypercapnia and hypocapnia elicit cerebral vasodilation and vasoconstriction, respectively [45]. It should be noted that cerebral arterioles are more responsive to changes in arterial CO2 than to sympathetic activation [46]. Thus, while measurements of arterial CO2 were not obtained in the present study, its effects on bulk cerebral blood flow may be inferred via changes in CVRi. On this point, the fall in CVRi during the struggle phase was relatively minimal compared with the rise in MAP (i.e., –1% v 18%). Thus, we propose that MAP was the primary mediator of the rise in bulk cerebral blood flow during this phase of breath holding.

The mounting arterial hypertension during the struggle phase could not be explained by a further increase in TPR – a finding consistent with that previously reported from our laboratory [9], [16]. It is possible that the magnitude of peripheral tissue hypoxia and/or hypercapnia at this time directly opposed further vasoconstriction, despite increasing sympathetic vasoconstrictor drive [15], [43]. Accordingly, the rising MAP and, by extension, CBFV were primarily due to an increasing systemic CO. Previous investigators have hypothesised that involuntary respiratory contractions may, at least in part, explain the progressive augmentation of left ventricular output during the struggle phase [9], [16].

Due to the location of the heart within the thorax, changes in ITP are directly transmitted to the right atria [18], [19]. In the present study, the increasingly negative inspiratory ITP swings during the struggle phase most likely decreased right atrial pressure, favouring venous return to the right heart [17]. Certainly, recent evidence suggests that venous blood flow through the inferior vena cava abruptly increases at the onset of respiratory contractions during breath holding [16]. Our data are consistent with the above hypothesis in that 63% of the variance observed in ΔLVSV could be accounted for by the increasing negativity of inspiratory efforts during the struggle phase (i.e., ΔITPinsp). Stated in other words, the more negative that inspiratory ITP swings became during the struggle phase, the larger was the increase in LVSV (see Figure 3). To the best of our knowledge, this is the first study to clearly demonstrate that the augmenting CO during the struggle phase is related to a ‘respiratory-mediated’ facilitation of left ventricular output.

‘Respiratory’ Variability in Cerebral Blood Flow During the Struggle Phase

During tidal inspiration, the increasing negativity of ITP transiently impairs left ventricular filling (decreased preload) and emptying (increased afterload), thereby causing an abrupt fall in LVSV – the converse is true during tidal expiration [22], [23]. Such respiratory-synchronous modulation of left ventricular output produces respiratory variability in MAP and cerebral blood flow [24][27]. Given that respiratory effort directly modulates ITP, it is perhaps not surprising that , , abruptly increased at the beginning of the struggle phase. Moreover, Path analysis revealed that indirectly mediated via throughout the struggle phase. These results provide strong evidence that our values for represented the haemodynamic effects of involuntary respiratory contractions, and were not influenced by other physiological processes operating in the same passband of frequencies (i.e., 0.20 to 0.80 Hz).

Cerebral Blood Flow and Breath Holding Performance

It can be reasoned that respiratory variability in cerebral blood flow and, by extension, O2 delivery during this phase induces transient mismatching between cerebral tissue O2 supply and utilisation. Cerebral tissue oxygenation would therefore vary in phase with ITP, evidencing regional tissue hypoxia and hyperoxia during static inspiratory and expiratory efforts, respectively. Recent observations made by Dujic et al. [9] support this thesis. These investigators reported that phasic oscillations in cerebral tissue oxygenation spontaneously emerge at the onset of involuntary respiratory contractions during breath holding in elite apnoea divers.

One may speculate that repetitive periods of cerebral hypoperfusion may affect local hypoxia in the vicinity of central chemoreceptors, constituting a ventilatory stimulus that is eventually perceived as breathlessness [47], [48]. Accordingly, those individuals who display greater respiratory variability in cerebral blood flow may also experience a heightened sensation of breathlessness during the struggle phase. The augmented dyspnoea would directly oppose apnoeic effort in these individuals, impairing breath holding performance. An alternative explanation for the inverse relation between and struggle phase duration is that increasing respiratory pressure-development, as occurs during this phase of breath holding [2], [3], necessitates an increase in energetic expenditure of the respiratory muscles. In turn, the progressively augmenting respiratory muscle O2 consumption may therefore compound the fall in O2 delivery to the cerebral tissues during the struggle phase, impairing overall breath holding performance in a manner that is proportional to . Nevertheless, we demonstrate here, for the first time, that struggle phase duration is directly proportional to ΔCBFV yet inversely related to . Moreover, the semi-partial R2 values for ΔCBFV and were 26% and 19%, respectively, indicating that the independent effects of each parameter on struggle phase duration were similar in strength, though opposite in direction. In turn, it may be stated that while a greater rise in bulk cerebral blood flow during the struggle phase facilitates longer breath holding, a larger amount of respiratory variability in cerebral blood flow negatively impacts on breath holding performance in man. Further studies are required to validate the cause-and-effect nature of this relationship.

Methodological Considerations

Transcranial Doppler ultrasound of the right middle cerebral artery (MCA) was used in the present study to estimate bulk cerebral blood flow during maximal breath holding [49]. This technique does not measure volumetric blood flow per se, rather it quantifies the velocity of blood passing through the insonated artery [50]. Changes in the diameter of the insonated artery may therefore alter blood velocity, independent of volumetric flow through the MCA. However, the weight of available evidence suggests that MCA diameter is relatively stable over a wide range of MAP and arterial blood gases [51], [52]. To maximise the length of apnoea in our subjects, each breath-hold manoeuvre was initiated from total lung capacity. The attendant lung hyperinflation created difficulty in obtaining beat-by-beat measurements of LVSV via Doppler flowmetry of the ascending aorta. Therefore, we obtained LVSV measurements using finger photoplethysmography and the Modelflow algorithm [33]. It was not feasible to calibrate the Finometer against a ‘gold-standard’ measure of CO in our study and, consequently, the absolute magnitudes of the haemodynamic parameters reported herein may be inaccurate [53]. This limitation was circumvented by using a repeated-measures design, where the errors introduced from uncalibrated Modelflow simulations are accounted for on a subject-by-subject basis. In addition to the above concerns, the physical constraints imposed by the Doppler probe headset precluded the measurement of cerebral tissue oxygenation via near-infrared spectroscopy [9].

Finally, the present study examined the impact of involuntary respiratory contractions on the cerebral haemodynamic response to maximal apnoea in trained divers under “dry” laboratory conditions. Therefore, we caution that our findings may not readily generalise to untrained/naïve participants, or to the conditions imposed on subjects when breath holding under water at depth. On this latter point, it is currently not feasible to use transcranial Doppler ultrasound to measure CBFV under water. Future studies may overcome these limitations via computational modelling of the pressure-flow characteristics of peripheral and cerebral vessels [54][56].

Summary of Findings

The findings of our study demonstrate that involuntary respiratory contractions play an important, yet complex role in mediating the CBFV response to maximal breath holding. On the one hand, involuntary respiratory contractions seemingly enhance central and cerebral haemodynamics during the struggle phase. While on the other hand, these contractions affect marked respiratory variability in cerebral blood flow and, presumably, O2 delivery to the brain. Most importantly, our data indicates that larger amounts of respiratory variability in cerebral blood flow during the struggle phase negatively impacts on breath holding performance.

Acknowledgments

The Authors thank each subject for their participation in the present study. The Authors are also grateful for the guidance and advice provided by Drs Surendran Sabapathy and Sophie Lalande during the preparation of this manuscript.

Author Contributions

Conceived and designed the experiments: TJC BDJ ZD. Performed the experiments: TJC BDJ ZD TB ML PZM. Analyzed the data: TJC JJK. Contributed reagents/materials/analysis tools: TJC BDJ ZD. Wrote the paper: TJC JJK ZD BDJ.

References

  1. 1. Lin YC, Lally DA, Moore TO, Hong SK (1974) Physiological and conventional breath-hold breaking points. Journal of Applied Physiology 37: 291–296.
  2. 2. Cross TJ, Breskovic T, Sabapathy S, Zubin Maslov P, Johnson BD, et al.. (2012) Respiratory Muscle Pressure-Development during Breath Holding in Apnoea Divers. Medicine and science in sports and exercise.
  3. 3. Whitelaw WA, McBride B, Amar J, Corbet K (1981) Respiratory neuromuscular output during breath holding. Journal of Applied Physiology 50: 435–443.
  4. 4. Godfrey S, Campbell EJ (1968) The control of breath holding. Respiration Physiology 5: 385–400.
  5. 5. Banzett RB, Lansing RW, Brown R, Topulos GP, Yager D, et al. (1990) 'Air hunger' from increased PCO2 persists after complete neuromuscular block in humans. Respiration Physiology 81: 1–17.
  6. 6. Parkes MJ (2006) Breath-holding and its breakpoint. Experimental Physiology 91: 1–15.
  7. 7. Breskovic T, Lojpur M, Maslov PZ, Cross TJ, Kraljevic J, et al. (2012) The influence of varying inspired fractions of O2 and CO2 on the development of involuntary breathing movements during maximal apnoea. Respiratory Physiology & Neurobiology 181: 228–233.
  8. 8. Palada I, Obad A, Bakovic D, Valic Z, Ivancev V, et al. (2007) Cerebral and peripheral hemodynamics and oxygenation during maximal dry breath-holds. Respiratory Physiology & Neurobiology 157: 374–381.
  9. 9. Dujic Z, Uglesic L, Breskovic T, Valic Z, Heusser K, et al. (2009) Involuntary breathing movements improve cerebral oxygenation during apnea struggle phase in elite divers. Journal of Applied Physiology 107: 1840–1846.
  10. 10. Craig AB Jr (1976) Summary of 58 cases of loss of consciousness during underwater swimming and diving. Medicine and science in sports 8: 171–175.
  11. 11. Albert B, Craig J (1961) Causes of loss of consciousness during underwater swimming. Journal of Applied Physiology 16: 583–586.
  12. 12. Andersson JPA, Linér MH, Rünow E, Schagatay EKA (2002) Diving response and arterial oxygen saturation during apnea and exercise in breath-hold divers. Journal of Applied Physiology 93: 882–886.
  13. 13. Przyby<$>\scale60%\raster="rg1"<$>owski T, Bangash M-F, Reichmuth K, Morgan BJ, Skatrud JB, et al. (2003) Mechanisms of the cerebrovascular response to apnoea in humans. The Journal of Physiology 548: 323–332.
  14. 14. Kjeld T, Pott FC, Secher NH (2009) Facial immersion in cold water enhances cerebral blood velocity during breath-hold exercise in humans. Journal of Applied Physiology 106: 1243–1248.
  15. 15. Heusser K, Dzamonja G, Tank J, Palada I, Valic Z, et al. (2009) Cardiovascular Regulation During Apnea in Elite Divers. Hypertension 53: 719–724.
  16. 16. Palada I, Bakovic D, Valic Z, Obad A, Ivancev V, et al. (2008) Restoration of hemodynamics in apnea struggle phase in association with involuntary breathing movements. Respiratory Physiology & Neurobiology 161: 174–181.
  17. 17. Pinsky MR (1984) Determinants of pulmonary arterial flow variation during respiration. Journal of Applied Physiology 56: 1237–1245.
  18. 18. Wise RA, Robotham JL, Summer WR (1981) Effects of spontaneous ventilation on the circulation. Lung 159: 175–186.
  19. 19. Brecher GA, Hubay CA (1955) Pulmonary blood flow and venous return during spontaneous respiration. Circulation Research 3: 210–214.
  20. 20. Virolainen J, Ventila M, Turto H, Kupari M (1995) Influence of negative intrathoracic pressure on right atrial and systemic venous dynamics. European Heart Journal 16: 1293–1299.
  21. 21. Kilburn KH, Sieker HO (1960) Hemodynamic effects of continuous positive and negative pressure breathing in normal man. Circulation Research 8: 660–669.
  22. 22. Innes JA, De Cort SC, Kox W, Guz A (1993) Within-breath modulation of left ventricular function during normal breathing and positive-pressure ventilation in man. The Journal of Physiology 460: 487–502.
  23. 23. Guz A, Innes JA, Murphy K (1987) Respiratory modulation of left ventricular stroke volume in man measured using pulsed Doppler ultrasound. The Journal of Physiology 393: 499–512.
  24. 24. Toska K, Eriksen M (1993) Respiration-synchronous fluctuations in stroke volume, heart rate and arterial pressure in humans. The Journal of Physiology 472: 501–512.
  25. 25. Elstad M, Toska K, Chon KH, Raeder EA, Cohen RJ (2001) Respiratory sinus arrhythmia: opposite effects on systolic and mean arterial pressure in supine humans. The Journal of Physiology 536: 251–259.
  26. 26. Coughtrey H, Rennie JM, Evans DH, Cole TJ (1993) Factors associated with respiration induced variability in cerebral blood flow velocity. Archives of Disease in Childhood 68: 312–316.
  27. 27. Zhang R, Zuckerman JH, Giller CA, Levine BD (1998) Transfer function analysis of dynamic cerebral autoregulation in humans. American Journal of Physiology - Heart and Circulatory Physiology 274: H233–H241.
  28. 28. Baković D, Valic Z, Eterović D, Vuković I, Obad A, et al. (2003) Spleen volume and blood flow response to repeated breath-hold apneas. Journal of Applied Physiology 95: 1460–1466.
  29. 29. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, et al. (2005) Standardisation of spirometry. European Respiratory Journal 26: 319–338.
  30. 30. Milic-Emili J, Mead J, Turner JM, Glauser EM (1964) Improved Technique for Estimating Pleural Pressure from Esophageal Balloons. Journal of Applied Physiology 19: 207–211.
  31. 31. Baydur A, Cha EJ, Sassoon CS (1987) Validation of esophageal balloon technique at different lung volumes and postures. Journal of Applied Physiology 62: 315–321.
  32. 32. O'Leary DD, Shoemaker JK, Edwards MR, Hughson RL (2004) Spontaneous beat-by-beat fluctuations of total peripheral and cerebrovascular resistance in response to tilt. American Journal of Physiology - Regulatory, Integrative and Comparative Physiology 287: R670–R679.
  33. 33. Jellema WT, Wesseling KH, Groeneveld JAB, Stoutenbeek CP, Thijs LG, et al. (1999) Continuous Cardiac Output in Septic Shock by Simulating a Model of the Aortic Input Impedance: A Comparison with Bolus Injection Thermodilution. Anesthesiology 90: 1317–1328.
  34. 34. Malliani A, Pagani M, Lombardi F, Cerutti S (1991) Cardiovascular neural regulation explored in the frequency domain. Circulation 84: 482–492.
  35. 35. Boker SM, McArdle JJ (2005) Path Analysis and Path Diagrams. Encyclopedia of Statistics in Behavioral Science: John Wiley & Sons, Ltd.
  36. 36. Blunch N (2008) Introduction to Structural Equation Modelling Using SPSS and AMOS. London, England: SAGE Publications, Ltd.
  37. 37. Shipley B (2000) Cause & Correlation in Biology : A User's Guide to Path Analysis, Structural Equations & Causal Inference. Port Chester, NY, USA: Cambridge University Press.
  38. 38. Rudikoff MT, Maughan WL, Effron M, Freund P, Weisfeldt ML (1980) Mechanisms of blood flow during cardiopulmonary resuscitation. Circulation 61: 345–352.
  39. 39. Tacq J (2004) Part Correlation. The SAGE Encyclopedia of Social Science Research Methods: Sage Publications, Inc.
  40. 40. Guzman PA, Maughan WL, Yin FC, Eaton LW, Brinker JA, et al. (1981) Transseptal pressure gradient with leftward septal displacement during the Mueller manoeuvre in man. British Heart Journal 46: 657–662.
  41. 41. Rahn H, Otis AB, Chadwick LE, Fenn WO (1946) The pressure-volume diagram of the thorax and lung. American Journal of Physiology 146: 161–178.
  42. 42. Ferrigno M, Hickey DD, Liner MH, Lundgren CE (1986) Cardiac performance in humans during breath holding. Journal of Applied Physiology 60: 1871–1877.
  43. 43. Breskovic T, Steinback CD, Salmanpour A, Shoemaker JK, Dujic Z (2011) Recruitment pattern of sympathetic neurons during breath-holding at different lung volumes in apnea divers and controls. Autonomic Neuroscience 164: 74–81.
  44. 44. Pan AW, He J, Kinouchi Y, Yamaguchi H, Miyamoto H (1997) Blood flow in the carotid artery during breath-holding in relation to diving bradycardia. European Journal of Applied Physiology and Occupational Physiology 75: 388–395.
  45. 45. Kety SS, Schmidt CF (1948) The effects of altered arterial tensions of carbon dioxide and oxygen on cerebral blood flow and cerebral oxygen consumption of normal young men. Journal of Clinical Investigation 27: 484–492.
  46. 46. Ainslie PN, Ashmead JC, Ide K, Morgan BJ, Poulin MJ (2005) Differential responses to CO2 and sympathetic stimulation in the cerebral and femoral circulations in humans. The Journal of Physiology 566: 613–624.
  47. 47. von Leupoldt A, Dahme B (2005) Cortical Substrates for the Perception of Dyspnea. Chest 128: 345–354.
  48. 48. Moosavi SH, Golestanian E, Binks AP, Lansing RW, Brown R, et al. (2003) Hypoxic and hypercapnic drives to breathe generate equivalent levels of air hunger in humans. Journal of Applied Physiology 94: 141–154.
  49. 49. Dahl A, Lindegaard KF, Russell D, Nyberg-Hansen R, Rootwelt K, et al. (1992) A comparison of transcranial Doppler and cerebral blood flow studies to assess cerebral vasoreactivity. Stroke 23: 15–19.
  50. 50. Kontos HA (1989) Validity of cerebral arterial blood flow calculations from velocity measurements. Stroke 20: 1–3.
  51. 51. Poulin MJ, Robbins PA (1996) Indexes of Flow and Cross-sectional Area of the Middle Cerebral Artery Using Doppler Ultrasound During Hypoxia and Hypercapnia in Humans. Stroke 27: 2244–2250.
  52. 52. Serrador JM, Picot PA, Rutt BK, Shoemaker JK, Bondar RL (2000) MRI Measures of Middle Cerebral Artery Diameter in Conscious Humans During Simulated Orthostasis. Stroke 31: 1672–1678.
  53. 53. Harms MP, Wesseling KH, Pott F, Jenstrup M, Van Goudoever J, et al. (1999) Continuous stroke volume monitoring by modelling flow from non-invasive measurement of arterial pressure in humans under orthostatic stress. Clinical Science (London) 97: 291–301.
  54. 54. Wong KKL, Mazumdar J, Pincombe B, Worthley SG, Sanders P, et al. (2006) Theoretical modeling of micro-scale biological phenomena in human coronary arteries. Medical & biological engineering & computing 44: 971–982.
  55. 55. Wong KKL, Tu J, Mazumdar J, Abbott D (2010) Modelling of blood flow resistance for an atherosclerotic artery with multiple stenoses and poststenotic dilatations. ANZIAM Journal 51.
  56. 56. Wong KKL, Tu J, Sun Z, Dissanayake D (2013) Methods in Research and Development of Biomedical Devices. Singapore: World Scientific Publishing Co.