# A new analysis of heart rate variability in the assessment of fetal parasympathetic activity: an experimental study in a fetal sheep model

Charles Garabedian, C. Champion, Elie Servan-Schreiber, Laura Butruille, Estelle Aubry, Dyuti Sharma, Regis Logier, Philippe Deruelle, Laurent Storme, Veronique Debarge, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Charles Garabedian, C. Champion, Elie Servan-Schreiber, Laura Butruille, Estelle Aubry, et al. A new analysis of heart rate variability in the assessment of fetal parasympathetic activity: an experimental study in a fetal sheep model. PLoS ONE, 2017, 12 (7), pp.e0180653. 10.1371/journal.pone.0180653. hal-02519669

HAL Id: hal-02519669

https://hal.univ-lille.fr/hal-02519669

Submitted on 26 Mar 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# A new analysis of heart rate variability in the assessment of fetal parasympathetic activity: An experimental study in a fetal sheep model

C. Garabedian¹,²*, C. Champion¹,³, E. Servan-Schreiber¹,², L. Butruille¹, E. Aubry¹,³, D. Sharma¹,³, R. Logier⁴, P. Deruelle¹,², L. Storme¹,⁵, V. Houfflin-Debarge¹,², J. De Jonckheere¹,⁴

1 Univ. Lille, EA 4489 - Perinatal Environment and Health, Lille, France, 2 CHU Lille, Department of Obstetrics, Lille, France, 3 CHU Lille, Department of Pediatric Surgery, Lille, France, 4 CHU Lille, CIC-IT 1403, MRRC, Lille, France, 5 CHU Lille, Department of Neonatology, Lille, France

* charles.garabedian@chru-lille.fr

<figure><img src="image_3.png" /></figure>

## Abstract

Analysis of heart rate variability (HRV) is a recognized tool in the assessment of autonomic nervous system (ANS) activity. Indeed, both time and spectral analysis techniques enable us to obtain indexes that are related to the way the ANS regulates the heart rate. However, these techniques are limited in terms of the lack of thresholds of the numerical indexes, which is primarily due to high inter-subject variability. We proposed a new fetal HRV analysis method related to the parasympathetic activity of the ANS. The aim of this study was to evaluate the performance of our method compared to commonly used HRV analysis, with regard to i) the ability to detect changes in ANS activity and ii) inter-subject variability. This study was performed in seven sheep fetuses. In order to evaluate the sensitivity and specificity of our index in evaluating parasympathetic activity, we directly administered 2.5 mg intravenous atropine, to inhibit parasympathetic tone, and 5 mg propranolol to block sympathetic activity. Our index, as well as time analysis (root mean square of the successive differences; RMSSD) and spectral analysis (high frequency (HF) and low frequency (LF) spectral components obtained via fast Fourier transform), were measured before and after injection. Inter-subject variability was estimated by the coefficient of variance (%CV). In order to evaluate the ability of HRV parameters to detect fetal parasympathetic decrease, we also estimated the effect size for each HRV parameter before and after injections. As expected, our index, the HF spectral component, and the RMSSD were reduced after the atropine injection. Moreover, our index presented a higher effect size. The %CV was far lower for our index than for RMSSD, HF, and LF. Although LF decreased after propranolol administration, fetal stress index, RMSSD, and HF were not significantly different, confirming the fact that those indexes are specific to the parasympathetic nervous system. In conclusion, our method appeared to be effective in detecting parasympathetic inhibition. Moreover, inter-subject variability was much lower, and effect size higher, with our method compared to other HRV analysis methods.

<figure><img src="image_4.png" /></figure>

### OPEN ACCESS

**Citation:** Garabedian C, Champion C, Servan-Schreiber E, Butruille L, Aubry E, Sharma D, et al. (2017) A new analysis of heart rate variability in the assessment of fetal parasympathetic activity: An experimental study in a fetal sheep model. PLoS ONE 12(7): e0180653. https://doi.org/10.1371/journal.pone.0180653

**Editor:** Elena Tolkacheva, University of Minnesota, UNITED STATES

**Received:** January 25, 2017

**Accepted:** June 19, 2017

**Published:** July 10, 2017

**Copyright:** © 2017 Garabedian et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The authors received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.
## Introduction

Under hypoxia, the fetus develops several coping mechanisms, including an increase in $O_2$ extraction by tissues and the redistribution of blood flow toward the vital organs. These adaptive mechanisms are regulated by the autonomic nervous system (ANS) via sympathetic and parasympathetic activities. The ANS is highly sensitive to hypoxia and acidosis, through various well-described mechanisms, including the stimulation of chemoreceptors and baroreceptors, or the direct actions of hypoxia and hypercapnia on the brainstem [1]. Therefore, an analysis of ANS activity could be an efficient tool in evaluating fetal well-being [2].

Analysis of heart rate variability (HRV) is a recognized non-invasive tool that is used to assess ANS regulation, and a variety of time and spectral analysis methods have been proposed to evaluate fetal HRV [3]. One of the main methods used to perform time analysis is the root mean square of the successive differences (RMSSD) of R—R intervals in the electrocardiogram (ECG) [4], which corresponds to an estimation of the short-term variations analysis. The spectral analysis of the HRV can be quantified by the use of signal processing algorithms allowing us to obtain the heart rate frequency components of the following frequency ranges: very low frequency (VLF), a range of between 0 and 0.04 Hz; low frequency (LF), a range of between 0.04 and 0.15 Hz; and high frequency (HF), above 0.15 Hz [3]. The VLF range is related to thermoregulatory and endocrine activities; the LF range is related to the modulation of both the sympathetic and parasympathetic nervous systems, and is dependent on the feedback activity of baroreceptors; the HF range is centered on the respiratory frequency and is related to parasympathetic modulation only [5]. Using spectral analysis of fetal HRV, several previous studies have shown significant changes in the HF range in hypoxic fetuses, suggesting variations in parasympathetic activity [6–8].

However, existing signal processing methods are limited in terms of the lack of thresholds of the numerical indexes, which is primarily due to high inter-subject variability. This lack of thresholds to enable the prediction of ANS changes means that the methods cannot be used in obstetrical clinical practice.

We proposed a new original method that allowed us to obtain a numerical index relative to the parasympathetic activity of the ANS. Although most of the published methods are based on spectral or time analysis of HRV, this numerical index uses spectral analysis to filter the signal in order to keep only HF oscillations, and then computes the magnitude of these oscillations in the time domain.

The aim of this study was to evaluate the performance of our method compared to commonly used HRV analysis, regarding i) the ability to detect the parasympathetic activity of the fetus and ii) inter-subject variability.

## Materials and methods

### Surgical preparation

Near-term pregnant Charmoise sheep (INRA, Leudeville, France) with a gestational age of 123 +/- 2 days, term = 145 days, underwent the surgical procedure. The anesthetic and surgical technique protocols followed those previously described by our team or in previous studies [9,10]. The sheep fasted for 24 h before general anesthesia and surgery. Before surgery, sheep were placed in a supine position, anesthetized with an intravenous injection of xylazine (Sedaxylan®, CEVA Santé Animale, France), intubated, and maintained with isoflurane 2% (Aerrane®, Baxter, France). The uterus was exteriorized via a maternal midline laparotomy. Following hysterotomy, a 4Fr diameter catheter (Arrow®, USA) was inserted into the femoral artery and vein until reaching the abdominal aorta and the inferior vena cava, respectively,
through the femoral approach. Four precordial electrodes (MYWIRE 10®, MAQUET, Germany) were placed subcutaneously in order to detect an ECG signal. A 5F5-diameter catheter (Arrow®, USA) was placed into the amniotic cavity for measuring baseline pressure (intra-amniotic pressure; IAP). The fetal arterial catheter and intra-amniotic catheter were connected to pressure sensors (Pressure Monitoring Kit®, Baxter, France) that were connected to a blood pressure monitor (monitor Merlin, Hewlett Packard, Palo Alto, CA, USA). The mean arterial pressure (MAP) was measured from the blood pressure phasic signal and referenced to the IAP (calculated MAP = observed MAP—observed IAP). All hemodynamic data were recorded with the Physiotrace® data acquisition board (Estaris monitoring, Lille, France).

## Experimental procedure

The experiments began only after the sheep had had 48 hours of rest. Prior to injection (one per 24h per animal), an initial 30-minute period, called the stability period, was recorded. Fetal hemodynamic parameters (heart rate, blood pressure, and intra-amniotic pressure) were recorded every 5 minutes for the duration of the experiment (30 min). Recording of the electrocardiographic signal was continuous over the same period. Fetal arterial and venous blood gases were analyzed at the beginning and end of each protocol to ensure that the animals were healthy and representative of the model.

Two drugs were tested according to their effects on the autonomic nervous system: 2.5 mg atropine (Atropine, Aguettant, France), to evaluate the effect of parasympathetic nervous system inhibition, and 5 mg propranolol (Karnodyl®, Primius Lab, France), a non-selective beta-blocker with a sympatholytic effect [11], were intravenously injected into the femoral vein.

Following the recovery period, the ewe and the fetus were sacrificed via an injection of embutramide(T61®, Intervet, France).

## ECG signal preprocessing

The fetal ECG signal was recorded using the Physiotrace™ data acquisition board (Physiotrace™, Estaris Monitoring, Lille, France) [12], and ECG analysis was conducted off-line using an automatic R wave detection algorithm. R waves that missed detection, false detection, and ectopic beats were identified and replaced by a linear interpolation using a specific artifact removal algorithm [13]. As recommended, the resulting R-R series were re-sampled at an 8 Hz frequency using a linear interpolation algorithm [3].

## Fetal stress index analysis

Jeanne et al. studied HRV in a simulated environment and confirmed that HRV spectral analysis was sensitive to i) the basal heart rate, ii) the heart rate variations magnitudes and iii) the heart rate oscillations frequency [14]. In order to reduce inter-subject variability we therefore developed the algorithm described below;

The re-sampled R-R series was isolated in a 64-second moving window. In a first step, the signal is mean centered to eliminate the influence of the basal heart rate value.

The mean value (M) is computed as:

$$M = \frac{1}{N} \sum_{i=1}^{N} (RR_{i}),$$

Where RRᵢ represents the RR samples values and N represent the number of samples in the window. M is then subtracted from each sample of the window as: RR'ᵢ = (RRᵢ-M).
In a 2nd step, the signal is normalized to reduce the influence of the heart rate variations magnitude. The norm value (S) is then calculated as:

$$S = \sqrt{\sum_{i=1}^{N} (RR'_{i})^{2}}.$$

and each RR'$_{i}$ is divided by S: RR"$_{i}$ = RR'$_{i}$ / S.

The normalized RR" series is then high-pass filtered above 0.15 Hz, using a 4 coefficient Daubeuchies wavelet-based filter in order to keep only high frequency oscillations.

In order to eliminate the influence on the resulting oscillations frequency, we estimate the area under the envelope of the mean, normalized and filtered R-R series. The local maximum and minimum are detected via this filtered signal, and the lower and upper envelopes were plotted.

Jeanne et al. demonstrated that this method was not sensitive to basal heart rate, heart rate variations magnitude or heart rate HF oscillations frequency while keeping a strong correlation with commonly used HRV spectral measurements [14].

In order to increase the time-sensitivity of the method and to allow detecting rapid transient drops in the parasympathetic activity, we divided the envelope in four 16-seconds sub-areas A1, A2, A3, and A4 (Fig 1). AUC$_{min}$ was determined as the minimum value of A1, A2, A3, and A4 and fetal stress index (FSI) was defined as: FSI = 100 x (5.1 x AUC$_{min}$ + 1.2) / 12.8, where a = 5.1 and b = 1.2 were empirically determined in a dataset of more than 200 R-R records analysis in order to maintain good coherence between the visual aspect of the filtered R-R series (Fig 1) and a quantitative measurement ranging from 0 to 100.

Continuous computation of the FSI was assumed by sliding the moving window with a 1 second moving period. Mean values of the FSI were analyzed over two predefined 5-minute periods: 1) between 10 and 5 minutes before atropine or propranolol administration, and 2) between 5 and 10 minutes after atropine or propranolol administration.

The period around the time of injection was not analyzed in order to avoid any mathematical artifacts due to a rapid increase or decrease of fetal heart rate.

<figure><img src="image_5.png" /><figcaption>Fig 1. A normalized and filtered RR series (black curve, in normalized unit). Areas A1, A2, A3, and A4 are computed between the lower and upper limits (area in grey). The lowest area (AUCmin) is then selected (A4) and Fetal Stress Index is computed.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0180653.g001
## Time and spectral analysis

HF and LF spectral components were computed using a Fast Fourier Transform (FFT) performed in a 256-seconds moving window with an 8-Hz sampling frequency (3). The moving window was mean-centered and passed through a Bartlett ponderation window prior to the FFT computation. Continuous computation of HF and LF was assumed by sliding the moving window with a 1 second moving period. Mean values of HF and LF were analyzed over the two previously defined 5-minute periods. Normalized HF was also computing as HF/ (HF +LF). RMSSD time analysis was also computed for these two periods.

## Ethics

All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. The anesthesia, surgery, and experimentation protocols were performed as recommended by the Ministry of Higher Education and Research, and the study was approved by the Animal Experimentation Ethics Committee (APAFiS #5160).

## Statistical analysis

Given the small sample size, a Wilcoxon non parametrical test was used to compare data obtained before and after injections. Statistical significance was defined as $p < 0.05$. In order to evaluate the ability of HRV parameters with regard to differentiating the periods before and after injection, we also estimated the effect size (i.e. the magnitude of the difference) for each parameter, before and after drug injection.

Effect size was evaluated through the Cohen's d coefficient, defined as:

$$d = \frac{\bar{D}}{\sqrt{\frac{SS_D}{n-1}}}$$

Where $\bar{D}$ represents the mean difference between data before and after drug injection and $SS_D$ represents the sum of squared deviations (i.e. the sum of scores of deviations from the mean difference scores). The inter-subject variability was estimated by computing the coefficient of variance (%CV), defined as: $\%CV = 100 * SD/M$, where SD represents the standard deviation and M represents the mean value. %CV was computed for FSI, LF, HF, and RMSSD prior to atropine and propranolol injection. Hemodynamic, gazometric, and HRV parameters were expressed as median values (1st quartile– 3rd quartile). We used Statistical Package for the Social Sciences, version 20.0, software (IBM, Armonk, NY, USA).

## Results

Seven near-term pregnant sheep were chronically instrumented.

### Atropine

Seven experimental procedures were performed. Table 1 shows the data obtained before and after atropine injection. Following direct intravenous injection of atropine, HR significantly increased from 156(140.2–168.2) to 190(181.2–199.2) bpm (p = 0.028), while the mean arterial blood pressure remained the same (47(45–59)–52(49–59), p = 0.225); HF significantly decreased from 52.9(37.2–75) to 25.4(21.1–39.1) (p = 0.043); RMSSD significantly decreased from 8.7(6.3–9.2) to 4(3–6.5) (p = 0.018); HF/(HF+LF) significantly decreased from 0.68(0.61–0.78) to 0.53(0.47–0.65) (p = 0.018); and FSI significantly decreased from 57.3(52.4–61.2) to 40.4(36–46.2) (p = 0.018). %CV was important for HF, LF and RMSSD (respectively 33.4%,
<table><caption>**Table 1.** Spectral analysis (HF, LF and normalized HF), RMSSD and fetal stress index values before and after direct intravenous administration of 2.5 mg atropine. The results are presented as median values (1st–3rd quartiles). A Wilcoxon test was used for analysis. Statistical significance was set as p <0.05.</caption><thead><tr><th></th><th>Before</th><th>%CV before</th><th>After</th><th>%CV after</th><th>p</th><th>Effect size</th></tr></thead><tbody><tr><td>HF</td><td>52.9(37.2–75)</td><td>33.4</td><td>25.4(21.1–39.1)</td><td>30.2</td><td>0.043</td><td>1.13</td></tr><tr><td>LF</td><td>22.3(16–31.1)</td><td>35.2</td><td>18.1(14–41.8)</td><td>51.1</td><td>0.398</td><td>0.35</td></tr><tr><td>HF/(HF+LF)</td><td>0.68(0.61–0.78)</td><td>11.6</td><td>0.53(0.47–0.65)</td><td>18.3</td><td>0.018</td><td>0.92</td></tr><tr><td>RMSSD</td><td>8.7(6.3–9.2)</td><td>49.3</td><td>4(3–6.5)</td><td>37.6</td><td>0.018</td><td>1.13</td></tr><tr><td>FSI</td><td>57.3(52.4–61.2)</td><td>9.9</td><td>40.4(36–46.2)</td><td>17.3</td><td>0.018</td><td>1.49</td></tr></tbody></table>

HF = high frequency; LF = low frequency, RMSSD = root mean square of the successive differences;

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

35.2% and 49.3% before injection). %CV was far lower for HF(HF+LF) (%CV = 11.6 before injection) and the FSI (%CV = 9.9 before injection). Although the RMSSD, HF(HF+LF) and FSI p values were similar (p = 0.018), the effect size was greater for FSI. Indeed, as Gail *et al.* found, HF, HF/(HF+LF) and RMSSD were associated with a large effect size (>0.8) and FSI with a very large effect size (>1.3) [15].

Fig 2 shows the box plot of HF, LF, RMSSD and FSI variations before and after direct intravenous atropine administration.

## Propranolol

Seven experimental procedures were performed; one was excluded due to poor ECG quality. Table 2 shows the data obtained before and after propranolol injection. Following propranolol injection, ABP remained stable (40(36–44) vs (42(39–45), p = 0.092); HR decreased (161(155–165) vs 147(143–155), p = 0.046, as expected; HF, HF/(HF+LF), RMSSD and FSI did not vary significantly; and LF was lower (28(20.9–33.6) vs 12.7(6–16.2), p = 0.046).

Gazometric parameters were comparable before and after the two protocols (Table 3).

## Discussion

A computerized analysis of FHR variability could provide promising markers of fetal distress, monitored via analysis of the ANS [6,16]. We developed an original HRV analysis algorithm toobtain a numerical index relative to the parasympathetic activity of the ANS. We hypothesized that this algorithm could be adapted for fetal HRV analysis in order to obtain an FSI related to fetal parasympathetic activity, and we compared this to common methods of HRV analysis.

In order to evaluate the effect of inhibition of the parasympathetic nervous system on HRV, we used direct administration of 2.5 mg intravenous atropine, which is often used in this manner [4,11]. Frasch et al. observed an FHR increase after atropine injection in a sheep model (181 vs 236 bpm, p < .05) [4]. Moreover, the spectralpower of the RMSSD and HF bands decreased, suggesting a decrease in vagal modulationof FHR. The spectral power of the LF bandincreased, suggesting vagal-mediated disinhibition ofsympathetic modulation of FHR (4). In our study, the FSI, HF, HF/(HF+LF), and RMMSD dropped after the atropine injection, as expected, indicating that the FSI is highly related to parasympathetic activity.

We then evaluated the inhibition of the sympathetic system by the use of the beta-blocker, propranolol; Frasch et al. had previously described the blocking effect of this system via propranolol use in seven sheep [11]. As expected, we did not observe any change in the FSI, RMMSD, and HF following propranolol administration. LF decreased, suggesting vagal
<figure><img src="image_6.png" /><figcaption>Fig 2. High frequency, low frequency, root mean square of the successive differences and fetal stress index before and after administration of 2.5 mg atropine. * → p<0.05 versus before atropine injection.</figcaption></figure>

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

blockade. Thus, these experimental results, correlated with the findings of previous studies, allowed us to confirm the specificity of our index for parasympathetic tone.

The interest in HRV in the present study remains in the fact that it could be a new tool to predict fetal acidosis during labor [2]. In a model of repeated cord occlusion in sheep, an increase of RMSSD was observed in association with FHR decelerations, reflecting initial vagal activation during fetal asphyxia, and the authors concluded that RMSSD mayaccurately identify hypoxic fetuses at an early stage [4]. Chung et al.investigated HRV at very low, low, and

<table><caption>**Table 2.** Spectral analysis (HF, LF and normalized HF), RMSSD and fetal stress index values before and after direct intravenous administration of propranolol. The results are presented as median values (1st–3rd quartiles). A Wilcoxon test was used for analysis. Statistical significance was set as p <0.05.</caption><thead><tr><td></td><td>Before</td><td>%CV before</td><td>After</td><td>%CV after</td><td>p</td><td>Effect size</td></tr></thead><tbody><tr><td>HF</td><td>52.3(34.8–78.1)</td><td>53.5</td><td>25.4(17.5–42.6)</td><td>43.1</td><td>0.075</td><td>0.81</td></tr><tr><td>LF</td><td>28(20.9–33.6)</td><td>35.4</td><td>12.7(6–16.2)</td><td>51.8</td><td>0.046</td><td>1.25</td></tr><tr><td>HF/(HF+LF)</td><td>0.66(0.64–0.76)</td><td>9.22</td><td>0.70(0.65–0.79)</td><td>10.9</td><td>0.345</td><td>0.38</td></tr><tr><td>RMSSD</td><td>6.7(3.5–13.5)</td><td>73.6</td><td>6.4(4.1–18.6)</td><td>91.7</td><td>0.917</td><td>0.39</td></tr><tr><td>FSI</td><td>56.5(48.25–59.25)</td><td>10.9</td><td>61.5(51–76)</td><td>21.9</td><td>0.599</td><td>0.45</td></tr></tbody></table>

HF = high frequency; LF = low frequency, RMSSD = root mean square of the successive differences;

https://doi.org/10.1371/journal.pone.0180653.t002
<table><caption>**Table 3.** Blood gases before and after injection of atropine and propranolol, respectively. The results are presented as median values (1st–3rd quartiles).</caption><thead><tr><th rowspan="2"></th><th colspan="3">Atropine</th><th colspan="3">Propranolol</th></tr><tr><th>Before</th><th>After</th><th>p</th><th>Before</th><th>After</th><th>p</th></tr></thead><tbody><tr><td>pH</td><td>7.39 (7.37–7.41)</td><td>7.40 (7.38–7.41)</td><td>0.39</td><td>7.38(7.26–7.40)</td><td>7.37(7.27–7.39)</td><td>0.22</td></tr><tr><td>pO2(mmHg)</td><td>15.0 (15.0–17.0)</td><td>16.00 (15.0–16.0)</td><td>0.26</td><td>16.0 (14.5–17.5)</td><td>16.0 (14.7–17.2)</td><td>1</td></tr><tr><td>pCO2(mmHg)</td><td>45.1 (42.5–47.9)</td><td>46.5(41.4–47.1)</td><td>0.87</td><td>41.0(39.3–45.7)</td><td>40.7(39.1–43.9)</td><td>0.67</td></tr><tr><td>Lactates(mmol/L)</td><td>1.2(1.1–1.3)</td><td>1.2(1.1–1.3)</td><td>0.25</td><td>1.1(1.0–1.5)</td><td>1.1(0.5–1.4)</td><td>0.34</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0180653.t003

high frequencies during the last 2 hours before delivery [8], and observed a significant decrease in VLF, LF, and HF in acidotic fetuses. Siiraet al. also studied FHR variability during the last hour of labor according to the pH at birth [16]. Before delivery, they observed a significant decrease in HF variations in the acidotic group. Van Laar et al. found that acidotic fetuses had a significantly decreased in HF/(HF+LF) during the last 30 minutes of labor, compared with non-acidotic fetuses [6]. In all these studies, authors observed a drop in parasympathetic activity, as evaluated by the HF spectral component in acidotic fetuses. However, non-routine clinical use is proposed due to the lack of thresholds allowing the prediction of ANS changes related to fetal acidosis. In the present study, we compared common HRV analysis methods to the FSI and we showed that the FSI was associated with a higher effect size, showing that this parameter presented a higher power to detect parasympathetic inhibition. This very large effect size could help in the determination of thresholds for fetal acidosis prediction. Moreover, we compared dispersion of the data defined as the coefficient of variance %CV. The results highlighted the fact that inter-subject variability was far lower with the FSI than with spectral analysis (HF and LF) and RMSSD.

This experimental model is associated with some limitations. The surgery performed to position electrodes and catheters may induce fetal pain, which is known to be responsible for a decrease in parasympathetic activity. However, all animals underwent the same surgical procedure, and a rest period of 48 hours was observed between surgery and experiment.

## Conclusion

Although most of the published methods were based on a spectral or time analysis of HRV, our computerized index results from a filtering method, followed by a continuous analysis of the filtered signal magnitude in the time domain. In this study, we demonstrated that the FSI was associated with a lower inter-subject variability and a higher effect size, suggesting that this index could better detect a parasympathetic decrease than commonly used HRV analysis. However, our results are preliminary and must be confirmed in a genuine clinical context (i.e. fetal acidosis prediction).

## Supporting information

S1 File. RMSSD, spectral analysis and FSI datas.
(PDF)

## Acknowledgments

We thank Dr. Thomas Hubert and all the staff of the Research Experimental Department for their veterinary care and their expert assistance with sheep surgery.
## Author Contributions

**Conceptualization:** C. Garabedian, R. Logier, P. Deruelle, L. Storme, V. Houfflin-Debarge, J. De Jonckheere.

**Data curation:** C. Champion, E. Servan-Schreiber, E. Aubry, D. Sharma.

**Formal analysis:** L. Butruille.

**Investigation:** C. Garabedian, C. Champion, E. Servan-Schreiber.

**Methodology:** C. Garabedian, E. Aubry, D. Sharma, P. Deruelle, L. Storme, V. Houfflin-Debarge, J. De Jonckheere.

**Project administration:** P. Deruelle, L. Storme.

**Resources:** P. Deruelle, L. Storme.

**Software:** L. Butruille, J. De Jonckheere.

**Supervision:** P. Deruelle, L. Storme, V. Houfflin-Debarge, J. De Jonckheere.

**Validation:** C. Garabedian, J. De Jonckheere.

**Visualization:** C. Garabedian, C. Champion, E. Servan-Schreiber.

**Writing – original draft:** C. Garabedian.

**Writing – review & editing:** V. Houfflin-Debarge, J. De Jonckheere.

## References

1. Ugwumadu A. Understanding cardiotocographic patterns associated with intrapartum fetal hypoxia and neurologic injury. Best Pract Res Clin Obstet Gynaecol. 2013 Aug; 27(4):509–36. https://doi.org/10.1016/j.bpobgyn.2013.04.002 PMID: 23702579

2. Van Laar Jo, Porath Mm, Peters Ch, Oei Sg. Spectral analysis of fetal heart rate variability for fetal surveillance: review of the literature. Acta Obstet Gynecol Scand. 2008 Mar 1; 87(3):300–6. https://doi.org/10.1080/00016340801898950 PMID: 18307069

3. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electro-physiology. Eur Heart J. 1996 Mar; 17(3):354–81. PMID: 8737210

4. Frasch MG, Müller T, Weiss C, Schwab K, Schubert H, Schwab M. Heart rate variability analysis allows early asphyxia detection in ovine fetus. Reprod Sci Thousand Oaks Calif. 2009 May; 16(5):509–17.

5. Akselrod S, Amitayt Y, Lang RM, Mor-Avi V, Keselbrener L. Spectral analysis of left ventricular area variability as a tool to improve the understanding of cardiac autonomic control. Physiol Meas. 2000 May; 21 (2):319–31. PMID: 10847198

6. Van Laar J, Peters C, Vullings R, Houterman S, Bergmans J, Oei S. Fetal autonomic response to severe acidaemia during labour. BJOG Int J Obstet Gynaecol. 2010 Mar 1; 117(4):429–37.

7. Siira SM, Ojala TH, Vahlberg TJ, Rosén KG, Ekholm EM. Do spectral bands of fetal heart rate variability associate with concomitant fetal scalp pH? Early Hum Dev. 2013 Sep; 89(9):739–42. https://doi.org/10.1016/j.earlhumdev.2013.05.007 PMID: 23809772

8. Chung DY, Sim YB, Park KT, Yi SH, Shin JC, Kim SP. Spectral analysis of fetal heart rate variability as a predictor of intrapartum fetal distress. Int J Gynaecol Obstet Off Organ Int Fed Gynaecol Obstet. 2001 May; 73(2):109–16.

9. Houfflin-Debarge V, Sabbah-Briffaut E, Aubry E, Deruelle P, Alexandre C, Storme L. Effects of environmental tobacco smoke on the pulmonary circulation in the ovine fetus. Am J Obstet Gynecol. 2011 May; 204(5):450.e8–450.e14.

10. Aubry E, Fayoux P, Jani J, Deprest J, Deruelle P, Houfflin-Debarge V, Storme L. Tracheal occlusion alters pulmonary circulation in the fetal lamb with normally developing lungs. J Pediatr Surg. 2013 Mar; 48(3):481–7. https://doi.org/10.1016/j.jpedsurg.2012.08.024 PMID: 23480900

11. Frasch MG, Müller T, Hoyer D, Weiss C, Schubert H, Schwab M. Nonlinear properties of vagal and sympathetic modulations of heart rate variability in ovine fetus near term. Am J Physiol Regul
Integr Comp Physiol. 2009 Mar; 296(3):R702-7. https://doi.org/10.1152/ajpregu.90474.2008 PMID: 19109371

12. De Jonckheere J, Logier R, Dassonneville A, Delmar G, Vasseur C. PhysioTrace: An efficient toolkit for biomedical signal processing. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Annu Conf. 2005; 7:6739-41.

13. Logier R, De Jonckheere J, Dassonneville A. An efficient algorithm for R-R intervals series filtering. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Annu Conf. 2004; 6:3937-40.

14. Jeanne M, Logier R, De Jonckheere J, Tavernier B. Validation of a graphic measurement of heart rate variability to assess analgesia/nociception balance during general anesthesia. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Annu Conf. 2009; 2009:1840-3.

15. Sullivan GM, Feinn R. Using Effect Size-or Why the P Value Is Not Enough. J Grad Med Educ. 2012 Sep; 4(3):279-82. https://doi.org/10.4300/JGME-D-12-00156.1 PMID: 23997866

16. Siira SM, Ojala TH, Vahlberg TJ, Jalonen JO, Välimäki IA, Rosén KG, Ekholm EM. Marked fetal acidosis and specific changes in power spectrum analysis of fetal heart rate variability recorded during the last hour of labour. BJOG Int J Obstet Gynaecol. 2005 Apr; 112(4):418-23.