Construction and Verification of Adverse Pregnancy Outcome Prediction Model for Pregnancy Complicated with Autoimmune Diseases
Objective: To construct and verify an individualized prediction model for adverse pregnancy outcomes in pregnant patients with autoimmune diseases. Methods: The data of patients admitted to the Department of Obstetrics, Affiliated Hospital of Qingdao University from January 2016 to December 2022 who were diagnosed with pregnancy combined with definite autoimmune diseases were retrospectively analyzed. Using univariate analysis, multivariate Logistic regression analysis and LASSO to screen independent risk factors for adverse pregnancy outcomes, a nomogram prediction model was established. Test set was used to evaluate the prediction efficiency, ROC curve and Bootstrop self-sampling method were used to verify the model, and the distinguishing ability and calibration of the model were further evaluated. Results: 1) Through univariate, multivariate and LASSO regression analysis, the optimal predictors of adverse pregnancy outcomes with autoimmune diseases were selected: hypertension, number of other complications, hemoglobin, N/L ratio, D-dimer. 2) The application of the above optimal predictors to establish a predictive model for adverse pregnancy outcomes of pregnancy combined with autoimmune diseases has a good performance, with a sensitivity of 74.3%, specificity of 82.5%, and AUC of 0.843 (95% CI: The sensitivity and specificity of the validation set model were 51.3%, 93.1%, and AUC was 0.75 (95% CI: 0.656~0.845). 3) The calibration curve was obtained by Bootstrop repeated self-sampling method, and the simulated curve had a strong consistency with the actual curve. Conclusion: In this study, a prediction model for adverse pregnancy outcomes of pregnancy with autoimmune diseases was established, which can be used as a relatively simple tool to predict the outcome of pregnancy with autoimmune diseases.
Autoimmune Disease
自身免疫性疾病(autoimmune disease, AID)是一类在遗传因素和外界环境因素综合作用下,破坏了机体自身平衡状态,免疫系统发生紊乱导致的自身组织器官受到免疫系统攻击的疾病,其发病机制不清,发病过程复杂多变,临床表现不一,结局多样。自身免疫性疾病的发病率女性明显高于男性,孕龄妇女多见,并且多在妊娠期间首次发现
回顾性收集2016年01月~2022年12月于青岛大学附属医院入院诊断为妊娠合并某种自身免疫性疾病的614名患者的临床资料,按7:3的比例完全随机化抽样的方式随机分为实验组(435人)和验证组(179人),实验组根据是否发生出现不良妊娠结局,将其分为不良结局组(92人)和非不良结局组(343人)。本研究经过青岛大学附属医院伦理委员会审核通过。
本研究为复合不良结局,只要发生其中任何一种不良妊娠结局,即认为患者发生了不良妊娠结局。孕产妇的不良结局包括:流产、胎停、胎死宫内、早产、出生低体重儿、新生儿窒息、产后出血等。
年龄、体重指数(BMI),初/经产妇,受孕方式、不良孕产史次数、治疗性流产次数、是否有家族史等;合并症;并发症及相关实验室检查包括血常规 + CRP、尿液检查、凝血检查、D-二聚体、红细胞沉降率、生化检查、补体检查等。
本研究采用R进行统计学分析,将两组资料具有统计学意义的变量纳入套索算法(Least Absolute Shrinkage and Selection Operator, LASSO)进行变量选择,将LASSO筛选出的这些变量再纳入多因素Logistic回归分析,得到最优预测变量。根据是否发生的不良妊娠结局作为因变量,在R软件中建立妊娠合并自身免疫性疾病不良妊娠结局的个性化nomogram预测模型,然后对模型进行验证及效能评价。通过Bootstrap重复自抽样法对列线图模型进行内部验证,经过Bootstrap自抽样1000次获得校准曲线进行评估,在软件中绘制出独立指标及列线图模型对不良结局预测价值对应的ROC曲线,并计算相应的AUC及95% CI,双侧P < 0.05时,认为有统计学意义。
选取2016年01月~2022年12月在青岛大学附属医院入院诊断为妊娠合并自身免疫性疾病且符合纳入标准及排除标准的共计614人,以7:3的比例完全随机化抽样分为实验组(435人)与验证组(179人),实验组根据患者是否发生不良结局分为不良结局组和非不良结局组。
根据是否发生不良结局,分为不良结局组(用1表示)和非不良结局组(用0表示)。
0 |
1 |
P |
|
n |
343 |
92 |
|
age (year) |
31.95 (4.51) |
31.17 (4.23) |
0.075 |
Pregnancy History = 1 (times) |
178 (37.0) |
47 (35.3) |
0.801 |
BMI (kg/m2) |
27.28 (3.75) |
26.85 (4.19) |
0.25 |
tube = 1 (times) |
58 (12.1) |
14 (10.5) |
0.739 |
Family History = 1 (times) |
12 (2.5) |
7 (5.3) |
0.177 |
Abnormal Pregnancy (times) |
0.98 (1.28) |
1.20 (1.44) |
0.095 |
hemoglobin (g/L) |
123.04 (13.31) |
118.39 (17.62) |
0.001 |
N/L |
4.43 (2.46) |
5.41 (4.94) |
0.002 |
bloodplatelet (109/L) |
187.59 (60.57) |
179.20 (77.38) |
0.185 |
CRP |
2.03 (5.19) |
3.09 (7.57) |
0.063 |
ESR (mm/L) |
28.54 (14.27) |
33.20 (17.82) |
0.002 |
APTT (s) |
31.04 (16.70) |
31.04 (5.52) |
0.996 |
INR |
0.91 (0.17) |
0.89 (0.08) |
0.162 |
DD (mg/L) |
1038.79 (709.99) |
856.39 (601.89) |
0.007 |
Total Protein (g/L) |
61.46 (5.88) |
59.28 (6.66) |
<0.001 |
albumin (g/L) |
40.38 (147.97) |
31.61 (4.14) |
0.495 |
Total Bilirubin (umol/L) |
9.30 (3.69) |
7.84 (3.48) |
<0.001 |
AST (U/L) |
17.14 (7.62) |
19.50 (13.14) |
0.008 |
ALT (U/L) |
14.56 (14.78) |
19.57 (23.24) |
0.003 |
Uric Acid (mmol/L) |
278.25 (80.11) |
323.53 (116.17) |
<0.001 |
Urea Nitrogen (mmol/L) |
3.43 (1.19) |
4.16 (2.01) |
<0.001 |
creatinine (umol/L) |
57.67 (16.56) |
55.62 (18.80) |
0.22 |
C1q (mg/L) |
158.92 (55.04) |
147.05 (33.49) |
0.018 |
Urine Occult Blood |
481 |
133 |
0.052 |
0 |
432 (89.8) |
101 (75.9) |
|
1 |
34 (7.1) |
18 (13.5) |
|
2 |
12 (2.5) |
11 (8.3) |
|
3 |
3 (0.6) |
3 (2.3) |
|
Urine Protein (%) |
<0.001 |
||
0 |
354 (73.6) |
61 (45.9) |
|
1 |
113 (23.5) |
34 (25.6) |
|
2 |
10 (2.1) |
17 (12.8) |
|
3 |
3 (0.6) |
19 (14.3) |
|
4 |
1 (0.2) |
2 (1.5) |
|
C3 (g/L) |
1.11 (0.22) |
1.25 (3.23) |
0.335 |
C4 (g/L) |
0.22 (0.13) |
1.46 (14.55) |
0.061 |
cesarean = 1 (times) |
370 (76.9) |
105 (78.9) |
0.706 |
others (times) |
0.5 (0.63) |
1.04 (0.8) |
<0.001 |
Geweeks (w) |
38.33 (1.02) |
32.48 (4.84) |
<0.001 |
Informal examination (times) |
18 (5.2) |
25 (27.2) |
<0.001 |
备注:1) tube:试管婴儿。2) 尿蛋白 < 0.15 g:阴性,用0表示;尿蛋白0.2~1 g/L:+,用1表示;尿蛋白1~2 g/L:++,用2表示;尿蛋白2~4 g/L:+++,用3表示;尿蛋白 > 4 g/L:++++,用4表示。3) 尿潜血阴性,用0表示;尿潜血(+),用1表示,尿潜血(++),用2表示;尿潜血(+++),用3表示。4) others:指与胎盘功能有关的并发症数量,如羊水过少、胎儿窘迫、B超示子宫血流异常等。
由
不良结局组较非不良结局组相比:是否正规产检、入院孕周、高血压、其它并发症数量,两组具有统计学意义(P < 0.05)。
通过
通过
将不良结局组和非不良结局组两组所有相关因素进行单因素Logistic回归分析,结果显示:高血压、其它并发症数量、Hb、N/L比值、ESR、D-D二聚体、总蛋白、胆红素、AST、ALT、尿酸、尿素氮、C1q、尿蛋白、入院孕周、未规范产检等16个指标认为有统计学意义,是发生不良妊娠结局的危险因素。
将单因素分析中具有统计学意义的因素纳入LASSO中进行变量筛选:如
将单因素分析中P < 0.05的变量经LASSO降维处理及多因素Logistic回归分析筛选出最终纳入的变量,因变量为是否出现不良妊娠结局,多因素Logistic回归结果:Hb (OR = 0.352, 95% CI = 0.2579~0.4816)、高血压(OR = 3.096, 95% CI = 1.050~9.126)、N/L (OR = 1.916, 95% CI = 1.215~3.028)、其它并发症数量(OR = 1.693, 95% CI = 1.047~2.738)、D-二聚体(OR = 0.996, 95% CI = 0.995~0.998);这些均是妊娠合并自身免疫性疾病不良结局得独立危险因素(P < 0.05)。进一步判断上述独立危险因素有无多重共线性,方差膨胀因子(VIF)分别是1.5、1.2、1.4、1.9、1.1,说明变量间不存在共线性。
将多因素筛选出来的变量,纳入二分类Logistic回归分析,以建立预测模型(如
OR |
SE |
Z |
P > IZI |
95% CI |
|
hemoglobin |
0.352 |
0.056 |
−6.54 |
0.000 |
0.258~0.482 |
hypertention |
3.096 |
1.708 |
2.05 |
0.04 |
1.050~9.126 |
N/L |
1.916 |
0.001 |
−6.24 |
0.000 |
1.215~3.028 |
others |
1.693 |
0.415 |
2.15 |
0.032 |
1.047~2.738 |
D-D |
0.996 |
0.060 |
2.084 |
0.037 |
0.995~0.998 |
基于多因素Logistic回归结果,应用R软件在预测模型中列入五个独立危险因素,然后建立妊娠合并自身免疫性疾病不良妊娠结局的个性化列线图预测模型。列线图如下
例如:一名35岁的妊娠合并自身免疫性疾病患者,入院时高血压(47分),有其它两个并发症(40分),Hb 80 g/L (65分),N/L5 (8分),D-二聚体1000 (50分),各预测指标累计得分约47 + 40 + 65 + 8 + 50 = 210分,相应预测风险 > 70%,表明该患者发生不良妊娠结局为高风险。
由
在实验组中计算并绘制校准曲线(Calibration curve)。将模型预测概率所对应的实际发生概率绘制成一条曲线即为校准曲线,可以用来预测概率与实际概率的一致程度。校准曲线距离对角线越远,说明预测的误差越大。通过Bootstrap重复自抽样法对列线图模型进行内部验证,经过Bootstrap自抽样1000次获得校准曲线,模拟曲线与实际曲线绝对误差是0.016,如
利用模型得出的结果,绘制出临床决策曲线,如
以上结果显示,应用Logistic回归建立的列线图预测模型能够较好的进行妊娠合并自身免疫性疾病不良结局的预测。
自身免疫性疾病(autoimmune disease, AID)是一种慢性全身性疾病,通常影响生育年龄妇女,出现大量的不良妊娠结局。近些年通过对自身免疫性疾病的研究及治疗的进展,使绝大多数病情稳定的自身免疫性疾病患者能够正常妊娠并安全度过妊娠期,获得良好的妊娠结局。
自身免疫性疾病是发生高血压的危险因素
SLE、APS等自身免疫性疾病存在较明显凝血功能异常和易栓倾向。D-二聚体是一种近几年来研究和应用较为广泛的凝血指标
血栓性疾病与炎症反应是密切相关的
本实验并未发现体重指数、胆固醇对于妊娠结局的影响
妊娠合并自身免疫性疾病属于高危妊娠,导致母胎不良妊娠结局显著增高
本研究选用的预测模型的数据来源于同一的三级甲等医疗机构,不良结局的样本量相对较少,且缺少中小医院样本(目前中小医院对于自身免疫性疾病认识不足诊疗相对较低,病例少)进行外部验证。但本研究为回顾性研究,不可避免的出现选择偏倚。除此外,患者的CRP及免疫性指标不是我们入院的常规检查中,且有明确意义激素用量等信息存在严重数据缺失,因此这些数据未被纳入本研究。
应用筛选出的5个最优预测变量建立妊娠合并自身免疫性疾病不良妊娠结局的预测模型的性能较好,经内部验证显示该模型准确性高,有较高的诊断价值。
有利于临床医师早期识别妊娠合并自身免疫性疾病患者的高危患者,尽早采取必要的干预措施,改善母婴不良结局。
*通讯作者。