基因组与生态学系统思维框架下的波多黎各致病性钩端螺旋体研究

《Frontiers in Public Health》：Genomic and ecological systems-thinking framework for pathogenic Leptospira in Puerto Rico

　　

摘要
引言：钩端螺旋体病是一种复杂的人畜共患病，需要高分辨率监测。研究人员采用系统思维框架，整合基因组与生态学数据，绘制了波多黎各共循环致病性钩端螺旋体谱系的地理和宿主结构。
方法：分析了来自人类、家养动物和野生动物宿主的问号钩端螺旋体（L. interrogans）、博氏钩端螺旋体（L. borgpetersenii）和基施纳钩端螺旋体（L. kirschneri）共44个核心基因组。利用根-尖回归、距离隔离分析和校准的单核苷酸多态性（SNP）阈值（≤1、≤5和≤10 SNPs）整合时空和景观元数据，以定义传播簇。
结果：钩端螺旋体物种表现出由地理分隔的不同生态途径，解释了问号钩端螺旋体基因组方差的56%和博氏钩端螺旋体的91%（PERMANOVA）。问号钩端螺旋体显示出跨多宿主的高景观连通性，形成局部网络（≤1至≤10 SNPs），捕获了活跃的溢出事件（≤1 SNP下的人-鼠关联），并解析为啮齿动物宿主特异性谱系（R2 = 0.34）。相反，博氏钩端螺旋体表现出空间和时间的基因组同质性，在以小家鼠（Mus musculus）为主的未分割传播库中缺乏宿主相关结构。因此，固定基因组阈值产生了不同结果：问号钩端螺旋体解析为4至5个离散的扩展簇，而博氏钩端螺旋体在≤10 SNP阈值下归为单一均匀种群。
结论：共循环的致病性钩端螺旋体占据由不同宿主限制和环境持久性塑造的独特生态位。固定基因组阈值缺乏普遍适用性；有效的基因组流行病学监测必须采用物种特异性阈值校准，以准确绘制传播路径。

基于基因组与生态学系统思维框架的波多黎各致病性钩端螺旋体研究解读

**研究背景与意义**

钩端螺旋体病是由钩端螺旋体属（Leptospira）螺旋体引起的一种全球重要人畜共患病，每年导致超过一百万临床病例和约58,900例死亡。该病的传播生态受环境因素、多种动物宿主以及钩端螺旋体物种基因组多样性的复杂相互作用驱动。钩端螺旋体属分为腐生型（S）和致病型（P）两个主要进化枝，其中致病型P1亚枝包含高毒力（P1+）和低毒力（P1?）病原体。钩端螺旋体的生态涉及一系列维持宿主（主要是啮齿动物，其维持慢性无症状感染）和意外宿主（包括人类，通常发展急性感染而不长期排菌）。历史上，确定感染源的方法依赖血清型宿主适应性的假设，但跨物种血清型共享和费力的分离株血清分型使传播途径的推断复杂化。全基因组测序（WGS）通过利用核心基因组单核苷酸多态性（SNP）比对来描绘高分辨率传播簇和跨物种溢出事件，有助于克服这些障碍。然而，由于元数据不可用或不完整、特定地点和时间的序列数量少以及固有的诊断挑战，钩端螺旋体的基因组流行病学研究仍然有限。此前的研究调查了城市温带气候中的单一宿主物种或特定暴发事件，而热带流行地区共循环钩端螺旋体物种的基因组多样性和时空聚类仍未得到充分探索。

**研究目的与方法**

本研究采用系统思维框架，分析波多黎各的致病性钩端螺旋体基因组。研究认识到监测缺口常常限制传统线性传播链的重建，因此将核心基因组结构解释为由宿主分布和景观组织塑造的生态连通性的直接体现。研究具体目标包括：（1）绘制钩端螺旋体多样性的多物种景观分区和宿主结构；（2）使用校准基因组阈值识别时空关联的传播簇；（3）弥合基因组连通性与景观尺度生态传播之间的差距。

**关键技术与方法**

研究人员使用了以下主要技术方法：从美国国家生物技术信息中心（NCBI）数据库筛选67个本地分离株的WGS元数据，排除宏基因组（n = 11）和重复记录（n = 2）后，对剩余54个分离株的原始读数通过Bactopia v3.1.0流程处理，使用FastQC和Trimmomatic进行质量控制和接头修剪，采用Shovill v1.1.0进行基因组组装，并依据五项质量标准评估。最终数据集包含44个基因组：问号钩端螺旋体（n = 24）、博氏钩端螺旋体（n = 15）和基施纳钩端螺旋体（n = 5）。核心基因组SNP比对通过Bactopia Snippy流程生成，使用Gubbins去除重组区域，IQ-Tree推断最大似然系统发育。时空和景观元数据整合采用根-尖回归、距离隔离分析（Mantel检验）、置换多元方差分析（PERMANOVA）、Hodges-Lehmann（HL）中位数差异对比以及固定SNP阈值（≤1、≤5、≤10 SNPs）定义传播簇。样本采集地点覆盖13个地理采样点，宿主包括小家鼠（n = 19）、褐家鼠（n = 10）、黑家鼠（n = 7）、人类临床病例（n = 4）、猫鼬（n = 1）和狗（n = 1）。该研究发表于《Frontiers in Public Health》。

**研究结果**

**3.1 数据集基线与同域种群结构**
44个致病性钩端螺旋体序列分布于13个地理采样位点，揭示了同域宿主重叠、单个ZCTA边界内的多物种共循环以及人类溢出与当地啮齿动物储库之间的共现。样本分为问号钩端螺旋体（24株）、博氏钩端螺旋体（15株，包括两个分歧的牛源谱系）和基施纳钩端螺旋体（5株）。

**3.2 根-尖回归分析与进化动态**
根-尖回归优化揭示了三种共循环钩端螺旋体物种的不同时间信号和对比微进化路径。问号钩端螺旋体在宏观时间线上呈线性时间信号（R2 = 0.23，替代率 = 3.466 × 10?2 替换/位点/年），在当代微观时间线上关联增强（R2 = 0.30），分子钟速率加快（5.55 × 10?1 替换/位点/年）。相反，博氏钩端螺旋体在宏观时间线上呈线性进展（R2 = 0.24），但微观时间线敏感性分析显示当代时间信号崩溃（R2 = 0.00），证实其局部宿主生态位内的短期时间基因组同质性。基施纳钩端螺旋体在当代窗口内表现出同步进化钟（R2 = 0.62）。

**3.3 问号钩端螺旋体基因组结果**
24个问号钩端螺旋体基因组的中位成对距离为24个SNPs（IQR: 15–28，范围: 0–37）。PERMANOVA确认采样ZCTA解释了总基因组方差的56%（F = 2.39; p = 0.002），Mantel检验确认了距离隔离模式（r = 0.19, p = 0.01）。地理分层评估显示空间衰减梯度：局部对（<5 km）中位数为19个SNPs，区域对（5–20 km）为20个SNPs，全岛对（>20 km）为24个SNPs。HL对比证实局部分离株相对于全岛对更接近中位数5个SNPs。固定阈值≤1、≤5和≤10 SNPs分别解析出4、5和5个不同簇。在≤1 SNP阈值下，一个人类临床病例（11-H10）与两个当地褐家鼠分离株在同一地点以0–1个SNP差异聚类，表明活跃的溢出事件。宿主物种来源解释了问号钩端螺旋体核心基因组SNP多样性的38%方差（PERMANOVA: R2 = 0.38, F = 2.19, p = 0.008）。在啮齿动物子集中，基因组多样性由特定啮齿动物物种结构化（PERMANOVA: R2 = 0.34; F = 4.32; p = 0.002），种内对比种间对更相关（HL位移 = 4 SNPs）。小家鼠分离株比褐家鼠（HL位移 = 7 SNPs）和黑家鼠（HL位移 = 10 SNPs）更保守。小家鼠内的SNP距离呈双峰分布，与宿主年龄相关：成年鼠对接近相同，而幼年鼠显示更高基线变异（HL位移 = 14 SNPs）。

**3.4 博氏钩端螺旋体基因组结果**
排除两个牛源分离株（分歧>40,000 SNPs）后，剩余13个基因组的中位成对距离为9个SNPs（IQR: 7–12）。ZCTA边界解释了总基因组方差的91%（F = 7.24; p < 0.001）。Mantel检验显示空间相关性（r = 0.22, p = 0.047），但刀切敏感性分析显示该连续信号脆弱。地理尺度揭示了二元分割而非渐进衰减梯度：局部对中位数为3个SNPs，全岛对为10个SNPs。固定阈值≤1、≤5和≤10 SNPs分别解析出2、2和1个簇。在≤5 SNP阈值下，2020年人类临床病例在扩大队列中连接，空间分类为“Distant”。在≤10 SNP阈值下，所有分离株包含在单个簇中。宿主物种身份占核心基因组SNP方差的26%，但关联不显著（PERMANOVA: F = 1.75, p = 0.13），种内对与种间对的基因组距离相当（HL位移 = 1 SNP）。小家鼠分离株对维持9个SNPs的基线中位距离，确认缺乏宿主物种间结构分离。

**3.5 基施纳钩端螺旋体基因组结果**
5个基施纳钩端螺旋体基因组来自四个采样点的两种宿主物种，中位成对距离为23个SNPs（IQR: 21–25，范围: 12–28）。在S11位点共定位的两个人源和鼠源分离株相差12个SNPs，表明这些啮齿动物可能是意外宿主或共享未采样主要储库的指示物。

**讨论与结论**

**讨论部分总结**：比较根-尖时间动态显示，问号钩端螺旋体维持活跃的当代分子钟（R2 = 0.30），而博氏钩端螺旋体表现出时间基因组同质性（R2 = 0.00）。博氏钩端螺旋体时间信号在移除2015年历史锚定分离株后完全崩溃，表明标准分子钟模型可能因杠杆伪影产生假阳性线性进化速率。景观连通性分析显示，ZCTA边界对两种物种基因组方差的解释力差异显著（56%对91%），反映问号钩端螺旋体的高连通性种群与博氏钩端螺旋体的碎片化地理亚群结构。固定基因组阈值缺乏普适性：对问号钩端螺旋体，≤10 SNP的放宽阈值可捕获局部时空关联网络；但对博氏钩端螺旋体，整个种群在≤10 SNP阈值下归为单一未分化簇，可能掩盖活跃传播动态。文献记载博氏钩端螺旋体经历插入序列介导的基因组缩减，失去环境感知和营养转运相关基因，促使其依赖直接宿主间传播循环；而问号钩端螺旋体保留大型结构完整基因组，支持环境持久性和多宿主地理传播。宿主物种与种群统计关联显示，问号钩端螺旋体的啮齿动物种群结构与基因组分配相关，种内传播模式与宿主行为生态学一致；博氏钩端螺旋体则缺乏宿主物种分配，作为由小家鼠主导的共享库运作。

**研究结论翻译**：本研究表明，波多黎各共循环的致病性钩端螺旋体物种占据独特的生态和进化生态位。问号钩端螺旋体的传播受高环境持久性、结构景观连通性和清晰宿主物种分配的影响。相反，博氏钩端螺旋体依赖高度均匀、宿主适应的克隆库，其特征是小哺乳动物种群内的短期基因组同质性。基因组监测框架不能依赖单一的通用属级传播阈值；必须进行生态分层和物种特异性校准，以准确绘制时空传播图并减轻人畜共患病风险。