《European Journal of Clinical Microbiology & Infectious Diseases》：Genomic landscape of 340 virulent Acinetobacter bacteriophages reveals anti-CRISPR–enriched candidates for therapeutic prioritization
编辑推荐：
目的：碳青霉烯耐药鲍曼不动杆菌（CRAB）是全球健康重大威胁，现有抗生素愈发不足。本研究旨在建立一个全面的基因组框架，用于合理优先筛选烈性不动杆菌噬菌体作为治疗候选。方法：研究人员对340种烈性不动杆菌噬菌体进行了大规模比较基因组分析，整合了系统发育重建、泛基
目的：碳青霉烯耐药鲍曼不动杆菌（CRAB）是全球健康重大威胁，现有抗生素愈发不足。本研究旨在建立一个全面的基因组框架，用于合理优先筛选烈性不动杆菌噬菌体作为治疗候选。方法：研究人员对340种烈性不动杆菌噬菌体进行了大规模比较基因组分析，整合了系统发育重建、泛基因组分析、基于CRISPR间隔子（CRISPR spacer）的宿主相互作用映射、抗CRISPR（Anti-CRISPR）蛋白鉴定以及系统性抗菌药物耐药（AMR）基因筛查。结果：基因组大小跨度接近20倍，基因组大小与GC含量呈显著负相关（R2=0.139，ρ=-0.630）。系统发育分析揭示多个谱系高度分化，无优势进化枝。泛基因组分析识别出20,982个独特蛋白家族，其中76.2%为云基因，证实了高度开放的基因组架构。CRISPR间隔子匹配在100个噬菌体基因组中获得1,480个高置信度匹配，提供了广泛历史感染性的分子证据。抗CRISPR分析鉴定出不动杆菌噬菌体XC1为一个卓越的治疗候选，它携带55个预测的抗CRISPR蛋白，具有典型调控位点架构。AMR筛查检测到21种不同的AMR基因同源物（宽松RGI命中，氨基酸同一性22.5%-47.1%），在数据集中异质性分布，证实存在大量治疗上干净的候选物，同时标记出部分需在使用前进一步审慎评估的噬菌体。结论：这些发现为针对多重耐药不动杆菌感染的合理噬菌体候选物优先筛选提供了多标准基因组框架，对循证噬菌体疗法开发具有直接影响。
**背景与问题**
抗菌药物耐药（AMR）是21世纪最严峻的全球健康挑战之一，其中碳青霉烯耐药鲍曼不动杆菌（CRAB）被WHO列为优先病原体。由于现有抗生素愈发不足，噬菌体疗法作为替代策略重新受到关注。然而，针对不动杆菌噬菌体的基因组多样性、宿主相互作用历史及治疗安全性的大规模系统评估尚不充分，缺乏多标准整合框架来指导治疗候选物的合理优先筛选。
**研究内容与意义**
研究人员对340种烈性不动杆菌噬菌体进行了大规模比较基因组分析，整合了系统发育重建、泛基因组分析、CRISPR间隔子（spacer）宿主相互作用映射、抗CRISPR（Anti-CRISPR）蛋白鉴定及抗菌药物耐药（AMR）基因筛查。通过揭示基因组多样性、历史感染模式以及AMR基因分布，建立了多标准基因组框架用于合理优先筛选治疗候选噬菌体。论文发表于《European Journal of Clinical Microbiology》。
**主要关键技术方法**
1. **系统发育重建**：基于k-mer频率（k=13）的无比对全基因组距离计算，使用非加权配对算术平均法（UPGMA）构建系统发育树。
2. **泛基因组分析**：利用Prokka对340个噬菌体基因组进行蛋白编码序列注释，基于精确氨基酸序列同一性聚类为非冗余蛋白家族，并分类为云基因、壳基因和持久基因。
3. **CRISPR间隔子匹配**：从350个完整鲍曼不动杆菌基因组（NCBI GenBank来源）中提取CRISPR间隔子，使用BLASTn（最小95%同一性，25 bp比对长度）比对至噬菌体基因组数据库。
4. **抗CRISPR蛋白鉴定**：采用 guilt-by-association 策略，基于蛋白大小（典型Acr范围）、邻近Aca样调节基因及跨基因组重现性计算复合支持分数（最高3.0）。
5. **AMR基因筛查**：使用耐药基因识别（RGI）工具比对至综合抗生素耐药数据库（CARD 3.2.7），记录Loose hits的氨基酸同一性，并对dfrA阳性噬菌体进行整合子结构分析（IntegronFinder）。
**研究结果**
**数据集组成与基因组大小分布**
通过PhageAI预测感染生活方式，340个基因组被分类为烈性（86.51%），长度范围15.8 kb–286.4 kb（中位数45.2 kb），GC含量30.1%–51.6%（中位数39.4%），呈现近20倍跨度，表明显著谱系间基因组策略差异。
**基因组大小与GC含量显著负相关**
线性回归显示基因组大小与GC含量显著负相关（R2=0.139，p=1.16×10?12），斯皮尔曼秩相关确认强单调负关联（ρ=-0.630，p=4.97×10?3?）。较小噬菌体基因组GC含量中位数较高且变异更大，较大基因组（＞100 kb）GC含量较低且分布更窄。
**系统发育分析揭示广泛多样性**
全基因组UPGMA树显示高度异质拓扑结构，多个分离良好的谱系，无单一优势进化枝。紧密相关的噬菌体形成紧凑子进化枝，而大部分基因组占据孤立系统发育位置（长末端分支），表明独立进化轨迹和谱系特异性分歧。
**泛基因组分析揭示开放基因组架构**
共鉴定20,982个独特蛋白家族，76.2%为云基因（单基因组检测），无蛋白家族在所有340个基因组中100%保守。泛基因组累积曲线连续上升无平台，证实高度开放架构，水平基因转移和模块化重组是主要驱动力。
**频率分析识别保守的基板和裂解相关蛋白**
基板相关组分（楔形蛋白、轮毂组装蛋白、尾管连接蛋白、穿刺装置）在48–50个噬菌体基因组中最常见，反映其在感染中不可或缺的作用。内溶素（endolysin）是裂解中最常检测的蛋白（6个基因组），尾刺蛋白、肽聚糖转糖苷酶gp16等仅出现于少数基因组。
**CRISPR间隔子匹配识别具有广泛历史感染证据的噬菌体**
340个烈性噬菌体基因组中100个（29.4%）获得1,480个高置信度CRISPR间隔子匹配。多个噬菌体（如vB_AbaM_IME284、P11B、HN_Aba_01）各自积累≥28个匹配，表明在临床菌株中有反复广泛的历史感染事件。
**抗CRISPR富集噬菌体的治疗优先筛选**
Guilt-by-association分析发现大多数噬菌体编码1–4个预测Acr候选，而不动杆菌噬菌体XC1携带55个预测Acr候选，远超其他噬菌体。其候选蛋白（如DMIAEELE_00254）获得最高支持分数（3.0），聚集在Gp33样转录调节因子旁，符合典型Acr调控位点架构。
**AMR基因在噬菌体基因组中的分布**
系统筛选检测到21种AMR基因同源物（Loose RGI hits，氨基酸同一性22.5%-47.1%），涉及喹诺酮耐药（gyrA/gyrB/parC/parE）、叶酸途径（dfrA家族）和多重耐药或消毒剂耐药（emrE、KpnF、qac家族）。71个携带dfrA同源物的基因组经整合子分析未发现任何整合酶基因或完整整合子结构。多数噬菌体携带少量或无AMR基因，确认存在大量治疗上干净的候选物。
**讨论与结论**
讨论部分强调：基因组大小与GC含量负相关中仅有约14%变异被解释，表明其他进化压力（如宿主密码子偏好、水平基因转移）独立影响核苷酸组成；系统发育异质性符合噬菌体镶嵌模型，支持多谱系组合用于鸡尾酒设计；开放泛基因组中大量云蛋白可能编码新抗防御机制；CRISPR间隔子匹配为历史感染提供分子证据，但约70%无匹配的噬菌体仍可能感染CRISPR缺陷菌株或编码高效Acr机制；Acr富集的XC1噬菌体可作为克服CRISPR免疫的载体和工程平台；AMR筛查中低同源性同源物不具功能耐药性，但基因组安全筛查必不可少。
**研究结论**
340种烈性不动杆菌噬菌体的基因组分析揭示了基因组组成、系统发育架构和基因含量组织的广泛多样性。噬菌体群体呈现以附属基因含量为主的开放泛基因组，在序列水平上无可检测的保守核心基因组，反映了水平基因转移和模块化重组是塑造该噬菌体集合的主要进化力量。CRISPR间隔子匹配与抗CRISPR分析共同鉴定出具有广泛历史感染性和强抗防御能力的噬菌体，其中不动杆菌噬菌体XC1因其异常庞大的预测抗CRISPR蛋白库而成为未来实验评估的高优先级候选。系统性AMR筛查确认了基因组干净的候选治疗噬菌体的可用性，同时标记出部分携带AMR基因同源物的噬菌体，这些噬菌体在排除临床使用前需进一步验证。总体上，这些发现增进了对不动杆菌噬菌体多样性和噬菌体-宿主相互作用动力学的理解，并为针对碳青霉烯耐药鲍曼不动杆菌的噬菌体候选物合理优先筛选建立了基于原则的多标准基因组框架。对已鉴定高优先级候选物的进一步实验验证，有望将这些基因组见解转化为有效针对多重耐药革兰阴性感染的噬菌体干预手段。
