《Applied and Environmental Microbiology》：Genomic and functional divergence analysis of closely related Teseptimavirus phages reveals dual HNH endonuclease positioning and its role in plaque formation
编辑推荐：
抗生素耐药性重新激发了人们对噬菌体作为抗菌剂的兴趣，然而噬菌体-宿主相互作用的复杂性往往导致不一致的感染结果，这凸显出解析密切相关噬菌体之间表型多样性分子基础的必要性。本研究比较了两株密切相关的T7样噬菌体——多价噬菌体C5和N30，二者均属于Teseptim
抗生素耐药性重新激发了人们对噬菌体作为抗菌剂的兴趣，然而噬菌体-宿主相互作用的复杂性往往导致不一致的感染结果，这凸显出解析密切相关噬菌体之间表型多样性分子基础的必要性。本研究比较了两株密切相关的T7样噬菌体——多价噬菌体C5和N30，二者均属于Teseptimavirus属，共享89.67%的基因组一致性且具有相同的宿主范围。尽管C5表现出更高的吸附率、更大的噬菌斑尺寸和更高的铺板效率（efficiency of plating, EOP），N30却显示出更强的抑菌活性和更大的爆裂量（burst size），这强调了全面评估噬菌体感染特性对治疗或生物防治应用的重要性。比较基因组学揭示了噬菌体C5和N30之间整体高度相似，但在一个假定蛋白和一个组氨酸-天冬酰胺-组氨酸（histidine-asparagine-histidine, HNH）内切酶上存在显著差异。对Teseptimavirus属的更广泛分析进一步揭示了HNH内切酶的两种不同基因组定位：一种位于基因组中段紧邻病毒粒子装配蛋白，另一种位于基因组末端紧邻终止酶大亚基（terminase large subunit），提示其在衣壳装配和DNA包装中的功能特化。三级结构预测与比较结果显示，尽管序列分歧中等，两种HNH内切酶仍表现出中度的结构相似性。此外，N30的HNH基因中一个单核苷酸替换显著提高了其在一种受试宿主菌株上的铺板效率，证明与终止酶相关的HNH内切酶在噬菌斑形成中发挥关键作用。综上所述，这些发现首次提供了Teseptimavirus属内HNH内切酶双基因组定位的证据，并将HNH内切酶确定为噬菌斑形成的关键遗传决定因素，为噬菌体感染变异性和宿主适应的分子机制提供了新见解。
**一、研究背景与科学问题**
细菌噬菌体（噬菌体）是生物圈中最丰富的生物实体，是驱动多种生态系统中细菌死亡和种群更替的关键因素。全球日益严峻的抗生素耐药性危机重新激发了人们对噬菌体疗法作为替代抗菌策略的关注。与广谱抗生素不同，噬菌体通过高度特异的噬菌体-宿主相互作用感染细菌，二者之间持续进行着协同进化竞赛。这种持续的共同进化常导致治疗挑战，包括宿主范围改变和感染效率差异。然而，即使亲缘关系很近的噬菌体，在相同细菌宿主上也可能表现出显著不同的感染结果，这反映了噬菌体-宿主相互作用的复杂性，并限制了噬菌体治疗的可预测性。因此，阐明这些表型差异背后的分子机制对于提高噬菌体应用的预测性和稳健性至关重要。
噬菌体感染是一个多步骤过程，包括宿主吸附、基因组注射、胞内复制、病毒粒子装配和细胞裂解。任何阶段的变异都可能导致感染表型的可测量差异。研究人员此前从猪肠道中分离了两株多价噬菌体C5和N30，它们感染大肠杆菌（*Escherichia coli*）和弗氏柠檬酸杆菌（*Citrobacter freundii*）菌株。尽管共享89.67%的基因组相似性并在五种受试菌株上表现出相同的宿主范围，C5始终展现出比N30更高的铺板效率（efficiency of plating, EOP）。这种亲缘关系密切的噬菌体之间的表型分歧为研究感染效率变异的遗传决定因素提供了契机。本研究的目的是系统比较C5和N30的关键感染表型并开展比较基因组分析，以鉴定可能与感染力差异相关的候选基因。
**二、研究人员开展的研究与主要发现**
研究人员系统比较了两株噬菌体在形态学、噬菌斑尺寸、吸附率、EOP值、抑菌活性和一步生长曲线等多个参数上的差异，发现C5在吸附率、噬菌斑尺寸和EOP值方面普遍占优，但N30在抑菌活性和爆裂量方面呈现更强的表现。比较基因组学分析表明，两株噬菌体基因组整体高度相似，但C5_ORF28和N30_ORF44无法配对，它们编码的蛋白均含有HNH内切酶结构域但序列一致性仅为46.25%。通过五轮盲传筛选，研究人员获得了N30的单核苷酸突变体N30-5，其在*E. coli* 94上的EOP显著提高，而这一突变仅发生在N30_ORF44（编码HNH内切酶）的G38510A位点，导致A135T氨基酸替换。进一步分析Teseptimavirus属内已分类噬菌体后发现，HNH内切酶基因存在两种不同的基因组定位：一种位于基因组中段紧邻病毒粒子装配蛋白基因，另一种位于基因组末端紧邻终止酶大亚基基因。三级结构预测发现两种HNH内切酶具有中等结构相似性（TM评分0.53，RMSD 2.21 ?），且N30_ORF44在S45位置存在一个额外残基。该研究的重要性在于首次发现了HNH内切酶在Teseptimavirus属中的双定位现象，并鉴定其为噬菌斑形成的关键遗传决定因素，为理解噬菌体感染变异性和宿主适应性的分子机制提供了新视角。
**研究方法主要包括**：利用透射电子显微镜进行噬菌体形态观察（样本来自本实验室分离的两株噬菌体C5和N30）；噬菌斑形态观察采用双层琼脂法加定量图像分析；EOP测定采用双层琼脂平板法，以*C. freundii* 113ju作为参考菌株（EOP=1）；吸附率测定在MOI为0.01条件下混育10 min后离心去上清，再通过双层琼脂法测定吸附后的噬菌体滴度；抑菌试验和一步生长曲线均在MOI 0.01条件下进行，连续监测7 h的OD
600变化和300 min内噬菌体滴度变化。基因组测序采用Illumina HiSeq平台，使用Trimmomatic v0.33和SPAdes v3.13.11进行质控和组装，GeneMarkS用于ORF预测，BLASTp用于蛋白比对，Clinker v0.0.32用于基因组可视化。蛋白序列比对使用ClustalW Omega，保守结构域鉴定通过Pfam数据库，三维结构预测用AlphaFold2，结构比对用TM-align，配体结合位点预测通过COACH服务器。
**三、研究结果**
**（一）噬菌体C5与N30感染特性的比较**：透射电镜显示两者均为典型短尾噬菌体科（Podoviridae）形态，具有二十面体衣壳和短而无收缩性的尾部。C5在五种细菌中的两种上产生明显更大的噬菌斑（如*C. freundii* 113ju，*E. coli* 94），在三种宿主上表现出显著更高的吸附率，并在所有五种细菌上具有更高的EOP值。相比之下，N30在*C. freundii* 113ju、*E. coli* 235和*E. coli* 94三种代表性宿主上均表现出明显更强的抑菌活性，且N30在*C. freundii* 113ju和*E. coli* 235中的爆裂量显著大于C5。
**（二）C5和N30的比较基因组学与蛋白序列比对**：C5基因组全长为39,551 bp，N30为39,266 bp，两者GC含量分别为48.82%和48.61%。C5和N30各有44个开放阅读框（open reading frame, ORF），其中C5_ORF28和N30_ORF44因序列一致性低（46.25%）而无法配对。C5_ORF28与T7_gp7.7的蛋白存在71.54%的一致性，而N30_ORF44仅有45%一致性。PfAM结构域分析显示两者均含HNH内切酶结构域，位于其一级序列N端56%区域内。另一对低一致性ORF为C5_ORF40与N30_ORF39（46%一致性），其BLASTp最佳匹配均来自Studiervirinae亚科噬菌体但全为未表征蛋白，其中超过半数的同源序列属于Teetrevirus属。
**（三）导致噬菌体EOP显著变化的突变鉴定**：经在*E. coli* 94上进行五轮盲传，获得N30-5突变株，其在*E. coli* 94上的EOP较N30显著提高（0.041 vs 0.0037），但仍显著低于C5（0.75）。N30-5与N30在噬菌斑尺寸、吸附率和生长曲线方面无显著差异。全基因组序列比较显示N30-5唯一突变位于N30-5_ORF44的G38510A位点，编码蛋白中产生A135T替换。N30_ORF44编码HNH内切酶，突变位于非保守HNH结构域之外，靠近蛋白C端。
**（四）HNH内切酶的序列与三维结构比较**：BLASTp分析显示两种HNH蛋白共享54%覆盖率和46%序列一致性，Pfam预测两者均含HNH内切酶结构域且属于His-Me_finger家族。AlphaFold2三维结构预测揭示两种蛋白N端区域存在大量重叠，超出保守结构域范围，表明其三级结构较一级序列更为保守。N30_ORF44比C5_ORF28多一个氨基酸残基（S45）。三维结构比对得出RMSD为2.21 ?，TM评分为0.53，中等结构相似性。COACH服务器预测C5_ORF28的潜在催化残基为C50、C55、H60，形成一个Cys-Cys-His锌离子配位基序。
**（五）Teseptimavirus属内HNH基因位置分析**：对在ICTV 2023版总物种列表中的Teseptimavirus属噬菌体进行比较分析发现，与N30_ORF44高度同源（≥90%覆盖率）的序列均位于基因组末端紧邻终止酶大亚基基因位点，而与C5_ORF28高度同源的序列则位于基因组中段中紧邻病毒粒子装配蛋白基因位点。其中仅有*E. coli* HZ2R8噬菌体同时含有可匹配至C5_ORF28和N30_ORF44的序列。
**四、讨论与结论总结**
讨论部分指出，C5和N30在感染表型上的多面性差异不能以单一指标（如噬菌斑大小或EOP）进行评价。噬菌斑尺寸和EOP强烈受扩散速率和潜伏期等因素影响，可能无法反映噬菌体在液体培养中的复制动态。C5虽具有较高吸附率，但过快或高频率的吸附可能导致裂解失控空间中子代噬菌体分散受限，限制有效感染。N30较强的抑菌活性与其较大的爆裂量和更高的子代产量相一致，提示胞内复制效率较初始吸附在决定杀灭能力方面可能发挥更关键作用。
HNH内切酶在Teseptimavirus属中呈现双基因结构、双基因组位置的现象系首次报道，一个位于基因组中段参与病毒粒子形态发生和装配，另一个位于基因组末端辅助DNA包装。HNH与EOP变化之间的关联虽初步获证，但突变仅在*E. coli* 94中效应显著，提示其作用具有宿主依赖性。研究人员明确指出，A135T突变通过较小EOP值改变（EOP提高约9倍）但未显著高于C5的水平，HNH内切酶的详细分子机制仍需进一步验证；同时不能排除其他因素如宿主细胞内限制修饰系统的影响。此外，C5_ORF40与N30_ORF39之间低序列相似性也指向两株噬菌体间潜在的感染动态相关基因变异。
研究结论：本研究首次在Teseptimavirus属中鉴定出HNH内切酶的双基因组定位现象，HNH内切酶或与病毒粒子装配模块相关或与DNA包装模块相关，提示其功能特化。同时证明单个HNH基因核苷酸替换即可显著改变噬菌体EOP，表明HNH内切酶是决定噬菌体感染力形成的关键因子，其基因组上下文可能塑造功能结局。本工作为密切相关的噬菌体之间微小遗传差异驱动表型分歧提供了新认识，并指出HNH内切酶可作为噬菌体理性改造的靶点，应用于治疗和生物防治目的。
（注：本文解读严格基于原文内容，未作任何推测性扩增。）
