CN109856384B - 一种应用代谢组学技术评价酵母培养物有效物组的方法 - Google Patents

一种应用代谢组学技术评价酵母培养物有效物组的方法 Download PDF

Info

Publication number
CN109856384B
CN109856384B CN201811470093.9A CN201811470093A CN109856384B CN 109856384 B CN109856384 B CN 109856384B CN 201811470093 A CN201811470093 A CN 201811470093A CN 109856384 B CN109856384 B CN 109856384B
Authority
CN
China
Prior art keywords
group
weight
effective substance
yeast culture
index
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811470093.9A
Other languages
English (en)
Other versions
CN109856384A (zh
Inventor
甄玉国
孙喆
王涛
张学峰
赵巍
陈雪
赵小丽
郑艳秋
张玲玲
张维刚
谭胜男
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Borui Technology Co ltd
Original Assignee
Changchun Borui Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changchun Borui Technology Co ltd filed Critical Changchun Borui Technology Co ltd
Priority to CN201811470093.9A priority Critical patent/CN109856384B/zh
Publication of CN109856384A publication Critical patent/CN109856384A/zh
Application granted granted Critical
Publication of CN109856384B publication Critical patent/CN109856384B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明涉及微生物有效物组评价领域,具体地涉及一种应用代谢组学技术评价酵母培养物有效物组的方法,利用代谢组学的检测技术对酵母培养物有效物组进行解析,并将多变量降维的数学思路运用到挖掘有效物组的组成中,参照筛选解析的有效物组成分,采用化学单品模拟有效物组的构成关系进行配伍组合,再通过动物饲养试验对有效物组的功效进行验证。得出以下结论:酵母培养物有效物组作为日粮添加剂能显著提高肉仔鸡生长性能、免疫功能、以及抗氧化能力和屠宰性能。

Description

一种应用代谢组学技术评价酵母培养物有效物组的方法
技术领域
本发明涉及微生物有效物组评价领域,特别是涉及一种应用代谢组学技术评价酵母培养物有效物组的方法。
背景技术
代谢组学是系统生物学的一门新兴科学,属于是系统生物学的重要组成部分,是构成系统生物学的基因组学、蛋白质组及转录组学之后发展起来的一门科学。代谢组学的概念的提出,最早是由1971年Devaux等提出的代谢谱分析。他们通过气相色谱质谱联用技术(GC-MS)对试验者体液的代谢物进行定性、定量测定,此方法主要用来对某些疾病进行筛选和诊断,至今仍在应用。
代谢组学是利用现代生物分析手段方法,来研究生物体或细胞在受到外界因素刺激或扰动之后,处于特定生理时期内的,并且与代谢途径相关的小分子代谢物(分子质量<1000)(如氨基酸、脂肪酸、糖类、维生素和脂肪等)随时空的动态变化规律,来探明生物体代谢途径的一种技术。它分析对象主要是血浆或血清、唾液、尿液以及细胞和组织的提取液等。
自20世纪90年代中期,代谢组学快速发展起来,被广泛地应用在疾病诊断药物开发等方面。有关植物代谢组学、微生物代谢组学等领域的报道很多。而代谢组学技术在动物营养和饲料领域的应用还处于起步阶段。营养需要量和饲料营养价值评定一直是动物营养研究的两大主题。而代谢组学研究技术日益成熟,扩展了动物营养学研究的范围,使其不仅从宏观扩展还是到微观深入都提供了强大的技术支撑。
酵母培养物作为一种极具替代抗生素潜力的产品,具有提高动物生产性能,调节免疫功能、改善胃肠道微生态平衡,保障肠道健康等多种营养保健的功效,在畜牧生产中具有广泛的应用前景。而酵母培养物是特定的培养基经酵母菌发酵技术制得的一种营养丰富、成分复杂微生态发酵产品。多年研究普遍认为,其营养保健作用是通过其发酵代谢产物来实现的,然而酵母培养物所包含的代谢产物种类繁多,除少数已知的营养物质外,多数是以各种化合物形式存在的安全性未知成分,其主要的或特定的有效作用成分未有确切描述和解释,因此酵母培养物成分中发挥功能的主要活性成分有及其之间的协同作用,是需要进一步研究确证的问题。
在饲料研发领域,利用代谢组学可以对饲料品质等方面作评价,目前应用最多的是评估饲用微生物在肠道定植和繁殖成优势菌群情况,以便直接研究益生菌的作用机制,有助于益生菌类饲料的开发。另外利用代谢组组学快速取样技术并结合其他分析技术,可以实现快速、高频率定量对细胞内代谢物进行测定,用于监测发酵动态过程。但至今并未有应用代谢组学技术评价酵母培养物有效物组的方法报道。
发明内容
本发明基于以往有关酵母培养物的代谢组学分析方法的建立、多变量统计分析如主成分分析等研究成果研究思想方法、在建立的气相色谱质谱联用的 GC-MS检测方法基础上,针对酵母培养物“有效物组”的探究,开展的研发。从而为实现客观、准确、快捷地评价酵母培养物产品品质,科学认识其功能实质提供了重要的研究参考依据。因为酵母培养物的代谢产物很多,本发明专利的目的就是找出起到关键作用的物质,本发明称之为有效物组,这个组可以是一个物质也可以是多个物质的组合。
本发明提供了一种应用代谢组学技术评价酵母培养物有效物组的方法,得出以下结论:酵母培养物有效物组作为日粮添加剂能显著提高肉仔鸡生长性能、免疫功能、对抗氧化能力和屠宰性能。
解决上述技术问题的技术方案如下:
本发明设计了一种应用代谢组学技术评价酵母培养物有效物组的方法,是利用代谢组学的检测技术对酵母培养物有效物组进行解析,并将多变量降维的数学思路运用到挖掘有效物组的组成中,参照筛选解析的有效物组成分,采用化学单品模拟有效物组的构成关系进行配伍组合,再通过动物饲养试验对有效物组的功效进行验证。
进一步地说,代谢组学的检测技术包括样品采集和处理、代谢物的鉴定和多元统计分析。
样品采集和处理是选用肉鸡为样品,采集和处理方法为肉鸡于21、42日龄进行屠宰实验,每个重复选2只宰杀并解剖,解剖后将一侧盲肠结扎,置于干燥的灭菌试管中,保存于-20℃保存,待用于菌群高通量测序,每组取肉仔鸡盲肠内容物0.5g;
更进一步地说,采用化学单品模拟有效物组的构成关系进行配伍组合具体表现为回归模型的建立与应用,所述的模型建立对生产性能中的日增重指标与筛选到的酵母培养物有效物组之间关系进行研究,利用多元统计分析方法,建立生产性能与酵母培养物有效物组的函数关系,从而进一步揭示有效物组的组合配伍与生产性能的关系。
AA肉仔鸡具有生长快,饲养周期短及饲料转化率高等优点,被广泛应用在肉鸡饲养生产中。AA肉仔鸡的生产性能的高低是直接影响到生产水平和经济效益的一个重要指标。而饲料的配方、鸡舍的温度、光照等因素影响着肉仔鸡的生产性能。
通过肉仔鸡饲养试验,筛选出生产性能最佳的发酵产品为发酵24h的酵母培养物,生产性能较差的是发酵60h的酵母培养物,为了找到造成这种差异的物质,因此在这两组中寻找潜在生物标志物即有效物组的组成。OPLS-DA使用 SIMCA-P-13.0(Umetrics AB,
Figure BDA0001890708420000032
Sweden)软件完成,模型的拟合度及解释率都较高,其中R2X(cum)和R2Y(cum)分别表示模型对X和Y矩阵的解释率,Q2(cum)表示模型的可预测度。R2X(cum)为0.963,R2Y(cum)为0.998, Q2(cum)为0.653,R2和Q2值通常高于0.4说明模型尚可,越接近越说明模型越稳定可靠。(具体操作见实施例部分详细说明)
1)模型建立定义
取各试验组且不同饲养时期的日增重数据,将日增重设定为Y;将酵母培养物有效物组筛选得到的其组成即生物标志物甘氨酸、果糖、肌醇、吡喃(型) 半乳糖、蔗糖相对应地分别标注为X1、X2、X3、X4、X5。X表示各组中某物质的采食量,Xn为Xn峰面积/Xn所在产品组别总峰面积*采食量。利用SPSS进行回归分析,建立了日增重与有效物组中各物质的回归方程Y=f(x)。
2)回归方程建立
第一步在SPSS中输入相关数据,通过酵母培养物有效物组筛选得到的其组成成分X1、X2、X3、X4、X5及相对应的日增重y。计算数据的基本统计量,结果如下表所示:
Figure BDA0001890708420000031
第二步是计算X1、X2、X3、X4、X5与y的偏相关性,结果如下表所示,经 spss计算,从表中结果可以看出X1、X2、X3、X4、X5与y之间均存在较强相关性,存在组建回归方程的可能。
Figure BDA0001890708420000041
第三步,根据数据绘制散点图;在相关性分析的基础之上,根据散点图特点,结合曲线估计可确认各成分与生产性能之间关系大致为立方函数。通过非线性回归分析得出回归方程基本结构:
Y=m+n*x1*x2*x3*x4*x5+a1*x1*x1*x1+a2*x1*x1+a3*x1+b1*x2*x2*x2+b2*x 2*x2+b3*x2+c1*x3*x3*x3+c2*x3*x3+c3*x3+d1*x4*x4*x4+d2*x4*x4+d3*x4+e1* x5*x5*x5+e2*x5*x5+e3*x5
指定各参数系数的初始数值。经spss计算,得到回归方程中各参数系数(即 X1、X2、X3、X4、X5与y在回归方程中的系数)结果如下:
Figure BDA0001890708420000042
Figure BDA0001890708420000051
由结果可知,R2=0.79即回归方程拟合度约为80%,拟合度较为理想,此回归函数可以信任。代入各参数系数,得到回归方程如下:
Y=32.994-1.458*x1*x2*x3*x4*x5-0.612*x1*x1*x1+7.506*x1*x1-22.364*x 1-0.819*x2*x2*x2-2.756*x2*x2+24.360*x2+1.848*x3*x3*x3-13.526*x3*x3+1 9.291*x3+7.630*x4*x4*x4-32.987*x4*x4+38.317*x4-0.134*x5*x5*x5+2.059* x5*x5-6.071*x5。
为从数学统计的角度验证选取的五种成分是否为决定生产性能的主要因素,根据已获得的X1、X2、X3、X4、X5与y之间的回归方程,进行如下计算:
(1)选取y值。根据实验数据,选取在现有条件下,添加全部成分时所能获取的最佳生产性能值。把该值作为下一步计算的y值。
(2)把上一步得到的y值带入已知的X1、X2、X3、X4、X5与y之间回归方程,同时根据现有试验条件指定X1、X2、X3、X4、X5的取值范围。编写vb 计算程序,计算在指定范围内是否存在能通过回归方程得到指定y值的X1、X2、 X3、X4、X5。经计算得到符合要求的X1、X2、X3、X4、X5。
(3)计算上一步得到的X1、X2、X3、X4、X5间的比值。在实践中按得到的比值调配X1、X2、X3、X4、X5五种成分的比例,经验证得到的生产性能值和指定的y值不存在显著性差异,继而证明,只选取五种成分同样可以达到选取全部成分所能达到的生产性能值。
(4)本数学模型仅适用于在正常肉仔鸡实践生产范围内,有一定的局限性。
该方法具体包括以下步骤:
(1)采用单因素完全随机分组设计,选取336羽1日龄健康体重相近的 AA(爱拔益加)肉仔鸡,公母各半随机分为7组,每组6个重复,每个重复8 只鸡;试验期为42d;
酵母培养物有效物组筛选分组:对照组饲喂基础日粮,在基础日粮中添加5 个时间差异化培养物的酿酒酵母培养物,各处理组添加量为0.192%(比例为发酵液干物质质量/饲料质量),(即发酵时间分别为YC-12h、YC-24h、YC-36h、 YC-48h、YC-60h),以及达农威益康VXP产品。
DZ为对照组饲喂基础日粮,ABC组为酵母培养物有效物组,即分别饲喂在基础日粮基础上添加有效物组3个配伍组合,D组饲喂发酵24h的酵母培养物 (YC24h),E、F组为XP组,在基础日粮基础上分别添加0.1%、1%“XP”饲料添加剂;整个试验周期饲粮中不添加抗生素与其他促生长剂;
所述的有效物组为甘氨酸、果糖、肌醇、吡喃(型)半乳糖、蔗糖;试验日粮分配表如下表所示:
Figure BDA0001890708420000061
ABC组在YC中所占比例配伍情况如下表所示:
Figure BDA0001890708420000062
各物质配伍说明:取生产性能为发酵24hYC产品的日增重,得到相应的X 值,即配伍1(有效物组A);在此基础上调整各物质的比例,所述的比例为各物质±10%,配置各物质的含量,为配伍2和3(有效物组B、C组);
(2)血清样品采集:在21日龄和42日龄时,随机选择2只健康和体重相近的鸡从颈静脉采血,采血前12h对试验鸡停料不停水;采血收集完成后,于室温放置待血清析出后将其分装于离心管中,3000r/min离心10min,取其上清分装于Eppendorf试管中,并立即浸入液氮,置-80℃低温保存,用于生化分析用;
器官组织样品采集:在21日龄和42日龄时,随机选择2只健康和体重相近的试验鸡,空腹12h后进行屠宰。取胸腺、脾脏、肌胃、腺胃和法氏囊用滤纸吸去血渍,剔除脂肪后称鲜重并记录,样品保存在-80℃待用;
肌肉样品的采集:于试验第42日结束时每个重复随机选取与该重复平均体重相近的肉鸡2只,颈静脉放血致死,取胸肌和腿肌样品(去除脂肪和结缔组织),称重并记录;
(3)生产性能:试验在21和42日龄时,以重复为单位,分别于试验开始的1、7、14、21、28、35、42d早晨,记录鸡初始重和末重,并记录整个试验期供料量、剩余料量和损失料量,计算平均日采食量(ADFI)、平均日增重(ADG) 和料重比(F/G);
免疫器官指数:在试验的21d和42d,从每个重复中随机选取2只鸡,空腹称体量,屠宰后,分离胸腺、法氏囊和脾脏等免疫器官并称重、记录,器官指数=器官质量(mg)/活体质量(g);
血清生化指标:使用全自动血液生化分析仪测定血中白蛋白,球蛋白,血清尿素氮与碱性磷酸酶;
血清免疫指标:在第21、42日龄进行屠宰,采血管收集血液,3500r/m, 4℃离心10min,取血清于1mL离心管,至于-20℃保存,用于血清总蛋白(TP)、白蛋白(ALB)、尿素氮(UREA)、葡萄糖(GLU)等测定,检测使用全自动生化分析仪;分泌性免疫球蛋白A(IgA)、分泌性免疫球蛋白G(IgG)用ELISA 试剂盒(南京建成生物工程研究所)测定;
血清抗氧化指标:血清中超氧化物歧化酶(SOD)活性、谷胱甘肽过氧化物酶(GPX)活性和丙二醛(MDA)含量使用南京建成生物工程研究所生产的试剂盒测定;
屠宰性能测定:在试验42d时,每个重复各随机取2只鸡,进行屠宰测定;活体重:禁食12h,宰前称重;屠体重:屠宰后,颈静脉放血并去毛后再次称得的重量;全净膛重为屠体重去除气管、食道、嗦囊、肠、脾脏、胰腺、胆囊和生殖器官心、肝、肾、腺胃、肌胃、腹脂、头、脚的后的重量。分割左侧胸肌、腿肌并称重,分别为胸肉重和腿肉重;腹脂重主要是腹部板油及肌胃周围的脂肪的重量。计算公式:腹脂率=腹脂重/全净膛重×100%;胸肌率=胸肌重/全膛净重×100%腿肌率=腿肌重/全膛净重×100%;全净膛率=全膛净重/活重×100%;
(4)DNA提取
称取21日龄肉鸡的盲肠内容物0.3g于无菌PBS中,涡轮震荡30S后,10000 rpm转速下离心10min后取上清溶于DNA缓冲液中;肠道内容物DNA提取方法参照试剂盒QIAamp FastDNA Stool Mini Kit操作步骤(QIAGEN,German);
16S r RNA基因V3~V4区的PCR扩增
将样品DNA稀释至20ng/μL作为模板,扩增细菌16S rRNA V3-4区。使用常规引物338F(5′-ACT CCTAC G GGAGGCAGCA-3′)和806R(5′-GGA CTAC H VGGGTWTCTAAT-3′)。
PCR扩增体系为25μL:5×reaction buffer 5μL,5×GC buffer 5μL, dNTP(2.5mM)2μL,正向引物(10uM)1μL,反向引物(10uM)1μL,模板DNA 2μL,ddH2O 8.75μL,Q5DNA酶0.25μL。PCR扩增条件:98℃预变性2min;然后98℃变性15s,55℃退火30s,72℃延伸30s,25-30个循环;最后72℃延伸5min,10℃保持待测定。PCR扩增产物通过2%琼脂糖凝胶电泳进行检测,采用AXYGEN公司的凝胶回收试剂盒切胶回收。将PCR扩增回收产物通过荧光试剂为Quant-iT PicoGreen dsDNA Assay Kit进行定量;
(5)数据处理与分析:
生产应用试验数据经Excel 2007初步整理后,采用SPSS 23.0统计软件中One-WayANOVA进行方差分析,多重比较采用Duncan's法;P<0.05为差异显著,P<0.01为差异极显著,数据用“平均数±标准差(Mean±SD)”表示;
运用QIIME软件识别疑问序列,并剔除:5’端引物错配碱基数>1的序列;含有连续相同碱基数>8的序列;再调用USEARCH检查并剔除嵌合体序列;
Alpha多样性指数:用于比较不同样本的多样性,首先对序列量拉平处理,从而校正测序深度引起的多样性差异;再使用QIIME软件分别对每个样本计算侧重于体现群落丰富度Chao1指数和ACE指数,以及兼顾群落均匀度的Shannon 指数和Simpson指数;
进一步地说,步骤(1)中所述的基础日粮的组成由下表所示:
Figure BDA0001890708420000081
Figure BDA0001890708420000091
其中,每千克预混料包含:维生素A 5000IU;维生素D3 1500IU;维生素E 15IU;维生素K3 0.8mg;维生素B12 0.01mg;叶酸0.5mg;烟酸50 mg;泛酸8mg;生物素0.1mg;吡哆醇2.2mg;核黄素4.4mg;硫胺素1.6 mg。
本发明的有益效果是:
本发明设计思路清晰,得出以下结论:酵母培养物有效物组作为日粮添加剂能显著提高肉仔鸡生长性能、免疫功能、对抗氧化能力和屠宰性能。日粮中添加不同处理组对肉仔鸡21日龄及42日龄时脾脏指数、胸腺指数和法氏囊指数的影响差异不显著。在42日龄时,各添加处理组胸腺指数较对照组呈增加趋势,但并不显著;添加酵母培养物有效物组显著提高了21日龄与42日龄肉鸡的体液免疫指标;酵母培养物有效物组血清中各抗氧化指标较对照组有增加趋势,但差异不显著。而在21日龄时有效物组B与YC组D血清中MDA的含量低于对照组与XP组,其中有效物组B血清中的MDA含量最低,比对照组降低了 62.27%。各处理组对42d肉仔鸡屠宰率、全净膛屠宰率、胸肌率、腿肌率无显著影响,但有一定的升高趋势。以97%作为相似性阐值统计得到8000多个OTU,且不同组别OUT数量不相同,酵母培养物及其配伍处理组较对照组数量有所减少,说明酵母培养物的添加降低了盲肠菌群的多样性,但添加处理组间并没有显著性变化。
下面结合附图和具体实施方式对本发明作进一步详细的说明。
附图说明:
图1为酵母培养物对肉仔鸡生长性能的影响;
图2是2个发酵时间点样品的OPLS-DA得分图结果;
图3是基于OPLS-DA得分图模式下的S-Plot载荷图;
图4为散点图;
具体实施方式
实施例1:
准备原料:AA肉仔鸡购自吉林德详公司;酿酒酵母培养物由吉林农业大学吉农博瑞研发中心提供;达农威益康VXP是购于达农威中国有限公司;
N,O-双(三甲基硅烷)三氟乙酰胺(BSTFA含1%TMCS)、三甲基氯硅烷是美国Sigma-Aldrich公司的TMCS;上海阿拉丁生化科技股份有限公司的色谱纯吡啶;
美国TEDIA公司的正己烷、美国Sigma-Aldrich公司的甲氧胺盐酸盐。
美国Agilent公司的Agilent7890A/5975C气相色谱-质谱联用仪,配有 Agilent化学工作站:MSD Chem station(E.02.01.1177,Agilent),自动进样器(Agilent7683),Agilent HP-5MS(30m×0.25mm×0.25μm)石英毛细管柱。德国Hettich公司的MIKRO-22R台式低温冷冻离心机,北京博医康有限公司的FD-1冷冻干燥机。
(1)采用单因素完全随机分组设计,选取336羽1日龄健康体重相近的 AA(爱拔益加)肉仔鸡,公母各半随机分为7组,每组6个重复,每个重复8 只鸡;试验期为42d;
酵母培养物有效物组筛选分组:对照组饲喂基础日粮,在基础日粮中添加5 个时间差异化培养物的酿酒酵母培养物,各处理组添加量为0.192%(比例为发酵液干物质质量/饲料质量),(即发酵时间分别为YC-12h、YC-24h、YC-36h、 YC-48h、YC-60h),以及达农威益康VXP产品,上述材料均由吉林农业大学吉农博瑞奶牛饲料研发中心提供。
酵母培养物有效物组验证分组:DZ为对照组饲喂基础日粮,ABC组为酵母培养物有效物组,即分别饲喂在基础日粮基础上添加有效物组3个配伍组合,D 组饲喂发酵24h的酵母培养物(YC24h),E、F组为XP组,在基础日粮基础上分别添加0.1%、1%“XP”饲料添加剂;整个试验周期饲粮中不添加抗生素与其他促生长剂;
所述的有效物组为甘氨酸、果糖、肌醇、吡喃(型)半乳糖、蔗糖;试验日粮分配表如下表所示:
试验组 DZ组 A组 B组 C组 D组 E组 F组
基础日粮
试验日粮 0对照 配伍1 配伍2 配伍3 24hYC 0.1%XP 1%XP
ABC组在YC中所占比例配伍情况如下表所示:
Figure BDA0001890708420000101
Figure BDA0001890708420000111
各物质配伍说明:取生产性能为发酵24hYC产品的日增重,得到相应的X 值,即配伍1(有效物组A);在此基础上调整各物质的比例,所述的比例为各物质±10%,配置各物质的含量,为配伍2和3(有效物组B、C组);
(2)血清样品采集:在21日龄和42日龄时,随机选择2只健康和体重相近的鸡从颈静脉采血,采血前12h对试验鸡停料不停水;采血收集完成后,于室温放置待血清析出后将其分装于离心管中,3000r/min离心10min,取其上清分装于Eppendorf试管中,并立即浸入液氮,置-80℃低温保存,用于生化分析用;
器官组织样品采集:在21日龄和42日龄时,随机选择2只健康和体重相近的试验鸡,空腹12h后进行屠宰。取胸腺、脾脏、肌胃、腺胃和法氏囊用滤纸吸去血渍,剔除脂肪后称鲜重并记录,样品保存在-80℃待用;
肌肉样品的采集:于试验第42日结束时每个重复随机选取与该重复平均体重相近的肉鸡2只,颈静脉放血致死,取胸肌和腿肌样品(去除脂肪和结缔组织),称重并记录;
(3)生产性能:试验在21和42日龄时,以重复为单位,分别于试验开始的1、7、14、21、28、35、42d早晨,记录鸡初始重和末重,并记录整个试验期供料量、剩余料量和损失料量,计算平均日采食量(ADFI)、平均日增重(ADG) 和料重比(F/G);
免疫器官指数:在试验的21d和42d,从每个重复中随机选取2只鸡,空腹称体量,屠宰后,分离胸腺、法氏囊和脾脏等免疫器官并称重、记录,器官指数=器官质量(mg)/活体质量(g);
血清生化指标:使用全自动血液生化分析仪测定血中白蛋白,球蛋白,血清尿素氮与碱性磷酸酶;
血清免疫指标:在第21、42日龄进行屠宰,采血管收集血液,3500r/m, 4℃离心10min,取血清于1mL离心管,至于-20℃保存,用于血清总蛋白(TP)、白蛋白(ALB)、尿素氮(UREA)、葡萄糖(GLU)等测定,检测使用全自动生化分析仪;分泌性免疫球蛋白A(IgA)、分泌性免疫球蛋白G(IgG)用ELISA 试剂盒(南京建成生物工程研究所)测定;
血清抗氧化指标:血清中超氧化物歧化酶(SOD)活性、谷胱甘肽过氧化物酶(GPX)活性和丙二醛(MDA)含量使用南京建成生物工程研究所生产的试剂盒测定;
屠宰性能测定:在试验42d时,每个重复各随机取2只鸡,进行屠宰测定;活体重:禁食12h,宰前称重;屠体重:屠宰后,颈静脉放血并去毛后再次称得的重量;全净膛重为屠体重去除气管、食道、嗦囊、肠、脾脏、胰腺、胆囊和生殖器官心、肝、肾、腺胃、肌胃、腹脂、头、脚的后的重量。分割左侧胸肌、腿肌并称重,分别为胸肉重和腿肉重;腹脂重主要是腹部板油及肌胃周围的脂肪的重量。计算公式:腹脂率=腹脂重/全净膛重×100%;胸肌率=胸肌重/全膛净重×100%腿肌率=腿肌重/全膛净重×100%;全净膛率=全膛净重/活重×100%;
(4)DNA提取
称取21日龄肉鸡的盲肠内容物0.3gg于无菌PBS中,涡轮震荡30S后,10000 rpm转速下离心10min后取上清溶于DNA缓冲液中;肠道内容物DNA提取方法参照试剂盒QIAampFast DNA Stool Mini Kit操作步躁(QIAGEN,German), 样品16S rRNA基因V3-V4区的PCR扩增由上海派森诺科技股份有限公司完成;
16S r RNA基因V3~V4区的PCR扩增
将样品DNA稀释至20ng/μL作为模板,扩增细菌16S rRNA V3-4区。使用常规引物338F(5′-ACT CCTAC G GGAGGCAGCA-3′)和806R(5′-GGA CTAC H VGGGTWTCTAAT-3′),MiSeq测序平台检测。
PCR扩增体系为25μL:5×reaction buffer 5μL,5×GC buffer 5μL, dNTP(2.5mM)2μL,正向引物(10uM)1μL,反向引物(10uM)1μL,模板DNA 2μL,ddH2O 8.75μL,Q5DNA酶0.25μL。PCR扩增条件:98℃预变性2min;然后98℃变性15s,55℃退火30s,72℃延伸30s,25-30个循环;最后72℃延伸5min,10℃保持待测定。PCR扩增产物通过2%琼脂糖凝胶电泳进行检测,采用AXYGEN公司的凝胶回收试剂盒切胶回收。将PCR扩增回收产物通过荧光试剂为Quant-iT PicoGreen dsDNA Assay Kit进行定量;
(5)数据处理与分析:
生产应用试验数据经Excel 2007初步整理后,采用SPSS 23.0统计软件中One-WayANOVA进行方差分析,多重比较采用Duncan's法;P<0.05为差异显著,P<0.01为差异极显著,数据用“平均数±标准差(Mean±SD)”表示;
采用Illumina MiSeq平台对群落DNA片段进行双端(Paired-end)测序;首先对FASTQ格式的双端序列逐一作质量筛查,要求截断后的序列长度≥150 bp,且不允许存在模糊碱基;再利用FLASH软件,对通过质量初筛的双端序列进行配对连接;最后,根据样本所对应的Index信息,将连接后的序列识别分配入对应样本,获得样本的有效序列;
运用QIIME软件识别疑问序列,并剔除:5’端引物错配碱基数>1的序列;含有连续相同碱基数>8的序列;再调用USEARCH检查并剔除嵌合体序列;
Alpha多样性指数:用于比较不同样本的多样性,首先对序列量拉平处理,从而校正测序深度引起的多样性差异;再使用QIIME软件分别对每个样本计算侧重于体现群落丰富度Chao1指数和ACE指数,以及兼顾群落均匀度的Shannon 指数和Simpson指数;
进一步地说,步骤(1)中所述的基础日粮的组成由下表所示:
日粮组成 1-21d 22-42d 营养成分 1-21d 22-42d
玉米 52.07 55.57 代谢能(MJ/kg) 12.746 12.719
大豆粕 35.50 34.00 粗蛋白(%) 20.722 19.609
鱼粉 5.00 3.60 钙(%) 0.910 0.865
大豆油 4.00 3.5 总磷(%) 0.676 0.584
磷酸氢钙 1.00 0.70 有效磷(%) 0.447 0.356
石粉 1.00 1.20 赖氨酸(%) 1.347 0.459
食盐 0.30 0.30 蛋氨酸(%) 1.248 0.432
赖氨酸 0.04 0.04
蛋氨酸 0.09 0.09
1%预混料 1.00 1.00
合计 100 100
其中,每千克预混料包含:维生素A 5000IU;维生素D3 1500IU;维生素E 15IU;维生素K3 0.8mg;维生素B12 0.01mg;叶酸0.5mg;烟酸50 mg;泛酸8mg;生物素0.1mg;吡哆醇2.2mg;核黄素4.4mg;硫胺素1.6 mg。
由图1可知,在整个试验期7-42d,试验组料重比(F/G)均优于对照组,其中添加酿酒酵母培养物组对肉仔鸡的料重比(F/G)都有不同程度的降低趋势,其中YC24组显著低于对照组。
酵母培养物有效物组的筛选
根据肉仔鸡饲养试验,筛选出生产性能最佳的发酵产品为发酵24h的酵母培养物(YC24组),生产性能较差的是发酵60h的酵母培养物(YC60组),为了找到造成这种差异的物质,因此在这2组中寻找潜在生物标志物即有效物组的组成。OPLS-DA使用SIMCA-P-13.0(Umetrics AB,
Figure BDA0001890708420000131
Sweden)软件完成,模型的拟合度及解释率都较高,其中R2X(cum)和R2Y(cum)分别表示模型对X和Y矩阵的解释率,Q2(cum)表示模型的可预测度。R2X(cum)为0.963, R2Y(cum)为0.998,Q2(cum)为0.653,R2和Q2值通常高于0.4说明模型尚可,越接近与1说明模型越稳定可靠。2个发酵时间点样品的OPLS-DA得分图结果如图2,可以直观看出发酵24h组与发酵60h组有着较好的分离。图3是基于OPLS-DA得分图模式下的S-Plot载荷图。在载荷图中可变量,数据点距离原点越远据点的贡献是正性的,反之则为负性的关性,数据点距离原点越远则为负相关。所以,在S型曲线两端的数据对分别代表了每个样本组中可信度最高的化合物,即特征代谢物。首先与NIST质谱数据库的比对,这2组的差异代谢物质整理见下表。
Figure BDA0001890708420000141
注:上述物质都被包含在2组中。
筛选标志物的标准是VIP值>1且P值<0.05。在差异物质中筛选符合VIP 值>1且P值<0.05的物质。筛选到VIP值大于1的代谢产物后,对B、E组分别行独立样本t检验,在这些代谢产物中以p<0.05为筛选标准,最终得出酵母培养物的生物标志物即有效物组的组成。特征代谢物见下表。由下表可知,筛选到10种生物标志物,酵母培养物的生物标志物即有效物组的组成主要为氨基酸和糖类。其中甘氨酸、果糖、肌醇与酵母培养物产品达农威VXP报道的主要成一致。
Figure BDA0001890708420000151
生物标志物的确定:为了确证这些物质的存在,我们进一步定量定性表中物质,最终筛选并确定以下5种物质为酵母培养物的生物标志物即有效物组成分见下表。
Figure BDA0001890708420000152
回归模型的建立与应用
AA肉仔鸡具有生长快,饲养周期短及饲料转化率高等优点,被广泛应用在肉鸡饲养生产中。AA肉仔鸡的生产性能的高低是直接影响到生产水平和经济效益的一个重要指标。而饲料的配方、鸡舍的温度、光照等因素影响着肉仔鸡的生产性能。本模型建立对生产性能中的日增重指标与筛选到的酵母培养物有效物组之间关系进行研究,利用多元统计分析方法,建立生产性能与酵母培养物有效物组的函数关系,从而进一步揭示有效物组的组合配伍与生产性能的关系,为酵母培养物有效物组的确证提供参考依据。
模型建立定义
在试验中,取本发明各试验组且不同饲养时期的日增重数据,将日增重设定为Y;将酵母培养物有效物组筛选得到的其组成即生物标志物甘氨酸、果糖、肌醇、吡喃(型)半乳糖、蔗糖相对应地分别标注为X1、X2、X3、X4、X5。X 表示各组中某物质的采食量,Xn为Xn峰面积/Xn所在产品组别总峰面积*采食量。利用SPSS进行回归分析,建立了日增重与有效物组中各物质的回归方程Y=f (x)。
回归方程建立
第一步在SPSS中输入相关数据,通过酵母培养物有效物组筛选得到的其组成成分X1、X2、X3、X4、X5及相对应的日增重y。(注:根据统计学原理,样本数据越充分,即得到的X1、X2、X3、X4、X5和y的数据量越大,结论越具有说服力。)计算数据的基本统计量,结果如下表所示:
Figure BDA0001890708420000161
第二步计算X1、X2、X3、X4、X5与y的偏相关性,结果如下表,经spss 计算,从表中结果可以看出X1、X2、X3、X4、X5与y之间均存在较强相关性,
Figure BDA0001890708420000162
Figure BDA0001890708420000171
第三步,根据数据绘制散点图,如图4所示:
在相关性分析的基础之上,根据散点图特点(主要看第一行的散点趋势),结合曲线估计可确认各成分与生产性能之间关系大致为立方函数。通过非线性回归分析得出回归方程基本结构:
Y=m+n*x1*x2*x3*x4*x5+a1*x1*x1*x1+a2*x1*x1+a3*x1+b1*x2*x2*x2+b2*x 2*x2+b3*x2+c1*x3*x3*x3+c2*x3*x3+c3*x3+d1*x4*x4*x4+d2*x4*x4+d3*x4+e1* x5*x5*x5+e2*x5*x5+e3*x5
指定各参数系数的初始数值。经spss计算,得到回归方程中各参数系数(即 X1、X2、X3、X4、X5与y在回归方程中的系数)结果如下:
Figure BDA0001890708420000172
Figure BDA0001890708420000181
由结果可知,R2=0.79即回归方程拟合度约为80%,拟合度较为理想,此回归函数可以信任。带入各参数系数,得到回归方程如下:
Y=32.994-1.458*x1*x2*x3*x4*x5-0.612*x1*x1*x1+7.506*x1*x1-22.364*x 1-0.819*x2*x2*x2-2.756*x2*x2+24.360*x2+1.848*x3*x3*x3-13.526*x3*x3+1 9.291*x3+7.630*x4*x4*x4-32.987*x4*x4+38.317*x4-0.134*x5*x5*x5+2.059* x5*x5-6.071*x5
回归方程应用
为从数学统计的角度验证选取的五种成分是否为决定生产性能的主要因素,根据已获得的X1、X2、X3、X4、X5与y之间的回归方程,进行如下计算:
(1)选取y值。根据实验数据,选取在现有条件下,添加全部成分时所能获取的最佳生产性能值。把该值作为下一步计算的y值。
(2)把上一步得到的y值带入已知的X1、X2、X3、X4、X5与y之间回归方程,同时根据现有试验条件指定X1、X2、X3、X4、X5的取值范围。编写vb 计算程序,计算在指定范围内是否存在能通过回归方程得到指定y值的X1、X2、 X3、X4、X5。经计算得到符合要求的X1、X2、X3、X4、X5。
(3)计算上一步得到的X1、X2、X3、X4、X5间的比值。在实践中按得到的比值调配X1、X2、X3、X4、X5五种成分的比例,经验证得到的生产性能值和指定的y值不存在显著性差异,继而证明,只选取五种成分同样可以达到选取全部成分所能达到的生产性能值。
(4)本数学模型仅适用于在正常肉仔鸡实践生产范围内,有一定的局限性。
日粮中添加不同配伍的酵母培养物有效物组对肉仔鸡生产性能的影响结果由下表所示:
Figure BDA0001890708420000191
注:同行数据肩表字母全部不相同为差异显著(P<0.05),无肩表或肩标有一个字母相同者即为差异不显著(P>0.05)。
由上表可知,在7~21日龄,与对照组相比,各添加处理组可显著提高肉仔鸡日增重、显著降低料重比(P<0.05),对日采食量也有显著影响(P<0.05)。其中有效物组(ABC)日增重均大于D组(YC24h)(P>0.05),较D组(YC24h) 分别提高1.71%、4.96%和2.69%,XP(0.1%)组日增重略高于有效物组。而各处理组体均显著低于对照组(P<0.05),其中以有效物组C的料重比最低。除有效物组A,对照组的采食量显著高于其他各处理组(P<0.05)。
在22~42日龄,各添加处理组对肉仔鸡采食量、料重比有显著影响 (P<0.05),对日增重影响不显著(P>0.05)。其中对照组的采食量显著高于其他各添加处理组(P<0.05)。与对照组相比,各添加处理组的日增重呈增加趋势,但差异不显著(P>0.05)。而各处理组的料重比较对照组显著降低(P<0.05),但各处理组间的料重比差异不显著(P>0.05)。
在7~42日龄,各添加处理组对肉仔鸡日增重、采食量、料重比均有显著影响(P<0.05),与对照组相比,各添加组的日增重都显著高于对照组(P<0.05),采食量、料重比均显著低于对照组(P<0.05),但各添加组间差异不显著 (P>0.05)。各个处理组间,有效物组C料重比最低,与D组(YC24h)组基本相同。
从多方面生产指标比较酵母培养物有效物组与发酵24h产品的饲养效果比较,有效物组的配伍组合的效果能达到甚至超过发酵产品的应用效果,说明在本试验条件下,筛选得到的生物标志物可认为是酵母培养物有效物组的主要成分组合。
日粮中添加不同处理的酵母培养物对肉仔鸡血清生化指标的影响
饲粮中添加不同处理的酵母培养物对肉仔鸡血清生化指标的影响见下表。结果表明,饲粮中添加不同添加处理组对21日龄和42日龄肉鸡血清总蛋白质、白蛋白、谷草转氨酶、碱性磷酸酶、血清尿素氮及血清球蛋白均有不同程度的影响。
21日龄时,对于血清总蛋白、白蛋白,血清球蛋白各试验处理组含量较对照组均有所提高,但差异并不显著(P<0.05)。对于谷草转氨酶,配伍组AB含量显著小于对照组(P<0.05),与对照组相比,有效物物组A、B组含量分别降低13.71%(P<0.05)、13.39%(P<0.05)。而各试验处理组间无显著差异(P >0.05)。对于碱性磷酸酶,各试验处理组含量均高于对照组,其中有效物组C 含量显著高于对照组与XP组(P<0.05)。对于血清尿素氮,各试验处理组含量均显著低于对照组(P<0.05),其中以E、F组(XP组)含量最低,分别较对照组降低42.33%(P<0.05)、48.17%(P<0.05)。
42日龄时,对于血清总蛋白、谷草转氨酶、碱性磷酸酶及血清球蛋白,各试验处理组含量分别较对照组有所提高,但与对照组相比差异不显著(P> 0.05)。对于白蛋白,D组(YC组)含量显著高于有效物组A(P<0.05),但与其他各组差异不显著(P>0.05)。对于血清尿素氮,有效物组C与E、F组(XP) 含量显著低于对照组(P<0.05),其中以E组最低,C组次之。有效物组ABC与 D组(YC24h)比较无显著差异(P>0.05)。
Figure BDA0001890708420000201
Figure BDA0001890708420000211
注:同行数据肩表字母全部不相同为差异显著(P<0.05),无肩表或肩标有一个字母相同者即为差异不显著(P>0.05)。
日粮中不同添加处理酵母培养物对肉仔鸡免疫器官指数的影响见下表。
由表可知日粮中添加不同处理组对肉仔鸡21日龄及42日龄时脾脏指数、胸腺指数和法氏囊指数的影响差异不显著(P>0.05)。在42日龄时,各添加处理组胸腺指数较对照组呈增加趋势,但并不显著(P>0.05)。
Figure BDA0001890708420000212
Figure BDA0001890708420000221
注:同列数据肩标相同字母者表示差异不显(P>0.05),字母不同者表示差异显著(P<0.05)。
日粮中添加不同配伍的酵母培养物有效物组对肉仔鸡IgA、IgG的影响如下表所示。
由表可知,21与42日龄,与对照组相比,各试验处理组IgA、IgG含量均有所升高。其中21日龄时,F组IgA的含量显著高于对照组(P<0.05),而其余各组间无显著差异(P>0.05)。有效物组C的IgG含量显著高于对照组(P<0.05)。 42日龄时,有效物组B中IgA、IgG的含量均显著高于对照组(P<0.05)。
Figure BDA0001890708420000222
注:同列数据肩标相同字母者表示差异不显(P>0.05),字母不同者表示差异显著(P<0.05)。
日粮中添加不同配伍的酵母培养物有效物组对肉仔鸡血清抗氧化指标的影响如下表所示:
从表可以看出,对于21日龄肉仔鸡,与对照组相比,日粮中添加XP产品的E、F组能显著性的提高血清中GPX,SOD和CAT的活性(P<0.05)。其他添加处理组血清中各抗氧化指标较对照组有增加趋势,但差异不显著(P>0.05)。有效物组B与D组血清中MDA的含量低于对照组与XP组,其中有效物组B血清中的MDA含量最低,比对照组降低了62.27%(P>0.05)。对于42日龄肉仔鸡,与对照组相比,各添加处理组血清中GPX,SOD和CAT活性均有不同程度的提高,而血清中MDA的含量较对照组略有升高,但影响作用不明显(P>0.05)。
Figure BDA0001890708420000231
注:同列数据肩标相同字母者表示差异不显(P>0.05),字母不同者表示差异显著(P<0.05)。
日粮中添加不同配伍的酵母培养物有效物组对肉仔鸡屠宰性能的影响如下表所示:
从表中可知,与对照组相比,各处理组对42d肉仔鸡屠宰率、全净膛屠宰率、胸肌率、腿肌率无显著影响(P>0.05),但有一定的升高趋势。其中,添加处理组与对照组比全净膛率分别提高了5.92%、0.71%、5.41%、8.06%、6.26%、 7.42%。各处理组与对照组比胸肌率分别提高了6.79%、9.81%、16.07%、15.68%、 15.36%、13.98%。而腹脂率DEF组较对照组有所下降(P>0.05),有效物组B的腹脂率显著高于E组(XP组),而各处理组间差异不显著(P>0.05)。
Figure BDA0001890708420000241
注:同列数据肩标相同字母者表示差异不显(P>0.05),字母不同者表示差异显著(P<0.05)。
样品reads数及OTUs统计
对每个样品的reads数进行统计,总共得到1246374个序列(如下表)。用Uclstvl.2.22软件根据序列相似性进行聚类得到OTU时选择97%作为相似性阔值(Edgar等,2010)。
Figure BDA0001890708420000242
Figure BDA0001890708420000251
注:DZ1、2、3、4、5分别为21曰龄对照组肉鸡盲肠内容物;ABC Z1、2、3、 4、5分别为21曰龄配伍组肉鸡盲肠内容物;D-12345分别为21曰龄YC(24h) 肉鸡盲肠内容物;E、F-12345分别为21曰龄XP组肉鸡盲肠内容物(下同)
由结果可知每个样本的Reads数在30,000以上,最高可达40,000多,可满足后续试验需要。各组肉鸡盲肠微生物的OTUs的数量范围在4000~8000之间,从不同日粮处理组来看,有效物组(ABC)盲肠内容物中菌群的多样性高于 XP组,与D组(YC24h)多样性水平相当。
物种多样性指数分析
Alpha多样性常用的度量指数主要包括Chao1指数、ACE指数,Shannon指数和Simpson指数。一般而言,Chao1或ACE指数越大,表明群落的丰富度越高。 Shannon多样性指数综合考虑了群落的丰富度和均匀度。Shannon指数值越高,表明群落的多样性越高。Simpson指数值越高,表明群落多样性越高。
样品名 Simpson Chao1 ACE Shannon
DZ 0.977655 1955.33 2079.972 7.956
A 0.956427 1587.318 1658.108 7.578
B 0.956427 1587.318 1658.108 7.578
C 0.983093 1552.804 1595.392 8.214
D 0.971534 1571.862 1579.36 7.735
E 0.980181 1565.542 1588.652 8.046
F 0.978281 1669.62 1770.232 8.004
由表可知,Chao1、ACE指数,对照组均高于试验处理组,但差异不显著。说明对照组微生物的数目要多于其他实验组,有效物组(ABC)与D组(YC24h) 丰度基本相同。对于Simpson、Shannon指数,有效物组C高于其他各组,但各组间差异不显著,说明配伍组C微生物多样性要丰富一些,并且该组中各物种之间个体分配越均匀。同时说明有效物组C各物质比例调配能够增加盲肠的菌群多样性。而添加XP组,也同样较对照组丰富了菌群的丰度及多样性。
本发明思路清晰,通过试验后得出结论:酵母培养物有效物组对肉仔鸡生长性能明显提高;日粮中添加不同处理组对肉仔鸡21日龄及42日龄时脾脏指数、胸腺指数和法氏囊指数的影响差异不显著(P>0.05)。在42日龄时,各添加处理组胸腺指数较对照组呈增加趋势,但并不显著(P>0.05);添加酵母培养物有效物组显著提高了21日龄与42日龄肉鸡的体液免疫指标(P<0.05);酵母培养物有效物组血清中各抗氧化指标较对照组有增加趋势,但差异不显著(P >0.05)。而在21日龄时有效物组B与YC组D血清中MDA的含量低于对照组与 XP组,其中有效物组B血清中的MDA含量最低,比对照组降低了62.27%(P >0.05)。各处理组对42d肉仔鸡屠宰率、全净膛屠宰率、胸肌率、腿肌率无显著影响(P>0.05),但有一定的升高趋势。以97%作为相似性阐值统计得到8000 多个OTU,且不同组别OUT数量不相同,酵母培养物及其配伍处理组较对照组数量有所减少,说明酵母培养物的添加降低了盲肠菌群的多样性,但添加处理组间并没有显著性变化。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质上对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。

Claims (3)

1.一种应用代谢组学技术评价酵母培养物有效物组的方法,其特征在于,利用代谢组学的检测技术对酵母培养物有效物组进行解析,并将多变量降维的数学思路运用到挖掘有效物组的组成中,参照筛选解析的有效物组成分,采用化学单品模拟有效物组的构成关系进行配伍组合,再通过动物饲养试验对有效物组的功效进行验证;
所述的代谢组学的检测技术包括样品采集和处理、代谢物的鉴定和多元统计分析;
所述的采用化学单品模拟有效物组的构成关系进行配伍组合具体表现为回归模型的建立与应用,所述的模型建立对生产性能中的日增重指标与筛选到的酵母培养物有效物组之间关系进行研究,利用多元统计分析方法,建立生产性能与酵母培养物有效物组的函数关系,从而进一步揭示有效物组的组合配伍与生产性能的关系;
包括以下步骤:
(1)筛选酵母培养物有效物组
采用单因素完全随机分组设计,选取336羽1日龄健康体重相近的AA爱拔益加肉仔鸡,公母各半随机分为7组,每组6个重复,每个重复8只鸡;试验期为42d;
酵母培养物有效物组筛选分组:对照组饲喂基础日粮,在基础日粮中添加5个时间差异化培养物的酿酒酵母培养物,各处理组添加量为0.192%,以及达农威益康V XP;
DZ为对照组饲喂基础日粮,ABC组为酵母培养物有效物组,即分别饲喂在基础日粮基础上添加有效物组3个配伍组合,D组饲喂发酵24h的酵母培养物24hYC,E、F组为V XP组,在基础日粮基础上分别添加0.1%、1%“V XP”饲料添加剂;整个试验周期饲粮中不添加抗生素与其他促生长剂;
所述的有效物组为甘氨酸、果糖、肌醇、吡喃型半乳糖、蔗糖;试验日粮分配表如下表所示:
试验组 DZ组 A组 B组 C组 D组 E组 F组 基础日粮 试验日粮 0对照 配伍1 配伍2 配伍3 24hYC 0.1%V XP 1%V XP
ABC组在YC中所占比例配伍情况如下表所示:
Figure FDA0003532851050000021
各物质配伍说明:取生产性能为发酵24hYC产品的日增重,得到相应的X值,即配伍1有效物组A;在此基础上调整各物质的比例,所述的比例为各物质±10%,配置各物质的含量,为配伍2和3有效物组B、C组;
(2)血清样品采集:采样对象为21日龄和42日龄时,随机选择2只健康和体重相近的鸡;采集到的血清样品于室温放置,待血清析出后将其分装于离心管中,3000r/min离心10min,取其上清分装于Eppendorf试管中,并立即浸入液氮,置-80℃低温保存,用于生化分析用;
器官组织样品采集:采集对象为21日龄和42日龄时,随机选择2只健康和体重相近的试验鸡,采集到的器官组织样品,剔除脂肪后称鲜重并记录,样品保存在-80℃待用;
肌肉样品的采集:于试验第42日结束时每个重复随机选取与该重复平均体重相近的肉鸡2只,采集胸肌和腿肌样品,去除脂肪和结缔组织,称重并记录;
(3)生产性能:试验在21和42日龄时,以重复为单位,分别于试验开始的1、7、14、21、28、35、42d早晨,记录鸡初始重和末重,并记录整个试验期供料量、剩余料量和损失料量,计算平均日采食量ADFI、平均日增重ADG和料重比F/G;
免疫器官指数:在试验的21d和42d,采集免疫器官并称重、记录,器官指数=器官质量mg/活体质量g;
血清生化指标:使用全自动血液生化分析仪测定血中白蛋白,球蛋白,血清尿素氮与碱性磷酸酶;
血清免疫指标:采集对象为21、42日龄鸡,收集血液,3500r/m,4℃离心10min,取血清于1mL离心管,至于-20℃保存,用于血清总蛋白TP、白蛋白ALB、尿素氮UREA、葡萄糖GLU等测定,检测使用全自动生化分析仪;分泌性免疫球蛋白A IgA、分泌性免疫球蛋白G IgG用ELISA试剂盒南京建成生物工程研究所测定;
血清抗氧化指标:血清中超氧化物歧化酶SOD活性、谷胱甘肽过氧化物酶GPX活性和丙二醛MDA含量使用南京建成生物工程研究所生产的试剂盒测定;
屠宰性能测定:在试验42d时,每个重复各随机取2只鸡;活体重:禁食12h,测试前后均称重;全净膛重为屠体重去除气管、食道、嗦囊、肠、脾脏、胰腺、胆囊和生殖器官心、肝、肾、腺胃、肌胃、腹脂、头、脚的后的重量;分割左侧胸肌、腿肌并称重,分别为胸肌重和腿肌重;腹脂重主要是腹部板油及肌胃周围的脂肪的重量;计算公式:腹脂率=腹脂重/全净膛重×100%;胸肌率=胸肌重/全膛净重×100%腿肌率=腿肌重/全膛净重×100%;全净膛率=全膛净重/活体重×100%;
(4)DNA提取
称取21日龄肉鸡的盲肠内容物0.3g于无菌PBS中,涡轮震荡30S后,10000rpm转速下离心10min后取上清溶于DNA缓冲液中;肠道内容物DNA提取方法参照试剂盒QIAamp Fast DNAStool Mini Kit操作步骤QIAGEN,German,样品16S rRNA基因V3-V4区的PCR扩增及Illumina高通量测序由上海派森诺科技股份有限公司完成;
16Sr RNA基因V3~V4区的PCR扩增
将样品DNA稀释至20ng/μL作为模板,扩增细菌16S rRNA V3-4区;使用常规引物338F5′-ACT CCTAC G GGAGGCAGCA-3′和806R 5′-GGA CTAC H VGGGTWTCTAAT-3′,MiSeq测序平台检测;
PCR扩增体系为25μL:5×reaction buffer 5μL,5×GC buffer 5μL,dNTP 2.5mM 2μL,正向引物10μM 1μL,反向引物10μM 1μL,模板DNA 2μL,ddH2O 8.75μL,Q5 DNA酶0.25μL;PCR扩增条件:98℃预变性2min;然后98℃变性15s,55℃退火30s,72℃延伸30s,25-30个循环;最后72℃延伸5min,10℃保持待测定,PCR扩增产物通过2%琼脂糖凝胶电泳进行检测,采用AXYGEN公司的凝胶回收试剂盒切胶回收;将PCR扩增回收产物通过荧光试剂为Quant-iTPicoGreen dsDNA Assay Kit进行定量;
高通量测序
采用MiSeq测序仪进行2×300bp的双端测序,相应试剂为MiSeq Reagent Kit V3600cycles;盲肠菌群多样性检测由派森诺生物上海完成;
(5)数据处理与分析:
生产应用试验数据经Excel 2007初步整理后,采用SPSS 23.0统计软件中One-WayANOVA进行方差分析,多重比较采用Duncan's法;P<0.05为差异显著,P<0.01为差异极显著,数据用“平均数±标准差Mean±SD”表示;
采用Illumina MiSeq平台对群落DNA片段进行双端Paired-end测序;首先对FASTQ格式的双端序列逐一作质量筛查,要求截断后的序列长度≥150bp,且不允许存在模糊碱基;再利用FLASH软件,对通过质量初筛的双端序列进行配对连接;最后,根据样本所对应的Index信息,将连接后的序列识别分配入对应样本,获得样本的有效序列;
运用QIIME软件识别疑问序列,并剔除:5’端引物错配碱基数>1的序列;含有连续相同碱基数>8的序列;再调用USEARCH检查并剔除嵌合体序列;
Alpha多样性指数:用于比较不同样本的多样性,首先对序列量拉平处理,从而校正测序深度引起的多样性差异;再使用QIIME软件分别对每个样本计算侧重于体现群落丰富度Chao1指数和ACE指数,以及兼顾群落均匀度的Shannon指数和Simpson指数。
2.根据权利要求1所述的应用代谢组学技术评价酵母培养物有效物组的方法,其特征在于,步骤(1)中所述的基础日粮的组成由下表所示:
日粮组成 1-21d 22-42d 营养成分 1-21d 22-42d 玉米 52.07 55.57 代谢能MJ/kg 12.746 12.719 大豆粕 35.50 34.00 粗蛋白% 20.722 19.609 鱼粉 5.00 3.60 钙% 0.910 0.865 大豆油 4.00 3.5 总磷% 0.676 0.584 磷酸氢钙 1.00 0.70 有效磷% 0.447 0.356 石粉 1.00 1.20 赖氨酸% 1.347 0.459 食盐 0.30 0.30 蛋氨酸% 1.248 0.432 赖氨酸 0.04 0.04 蛋氨酸 0.09 0.09 1%预混料 1.00 1.00 合计 100 100
其中,每千克预混料包含:维生素A 5000IU;维生素D3 1500IU;维生素E15IU;维生素K30.8mg;维生素B12 0.01mg;叶酸0.5mg;烟酸50mg;泛酸8mg;生物素0.1mg;吡哆醇2.2mg;核黄素4.4mg;硫胺素1.6mg。
3.根据权利要求1所述的应用代谢组学技术评价酵母培养物有效物组的方法,其特征在于,步骤(1)中所述的肉仔鸡购自吉林德详公司,酵母培养物由吉林农业大学吉农博瑞研发中心提供,所述的V XP为达农威中国有限公司的益康V XP。
CN201811470093.9A 2018-11-28 2018-11-28 一种应用代谢组学技术评价酵母培养物有效物组的方法 Active CN109856384B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811470093.9A CN109856384B (zh) 2018-11-28 2018-11-28 一种应用代谢组学技术评价酵母培养物有效物组的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811470093.9A CN109856384B (zh) 2018-11-28 2018-11-28 一种应用代谢组学技术评价酵母培养物有效物组的方法

Publications (2)

Publication Number Publication Date
CN109856384A CN109856384A (zh) 2019-06-07
CN109856384B true CN109856384B (zh) 2022-06-07

Family

ID=66890542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811470093.9A Active CN109856384B (zh) 2018-11-28 2018-11-28 一种应用代谢组学技术评价酵母培养物有效物组的方法

Country Status (1)

Country Link
CN (1) CN109856384B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113487039B (zh) * 2021-06-29 2023-08-22 山东大学 基于深度强化学习的智能体自适应决策生成方法及系统
CN114216982A (zh) * 2021-12-14 2022-03-22 黑龙江飞鹤乳业有限公司 维生素b12的检测方法
CN114720618A (zh) * 2022-03-30 2022-07-08 集美大学 一种法夫酵母代谢产物的分析方法
CN117757891B (zh) * 2024-02-22 2024-05-31 潍坊华卓生物科技有限公司 预防h9亚型禽流感病毒的功能益生菌反向筛选方法及应用

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149280A (zh) * 2011-12-07 2013-06-12 中国农业大学 一种应用代谢组学技术评价动物个体营养状况的方法
CN103245747A (zh) * 2013-05-23 2013-08-14 中国农业大学 一种应用代谢组学技术预测猪日粮代谢能水平的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149280A (zh) * 2011-12-07 2013-06-12 中国农业大学 一种应用代谢组学技术评价动物个体营养状况的方法
CN103245747A (zh) * 2013-05-23 2013-08-14 中国农业大学 一种应用代谢组学技术预测猪日粮代谢能水平的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
酵母培养物对肉仔鸡生长性能及盲肠菌群的影响;甄玉国等;《中国畜牧杂志》;20170610(第06期);全文 *
酿酒酵母培养物对肉仔鸡生长性能及免疫功能的影响;孙喆等;《中国兽医学报》;20180630;第38卷(第6期);第1222-1227页 *
高通量测序分析酵母培养物对肉仔鸡盲肠菌群的影响;孙喆等;《中国畜牧杂志》;20171210;第53卷(第12期);全文 *

Also Published As

Publication number Publication date
CN109856384A (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
CN109856384B (zh) 一种应用代谢组学技术评价酵母培养物有效物组的方法
Gong et al. Effects of zinc bacitracin, bird age and access to range on bacterial microbiota in the ileum and caeca of broiler chickens
Kasparovska et al. Effects of isoflavone-enriched feed on the rumen microbiota in dairy cows
de Melo Pereira et al. An updated review on bacterial community composition of traditional fermented milk products: what next-generation sequencing has revealed so far?
US20060276973A1 (en) Predicting animal performance
TW201713775A (zh) 用於自複雜的異源群集中分析微生物品系、預測及識別其機能性關係及交互作用,且依據前述選擇並合成微生物系集之方法、裝置及系統
Neumann et al. Differences in major bacterial populations in the intestines of mature broilers after feeding virginiamycin or bacitracin methylene disalicylate
Sales et al. Heat stress influence the microbiota and organic acids concentration in beef cattle rumen
Wang et al. Effects of Bacillus subtilis and coccidial vaccination on cecal microbial diversity and composition of Eimeria-challenged male broilers
JP2020503048A (ja) トレーサー分析論による複雑な不均一コミュニティの微生物株の解析、その機能的関連性及び相互作用の決定、ならびに微生物アンサンブル(投与される微生物アンサンブル及び接種される微生物アンサンブルを含む)の合成、のための方法、装置、及びシステム
Betiku et al. Water system is a controlling variable modulating bacterial diversity of gastrointestinal tract and performance in rainbow trout
He et al. Determination of the available energy, standardized ileal digestibility of amino acids of fermented corn germ meal replacing soybean meal in growing pig diets
Olson et al. Application of microbial analyses to feeds and potential implications for poultry nutrition
Bao et al. Diversity analysis of bacterial and function prediction in hurunge from mongolia
Bao et al. Diversity analysis and gene function prediction of bacteria and fungi of Bactrian camel milk and naturally fermented camel milk from Alxa in Inner Mongolia
Lin et al. Comparison of growth performance and rumen metabolic pathways in sheep and goats under the same feeding pattern
Seidel et al. Impact of concentrations of camphor on the in vitro mixed ruminal microorganism fermentation from goats selected for consumption of low and high levels of Juniperus spp.-2
Milanov et al. Toxin genotypes of Clostridium perfringens in animal feed and their role in the ethiology of enterotoxemia in domestic animals
Rzeznitzeck et al. Morphology, microbiota, and metabolome along the intestinal tract of female turkeys
Zened et al. Milk replacer feeding once or twice a day did not change the ruminal metabolomic profile and the microbial diversity of dairy calves from birth to weaning
Chai et al. Effect of weaning time on growth performance and rumen development of Hu lambs
KR102169139B1 (ko) 삼채 잎의 분말을 유효성분으로 포함하는 가금류의 장내 균총 개선용 조성물
He et al. Vitamin-B12 enrichment in tempeh by co-culture with Propionibacterium freudenreichii during fermentation
US20220218764A1 (en) Methods for improving treatment of equine colic by administration of a synthetic bioensemble or purified strains thereof
RU2694256C2 (ru) Кормовая добавка для профилактики бактерионосительства микроорганизмов рода Salmonella у сельскохозяйственной птицы и способ ее применения

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 130000 Lanjia Street, Kuancheng District, Changchun City, Jilin Province, 1777

Applicant after: Changchun Borui Technology Co.,Ltd.

Address before: 130000 Section 161, floor 5, No. 57, Changjiang Road, Changchun Economic Development Zone, Jilin Province

Applicant before: CHANGCHUN BORUI AGRICULTURE AND ANIMAL HUSBANDRY GROUP CO.,LTD.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A method for evaluating the effective composition of yeast culture using metabolomics techniques

Granted publication date: 20220607

Pledgee: China Construction Bank Co.,Ltd. Changchun Science and Technology Sub branch

Pledgor: Changchun Borui Technology Co.,Ltd.

Registration number: Y2024220000009

PE01 Entry into force of the registration of the contract for pledge of patent right