宏基因组分析杜洛克×长白×大白猪不同生长阶段肠道微生物组成与抗生素抗性基因关联

《Environmental Technology & Innovation》：Metagenomic analysis of microbial composition and antibiotic resistance gene associations in Duroc × Landrace × Yorkshire pigs’ intestinal microbiota across different growth stages

　　

背景：随着集约化养猪业的快速发展和抗生素的广泛使用，猪的肠道微生物群落和抗生素抗性基因（antibiotic resistance genes，ARGs）因其与动物健康和环境安全的相关性而受到越来越多的关注。表征肠道微生物群和ARG谱的阶段性差异对于理解牲畜粪便中ARG对环境的潜在负荷以及为商业化猪生产提供基于证据的管理策略至关重要。
方法：本研究调查了来自集约化养殖场的不同生长阶段的杜洛克×长白×大白（Duroc × Landrace × Yorkshire，DLY）杂交猪，包括仔猪组（n=13）、育肥猪组（n=19）和母猪组（n=7）。收集粪便样本，并采用宏基因组测序结合生物信息学方法进行分析，以检查不同生长阶段肠道微生物组成和抗性基因的结构、功能和宿主归属。
结果：与仔猪和母猪相比，育肥猪表现出显著更高的基因数量和物种丰富度（P < 0.05）。厚壁菌门（Firmicutes）和拟杆菌门（Bacteroidetes）分别是仔猪组和育肥猪组中的优势菌门，而母猪组的变形菌门（Proteobacteria）显著增加，拟杆菌门丰度显著降低（P < 0.05）。Beta多样性分析证实了各组间存在统计学显著的结构分化（PERMANOVA，Adonis R² = 0.165，P = 0.008）。ARGs在育肥猪组中表现出最高丰度，tetQ和tetO分别是仔猪组和育肥猪组中的主要抗性基因，且tetO丰度显著高于其他组（P < 0.05）。ARGs主要归属于厚壁菌门、拟杆菌门和变形菌门，其中普雷沃氏菌属（Prevotella）是ARG相关序列中主要注释的属。相关性分析显示，大多数ARGs（tetQ、sul1和floR）与变形菌门和拟杆菌门类群呈正相关，而tetL抗性基因与链球菌属（Streptococcus）显示出高度特异性的正相关（P < 0.001），表明不同抗性基因与特定微生物类群之间存在推定的分类群-ARG关联。
结论：这项横断面研究表征了DLY猪在三个代表性生产阶段粪便微生物群和ARGs的组成及共丰度关联，为理解阶段性粪便ARG输出提供了基线数据，并为制定有针对性的管理策略以减少集约化养猪场ARG向环境扩散提供了依据。

# 论文解读：不同生长阶段杜洛克×长白×大白猪肠道微生物组成与抗生素抗性基因关联的宏基因组分析

## 研究背景与问题

随着集约化养猪业的快速发展，抗生素曾广泛用于促进生长和预防疾病。然而，抗生素的滥用导致抗生素抗性基因（antibiotic resistance genes，ARGs）在环境中大量积累和传播，成为全球关注的公共卫生和生态安全问题。猪粪作为ARGs的重要储库，施用于农田后可能污染土壤和作物，威胁人类健康。即使未接受抗生素治疗的健康动物粪便中也含有丰富的ARGs，表明ARG污染的普遍性和复杂性。猪肠道微生物群落随生长阶段呈现动态变化，受到营养需求、生理成熟和饲养条件等因素影响。先前研究多集中于特定生长阶段、地方品种或一般杂交猪，而对中国主导商业品种杜洛克×长白×大白（Duroc × Landrace × Yorkshire，DLY）杂交猪不同生产阶段的系统性横断面比较仍然缺乏，限制了对阶段性粪便ARG谱及其潜在分类学关联的理解。为此，研究人员开展了此项研究。

## 研究内容与意义

研究人员从一家大型集约化商业农场采集了DLY猪三个代表性生产阶段（仔猪13头、育肥猪19头、母猪7头）的粪便样本，采用宏基因组测序（shotgun metagenomic sequencing）结合生物信息学方法，全面分析了不同生长阶段肠道微生物群落的结构、功能以及ARG的分布特征与宿主归属。研究揭示了生长阶段对微生物组成和ARG谱的显著影响，为评估阶段性环境ARG负荷和制定阶段特异性粪便管理策略提供了基线数据。该论文发表于《Environmental Technology》。

## 关键技术方法

研究人员采用宏基因组测序方法对来自单一集约化农场的DLY猪三个生长阶段（仔猪、育肥猪、母猪）的粪便样本进行分析。样本采集后，使用QIAamp DNA Stool Mini Kit提取总DNA，经文库构建后在Illumina HiSeq平台进行双端150 bp测序。生物信息学分析包括：使用Readfq进行质量过滤，Bowtie2去除宿主基因组序列，MEGAHIT进行个体和混合组装，MetaGeneMark预测开放阅读框，CD-HIT构建非冗余基因集。分类注释通过DIAMOND比对NCBI NR数据库并结合LCA算法实现；ARG注释通过RGI比对CARD数据库（仅保留“Perfect”或“Strict”命中）。采用α多样性（Shannon、Chao1等）、β多样性（PERMANOVA）和LEfSe分析微生物群落差异，用Spearman相关分析评估ARG与微生物类群的关联。

## 研究结果

### 3.1 不同生长阶段DLY猪粪便微生物群落基因组组成、丰度和多样性分析
通过Venn图、α多样性和β多样性分析发现，育肥猪的基因总数和物种丰富度显著高于仔猪和母猪（P < 0.05）。厚壁菌门（Firmicutes）和拟杆菌门（Bacteroidetes）分别是仔猪和育肥猪的优势菌门；母猪组变形菌门（Proteobacteria）丰度显著升高，拟杆菌门丰度显著降低（P < 0.05）。β多样性分析（PERMANOVA，Adonis R² = 0.165，P = 0.008）确认了组间显著的结构分化。LEfSe分析显示拟杆菌相关类群和乳杆菌目在仔猪和育肥猪中富集，肠杆菌科等变形菌类群在母猪中更突出。

### 3.2 不同生长阶段DLY猪各分类水平细菌类群的聚类和差异丰度分析
通过热图聚类和Student's t检验，在多个分类水平识别出显著差异的细菌类群。门水平上，放线菌门在仔猪组中显著高于育肥猪（P < 0.01）；拟杆菌门在母猪组中显著低于育肥猪（P < 0.01）。属水平上，链球菌属在育肥猪中丰度显著高于其他两组（P < 0.05）；乳酸杆菌属和普雷沃氏菌属在母猪组中显著降低（P < 0.05）。种水平上，罗伊氏乳杆菌和猪肠乳杆菌在母猪组中显著低于其他阶段（P < 0.05）。

### 3.3 不同生长阶段DLY猪粪便微生物群落基因功能分析
通过比对KEGG、eggNOG和CAZy数据库进行功能注释，结果显示所有样本中代谢相关功能（KEGG Level 1）占主导，尤其是碳水化合物代谢和氨基酸代谢（Level 2）。CAZy酶家族中糖苷水解酶和糖基转移酶丰度最高。eggNOG注释中，除未知功能类别外，复制、重组和修复相关基因相对丰度最高。三个生长阶段的功能谱整体上大体相似。

### 3.4 不同生长阶段DLY猪粪便中抗生素抗性基因的分布特征
通过比对CARD数据库，发现ARG序列总数在育肥猪和母猪中显著高于仔猪（P < 0.05），但ARG类型数（ARO类别）组间无显著差异。ARG类别以四环素类、大环内酯类-林可酰胺类-链阳菌素类（macrolide-lincosamide-streptogramin，M-L-S）和氨基糖苷类为主。核心ARGs（≥80%样本检出）占总ARGs的60%以上。丰度前五的ARG亚型为tetQ、aph(3')-IIIa、tetW/N/W、tetO和tetW。tetO在育肥猪中相对丰度显著高于其他两组（P < 0.05）；aph(3')-IIIa在仔猪中最高。PCA和PCoA显示各组ARG谱部分分离但重叠显著。

### 3.5 不同生长阶段DLY猪ARGs的聚类和差异丰度分析
通过热图聚类和t检验，发现特定ARG亚型在生长阶段间呈现显著差异。tetL和tetO在育肥猪中显著高于仔猪和母猪（P < 0.05）；aph(3')-IIIa、tetQ、fexA等ARGs在仔猪中显著高于母猪（P < 0.05）；cmx在仔猪中显著高于育肥猪（P < 0.05）。这些结果表明ARG分布具有阶段性特征。

### 3.6 不同生长阶段DLY猪ARGs的微生物宿主归属和分类关系
通过将ARG注释序列映射到分类注释，发现ARGs主要归属于厚壁菌门、变形菌门和拟杆菌门，且各门比例在生长阶段间呈现差异：变形菌门相关ARGs比例逐渐升高，拟杆菌门相关ARGs在母猪组中最低。从属级看，普雷沃氏菌属、梭菌属和拟杆菌属是主要注释属。

### 3.7 不同生长阶段DLY猪ARGs与微生物群落的Spearman相关性分析
Spearman相关分析显示，多个ARGs（包括AAC3-Ib、cmx、floR、mphD、sul1和APH6-Id）与变形菌门及其下属类群（γ-变形菌纲、假单胞菌科）呈显著正相关，与链球菌科、链球菌属及普雷沃氏菌科呈负相关。tetQ与拟杆菌门、普雷沃氏菌科、普雷沃氏菌属呈显著正相关（P < 0.001）。tetL抗性基因与链球菌科及链球菌属呈高度特异性的正相关（P < 0.001），与其他微生物类群均负相关。

## 讨论总结与结论翻译

讨论部分指出，本研究采用横断面设计，观察到的差异反映的是与生产阶段相关的组间变异，而非个体纵向变化。若干与生长阶段相混淆的因素（如日粮配方、体重、繁殖状态、饲养条件及既往抗生素暴露史）无法完全剥离。研究结果与先前报道一致，如拟杆菌门丰度在母猪组中最低，物种丰富度随生产阶段升高；功能分析显示代谢相关功能占主导，且在三个阶段间宏观相似。ARG方面，检测到的主要类别（四环素、MLS、氨基糖苷类）和优势亚型（tetQ、aph(3')-IIIa等）与先前猪粪抗性组研究一致。然而，与Gaire等人纵向研究中ARG丰度随年龄下降的趋势相反，提示研究设计（横断面对比纵向）和管理背景对结果解释至关重要。研究还讨论了分类注释和相关性分析的局限性：基于短读长的注释仅提供推定宿主关系；Spearman相关性不构成基因组共定位或水平转移证据。功能分析仅作为描述性基线，未进行正式差异富集检验。未来需结合纵向队列设计、长读长测序、可移动遗传元件分析和环境配对采样来验证和扩展这些发现。

**结论翻译**：
这项横断面宏基因组研究评估了杜洛克×长白×大白（DLY）猪不同生产阶段粪便微生物群落和抗生素抗性基因（ARGs）的谱系。结果表明，育肥猪的估算物种丰富度在统计学上高于仔猪。此外，厚壁菌门和拟杆菌门在仔猪和育肥猪的群落结构中占主导，而母猪组呈现明显的组成性变化，表现为变形菌门相对丰度增加和拟杆菌门降低。同时，ARGs的总序列计数在评估的阶段间存在显著差异，年龄较大的队列通常表现出更高水平。核心ARGs在所有组中均被检测到，占总抗性基因序列的60%以上。特定抗性谱也表现出可测量的阶段性变异。例如，tetQ基因在仔猪和育肥猪中保持最高相对丰度，而tetO比例在育肥猪队列中显著升高。ARG携带序列的分类注释主要归类为厚壁菌门、拟杆菌门和变形菌门，其中普雷沃氏菌属是主要注释属。随后的相关性分析确定了抗性基因与微生物类群之间不同的统计共丰度模式，尤其突出了tetL基因与链球菌属之间的排他性共现。总体而言，这些描述性结果记录了不同生长阶段肠道微生物组成和抗性组谱的结构变异，为未来对猪生产过程中抗菌素耐药性动态的纵向研究提供了基线观察数据。