脓毒症是由感染引发的,以机体免疫反应失调为特征的一类疾病。脓毒症相关死亡约占全球总死亡人数的20%
[1, 2]。早期识别高危患者并实施针对性干预,对于降低病死率具有重要的临床意义。目前,降钙素原(PCT)、C反应蛋白(CRP)和中性粒细胞表面抗原CD64等生物标志物已被证实可预测脓毒症早期预后
[3, 4],但其在临床应用中仍存在局限性。如PCT在非感染性炎症时亦可升高,特异性不足;CRP的敏感性和特异性均较差,难以准确区分脓毒症与其他炎症疾病;CD64虽敏感性和特异性均较高,但检测成本高、需特殊设备,在临床难以广泛应用
[5, 6]。因此,寻找一种新生物标志物,既具有较好的敏感性和特异性,又易于临床推广,将有助于早期预测脓毒症相关死亡。
近年来,红细胞分布宽度 (RDW)的升高被证实是脓毒症死亡的独立危险因素
[7-9]。RDW反映红细胞体积异质性
[10],其升高提示红细胞稳态失调,可能与端粒缩短、氧化应激、炎症、以及促红细胞生成素异常导致的生成受损或代谢与存活异常有关。同时,红细胞压积 (HCT) 作为浓缩红细胞相对于全血红细胞体积的指标,既往研究显示低HCT是脓毒症30 d内死亡率增加的独立危险因素
[11]。然而,RDW或HCT水平易受营养状况、液体复苏以及红细胞输注等因素的干扰
[12-14]。
既往研究发现复合指标在预测脓毒症预后较单一指标稳定性更好。如一项研究发现乳酸/血细胞比值可预测脓毒症死亡
[15],但未对混杂因素进行校正。也有研究发现中性粒细胞/淋巴细胞比值以及淋巴细胞/单核细胞比值可用于预测脓毒症预后,但研究样本量小,仅47例
[16]。相比之下,RDW/HCT可能更能综合反映生理状况,降低单项指标的个体变异干扰。原因在于:脓毒症相关的全身炎症反应和氧化应激可抑制骨髓造血功能、缩短红细胞寿命,促进未成熟红细胞进入外周循环,导致RDW升高
[8];另一方面血液稀释、贫血以及红细胞破坏增加可致HCT下降,降低机体携氧能力,加重组织缺氧和多器官功能损伤,影响患者预后
[11]。基于此,RDW/HCT 比值的变化可综合反映脓毒症患者红细胞生成、炎症反应和氧化应激等病理生理过程。然而,目前仅有慢性阻塞性肺部疾病急性加重期的队列研究报道RDW/HCT升高可预测28 d全因死亡
[17],该比值在预测脓毒症预后中的作用尚不明确。
据此,本研究拟基于美国重症监护医学信息数据IV 3.1(MIMIC-IV)数据库,采用回顾性队列研究的方法,利用Kaplan-Meier(K-M)生存曲线、多变量Cox回归和限制性立方样条回归(RCS)分析RDW/HCT比值与脓毒症患者28 d全因死亡率的相关性,并通过亚组分析评估结果的稳健性。同时,基于电子重症监护病房协作研究数据库2.0(eICU-CRD)进行外部验证,以期为脓毒症早期死亡识别提供简单有效的评估指标。
1 资料和方法
1.1 数据来源
本研究采用回顾性队列研究方法,从MIMIC-IV 3.1纳入27 335例成年脓毒症患者。验证队列为电子重症监护病房协作研究数据库eICU-CRD,该数据库收集了2014~2015年来自美国208家医院的患者信息,且与MIMIC-IV由不同医院组成,彼此相互独立。本研究从eICU-CRD 2.0共纳入17 406例成年脓毒症患者。两数据库均已对患者信息进行去隐私化处理,符合医学伦理学标准。研究已获得MIMIC-IV 3.1及eICU-CRD 2.0数据库使用授权 (证书编号:63806772)。纳入标准 :符合脓毒症3.0诊断标准且首次入重症监护病房(ICU)的患者。排除标准:年龄<18岁;入住ICU(<24 h);RDW或HCT缺失。
1.2 数据收集
使用PostgresSQL(v16.3.2)和Navicate Premium(v16)提取以下数据:一般资料:年龄、性别、种族;合并症:心梗、充血性心力衰竭、脑血管疾病、恶性肿瘤、糖尿病、慢性肾脏疾病、急性肾损伤;入住ICU 24 h之内的实验室化验结果:RDW、HCT、血小板(PLT)、红细胞(RBC)、白细胞(WBC)、凝血酶原时间(PT)、部分促凝血酶原时间(PTT)、国际标准化比值(INR)、血清肌酐(Scr)、血尿素氮(BUN)、血糖以及电解质(钾、钠);生命体征:心率、平均动脉压、呼吸频率;临床治疗措施:血管活性药物、激素类药物、连续性肾脏替代治疗(CRRT),机械通气;入住ICU 24 h之内的评估:查尔森合并症指数(CCI)、序贯器官衰竭评分(SOFA)、格拉斯哥昏迷评分法(GCS)。
1.3 结局指标
28 d全因死亡率:定义为从入ICU当日起28 d内的任何原因死亡。
1.4 统计学分析
本研究采用R 4.3.3 软件行统计分析。连续变量的正态性采用 Shapiro-Wilk检验进行评估,并结合直方图及Q–Q图对判断数据分布,符合正态分布的连续性变量用均数±标准差表示,不符合正态分布的连续性变量用中位数(四分位数间距)表示。分类变量以频率(百分比)表示。基线资料组间比较采用单因素方差分析、K-W检验和卡方检验。对缺失比例小于20%的变量采用多重插补法,通过多重插补链式方程(MICE)生成5个插补数据集,经一致性评价后,选取其中一套插补数据进行后续分析。本研究遵循加强流行病学中观察性研究报告质量声明(STROBE)声明, 构建多因素Cox比例风险模型分析RDW/HCT比值对预后的独立影响。模型1:未校正任何混杂因素;模型2: 校正年龄、性别、种族;MIMIC-IV的模型3:在模型2基础上加入CCI、SOFA、GCS、RBC、WBC、PT、PTT、INR、Scr、BUN、血糖、血钠、血钾、既往病史(心梗、脑血管疾病、充血性心力衰竭、恶性肿瘤、慢性肾脏疾病、急性肾损伤)、心率、平均动脉压、呼吸频率、血管活性药物、激素类药物、CRRT;eICU-CRD的模型3校正因素同MIMIC-IV的模型3的校正因素,除外糖尿病和CRRT。
通过RCS分析,设置4个节点(第5、35、65和95百分位),构建样条函数,拟合RDW/HCT与死亡率的关系。若存在非线性关联,通过递归法及最大似然法确定拐点,并利用似然比检验对比分段Cox回归和普通Cox回归模型的拟合优度,检验是否存在阈值或饱和效应,P<0.05说明存在显著效应,两者关系为非线性关系。
利用亚组分析探讨不同亚组中连续RDW/HCT与28 d死亡率之间的关系。亚组分析采用分层Cox回归,校正混杂因素同模型3,使用似然比检验进行交互作用评价,并绘制亚组森林图。所有检验均为双侧检验,P<0.05被认为差异具有统计学意义。
2 结果
2.1 纳入人群总体特征
本研究纳入MIMIC-IV 3.1中的27 335例脓毒症患者,中位年龄为68岁,其中58.1%为男性;28 d内死亡率为19.3% (n=5276)。RDW/HCT中位数为0.46 (IQR: 0.38~0.58)。根据RDW/HCT比值的四分位数将患者分为4组:Q1(RDW/HCT≤0.384;n=6834)、Q2(0.384<RDW/HCT≤0.463;n=6834)、Q3(0.463<RDW/HCT≤0.576;n=6833)和Q4(RDW/HCT>0.576;n=6834)。4组患者28 d内的死亡率分别是16.8%、16.6%、18.0%、25.9%。
脓毒症患者基线特征分析显示(
表1)。RDW/HCT比值较高组死亡率高于较低组,差异有统计学意义(
P<0.001)。与RDW/HCT比值较低组相比,比值较高组患者危重程度更重(SOFA评分,GCS评分),合并心衰、恶性肿瘤、糖尿病、肾病的比例较高,PT、PTT、INR、Scr、BUN较高,但HCT、RBC、WBC、血糖水平较低。
2.2 K-M生存曲线分析
K-M生存曲线(
图1)显示,4组脓毒症患者的28 d生存率有统计学差异(
P<0.001)。与RDW/HCT比值较低组相比(Q1、Q2),比值较高组(Q3、Q4)患者的生存率显著降低,其中Q4组生存率最低.
2.3 多因素Cox回归分析
根据28 d预后将患者分为死亡组与未死亡组。与未死亡组患者相比,死亡组患者入ICU时SOFA评分更高,RDW/HCT更高(
表2)。两组在年龄、性别、种族、疾病危重程度(SOFA评分,GCS评分)、RDW、RBC、WBC、PT、PTT、INR、Scr、BUN、血糖、钾、钠、合并心梗、脑血管疾病、充血性心力衰竭、恶性肿瘤、肾脏疾病等疾病的比例,以及心率、呼吸频率、血压、使用血管活性药物、激素类药物、CRRT的比例均有显著差异(
P<0.05,
表2)。
Cox比例风险回归分析显示(
表3),RDW/HCT作为连续变量与28 d死亡风险呈正相关。模型1中,HR=3.29,95%
CI:2.87~3.78,
P<0.001;模型2中,HR=3.39,95%
CI:2.94-3.91,
P<0.001;模型3中,HR=1.89,95%
CI:1.52~2.36,
P<0.001。进一步以Q2组作为参照组。结果显示,Q4组在所有模型中均显示与28 d死亡风险显著正相关(
P<0.001)。其中,模型3显示,与Q2组相比,Q4组死亡风险增加31%(HR=1.31,95%
CI:1.19~1.44,
P<0.001),但Q1与Q3组与Q2组相比,死亡风险无统计学意义。趋势性检验提示RDW/HCT与28 d死亡率之间可能存在非线性关系(
P<0.001)。
2.4 RCS分析和饱和阈值效应分析
进一步行RCS分析,显示RDW/HCT比值和患者28 d死亡风险呈非线性“U型”关系(
P<0.001,
图2)。阈值效应分析显示,28 d死亡风险的拐点为0.35(
表4)。RDW/HCT<0.35,比值与28 d死亡率呈显著负相关(HR=0.022,95%
CI:0.003~0.162,
P<0.001);RDW/HCT≥0.35,28 d死亡风险随比值升高而显著增高(HR=2.034,95%
CI:1.628~2.541,
P<0.001)。进一步比较标准线性模型与分段线性模型的拟合优度,对数似然比检验证实分段线性模型显著优于标准线性模型(对数似然比
P<0.001),提示阈值效应的存在。
2.5 RDW/HCT的亚组分析结果
进一步分析不同亚组之间RDW/HCT与脓毒症患者28 d死亡率之间的关系。根据性别、年龄、心肌梗死、充血性心力衰竭、慢性肾脏疾病、恶性肿瘤、是否使用激素类药物进行亚组分析,结果显示:RDW/HCT在各亚组之间均无显著交互作用(P>0.05),RDW/HCT与脓毒症患者28 d死亡率之间的关系具有稳健性。
2.6 eICU-CRD外部验证
2.6.1 eICU-CRD人群总体特征
本研究在eICU-CRD数据库纳入17 406例脓毒症患者,中位年龄67岁,48.1%为男性;28 d内死亡率为16.0%(2785例)。RDW/HCT的中位数为0.48 (IQR:0.40~0.58)。根据RDW/HCT比值4分位数分组: Q1(RDW/HCT≤0.400;
n=3813)、Q2(0.400<RDW/HCT≤0.480;
n=4551)、Q3(0.480<RDW/HCT≤0.580;
n=4499)和Q4 (RDW/HCT>0.580;
n=4543)。四组患者28 d内的死亡率分别是11.1%、13.2%、16.4%和22.5%。
表5列出了4组脓毒症患者的基线特征。RDW/HCT比值较高组死亡率显著高于比值较低组(
P<0.001)。与RDW/HCT比值较低组相比,比值较高组的患者SOFA评分更高,合并恶性肿瘤及慢性肾脏疾病的比例较高,RDW、WBC、PT、PTT、INR、Scr、BUN、血钾浓度较高,但HCT、RBC、血糖水平较低。
2.6.2 Cox分析
进一步根据患者是否死亡分为两组。与未死亡组患者相比,死亡组入住ICU时SOFA评分及RDW/HCT比值更高(
表6)。两组年龄、种族、危重程度(SOFA评分,GCS评分)、RDW、HCT、RBC、WBC、PT、PTT、INR、Scr、BUN、钾、钠,合并心梗、脑血管疾病、心力衰竭、恶性肿瘤、肾脏疾病等疾病的比例,以及心率、呼吸频率、血压、血管活性药物、糖皮质激素的比例均有显著差异(均
P<0.05)。
在eICU-CRD队列(
表3),模型1显示:RDW/HCT作为连续变量与28 d死亡率呈显著正相关(HR=3.10, 95%
CI:2.45~3.94,
P<0.001),但在模型2和模型3中,多变量充分校正后,两者无显著相关性。
根据RDW/HCT比值分为4组,模型3发现:与Q2组相比,Q1组28 d死亡风险显著升高 (HR=1.23,95% CI:1.08~1.40,P=0.002),Q3组风险显著降低(HR=0.88, 95% CI:0.79~0.99,P=0.030),Q4组则无显著性,提示RDW/HCT与28 d死亡率之间可能存在非线性关系。
2.6.3 K-M生存曲线分析
K-M生存曲线(
图4)显示,RDW/HCT比值与28 d的生存率差异有统计学意义(
P<0.001)。低水平组(Q1、Q2)生存率高于高水平组(Q3、Q4)。
2.6.4 RCS分析和饱和阈值效应分析
RCS分析显示RDW/HCT比值与28 d死亡率存在非线性关系(
P<0.001),其趋势与MIMIC-IV队列一致,均呈“U型”关系。饱和阈值效应分析显示死亡风险拐点为0.48(
表4),RDW/HCT<0.48时,比值与28 d死亡率呈显著负相关 (HR=0.031, 95%
CI:0.011~0.089,
P<0.001),RDW/HCT≥0.48,28 d死亡风险随比值升高而显著增加 (HR=1.969,95%
CI:1.402~2.766,
P<0.001),且分段线性模型显著优于标准线性模型(对数似然比
P<0.001)。
3 讨论
本研究基于MIMIC-IV和eICU-CRD数据库的大样本分析发现,RDW/HCT比值与脓毒症患者28 d全因死亡呈非线性“U型”关系。且均存在显著的阈值效应。这一结果在不同数据库中表现出良好的稳健性。尽管两数据库的风险拐点存在差异(MIMIC-IV为0.35,eICU-CRD为0.48),这可能与不同来源患者的异质性、实验室检测方法的差异以及人群分布特征有关,但整体非线性趋势的高度一致性提示RDW/HCT作为脓毒症预后评估指标的临床价值。
相对于单独使用RDW或HCT,RDW/HCT比值综合了红细胞数量和红细胞体积的变化,能更全面地反映红细胞在生成与功能方面的双重异常。既往研究关于RDW/HCT比值与疾病预后的相关性研究较少。Li 等
[18]研究发现RDW能较好预测危重心血管病人的全因死亡率。两者存在非线性关系,表现为死亡率随着RDW的降低而下降;直至拐点,死亡率随着RDW的升高而上升,随后进入平台期。也有研究发现HCT与急性缺血性脑卒中男性患者住院期间死亡率存在“U型”相关,过高的HCT值与死亡风险升高相关
[19]。本研究为避免RDW或HCT受患者自身的慢性疾病、营养状况和输血补液的影响,采用RDW/HCT,发现在MIMIC-IV脓毒症人群中,当RDW/HCT比值≥0.35,RDW/HCT比值每增加1,脓毒症患者的28 d死亡风险增加2.034倍,这提示该复合指标识别高死亡风险的脓毒症人群具有较高的敏感性。
然而,RDW/HCT比值增高与脓毒症患者死亡风险增加的机制尚不清楚,可能与以下几个因素有关:氧供与微循环障碍:低HCT意味着血液中红细胞数量减少,血液携氧能力下降。Luo等
[11]研究显示,与常规HCT组相比,低HCT组脓毒症患者的30 d死亡率增加了1.589倍(HR=1.589,95%
CI:1.009~2.979)。在此基础上,RDW升高表明红细胞变异性增大,可能增加微循环中的血流阻力,影响组织灌注,从而加重器官缺氧与缺血性损伤
[20]。红细胞稳态受损:RDW升高与红细胞生成紊乱相关,且与心血管疾病、慢性肺部疾病、急性脑梗死等多种疾病的预后有关
[21-24]。红细胞数量与质量的异常可能反映了脓毒症引起的红细胞生成功能的衰退与寿命的缩短。研究表明,RDW对脓毒症死亡率的预测能力与SOFA评分和APACHE-II评分相当
[25]。血流动力学不稳定:RDW升高和HCT下降可能反映了脓毒症患者血容量变化或贫血加重,引发血压不稳定、低血压和心输出量下降等,增加脓毒症死亡风险
[26, 27]。炎症与氧化应激叠加:高RDW伴随炎症级联反应的放大与氧化应激的加剧,产生的活性氧还可致红细胞变形性降低;若HCT同时降低,将进一步加重机体缺氧状况,加速重要脏器衰竭
[13, 28-30]。
本研究存在一定局限性。首先,作为回顾性研究,无法排除选择偏倚和数据收集过程中潜在的误差。其次,MIMIC-IV数据库与eICU-CRD的RDW/HCT的拐点存在差异,但整体趋势的一致性表明该非线性关系具有稳健性,未来需在更大规模、多中心及更具一致性的人群中开展前瞻性临床研究,明确其临床指导意义。
综上所述,RDW/HCT比值升高与脓毒症患者28 d全因死亡率有关。作为一种简便、高效的生物标志物,该指标有助于脓毒症患者早期死亡风险分层及预后评估。
国家自然科学基金(82002084)
国家自然科学基金(81971806)
四川省科技计划项目(2024NSFC1535)
四川省科技计划项目(2023YFH0074)