Spectral analysis of antepartum fetal heart rate variability from fetal magnetocardiogram recordings

The magnetocardiogram

The magnetic field produced by the heart's electrical activity is called the magnetocardiogram (MCG). The first 20 years of MCG research established most of the concepts, instrumentation, and computational algorithms in the field. Additional insights into fundamental mechanisms of biomagnetism were gained by studying isolated hearts or even isolated pieces of cardiac tissue. Much effort has gone into calculating the MCG using computer models, including solving the inverse problem of deducing the bioelectric sources from biomagnetic measurements. Recently, most magnetocardiographic research has focused on clinical applications, driven in part by new technologies to measure weak biomagnetic fields.

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

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.

A computer-aided approach to detect the fetal behavioral states using multi-sensor Magnetocardiographic recordings

We propose a novel computational approach to automatically identify the fetal heart rate patterns (fHRPs), which are reflective of sleep/awake states. By combining these patterns with presence or absence of movements, a fetal behavioral state (fBS) was determined. The expert scores were used as the gold standard and objective thresholds for the detection procedure were obtained using Receiver Operating Characteristics (ROC) analysis. To assess the performance, intraclass correlation was computed between the proposed approach and the mutually agreed expert scores. The detected fHRPs were then associated to their corresponding fBS based on the fetal movement obtained from fetal magnetocardiogaphic (fMCG) signals. This approach may aid clinicians in objectively assessing the fBS and monitoring fetal wellbeing.

Differences in the sleep states of IUGR and low-risk fetuses: An MCG study

BACKGROUND

Intrauterine growth restriction (IUGR) is a fetal condition characterized by growth-rate reduction. Afflicted fetuses tend to display abnormalities in heart rate.

OBJECTIVE

To study the differences in the heart-rate variability of low-risk fetuses and IUGR fetuses during different behavioral states.

METHODS

A total of 40 fetal magnetocardiograms were analyzed from 20 low-risk and 20 IUGR fetuses recorded using a 151-sensor SQUID-array system. The maternal cardiac signals were attenuated using signal-space projection. Fetal R waves were identified using an adaptive Hilbert transform approach and fetal heart rate was calculated. In each three-minute window, the heart rate was classified into patterns reflective of quiet sleep (pattern A) and active sleep (pattern B) using the criteria of Nijhuis. Two adjacent 3-min windows exhibiting the same pattern were selected for analysis from every dataset. Heart-rate variability in that 6-min window was characterized using three measures, standard deviation of normal to normal (SDNN), root mean square of successive differences (RMSSD) and phase plane area (PPA).

RESULTS

All three measures tended to be lower in the IUGR group compared to the low-risk group. However, when the measures were analyzed in patterns, only PPA showed significant difference between the risk groups in pattern A, whereas both PPA and SDNN showed highly significant risk-group differences in pattern B. RMSSD did not show any significant risk-group difference.

CONCLUSION

The result signifies that the heart-rate variability of IUGR fetuses is different from that of low-risk fetuses, and only PPA was able to capture the HRV differences in both quiet and active states. The difference between these two groups of fetuses shows that the fetal-activity states are potential confounders when characterizing heart-rate variability.

Effects of fetal respiratory movements on the short-term fractal properties of heart rate variability

We evaluated the effect of fetal respiratory movements (RM) on the heart rate (HR) fractal dynamics.Abdominal ECG recordings were collected from low-middle-risk pregnant woman at rest. Mean gestational age was 34.8 ± 3.7 weeks. Ultrasound images were simultaneously acquired determining if RM were exhibited by fetuses. 13 pairs of HR series were compared. Each pair included 5 min of data from the same fetus either during the manifestation of RM or when there was no persistent indication of them. Detrended fluctuation analysis was applied to these series for obtaining the scaling exponent α1. HR series were also assessed using the conventional parameters RMSSD and HF power.The main findings of this contribution were the lack of significant changes in the scaling exponent α1 of fetal HR fluctuations as a result of RM. By contrast, HF power and RMSSD did show significant changes associated with the manifestation of fetal RM (p < 0.001 and p < 0.05, respectively). Yet the scaling exponent was the only parameter showing a significant relationship with the particular frequency of fetal RM (r s  = 0.6, p < 0.03). Given the invariability of α1 regarding the manifestation of fetal RM, we consider that the HR short-term fractal properties are convenient for assessing the cardiovascular prenatal regulation.

Change of Spectral Analysis of Fetal Heart Rate During Clinical Hypnosis: a Prospective Randomised Trial from the 20th Week of Gestation Till Term

Objective: To investigate the functional adaptive process of the fetal autonomic nervous system during hypnosis from the 20th week of gestation till term. Are there changes in the power spectrum analysis of fetal heart rate when the mother is having a clinical hypnosis or control period? Study Design: Fourty-nine FHR recordings were analysed. Included recordings were from singletons and abdominal fetal ECG-monitored pregnancies. All women were randomised to receive clinical hypnosis followed by a period with no intervention or vice versa. Statistical analyses were performed with the Wilcoxon signed ranks and Spearman rho correlation tests. Results: There was a significant difference found between fetal heart rate at baseline (144.3 ± 6.0) and hypnosis (142.1 ± 6.4). A difference was also detected between the standard deviation of the heart rate between baseline (6.7 ± 1.9) and hypnosis (6.8 ± 3.5). LFnu was smaller during baseline (80.2 ± 5.3) than during hypnosis (82.1 ± 5.7), whereas HFnu was significantly larger (19.8 ± 5.3 vs. 17.9 ± 5.7). There was no correlation between the gestation age and the change in LFnu, HFnu or ratio LF/HF due to the hypnosis intervention. Conclusion: The functional adaptive process of the fetal autonomic system during hypnosis is reflected by a sympathovagal shift towards increased sympathetic modulation.

Characterization of the fetal diaphragmatic magnetomyogram and the effect of breathing movements on cardiac metrics of rate and variability

Breathing movements are one of the earliest fetal motor behaviors to emerge and are a hallmark of fetal well-being. Fetal respiratory sinus arrhythmia (RSA) has been documented but efforts to quantify the influence of breathing on heart rate (HR) and heart rate variability (HRV) are difficult due to the episodic nature of fetal breathing activity. We used a dedicated fetal biomagnetometer to acquire the magnetocardiogram (MCG) between 36 and 38 weeks gestational age (GA). We identified and characterized a waveform observed in the raw data and independent component decomposition that we attribute to fetal diaphragmatic movements during breathing episodes. RSA and increased high frequency power in a time-frequency analysis of the IBI time-series was observed during fetal breathing periods. Using the diaphragmatic magnetomyogram (dMMG) as a marker, we compared time and frequency domain metrics of heart rate and heart rate variability between breathing and non-breathing epochs. Fetal breathing activity resulted in significantly lower HR, increased high frequency power, greater sympathovagal balance, increased short-term HRV and greater parasympathetic input relative to non-breathing episodes confirming the specificity of fetal breathing movements on parasympathetic cardiac influence. No significant differences between breathing and non-breathing epochs were found in two metrics reflecting total HRV or very low, low and intermediate frequency bands. Using the fetal dMMG as a marker, biomagnetometry can help to elucidate the electrophysiologic mechanisms associated with diaphragmatic motor function and may be used to study the longitudinal development of human fetal cardiac autonomic control and breathing activity.

Automatic identification of fetal breathing movements in fetal RR interval time series

Fetal breathing movements are associated with respiratory sinus arrhythmia (RSA). We present an algorithm which processes RR interval time series in the time and frequency domain, identifying spectral peaks with characteristics consistent with fetal RSA. Tested on 50 data sets from the second and third trimester, the algorithm had a sensitivity of 96.1%, false positive rate 35.7%, false negative rate 3.9%. The characteristics of automatically and visually identified episodes were very similar and corresponded the expected changes over gestation. The method is suited for easy and reliable identification of fetal breathing movements.

Aerobic exercise during pregnancy influences fetal cardiac autonomic control of heart rate and heart rate variability

BACKGROUND

Previous studies using ultrasound technology showed that fetal heart rate (HR) may be responsive to maternal aerobic exercise. Although it is recognized that cardiac autonomic control may be influenced by the intrauterine environment, little is known about how maternal exercise affects fetal heart development.

AIMS

This study tested the hypothesis that regular maternal exercise throughout gestation influences fetal cardiac autonomic control of HR and heart rate variability (HRV) when compared to fetuses of non-exercising women.

STUDY DESIGN

Magnetocardiograms (MCGs) were recorded using a dedicated fetal biomagnetometer at 28, 32 and 36 weeks gestational age (GA) from 26 regularly exercising (>30 min of aerobic exercise, 3x per week) and 35 healthy, non-exercising pregnant women. Fetal MCG was isolated and normal R-peaks were marked to derive fetal HR and HRV in the time and frequency domains. We applied a mixed-effects model to investigate the effects of exercise, GA and fetal activity state.

RESULTS

At 36 weeks GA, during the active fetal state, fetal HR was significantly lower in the exercise group (p=<0.0006). Post-hoc comparisons showed significantly increased HRV in the exercise group during the active fetal state at 36 weeks GA for both time and frequency domain measures.

CONCLUSION

These results indicate that regular maternal exercise throughout gestation results in significantly lower fetal HR and increased HRV.

Fetal heart rate variability reveals differential dynamics in the intrauterine development of the sympathetic and parasympathetic branches of the autonomic nervous system

The aim of this study was to investigate the hypothesis that fetal beat-to-beat heart rate variability (fHRV) displays the different time scales of sympatho-vagal development prior to and after 32 weeks of gestation (wks GA). Ninety-two magnetocardiograms of singletons with normal courses of pregnancy between 24 + 1 and 41 + 6 wks GA were studied. Heart rate patterns were either quiet/non-accelerative (fHRP I) or active/accelerative (fHRP II) and recording quality sufficient for fHRV. The sample was divided into the GA groups <32 wks GA/>32 wks GA. Linear parameters of fHRV were calculated: mean heart rate (mHR), SDNN and RMSSD of normal-to-normal interbeat intervals, power in the low (0.04-0.15 Hz) and high frequency range (0.15-0.4 Hz) and the ratios SDNN/RMSSD and LF/HF as markers for sympatho-vagal balance. fHRP I is characterized by decreasing SDNN/RMSSD, LF/HF and mHR. The decrease is more pronounced <32 wks GA. Beyond that GA SDNN/RMSSD is predominantly determined by RMSSD during fHRP I and by SDNN during fHRP II. In contrast to fHRP I, during fHRP II, mHR is positively correlated to SDNN/RMSSD instead of SDNN >32 wks GA. LF/HF increases in fHRP II during the first half of the third trimester. Non-accelerative fHRP are indicative of parasympathetic dominance >32 wks GA. In contrast, the sympathetic accentuation during accelerative fHRP is displayed in the interrelations between mHR, SDNN and SDNN/RMSSD. Prior to 32 wks GA, fHRV reveals the increasing activity of the respective branches of the autonomic nervous system differentiating the types of fHRP.

Fetal heart rate variability in growth restricted fetuses

Intrauterine growth restriction (IUGR) remains a major problem in perinatal medicine because of the variety of its underlying causes and the prediction of its outcome. Characteristics of heartbeat interval patterns are associated with neuro-vegetative and humoral regulatory processes. Fetal magnetocardiography allows non-invasive assessment of these processes with high precision throughout the second half of gestation. The aim of our study was the analysis of linear and non-linear parameters of fetal heart rate fluctuations to distinguish between IUGR fetuses and a cohort of normal subjects, both pre-selected from heart-rate traces representing a quiet state of activity in the third trimester of gestation.

Permutation entropy improves fetal behavioural state classification based on heart rate analysis from biomagnetic recordings in near term fetuses

The relevance of the complexity of fetal heart rate fluctuations with regard to the classification of fetal behavioural states has not been satisfyingly clarified so far. Because of the short behavioural states, the permutation entropy provides an advantageous complexity estimation leading to the Kullback-Leibler entropy (KLE). We test the hypothesis that parameters derived from KLE can improve the classification of fetal behaviour states based on classical heart rate fluctuation parameters (SDNN, RMSSD, ln(LF), ln(HF)). From measured heartbeat sequences (35 healthy fetuses at a gestational age between 35 and 40 completed weeks) representative intervals of 256 heartbeats were visually preclassified into fetal behavioural states. Employing discriminant analysis to separate the states 1F, 2F and 4F, the best classification result by classical parameters was 80.0% (SDNN). After additionally considering KLE parameters it was improved significantly (p<0.0005) to 94.3% (ln(LF), KLE_Mean). It could be confirmed that KLE can improve the state classification. This might reflect the consideration of different physiological aspects by classical and complexity measures.

Fetal source extraction from magnetocardiographic recordings by dependent component analysis

Fetal magnetocardiography (fMCG) has been extensively reported in the literature as a non-invasive, prenatal technique that can be used to monitor various functions of the fetal heart. However, fMCG signals often have low signal-to-noise ratio (SNR) and are contaminated by strong interference from the mother's magnetocardiogram signal. A promising, efficient tool for extracting signals, even under low SNR conditions, is blind source separation (BSS), or independent component analysis (ICA). Herein we propose an algorithm based on a variation of ICA, where the signal of interest is extracted using a time delay obtained from an autocorrelation analysis. We model the system using autoregression, and identify the signal component of interest from the poles of the autocorrelation function. We show that the method is effective in removing the maternal signal, and is computationally efficient. We also compare our results to more established ICA methods, such as FastICA.

Changes in the frequency power spectrum of fetal heart rate in the course of pregnancy

OBJECTIVE

The aim of this study was to examine changes in the heart rate variability based on the frequency power spectrum of healthy fetuses during the second and third trimester of pregnancy.

METHODS

We analyzed 222 fetal magnetocardiograms recorded in 49 healthy singleton pregnancies between the 16th and 42nd week. Discrete Fourier transformation was performed on the time-based step function of the RR-intervals. Changes of spectral density in the frequency spectrum in various bands between 0.003 to 1 Hz, including low-frequency (LF: 0.04-0.15 Hz) and high-frequency (HF: 0.15-0.40 Hz) bands, were examined as a function of gestational age.

RESULTS

Spectral density between 0.003 to 1.0 Hz increased with gestational age with large changes, in particular, at lower frequencies. At approximately the 32nd week, the rate of increase in power slowed substantially. Prior to this time, the rates of change in power were different for the bands 0.003 to 0.40 Hz, 0.40 to 0.60 Hz and 0.60 to 1.0 Hz. LF and HF showed similar development, with HF increasing slightly more rapidly.

CONCLUSION

We conclude that characteristic spectral bands that increase in spectral density at different rates during the second and third trimester may be identified. They most likely reflect developmental changes and behavioral states during pregnancy.

Simultaneity of foetal heart rate acceleration and foetal trunk movement determined by foetal magnetocardiogram actocardiography

We describe a new method of actocardiography, foetal magnetocardiogram (fMCG) actocardiography, which is based on the high sensitivity of the fMCG to foetal trunk movements. We demonstrate the efficacy of the method by applying it to assess the simultaneity of onset of foetal heart rate (FHR) acceleration and foetal trunk movement. The analysis was restricted to events for which the onset of FHR accelerations and foetal movements could be determined accurately, i.e. when FHR was stable and near the average quiescent level just prior to the acceleration. We found that FHR accelerations coincided with or preceded foetal movements nearly all the time. This supports the hypothesis of coordinated control of FHR accelerations and foetal movements more strongly than prior studies, based on other techniques. We also found that beat-to-beat FHR variability often decreased at or near the start of FHR accelerations and that this occurrence was an accurate marker of foetal movement onset, even when foetal movement onset lagged FHR accelerations.

Dynamic analysis of beat-to-beat fetal heart rate variability recorded by SQUID magnetometer: quantification of sympatho-vagal balance

BACKGROUND

Quantitative analysis of fetal heart rate variability (HRV) can be used to investigate the neural control mechanisms of fetal cardiac activity. However, conventional power spectrum methods do not reveal the full complexity of the time-varying sympatho-vagal balance in the fetus.

AIM

This study was carried out to explore alternative digital signal processing methods of analysing fetal HRV in time domain (rather than frequency domain), in line with most types of physiological monitoring.

METHODS

The beat-to-beat fetal heart rate was obtained by Superconducting Quantum Interference Device (SQUID) magnetocardiographic recording. These data were filtered within appropriately selected frequency bands: high frequency (HF) f>0.2 Hz, low frequency (LF) 0.05<f<0.2 Hz and very low frequency (VLF) f<0.05 Hz. To quantify the dynamics of sympathetic-parasympathetic interaction, the integrated amplitude of the HF curve (within a moving 4-s time window) was compared with the LF component. Also, the interaction of the VLF component with both HF and LF components was analysed.

RESULTS

A high degree of correlation was observed between these components for extended periods of time, although the presence of a strong correlation was found to depend on the fetal behavioural state. It was further observed that the low-frequency oscillations lag high-frequency activity by 1-3 s.

CONCLUSIONS

We propose a numerical parameter, based on the ratio of the LF to HF components of the fetal HRV as a quantitative measure of the relative gain of the parasympathetic system. Our results show this parameter to decrease with gestational age of the fetus.

Human fetuses have nonlinear cardiac dynamics

Approximate entropy (ApEn) is a statistic that quantifies regularity in time series data, and this parameter has several features that make it attractive for analyzing physiological systems. In this study, ApEn was used to detect nonlinearities in the heart rate (HR) patterns of 12 low-risk human fetuses between 38 and 40 wk of gestation. The fetal cardiac electrical signal was sampled at a rate of 1,024 Hz by using Ag-AgCl electrodes positioned across the mother's abdomen, and fetal R waves were extracted by using adaptive signal processing techniques. To test for nonlinearity, ApEn for the original HR time series was compared with ApEn for three dynamic models: temporally uncorrelated noise, linearly correlated noise, and linearly correlated noise with nonlinear distortion. Each model had the same mean and SD in HR as the original time series, and one model also preserved the Fourier power spectrum. We estimated that noise accounted for 17.2-44.5% of the total between-fetus variance in ApEn. Nevertheless, ApEn for the original time series data still differed significantly from ApEn for the three dynamic models for both group comparisons and individual fetuses. We concluded that the HR time series, in low-risk human fetuses, could not be modeled as temporally uncorrelated noise, linearly correlated noise, or static filtering of linearly correlated noise.

Noninvasive in utero assessment of PR and QRS intervals from the fetal magnetocardiogram

Fetal cardiac waveform intervals were assessed from fetal magnetocardiogram (FMCG) recordings taken from 59 pregnant women at 17-41 weeks gestation. Beyond 27 weeks' gestation PR and QRS intervals, measured from averaged waveforms, could be obtained from all subjects; however, prior to 21 weeks' gestation the success rate was 50% or less due to low signal-to-noise ratio. QT interval could not be assessed accurately in most subjects. Weak, but statistically significant, correlations with gestational age were found for PR interval (n = 145, R2 = 0.033, P = 0.028) and QRS interval (n = 145, R2 = 0.140, P < 0.0005). Abnormal waveform morphology was documented in several patients with cardiac malformations.

Heart rate dynamics in low risk human fetuses

Evaluation of nonlinear heart rate (HR) dynamics has received considerable attention in the pediatric literature because such analyses not only provide insight into underlying control mechanisms, but may also help to differentiate between normal and abnormal infants. The purpose of this study was to determine, in eight low risk human fetuses, if nonlinear HR dynamics could be identified by analyzing the dispersion of interbeat intervals at slow (Ds) and fast (Df) HRs. The fetal cardiac electrical signal was captured transabdominally at a resolution of +/- 1 ms. To test the null hypothesis, that the time series is the result of a linear stochastic process, Ds and Df for the original time series were compared with the values calculated for three linear models. The linear models were constructed to preserve the major statistical properties of the original time series, including the mean, SD, and the Fourier power spectrum. For each fetus, there was no evidence of nonlinear cardiac dynamics based on analyses of Ds and Df. In contrast, the distribution of adjacent R-R intervals and the pattern of change across three successive interbeat intervals both revealed significant nonlinearities in HR control in each fetus. If the difference between normal and abnormal infants is the result of aberrant control of nonlinear processes, then our findings indicate that parameters which describe the nonlinearity may be more useful then Ds and Df in assigning a risk status.

[Fetal magnetocardiography]

Fetal magnetocardiography is a new, alternative method for prenatal surveillance. The fetal magnetocardiogram (FMCG) registers the magnetic field produced by conduction currents in the fetal heart. Compared to the fetal electrocardiogram, the propagation of magnetic fields is relatively undisturbed by surrounding tissue. The FMCG thus has the advantage of a higher signal-to-noise ratio and can be acquired earlier pregnancy. Also, the high temporal resolution of the signal permits a significantly more precise determination of fetal heart rate parameters than fetal ultrasound. FMCG registration using a biomagnetometer is noninvasive and can be performed as of the second trimeter. It can be used to examine signal morphology, cardiac time intervals, heart rate variability as well as cardiac magnetic fields. To date, arrhythmic activity has been observed in the form of supraventricular and ventricular ectopies as well as atrial flutter, atrio-ventricular block, atrial tachycardia and Torsades de Pointes tachycardia. We also report here on the presence of short episodes of bradycardia in the second trimester of normal pregnancy. Measurement of the magnetic field strength at various locations above the abdomen has allowed the reconstruction of the fetal cardiac magnetic field and the determination of its relation to the position of the fetus. Signal averaging has permitted the precise examination of signal amplitude and cardiac time intervals and has shown that they increase in the course of pregnancy. Heart rate variability could be quantified in the time and frequency domain as well as using parameters of nonlinear dynamics. The results demonstrated an increase of variability and complexity over gestational age. Furthermore spectral analysis of fetal heart arte data could be associated with sympathetic and parasympathetic activity as well as, with respiration. Although the studies presenting these results have involved only limited numbers of observations, they demonstrate the potential of the method in the examination of the fetal conductive system, arrhythmias, congential defects, growth, development of the autonomic system, acidosis and distress. Furthermore, first results in pathological cases indicate that it may become a valuable tool in prenatal diagnostics. Improvements in instrumentation as well as prospective multicenter studies with larger numbers of appropriate subjects are required to determine whether magnetocardiography will establish itself as a new tool in clinical fetal, surveillance.

Foetal magnetocardiogram amplitude oscillations associated with respiratory sinus arrhythmia

We show that oscillations at foetal breathing frequencies observed in foetal heart rate (FHR) tracings obtained from foetal magnetocardiogram (FMCG) recordings are often accompanied by synchronous modulation of FMCG signal amplitude. This implies an association between the FHR oscillations and foetal chest wall movements, corroborating the hypothesis that the oscillations are due to respiratory sinus arrhythmia.