共生菌群与寄生性线虫捻转血矛线虫伊维菌素耐药性的关联

《Animal Diseases》：Symbiotic bacteria are associated with ivermectin resistance in the parasitic nematode Haemonchus contortus

　　

抗蠕虫药物（anthelmintic）耐药性在胃肠道线虫中的广泛传播对全球畜牧业生产及可持续性构成严重威胁。尽管寄生虫基因组层面的适应性改变已被广泛研究，但寄生虫相关微生物组（microbiome）在调控药物敏感性方面的作用仍有待探索。本研究旨在调查血液寄生性线虫捻转血矛线虫（Haemonchus contortus）的伊维菌素（ivermectin, IVM）耐药株与敏感株的共生微生物群，并评估其与伊维菌素耐药性的潜在关联。比较性微生物组分析揭示，嗜麦芽窄食单胞菌（Stenotrophomonas maltophilia）可作为潜在的耐药性生物标志物，其在耐药株中的相对丰度显著升高3.4倍。此外，功能预测分析表明，遗传信息处理及细胞运动相关通路在耐药株微生物群中显著富集。幼虫致死率实验证实，庆大霉素（gentamicin）处理显著增强幼虫对伊维菌素的敏感性。上述发现提示，共生菌群与捻转血矛线虫的伊维菌素耐药性存在相关性，且该作用可能与寄生虫遗传背景因素协同发挥效应。

## 研究背景与问题提出

伊维菌素（ivermectin, IVM）作为大环内酯类（macrocyclic lactone）抗寄生虫药物，自1981年问世以来一直是 livestock 及人类寄生虫病防控的基石性药物，其广谱高效特性彻底革新了线虫感染的管理策略。然而，该药的密集及不规范使用施加了巨大的选择压力，导致抗蠕虫药物耐药性（anthelmintic resistance, AR）在全球范围内广泛涌现。捻转血矛线虫（Haemonchus contortus）作为小反刍动物最具致病性的胃肠道线虫，同时也是研究抗蠕虫药物耐药性的核心模式生物，其伊维菌素耐药性已在全球24个国家被报道，严重威胁绵羊及山羊产业的可持续发展。

数十年来，研究工作主要集中于解析寄生虫基因组内部的耐药性分子机制，包括谷氨酸门控氯离子通道（glutamate-gated chloride channel）靶点突变、P-糖蛋白（P-glycoprotein）等细胞外排泵上调，以及细胞色素P450（cytochrome P450, CYP）和谷胱甘肽-S-转移酶（glutathione S-transferase）等解毒酶类的改变。尽管上述基因组机制的研究取得显著进展，但现有机制尚无法完全解释田间种群中观察到的耐药性表型多样性，提示可能存在其他协同作用因素。随着药物微生物组学（pharmacomicrobiomics）的兴起，共生菌群在宿主药物代谢中的作用日益受到关注。研究表明共生菌可直接参与杀虫剂的代谢降解，在秀丽隐杆线虫（Caenorhabditis elegans）等模式生物中，共生菌参与营养代谢、寿命调控及免疫防御等多种宿主生理过程。捻转血矛线虫本身亦携带稳定的细菌群落，且抗菌处理显著影响幼虫存活，暗示共生菌在寄生虫生理中具有重要功能，但其是否参与伊维菌素耐药性仍缺乏证据。基于此，研究人员提出科学假设：捻转血矛线虫的共生菌群与其伊维菌素耐药性存在关联。

## 关键技术方法

本研究采用体内山羊感染模型获取实验材料，涉及两株伊维菌素敏感株（AS、ZJ）和三株耐药株（WS、CY、XJR），每株设5个混合虫体样本。微生物群落分析采用 PacBio 单分子实时测序（single-molecule real-time, SMRT）技术对16S rRNA基因全长（V1-V9区）进行测序；功能预测采用 Tax4Fun 软件基于 Kyoto Encyclopedia of Genes and Genomes（KEGG）数据库进行。差异物种鉴定采用线性判别分析效应量（linear discriminant analysis effect size, LEfSe）方法。功能验证实验选取耐药株XJR的第四期幼虫（L4），经亚致死浓度庆大霉素处理降解共生菌后，进行伊维菌素敏感性测定，并通过定量PCR（qPCR）检测嗜麦芽窄食单胞菌丰度变化。

## 研究结果

### 样本采集与16S rRNA基因测序统计

研究人员利用体内山羊感染模型采集伊维菌素敏感株（S组）和耐药株（R组）样本，共纳入5个样本组。经环化共有序列（circular consensus sequencing, CCS）读长生成、引物修剪及质量过滤后，成功保留334,128条序列，每样本平均13,365条。物种累积曲线趋于平缓，表明测序深度足以捕获微生物多样性。

### 捻转血矛线虫共生菌群的微生物多样性

全长16S rRNA基因测序显示，敏感组含2,109个操作分类单元（operational taxonomic unit, OTU），耐药组含1,756个OTU。α多样性指数（观测物种数、Shannon指数及Good's覆盖度）在两组间无显著差异，表明伊维菌素耐药性不改变整体微生物群落丰富度或均匀性。基于Bray-Curtis相异度矩阵的主坐标分析（principal coordinate analysis, PCoA）显示两组部分重叠，但相似性分析（analysis of similarities, ANOSIM）证实组间群落结构存在中等程度但统计学显著的差异（R=0.35, P=0.001）。门水平聚类分析显示，地理来源相近的耐药株WS和CY聚类最为紧密，提示地理因素可能影响微生物群落结构。

### 伊维菌素敏感株与耐药株的微生物组成差异

在门水平，两组共享变形菌门（Proteobacteria）、厚壁菌门（Firmicutes）和放线菌门（Actinobacteria）等优势类群。科水平分析显示，敏感组富集伯克氏菌科（Burkholderiaceae）、假单胞菌科（Pseudomonadaceae）和瘤胃球菌科（Ruminococcaceae），而耐药组以伯克氏菌科、黄单胞菌科（Xanthomonadaceae）和瘤胃球菌科为主。统计学分析显示黄单胞菌科在耐药组中显著富集。属水平上，嗜麦芽窄食单胞菌属（Stenotrophomonas）在耐药组显著富集；种水平上，嗜麦芽窄食单胞菌（S. maltophilia）、皮氏罗尔斯顿氏菌（Ralstonia pickettii）和约氏乳酸杆菌（Lactobacillus johnsonii）在耐药组相对丰度显著更高，其中S. maltophilia丰度升高3.4倍。

### 耐药性相关共生菌的鉴定

LEfSe分析鉴定出稳健的分辨生物标志物（线性判别分析得分>4.0）。敏感组显著富集拟杆菌门（Bacteroidetes），而耐药组特征性富集黄单胞菌科。属水平上，嗜麦芽窄食单胞菌属和罗尔斯顿氏菌属在耐药组显著富集；种水平上，S. maltophilia被确定为生物标志物。

### 伊维菌素敏感株与耐药株共生菌的功能预测

采用Tax4Fun基于16S rRNA基因测序数据进行功能预测。KEGG（Kyoto Encyclopedia of Genes and Genomes）一级水平分析显示两组主要功能组成相似，以代谢、遗传信息处理和环境信息处理为主。二级水平分析表明，耐药组"遗传信息处理"和"细胞运动"相关通路相对丰度显著增加，而"聚糖生物合成与代谢"通路显著下调。三级水平分析进一步揭示了具体代谢通路的组成差异。针对"细胞色素P450介导的外源物质代谢"通路的评估显示，尽管耐药组均值略高，但差异未达统计学显著水平。

### 抗菌处理显著增强耐药幼虫对伊维菌素的敏感性

为功能验证共生菌对伊维菌素耐药性的贡献，研究人员对耐药株XJR的L4幼虫进行庆大霉素处理以降解共生菌。毒性预实验确定300 μg/mL为最佳亚致死浓度，该剂量对幼虫活力无显著影响。qPCR定量证实庆大霉素处理有效降低幼虫体内S. maltophilia的相对丰度。低浓度伊维菌素（50 nM）联合处理仅导致死亡率非显著性增加，但在判别剂量100 nM伊维菌素条件下，庆大霉素预处理组幼虫死亡率较单独伊维菌素组显著升高（P=0.0070），表明抗菌处理显著增强耐药幼虫对伊维菌素的敏感性。

## 讨论总结

研究人员就共生菌群变异与伊维菌素耐药性的关联展开深入讨论。传统上抗蠕虫药物耐药性主要归因于寄生虫基因组变异，但这些经典机制无法完全解释田间耐药表型的多样性。本研究通过全长16S rRNA基因测序发现伊维菌素耐药株与敏感株间共生菌群存在显著差异，且抗菌介导的菌群清除显著增强幼虫对伊维菌素的敏感性，提示共生菌群与耐药性表型相关。研究人员提出，菌株间的遗传差异构成耐药性基线，而共生菌群的改变与之相关，可能协同促成表型复杂性，从而形成更为综合的耐药模型。

尽管α多样性无显著差异且PCoA存在重叠，ANOSIM分析证实了组间微生物群落结构的统计学显著差异，表明特定分类群的定向驱动可能是造成组间相异度的关键。功能预测分析提示遗传信息处理和细胞运动等广泛通路存在功能分化，但细胞色素P450介导的外源物质代谢等特异性降解通路未显著富集，这可能源于微生物功能冗余及16S rRNA预测工具在精细功能基因解析上的固有限制。

关于地理因素的影响，研究人员指出不同地理隔离菌株的共生菌组成差异显著，但地理因素并不能完全掩盖耐药性表型本身的作用效应。值得注意的是，尽管耐药株XJR与敏感株AS在实验室内共同保存数年，二者在聚类分析中距离较远，这降低了近期环境趋同塑造群落结构的可能性，进一步凸显耐药性表型与共生细菌组成之间的关联性。

在技术方法层面，本研究采用PacBio SMRT测序靶向16S rRNA基因V1-V9高变区，相较于传统仅限于V3-V4区的二代测序，提供了种水平的增强分类学分辨率。然而，该技术对低丰度物种鉴定可能存在影响。

关于共生菌获得与发育阶段动态，研究人员讨论了线虫从自由生活阶段到寄生阶段过程中共生菌的潜在获得途径，包括水平摄取和垂直传播等方式。但由于本研究微生物组分析主要局限于成虫样本，尚无法完全捕获不同发育阶段共生菌的组成和功能转变。

针对S. maltophilia的潜在耐药贡献，研究人员引用前期研究指出该菌可降解阿维菌素（伊维菌素的结构类似物），推测其可能直接参与伊维菌素降解；或通过激活寄生虫内在解毒系统、增强一般胁迫耐受性等间接途径发挥作用。尽管其他差异富集物种亦被观察到，但S. maltophilia在种水平上相对丰度最高，且LEfSe分析在所有耐药株中呈现一致的富集趋势，故被优先选择进行深入研究。

在研究局限性方面，研究人员坦承当前证据尚属相关性层面，缺乏严谨的"回补"实验以确立因果关系；地理混杂因素尚无法完全排除；功能探索主要基于16S rRNA基因序列的生物信息学推断，有待转录组学、代谢组学或体外降解实验验证；样本量限制及基于单一判别剂量的二元耐药分类阻碍构建细菌丰度与耐药水平的定量回归模型；该共生菌介导的耐药机制在其他地理来源野生型菌株或其他胃肠道线虫中的普遍存在性亦有待系统调查。

## 研究结论

本研究提供了支持共生菌群与捻转血矛线虫伊维菌素耐药性存在潜在关联的证据。研究人员鉴定出嗜麦芽窄食单胞菌（Stenotrophomonas maltophilia）作为关键的耐药性相关生物标志物。抗菌介导的细菌清除显著增强伊维菌素敏感性。上述发现凸显寄生虫相关共生菌群作为开发整合性抗蠕虫药物耐药性管理策略的潜在靶标价值。