Skip to main content

A hierarchical procedure to select intrauterine and extrauterine factors for methodological validation of preterm birth risk estimation

Abstract

Background

Etiopathogenesis of preterm birth (PTB) is multifactorial, with a universe of risk factors interplaying between the mother and the environment. It is of utmost importance to identify the most informative factors in order to estimate the degree of PTB risk and trace an individualized profile. The aims of the present study were: 1) to identify all acknowledged risk factors for PTB and to select the most informative ones for defining an accurate model of risk prediction; 2) to verify predictive accuracy of the model and 3) to identify group profiles according to the degree of PTB risk based on the most informative factors.

Methods

The Maternal Frailty Inventory (MaFra) was created based on a systematic review of the literature including 174 identified intrauterine (IU) and extrauterine (EU) factors. A sample of 111 pregnant women previously categorized in low or high risk for PTB below 37 weeks, according to ACOG guidelines, underwent the MaFra Inventory. First, univariate logistic regression enabled p-value ordering and the Akaike Information Criterion (AIC) selected the model including the most informative MaFra factors. Second, random forest classifier verified the overall predictive accuracy of the model. Third, fuzzy c-means clustering assigned group membership based on the most informative MaFra factors.

Results

The most informative and parsimonious model selected through AIC included Placenta Previa, Pregnancy Induced Hypertension, Antibiotics, Cervix Length, Physical Exercise, Fetal Growth, Maternal Anxiety, Preeclampsia, Antihypertensives. The random forest classifier including only the most informative IU and EU factors achieved an overall accuracy of 81.08% and an AUC of 0.8122. The cluster analysis identified three groups of typical pregnant women, profiled on the basis of the most informative IU and EU risk factors from a lower to a higher degree of PTB risk, which paralleled time of birth delivery.

Conclusions

This study establishes a generalized methodology for building-up an evidence-based holistic risk assessment for PTB to be used in clinical practice. Relevant and essential factors were selected and were able to provide an accurate estimation of degree of PTB risk based on the most informative constellation of IU and EU factors.

Peer Review reports

Background

Etiologic vs generalized prediction of PTB

Preterm Birth (PTB) affects 5 to 18% of pregnancies and it is the leading cause of neonatal death [1]. The major etiological dichotomy of PTB involves primarily iatrogenic (iPTB) or spontaneous (sPTB) components, arising generally from different pathophysiology. However, differentiation is often difficult and clear phenotypic and etiologic classification of PTB remains a controversial issue [2]. Etiopathogenesis of PTB is multifactorial and factors possibly contributing to but not completely explaining PTB include: obstetrical and gynecological history (i.e. prior PTBs, short cervical length), chronic medical conditions (i.e. arterial hypertension, diabetes mellitus) and pregnancy complications (i.e. preeclampsia, urinary tracts infections) [3]. A number of other risk factors for PTB have been globally identified and their relevance is broadly acknowledged including: a history of pregnancy complications [4], placental abnormalities or abruption [5], preeclampsia [6], cigarette smoking [7], infections [8] and fetal growth restriction [911]. However, little is known about the interplay of these components along with other environmental and social factors. Previous studies produced solid evidence for Bayesian PTB risk estimation based upon history, specific clinical features or biomarkers, for both iPTB or sPTB [1214]). However, these prediction models, despite remarkable performances are limited to the capacity of prediction within the mechanisms for which they were defined and tested. Given the etiologic diversity of PTB it has been referred to as the “great obstetrical syndrome” or “Preterm Birth Syndrome” [15], characterized by both uterine and extrauterine components. There is an urgent need of holistic-generalized prediction models of PTB, capable of encompassing all or most etiologic mechanisms of PTB, in order to provide a useful tool at the bedside for clinical decision purposes. Therefore, we describe here the method of definition of a multi-layer prediction model, suitable for assisting clinicians in planning patients monitoring or treatment(s). Quantification and assessment of risk factors is a clinically established method to characterize the risk of PTB as defined by ACOG [16, 17]. However, it is currently difficult to translate this information into a comprehensive risk assessment capable of attributing to different factors or clusters specific weights, thus guiding a patient-specific approach. This risk individualization would lead to definite managements and treatments within a precision medicine approach.

Study aims and objectives

The primary aim of this study was to identify all available risk factors for PTB and to select the most relevant among the overall group for the scope of defining accurate risk prediction. This primary aim was pursued first by quantifying the risk of PTB associated with intrauterine (IU) and extrauterine (EU) factors separately to evaluate the contribution of factors in each gestational dimension, and then by searching for the most informative risk factors and select the best model including both IU and EU factors with the greatest contribution for the degree of estimated PTB risk. The second purpose of this study was to verify the predictive accuracy of the model including solely the most informative IU and EU risk factors of PTB birth identified at the previous step, as compared to a model including all the IU and EU risk factors considered in the present study in order to pinpoint the PTB constellation of risk factors according to which clinicians may attribute the degree of estimated PTB risk and increase confidence in managing potential PTB. The ultimate objective of this study was to overcome the high/low risk dichotomy of PTB risk through a clustering approach in order to verify if our sample of pregnant women may distribute in more than 2 groups based on the most informative IU and EU risk factors. Such an approach may allow us to assign to each pregnant woman a group membership probability value based on relevant risk factors able to define potential group boundaries along a continuum that moves from a lower to a higher degree of PTB risk. Such a gestational constellation of both IU and EU factors defining PTB risk, may set the basis for future creation of a clinical assessment tool reliably evaluating and quantifying the patient-specific risk for PTB that could indicate in each specific clinical case the most effective clinical management including: monitoring, surveillance, prophylaxis or therapy.

Materials and methods

Participants

This study was based on data collected from pregnant women (n = 116) enrolled in the Obstetrics Gynecology Department of the San Raffaele Hospital in Milan, between May 2019 and December 2019. The patients who met the inclusion criteria were informed on the aim of the study, the Maternal Frailty Inventory (see below for details) administration and the time commitment necessary and were asked if they wanted to participate. The criteria for inclusion were: mothers with a single fetus, gestational age below 37 weeks. The exclusion criteria were: fetuses with prenatal diagnosis of major malformation with subsequent medical abortion, twin pregnancies, pregnant women who did not sign the informed consent. All pregnant women gave written informed consent before inclusion in the study and research was performed in compliance with the principles of the Helsinki Declaration, and was approved by the Ethics Committee of the San Raffaele Hospital, Milan. STROBE guidelines for cross-sectional studies were followed [18]. A final sample of 111 pregnant women, ranging in age from 16 to 45 years with a mean age of 32.5 years (standard deviation 6.05) was included in the present study.

Risk factors definitions

We classified the sample of 111 pregnant women in low (i.e. n=68) or high risk (i.e. n=43) for PTB, according to the 2012 ACOG clinical management guidelines [17], considering risk factors related to obstetrical and gynecological history and pregnancy complications. Each risk factor has been labeled in relation to the environment from which it originates: uterus (U), placenta (P), fetus (F), or extrauterine indicating those clinical conditions that do not originate from the utero-placenta-fetal system. Both mother age (p=.075) and gestational age (p=.71) did not significantly differ between low and high risk subsamples of pregnant women (see Table 1 and Fig. 1). PTB was defined as any delivery below 37 weeks and 0 days of gestational age, regardless of the etiology (spontaneous or iatrogenic). Dating of pregnancy was confirmed in all cases included in the study with ultrasound crown-rump length measurements obtained in the first trimester. Fetal growth restriction was defined in presence of fetal biometry below the 3rd centile for gestational age or, with a drop in abdominal circumference below the 10th centile or EFW of >2 quartiles or >50 percentiles and/or abnormalities of fetal or uterine Doppler [11]. Pregnancy induced hypertension was defined as systolic blood pressure ≥140 mmHg and/or diastolic blood pressure ≥90 mmHg developing after 20 weeks of gestation, whereas hypertension that antecedes pregnancy or was present on at least two occasions before the 20th week of gestation or persists longer than 12 weeks postpartum was defined as chronic hypertension. Preeclampsia was defined as a new onset of hypertension and proteinuria or a new onset of hypertension and significant end-organ dysfunction with or without proteinuria after 20 weeks of gestation or postpartum in a previously normotensive woman [23]. Short cervix was defined in presence of transvaginal ultrasound measurement cervical length below 15 mm at 20-24 weeks [28]. Polydramnios was defined as a deepest amniotic fluid above 8 cm, oligohydramnios as an amniotic fluid index below 5 cm. Placenta previa was defined when the lower edge of the placental insertion was distant from the internal os ≤15 mm at time of delivery, and abruption placentae as any hemorrhage above 20 weeks trimester resulting in the evidence of bleeding from the placental bed.

Fig. 1
figure 1

Risk factors distribution. We show the distribution of risk factors in the sample of pregnant women according to maternal age and fetal gestational age

Table 1 Risk Factors preterm birth (PTB). Risk factors for PTB used for sample classification in low/high PTB risk and their definitions

Maternal frailty inventory of risk factors (MaFra)

In order to identify a maternal frailty for PTB, we developed the Maternal Frailty (MaFra) Inventory, identifying risk factors as IU or EU [29]. IU factors encompass all aspects of the uterine environment (i.e. exclusively uterus and its contents including the utero-placenta-fetal system) whereas EU factors involve all other aspects arising from extrinsic factors or different anatomical regions (e.g. systemic infections, autoimmune syndromes, medi- cations, stress or anxiety, etc). A systematic search of PubMed, CINAHL, Google Scholar and PsycINFO was performed in accordance with a detailed search strategy in order to pinpoint all relevant IU and EU PTB risk factors. All published citations were cross-referenced from other relevant studies. Studies that investigated the effects or the association between maternal gynecological, obstetrical, environment, mental state, lifestyle and any other socio-demographical factor associated with PTB were eligible for inclusion. Two authors (MC and FT) independently reviewed all studies to identify eligible studies for each risk factor included in the MaFra inventory. We included studies published in English only, from January 1990 through September 2019. We included systematic reviews, meta-analysis, randomized controlled trials, cohort, case-control and cross-sectional studies in which PTB is the outcome of interest. The reference lists of review articles and relevant meta-analyses were checked for additional references. Only peer-reviewed literature was included. Studies from high-income, middle-income and developing countries were all included. In addition, reported epidemiological data by national and international organizations [3033] was accounted for risk factor inclusion. Finally, The MaFra Inventory comprises 174 variables encompassing all intrauterine and extra-uterine factors included and referenced in (Table 2).

Table 2 Maternal Frailty (MaFra) Inventory. Intrauterine and extrauterine preterm birth risk factors included in the MaFra Inventory. 150 items assess the 71 listed factors while 24 Items collect more general sociodemographic, anamnestic and basic pregnancy history information for a total of 174 items included in the MaFra Inventory
Table 2 Maternal Frailty (MaFra) Inventory. Intrauterine and extrauterine preterm birth risk factors included in the MaFra Inventory. 150 items assess the 71 listed factors while 24 Items collect more general sociodemographic, anamnestic and basic pregnancy history information for a total of 174 items included in the MaFra Inventory (Continued)

Data preprocessing

After the data collection process, data cleaning procedures implied that all the variables with a) more than 50% of missing values or b) less than 3 positive occurrences for any dichotomic categorical variable (e.g. no subject has declared use of bronchodilators or HIV infection) were not taken in consideration for further analyses steps with more than half subjects with missing values. Along the same lines, we collapsed nested categorical variables (i.e. that depend on a previous positive answer to another question) and aggregated levels whenever a) or b) applied. After data cleaning, a subset of 86 variables of the original 174 were retained for all analyses steps. The median non-response rate for each variable, considering also missing values, was 0.9% while the average is 6%. To constitute the final dataset, we use random forest imputation (see e.g. [89] for an overview) to avoid slicing effects (i.e. NAs in different places for different variables that inevitably reduce the sample size) at the modeling step of the explanatory phase. This procedure is commonly considered robust especially in a low/moderate non response framework as the one of our study.

Data analyses

Explanatory phase

We performed a univariate logistic regression (see ch. 6 in [90] for a detailed overview) for each variable and ordered the risk factors based upon the associated p-values, from the smallest to the biggest. To avoid the potential imputation bias on the ordering procedure, we make use of the unimputed Mafra dataset. In order to deal with perfect separation (see e.g. [91]) (i.e. all subjects showing this risk factor labelled as high PTB risk) we used a bias reduction technique [92] leading to increased accuracy of all the estimates. To select the model including the most informative MaFra factors, we use AIC (i.e. Akaike Information Criterion see [93]) which is a well-known model selection criterion (see e.g. [94] ch. 2 and 3 for a detailed discussion). Based on the p-value ordering, relevant factors were added in a model one at a time and for each successive logistic regression model an AIC value was calculated. The model search ended when the sequence of AIC values reached its minimum, since further addition of other risk factors to the model increased the AIC. Model selection was performed including a) IU factors only b) EU factors only and c) including both IU and EU factors. At this modeling stage, to avoid slicing effects, we use the imputed dataset as described in the preprocessing section. To summarize, the aim of the explanatory phase was to select the most informative risk factors, representative of the Mafra inventory. Thus, p-value ordering was coupled with a model selection criterion to stop the search in order to select a group of variables which best explain PTB risk in our sample. Finally, we analyzed two models, based on the distinction between IU and EU risk factors, to understand the specific contribution of the two gestational dimensions and a third model which combined both IU and EU risk factors to identify the most informative factors irrespective of gestational dimension. For this combined evaluation, AIC differences were also used to assess the level of empirical support of a given model with respect to the best model (i.e. the one that reaches minimum AIC). The i-th AIC difference, between the i-th model and the best model, was calculated as Δi=AICiAICmin. Starting from values of Δi>4, the level of empirical support for the i-th model, with respect to the best model is considerably less (see [94] ch. 2 for a detailed explanation).

Predictive phase

We tested the ability of the most informative IU and EU factors selected in the explanatory phase to summarize the MaFra Dataset in a classification task. In particular we compared the out of sample classification error of a random forest [95] trained with all the MaFra risk factors against a random forest trained only with the subset of risk factors determined by AIC. We estimated the classification error by means of leave-one-out-cross-validation (i.e. LOOCV see ch. 7 in [96] for a detailed explanation). In our case, it consisted of leaving one mother out as a test set and predicting her degree of PTB risk with random forest trained on the other 110 remaining mothers. We did this for every subject to obtain the desired estimate of the out of sample classification error for both competitors. The aim of this predictive phase was to understand if the subset of IU and EU factors, selected by AIC in the explanatory phase, reached a comparable or higher overall classification accuracy than the entire pool of both IU and EU risk factors included in the MaFra Inventory.

Clustering phase

We evaluated the degree of PTB risk with fuzzy C-means (i.e. FCM [97]) clustering of the most informative IU and EU factors selected in the explanatory phase and used for classification in the predictive phase. The main steps in the FCM algorithm included: (1) fix a number of clusters and assign coefficients randomly to each data point for being in one of the clusters; (2) compute the centroids (i.e. mean of all points, weighted by their degree of being into the cluster) and, for each data point, update the coefficients; (3) repetition until convergence.

Results

Explanatory phase

The ordered univariate p-values of the selected MaFra factors are reported in Table 3. For the model including only IU variables, AIC reached global minimum (i.e. 107.62) with the inclusion of Placenta Previa (PP), Pregnancy Induced Hypertension (PH), Cervix Length (CL), Fetal Growth Restriction (FG), Preeclampsia (PC) and Fetal Sex (FS) (see Fig. 2); for the model including only EU variables, AIC reached global minimum (i.e. 121.74) with the inclusion of Antibiotics (AB), Physical Exercise (PE), Maternal Anxiety (AX), Antihypertensives (AH), Depression Level (DL) and Hormones Medication (HO) (see Fig. 2); for the model including both IU and EU variables, AIC reached a local minimum (i.e. 104.47) including 9 risk factors: PP, PH, AB, CL, PE, FG, AX, PC, AH. On the other hand, a global minimum (i.e. 103.71) is reached with the addition of two more factors, namely DL and FS (see Fig. 3). We compared the two models in terms of AIC differences (Δi=AICiAICmin see [94] ch.2 for a detailed explanation). There was substantial empirical support (i.e. Δi<4) for the models with respectively 9 and 11 risk factors. There was also evidence (see Table 4) that adding the risk factors DL and FS does not bring additional explanatory power since the models are almost equivalent in terms of AIC. Thus we focused on the model with 9 risk factors, 5 IU and 4 EU, which is the most parsimonious. Since the role EU factors in explaining and predicting PTB risk is not clear-cut through the literature and thus requires validation and replicability from further studies, we also present a detailed overview of the 4 EU factors selected by AIC in the Supplementary material.

Fig. 2
figure 2

Akaike Information Criterion (AIC) Model Selection for intrauterine (IU) and extrauterine (EU) risk factors. Forward selection of both IU and EU most representative risk factors by AIC. On the x-axis we display the labels of the risk factors entered with respect to the order of Table 3; on the y-axis the AIC values. Each point on the graph represents the AIC value scored by a model composed by all the risk factors on the left of the point (including the label above the point). The green color represents a decrease in the expected preterm birth risk whenever the specific risk factor augments holding the other variables constant at a certain value (e.g. increasing the cervical length IU_CL leads to a lower expected risk of preterm birth). For dichotomous variables, the variation is from absence to presence (e.g. for EU_PE, doing physical exercise decrease the expected preterm birth risk). In the specific case of IU_FS (Fetal Sex), the variation is from male to female. On the other hand, the red color represents an increase in the expected preterm birth risk whenever a specific dichotomous risk factor becomes active holding the other variables constant at a certain value (e.g. for EU_AB, giving antibiotics leads to an expected increase of the preterm birth risk). Both colors directly reflect the signs of the beta coefficients associated to each risk factor in the specific logistic regression model. Finally the grey colored points represent other possible risk factors that were not added because the criterion had already reached its minimum

Fig. 3
figure 3

Akaike Information Criterion (AIC) Model Selection for all Maternal Frailty (MaFra) risk factors. Forward selection of the most representative risk factors by AIC in the MaFra dataset. On the x-axis we display the labels of the risk factors entered with respect to the order of Table 3; on the y-axis the AIC values. Each point on the graph represents the AIC value scored by a model composed by all the risk factors on the left of the point (including the label above the point). The green color represents a decrease in the expected preterm birth risk whenever the specific risk factor augments holding the other variables constant at a certain value (e.g. increasing the cervical length IU_CL leads to a lower expected risk of preterm birth). For dichotomous variables, the variation is from absence to presence (e.g. for EU_PE, doing physical exercise decrease the preterm birth risk). In the specific case of IU_FS (Fetal Sex), the variation is from male to female. On the other hand, the red color represents an increase in the expected preterm birth risk whenever a specific dichotomous risk factor becomes active holding the other variables constant at a certain value (e.g. for EU_AB, giving antibiotics leads to an increase of the expected preterm birth risk). Both colors directly reflect the signs of the beta coefficients associated to each risk factor in the specific logistic regression model. Finally the grey colored points represent other possible risk factors that were not added because the criterion had already reached its minimum

Table 3 Ordered univariate logistic regression p-values of the risk factors below the 0.05 threshold
Table 4 Akaike information criterion (AIC) differences among the candidate nested models presented in Fig. 3

Predictive phase

In the LOOCV analysis including only the constellation of IU and EU variables identified in the explanatory phase, the random forest classifier achieved a best accuracy of 81.08% and an area under the curve (AUC) of 0.8122. On the other hand, for the LOOCV analysis including the full MaFra Inventory, the random forest classifier yielded a best accuracy of 80.18% and an AUC of 0.7374. We show in Fig. 4 the ROC curves of both classifiers. LOOCV variable importance is summarized in a variable importance network. The size of the vertices are relative to their average permutation score in Random Forest and the edges are related to the correlation structure of these variables (see Fig. 5). In particular, correlation among variables were computed with the point biserial (between continuous and categorical variables) and phi correlation coefficients (between categorical variables) which, in our particular case, are both equivalent to the Pearson correlation coefficient. In addition, correlation coefficients were considered significant at a p<0.05 threshold and colored accordingly (see Fig. 5).

Fig. 4
figure 4

ROC curve comparison between Akaike Information Criterion (AIC) selected and all Maternal Frailty (MaFra) risk factors. We compare the ROC curves of a random forest [95] trained only on the subset of 9 risk factors determined by AIC in the explanatory phase against a random forest trained on all the MaFra risk factors. We estimate the classification error for each competitor by means of leave-one-out-cross-validation (i.e. LOOCV see ch. 7 in [96] for a detailed explanation). We present also the area under the curve specific to each competitor in different colors (i.e. orange for the 9 risk factors selected by AIC and blue for the whole set of MaFra risk factors)

Fig. 5
figure 5

Variable importance network for the risk factors selected by Akaike Information Criterion (AIC). We build a network from the variable importance ([96] ch.15) of a random forest trained on the subset of 9 risk factors determined by AIC in the explanatory phase. The size of each vertex is proportional to the average variable importance while each edge represents the Pearson correlation coefficient between the two risk factors. In addition, the thickness of each edge is proportional to the strength of the correlation and colored edges imply a significant (i.e. pvalue<0.05) relationship either negative (red edge) or positive (green edge)

Clustering phase

The selected subset of relevant IU or EU factors that model PTB risk identified three clusters of pregnant women (see Fig. 6). In Cluster 1 (green), the typical woman has a longer than average CL (3.36 cm; average CL: 3.01 cm), has not taken AB, exercises more, is not anxious and has been treated with AH. We refer to this group of women as the low PTB risk cluster. In Cluster 2 (purple) instead the typical woman has taken AB, exercises less and is more anxious than low PTB risk women, however with a longer than average CL (3.24 cm), thus we refer to them as the middle PTB risk cluster. In Cluster 3 (orange) the typical woman presents a shorter than average CL (1.37 cm) than both low and middle risk mothers and PP. Moreover, these women take more AB, are more anxious and exercise less than low PTB risk women. Therefore we refer to them as the high PTB risk cluster. The matrix of centroids derived from fuzzy C-means clustering on the risk factors selected by AIC in the explanatory phase (see Supplementary material) was used to identify means or proportions in each cluster and spot the informative factors which enable grouping interpretation. Furthermore, in order to validate the clustering procedure, we focused on a subset of hard clustered subjects (Cluster 1: n=6; Cluster 2: n=8; Cluster 3: n=9) with coefficients of cluster membership above 0.7, showing the characteristics outlined above for each PTB degree of risk cluster. Namely, we assessed the time of delivery at birth for these women in order to parallel PTB risk cluster membership with birth gestational week. Time delivery of birth differed on average for Clusters 1, 2 and 3 (see Table 5). Namely, gestational week at birth decreases on average with degree of PTB risk cluster membership (i.e. Cluster 3 (High PTB Risk) < Cluster 2 (Middle PTB risk) < Cluster 1 (Low PTB risk) (see Table 5 for all cluster descriptive values). Preterm birth (<37 weeks of gestation) for hard clustered subjects included in the cluster validation procedure was for the larger part spontaneous (sPTB) (n=7) rather than elective (iPTB) (n=2).

Fig. 6
figure 6

Fuzzy C-means clustering for the risk factors selected by Akaike Information Criterion (AIC) and cluster centroids description. We present three clusters obtained by fuzzy C-means method on the 9 most informative variables selected by AIC at the explanatory phase. We split mothers based on the high/low preterm birth risk classification and we show them in ascending order with respect to cluster 1 membership score (i.e. 0 = no membership; 1 = total membership). We also show the informative preterm birth risk factors for Cluster 1, 2 and 3. The centroids matrix is also provided as a Supplementary material

Table 5 Summary statistics of the gestational week (GW) related to a subset (Cluster 1: n=6; Cluster 2: n=8; Cluster 3: n=9) of hard clustered subjects (i.e. with membership score greater or equal than 0.7)

Discussion

Explanatory phase

First of all, this study pinpoints relevant clinical risk factors related to prematurity within a comprehensive list of available factors described in the Maternal Frailty inventory classified as IU or EU factors.

Intrauterine factors

Within the intrauterine constellation of risk factors, our model found that the elements with the highest probability to predict high-risk of PTB to be: placenta previa (PP), pregnancy induced hypertension (PH), cervix length (CL), fetal growth restriction (FG), preeclampsia (PC) and fetal sex (FS). In particular, if PC, PH and/or FG arise during pregnancy, there is an increased likelihood of PTB since all of these clinical conditions share the same placental etiopathogenetic substrate representing a continuum, inherently characterizing risks for PTB [98, 99]. A connection between FG and PTB has been demonstrated via two case-controlled studies, one conducted in Europe [100] and the other in the USA [101]. In the EUROPOP study, 23% of preterm infants were found to be below the 10th percentile of fetal growth standards [100]. Whereas, Bukowski et al. [101] found that 25% of the fetuses delivered at or before 34 weeks did not reach the 5th percentile, 33% did not reach the 10th percentile, and 60% reached the 25th percentile of their individual growth potential. The association between PTB and FG is well established and recent literature confirms this concept [102, 103]. Independent studies described cervical shortening as a strong risk factor for sPTB both in asymptomatic [28, 104107] and symptomatic pregnant woman [108]. Moreover, a clear continuous inverse relationship was identified between CL at mid trimester ultrasound assessment and the probability of sPTB, with different cut-offs points proposed as a threshold (from 15 and 25 mm). Some authors spotlight a strong correlation of the association PP-short cervix with PTB risk [24]. Stafford et al. [109] reported that women with PP and CL of 30 mm or less were more likely to deliver before 37 weeks than those with a longer cervix, 69% compared with 21% respectively. Male sex is a risk factor for PTB [43], with noted increased association with hormonal differences between the sexes, relatively higher birth weight, and infection-related pathways [110, 111]. A possible explanation may be the higher presence of inflammatory markers found in male placentas compared with female placentas [112, 113], possibly acting synergistically with placental corticotropin-releasing hormone (CRH) production leading to the initiation of labor.

Extrauterine factors

Antibiotics (AB), physical exercise (PE), anxiety (AX), antihypertensive (AH) administration, depression (DL), hormones (HO) administration emerged as relevant elements of the EU constellation of risk factors predicting high/low risk of PTB in clinical context. AB administration can be linked to the presence of an infection during pregnancy [114]. However, the ORACLE trial provided evidence that AB should not be routinely prescribed for women in spontaneous preterm labour without evidence of clinical infection [115]. In our sample, 43 women consumed AB during pregnancy however, 27 of them did not actually have an infection. Our hypothesis is that AB administration during pregnancy likely leads to alterations in the maternal microbiome and previous studies showed that AB treatment to treat maternal infection did not reduce the incidence of PTB [61, 116, 117]. Furthermore, there is increasing evidence that some vaginal microbiomes are associated with an increased risk for PTB [118]and it is feasible that the microbiome alteration caused by AB consumption disrupts the EU environment predisposing for PTB. Moreover, specific components of the feto-placental microbiome have been linked to particular pregnancy complications, some that are strongly associated with PTB: pPROM, FG and PC [114]. PE may instead reduce the risk of PTB by increasing placental perfusion. This decreased risk is presumed to reduce oxidative stress [119], inhibiting the production of maternal and fetal cortisol and reducing the placental production of CRH. This inhibition of CRH prevents onset of uterine contractions and labor. DL has been found to be tightly related to adverse birth outcomes. A meta-analysis of nine studies (with 5,540 women) indicated that the risk of PTB (2.41; 1.47–3.56) was higher in depressed mothers compared to mothers without depression [120]. In general, DL, stress, and AX are associated with an increase in hypothalamic corticotropin-releasing factor (CRF) release and plasma cortisol concentrations. The elevation of placental CRF could likely initiate uterine contractions and cervical ripening [83]. This mechanism in normal, low risk pregnancy may determine a risk increase as it was described. However, AX plays a major role in adherence to treatments and compliance in high-risk pregnancies and may contribute to improve the outcome favoring prevention of noxae of greater magnitude. The relationship between AH administration during pregnancy and PTB risk may be related to possible underlying maternal hypertensive disorders, either before or during pregnancy. Certain AH medication, commonly used during pregnancy, such as beta-blockers, may diminish placental blood flow due to the subsequent vasoconstriction of placental vessels leading to a sort of placental haemodynamic insufficiency consequently leading to potential FG and increasing PTB risk [121]. Since sPTB is likely the final common pathway of several pathogenic processes, a single intervention such as progesterone supplementation is unlikely to benefit all women at risk. This concept was also shown in a major randomized controlled trial for prevention of PTB in women with short cervix in which the beneficial effect of natural progesterone was different according to the degree of cervical shortening [122]. A logistic regression analysis demonstrated significant treatment-genotype interactions, which could result in either a beneficial or harmful treatment response [123]. Studies have also found that women with certain characteristics, such as vaginal bleeding, gonorrhea, or chlamydia in the current pregnancy, a late PTB in a past pregnancy or penultimate PTB, are less likely to have a significant risk reduction [124]. We speculate that, because high-risk PTB women were more likely to deliver preterm, they were also more likely to receive treatment to attempt to delay delivery, such as progesterone supplementation. Finally, medically assisted conception were recently related to an increased risk of PTB both of spontaneous and iatrogenic etiology [22, 125].

Predictive phase

This study highlighted that a model including only the most informative and relevant selected IU and EU factors is able to classify pregnant women according to the PTB risk dichotomy with a higher accuracy as compared to a model comprising the full list of factors somehow associated to PTB and assessed through the MaFra. A precise constellation of the most informative IU and EU factors significantly impacting on prediction of PTB degree of risk may help guiding clinicians through an otherwise very extensive universe of factors. On one side the IU dimension is a mapped risk territory in the sense that most IU factors are already considered for PTB risk evaluation, on the other hand the EU dimension is populated by many factors related to maternal environment or lifestyle which may only somehow be associated with PTB risk. However, it is extremely important to pinpoint only those EU factors which may considerably add relevant information for the quantification of the degree of risk evaluated only upon clinical groundings in order to increase confidence for attributing PTB risk. In our sample we identified a specific constellation of informative IU and EU factors, however such factor configuration could be on the one hand externally validated on a new sample of mothers or on the other hand alternative configurations may stand out for different populations of pregnant women, therefore we believe that a hierarchical procedure which “explains” the entire gestational dimension on the basis of relevant emerging factors in each dimension for specific populations and able to inform a flexible and dynamic “predictive” model of PTB risk, can guide clinicians through the immense universe of potential factors without losing the focus on PTB risk salient influencers. An increasing number of risk factors are thought to interact to cause a transition from uterine quiescence toward preterm labour and sPTB [29] as well as that from diagnosis to decision of elective iPTB. These risk factors have been found to interact in various ways in the definition of phenotype and etiology of PTB. These factors perhaps aggregate mimicking constellations, thus escalating the potential for a woman to experience PTB. This universe of risk factors, if successfully mapped, could be visualized as a multi-layer model covering the expanse of IU and EU risk factors. This model could provide a more accurate assessment of the pregnancy and evolve into a clinical tool used to calculate the patient-specific risk and their risk specific management or treatment.

Clustering phase

We have identified three clusters related to degree of PTB risk ranging from high to middle and low, however a probability of group membership is assigned to each pregnant woman in a continuous fashion from 0 to 1, meaning that some women can be almost uniquely assigned to one cluster (i.e. >0.7), however for other women cluster probability membership may be slightly spread among the different clusters. The clustering phase outlined the definition of a frailty profile as a combination of informative IU and EU factors encompassed within the gestational environment (see Fig. 6).

Furthermore, identified clusters mirror the targeting of the delivery time of birth (i.e. gestational week) according to the probability (i.e. membership scores) of falling within the low, medium and high portions of the PTB risk constellation. The high PTB risk group of women clustered around a combination of biological and environmental factors (i.e. IU and EU factors). The middle-risk group clustered around mainly environmental factors (i.e. solely EU factors) without the influence of a pre-existent clinical frailty. It therefore seems that extrauterine factors may have an exclusive weight for the definition of a middle degree of risk and may also have potential impact on childbirth evidenced by the fact that mothers belonging to the middle risk cluster give birth earlier. On the one hand, clusters are clearly characterized by the probability of coexistence of certain risk factors, on the other hand, the effect of the interplay between risk factors more rarely associated should not be underestimated. The specific magnitude of absolute risk for each specific patient will actually be determined by the exploration of these specific associations, which we are willing to assess in future larger studies. For instance, we may claim that short CL added to PP could be informative of a higher degree of risk of PTB as compared to short CL added to maternal AX, so on and so forth. There will be room for future research in exploring these complex constellations of risk factors through Bayesian models of risk predictions or supervised machine learning techniques, the former with an a priori background risk model, the latter with classifiers on a degree of PTB risk continuum. Thus, this study sets the groundings for exploring a more extended universe of risk factors, including for instance the links between abnormal cardiovascular or placental function and PTB and this concept may be followed by the recognition of novel risk factors or potential biomarkers for prediction of PTB. It is known that abnormal cardiovascular adaptation in pregnancy leads to increased risk of preeclampsia, fetal growth restriction and potentially iatrogenic PTB [126]. Recent evidences show that abnormally low uterine arteries pulsatility index is described in pregnancies after in vitro fertilization and frozen blastocyst transfer as well as oocyte donation [22, 125]. In IVF/ICSI group there is on average a higher maternal age and high-risk of both preeclampsia and PTB (spontaneous or iatrogenic). Therefore, novel Doppler observations may help in the definition of specific phenotypes and novel risk factors with potential for generalization to all patients. From a more clinical perspective, frailty profiles may allow to differentiate resources allocation for patients according to their risk composition of IU and EU factors in relation to number of prenatal visits, lab and instrumental investigations such as ultrasound. A subgroup at very low risk would require minimal monitoring and intervention and a very high risk, would require intensive monitoring and treatment. Between these two extremes, at a middle risk PTB risk point, there would be room for continuous patient-specific tailoring of monitoring timing, with a criterion of proportionality according to cluster membership probability values and frailty profiles. For instance, clinicians could appropriately and timely indicate the use of prophylactic betametasone and magnesium sulphate for cases with a specific probability of PTB risk in line with a precise frailty profile, avoiding these treatments for cases below predefined membership values and with a non-overlapping frailty profile. This opens room for the creation of patient-specific risk assessment of PTB and precision medicine approach in this field.

Clinical translation of the model and future research directions

Obstetrics is going towards early individualization of risks as well as targeted monitoring and interventions [127]. The model described allows longitudinal risk reassessment on repeated visits of individual patients in which the risk background of each new visit would be the posterior risk found in the previous one. Repeated assessment of patients from pre-gestational period to the first trimester may favor preventive strategies for primary prevention. Reassessment at second and third trimester would permit minimizing medicalization in the lower risk groups and increasing surveillance and preparing therapeutic intervention in higher risk groups, proportionally to the risk magnitude of PTB. Thus, we believe that clinical management may be determined by the output of our risk assessment with reproducible and exact quantification. This study opens room to future research on this topic with the aim of accurate patient-specific risk estimation of PTB. This method would also help creating homogeneous study groups on which testing specific prophylactic or therapeutic interventions.

Conclusions

This study identified with a formal and reproducible methodological framework the selection of relevant and essential IU and EU risk factors able to provide an accurate definition of PTB risk. The first two phase (i.e. explanatory and predictive) identify and cross-validate significant predictors of PTB risk in our sample among both intrauterine factors and extrauterine factors. The third phase (i.e. clustering) spotlights the presence of specific factors which increase (i.e. anxiety) or reduce (i.e. physical exercise) the degree of risk of PTB, suggesting that PTB risk dimensionality is more complex than a binary universe and should be observed through the lens of a preterm birth syndrome in terms of intrauterine-extrauterine interactions that may predispose more or less to preterm birth risk. Founding on the background of IU and EU risk factors, this study establishes a generalized methodology for building-up an evidence-based holistic risk assessment for PTB to be used in clinical practice. Our proposed approach would contribute to increased reproducibility, exact quantification of risk assessment and would help in delivering patient-specific interventions.

Availability of data and materials

The dataset generated and analysed during the current study is not publicly available due to ongoing data analysis beyond the scope of this study but is available from cesare.miglioli@unige.ch on reasonable request through a private repository on github. The Maternal Frailty Inventory is available from dellarosa.pasquale@hsr.it on reasonable request.

References

  1. Liu L, Johnson HL, Cousens S, Perin J, Scott S, Lawn JE, et al. Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since 2000. Lancet. 2012; 379(9832):2151–61.

    Article  PubMed  Google Scholar 

  2. Goldenberg RL, Gravett MG, Iams J, Papageorghiou AT, Waller SA, Kramer M, et al. The preterm birth syndrome: issues to consider in creating a classification system. Am J Obstet Gynecol. 2012; 206(2):113–8.

    Article  PubMed  Google Scholar 

  3. Committee on Practice Bulletins–Obstetrics T, et al. Practice bulletin no. 130: prediction and prevention of preterm birth. Obstet Gynecol. 2012; 120(4):964.

    Article  Google Scholar 

  4. Phillips C, Velji Z, Hanly C, Metcalfe A. Risk of recurrent spontaneous preterm birth: a systematic review and meta-analysis. BMJ Open. 2017; 7(6):e015402.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Downes KL, Grantz KL, Shenassa ED. Maternal, labor, delivery, and perinatal outcomes associated with placental abruption: a systematic review. Am J Perinatol. 2017; 34(10):935.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lisonkova S, Joseph K. Incidence of preeclampsia: risk factors and outcomes associated with early-versus late-onset disease. Am J Obstet Gynecol. 2013; 209(6):544–e1.

    Article  PubMed  Google Scholar 

  7. Lumley J, Chamberlain C, Dowswell T, Oliver S, Oakley L, Watson L. Interventions for promoting smoking cessation during pregnancy. Cochrane Database Syst Rev. 2009; 3:CD001055.

    Google Scholar 

  8. Wingert A, Pillay J, Sebastianski M, Gates M, Featherstone R, Shave K, et al. Asymptomatic bacteriuria in pregnancy: systematic reviews of screening and treatment effectiveness and patient preferences. BMJ Open. 2019; 9(3):e021347.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Gardosi JO. Prematurity and fetal growth restriction. Early Hum Dev. 2005; 81(1):43–9.

    Article  PubMed  Google Scholar 

  10. Figueras F, Gratacos E. Update on the diagnosis and classification of fetal growth restriction and proposal of a stage-based management protocol. Fetal Diagn Ther. 2014; 36(2):86–98.

    Article  PubMed  Google Scholar 

  11. Lees C, Stampalija T, Baschat A, da Silva Costa F, Ferrazzi E, Figueras F, et al. ISUOG Practice Guidelines: diagnosis and management of small-for-gestational-age fetus and fetal growth restriction. Ultrasound Obstet Gynecol. 2020; 56(2):298–312.

    Article  CAS  PubMed  Google Scholar 

  12. O’Gorman N, Wright D, Poon L, Rolnik DL, Syngelaki A, de Alvarado M, et al. Multicenter screening for pre-eclampsia by maternal factors and biomarkers at 11–13 weeks’ gestation: comparison with NICE guidelines and ACOG recommendations. Ultrasound Obstet Gynecol. 2017; 49(6):756–60.

    Article  PubMed  Google Scholar 

  13. Khalil A, Elkhouli M, Garcia-Mandujano R, Chiriac R, Nicolaides K. Maternal hemodynamics at 11–13 weeks of gestation and preterm birth. Ultrasound Obstet Gynecol. 2012; 40(1):35–9.

    Article  CAS  PubMed  Google Scholar 

  14. Greco E, Gupta R, Syngelaki A, Poon LC, Nicolaides KH. First-trimester screening for spontaneous preterm delivery with maternal characteristics and cervical length. Fetal Diagn Ther. 2012; 31(3):154–61.

    Article  PubMed  Google Scholar 

  15. Romero R, Espinoza J, Kusanovic JP, Gotsch F, Hassan S, Erez O, et al. The preterm parturition syndrome. BJOGInt J Obstet Gynaecol. 2006; 113:17–42.

    Article  CAS  Google Scholar 

  16. American College of Obstetricians Gynecologists. Assessment of risk factors for preterm birth: ACOG practice bulletin No. 31. Obstet Gynecol. 2001; 98(4):709–16.

    Article  Google Scholar 

  17. American College of Obstetricians and Gynecologists, et al. ACOG practice bulletin no. 127: Management of preterm labor. Obstet Gynecol. 2012; 119(6):1308–17.

    Article  Google Scholar 

  18. Vandenbroucke JP, Von Elm E, Altman DG, Gøtzsche PC, Mulrow CD, Pocock SJ, et al. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): explanation and elaboration. PLoS Med. 2007; 4(10):e297.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Barros-Silva J, Pedrosa AC, Matias A. Sonographic measurement of cervical length as a predictor of preterm delivery: a systematic review. J Perinat Med. 2014; 42(3):281–93.

    Article  PubMed  Google Scholar 

  20. Pri-Paz S, Khalek N, Fuchs K, Simpson L. Maximal amniotic fluid index as a prognostic factor in pregnancies complicated by polyhydramnios. Ultrasound Obstet Gynecol. 2012; 39(6):648–53.

    Article  CAS  PubMed  Google Scholar 

  21. Mercer BM. Preterm premature rupture of the membranes. Obstet Gynecol. 2003; 101(1):178–93.

    PubMed  Google Scholar 

  22. Cavoretto P, Candiani M, Giorgione V, Inversetti A, Abu-Saba M, Tiberio F, et al. Risk of spontaneous preterm birth in singleton pregnancies conceived after IVF/ICSI treatment: meta-analysis of cohort studies. Ultrasound Obstet Gynecol. 2018; 51(1):43–53.

    Article  CAS  PubMed  Google Scholar 

  23. American College of Obstetricians and Gynecologists, et al. Gestational hypertension and preeclampsia: ACOG Practice Bulletin, number 222. Obstet Gynecol. 2020; 135(6):e237–60.

    Article  Google Scholar 

  24. Fishman SG, Chasen ST, Maheshwari B. Risk factors for preterm delivery with placenta previa. J Perinat Med. 2012; 40(1):39–42.

    Article  Google Scholar 

  25. American College of Obstetricians and Gynecologists, et al. ACOG Practice Bulletin No. 204: fetal growth restriction. Obstet Gynecol. 2019; 133(2):e97–109.

    Article  Google Scholar 

  26. Kalinderi K, Delkos D, Kalinderis M, Athanasiadis A, Kalogiannidis I. Urinary tract infection during pregnancy: current concepts on a common multifaceted problem. J Obstet Gynaecol. 2018; 38(4):448–53.

    Article  PubMed  Google Scholar 

  27. Kolstad KD, Mayo JA, Chung L, Chaichian Y, Kelly VM, Druzin M, et al. Preterm birth phenotypes in women with autoimmune rheumatic diseases: a population-based cohort study. BJOG Int J Obstet Gynaecol. 2020; 127(1):70–8.

    Article  CAS  Google Scholar 

  28. Celik E, To M, Gajewska K, Smith G, Nicolaides K. Cervical length and obstetric history predict spontaneous preterm birth: development and validation of a model to provide individualized risk assessment. Ultrasound Obstet Gynecol. 2008; 31(5):549–54.

    Article  CAS  PubMed  Google Scholar 

  29. Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. The Lancet. 2008; 371(9606):75–84.

    Article  Google Scholar 

  30. FIGO Working Group on Good Clinical Practice in Maternal–Fetal Medicine, Di Renzo GC, Fonseca E, Gratacos E, Hassan S, Kurtser M, et al. Good clinical practice advice: Prediction of preterm labor and preterm premature rupture of membranes. Int J Gynecol Obstetr. 2019; 144(3):340–6.

    Article  Google Scholar 

  31. Poon LC, Shennan A, Hyett JA, Kapur A, Hadar E, Divakar H, et al. The International Federation of Gynecology and Obstetrics (FIGO) initiative on pre-eclampsia: A pragmatic guide for first-trimester screening and prevention. Int J Gynecol Obstet. 2019; 145:1–33.

    Article  Google Scholar 

  32. Stoll BJ, Hansen NI, Bell EF, Shankaran S, Laptook AR, Walsh MC, et al. Neonatal outcomes of extremely preterm infants from the NICHD Neonatal Research Network. Pediatrics. 2010; 126(3):443–56.

    Article  PubMed  Google Scholar 

  33. Raju TN, Mercer BM, Burchfield DJ, Joseph Jr GF. Periviable birth: executive summary of a joint workshop by the Eunice Kennedy Shriver National Institute of child health and human development, Society for Maternal-Fetal medicine, American Academy of pediatrics, and American College of obstetricians and Gynecologists. Am J Obstet Gynecol. 2014; 210(5):406–17.

    Article  PubMed  Google Scholar 

  34. Wagura P, Wasunna A, Laving A, Wamalwa D, et al. Prevalence and factors associated with preterm birth at kenyatta national hospital. BMC Pregnancy Childbirth. 2018; 18(1):1–8.

    Article  Google Scholar 

  35. Zhang Y, Zhou J, Ma Y, Liu L, Xia Q, Fan D, et al. Mode of delivery and preterm birth in subsequent births: A systematic review and meta-analysis. PloS ONE. 2019; 14(3):e0213784.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Molitoris J. Breast-feeding During Pregnancy and the Risk of Miscarriage. Perspect Sex Reprod Health. 2019; 51(3):153–63.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Rafael TJ, Hoffman MK, Leiby BE, Berghella V. Gestational age of previous twin preterm birth as a predictor for subsequent singleton preterm birth. Am J Obstet Gynecol. 2012; 206(2):156–e1.

    Article  PubMed  Google Scholar 

  38. Saccone G, Perriera L, Berghella V. Prior uterine evacuation of pregnancy as independent risk factor for preterm birth: a systematic review and metaanalysis. Am J Obstet Gynecol. 2016; 214(5):572–91.

    Article  PubMed  Google Scholar 

  39. Klatsky PC, Tran ND, Caughey AB, Fujimoto VY. Fibroids and reproductive outcomes: a systematic literature review from conception to delivery. Am J Obstet Gynecol. 2008; 198(4):357–66.

    Article  PubMed  Google Scholar 

  40. Granese R, Mantegna S, Mondello S, Amadore D, Imbesi G, Calagna G, et al. Preterm birth: incidence, risk factors and second trimester cervical length in a single center population. A two-year retrospective study. Eur Rev Med Pharmacol Sci. 2017; 21(19):4270–7.

    CAS  PubMed  Google Scholar 

  41. Goldenberg RL, Iams JD, Das A, Mercer BM, Meis PJ, Moawad AH, et al. The preterm prediction study: sequential cervical length and fetal fibronectin testing for the prediction of spontaneous preterm birth. Am J Obstet Gynecol. 2000; 182(3):636–43.

    Article  CAS  PubMed  Google Scholar 

  42. Hagberg H, Mallard C, Jacobsson B. Role of cytokines in preterm labour and brain injury. BJOG Int J Obstet Gynaecol. 2005; 112:16–8.

    Article  CAS  Google Scholar 

  43. Peelen MJ, Kazemier BM, Ravelli AC, De Groot CJ, Van Der Post JA, Mol BW, et al. Impact of fetal gender on the risk of preterm birth, a national cohort study. Acta Obstet Gynecol Scand. 2016; 95(9):1034–41.

    Article  PubMed  Google Scholar 

  44. Razak A, Florendo-Chin A, Banfield L, Wahab MA, McDonald S, Shah P, et al. Pregnancy-induced hypertension and neonatal outcomes: a systematic review and meta-analysis. J Perinatol. 2018; 38(1):46–53.

    Article  CAS  PubMed  Google Scholar 

  45. Sibai BM. Preeclampsia as a cause of preterm and late preterm (near-term) births. Semin Perinatol. 2006; 30(1):16–19.

    Article  PubMed  Google Scholar 

  46. Turienzo CF, Sandall J, Peacock JL. Models of antenatal care to reduce and prevent preterm birth: a systematic review and meta-analysis. BMJ Open. 2016; 6(1):e009044.

    Article  Google Scholar 

  47. McDonald SD, Han Z, Mulla S, Beyene J. Overweight and obesity in mothers and risk of preterm birth and low birth weight infants: systematic review and meta-analyses. BMJ. 2010; 341:c3428.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Rosenberg TJ, Garbers S, Lipkind H, Chiasson MA. Maternal obesity and diabetes as risk factors for adverse pregnancy outcomes: differences among 4 racial/ethnic groups. Am J Public Health. 2005; 95(9):1545–51.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Sheehan PM, Nankervis A, Araujo Júnior E, Da Silva Costa F. Maternal thyroid disease and preterm birth: systematic review and meta-analysis. J Clin Endocrinol Metab. 2015; 100(11):4325–31.

    Article  CAS  PubMed  Google Scholar 

  50. Dreier JW, Andersen AMN, Berg-Beckhoff G. Systematic review and meta-analyses: fever in pregnancy and health impacts in the offspring. Pediatrics. 2014; 133(3):e674–88.

    Article  PubMed  Google Scholar 

  51. Lao TT, Hui AS, Sahota DS. Prior pregnancy and antenatal rubella sero-negativity–evidence of persistent maternal immunologic alteration?Am J Reprod Immunol. 2017; 78(3):e12727.

    Article  Google Scholar 

  52. Romero R, Avila C, Brekus CA, Morotti R. The role of systemic and intrauterine infection in preterm parturition. Ann N Y Acad Sci. 1991; 622(1):355–75.

    Article  CAS  PubMed  Google Scholar 

  53. Rahimi R, Nikfar S, Abdollahi M. Meta-analysis finds use of inhaled corticosteroids during pregnancy safe: a systematic meta-analysis review. Hum Exp Toxicol. 2006; 25(8):447–52.

    Article  CAS  PubMed  Google Scholar 

  54. Sujan AC, Quinn PD, Rickert ME, Wiggs KK, Lichtenstein P, Larsson H, et al. Maternal prescribed opioid analgesic use during pregnancy and associations with adverse birth outcomes: A population-based study. PLoS Med. 2019; 16(12):e1002980.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Su CY, Lin HC, Cheng HC, Yen AMF, Chen YH, Kao S. Pregnancy outcomes of anti-hypertensives for women with chronic hypertension: a population-based study. PLoS ONE. 2013; 8(2):e53844.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Asker C, Wikner BN, Källén B. Use of antiemetic drugs during pregnancy in Sweden. Eur J Clin Pharmacol. 2005; 61(12):899–906.

    Article  PubMed  Google Scholar 

  57. Aldridge TD, Hartmann KE, Michels KA, Velez Edwards DR. First-trimester antihistamine exposure and risk of spontaneous abortion or preterm birth. Pharmacoepidemiol Drug Saf. 2014; 23(10):1043–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Nielsen GL, Sorensen HT, Larsen H, Pedersen L. Risk of adverse birth outcome and miscarriage in pregnant users of non-steroidal anti-inflammatory drugs: population based observational study and case-control study. BMJ. 2001; 322(7281):266–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Cetingoz E, Cam C, Sakallı M, Karateke A, Celik C, Sancak A. Progesterone effects on preterm birth in high-risk pregnancies: a randomized placebo-controlled trial. Arch Gynecol Obstet. 2011; 283(3):423–9.

    Article  CAS  PubMed  Google Scholar 

  60. Fell D, Platt R, Lanes A, Wilson K, Kaufman J, Basso O, et al. Fetal death and preterm birth associated with maternal influenza vaccination: systematic review. BJOG Int J Obstet Gynaecol. 2015; 122(1):17–26.

    Article  CAS  Google Scholar 

  61. Smaill FM, Vazquez JC. Antibiotics for asymptomatic bacteriuria in pregnancy. Cochrane Database Syst Rev. 2019; 11:CD000490.

    Google Scholar 

  62. Mantovani E, Filippini F, Bortolus R, Franchi M. Folic acid supplementation and preterm birth: results from observational studies. BioMed Res Int. 2014; 2014:481914.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Charlton BM, Mølgaard-Nielsen D, Svanström H, Wohlfahrt J, Pasternak B, Melbye M. Maternal use of oral contraceptives and risk of birth defects in Denmark: prospective, nationwide cohort study. BMJ. 2016; 352:h6712.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  64. Sibai BM, Caritis SN, Hauth JC, MacPherson C, VanDorsten JP, Klebanoff M, et al. Preterm delivery in women with pregestational diabetes mellitus or chronic hypertension relative to women with uncomplicated pregnancies. Am J Obstet Gynecol. 2000; 183(6):1520–4.

    Article  CAS  PubMed  Google Scholar 

  65. Shah S, Banh ET, Koury K, Bhatia G, Nandi R, Gulur P. Pain management in pregnancy: multimodal approaches. Pain Res Treat. 2015; 2015:987483.

    PubMed  PubMed Central  Google Scholar 

  66. Nykjaer C, Alwan NA, Greenwood DC, Simpson NA, Hay AW, White KL, et al. Maternal alcohol intake prior to and during pregnancy and risk of adverse birth outcomes: evidence from a British cohort. J Epidemiol Community Health. 2014; 68(6):542–9.

    Article  PubMed  Google Scholar 

  67. Ludlow JP, Evans SF, Hulse G. Obstetric and perinatal outcomes in pregnancies associated with illicit substance abuse. Aust N Z J Obstet Gynaecol. 2004; 44(4):302–6.

    Article  PubMed  Google Scholar 

  68. Kobayashi S, Sata F, Murata K, Saijo Y, Araki A, Miyashita C, et al. Dose-dependent associations between prenatal caffeine consumption and small for gestational age, preterm birth, and reduced birthweight in the Japan Environment and Children’s Study. Paediatr Perinat Epidemiol. 2019; 33(3):185–94.

    Article  PubMed  Google Scholar 

  69. Kramer MS, Lydon J, Séguin L, Goulet L, Kahn SR, McNamara H, et al. Stress pathways to spontaneous preterm birth: the role of stressors, psychological distress, and stress hormones. Am J Epidemiol. 2009; 169(11):1319–26.

    Article  PubMed  Google Scholar 

  70. Cnattingius S, Villamor E, Johansson S, Bonamy AKE, Persson M, Wikström AK, et al. Maternal obesity and risk of preterm delivery. Jama. 2013; 309(22):2362–70.

    Article  CAS  PubMed  Google Scholar 

  71. Fuchs F, Monet B, Ducruet T, Chaillet N, Audibert F. Effect of maternal age on the risk of preterm birth: A large cohort study. PloS ONE. 2018; 13(1):e0191002.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  72. Murphy DJ. Epidemiology and environmental factors in preterm labour. Best Pract Res Clin Obstet Gynaecol. 2007; 21(5):773–89.

    Article  PubMed  Google Scholar 

  73. Srinivasjois RM, Shah S, Shah PS, Knowledge Synthesis Group on Determinants of Preterm/LBW Births. Biracial couples and adverse birth outcomes: a systematic review and meta-analyses. Acta Obstet Gynecol Scand. 2012; 91(10):1134–46.

    Article  PubMed  Google Scholar 

  74. Vatten LJ, Skjærven R. Effects on pregnancy outcome of changing partner between first two births: prospective population study. BMJ. 2003; 327(7424):1138.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Bilgin A, Mendonca M, Wolke D. Preterm birth/low birth weight and markers reflective of wealth in adulthood: a meta-analysis. Pediatrics. 2018; 142(1):e20173625.

    Article  PubMed  Google Scholar 

  76. Mozurkewich EL, Luke B, Avni M, Wolf FM. Working conditions and adverse pregnancy outcome: a meta-analysis. Obstet Gynecol. 2000; 95(4):623–35.

    CAS  PubMed  Google Scholar 

  77. Li X, Huang S, Jiao A, Yang X, Yun J, Wang Y, et al. Association between ambient fine particulate matter and preterm birth or term low birth weight: an updated systematic review and meta-analysis. Environ Pollut. 2017; 227:596–605.

    Article  CAS  PubMed  Google Scholar 

  78. Ikehara S, Kimura T, Kakigano A, Sato T, Iso H, Group JECS, et al. Association between maternal alcohol consumption during pregnancy and risk of preterm delivery: the Japan Environment and Children’s Study. BJOG Int J Obstet Gynaecol. 2019; 126(12):1448–54.

    Article  CAS  Google Scholar 

  79. Greenwood DC, Thatcher NJ, Ye J, Garrard L, Keogh G, King LG, et al. Caffeine intake during pregnancy and adverse birth outcomes: a systematic review and dose-response meta-analysis. Eur J Epidemiol. 2014; 29(10):725–34.

    Article  CAS  PubMed  Google Scholar 

  80. Imdad A, Bhutta ZA. Maternal nutrition and birth outcomes: Effect of balanced protein-energy supplementation. Paediatr Perinat Epidemiol. 2012; 26:178–90.

    Article  PubMed  Google Scholar 

  81. Col-Araz N. Evaluation of factors affecting birth weight and preterm birth in southern Turkey. J Pak Med Assoc. 2013; 63(4):459–62.

    PubMed  Google Scholar 

  82. Di Mascio D, Magro-Malosso ER, Saccone G, Marhefka GD, Berghella V. Exercise during pregnancy in normal-weight women and risk of preterm birth: a systematic review and meta-analysis of randomized controlled trials. Am J Obstet Gynecol. 2016; 215(5):561–71.

    Article  PubMed  Google Scholar 

  83. Dayan J, Creveuil C, Herlicoviez M, Herbel C, Baranger E, Savoye C, et al. Role of anxiety and depression in the onset of spontaneous preterm labor. Am J Epidemiol. 2002; 155(4):293–301.

    Article  CAS  PubMed  Google Scholar 

  84. Staneva A, Bogossian F, Pritchard M, Wittkowski A. The effects of maternal depression, anxiety, and perceived stress during pregnancy on preterm birth: A systematic review. Women Birth. 2015; 28(3):179–93.

    Article  PubMed  Google Scholar 

  85. Dole N, Savitz DA, Hertz-Picciotto I, Siega-Riz AM, McMahon MJ, Buekens P. Maternal stress and preterm birth. Am J Epidemiol. 2003; 157(1):14–24.

    Article  CAS  PubMed  Google Scholar 

  86. Zhang JJ, Ma XX, Hao L, Liu LJ, Lv JC, Zhang H. A systematic review and meta-analysis of outcomes of pregnancy in CKD and CKD outcomes in pregnancy. Clin J Am Soc Nephrol. 2015; 10(11):1964–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Crider KS, Whitehead N, Buus RM. Genetic variation associated with preterm birth: a HuGE review. Genet Med. 2005; 7(9):593–604.

    Article  CAS  PubMed  Google Scholar 

  88. Short CES, Taylor GP. Antiretroviral therapy and preterm birth in HIV-infected women. Expert Rev Anti-Infect Ther. 2014; 12(3):293–306.

    Article  CAS  PubMed  Google Scholar 

  89. Tang F, Ishwaran H. Random forest missing data algorithms. Stat Anal Data Mining ASA Data Sci J. 2017; 10(6):363–77.

    Article  Google Scholar 

  90. Vittinghoff E, Glidden DV, Shiboski SC, McCulloch CE. Regression methods in biostatistics: linear, logistic, survival, and repeated measures models.. Springer Science & Business Media: 2011. p. 139–202.

  91. Heinze G. A comparative investigation of methods for logistic regression with separated or nearly separated data. Stat Med. 2006; 25(24):4216–26.

    Article  PubMed  Google Scholar 

  92. Kosmidis I, Firth D, et al. A generic algorithm for reducing bias in parametric estimation. Elec J Stat. 2010; 4:1097–112.

    Article  Google Scholar 

  93. Akaike H. Information theory and an extension of the maximum likelihood principle. In: Selected papers of hirotugu akaike. New York: Springer: 1998. p. 199–213.

    Google Scholar 

  94. Burnham KP, Anderson DR, Vol. 2. A practical information-theoretic approach. Model selection and multimodel inference, 2nd ed. New York: Springer; 2002.

    Google Scholar 

  95. Breiman L. Random forests. Mach Learn. 2001; 45(1):5–32.

    Article  Google Scholar 

  96. Hastie T, Tibshirani R, Friedman J. The elements of statistical learning: data mining, inference, and prediction.Springer Science & Business Media; 2009.

  97. Bezdek JC, Ehrlich R, Full W. FCM: The fuzzy c-means clustering algorithm. Comput Geosci. 1984; 10(2-3):191–203.

    Article  Google Scholar 

  98. Morgan TK. Role of the placenta in preterm birth: a review. Am J Perinatol. 2016; 33(03):258–66.

    Article  PubMed  Google Scholar 

  99. Woelkers D, Barton J, von Dadelszen P, Sibai B. [71-OR]: the revised 2013 ACOG definitions of hypertensive disorders of pregnancy significantly increase the diagnostic prevalence of preeclampsia. Pregnancy Hypertens Int J Women’s Cardiovasc Health. 2015; 5(1):38.

    Google Scholar 

  100. Zeitlin J, Ancel P, Saurel-Cubizolles M, Papiernik E. The relationship between intrauterine growth restriction and preterm delivery: an empirical approach using data from a European case-control study. BJOG Int J Obstet Gynaecol. 2000; 107(6):750–8.

    Article  CAS  Google Scholar 

  101. Bukowski R, Gahn D, Denning J, Saade G. Impairment of growth in fetuses destined to deliver preterm. Am J Obstet Gynecol. 2001; 185(2):463–7.

    Article  CAS  PubMed  Google Scholar 

  102. Kabiri D, Romero R, Gudicha DW, Hernandez-Andrade E, Pacora P, Benshalom-Tirosh N, et al. Prediction of adverse perinatal outcome by fetal biometry: comparison of customized and population-based standards. Ultrasound Obstet Gynecol. 2020; 55(2):177–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  103. Ganzevoort W, Thornton JG, Marlow N, Thilaganathan B, Arabin B, Prefumo F, et al. Comparative analysis of 2-year outcomes in GRIT and TRUFFLE trials. Ultrasound Obstet Gynecol. 2020; 55(1):68–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Rozenberg P, Gillet A, Ville Y. Transvaginal sonographic examination of the cervix in asymptomatic pregnant women: review of the literature. Ultrasound Obstet Gynecol. 2002; 19(3):302–11.

    Article  CAS  PubMed  Google Scholar 

  105. To MS, Fonseca EB, Molina FS, Cacho AM, Nicolaides KH. Maternal characteristics and cervical length in the prediction of spontaneous early preterm delivery in twins. Am J Obstet Gynecol. 2006; 194(5):1360–5.

    Article  PubMed  Google Scholar 

  106. Crane J, Hutchens D. Transvaginal sonographic measurement of cervical length to predict preterm birth in asymptomatic women at increased risk: a systematic review. Ultrasound Obstet Gynecol. 2008; 31(5):579–87.

    Article  CAS  PubMed  Google Scholar 

  107. Greco E, Lange A, Ushakov F, Calvo JR, Nicolaides KH. Prediction of spontaneous preterm delivery from endocervical length at 11 to 13 weeks. Prenat Diagn. 2011; 31(1):84–9.

    Article  PubMed  Google Scholar 

  108. Berghella V, Palacio M, Ness A, Alfirevic Z, Nicolaides K, Saccone G. Cervical length screening for prevention of preterm birth in singleton pregnancy with threatened preterm labor: systematic review and meta-analysis of randomized controlled trials using individual patient-level data. Ultrasound Obstet Gynecol. 2017; 49(3):322–9.

    Article  CAS  PubMed  Google Scholar 

  109. Stafford IA, Dashe JS, Shivvers SA, Alexander JM, McIntire DD, Leveno KJ. Ultrasonographic cervical length and risk of hemorrhage in pregnancies with placenta previa. Obstet Gynecol. 2010; 116(3):595–600.

    Article  PubMed  Google Scholar 

  110. Hall CS, Domhoff GW, Blick KA, Weesner KE. The dreams of college men and women in 1950 and 1980: A comparison of dream contents and sex differences. Sleep. 1982; 5(2):188–94.

    Article  CAS  PubMed  Google Scholar 

  111. Pararas M, Skevaki C, Kafetzis D. Preterm birth due to maternal infection: causative pathogens and modes of prevention. Eur J Clin Microbiol Infect Dis. 2006; 25(9):562–9.

    Article  CAS  PubMed  Google Scholar 

  112. Ghidini A, Salafia CM. Gender differences of placental dysfunction in severe prematurity. BJOG Int J Obstet Gynaecol. 2005; 112(2):140–4.

    Article  Google Scholar 

  113. Clifton V. Sex and the human placenta: mediating differential strategies of fetal growth and survival. Placenta. 2010; 31:S33–9.

    Article  PubMed  CAS  Google Scholar 

  114. Kuperman AA, Koren O. Antibiotic use during pregnancy: how bad is it?BMC Medicine. 2016; 14(1):91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  115. Kenyon S, Taylor D, Tarnow-Mordi WF, et al. Broad-spectrum antibiotics for spontaneous preterm labour: the ORACLE II randomised trial. Lancet. 2001; 357(9261):989–94.

    Article  CAS  PubMed  Google Scholar 

  116. Lamont RF, Nhan-Chang CL, Sobel JD, Workowski K, Conde-Agudelo A, Romero R. Treatment of abnormal vaginal flora in early pregnancy with clindamycin for the prevention of spontaneous preterm birth: a systematic review and metaanalysis. Am J Obstet Gynecol. 2011; 205(3):177–90.

    Article  PubMed  PubMed Central  Google Scholar 

  117. Schneeberger C, Geerlings SE, Middleton P, Crowther CA. Interventions for preventing recurrent urinary tract infection during pregnancy. Cochrane Database Syst Rev. 2012; 11:CD009279.

    PubMed  Google Scholar 

  118. Fettweis JM, Serrano MG, Brooks JP, Edwards DJ, Girerd PH, Parikh HI, et al. The vaginal microbiome and preterm birth. Nat Med. 2019; 25(6):1012–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  119. Armson AB. Preterm labor. Obstet Evid Based Guidelines. 2017; 3:213–25.

    Google Scholar 

  120. Fekadu Dadi A, Miller ER, Mwanri L. Antenatal depression and its association with adverse birth outcomes in low and middle-income countries: A systematic review and meta-analysis. PloS ONE. 2020; 15(1):e0227323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  121. Petersen KM, Jimenez-Solem E, Andersen JT, Petersen M, Brødbæk K, Køber L, et al. β-Blocker treatment during pregnancy and adverse pregnancy outcomes: a nationwide population-based cohort study. BMJ Open. 2012; 2(4):e001185.

    Article  PubMed  PubMed Central  Google Scholar 

  122. Fonseca E, Celik E, Parra M, Singh M, Nicolaides K. Fetal Medicine Foundation Second Trimester Screening G. Progesterone and the risk of preterm birth among women with a short cervix. N Engl J Med. 2007; 357(5):462–9.

    Article  CAS  PubMed  Google Scholar 

  123. Manuck TA, Lai Y, Meis PJ, Dombrowski MP, Sibai B, Spong CY, et al. Progesterone receptor polymorphisms and clinical response to 17-alpha-hydroxyprogesterone caproate. Am J Obstet Gynecol. 2011; 205(2):135–e1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  124. Manuck TA, Stoddard GJ, Fry RC, Esplin MS, Varner MW. Nonresponse to 17-alpha hydroxyprogesterone caproate for recurrent spontaneous preterm birth prevention: clinical prediction and generation of a risk scoring system. Am J Obstet Gynecol. 2016; 215(5):622–e1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  125. Cavoretto PI, Giorgione V, Sotiriadis A, Viganò P, Papaleo E, Galdini A, Gaeta G, Candiani M. IVF/ICSI treatment and the risk of iatrogenic preterm birth in singleton pregnancies: systematic review and meta-analysis of cohort studies. J Matern Fetal Neonatal Med. 2020; 17(4):1–10.

    Google Scholar 

  126. Thilaganathan B, Kalafat E. Cardiovascular system in preeclampsia and beyond. Hypertension. 2019; 73(3):522–31.

    Article  CAS  PubMed  Google Scholar 

  127. Nicolaides KH. Turning the pyramid of prenatal care. Fetal Diagn Ther. 2011; 29(3):183–96.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all the pregnant women for their participation and their motivation. They would also like to thank Dr. Mirko Pozzoni and Dr. Silvia Lina Spinillo at the San Raffaele Hospital for their contribution in subject recruitment and for their assistance in MaFra administration and scoring. Edoardo Vignotto acknowledges funding from the Swiss National Science Foundation (Doc.Mobility Grant 188229).

Funding

This work was supported by the Italian Ministry of Health’s “Ricerca Finalizzata 2016” (grant number RF-2016-02364081; Principal Investigator: Dr. Pasquale Anthony Della Rosa).

Author information

Authors and Affiliations

Authors

Contributions

Pasquale Anthony Della Rosa: Conceptualization, Methodology, Formal analysis, Visualization, Writing - Original Draft, Writing - Review and Editing.

Cesare Miglioli: Conceptualization, Methodology, Formal analysis, Visualization, Writing - Original Draft, Writing - Review and Editing.

Martina Caglioni: Writing - Original Draft, Writing - Review and Editing.

Francesca Tiberio: Writing - Original Draft, Writing - Review and Editing.

Kelsey H.H. Mosser: Writing - Review and Editing.

Edoardo Vignotto: Methodology.

Matteo Canini: Visualization, Writing - Review and Editing.

Cristina Baldoli: Investigation, Supervision.

Andrea Falini: Supervision, Project administration.

Massimo Candiani: Supervision, Project administration.

Paolo Cavoretto: Conceptualization, Writing - Original Draft, Writing - Review and Editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Paolo Cavoretto.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the San Raffaele Hospital, Milan (39-OSR/2018) and was conducted in compliance with the principles of the Helsinki Declaration. All pregnant women received documentation and guidance describing all aspects of the study (e.g. Maternal Frailty Inventory content, benefits, scope). All pregnant women included in this study (n=111) provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

Detailed overview of the extrauterine (EU) factors selected by Akaike Information Criterion (AIC) in the explanatory phase and the centroids matrix derived from fuzzy C-means clustering on the risk factors selected by AIC in the explanatory phase.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Della Rosa, P.A., Miglioli, C., Caglioni, M. et al. A hierarchical procedure to select intrauterine and extrauterine factors for methodological validation of preterm birth risk estimation. BMC Pregnancy Childbirth 21, 306 (2021). https://doi.org/10.1186/s12884-021-03654-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12884-021-03654-3

Keywords