Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 30 May 2022
Sec. Family Medicine and Primary Care

Development and Validation of Multi-Stage Prediction Models for Pre-eclampsia: A Retrospective Cohort Study on Chinese Women

\nZeyu TangZeyu Tang1Yuelong JiYuelong Ji2Shuang ZhouShuang Zhou2Tao SuTao Su3Zhichao YuanZhichao Yuan2Na HanNa Han3Jinzhu Jia,
Jinzhu Jia1,4*Haijun Wang,
Haijun Wang2,5*
  • 1Department of Biostatistics, School of Public Health, Peking University, Beijing, China
  • 2Department of Maternal and Child Health, School of Public Health, Peking University, Beijing, China
  • 3Maternal and Child Health Care Hospital of Tongzhou District, Beijing, China
  • 4Center for Statistical Science, Peking University, Beijing, China
  • 5National Health Commission Key Laboratory of Reproductive Health, Beijing, China

Objective: The aim of this study is to develop multistage prediction models for pre-eclampsia (PE) covering almost the entire pregnancy period based on routine antenatal measurements and to propose a risk screening strategy.

Methods: This was a retrospective cohort study that included 20582 singleton pregnant women with the last menstruation between January 1, 2013 and December 31, 2019. Of the 20582 women, 717 (3.48%) developed pre-eclampsia, including 46 (0.22%) with early-onset pre-eclampsia and 119 (0.58%) preterm pre-eclampsia. We randomly divided the dataset into the training set (N = 15665), the testing set (N = 3917), and the validation set (N = 1000). Least Absolute Shrinkage And Selection Operator (LASSO) was used to do variable selection from demographic characteristics, blood pressure, blood routine examination and biochemical tests. Logistic regression was used to develop prediction models at eight periods: 5–10 weeks, 11–13 weeks, 14–18 weeks, 19–23 weeks, 24–27 weeks, 28–31 weeks, 32–35 weeks, and 36–39 weeks of gestation. We calculated the AUROC (Area Under the Receiver Operating Characteristic Curve) on the test set and validated the screening strategy on the validation set.

Results: We found that uric acid tested from 5–10 weeks of gestation, platelets tested at 18–23 and 24–31 weeks of gestation, and alkaline phosphatase tested at 28–31, 32–35 and 36–39 weeks of gestation can further improve the prediction performance of models. The AUROC of the optimal prediction models on the test set gradually increased from 0.71 at 5–10 weeks to 0.80 at 24–27 weeks, and then gradually increased to 0.95 at 36–39 weeks of gestation. At sensitivity level of 0.98, our screening strategy can identify about 94.8% of women who will develop pre-eclampsia and reduce about 40% of the healthy women to be screened by 28–31 weeks of pregnancy.

Conclusion: We developed multistage prediction models and a risk screening strategy, biomarkers of which were part of routine test items and did not need extra costs. The prediction window has been advanced to 5–10 weeks, which has allowed time for aspirin intervention and other means for PE high-risk groups.

Introduction

Pre-eclampsia is a pregnancy related syndrome defined as newly occurred hypertension at or after 20 weeks of gestation, accompanied by proteinuria or other organs damage (1). The incidence of pre-eclampsia worldwide is 0.2-9.2% (2). A study including 112,386 pregnant women in China showed that the incidence of pre-eclampsia in China was approximately 2.87% in 2011 (3). Every year, approximately 76 thousands women and half million infants died from pre-eclampsia worldwide. Pre-eclampsia can have adverse effects on pregnant women, for instance, causing damage to the liver and kidney systems (4). If left untreated, it can lead to pulmonary edema, eclampsia, brain damage, and even maternal death (57). Pregnant women and their children affected by pre-eclampsia are at increased risk for long-term cardiovascular and chronic diseases, including chronic hypertension, stroke, metabolic syndrome, and cognitive impairment (816).

In the first trimester, taking pharmacologic interventions (e.g., aspirin) for high-risk pregnant women can reduce the risk of early-onset and preterm pre-eclampsia (17, 18). It can reduce the incidence of adverse perinatal outcomes by intensive monitoring and selecting the appropriate time of delivery during the second or third trimester (19). Early identification of high-risk groups will help to take interventions in advance. Therefore, it is of great significance to develop risk prediction models for pre-eclampsia.

However, pre-eclampsia related prediction models have been developed mainly in developed countries (20). In recent years, there were some studies developing prediction models of pre-eclampsia based on Chinese population (2126). These studies were mainly based on the hospital electronic medical data system and carried out in eastern China (three in Shanghai and one in Tianjin), screening pregnant women in the first and second trimester of pregnancy. However, these studies primarily focused on specific high-risk groups or used expensive biomarkers beyond the scope of routine testing, which limits their generalizability. In addition, these studies did not consider pre-eclampsia subtypes (2326), possible bias caused by the process of variable selection (23, 24, 26), and insufficient number of outcome events (2123, 25).

In this study, based on 20582 pregnant women in Beijing, we aimed to develop multistage prediction models covering almost the entire pregnancy period by selecting valuable predictors from routine antenatal measurements, and a risk screening strategy based on the optimal models.

Methods

Study Population

This was a Peking University Retrospective Birth Cohort in Tongzhou based on the hospital information system, including singleton pregnant women having prenatal care, delivery and outcome records, with last menstruation between 1 January 2013 and 31 December 2019 and with delivery gestational weeks no less than 28 weeks in Tongzhou Maternal and Child Health Care Hospital of Beijing. We further selected pregnant women with the latest record of deliveries, and with at least one, two and two antenatal examination records in the first, second, and third trimester, respectively. We excluded pregnant women using assisted reproductive technology, having systemic lupus erythematosus, or chronic hypertension, or gestational hypertension without pre-eclampsia. Also, we excluded pregnant women who lacked blood pressure measurements, blood routine examination and biochemical tests records at 5–10 weeks of gestation. Finally, we included 20,582 women in our study. The inclusion and exclusion criteria can be seen in Figure 1. The studies involving human participants were reviewed and approved by Institutional Review Board of Peking University Health Science Center (No. IRB00001052-21023).

FIGURE 1
www.frontiersin.org

Figure 1. The inclusion and exclusion criteria for the study population.

Maternal Characteristics, Medical History and Biomarkers

We extracted maternal characteristics, medical history, blood pressure measurements, blood routine examination and biochemical tests from the electronic data system of Tongzhou Maternal and Child Health Hospital of Beijing. Maternal characteristics included maternal height (cm), pre-pregnancy weight (kg), pre-pregnancy BMI, pre-gestational diabetes mellitus, ethnicity, parity, gravidity, abortion history, family history of hypertension, family history of diabetes, maternal age, husband age. Blood pressure measurement records included systolic blood pressure (SBP) and diastolic blood pressure (DBP) measurement records and we calculated mean arterial pressure (MAP) by MAP = (SBP+2 × DBP)/3. Blood routine examination and biochemical tests included 48 biomarkers,including hemoglobin, mean red blood cell volume, platelet, red blood cell, white blood cell, alanine aminotransferase, aspartate aminotransferase, urea nitrogen, and calcium, etc. Lists of all biomarkers were shown in Supplementary Tables 1a,b.

Maternal Outcomes

Based on the Guidelines on Diagnoses and Treatments of Hypertensive Disorders in Pregnancy by Chinese Society of Obstetrics and Gynecology, pre-eclampsia is defined as hypertension first appeared after 20 weeks of pregnancy, systolic blood pressure ≥140 and/or diastolic blood pressure ≥90 mmHg, with proteinuria or with any of the following organs or systems involved: heart, lung, liver, kidney and other important organs, or abnormal changes in blood system, digestive system and nervous system, placenta and fetus involved, etc (27). According to the time of diagnosis, we divided pre-eclampsia into early-onset pre-eclampsia (<34+0 weeks of gestation), preterm pre-eclampsia (<37+0 weeks of gestation), late-onset pre-eclampsia (≥34+0 weeks of gestation) and term pre-eclampsia (≥37+0 weeks of gestation).

Statistics

Division of Gestational Weeks

The first trimester, the second trimester and the third trimester were defined as < = 13 weeks, 14–27 weeks and > = 28 weeks, respectively (24, 28). We divided the whole pregnancy into eight periods: 5–10 weeks, 11–13 weeks, 14–18 weeks, 19–23 weeks, 24–27 weeks, 28–31 weeks, 32–35 weeks and 36–39 weeks. Therefore, the first trimester, the second trimester and the third trimester of pregnancy included 2, 3 and 3 time periods, respectively.

Preprocessing Variables and Imputing Missing Values

The pregnant women's height, pre pregnancy weight, pre pregnancy BMI, pregnant women's age and husband's age were used as continuous variables. We divided pre-gestational diabetes mellitus, family history of hypertension, abortion history, family history of diabetes into yes or no and we divided ethnicity into Han nationality or other nationalities. Parity was divided into primiparous or multiparous, and gravidity was divided into the first pregnancy, the second pregnancy, the third and more pregnancies. We divided the pregnancy season into spring, summer, autumn and winter according to the date of the last menstruation after verification.

We also kept blood pressure, routine blood and biochemical tests as continuous variables. If pregnant women had more than one measurement record during a certain period, we took the average value of these records. Therefore, regardless of missing values, each pregnant woman should have eight records of mean arterial pressure and biomarker measurements throughout pregnancy. There was one measurement record in each of the eight time periods. In particular, pregnant women rarely did biochemical tests in the second trimester of pregnancy. Therefore, there was no such test value in the 14–18 weeks, 19–23 weeks, and 24–27 weeks, so the variable was marked as missing.

The missing rate of demographic characteristic variables was <2%. We used the median to impute the missing values of continuous variables and the mode to impute the missing values of categorical variables. For the mean arterial pressure measurements, blood routine examination and biochemical tests (except the second trimester of pregnancy), the missing values were imputed by the values in the previous period. For example, missing values at 32–35 weeks of gestation were imputed by the measurements at 28–31 weeks of gestation. A pregnant woman has a prenatal examination at 32–35 weeks of pregnancy, but for some reason, she lacked a blood routine examination. Then, this missing value can be replaced by the value tested during 28–31 weeks. In clinical application, this method is easy for health professionals to understand and use it. If a pregnant woman is unable to get biochemical indicators for some reasons, this method can be used to quickly impute the missing value. In addition, by using previous measurements, this method can avoid reverse causality. The missing data imputed by the last observation carried forward (LOCF) method was used in many clinical antidepressant trials and other high level clinical studies (2931).

Development and Validation of Prediction Model

We sampled from women with and without pre-eclampsia to form a training set (N = 15665), a test set (N = 3917) and a validation set (N = 1000). The incidence rate of pre-eclampsia was almost the same in each data set (3.5%). We developed prediction models on the training set and verify the model on the test set. Based on the calculation results from the training set and test set, a screening strategy was formed and tested on the validation set. We did univariate analysis on demographic characteristics to compare the differences between each group and the control group (women without any types of pre-eclampsia), for example, comparison of women with early-onset pre-eclampsia and without any types of pre-eclampsia. We used Chi-test for categorical variables (Fisher exact test when the counts in some cells was fewer than five) and Kruskal-Wallis Rank Sum test for continuous variables.

We developed prediction models by using Logistic Regression on the imputed training set and calculated the area under receiver operator curve (AUROC) on the test set to reflect the ability of prediction. We used LASSO technique to select the variables with the highest priority and the second priority to join in the model according to the order in which the coefficients enter into the model along the solution path.

Least Absolute Shrinkage And Selection Operator (LASSO) is a shrinking technique used to do variable selection. It has been used to do variable selection and avoid overfitting for development of pre-eclampsia prediction models in previous studies (25, 32). LASSO was described as follows. For more details, we recommended readers to related materials (33, 34).

The objective of LASSO is to find the proper coefficients to minimize the loss function:

min - 1NΣi = 1Nl(yi|Xi)+λΣj = 1J|βj|,    (1)

where Xi = (xi1, xi2, ……, xiN) is the vector of observed variables of the subject i, yi is the outcome of the subject i, yiϵ{0, 1}, yi = 1 if the subject i developed pre-eclampsia and yi = 0 if the subject i did not develop pre-eclampsia. |βj| is the absolute value of the coefficient of the variable xij. λ is a tuning parameter: when it is big enough, all coefficients are zeros and the coefficients gradually turn to non-zero with the decrease of λ. l(yi|Xi) is the log-likelihood function of the logistic regression, given by the following equation:

l(yi|Xi)=yi·(β0+Σj=1Jβjxij)                   log[ 1+exp(β0+Σj=1Jβjxij) ].    (2)

We determined a basic variable set, which was forced into models, and then selected the variables that can improve AUROC on the test set from the demographical variables, mean arterial pressure and biochemical markers in proper order in four stages. We developed predictions on each period (5–10 weeks, 11–13 weeks, 14–18 weeks, 19–23 weeks, 24–27 weeks, 28–31 weeks, 32–35 weeks and 36–39 weeks). The four stages of developing models for each period were described as follows.

Stage 1

As reported in previous studies, pre-pregnancy BMI, pre-gestational diabetes mellitus (PGDM), parity, family history of hypertension and maternal age were highly associated with pre-eclampsia, thus these five variables were forced into models and denoted as “basic variables” (1, 35, 36). Then, we added other demographical variables and medical history into the model, including maternal height, pre-pregnancy weight (kg), ethnicity, gravidity, abortion history, husband age (years), and pregnancy season. Last, we selected the optimal set of variables with the highest AUROC denoted as Stage1-Optimal Variables.

Stage 2

Based on Stage1-Optimal Variables, we added MAPs measured at baseline (5–10 weeks) and the current period and then tested whether MAPs could improve the ability of prediction. For example, for women took prenatal examinations during 28–31 weeks, we added MAPs measured at 5–10 weeks and 28–31 weeks into the model. The optimal set of variables with the highest AUROC was denoted as Stage2-Optimal Variables.

Stage 3

Based on Stage2-Optimal Variables, we added biomarkers measured at baseline (5–10 weeks) and the current period. By using LASSO, we selected the optimal set of variables with the highest AUROC denoted as Stage3-Optimal Variables. As mentioned before, we did not add biochemical tests measured in the second trimester, because it was rarely tested in the hospital.

Stage 4

At the current period, we further added all MAPs measured in previous periods based on Stage2-Optimal Variables and added all biomarkers measured in previous periods based on Stage3-Optimal Variables to check whether adding other variables could improve the AUROC.

We excluded women who had used aspirin during pregnancy and conducted sensitivity analysis by redeveloping optimal prediction models to test whether the valuable biomarkers we found could still improve the predictive ability of models (37).

After developing the best prediction model at each stage, we selected the appropriate risk cutoff values to make the sensitivity of the prediction models reach 0.95, 0.96, 0.97, 0.98, and 0.99 in each period. Then, we screened the population in the validation set with different sensitivities to form the best screening strategy.

We used R software (38) to do all calculations. The “glmnet” package in R was used to fit logistic regressions via LASSO (39). Details about how to perform it in R are available at vignettes of the “glmnet” package.

Results

Comparison of Characteristics Between Women With and Without Pre-eclampsia

Table 1 showed demographical characteristics of women with and without pre-eclampsia. Compared with those without pre-eclampsia, women with pre-eclampsia had higher body weight and body mass index (BMI), and had higher proportion of pre-gestational diabetes, primipara, first pregnancy, family history of hypertension, family history of diabetes, pregnancy in spring and winter, and gestational diabetes. The differences between the two groups were statistically significant. Similar results were seen in women with early-onset and preterm pre-eclampsia. However, there were not obvious difference in the proportion of pre-gestational diabetes mellitus, family history of hypertension, family history of diabetes, and pregnancy season between women with and without early-onset pre-eclampsia, between women with and without preterm pre-eclampsia.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of pregnancies in imputed datasets.

Results of Variables Selection

Results of Stage 1

Supplementary Tables 2a,b, showed that before 20 weeks of gestation, the AUROC of the model on basic variables was 0.68 for all pre-eclampsia, 0.73 and 0.74 for early-onset and preterm pre-eclampsia, respectively. The AUROC of late-onset pre-eclampsia and term preeclampsia were 0.68 and 0.67, respectively. The model developed during 24–27 weeks had similar results. However, additional demographical characteristics before 20 weeks of gestation and gestational diabetes mellitus at 24–27 weeks of gestation did not increase AUROC on the test set. The Stage1-Optimal Variables included pre-pregnancy BMI, pre-gestational diabetes mellitus (PGDM), parity, family history of hypertension and maternal age.

Results of Stage 2

Tables 2A,B showed the results of Stage 2. For all pre-eclampsia, the AUROC increased from 0.68 to 0.70 with the addition of mean arterial pressure (MAP) measured at 5–10 weeks. With the increase of gestational weeks, the AUROC gradually increased to 0.78 with addition of MAPs measured at 5–10 and 28–31 weeks, and remained stable at 32–35 weeks and 36–39 weeks. Similar trends were seen in the late-onset and term pre-eclampsia, where the AUROC gradually increased from around 0.68 to 0.77.

TABLE 2A
www.frontiersin.org

Table 2A. Performance of prediction model based on characteristics, medical history, mean arterial pressure, and biomarkers for all, early-onset and preterm pre-eclampsia.

TABLE 2B
www.frontiersin.org

Table 2B. Performance of prediction model based on characteristics, medical history, mean arterial pressure, and biomarkers for late-onset and term pre-eclampsia.

For early-onset and preterm pre-eclampsia, AUROC gradually increased from about 0.73 to 0.89 and 0.80 at 28–31 weeks of gestation, respectively. Generally, based on the Stage 1-Optimal Variables, the addition of MAPs increased AUROC for all pre-eclampsia and its subtypes, and the AUROC was higher for early-onset and preterm preeclampsia than that for late-onset and term pre-eclampsia. The Stage 2-Optimal Variables included Stage 1-Optimal Variables and MAPs measured at baseline (5–10 weeks) and at the current prenatal examination period.

Results of Stage 3

As can be seen in Tables 2A,B, based on Stage 2-Optimal Variables, we found that the AUROC of all types of pre-eclampsia can be slightly improved by adding uric acid test in 5–10 weeks and 11–13 weeks of gestation. At 18–23 weeks and 24–27 weeks of gestation, the addition of platelets tested at baseline and the current period can improve AUROC for all types of pre-eclampsia (0.79–0.86). Also, we found that adding uric acid and alkaline phosphatase tested in the current period at 28-31 weeks, 32–35 weeks and 36–39 weeks of pregnancy can significantly improve the prediction ability for all types of pre-eclampsia, especially for pre-eclampsia occurred after 31 weeks (including late-onset and term pre-eclampsia). For all pre-eclampsia, the AUROC increased from approximately 0.78 to 0.86 (28–31 weeks), 0.89 (32–35 weeks) and 0.95 (36–39 weeks), respectively. Stage3-Optimal Variables are as follows: (i) Stage2-Optimal Variables plus uric acid tested in the current period in 5–10 weeks and 11–13 weeks; (ii) Stage 2-Optimal Variables without adding biomarkers in 14–18 weeks; (iii) Stage 2-Optimal Variables plus platelets tested at baseline and the current period in 18–23 weeks and 24–27 weeks; (iv) Stage 2-Optimal Variables plus uric acid and alkaline phosphatase tested in the current period at 28–31 weeks, 32–35 weeks and 36–39 weeks.

Results of Stage 4

As can be seen in Supplementary Tables 3a,b, we found that adding the MAPs measured in all previous periods to Stage 2-Optimal Variables and adding biomarkers tested in all previous periods to Stage 3-Optimal Variables almost did not improve AUROC. In order to prevent over fitting, we chose Stage 3-Optimal Variables to develop our final predictions models.

Results of Sensitivity Analysis

There were 104 and 28 women using aspirin in the training set (N = 15665) and test set (N = 3917), respectively. For the optimal prediction models in each stage, the results were very similar in the datasets with and without women who had used aspirin, which were shown in Supplementary Tables 4a–e. Although a few women had used aspirin during pregnancy, it was less likely to affect our results.

Best Prediction Model in Each Period

The Best Prediction Models in Each Period Are Shown as Below

The risk score for pre-eclampsia of each woman is calculated by odds/(1+odds), where odds = eY.

Period 1 (5–10 Weeks)

Y = (−13.40) + (0.102*BMI) + (−0.920 if Multiparous) + (0.0445* Maternal Age) + (0.433* if Pre–gestational diabetes mellitus) + (0.320* if Family history of hypertension) + (0.0551*).

Period 2 (11–13 Weeks)

Y = (−15.10) + (0.0890*BMI) + (−0.888 if Multiparous) + (0.0479* Maternal Age) + (0.407* if Pre–gestational diabetes mellitus) + (0.319* if Family history of hypertension) + (0.0233* MAP measured at 5–10 weeks) + (0.0503* MAP measured at 11–13 weeks) + (0.00609* uric acid test at 11–13 weeks).

Period 3 (14–18 Weeks)

Y = (−15.70) + (0.110*BMI) + (−0.903 if Multiparous) + (0.0460* Maternal Age) + (0.499* if Pre-gestational diabetes mellitus) + (0.347* if Family history of hypertension) + (0.0182* MAP measured at 5–10 weeks) + (0.0704* MAP measured at 14–18 weeks).

Period 4 (19–23 Weeks)

Y = (−15.80) + (0.110*BMI) + (−0.917 if Multiparous) + (0.0531* Maternal Age) + (0.453* if Pre-gestational diabetes mellitus) + (0.411* if Family history of hypertension) + (0.0253* MAP measured at 5–10 weeks) + (0.0648* MAP measured at 19–23 weeks) + (−0.0237* platelets test at 5–10 weeks) + (0.0231* platelets test at 19–23 weeks).

Period 5 (24–27 Weeks)

Y = (−16.10) + (0.110*BMI) + (−0.929 if Multiparous) + (0.0531* Maternal Age) + (0.400* if Pre–gestational diabetes mellitus) + (0.408* if Family history of hypertension) + (0.0274* MAP measured at 5–10 weeks) + (0.0664* MAP measured at 24–27 weeks) + (−0.0222* platelets test at 5–10 weeks) + (0.0217* platelets test at 24–27 weeks).

Period 6 (28–31 Weeks)

Y = (−15.60) + (0.0778*BMI) + (−0.941 if Multiparous) + (0.0611* Maternal Age) + (0.410* if Pre-gestational diabetes mellitus) + (0.375* if Family history of hypertension) + (0.0262* MAP measured at 5–10 weeks) + (0.0778* MAP measured at 28–31 weeks) + (0.00974* uric acid test at 28–31 weeks) + (−0.0578* alkaline phosphatase test at 28–31 weeks).

Period 7 (32–35 Weeks)

Y = (−16.00) + (0.0586*BMI) + (−0.968 if Multiparous) + (0.0445* Maternal Age) + (0.509* if Pre–gestational diabetes mellitus) + (0.455* if Family history of hypertension) + (0.0203* MAP measured at 5–10 weeks) + (0.110* MAP measured at 32–35 weeks) + (0.00938* uric acid test at 32–35 weeks) + (−0.0751* alkaline phosphatase test at 32–35 weeks).

Period 8 (36–39 Weeks)

Y = (−14.90) + (0.0252*BMI) + (−1.430 if Multiparous) + (0.0357* Maternal Age) + (0.617* if Pre-gestational diabetes mellitus) + (0.332* if Family history of hypertension) + (−0.0144* MAP measured at 5–10 weeks) + (0.129* MAP measured at 36–39 weeks) + (0.0122* uric acid test at 36–39 weeks) + (−0.0600* alkaline phosphatase test at 36–39 weeks).

Screening Strategy

The screening strategy was shown in Table 3 and Figure 2. We selected the 0.01, 0.011, 0.008, 0.009, 0.01, 0.003, 0.001, and 0.002 as risk cutoff values for each stage to make the sensitivity reach 0.98 (Table 3). Figure 2 showed the screening results of the validation set using the optimal prediction models with the sensitivity of 0.98. From 5–10 weeks of gestation, 1000 women participated in prenatal examination. Risk scores were calculated by the optimal prediction models for each pregnant woman participating in prenatal examination. Pregnant women exceeding the risk cutoff value were classified as high-risk group, and the risk assessment will be carried out in the next prenatal examination period; Otherwise, they were classified as the low-risk group and will not be examined in subsequent prenatal examinations. According to this method, after the prenatal examination of 14–18 weeks of pregnancy, 796 women were divided into the high-risk group, of which 35 will develop pre-eclampsia, and 204 women were divided into the low-risk group, of which 0 will develop pre-eclampsia. The number of women who need to continue screening decreased by about 20%. After the prenatal examination at 32–35 weeks of gestation, 556 women were divided into high-risk group, of which 28 will develop pre-eclampsia; 434 women were divided into the low-risk group, of which 3 will develop pre-eclampsia. The number of women who need continue screening decreased by about 45%.

TABLE 3
www.frontiersin.org

Table 3. The optimal multi–stage prediction models with the sensitivity of 0.98.

FIGURE 2
www.frontiersin.org

Figure 2. Screening results of the validation set using the optimal prediction models with the sensitivity of 0.98. From 5–10 weeks of gestation, 1000 women participated in prenatal examination.

Discussion

Summary

In this study, we established multi-stage pre-eclampsia risk prediction models throughout pregnancy based on 20582 pregnant women in China. We sequentially select valuable variables from demographical characteristics, mean arterial pressure, blood routine examination, and biochemical biomarkers to the prediction models. We found that uric acid tested from 5–10 weeks of gestation, platelets tested at 18-23 and 24–31 weeks of gestation, and alkaline phosphatase tested at 28–31, 32–35 and 36–39 weeks of gestation can further improve the prediction performance of models. The AUROC of the optimal prediction models on the test set gradually increased from 0.71 at 5–10 weeks to 0.80 at 24–27 weeks, and then gradually increased to 0.95 at 36–39 weeks of gestation. Based on the optimal prediction models, we established a multi-stage screening strategy from 5–10 weeks of pregnancy, which can add about 94.3% of women who will develop pre-eclampsia and reduce about 40% of the healthy women to be screened by 28–31 weeks of pregnancy.

Comparison With Previous Studies

Some studies used demographic characteristics and medical history to build prediction models. The AUROC of prediction models developed by David Wright et al. for all pre-eclampsia, preterm pre-eclampsia, and early-onset pre-eclampsia were 0.76, 0.79, and 0.81, respectively, which were higher than 0.68, 0.74, and 0.73 in our study (40). Similar results were seen in models given by LCY Poon et al. (41). This may be because more high-risk groups were included in these two studies, for example, the model established by David Wright et al. included high-risk groups such as chronic hypertension, systemic lupus erythematosus or antiphospholipid syndrome, and assisted reproduction. The prediction model may perform better on high-risk groups.

We found that on the basis of demographic characteristics and medical history, for all pre-eclampsia, the addition of mean arterial pressure can improve the AUROC of the model. The later the gestational week for the measurement of MAPs, the higher the AUROC. It gradually increased from 0.68 to 0.78 at 28–31 weeks of gestation, and then remained stable at 32–35 weeks and 36–39 weeks of gestation. This finding is basically consistent with previous studies (4244). Based on 12996 British women, Wallis et al. found that the AUROC increased from 0.79 to 0.80, 0.80, 0.84, 0.84, 0.84, 0.86, and 0.88 after adding the mean arterial pressure measured at 20, 25, 28, 31, 34 and 36 weeks to models with demographic characteristics, medical history, and MAP measured at baseline (43). Tayyar A. et al. also found the predictive value of adding MAPs at 12, 22, 32 and 36 weeks of gestation (44).

We found that uric acid tested from 5–10 weeks of gestation, platelets tested at 18–23 and 24–31 weeks of gestation, and alkaline phosphatase tested at 28–31, 32–35 and 36–39 weeks of gestation can further improve the prediction performance of models. There are some potential explanations for uric acid, platelets, and alkaline phosphatase to improve prediction performance. Platelets and uric acid have predictive values for pre-eclampsia. A systematic review of 69 studies showed that pregnant women with pre-eclampsia had a higher average platelet volume than the normal (45). Placental related diseases are associated with the transitional activation of maternal platelets, such as preeclampsia (46). Therefore, the detection of platelet function and over activation has a certain predictive value for preeclampsia (47). There was a strong association between uric acid and pre-eclampsia (48). The increased blood pressure can lead to organ damage, such as liver or kidney damage. Renal injury may increase the level of serum uric acid, and then increase the risk of pre-eclampsia (49, 50). It was reported that pregnant women with subsequent preeclampsia had elevated uric acid levels as early as 10 weeks of gestation (51). This phenomenon is basically consistent with the findings of this study. For alkaline phosphatase, some studies showed that the level of alkaline phosphatase was higher in women with pre-eclampsia than that in normal women (52, 53). Alkaline phosphatase may be related to pre-eclampsia, but it is necessary to further study its predictive value.

Comparison With Chinese Population-Based Research

In recent years, there were some studies developing prediction models of pre-eclampsia based on Chinese population (2126). These studies were mainly based on the hospital electronic medical record data system and mainly carried out in eastern China (three in Shanghai, one in Tianjin and one in Shanxi), screening pregnant women in the first and second trimester of pregnancy. Two studies predicted early-onset pre-eclampsia where the detection rate was between 40.7%−73.2% with the false positive rate of 10%. Four studies predicted all pre-eclampsia (undifferentiated subtypes), with AUROC from 0.86 to 0.98.

However, the populations and predictors used in these studies were quite different. Jiang et al. used demographic characteristics and biochemical markers (a complete blood count, serum albumin, serum uric acid, 24-h urinary protein, antinuclear antibodies, anti–double-stranded DNA antibody, antiphospholipid antibodies, etc.) to achieve an AUROC of 0.975 for pregnant women with systemic lupus erythematosus (23). Chen et al. predicted the risk of early-onset preeclampsia for pregnant women with twin pregnancy, where the AUROC was 0.82 (95% CI: 0.76–0.88), and the detection rate was 40.7% (false positive 10%) (22). Wang et al. developed a prediction model for pre-eclampsia using 31 blood flow related parameters such as vascularization index (VI), blood flow index (FI) and vascularization blood flow index (VFI) related to uterus and placenta, and the AUROC reached 0.877 (25). However, these studies aimed at specific high-risk groups, or used expensive indicators that were not within the scope of routine testing, so it may not be conducive to promotion. In addition, these studies also have other limitations: not consider pre-eclampsia subtypes (2326), possible bias caused by variable screening process (23, 24, 26), and insufficient number of outcome events (2123, 25).

Compared With the Screening Strategies of Previous Studies

Wallis et al. proposed a multi-stage screening strategy based on the British women: the first screening was conducted at < 18 weeks of gestation, and the second to fifth screening was conducted by using mean arterial pressure at 20, 25, 28, and 31 weeks of gestation. We and Wallis et al. tested the screening strategy based on the same number of pregnant women. Compared with the screening strategy proposed by Wallis et al., we advanced the prediction window as early as 5–10 weeks. After completing the screening at 28–31 weeks, we reduced the number of people to be screened by about 45%, while Wallis et al. reduced by about 35%. In particular, we extended the screening period to 36–39 weeks, which may effectively screen term pre-eclampsia and even later pre-eclampsia. Although the two strategies target different populations and it is necessary to compare the screening strategies of the two sides in the same population in the future, our screening strategy has the potential for improvement.

Advantages, Clinical Value and Limitations of the Study

We have several advantages in our study. First, our sample size and the number of outcomes were larger than that in previous prediction models based on Chinese population. Second, compared with the stepwise regression used in many studies, we used LASSO for variable selection, so it is less likely to over fitting, and we tested our model on the randomly divided test set and validation set. Third, we found that on the basis of demographic characteristics, medical history, and mean arterial pressure, the addition of uric acid can improve the prediction ability of the model from the first trimester; the addition of platelet and alkaline phosphatase can improve the prediction ability of the model in the second and third trimester, respectively. We did not test the predictive value of uric acid in the second trimester, because pregnant women rarely did biochemical tests in this period. We suggest that pregnant women in the second trimester do additional uric acid tests in the future. In particular, at 28–31 weeks, 32–35 weeks and 36–39 weeks of pregnancy, the AUROC of the prediction model for all pre-eclampsia reached 0.86, 0.89 and 0.95, respectively, which has the value in predicting late-onset and term pre-eclampsia. Fourth, to our knowledge, our study is the first study to develop multi-stage prediction models and propose a screening scheme based on the Chinese women. The screening time covers almost the entire pregnancy period: the earliest to 5–10 weeks of pregnancy and the latest to 36–39 weeks of pregnancy.

Several clinical values can be seen in our study. First, we found that uric acid, platelets and alkaline phosphatase can improve the predictive ability of the model, which are part of the routine test items and it does not need extra costs. Second, a screening system has been developed for the Chinese population and policy environment, and the prediction window has been advanced to 5–10 weeks, which allows time for aspirin intervention and other means for high-risk groups. There were strong evidences supporting that aspirin is the only drug preventing pre-eclampsia (54). Some guidelines recommended women with high or moderate clinical risk factors to use aspirin starting at the first or second trimester (1, 35, 36, 55). ACOG recommended using aspirin starting between 12 and 28 weeks of gestation, ideally before 16 weeks, and ISSHP also supported it was ideally used before 16 weeks (1, 55). A clinical trial involving 14361 women suggested that using low-dose aspirin starting at 6–13 weeks of gestation could reduce risk of preterm pre-eclampsia and perinatal mortality in low-income and middle-income countries (56). The study advanced the recommended time of aspirin use to 6 weeks of gestation. Third, for women identified at high risks during the second or third trimester, and missing the best window for using aspirin, clinicians could take intensive surveillance or hospitalization, and carefully select the best delivery time for pregnant women. Some previous studies developed prediction models in the second or third trimester (43, 5760). It is useful to develop multi-stage prediction models covering the first, second, and third trimester. Forth, in the third trimester of pregnancy, the prediction ability of late-onset and term pre-eclampsia was improved by adding uric acid and alkaline phosphatase. Although late-onset and term pre-eclampsia have a lower hazard than early-onset, preterm pre-eclampsia, the former two subtypes have a higher incidence rate. Doctors can treat pre-eclampsia by delivery, however, the early birth of newborns may still be detrimental to their future growth.

However, the study has several limitations. First, the study had a large sample size and a large number of outcomes, but the incidence rate of early-onset pre-eclampsia was small, because the incidence rate of early-onset pre-eclampsia was particularly low. We used LASSO for variable selection, so it was less likely to cause over fitting. Second, there were missing values in biomarkers, but it was in line with the actual clinical situation. We used a simple imputation method, that is, using the previous tested values to impute the missing values. This method is easy for health professionals to use. Third, we did not conduct external validation. But we conducted internal validation on the test set and validation set by randomly dividing the data set. Moreover, we used LASSO to select variables, which reduced the possibility of over fitting. In the future, our model awaits validation on other data sets.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions the data that support the findings of this study are available from Tongzhou Maternal and Child Health Care Hospital of Beijing but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of Tongzhou Maternal and Child Health Care Hospital of Beijing. Requests to access these datasets should be directed to Haijun Wang, whjun@pku.edu.cn.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Peking University Health Science Center (No. IRB00001052-21023). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

ZT, YJ, SZ, TS, ZY, and NH: data curation. ZT, YJ, JJ, and HW: methodology. ZT and JJ: software. ZT: writing—original draft. YJ, JJ, and HW: writing—review & editing. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by Beijing Natural Science Foundation (7212144) and National Natural Science Foundation of China (92046019).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We sincerely thank the research group of Peking University Retrospective Birth Cohort in Tongzhou based on the hospital information system. We appreciated the health professionals in Tongzhou Maternal and Child Health Care Hospital of Beijing for data collection and management.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.911975/full#supplementary-material

References

1. Brown MA, Magee LA, Kenny LC, Karumanchi SA, Mccarthy FP, Saito S, et al. The hypertensive disorders of pregnancy: ISSHP classification, diagnosis & management recommendations for international practice. Pregnancy Hypertens. (2018) 13:291–310. doi: 10.1161/HYPERTENSIONAHA.117.10803

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Umesawa M, Kobashi G. Epidemiology of hypertensive disorders in pregnancy: prevalence, risk factors, predictors and prognosis. Hypertens Res. (2017) 40:213–20. doi: 10.1038/hr.2016.126

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Obukhov AG, Ye C, Ruan Y, Zou L, Li G, Li C, et al. The 2011 Survey on Hypertensive Disorders of Pregnancy (HDP) in China: prevalence, risk factors, complications, pregnancy and perinatal outcomes. PLoS ONE. (2014) 9:e100180. doi: 10.1371/journal.pone.0100180

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ackerman CM, Platner MH, Spatz ES, Illuzzi JL, Xu X, Campbell KH, et al., Severe cardiovascular morbidity in women with hypertensive diseases during delivery hospitalization. Am J Obstet Gynecol. (2019) 220:582 e1-582 e11. doi: 10.1016/j.ajog.2019.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Stevens W, Shih T, Incerti D, Ton TGN, Lee HC, Peneva D, et al. Short-term costs of preeclampsia to the United States health care system. Am J Obstet Gynecol. (2017) 217:237–48 e16. doi: 10.1016/j.ajog.2017.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Hitti J Sienas L Walker S Benedetti TJ„ Easterling T Contribution Contribution of hypertension to severe maternal morbidity. Am J Obstet Gynecol. (2018) 219:405 e1–405 e7. doi: 10.1016/j.ajog.2018.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Theilen LH, Meeks H, Fraser A, Esplin MS, Smith KR, Varner MW. Long-term mortality risk and life expectancy following recurrent hypertensive disease of pregnancy. Am J Obstet Gynecol. (2018) 219:107 e1–107 e6. doi: 10.1016/j.ajog.2018.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Staff AC. Long-term cardiovascular health after stopping pre-eclampsia. Lancet. (2019) 394:1120–1. doi: 10.1016/S0140-6736(19)31993-2

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Callaway LK, Lawlor DA, O'callaghan M, Williams GM, Najman JM Mcintyre HD. Diabetes mellitus in the 21 years after a pregnancy that was complicated by hypertension: findings from a prospective cohort study. Am J Obstet Gynecol. (2007). 197:492 e1–7. doi: 10.1016/j.ajog.2007.03.033

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Khashan AS, Evans M, Kublickas M, Mccarthy FP, Kenny LC, Stenvinkel P, et al. Preeclampsia and risk of end stage kidney disease: A Swedish nationwide cohort study. PLoS Med. (2019) 16:e1002875. doi: 10.1371/journal.pmed.1002875

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Carr DB, Newton KM, Utzschneider KM, Tong J, Gerchman F, Kahn SE, et al. Preeclampsia and risk of developing subsequent diabetes. Hypertens Pregnancy. (2009) 28:435–47. doi: 10.3109/10641950802629675

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Magnussen Eb VL. Smith Gd, Romundstad Pr, Hypertensive disorders in pregnancy and subsequently measured cardiovascular risk factors. Obstet Gynecol. (2009) 114:961–70. doi: 10.1097/AOG.0b013e3181bb0dfc

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mcdonald SD, Han Z, Walsh MW, Gerstein HC., Devereaux P.J. Kidney disease after preeclampsia: a systematic review and meta-analysis. Am J Kidney Dis. (2010) 55:1026–39. doi: 10.1053/j.ajkd.2009.12.036

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Brown MC, Best KE, Pearce MS, Waugh J, Robson SC., Bell R. Cardiovascular disease risk in women with pre-eclampsia: systematic review and meta-analysis. Eur J Epidemiol. (2013) 28:1–19. doi: 10.1007/s10654-013-9762-6

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Bokslag A, Teunissen PW, Franssen C, Van Kesteren F, Kamp O, Ganzevoort W, et al. Effect of early-onset preeclampsia on cardiovascular risk in the fifth decade of life. Am J Obstet Gynecol. (2017) 216:523 e1-523 e7. doi: 10.1016/j.ajog.2017.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Fields JA, Garovic VD, Mielke MM, Kantarci K, Jayachandran M, White WM, et al. Preeclampsia and cognitive impairment later in life. Am J Obstet Gynecol. (2017) 217:74 e1-74 e11. doi: 10.1016/j.ajog.2017.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Roberge S, Villa P, Nicolaides K, Giguere Y, Vainio M, Bakthi A, et al. Early administration of low-dose aspirin for the prevention of preterm and term preeclampsia: a systematic review and meta-analysis. Fetal Diagn Ther. (2012) 31:141–6. doi: 10.1159/000336662

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bujold E, Roberge S, Lacasse Y, Bureau M, Audibert F, Marcoux S, et al. Prevention of preeclampsia and intrauterine growth restriction with aspirin started in early pregnancy: a meta-analysis. Obstet Gynecol. (2010) 116:402–14. doi: 10.1097/AOG.0b013e3181e9322a

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Koopmans CM, Bijlenga D, Groen H, Vijgen SMC, Aarnoudse JG, Bekedam DJ, et al. Induction of labour versus expectant monitoring for gestational hypertension or mild pre-eclampsia after 36 weeks' gestation (HYPITAT): a multicentre, open-label randomised controlled trial. Lancet. (2009) 374:979–88. doi: 10.1016/S0140-6736(09)60736-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Antwi E, Amoakoh-Coleman M, Vieira DL, Madhavaram S, Koram KA, Grobbee DE, et al. Systematic review of prediction models for gestational hypertension and preeclampsia. PLoS ONE. (2020) 15:e0230955. doi: 10.1371/journal.pone.0230955

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Chang Y, Chen X, Cui HY, Li X., Xu YL. New Predictive Model at 11(+0) to 13(+6) Gestational Weeks for Early-Onset Preeclampsia With Fetal Growth Restriction Reprod Sci. (2017) 24:783–9. doi: 10.1177/1933719116669053

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Chen J, Zhao D, Liu Y, Zhou J, Zou G, Zhang Y, et al. Screening for preeclampsia in low-risk twin pregnancies at early gestation. Acta Obstet Gynecol Scand. (2020) 99:1346–53. doi: 10.1111/aogs.13890

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Jiang M, Wang Y, Fu Q., Lin S., Wu J„, Di W, Preeclampsia risk prediction model for chinese pregnant patients with systemic lupus erythematosus Arthritis Care Res. (2020) 72:1602–10. doi: 10.1002/acr.24265

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Mi B, Wen X, Li S, Liu D, Lei F, Liu X, et al. Parameterization of the mid-trimester drop in blood pressure trajectory during pregnancy and its utility for predicting preeclampsia. J Hypertens. (2020) 38:1355–66. doi: 10.1097/HJH.0000000000002395

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wang W, Wang Y, Yuan T, Zhang H, Li C, Li X, et al. Nomogram-based prediction of pre-eclampsia in the first trimester of gestation. Pregnancy Hypertens. (2020) 21:145–51. doi: 10.1016/j.preghy.2020.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yue CY, Gao JP, Zhang CY. Ni YH„ Ying CM. Development and validation of a nomogram for the early prediction of preeclampsia in pregnant Chinese women. Hypertens Res. (2020) 44:417–25. doi: 10.1038/s41440-020-00558-1

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Diagnosis and treatment guideline of hypertensive disorders in pregnancy 2020. Zhonghua fu chan ke za zhi. (2020) 50:227–38.

Google Scholar

28. World Health Organization. WHO Recommendations on Antenatal Care for a Positive Pregnancy Experience. Geneva: World Health Organization (2016).

Google Scholar

29. Abe K, Aoki M, Tsuji S, Itoyama Y, Sobue G, Togo M, et al. Safety and efficacy of edaravone in well defined patients with amyotrophic lateral sclerosis: a randomised, double-blind, placebo-controlled trial. Lancet Neurol. (2017) 16:505–12. doi: 10.1016/S1474-4422(17)30115-1

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Behr J, Prasse A, Kreuter M, Johow J, Rabe KF, Bonella F, et al. Pirfenidone in patients with progressive fibrotic interstitial lung diseases other than idiopathic pulmonary fibrosis (RELIEF): a double-blind, randomised, placebo-controlled, phase 2b trial. Lancet Respir Med. (2021) 9:476–86. doi: 10.1016/S2213-2600(20)30554-3

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Woolley SB, Cardoni AA. Goethe JW. Last-observation-carried-forward imputation method in clinical efficacy trials: review of 352 antidepressant studies. Pharmacotherapy. (2009) 29:1408–16. doi: 10.1592/phco.29.12.1408

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Baschat AA, Magder LS., Doyle LE, Atlas RO, Jenkins CB, Blitzer MG Prediction of preeclampsia utilizing the first trimester screening examination. Am J Obstet Gynecol. (2014) 211:514 e1–7 doi: 10.1016/j.ajog.2014.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

33. James G, Witten D, Hastie T„ Tibshirani R An Introduction to Statistical Learning. Vol. 112. Manhattan, NY: Springer (2013). doi: 10.1007/978-1-4614-7138-7

CrossRef Full Text | Google Scholar

34. Hastie T, Tibshirani R, Friedman JH„ Friedman JH. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. Vol. 2. Manhattan, NY: Springer (2009). doi: 10.1007/978-0-387-84858-7

PubMed Abstract | CrossRef Full Text | Google Scholar

35. 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 Gynaecol Obstet. (2019) 145 Suppl 1:1–33. doi: 10.1002/ijgo.12802

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Webster K, Fishburn S, Maresh M, Findlay SC, Chappell LC, Guideline C. Diagnosis and management of hypertension in pregnancy: summary of updated NICE guidance. BMJ. (2019) 366:l5119. doi: 10.1136/bmj.l5119

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cheong-See F, Allotey J, Marlin N, Mol BW, Schuit E, Ter Riet G, et al. Prediction models in obstetrics: understanding the treatment paradox and potential solutions to the threat it poses. BJOG. (2016) 123:1060–4. doi: 10.1111/1471-0528.13859

PubMed Abstract | CrossRef Full Text | Google Scholar

38. R Core Team. R. A language environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. (2019). Available online at: https://www.R-project.org/

Google Scholar

39. Noah Simon, Jerome Friedman, Trevor Hastie, Rob Tibshirani. Regularization paths for Cox's proportional hazards model via coordinate descent. J Stat Softw. (2011) 39, 1–13. https://www.jstatsoft.org/v39/i05/

PubMed Abstract | Google Scholar

40. Wright D, Syngelaki A, Akolekar R, Poon LC, Nicolaides KH. Competing risks model in screening for preeclampsia by maternal characteristics and medical history. Am J Obstet Gynecol. (2015) 213:62 e1–62 e10 doi: 10.1016/j.ajog.2015.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Poon LC, Kametas NA. Chelemen T. Leal A, Nicolaides KH, Maternal risk factors for hypertensive disorders in pregnancy: a multivariate approach. J Hum Hypertens. (2010) 24:104–10. doi: 10.1038/jhh.2009.45

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Cnossen JS, Vollebregt KC, De Vrieze N, Ter Riet G, Mol BW, Franx A, et al. Accuracy of mean arterial pressure and blood pressure measurements in predicting pre-eclampsia: systematic review and meta-analysis. BMJ. (2008) 336:1117–20. doi: 10.1136/bmj.39540.522049.BE

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Macdonald-Wallis C, Silverwood RJ, De Stavola BL, Inskip H, Cooper C, Godfrey KM, et al. Antenatal blood pressure for prediction of pre-eclampsia, preterm birth, and small for gestational age babies: development and validation in two general population cohorts. BMJ. (2015) 351:h5948. doi: 10.1136/bmj.h5948

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Tayyar A, Krithinakis K, Wright A. Wright D, Nicolaides KH. Mean arterial pressure at 12, 22, 32 and 36 weeks' gestation in screening for pre-eclampsia Ultrasound. Obstet Gynecol. (2016) 47:573–9. doi: 10.1002/uog.15815

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Jakobsen C, Larsen JB. Fuglsang J, Hvas A-M. Platelet function in preeclampsia - a systematic review and meta-analysis. Platelets. (2019) 30:549–62. doi: 10.1080/09537104.2019.1595561

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Moser G Guettler J Forstner D Gauster M Maternal Maternal Platelets—Friend or Foe of the Human Placenta? Int j mol sci. (2019). 20 doi: 10.3390/ijms20225639

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kazmi RS, Cooper AJ, Lwaleed BA. Platelet function in pre-eclampsia. Semin Thromb Hemost. (2011) 37:131–6. doi: 10.1055/s-0030-1270339

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Khaliq OP, Konoshita T. Moodley J, Naicker T, The Role of Uric Acid in Preeclampsia: Is Uric Acid a Causative Factor or a Sign of Preeclampsia? Curr Hypertens Rep. (2018) 20:80. doi: 10.1007/s11906-018-0878-7

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Hahn K, Kanbay M, Lanaspa MA, Johnson RJ, Ejaz AA. Serum uric acid and acute kidney injury: A mini review. J Adv Res. (2017) 8:529–36. doi: 10.1016/j.jare.2016.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Müller-Deile J., Schiffer M. Preeclampsia from a renal point of view: Insides into disease models, biomarkers and therapy. World J Nephrol. (2014) 3:169–81. doi: 10.5527/wjn.v3.i4.169

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wu Y, Xiong X, Fraser WD, Luo Z.-C. Association of uric acid with progression to preeclampsia and development of adverse conditions in gestational hypertensive pregnancies. Am J Hypertens. (2012) 25:711–7. doi: 10.1038/ajh.2012.18

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Duan Z, Li C, Leung WT, Wu J, Wang M, Ying C, et al. Alterations of Several Serum Parameters Are Associated with Preeclampsia and May Be Potential Markers for the Assessment of PE Severity. Dis Markers. (2020) 2020:7815214. doi: 10.1155/2020/7815214

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Chaparro A, Gaedechens D, Ramírez V, Zuñiga E, Kusanovic JP, Inostroza C, et al. Placental biomarkers and angiogenic factors in oral fluids of patients with preeclampsia. Prenat Diagn. (2016) 36:476–82. doi: 10.1002/pd.4811

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Chappell LC. Cluver CA. Kingdom J, Tong S, Pre-eclampsia. Lancet. (2021) 398:341–54. doi: 10.1016/S0140-6736(20)32335-7

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Gestational hypertension and preeclampsia: ACOG practice bulletin number 222. Obstet Gynecol. (2020) 135:e237–e260. doi: 10.1097/AOG.0000000000003891

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Hoffman MK, Goudar SS, Kodkany BS, Metgud M, Somannavar M, Okitawutshu J, et al. Low-dose aspirin for the prevention of preterm delivery in nulliparous women with a singleton pregnancy (ASPIRIN): a randomised, double-blind, placebo-controlled trial. Lancet. (2020) 395:285–93. doi: 10.1016/S0140-6736(19)32973-3

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Tsiakkas A, Saiid Y, Wright A, Wright D, Nicolaides KH. Competing risks model in screening for preeclampsia by maternal factors and biomarkers at 30-34 weeks' gestation. Am J Obstet Gynecol. (2016) 215:87 e1–87 e17 doi: 10.1016/j.ajog.2016.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Andrietti S, Silva M. Wright A. Wright D, Nicolaides KH, Competing-risks model in screening for pre-eclampsia by maternal factors and biomarkers at 35-37 weeks' gestation Ultrasound. Obstet Gynecol. (2016) 48:72–9. doi: 10.1002/uog.15812

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Tayyar A, Garcia-Tizon Larroca S. Poon LC. Wright D, Nicolaides KH, Competing risk model in screening for preeclampsia by mean arterial pressure and uterine artery pulsatility index at 30-33 weeks' gestation. Fetal Diagn Ther. (2014) 36:18–27. doi: 10.1159/000360792

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Ciobanu A, Wright A, Panaitescu A, Syngelaki A, Wright D, Nicolaides KH. Prediction of imminent preeclampsia at 35–37 weeks gestation. Am J Obstet Gynecol. (2019) 220:584 e1–584 e11 doi: 10.1016/j.ajog.2019.01.235

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: pre-eclampsia, multi-stage prediction model, screening strategy, pregnancy, logistic regression, LASSO

Citation: Tang Z, Ji Y, Zhou S, Su T, Yuan Z, Han N, Jia J and Wang H (2022) Development and Validation of Multi-Stage Prediction Models for Pre-eclampsia: A Retrospective Cohort Study on Chinese Women. Front. Public Health 10:911975. doi: 10.3389/fpubh.2022.911975

Received: 03 April 2022; Accepted: 04 May 2022;
Published: 30 May 2022.

Edited by:

Redhwan Ahmed Al-Naggar, National University of Malaysia, Malaysia

Reviewed by:

Mustafa Bayraktar, Atatürk University, Turkey
Doris Cristina Quintero-Lesmes, Fundación Cardiovascular de Colombia, Colombia
Jose Guida, State University of Campinas, Brazil

Copyright © 2022 Tang, Ji, Zhou, Su, Yuan, Han, Jia and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jinzhu Jia, jzjia@math.pku.edu.cn; Haijun Wang, whjun@pku.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.