小丝氨酸重组酶是发现抗噬菌体防御系统的标志物

《PLOS Biology》：Small serine recombinases are markers for antiphage defense system discovery

　　

噬菌体疗法的重新兴起凸显了理解细菌通过抗噬菌体防御系统抵御噬菌体感染的必要性。传统上，鉴定抗噬菌体防御系统的策略缺乏通量，或对防御系统研究不足的细菌物种存在局限性。在此，研究人员开发了一个生物信息学流程，利用小丝氨酸重组酶来识别已知和未知的抗噬菌体防御系统。使用该方法查询参考基因组和宏基因组，研究人员表明小丝氨酸重组酶基因与抗噬菌体防御系统存在遗传关联，并可作为在多样细菌门中发现这些系统的诱饵。通过共转录预测和蛋白质结构域丰度的统计学分析，研究人员通过实验验证了该生物信息学方法，发现KAP P-loop NTPase与假定的抗噬菌体结构域融合，并证实原核Schlafen蛋白是一类新的抗噬菌体防御系统。该研究表明，小丝氨酸重组酶是跨多样细菌门发现抗噬菌体防御的可靠遗传标志物。

**论文解读：基于小丝氨酸重组酶的RADS流程及其在抗噬菌体防御系统发现中的应用**

**一、研究背景与科学问题**

多重耐药（MDR）细菌感染已成为严重的公共卫生威胁，可用的治疗手段极为有限，因此噬菌体（噬菌体）疗法重新受到关注。然而，细菌通过抗噬菌体防御系统适应噬菌体捕食，直接威胁噬菌体疗法的疗效。因此，全面理解抗噬菌体防御系统的多样性、活性和演化对于开发成功的噬菌体疗法至关重要。传统上，抗噬菌体防御系统的发现依赖于两类策略：其一，利用防御系统倾向于聚集在“防御岛”的特性，通过已知防御系统邻近未知基因来鉴定新系统；其二，通过构建基因组文库并在异源宿主中筛选噬菌体抗性。这些方法存在明显的局限性：缺乏通量，且研究样本多集中于大肠杆菌（*Escherichia coli*）和枯草芽孢杆菌（*Bacillus subtilis*）等模式细菌。新近发现的多数防御系统主要来自假单胞菌门（Pseudomonadota），在其他细菌门中代表性不足，提示传统方法限制了新防御系统的发现。此外，近期研究表明抗噬菌体防御系统常定位于移动遗传元件（MGEs）边界，且支持DNA移动的基因（如重组酶、整合酶）常嵌入或邻近防御岛。研究人员此前在粪肠球菌（*Enterococcus faecalis*）中鉴定了一个与小型丝氨酸重组酶（EF_B0058）共转录的IV型限制酶（TIV-RE, EF_B0059），但该重组酶对噬菌体抗性并非必需，提示其可能与防御系统的遗传维持有关。基于此，研究人员提出假设：其他抗噬菌体防御系统可能位于这些小型丝氨酸重组酶同源物的侧翼基因组区域，并可能与之共转录。为验证该假设，研究人员开发了名为“重组酶相关防御搜索”（Recombinase-Associated Defense Search, RADS）的生物信息学流程，利用小型丝氨酸重组酶作为遗传标记，系统发掘多样细菌中的抗噬菌体防御系统。

**二、主要技术方法**

研究人员构建了RADS生物信息学管道，其核心流程包括：（1）以粪肠球菌丝氨酸重组酶EF_B0058为诱饵，通过Diamond BLASTp在目标基因组或宏基因组重叠群中搜索同源物（氨基酸一致度≥30%）；（2）使用Seqkit提取重组酶基因上下游各5 kb的侧翼序列作为RADS重叠群；（3）利用Prodigal预测开放阅读框（ORFs），并通过InterProScan注释蛋白质结构域；（4）使用DefenseFinder识别已知抗噬菌体防御系统；（5）通过共转录预测（下游ORF同链且距离≤100 bp）和基于二项分布的统计检验（Benjamini-Hochberg校正）优先排序候选ORF。研究还采用了AlphaFold3结构建模、Kraken2宏基因组分类、系统发育分析等方法。候选防御基因被克隆至组成型表达质粒pLZ12A，并在粪肠球菌OG1RF（无前噬菌体的模式菌株）中进行异源表达与噬菌体攻毒实验。数据来源包括RefSeq（Bacillota和Bacteroidota完整基因组，2026年3月访问）、人类微生物组计划DACC粪便宏基因组（2,441个样本，2024年7月访问）和土壤宏基因组数据集（3,304个宏基因组）。

**三、主要研究结果**

**1. RADS流程的开发及其在Bacillota门中鉴定已知防御系统**

研究人员以EF_B0058为诱饵对Bacillota门13,553个完整基因组运行RADS，鉴定出9,478个同源重组酶，发现率为0.215 hits/Mb。AlphaFold3建模显示，与EF_B0058同源性最低的9个蛋白与大肠杆菌噬菌体Mu的小型丝氨酸重组酶结构高度相似，平均RMSD值为0.939 ?，验证了这些蛋白确为小型丝氨酸重组酶。系统发育分析显示这些同源物与EF_B0058聚类，且独立于大丝氨酸和酪氨酸重组酶。DefenseFinder注释显示，RADS重叠群中已知防御系统的发现率为每个重叠群0.681个，涵盖流产感染、CBASS、核酸限制、毒素/抗毒素、反转录子、tRNA降解和未知机制等类别，分布于29个细菌科和4个目。该部分证明了RADS能够在之前未被充分研究的细菌类群中识别已知防御系统。

**2. RADS在宏基因组中鉴定多样细菌门的防御系统**

对人肠道宏基因组运行RADS，小丝氨酸重组酶重叠群的发现率为0.016%。输入与输出的菌门分布高度一致（Bacillota 69.41%、Bacteroidota 21.34%、Pseudomonadota 3.26%等），表明小丝氨酸重组酶在微生物菌门中分布广泛。DefenseFinder在9个细菌门中鉴定出高丰度的已知防御系统，发现率达30.2%。在土壤宏基因组中，尽管重叠群质量较低（平均长度249.6 bp），RADS仍鉴定出34个防御系统，发现率高达85%，且输入与输出的菌门分布相似。对Bacteroidota门2,035个完整基因组的分析获得278个RADS重叠群，发现率为0.035 hits/Mb，已知防御系统发现率为每个重叠群0.877个，分布跨越13个科和5个目。值得注意的是，酪氨酸重组酶作为诱饵时表现出菌门偏好性，而小丝氨酸重组酶则具有独特的跨菌门无偏性，凸显了其作为通用遗传标记的优势。

**3. 二项分布与共转录分析发现新防御基因**

研究人员对Bacillota数据集中与小丝氨酸重组酶共转录的ORFs进行统计分析。结果显示，这些共转录ORFs的防御分数（Defense Score）中位数为0.0，意味着传统基于防御岛的“共关联”研究策略很可能遗漏这些系统。通过二项分布检验（Benjamini-Hochberg校正P＜0.05）并结合共转录预测，研究人员选出6个显著富集结构域的ORFs进行功能验证，包括DUF6056、假定噬菌体头部形态发生、超感染排除、Toprim DNA拓扑异构酶、Schlafen和KAP P-loop NTPase结构域蛋白。在粪肠球菌OG1RF中异源表达并用5株基因组多样的粪肠球菌噬菌体进行攻毒实验，发现其中3个ORFs表现出强效抗噬菌体活性（效率蚀斑（EOP）降低达6个对数数量级）：KAP P-loop NTPase、Schlafen和超感染排除结构域蛋白；另有2个（Toprim和DUF6056）表现出弱至可忽略的防御活性。

**4. KAP P-loop NTPase结构域与抗噬菌体结构域的融合**

对粪肠球菌KAP NTPase ORF（EF_KAP）的深入研究表明，该蛋白由N端跨膜螺旋、KAP NTPase结构域（含Walker B基序）和322个氨基酸的未知功能区域（RUF）组成。Walker B基序双点突变（D246A、D247A）完全消除抗噬菌体防御活性，表明其酶活性（核苷酸水解）为防御所必需；缺失RUF同样使防御失活，证明两个区域协同发挥作用。进一步，研究人员从人类肠道宏基因组RADS重叠群中选取3个来自不同宿主菌属（Fac*calibacterium*、*Lachnospira*和未分类重叠群）且C端区域各异的KAP P-loop NTPase蛋白，均能在粪肠球菌中赋予噬菌体抗性（EOP降低达3—6个对数数量级）。这表明KAP P-loop NTPase结构域融合未知功能区域构成了一类此前未描述的抗噬菌体防御蛋白：KAP结构域作为调控模块，激活或调控C端效应区域执行防御功能，而C端区域的多样性暗示其代表了多种抗噬菌体机制。

**四、讨论与结论**

该研究表明，小丝氨酸重组酶是跨细菌门发现抗噬菌体防御系统的可靠遗传标记。RADS流程的发现率在不同菌门间存在差异（Bacillota为0.681 hits/重叠群，Bacteroidota为0.877 hits/重叠群），这可能归因于DefenseFinder本身对假单胞菌门的偏向性及同源搜索在远缘物种中的局限性，凸显了在更多菌门中开展防御系统发现和功能解析的重要性。毒素/抗毒素系统是RADS重叠群中最常见的防御类别，这可能与其功能上的“混杂性”（可限制多种噬菌体）及移动适应性强有关；CRISPR-Cas系统在RADS重叠群中极为罕见，因其对靶噬菌体高度特异且需大型CRISPR阵列。研究人员推测，小丝氨酸重组酶可能通过促进水平基因转移、介导DNA倒位调控防御基因表达、或在噬菌体捕食压力下通过快速重组与基因转移增强群体遗传多样性等方式，与防御系统协同演化。对于KAP P-loop NTPase，尽管二十年前已有其参与抗噬菌体防御的假说，本研究首次提供了实验证据：该类蛋白代表一种新类别，其KAP结构域融合于多样的抗噬菌体效应区域，通过酶活性调控防御功能。研究还独立证实了原核Schlafen蛋白的抗噬菌体防御功能，其C端缺乏已知结构域却仍具有广谱防御活性，提示存在尚未阐明的噬菌体感知机制。综上，RADS为系统发现抗噬菌体防御系统提供了新的范式，未来可扩展至任何目的基因的侧翼基因组分析，有助于更均衡地评估防御系统在细菌域中的分布与演化。