基于16S rRNA测序探究杜洛克猪饲料利用率与肠道菌群的关系
【摘要】 为探究杜洛克猪饲料利用率与肠道菌群的关系，本研究计算了91头杜洛克公猪30~100kg的校正饲料转化率（FCR），选择FCR数值最高和最低各10头猪分别作为HFCR组和LFCR组，采用微生物16SrRNA基因测序技术分析两组间肠道微生物多样性、组成与结构差异，并对差异微生物功能进行预测。结果显示:两组间肠道微生物在总体Alpha多样性指标上无显著性差异，但在微生物组成上存在多个显著差异的菌属;Fimicutes和Bacteroidetes为杜洛克猪的优势菌门，Prevotella、norank_f_Ruminococcaceae和Lactobacillus是主要的优势菌属;在门水平上，HFCR组WPS-2的丰度显著高于LFCR组;在属水平上，HFCR组个体中unclassified_f_Ruminococcaceae、unclassified_o_Bacteroidales、norank_p_WPS-2和rc4-4的丰度显著高于LFCR组，Collinsella、norank_f_Enterobacteriaceae、Staphylococcus、1-68和norank_o__Streptophyta的丰度显著低于LFCR组。KEGG分析发现，LFCR组和HFCR组肠道菌群基因主要富集在Aminosugarandnucleotidesugarmetabolism、Purinemetabolism、Glycolysis/Gluconeogenesis和Pyrimidinemetabolism等功能通路，但在两组间并无显著差异。本研究结果为利用肠道微生物提高猪的饲料利用率提供了一定的参考和数据基础。关键词:杜洛克;饲料利用率;肠道微生物中图分类号:S828.2文献标识码:ADOI编号:10.19556/j.0258-7033.20200710-04猪在我国畜牧业养殖中占有重要地位，据《中华人民共和国2019年国民经济和社会发展统计公报》显示，2019年我国猪肉产量虽同比下降21.3%，但仍达4255万t。成本核算发现，饲料成本占整个养猪生产成本的60%以上，同时由于猪饲料是玉米-豆粕型饲料，这在一定程度上导致了猪用粮食增加，因此提高饲料利用率是生猪产业发展的必由之路。遗传、疾病、饮食和管理等诸多因素都会影响猪的饲料利用率[1-3]。近年来随着测序技术的发展以及对猪的肠道微生物组的研究，发现肠道菌群是食物消化的重要贡献者，影响宿主的整体生理生长、免疫反应和发病机制，对猪的饲料收稿日期:2020-07-14;修回日期:2020-08-11资助项目:国家现代农业产业技术体系（nycytcgccxtd-15-01）;广西创新驱动发展专项（桂科AA17204029、桂科AA18118051）作者简介:田威龙（1996-），男，安徽宿州人，硕士，主要从事猪遗传育种研究，E-mail:2103755661@qq.com*通讯作者:兰干球（1963-），男，广西来宾人，博士，主要从事猪遗传育种研究，E-mail:gqlan@gxu.edu.cn;梁晶（1989-），男，安徽淮北人，博士，主要从事猪遗传育种研究，E-mail:liangjing@gxu.edu.cn利用率产生重要影响[4-5]。Tan等[6]研究发现，高、低饲料转化率（FCR）长白猪肠道不同位置存在的差异微生物，揭示了与FCR相关的潜在肠道微生物标志物。Mccormack等[7]研究发现，饲料利用效率更高的杂交猪中有益细菌的相对丰度较高，尤其是梭菌和拟杆菌属的成员;而红球菌和红藻菌等潜在有害细菌的相对丰度较低，表明肠道菌群与饲料利用效率密切相关。因此，了解猪肠道菌群的组成和差异可为通过调节肠道微生物的群落、提高猪饲料利用效率提供一定策略。本研究以杜洛克为研究对象，探索高FCR和低FCR猪肠道微生物的组成和差异，为进一步挖掘与猪饲料利用率密切相关的肠道微生物提供理论基础。1材料与方法1.1动物饲养与样品收集试验所用91头美系纯种杜洛克公猪来自广西南宁市某核心育种场。仔猪断奶后，同时转移至相同保育条件下饲养，并在30日龄时转移至同一育肥舍，使用奥斯本自动喂料系统对猪的采食量-146-遗传育种·GeneticsandBreeding2021年第57卷第06期和体重进行测定。对所有试验猪提供相同的玉米-豆粕型基础日粮（NRC，2012），其组成及营养成分见表1。测定期间自由饮水与采食，未使用任何抗生素或药物。当达到（80±1.5）日龄时，采集个体新鲜粪便样品，放入液氮速冻后转至-80℃冰箱保存至使用。表1育肥猪基础饲粮配方及营养成分含量项目原料组成,%玉米麸皮豆粕花生粕磷酸氢钙食盐70%赖氨酸苏氨酸大豆磷脂预混料①石粉营养成分②消化能,MJ/kg粗蛋白质,%粗纤维,%粗灰质,%钙,%总磷,%赖氨酸,%蛋氨酸,%苏氨酸,%74.505.5011.504.001.000.350.700.150.31.001.0013.6813.358.0010.230.450.530.530.270.65注:①每千克全价饲料提供:铜15mg，铁100mg，锰85mg，锌80mg，碘0.5mg，硒0.3mg，维生素A8000IU，维生素D31700IU，维生素E13IU，维生素K1.5mg，维生素B11.5mg，维生素B25mg，维生素B62.5mg，维生素B120.01mg，烟酰胺23mg，叶酸1mg，泛酸9mg，生物素0.10mg。②除代谢能、赖氨酸、蛋氨酸、苏氨酸为计算值，其余均为为实测值。1.2饲料利用率计算和分组利用从奥斯本自动喂料系统中导出的测定期间所有猪只的采食量数据和在初测和结测时称量个体体重数据，依据《全国种猪遗传评估方案(试行)》校正公式计算出校正30~100kgFCR并对FCR进行排序，将FCR数值最低的10只猪作为LFCR组，将FCR数值最高的10只猪作为HFCR组。1.3粪便DNA的提取和测序参照FastDNASPINsoilkit（MPBiomedicals，SantaAna，CA，USA）试剂盒说明书提取粪便基因组DNA。提取的DNA使用微量紫外分光光度计（NanoDrop2000，美国）测量DNA浓度，使用1%琼脂糖凝胶电泳检测DNA质量。扩增细菌16SrRNA基V3-V4区，引物为通用引物（F:338F，5'-ACTCCTACGGGAGGCAGCA-3';806R，5'-GGACTACHVGGGTWTCTAAT-3'）。扩增条件:98℃初始变性1min，98℃变性10s，50℃退火30s，72℃延伸30s，30个循环数，72℃延伸5min。扩增产物经凝胶电泳检测后送上海美吉生物医药科技有限责任公司使用IlluminaMiseq平台测序。1.4生物信息学分析和统计利用QIIME（QuantitativeInsightsIntoMicrobialEcology）（version1.9.1）对16SrRNA原始数据进行质控，使用USEARCH软件软件（version7.0）按相似度≥97%的序列配给相同的操作分类单位（OTU）;RDPClassifier贝叶斯算法（version2.11）对OTUs序列分类注释;使用R语言对粪便菌群门和属水平的PCoA进行分析，并由STAMP软件（version2.1.3）可视化;用Mothur软件（version1.30.2）计算Alpha指数，包括Sobs、Shannon、ACE、Chao1;利用PICRUSt对OTUs丰度标准化，并将每个标准化的OTU对应到相应的Greengenes，获得相应的KEGG信息，并用STAMP对丰度化的KEGG进行差异分析。2结果2.1样本饲料利用率统计如表2所示，LFCR组的FCR为1.661~1.919，平均为1.830，说明其饲料利用效率高;HFCR组的FCR为2.397~2.621，平均为2.469。HFCR组和LFCR组猪的校正30~100kgFCR差异极显著。表2HFCR组和LFCR组猪的FCR表型指标LFCR组HFCR组P值30~100kgFCR1.830±0.2592.469±0.259<0.0012.2样本序列信息HFCR组获得441444个16SrRNA基因序列，LFCR组获得483951个16SrRNA基因序列，共925395个。基于每个样品的OTU指数得到所有样品的稀释性曲线（图1），Shannon指数和Sobs指数曲线上升后均逐步趋于平稳，说明本试验样品测序深度能覆盖肠道菌群的大部分菌种，结果可用于后续数据分析。2.3Alpha多样性分析本研究以细菌物种的多样性指数（Shannon、Simpson）、丰富度指数（Chao1、ACE）、测序覆盖率指数（Coverage）作为猪肠道微生物Alpha多样性的参数，结果显示，HFCR和LFCR组间Alpha多样性指数均无显著性差异（表3）。-147-2021年第57卷第06期ScienceandTechnology·科学技术HFCRLFCRHFCRLFCR5.554.543.532.521.510.50levelUYOnoxedninonnahSABlevelUYOnoxednisboS8007006005004003002001000020004000600080001000012000140001600018000200002200024000采样数量NumberofReadsSampled020004000600080001000012000140001600018000200002200024000采样数量NumberofReadsSampled图中一条曲线代表一个样品，横坐标为从样品中选取的序列数，纵坐标为从样品中观察到的OTU指数。图1Shannon指数曲线（A）Sobs指数曲线（B）分析表3杜洛克猪不同样本的Alpha多样性估计参数HFCR组LFCR组Shannon指数4.551±0.452244.3001±0.22806Simpson指数0.043494±0.0406190.042378±0.9347ACE指数Chao指数838.37±115.76774.81±63.822855.57±113.55789.28±66.775P值0.13470.93470.14580.129Coverage指数0.99512±0.00121360.99607±0.00089270.061032.4基于门和属水平的微生物群落结构组成在门水平上，共鉴定到14个门水平菌，其中HFCR和LFCR组共有门有13个（92.86%），LFCR组有1个独有菌门为p__Fusobacteria（图2-A）。2组均检测到2个优势门（图2-C），分别为Fimicutes（HFCR组占69.57%，LFCR组占73.01%）和Bacteroidetes（HFCR组占26.83%，LFCR组占22.51%），均占据了各组总数的95%以上。在属水平上（图2-B），共鉴定到114个属水平菌，其中HFCR和LFCR组共有的属有110个，HFCR组有1个特有的属为g_unclassified_f_Bacteroidaceae，LFCR组有3个特有的属为g_Fusobacterium、g_1-68、g_unclassified_f_Planococcaceae。2组均检测到4个优势属，分别为Prevotella（HFCR组占15.22%，LFCR组占14.20%）、norank_f_Ruminococcaceae（HFCR组占11.13%，LFCR组占9.82%）、Lactobacillus（HFCR组占10.97%，LFCR组占9.90%）和norank_f_Lachnospiraceae（HFCR组占7.42%，LFCR组占4.54%）（图2-D）。AHFCRLFCRHFCR0131LFCRHFCRLFCRHFCR11103LFCRHFCRLFCRGroupBCDlevelmuiyhPnoecnadnubaytinummocfotbacreP10.80.60.40.2010.8levelsuneGno0.6ecnadnuba0.4ytinummocfo0.2tbacreP0HFCRLFCRGroupPrevotellanorank_f__RuminococcaceaeLactobacillusnorank_f__Lachnospiraceaenorank_f__S24-7norank_f__ClostridiaceaeStreptococcusOscillospiraMegasphaeranorank_o__ClostridialesSMB53RuminococcusRoseburiaMegamonasCoprococcusBlautianorank_o__Bacteroidales[Prevotella]unclassified_o__ClostridialesFaecalibacteriumDoreanorank_f__Christensenellaceaeunclassified_f__Ruminococcaceaenorank_f__CoriobacteriaceaeTreponemaothersVenn图展示的HFCR和LFCR组共有的门（A）和属（B），条形图展示的是HFCR和LFCR组在门（C）和属（D）水平的物种组成。图2HFCR和LFCR组的微生物群落结构组成-148-遗传育种·GeneticsandBreeding2021年第57卷第06期2.5基于门和属水平的微生物组成差异利用主坐标（Principalco-ordinatesanalysis，PCoA）比较HFCR和LFCR组微生物在门和属水平上的组成，从PCoA结果显示，在门和属水平上微生物组成均未依据FCR表型明显分为2组（图3）。采用LEfSe分析比较HFCR和LFCR组之间菌群差异性，分支图显示了HFCR和LFCR组中在微生物的门、纲、目、科、属分类水平上最丰富的微生物群体（图4-A），共发现23个差异项，其中14个在LFCR组中富集，9个在HFCR组中富集，LFCR组中的差异项多于HFCR组。LDA分析（图4-B）结果显示，14个在LFCR组中富集的菌群分别为p_Actinobacteria、g_Collinsella、f_Enterobacteriaceae、o_Enterobacteriales、g_norank_f_Enterobacteriaceae、g_1-68、g_staphylococcus、g_Staphylococcus、g_norank_o_Streptophyta、c_Chloroplast、o_Streptophyta、f_norank_Streptophyta、f_tissierellaceae、o_bacillales、Streptococcus、Christensenellaceae、Ruminococcaceae，9个在HFCR组中富集的菌群分别为g_unclassified_o_Bacteroidales、f_unclassified_o_Bacteroidales、g_rc4-4、p_WPS-2、g_norank_p_WPS-2、o_norank_p_WPS-2、f_norank_p_WPS-2、c_norank_p_WPS-2、g_unclassified_f_Ruminococcaceae。APCoAonPhylumlevelD16HFCRLFCRD30D6D61D68D37D54D67D62D75D72D60D81D84D12D20D56D190.08）%76.3（2CP0.060.040.020-0.02-0.04-0.06-0.08D79D63-0.4-0.35-0.3-0.25-0.2-0.15-0.1-0.0500.050.10.150.20.25PC1（88.07%）PCoAonGenuslevelB）%99.412（CP0.30.250.20.150.10.050-0.05-0.1-0.15-0.2-0.25D61D6HFCRLFCRD63D79D54D19D60D30D12D68D16D84D37D67D56D62D81D75D72D20-0.6-0.5-0.4-0.3-0.2-0.100.10.20.3图3门水平（A）和属水平（B）的PCoA分析PC1（34.97%）在门、属水平比较HFCR和LFCR组前150种微生物类群存在显著差异的微生物，结果发现，在门水平（图5-A），WPS-2在HFCR组中的丰度显著高于LFCR组;在属水平上（图5-B），unclassified_f_Ruminococcaceae、unclassified_o_Bacteroidales、norank_p_WPS-2和rc4-4在HFCR组中的丰度显著高于LFCR组，Collinsella、norank_f_Enterobacteriaceae、Staphylococcus、1-68和norank_o_Streptophyta在LFCR组中的丰度显著高于HFCR组。2.6杜洛克猪肠道微生物组功能预测分析基于PICRUSt分析平台的16S功能预测将获得的OTU丰度表与Greengene数据库比对，获得相应的COG及KO功能信息及其丰度。结果发现，杜洛克猪肠道微生物组功能主要集中于Aminoacidtransportandmetabolism、Carbohydratetransportandmetabolism、Transcription、Cellwall/membrane/envelopebiogenesis、Generalfunctionpredictiononly等方面（图6-A、6-B），高低组无明显差异;KEGG预测分析发现，HFCR组和LFCR组无显著的差异，基因主要富集Aminosugarandnucleotidesugarmetabolism、purinemetabolism、Glycolysis/Gluconeogenesis、Pyrimidinemetabolism等功能通路（图6-C）。3讨论由于我国本地猪种基本都是脂肪型猪，具有耐粗饲和抗病能力强等优点，但存在瘦肉率低、生长速度慢和饲料利用率低等问题，因此为了适应现代生猪产业的要求，我国大量引进西方商业猪种，以改良本地品种或者直接进行繁育作为商品猪，其中杜洛克通常作为最重要的终端父本使用。杜洛克在配套系中作为父本主要是其具有生长速度快、瘦肉率高、饲料利用率高等优良的生长性状，这也是杜洛克在育种工作中的主要育种目标。饲料利用率受饮食、遗传等众多因素的影响，而这些因素恰恰也是肠道菌群的重要影响因素。新一代测序技术的出现极大扩展了人们对于肠道微生物的认识，猪的肠道微生物与猪健康和生长效率之间的关系也越来越受到关注。前人发现当仔猪出生后，暴露在各种环境条件下，大量细菌开始在体内定殖，但是其肠道菌群组成会随着时间的推移而不断变化，而不是一成不变[8];尤其是在断奶后，由于饮食结构发生显著变化，肠道菌群的组成会发生很大的改变[9];而肠道微生物群和粪便微生物群-149-2021年第57卷第06期ScienceandTechnology·科学技术ABvHFCRLFCRa:o__Enterobacterialesb:f__Enterobacteriaceaec:g__norank_f__Enterobacteriaceaed:f___Tissierellaceae_e:g__1_68f:g__rc4_4g:o__Bacillalesh:f__Staphylococcaceaei:g__Staphylococcusj:p__Actinobacteriak:g__Collinsellal:c__Chloroplastm:o__Streptophytan:f__norank_o__Streptophytao:g__norank_o__Streptophytap:f__unclassified_o__Bacteroidalesq:g__unclassified_o__Bacteroidalesr:g__unclassified_f__Bacteroidaceaes:p__WPS_2t:c__norank_p__WPS_2u:o__norank_p__WPS_2v:f__norank_p__WPS_2w:g__norank_p__WPS_2g_unclassified_o_Bacteroidalesf_unclassified_o_Bacteroidalesg_rc4_4f_norank_p_WPS_2g_norank_p_WPS_2p_WPS_2o_norank_p_WPS_2c_norank_p_WPS_2g_unclassified_f_BacteroidalesHFCRLFCRp_Actinobocteriag_Collinsellaf_Enterobacteriaceaeo_Enterobacteriaceaeg_norank_f_Enterobacteriaceaeg_1_68g_Staphylococcusf_Staphylococcaceaeg_norank_o_Streptophytao_Streptophytac_Chloroplastf_norank_o_Streptophytaf_Tissierellaceae_o_Bacillales4.03.53.02.52.01.51.00.50.00.51.01.52.02.53.03.54.0A为微生物LEfSe分析分枝图，B为用LDA分值展示LEfSe分析结果。图4HFCR和LFCR组微生物差异分析LDASCOER(log10)之间具有密切的关系，可以通过粪便微生物的群落变化反映肠道菌群的改变[10]。所以本研究利用16SrRNA测序技术探究粪便菌群与杜洛克饲料利用率的关系，为进定联系[7]。本研究中试验猪均饲养在同一育肥舍，饲喂相同的育肥日粮，采取相同的饲养方式，尽可能减少外界环境因素对于肠道菌群的影响。肠道菌群呈现动态发一步利用肠道菌群提高饲料利用率提供参考依据。有研究表明，肠道微生物与猪的生长性状具有一育过程，直至育肥期才趋于稳定，故本研究测定育肥期30~100kg杜洛克FCR，探究此阶段高低FCR个体间的-150-遗传育种·GeneticsandBreeding2021年第57卷第06期AWPS-295%confidenceintervals0.0000.0100.0200.0300.040-5-4-3-2-1012345Proportions(%)Differencebetweenproportions(%)95%confidenceintervalsBCollinsellaunclassified_o_Bacteroidalesnorank_f_Enterobacteriaceaenorank_f_WPS-2unclassified_F_BacteroidalesStaphylococcusnorank_o_Streptophytarc4-41-66HFCRLFCRHFCRLFCR0.01465Pvalue0.020840.081030.025360.014650.034980.025810.03210.0020480.03484Pvalue0.000.100.200.300.400.50Proportions(%)-5-4-3-2-1012345Differencebetweenproportions(%)图5在门水平（A）和属水平（B）的差异菌群微生物多样性差异对猪的生长性能的影响具有合理性。本研究发现在门水平上，HFCR组和LFCR组中共鉴定到14个门水平菌，其中，2组中检测到2个优势菌门，分别为Fimicutes和Bacteroidetes，这与之前的报道结果类似[11];在属水平，2组共鉴定到114个属水平菌，其中Prevotella、norank_f_Ruminococcaceae、Lactobacillus是2组中占据前三位的优势菌属。Pajarillo等[12]对杜洛克进行16SrRNA高通量测序，发现普氏菌属（Prevotella）一直是其优势菌属，与本研究结果一致。研究发现普氏菌属有利于降解黏蛋白和植物性碳水化合物（如半纤维素和木聚糖），这在饲料消化中具有重要的作用[12]。本研究的优势菌属乳杆菌属（Lactobacillus）多被认为是益生菌。黄金秀等[14]研究发现，作为肠道内优势的乳杆菌之一的罗伊氏乳杆菌，可以通过重塑猪肠道菌群结构，进而影响肠道黏膜免疫功能，调节肠道生的门、纲、目、科、属分类水平上共检测到23个差异项，可作为潜在的生物标志物，其中g_unclassified_o_Bacteroidales等9个在HFCR组中富集，p_Actinobacteria等14个菌群在LFCR组中富集。前人报道，Actinobacteria与脂肪呈正相关，暗示着Actinobacteria的富集可能有利于脂肪的沉积[16];Mccormack等[7]研究也发现，高饲料效率猪回肠中Actinobacteria丰度是低饲料利用效率猪丰度的3倍，说明Actinobacteria可能是潜在标志饲料利用率的菌群之一，同时也证明了粪便微生物组可在一定程度上代表肠道微生物组。属水平进一步分析发现，LFCR组中的差异项显著多于HFCR组且Collinsella、norank_f_Enterobacteriaceae、Staphylococcus、1-68、norank_o_Streptophyta在LFCR组中的丰度显著高于HFCR组，与相关的报道高低组的差异菌群不符[6]，这可能是由性别、年龄、品种、饮食、态平衡。乳酸杆菌属可在均质发酵或异源发酵中将碳水环境、采样方式和保存方式等诸多因素引起的。本研究化合物转化为乳酸，或在异质发酵中将碳水化合物转化为乙酸，这2种酸产生的酸化环境可以抑制其他微生物的生长[15]，有利于维持胃肠道健康以及营养物质吸收，对于猪的生长性能也具有一定意义，也在一定程度上说明了动物的肠道菌群可以影响宿主的生长性状。从本研究PCOA的结果可以看出，在门和属水平上微生物组成都没有根据FCR分为2组，说明高低FCR杜洛克粪便微生物在菌群结构上具有一定的相似性，但是两组间仍表现出具体的差异性，这种差异性很大程度上是由于2组个体在饲料利用率性状上的差别造成的。LEfSe分析显示，HFCR组和LFCR组中在微生物是由粪便微生物代替肠道微生物，结果不同也说明粪便微生物虽然可以在一定程度上代表肠道微生物，但仍有一定差异，虽然已有研究表明肠道微生物与猪的FCR有潜在的联系，但相关菌群的功能并未研究的很透彻。此外，因为宿主本身FCR不同导致肠道微生物组成的不同，还是因为微生物组成不同才导致FCR不同，仍需进一步深入研究。功能预测发现杜洛克高低FCR组的肠道微生物组功能无显著差异，但主要集中于氨基酸运输和代谢、碳水化合物运输和代谢。小肠中的非消化性碳水化合物包括纤维素、木聚糖和抗性淀粉，可以被肠道中的微生物-151-2021年第57卷第06期ScienceandTechnology·科学技术15000001250000100000075000050000025000001500000125000010000007500005000002500000AecnadnubABecnadnubACCOGfunctionclassification:HFCRA:RNAprocessingandmodificationB:ChromatinstructureanddynamicsC:EnergyproductionandconversionD:Cellcyclecontrol,celldivision,chromosomepartitioningE:AminoacidtransportandmetabolismF:NucleotidetransportandmetabolismG:CarbohydratetransportandmetabolismH:CoenzymetransportandmetabolismI:LipidtransportandmetabolismJ:Translation,ribosomalstructureandbiogenesisK:TranscriptionL:Replication,recombinationandrepairM:Cellwall/membrane/envelopebiogenesisN:CellmotilityO:Posttranslationalmodification,proteinturnover,chaperonesP:InorganiciontransportandmetabolismQ:Secondarymetabolitesbiosynthesis,transportandcatabolismR:GeneralfunctionpredictiononlyS:FunctionunknownT:SignaltransductionmechanismsU:Intracellulartrafficking,secretion,andvesiculartransportV:DefensemechanismsW:ExtracellularstructuresZ:CytoskeletonABCDEFGHIJKLMNOPQRSTUVWZFunctionClassCOGfunctionclassification:LFCRABCDEFGHIJKLMNOPQRSTUVWZFunctionClassHeatmapofPathwayLevel3A:RNAprocessingandmodificationB:ChromatinstructureanddynamicsC:EnergyproductionandconversionD:Cellcyclecontrol,celldivision,chromosomepartitioningE:AminoacidtransportandmetabolismF:NucleotidetransportandmetabolismG:CarbohydratetransportandmetabolismH:CoenzymetransportandmetabolismI:LipidtransportandmetabolismJ:Translation,ribosomalstructureandbiogenesisK:TranscriptionL:Replication,recombinationandrepairM:Cellwall/membrane/envelopebiogenesisN:CellmotilityO:Posttranslationalmodification,proteinturnover,chaperonesP:InorganiciontransportandmetabolismQ:Secondarymetabolitesbiosynthesis,transportandcatabolismR:GeneralfunctionpredictiononlyS:FunctionunknownT:SignaltransductionmechanismsU:Intracellulartrafficking,secretion,andvesiculartransportV:DefensemechanismsW:ExtracellularstructuresZ:CytoskeletonBiosynthesisofaminoacidsCarbonmetabolismRibosomeABCtransportersPurinemetabolismPyrimidinemetabolismQuorumsensingTwo-componentsystemStarchandsucrosemetabolismAminosugarandnucleotidesugarmetabolismGlycolysis/GluconeogenesisAminoacyl-tRNAbiosynthesisPyruvatemetabolismCysteineandmethioninemetabolismHomologousrecombinationAlanine,aspartateandglutamatemetabolismOxidativephosphorylationGlycine,serineandthreoninemetabolismPeptidonglycanbilsynthesisCarbonfixationpathwaysinprokaryotes6.36.26.16.05.95.85.75.65.5D30D68D37D54D84D16D75D79D72柱状图展示了HFCR组（A）和LFCR组（B）GOE功能分类结果，热图（C）展示了KEGG三级代谢通路。D19D12D20D61D60D67D63D81D6D56D62图6HFCR和LFCR组猪粪便微生物功能代谢通路分析-152-遗传育种·GeneticsandBreeding2021年第57卷第06期发酵产生挥发性脂肪酸及其衍生物，从而为微生物和宿主生长提供所需的能量[17]，在育成猪配方日粮中主要成分包括富含纤维的玉米和高蛋白大豆，因此猪肠道微生物组的此种构成可能具有更高的纤维素消化利用能力;而发酵膳食纤维在肠道菌群的作用下可产生短链脂肪酸（SCFA），而SCFAs可以减少肠道炎症，从而提高肠道的吸收能力，并提高猪的FCR，推测以氨基酸运输和代谢、碳水化合物运输和代谢为主要功能的肠道菌群是环境和饲料等多重因素导致的，可能与FCR有关系，其机理还需要更深入研究。总之，本研究发现高低FCR组之间的粪便微生物组存在一定差异，支持了肠道微生物群的组成和潜在功能与FCR相联系的假说。但这些差异被认为是微妙的，因为在所有检测到的分类单元中，差异相对较小，并且大多数确实以较低的相对丰度（<2%）出现。不过这些分类单元仍可能会影响FCR，肠道菌群内部的复杂相互作用也可能会对宿主体内稳态和FCR产生最大影响，所以探索猪高、低FCR肠道微生物的组成和差异，挖掘与FCR相关的肠道菌群，进而通过饮食干预猪肠道微生物的组成，提高猪的FCR将是一种非常有潜力的策略和研究思路。4结论本研究采用16SrRNA技术分析杜洛克公猪高、低FCR组的肠道微生物，揭示了与FCR相关的复杂细菌群落。在整个菌群结构上，高低FCR组的优势菌在类别上相同，仅在微生物的相对比例上有所差别。在高和低FCR组中发现了与FCR相关的潜在微生物标记物。功能预测分析发现杜洛克猪肠道微生物主要参与氨基酸运输和代谢、碳水化合物运输和代谢等营养代谢通路。本研究阐述了与猪饲料利用率相关的差异菌群，为进一步研究杜洛克猪饲料利用率与肠道微生物群的关系提供了一定的理论基础。参考文献:[1]WoltmannMD,ClutterAC,BuchananDS,etal.Growthandcarcasscharacteristicsofpigsselectedforfastorslowgaininrelationtofeedintakeandefficiency[J].JAnimSci,1992,70:1049-1059.[2]ArmstrongTA,SpearsJW,CrenshawTD,etal.Boronsuppl-ementationofasemipurifieddietforweanlingpigsimprovesfeedefficiencyandbonestrengthcharacteristicsandaltersplasmalipidmetabolites[J].JNutr,2000,130:2575-2581.[3]SaintilanR,BrossardL,VautierB,etal.Phenotypicandgeneticrelationshipsbetweengrowthandfeedintakecurvesandfeedefficiencyandaminoacidrequirementsinthegrowingpig[J].Animal,2015,9:18-27.[4]VigorsS,SweeneyT,O'SheaCJ,etal.Pigsthataredivergentinfeedefficiency,differinintestinalenzymeandnutrienttransportergeneexpression,nutrientdigestibilityandmicrobialactivity[J].Animal,2016,10:1848-1855.[5]RichardsJD,GongJ,deLangeCFM.Thegastrointestinalmicrobiotaanditsroleinmonogastricnutritionandhealthwithanemphasisonpigs:Currentunderstanding,possiblemodulations,andnewtechnologiesforecologicalstudies[J].CanJAnimSci,2005,85:421-435.[6]TanZ,WangY,YangT,etal.DifferencesingutmicrobiotacompositioninfinishingLandracepigswithlowandhighfeedconversionratios[J].AntonieVanLeeuwenhoek,2018,111:1673-1685.[7]MccormackUM,CuriaoT,BuzoianuSG,etal.Exploringapossiblelinkbetweentheintestinalmicrobiotaandfeedefficiencyinpigs[J].ApplEnvironMicrobiol,2017,83(15):e00380-17.[8]PalmerC,BikEM,DiGiulioDB,etal.Developmentofthehumaninfantintestinalmicrobiota[J].PLoSBiol,2007,5(7):e177.[9]ParksBW,NamE,OrgE,etal.Geneticcontrolofobesityandgutmicrobiotacompositioninresponsetohigh-fat,high-sucrosedietinmice[J].CellMetab,2013,17(1):141-152.[10]StanleyD,GeierMS,ChenH,etal.Comparisonoffecalandcecalmicrobiotasrevealsqualitativesimilaritiesbutquanti-tativedifferences[J].BMCMicrobiol,2015,15(1):51.[11]QuanJ,WuZ,YeY,etal.Metagenomiccharacterizationofintestinalregionsinpigswithcontrastingfeedefficiency[J].FrontMicrobiol,2020,11:32.[12]PajarilloEA,ChaeJP,BalolongMP,etal.Characterizationofthefecalmicrobialcommunitiesofdurocpigsusing16SrRNAgenepyrosequencing[J].Asian-AustralasJAnimSci,2015,28(4):584-591.[13]LamendellaR,DomingoJW,GhoshS,etal.Comparativefecalmetagenomicsunveilsuniquefunctionalcapacityoftheswinegut[J].BMCMicrobiol,2011,11(1):103.[14]黄金秀,王琪,肖融,等.罗伊氏乳杆菌调节仔猪肠道黏膜免疫功能的研究进展[J/OL].动物营养学报,2020,32(6):2454-2459.[15]SamiM,YamashitaH,HironoT,etal.Hop-resistantLactob-acillusbreviscontainsanovelplasmidharboringamultidrugresistance-likegene[J].FermentBioeng,1997,84:1-6.[16]WuGD,ChenJ,HoffmannC,etal.Linkinglong-termdietarypatternswithgutmicrobialenterotypes[J].Science,2011,334(6052):105-108.[17]SinghKM,ShahT,DeshpandeS,etal.Highthroughput16SrRNAgene-basedpyrosequencinganalysisofthefecalmicrobiotaofhighFCRandlowFCRbroilergrowers[J].MolBiolRep,2012,39(12):10595-10602.（责任编辑:周会会）-153-2021年第57卷第06期ScienceandTechnology·科学技术
