IF50.5｜古欧亚人群病原体时空分布图谱揭示人畜共患病起源

1,313例古人类宏基因组筛查显示，人畜共患病约6,500年前才出现，与畜牧化同步。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

传染病在人类历史中持续造成重大死亡，但病原体的起源、传播动态及其与人类生活方式变迁的关系仍不清楚。农业革命后，定居社区扩大、动物驯化普及，理论上增加了人畜共患病风险，即所谓“第一次流行病学转变”。然而，直接分子证据长期缺乏，该假说存在争议。古DNA技术虽已成功重建鼠疫、结核、天花等个别病原体基因组，但缺乏对古代细菌、病毒和寄生虫的整体景观研究。

科学问题

卡在哪里： 以往古病原体研究多为单物种、小样本，无法系统评估不同病原体类群在时空上的分布差异，也难以检验畜牧化对疾病负担的影响。

本文要回答： 利用大规模古人类宏基因组数据，绘制欧亚大陆病原体时空分布图谱，并探究人畜共患病出现的时间及驱动因素。

技术路线

作者主要用了：宏基因组分类、参考基因组比对、aDNA损伤认证、贝叶斯时间序列、时空建模。

作者开发了一套高通量计算流程，对1,313例古人类（距今37,000年）的鸟枪法测序数据进行微生物DNA筛选。首先用KrakenUniq进行k-mer分类，再对136个细菌/原生动物属和1,356个病毒属进行参考基因组比对，并通过metaDMG估计aDNA损伤模式来认证古代来源。最终得到5,486个物种-样本命中，其中3,384个涉及已知人类病原体。该流程结合了分类、比对和损伤认证，有效区分了内源性感染与环境污染。

图 1　Fig. 1。图1展示了1,313个样本的地理分布和年代范围，覆盖欧亚大陆西部、中亚和北亚及东南亚，时间跨度37,000年。柱状图显示不同地区样本来源的遗址类型（如墓葬、聚落等），其中牙齿样本占86%。该图奠定了研究的数据基础，显示样本在时空上的广泛覆盖，为后续病原体分布分析提供背景。

核心亮点

亮点 1｜古代微生物DNA广泛存在

在1,005个样本中检测到5,486个古代微生物命中，涉及492个物种、136个属。其中3,384个命中属于214种已知人类病原体。口腔微生物如Actinomyces（28.5%样本）和Streptococcus（18.1%）检出率最高，土壤相关属如Clostridium（18.9%）和Pseudomonas（8.3%）也常见。通过ANI和多位点率分析，发现口腔微生物和环境微生物呈现混合来源，而病原体如Y. pestis和M. leprae则表现为单一菌株感染特征。

图 2　Fig. 2。图2a显示不同属的古代微生物命中总数和检出率，细菌属如Actinomyces、Streptococcus、Clostridium等占主导。图2b展示ANI分布，Yersinia属内所有物种ANI均高，而Mycobacterium属中只有M. leprae的ANI>99%。图2c显示reads富集倍数，M. leprae富集超过100倍，而B. recurrentis富集最低。该图揭示了不同属内物种的相似性和丰度差异，反映了内源性与外源性来源的混合。

亮点 2｜人畜共患病约6,500年前出现

将古代微生物按传播途径和宿主分类后，发现人畜共患病类群在约6,500年前才首次被检测到，并在约5,000年前达到高峰，与牲畜驯化的广泛普及时间吻合。贝叶斯时间序列分解显示，人畜共患病检出率从约6,000年前开始上升，并持续到中世纪。相比之下，人际传播病原体在更早时期就已存在，但检出率从约2,500年前才开始增加。这为“第一次流行病学转变”提供了直接分子证据。

图 3　Fig. 3。图3分别展示了Y. pestis、B. recurrentis、M. leprae、T. pallidum、Y. enterocolitica、Leptospira等病原体的地理分布和时间线。每个面板中，白色方块为所有样本，彩色圆圈为阳性病例。例如，Y. pestis病例广泛分布于欧亚大陆，时间从约5,700年前到近代；B. recurrentis在铁器时代和维京时代集中出现。该图直观呈现了各病原体的时空动态。

亮点 3｜鼠疫耶尔森菌的时空动态

新发现42例Y. pestis感染（检出率3%），将鼠疫的时空分布大幅扩展。最早3例距今约5,700-5,300年，分布在俄罗斯西部、中亚和贝加尔湖地区，挑战了早期鼠疫仅为孤立人畜共患溢出的观点。时间序列分析显示，鼠疫检出率在约5,000年前达到高峰，随后在约2,800年前突然下降，进入约800年的检测空白期，直到约2,000年前再次出现，与历史记载的鼠疫大流行吻合。这种模式与LNBA-菌株灭绝和 flea 传播适应性进化相关。

图 4　Fig. 4。图4展示了Y. pestis和B. recurrentis检出率的贝叶斯变化点检测结果。上方面板显示检出率趋势线和95%可信区间，中间面板显示变化点的概率分布，下方面板显示趋势斜率为正（红色）、负（蓝色）或零（白色）的概率。Y. pestis在约5,000年前达到高峰，约2,800年前突然下降，约2,000年前再次上升；B. recurrentis则在约2,800年前后从零星出现转为持续检出。该图量化了两种病原体的流行模式转变。

亮点 4｜回归热螺旋体的流行史

新报告34例Borrelia recurrentis感染（检出率2.5%），最早病例来自新石器时代斯堪的纳维亚（距今约5,600年），表明体虱在农业社会早期已作为疾病媒介。检出率在铁器时代和维京时代达到高峰，且从约2,800年前开始持续检出，与Iron Age进化支的出现时间重叠。与鼠疫的检测空白期形成互补，说明该时期鼠疫缺席并非采样偏差。

图 5　Fig. 5。图5a显示不同微生物类别的相对检出率随时间的变化，口腔微生物始终占比最高，环境类微生物稳定存在，而人畜共患病类在约6,500年前才出现。图5b显示人畜共患病和人际传播病原体检出率的贝叶斯趋势线，人畜共患病从约6,000年前开始上升，人际传播病原体从约2,500年前开始增加。该图支持了农业和畜牧化导致疾病负担增加的结论。

亮点 5｜多种病原体首次在古代样本中发现

除鼠疫和回归热外，还首次在古代样本中检测到Yersinia enterocolitica（12例）、Leptospira interrogans（20例）、Corynebacterium diphtheriae（2例）等。疟疾方面，发现9例感染，涉及间日疟原虫、三日疟原虫和恶性疟原虫。此外，在维京时代个体中发现麻风分枝杆菌与天花病毒共感染，以及麻风与钩端螺旋体共感染等15例共感染案例。

亮点 6｜贝叶斯时空建模揭示驱动因素

利用分层贝叶斯模型分析病原体检出率与气候、人类迁移、祖先成分等协变量的关系，发现样本年龄是感染类病原体检出率的最重要预测因子，且呈负相关。经度对鼠疫和回归热检出率有正向效应，表明欧亚草原东部地区感染率更高，可能与畜牧生活方式有关。牙齿样本比岩骨更能保存病原体DNA，尤其对于口腔微生物和血源性病原体。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 宏基因组分类、参考基因组比对、aDNA损伤认证、贝叶斯时间序列、时空建模
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事古DNA、病原体进化、流行病学或生物信息学的研究者。

带走一句

大规模古宏基因组筛查结合严格aDNA损伤认证，可系统重建病原体时空动态，为流行病学转变提供直接证据。