患者样本选择与深度测序
确保样本在感染时间上可比，并获取细菌序列数据以分析病毒群体
从2013-2016年在几内亚欧洲移动实验室（EMLab）收治的EVD患者中选取血液样本，根据症状出现时间和序列质量进行筛选，最终纳入38例存活和96例死亡患者。利用RNA-seq文库构建和HiSeq2500测序获得序列数据。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
埃博拉病毒病患者的优势病毒基因组序列周围的变异对病毒载量和结局的影响
Variation around the dominant viral genome sequence contributes to viral load and outcome in patients with Ebola virus disease
实验完整度
高包含患者转录组深度测序的群体遗传学分析、相关性和统计学检验，并在细胞培养中使用EBOV微型基因组系统进行功能验证（LSTOP和L3mut突变体的活性及VP35相互作用），形成多层次的分子和功能验证。
主要模型
急性EVD患者的血液样本 BSR-T7细胞（稳定表达T7聚合酶的BHK细胞系）
重点核对
患者分组：住院存活组和住院死亡组 患者样本症状出现后平均时间：存活组6.4天，死亡组5.9天 病毒载量：以1/Ct值表示 L蛋白突变体：LSTOP（Q986→stop）和L3mut（N572S，Q986R，F2061S） 微型基因组系统转染质粒用量：pMG 500ng，pNP 250ng，pVP30 125ng，pVP35 125ng，pL 125ng，pRLTK 250ng（每孔24孔板）
摘要
背景：病毒载量是埃博拉病毒病（EVD）患者结局的主要决定因素，高病毒载量往往导致死亡。2013-2016年埃博拉病毒（EBOV）暴发流行中的证据表明，不同的病毒基因型在患者中可表现出不同的表型。此外，由于病毒RNA合成在单个患者内的易错性，EBOV基因组以优势病毒基因组序列为中心存在。患者体内的低频率变异可能对病毒蛋白功能产生影响，从而影响总体表型。为了研究这些低频率变异的影响，对急性EVD患者的血液样本进行了深度测序。结果：我们比较了急性EVD存活患者与死亡患者之间的低频变异频率。发现了非同义变异，这些变异可能影响病毒蛋白功能。在L基因（编码病毒RNA依赖性RNA聚合酶（RdRp））的三个密码子位点发现了最高的替换频率。在病毒复制实验中重现这些替换，发现这些替换导致异常的病毒RNA合成，并与患者结局相关。结论：总之，这些发现支持以下观点：在EVD存活患者中，某些情况下病毒的遗传变异导致了影响病毒蛋白功能的有害突变，从而降低了病毒载量。此类突变也可能导致病毒持续存在，并与复发性感染相关。关键词：埃博拉病毒、埃博拉病毒病、病毒遗传学、进化、患者结局。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 埃博拉病毒（EBOV）在患者体内的低频率变异（minor variants）是否影响病毒载量和患者结局？
- 核心机制
- 低频率变异导致的L蛋白突变（N572S，Q986R，F2061S）以及Q986位点的提前终止密码子，削弱了L蛋白的RNA依赖性RNA聚合酶活性，并且可能通过与VP35等聚合酶辅因子结合而竞争性抑制功能，从而降低病毒RNA合成，导致病毒载量降低，有利于患者存活。
- 主要证据
- 对急性EVD患者血液样本的深度测序分析，发现L基因三个位点的低频变异与病毒载量负相关，且存活患者中更频繁；在细胞培养中利用EBOV微型基因组系统证实LSTOP无RdRp活性，L3mut活性显著降低，且LSTOP能与VP35相互作用。
- 研究意义
- 这些发现提示病毒基因组低频突变可能通过降低病毒复制能力而促进患者存活，并可能与病毒持续感染和复发性感染有关，为理解EVD的发病机制和抗病毒治疗（如favipiravir）的机制提供新见解。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
病毒基因组序列组装与变异分析
确定每个患者的优势病毒基因组序列，并量化基因组范围内的低频变异频率
将测序reads比对到人参考基因组和EBOV参考基因组，去除人源reads，再比对到患者特异的优势基因组序列，生成BAM文件，使用DiversiTools计算核苷酸和氨基酸位点的变异频率。
统计分析与核心位点关联鉴定
比较存活组和死亡组病毒载量及变异频率的差异，并鉴定与病毒载量相关的关键变异位点
采用t检验、Wilcoxon秩和检验、Spearman相关分析和广义线性模型（glm）分析病毒载量（1/Ct）与变异频率的关系。发现L蛋白572、986和2061位点的非同义变异与病毒载量负相关，且存活组中更频繁。
功能验证：微型基因组系统
验证L蛋白突变体（LSTOP和L3mut）对病毒RNA合成功能的影响
使用EBOV Makona株的微型基因组系统，将野生型L、LSTOP（Q986→stop）或L3mut（N572S，Q986R，F2061S）表达质粒转染BSR-T7细胞，检测荧光素酶活性。
蛋白质互作验证：共免疫沉淀
验证LSTOP突变体是否仍能与VP35相互作用，支持其可能通过竞争性结合抑制病毒RNA合成的模型
在BSR-T7细胞中表达VP35-eGFP和野生型L或LSTOP，用GFP-Trap琼脂糖珠子进行共免疫沉淀，通过Western blot检测相互作用。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
样本处理与测序文库构建
统计分析与数据可视化
功能验证实验
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| RNA-seq文库制备 | Epicentre ScriptSeq v2 RNA-Seq文库制备试剂盒 | Epicentre | -- |
| RNA-seq文库纯化 | AMPure XP磁珠 | -- | -- |
| 文库质量检测 | Qubit荧光计 | Thermo Fisher | -- |
| 安捷伦2100生物分析仪 | Agilent | -- | |
| 测序 | HiSeq2500测序仪 | Illumina | -- |
| 数据处理 | Cutadapt | -- | -- |
| Sickle | -- | -- | |
| Hisat2 | -- | -- | |
| Bowtie2 | -- | -- | |
| SAMtools | -- | -- | |
| Picard | -- | -- | |
| Quasirecomb | -- | -- | |
| DiversiTools | -- | -- | |
| 统计分析 | R语言 | -- | -- |
| ggplot2 | -- | -- | |
| 质粒构建 | pUC57载体 | -- | -- |
| pEGFP-N1载体 | -- | -- | |
| GeneArt基因合成服务 | Thermo Fisher Scientific | -- | |
| SURE感受态细胞 | Agilent | -- | |
| 亚克隆效率DH5α感受态细胞 | Invitrogen | -- | |
| QIAGEN质粒大提试剂盒 | Qiagen | -- | |
| NanoDrop ONE分光光度计 | Thermo Scientific | -- | |
| SacI-HF限制性内切酶 | New England Biolab | -- | |
| SalI-HF限制性内切酶 | New England Biolab | -- | |
| BglII-HF限制性内切酶 | New England Biolab | -- | |
| 细胞培养 | 杜氏改良Eagle培养基 | Sigma-Aldrich | -- |
| 胎牛血清 | Sigma-Aldrich | -- | |
| GlutaMAX添加剂 | Thermo Fisher | -- | |
| 青霉素-链霉素 | -- | -- | |
| 遗传霉素G418 | Thermo Fisher | -- | |
| Lipofectamine 2000转染试剂 | Invitrogen | -- | |
| 微型基因组系统 | 被动裂解缓冲液 | Promega | -- |
| 双荧光素酶报告基因检测 | Promega | -- | |
| GloMax Explorer板式发光检测仪 | Promega | -- | |
| 共免疫沉淀 | 磷酸钙 | -- | -- |
| HEPES缓冲盐水 (2x) | -- | -- | |
| Halt蛋白酶抑制剂混合物（无EDTA） | Thermo Scientific | -- | |
| GFP-Trap琼脂糖珠子 | ChromoTek | -- | |
| Western blot | Pierce BCA蛋白定量试剂盒 | Thermo Scientific | -- |
| Laemmli样品缓冲液 | Bio-Rad | -- | |
| 脱脂奶粉 | Marvel | -- | |
| 抗EBOV-L抗体 | IBT Bioservices | 0301-045 | |
| EBOV-NP抗体 | IBT Bioservices | 0301-012 | |
| EBOV-VP35抗体 | IBT Bioservices | 0301-040 | |
| 抗eGFP抗体 | Santa Cruz | sc-8334 | |
| 抗萤火虫荧光素酶抗体 | Abcam | ab185923 | |
| 抗GAPDH抗体 | Abcam | ab8245 | |
| 抗小鼠二抗 | Sigma | A4416 | |
| 抗兔二抗 | Sigma | A6154 | |
| Clarity Western ECL印迹底物 | Bio-Rad | -- | |
| ChemiDoc Touch凝胶成像系统 | Bio-Rad | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 样本选择与RNA-seq测序 | 患者症状出现到取样的平均时间，存活组和死亡组之间需无显著差异；序列质量过滤标准（组装长度>18800nt且无gap） 阅读提示：Results部分第一段 |
| 病毒基因组变异分析 | 参考基因组来源与版本（人类GRCh38和EBOV参考序列KY426690）；比对参数（--local -X 2000 --no-mixed）；变异调用流程（DiversiTools）和P值阈值 阅读提示：Methods部分：Dominant genome generation及Nucleotide and amino acid deviations |
| 病毒载量测定 | 病毒载量的表示方式（1/Ct值）及其与变异频率的相关性分析 阅读提示：Results部分（图1和图2） |
| 微型基因组功能性实验 | 细胞系（BSR-T7），转染试剂（Lipofectamine 2000），质粒用量（pMG 500ng/孔等），检测时间（转染后24小时），荧光素酶检测方法（Dual-Luciferase） 阅读提示：Methods部分：Cell culture及EBOV transcription/replication mini-genome system |
| 共免疫沉淀实验 | 转染方法（磷酸钙），质粒用量（pMG 17.5μg等），细胞数（1.5×10^6/盘），孵育时间（48小时），IP珠子（GFP-Trap）和Western blot抗体 阅读提示：Methods部分：Co-immunoprecipitation of VP35-eGFP |
| 统计分析 | 统计检验方法（t检验，Wilcoxon秩和检验，Spearman相关，glm），显著性阈值（P<0.05），glm模型族（Gamma） 阅读提示：Methods部分：Statistical analysis |