| 送交者: 苦难与荣耀[☆★★声望品衔12★★☆] 于 2021-07-02 21:24 已读 1175 次 2 赞 | 苦难与荣耀的个人频道 |
长标题：美国病毒学家Bloom新冠溯源研究恢复、依据的基础数据，是武汉大学“纳米孔靶向测序检测”技术研发、测试过程中测序、输出的新冠病毒基因序列片断，这些基因序列片断准确度欠佳，存在大量的核苷酸测序缺失。 6park.com对美国病毒学家Jesse Bloom的新冠溯源论文作了初步研究，发现了一些疑点、缺陷和错误。本文整理的是与论文的基础研究数据相关的一组问题。成文仓促，请各位网友批评、指正。 6park.com2021年6月22日，美国西雅图贺勤森癌症研究中心(Fred Hutchinson Cancer Research Center，福瑞德·哈金森癌症研究中心)的病毒学家杰西·布鲁姆（Jesse Bloom）在生物学预印本平台biorxiv发布了如下论文：Recovery of deleted deep sequencing data sheds more light on the early Wuhan SARS-CoV-2 epidemic（恢复删除的深度测序数据为新冠病毒在武汉的早期流行提供了更多信息） https://www.biorxiv.org/content/10.1101/2021.06.18.449051v1 https://www.biorxiv.org/content/10.1101/2021.06.18.449051v1.full 6park.com本文将说明以下6点： 1、Bloom恢复的基因序列，是武汉大学“纳米孔靶向测序检测”技术研发过程中产生的测试数据； 6park.com2、关于数据是如何删除的，还存在疑点，还不能下定论。一种可能是：武汉大学在完成相关项目后，向NIH请求撤回（删除或移走）了他们在研发、测试阶段输出的，提交到SRA保存的新冠病毒基因序列片断； 6park.com3、“纳米孔靶向测序检测”技术是qPCR病毒检测技术的替代技术，其核心功能是病毒检测，而非病毒测序。检测样本中是否存在目标病毒是这一技术的首要功能，对病毒进行部分测序并输出测序结果是其辅助功能。而且，这一技术测序、输出的不是新冠病毒的全基因组序列，只是与检测功能相关的全基因组序列的若干片断； 6park.com4、“纳米孔靶向测序检测”技术输出的基因序列，允许更大误差和更高错误率，不应将它们等同于专门测序得到的，可向国际生物数据库正式提交的，可用于严谨基因科学研究的，权威性的基因序列。即使这些基因序列片断仍存在科学研究价值，在将它们用于对准确性要求极高的基因科学研究，特别是新冠溯源研究时，也应当非常审慎，也应当对其中的测序错误可能产生的影响、误导进行充分的分析、评估。 6park.com5、Bloom恢复的序列片断证明，武汉大学“纳米孔靶向测序检测”技术的测序结果中存在大量的核苷酸测序缺失。 6park.com6、武汉大学“纳米孔靶向测序检测”技术输出的欠精确、可靠的研发、测试数据，被Bloom拿来充当了其溯源论文的基础数据。是否采取措施规避了基础数据中的测序错误？测序错误可能对研究过程及论文结论造成什么样的影响或误导？Bloom在论文中没有进行讨论、评估。 6park.com下面对上述6点予以展开。 6park.comBloom在论文中说，他通过谷歌云恢复了去年6月被删除的一些测序数据，他通过这些数据重建了13个新冠病毒的基因序列。论文提供了这13个基因序列的下载地址： https://github.com/jbloom/SARS-CoV-2_PRJNA612766/raw/main/results/consensus/consensus_seqs.csv 6park.com需要指出的是，Bloom重建的13个序列都不是完整的新冠病毒全基因组序列，每个序列都是新冠全基因组序列一大一小的两个片断。Bloom未得到全基因组序列，不是因为他恢复得不完全，而是因为他恢复的数据本来就是基因组片断。这些新冠基因组片断来自武汉大学“纳米孔靶向测序检测”（Nanopore Targeted Sequencing，简称NTS）项目，是这一项目研发、测试过程中产生、输出的。NTS技术不需要检测新冠病毒的全基因组，只需要比对、测序全基因组中的约12个区段。 6park.com后续内容将涉及一个生物信息数据库SRA。SRA，即Sequence Read Archive（序列读取档案）是NIH（National Institutes of Health，美国国立卫生研究院）管理、维护的两个生物学数据库之一，另一个是GenBank。 6park.comNTS项目使用SRA为研发的数据存储平台，它在SRA上的项目代号为Bio-Project PRJNA612766。Bloom论文中说，截止2020年3月30日，PRJNA612766项目共向SRA提交了282 份数据，其中241份数据的相关信息被一位名叫Carlos Farkas的科学家（等人）整理在一个Excel表格中，Bloom就是依据表格中的信息恢复了部分数据并重建了13个序列片断。Bloom恢复的数据只是PRJNA612766即武汉大学NTS项目提交数据的很小一部分。我粗略地浏览了该表格，NTS项目最早的数据提交时间应该是2020年1月15日。顺便附上Bloom论文提供的Farkas表格的下载地址。 https://web.archive.org/web/20210502130356/https://dfzljdn9uc3pi.cloudfront.net/2020/9255/1/Supplementary_Table_1.xlsx. 6park.comBloom重建的13个序列片断的对象时间戳都是2020年2月15日。Bloom推测，对象时间戳可能是指数据上传到SRA的时间。也就是说，Bloom恢复的数据应该是同一批数据，它们都是2020年2月15日这天上传到SRA的。 6park.comBloom自已通过谷歌云恢复、重建了13个序列片断。其实，他可以通过NIH，利用SRA的数据库备份系统恢复所有相关数据。论文没有提及这一数据恢复途径，没有解释Bloom为什么不通过这一途径完整地、更可靠地恢复数据。 6park.comBloom论文引用的武汉大学NTS（纳米孔靶向测序检测）研究的相关论文为：Nanopore target sequencing for accurate and comprehensive detection of SARS-CoV-2 and other respiratory viruses（纳米孔靶向测序可准确、全面检测 SARS-CoV-2 和其他呼吸道病毒） https://www.medrxiv.org/content/10.1101/2020.03.04.20029538v1.full-text 6park.com这一论文投稿到预印版平台medRxiv的时间是2020年2月29日，正式发布时间是2020年3月6日。 6park.com2020年3月4日，中国人民日报和新华网报道了这一新技术： 武汉大学研发纳米孔靶向测序检测方法 http://www.xinhuanet.com/science/2020-03/04/c_138841605.htm 6park.com概要地介绍一下武汉大学“纳米孔靶向测序检测”技术（NTS）的功能特点，这将有助于我们作出进一步的正确判断： 1、相比传统的qPCR（quantitative polymerase chain reaction，实时荧光定量PCR，或实时定量聚合酶链反应）检测30%~50%的阳性检出率，NTS将阳性检出率提升了43.8%，达到约75%~94%； 2、对于高浓度病毒样本，NTS仅需测序10分钟即可检测阳性，即使极低浓度病毒样本，也仅需测序4小时完成检测，从收到样本到出具结果，全程控制在6-10小时； 3、NTS可在测序后4小时内高敏感性、高准确性地同时检测SARS-CoV-2和其他10大类、40余种呼吸道病毒； 4、NTS最低检测敏感度是广泛使用的传统qPCR的100倍； 5、NTS还可输出检测样本中病原体（如新冠病毒）的基因组序列片断，可用于考察新冠病毒基因组的变异情况，监控病毒变异引起的毒性与传播能力改变。这是qPCR所没有的功能，qPCR病毒检测只作基因比对，不作基因测序，不记录、输出任何基因组序列。 6、NTS所需的纳米孔测序平台对实验室要求不高，其中最小测序仪MinION是便携式的，因此NTS也适合不同级别的医院使用。 6park.com由上述功能细节可知，NTS技术是面向、服务于医院和普通实验室的；NTS技术对实验室条件要求不高，这应该暗示着它并不能胜任有着极高精度的专业测序；第5组功能是锦上添花性质的，该组部分功能很可能有夸大其词的噱头成份。一个仓促研发的产品，如果号称拥有众多强大功能，那很可能意味着，至少，它的相当一部分功能是不尽善尽美的。 6park.com下面讨论数据撤回、删除的有关情况。 6park.comBloom在论文中说，当他按照Carlos Farkas论文和表格的指引去查阅武汉大学提交到SRA的数据时，在 NCBI-SRA系统中已经搜索不到相关项目PRJNA612766了，用Farkas提供的Accession ID直接搜索该项目下的相关序列数据，搜索结果提示：序列数据已被删除。 6park.comBloom说，SRA 被设计为深度测序数据的永久存档，上传数据到SRA后，只能通过向 SRA 工作人员发送电子邮件来删除数据。为说明可以通过发送邮件来请求删除数据，Bloom举了一个例子，他提供了一个Xiao姓科学家向SRA工作人员发送email，请求删除数据的电子邮件截图。 6park.com6park.comBloom说，这位xiao姓科学家是一篇穿山甲冠状病毒论文的lead author。Bloom论文的“Literature Cited”（引用的文献）中指示了xiao姓科学家的论文： https://www.nature.com/articles/s41586-020-2313-x 6park.com由论文作者列表可知，这个xiao姓的lead author或者是论文的通讯作者之一，华南农业大学特聘教授肖立华，或者是论文的第一作者Kangpeng Xiao。 6park.comBloom提供Xiao姓科学家的email截图是为了例证：SRA数据可通过email请求撤回、删除。由截图可见，Xiao姓科学家删除的项目与武汉大学的NTS项目二者的BioProject ID不同，它们是两个不同的项目。 6park.comBloom是怎么得到Xiao姓科学家邮件截图的？Bloom也未作说明。邮件截图是NIH提供的吗？如果是的话，NIH为什么不直接提供武汉大学PRJNA612766项目（即“纳米孔靶向测序检测”项目）相关的邮件，而要提供一个无关项目的邮件呢？如果邮件不是NIH提供的，那么，如此私密的信息，Bloom又是怎么得到的呢？Bloom和xiao姓科学家有什么直接或间接的关系吗？ 6park.com这些疑点暂时无法得到澄清。至少有以下两种可能性： 1、xiao姓科学家的邮件是NIH提供的，但NIH未提供或未能提供武汉大学PRJNA612766项目组人员请求删除数据的邮件； 2、xiao姓科学家的邮件是Bloom通过私人渠道获得的。在2018回国工作，被华南农业大学特聘前，肖立华在美国学习、工作了至少18年，在美国CDC工作了14年。 6park.com小结一下上述几段内容。Bloom称数据被删除了，称SRA数据只能通过发送email请求删除，他还举了一个其它项目通过这一途径请求删除数据的例子，但他未提供武汉大学项目组人员要求删除数据的直接证据。Bloom论文提供的信息，让我无法肯定，相关数据确实是武汉大学项目组人员要求删除的，我不能排除数据删除的其它可能性。 6park.com华尔街日报（The Wall Street Journal）6月28日的一篇文章含糊其辞、语焉不详地声称：NIH证实，应一名中国研究人员的申请删除了这些序列。这篇新闻的标题为：美国应中方要求删除新冠基因序列，病毒溯源难度加大 https://cn.wsj.com/articles/%E7%BE%8E%E5%9B%BD%E5%BA%94%E4%B8%AD%E6%96%B9%E8%A6%81%E6%B1%82%E5%88%A0%E9%99%A4%E6%96%B0%E5%86%A0%E5%9F%BA%E5%9B%A0%E5%BA%8F%E5%88%97%EF%BC%8C%E7%97%85%E6%AF%92%E6%BA%AF%E6%BA%90%E9%9A%BE%E5%BA%A6%E5%8A%A0%E5%A4%A7-11624509310 6park.com令人不解的是，对如此重要的事情，华尔街日报的文章竟不说明，NIH何时证实了相关信息，证实相关信息的是NIH的哪一位工作人员，其title是什么，证实信息的有关声明发表于何处。无法判定华尔街日报的做法是有意，抑或仅仅是大意疏忽。 6park.com尽管至今仍疑点重重，但我不想过多纠缠数据到底是怎么删除的，这不是第一位的问题。以下，我将假定相关数据确实被删除了，而且确实是武汉大学项目组人员请求删除的。 6park.comBloom提到了一个叫Aisu Fu的人，论文说，武汉大学PRJNA612766项目的病毒样本是Aisu Fu和武汉大学人民医院搜集的。论文没有提供此人的更多信息，事实上，Aisu Fu在论文中只出现了一次。 6park.comAisu Fu是谁呢？ 6park.comAisi Fu，中文名付爱思，是武汉臻熙医学检验实验室有限公司的总负责人，他与武汉大学药学院刘天罡教授，武汉大学人民医院李艳教授、余锂镭教授是NTS技术的共同研发者，这些信息可由以下新闻获得：武汉大学新闻网- 重磅！武汉大学联合团队开发纳米孔靶向测序 大幅提升新冠病毒阳性检出率 https://news.whu.edu.cn/info/1002/57753.htm 6park.com在前面提到过的武汉大学medRxiv预印本论文中，付爱思是第二作者，不过，在那篇论文中，他的署名不是Aisu Fu，而是Aisi Fu。 https://www.medrxiv.org/content/10.1101/2020.03.04.20029538v1.full 6park.com假定Bloom恢复的SRA数据确实是武汉大学请求删除的，那么，我认为，情况应该是：NTS（纳米孔靶向测序检测）项目研发完成后，武汉大学项目组人员向NIH-SRA发出申请，撤回、删除了他们提交到SRA的数据。这些数据本来就是项目研发、测试过程中产生的非正式数据。 6park.com这一撤回、删除非正式的研发、测试数据的做法有什么不正常吗？隐藏着什么不良动机吗？我看不出来。 6park.com学术论文应当专注于学术本身，基于事实，有一说一，有二说二，力求客观中立，避免被利用为政治工具，更不主动充当政治工具。妄加揣测，轻率贸然地陷人以罪，是有违科学精神，不道德，不负责任的行为。Bloom是怎么做的呢？在论文的Discussion部分，Bloom对中国科学家提出了如下指控：显然，对样本进行完全测序比偷偷删除部分序列更能提供科学信息。。。这些序列似乎很可能被删除以掩盖它们的存在。原文为：and it clearly would have been more scientifically informative to fully sequence the samples rather than surreptitiously delete the partial sequences。。。It therefore seems likely the sequences were deleted to obscure their existence. 6park.comBloom的这些言辞，是强词夺理、自相矛盾、逻辑错乱、无耻下作的造谣中伤。 6park.com首先，Bloom恢复的基因序列片断并非来自专门测序过程，它们是NTS（纳米孔靶向测序检测）技术研发、测试过程中产生的非正式数据，NTS主要用于医学检测、临床诊断，其输出的基因序列不是为严肃基因研究准备的。 6park.com第二，武汉大学在NTS研究中没有进行完全测序，不是他们故意不进行完全测序，而是因为NTS不需要进行全基因组测序，对基因组的某些重要片断进行测序就很充分了。事实上，NTS的序列检测范围，已经大大超过传统qPCR的序列比对范围，“相当于撒下了十几张大网”，同时捕捉病毒样本中的可疑基因片断。 这是NTS阳性检出率大大提高的根本原因。 6park.com 6park.com第三，Bloom的溯源研究论文基于一个最基本的假设：新冠病毒是自然演化产生的。如果这一假设不成立，Bloom的论文就崩溃了。Bloom一方面以武汉大学的数据作为自己自然演化理论的基础数据，基本依据，一方面又指控武汉大学的科学家偷偷摸摸删除这些可支持其自然演化理论的数据。这是一种非常错乱的逻辑。武汉大学的科学家为什么要删除“新冠自然演化”的证据？删除这些证据对中国科学家，对中国政府有什么好处？删除“新冠自然演化”的证据，掩盖“新冠自然演化”的“真相”对中国科学家，对中国政府有什么好处？ 6park.com第四，2017年12月19日，Trump政府解除了奥巴马政府3年前颁布的“功能增益研究”(Gain-of-Function, G-o-F)禁令，允许美国科学家重新申请联邦经费，开展功能增益研究，在实验室中研发、制造更具致病能力或更具传播能力的病毒或其它病原体。奥巴马禁令颁布于2014年10月22日。美国的功能增益研究与新冠病毒的出现没有关系吗？奥巴马功能增益研究禁令解除两年后，新冠病毒就出现了，这只是一种巧合吗？从2017年12月解除功能增益研究禁令，到2019年11月前后出现新冠病毒，在这近两年的时间里，美国的病毒学家们没有功能增益-改造出任何一种可怕的病毒吗？ 对于Trump政府解禁功能增益研究，打开了潘多拉灾难之盒，对于美国科学家可以在联邦经费支持下，合法地在实验室功能增益-改造病毒，Bloom不置一词，装作没有这回事。事实上，对解除奥巴马禁令，重启功能增益研究，美国科学界讳莫如深，全体失声；Bllom不仅不做深刻的自我反思、自我检讨，反而扮出清白无辜的模样对中国科学家造谣中伤，栽赃嫁祸，其言行非常无耻下作。 6park.com下面讨论Bloom所恢复数据的精确性、可靠性问题。 6park.com由研发目的、用途及功能特点可判断，NTS检测技术是传统qPCR检测的替代技术，它的首要目的是病毒检测，而非病毒测序。即，确定样本中是否存在目标病毒是它的首要功能，测序并输出基因序列片断只是它的辅助功能。NTS技术与专门的测序明显不同，它的测序精度难以与专门测序相提并论，它的测序结果无法直接提交到国际基因数据库作为权威数据供严谨的基因研究使用。将NTS输出的基因序列片断等同于专门测序得出的权威基因组序列，并将之用于对精确度、可靠度要求极高的新冠溯源研究，我认为是不恰当的。即使这些数据确有溯源研究的价值，如线索价值，在使用时也应当非常审慎、小心，避免被数据中的错误误导。 6park.com相比qPCR，NTS技术的阳性检出率虽大大提高，达到了约75%~94%，但仍远远称不上高度精准，可以想见，它的测序功能的精准程度也是有限的；同时，NTS是疫情发生后短期内开发出来的，数据上传的2月15日，项目刚刚启动了约一个月，NTS技术尚在研发阶段，其输出的基因序列更可能存在误差、偏差甚至错误。 6park.comNTS测序功能精确度欠佳。这不是一个推测，而是一个事实。Bloom恢复、重建的基因序列片断就是这一事实的确凿证明。 6park.comBloom在论文中说：I aligned the recovered deep sequencing data to the SARS-CoV-2 genome using minimap2。。。 即：我使用minimap2比对了恢复的深度测序数据与新冠病毒基因组的一致性。minimap2是一种基因组序列比对工具。 6park.com稍作间隔后，Bloom提供了如下表格： NTS多区段比对、检测与qPCR有限位置比对对照图
6park.com表格第一列的新冠病毒sample一共有14个，其中13个对应Bloom从SRA恢复、重建的基因序列片断，另外一个来自某个2月住院的患者。 6park.com表格中的第二列应该是14个sample与proCov2的基因序列一致性比对结果。proCov2是天普大学（Temple University）科学家Sudhir Kumar提出的一个虚拟的新冠病毒的祖病毒。proCov2与最早发现的新冠病毒样本之一WuHan-hu-1只相差三个核苷酸，将WuHan-hu-1进行以下三个单核苷酸的更改：C8782T、C18060T和 T28144C，就得到了proCov2。 6park.com注1：C8782T代表：将基因序列中8782位点的胞嘧啶C对应的核苷酸（碱基对）改为胸腺嘧啶T对应的核苷酸。 6park.com注2：WuHan-hu-1的基因序列是上海复旦大学张永振团队2020年1月5日上传的，是第一个上传到国际生物信息数据库的新冠病毒全基因组序列。WuHan-hu-1的病毒样本由武汉市中心医院采集提供，采集时间是2019年12月30日或26日，采集自一名41岁的陈姓新冠早期患者，该患者是华南海鲜市场的一名个体经营者。 6park.comproCov2与Wuhan-hu-1只相差3个核苷酸，而二者基因组序列（核苷酸序列）长度均为29903（含近3万个核苷酸）。易知，二者基因组序列的差异度约为0.01%（万分之一），即一致性约为99.99%。所以，Table-1中各病毒sample与proCov2的基因序列一致性比对结果，可视为这些sample与Wuhan-hu-1的一致性比对结果。 6park.com要注意的是，由于Bloom恢复的基因序列，也就是武汉大学NTS技术输出的基因序列不是全基因组序列，而只是全基因组序列的片断，因些，Table-1显示的不是病毒间全基因组序列的一致性，而是对应的基因组片断的一致性。Table-1标题栏已标明，比对的序列区间是21570~29550，这一区间含近8000个核苷酸。 6park.com观察各个病毒sample与Wuhan-hu-1或proCov2的一致性，可知：Table-1所列14个病毒sample中，8个与Wuhan-hu-1的一致性差别超过1%，6个超过2%，4个超过3%，3个超过4%。 6park.com由于比对的基因片断含近8000个核苷酸，因此，如果序列一致性差别超过1%，那就意味着该病毒Sampe与Wuhan-hu-1至少有80个核苷酸不同； 如果一致性差别超过2%，就至少有160个核苷酸不同； 如果一致性差别超过3%，就至少有240个核苷酸不同； 如果一致性差别超过4%，就至少有320个核苷酸不同。 6park.com比对区间外还有2万多个核苷酸位点，如果考虑这些位点，那么这些病毒sample与Wuhan-hu-1的核苷酸差异可能会更多。 6park.com新冠病毒的变异速度是：一个病毒平均一年产生约25个核苷酸突变。产生上述规模的突变，正常情况下需要几年、十几年的时间。在疫情早期的2月15日前，武汉大学人民医院的13位患者，其体内的新冠病毒同时发生了如此显著的突变，这可能吗？ 6park.com对比一下迄今为止，新冠病毒的实际变异情况： 英国Alpha变种B.1.1.7的核苷酸变异位点约为28~32个； 南非Beta变种B.1.351的变异位点约为23个； 巴西Gamma变种P.1的变异位点约为17个； 印度Delta变种B.1.617.2的变异位点约为13~17个。 6park.com可见，疫情发展至今，四个最重要的新冠病毒变种，其变异位点数都小于40个核苷酸；相比之下，Bloom恢复的数据有8个呈现出了超过80个核苷酸位点的不同，有3个呈现出了超过320个核苷酸位点的不同。而且，这么大幅度的核苷酸差异发生在去年2月15日前发现的病毒样本中。 6park.com这些核苷酸位点的差异都是突变造成的吗？如果是的话，如此超常的突变，Bloom怎么会视而不见呢？这个困惑昨天未能解决。今天，我返回Bloom的论文，下载了几个Bloom恢复、重建的基因序列，将它们的基因序列与Wuhan-Hu-1的基因序列加以比对（使用NCBI Blast工具），而后发现：原来，绝大多数核苷酸差异对应着核苷酸缺失，与变异无关。Bloom重建序列中存在着大量的字母N，每一个字母N都代表着其所在位点的核苷酸缺失。 6park.comBloom重建序列中的大量核苷酸缺失，要么是Bloom重建序列时产生的，要么是武汉大学NTS测序错误造成的。 6park.com对所采用的基础数据的准确性、可靠性问题，以及它们对研究过程、论文结论可能造成的影响，Bloom在论文中没有进行讨论和评估。 6park.com（正文完） 6park.com 6park.com相关文章： 6park.com虚假伪劣、贼喊捉贼，打假美国之音溯源新发现 https://club.6parkbbs.com/bolun/index.php?app=forum&act=threadview&tid=15950371 6park.com中国科学家是否删除了新冠早期基因组数据？ https://club.6parkbbs.com/bolun/index.php?app=forum&act=threadview&tid=15949751 6park.com 附录：Bloom重建的基因序列片断示例 https://club.6parkbbs.com/pk/index.php?app=forum&act=threadview&tid=14492427 6park.com该示例对应Table-1中的C9，即第四个病毒Sample，序列中的N代表核苷酸缺失的位点。该示例是C9基因序列片断的一部分，未展示全部近8000个核苷酸（碱基对）。 6park.comBloom论文提供了13个重建新冠基因序列片断的下载地址： https://github.com/jbloom/SARS-CoV-2_PRJNA612766/raw/main/results/consensus/consensus_seqs.csv 6park.com 6park.com贴主:苦难与荣耀于2021_07_02 21:47:59编辑 6park.com贴主:苦难与荣耀于2021_07_03 0:00:45编辑 6park.com贴主:苦难与荣耀于2021_07_03 23:11:56编辑 6park.com贴主:苦难与荣耀于2021_07_04 4:45:00编辑
贴主:苦难与荣耀于2021_07_05 5:24:54编辑
请点这里投票，“赞”助支持！