MPB：华大孙海汐等-从细菌基因组中预测活性前噬菌体工具Prophage Hunter的使用流程和常见问题

为进一步提高《微生物组实验手册》稿件质量，本项目新增大众评审环节。文章在通过同行评审后，采用公众号推送方式分享全文，任何人均可在线提交修改意见。公众号格式显示略有问题，建议电脑端点击文末阅读原文下载PDF审稿。在线文档(https://kdocs.cn/l/cL8RRqHIL)大众评审页面登记姓名、单位和行号索引的修改建议。修改意见的征集截止时间为推文发布后的72小时，文章将会结合有建设性的修改意见进一步修改后获得DOI在线发表，同时根据贡献程度列为审稿人或致谢。感谢广大同行提出宝贵意见。

从细菌基因组中预测活性前噬菌体工具Prophage Hunter的使用流程和常见问题

The Usage Process and Common Problems of Prophage Hunter, a Tool to Capture Active Phage from the Bacterial Genome

孙海汐^{1, 2, #, *}，李敏^{1, 2, 3, #}，宋文琛^{1, 2, #}，肖敏凤^{1, 2, 3, *}

1.深圳华大生命科学研究院,深圳；2.深圳市未知病原体应急检测重点实验室，深圳；3.中国科学院大学华大教育中心，深圳

*通讯作者邮箱: [e-mail]; [e-mail]

#共同第一作者/同等贡献

摘要：识别具有活性的前噬菌体，对于研究噬菌体和细菌的协同进化、噬菌体的生理生化以及工程化设计噬菌体等多种用途至关重要。这里，我们介绍了Prophage Hunter，该工具旨在从细菌的全基因组序列中寻找具有活性的前噬菌体。结合序列相似性的匹配及基于遗传特征的机器学习分类模型，我们开发了一种新颖的评分系统，该系统在验证数据集上预测活性前噬菌体中表现出比当前工具更高的准确性。该工具也可以选择跳过序列相似性匹配，这样有更多可能发现新颖的前噬菌体。Prophage Hunter提供一站式网站服务，包括从细菌基因组中提取前噬菌体基因组、评估前噬菌体的活性、鉴定系统发育相关的噬菌体、注释噬菌体蛋白的功能及可视化前噬菌体基因组位置信息等。Prophage Hunter可在https://pro-hunter.genomics.cn/免费使用。

关键词：前噬菌体，细菌，机器学习，注释，一站式分析

仪器设备

1.个人电脑：安装主流浏览器 (Chrome/Safari) 即可

实验步骤

1.准备输入数据：细菌基因组序列,可包含一条或多条序列 (FASTA格式，图1)。

图1. FASTA格式序列

2.在浏览器中输入https://pro-hunter.bgi.com/，进入网站主页 (图2)。

图2. Prophage Hunter网站主页

3.在网页的导航栏上单击"Start Hunting (开始狩猎)"按钮启动Prophage hunter程序 (图3)。

图3. 单击红色框中的Start Hunting (开始狩猎) 按钮以启动程序

4.在主页上，单击"Browse (浏览)"将一个或多个核苷酸序列以FASTA文件格式上传到网站 (图4a)。默认情况下，Prophage Hunter使用相似性搜索策略来标识初始前噬菌体区域。也可以通过勾选"Skip similarity matching (跳过相似性匹配)"框 (图4b)，用户可以跳过此过程以识别新型噬菌体。本示例以地衣芽孢杆菌Bacillus licheniformis DSM 13 (Accession Num.:NC_006322.1) 为输入。

图4. 单击"Browse (浏览)"上传FASTA文件

5.可选择输入电子邮箱地址以接收指向分析报告的超链接 (图5)。请注意，此超链接将在一周后过期。若不输入，则需在提交序列后跳转的分析页面等待分析完成，或自行将提交序列后跳转的分析页面链接复制保存，以便查看结果（详见步骤7、8）。

图5. 输入电子邮件地址以接收分析报告

6.勾选以加入"Join User Experience Improvement Program (用户体验改善计划)" (图6a)。有关更多详细信息，请单击"Join User Experience Improvement Program (用户体验改善计划)" (图6b)。

图6a. 勾选加入用户体验改善计划

图6b. 用户体验改善计划具体信息

7.点击"START HUNTING (开始狩猎)"按钮开始分析。

图7. 开始分析

8.分析完成需要等待5-15分钟，分析完成后点击蓝色区域链接即可到达结果页面 (图8)。

图8. 等待分析完成

结果

1.在结果页面的顶部是一个基因组浏览器，显示了每个预测的前噬菌体区域 (图9a)。活性前噬菌体区域（Category 为Active）和模糊区域（Category 为Ambiguous，即难以判断活性的前噬菌体区域）分别以天蓝色和灰色着色 (图9b)。若提交文件中含有多个序列，用户可以通过单击左上方的下拉菜单切换到基因组的其他染色体或Scaffold序列 (图9c)。

图9. 基因组浏览器显示预测的前噬菌体区域的基因座和活性类别

2.在结果页面底部的表格显示了每个预测的前噬菌体基因组区域的详细信息 (图10a)：

- Candidate ID (候选ID)：预测的前噬菌体基因组区域的ID；

- Sequence ID (序列ID)：输入文件中显示的细菌基因组序列ID；

- Start (开始)：预测的前噬菌体基因组区域的开始位置；

- End (结束)：预测的前噬菌体基因组区域的结束位置；

- Length (长度)：预测的前噬菌体基因组区域的长度；

- Category (类别)：预测的前噬菌体基因组活性类别，Active表示预测为具有活性(Score > 0.8)，Ambiguous表示该预测区域活性难以判断 (0.5 < Score < 0.8)；

- Score (分数)：预测的前噬菌体基因组活性分数，分数越高，支持预测区域的活性证据越强；

- Closest phage (最近源的噬菌体)：与现有噬菌体库相比，预测的前噬菌体区域最近源的噬菌体，" N/A"表示未在噬菌体库中找到近源噬菌体；

- Gene number (基因数目)：预测的前噬菌体区域内基因的数目。

图10. 每个预测的前噬菌体区域详细信息。

用户可以单击相应的按钮以下载对应预测的前噬菌体区域的基因组DNA序列、CDS序列或蛋白质序列 (图10b) 或所有分析结果 (图10c)。用户还可以单击最接近的噬菌体的分类名称以查看详细信息 (图10d)。

3.在结果页面顶部的基因组浏览器或下方的表格中单击预测的前噬菌体区域，可以查看预测的前噬菌体区域每个基因的注释情况 (图11a) 和与该区域同源的前5个最接近的噬菌体 (图11b)。本示例为点击表格中 Candidate 22 。

    图11a. 预测的前噬菌体中基因在数据库中 (NCBI NR、Pfam和InterPro) 的注释情况

图11b. 与预测区域同源的前5个最接近的噬菌体

致谢

本项目由国家重点研发计划项目 (2020YFA0908700)、深圳市孔雀团队项目(KQTD2015033117210153)支持。

参考文献

1.Song, W., Sun, H. X., Zhang, C., Cheng, L., Peng, Y., Deng, Z., Wang, D., Wang, Y., Hu, M., Liu, W., Yang, H., Shen, Y., Li, J., You, L. and Xiao, M. (2019). Prophage Hunter: an integrative hunting tool for active prophages. Nucleic Acids Res  47(W1):W74-W80.

猜你喜欢

系列教程：微生物组入门 Biostar 微生物组  宏基因组

专业技能：学术图表 高分文章 生信宝典 不可或缺的人

一文读懂：宏基因组 寄生虫益处 进化树

必备技能：提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

扩增子分析：图表解读 分析流程 统计绘图

在线工具：16S预测培养基 生信绘图

科研经验：云笔记  云协作 公众号

编程模板: Shell  R Perl

生物科普:  肠道细菌 人体上的生命 生命大跃进  细胞暗战 人体奥秘  

写在后面

为鼓励读者交流、快速解决科研困难，我们建立了“宏基因组”专业讨论群，目前己有国内外5000+ 一线科研人员加入。参与讨论，获得专业解答，欢迎分享此文至朋友圈，并扫码加主编好友带你入群，务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份，另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助，首先阅读《如何优雅的提问》学习解决问题思路，仍未解决群内讨论，问题不私聊，帮助同行。

学习16S扩增子、宏基因组科研思路和分析实战，关注“宏基因组”

点击阅读原文下载PDF审稿，或浏览器直接访问下载链接：http://210.75.224.110/github/MicrobiomeProtocol/04Review/210226/2003952HaixiSun970936/Protocol2003952.pdf