《Nature Communications》：Reply to: Evidence of widespread sequencing errors in global avian influenza study
编辑推荐：
研究人员针对全球禽流感（Avian Influenza Virus, AIV）监测中出现的序列异常问题作出回应：异常源于向GISAID数据库批量提交过程中的操作失误。依GISAID建议，研究人员使用命令行工具FluCLI上传数据；在文件重构时，本用于系统发育树
研究人员针对全球禽流感（Avian Influenza Virus, AIV）监测中出现的序列异常问题作出回应：异常源于向GISAID数据库批量提交过程中的操作失误。依GISAID建议，研究人员使用命令行工具FluCLI上传数据；在文件重构时，本用于系统发育树构建、混合了研究序列与参考序列的分析文件被误传，且样本标识符被错误格式化，导致研究序列与参考序列重复、混杂。原始分析流程中并未出现此类异常。误传序列已在GISAID与GenBank按原 accession 号更正（版本2）。基于校正后序列，血凝素（Hemagglutinin, HA）与神经氨酸酶（Neuraminidase, NA）不再呈现跨地域、跨时段的100%核苷酸一致性；参照AIV substitution rate（HA约3.92×10-3 每位点每年，NA约3.61×10-3 每位点每年），校正序列符合流感病毒预期进化速率。TempEst与最大似然法（General Time Reversible模型）系统发育分析显示，校正后序列无离群节点，且与澳大利亚2020年报道序列非完全同源。研究不质疑澳大利亚无禽流感状态，亦不支持将任一国家认定为病毒起源地。巴布亚新几内亚（Papua New Guinea）地处东亚–澳大拉西亚迁飞廊道（East Asian–Australasian Flyway），可经西伯利亚、阿拉斯加等路径与其他大陆廊道间接相连；气候胁迫、监测空白与小型散养禽业导致病毒可沉默循环。H3分支注释错误源自NextClade自动归类，研究人员已将相关序列订正为“unassigned”。样本为环境排泄物，经DNA/RNA Shield灭活保存，不涉及活体操作与病毒培养，伦理获豁免。综上，公开记录异常为提交错误所致，研究结论基于真实数据成立。
研究背景与问题提出
禽流感（Avian Influenza Virus, AIV）作为具有大流行潜力的RNA病毒，其真实传播格局长期受限于全球监测不均：高收入国家具备较完整基因组监测体系，而全球南方（Global South）偏远岛屿、热带海岸与冲突地区长期缺乏系统采样。野生鸟类粪便作为环境中可被动收集的无创样本，理论上可弥补传统活体采样难以覆盖的生态盲区。然而，当研究人员试图基于大范围环境废物宏基因组揭示AIV空间扩散时，公开数据库中出现“跨大洲、跨年份100%一致”的异常序列，引发同行对数据可靠性的质疑。争议核心并非仅限技术失误，更触及迁飞生态学解释、国家溯源污名化与遗传资源合规边界。该回应文章发表于《Nature Communications》，是对Matters Arising质疑的正式答复与记录校正。
主要技术方法概览
研究人员依托2018年启动的环境废水与废弃物的监测计划，在亚洲、太平洋、非洲及中东多国偏远区域收集野生鸟类粪便环境样本，样本置入DNA/RNA Shield保存液以灭活病毒并保留核酸；核心分析包括：基于GISAID/GenBank的序列下载与BLAST比对、最大似然法（General Time Reversible模型）系统发育重建、TempEst时间信号检验、NextClade自动分型与人工校正，以及迁飞廊道（flyway）生态学交叉验证。研究不涉及病毒培养、活体接触或质粒操作，仅以环境核酸测序与生物信息学推断为技术主干。
研究结果解读
序列误传导致100%同源性假象
研究人员确认，异常100%重复来自提交环节：用于系统树预分析、含参考序列的混合文件被FluCLI误传，样本编号重排使研究序列与公共参考序列错位叠加。以校正后序列重新计算，HA与NA区段突变负荷符合Chen与Holmes（2006）所估AIV进化速率（HA约6.7个突变/年，NA约5.1个突变/年），跨区静止演化假象消失。
TempEst与时间系统发育再分析
下载1996–2024年亚洲NA全长序列817条，以最大似然法建树，校正序列未产生偏离时间梯度的离群分支，说明原结论所依数据在本地分析阶段并未嵌入错误参考，异常仅存于对外发布库。
澳大利亚状态澄清
系统树显示校正序列与澳洲2020年两条序列不完全等同。研究人员明确：原论文未声称澳洲为起源地，亦不质疑其无疫状态；流感病毒共享变体可因监测缺口在多地“沉默循环”，溯源归因需谨慎。
迁飞廊道与全球循环
巴布亚新几内亚位于东亚–澳大拉西亚迁飞廊道，并可经北西伯利亚、阿拉斯加链接中亚、西非–东非廊道；太平洋与美洲廊道在南北半球存在重叠。研究区含大量无人监测岛礁，叠加气旋、干旱、海平面上升与COVID-19期间监测转移，使传统路线扰动、病毒隐匿扩散具备生态学合理性。
H3分支注释问题
NextClade对禽源H3偶发归入人源分支，研究人员承担校验责任，将相关条目在作者勘误中改为“unassigned”并保留原始进化关系不确定性。
伦理与监管合规
粪便属公共环境废弃物，不含活体、不培养、不涉CITES标本或Nagoya Protocol遗传资源转移；泰国及多数采样国未批准名古屋议定书，本地权威确认无需特殊许可，机构伦理委员会给予豁免。
讨论与结论翻译
讨论指出，全球流感认知受采样偏见主宰：北方高纬度与富裕地区数据稠密，而岛弧、撒哈拉以南与南亚散户禽业近乎黑箱。校正后数据支持原研究主旨——偏远生态节点可能是AIV无声储存库，但绝不指向单一国家起源。研究人员强调，在缺少全物种、全地域覆盖前，任何“不可能传播”断言都缺乏稳健性。结论部分可译为：所识别异常皆因GISAID误传所致；研究人员承担全部责任，已撤回并替换受影响记录，强化内部提交流程，原 manusc 结论基于正确数据有效；禽流感全球循环之复杂性要求更公平、更密集、跨政治边界的环境监测网络，且学术表达须避免将遗传相似性强行等同于地理起源。
现有解读已完成，接下来可以选择一个方向继续深化：
把校正前后数据差异整理成评审要点 评估该回应的学术严谨性与遗留疑点 把校正前后数据差异整理成评审要点 评估该回应的学术严谨性与遗留疑点
