以下是 CheckM 结果表中各列的中文解释:
-
Bin Id (基因组分箱编号)
- 这是每个组装出的基因组(MAG)的唯一标识符。在你的数据中,它们以 “RKP” 开头(例如 RKP53, RKP5)。这对应于你输入文件中的样本名称。
-
Marker Lineage (标记谱系/参考进化分支)
- 这表示 CheckM 用来评估该基因组质量的参考进化树分支。
- CheckM 不会用同一套标准去衡量所有细菌,而是根据该基因组最可能属于的分类群(如“根节点 root”、“细菌界 kBacteria”、“葡萄球菌属 gStaphylococcus”等),选择一套特定的单拷贝标记基因集来进行比对。
- 注意: 如果显示为
root (UID1)或k__Bacteria,通常意味着该基因组质量较差或太短,CheckM 无法将其精确归类到更具体的分类层级,只能用最通用的标记集来估算,因此其结果的参考价值较低。
-
# Genomes (# 参考基因组数)
- 在选定的“标记谱系”中,包含的参考基因组数量。
- 这个数字越大,说明该分类群在数据库中越丰富,CheckM 评估时使用的统计模型就越可靠。例如
g__Staphylococcus有 60 个参考基因组,而root有 5656 个。
-
# Markers (# 标记基因数)
- 在该特定谱系中,用于评估的单拷贝标记基因的总数。
- CheckM 通过查找这些保守基因是否存在、是否唯一来判断组装质量。不同分类群的标记基因数量不同(例如葡萄球菌属有 773 个标记,而根节点只有 56 个)。
-
Completeness (%) (完整度)
- 估计该基因组包含了多少比例的预期标记基因。
- 含义: 数值越高越好。100% 表示找到了所有预期的单拷贝标记基因。
- 警示: 高完整度并不总是代表高质量。如果 contamination(污染率)也很高(如你的 RKP53 达到 100% 完整度但 104% 污染),说明这个 bin 里混入了太多其他物种的序列,导致标记基因重复出现,从而虚高了完整度。
-
Contamination (%) (污染率)
- 估计该基因组中多余/重复的标记基因比例。
- 含义: 数值越低越好。理想情况下应为 0%。
- 计算逻辑: 如果某个本该是“单拷贝”的基因在你的 bin 里出现了 2 次或更多,CheckM 就会认为这是污染(即混入了其他菌株或物种的 DNA)。
- 你的数据情况: 很多样本(如 RKP53, RKP43, RKP31)污染率超过 100%,这意味着平均每个标记基因都出现了两次以上,这是严重的嵌合体(Chimeric bin),必须剔除。
-
Strain Heterogeneity (%) (菌株异质性)
- 估计该基因组内部是否存在多个近缘菌株混合的情况。
- 含义: 数值越低越好。
- 原理: 当存在多个相似菌株时,某些标记基因可能会出现轻微的序列变异(多态性),或者部分标记基因出现双份。CheckM 通过分析这些信号来估算异质性。
- 影响: 高异质性(如 RKP53 的 82.76%)意味着这个 bin 不是单一纯培养物,而是多个相似菌株的混合物。这对于构建精确的系统发育树或进行 SNP 分析是非常不利的。
Table 1: CheckM Quality Assessment Results
| Bin Id | Marker Lineage | # Genomes | # Markers | Completeness (%) | Contamination (%) | Strain Heterogeneity (%) |
|---|---|---|---|---|---|---|
| RKP53 | root (UID1) | 5656 | 56 | 100.00 | 104.17 | 82.76 |
| RKP5 | c__Gammaproteobacteria (UID4445) | 228 | 583 | 100.00 | 0.14 | 0.00 |
| RKP43 | root (UID1) | 5656 | 56 | 100.00 | 104.17 | 0.00 |
| RKP42 | k__Bacteria (UID203) | 5449 | 104 | 100.00 | 99.84 | 0.00 |
| RKP37 | k__Bacteria (UID203) | 5449 | 104 | 100.00 | 99.84 | 0.00 |
| RKP33 | k__Bacteria (UID203) | 5449 | 104 | 100.00 | 101.57 | 0.00 |
| RKP31 | root (UID1) | 5656 | 56 | 100.00 | 104.92 | 79.03 |
| RKP1 | root (UID1) | 5656 | 56 | 100.00 | 100.00 | 0.00 |
| RKP47 | g__Burkholderia (UID4006) | 64 | 769 | 99.87 | 0.00 | 0.00 |
| RKP38 | g__Staphylococcus (UID294) | 60 | 773 | 99.81 | 0.00 | 0.00 |
| RKP30 | g__Staphylococcus (UID294) | 60 | 773 | 99.81 | 2.61 | 7.14 |
| RKP54 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.00 | 0.00 |
| RKP25 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.28 | 0.00 |
| RKP24 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.28 | 0.00 |
| RKP19 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.00 | 0.00 |
| RKP18 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.00 | 0.00 |
| RKP17 | g__Staphylococcus (UID294) | 60 | 773 | 99.67 | 0.00 | 0.00 |
| RKP26 | o__Lactobacillales (UID544) | 293 | 475 | 99.63 | 0.00 | 0.00 |
| RKP15 | o__Lactobacillales (UID544) | 293 | 475 | 99.63 | 0.56 | 0.00 |
| RKP7 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP56 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP50 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 25.90 | 24.64 |
| RKP49 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP45 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP40 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP39 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP3 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP29 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP22 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP21 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP13 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP12 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP11 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP10 | g__Staphylococcus (UID298) | 56 | 805 | 99.62 | 0.00 | 0.00 |
| RKP36 | g__Staphylococcus (UID298) | 56 | 805 | 99.61 | 3.41 | 30.56 |
| RKP35 | g__Staphylococcus (UID301) | 45 | 940 | 99.51 | 0.08 | 0.00 |
| RKP34 | g__Staphylococcus (UID301) | 45 | 940 | 99.51 | 0.08 | 0.00 |
| RKP6 | g__Staphylococcus (UID298) | 56 | 805 | 99.51 | 1.79 | 57.14 |
| RKP4 | g__Staphylococcus (UID298) | 56 | 805 | 99.51 | 0.17 | 100.00 |
| RKP9 | g__Staphylococcus (UID298) | 56 | 805 | 99.48 | 0.00 | 0.00 |
| RKP8 | g__Staphylococcus (UID298) | 56 | 805 | 99.48 | 0.00 | 0.00 |
| RKP41 | c__Bacilli (UID285) | 586 | 325 | 99.45 | 0.55 | 0.00 |
| RKP2 | g__Staphylococcus (UID301) | 45 | 940 | 99.40 | 0.10 | 0.00 |
| RKP48 | g__Staphylococcus (UID298) | 56 | 805 | 99.38 | 1.14 | 50.00 |
| RKP46 | g__Staphylococcus (UID298) | 56 | 805 | 99.38 | 1.14 | 50.00 |
| RKP16 | g__Staphylococcus (UID298) | 56 | 805 | 99.38 | 0.59 | 0.00 |
| RKP14 | g__Staphylococcus (UID298) | 56 | 805 | 99.38 | 0.57 | 0.00 |
| RKP55 | g__Staphylococcus (UID298) | 56 | 805 | 99.33 | 0.57 | 0.00 |
| RKP52 | g__Staphylococcus (UID298) | 56 | 805 | 99.33 | 0.57 | 0.00 |
| RKP51 | g__Staphylococcus (UID298) | 56 | 805 | 99.33 | 0.57 | 0.00 |
| RKP32 | g__Staphylococcus (UID294) | 60 | 773 | 99.25 | 0.00 | 0.00 |
| RKP28 | f__Micrococcaceae (UID1623) | 39 | 457 | 98.70 | 0.23 | 0.00 |
| RKP23 | k__Bacteria (UID203) | 5449 | 104 | 98.28 | 37.93 | 85.48 |
| RKP44 | g__Staphylococcus (UID298) | 56 | 805 | 97.34 | 5.16 | 37.14 |
Explain: RKP4 g__Staphylococcus (UID298) 0.17 100.00
这是一个非常典型且有趣的现象。在 CheckM 的结果中,低污染(Contamination = 0.17%) 但 极高菌株异质性(Strain Heterogeneity = 100.00%) 通常意味着:你的这个 Bin 里实际上混合了多个高度相似的菌株(Strains),而不是多个不同的物种(Species)。
以下是详细的解读:
1. 为什么污染率这么低 (0.17%)?
- 定义: CheckM 的“污染”主要检测的是不同物种之间的标记基因重复。例如,如果 Bin 里同时有 Staphylococcus aureus 和 Escherichia coli,它们的核心基因完全不同,CheckM 会发现大量非预期的标记基因,从而报出高污染。
- 现状: RKP4 被 GTDB-Tk 鉴定为 Staphylococcus haemolyticus。Bin 里的所有序列都属于这个物种(或极近缘物种)。因为物种层面的核心基因集是统一的,所以没有触发“跨物种污染”的警报。
2. 为什么菌株异质性高达 100%?
- 定义: “菌株异质性”检测的是同一物种内不同菌株之间的序列差异。
- 机制: 当样本中存在两个或多个高度相似的 S. haemolyticus 菌株时,组装软件(Assembler)往往无法将它们完全分开,而是会将它们组装成一个个“嵌合”的 Contig,或者将属于不同菌株的相似区域合并在一起。
- CheckM 的信号: CheckM 在比对标记基因时,发现某些本该是“单拷贝”的基因位点上,存在显著的序列多态性(SNPs/Indels),或者同一个基因出现了两个高度相似但不完全相同的副本(Paralogs due to strain variation)。
- 100% 的含义: 这表示在该 Bin 的所有可评估标记基因中,几乎每一个都检测到了这种“多菌株混合”的信号。这说明该 Bin 不是由一个纯克隆(Pure Clone)组成的,而是一个种群(Population)的混合体。
3. 结合 GTDB-Tk 结果看 RKP4
- 分类: Staphylococcus haemolyticus
- ANI: 98.9% (vs GCF_006094395.1)
- AF (Alignment Fraction): 0.876
- 注意: AF 为 0.876 意味着只有约 87.6% 的参考基因组与你的 Bin 比对上了。对于高质量 MAG,AF 通常 >0.9 甚至 >0.95。
- 关联: 较低的 AF 加上 100% 的异质性,强烈暗示你的 Bin 可能是一个泛基因组(Pangenome)级别的混合,或者包含了该物种中某个参考基因组所没有的独特区域(可能是另一个菌株特有的基因岛)。
4. 对下游分析的影响与建议
| 分析类型 | 影响程度 | 建议 |
|---|---|---|
| 物种注释/功能概况 | 🟢 低 | 基本可用。因为都是 S. haemolyticus,物种层面的功能预测(如代谢通路)通常是准确的。 |
| 系统发育树 (Phylogeny) | 🟡 中 | 需谨慎。 如果建树使用的是单拷贝标记基因(Single-copy markers),混合菌株会导致分支长度异常或拓扑结构错误。建议在建树前使用工具(如 phyloFlash 或手动筛选)提取一致性序列,或者直接剔除该样本以保证树的稳健性。 |
| SNP calling / 进化分析 | 高 | 不可用。 100% 异质性意味着你无法区分哪些变异是测序错误、组装错误,还是真实的菌株间差异。强行做 SNP 分析会得到错误的进化距离。 |
| 毒力因子/耐药基因 | 🟡 中 | 需注意基因拷贝数。混合菌株可能导致某些耐药基因看起来有“多个拷贝”,但实际上可能分属不同菌株。 |
总结
RKP4 是一个“干净但混杂”的 Bin。 它没有混入其他物种(低污染),但它内部包含了多个 S. haemolyticus 菌株的序列(高异质性)。
- 如果你的目标是构建高精度的系统发育树: 听从你之前的建议,删除 RKP4 是最稳妥的选择,因为它会引入噪音。
- 如果你只想看群落组成或功能潜力: 可以保留,但需知晓它代表的是一个菌株群体而非单一基因组。