Developmental milestones of the autonomic nervous system revealed via longitudinal monitoring of fetal heart rate variability

Abstract

Background

Fetal heart rate variability (fHRV) of normal-to-normal (NN) beat intervals provides high-temporal resolution access to assess the functioning of the autonomic nervous system (ANS).

Aim

To determine critical periods of fetal autonomic maturation. The developmental pace is hypothesized to change with gestational age (GA).

Study design

Prospective longitudinal observational study.

Subjects

60 healthy singleton fetuses were followed up by fetal magnetocardiographic heart rate monitoring 4–11 times (median 6) during the second half of gestation.

Outcome measure

FHRV parameters, accounting for differential aspects of the ANS, were studied applying linear mixed models over four predefined pregnancy segments of interest (SoI: <27; 27+0–31+0; 31+1–35+0; >35+1 weeks GA). Periods of fetal active sleep and quiescence were accounted for separately.

Results

Skewness of the NN interval distribution VLF/LF band power ratio and complexity describe a saturation function throughout the period of interest. A decreasing LF/HF ratio and an increase in pNN5 indicate a concurrent shift in sympathovagal balance. Fluctuation amplitude and parameters of short-term variability (RMSSD, HF band) mark a second acceleration towards term. In contrast, fetal quiescence is characterized by sequential, but low-margin transformations; ascending overall variability followed by an increase of complexity and superseded by fluctuation amplitude.

Conclusions

An increase in sympathetic activation, connected with by a higher ability of parasympathetic modulation and baseline stabilization, is reached during the transition from the late 2nd into the early 3rd trimester. Pattern characteristics indicating fetal well-being saturate at 35 weeks GA. Pronounced fetal breathing efforts near-term mirror in fHRV as respiratory sinus arrhythmia.

About the Authors:

Uwe Schneider

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Obstetrics, Division of Prenatal Diagnostics and Fetal Physiology, Jena University Hospital, Jena, Germany

ORCID logo http://orcid.org/0000-0002-7920-4350

Franziska Bode

Roles Data curation, Formal analysis, Visualization, Writing – original draft

Affiliation: Department of Obstetrics, Division of Prenatal Diagnostics and Fetal Physiology, Jena University Hospital, Jena, Germany

Alexander Schmidt

Roles Formal analysis, Methodology, Validation

Affiliation: Hans Berger Clinic of Neurology, Biomagnetic Center, Jena University Hospital, Jena, Germany

Samuel Nowack

Roles Data curation, Formal analysis, Methodology

Affiliation: Hans Berger Clinic of Neurology, Biomagnetic Center, Jena University Hospital, Jena, Germany

Anja Rudolph

Roles Conceptualization, Data curation, Investigation, Methodology

Affiliation: Department of Obstetrics, Division of Prenatal Diagnostics and Fetal Physiology, Jena University Hospital, Jena, Germany

Eva-Maria Doelcker

Roles Formal analysis, Methodology, Validation, Writing – review & editing

Affiliations Hans Berger Clinic of Neurology, Biomagnetic Center, Jena University Hospital, Jena, Germany, Institute of Biomedical Engineering and Informatics, Technical University, Ilmenau, Germany

Peter Schlattmann

Roles Software, Supervision, Validation

Affiliation: Institute for Medical Statistics, Computer and Data Sciences, Jena University Hospital, Jena, Germany

Theresa Götz

Roles Formal analysis, Methodology, Software, Validation, Writing – review & editing

Affiliations Hans Berger Clinic of Neurology, Biomagnetic Center, Jena University Hospital, Jena, Germany, Institute for Medical Statistics, Computer and Data Sciences, Jena University Hospital, Jena, Germany

Dirk Hoyer

Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliation: Hans Berger Clinic of Neurology, Biomagnetic Center, Jena University Hospital, Jena, Germany

Introduction

Developing the autonomic ability to adapt to varying amounts of supply and demand in the organism is one of the key necessities during fetal maturation. This process of autonomic regulation in utero reaches a pace that may never be observed again later in life. With regard to the concept of ‘developmental origins of adult disease’, maturation in utero may be highly susceptible to acute and chronic influences with long-lasting consequences [1,2].

Autonomic capacity is mirrored in cardiovascular regulation. Since the autonomic nervous system (ANS) is involved in the regulation of nearly all organs, cardiac autonomic control, which is accessible from fetal heart rate patterns (fHRP), provides relevant diagnostic and prognostic information, not only on current fetal well-being, but also on functional autonomic brain age [3,4].

On the basis of neuroanatomical studies, a sequential development of the ANS has been suggested [5]: Both the differentiation of the lateral zone of the hypothalamus, as well as an increasing myelination of the vagal nerve, indicate a first propelling maturation of the parasympathetic branch at the end of the second trimester [6,7]. From 32 weeks of gestation (WGA) onwards, both baroreflex responsiveness and the observation of increasing respiratory sinus arrhythmia indicate advancing modulatory capacities of the vagal nerve. Simultaneously, the coordination of fetal movements and requirement driven heart rate accelerations can be seen as the expression of advancing sympathetic responsiveness [8]. Therefore, different essential periods of ANS developmental dynamics such as (a) the transitional period from the late second into the early third trimester, (b) the time around 30 to 32 WGA, and (c) the near-term pregnancy beyond 35 WGA, can be expected.

HRV analysis is nowadays an internationally standardized procedure to assess the autonomic regulation of the cardiovascular system [9]. It involves the extraction of linear and non-linear (complexity) parameters the train of cardiac beat-to-beat cycles. Interpretation requires a priori knowledge of the physiological background of the heart beat sequences Linear analyses have in previous studies been proven viable in the fetus, both in time and frequency domains.) [9–13]. The temporal resolution of the cardiac inter-beat intervals considerably influences the results. Therefore, an electrophysiological method of recording, i.e. fetal magnetocardiography (fMCG), bears methodical advantages [10]. Recently, HRV indices that reflect general principles of maturation, adaptation and self-organization, such as increasing fluctuation amplitude, increasing complexity and the formation of characteristic patterns, have been investigated. [4,14, 15].

In addition to the physiological changes with increasing WGA, the concurrent fetal neuro-behavioral state (state of fetal activity) must be taken into account when interpreting fHRP. Rest/activity cycles have been described from about 23 WGA onwards with progressing discernibility between 26 and 32 WGA. Synchronization of neuro-behavioral variables like fHRP, body movements and eye-movements are markers of the developmental integrity of the fetus [15–17]. From 32 WGA onwards state synchronization was observed in up to 80% of observation times.

Indices of both sympathetic activation and vagal modulation increase with growing fetal age [4,18–20]. The transitional period between the second and the third trimesters of pregnancy is characterized by the appearance of physiological decelerations, a steeply enhancing vagal regulation, followed by an increase in overall variability and complexity [10,21,22]. Beyond 30 weeks GA development is characterized by a widespread distribution of inter-individual characteristics of fHRV. This distribution can in part be related to the concurrent fetal neurobehavioral states [18]. There is a strong interconnection between fHR accelerations and intentional body movements [23].

Therefore, several gestational segments of different maturational pace are proposed: a developmental surge of the sympathetic nervous system from the late second into the early third trimester, followed by a transitional period around 30 WGA [21]. The late third trimester is characterized by (1) an increase in heart rate patterns indicative of respiratory sinus arrhythmia in association with fetal thoracic movements and (2) the synchronization of accelerative fHRP in association with fetal activity. A longitudinal study overcoming the inter-individual variance and taking these different stages of maturation into account is pending. The aim of the present study is therefore, to augment the evidence of these critical periods of normal autonomic development in utero from fHRV analysis.

Materials and methods

Study population

This is a prospective longitudinal observational study in a cohort of healthy pregnant women and their singleton fetuses, with normal course of pregnancies according to standard maternity care in Germany and normal maternal and perinatal outcomes. All subjects voluntarily participated and gave their written consent after explicit information on the character of the investigation. The study was approved by the local Ethics Committee of the Medical Faculty of the Friedrich Schiller University of Jena (Reg. Nr 1104-04/03 and amendments).

From the initially recruited 61 subjects one had to be excluded retrospectively because of a prenatally unknown congenital heart defect. All the 60 neonates were reported to be healthy and mature at birth. As a priori exclusion criteria we considered: maternal age < 18 y, multiplets, active preterm labor, cardiovascular disease, administration of any medication with known cardiovascular effects, diseases from the gestational hypertensive spectrum, maternal diabetes, substance abuse (nicotine, alcohol, drugs), hints of fetal distress from conventional fetal monitoring (CTG, ultrasound), intrauterine growth restriction/abnormal findings in Doppler ultrasound, known chromosomal or other congenital abnormalities, as well as previous exposure to synthetic steroids to enhance fetal lung maturation.

312 monitoring sessions in non-stress situations from 60 fetuses between 19 and 39 completed WGA were included in our analysis. Initially, a biweekly pattern of follow-up was intended. The numbers of repetitive sessions ranged from 4–11 over the second half of pregnancy (mean 6 sessions/subject).

Therefore, on the basis of our presuppositions, the gestational period of interest was subdivided into four developmental segments of interest (SoI {1–4}): <27 WGA; 27+0–31+0 WGA; 31+1–35+0 WGA; >35+1 WGA, constituting respective developmental milestones of interest. This resulted in at least one monitoring session in each of the 4 SoI per subject.

Methods

All magnetocardiographic (fMCG) monitoring sessions were performed over 30 min, sampled at 1024 Hz in a magnetically shielded room at the Biomagnetic Center, Department of Neurology, Jena University Hospital using the vector-magnetograph ARGOS 200 (ATB, Chieti, Italy) during normal daytime working hours and according to published standards of procedure [4,24–26]. The pregnant women were positioned supine or with a slight twist to either side to prevent compression of the inferior vena cava by the pregnant uterus. The Dewar containing the magnetometers was positioned with its curvature above the fetal heart after sonographic localization and as close to the maternal abdominal wall as possible without direct contact [27]. A maternal 3-lead-ECG was additionally used to distinguish between fetal and maternal cardiac representation in the multichannel magnetic signal in cases of equivocalness.

After removal of maternal cardiac signals by independent component analysis [23], the fetal heartbeats were automatically detected and normal-to-normal (NN) beat intervals series were calculated. The NN series were screened for artifacts, arrhythmias and nonstationarities both automatically and visually. The rate of interpolated artifacts was below 5% in all cases [9]. The following distribution of recordings per pregnancy SoI served for further analysis: SoI {1} 84; SoI {2} 82; SoI {3} 80; SoI {4} 66.

Fetal heart rate variability analysis (fHRV) was performed using the parameters as listed in Table 1 [4,9,13,14,18,22]. Frequency analysis was based on the frequency bands as described by David et al. [13]. Parameters were grouped with regard to their representation of developmental and regulatory aspects of the ANS (see Table 1).

[Figure omitted. See PDF.]

Table 1. Parameters of fetal heart rate variability.

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

In addition to the analysis of the entire 30 min recordings (unclassified data), 10 min sections of active and quiet fetal neurobehavioral states (quiescence/1F and active sleep/2F according to [15]) were selected from visual inspection of the heart rate pattern printout after a consensus decision by three independent obstetricians (Table 2). Episodes of active sleep could be found in 263 of the 312 recorded datasets (84.3%) and fetal quiescence of at least 10 min length in 65 traces (20.8%; SoI{1} 20, SoI{2} 15, SoI{3} 13, SoI{4} 17 recordings). Generally, a single episode per state and dataset was considered. This standardized procedure was performed blinded to the heart rate analysis, according to the standard criteria previously described and extended to gestational ages prior to 32 WGA [4,15,18,20]. Episodes of active awakeness were not considered separately due to their rare occurrence.

[Figure omitted. See PDF.]

Table 2. Fetal state of activity.

https://doi.org/10.1371/journal.pone.0200799.t002

The original 312 recordings (RR interval series over 30 min) and the supporting information file, including the recording IDs and gestational age at each recording, constitute the minimal data set and will be made available by a download directory of the Biomagnetic Center, Department of Neurology, Jena University Hospital (please, contact [email protected] using the keyword biomagnetic data for details)

Statistical analysis

Our described longitudinal design of non-independent recordings at different time points and with inconsistent representation of the single subjects required the development of a special statistical model that sufficiently takes all these aspects into consideration [28]. We therefore applied a linear mixed model; the so-called PRO MIXED (SAS 3.4. Basic Edition) [29,30]. The PRO MIXED model offers the opportunity to perform both the restricted maximum likelihood (REML) and the maximum likelihood methods to account for inter- and intra-correlation between the datasets. In the described investigation the REML method was used for the smaller estimation error and combined with the between-within-method to calculate the degrees of freedom in non-independent data sets from one subject [29]. Within the 4 separate SoI the ‘Estimator’, as a measure of developmental dynamics, was calculated for each parameter from both the state independent (30 min) and pre-selected state-dependent NN interval segments (10 min). To account for differences between the predefined SoI the linear ‘Contrast’ was applied. The Quality of the model was confirmed by zero-model likelihood ratio, significance was tested by χ2 test (p<0.05, trend p<0.10).

Results

The parameters indicating fluctuation amplitude/sympathetic activation, pattern formation and complexity were consistent, when comparing the results from unclassified 30-min-recordings and those segments representing active fetal sleep (Figs 1 and 2). An overview of the numerical data is given in Table 3.

[Figure omitted. See PDF.]

Fig 1.

Developmental course of ACTAMP20 (a-c) and SDNN (d-f) in the linear mixed model. (a, d) unclassified datasets, (b, e) fetal active sleep, (c, f) fetal quiescence. x–Segment of Interest {< 27 WGA; 27–30 WGA; 31–35 WGA; >35 WGA}; completed weeks of GA; y–Estimator in [ms] units. Numerical data see Table 2. Significant (p<0.05) contrasts are marked by *, trends (p<0.10) by +.

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

[Figure omitted. See PDF.]

Fig 2.

Developmental course of gMSE(3) (a, b) and skewness (c, d) in the linear mixed model. (a) unclassified datasets: SoI {1–2} χ2 = 11,81; p<0.001, SoI {2–3} χ2 = 10.04; p = 0.0013; (b) fetal active sleep SoI {1–2} χ2 = 6.7, p = 0.0096, SoI {2–3} χ2 = 12.17, p<0.001; (c) unclassified datasets: SoI {1–2} χ2 = 70.31, p<0.0001, SoI {2–3} χ2 = 21.29, p<0.001, SoI {3–4} χ2 = 3.16, p = 0.076; (d) fetal active sleep: SoI {1–2} χ2 = 35.50, p<0.0001, SoI {2–3} χ2 = 21.13, p<0.0001. x–Segment of Interest {< 27 WGA; 27–30 WGA; 31–35 WGA;>35 WGA}; completed weeks of GA; y–Estimator in [bitnorm] Significant (p<0.05) contrasts are marked by *, trends (p<0.10) by +. Results for ‘fetal quiescence’ are not shown because they did not reach significant levels (gMSE(3): SoI {2–3}, χ2 = 3.41; p = 0.07).

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

ACTAMP20 and SDNN display a two-step-increase that is pronounced between SoI {1–2} and SoI {3–4} (Fig 1, Table 3). The observed increase in overall fluctuation amplitude is flanked by a steady increment of short term variability throughout the gestational period of interest (pNN5, Fig 3, Table 3). In addition, RMSSD (unclassified χ2 = 3.59, p = 0.06; active sleep χ2 = 7.82; p = 0.005) and HF (Table 3) show their most remarkable increase between SoI {3–4}. In contrast, both skewness and complexity follow a saturation function towards term. VLF/LF (unclassified, increase, p = 0.0098) and LF/HF (active sleep, decrease, p = 0.0624) ratios changed between SoI {1–2} (Table 2).

[Figure omitted. See PDF.]

Fig 3. Developmental course of pNN5 in the linear mixed model.

(a) unclassified datasets, (b) fetal active sleep, (c) fetal quiescence. x–Segment of Interest {< 27 WGA; 27–30 WGA; 31–35 WGA; >35 WGA}; completed weeks of GA; y–Estimator in [%]. Numerical data see Table 2. Significant (p<0.05) contrasts are marked by *, trends (p<0.10) by +.

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

[Figure omitted. See PDF.]

Table 3. Contrast analysis of fHRV parameters between pregnancy segments of interest.

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

With restriction to the low number of observed periods of fetal quiescence, fHRV changes between SoI could only be demonstrated in a few cases as trends: SDNN between SoI {1–2} (χ2 = 2.95; p = 0.09, Fig 2C) and {3–4} (χ2 = 2.18; p = 0.14; Fig 2C), ACTAMP20 between SoI {3–4} (χ2 = 3.14; p = 0.08; Fig 2F), VLF between SoI {1–2} (χ2 = 5.07; p = 0.024) and gMSE(3) between SoI {2–3} (χ2 = 3.41; p = 0.07). The increase in pNN5 between SoI {2–3} (χ2 = 2.57, Fig 1C) did not reach an appropriate level of significance (p = 0.11).

Discussion

In this longitudinal intra-individual monitoring study we were able to demonstrate several key features of normal fetal autonomic maturation during the second half of gestation (Fig 4): The transitional period from the late second into the early third trimester is characterized by an acceleration of autonomic maturation. The general principles of maturation as described earlier are fulfilled by the increase in fluctuation amplitude (ACTAMP20, SDNN), complexity (gMSE(3)), the formation of characteristic patterns, like a shift from decelerating towards accelerative heart rate patterns (skewness), and a stabilization of baseline heart rate (VLF/LF). These characteristics are accompanied by an increase in parasympathetic modulation (pNN5, LF/HF ratio) that is more pronounced during periods of active sleep when higher sympathetic activation is self-evident (ACTAMP20, SDNN) [18]. Additionally, the period close to term is characterized by a steep increase in high frequency variability (HF, RMSSD) that most likely represents respiratory sinus arrhythmia (RSA) [31,32]. RSA is associated with breathing efforts of the maturing fetus. These thoracic movements are characteristically not confined to either of the fetal neurobehavioral states but occur more frequently during periods of active sleep [33].

[Figure omitted. See PDF.]

Fig 4. Autonomic maturation as reflected by heart rate regulation.

Schematic summary of the findings for illustration purposes. The scales of the individual components contributing to heart rate regulation are arbitrary and do not reflect absolute values, but coincidental qualitative changes throughout the second half of gestation.

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

Our results support the concept, that increasing sympathetic activation is accompanied by increased vagal modulation in states of preserved sympatho-vagal balance, as studied here [18]. The previous theory of ‘polyvagal development’ [5] which sees a sequential development of vagal, followed by sympathetic and again vagal maturation, can only in parts be confirmed. Our results nicely fit the two- step parasympathetic development as described from neuroanatomic studies. Cheng et al. and Porges et al. describe the maturation of primitive neurons in the parasympathetic brain stem nuclei and the development of the myelin sheaths of the vagal nerve reaching their maximum pace around 30 WGA [5,6,34,35]. In a second step, RSA is coordinated further beyond 30 WGA. Suess et al. were able to demonstrate a lack of coordination in the swallowing process to RSA in preterm infants prior to 30 WGA, in comparison to a cohort beyond 32 weeks of developmental age [36].

With focus on periods of fetal quiescence very little changes were observed during the gestational period of interest, confirming previous observations from cross-sectional studies [18,22]. Both vagal modulation (pNN5) and complexity (gMSE(3)) did solely increase around 32 WGA at the transition from SoI {2–3}. Hence, characterizing the apparent fetal state of neurobehavioral activity does essentially influence analysis of autonomic function. In previous cross-sectional studies, we were able to demonstrate an increase in beat-to-beat complexity during fetal quiescence with advancing GA [22]. Even though this is, to our knowledge, the largest prospectively conceived longitudinal study performed on fMCG heart rate monitoring so far, we finally have to contemplate, that episodes of fetal quiescence might have been underrepresented in the study population to reach significant results.

In previous work we could demonstrate, that the use of parameters like pNN5, ACTAMP20, skewness or gMSE(3) more consistently described fetal autonomic maturation than parameters primarily proposed in adult HRV (classical according to Task Force 1996) [4,9,22]. With respect to the frequency domain, a characteristic shift towards higher frequency variability from the LF to the HF spectrum has previously been described. Hence, a loss of additional information in the VLF/LF ratio from about 32 WGA onwards is to be expected [11,37].

The different HRV indices mainly address different physiological aspects. There are some known redundancies concerning vagal activity (RMSSD, pNN5, HF) and overall fluctuation range (SDNN, ACTAMP20). Therefore, we refrained from performing multiple-testing adjustments for different physiological aspects. The linear mixed model statistics inherently consider all effects and repetitions for the chosen HRV indices.

Obviously, visual analysis of the fetal neurobehavioral state of activity, solely from heart rate patterns, doesn’t represent the gold standard [15,16]. Nevertheless, the described method has been consistently used over the last years, was able to cluster fHRV results accordingly and has been adapted by research colleagues worldwide [4,14,18–22,37,38]. More recent efforts have been made to include fetal movements into analysis [23,39]. The resemblance of results between unclassified data and those, attributed to active fetal sleep in our present study, makes it unlikely that additional information on state appearance will cast doubt on the general principles elaborated here.

Continuous antenatal heart rate tracing is based on two physical principles: (i) Doppler ultrasound-based cardiotocography (CTG), although not able to identify each individual heart beat delivers a reasonable temporal resolution for daily clinical practice and algorithms for computerized analysis, i.e. according to the Dawes/Redman criteria, aid in reducing its shortcomings [40–42]. (ii) In comparison to the CTG, electrophysiological methods such as fetal electrocardiography (fECG) or fetal magnetocardiography (fMCG) allow the precise QRS complex detection.

Magnetocardiography, though sophisticated and providing optimum temporal resolution, failed to reach worldwide significance throughout its recent 40 years of application [43]. The method is far too costly and specialist-driven to be widely distributed. Moreover, fECG is burdened with severe signal attenuation due to the vernix caseosa, making it difficult to apply between 28 and 34 WGA [44–45]. Therefore, efforts should be encouraged to determine how much of the gathered information might be detectable by routine clinical applications.

Conclusion

In conclusion, our study demonstrates changes in the developmental pace of the ANS during intrauterine life that can be monitored by selectively applying parameters of fHRV. An increase in sympathetic activation is accompanied by a higher ability of parasympathetic modulation and this maturational milestone is reached during the transitional period from the late second into the early third trimester. The near term period is characterized by pronounced periods of fetal breathing efforts that mirror in fetal heart rate regulation. The pattern characteristics of fetal active sleep are stable from 35 WGA onwards.

Acknowledgments

The recordings, preprocessing, and technical service were done Franziska Jaenicke, Esther Heinicke, Ulrike Wallwitz, Isabelle Kynass, Janine Zöllkau, Kathrin Kumm, Sophia Leibl, Adelina Pytlik, Stefan Claus, Tina Radtke, Ralph Huonker. English editing by a native speaker was performed by Anna K Multhaup.

Citation: Schneider U, Bode F, Schmidt A, Nowack S, Rudolph A, Doelcker E-M, et al. (2018) Developmental milestones of the autonomic nervous system revealed via longitudinal monitoring of fetal heart rate variability. PLoS ONE 13(7): e0200799. https://doi.org/10.1371/journal.pone.0200799

References

1. Barker DJ. In utero programming of chronic disease. Clinical Science. 1998; 95: 115–28. pmid:9680492

2. Van den Bergh BR, Mulder EJ, Mennes M, Glover V. Antenatal maternal anxiety and stress and the neurobehavioural development of the fetus and child: links and possible mechanisms. A review. Neurosci Biobehav Rev. 2005; 29: 237–258. pmid:15811496

3. FIGO 2011 (Federation Internationale de Gynecologie et d'Obstetrics), Maternity Guidelines, Federal Gazette, Institute for Quality and Efficiency in Health Care, Germany

4. Hoyer D, Tetschke F, Jaekel S, Nowack S, Witte OW, Schleussner E, Schneider U. Fetal functional brain age assessed from universal developmental indices obtained from neuro-vegetative activity patterns. PloS One. 2013; 8: e74431. pmid:24058564

5. Porges SW, Furman SA. The Early Development of the Autonomic Nervous System Provides a Neural Platform for Social Behavior: A Polyvagal Perspective. Infant Child Dev. 2011; 20: 106–118. pmid:21516219

6. Cheng G, Zhou X, Qu J, Ashwell KW, Paxinos G. Central vagal sensory and motor connections: human embryonic and fetal development. Auton Neurosci. 2004; 114: 83–96. pmid:15331048

7. Koutcherov Y, Mai JK, Paxinos G. Hypothalamus of the human fetus. J Chem Neuroanat. 2003; 26: 253–270. pmid:14729128

8. DiPietro JA, Hodgson DM, Costigan KA, Hilton SC, Johnson TR. Fetal neurobehavioral development. Child Dev. 1996; 67: 2553–2567. pmid:9022256

9. Force Task. 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 Electrophysiology. European Heart Journal. 1996; 17: 354–381. pmid:8737210

10. Van Leeuwen P, Lange S, Bettermann H, Gronemeyer D, Hatzmann W. Fetal heart rate variability and complexity in the course of pregnancy. Early Hum Dev. 1999; 54: 259–269. pmid:10321792

11. Zhuravlev YE, Rassi D, Mishin AA, Emery SJ. Dynamic analysis of beat-to-beat fetal heart rate variability recorded by SQUID magnetometer: quantification of sympatho-vagal balance. Early Hum Dev. 2002; 66: 1–10. pmid:11834343

12. Wakai RT, Assessment of fetal neurodevelopment via fetal magnetocardiography, Exp Neurol. 2004; 190, Suppl. 1: S65–71.

13. David M, Hirsch M, Karin J, Toledo E, Akselrod S. An estimate of fetal autonomic state by time-frequency analysis of fetal heart rate variability. Journal Appl Physiol. 2007; 102: 1057–1064.

14. Hoyer D, Kowalski EM, Schmidt A, Tetschke F, Nowack S, Rudolph A, et al. Fetal autonomic brain age scores, segmented heart rate variability analysis, and traditional short term variability. Front Hum Neurosci. 2014; 8: 948. pmid:25505399

15. Hoyer D, Schneider U, Kowalski EM, Schmidt A, Witte OW, Schleussner E et al. Validation of functional fetal autonomic brain age score fABAS in 5 min short recordings. Physiol Meas. 2015; 36: 2369–2378. pmid:26489779

16. Nijhuis J, Prechtl HF, Martin CB Jr., Bots RS. Are there behavioural states in the human fetus? Early Hum Dev. 1982; 6: 177–195. pmid:7094856

17. Pillai M, James D. The development of fetal heart rate patterns during normal pregnancy. Obstet Gynecol. 1990; 76: 812–816. pmid:2216230

18. Schneider U, Frank B, Fiedler A, Kaehler C, Hoyer D, Liehr M, Haueisen J, Schleussner E. Human fetal heart rate variability-characteristics of autonomic regulation in the third trimester of gestation. Journal Perinat Med.2008; 36: 433–441.

19. Lange S, Van Leeuwen P, Schneider U, Frank B, Hoyer D, Geue D et al. Heart rate features in fetal behavioural states. Early Hum Dev. 2009; 85: 131–135. pmid:18757143

20. Brändle J, Preissl H, Draganova R, Ortiz E, Kagan KO, Abele H et al. Heart rate variability parameters and fetal movement complement fetal behavioral states detection via magnetography to monitor neurovegetative development. Front Hum Neurosci. 2015; 9: 147. pmid:25904855

21. Van Leeuwen P, Geue D, Lange S, Hatzmann W, Gronemeyer D. Changes in the frequency power spectrum of fetal heart rate in the course of pregnancy. Prenat Diagn. 2003; 23: 909–916. pmid:14634977

22. Hoyer D, Heinicke E, Jaekel S, Tetschke F, DiPietro Paolo D, Haueisen J et al. Indices of fetal development derived from heart rate patterns. Early Hum Dev. 2009; 85: 379–386. pmid:19188032

23. Schmidt A, Schneider S, Witte OW, Schleussner E, Hoyer D. Developing fetal motor-cardiovascular coordination analysed from multi-channel magnetocardiography. Physiol Meas. 2014; 35: 1943–59. pmid:25229562

24. Kiefer-Schmidt I, Lim M, Wacker-Gussmann A, Ortiz E, Abele H, Kagan KO et al. Fetal magnetocardiography (fMCG): moving forward in the establishment of clinical reference data by advanced biomagnetic instrumentation and analysis. J Perinat Med. 2012; 40: 277–86. pmid:22505507

25. Lange S, Van Leeuwen P, Geue D, Hatzmann W, Gronemeyer D. Influence of gestational age, heart rate, gender and time of day on fetal heart rate variability. Med Biol Eng Comput. 2005; 43: 481–486. pmid:16255430

26. Grimm B, Haueisen J, Huotilainen M, Lange S, Van Leeuwen P, Menendez T et al. Recommended standards for fetal magnetocardiography, Pacing Clin Electrophysiol. 2003; 26: 2121–2126. pmid:14622314

27. Sturm R, Muller HP, Pasquarelli A, Demelis M, Erne SN, Terinde R et al. Multi-channel magnetocardiography for detecting beat morphology variations in fetal arrhythmias. Prenat Diagn. 2004; 24: 1–9. pmid:14755401

28. Davis CS, Statistical Methods for the Analysis of Repeated Measurements, 1. Aufl. New York: Springer-Verlag. 2002: 416.

29. Ramon C, Littell GAM, Stroup WW, Wolfinger R. SAS System for Mixed Models. North Carolina, USA: SAS Institute Inc. 2006.

30. Wright SP and Russel D,. Repeated Measures Analysis Using Mixed Models: Some Simulation Results. Modelling Longitudinal and Spatially Correlated Data. New York: Springer Verlag. 1997: 147–157.

31. Gustafson KM, May LE, Yeh HW, Million SK, Allen JJ. Fetal cardiac autonomic control during breathing and non-breathing epochs: the effect of maternal exercise. Early Hum Dev. 2012; 88: 539–46. pmid:22264436

32. Wakai RT, Wang M, Pedron SL, Reid DL, Martin CB Jr.. Spectral analysis of antepartum fetal heart rate variability from fetal magnetocardiogram recordings. Early Hum Dev. 1993; 35: 15–24. pmid:8293713

33. Mulder EJH, Boersma M, Meeuse M, Van der Waal M, Van de Werd , Visser GHA. Patterns of breathing movements in the near-term human fetus: relationship to behavioural states. Early Hum Dev. 1994; 36: 127–35. pmid:8200321

34. Dubois J, Dehaene-Lambertz G, Kulikova S, Poupon C, Huppi PS, Hertz-Pannier L. The early development of brain white matter: a review of imaging studies in fetuses, newborns and infants. Neuroscience. 2014; 276: 48–71. pmid:24378955

35. Holland MA, Miller KE, Kuhl E. Emerging Brain Morphologies from Axonal Elongation. Ann Biomed Eng. 2015; 43: 1640–1653. pmid:25824370

36. Suess PE, Alpan G, Dulkerian SJ, Doussard-Roosevelt J, Porges SW et al. Respiratory sinus arrhythmia during feeding: a measure of vagal regulation of metabolism, ingestion, and digestion in preterm infants. Dev Med Child Neurol. 2000; 42: 169–173. pmid:10755456

37. Schneider U, Schleussner E, Fiedler A, Jaekel S, Liehr M, Haueisen J et al. Fetal heart rate variability reveals differential dynamics in the intrauterine development of the sympathetic and parasympathetic branches of the autonomic nervous system. Physiol Meas. 2009; 30: 215–226. pmid:19179746

38. Moraes ER, Murta LO, Baffa O, Wakai RT, Comani S. Linear and nonlinear measures of fetal heart rate patterns evaluated on very short fetal magnetocardiograms. Physiol Meas. 2012; 33: 1563–83. pmid:22945491

39. Vairavan S, Ulusar UD, Eswaran H, Preissl H, Wilson JD, Mckelvey SS et al. A computer-aided approach to detect the fetal behavioral states using multi-sensor Magnetocardiographic recordings. Comput Biol Med. 2016; 69: 44–51. pmid:26717240

40. Pardey J, Moulden M, Redman CW. A computer system for the numerical analysis of nonstress tests. Am J Obstet Gynecol. 2002; 186: 1095–103. pmid:12015543

41. Serra V, Bellver J, Moulden M, Redman CW. Computerized analysis of normal fetal heart rate pattern throughout gestation. Ultrasound Obstet Gynecol. 2009; 34: 74–9. pmid:19489020

42. Serra V, Moulden M, Bellver J, Redman CW. The value of the short-term fetal heart rate variation for timing the delivery of growth-retarded fetuses. Br Journal Obstet Gynecol. 2008; 115: 1101–7.

43. Kariniemi V, Ahopelto J, Karp PJ, Katila TE. The fetal magnetocardiogram. J Perinat Med. 1974; 2: 214–216. pmid:4468305

44. Stinstra JG, Peters MJ. The influence of fetoabdominal tissues on fetal ECGs and MCGs. Arch Physiol Biochem. 2002; 10: 165–176.

45. Andreotti F, Riedl M, Himmelsbach T, Wedekind D, Wessel N, Stepan H et al. Robust fetal ECG extraction and detection from abdominal leads Physiol Meas. 2014; 35: 1551–1567. pmid:25071095