耐药基因检测与结果解读
耐药基因检测是在已测序的细菌分离株或宏基因组数据中，识别已知的抗微生物耐药基因（ARG）和耐药相关点突变，说明样本携带哪些耐药决定因子、可能涉及哪些药物类别。输入通常是分离株的基因组组装、分离株原始读段，或样本的宏基因组读段；分析靠专门的参考数据库和配套工具完成，常用的有 CARD 配 RGI、AMRFinderPlus 和 ResFinder。得到的基因清单只是基因型，要与实验室药敏结果比对，才能知道它对表型的解释力有多强。结果多用于临床微生物学研究、医院感染控制和耐药监测。本页介绍的是科研与监测分析方法，不替代临床药敏报告与治疗决策。
什么时候用
- 描述一批临床分离株的耐药基因谱，比较不同科室或时间段
- 医院感染暴发调查中核对可疑菌株的耐药决定因子是否一致
- 耐药监测项目中统计碳青霉烯酶等重点基因的检出率与趋势
- 药敏显示耐药但机制不明时，从基因组中找可解释的基因或突变
- 在宏基因组样本中评估耐药基因的总负荷与组成
- 判断重点基因位于染色体还是质粒上，用于感染控制评估
不适合的情况
- 直接出临床报告：治疗决策仍以标准药敏试验结果为准
- 发现全新机制：基于已知数据库只能找到已收录的基因
- 只有 16S 扩增子数据：扩增子里没有耐药基因的序列信息
需要准备的数据
- 基因组组装或 contig
- 分离株建议先组装再检索，需记录组装质量（contig 数、N50）和测序平台。
- 原始读段 FASTQ
- 双端读段，已去接头和低质量碱基，用于读段级检索与点突变判读。
- 宏基因组读段
- 需记录测序深度和宿主序列去除情况；同批样本深度尽量接近，便于归一化。
- 参考数据库
- CARD、NCBI Reference Gene Catalog 或 ResFinder 数据库，记录版本与下载日期。
- 菌株与药敏信息
- 物种鉴定结果、MIC 或纸片法结果及判读标准（CLSI 或 EUCAST）、采样时间和来源。
分析步骤
- 1
明确问题与路线
先确定要回答的是单株菌携带哪些耐药基因，还是一批样本的耐药基因负荷。分离株通常走组装后检索，宏基因组和需要检出低丰度基因的场景走读段级检索。两条路线的阈值设定和定量口径不同，中途换路线会让结果无法合并。
- 2
质控与组装
先对读段做质控和修剪，分离株再做从头组装，并确认物种鉴定结果与是否存在污染。组装过于碎片化、contig 数异常多时，基因容易被打断，后面会出现大量只检出一部分序列的命中。
- 3
选择数据库与工具
按目的选库：CARD 的 ARO 本体便于按机制归类，NCBI Reference Gene Catalog 配 AMRFinderPlus 命名统一并支持物种特异位点，ResFinder 侧重获得性基因。记录库与工具版本，监测项目应全程固定版本。
- 4
比对与阈值过滤
运行检索后按一致性和覆盖度过滤命中。RGI 的 perfect 与 strict 命中可直接采纳，loose 命中需逐条人工核对。覆盖度明显偏低或位于 contig 末端的命中标注为不完整，不要当成完整基因上报。
- 5
点突变判读
喹诺酮、利福平等药物的耐药常由靶基因的特定密码子替换决定，必须比对到具体位置，而不是看基因有无。AMRFinderPlus 的 --organism 参数和 ResFinder 配套的 PointFinder 提供物种特异的位点清单；没有对应物种时只报告获得性基因。
- 6
基因组背景与定位
查看重点基因所在 contig 的上下文，包括相邻的插入序列、整合子结构和质粒复制子（可用 PlasmidFinder 判断），据此推测它位于质粒还是染色体。短读长组装常把质粒打碎，必要时补长读长做混合组装。
- 7
与药敏结果比对
按药物逐一统计基因型预测与药敏表型的一致性，并看不一致的方向。表型耐药而基因型阴性，提示机制未被收录，或由外排泵、孔蛋白缺失、表达水平变化导致；基因型阳性而表型敏感，提示基因未表达或只带来低水平耐药。
- 8
汇总与报告
输出每株或每样本的基因清单、药物类别汇总和不一致病例列表，写明数据库版本、工具版本、阈值和药敏判读标准。监测用途还要给出重点基因的检出率和时间趋势，并说明纳入菌株的来源构成。
关键参数与决策
一致性与覆盖度阈值
获得性耐药基因常从 90% 一致性、60% 覆盖度起步（经验值）。阈值放低会把同家族中不介导耐药的同源基因算进来，放得过高则漏掉新等位基因。重点基因建议再看比对细节和功能域是否完整。
RGI 命中标准
perfect 指与参考序列完全一致，strict 指落在模型的分值阈值内，loose 指低于阈值的远同源命中。报告通常只采纳 perfect 和 strict，loose 仅作探索用途，必须人工确认后才能提及。
物种参数
点突变判读依赖物种特异的位点清单，物种选错会直接漏判或误判。先确认鉴定结果再选 --organism 或 PointFinder 的物种；工具不支持该物种时，只报告获得性基因并在方法中说明。
核酸输入与蛋白输入
蛋白序列检索对远同源更敏感，但依赖基因预测是否正确；核酸检索能发现假基因和截断，却更容易受组装错误影响。两种输入结合使用可以减少漏检，差异处逐条核对。
宏基因组归一化
常用每百万读段的比对数（RPKM）或每个基因组当量的拷贝数。按基因组当量归一化更接近「携带该基因的细菌占多少比例」，适合跨样本比较；只用原始计数会被测序深度主导。
数据库版本
基因命名和收录范围会随版本变化，同一基因在不同版本里可能换名或拆分。同一项目的所有样本用同一版本重跑，报告中写明版本和日期，跨年度趋势分析尤其要注意。
常见错误
- 把检出基因直接当成耐药：基因可能不表达、被插入序列打断，或只带来低水平耐药，必须与药敏结果核对。
- 采纳 loose 或低一致性命中就报阳性：会把同家族中不介导耐药的同源基因计入，假阳性明显增多。
- 只看 gyrA、parC 有无来判断喹诺酮耐药：这些基因普遍存在，只有特定密码子替换才有解释意义。
- 组装不完整却按完整基因上报：contig 断裂会让基因只剩一段，覆盖度字段必须和基因名一起看。
- 宏基因组直接比较原始检出计数：测序深度和菌量不同，需按基因组当量或 16S 拷贝归一化后再比。
- 用未校正的长读长读段判点突变：碱基错误和移码会造成假的氨基酸替换，应先得到校正后的共识序列。
结果怎么看
耐药基因检出表
每行一个命中，含基因名、一致性、覆盖度、所在 contig 与位置、对应药物类别。优先看覆盖度接近完整的命中；若大量命中覆盖度只有一半左右，先回头检查组装质量而不是直接解读。
耐药基因谱热图
行为菌株、列为基因或药物类别，用于比较一批分离株。同一科室、同一时间段的菌株呈现几乎相同的基因谱时，提示可能存在同一克隆的播散，需结合菌株追踪的 SNP 距离确认。
基因型与表型一致性表
按药物列出基因型预测与药敏结果的四格分布和一致率。某个药物一致率明显偏低，通常说明该药的耐药主要由外排、孔蛋白或表达调控决定，基因清单对它的解释力有限，报告中要写明。
基因组背景图
展示重点基因所在 contig 的邻近结构，如插入序列和整合子。基因落在带质粒复制子的 contig 上，提示它处于可移动元件背景中，在感染控制评估里权重更高，也更值得持续监测。
宏基因组 ARG 丰度表
归一化后的各基因或药物类别丰度。若总丰度与测序深度强相关，说明归一化没做好；若丰度集中在少数常见共生菌携带的基因，解读时要说明来源，不要直接当成临床耐药风险。
和相近方法怎么选
CARD/RGI、AMRFinderPlus 与 ResFinder
CARD 带 ARO 本体和机制分类，便于按作用机理归类，并区分同源模型与变异模型；AMRFinderPlus 基于 NCBI 的参考基因目录，命名统一、带 HMM，并支持物种特异的点突变判读；ResFinder 侧重获得性基因，配 PointFinder 处理染色体突变。实践中常跑两套交叉比对，差异处人工核对。
基于组装与基于读段
组装后检索能给出基因所在的序列背景，便于判断完整性和是否位于可移动元件上；读段级检索不依赖组装，对低丰度基因更灵敏，适合宏基因组，但拿不到上下文，且高度同源的基因之间容易误分配读段。两种结果的数值口径不同，不要直接混用。
基因型检测与表型药敏试验
药敏试验直接测量抑制浓度，是判读耐药的标准方法；基因型检测给出机制层面的解释和可用于流行病学追踪的标记。对部分药物与物种，两者一致性较好；对以外排和表达调控为主的组合则较差。二者互补，不能互相替代。
耐药基因检测与菌株追踪
耐药基因检测回答样本携带哪些耐药决定因子，菌株追踪回答两个样本的细菌是否同源。暴发调查通常两者同做：基因谱相近但核心基因组 SNP 距离很远，更像是相同耐药元件分别存在于不同菌株，而不是同一克隆播散。
常见问题
检出耐药基因就一定耐药吗？
不一定。基因可能不表达、拷贝数低、被插入序列打断，或只带来低水平耐药而未超过判读临界值。反过来，药敏耐药也可能找不到对应基因，因为外排泵过表达、孔蛋白缺失和调控区变化不一定被数据库收录。基因型结果应与药敏结果一起解读。
CARD 和 AMRFinderPlus 该选哪个？
看目的。需要按机制和药物本体归类、或分析远同源与点突变模型时，CARD 配 RGI 更方便；需要与公共监测数据命名一致、并做物种特异的点突变判读时，AMRFinderPlus 更合适。两者收录范围和命名不完全相同，重点基因建议都跑一遍并核对差异。
一致性和覆盖度阈值设多少合适？
获得性基因常从 90% 一致性、60% 覆盖度起步（经验值），并把低覆盖度命中单独标为不完整。想找新等位基因可把一致性降到 80% 左右，但要逐条核对比对区域和功能域是否完整。阈值必须写进方法，不同阈值下的结果不可直接比较。
宏基因组里怎么知道耐药基因属于哪个细菌？
短读长数据很难可靠回答。较稳妥的做法是先组装分箱得到 MAG，看基因是否落在某个 MAG 内；或用长读长，把耐药基因与同一条读段上的物种标记基因联系起来。仅凭样本里同时存在某个基因和某个物种，不足以把两者对应起来。
基因在质粒上还是染色体上，需要区分吗？
在监测和感染控制里有意义：位于质粒等可移动元件背景中的基因，更容易在不同菌株里反复出现，监测时需要关注同类元件的分布范围。判断依赖序列背景，短读长组装常把质粒打碎，通常要靠长读长或混合组装才能较好闭合。
在 Codar 里怎么做
在 Codar 中提供分离株组装或宏基因组读段、物种鉴定结果和已有药敏信息后，会先确认检索路线、参考数据库以及一致性与覆盖度阈值，再进行抗微生物耐药基因检测，给出基因检出表、药物类别汇总、耐药基因谱热图和重点基因所在 contig 的背景信息；提供药敏结果时同时统计基因型与表型的一致性。最后交付结果表、图、分析报告和可复现脚本。
我有 40 株肺炎克雷伯菌分离株的组装基因组和纸片法药敏结果，想查耐药基因、看碳青霉烯酶分布，并统计基因型和表型是否一致。