This paper presents an observational prospective study carried out on a population of pregnant women monitored by means of the three major fetal surveillance tests reported in the Background section. Data were collected at the Department of Mother and Child at the University Hospital of Naples Federico II, in collaboration with the Polytechnic of Milan and the University of Pavia. The population consisted of a total of 152 pregnant women: 25 had a diagnosis of pregestational diabetes mellitus (PGDM), 61 were diagnosed with gestational diabetes mellitus (GDM) and 66 were controls.
All pregnant women were admitted, followed up and gave birth in the same Obstetric Unit. Pregnant women were enrolled starting from the third trimester of gestation until the end of gestation. All enrolled subjects read and signed informed consent for participation in the study.
Inclusion criteria for the entire population were defined as follows:
- singleton pregnancy;
- certain determination of pregnancy dating calculated from the first day of the last menstrual period and confirmed by ultrasound measurements according to the population nomograms (8);
- gestational age at or beyond 26 weeks of gestation.
- no known pathologies which are not related to diabetes.
Healthy Controls were selected among those who matched the previous three inclusion criteria by adding the condition that they were defined as euglycemic during the screening examinations, i.e., 2-hour oral glucose tolerance test (OGTT) and had normal fasting blood glucose levels for all three pregnancy trimesters.
GDMs were selected on the basis of the OGTT test generally performed at 24-28 weeks of gestation. In a few cases, it was conducted earlier (i.e., at 16-18 weeks of gestation) based on additional risk factors for gestational diabetes. The 2-hour diagnostic OGTT consisted of a three-step approach, with the determination of fasting, 1-hour and 2-hour glucose levels after administration of a 75 g oral glucose solution (9).
PGDM was in most cases a known condition that preceded pregnancy. In a few cases only, PGDM was diagnosed in the first trimester or early second trimester with standard diagnostic criteria, defined as hemoglobin A1C (HbA1C) of 6.5% or higher, fasting plasma glucose of 126 mg/dL or higher, or 2-hour glucose of 200 mg/dL or greater during a 75-g oral glucose tolerance test (10).
Fasting glucose levels at the first and second trimester and pre-pregnancy BMI were retrospectively obtained from hospital records. Population characteristics are summarized in Table 1.
Table 1 Maternal characteristics of the study population.
|
Controls
|
GDMs
|
PGDMs
|
N
|
66
|
61
|
25
|
Pre-pregnancy BMI (Kg/m2)
|
25.4±4.9
|
29.8±7.7
|
27.9±5.5
|
Age (years)
|
31.7±5.2
|
35.0±5.3
|
31.6±5.8
|
Number of pregnancies
|
2 (1-4)
|
2 (1-9)
|
2 (1-4)
|
Fasting glucose 1st trim
|
74 (69-82)
|
86 (76-92)
|
123 (114.7-143.7)
|
Fasting glucose 2nd trim
|
78 (72-85)
|
89 (82.5-99)
|
133 (122.5-154.5)
|
Fasting glucose 3rd trim*
|
83 (76-88)
|
90 (82.0-99.2)
|
142 (114.0-174.2)
|
1-h Plasma glucose
|
142 (132-153)
|
184.5 (169-201)
|
-
|
2-h Plasma glucose
|
123 (132-118)
|
141 (118.5-159)
|
-
|
urinary proteinuria*
|
0 (0%)
|
5 (8.2%)
|
0 (0%)
|
urinary glycosuria*
|
0 (0%)
|
5 (8.2%)
|
11 (44%)
|
gestational hypertension*
|
0 (0%)
|
11 (18%)
|
5 (20%)
|
chronic hypertension
|
0 (0%)
|
4 (6.6%)
|
1 (4%)
|
preeclampsia*
|
0 (0%)
|
3 (4.9%)
|
0 (0%)
|
diabetic retinopathy
|
0 (0%)
|
0 (0%)
|
1 (4%)
|
chronic renal failure
|
0 (0%)
|
0 (0%)
|
1 (4%)
|
Asterisks (*) indicate values or diagnoses obtained during the course of the study
2.1 Experimental Protocol
Fetal monitoring included fetal biometry and amniotic fluid evaluation, Doppler velocimetry of the UA and MCA and antepartum cCTG monitoring. Moreover, to assess maternal glycemic control, fasting blood glucose levels were measured during the third trimester of pregnancy in all groups.
2.1.1 Amniotic fluid
Concerning amniotic fluid evaluation, the deepest vertical amniotic fluid pocket (DVP) was calculated in all pregnancies through the identification and measurement of the largest vertical diameter within an amniotic fluid pocket. Values lower than 2 cm identify oligohydramnios, while values over 8 cm are considered polyhydramnios (11).
2.1.2 Doppler velocimetry
The antepartum Doppler velocimetry was performed using a Voluson E8 (General Electric Healthcare Technologies, Milwaukee, WI, USA) ultrasound machine equipped with a 2–8MHz transabdominal transducer and color Doppler imaging was used to identify the UA and MCA vessels. Doppler flow spectra were obtained from the UA at the midsection of the umbilical cord and the middle portion of the MCA, as described in (12). The indices considered for the analyses are UA and MCA Pulsatility indices (UA-PI, MCA-PI) and the analysis was performed at the 26th, 32nd, and 36th weeks of gestation.
2.1.3 Antepartum CTG
The antepartum computerized CTG (cCTG) monitoring was performed in a controlled clinical environment with the patient lying on an armchair. The cCTG records were obtained using Philips Avalon FM30, which is equipped with an ultrasound transducer and a transabdominal tocodynamometer. Recordings lasted at least 40 minutes.
In our study, the FHR signal was collected simultaneously employing two computerized systems, namely the Philips IntelliSpace Perinatal (www.philips.com/IntelliSpacePerinatal) and the 2CTG2 system (13).
Although in pregnancies complicated by PGDM CTG monitoring is considered appropriate starting from the 32nd week of gestation, it can be performed at earlier gestational ages in case of additional complications. In our study, we decided to start CTG monitoring from the 31st week of gestation with approximately weekly frequency. Recordings include the fetal heart rate (FHR) and uterine contractions signals, and the fetal movements series (FM), which is set to 1 when the mother presses a button to indicate a perceived movement and 0 when she does not. A collection of parameters, both linear and nonlinear, were calculated from the CTG traces stored by the 2CTG2 software.
They include the average of the FHR baseline (MeanFHR) calculated using the Mantel algorithm (14), the number of accelerations per hour (#Acc), and the Short Time Variation (STV) (15). All these are commonly used CTG time-domain parameters (6,15).
Moreover, the distribution of power in LF (0.03-0.15 Hz), MF (0.15- 0.5 Hz), and HF (0.5-1 Hz) was calculated and expressed as the percentage of the total signal power (5). These features were computed on windows of three minutes which were then averaged within the same tracing.
Other parameters which are not used in clinical practice were also examined to provide a more comprehensive analysis of the FHR. In particular, two features that estimate signal complexity were included: the Sample Entropy (SampEn), which was computed setting x=1, m=0.1, and r=1 (16), and the ternary Lempel-Ziv complexity (LZ3), which was computed as reported in (17) setting p to 0.01. To compute this parameter the FHR signal is re-codified to a trace composed of 3 symbols which represent signal increase, decrease or stationary state. LZ3 quantifies the emergence of new patterns in this trace. Values close to 1 indicate high irregularity of the signal while lower values indicate the presence of repetitive patterns. Both SampEn and LZ3 were computed over three-minute windows, which were then averaged within the same recording.
Finally, we applied the Phase Rectified Signal Averaging (PRSA) technique to highlight hidden oscillations in the signal (18). The value of L was set to 40 FHR samples, T to 1 and s to 2. From the PRSA curve we computed the Deceleration Capacity (DC), the Deceleration Reserve (DR) (19) and, by applying the CWT, we estimated the power in MF and HF at the anchor point (20). The values obtained (MFprsa and HFprsa) were then normalized by the total power of the PRSA signal.
Fetal movements (MovF%) were estimated by measuring the percentage of time mothers press the dedicated button over the entire CTG recording.
These measurements were done for routine follow-up in cases and only for study purposes in the controls.
2.1.4 Pregnancy outcomes
Neonatal data were collected for all newborns. They include week and modality of delivery, birthweight, Apgar score at 1 and 5 minutes, UA blood gas values at birth, and access to the neonatal intensive care unit (NICU). UA blood gas values collected are: pH, O2 and CO2 partial pressure (pO2, pCO2), cord arterial base excess (BECF), lactates and bicarbonate (HCO3).
2.2 Statistical Analysis
UA-PI and MCA-PI, which were assessed at three different time points during the last trimester, were compared among groups using repeated measure analysis of variance (ANOVA). In particular, we looked for between-subjects effects in a repeated measures model. When performing pairwise comparisons, p-values were corrected with Bonferroni correction. Concerning the occurrence of oligohydramnios and polyhydramnios, differences between groups were evaluated using Fisher’s exact test.
The comparison of cCTG parameters between the three groups was made through the non-parametric Kruskal-Wallis test since normality was not respected for several of the analyzed features. The test was applied independently for each group of weeks. Additionally, Linear Mixed Effect models (LME) were also employed to make use of all the available recordings and account for repeated measures. For each of the CTG features discussed, we fitted an LME model. The fixed independent variables are the group (i.e., Controls, GDMs, PGDMs) and gestational age, while the subject is modeled as the random factor to account for repeated measures. Before fitting the models, variables were log-transformed when necessary to ensure that they were approximately normally distributed, and all were normalized using the standard scaler.
Umbilical cord gas values were compared through 2-way ANOVA, considering the mode of delivery (spontaneous/cesarean) and the groups as factors. Post-hoc comparisons between groups were corrected using the Bonferroni procedure. The prevalence of Apgar scores less than 7, neonatal intensive care admission, small for gestational age and large for gestational age were compared across groups using Fisher's exact test.
Analyses were conducted in MATLAB R2022b (The Math Works, Inc.).