《Frontiers in Immunology》：Genome-resolved gut metagenomics identifies an Escherichia coli-Collinsella signature associated with Wagner 4 gangrenous diabetic foot ulcers
编辑推荐：
糖尿病足溃疡（DFU）是2型糖尿病（T2DM）的主要并发症，但肠道微生物组是否捕获与晚期溃疡严重程度相关的系统性微生物特征仍不清楚。研究人员对43例伴活动性DFU的T2DM患者粪便样本进行鸟枪法宏基因组测序，比较Wagner 1–3级（n=30）与Wagner
糖尿病足溃疡（DFU）是2型糖尿病（T2DM）的主要并发症，但肠道微生物组是否捕获与晚期溃疡严重程度相关的系统性微生物特征仍不清楚。研究人员对43例伴活动性DFU的T2DM患者粪便样本进行鸟枪法宏基因组测序，比较Wagner 1–3级（n=30）与Wagner 4级坏疽性疾病（n=13）。从头组装与分箱回收440个去重复且达中等质量或高完整度/低污染阈值的宏基因组组装基因组（MAG）。群落水平多样性与优势分类组成未能区分Wagner 4与Wagner 1–3，表明晚期疾病未表现为广泛生态重构。特征水平分析识别出基因组解析的MAG谱。为优先筛选稳健候选，研究人员联合两种互补方法：随机森林（RF）稳定性选择识别出24个在重采样折叠中可重复高分类重要性的MAG，以及协变量校正MaAsLin2差异丰度检验。两法交集优先出三个各法均支持的MAG：一个在Wagner 4富集的Escherichia coli MAG与两个在Wagner 4耗竭的Collinsella MAG。此三MAG特征（袋外AUC=0.703）保留更宽24-MAG RF分类器大部分判别信息，分类器解释、差异丰度检验与单MAG丰度分布方向一致。功能注释进一步分离Wagner 4富集的Escherichia coli与Collinsella MAGs：Escherichia coli MAG携带抗生素耐药与毒力因子信号并具有呼吸代谢能力，而两Collinsella MAGs未检出耐药与毒力命中且代谢谱紧凑。探索性临床关联分析将该Escherichia coli–Collinsella丰度评分与更长DFU病程联系，符合慢性或晚期疾病负担的肠道微生物相关物。这些发现支持纵向肠道宏基因组验证，以确定该信号是否追踪DFU进展、治疗反应或恢复。
研究背景与立项依据：糖尿病足溃疡（DFU）是2型糖尿病（T2DM）后果严重的并发症，影响全球约6%糖尿病患者，伴随高复发、感染、住院与截肢风险。现有微生物学研究多聚焦溃疡局部，伤口宏基因组显示菌群随溃疡深度、血糖控制与临床结局变化；而“肠-皮肤轴”假说提出肠道微生物及其代谢物可经循环免疫与代谢信号影响远端皮肤与伤口生物学。既往16S rRNA研究仅达属级分辨率，无法解析单基因组功能。因此研究人员开展本研究，以基因组解析鸟枪法宏基因组学询问肠道MAG是否关联DFU严重程度（Wagner 1–3 vs Wagner 4），论文发表于《Frontiers in Immunology》。
主要技术方法：研究人员前瞻性纳入广西医科大学第一附属医院43例T2DM伴活动性DFU住院患者，按入院Wagner分级分Wagner 1–3级（n=30）与Wagner 4级（n=13），排除近4周全身抗生素暴露者；收集入院粪便标本于DNA/RNA Shield缓冲液?80℃保存。对样本行鸟枪法宏基因组测序（DNBSEQ-T7，双端150 bp，约10 Gb/样），经KneadData质控与GRCh38宿主去污染后，MEGAHIT从头组装、MetaBAT2分箱、CheckM质量评估、dRep去重复得440个MAG；CoverM定量，GTDB-Tk注译。统计上并行随机森林（RF）五折稳定性选择（Boruta式Gini重要性，≥3折入选）与协变量校正MaAsLin2差异丰度检验，交集定三MAG；SHAP解释方向，ROC/DeLong评AUC；RGI/CARD注译抗生素耐药基因（ARG）、ABRicate/VFDB注译毒力因子（VF）、DRAM注译代谢模块；Spearman与相关线性模型析临床关联。
研究结果：
队列特征显示Wagner 4炎症负担更重：两组年龄、BMI、HbA1c、糖尿病病程相似（P>0.4）；Wagner 4组C反应蛋白（CRP）中位数44.4 mg/L（vs 5.0，P=0.005）、红细胞沉降率（ESR）92 mm/h（vs 58，P=0.003）更高，血红蛋白更低（86.0 g/L vs 116.5，P=0.022），且男性比例更高（92.3% vs 46.7%，P=0.013），性别作为协变量纳入后续模型。
群落水平肠道微生物组未分离Wagner 4与Wagner 1–3：440 MAG跨11门以Bacillota（69.8%）、Bacteroidota（12.5%）、Actinomycetota（9.5%）、Pseudomonadota（3.9%）为主；Inverse Simpson与Shannon α多样性无组间差（P=0.47、0.20），Bray-Curtis β多样性PERMANOVA无全局结构（R2=0.022，P=0.492），前十丰度MAG仅见Escherichia coli在Wagner 4略扩张，说明严重度变异不在总多样性或优势分类结构。
稳定肠道MAG特征区分Wagner 4：RF五折稳定性选择得24个RF-selected MAG（含Escherichia coli、Klebsiella pneumoniae、Proteus mirabilis及Collinsella、Gemmiger等），SHAP显示高Escherichia coli丰度推预测向Wagner 4，高Collinsella数种推向Wagner 1–3。
Wagner 4在分类各级均现Escherichia富集与Collinsella耗竭：MaAsLin2校正年龄/性别/BMI/HbA1c后，属级Escherichia（+2.98，P=0.043）与Duodenibacillus（+3.02，P=0.028）富集；种/MAG级Escherichia coli（+2.93，P=0.045）富集，9个Collinsella MAG中8个系数为负、5个P<0.10，其中Collinsella sp. [S015_bin.13]（P=0.033）、[S028_bin.32]（P=0.035）、[S033_bin.6]（P=0.049）耗竭。
三MAG信号（Escherichia coli富集＋Collinsella sp905214615与Collinsella sp. [S028_bin.32]耗竭）捕获Wagner 4：RF与MaAsLin2交集得此三MAG（完整度94.3%–97.8%、污染≤0.81%）；三MAG RF模型袋外AUC=0.703（五折CV 0.733±0.177），与24-MAG模型（0.728）无统计差（DeLong P=0.78）。三MAG在两组患病率均近100%，但丰度重分布：Escherichia coli中位5.56% vs 0.46%（P=0.083），两Collinsella MAGs在Wagner 4显著更低（P=0.019、0.046），即信号为丰度再分配而非有无之别。
功能注释将Wagner 4富集Escherichia coli MAG联于毒力富集与呼吸代谢能力：该Escherichia coli MAG携3个RGI Perfect/Strict ARG命中覆盖头孢菌素、青霉素β-内酰胺、氟喹诺酮、四环素、甘氨环素、酚醇、利福霉素及消毒 antiseptics（另注碳青霉烯与单环β-内酰胺），含99个VFDB毒力因子（运动30、黏附23、效应物递送17）；DRAM示糖酵解完整度1.00、磷酸戊糖0.86、复合体II 1.00、TCA 0.88、复合体IV部分0.33，具兼性厌氧呼吸能力。两Collinsella MAGs无ARG/VF命中，DRAM仅含糖酵解（0.78）、磷酸戊糖（0.57–0.71）、F型ATP合酶（0.75），缺呼吸复合体II/IV、TCA不全（0.125）、无SCFA乙酸合成，呈代谢紧凑。
三MAG丰度谱关联DFU病程：复合Escherichia coli–Collinsella评分（Escherichia coli z分减两Collinsella z分均值）与DFU病程正相关（Spearman ρ=0.50，P=5.76×10-4，BH q=0.013），校正年龄/性别/BMI/HbA1c后β=0.46（P=0.003）；单MAG中Escherichia coli ρ=0.37（P=0.016），两Collinsella分别ρ=?0.39（P=0.010）、?0.31（P=0.044），提示该肠道信号为慢性/晚期疾病负担相关物。
讨论与结论翻译：在本DFU基因组解析肠道宏基因组研究中，最清晰微生物信号区分Wagner 4与Wagner 1–3，而非追踪全局群落多样性；该特征应首要解释为晚期、坏疽性或缺血性DFU的候选特征，非全部Wagner级间平滑梯度。信号经RF稳定性选择、SHAP方向性、MaAsLin2差异丰度、跨法交集、功能注释与探索性临床关联一致支持。群落α/β多样性不分离与MAG级对立丰度不冲突，因少数基因组解析谱系反向移动可被总体生态指标掩盖。Wagner 4富集Escherichia coli MAG具pathobiont样特征，伴CRP/ESR升高与贫血的背景炎症中，肠炎症促Enterobacteriaceae过度生长、硝酸盐呼吸优势机制使其生物学合理；Collinsella广泛耗竭（8/9 MAG负系数）反映该属内代谢紧凑糖解谱系收缩，而非单纯SCFA生产者丢失。复合评分与DFU病程正相关支持“晚期炎症疾病肠道印记”解读。研究局限为单中心横断面、n=43（Wagner 4仅13）、性别失衡（已校正与男性亚组敏感分析方向一致）、不能定因果、未控二甲双胍与住院暴露；三MAG签名为假设生成级优先候选，非验证诊断标志物。结论：基因组解析粪便宏基因组在T2DM伴DFU患者中识别出Escherichia coli富集与多Collinsella MAG协调耗竭的三MAG签名，关联Wagner 4坏疽性严重度与更长溃疡病程，功能上对应耐药/毒力/呼吸型Escherichia coli对代谢紧凑型Collinsella的反向格局，需纵向多中心队列验证其是否追踪进展、治疗反应与恢复。
