肺炎支原体（Mycoplasma pneumoniae, MP）是儿童社区获得性肺炎的重要病原之一^[1]。据统计，在住院儿童社区获得性肺炎中，肺炎支原体肺炎（Mycoplasma pneumoniae pneumonia, MPP）约占10%~40%^[2-4]。由于MP耐药、免疫功能紊乱、混合感染和MP载体量过多等原因，导致难治性肺炎支原体肺炎（refractory Mycoplasma pneumoniae pneumonia, RMPP）的发病率逐年升高^[5]。Tamura等^[6]将经大环内酯类抗菌药物正规治疗7 d及以上，临床征象加重、持续发热、肺部影像学加重者，定义为RMPP。若不及时治疗，易留有肺不张、肺实变、胸腔积液等并发症，但目前对于RMPP的早期诊断尚无有效监测指标^[7]。

乳酸脱氢酶（lactate dehydrogenase, LDH）是糖酵解途径中一种重要的酶，广泛存在于组织细胞中，其中以肝脏中活性最高，其次为心脏、骨骼肌、肾脏。即便只有很少的肺组织破坏所释放的酶就能引起血清中LDH浓度变化，所以当肺组织缺氧坏死时，细胞膜通透性增加，酶从细胞内释放至血液，使血液中LDH升高^[8]。肺炎时由于肺部炎症及缺氧，血清LDH升高，其水平变化可以反映疾病严重程度，是一项重要的血清学指标^[9]。近年来，许多文献报道LDH在RMPP患者中明显升高，通过普通logistic回归分析发现LDH是早期识别RMPP的一个重要指标^[10-11]。但目前临床中LDH值和RMPP的发生并非为简单的线性关系，单纯通过logistic回归分析求得比值比（OR值），并不能精确地评估LDH值的变化与RMPP发生的关系，因此本研究通过限制性立方样条函数探究LDH连续性变化与RMPP的关联强度，并采用决策曲线分析来评估其临床预测价值。

1 资料与方法 1.1 研究对象

以2013年2月至2019年4月在重庆医科大学附属儿童医院住院的3 764例MPP患儿为研究对象。纳入标准：符合MPP的诊断依据，即临床上有肺炎的表现和/或影像学改变，结合MP病原学检查即可诊断为MPP^[12]。MP病原学检查：单次MP抗体滴度≥1 : 160可作为诊断MP近期感染或急性感染的参考；核酸诊断技术特异性强、敏感、快速，可用于早期诊断^[13-14]。排除标准:既往有慢性肺疾病、免疫缺陷病、结缔组织疾病、血液系统疾病者；肺炎恢复期住院者；病史资料不全者。其中，符合RMPP定义^[12]，即经大环内酯类抗菌药物正规治疗7 d及以上，临床征象加重、持续发热、肺部影像学加重者有73例，为难治组。在余下3 691例MPP患儿中采用倾向性评分法，按1 : 2的比例进行匹配，即每一个难治组个体与2个有最相似倾向性评分值的非RMPP个体匹配，最终纳入146例非RMPP患儿为普通组。

1.2 数据收集

收集所有患儿的基线资料：年龄、性别、入院时的病程、体重，利用倾向性评分匹配法对两组患儿进行匹配，得到组间协变量均衡的样本，再收集样本的临床表现、实验室及影像学检查结果。临床表现资料主要包括是否有高热、气促、细湿啰音、哮鸣音及皮肤损害。实验室检查主要包括白细胞（WBC）计数、血小板（PLT）计数、中性粒细胞百分比（N%）、血红蛋白（Hb）、C反应蛋白（CRP）、降钙素原（PCT）、白蛋白、谷氨酸-丙酮酸氨基转移酶（GPT）、天门冬氨酸氨基转移酶（AST）、LDH含量及鼻咽抽吸物MP-DNA载量等。影像学检查结果主要包括有无胸腔积液、肺实变、肺不张及肺坏死。

1.3 统计学分析

采用SPSS 24.0和R 3.60软件进行统计学分析。正态分布计量资料以均数±标准差（x±s）表示，组间比较采用独立样本t检验；偏态分布计量资料以中位数（四分位数间距）[M（IQR）]表示，组间比较采用Mann-Whitney U检验。计数资料以百分率（%）表示，组间比较采用卡方检验。先比较难治组与普通组不同指标的差异，再进一步进行logistic回归分析，若单因素logistic回归分析结果中自变量P < 0.1，则纳入该变量行多元logistic回归分析，以向前法（forward）筛选自变量，以P < 0.05为差异有统计学意义。通过logistic回归模型绘制LDH与RMPP发生的概率曲线。再以LDH为预测指标，敏感性为纵坐标，1-特异性为横坐标，绘制ROC曲线，同时通过约登指数法来寻找LDH最佳截断值。本研究通过限制性立方样条函数评估LDH和RMPP发病OR值间的关系。限制性立方样条图是应用样条函数与logistic回归相结合的限制性立方样条法分析自变量与应变量之间的非线性关系，以LDH最佳截断值为参考值，应用限制性立方样条法分析LDH与RMPP的剂量反应关系，其中LDH为连续变量，根据其百分位数的分布取4个节点（P₅、P₃₅、P₆₅、P₉₅），在R 3.60软件中绘制限制性立方样条图。最后对LDH进行决策曲线分析。决策曲线分析是一种评估临床预测模型、诊断实验和分子标记物的统计学方法，以净获益率为纵坐标，高风险阈值为横坐标，绘制DCA曲线^[15]，整体分析LDH这个指标在临床应用中在不同的概率下做决策，能否产生有效的临床价值。

2 结果 2.1 一般情况

普通组和难治组两组患儿年龄、性别、体重及入院时的病程的比较差异无统计学意义（P > 0.05），见表 1。

2.2 两组相关临床因素的比较

普通组和难治组两组高热发生率、WBC计数、PLT计数、N%及血清CRP、PCT、Hb、白蛋白、GPT、AST、LDH含量的比较差异有统计学意义（P < 0.05）；两组鼻咽抽吸物MP-DNA载量及胸腔积液、肺实变、肺不张、气促、皮肤损害发生率的比较差异有统计学意义（P < 0.05）。见表 1。

2.3 LDH值和RMPP发生的关系

先以LDH值为横坐标，RMPP发生概率为纵坐标，通过单因素logistic回归绘制LDH与RMPP发生的概率曲线如图 1，可见不同的LDH值均有相对应的RMPP发生概率，不同的RMPP发生概率也有相对应的LDH值。再纳入经单因素分析具有统计学意义的指标作为自变量，是否发生RMPP作为因变量，行多因素logistic回归分析，结果显示：高热、Hb值、LDH值、肺实变是RMPP发生的独立预测因素（均P < 0.05），见表 2。

2.4 LDH对RMPP发生的预测价值

以LDH为预测指标，敏感性为纵坐标，1-特异性为横坐标，绘制ROC曲线如图 2，计算ROC曲线下面积为0.836（95%CI：0.778~0.894）。当LDH为384.5 U/L时，约登指数最高，为0.599，故将LDH为384.5 U/L设置为最佳截断值，预测RMPP的敏感性和特异性分别为78.1%和80.8%。

2.5 LDH与RMPP关联强度的剂量反应分析

以LDH最佳截断值384.5 U/L为参考值，应用样条函数与logistic回归相结合的限制性立方样条法分析LDH与RMPP发生的剂量反应关系，如图 3，图中横坐标为LDH的连续变化，纵坐标为相对应的预测值（OR），阴影部分表示95%CI。由图可直观看出：LDH连续变化与RMPP的关联强度呈非线性剂量反应关系（P < 0.01）。结果显示，LDH与RMPP发生大致呈显著的正相关，特别当LDH > 384.5 U/L时，随着LDH升高，RMPP患病的风险显著上升。

2.6 LDH与RMPP的决策曲线分析

敏感性、特异性和ROC曲线下面积仅测量该模型的诊断准确性，未能考虑特定模型的临床效应，而决策曲线分析的优势在于它将患者或决策者的偏好整合到分析中^[15]。以净获益率为纵坐标，高风险阈值为横坐标，绘制DCA曲线（图 4）。其中高风险阈值设为（0，1）。由图 4可看出，在高风险阈值为0.1~0.7时（由图 1可获得LDH值相对应的发生风险值），净获益率大于0，有临床意义；且高风险阈值在0.1~0.7之间取值越小，净获益率越高，临床意义越大。

3 讨论

RMMP常合并许多严重肺内外并发症，病情进展迅速，危及生命^[12]。因此，尽早识别、及时治疗在临床诊疗中特别重要。有研究显示，MPP可导致多种炎性细胞浸润和细胞因子释放增加。因此，可将血清中一些炎性因子或细胞因子作为早期识别MPP的生物标志物^[4]。有多篇文献报道，血清LDH可用于预测儿童RMPP的发生^[16-17]。

本研究发现，LDH是儿童RMPP发生的独立预测因素，其ROC曲线下面积为0.836，表明LDH预测儿童RMPP的准确性较高，与国内外研究结果一致^[18-19]。本研究中LDH的最佳截断值为384.5 U/L，敏感性为78.1%，特异性81.8%。在Lu等^[16]的研究中，LDH的截断值为379 U/L，敏感性为48%，特异性85.8%。在Inamura等^[19]的研究中，LDH的截断值为410 U/L。这些研究中LDH截断值与本研究略有不同，可能与样本量、地区、样本人群不同等有关。既往的研究是在传统的logistic模型中计算LDH与RMPP发生的关联强度，忽略了LDH连续性变化与RMPP关联强度的变化轨迹。人为地将LDH进行分段研究不仅可能会损失信息量，还可能得到不准确的结果。所以不能只强调带有一定主观性的某个LDH分段的RMPP发生的风险，更需着眼于描述LDH实际发生细微变化时的剂量-反应曲线，以尽可能客观、清楚地展现二者的关联性。而本研究通过限制性立方样条模型将LDH定量数据与RMPP的发生结合，实现了非线性剂量-反应关系的连续性呈现，并直观地刻画了关联强度的变化轨迹。同时通过logistic回归模型构建LDH与RMPP发生概率的关系，显示不同的LDH值均有相对应的RMPP发生概率，不同的RMPP发生概率也有相对应的LDH值；而在决策曲线分析中，横坐标为不同的LDH值相对应的RMPP发生风险，纵坐标为净获益率，显示在发病风险阈值为0.1~0.7时，净获益率大于0，且发病风险阈值在0.1~0.7之间取值越小，净获益率越高，提示临床上当LDH在此范围（RMPP发生风险阈值为0.1~0.7相对应的LDH值范围）时，净获益率大于0，提示预测RMPP的发生有临床意义，且LDH值在此范围取值越小，净获益率越高，提示预测RMPP的发生越有临床意义。目前的临床研究大部分着重于利用ROC曲线对预测模型的准确性进行评价，对模型造成的临床后果很少涉及。但只关注准确性是不全面的，还需了解预测模型能否对临床实践有所帮助。而决策曲线分析就是这样一种简洁的决策分析模型，更倾向于分析决策是否有临床获益。

本研究的优势在于利用倾向性评分匹配法对患儿基线数据进行了匹配，避免了因基线数据差异而出现偏倚。通过多元logistic回归分析发现LDH是RMPP发生的独立预测因素，用剂量反应分析观察不同LDH的变化对RMPP发生风险的预测，并采用决策曲线分析来评估其预测价值。

本研究也存在一定的局限性，即缺少外对照，无法对建立的模型进行扩展性的验证。因此，后续尚需对当前建立的指标进行前瞻性试验，了解该模型的稳定性。

总之，本研究显示LDH是儿童RMPP发生的独立预测因素，对RMPP的预测有重要价值，与RMPP发生大致呈显著的非线性剂量反应关系，特别当LDH > 384.5 U/L时，随着LDH的升高，RMPP发生的风险显著上升。