首次证实！华南海鲜市场不是病毒发源地！只是加速了传播，11月下旬或已开始人传人，12.8日和1月6日有2次大规模感染

全基因组测序数据首次证实：华南海鲜市场的新冠病毒来自其他地方，海鲜市场不是新冠发源地！

11月下旬可能就已经发生人传人，华南海鲜市场加速了病毒传播，可能在12月8日和1月6日发生过2次大规模新冠病毒传染。

新冠病毒爆发已两个多月，确定华南海鲜市场是不是唯一的发源地，对于寻找病毒的来源，以及确定中间宿主，对疫情的控制和避免再次爆发具有至关重要的意义。

中科院西双版纳热带植物园郁文彬教授联合华南农大和北京脑科中心的科研人员收集了来自全球四大洲12个国家的93个新型冠状病毒样本的基因组数据（截止2月12日），通过全基因组数据解析，分析后发现：

1.华南海鲜市场的新冠病毒是从其他地方传入进来的，并不是发源地；

2.新冠病毒可能在12月初，甚至11月下旬已经开始人传人，随后在华南海鲜市场加快了人际传播，在市场中发生快速传播，并蔓延到市场之外；

3.12月8日和1月6日发生过2次大规模新冠病毒传染；

4.现扩散的病例至少来自于3个途经，美国包括了5个来源。

5.新冠病毒基因组尚未发生重组，也就是说不会产生新的全新的病毒毒株。

研究于2 月 19 日发表在中国科学院科技论文预发布平台 ChinaXiv 上。以下报告解读来自西双版纳热带植物园（ID:XTBGly）。

研究数据来自全球各领域上传到GISAID EpiFluTM数据库中的共享数据，通过全基因组数据解析，追溯传染源及扩散路径。研究发现，收到的93个样本包含58种单倍型，可以归纳为五组（图1），包括3个古老超级传播者单倍型（H1, H3和H13）和2个新的超级传播者单倍型（H56和mv2）；华南海鲜市场的新型冠状病毒是从其他地方传入进来，在市场中发生快速传播蔓延到市场之外；同时，现扩散的病例至少来自于3个途径。新型冠状病毒在2月12日之前发生过2次明显的种群扩张（分别是12月8日和1月6日）。

图1. 新型冠状病毒58种单倍型的演化关系和地理分布格局（A，B），单倍型之间的可能演化关系（C），以及新型冠状病毒的可能传播和扩散路线(D)。A和B圆圈中的数据是样本数量。

华南海鲜市场的新冠病毒是从其他地方传入的

基于120个变异位点得到58种单倍型（基因类型），单倍型演化关系显示，单倍型H13和H38是比较“古老的”单倍型，通过一个中间载体（mv1，可能是一个祖先单倍型，可能是来自中间宿主或者“零号病人”）与蝙蝠冠状病毒RaTG13关联，并通过单倍型H3衍生出了单倍型H1。与华南海鲜市场有关联的患者样品单倍型都是H1及其衍生的单倍型H2，H8-H12（图1, A），而一份武汉样品单倍型H3与华南海鲜市场无关。可见，华南海鲜市场的新型冠状病毒是从其他地方传入进来，在市场中发生快速传播蔓延到市场之外。

另外，根据病患发病时间记录和种群扩张时间推断，也印证了之前《柳叶刀》的一篇研究华南海鲜市场不是病毒发源地的推论。

2020年1月24日，北京中日友好医院曹彬及北京协和医学院王建伟共同通讯在国际顶级医学期刊Lancet 《柳叶刀》在线发表题为“Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China”的研究论文，该研究在追踪第一批 41 个确诊病例发现，这批最早在 12 月 1 日报告的病例中，总共有 13 个病例与海鲜市场无接触史。也就是说，病毒在进入海鲜市场之前就已经潜入武汉居民的生活中。

点击图片查看详情

对“古老的”单倍型H13和H38的病毒样品溯源发现分别是来自深圳的病患（广东首例）和美国华盛顿州的病患（美国首例）。他们的旅行记录表明应该都是2019年12月底至2020年1月初在武汉探亲期间被感染的。现有武汉样本中没有检测到H13和H38单倍型，可能是因为现有样品主要采自几家定点医院，而且样品采集时间局限于2019年12月24日和2020年1月5日。如果能在武汉其他医院早期的病患检测到这两种单倍型，将对于寻找病毒来源非常有帮助。

图 | 新型冠状病毒单倍型的样品采集时间情况。红色圈的样品是确认与华南海鲜市场有关；蓝色圈的样本确认与海鲜市场无关。

追踪“零号病人”：可能是接触蝙蝠的人，理论上是能够追溯的

DeepTech深科技（ID:deeptechchina）采访了该论文作者郁文彬。据DeepTech报道：

华中农业大学一位匿名教授评价说，这是一项基于现有新型冠状病毒基因组序列的研究。理论上根据其变异的情况，是能够溯源各个样本之间的关系，到底是父子关系还是兄弟姐妹关系，到底复制了多少代。那么就可以最终追溯出最早的新冠病毒从哪里来。当然这里有一个假设是病毒是从动物传播到人。

在华南海鲜市场之外发现单倍型 H3 也就意味着，H3 携带者传染了去华南海鲜市场的人，结果将人流拥挤的华南海鲜市场变成了病毒传播爆发点。

他认为，溯源很重要，这个研究是通过大数据分析的，非常有意义。最核心的调查是最早感染新冠病毒的那几个人，这尤其需要医院配合，但是医院可能现在没心思、没精力来做这件事，那些医生可能也没有进化基因组学的背景，所以他们也可能没有认识到这件事的价值。

他说，至于这个最原始的携带者，也就是“零号病人”，可能是接触蝙蝠的人。理论上还是能够追溯的，即使零号病人已经去世，在某个医院应该还保留了他的血液样本。

郁文彬认为，最原始的单倍型不太可能存在现在的患者身上，因为病毒的变异太快。他们现在能拿到的数据已经是最早是 12 月 24 日，因为没有早期的样本就很难去溯源，“我们猜测国内医院可能没有保存早期患者的组织材料，当然我们非常希望有早期的样品存在”。

郁文彬表示，这个研究有个不足，那就是 93 个样本中，武汉样本主要是早期的，如果有更多武汉样本进行基因组测序的话，可能在溯源方面可以找到更多的证据，比如说找到 H13 和 H38 单倍型，可能可以帮助找到病毒来源。

新冠病毒2月12日之前发生过2次明显的种群扩张分别是12月8日和1月6日

根据新型冠状病毒基因组发生重组时间推算1月之前的种群扩张发生时间是12月8日，该结果暗示病毒可能在12月初，甚至11月下旬即已经开始有人际传播，随后在华南海鲜市场加快了人际传播（图2）。研究推算2月份之前的种群扩张时间在1月6日，这个可能与元旦假期有关联。需要指出，这一天国家疾控中心发布了2级应急响应。当时的预警起到了一些警示作用，公众活动和出行都有所减少。如果当时的警示能引起大众更广泛的重视，那么1月份中下旬向全国和全球蔓延的病例会有所降低。研究人员进一步确认我国其他9个省区和其他11个国家的感染病例基本都是从武汉直接或者间接输入而来。

现扩散的病例至少来自于3个途径美国包括了5个来源

为了能够细分来源，研究人员将58种单倍型分成了五组（图1），采用标准是3个中心（古老超级传播者）单倍型（H1, H3和H13）和2个新的超级传播者单倍型（H56和mv2）。以此鉴别出广东的病毒可能有三个来源，重庆和台湾的病毒有两个来源。其中，广东深圳一家人在早期就通过人传人进行了传播。有较多样本的澳大利亚、法国、日本和美国，他们的患者感染源至少有两个，尤其是美国包括了五个来源。非常值得关注的是H56这个超级传播者单倍型，它同时是澳大利亚、法国和美国，以及我国台湾患者的传染源。其他国家患者因为样品比较少，大多数的来源比较单一，他们除了是武汉旅游输入或在武汉感染外，有一些人可能是在广东、新加坡等地被感染。

图3. 120个变异位点在8个编码区的分布情况。统计类型包括替换或替代（左上角），密码子位置1-3（右上角），同义突变或非同义突变（左下角），和氨基酸性质（右下角）。

新冠病毒基因组尚未发生重组事件

研究人员发现新型冠状病毒基因组没有发生重组事件，93个基因组之间有120核苷酸发生了突变（0.41%序列长度），并均匀分散在10个编码区（χ2=1.958, df=9, P=0.99）。120个突变的核苷酸关联了119个氨基酸密码子，其中79个密码子 (65.83%)改变了氨基酸类型，并有42个(53.17%)氨基酸理化性质都被改变（图3）。这些氨基酸类型以及理化性质改变是否会影响新型冠状病毒的活性暂不清楚，需要其他蛋白组学和结构生物学方面的专业人士进行验证。本研究是版纳植物园综合保护中心生物多样性研究组的科研人员利用其在系统与演化领域的专长开展的，本研究提到单倍型演化关系分析方法可以结合到传染病学研究中，对于寻找传染源，以及精确的传播和扩散方向能提供非常重要的信息。

“新冠病毒基因组尚未发生重组事件”的结论是一个好消息，也就是说不会产生新的全新的病毒毒株。

相关研究结果已提交到中国科学院科技论文预发布平台

http://www.chinaxiv.org/abs/202002.00033，

可供疾控防疫相关部门参考，并接受专家学者的交流和讨论。

版权声明

本文由科研大匠综合自西双版纳热带植物园（ID:XTBGly）、DeepTech深科技，科研大匠，仅作分享。

更多阅读：

• 新冠肺炎，究竟是哪个国家的阴谋？