温带和丝状噬菌体作为石珊瑚组织损失病中细菌毒力的储存库

《Coral Reefs》：Temperate and filamentous bacteriophages as reservoirs of bacterial virulence in stony coral tissue loss disease

　　

溶原转换（Lysogenic conversion），即噬菌体（bacteriophages，感染细菌的病毒）整合到细菌基因组中并赋予宿主新的表型特征，是病原体出现的一个公认机制。虽然这一过程是人类和其他动物许多细菌性疾病的基础，但它在珊瑚疾病中的作用仍未得到充分探索。在此，研究人员调查了溶原转换在石珊瑚组织损失病（stony coral tissue loss disease，SCTLD）中的潜在作用，该病自2014年以来已导致加勒比珊瑚前所未有的死亡率，且其病原体仍未知。通过对27份原始珊瑚样本和来自佛罗里达礁脉（Florida Reef Tract）的174个公开宏基因组（metagenomes）（涵盖七个珊瑚物种），研究人员在患病珊瑚或患病珊瑚的视觉健康部分（分别为diseased tissue directly adjacent to lesion，DD和visually healthy tissue on diseased coral，HD）中检测到比视觉健康珊瑚（visually healthy coral，VH）更多携带与溶原性（lysogeny）相关基因组特征的独特病毒基因组。SCTLD相关的噬菌体主要预测感染先前与该病相关的细菌类群，如Vibrionales、Rhodobacterales和Flavobacteriia，并携带丰富且广泛分布的毒力因子（virulence factors），这些因子具有增强细菌定植、竞争或直接宿主损伤的潜力，包括Tse2、Zot、RTX、pneumolysin和cytolytic delta-endotoxin的同源物（homologs）。尽管因果关系尚未解决，研究人员的发现表明噬菌体具有在受SCTLD影响的珊瑚中横向转移（laterally transfer）细菌毒力基因的基因组能力。通过溶原转换获得基因可能有助于细菌毒力，同时维持群落分类谱，有助于解释先前的群落分析观察结果，并为疾病发病机制（pathogenesis）提供机制框架。

石珊瑚组织损失病（stony coral tissue loss disease，SCTLD）自2014年以来对加勒比和佛罗里达礁脉造成毁灭性影响，累及30多种造礁珊瑚，最易感物种死亡率达40%–100%。尽管传播模式（如沿佛罗里达礁脉的聚集热点）提示存在传染性病原体，但病因仍未知。先前研究显示，抗生素可减缓病变进展，但基于16S rRNA测序的微生物组研究未发现单一、一致的细菌病原体，且培养分离的细菌在实验中对组织损失的重现性不一致。同时，病毒（尤其是噬菌体）的作用受到关注：溶原转换（lysogenic conversion）——温带噬菌体（temperate phages）或丝状噬菌体（filamentous bacteriophages）整合入细菌基因组并改变宿主表型——是人类和动物病原体出现的已知机制（如霍乱弧菌中的霍乱毒素基因由丝状前噬菌体编码），但在珊瑚疾病中尚未被探索。本研究旨在调查溶原转换在SCTLD发病机制中的潜在作用，即噬菌体是否通过横向转移毒力基因促进细菌致病性。

研究人员收集了来自美国佛罗里达群岛（Florida Keys）的三个地点、三种珊瑚物种的27份原始样本（包括病变邻近组织DD和同一患病珊瑚的视觉健康组织HD），并整合了174个来自佛罗里达礁脉的公开宏基因组（涵盖七个珊瑚物种），共计201个宏基因组（DD 117份、HD 34份、VH 50份）。主要关键技术方法包括：（1）病毒检测与溶原性识别：使用geNomad、VIBRANT和CheckV鉴定病毒基因组，并通过五条互补策略筛选具有溶原潜力的候选病毒（前噬菌体proviruses、推定的温带噬菌体、丝状噬菌体Inoviridae），最终获得6551个目标噬菌体，去冗余后1384个代表性病毒基因组。（2）基因组分析：用BLASTn进行宿主预测（结合iPHoP和前噬菌体侧翼序列比对），用MetaCerberus和自定义VFDB隐马尔可夫模型进行毒力因子注释，通过vConTACT2构建基因共享网络。（3）统计比较：基于相对片段丰度（fractional abundance）计算Bray-Curtis距离并进行PERMANOVA检验，Kruskal-Wallis检验比较多样性指标和毒力基因丰度差异，指示物种分析（Indicator Species Analysis）识别疾病指示病毒。

**结果**

**Viral signatures of lysogeny in SCTLD metagenomes（SCTLD宏基因组中的溶原性病毒特征）**
通过多方法鉴定，共发现117个独特前噬菌体簇代表，其中112个来自患病珊瑚（DD或HD）。在1757个溶原转换候选病毒（包括前噬菌体、温带病毒和Inoviridae）中，137个仅存在于DD样本，99个同时存在于DD和HD样本，仅7个仅存在于VH样本。指示物种分析识别出9个与患病样本相关的指示病毒，其中2个为DD唯一指示病毒（如SINT_DD_R2S43C_NODE_111，特异性96%，保真度64%），7个为DD和HD共有指示病毒。

**Community profiles of lysogenic conversion candidates（溶原转换候选病毒的群落特征）**
非度量多维尺度分析（NMDS）显示，三组（DD、HD、VH）间病毒群落组成有显著差异（PERMANOVA，R²=0.041，p=0.033），HD组与DD和VH均显著不同。Alpha多样性指标中，HD样本的Shannon指数（4.89±0.12）和Simpson指数（0.97±0.005）显著高于DD和VH，均匀度也更高，但丰富度较低。科水平分析发现，Autographiviridae和Inoviridae在DD样本中丰度和流行率更高，Zobellviridae在HD样本中显著富集，而部分Caudoviricetes科仅在DD样本中检出（流行率低）。

**Virulence genes in lysogenic conversion candidate viruses（溶原转换候选病毒中的毒力基因）**
DD样本中检测到40种毒力基因（VFs），高于HD（38种）和VH（34种）。部分毒力因子同源物（如Zot、RTX蛋白FrpC、胆固醇依赖性溶细胞素Pneumolysin、羧酸碳储存调节因子CsrA）仅存在于患病珊瑚样本（DD/HD），尽管流行率较低（0.85%–6.84%）。丰度差异分析显示，LPS相关基因、极鞭毛基因和I型菌毛基因在DD/HD中显著高于VH，而VI型分泌系统（T6SS）基因在DD/HD中显著降低。

**Genome similarity network of lysogenic conversion candidates（溶原转换候选病毒的基因组相似性网络）**
经过对231个SCTLD相关溶原转换候选病毒（包括毒力基因编码病毒和疾病指示病毒）的基因共享网络分析，发现至少9个噬菌体科，主要宿主包括Alphaproteobacteria（40个病毒连接，其中58%编码毒力因子）、Clostridia（32个连接，50%编码毒力因子）和Bacteroidia等。三个Inoviridae病毒均仅存在于患病样本并携带多个毒力基因。网络中的病毒宿主预测包括Vibrionales、Rhodobacterales、Flavobacteriia等先前与SCTLD相关的类群。

**Virulence gene-encoding bacteriophages associated with SCTLD（与SCTLD相关的毒力基因编码噬菌体）**
在65个SCTLD致病性候选病毒中，13个编码可能直接损伤真核细胞的毒力因子，包括Tse2（NAD依赖性细胞毒素）、Exotoxin type G、细胞溶解δ-内毒素（cyt2Bb1）、热稳定肠毒素A（ystA）、Zonula occludens toxin（Zot）、辅助霍乱肠毒素（Ace）、铁调控RTX毒素（frpA）和硫醇激活的溶细胞素（Pneumolysin/Seeligeriolysin同源物）。例如，SCTLD最显著指示病毒（SINT_DD_R2S43C_NODE_111，感染Vibrio的Caudoviricetes）编码Tse2同源物；一个感染Flavobacteriia的前噬菌体（MCAV_DD_McD6_NODE_117）编码外毒素G和细胞溶解δ-内毒素；四个Inoviridae丝状噬菌体编码Zot和Ace，可能破坏珊瑚上皮屏障。

**讨论与结论**
讨论部分指出，温带病毒在患病珊瑚中富集，而非仅普遍存在（85%的温带病毒在健康和患病珊瑚中均有，但近15%仅见于患病样本），表明疾病状态改变了溶原病毒群落。HD样本的病毒群落Alpha多样性升高和组成差异可能反映病变形成前的早期微生物扰动或病变邻近效应。Autographiviridae（通常为裂解型，但部分携带溶原模块）和Inoviridae（慢性感染型）在DD样本中富集，提示它们可能与致病细菌互作。基因共享网络显示，Alphaproteobacteria（特别是Rhodobacterales）和Clostridia是SCTLD相关溶原转换候选病毒的主要宿主，且高比例携带毒力基因，暗示噬菌体介导的基因转移可能促成这些类群的致病性转变。毒力因子分析发现，DD和HD样本中多种潜在真核细胞损伤毒素（如Tse2、Zot、RTX、溶细胞素）的噬菌体编码，与已知病原体（如霍乱弧菌）中的溶原转换机制相似，说明温带和丝状噬菌体可作为细菌毒力的移动储库。

结论部分翻译如下：总体而言，研究结果突出了溶原转换在SCTLD发病机制中的潜在作用。在患病珊瑚中富集的温带病毒携带多种推定的毒力因子，包括已知破坏真核细胞膜、干扰组织屏障和调节宿主-微生物互作基因的同源物。这一模式提出了一个概念模型，即温带噬菌体的溶原转换导致疾病结局的异质性，尽管因果关联无法从这些数据中建立。根据特定前噬菌体的有无、诱导状态及水平基因转移事件，细菌共生体可能获得或失去细胞毒性能力。在此框架下，噬菌体可能导致珊瑚种内和种间的差异易感性，以及毒力和死亡率连续谱（类似于产毒与非产毒霍乱弧菌）。本研究中识别的特定病毒和毒力因子为实验验证提供了具体靶点，例如基于CRISPR的前噬菌体编辑、诱导试验或转录组调查。未来工作还可探索这些SCTLD相关细菌宿主的丝状病毒，以识别选择性靶向宿主细菌的基因（如PftP4样多肽），为环境机会病原体的管理提供新工具。