支气管肺发育不良(bronchopulmonary dysplasia, BPD)是早产儿最常见的慢性肺疾病,也是导致神经发育障碍、长期呼吸系统并发症及死亡的主要原因
[1]。尽管围产医学和新生儿重症监护技术不断进步,但BPD的发生率仍居高不下,尤其在极早产儿及超早产儿中更为显著
[2-3]。目前,最有效的预防策略是避免极端早产
[4],但当早产不可避免时,在生后关键窗口期(如生后7~14 d)进行精准干预至关重要
[5-6]。尽管已有许多干预措施在临床试验中进行评估,但治疗效果有限
[7]。部分原因在于缺乏可靠的早期预测工具,导致无法精准识别高危患儿,甚至可能造成无效或有害的干预
[8-9]。因此,开发适用于胎龄<32周早产儿的BPD动态预测模型,对个体化风险评估、优化早期干预具有重要意义。
理想的BPD预测模型应具备易获取、高准确性、可预测远期结局、反映病理机制及多中心验证等特点
[10]。此类模型不仅有助于早期识别高风险患儿,降低医疗成本,还能为个体化干预提供依据
[11]。目前大部分预测模型,如2011年美国国家儿童健康与人类发展研究所(National Institute of Child Health and Human Development, NICHD)团队构建的多时间点预测模型
[12]、2019年Valenzuela-Stutman等
[13]构建的动态预测模型均基于过时的临床数据,且鲜有适用于床旁决策的工具
[14]。本研究基于NICHD 2018年BPD诊断标准
[6],通过整合围产期特征、疾病严重度指标[如新生儿急性生理学评分Ⅱ(Score for Neonatal Acute Physiology Ⅱ, SNAP-Ⅱ)
[15]],以及呼吸支持、营养支持等数据,构建生后第1、7、14、21及28天的多时间点动态预测模型及相应的列线图,为临床提供动态的可视化风险量化工具,并为BPD的早期诊断和个体化干预提供依据。
1 资料与方法
1.1 研究对象
回顾性收集2016年1月—2020年11月期间在中南大学湘雅二医院出生并转入新生儿重症监护室(neonatal intensive care unit, NICU)的早产儿的临床资料。研究对象的纳入流程见
图1。
纳入标准:(1)在中南大学湘雅二医院产科分娩并于生后24 h内入住NICU;(2)出生胎龄<32周的早产儿。
排除标准:(1)严重先天畸形(如复杂紫绀型先天性心脏病、膈疝等);(2)先天性遗传病或代谢性疾病(如21三体综合征、苯丙酮尿症等);(3)肺或肺血管先天性发育不良。
剔除标准:(1)14 d内死亡,或在14 d至校正胎龄36周之前因非持续性实质性肺部疾病和呼吸衰竭以外的原因死亡;(2)因其他原因放弃治疗的意外出院;(3)关键变量缺失>20%。
本研究已获中南大学湘雅二医院医学伦理委员会批准,批准文号:(2020)伦审[临研]第076号,并取得患儿监护人的知情同意。
1.2 分组
根据2018年的NICHD BPD定义
[6],将研究对象分为BPD组和非BPD组。2016年1月—2019年12月入住NICU的新生儿为模型构建组;2020年1—11月住院的新生儿为模型验证组。
1.3 临床资料收集
收集研究对象的临床数据:基本资料(胎龄、出生体重、SNAP-Ⅱ评分等);母亲围产期资料(产前使用糖皮质激素、绒毛膜羊膜炎、胎膜早破等);营养支持情况(开奶时间、禁食时间、平均加奶速度、热量摄入等);呼吸支持情况[用氧总时间、机械通气时间、吸入氧浓度(fraction of inspired oxygen, FiO2)等];治疗干预措施[外源性肺表面活性物质(pulmonary surfactant, PS)干预、静脉使用地塞米松、物理治疗开始日龄等];实验室检查结果及住院期间各系统并发症[痰培养等结果,以及肺炎、败血症、呼吸窘迫综合征(respiratory distress syndrome, RDS)、有血流动力学意义的动脉导管未闭(hemodynamically significant patent ductus arteriosus, hsPDA)、坏死性小肠结肠炎(necrotizing enterocolitis, NEC)、低钙血症等并发症]。
1.4 BPD多时间点临床预测模型的构建
通过单因素logistic回归分析和联合受试者操作特征曲线(receiver operating characteristic curve, ROC曲线),筛选出患儿生后第1、7、14、21、28天5个时间点的危险因素及有诊断价值的指标。将单因素分析中P值小于0.05的变量纳入多因素logistic回归分析,最终根据ROC曲线下面积(area under the curve, AUC)选定对预测BPD贡献度最大的5个变量进入最终预测模型。
1.5 统计学分析
采用SPSS 25.0进行统计学分析。计量资料以均数±标准差或中位数(四分位数间距)表示,组间比较使用两样本t检验或Wilcoxon秩和检验。计数资料用例数和百分率(%)表示,组间比较采用卡方检验。将单因素分析中有统计学意义的因素纳入多因素logistic回归分析,构建预测模型。采用ROC曲线评价模型的预测能力,采用Hosmer-Lemeshow拟合优度检验评价模型的校准能力。P<0.05为差异具有统计学意义。
2 结果
本研究最终纳入472例早产儿,其中非BPD组325例,BPD组147例。两组的基本情况见
表1。两组性别、胎龄、出生体重及极低出生体重儿和超低出生体重儿的比例的比较,差异均有统计学意义(
P<0.05)。模型构建组392例早产儿中,非BPD 272例,BPD 120例;模型验证组80例早产儿中,非BPD 53例,BPD 27例。
2.1 胎龄<32周早产儿BPD发生的单因素分析
单因素logistic分析显示,SNAP-Ⅱ分数越高、胎龄越小、出生体重越低、Apgar评分越低、恢复出生体重日龄越长、平均体重增长越慢,发生BPD的风险越高(
P<0.05)。母亲产前使用糖皮质激素可使BPD风险降低,痰培养阳性与BPD风险增加有关。当存在RDS、早产儿呼吸暂停(apnea of prematurity, AOP)、肺炎、肺出血、hsPDA、NEC、Ⅲ~Ⅳ级脑室内出血(intraventricular hemorrhage, IVH)、早产儿视网膜病变(retinopathy of prematurity, ROP)、低钙血症时,患儿发生BPD的风险增加(
P<0.05)。禁食时间及微量喂养时间越长、开奶时间越晚、平均加奶速度越慢,以及热量的摄入不足均与BPD的发生风险增加有关(
P<0.05)。用氧时间越长、机械通气时间越长,BPD的发生风险越高(
P<0.05)。生后第1天FiO
2越高,BPD发生的风险越大(
P<0.05)。需PS干预、静脉使用地塞米松与BPD的发生有关(
P<0.05);物理治疗开始时间越晚,BPD的发生风险越高(
P<0.05)。见
表2。
2.2 胎龄<32周早产儿不同时间点BPD发生的多因素logistic回归分析
将上述单因素分析中差异具有统计学意义的指标纳入多因素logistic回归分析。分析过程中将连续变量(禁食时间、开奶时间、微量喂养时间、平均加奶速度、恢复出生体重日龄、平均体重增长、机械通气时长、物理治疗开始日龄)转变为分类变量,分类标准根据ROC曲线及约登指数取最佳截断值。最终筛选出胎龄<32周早产儿在生后第1、7、14、21、28天5个时间点的BPD发生的独立预测因素,分别为5个、9个、7个、8个、9个(
表3)。
2.3 胎龄<32周早产儿BPD多时间点预测模型的构建
根据ROC曲线分析,每个时间点(生后第1、7、14、21、28天)各选定了5个对预测BPD贡献度最大的变量进入最终预测模型,5个时间点的最终模型ROC曲线见
图2。第1、7、14、21、28天模型预测BPD的AUC分别为0.917、0.927、0.939、0.944、0.968。
2.4 胎龄<32周早产儿BPD多时间点预测模型的列线图
通过构建多因素logistic回归模型,根据回归系数给每个影响因素的取值水平进行赋分,将各评分相加得到分数总和,通过总分与BPD发生概率之间的函数转换计算出BPD的预测值。预测患儿在生后第1、7、14、21、28天发生BPD概率的列线图见
图3。
2.5 胎龄<32周早产儿BPD多时间点预测模型的内部验证
将验证组80例早产儿(非BPD组53例,BPD组27例)的相关数据代入5个时间点的预测模型中,绘制ROC曲线进行验证。结果显示,各时间点模型在验证组中仍保持良好的判别能力,AUC值分别为0.899(第1天)、0.917(第7天)、0.908(第14天)、0.958(第21天)和0.912(第28天),见
图4。
进一步采用Hosmer-Lemeshow拟合优度检验评估模型的校准能力,显示生后第1、7、14、21、28天的
值分别为9.64、5.83、6.71、8.24、9.99;
P值分别为0.291、0.666、0.568、0.410、0.266。5个时间点的拟合优度校准曲线见
图5。
3 讨论
BPD是早产儿中一种常见且严重的肺部疾病,其病因复杂,涉及遗传、环境、围产期和新生儿期的多种因素
[3,5]。本研究探讨了胎龄<32周的早产儿发生BPD的预测因素,研究结果揭示了出生体重、胎龄、SNAP-Ⅱ评分、机械通气持续时间、FiO
2等多个影响BPD发生的关键因素,并根据这些因素构建了具有较高预测效能的多时间点预测模型。
患儿的胎龄、出生体重是BPD发生的重要影响因素
[16]。本研究结果显示,胎龄越小、出生体重越低,发生BPD的风险越高,这一结果印证了极低出生体重儿在BPD易感性中的核心地位
[16-17]。既往研究显示,在胎龄22~24周的早产儿中,有80%被诊断为BPD,而在胎龄大于28周的婴儿中仅为20%
[18]。在患有BPD的早产儿中,95%是极低出生体重儿
[19],5%是超低出生体重儿
[20]。在预测模型层面,出生体重与胎龄分别参与了5个和4个模型,进一步证明了二者在BPD发生机制中的重要作用。在疾病严重程度指标方面,SNAP-Ⅱ评分是BPD的独立危险因素(
OR=1.058~1.122),在第1天至第14天的模型中均展现了较高的贡献度,这与基于生理储备耗竭理论的预测模型构建原则
[21]相符。在呼吸支持方面,是否进行机械通气、FiO
2、有创机械通气时间及机械通气时间均是BPD发生的重要预测因素。长时间的机械通气和高浓度氧疗可能导致肺损伤,进而引发BPD,因此临床实践中应尽量减少有创通气的应用,采取保护性肺通气策略
[22]。此外,早期识别和干预呼吸支持需求,缩短机械通气时间,有助于降低BPD的发生风险
[23]。胎龄<32周早产儿的营养状况对BPD的发生有显著影响。本研究结果提示,开奶时间越晚、禁食时间越长、平均加奶速度越慢均会增加BPD的风险(
OR=2.090~3.848),营养支持情况在第7天起参与模型贡献。营养不良会导致早产儿生长发育迟缓,影响肺部修复和发育
[24]。优化早产儿的营养支持策略,对预防BPD至关重要。早期充足的营养支持不仅有助于改善肺功能,还可以减少其他并发症的发生,进一步降低BPD的风险
[25]。
本研究通过整合多时间点动态监测数据,构建了具有阶段特异性的多时间点独立的BPD临床预测模型。结果显示,各日龄节点(第1、7、14、21、28天)预测模型的AUC分别为0.917、0.927、0.939、0.944、0.968,这说明随着早产儿肺发育进程的推进及临床表型的逐步显现
[3],动态整合时序性临床数据可提升预测模型的判别能力,这与既往研究中预测模型随观察时间延长而优化的结论
[12-13]形成互证。NICHD研究团队
[12]于2011年开发的6个时间点预测模型(AUC为0.793~0.854)虽证实了预测效能的时序增强现象,但其各时间点均采用固定变量,可能忽视了不同病理阶段的特异性风险要素;Valenzuela-Stutman等
[13]于2019年建立的动态预测模型(AUC为0.788~0.894)虽实现了变量筛选的时序优化,在生后各日龄节点筛选出不同的预测因子,但其基于传统氧依赖28 d的诊断标准,未能充分反映现代围产医学进展下“新型BPD”的疾病诊疗特征
[14]。相较而言,本研究模型体系采用NICHD 2018年BPD诊断标准,通过融合多阶段不同时间点特异性预测因子,可能更具临床适用性。
涵盖多个时间点的预测模型能够在早产儿的不同发育阶段动态评估BPD的风险,这种多时间点的预测模型能够整合临床特征的时序性变化
[3],提高预测效能,从而可为临床干预提供更为精准的参考。多时间点的风险预测形式使临床医生能够依据患儿日龄选择适配的预测因子和干预策略,例如,早期模型(第1~7天)中纳入的出生体重、胎龄及SNAP-Ⅱ评分等围产期参数,可为产房内肺保护性通气策略的制订提供依据
[26];中期模型(第7~14天)中机械通气时长等呼吸支持参数的引入,提示肺损伤的累积效应开始显现
[3],为呼吸机模式选择及撤机时机的选择提供量化参考
[27];后期模型(第14~28天)中营养指标(开奶时间)及治疗措施(地塞米松使用)的贡献度提升,则与BPD病理进程中炎症修复与组织重塑的动态平衡机制相契合
[5,28]。这种预测因子的阶段特异性演变与既往研究
[3]的发现一致,可为临床实施更精准的时序干预提供理论依据。在本研究中,模型AUC随日龄增加而增加,从第1天的0.917升至第28天的0.968,这可能源于随着临床信息的累积,中后期模型可整合机械通气暴露时长、氧疗依赖性等更具判别力的时序数据从而提高模型预测效能,提示多时间点动态风险评估较传统单时间点预测具有显著优越性,这与既往Meta分析中动态模型较单时间点模型更具优势的结论
[10,14]一致。值得注意的是,针对BPD干预的关键时间窗(生后第7~14 天)
[8],本研究第7天(AUC=0.927)及第14天(AUC=0.939)模型均展现出优异的判别能力,其AUC值均高于Laughon等
[12]基于固定变量的预测模型(第7天AUC=0.818;第14天AUC=0.827)。基于上述多时间点动态模型,构建床旁可用的列线图工具,通过将多变量回归方程转化为可视化评分系统,使临床医生可以根据患儿的具体情况快速得出预测概率,从而实现个性化的风险管理。本研究内部验证结果显示5个时间点的模型对胎龄<32周早产儿是否发生BPD均具有良好判别能力,Hosmer-Lemeshow拟合优度检验亦显示模型预测概率与实际发生概率具有良好的一致性。
本研究的局限性如下:首先,单中心回顾性设计可能导致呼吸管理方案同质化对预测因子的筛选产生偏倚,需通过多中心外部验证确认模型的普适性;其次,当前纳入的预测变量以常规临床指标为主,未来需整合表观遗传标记物、肺部超声参数等新型生物标志物以提升模型的前沿性;最后,本研究样本量相对较少,扩大样本量可能揭示更多潜在预测因子,这需要通过多中心协作研究进一步探索。
总之,本研究成功构建了胎龄<32周早产儿的BPD多时间点动态预测模型,生后第1~28天的5个独立模型均展现出良好的判别能力,预测效能随日龄递增。模型揭示了BPD风险预测因素的动态演变规律,从早期围产期指标到后期治疗因素的时序性变化,可为临床提供精准干预的靶点和可视化决策工具。
中南大学湘雅早产儿临床大数据库系统建设(056)