武汉新型冠状病毒2019-nCoV基因组分析
武汉新型冠状病毒2019-nCoV基因组分析
写在前面
2020 年 1 月 10 日，武汉发现的新型冠状病毒的基因组被破译。该病毒测序由复旦大学生物医学研究院张永振教授领导的协作团队完成。悉尼大学的 Edward C. Holmes 代表团队上传结果，并最先在 Twitter 公布了这一消息。目前该病毒序列在 virological.org 网站中可以查询，也已经上传到了 GenBank 中（但未公开）。序列下载地址如下：
http://virological.org/t/initial-genome-release-of-novel-coronavirus/319
研究团队同时声明：“请随时下载、共享、使用和分析此数据。如果您希望在杂志上发表使用这些数据得到的研究结果，请与我们联系。如果您还有其他疑问，请直接与我们联系。”
据1月12日世界卫生组织官微报道，中国已与世卫组织分享新型冠状病毒（2019-nCoV）全基因组序列，可以看出该新型冠状病毒目前已被命名为2019-nCoV（2019新型冠状病毒）。
上一篇文章发表后，大家关注度很高，为避免简单粗暴的以讹传讹，本次进行后续分析报道。
基因组注释
利用NCBI上的ORFfinder工具进行该病毒的基因组注释分析。
软件给出了14个主要的开放阅读框（可以简单理解为具有编码基因的能力）。下图中标明了orf1ab和四个主要结构蛋白的位置。
结构基因S的比较
武汉毒株的S基因与bat-SL-CoVZC45的核苷酸相似性只有83.45%（这个数值应该取得是下图range1的结果，整个S基因的相似性应该低于这个数值）。
下图是核苷酸比对的情况，可以看到S基因中间有一段差异很大。
其他结构基因的比较
E基因的核苷酸序列相似性为98.68%，M基因的核苷酸序列相似性为93.42%，N基因的核苷酸序列相似性为91.03%。氨基酸序列相似性暂未分析，但结果应该跟核苷酸序列相似性相差不大。
系统进化树的构建
通过下载一系列相关的SARS-CoV、bat-SARS-CoV以及一些蝙蝠来源的α、β冠状病毒全基因组序列，通过MEGAX软件进行多重序列比对并构建了系统进化树。可以看到武汉病毒跟bat-SL-CoVZC45属于同一分支，但距离人源的SARS-CoV病毒较远。
武汉病毒与SARS-CoV的关系
画重点，武汉病毒不是SARS-CoV病毒，不是SARS-CoV病毒，不是SARS-CoV病毒，二者不能划等号。
猜你喜欢
10000+：菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature Cell专刊 肠道指挥大脑
文献阅读 热心肠 SemanticScholar Geenmedical
16S功能预测 PICRUSt FAPROTAX Bugbase Tax4Fun
生物科普: 肠道细菌 人体上的生命 生命大跃进 细胞暗战 人体奥秘
写在后面
为鼓励读者交流、快速解决科研困难，我们建立了“宏基因组”专业讨论群，目前己有国内外5000+ 一线科研人员加入。参与讨论，获得专业解答，欢迎分享此文至朋友圈，并扫码加主编好友带你入群，务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份，另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助，首先阅读《如何优雅的提问》学习解决问题思路，仍未解决群内讨论，问题不私聊，帮助同行。
学习16S扩增子、宏基因组科研思路和分析实战，关注“宏基因组”
点击阅读原文，跳转最新文章目录阅读