一种基于单菌体物理形态特征的抗菌药物敏感性检测方法
一、引言
细菌感染是全球高死亡负担的重要病因，抗生素耐药性（AMR）进一步加剧了感染诊疗难度。抗菌药物敏感性测试（AST）是指导精准用药、遏制耐药传播的关键环节。然而，临床主流方法——琼脂扩散法、肉汤稀释法及其衍生技术——依赖细菌群体在抗生素作用下的生长抑制结果，通常需8–24小时甚至更长的体外培养周期，难以满足感染早期精准用药的时效需求。此外，群体水平的检测本质上是菌群平均响应，无法揭示单个细菌在药物作用早期的异质性变化，可能掩盖个体差异，影响敏感性与稳定性。
近年来，基于显微成像与人工智能的单细菌快速药敏检测成为研究热点。以公开号CN121249842A为代表的一类方案，采用“多抗生素浓度+同一细菌连续多帧序列图像+ROI序列输入端到端深度学习模型”的技术路线。该类方案虽能实现单细菌存活或易感状态的判别，但仍存在显著不足：其一，依赖连续时间序列图像采集与序列级建模，数据获取与计算复杂度高；其二，端到端深度学习模型直接输出分类结果，特征表达缺乏可解释性，难以从形态学层面追溯药物作用机制；其三，模型训练基于多浓度条件，侧重于单细菌状态判别，与临床AST中“特定菌种-特定抗生素-标准判读体系”的快速药敏需求存在差距，且无法明确不同菌药组合达到稳定判别能力的最早时间窗口。
因此，亟需发展一种不同于连续序列成像深度学习方案的快速药敏检测方法：在临床相关药物作用条件下，基于单菌体终点静态图像的可解释形态学定量特征，构建面向特定菌种-抗生素组合的快速药敏判别模型，在缩短检测时间的同时提升方法的可解释性、标准化潜力及临床转化适用性。
二、核心框架：“菌种-抗生素-处理时长”独立检测单元
本方法的第一性原理是将每次药敏检测组织为独立的分析单元（Independent Detection Unit, IDU），其定义为：
IDU = 目标革兰阴性菌 × 抗生素种类 × 处理时长
菌种选取范围涵盖临床常见革兰阴性菌：大肠杆菌、肺炎克雷伯菌、铜绿假单胞菌。抗生素种类选自氨曲南、头孢他啶、美罗培南、环丙沙星。处理时长梯度设为15分钟、30分钟、60分钟、90分钟。以大肠杆菌为例，4种抗生素×4个时间点共构成16个IDU。每个IDU内的所有单菌图像继承该菌株-抗生素对应的金标准药敏标签（敏感S/非敏感non-S），不跨浓度、不跨菌种混合训练。
这一设计的核心价值在于：将药敏判别从“判断细菌存活与否”转换为“量化特定药物作用下特定菌种在特定时间点的形态偏移是否达到可判别阈值”。它同时解决了三个关键问题：一是对接临床折点（采用CLSI敏感折点浓度，而非自设多浓度筛选IC₅₀）；二是给出时间窗口（通过横向比较15、30、60、90分钟的AUC，取首个达标时长作为该组合的最快检测时间）；三是保留可解释性（形态指标为人可定义的几何量，而非神经网络黑盒激活值）。
三、实验流程与方法
3.1 药物处理
针对单个IDU，从血平板上挑取目标菌临床分离的新鲜菌落，配制菌悬液并调整至预定处理浓度（0.01 OD₆₀₀）。将设定浓度的抗生素溶液稀释加入米勒-欣顿肉汤（MHB）培养基中混匀，再将调制好的菌悬液接种。抗生素浓度依据菌种分类设定：肠杆菌（大肠杆菌、肺炎克雷伯菌）采用头孢他啶4 μg/mL、氨曲南4 μg/mL、美罗培南1 μg/mL、环丙沙星0.25 μg/mL；非发酵菌（铜绿假单胞菌）采用头孢他啶8 μg/mL、氨曲南8 μg/mL、美罗培南2 μg/mL、环丙沙星0.5 μg/mL。将培养体系置于37℃、160 rpm恒温摇床，孵育至该IDU设定的处理时长（15/30/60/90分钟）即终止。
3.2 单菌体成像
孵育结束后，对菌体样本进行固定（4%多聚甲醛，室温5分钟）与细胞核染色（Hoechst 33342，避光5分钟）。取1.75 μL处理后的菌体样本滴片，覆盖24 mm×24 mm盖玻片，采用倒置共聚焦激光扫描显微镜（Zeiss LSM880，405 nm激发，100×油镜）采集明场与荧光双通道图像，分辨率为1024×1024像素。每个IDU记录不少于3个视野，图像保存为.tiff格式。此处需特别强调：本方法采用终点静态采像，即在预设处理时长到达后对固定染色样本进行单帧成像，而非对同一细菌进行连续多帧序列图像采集。
3.3 单菌体分割与多维特征提取
使用Cellpose实例分割算法从每个单菌体图像中识别并分割出单个菌体，生成掩膜（mask），将单个菌体从背景及相邻目标中分离。基于分割结果，提取单个菌体的六项物理形态学指标：
• 长轴长度（Major Length, μm）：菌体最大延展尺寸；
• 短轴长度（Minor Length, μm）：垂直于长轴方向的最大尺寸；
• 实心度（Solidity）：mask面积与凸包面积的比值，表征轮廓缺损程度；
• 圆度（Roundness）：4π·面积/周长²，趋近1表示接近圆形；
• 骨架长度（Skeleton Length, μm）：形态中心骨架化后的路径长度；
• 核/胞比（Nucleus/Cell Ratio）：Hoechst荧光区域面积与菌体总mask面积的比值。
前两项与骨架长度刻画菌体长度及延展特征，实心度与圆度表征轮廓规则性和形态完整性，核/胞比反映菌体内核酸凝聚态变化。所有指标经图像标尺统一校准，量纲为μm或无量纲。
四、单一指标性能评估与联合建模
4.1 单指标评估
在每个IDU内部，对六项指标分别进行两项分析：首先采用Mann–Whitney U检验比较敏感菌与非敏感菌之间的指标分布差异，验证各指标在不同药敏表型间是否存在显著差异；其次，针对单一指标构建受试者工作特征（ROC）曲线，计算曲线下面积（AUC），量化其对敏感/非敏感状态的区分能力。单指标AUC值越接近1，表明该指标单独即可具备良好的药敏判别能力。此步旨在明确哪些形态指标具有独立判别效能，为后续联合建模提供特征筛选依据。
4.2 联合建模与最早时间窗口确定
将不同IDU的六项形态学指标作为联合输入变量，导入多种传统机器学习模型，包括逻辑回归（L1/L2正则化）、随机森林、直方图梯度提升（HistGB）、极端梯度提升（XGBoost）、轻量梯度提升（LightGBM）。模型输入为x ∈ ℝ⁶，输出为药物敏感/非敏感二分类标签。数据集按8:2比例划分为训练集与测试集，采用分层五折交叉验证进行训练与性能评估，评价指标为AUC与分类准确率。
通过比较不同处理时长下模型的AUC与准确率，确定各菌药组合达到预设性能阈值（AUC≥0.80且准确率≥0.80，优选≥0.90）的最短处理时长，即为该组合的最早可判别时间窗口。这一分析目标与仅基于连续序列图像进行单次分类的方案存在本质区别：本方法不仅实现药敏分类，更建立了不同菌药组合的时间响应规律。
五、与序列端到端方案的本质区别
本方法与现有基于连续序列图像采集并直接采用端到端深度学习模型的方案在技术路线上存在实质性差异，具体体现在以下几个方面：
维度 序列端到端方案 本IDU形态方案
输入数据 同菌多帧ROI序列 终点单帧静态图像
药物浓度 多浓度扫描 单临床折点浓度
特征构建 CNN自动提取 六项显式几何/荧光量
标签来源 存活/耐药状态 菌种-抗生素金标准S/non-S
时间处理 序列时序建模 离散时间点比选
输出内容 单次分类结果 分类结果+最早判别窗口
可解释性 低 高（指标可回溯至药理效应）
本方法不依赖连续动态视频信息，而是在临床敏感折点浓度条件下，通过预设时间点终点静态显微成像、单菌分割及显式形态学特征提取，结合单指标统计学评价与多维特征联合建模，实现对细菌敏感/非敏感状态的快速判别，并进一步确定不同菌种与抗生素组合达到稳定判别能力的最早时间窗口。这一路线在数据采集复杂度、特征可解释性、临床判读逻辑匹配度等方面均优于现有序列端到端方案。
六、工程化映射：丹麦Biosense oCelloScope系统
上述实验流程若要进入临床实践，必须将“固定染核+共聚焦成像+Cellpose分割+Python特征提取+机器学习建模”这一链条封装为一体化、标准化的检测系统。丹麦Biosense公司推出的oCelloScope微生物快速药敏监测系统在架构上与本文提出的IDU框架高度契合。
oCelloScope系统的核心特征包括：
• 以单菌体显微图像为入口，不依赖连续视频追帧，而是采集药物作用后预设时间点的静态视野；
• 采用自动化明场/荧光显微扫描替代手动滴片共聚焦，通量适配多孔板；
• 分析路径遵循“单菌分割→形态量化→统计/机器学习判别”的显式流程，而非端到端视频CNN；
• 判读逻辑按“菌种-抗生素-作用时长”组织，输出敏感/非敏感结论及可解释的形态偏移报告；
• 目标是将传统8–24小时的培养法压缩至小时级甚至亚小时级，同时保留每项形态指标的统计证据链。
换言之，本文披露的IDU框架、六项特征集、单指标AUC前置评估、最早时间窗口判定规则，恰好构成了oCelloScope类系统的算法与实验设计规范；而oCelloScope则将这套规范从“论文流程”转化为“开机选菌选药选时间自动出报告”的仪器行为。两者之间的关系是：前者提供可审计的方法学骨架，后者提供可复现的工程载体。
七、实证数据摘要
7.1 大肠杆菌-美罗培南组合
使用美罗培南（1 μg/mL）处理，60分钟时长下长轴长度、实心度、骨架长度在敏感与非敏感菌株间差异最为显著（P<0.0001），单指标AUC在该时间点达到峰值。联合随机森林模型在60分钟时AUC可达0.9以上。
7.2 肺炎克雷伯菌-头孢他啶组合
使用头孢他啶（4 μg/mL）处理，30分钟时长下圆度与实心度指标在敏感与非敏感菌株间区分明显，单指标AUC在30分钟时开始显著上升，联合模型在30分钟即可达到预设性能阈值。
7.3 铜绿假单胞菌-氨曲南组合
使用氨曲南（8 μg/mL）处理，15–30分钟内圆度与核/胞比指标即呈现显著差异，单指标AUC在30分钟时接近平台，联合LightGBM模型在30分钟达标。这是三例中最早的可判别时间窗口。
上述结果表明，不同菌种与抗生素组合的形态学响应速度存在差异：铜绿假单胞菌在氨曲南作用后15–30分钟即可形成稳定判别特征，而部分肠杆菌在β-内酰胺类药物作用下需更长时间（60–90分钟）才能达到稳定判别状态。这一发现为建立不同菌药组合的快速药敏检测时间窗口提供了实证基础。
八、总结
本文提出了一种基于单菌体物理形态特征的抗菌药物敏感性检测方法，其核心创新在于：
• 以“菌种-抗生素-处理时长”建立独立检测单元，在临床敏感折点浓度下进行短时药物处理与终点静态显微成像；
• 提取长轴长度、短轴长度、实心度、圆度、骨架长度、核/胞比六项显式形态学指标，并开展单指标统计学评价与多指标联合机器学习建模；
• 通过比较不同处理时长下的模型性能，确定各菌药组合达到稳定判别能力的最早时间窗口。
该方法在技术路线、输入数据形式、特征构建方式、标签定义逻辑及最终判读目标方面，均与现有基于连续序列图像采集并直接采用端到端深度学习模型的方案存在本质区别。同时，该方法与丹麦Biosense oCelloScope系统的工程化架构高度兼容，具备良好的标准化潜力和临床转化前景。总体而言，本方法将药敏检测周期从8–24小时缩短至90分钟以内，在保证判别准确性的同时提升了方法的可解释性与临床适用性，为快速药敏检测技术提供了新的实现路径。
相关新闻推荐
1、低盐发酵甜瓣子体系中产生物胺的菌株与降解菌株能力研究【上】
2、实时监测蓝藻细胞内用于外源生物合成的苯丙氨酸相对含量的方法