全基因组测序与质量控制
获取并筛选用于系统发育分析的Mab临床分离株数据。
从公开数据库及文献收集Mab WGS数据，进行质量控制和亚种分配。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
脓肿分枝杆菌临床优势克隆群的突变率与适应性变异
Mutation rates and adaptive variation among the clinically dominant clusters of Mycobacterium abscessus
实验完整度
高结合系统发育分析、溯祖分析、祖先模拟和体外突变率测定等多层级验证，支持突变率减慢假说。
主要模型
M. a. abscessus 临床分离株 M. a. massiliense 临床分离株 M. bolletii 临床分离株
重点核对
483例患者分离株的全基因组序列 使用Treemmer修剪系统发育树以保留95%的多样性 BEAST溯祖分析估计分子钟速率 fluctuation assay测定自发突变率 UvrD/Rep解旋酶基因的突变富集
摘要
脓肿分枝杆菌（Mab）是一种多重耐药病原体，越来越多地导致严重的肺部感染。对Mab全基因组序列（WGS）的分析表明，来自不同地理位置的临床分离株形成密集的遗传聚类。这被解释为支持患者间传播，但流行病学研究反驳了这一解释。在这里，我们提供了与系统发育聚类出现同时发生的Mab分子钟速率减慢的证据。我们使用公开可用的483个Mab患者分离株的WGS进行了系统发育推断。我们采用子采样方法结合溯祖分析来估计树的长内部枝上的分子钟速率，表明与系统发育聚类内的枝相比，长期分子钟速率更快。我们使用祖先模拟来预测时钟速率变异对系统发育聚类的影响，发现观察到的系统发育中的聚类程度更容易用时钟速率减慢来解释，而不是用传播来解释。我们还发现，系统发育聚类富含影响DNA修复机制的突变，并报告聚类分离株在体外的自发突变率较低。我们提出，Mab通过DNA修复基因的变异适应宿主环境，影响其突变率，并表现为系统发育聚类。这些结果挑战了Mab中的系统发育聚类可归因于人际传播的模型，并增进了我们对新兴兼性病原体传播推断的理解。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 脓肿分枝杆菌的系统发育聚类是否由突变率变化而非人际传播驱动？
- 核心机制
- DNA修复基因（如UvrD/Rep家族解旋酶）的突变导致突变率降低，从而在系统发育上表现为密集聚类。
- 主要证据
- 基于483个分离株的系统发育分析、溯祖分析估计的分子钟速率、祖先模拟以及体外突变率测定，显示聚类分离株的突变率较低。
- 研究意义
- 挑战了系统发育聚类作为人际传播证据的模型，有助于理解新兴兼性病原体的传播推断。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
系统发育树构建
确定Mab分离株的进化关系以及密集聚类的存在。
使用核心基因组SNP比对，通过IQ-TREE构建最大似然树，并鉴定优势克隆群(DCs)。
分子钟速率估计
估计长期（内部长枝）与近期（聚类内）的进化速率差异。
通过Treemmer修剪树，利用BEAST进行溯祖分析估计长枝的分子钟速率；通过患者内纵向分离株对估计聚类内的突变率。
祖先模拟检验聚类驱动因素
评估种群规模效应和突变率变化对系统发育聚类的解释能力。
使用msprime模拟祖先历史，设置聚类亚群与非聚类亚群，比较不同有效种群大小和突变率变化下的聚类程度。
鉴定富集突变与突变率测定
识别与聚类相关的遗传变异，并验证聚类菌株的突变率降低。
使用SNPPar识别独立突变，通过Fisher精确检验富集；选择代表性菌株进行fluctuation assay测定自发突变率。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 全基因组测序处理 | SRA工具包 | -- | v2.9.6 |
| 序列处理与质检 | Trimmomatic | -- | v0.36 |
| 质量控制 | FastQC | -- | v0.11.8 |
| 基因组组装 | SPAdes | -- | v3.11.1 |
| 亚种分配 | fastANI | -- | v1.2 |
| 序列比对 | BWA MEM | -- | v0.7.17 |
| 变异检测 | Pilon | -- | v1.23 |
| 变异过滤 | bcftools | -- | v1.10.2 |
| 核心基因组推断 | Roary | -- | -- |
| 可定位性计算 | GenMap | -- | v1.1 |
| 重组去除 | Gubbins | -- | v2.4.1 |
| 系统发育推断 | IQ-TREE | -- | v1.6.12 |
| 树木修剪 | Treemmer | -- | v0.3 |
| 溯祖分析 | BEAST | -- | v2.6.3 |
| 模型选择 | ModelTest-NG | -- | v0.1.6 |
| 进化分析 | MegaX | -- | -- |
| 重组推断 | ClonalFrameML | -- | -- |
| 祖先模拟 | msprime | -- | v1.0.2 |
| 数据处理 | Python statsmodels | -- | v0.12.2 |
| pingouin | -- | -- | |
| Biopython | -- | v1.79 | |
| 细菌培养 | Middlebrook 7H9肉汤 | BD Diagnostics | 271310 |
| 吐温-80 | Sigma-Aldrich | P1754 | |
| OADC添加剂 | VWR | 90000-614 | |
| Middlebrook 7H10琼脂 | BD Diagnostics | 262710 | |
| 突变率测定 | 硫酸阿米卡星 | Sigma-Aldrich | A2324 |
| 培养板 | 24孔培养板 | VWR | 10861-558 |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 全基因组测序与质量控制 | 样本来源、测序平台、覆盖度阈值（≥80% at ≥20×）、平均核苷酸一致性阈值（≥98%） 阅读提示：Methods: Data Sources and Description, Quality Control of Illumina Reads |
| 系统发育树构建 | 核心基因组定义、重组去除方法（Gubbins）、每位患者仅取一个分离株、聚类定义标准 阅读提示：Methods: Variant Calling and SNP Filtering, Core Genome Inference, Phylogenetic Tree Inference; Results: Phylogenetic Analysis Confirms... |
| 分子钟速率估计 | 修剪策略（95%多样性保留）、BEAST模型设置（松弛分子钟、树先验）、样本量及链长 阅读提示：Methods: Coalescent Analysis; Table 1 |
| 体外突变率测定 | 菌株选择（是否属于优势克隆群）、培养基组成、培养时间（6天）、阿米卡星浓度（300 μg/mL）、起始菌量（5000 cfu/孔）、重复数（n=3）、计算方法（Shinyflan/ML估计） 阅读提示：Methods: M. abscessus Culture, Fluctuation Assay |
| 祖先模拟 | 有效种群大小、突变率变化范围、采样时间分布、重复次数（1000） 阅读提示：Methods: Ancestry Simulations |