基于宏基因组学解析巴基斯坦开伯尔-普赫图赫瓦省奶牛粪便微生物组、耐药组与质粒组及其“同一个健康”意义

《Veterinary Sciences》：Shotgun Metagenomics Reveals Microbial Diversity, Resistome, and Plasmidome in Dairy Cattle Feces Shehla Shehla, Muhammad Kashif Obaid, Sadaf Niaz, Munir Ahmad Khan, Anum Ali Ahmad, Mostafa A. Abdel-Maksoud, Abdulaziz Alamri, Salman Alrokayan, Muhammad Shoaib and Qiaoyun Ren + 1 author

为了应对牲畜粪便作为抗生素耐药性（AMR）传播热点和新兴“同一个健康”威胁的挑战，研究人员对巴基斯坦三个地区的健康奶牛粪便进行了宏基因组学研究。该研究首次系统揭示了该地区奶牛粪便的微生物组成、多样性，检测到大量且多样化的抗菌素耐药基因（ARGs）和质粒，并证实了包括能抵抗最后防线抗生素（如tigecycline）的基因tet(X)在内的关键ARGs的存在。这项研究强调了奶牛肠道是重要的耐药基因储存库和水平基因转移（HGT）热点，为评估区域性的公共卫生与环境风险提供了重要的基线数据。

　　

奶牛不仅是重要的经济动物，其肠道和排出的粪便更是一个复杂而活跃的微生物“宇宙”。这个生态系统中居住着数以万亿计的细菌、古菌、病毒和真核微生物，它们不仅影响着奶牛的健康与生产力，还可能成为威胁人类和动物健康的“特洛伊木马”——抗生素耐药基因（ARGs）的储存库和传播载体。随着“同一个健康”理念的深入人心，人们越来越认识到，动物、人类和环境的健康紧密交织。奶牛粪便作为重要的有机肥来源和环境污染物，其中潜藏的耐药基因和机会性病原体，可能通过食物链、水体和直接接触等途径扩散，对公共卫生构成严峻挑战。然而，在巴基斯坦这样的农业国，关于奶牛粪便微生物群落及其耐药性潜力的系统性研究仍很匮乏，这限制了我们评估和应对这一区域性健康风险的能力。为了填补这一知识空白，一项名为《Shotgun Metagenomics Reveals Microbial Diversity, Resistome, and Plasmidome in Dairy Cattle Fees》的研究在《Veterinary Sciences》期刊上发表。该研究首次对巴基斯坦开伯尔-普赫图赫瓦省三个不同生态区的奶牛粪便进行了全面的宏基因组学“人口普查”，不仅描绘了其微生物群的“居民”图谱，更深入探查了其中暗藏的“军火库”（耐药基因）和“运载工具”（质粒），揭示了该地区作为抗生素耐药性热点区域的潜在风险。

为了开展这项研究，研究人员从巴基斯坦开伯尔-普赫图赫瓦省的三个地区（马尔丹、白沙瓦和德拉伊斯梅尔汗）的健康奶牛中，共收集了150份粪便样本。由于预算限制，每个地区的50份样本被等量混合，最终形成三个复合样本进行后续分析。研究采用了几项关键技术：首先，使用Illumina平台进行鸟枪法宏基因组测序，以获得样本中所有微生物的DNA序列信息。接着，利用生物信息学流程对测序数据进行质量控制、组装和基因预测。然后，通过比对微生物非冗余蛋白数据库进行物种注释，并使用ResFinder和PlasmidFinder工具分别鉴定抗菌素耐药基因和质粒。此外，还通过比对CAZy数据库分析了碳水化合物活性酶（CAZymes）的功能谱。这些方法共同构成了从样本到数据的完整分析链条。

研究结果

3.1. 测序数据概览

研究共产生了16,300 Mbp的原始数据，经过质量控制后获得16,030 Mbp的干净数据，平均有效数据率为98.39%。组装和基因预测后，共识别出797,541个开放阅读框（ORFs）。

3.2. 核心-PAN基因

核心基因分析显示，三个样本（FC1, FC2, FC3）共享248,137个基因，同时每个样本也有其特有的基因库，反映了微生物组成的共性与地区特异性。

3.3. 粪便微生物组成

微生物组成分析显示，在界水平，细菌占绝对主导（84.00%–91.00%），其次是病毒、古菌和真核生物。在门水平，假单胞菌门、芽孢杆菌门和拟杆菌门是主要类群。在物种水平，宽厄菲菌、大肠杆菌、普罗威登斯菌和斯氏弓形杆菌等是优势或高丰度物种。其中，大肠杆菌、普罗威登斯菌和斯氏弓形杆菌等机会性病原体的较高丰度，可能暗示了肠道内存在有利于其定植的亚临床环境。从地区差异看，白沙瓦（FC2）样本观察到的物种数量最多，这可能反映了农场管理、当地抗生素使用、饮食组成或环境暴露的潜在差异。

3.4. 粪便微生物多样性

Alpha多样性分析（香农指数、辛普森指数、ACE丰富度估计量）表明，从高级到低级的分类水平，微生物多样性和丰富度逐渐增加。Bray-Curtis相异度聚类树显示，FC1和FC2样本的微生物组成更为相似，而与FC3差异较大，这表明不同地区奶牛粪便的微生物群落结构存在差异。

3.5. 抗菌素耐药基因与质粒

耐药组分析在所有样本中均检测到大量耐药基因（40-49个），涵盖了β-内酰胺类、氨基糖苷类、磺胺类、甲氧苄啶、氟喹诺酮类、酚类等多种抗生素类别。尤为关键的是，所有样本均检出了可导致对最后防线抗生素替加环素耐药的tet(X)基因，以及重要的超广谱β-内酰胺酶基因如bla_CTX-M-15。质粒组分析显示，样本主要携带Col和IncF家族质粒，其中在FC3样本的一个Col440I质粒相关重叠群中还发现了qnrB19基因，这为耐药基因通过质粒进行水平基因转移（HGT）提供了直接证据。这些发现表明，奶牛肠道是活跃的耐药基因储存和交换热点。

3.6. CAZyme鉴定

功能分析显示，糖苷水解酶（GHs）和糖基转移酶（GTs）是粪便微生物群中最主要的碳水化合物活性酶（CAZymes），其次是碳水化合物结合模块（CBMs）。具体而言，GT2、GT4、GH13和CBM50的丰度在所有样本中均较高。这些酶在复杂碳水化合物的降解、微生物细胞壁构建和代谢中起着关键作用，反映了粪便微生物群落强大的碳水化合物代谢潜能。

研究结论与重要意义

本研究首次为巴基斯坦开伯尔-普赫图赫瓦省三个地区的奶牛粪便微生物组建立了基线图谱。结论明确指出，该地区奶牛的粪便微生物组并非简单的废物，而是一个复杂且具有重要“同一个健康”意义的生态系统。研究不仅揭示了高微生物多样性及可能受地理和环境因素影响的组成差异，更关键的是，发现了广泛且多样化的耐药组和质粒组。机会性病原体在粪便群落中的突出地位，提示了牛群中可能存在需要进一步调查的亚临床胃肠道健康问题。而大量耐药基因和质粒的检出，特别是关键基因如tet(X)和bla_CTX-M-15，以及质粒上qnrB19基因的发现，强有力地证实了这些地区的奶牛是能够通过水平基因转移进行传播的抗生素耐药基因的重要储存库。

这项研究的意义重大。首先，它填补了巴基斯坦特定地区奶牛粪便微生物组和耐药组数据的空白，为后续研究提供了宝贵的比较基线。其次，它将奶牛肠道定位为一个动态的“生物反应器”和耐药基因交换“热点”，提升了我们对农业环境中抗生素耐药性（AMR）产生和传播机制的认识。最重要的是，研究结果发出了明确的公共卫生警报：即使在未直接施压抗生素的情况下，奶牛粪便中也可能携带并传播对临床最后防线药物构成威胁的耐药基因，这凸显了对畜牧业进行持续的宏基因组学监测、加强动物粪便管理以及制定基于“同一个健康”理念的AMR遏制策略的紧迫性。未来研究应扩大样本量、进行个体水平测序，并探索微生物特征与宿主健康、生产数据及具体管理实践之间的关联，从而将描述性发现转化为可应用的干预和管理策略。