《Biology》：Dynamics of Antibiotic Resistance Gene Profiles in Captive Forest Musk Deer (Moschus berezovskii) Along a Breeding Duration Gradient
编辑推荐：
背景：为保护野生种群并确保麝香可持续供应，中国启动了林麝圈养繁殖。沿育种时长梯度的圈养林麝肠道抗生素抗性基因（ARG）谱的时间动态仍鲜有表征。方法：研究人员采用宏基因组测序系统表征ARG谱及潜在迁移性。样本按育种时长分为短期、中期和长期组。结果：所有样本中总计
背景：为保护野生种群并确保麝香可持续供应，中国启动了林麝圈养繁殖。沿育种时长梯度的圈养林麝肠道抗生素抗性基因（ARG）谱的时间动态仍鲜有表征。方法：研究人员采用宏基因组测序系统表征ARG谱及潜在迁移性。样本按育种时长分为短期、中期和长期组。结果：所有样本中总计注释331个ARG亚型和71个移动遗传元件（MGE）亚型。ARG Shannon多样性各组间总体差异显著（Kruskal-Wallis，p = 0.03）；Bonferroni校正Dunn检验未显示显著两两差异。Bray–Curtis PCoA显示短期组明显分离（p = 0.002），共享核心ARG亚型随育种年限延长逐渐增加。ARG与MGE丰度强正相关（r = 0.85，p = 0.0001）。共恢复63条携带ARG-MGE共定位复合体的contig，短期组此类contig比例最高，且在1–3 kb基因组区间内呈现紧密物理连锁。结论：育种时长与圈养林麝肠道ARG特征相关，短期圈养ARG-MGE共定位更高，提示迁移可能性更大。
研究背景与立项依据
林麝（Moschus berezovskii）为中国一级保护动物，其雄性麝香是传统中药与高档香料的重要原料。自1958年起，中国以圈养繁殖作为迁地保护核心策略，但集约化养殖使林麝易患消化道与呼吸道疾病，经验性抗菌治疗普遍，β-内酰胺类、氟喹诺酮类与氨基糖苷类常用。既往研究多聚焦野生与圈养静态比较，而圈养本身非均一条件，管理强度、种群密度、膳食结构与抗菌药暴露均随育种年限变化，形成持续选择压。沿育种时长梯度刻画肠道抗生素抗性基因（ARG，antibiotic resistance gene）与移动遗传元件（MGE，mobile genetic element）动态，对One Health框架下耐药污染管控具关键意义。该研究发表于《Biology》，通过宏基因组学揭示时间维度耐药组演化轨迹。
主要技术方法概览
研究人员按育种年限从15个独立鹿场选样，分短期（ST，<5年）、中期（MT，5–10年）、长期（LT，>10年）三组各5场；每场取10头健康成年林麝新鲜粪便等量混合成复合样本，液氮速冻后?80 ℃保存。提取总DNA后于MGI-T7平台行PE150宏基因组鸟枪测序；Fastp质控、BWA比对林麝参考基因组去除宿主序列得clean reads；MEGAHIT组装contig，METAProdigal预测ORF，CD-HIT聚类去冗；ARG以DIAMOND比对SARG数据库（E≤1×10?5，氨基酸一致性≥70%，覆盖≥40%），MGE以BLASTn比对mMGE数据库（E≤1×10?5，核苷酸一致性≥70%，覆盖≥40%）；Salmon算TPM丰度并log10转换。同contig内ARG-MGE视为共定位，统计物理间距。α多样性用Shannon指数，β多样性用Bray–Curtis PCoA，组间比较用Kruskal–Wallis＋Bonferroni校正Dunn检验，相关性用Spearman秩相关。
研究结果
3.1 ARG谱表征：15样本共获23型331亚型，ST组总丰度最高、LT最低但组间无显著差异（p=0.10）；四环素类、MLS（大环内酯-林可酰胺-链阳菌素）与多重耐药基因为主。Shannon多样性总体差异显著（p=0.03）但两两不显著，ST偏高；PCoA显示ST组与另两组分离（p=0.002）。
3.2 共享特征与高风险ARG：共享亚型占比从ST的67.1%升至MT 83.0%、LT 88.6%；ST独有56亚型，MT、LT仅7和6。lnu(C)、rpoB、CfxA6、APH(3′)-IIIa、ugd丰度跨组差异显著。风险分级检得31个Rank I与16个Rank II亚型，erm(B)/erm(F)/APH(3′)-IIIa为优势高危ARG；MT高危丰度最高但组间无统计差异（p=0.50）。
3.3 MGE丰度与组成：71亚型分转座酶、整合酶、质粒、IS元件、ist五类。MGE总丰度组间差异显著（p=0.04），组成差异显著（p=0.008）；转座酶与IS元件占主导，ST组IS比例偏高、整合酶偏低，tnpA与IS26丰度变化显著。
3.3（原文3.4）ARG-MGE相关性：总丰度Spearman相关r=0.85（p=0.0001）。亚型层面，ST中ANT(2″)-Ia与repUS16完全相关，APH(3′)-Ia、tet(X3)与Int-Tn916/IS91/tnpA强相关；MT中sul1-repUS16/tniB、tet(U)-rep14强相关；LT中ANT(2″)-Ia、erm(C)与repUS16强相关。
3.5 ARG-MGE共定位：63条共定位contig涵盖49 ARG与31 MGE亚型，含16高危ARG（如aadA5、CARB-4、tet(A/M)等）。ST组共定位contig比例最高；累积共现频率随物理距离上升，ST在1–3 kb区间积累最快、终频最高（6.21%@3732 bp），LT最低（4.55%）。APH(3′)-Ia-tnpA、rep14-tet(U)、repUS16-erm(C)等呈<1 kb紧密连锁。
讨论浓缩
圈养林麝肠道为重要ARG储库，优势类群与集约养殖系统一致。ARG丰度随育种时长呈非显著下降，可能与短期多源输入（土壤等环境储库）及后期核心亚型收敛有关；MT/LT组α多样性下降、共享核心上升，符合持续药压下耐药组收敛模型。ARG-MGE强相关与共定位contig提供潜在水平基因转移（HGT）基因组线索，但池化样本无法区分个体携带率，亦非活性转移直接证据。ST组虽丰度未最高，但共定位比例与中短距累积最强，提示短期圈养迁移潜力最大；LT组耐药组稳定但高危ARG仍与tnpA/IS91等连锁。高危ARG（APH(3″)-Ib、APH(6)-Id、mph(A)、tet(A)等）可经Tn3/IS26/Inc型质粒在人-畜-环境间共享，鹿场或成人与动物界面交换节点，但本研究未证实实际交换。
结论翻译
综上，育种时长与林麝肠道ARG特征相关：长期圈养关联更稳定、收敛的耐药组；短期圈养呈现更高多样性与更强ARG-MGE共定位，提示更高迁移潜力。鉴于观察性设计且缺抗菌药暴露记录，上述关联应谨慎解读。各时长组高危ARG与MGE的物理连锁凸显定期监测紧迫性。未来需整合用药数据、管理调查与接合等功能实验以确立因果并验证实际HGT频率。
