畜牧兽医学报  2024, Vol. 55 Issue (4): 1605-1614. DOI: 10.11843/j.issn.0366-6964.2024.04.023    PDF    
基于16S rRNA测序分析高低饲料转化率猪粪便微生物的组成差异
张帅1, 陈奎蓉1, 许迪1, 江山1, 王梦影1, 张坤1, 徐玉培1, 雷国凤1, 张志程1, 郭猛2, 赵云翔1,2, 兰干球1, 梁晶1     
1. 广西大学动物科学技术学院 广西畜禽繁育与疾病防控重点实验室,南宁 530004;
2. 广西扬翔股份有限公司,贵港 537100
摘要:提高猪饲料转化率是目前猪育种工作的重点,本研究对384头三元杂交母猪的FCR值进行排序,从中选取两端高饲料转化率(HFCR组)和低饲料转化率(LFCR组)猪各10头,采集粪便进行微生物16S rRNA基因测序分析。结果表明,LFCR组的Richness指数、Chao1指数、ACE指数均极显著高于HFCR组(P<0.01)。通过LEfSe分析,在属水平上筛选两组间具有显著差异的微生物,其中,在LFCR组筛选LDA>3的微生物主要是Prevotella_1、Treponema_2、Prevotellaceae_NK3B31_groupEubacterium_coprostanoligenes_groupp-1088-a5_gut_groupPrevotella_9、阿克曼菌属、Prevotellaceae_UCG_009、Ruminococcaceae_UCG-014、Ruminococcaceae_NK4A214_group,在HFCR组LDA>3的微生物主要是巨球藻属、小类杆菌属、链型杆菌属、假丁酸弧菌属、柯林斯氏菌、Solobacterium、罕见小球菌属、链球菌属。Spearman相关性分析结果表明,与FCR显著负相关的微生物大多与短链脂肪酸的产生有关,而与FCR显著正相关的微生物多数为潜在的致病性细菌。本研究通过筛选高、低饲料转化率猪粪便中的差异微生物,并对差异微生物丰度与饲料转化率性状进行相关性分析,最终筛选到与饲料转化率显著关联的粪便微生物,为猪饲料转化率性状的肠道微生物育种标记筛选提供一定参考。
关键词饲料转化率    16S rRNA    粪便微生物    
Analysis of Microbial Composition Differences in High and Low Feed Conversion Rates Pig Feces based on 16S rRNA Sequencing
ZHANG Shuai1, CHEN Kuirong1, XU Di1, JIANG Shan1, WANG Mengying1, ZHANG Kun1, XU Yupei1, LEI Guofeng1, ZHANG Zhicheng1, GUO Meng2, ZHAO Yunxiang1,2, LAN Ganqiu1, LIANG Jing1     
1. Guangxi Key Laboratory of Animal Breeding, Disease Control and Prevention, College of Animal Science and Technology, Guangxi University, Nanning 530004, China;
2. Guangxi Yangxiang Co., Ltd, Guigang 537100, China
Abstract: Improving feed conversion efficiency in pigs is currently a key focus of pig breeding efforts. In this study, FCR values of 384 crossbred sows were ranked, and 20 pigs with extreme feed conversion ratio (10 with the highest FCR and 10 with the lowest FCR) were selected to form the HFCR group and LFCR group, respectively. Fecal samples were collected for 16S rRNA gene sequencing analysis. The results showed that the Richness index, Chao1 index, and ACE index of the LFCR group are significantly higher than those of the HFCR group (P < 0.01). The LEfSe analysis was used to identify microorganisms at the genus level with significant differences between the two groups. In the LFCR group, microorganisms with LDA>3 were mainly included Prevotella_1, Treponema_2, Prevotellaceae_NK3B31_group, Eubacterium_coprostanoligenes_group, p-1088-a5_gut_group, Prevotella_9, Akkermansia, Prevotellaceae_UCG_009, Ruminococcaceae_UCG-014, and Ruminococcaceae_NK4A214_group. In the HFCR group, microorganisms with LDA>3 were mainly included Megasphaera, Dialister, Catenibacterium, Pseudobutyrivibrio, Collinsella, Solobacterium, Subdoligranulum, and Streptococcus. The results of Spearman correlation analysis indicated that microorganisms significantly negatively correlated with FCR are mostly associated with short-chain fatty acid production, while microorganisms significantly positively correlated with FCR were mostly potential pathogenic bacteria. This study screened differential microorganisms in feces from pigs with high and low feed conversion rates, also conducted correlation analysis between the abundance of differential microorganisms and feed conversion rate traits, and ultimately identified fecal microorganisms significantly associated with feed conversion rate, thus provided a reference for the selection of intestinal microbial breeding markers for pig feed conversion rate traits.
Key words: feed conversion rate    16S rRNA    fecal microbiota    

畜牧业是国家发展的重要基石,在西方发达国家农业经济中,畜牧业发展占总体70%左右。猪是人类最早驯化的家养动物,同时也是肉类产品的主要来源[1]。截止到2019年,全球猪肉消费量达到新高,占总肉类消费的33%,据推测, 到2060年,全球人口将急剧增长,这也意味着人类对猪肉需求量会更高。如何快速提高猪肉生产,满足人类需求已经成为全世界亟待解决的问题[2]

在整个养猪过程当中,饲料成本约占养猪成本的60%,猪的饲料利用效率在很大程度上影响着养殖生产和总体经济效益[3],因此提高饲料转化率是目前猪育种工作的一项重要任务。目前,衡量饲料利用效率的指标有两个:剩余采食量(residual feedintake,RFI)和饲料转化率(feed conversion rate,FCR)[4],RFI是实际采食量与预期采食量之间的差值[5];FCR在猪饲养中统称为“料肉比”,即平均日采食量与平均日增重的比值(ADFI/ADG)[6]。随着自动喂料系统的发展,猪饲料转化率记录更加准确,使得猪育种工作更加方便快捷[7]

饲料转化率受多个因素影响。研究表明,猪饲料转化率受到肠道微生物活性和组成的影响[8],特别是肠道微生物能够利用宿主不能消化的碳水化合物。例如纤维素和抗性多糖,其先被微生物降解为单糖或二聚糖,再经过发酵形成短链脂肪酸(SCFAs)[9],这些代谢产物能被宿主吸收,并提供能量[10]。基于16S rRNA发现梭状杆菌(Cloacibacillus)、幽门螺杆菌(Helicobacter)和振荡杆菌(Oscillibacter)在低FCR的鸡中更为丰富[11],Tan等[12]研究发现,高FCR猪的普雷沃氏菌CAG:604在盲肠显著富集;McCormack等[13]研究发现,低RFI猪肠道微生物某些特定菌增加与健康和饲料效率有关。寻找到更多的与FCR相关的肠道微生物有助于揭示微生物功能与宿主的生理和代谢之间的关系,使得FCR表征更加全面。

目前仍然不清楚,肠道微生物如何影响猪饲料转化率。本研究通过微生物16S rRNA基因测序技术,对高、低饲料转化率猪粪便微生物的组成结构差异进行对比分析,旨在探究猪粪便中的微生物与猪饲料转化率之间可能存在的关联。这些微生物可能作为一种生物标记物,为在育种中应用微生物标记提高猪饲料转化率的实践提供一定的理论基础。

1 材料与方法 1.1 动物和试验设计

本研究中使用的猪均来自广西扬翔股份有限公司。选择健康状况良好、断奶批次相同、体重接近(23~27 kg)的三元杂交母猪384头,于断奶后转至同一栋大型猪舍中饲养,每栏饲养12~15头,饲喂不含抗生素或药物的标准玉米-豆粕型饲料,自由采食和饮水。利用自动喂料系统(HAMOER Technology Co.;Ltd., Tianjin)在猪进食时识别其电子耳标并测定猪的采食量和体重,计算90~165日龄期间的FCR值。将试验猪的FCR值从高到低进行排序,选择FCR极端高、低个体各10头,分别定义为HFCR组和LFCR组。为减少环境因素对肠道微生物测序结果的影响,试验猪达到170日龄时在测定栏位采集粪便,时间限定在上午8:00—10:30之间。为防止采样时猪产生应激情况,使用一次性PE手套接取肛门内的粪便样品,清理表面后置于冻存管内并放入液氮中保存,后转入实验室-80 ℃冰箱长期保存用于后续分析。

1.2 粪便微生物分析

1.2.1 微生物总DNA提取和PCR反应   使用QIAamp Fast DNA粪便提取试剂盒(51604;Qiagen,德国)从每头猪粪便样本中分离并收集总细菌DNA。提取完成后,用Qubit检测DNA浓度,用1%琼脂糖凝胶电泳检测抽提DNA的完整度。利用带有barcode的特异引物341F(5′-CCTACGGG-NGGCWGCAG-3′)和806R(5′-GGACTACHVGG-GTWTCTAAT-3′)进行PCR扩增得到16S V3-V4扩增子区域。PCR条件为:95 ℃预变性3 min;95 ℃ 30 s,55 ℃ 30 s,72 ℃ 15 s,从第二步开始25个循环;72 ℃延伸5 min,4 ℃保存。PCR扩增完成后,用磁珠对产物进行纯化,用Qubit检测产物浓度, 并用1.5%琼脂糖凝胶电泳检测PCR产物大小是否与预期相符,PCR产物目的条带大小正确,浓度≥2 ng·μL-1,可进行后续试验。

1.2.2 文库构建和数据预处理   对上述PCR产物进行混样建库,加上测序接头,进行PCR扩增后用磁珠分选去掉引物二聚体等小片段,并记录文库index信息。反应条件:98 ℃预变性45 s;98 ℃ 15 s,60 ℃ 30 s,72 ℃ 30 s,8个循环;72 ℃延伸10 min,4 ℃保存。使用Qubit对文库浓度进行定量并计算浓度,文库浓度≥4 nmol·L-1,可正常上机。文库检测合格后,按照有效浓度要求将不同文库pooling后上机测序,扩增子测序平台为Illumina Miseq(MiSeq;Illumina Inc.,美国),测序策略为PE250。测序数据Q30(%)>85%满足下机要求,下机数据根据双端barcode拆分数据,并切除barcode信息,切除双端引物后的测序序列进行扩增子流程分析。

1.2.3 扩增子流程分析   扩增子生物信息学采用EasyAmplicon v 1.0进行分析[14]。原始的Illumina fastq文件经过质控(错误率小于1%)后,使用VSEARCH去冗余设置参数为min_unique_size为59,采用Usearch对AVSs进行100%相似度比较,然后去掉其中的质体和非细菌,设置参数为:sintax_cutoff 0.1。用Vegan包过滤低丰度ASVs后将所有样本以最小序列条数抽平获得ASVs丰度表。根据SILVA数据库(https://www.arb-silva.de/)进行物种注释,物种注释结果用于扩增子分析[15]

α多样性指数计算采用Vegan包,使用抽平后的数据得出Richness、Chao1、ACE、Shannon、Simpson和Invsimpson指数。基于bray_curtis距离绘制主坐标分析(PCoA)图。用柱状图展示top 5门水平高丰度菌、top 11科水平高丰度菌。使用LEfSe(Linear discriminant analysis Effect Size,线性判别分析)分析筛选LDA>3的属作为两组差异菌属进行讨论。使用PICRUSt对两组粪便微生物菌群功能进行预测[16],STAMP软件展示两组微生物功能有何差异。

1.3 统计分析

采用SPSS(version 26.0; SPSS Inc., USA)分析20头猪表型差异,结果均用“平均数±标准差(mean±SD)”方式表示。采用T-test检验计算P值,P<0.05表示差异显著,P<0.01表示差异极显著。使用SPSS软件中的Spearman相关检验计算表型与差异菌群之间相关性。

2 结果 2.1 高、低饲料转化率猪表型差异

高、低饲料转化率猪生长性状见表 1。与HFCR组相比,LFCR组的ADG显著提高(P<0.01),而ADFI和FCR都显著低于HFCR组(P<0.01)。

表 1 高、低饲料转化率猪的生长性能 Table 1 Growth performance of pigs with high and low feed conversion rates
2.2 高、低饲料转化率猪粪便微生物α多样性分析

α多样性主要用于表示物种丰富度和均匀度的指标,结果表明(表 2),LFCR组的Richness、Chao1和ACE指数都显著高于HFCR组(P<0.01)。而LFCR组的Shannon指数高于HFCR组,但没有达到显著性差异。

表 2 高、低饲料转化率猪的粪便微生物α多样性指数 Table 2 Feces microbial α-diversity index of pig with high and low feed conversion rates
2.3 高、低饲料转化率猪粪便微生物β多样性分析

基于Bray_Curtis距离计算出高低饲料转化率猪粪便微生物的结构相似性情况。PCo1的解释率为28.48%,PCo2的解释率为18.38%,高、低饲料转化率猪粪便微生物分布呈现左右分离状态,且解释率之和大于40%,表明高低饲料转化率分组间粪便微生物的整体组成结构差异较大(图 1)。

图 1 高、低饲料转化率猪的粪便微生物β多样性分析 Fig. 1 Feces microbial β diversity analysis of pig with high and low feed conversion rates
2.4 高、低饲料转化率猪粪便微生物物种组成分析

物种堆叠图主要用于展示不同组之间物种组成的差异。在门水平上(图 2)平均相对丰度排名前5的分别是:厚壁菌门(Firmicutes), 各组相对丰度为:LFCR=45.36%,HFCR=60.65%;拟杆菌门(Bacteroideta),各组相对丰度为:LFCR=30.09%,HFCR=20.34%;广古菌门(Euryarchaeota),各组相对丰度:LFCR=10.55%,HFCR=9.31%;螺旋体菌门(Spirochaetae),各组相对丰度为:LFCR=5.38%,HFCR=2.80%;放线菌门(Actinobacteria),各组相对丰度为:LFCR=1.86%,HFCR=2.49%。

图 2 高、低饲料转化率猪粪便微生物的门水平物种组成 Fig. 2 Phylum level species composition of microorganisms in feces of pigs with high and low feed conversion rates

在科水平上(图 3),平均相对丰度前11的菌群主要为以下几种,瘤胃球菌科(Ruminococcaceae),各组相对丰度为:LFCR=14.96%,HFCR=12.67%。拟杆菌_S24-7(Bacteroidales_S24-7),各组相对丰度为:LFCR=11.91%,HFCR=12.64%。链球菌科(Streptococcaceae),各组相对丰度:LFCR=7.32%,HFCR=13.21%。甲烷杆菌科(Methanobacteriaceae),各组相对丰度:LFCR=10.54%,HFCR=9.31%。毛螺菌科(Lachnospiraceae),各组相对丰度:LFCR=7.74%,HFCR=11.38%。普氏菌科(Prevotellaceae),各组相对丰度:LFCR=13.88%,HFCR=4.22%。乳酸杆菌科(Lactobacillaceae),各组相对丰度:LFCR=5.16%,HFCR=8.12%。螺旋体科(Spirochaetaceae),各组相对丰度:LFCR=5.38%,HFCR=2.80%。丹毒丝菌科(Erysipelotrichaceae),各组相对丰度:LFCR=1.86%,HFCR=5.26%。氨基酸球菌科(Acidaminococcaceae),各组相对丰度:LFCR=2.82%,HFCR=3.33%。理研菌科(Rikenellaceae),各组相对丰度:LFCR=2.40%,HFCR=2.28%。

图 3 高、低饲料转化率猪粪便微生物的科水平物种组成 Fig. 3 Family level species composition of microorganisms in feces of pigs with high and low feed conversion rates
2.5 LEfSe分析筛选高、低饲料转化率猪粪便微生物属水平差异

LEfSe分析主要可以对两组或多组进行比较,从而找出组间在丰度上具有显著差异的物种。在属水平上,筛选平均相对丰度大于0.5%的属进行LEfSe分析,结果如图所示(图 4)。在LFCR组LDA值大于3的属有:Prevotella_1、Treponema_2、Prevotellaceae_NK3B31_groupEubacterium_coprostanoligenes_groupp-1088-a5_gut_groupPrevotella_9、阿克曼菌属、Prevotellaceae_UCG_003、Ruminococcaceae_UCG-014、Ruminococcaceae_NK4A214_group。在HFCR组LDA值大于3的属有:巨球藻属、小类杆菌属、链型杆菌属、假丁酸弧菌属、柯林斯氏菌、Solobacterium、罕见小球菌属、链球菌属。

图 4 利用LEfSe分析筛选高、低饲料转化率猪粪便中的差异微生物属 Fig. 4 Using LEfSe analysis to screen differential microbial genera in feces of pigs with high and low feed conversion rates
2.6 PICRUSt预测高、低饲料转化率猪粪便微生物功能差异

PICRUSt可以用于细菌功能预测。在LFCR组的铁载体组非核糖体肽的生物合成、氮素代谢、脂质生物合成蛋白、转录机制、PPAR信号通路、蛋白质消化吸收、鞭毛组装、细菌运动蛋白、乙苯基降解、苯丙氨酸代谢、萜、哌啶和吡啶生物碱的生物合成、叶酸合成、霍乱弧菌致病周期、维生素B6代谢、谷氨酸能突触、脂多糖生物合成蛋白、脂多糖生物合成和叶绿醇降解17个功能较强。而HFCR组的硫胺素代谢、磷酸转移酶系统、糖酵解/异生、原发性免疫缺陷4个功能较强(图 5)。

图 5 利用PICRUSt预测高、低饲料转化率猪粪便微生物的功能差异 Fig. 5 Using PICRUSt to predict microbial functional differences in feces of pigs with high and low feed conversion rates
2.7 高、低饲料转化率猪粪便微生物差异菌属与生长性状关联分析

将两组猪生长性状和差异微生物进行Spearman相关性分析用于寻找影响饲料转化率的微生物(图 6)。与FCR值呈显著负相关的微生物有:Prevotella_1(r=-0.789,P<0.001)、Treponema_2(r=-0.548,P=0.012)、Prevotellaceae_NK3B31_group(r=-0.681,P<0.001)、Eubacterium_coprostanoligenes_group(r=-0.848,P<0.001)、p-1088-a5_gut_group(r=-0.54,P=0.014)、Prevotella_9(r=-0.503,P=0.024)、阿克曼菌属(r=-0.561,P=0.015)、Prevotellaceae_UCG_003(r=0.685,P<0.001)、Ruminococcaceae_UCG-014(r=-0.55,P=0.012)、Ruminococcaceae_NK4A214_group(r=-0.59,P<0.01)。与FCR值呈显著正相关的微生物有:巨球藻属(r=0.52,P=0.019)、柯林斯氏菌(r=0.74,P<0.001)、Solobacterium(r=0.728,P<0.001)、罕见小球菌属(r=0.5,P=0.025)、链球菌属(r=0.628,P<0.01)。

红色代表正相关,蓝色代表负相关。*代表P<0.05,**代表P<0.01,***代表P<0.001 Red represents positive correlation, while blue represents negative correlation。* Represents P < 0.05, ** represents P < 0.01, and*** represents P < 0.001 图 6 粪便微生物与生长表型的相关性分析 Fig. 6 Correlation analysis between fecal microorganisms and growth phenotype
3 讨论

LFCR组猪的日采食量较低而日增重较高,通俗讲为“吃的少,长得快”的猪。这类猪不仅对饲料的需求少同时生长速率较快。因为这类猪有利于降低企业生产成本增加效率,故常在猪育种中被留种。探究影响猪饲料转化率的潜在因素,是目前猪育种工作关注的重心。

3.1 粪便微生物多样性与猪饲料转化率的关系

微生物α多样性指数用于反映微生物群落的丰富度和均匀性。有研究表明,微生物多样性与猪生长性能存在联系[17]。Tan等[18]研究表明,LFCR猪的肠道微生物α多样性比HFCR猪更为丰富,这与本研究一致。Metzler-Zebeli等[19]发现,低RFI猪的盲肠微生物多样性也高于高RFI猪。本研究中,LFCR组的Richness、Chao1和ACE指数显著高于HFCR组。这表明LFCR组的粪便中存在更为多样化的微生物群落和有益菌,可能促进了猪的健康状态或提高了代谢水平,从而有助于提高饲料转化率。

3.2 LEfSe分析揭示高、低饲料转化率猪粪便微生物的组成及潜在功能差异

本研究通过LEfSe分析了高、低饲料转化率猪粪便微生物组成差异。Eubacterium_coprostanogenes_group属于真菌属,有研究表明,其可以产生短链脂肪酸作为肠上皮中的特殊营养物质和能量成分,保护肠黏膜屏障,降低人体炎症水平,增强胃肠动力[20]Prevotella_1、Prevotella_9、Prevotelaceae_NK3B31_groupPrevotellaceae_UCG-003属于普雷沃氏菌科,是猪肠中纤维降解菌的主要类型之一,参与果胶和蛋白质的代谢[21]Prevotella_1可以降解半纤维素或木聚糖,为身体提供能量[22]Prevotellaceae_UCG-003具有调节肠道功能的潜力,可降解纤维多糖产生短链脂肪酸[21]。然而,Xue等[23]发现,高RFI组绵羊瘤胃中Prevotella_1和Prevotellacea-UCG-003的相对丰度较高,这与本试验结果相反。可能的原因是Prevotell_1在反刍动物和单胃动物中发挥作用不同所导致。Jiang等[24]发现Prevotelaceae_NK3B31_group在断奶仔猪粪便中的丰度较高,并且与乙酸、丙酸和总SCFAs的浓度呈正相关,表明Prevotelaceae_NK3B31_group可以产生SCFAs。

Ruminococcaceae_UCG-014和Ruminococcaceae_NK4A214属于瘤胃球菌科。瘤胃球菌科成员具有纤维素降解和半纤维素降解活性,可产生乙酸盐、丁酸盐、甲酸盐和氢气。在一项研究中发现,Ruminococcaceae_UCG-014在低RFI牛瘤胃中丰度较高,这与本研究结果相似[21]。有研究证实,Ruminococcaceae_UCG-014可以产生丁酸[25]。阿克曼菌属已被确定为与体重增加相关的重要粪便生物标志物[26]。在绝经后妇女中,阿克曼菌属与胰岛素抵抗和血脂异常呈负相关[27]。肥胖个体的粪便中阿克曼菌属显著减少[28]。也有研究表明,阿克曼菌属可以降解黏蛋白,也是短链脂肪酸的生产者[29]。这些真菌属、普雷沃氏菌和瘤胃球菌在LFCR组粪便中大量存在,可能是它们可以产生大量短链脂肪酸为机体供给能量所以提高了猪的饲料转化率。

在人类研究中发现,柯林斯氏菌与肥胖和动脉粥样硬化有关[30]。同样,在高胆固醇血症患者的粪便中也发现了较多的柯林斯氏菌[31]。Lahti等[32]研究表明柯林斯氏菌和甘油三酯之间存在显著的正相关性。Solobacterium属于丹毒丝菌科,研究表明丹毒科在肥胖患者和炎症性肠病患者中发挥着重要作用[33];另外一项研究揭示,Solobacterium是结直肠癌癌症的标志性微生物群,在结直肠癌癌症患者的粪便中具有更高的丰度。链球菌是猪的常见病原体,不仅被认为是人类的人畜共患病原体,也是抗生素耐药性的主要宿主,可以将耐药性基因传递给其他致病性链球菌物种[34]。研究指出,罕见小球菌属与抑郁症患病概率之间存在显著的正相关[35]。这些菌在HFCR组显著富集,提示HFCR组可能存在大量致病菌,这些致病菌破坏了机体健康所以造成了饲料转化率低。假丁酸弧菌属在大肠后部显著富集,并与短链脂肪酸含量呈正相关[36]。在一项研究中,发现HFCR长白猪盲肠中的假丁酸弧菌属和巨球藻属显著富集,因为它们可以代谢乳酸产生游离脂肪酸[37],令人感兴趣的是,为什么这些产短链脂肪酸的细菌会在HFCR组显著富集。

3.3 与猪饲料转化率相关的粪便微生物及其在育种中的应用前景

Eubacterium_coprostanogenes_groupPrevotella_1、Prevotella_9、Prevotelaceae_NK3B31_group、普Prevotellaceae_UCG-003、Ruminococcaceae-UCG-014、Ruminococcaceae_NK4A214和阿克曼菌属,这些菌都具有产短链脂肪酸的能力与FCR值呈显著负相关,可能产生的短链脂肪酸供给机体所以LFCR饲料转化率高。而柯林斯氏菌、Solobacterium、罕见小球菌属和链球菌都与疾病发生有关且与FCR值呈显著正相关,可能是因为这些致病菌在肠道大量存在,破坏了肠道内环境,损坏机体健康因此降低了饲料转化率。有研究表明,肠道微生物和粪便微生物具有显著差异[38],肠道菌群更能反映宿主的系统发生关系,而粪便微生物只能反映宿主的饮食代谢情况。肠道菌群能够影响的不仅有营养物的消化吸收,还会影响大脑活动。肠道菌群失衡会破坏人体的诸多系统,还会导致肥胖、营养不良甚至癌症[39]。因此,用粪便微生物与生长表型进行相关性分析,也只能说明两者之间存在显著相关性,不能直接说明其存在因果关系。本研究筛选到的粪便微生物如何调控饲料转化率,仍需要进一步的试验探究。而粪便微生物有望作为一种易于检测的分子标志物,为动物育种时辅助筛选优良性状个体提供一定参考。

4 结论

本研究通过采集20头高、低饲料转化率的三元杂交母猪的粪便样品,利用16S rRNA测序技术分析两组间粪便微生物的组成与结构差异。结果显示,低饲料转化率组的粪便微生物α多样性显著高于高饲料转化率组。进一步分析发现,低饲料转化率组富集了大量产生短链脂肪酸的微生物,而这些微生物与饲料转化率值呈显著负相关。相反,高饲料转化率组的粪便中富集了大量的致病菌,这些菌与饲料转化率值呈显著正相关。最终,本研究鉴定了一些与饲料转化率相关的粪便微生物,这些微生物有望成为分子标记物,为促进猪育种进展、改善猪饲料转化率提供一定的理论基础。

参考文献
[1]
贾云飞, 郑伟程, 何泽军, 等. 我国畜牧业绿色发展: 现状、制约与战略应对[J]. 中国畜牧杂志, 2023, 59(1): 323-329.
JIA Y F, ZHENG W C, HE Z J, et al. Green development of animal husbandry in China: current situation, constraints, and strategic responses[J]. Chinese Journal of Animal Science, 2023, 59(1): 323-329. (in Chinese)
[2]
谢铿铮, 崔昊, 王楠, 等. 我国猪肉及种猪进出口现状分析[J]. 中国猪业, 2022, 17(3): 13-16, 23.
XIE K Z, CUI H, WANG N, et al. Analysis of the current situation of import and export of pork and breeding pigs in China[J]. China Swine Industry, 2022, 17(3): 13-16, 23. (in Chinese)
[3]
SI J L, FENG L L, GAO J Y, et al. Evaluating the association between feed efficiency and the fecal microbiota of early-life Duroc pigs using 16S rRNA sequencing[J]. AMB Express, 2020, 10(1): 115. DOI:10.1186/s13568-020-01050-2
[4]
NKRUMAH J D, BASARAB J A, PRICE M A, et al. Different measures of energetic efficiency and their phenotypic relationships with growth, feed intake, and ultrasound and carcass merit in hybrid cattle[J]. J Anim Sci, 2004, 82(8): 2451-2459. DOI:10.2527/2004.8282451x
[5]
BODDICKER N, GABLER N K, SPURLOCK M E, et al. Effects of ad libitum and restricted feeding on early production performance and body composition of Yorkshire pigs selected for reduced residual feed intake[J]. Animal, 2011, 5(9): 1344-1353. DOI:10.1017/S175173111100036X
[6]
MAURER R R, GREGORY K E. Contributions of ovum cytoplasm, uterine environment and postnatal environment to maternal effects in beef cattle[J]. J Anim Sci, 1990, 68(8): 2319-2332. DOI:10.2527/1990.6882319x
[7]
GILBERT H, BILLON Y, BROSSARD L, et al. Review: divergent selection for residual feed intake in the growing pig[J]. Animal, 2017, 11(9): 1427-1439. DOI:10.1017/S175173111600286X
[8]
TURNBAUGH P J, LEY R E, MAHOWALD M A, et al. An obesity-associated gut microbiome with increased capacity for energy harvest[J]. Nature, 2006, 444(7122): 1027-1031. DOI:10.1038/nature05414
[9]
CHANG E B, MARTINEZ-GURYN K. Small intestinal microbiota: the neglected stepchild needed for fat digestion and absorption[J]. Gut Microbes, 2019, 10(2): 235-240. DOI:10.1080/19490976.2018.1502539
[10]
KOH A, DE VADDER F, KOVATCHEVA-DATCHARY P, et al. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites[J]. Cell, 2016, 165(6): 1332-1345. DOI:10.1016/j.cell.2016.05.041
[11]
SINGH K M, SHAH T, DESHPANDE S, et al. High through put 16S rRNA gene-based pyrosequencing analysis of the fecal microbiota of high FCR and low FCR broiler growers[J]. Mol Biol Rep, 2012, 39(12): 10595-10602. DOI:10.1007/s11033-012-1947-7
[12]
TAN Z, YANG T, WANG Y, et al. Metagenomic analysis of cecal microbiome identified microbiota and functional capacities associated with feed efficiency in landrace finishing pigs[J]. Front Microbiol, 2017, 8: 1546. DOI:10.3389/fmicb.2017.01546
[13]
MCCORMACK U M, CURIÃO T, BUZOIANU S G, et al. Exploring a Possible Link between the Intestinal Microbiota and Feed Efficiency in Pigs[J]. Appl Environ Microbiol, 2017, 83(15): e00380-17.
[14]
ROGNES T, FLOURI T, NICHOLS B, et al. VSEARCH: a versatile open source tool for metagenomics[J]. PeerJ, 2016, 4: e2584. DOI:10.7717/peerj.2584
[15]
LI B, ZHANG H B, SHI L L, et al. Saccharomyces boulardii alleviates DSS-induced intestinal barrier dysfunction and inflammation in humanized mice[J]. Food Funct, 2022, 13(1): 102-112. DOI:10.1039/D1FO02752B
[16]
DOUGLAS G M, MAFFEI V J, ZANEVELD J R, et al. PICRUSt2 for prediction of metagenome functions[J]. Nat Biotechnol, 2020, 38(6): 685-688. DOI:10.1038/s41587-020-0548-6
[17]
郑先瑞, 卓明雪, 纪金丽, 等. 皖南黑猪不同生长阶段血清免疫指标及肠道菌群的特征分析[J]. 畜牧兽医学报, 2023, 54(9): 3770-3783.
ZHENG X R, ZHUO M X, JI J L, et al. Characteristics of serum immune indices and intestinal microbiota of Wannan black pigs at different growth stages[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(9): 3770-3783. (in Chinese)
[18]
TAN Z, WANG Y, YANG T, et al. Differences in gut microbiota composition in finishing Landrace pigs with low and high feed conversion ratios[J]. Antonie van Leeuwenhoek, 2018, 111(9): 1673-1685. DOI:10.1007/s10482-018-1057-1
[19]
METZLER-ZEBELI B U, LAWLOR P G, MAGOWAN E, et al. Interactions between metabolically active bacteria and host gene expression at the cecal mucosa in pigs of diverging feed efficiency[J]. J Anim Sci, 2018, 96(6): 2249-2264. DOI:10.1093/jas/sky118
[20]
WEI W, JIANG W B, TIAN Z, et al. Fecal g. Streptococcus and g. Eubacterium_coprostanoligenes_group combined with sphingosine to modulate the serum dyslipidemia in high-fat diet mice[J]. Clin Nutr, 2021, 40(6): 4234-4245. DOI:10.1016/j.clnu.2021.01.031
[21]
LIU Y, WU H, CHEN W B, et al. Rumen microbiome and metabolome of high and low residual feed intake angus heifers[J]. Front Vet Sci, 2022, 9: 812861. DOI:10.3389/fvets.2022.812861
[22]
WANG B, MA M P, DIAO Q Y, et al. Saponin-induced shifts in the rumen microbiome and metabolome of young cattle[J]. Front Microbiol, 2019, 10: 356. DOI:10.3389/fmicb.2019.00356
[23]
XUE M Y, SUN H Z, WU X H, et al. Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance[J]. Microbiome, 2020, 8(1): 64. DOI:10.1186/s40168-020-00819-8
[24]
JIANG X Y, LU N S, ZHAO H C, et al. The microbiome-metabolome response in the colon of piglets under the status of weaning stress[J]. Front Microbiol, 2020, 11: 2055. DOI:10.3389/fmicb.2020.02055
[25]
HUANG P, JIANG A Q, WANG X X, et al. NMN maintains intestinal homeostasis by regulating the gut microbiota[J]. Front Nutr, 2021, 8: 714604. DOI:10.3389/fnut.2021.714604
[26]
CHEN K L A, LIU X J, ZHAO Y C, et al. Long-term administration of conjugated estrogen and bazedoxifene decreased murine fecal β-glucuronidase activity without impacting overall microbiome community[J]. Sci Rep, 2018, 8(1): 8166. DOI:10.1038/s41598-018-26506-1
[27]
BRAHE L K, LE CHATELIER E, PRIFTI E, et al. Specific gut microbiota features and metabolic markers in postmenopausal women with obesity[J]. Nutr Diabetes, 2015, 5(6): e159. DOI:10.1038/nutd.2015.9
[28]
DERRIEN M, VAUGHAN E E, PLUGGE C M, et al. Akkermansia muciniphila gen. nov., sp. nov., a human intestinal mucin-degrading bacterium[J]. Int J Syst Evol Microbiol, 2004, 54(Pt 5): 1469-1476.
[29]
DERRIEN M, BELZER C, DE VOS W M. Akkermansia muciniphila and its role in regulating host functions[J]. Microb Pathog, 2017, 106: 171-181. DOI:10.1016/j.micpath.2016.02.005
[30]
CANDELA M, BIAGI E, SOVERINI M, et al. Modulation of gut microbiota dysbioses in type 2 diabetic patients by macrobiotic Ma-Pi 2 diet[J]. Br J Nutr, 2016, 116(1): 80-93. DOI:10.1017/S0007114516001045
[31]
KHAN T J, AHMED Y M, ZAMZAMI M A, et al. Atorvastatin treatment modulates the gut microbiota of the hypercholesterolemic patients[J]. OMICS A J Integrat Biol, 2018, 22(2): 154-163. DOI:10.1089/omi.2017.0130
[32]
LAHTI L, SALONEN A, KEKKONEN R A, et al. Associations between the human intestinal microbiota, Lactobacillus rhamnosus GG and serum lipids indicated by integrated analysis of high-throughput profiling data[J]. PeerJ, 2013, 1: e32. DOI:10.7717/peerj.32
[33]
STANLEY D, HUGHES R J, GEIER M S, et al. Bacteria within the gastrointestinal tract microbiota correlated with improved growth and feed conversion: challenges presented for the identification of performance enhancing probiotic bacteria[J]. Front Microbiol, 2016, 7: 187.
[34]
LAI L Y, DAI J, TANG H Y, et al. Streptococcus suis serotype 9 strain GZ0565 contains a type Ⅶ secretion system putative substrate EsxA that contributes to bacterial virulence and a vanZ-like gene that confers resistance to teicoplanin and dalbavancin in Streptococcus agalactiae[J]. Vet Microbiol, 2017, 205: 26-33. DOI:10.1016/j.vetmic.2017.04.030
[35]
RADJABZADEH D, BOSCH J A, UITTERLINDEN A G, et al. Gut microbiome-wide association study of depressive symptoms[J]. Nat Commun, 2022, 13(1): 7128. DOI:10.1038/s41467-022-34502-3
[36]
SONG Y Y, CHEN K, LV L, et al. Uncovering the biogeography of the microbial commmunity and its association with nutrient metabolism in the intestinal tract using a pig model[J]. Front Nutr, 2022, 9: 1003763. DOI:10.3389/fnut.2022.1003763
[37]
REYER H, OSTER M, MCCORMACK U M, et al. Host-microbiota interactions in ileum and caecum of pigs divergent in feed efficiency contribute to nutrient utilization[J]. Microorganisms, 2020, 8(4): 563. DOI:10.3390/microorganisms8040563
[38]
INGALA M R, SIMMONS N B, WULTSCH C, et al. Comparing microbiome sampling methods in a wild mammal: fecal and intestinal samples record different signals of host ecology, evolution[J]. Front Microbiol, 2018, 9: 803. DOI:10.3389/fmicb.2018.00803
[39]
胡希怡, 王慧, 李富宽, 等. 肠道菌群调控母幼分离子代行为异常的研究进展[J]. 畜牧兽医学报, 2023, 54(11): 4514-4525.
HU X Y, WANG H, LI F K, et al. Research progress on the regulation of intestinal flora on abnormal behavior of maternal separation in offspring[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(11): 4514-4525. (in Chinese)

(编辑   范子娟)