《Microorganisms》：Exploratory Shotgun Metagenomic Analysis of the Fecal Microbiome and Antimicrobial-Resistance Potential in Simmental Calves with Different Diarrheal Statuses
编辑推荐：
摘要：犊牛腹泻常伴随胃肠道微生物群的紊乱，但临床状态与微生物群落组成、基因含量功能潜力以及抗生素耐药性相关特征之间的关系仍不完全清楚。本探索性研究旨在利用鸟枪法宏基因组学（shotgun metagenomics）表征不同腹泻状态犊牛的这些特征。研究对象包括来
摘要：犊牛腹泻常伴随胃肠道微生物群的紊乱，但临床状态与微生物群落组成、基因含量功能潜力以及抗生素耐药性相关特征之间的关系仍不完全清楚。本探索性研究旨在利用鸟枪法宏基因组学（shotgun metagenomics）表征不同腹泻状态犊牛的这些特征。研究对象包括来自内蒙古察哈尔右翼中旗的六头约2月龄的西门塔尔犊牛，其中两头健康犊牛、两头轻度腹泻犊牛和两头重度腹泻犊牛。研究人员将基于读长的分类学分析与基于基因丰度加权的分类注释相结合，并使用eggNOG/COG、KEGG和CAZy注释以及经过保守筛选的CARD关联序列同源候选物来表征微生物的功能潜力。两个重度腹泻样本表现出较低的细菌丰富度和香农多样性（Shannon diversity），尽管主坐标分析（PCoA）并未显示根据临床严重程度的清晰分离。这两个样本显示出变形菌门（Proteobacteria）和肠杆菌科（Enterobacteriaceae）相关分类单元的高代表性，特别是埃希氏菌属（Escherichia）相关的分类单元，同时瘤胃球菌科（Ruminococcaceae）、粪杆菌属（Faecalibacterium）和普雷沃氏菌属（Prevotella）的代表性较低。基因含量谱显示，与运输、环境感知、运动性、分泌和细菌细胞包膜过程相关的选定功能代表性较高，而与聚糖和复杂碳水化合物相关的若干功能代表性较低。选定的CARD关联候选物在重度腹泻样本中也更为丰富。在这个小型探索性队列中，两个重度腹泻样本因此显示出一致的分类学和基因含量差异；然而，这些观察结果不应被解释为通用的严重程度特异性微生物组标志物，并且需要在更大的纵向队列中进行验证。
论文解读：探索性鸟枪法宏基因组学揭示西门塔尔犊牛不同腹泻状态下粪便微生物组与抗生素耐药性潜能的协同变化
1. 研究背景、现存问题与研究目的
犊牛腹泻是断奶前阶段导致发病率、死亡率、生长迟缓和经济损失的主要原因。其病因复杂，涉及肠道病原体、宿主免疫、营养、管理和环境条件等多种因素的相互作用。除了已知的传染性病原体外，胃肠道微生物群日益受到关注，因为生命早期的微生物定殖有助于营养利用、黏膜免疫发育、定殖抵抗力和胃肠道成熟。犊牛肠道微生物群在出生后迅速发育，并在断奶前后经历显著的年龄依赖性变化，个体间差异巨大。这使得早期断奶前成为一个特别动态的阶段，在此阶段微生物群落发育的紊乱可能与胃肠道疾病相关。
现有证据表明，犊牛腹泻伴随着粪便微生物多样性和群落组成的显著改变。许多研究报道了腹泻犊牛中细菌多样性降低，以及与粪杆菌属、瘤胃球菌科、普雷沃氏菌属、肠杆菌科和埃希氏菌/志贺氏菌属相关分类单元的变化。然而，不同研究和农场之间存在相当大的异质性，表明腹泻不太可能对应单一的普遍保守的分类学特征。此外，大多数早期研究依赖于16S rRNA基因测序，主要集中在细菌群落组成上。更近期的鸟枪法宏基因组学和多组学研究显示，腹泻相关的菌群失调也可能涉及微生物基因含量和代谢潜力的实质性变化。然而，在同一粪便宏基因组中，分类学组成、基因含量功能潜力、碳水化合物活性酶谱和抗生素耐药性相关特征之间的关系仍未得到充分表征，尤其是在单个样本水平上的整合分析。本研究旨在通过对健康、轻度和重度腹泻的西门塔尔犊牛进行鸟枪法宏基因组测序，综合评估这些方面，以确定不同腹泻状态的犊牛是否在粪便微生物群落组成和基因含量功能潜力上表现出相应的差异。
2. 主要关键技术方法
研究人员从内蒙古乌兰察布市察哈尔右翼中旗放牧的西门塔尔犊牛中采集粪便样本，共纳入6头约2月龄、体重相近的犊牛，分为健康组、轻度腹泻组和重度腹泻组。采用E.Z.N.A. stool DNA Kit提取粪便DNA，构建配对末端测序文库后，在Illumina NovaSeq 6000平台上进行鸟枪法宏基因组测序。测序数据经Trimmomatic v0.39进行质量控制，随后使用Kraken2 v2.0.6进行基于读长的分类学分析。同时，使用MEGAHIT v1.1.1进行从头组装，MetaProdigal v2.6.3预测开放阅读框，CD-HIT v4.8.1构建非冗余基因目录，并使用Salmon v1.1.0估计基因丰度。对预测蛋白序列进行了多重数据库注释：使用DIAMOND v0.9.22.123比对NCBI非冗余蛋白数据库进行基因中心分类注释；使用eggNOG-mapper v5.0进行eggNOG/COG功能注释；使用KofamScan v1.2.0进行KEGG注释；使用HMMER hmmscan v3.3.2比对dbCAN-HMMdb-V8进行CAZy注释；并使用BLASTP v2.9.0+比对CARD数据库，经保守筛选后鉴定抗生素耐药性相关序列同源候选物。
3. 研究结果
3.1 测序与质量控制
鸟枪法宏基因组测序共产生约86.24 Gb原始数据，经质量过滤后保留约59.58 Gb清洁数据。各样本的原始数据量和清洁数据量存在较大差异。宿主过滤后，非宿主读段比例在不同样本间变化显著。进入下游分类分析的读段中，细菌分配比例在两个重度腹泻样本中最高，分别达到80.03%和75.09%，而未分类读段比例则相应最低。
3.2 粪便微生物多样性与群落结构
通过对细菌属水平丰度谱进行稀疏化处理，发现两个重度腹泻样本的观测属丰富度和香农多样性均低于健康组和轻度腹泻组样本。基于Bray-Curtis相异度的主坐标分析未显示根据临床严重程度的清晰分离。虽然两个重度腹泻样本在排序空间中相对接近，但一个健康样本与它们的距离甚至更近，表明个体间存在巨大变异。
3.3 粪便细菌群落分类学组成
在门水平，变形菌门在两个重度腹泻样本中占主导地位，分别占细菌分类读段的86.76%和84.25%，而厚壁菌门和拟杆菌门的比例极低。在属水平，两个重度腹泻样本均以埃希氏菌属为主导，分别占77.29%和61.28%。相比之下，健康和轻度腹泻组的分类学组成呈现出高度的个体间异质性。
3.4 重度腹泻中肠杆菌科相关类群的扩张与特定厌氧类群的耗竭
进一步分析显示，重度腹泻样本中的变形菌门优势主要由肠杆菌科相关类群贡献，其在S1和S2中分别占84.75%和79.35%。与此同时，多个在非重度样本中常见的厌氧菌群，如瘤胃球菌科、粪杆菌属和普雷沃氏菌属，在两个重度腹泻样本中均显著减少。基于基因丰度加权的NR数据库分类分析也支持了这一总体分类模式。
3.5 eggNOG/COG和KEGG注释揭示的重度腹泻微生物组功能谱
eggNOG/COG注释结果显示，与细胞运动性、无机离子转运和代谢、次级代谢产物生物合成以及细胞内运输、分泌和囊泡运输相关的基因在两个重度腹泻样本中一致性地高于所有四个非重度样本。KEGG Level 2分析表明，重度腹泻样本中与膜运输、信号转导、细胞运动性和细胞群落-原核生物相关的基因丰度更高。在KEGG Level 3通路水平，ABC转运蛋白、磷酸转移酶系统、双组分系统、鞭毛组装、细菌趋化性、细菌分泌系统、脂多糖生物合成和铁载体群非核糖体肽的生物合成等八个通路在重度腹泻样本中丰度更高，而其他聚糖降解通路则呈现相反趋势。
3.6 CAZyme家族谱揭示的碳水化合物和细胞包膜相关功能潜力
CAZy注释显示，糖苷水解酶是最主要的CAZy类别。在家族水平上，与细菌细胞包膜生物合成或重塑相关的多个CAZy家族，如GT73、GT9、GH23、GH102和GH103，在重度腹泻样本中丰度更高。相反，与复杂碳水化合物利用相关的家族，如GH5、GH29、CBM37，以及纤维小体相关特征SLH和cohesin，则在重度腹泻样本中丰度较低。
3.7 六个粪便宏基因组中CARD关联序列同源候选物
经保守筛选后，共获得133个独特的目录查询基因作为CARD关联序列同源候选物。两个重度腹泻样本中这些候选物的总丰度最高。其中，CTX-M-65、TEM-1、APH(6)-Id、APH(3'')-Ib、aadA5、sul1、floR和mphA等八个CARD参考标签在两个重度腹泻样本中的丰度均高于所有四个非重度样本。
4. 讨论与结论
讨论总结：研究表明，重度腹泻犊牛的粪便微生物群落结构和宏基因组功能潜力发生了协同变化，其最显著的特征并非单一分类单元或通路，而是细菌多样性降低、肠杆菌科主导、多种厌氧菌群耗竭、与环境感知和运动等功能相关的基因含量增加、碳水化合物活性酶谱重塑以及特定抗生素耐药性基因丰度升高这几种特征的共存。尽管两个重度腹泻样本在这些互补的分析层面展现出相对一致的图谱，但主坐标分析并未支持一个由临床严重程度定义的独特群落簇，强调了该小型队列中的巨大个体间变异。研究结果与先前关于犊牛腹泻相关菌群失调的报告基本一致，尤其是肠杆菌科的扩张与某些厌氧菌群的减少。这种分类学变化与功能潜力的改变在生物学上是合理的，例如革兰氏阴性菌扩张可解释与细胞包膜、运动和铁摄取等相关功能的增强。同时，与复杂碳水化合物利用相关的CAZy家族丰度降低，与瘤胃球菌科等厌氧碳水化合物利用菌群的减少相符。对于抗生素耐药性潜能，重度腹泻样本中更高的耐药基因丰度可能与其肠杆菌科的优势有关，但由于缺乏详细的治疗记录，无法区分治疗、临床状态和微生物群落组成各自的贡献。
研究结论：本研究的主要贡献在于对同一小型队列中的多个宏基因组层进行了整合表征。研究人员将基于读长和基于基因丰度加权的分类学图谱与COG/KEGG、CAZy和CARD关联基因含量模式结合在同一宏基因组数据集中进行分析。在两个重度腹泻样本中观察到的一致跨层模式提供了假设生成的证据，表明群落组成、编码功能潜力和抗生素耐药性相关遗传潜力之间存在关联，这些关联可以在更大的纵向队列中进行检验。
