《mSystems》：Gain and loss of plasmid-borne antibiotic resistance genes are associated with chromosomal resistance presence in Enterobacteriaceae
编辑推荐：
摘要：质粒是抗生素抗性基因（antibiotic resistance genes, ARGs）传播的核心载体，其作为最具流动性和可进化性的遗传元件，具有广泛的宿主范围和高速率的基因更替，能极有效地在不同细菌谱系间扩散抗生素抗性。研究人员应用基于系统发育感知(
摘要：质粒是抗生素抗性基因（antibiotic resistance genes, ARGs）传播的核心载体，其作为最具流动性和可进化性的遗传元件，具有广泛的宿主范围和高速率的基因更替，能极有效地在不同细菌谱系间扩散抗生素抗性。研究人员应用基于系统发育感知(phylogeny-aware)的基因获得(gain)与丢失(loss)模型，对6895株肠杆菌科(Enterobacteriaceae)基因组中的质粒携带基因进行了量化，分析了基因家族获得、丢失、扩增(expansion)和缩减(reduction)四类进化过程。研究发现，总体而言，质粒携带ARG(pARG)与非抗性质粒基因的获得率相似，但其扩增率和缩减率显著更高。上述四种过程均呈强烈的种属依赖性(species-dependent)，抗生素类别的影响较小。此外，携带特定染色体ARG(chromosomal ARG, cARG)的细菌分支(clade)相较于不携带对应cARG的姐妹分支，表现出显著更高的质粒抗性基因获得率和更低的丢失率。研究人员还发现IncQ2质粒骨架与qnrS2相关联且专属分布于勒克氏菌(Leclercia adecarboxylata)，而携带mprF的Col(VCM04)质粒主要(71.4%)分布于柠檬酸杆菌属(Citrobacter)。综上，质粒介导的抗性主要受种属驱动，cARG可有效标记具有高质粒抗性基因获得能力的细菌谱系。
论文解读：《肠杆菌科中质粒携带抗生素抗性基因的获得与丢失与染色体抗性基因存在与否相关联》
该研究发表于《mSystems》。抗生素耐药性(antimicrobial resistance, AMR)是医学与进化生物学领域的重大挑战，细菌基因组进化中通过水平基因转移(horizontal gene transfer, HGT)获得的基因获得及选择作用下的基因丢失是塑造基因组 repertoire 的主要力量，其中可接合质粒(conjugative plasmids)等可移动遗传元件(mobile genetic elements, MGEs)是HGT的核心载体，对抗生素抗性基因(antibiotic resistance genes, ARGs)的传播至关重要。已有比较基因组学显示携带ARG的质粒分类单元具有更高流动性、更广宿主范围和更快的基因更替率，其可塑性由转座元件和重组系统增强，复制调控突变还可提高质粒拷贝数从而放大抗性水平，因此质粒是具频繁水平转移、快速更替和广泛遗传可塑性的动态复制子。肠杆菌科(Enterobacteriaceae)是临床重要耐药菌的主要储存库，包含大肠埃希菌(Escherichia)、克雷伯菌(Klebsiella)、肠杆菌(Enterobacter)和沙门菌(Salmonella)等，虽已有大量临床研究记录特定耐药质粒传播或单个ARG流行病学，但质粒携带ARG(plasmid-borne ARGs, pARGs)如何跨谱系维持与扩增、其进化动态是否区别于其他质粒功能基因尚不清楚。本研究通过对近7000株肠杆菌科基因组应用系统发育感知的出生–死亡(birth-death)模型，系统探究质粒介导ARG的更替(turnover)动态，量化各谱系的基因家族获得(gain)、丢失(loss)、扩增(expansion)和缩减(reduction)平衡，并评估染色体ARG(chromosomal ARGs, cARGs)与质粒抗性组(repertoire)的关联。
主要关键技术方法：
研究人员从GTDB r220筛选肠杆菌科高质量完整或近完整基因组（满足MIMAG标准：完整度>99%、污染<1%、平均contig长>5000 bp、contig数<500），每物种最多取300株并二次抽样，最终纳入6895个基因组涵盖159个种，其中5335个含鉴定出的质粒序列。用geNomad v1.11.0与Platon v1.7交叉确认质粒及其编码基因；用Resistance Gene Identifier (RGI) v6.0.5比对CARD数据库v4.0.1注释ARG；OrthoFinder v2.5.5构建各物种泛基因组并进行直系同源基因家族聚类，定义质粒关联基因家族为家族内≥70%成员为质粒来源基因；用Count v10.04基于物种树对各物种质粒直系同源基因家族拟合系统发育出生–死亡模型估计基因获得(κ)、丢失(μ)及基因重复/扩增(λ)参数，经最大似然迭代优化后用Posteriors模块计算各分支获得、丢失、扩增、缩减事件的后验概率，按分支总长标准化得进化速率；种间比较采用配对Wilcoxon检验，方差分解用普通最小二乘(OLS)模型；姐妹分支(sister clade)分析筛选含特定cARG家族≥80%基因组的分支与<20%的姐妹分支配对，用线性混合效应模型(LMM)评估cARG存在对pARG进化速率影响；用MOB-suite v3.1.9重构质粒单位，Mash距离建质粒相似性网络并做复制子(replicon)分型。
研究结果：
质粒携带抗生素抗性基因与非抗性质粒基因相比具有更高扩增与缩减率但获得动态相似
研究人员应用系统发育出生–死亡模型量化5335个含质粒基因组中质粒基因家族四类过程。比对80个物种中pARG与其他质粒基因(non-ARG)整体速率发现：pARG与non-ARG的获得率(gain rate)无显著差异（校正P=0.257），表明抗性基因家族整体新获过程遵循普通质粒进化途径；而pARG丢失率(loss rate)显著升高（校正P=1.23×10?2），扩增率(expansion rate)（校正P=3.34×10??）与缩减率(reduction rate)（校正P=4.03×10??）均显著更高，其中扩增中位log?比为1.28、缩减为0.83。表明pARG在基因家族水平获得无特异性信号，但其拷贝数动态（扩增与缩减）受种属水平驱动明显强于其他质粒基因。
抗生素抗性基因更替动态主要依赖菌种（种属依赖性）
将直系同源基因家族按CARD定义药物类别注释后分析发现，四类事件速率均强烈依赖于菌种（Kruskal–Wallis检验gain H=536.9 P=2.9×10??? η2=0.76；loss H=506.3 P=1.4×10??3 η2=0.71；expansion H=312.5 P=1.6×10?2? η2=0.39；reduction H=436.9 P=2.4×10??1 η2=0.60），OLS方差分解显示菌种身份解释获得、丢失、扩增、缩减总方差分别达66.8%、63.1%、42.1%、51.5%；而药物类别效应弱且不显著（gain P=0.74, loss P=0.58）或仅为小效应（expansion P=5.2×10?? η2=0.09；reduction P=1.1×10?3 η2=0.04），解释方差均<7%。此外高扩增率菌种对同一药物类别通常也呈高缩减率（Spearman ρ=0.792, 置换P≈0.0002），而高获得率菌种未必对应高丢失率（ρ=?0.196, P≈0.079），说明扩增与缩减耦合紧密，获得与丢失相对独立。
染色体抗生素抗性基因与pARG获得率升高相关联
经姐妹分支分析筛选464对分支（含80个种、84个RGI定义ARG家族），其中一支≥80%基因组含特定cARG家族而姐妹支<20%。LMM结果显示：含cARG分支的pARG获得率约为姐妹支2.47倍（估计值0.91对数尺度，P=6.75×10??），pARG丢失率显著降低（fold-change=0.138, P=6.10×10?2?）；而其他质粒基因在cARG分支获得率反而低于姐妹支（约3.7倍低，P=2.97×10??），丢失率也降低（P=1.44×10??）。交互项显示cARG存在对不同基因类型获得与丢失影响差异显著（gain交互P=2.2×10??, loss交互P=5.6×10??）。阈值敏感性测试中70/30至90/10各梯度均保持正向关联，cARG定义越严格pARG获得偏向越大（70/30阈值fold-change=1.65升至90/10阈值=4.91），pARG丢失抑制在各阈值均稳定显著（fold-change 0.183~0.125, 所有P<0.001）。cARG分支 median保有2.64个独特质粒簇 vs 姐妹支1.66（P=2.42×10?21），且pARG关联复制子类型在两姐妹支间median共享比率为0%，提示这些质粒多为分支分化后独立获得而非共祖继承，结合更低质粒基因丢失率说明cARG谱系更擅长获取并维持环境可动组(mobilome)。部分ARG家族如sul与dfr也显示本家族更高获得与更低丢失。
质粒携带抗生素抗性基因与质粒复制子类型的共现
重构质粒并做复制子分型及Mash距离网络分析发现特定质粒类型与ARG强关联：ColRNAI_rep_cluster_1987主要携带消毒抗性基因qacG；Col(VCM04)主要携带肽类抗性基因mprF（56个中40个即71.4%源自柠檬酸杆菌属Citrobacter）；IncQ2主要关联氟喹诺酮抗性基因qnrS2且仅发现于Leclercia adecarboxylata。广谱多药耐药相关质粒（如IncC、IncHI2A、IncN）则呈复杂异质ARG组成，少有明确单基因–复制子专一关联。
讨论与结论总结：
研究人员指出pARG与其他质粒基因在种属水平获得率相当，但扩增与缩减率显著更高且具强种属依赖，可能反映IS26转位单元、整合子盒重排及选择性扩增等基因水平可动元件活动，复制调控突变也可提升质粒拷贝数。cARG谱系展现更高pARG获得与更低丢失，此关联可由两种互不包含的进化情景解释——染色体整合为先导背景促进后续质粒摄取，或高强度质粒流下高频转座/整合致cARG作为高动态抗性可动组后果；无论方向cARG可作谱系倾向积累质粒介导ARG的基因组标记。cARG分支更高质粒负荷与独立后获性质粒、低丢失率共同支持其高效获取–维持环境可动组，多质粒共存利于不同pARG物理接近促重组与新耐药岛形成，再经低丢失率稳定下来推动多重耐药传播。少数Count推断高后验事件匹配已知暴发（如Enterobacter quasihormaechei中blaCTX-M获得对应法国ST873暴发、Escherichia fergusonii中mcr家族获得对应中国mcr-1流行等）验证模型可靠性。研究局限含未纳入ICE/IME及噬菌体、缺生态抗生素暴露史。
结论（翻译）：
研究发现肠杆菌科中质粒携带ARG(pARG)在整体基因家族水平遵循与其他质粒基因相同的更替动态，表明抗性扩散主要反映质粒水平过程而非ARG特有规则；在拷贝数变化层面pARG显示受谱系驱动且药物类别影响微弱的升高扩增与缩减，暗示反复扩增–缩减循环是一种灵活耐药策略；此外染色体抗性内容可作为进化时间尺度上细菌谱系具质粒介导ARG积累潜能的基因组标记。
