支气管肺发育不良(bronchopulmonary dysplasia, BPD)是一种常见于早产儿的慢性肺部疾病。BPD的主要病理特征包括肺泡发育不全、肺微血管发育障碍,以及肺泡数量和肺毛细血管床的减少,易导致不可逆的气道上皮损伤,即使经过充分治疗,仍可能对患儿的长期生长和发育产生不利影响
[1]。鉴于BPD的复杂病理过程及其深远影响,探索其潜在机制及治疗靶点尤为重要。微小核糖核酸(microRNA, miRNA)通过促进mRNA降解或抑制其翻译来调控基因表达,参与细胞增殖、分化和血管生成等多种生物学过程,并与肺部疾病的发生和进展密切相关。已有研究表明,miRNAs-15b(miR-15b)能够调节血管细胞的增殖、分化及凋亡,其过表达显著抑制新生血管和侧支动脉的形成
[2]。血管内皮生长因子(vascular endothelial growth factor, VEGF)作为一种重要的促血管生成因子,在肺部发育和修复过程中发挥着关键作用。研究发现,VEGF水平可作为新生儿肺损伤的生物标志物
[3]。尽管VEGF在BPD发病机制中的作用已有一定的研究基础
[4],但其与miR-15b的相互作用及其在BPD早产儿神经系统发育中的具体机制仍需进一步探讨。本研究旨在探索血清miR-15b和VEGF在BPD早产儿中的动态变化及其与近期神经系统发育的关系,这对于理解BPD的病理机制及寻找新的治疗策略具有重要意义。
1 资料与方法
1.1 临床资料
回顾性收集2020年1月—2023年2月我院新生儿重症监护病房(neonatal intensive care unit, NICU)收治的BPD早产儿170例为研究对象。纳入标准:(1)本院出生且出生后即转入NICU;(2)28周≤胎龄≤32周;(3)住院时间>28 d;(4)出生体重<1 500 g;(5)在研究数据收集及分析时,已确诊符合NICHD于2001年制定的BPD诊断及分度标准
[5]。排除标准:(1)新生儿重度窒息;(2)合并中枢神经系统感染;(3)颅内出血程度>Ⅲ度;(4)合并先天发育畸形。研究期间自动放弃3例,死亡1例,数据不完整者10例,最终纳入156例。本研究通过我院伦理委员会批准(批准号:2024-034),且获得患儿监护人的知情同意。
1.2 BPD诊断标准及分组
依据NICHD于2001年制定的BPD诊断标准
[6]进行诊断及分组。任何氧依赖(氧浓度>21%)≥28 d的新生儿被诊断为BPD。
对于出生胎龄<32周的BPD患儿,在校正胎龄达到36周或出院时,依据需氧浓度进行BPD分度:重度(氧浓度≥30%或机械通气)、中度(氧浓度<30%)和轻度(未用氧)。依据BPD分度,轻度88例,中度55例,重度13例,即轻度组88例与中重度组68例。
1.3 观察指标
通过查阅病历获取患儿出生后第1、7、14、28天静脉血检测结果,选择这些时间点检测血清miR-5b和VEGF的原因如下。生后早期是BPD关键期,第1天检测可获得基础血清标志物水平,为后续动态监测提供起始参考,后续时间点能捕捉到BPD发展过程中肺泡和肺血管系统损伤修复、炎症反应等变化对血清标志物水平的影响,有助于了解疾病进展。对于神经系统发育监测,早产儿前28 d神经发育关键事件频发,这些时间点的检测可观察标志物与神经发育的关联。虽然BPD诊断以持续通气28 d为标准,但血清最后采集时间(第28天)可从血清学角度反映此时机体状态,与BPD诊断时间互补,共同为评估患儿病情及神经系统发育提供依据。
神经发育评估(developmental screening test, DST)
[7]:所有受试早产儿在矫正月龄达到6月龄时,于我院儿童保健门诊由接受过专业培训的医生进行DST。DST主要用于0~72月龄儿童的评估。DST测试结果以智力指数(mental index, MI)和发育商(development quotient, DQ)表示。MI反映智力能区的发育水平,而DQ反映3个能区(社会适应能力、运动能力和智力能力)的总体发育水平。对于3岁以下儿童,根据DQ将其分为3个等级:DQ<70分为异常,70分≤DQ<85分为可疑,DQ≥85分为正常。本研究将DQ<70分的患儿纳入DST异常组,而未将DQ介于70分到85分之间的患儿(即DST可疑患者)纳入本研究,且将DQ≥85分的患儿纳入DST正常组。
1.4 统计学分析
通过SPSS 22.0软件对数据进行统计学分析,连续变量的正态性检验采用Shapiro-Wilk法,方差齐性检验采用Levene检验。正态分布的计量资料用均数±标准差()表示,两组间比较采用两样本t检验(方差齐时)或校正t检验(方差不齐时);非正态分布计量资料以中位数(四分位数间距)[M(P25,P75)]表示,组间比较采用Mann-Whitney U检验。计数资料以例数和百分率(%)表示,组间比较采用检验或Fisher确切概率法。使用重复测量方差分析探讨血清miR-15b、VEGF水平随时间变化的趋势及组间差异。为保证方差分析的有效性并控制Ⅰ类错误率,进行Mauchly's球形假设(Sphericity)检验;当时间×组间交互效应显著时,在各时间点上进行事后两两比较。使用logistic回归逐步校正多个潜在的混杂因素,分析影响神经系统发育不良的临床因素。使用TargetScanHuman 7.1预测miR-15b和VEGF的结合位点。使用逐步回归法和Bootstrap法检验分析VEGF在miR-15b和近期神经系统发育之间的中介效应。P<0.05为差异具有统计学意义。
2 结果
2.1 BPD早产儿出生情况及母亲孕期情况比较
与轻度组比较,中重度组早产儿小于胎龄儿比例、1 min Apgar评分≤7分比例、出生2周仍需机械通气比例、机械通气时间、总氧疗时间、抗生素使用时间、使用枸橼酸咖啡因比例、吸氧浓度、动脉血二氧化碳分压明显较高,母亲的先兆早产比例、重度子痫前期比例明显较高,早产儿胎龄、出生体重、身长、头围、动脉血pH、动脉血氧分压明显较低,差异均有统计学意义(均
P<0.05)。见
表1。
2.2 BPD早产儿血清miR-15b、VEGF动态变化
miR-15b和VEGF数据均符合球形假设(miR-15b:=6.342,df=5,P=0.273;VEGF:=8.216,df=5,P=0.145)。
重复测量方差分析显示,两组miR-15b和VEGF均存在组间差异、时间差异(
P<0.05),且时间与组间均有交叉效应(
P<0.05)。进一步行两因素的简单效应分析显示,出生第1天时,两组间血清miR-15b和VEGF水平差异无统计学意义(
P>0.05)。第7~28天,与轻度组相比,中重度组的miR-15b水平明显升高,而VEGF水平明显降低(
P<0.05)。轻度组的miR-15b水平随时间逐渐下降,VEGF水平上升;相反,中重度组的miR-15b水平随时间逐渐上升,VEGF水平下降;事后两两比较显示,差异具有统计学意义(
P<0.05)。见
表2、
图1。
2.3 BPD早产儿近期神经系统发育不良的临床因素分析
在不同严重程度BPD早产儿中,DST异常组的1 min Apgar评分≤7分比例、动脉血二氧化碳分压水平、机械通气时间、总氧疗时间、吸氧浓度、抗生素使用时间、使用枸橼酸咖啡因比例均明显高于DST正常组。此外,中重度组中,DST异常组的PaO
2水平明显低于DST正常组,出生2周仍需机械通气比例、先兆早产比例明显高于DST正常组(
P<0.05);轻度组中,DST异常组的重度子痫前期比例明显高于DST正常组(
P<0.05)。见
表3。
2.4 不同严重程度BPD早产儿血清miR-15b、VEGF与近期神经系统发育不良的关系分析
分别以不同严重程度BPD早产儿出生后第1~28天的miR-15b、VEGF水平为自变量,以近期神经系统发育情况(DST正常=0,DST异常=1)为因变量,并选择轻度组为参考组,进行logistic回归分析,模型1为调整1 min Apgar评分≤7分比例、PaCO
2水平、机械通气时间、总氧疗时间、吸氧浓度、抗生素使用时间;模型2在模型1的基础上调整使用枸橼酸咖啡因比例、动脉血氧分压水平及出生2周仍需机械通气、先兆早产、重度子痫前期比例等指标。结果显示,在未校正模型及逐步校正混杂因素后的模型1和模型2中,从第7天开始,与轻度组相比,其
OR值逐渐降低,提示VEGF水平越低,DST异常的风险越高(
P<0.001)。相反,miR-15b水平在第7天后,中重度组的
OR值逐渐增加,提示miR-15b水平越高,DST异常的风险越高(
P<0.001)。见
表4。
2.5 生物信息学预测
预测软件Target Scan Human 7.1预测miR-15b的靶基因,结果显示,miR-15b在VEGF 3'UTR 276-283bp间有1个保守结合位点,见
图3。提示miR-15b通过直接结合VEGF的3'UTR调控其表达。
2.6 VEGF对miR-15b和近期神经系统发育之间关系的中介作用模型的回归分析
分别统计早产儿出生后第1~28天的VEGF、miR-15b整体平均水平,并使用逐步回归法和Bootstrap法检验分析VEGF在miR-15b和近期神经系统发育之间的中介效应。模型1分析miR-15b与近期神经系统发育的关系,结果显示,调整混杂因素前后,miR-15b对近期神经系统发育均有明显负向影响(未调整:
β=-1.253,
t=-4.391,
P<0.001;调整后:
β=-1.230,
t=-4.231,
P<0.001),表明miR-15b直接抑制近期神经系统发育。模型2分析miR-15b与VEGF之间的关系,结果显示,调整混杂因素前后,miR-15b对VEGF有显著的负向影响(未调整:
β=-5.556,
t=-4.763,
P<0.001;调整后:
β=-5.452,
t=-4.633,
P<0.001),即miR-15b抑制VEGF表达;模型3同时考虑miR-15b和VEGF时发现,VEGF对近期神经系统发育有显著正向影响(未调整:
β=0.131,
t=6.269,
P<0.001;调整后:
β=0.131,
t=5.810,
P<0.001),而miR-15b的直接影响减弱(未调整:
β=-0.497,
t=-1.945,
P>0.05;调整后:
β=-0.592,
t=-1.984,
P>0.05)。见
表5。
进一步使用Bootstrap法(重复抽样5 000次)分析VEGF的中介效应,结果显示miR-15b通过抑制VEGF表达间接影响近期神经系统发育,间接效应值为-0.705(95%CI:-1.178~-0.372),表明这种间接路径具有统计学意义(P<0.05)。miR-15b对近期神经系统发育的总效应为-1.297(P<0.001),其中间接效应比例约为54.36%,说明miR-15b主要通过调节VEGF水平来影响早产儿的近期神经系统发育。
3 讨论
随着NICU中保护性肺通气策略的广泛应用及新型治疗手段的实施,新型BPD发病率呈上升趋势,其病理特征是肺泡及肺血管系统的发育受阻。本研究选择胎龄28~32周早产儿,原因是该胎龄段早产儿处于肺发育和神经系统发育的关键时期,且对生理病理因素敏感,BPD发生风险较高,有助于深入研究血清miR-15b、VEGF在BPD进程及其与近期神经系统发育的关系。VEGF为重要的促血管生成因子,在肺内主要由肺泡Ⅱ型上皮细胞于肺发育的肺泡期分泌,该时期分泌量相对少,影响肺泡发育。研究表明,孕妇出现绒毛膜羊膜炎、子痫前期等产前应激状况时,羊水及血液中VEGF拮抗剂水平上升,早产儿BPD风险增加,提示VEGF参与BPD的发生
[8]。越来越多的研究发现新型BPD发病机制中miRNA网络的复杂性,其中miR-15b在新型BPD的病理生理过程中起到一定作用。Xing等
[9]选取出生后12 h内的Wistar幼鼠,分为BPD组(暴露于高浓度氧60%~80%)和对照组(正常空气中);与对照组相比,BPD组中有19种miRNA,其中miR-15b明显上调。Sahni等
[10]采用高通量测序技术检测高氧诱导BPD小鼠模型中的非编码RNA,发现miR-15b表达上调。本研究发现,从BPD患儿出生第1天至第28天,轻度BPD早产儿血清VEGF水平呈明显上升趋势,而miR-15b则逐渐下降;而中重度BPD早产儿则相反。以上结果提示,不同严重程度BPD患儿血清VEGF及miR-15b的动态变化不同。在轻度BPD中,较低的miR-15b水平减弱了对VEGF表达的抑制作用,促进VEGF介导的修复过程,这可能是机体应对轻微损伤的自我保护机制。中重度BPD因肺部损伤重、炎症强,导致释放大量细胞因子,干扰VEGF的合成并加速其降解,且激活一些信号通路上调miR-15b,形成恶性循环,加重肺损伤。本研究通过预测软件Target Scan Human 7.1预测miR-15b的靶基因发现,miR-15b在VEGF 3'UTR 276-283bp间有1个保守结合位点,提示miR-15b可通过直接结合VEGF的3`UTR调控其表达,但缺乏具体的功能实验验证。已有研究表明,miR-15b可直接结合VEGF mRNA的3`UTR区域,抑制其翻译或促进mRNA降解,从而降低VEGF蛋白水平
[11]。
在BPD早产儿中,神经胶质细胞的分化及髓鞘化过程对外部因素高度敏感,长期氧疗可能导致神经行为发育迟缓乃至脑瘫等严重后遗症。VEGF作为神经发育的关键因子,可促进血管生成、神经细胞增殖、迁移及突触形成。有研究报道,miR-15b在阿尔茨海默病、帕金森病等神经系统疾病中表现出显著的调控作用
[12]。也有研究报道,miR-15b能负向调控VEGF的表达,影响神经胶质细胞的增殖、分化和髓鞘化
[13]。在BPD早产儿中,长期氧疗引发的氧化应激和炎症反应会打破miR-15b与VEGF的平衡,影响神经系统正常发育。有研究发现,BPD早产儿血清中miR-15b和VEGF水平异常变化与神经行为发育迟缓、认知功能障碍及脑性瘫痪等神经后遗症显著相关
[14]。本研究发现,两者是近期神经系统发育的独立预测因子。有研究采用体外实验模型(人胚肺成纤维细胞)
[15]及动物模型(BPD大鼠)
[16]发现,高氧条件下miR-15b水平上升,抑制VEGF生成。本研究进一步揭示了miR-15b通过VEGF间接影响早产儿近期神经系统发育的机制。结果显示,miR-15b不仅直接抑制近期神经系统发育(调整后:
β=-0.592,
t=-1.984,
P>0.05),尽管其直接影响有所减弱且不再具有统计显著性;还通过抑制VEGF的表达间接影响该过程[间接效应值为-0.705(95%
CI:-1.178~-0.372)]。miR-15b对近期神经系统发育的总效应为-1.297(95%
CI:-1.750~-0.844,
P<0.001),其中间接效应比例约为54.36%,这表明miR-15b主要通过调节VEGF水平来影响早产儿的近期神经系统发育。这与Yang等
[17]的研究结果一致。这提示VEGF在神经发育过程中扮演着关键角色,参与神经元迁移、分化和血管生成等过程。miR-15b通过抑制VEGF的表达,干扰这些过程,从而影响早期神经系统的正常发育。此外,该结果也为临床探索miR-15b和VEGF作为早产儿神经保护治疗的潜在靶点提供了理论依据。
在临床应用价值方面,本研究中血清miR-15b、VEGF水平均是DST的重要影响因素,存在独立相关性,趋势性检验存在统计学意义(P<0.001),提示二者在评估BPD早产儿近期神经系统发育方面具有潜在价值。在临床可操作性上,血清miR-15b和VEGF检测可在早产儿出生后的特定时间点(如第1、7、14、28天)进行常规采血检测,检测方法相对成熟,具有一定的可操作性。其可能的临床应用场景包括:在轻度BPD早产儿中,若发现miR-15b表达下降而VEGF表达上升,可早期提示患儿神经系统发育可能处于较好状态,但仍需密切监测,以便及时发现潜在问题并进行干预;而对于中重度BPD早产儿,若miR-15b持续升高、VEGF持续降低,可能预示更高的神经系统发育不良风险,此时可考虑采取更积极的干预措施,如神经保护治疗或加强神经系统发育监测等,从而为改善BPD早产儿的神经系统预后提供依据。
综上所述,BPD早产儿中,轻度与中重度患儿血清miR-15b、VEGF的变化趋势不同;miR-15b通过抑制VEGF表达间接影响近期神经系统的发育;miR-15b和VEGF在不同严重程度的BPD早产儿中具有不同的临床价值,特别是在轻度BPD中,miR-15b表达下降而VEGF表达上升,对DST的早期诊断和干预具有重要意义。本研究存在一定不足:(1)本研究为单中心研究,数据来源集中于一个特定的医疗单位,可能存在地域或机构特有的临床实践模式偏差;(2)本研究样本量有限,特别是中重度BPD患儿数量较少,可能限制了统计效力,导致某些潜在关联未能充分显现;(3)仅关注了短期内的神经系统发育评估,对于长期神经发育的影响尚不清楚;(4)未纳入非BPD早产儿作为对照组,无法明确这些生物标志物是否特异性地反映BPD病理过程;(5)未使用萤光素酶报告基因实验验证miR-15b与VEGF 3`UTR的结合能力,以及通过敲低或过表达miR-15b来观察VEGF表达的变化。未来本课题组将扩大样本量,涵盖更多中重度BPD患儿;通过多中心合作收集更广泛的数据,减少偏倚;进行长期跟踪随访;并纳入非BPD早产儿作为对照组,确保研究结论的可靠性和推广性;进行萤光素酶报告实验对VEGF和miR-15b进行验证。
山东省医药卫生科技发展计划项目(202406011275)