Oscillations of heart rate and respiration synchronize during poetry recitation

Dirk Cysarz, Dietrich von Bonin, Helmut Lackner, Peter Heusser, Maximilian Moser, Henrik Bettermann


The objective of this study was to investigate the synchronization between low-frequency breathing patterns and respiratory sinus arrhythmia (RSA) of heart rate during guided recitation of poetry, i.e., recitation of hexameter verse from ancient Greek literature performed in a therapeutic setting. Twenty healthy volunteers performed three different types of exercises with respect to a cross-sectional comparison: 1) recitation of hexameter verse, 2) controlled breathing, and 3) spontaneous breathing. Each exercise was divided into three successive measurements: a 15-min baseline measurement (S1), 20 min of exercise, and a 15-min effect measurement (S2). Breathing patterns and RSA were derived from respiratory traces and electrocardiograms, respectively, which were recorded simultaneously using an ambulatory device. The synchronization was then quantified by the index γ, which has been adopted from the analysis of weakly coupled chaotic oscillators. During recitation of hexameter verse, γ was high, indicating prominent cardiorespiratory synchronization. The controlled breathing exercise showed cardiorespiratory synchronization to a lesser extent and all resting periods (S1 and S2) had even fewer cardiorespiratory synchronization. During spontaneous breathing, cardiorespiratory synchronization was minimal and hardly observable. The results were largely determined by the extent of a low-frequency component in the breathing oscillations that emerged from the design of hexameter recitation. In conclusion, recitation of hexameter verse exerts a strong influence on RSA by a prominent low-frequency component in the breathing pattern, generating a strong cardiorespiratory synchronization.

  • creative arts therapy
  • cross-sectional study design
  • bivariate data analysis
  • heart rate variability

the effects of different breathing frequencies and patterns on cardiovascular regulation have been investigated extensively in recent years. In this context, various effects of poetry recitation on cardiovascular parameters, especially on heart rate oscillations, have been demonstrated (4, 9, 50). Bernardi et al. (4) found a frequency adjustment of breathing oscillations with endogenous blood pressure fluctuations (Mayer waves) and even cerebral blood flow oscillations during the recitation of the rosary and the “OM” mantra. This effect was attributed to the breathing frequency of ∼6 breaths/min induced by the metric of both religious verses. Furthermore, they noticed an increased arterial baroreflex sensitivity, which is a favorable long-term prognostic factor in cardiac patients (3). In another study, Bernardi et al. (5) observed a significant increase in arterial oxygen saturation (SaO2) during controlled breathing at frequencies of 15/6/3 breaths/min in patients with chronic heart failure and healthy controls. The strongest increase was found at a breathing frequency of 6 breaths/min. Thus recitation of specific poetry as a means to control breathing patterns was proposed and the rosary prayer to be “viewed as a health practice as well as a religious practice” (4).

In our own investigations on cardiovascular and cardiorespiratory regulation during and after recitation of poetry, we were led by the observation of therapists that creative arts greatly influence well-being in humans by various means. On a psychosomatic level, these therapies increase the salutogenetic potential in humans (2, 17). Furthermore, the autonomic regulation is affected to enhance the flexibility of regulatory processes to maintain stability and coherence between different functions. As a result, a temporal order of physiological functions appears (20). This seems to apply especially when therapeutic speech as rhythmic poetry is recited using different breathing modalities. These breathing modalities are used to activate or calm down the patient. For example, some exercises of Anthroposophic Therapeutic Speech (ATS) (14), which is originally based on the philosophy of R. Steiner, use the recitation of poetry to treat stress-related symptoms in the cardiorespiratory system (50).

The physiological influences of various breathing patterns on heart rate fluctuations are well known (1, 6, 15, 21). Thus, in a first study, we (50) examined heart rate variability (HRV) in healthy subjects during and after recitation of poetry. Two different examples of old poetry, hexameter and alliterative verse, were used. Although no pacemakers had been used, recitation of hexameter verse always modulated heart rate at a frequency of 12 and 6 beats/min (i.e., 2:1 frequency ratio). Furthermore, during 15 min of rest after the exercise, an increased high-frequency component of HRV (50) and a predominance of typical short “heart rate patterns” resulting from intermittent cardiorespiratory synchronization were observed (9). In contrast, these observations were not found during/after normal conversation, which was used as a control exercise. Thus recitation of poetry changed cardiorespiratory interaction, whereas normal conversation did not.

Many features of the cardiorespiratory control during recitation of poetry are still unknown. Recently, on the basis of simultaneous recordings of an electrocardiogram (ECG) and a respiratory trace, new techniques for the analysis of cardiorespiratory interaction have been developed (41, 42). They unambiguously revealed that heart rate and respiration may intermittently synchronize. The application of these techniques promises new information about the cardiorespiratory interaction, e.g., after myocardial infarction (22, 25).

In this study, we investigated the cardiorespiratory synchronization in healthy subjects during recitation of hexameter verse because our first results were most obvious with this exercise. Cardiorespiratory synchronization was analyzed with respect to a respiratory trace and the oscillations in heart rate induced by respiration, i.e., respiratory sinus arrhythmia (RSA). For this purpose, we adapted recently developed techniques because they offer to define a continuous phase that may easily be analyzed (23, 43). We show that this technique is able to capture important information. Three different exercises were compared using a cross-sectional study design: 1) recitation of hexameter verse, 2) controlled breathing, and 3) spontaneous breathing. The results of this study may improve our understanding of regulatory processes that maintain stability and coherence between different physiological functions because cardiorespiratory interaction seems to play a crucial role in this context.



Twenty-three healthy subjects without experience in ATS and without prior knowledge of the hexameter text used for the recitation were enrolled in the study. After an initial check, three subjects had to be excluded due to frequent ectopic heartbeats. The remaining 20 subjects (10 females; age: 43 ± 6.6, average ± SD; 3 smokers) were included and had no history of cardiovascular diseases, especially no hypo- or hypertension or antiarrhythmical therapy. All subjects gave their informed written consent to take part in the study. The study protocol was approved by the Ethics Commitee of the University of Berne, Switzerland.

Experimental Procedures

All subjects were invited individually three times to the therapy center at the same time of day. In each of the three sessions, the subjects performed a different exercise (in random order): hexameter recitation (H), controlled breathing (C), and spontaneous breathing (S), see Table 1. During each session, an ECG and the nasal/oral airflow were recorded simultaneously (see Data Acquisition and Preprocessing). The overall duration of each session was 50–60 min, divided into three successive measurements: 15-min quiet rest in a resting chair (baseline measurement S1; see Table 1), 20-min exercise measurement (H, C, S), and 15-min quiet rest in a resting chair (effect measurement S2). During S1 and S2, the subjects were allowed to breathe spontaneously. This procedure resulted in nine different measurements of each subject. To ensure comparable levels of physical activity during the three types of exercises, the subjects walked through the room at a pace of 50 steps/min (given by an electric metronome). The three experiments had to be at least 24 h apart but within 14 days.

View this table:
Table 1.

Experimental protocol

Hexameter recitation.

Hexameter is the oldest rhythmic verse in ancient Greece, where it is found in the two of the largest pieces of epic poetry yet known, the Iliad and the Odyssey. A line of hexameter verse contains six feet (a foot is the basic rhythmic unit), usually dactyls (a single dactyl comprises one long and two short parts). The line is mostly split in two parts by a break (caesura), which in the most archetypal form (Iliad) mainly takes place after the long part of the third dactyl. This results in a convenient and regular breathing pattern for the speaker (24). We used a piece from Homers Odyssey in a German translation that did not alter the rhythmic scheme of the verse1.

Hexameter recitation was carried out with the mentioned breathing pattern. The therapist walked at the same pace as the subject, speaking the first part of the hexameter line. The sequence of steps was divided as follows: inspiration always took one step, speaking three steps. To support a full inspiration, recitation was accompanied by lifting the stretched arms to the level of the shoulders during inspiration and lowering the arms during recitation. The subject listened to the text recited by the therapist without lifting the arms (but continued walking) and subsequently repeated it in the therapist's fashion. This alternating procedure was repeated for 20 min. The walking pace, the movement of the arms, and the alternating fashion of listening and reciting ensured a high temporal coordination between these tasks and made the recitation more comfortable.

Controlled breathing.

Controlled breathing consisted of breathing in the “hexameter pattern,” i.e., inspiration during one step, expiration during three steps. The duration of expiration (3 steps) was maintained by control of breath with the lips. Similar to hexameter recitation, the therapist first demonstrated the specific type of breathing. Again, the breathing was accompanied by lifting the stretched arms to the level of the shoulders during inspiration and lowering the arms during expiration. The subject “listened” to the breathing by the therapist without lifting the arms and subsequently imitated it. This alternating procedure ensured a comparable breathing pattern with respect to hexameter recitation.

Spontaneous breathing.

During spontaneous breathing, the subjects walked at the same pace as in the two other exercises, but they were instructed to breathe spontaneously. Furthermore, there were no restrictions with respect to the movement of the arms.

Data Acquisition and Preprocessing

The ECG (standard lead) and the uncalibrated nasal/oral airflow (derived by 3 thermistors that were placed next to the nostrils and in front of the mouth) were recorded simultaneously in all subjects using solid-state recorders (Medikorder MK2, Tom-Signal, Graz, Austria). The sampling rates of the ECG and the nasal/oral airflow were 3,000 and 50 Hz, respectively. This ensured an accuracy <1 ms for the times of the identified R peaks. The inspiratory and expiratory onsets were defined as the local minima and the local maxima, respectively, of the nasal/oral airflow because they were due to the change from exhaled warm air (warmed by the respiratory tract) to inhaled air at the temperature of the environment (and vice versa). For further analysis, the acquired data were saved to a file and were further processed using Matlab (The Mathworks, Natick, MA) and C routines. All automatically identified R peaks were visually controlled and edited where necessary. The manually edited R peaks (∼0.1% of all R peaks) had an accuracy of 5 ms because the ECG was recorded with a sampling rate of 200 Hz.

To obtain a heart rate time series with equidistant time steps, the times of successive R peaks were first converted to a R-R tachogram, i.e., the sequence of times between successive R peaks. Next, the resulting R-R tachogram was resampled at a rate of 5 Hz using linearly interpolated values. To get a time series for the nasal/oral airflow at corresponding sampling times, each 10th sample was used. These two time series share a common time axis and served as the basis for further calculations.

Figure 1A shows a short trace of the ECG and the simultaneous nasal/oral airflow during recitation of hexameter verse. The breathing oscillations show two different frequencies. A high-frequency component at ∼12 breaths/min resulted from the design of the hexameter exercise (duration of a full respiratory cycle: ∼4.8 s, see Experimental Procedures). Furthermore, a low-frequency component appears at ∼6 cycles/min. This low-frequency oscillation reflects the alternation of a shallow breathing cycle during listening and a deeper breathing cycle during recitation. The heart rate time series shows an oscillation that is obviously synchronized to the low-frequency oscillations of the breathing pattern. Thus, to obtain a nasal/oral time series that contains merely the low-frequency component of the breathing pattern, the time series was band-pass filtered. Because the filtering must leave the phase of the time series unaltered, the technique used in this study consisted of two consecutive steps. First, a moving average with two windows of different lengths (20 and 60 samples, i.e., 4 and 12 s) was calculated. Next, the filtered time series was obtained by subtraction of the moving average time series with the longer window length from the moving average time series with the shorter window length, see Fig. 1B. The changes in the amplitude of the oscillations caused by this filtering technique are not disadvantageous in the context of this study because the amplitudes are irrelevant for the analysis of synchronization. If higher moments of the time series have to be preserved, more sophisticated filtering techniques should be used, e.g., the Savitzky-Golay filter (34).

Fig. 1.

A: heart rate time series derived from the interpolated R-R tachogram (RRintp.; top) and nasal/oral airflow [in arbitrary units (a.u.); bottom] during recitation of hexameter before filtering. Obviously, the heart rate time series mainly shows a low-frequency oscillation, whereas the nasal/oral airflow contains a low- and a high-frequency oscillation. B: heart rate time series and nasal/oral airflow after the filtering procedure. In the nasal/oral airflow, only the low-frequency oscillation remains. The synchronization is obvious.

Analysis of Cardiorespiratory Phase Synchronization

The median heart rate and the median breathing frequency of each measurement served as the basic parameters. Because the breathing patterns of the three exercises modulate the heart rate differently, the extent of RSA of each measurement is expressed as the median of the longest R-R interval minus the shortest R-R interval of each respiratory cycle. In the following, two different techniques are used to quantify cardiorespiratory synchronization. The first one is based on the Fourier transformation of both filtered time series xheart and xresp, and their corresponding cross spectrum Chr(f) = (F xheart)(f)(F xresp)*(f) (where F denotes the Fourier transformation and * denotes complex conjugate). The normalized amplitude of Chr(f) is called the coherence function and may be used as a simple measure of synchronization (35) Math

The brackets 〈. . .〉 denote an average using epochs containing 300 samples, i.e., 1 min, and an overlap of 150 samples between successive epochs. In this study, ξ = Γhr(fresp) at the centre frequency fresp of the nasal/oral airflow serves as a simple means to quantify cardiorespiratory synchronization.

The second technique has its origin in the analysis of two weakly coupled chaotic systems (38, 39). First of all, because each oscillator is completely described by its amplitude and its phase as a function of time, a suitable phase has to be defined for a time series that merely contains oscillations in a relatively narrow frequency band. For the purposes of this study, the amplitude is not required and, hence, not calculated. In recent studies, the phase of such a time series was defined using the concept of analytic signals (37, 40, 43). They made use of the Hilbert-transform (ti) of the band-pass filtered time series x(ti). The Hilbert-transform is defined as the convolution of x(ti) and the function 1/(πti) (11) Math

Notice that x denotes a convolution. Practically, the convolution may be calculated more easily in the frequency domain using the Fourier-Transform of this equation because the convolution is then turned into a simple multiplication. With the definition of the Hilbert-transform, the phase φ(ti) of a time series is obtained by calculation of Math

With the use of this definition, the phase is calculated for the interpolated heart rate times series and the filtered nasal/oral airflow, yielding a phase φheart(ti) and a phase φresp(ti). Next, the difference φ(ti) between these two phases is Math

The heart rate time series and the nasal/oral time series, i.e., RSA and breathing pattern, are synchronized if this phase difference is constant, i.e., Math with δ being a constant offset [the phase difference φ(ti) needs not necessarily to be around zero]. In the case of synchronization, the constant phase difference φ(ti) is proportional to the time delay between both time series because neither the definition of the phases φheart(ti) and φresp(ti) nor the preceding filtering procedure introduced any additional constant or time-dependent phase shift.

Unfortunately, noise and other sources of interference in both time series lead to random-like “phase jumps” of ±2π in the sequence of φheart(ti) and φresp(ti). Thus, even in a synchronized state, the phase difference φ(ti) is not constant anymore. Hence, the above mentioned condition for a synchronized state is not suitable anymore. Instead, the synchronized state is characterized by a statistical preference of some values of Math i.e., a preference of some values of φ(ti) ± 2 (with n being an integer). Thus, in the case of synchronization of two systems that are disturbed by noise, the distribution of Ψ(ti) shows one unambiguous maximum. In this study, the distribution of Ψ(ti) is quantified by Math where brackets 〈. . .〉 denote an average (27, 46). Theoretically, if γ = 1, Ψ(ti) is constant because both time series are completely synchronized in a statistical sense. In this case, the location of the maximum Ψ(ti) is the preferred phase difference between both time series. If γ = 0, both time series are completely desynchronized because the values of Ψ(ti) are equally distributed in the range [−π, π], i.e., no preference of any phase at all. For further details, see the literature.

For real-world data, the lower bounds of γ and ξ have to be estimated because even in the absence of any coupling synchronized patterns may appear by chance. To accomplish this task, the concept of surrogate data for bivariate data is used (28). The bivariate surrogate data were created as follows. The nasal/oral airflow was left unchanged, whereas the sequence of the original R-R tachogram was randomized. Subsequently, the heart rate time series was constructed as described above. This procedure maintains the distribution of the R-R tachogram, i.e., the mean and standard deviation of the R-R distances are the same as in the original R-R tachogram, whereas the temporal structure is completely destroyed. Hence, any cardiorespiratory synchronization due to coupling is also destroyed. In the surrogate data, spurious synchronization may occur due to fluctuations of the R-R distances and the subsequent filtering procedure used for the synchronization analysis. Without going into detail, the 95% percentile of the surrogate data of all 180 measurements (one realization per measurement) yields γ = 0.14 (analysis of phase difference) and ξ = 0.23 (coherence analysis). These estimations serve as a lower bounds of γ and ξ. It shows that even in the absence of synchronization, the fluctuations in the randomized R-R tachogram result in low-frequency oscillations that may be spuriously synchronized to the nasal/oral airflow. In the following, γ and ξ serve as two different indexes of cardiorespiratory synchronization.

An example of cardiorespiratory phase synchronization analysis is illustrated in Fig. 2. The phases φheart(ti) and φresp(ti) have been calculated as described above. First, during the baseline measurement, cardiorespiratory interaction is not synchronized. As a result, the phase difference φ(ti) increases with time, see Fig. 2A. Thus the values of Ψ(ti) are almost equally distributed resulting in a low γ value (γ = 0.14) that is equal to the estimated lower bound. Hence, this example contains hardly any cardiorespiratory synchronization. Recitation of hexameter verse synchronizes both systems, which is reflected in both time series, see for example Fig. 1B. Thus the phase difference φ(ti) is approximately constant for some period of time and in the phase difference certain plateaus appear (Fig. 2B). Furthermore, the distribution of Ψ(ti) shows one distinct peak that leads to a much higher γ value (γ = 0.78). In this example, the peak of Ψ(ti) is located at Ψ(ti) = 1.58, which is approximately π/2. Hence, the phase of the respiratory time series advances the phase of the heart rate time series by approximately a quarter of an oscillation. Notice that the analysis of the coherence as a measure of cardiorespiratory synchronization yields ξ = 0.59 and ξ = 0.68, respectively, for the examples in Fig. 2, A and B.

Fig. 2.

A: increase of the phase difference φresp(ti) − φheart(ti) during baseline measurement S1 denotes a desynchronized state between the heart rate time series and the low-frequency component of the nasal/oral time series (γ = 0.14). B: cardiorespiratory synchronization during hexameter recitation: the phase difference shows plateaus at different levels and an obvious maximum in the distribution of Ψ(ti), which is also reflected by γ = 0.78.


The objective of this study was to assess the effects of hexameter recitation on cardiorespiratory synchronization in contrast to controlled breathing and spontaneous breathing. To this end, descriptive methods are used. Because the number of subjects is small (n = 20) and the distribution of the γ values is not known, the median instead of the mean is used to quantify the distributions. Median heart rate, median frequency of the low-frequency breathing oscillations, extent of RSA, the γ value, and the ξ value were calculated for each measurement of each subject. Effects of transitions at the beginning of each measurement were reduced by omitting the first 2 min of the recording. To avoid a bias due to slightly different durations of each measurement, the subsequent 13 min were analyzed. Seven of 180 measurements (20 subjects × 3 exercises × 3 measurements per exercise) were slightly shorter than 13 min. Subsequently, the three measurements of each exercise (S1, exercise, S2) were characterized by the median, lower, and upper quartile of the subjects' median heart rate, median frequency of low-frequency breathing oscillations, and extent of RSA. Box and whisker plots were used for visualization of the distribution of the γ and ξ values.

The nonparametric Friedman test was used to calculate the probability of equality between the three measurements of each session (S1, exercise, and S2) and between the three exercise measurements (hexameter, controlled breathing, and spontaneous breathing). Subsequently, an appropriate post hoc test for multiple comparisons was used to calculate the probability of equality between two different measurements (10). A PFriedman value near zero indicates a high probability of differences between the three measurements with respect to the analyzed parameter. For low values of PFriedman, low P values of the post hoc test indicate the pair of measurements with likely differences.


In total, 180 measurements were analyzed with respect to cardiorespiratory synchronization. All exercises were associated with an increase of heart rate compared with baseline (S1) and effect measurement (S2) (Table 2). The hexameter exercise showed the highest heart rate (82.9 beats/min). Furthermore, in all three exercises, the heart rate of S2 was decreased compared with S1. The frequency of the low-frequency breathing oscillations decreased during hexameter exercise and controlled breathing exercise compared with S1 and S2 (Table 3). During the spontaneous breathing exercise, this frequency increased compared with S1 and S2. The frequency of the low-frequency breathing oscillations was lowest for the hexameter exercise (6.4 breaths/min) and highest for the spontaneous breathing exercise (12.5 breaths/min). The measurements of S1 and S2 showed an intermediate frequency of the low-frequency component at ∼8.5 breaths/min. Although the frequency of the low-frequency breathing oscillations decreased during the hexameter exercise and during controlled breathing, the RSA did not change noticeably compared with S1 and S2 (Table 4). On the contrary, the increase of this frequency during spontaneous breathing decreased the extent of RSA compared with S1 and S2.

View this table:
Table 2.

Median, lower, and upper quartile of the subjects' median heart rate of each measurement

View this table:
Table 3.

Median, lower, and upper quartile of the subjects' median frequency of the low-frequency breathing oscillations of each measurement

View this table:
Table 4.

Total median, lower, and upper quartile of the subjects' extent of RSA of each measurement

The results of the γ values as a quantitative index of cardiorespiratory phase synchronization are shown in Fig. 3, A–C. Clearly, the recitation of hexameter and controlled breathing led to an increase of the γ values compared with S1 and S2. The increase was largest for the recitation of hexameter verse (hexameter exercise: median γ = 0.70, controlled breathing exercise: median γ = 0.57). For both exercises, the comparisons S1 vs. exercise and exercise vs. S2 yielded low P values. Thus, during these two exercises, the cardiorespiratory interaction was more synchronized compared with S1 and S2. In contrast to these findings, spontaneous breathing led to a decrease of the γ values compared with S1 and S2 (median γ = 0.15). This is very close to the lower bound of the γ values and indicates a high degree of desynchronization. Again, the comparison S1 vs. exercise and exercise vs. S2 yielded low P values. Hence, during the spontaneous breathing exercise, the cardiorespiratory interaction was less synchronized compared with S1 and S2. Comparing S1 and S2 of all exercises, the γ values always were γ ≈ 0.3.

Fig. 3.

Box and whisker plots of the γ values quantifying cardiorespiratory phase synchronization. A–C: comparison of the 3 measurements of each exercise. D: comparison of the 3 exercise measurements [hexameter (H), controlled breathing (C), spontaneous breathing (S)]. In all diagrams, a low value of PFriedman indicates likely differences between the γ values of the 3 measurements. The probability of similar γ values between 2 measurements is indicated by the P values above the box and whisker plot. The box plots show median and quartiles (horizontal lines), mean value (*), and maximum and minimum values (whiskers).

In Fig. 3D, the comparison of the γ values during the three different exercises is shown. This diagram shows that the difference of the γ values between the hexameter exercise and the controlled breathing exercise was large enough to yield a low P value. Thus reciting hexameter verse results more often in a synchronized cardiorespiratory interaction than does controlled breathing. Almost trivially, the low γ values of the spontaneous breathing exercise produced low P values between the other two exercises. As these γ values are also lower than before or after the exercise, the spontaneous breathing exercise desynchronizes the oscillations of heart rate and respiration.

The results of the distributions of ξ values derived from the coherence analysis are shown in Fig. 4, A–C. Although the hexameter recitation leads to an increase of the ξ values compared with S1 and S2, only the difference between hexameter and S2 is likely as indicated by the low P value. In contrast to the analysis of phase differences, the controlled breathing exercise does not show any difference between exercise, S1, and S2. And during the spontaneous breathing exercise, the ξ values are slightly lower compared with S1, but compared with S2, no difference is observable. The comparison of the three exercises in Fig. 4D shows a similar result compared with the γ values. Notice that all PFriedman values are increased compared with the analysis of phase synchronization, indicating that the differences between the measurements and exercises are smaller.

Fig. 4.

Box and whisker plots of the ξ values quantifying cardiorespiratory coherence as a simple measure of cardiorespiratory synchronization (for further details, see Fig. 3).

The synchronization of oscillations in heart rate and respiration during hexameter recitation and for controlled breathing allowed the calculation of the maximum of the distribution of Ψ(ti) as the preferred phase difference between both time series. The preferred phase difference differed largely in each exercise and thus could not easily be condensed in one number (see Table 5). Furthermore, most subjects showed a different preferred phase difference for recitation of hexameter verse than for controlled breathing. Remarkably, during both exercises, many subjects showed a preferred phase difference of ∼3/4π ≈ 2.4. In this case, the low-frequency breathing oscillations preceded the heart rate oscillations by this phase difference, i.e., less than half an oscillation.

View this table:
Table 5.

Preferred phase difference Ψ(ti) between heart rate and breathing oscillations during hexameter and controlled breathing exercise of each subject


In this study, cardiorespiratory interaction during recitation of hexameter verse has been investigated. In ATS, this exercise is often used to “harmonize” and “strengthen” rhythmic processes in humans. This is thought to be achieved by the recitation of speech that is supported by the therapeutic method the recitation is carried out: the breathing pattern is coordinated to the walking pattern and to the rhythm of the poetry. With respect to cardiorespiratory interaction, the results of the analysis of the phase difference and the coherence analysis revealed the following classification. First, during recitation of hexameter verse, the low-frequency oscillations of the breathing pattern were synchronized to a large extent with the heart rate oscillations. Second, the cardiorespiratory interaction was also synchronized during the controlled breathing exercise, but to a slightly lesser extent. Third, the resting periods before and after the exercises showed a further reduction of cardiorespiratory synchronization. Finally, during the spontaneous breathing exercise, the cardiorespiratory interaction was almost completely desynchronized. Rhythmic speech thus has the strongest impact on synchronization of low-frequency breathing oscillations and heart rate fluctuations, whereas cardiorespiratory interaction during everyday activities is rarely synchronized.

This kind of cardiorespiratory interaction has its origin in the modulation of heart rate by respiration, i.e., RSA (1). For breathing frequencies above 0.15 Hz, the magnitude of RSA serves as an index of parasympathetic activity of the autonomic nervous system (15) because the latencies of heart rate responses and the corresponding decay constants of the vagal cardiac effectors are far shorter compared with the sympathetic branch (7). However, this notion has been challenged recently by the finding that sympathetic activity also influences RSA in the frequency range of breathing oscillations (47). RSA is a frequency- and amplitude-dependent phenomenon. If the breathing frequency is constant, the magnitude of RSA increases as the tidal volume increases (21). For a given tidal volume, many studies consistently showed that breathing oscillations modulate heart rate strongest for frequencies below 0.15 Hz (6, 12, 18, 21, 32, 44). Mental stress decreases RSA, whereas reading aloud shifts RSA fluctuations into the low-frequency range of the breathing frequency and thus increases RSA (6). Furthermore, the magnitude of RSA is larger if a short, rapid inspiration is followed by a long expiration than vice versa (45). Remarkably, during speech production, the physiological demands regarding pulmonary gas exchange may be overridden within wide limits (13, 30). This effect is at least partly counterbalanced by an improved pulmonary gas exchange resulting from RSA (16, 19). This short and by no means complete description may give an impression of the complexity of interactions that affect the magnitude of RSA. Although the phase relationship between respiratory oscillations and modulations of heart rate has been investigated in experimental situations (1, 15), to the best of our knowledge cardiorespiratory synchronization between breathing oscillations and RSA during speech production has not been systematically investigated.

During hexameter recitation, most of the mentioned influences have to be taken into account. A slight mental activity was required to recite the text properly. The adjusted walking scheme during the exercise needed mild physical effort. Both mental and physical activity increased the heart rate. On the other hand, the recitation produced a low-frequency component in the breathing oscillations at ∼6 cycles/min, i.e., one-half of the actual breathing frequency. This oscillation appeared due to a longer duration and a larger amplitude of the breathing cycle during recitation compared with listening (Fig. 1). In contrast to the mental and physical activities, these properties lead to an increase of the magnitude of RSA. In summary, the reduction of the RSA due to an increase of the heart rate during recitation and controlled breathing was compensated by the increase of RSA due to the alternation of recitation and listening. Furthermore, the results of the analysis of cardiorespiratory synchronization revealed that during hexameter recitation RSA is synchronized with the low-frequency component of the respiratory oscillations. The phase difference between the phase of the low-frequency oscillations of the respiratory signal and the phase of the heart rate time series is 3/4π ≈ 2.4 for many subjects, which is in agreement with the results of an experimental study at low breathing frequencies (15).

The emergence of cardiorespiratory synchronization during a regular breathing pattern at a low (and almost constant) frequency may be explained by the regular excitatory and inhibitory effects of the central respiratory generators on vagal and sympathetic outflow. Surprisingly, the high-frequency breathing oscillations during hexameter recitation (unfiltered heart rate time series; Fig. 1A) modulated the heart rate only to a minor extent. Thus the low-frequency modulations of vagal and sympathetic cardiac effectors seem to override the modulations at higher frequencies. This effect may be due to the local maximum of the cardiorespiratory transfer function at low frequencies (∼0.1 Hz) (7). However, the complete physiological origin of the phase synchronization is difficult to explain, because the contributions of the central and peripheral mechanisms to the generation of RSA are not yet fully understood (7, 31). Other mechanisms, such as the optimization of pulmonary gas exchange and the increase in arterial oxygen saturation (5, 19), may also be of importance for the emergence of cardiorespiratory synchronization.

In principle, the same holds for the controlled breathing exercise except that this exercise showed synchronization to a lesser degree. A reason could be a slightly lower extent of RSA compared with the hexameter recitation. Although we did not control for it, the differences in the extent of synchronization and RSA might be due to a slightly larger tidal volume during recitation compared with controlled breathing. During spontaneous breathing, the breathing oscillations did not contain a low-frequency component. Hence, during this exercise, the reduction of the magnitude of RSA by the physical activity was not compensated and, as a consequence, cardiorespiratory synchronization hardly ever occurred. Although there were no restrictions on breathing behavior during the resting periods, cardiorespiratory synchronization then occurred more often than during spontaneous breathing but less often than during the other two exercises. This result may be explained by the reduced physical activity that decreased the heart rate and, in turn, increased the magnitude of RSA compared with the spontaneous breathing exercise. In principle, these results are in accordance with the recent observation that the extent of cardiorespiratory synchronization increases if the breathing frequency decreases (33). However, in the present study, the degree of synchronization depended on different breathing patterns and on physical activity, whereas in the cited study only the breathing frequency varied. A study where both, breathing frequency and physical activity, were varied with respect to cardiorespiratory synchronization has not been carried out yet.

Unlike the analysis of special heart rate patterns during hexameter recitation in our previous investigation that incorporated a combination of a cross-sectional and a longitudinal study design with a low number of subjects (9), we did not find any immediate effect after the exercises. In other words, hexameter recitation did not affect cardiorespiratory synchronization straight after the exercise. However, the analysis of the previous data only showed immediate effects after the subjects were familiar with the exercise, i.e., after 2–3 wk of exercise. Additionally, our previous investigation dealt with an analysis of specific “heart rate patterns” (8) that is based on cardiorespiratory interaction at breathing frequencies well above 0.1 Hz. In contrast, the present study is based on a cross-sectional study design analyzing cardiorespiratory interaction at breathing patterns of ∼0.1 Hz. Thus persistent effects of hexameter recitation on HRV and presumably also on cardiorespiratory interaction need a certain familiarization with the exercise.

Compared with established techniques to analyze bivariate data, such as the calculation of the coherence, the analysis of phase differences yields more consistent results with respect to cardiorespiratory interaction, i.e., to the interplay of respiratory oscillations with heart rate fluctuations. The calculation of the coherence separated the different exercises to a lesser degree, particularly the differences between S1, exercise, and S2 were less obvious. Thus especially during uncontrolled breathing, the coherence analysis seems to be less powerful. A reason is that the coherence mainly puts emphasis on the appearance of certain frequencies in both time series regardless of the phase difference between the oscillations. On the contrary, the analysis of phase differences is explicitly based on the interaction between both oscillations. A more detailed comparison and an analysis of the features and restrictions of each method are beyond the scope of this study and will be deferred to another study. Using the method presented here, we also expect cardiorespiratory synchronization for other techniques of artistic or religious speech creation, e.g., the “OM” meditation, because they also use breathing patterns at low frequencies (26, 29, 48, 49). Furthermore, data from many other studies investigating influences of low-frequency breathing patterns on HRV likely contain cardiorespiratory synchronization. More generally, because speech production is accompanied by an increase of the magnitude of RSA (36), we do not preclude the possibility that cardiorespiratory synchronization may also emerge to a certain extent during normal speech production.

In conclusion, the special breathing pattern used for the recitation of hexameter verse produced a strong cardiorespiratory synchronization with respect to low-frequency breathing oscillations and heart rate variations. Controlled breathing showed cardiorespiratory synchronization to a lesser extent. The physiological origin of this kind of cardiorespiratory interaction still needs to be explored in more detail. Furthermore, it remains to be shown if and how cardiorespiratory synchronization affects arterial oxygen saturation at low breathing frequencies (5). A possible explanation would be that cardiorespiratory synchronization is more advantageous with respect to oxygen uptake because it takes advantage of synergistic effects.


This work was financially supported by Weleda, Germany (to D. Bonin, D. Cysarz, and H. Bettermann) and Wala, Germany (to D. Bonin).


D. Auerbach helped in the preparation of the final version of the manuscript.


  • 1 The enormous length of this narrative poetry (the Iliad consists of 15,693 lines of hexameter), which has been recited by special “singers” (aoides) in full length (24), seems to suggest at least a harmonious, beneficial, health-maintaining influence not only on the speaker but also on the audience. However, there is a certain disagreement regarding the purpose of the caesura. Is it mainly introduced to take a breath or to put emphasis on certain words or parts of a phrase? We tend to prefer the former, although both interpretations might apply. In ATS, the rhythmic function of the caesura is of importance. The caesura is regarded to take the duration of another dactyl. Thus taking into account that only half a line is recited each time, each breath cycle in rhythmic hexameter recitation covers the length of 4 feet (recitation: 3 feet, caesura or breath taking: 1 foot). This pattern is experienced as especially harmonious.

  • 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.


View Abstract