流行病学监测与地理脱钩分析
明确中国SFTS发病率与病死率的地理分布模式。
收集2010-2023年全国27省27457例报告病例和1326例死亡数据，计算累积发病率和CFR，比较内陆与沿海地区及河南与浙江省。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
东亚大别班达病毒的系统地理流行病学：分化的传播网络与基因型相关的临床严重性
Phylogeographic epidemiology of Dabie bandavirus in East Asia: divergent transmission networks and genotype‑linked clinical severity
实验完整度
高整合1820个病毒基因组、936例临床队列与27457例监测数据，含系统地理推断、重组检测和突变-表型关联分析。
主要模型
中国信阳市SFTS临床队列（805例） 浙江省沿海SFTS临床队列（131例） 1820个SFTSV全球基因组数据集 中国2010-2023年全国SFTS监测病例（27457例）
重点核对
病例定义：急性发热≥38.0°C，伴胃肠不适或出血等症状，血小板减少和白细胞减少，SFTSV RNA阳性或病毒培养分离 病毒RNA提取：200 μl血清使用QIAsymphony Virus/Bacteria Mini Kit，RT-PCR使用ZJ Bio-Tech试剂盒，ABI 7500系统，Ct≤37为阳性 基因组测序：巢式RT-PCR使用PrimeScript One Step RT-PCR Kit（RR055A）和Premix Taq Kit（RR902A），ABI 3730xl测序 贝叶斯系统地理分析：S片段207条序列，GTR+Γ+I模型，对数正态松弛分子钟，5次独立MCMC各1亿代，ESS>200 突变-表型关联：774个有结局注释的基因组，多变量logistic回归调整年龄、性别和发病至入院间隔，Benjamini-Hochberg校正
摘要
背景：由大别班达病毒（SFTSV）引起的严重发热伴血小板减少综合征（SFTS）在东亚呈现发病率与病死率地理脱钩的模式。我们旨在阐明导致这种内陆-沿海流行病学分化背后的不同生态驱动因素和系统地理动力学。方法：整合1820个高质量全球SFTSV基因组、特征明确的临床队列（936例患者）和中国全国监测数据（27457例病例），构建了综合分析框架。采用生态建模、贝叶斯系统地理学和基因型-表型关联分析，追踪病毒的进化轨迹和临床意义。结果：识别出显著的SFTS“内陆高发病率 vs. 沿海高病死率”模式。SFTS发病率对气象因素表现出不同的敏感性；内陆传播对温度波动敏感，而沿海动态受日照阈值（>200 h/月）约束。相比之下，临床严重性的空间分化与区域病毒遗传结构的分布相关。内陆地区主要共循环基因型A、C和D，而沿海地区以基因型B为主。浙江省被确定为遗传枢纽，其重组频率显著高于内陆地区（11.0% vs. 3.5%，P < 0.001）。贝叶斯系统地理推断表明中国浙江省与韩国和日本之间存在频繁的谱系交换。临床上，基因型B和D分别与沿海和内陆地区死亡率升高相关，提示沿海重症表型由其基因型B主导的结构所塑造。此外，RdRp-N828S突变成为致死结局的稳健分子相关因素，值得进一步功能验证。结论：不同的气象因素和可能的海洋传播网络可能构成SFTS地理脱钩流行病学的基础。这些发现强调风险评估必须超越发病率本身，并为高风险热点地区的靶向监测和基因型特异性干预提供了系统地理信息框架。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 阐明东亚SFTS发病率与病死率地理脱钩（内陆高发病率 vs. 沿海高病死率）背后的生态驱动因素和系统地理动力学。
- 核心机制
- 内陆发病率主要与温度等宏观气象因素相关，沿海严重性则与基因型B主导的病毒遗传结构及可能的海洋传播网络相关；基因型B和D分别与沿海和内陆死亡率升高相关，RdRp-N828S突变与致死结局独立相关。
- 主要证据
- 基于1820个SFTSV基因组、936例临床队列和27457例监测数据的分析显示：内陆与沿海发病率及CFR差异显著；GAM模型解释内陆82.1%、沿海48.7%的发病率变异；浙江省重组频率显著高于内陆（11.0% vs. 3.5%，P<0.001）；贝叶斯系统地理推断显示浙江与韩国、日本间频繁谱系交换；多变量logistic回归显示RdRp-N828S与致死风险增加相关（OR=2.52，校正P<0.05）。
- 研究意义
- 提示风险评估需超越发病率，为高风险热点地区的靶向监测和基因型特异性干预提供系统地理信息框架。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
气象因素与发病率关联建模
探究内陆和沿海地区气象因素对SFTS发病率的不同影响。
选取信阳市（内陆）和舟山市（沿海）为代表，收集2012-2018年月度平均气温、降水量和日照时数，采用广义加性模型（GAM）分析滞后一个月的气象因素与SFTS发病率的关系。
病毒基因组系统发育与基因型地理分布
表征SFTSV遗传多样性、基因型分类及内陆与沿海地区的基因型分布差异。
从GenBank获取1820个高质量SFTSV基因组，进行各片段多序列比对和最大似然系统发育树构建，采用MDS验证基因型分类，层次聚类分析区域基因型组成。
重组与重配检测
评估浙江省与内陆地区及韩国的重组和重配频率差异。
通过比较各毒株L、M、S片段基因型分配识别重配体，使用RDP4软件套件检测重组事件。
贝叶斯系统地理推断
推断浙江省在中国、韩国和日本间海洋传播网络中的作用及病毒扩散路径。
选择S片段进行离散系统地理分析，去除重组序列，对序列进行降采样，使用BEAST进行贝叶斯推断，重建MCC树和空间传播路径。
基因型-临床结局关联分析
检验病毒基因型与临床严重性（病死率）的关联。
比较河南（805例）和浙江（131例）临床队列特征，分析各基因型CFR，多变量logistic回归和Kaplan-Meier生存分析。
突变-表型关联与选择压力分析
鉴定基因型B特异性标记和与致死结局相关的突变，评估选择压力。
使用774个有结局注释的基因组，多变量logistic回归分析突变与致死结局的关联，HyPhy进行正选择分析和BUSTED-PH比较基因型B与其他基因型的选择压力差异。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
重组与重配检测
宿主-地点网络分析
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 病毒RNA提取 | QIAsymphony病毒/细菌迷你试剂盒 | Qiagen | -- |
| 实时RT-PCR | SFTSV逆转录聚合酶链反应试剂盒 | ZJ Bio-Tech | -- |
| ABI 7500实时PCR系统 | Applied Biosystems | -- | |
| 巢式RT-PCR | PrimeScript一步法RT-PCR试剂盒 | TaKaRa Bio Inc. | RR055A |
| Premix Taq™试剂盒 | TaKaRa Bio Inc. | RR902A | |
| Biometra TOne PCR系统 | Analytik Jena | -- | |
| PCR产物纯化 | QIAquick PCR纯化试剂盒 | Qiagen | -- |
| 测序 | ABI 3730xl DNA分析仪 | Applied Biosystems | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 临床病例定义与检测 | 病例定义：急性发热≥38.0°C，伴临床症状和血小板减少、白细胞减少，SFTSV RNA阳性或病毒培养分离；RNA提取200 μl血清；实时PCR Ct≤37为阳性 阅读提示：Methods 'SFTS clinical cohorts and SFTSV detection' |
| 基因组测序 | 巢式RT-PCR引物对（Table S1）；第一轮使用PrimeScript One Step RT-PCR Kit (RR055A)，第二轮使用Premix Taq Kit (RR902A)；循环条件：94°C 45s, 56°C 45s, 72°C 90s，30循环；ABI 3730xl双向测序 阅读提示：Methods 'SFTSV genome sequencing and phylogenetic inference'，Table S1 |
| 系统发育分析 | MAFFT v7.407比对，TrimAl v1.4修剪；IQ-TREE v1.6.12，GTR+F+R（L和M）、TVMe+R（S），1000次超快自举；参考序列GenBank NC_018136–NC_018138 阅读提示：Methods 'SFTSV genome sequencing and phylogenetic inference' |
| 重组检测 | RDP4软件，7种算法（RDP, GENECONV, Bootscan, Maxchi, Chimaera, Siscan, 3Seq），至少4种方法支持，P<10^-6 阅读提示：Methods 'Detection of reassortment and recombination' |
| 贝叶斯系统地理推断 | S片段207条序列；GTR+Γ+I替换模型，对数正态松弛分子钟，非对称替换模型+BSSVS；5次独立MCMC各1亿代，每10000步采样，10% burn-in；ESS>200；BF>3且PP>0.50 阅读提示：Methods 'Spatiotemporal phylogeographic inference'，Fig. S1, Fig. S2 |
| 突变-表型关联 | 774个有结局注释基因组；多变量logistic回归调整年龄、性别、发病至入院间隔；Benjamini-Hochberg FDR校正；基因型B特异性标记：>90%基因型B序列且<10%其他基因型序列 阅读提示：Methods 'Mutation-phenotype association and evolutionary pressure' |
| 选择压力分析 | Codon-MSA生成密码子感知比对；排除重组株和提前终止密码子；冗余序列折叠；最终数据集：RdRp 728，Gn 615，Gc 596，NSs 509，NP 393；HyPhy v2.5.24，FEL/MEME/SLAC/FUBAR，至少2种方法（P<0.05或后验概率>0.90）；BUSTED-PH和contrast-FEL 阅读提示：Methods 'Mutation-phenotype association and evolutionary pressure' |
| 气象因素与GAM建模 | 信阳市和舟山市；2012-2018年月度平均气温、降水量、日照时数；1个月滞后；GAM公式含年分类变量；薄板回归样条；诊断检查基维数k 阅读提示：Methods 'Meteorological data source and statistical analysis' |