高精料日粮下湖羊结肠的适应性机制：微生物、上皮与功能的多维解析

高精料日粮下湖羊结肠的适应性机制：微生物、上皮与功能的多维解析一、引言1.1研究背景与意义随着我国经济水平的提升，人们对肉类、奶类等反刍动物产品的需求日益增长。在反刍动物养殖领域，为满足其高营养需求并确保高产，通常采用高精料日粮进行喂养。高精料日粮富含大量易发酵的碳水化合物，能够显著提升动物的生产性能。在奶牛养殖中，高精料日粮可使奶牛的产奶量明显增加；在肉羊育肥阶段，能加快肉羊的生长速度，缩短育肥周期，提高养殖经济效益。然而，大量研究表明，长期饲喂高精料日粮会对反刍动物的健康产生诸多不利影响。当反刍动物摄入过量的高精料日粮时，瘤胃内会发生剧烈的发酵过程。瘤胃微生物迅速分解淀粉，产生大量的挥发性脂肪酸（VFA），如乙酸、丙酸和丁酸等，同时乳酸的生成量也大幅增加。这使得瘤胃液的pH值急剧下降，当pH值低于5.8且持续3小时以上，就会引发亚急性瘤胃酸中毒（SARA）。SARA会导致瘤胃代谢紊乱，影响瘤胃微生物的正常生长和繁殖，降低瘤胃对饲料的消化和利用效率。过量的过瘤胃淀粉进入小肠和大肠后，会引起后肠道发酵参数的改变，如pH值降低、VFA浓度变化等，进而损害肠道上皮结构和功能。湖羊作为我国特有的绵羊品种，具有早熟、四季发情、一胎多羔、繁殖力高、泌乳性能好、生长发育快、肉质好、耐湿热等优良特性。目前，湖羊的养殖规模不断扩大，在我国肉羊产业中占据着重要地位。结肠作为湖羊消化系统的重要组成部分，对维持其肠道健康和整体生理功能具有关键作用。结肠内存在着复杂而多样的微生物区系，这些微生物与湖羊形成了紧密的共生关系。它们参与饲料的发酵和消化，产生VFA等代谢产物，为湖羊提供额外的能量来源；同时，还能合成维生素K、维生素B族等营养物质，对湖羊的生长发育和免疫功能有着重要影响。结肠上皮形态和功能的正常与否，直接关系到肠道的屏障功能、营养物质的吸收以及免疫调节等过程。正常的结肠上皮能够有效地阻止病原体的入侵，维持肠道内环境的稳定；同时，还能高效地吸收VFA等营养物质，为湖羊的生命活动提供支持。然而，高精料日粮的饲喂可能会打破湖羊结肠内的微生物平衡，导致有益菌数量减少，有害菌大量繁殖，进而引发一系列肠道健康问题。高精料日粮还可能对结肠上皮的形态和功能造成损害，影响其屏障功能和营养物质吸收能力。因此，深入研究湖羊结肠微生物区系及上皮形态和功能对高精料日粮的适应性应答，对于揭示高精料日粮对湖羊肠道健康的影响机制，制定合理的饲养管理策略，提高湖羊的生产性能和健康水平具有重要的理论和实践意义。本研究旨在系统地探究湖羊结肠微生物区系及上皮形态和功能对高精料日粮的适应性应答规律。通过分析不同精料水平下湖羊结肠发酵参数、微生物区系组成和功能、上皮形态结构以及相关基因表达的变化，揭示高精料日粮对湖羊结肠健康的影响机制。研究结果将为优化湖羊日粮配方、改善肠道健康、提高养殖效益提供科学依据，同时也为反刍动物营养研究领域提供新的理论参考，推动反刍动物养殖产业的可持续发展。1.2国内外研究现状在反刍动物养殖领域，高精料日粮的应用十分广泛。大量研究聚焦于高精料日粮对反刍动物瘤胃健康的影响。在瘤胃内环境方面，高精料日粮会显著改变瘤胃发酵模式及代谢产物的种类和浓度。李洋对山羊的研究表明，高精料饲粮使山羊瘤胃内乙酸浓度、乙酸/丙酸下降，丙酸和丁酸浓度上升，并产生大量组胺、脂多糖（LPS）、乙醇胺和D-乳酸等异常代谢产物，进而引起瘤胃代谢紊乱。Zhang等通过整合代谢组学和相关网络分析，发现高精料饲粮增加了奶牛瘤胃液中的细菌降解产物（包括黄嘌呤、次黄嘌呤、尿嘧啶等）、有毒化合物（包括脂多糖、生物胺、乙醇胺等）和15种氨基酸（包括丙氨酸、亮氨酸、甘氨酸等）的浓度，这些代谢产物的变化可能增加亚急性瘤胃酸中毒（SARA）的发生风险。瘤胃微生物区系也会受到高精料日粮的显著影响。Plaizier等基于16SrRNA基因V3-V4区域高通量测序分析奶牛瘤胃液的菌群，共观察到11个不同门，90种细菌群落；在门水平上，高精料饲粮增加了瘤胃液中厚壁菌门和放线菌门的相对丰度，而降低了蓝细菌门和拟杆菌门的相对丰度。厚壁菌门具有消化代谢淀粉类碳水化合物的能力，其相对丰度的增加会加速淀粉的降解，不利于高精料饲喂时瘤胃内环境的健康；而拟杆菌门能降解纤维素，且与瘤胃液pH关系密切，其相对丰度的下降可能与瘤胃液pH下降有关。Fernando等使用末端限制性片段长度多态性（T-RFLP）分析肉牛瘤胃液16SrRNA基因，发现谷物和干草比为20∶80、40∶60与60∶40、80∶20的细菌种群结构存在明显差异，且精料比例越高，瘤胃微生物菌群的变化越明显；在门水平上，低精料饲粮中肉牛瘤胃中含有大量的纤维杆菌，随着精料的增加，琥珀酸纤维杆菌的相对丰度逐渐降低，而高精料饲粮组肉牛瘤胃中含有大量的厚壁菌和变形菌。高精料日粮对反刍动物瘤胃上皮也会造成损伤。Steele等用高精料饲粮诱导羔羊发生SARA时发现其瘤胃上皮出现过度角质化，且瘤胃角质层厚度、乳头高度及宽度都呈下降趋势，同时还会破坏上皮细胞间的紧密连接，使瘤胃上皮细胞间的间隙增大，最终削弱瘤胃上皮的屏障功能。在对奶牛的研究中也发现，高精料饲粮组奶牛的瘤胃颗粒层、基底层和棘突层厚度以及活检瘤胃乳头中的整个复层鳞状瘤胃上皮（SSE）厚度均明显降低；通过电子显微镜高放大倍数观察发现，高精料组瘤胃颗粒层的桥粒数目减少，细胞黏附程度降低，表明颗粒细胞没有形成紧密的屏障。相比瘤胃，关于高精料日粮对反刍动物后肠，尤其是结肠的研究相对较少。在结肠微生物区系方面，有研究表明高精料日粮会改变结肠内容物及黏膜微生物区系。林丽梅等把15只湖羊公羊随机分为3组，分别饲喂非颗粒低谷物饲粮（CON）、非颗粒高谷物饲粮（HG）和颗粒高谷物饲粮（HP），结果表明，HG和HP组的发酵参数高于CON组，尤其是乙酸浓度；与CON组相比，HG饲喂组改变了结肠中共生细菌的组成，包括乙酸生成菌（如Acetitomaculumspp.）、丁酸生成菌（如Coprococcusspp.和Subdoligranulumspp.）和淀粉降解菌（如Prevotellaspp.，Roseburiaspp.，和Oscillibacterspp.）的增加。但目前对于结肠微生物区系在高精料日粮作用下的动态变化规律以及微生物之间的相互作用关系，仍缺乏深入系统的研究。在结肠上皮形态及屏障功能方面，高精料日粮对结肠上皮形态结构和屏障功能会产生影响。有研究通过光镜、电镜观察奶山羊后段肠道上皮屏障的形态和紧密连接的情况，发现高精料日粮会对其产生影响，但对于湖羊结肠上皮在高精料日粮下的具体形态学变化，如绒毛高度、隐窝深度等指标的研究还不够细致，且对其屏障功能受损的分子机制研究较少。在结肠上皮挥发性脂肪酸（VFA）吸收和代谢功能方面，高精料日粮会改变结肠内VFA的产生和组成，进而影响结肠上皮对VFA的吸收和代谢。然而，目前对于高精料日粮影响湖羊结肠上皮VFA吸收和代谢的相关基因表达及信号通路的研究还不够深入，缺乏全面系统的认识。1.3研究目标与内容本研究旨在深入探究湖羊结肠微生物区系及上皮形态和功能对高精料日粮的适应性应答机制，为湖羊的科学饲养提供理论依据。具体研究内容如下：湖羊结肠发酵及上皮形态和功能对高精料日粮的适应性应答：选取健康、体重相近的湖羊若干只，随机分为不同组，分别饲喂不同精料水平的日粮。在试验期内，定期采集湖羊的结肠内容物和组织样品。测定结肠内容物的pH值、挥发性脂肪酸（VFA）浓度、氨态氮含量等发酵参数，分析精料水平对结肠发酵的影响。通过制作结肠组织切片，利用苏木精-伊红（HE）染色，在显微镜下观察结肠上皮的绒毛高度、隐窝深度等形态学指标的变化；采用免疫组织化学或Westernblot等方法检测结肠黏膜紧密连接蛋白（如ZO-1、Occludin等）、上皮细胞增殖与凋亡相关基因（如PCNA、Bax、Bcl-2等）以及VFA吸收与代谢相关基因（如MCT1、CD36等）的表达水平，探究高精料日粮对结肠上皮屏障功能、细胞增殖与凋亡以及VFA吸收和代谢功能的影响。湖羊结肠内容物细菌菌群对高精料日粮的适应性应答：对上述试验湖羊的结肠内容物样品进行DNA提取，利用IlluminaMiSeq高通量测序技术对16SrRNA基因的特定区域进行测序，分析结肠内容物细菌菌群的丰富度、多样性以及群落组成的变化。通过PICRUSt等软件进行菌群功能预测，探讨高精料日粮对结肠内容物细菌菌群功能的影响。采用实时荧光定量PCR（qPCR）技术测定结肠内容物中总菌、厚壁菌门和拟杆菌门的基因拷贝数，进一步了解细菌数量的变化情况。湖羊结肠黏膜微生物区系和炎症因子表达对高精料日粮的适应性应答：提取湖羊结肠黏膜样品的DNA，进行IlluminaMiSeq测序，分析结肠黏膜微生物菌群的丰度、多样性和群落组成的变化。同样利用PICRUSt软件进行菌群功能预测，研究高精料日粮对结肠黏膜微生物菌群功能的影响。通过qPCR技术测定结肠黏膜中总菌、厚壁菌门和拟杆菌门的基因拷贝数。提取结肠黏膜的RNA，反转录为cDNA后，利用qPCR检测炎症因子相关基因（如IL-1β、IL-6、TNF-α等）的表达水平，分析高精料日粮对结肠黏膜炎症反应的影响，并通过相关性分析探讨微生物区系与炎症因子表达之间的关系。二、高精料日粮对湖羊结肠发酵及上皮形态和功能的影响2.1材料与方法2.1.1实验动物与饲养管理选取18只健康状况良好、体重相近（约为[X]kg）的6月龄湖羊公羊作为实验动物。实验羊被安置在江苏省农业科学院畜牧研究所的标准化羊舍内，羊舍保持通风良好、光照充足，温度控制在[X]℃-[X]℃，相对湿度维持在[X]%-[X]%。实验羊采用单栏饲养，以避免相互干扰，每只羊均配备独立的采食槽和饮水器，确保自由采食和饮水。2.1.2试验设计本实验采用单因素完全随机设计，将18只湖羊随机分为3组，每组6只羊。对照组（CON组）饲喂低精料日粮，精料与粗料的比例为30:70；实验组1（HG1组）饲喂中等精料日粮，精料与粗料的比例为50:50；实验组2（HG2组）饲喂高精料日粮，精料与粗料的比例为70:30。实验所用的精料主要由玉米、豆粕、麦麸等组成，粗料为苜蓿干草和羊草按照1:1的比例混合而成。日粮的营养水平根据NRC（2007）绵羊营养需要标准进行配制，具体日粮组成及营养水平如表1所示：原料CON组HG1组HG2组玉米（%）[X][X][X]豆粕（%）[X][X][X]麦麸（%）[X][X][X]苜蓿干草（%）[X][X][X]羊草（%）[X][X][X]预混料（%）[X][X][X]营养水平干物质（%）[X][X][X]粗蛋白（%）[X][X][X]粗脂肪（%）[X][X][X]中性洗涤纤维（%）[X][X][X]酸性洗涤纤维（%）[X][X][X]代谢能（MJ/kg）[X][X][X]实验分为预试期和正试期两个阶段。预试期为14天，在此期间对实验羊进行驱虫、防疫等处理，并让其适应实验环境和基础日粮；正试期为56天，在正试期内按照实验设计给各组羊饲喂相应的日粮。每天分别于08:00和16:00进行两次等量饲喂，保证每只羊每天的采食量相同，同时记录每只羊的采食量和剩余饲料量。2.1.3样品采集在正试期结束后的清晨，羊只空腹状态下，每组随机选取3只羊，采用颈静脉放血的方法进行屠宰。迅速采集结肠内容物和结肠组织样品。结肠内容物采集约10g，置于无菌冻存管中，立即放入液氮中速冻，随后转移至-80℃冰箱保存，用于结肠发酵参数测定和微生物区系分析；结肠组织样品采集约1cm×1cm大小，一部分用4%多聚甲醛溶液固定，用于制作组织切片观察上皮形态；另一部分迅速放入液氮中速冻，然后转移至-80℃冰箱保存，用于提取RNA，检测相关基因的表达水平。2.1.4样品测定结肠发酵参数的测定：将结肠内容物解冻后，称取1g，加入9mL无菌生理盐水，涡旋振荡10min，使内容物充分混匀。然后将混合液在4℃、12000r/min条件下离心10min，取上清液用于测定pH值、挥发性脂肪酸（VFA）浓度和氨态氮含量。pH值使用pH计直接测定；VFA浓度采用气相色谱仪（GC-2014C，岛津公司，日本）进行测定，色谱柱为毛细管柱（HP-INNOWAX，30m×0.32mm×0.25μm），进样口温度为250℃，检测器温度为280℃，初始柱温为100℃，保持3min，然后以10℃/min的速率升温至230℃，保持5min；氨态氮含量采用苯酚-次氯酸钠比色法进行测定。结肠组织切片处理：将固定好的结肠组织依次经过乙醇梯度脱水（70%、80%、90%、95%、100%乙醇各处理30min）、二甲苯透明（二甲苯Ⅰ、二甲苯Ⅱ各处理20min）、石蜡包埋等步骤，制成石蜡切片。切片厚度为5μm，采用苏木精-伊红（HE）染色法进行染色，然后在光学显微镜（BX53，奥林巴斯公司，日本）下观察结肠上皮的绒毛高度、隐窝深度等形态学指标。绒毛高度从绒毛顶端到绒毛与隐窝交界处的垂直距离，隐窝深度从隐窝开口处到隐窝底部的垂直距离，每个切片随机选取10个视野进行测量，取平均值。黏膜RNA提取和RT-qPCR：使用Trizol试剂（Invitrogen公司，美国）提取结肠黏膜组织的总RNA，按照试剂盒说明书的操作步骤进行。提取的RNA用超微量分光光度计（NanoDrop2000，ThermoFisherScientific公司，美国）测定浓度和纯度，要求A260/A280比值在1.8-2.0之间。取1μg总RNA，利用反转录试剂盒（PrimeScriptRTreagentKitwithgDNAEraser，TaKaRa公司，日本）将其反转录为cDNA。以cDNA为模板，采用SYBRGreen荧光定量PCR试剂盒（TBGreenPremixExTaqII，TaKaRa公司，日本）在实时荧光定量PCR仪（CFX96Touch，Bio-Rad公司，美国）上进行qPCR反应。反应体系为20μL，包括10μLSYBRGreenPremixExTaqII、0.8μL上游引物、0.8μL下游引物、2μLcDNA模板和6.4μLddH2O。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3个重复，以β-actin作为内参基因，采用2-ΔΔCt法计算目的基因的相对表达量。本实验检测的目的基因包括紧密连接蛋白（ZO-1、Occludin）、上皮细胞增殖与凋亡相关基因（PCNA、Bax、Bcl-2）以及VFA吸收与代谢相关基因（MCT1、CD36），引物序列如表2所示：|基因名称|引物序列（5'-3'）||----|----||ZO-1|F:[X]R:[X]||Occludin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||基因名称|引物序列（5'-3'）||----|----||ZO-1|F:[X]R:[X]||Occludin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||----|----||ZO-1|F:[X]R:[X]||Occludin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||ZO-1|F:[X]R:[X]||Occludin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||Occludin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||PCNA|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||Bax|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||Bcl-2|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||MCT1|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||CD36|F:[X]R:[X]||β-actin|F:[X]R:[X]||β-actin|F:[X]R:[X]|2.1.5数据分析实验数据使用SPSS22.0统计软件进行分析，采用单因素方差分析（One-wayANOVA）对各组数据进行差异显著性检验，当P<0.05时认为差异显著，当P<0.01时认为差异极显著。数据结果以“平均值±标准差（Mean±SD）”表示。利用Origin2021软件绘制图表，直观展示数据变化趋势。2.2结果与分析2.2.1湖羊结肠发酵参数各组湖羊结肠内容物的发酵参数测定结果如表3所示。与CON组相比，HG1组和HG2组结肠内容物的pH值显著降低（P<0.05），且HG2组的pH值低于HG1组，差异显著（P<0.05）。这表明随着精料水平的升高，湖羊结肠内容物的酸性增强。在挥发性脂肪酸（VFA）浓度方面，HG1组和HG2组的总VFA、乙酸、丙酸和丁酸浓度均显著高于CON组（P<0.05）。其中，HG2组的总VFA浓度最高，分别比CON组和HG1组增加了[X]%和[X]%；乙酸浓度也显著高于HG1组（P<0.05），丙酸和丁酸浓度在两组间差异不显著（P>0.05）。这说明高精料日粮可促进湖羊结肠内碳水化合物的发酵，产生更多的VFA。在氨态氮含量上，HG1组和HG2组显著高于CON组（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）。表明高精料日粮使湖羊结肠内蛋白质和氨基酸的分解代谢增强，产生更多的氨态氮。组别pH值总VFA（mmol/L）乙酸（mmol/L）丙酸（mmol/L）丁酸（mmol/L）氨态氮（mg/dL）CON组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]注：同行数据肩标不同小写字母表示差异显著（P<0.05），相同或无字母表示差异不显著（P>0.05），下同。2.2.2湖羊结肠上皮形态学湖羊结肠上皮的形态学指标测定结果如表4所示。HG1组和HG2组的结肠绒毛高度显著低于CON组（P<0.05），且HG2组低于HG1组，差异显著（P<0.05）。隐窝深度方面，HG1组和HG2组显著高于CON组（P<0.05），HG2组也显著高于HG1组（P<0.05）。绒毛高度/隐窝深度的比值能反映肠道上皮的消化吸收能力，HG1组和HG2组该比值显著低于CON组（P<0.05），且HG2组低于HG1组，差异显著（P<0.05）。这些结果表明，高精料日粮会导致湖羊结肠绒毛变短，隐窝加深，使绒毛高度/隐窝深度比值降低，进而影响结肠上皮的消化吸收功能。组别绒毛高度（μm）隐窝深度（μm）绒毛高度/隐窝深度CON组[X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X]2.2.3高精料日粮对湖羊结肠黏膜紧密连接蛋白和上皮细胞增殖与凋亡相关基因的表达的影响各组湖羊结肠黏膜紧密连接蛋白和上皮细胞增殖与凋亡相关基因的相对表达量测定结果如表5所示。与CON组相比，HG1组和HG2组结肠黏膜中紧密连接蛋白ZO-1和Occludin的基因表达量显著降低（P<0.05），且HG2组低于HG1组，差异显著（P<0.05）。这表明高精料日粮会破坏湖羊结肠黏膜紧密连接蛋白的表达，导致肠道屏障功能受损。在上皮细胞增殖相关基因PCNA的表达上，HG1组和HG2组显著低于CON组（P<0.05），且HG2组低于HG1组，差异显著（P<0.05），说明高精料日粮抑制了湖羊结肠上皮细胞的增殖。凋亡相关基因Bax的表达量，HG1组和HG2组显著高于CON组（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；而抗凋亡基因Bcl-2的表达量，HG1组和HG2组显著低于CON组（P<0.05），且HG2组低于HG1组，差异显著（P<0.05）。Bax/Bcl-2比值能够反映细胞凋亡的趋势，HG1组和HG2组该比值显著高于CON组（P<0.05），且HG2组高于HG1组，差异显著（P<0.05），表明高精料日粮促进了湖羊结肠上皮细胞的凋亡。组别ZO-1OccludinPCNABaxBcl-2Bax/Bcl-2CON组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X][X]±[X][X]±[X][X]±[X]2.2.4湖羊结肠上皮细胞增殖与凋亡相关基因的表达量和VFA的相关性分析对湖羊结肠上皮细胞增殖与凋亡相关基因的表达量和VFA浓度进行相关性分析，结果如表6所示。PCNA的表达量与总VFA、乙酸、丙酸和丁酸浓度均呈显著负相关（P<0.05），说明VFA浓度的升高会抑制结肠上皮细胞的增殖。Bax的表达量与总VFA、乙酸、丙酸和丁酸浓度均呈显著正相关（P<0.05），Bcl-2的表达量与总VFA、乙酸、丙酸和丁酸浓度均呈显著负相关（P<0.05），Bax/Bcl-2比值与总VFA、乙酸、丙酸和丁酸浓度均呈显著正相关（P<0.05）。这表明VFA浓度的升高会促进结肠上皮细胞的凋亡，且VFA浓度与结肠上皮细胞增殖和凋亡相关基因的表达密切相关。指标总VFA乙酸丙酸丁酸PCNAr=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]Baxr=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]Bcl-2r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]Bax/Bcl-2r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]r=[X]，P=[X]2.2.5高精料日粮对湖羊结肠黏膜VFA吸收与代谢相关基因的表达的影响各组湖羊结肠黏膜VFA吸收与代谢相关基因的相对表达量测定结果如表7所示。与CON组相比，HG1组和HG2组结肠黏膜中VFA转运蛋白MCT1和脂肪酸转运蛋白CD36的基因表达量显著降低（P<0.05），且HG2组低于HG1组，差异显著（P<0.05）。这表明高精料日粮会抑制湖羊结肠黏膜VFA吸收与代谢相关基因的表达，从而影响结肠上皮对VFA的吸收和代谢功能。组别MCT1CD36CON组[X]±[X][X]±[X]HG1组[X]±[X][X]±[X]HG2组[X]±[X][X]±[X]2.3讨论本研究结果表明，高精料日粮对湖羊结肠发酵及上皮形态和功能产生了显著影响。随着精料水平的升高，湖羊结肠内容物的pH值显著降低，这与以往的研究结果一致。高精料日粮中富含大量易发酵的碳水化合物，进入结肠后被微生物迅速发酵，产生大量的挥发性脂肪酸（VFA），导致结肠内容物的酸性增强，pH值下降。研究发现，当瘤胃内pH值低于5.8时，就会引发亚急性瘤胃酸中毒（SARA），而本研究中高精料组湖羊结肠内容物的pH值也出现了明显下降，这可能会对结肠内环境的稳定产生不利影响。高精料日粮使湖羊结肠内VFA浓度显著升高，这表明高精料日粮促进了结肠内碳水化合物的发酵。VFA是反刍动物重要的能量来源，其中乙酸、丙酸和丁酸是主要的VFA成分。在本研究中，乙酸、丙酸和丁酸浓度均显著升高，尤其是乙酸浓度增加最为明显。这可能是因为高精料日粮中的淀粉等碳水化合物在结肠内被微生物发酵，首先产生大量的乙酸。有研究表明，VFA可以通过调节肠道上皮细胞的增殖、分化和凋亡，影响肠道的生长和发育。然而，过高的VFA浓度也可能对肠道上皮细胞产生损伤作用，如破坏细胞的紧密连接，增加肠道通透性。氨态氮含量在高精料组显著升高，这说明高精料日粮使湖羊结肠内蛋白质和氨基酸的分解代谢增强。氨态氮是蛋白质和氨基酸分解的产物，其含量的增加可能与高精料日粮中蛋白质含量较高以及结肠微生物对蛋白质的分解能力增强有关。过高的氨态氮含量可能会对结肠黏膜产生刺激，影响肠道的正常功能。在结肠上皮形态方面，高精料日粮导致湖羊结肠绒毛高度显著降低，隐窝深度显著增加，绒毛高度/隐窝深度比值显著降低。绒毛高度的降低会减少肠道上皮的吸收面积，从而影响营养物质的吸收；隐窝深度的增加则可能反映出上皮细胞的更新速度加快，但也可能是肠道受到损伤后的一种代偿性反应。研究表明，绒毛高度/隐窝深度比值是反映肠道消化吸收功能的重要指标，该比值的降低表明高精料日粮损害了湖羊结肠上皮的消化吸收功能。高精料日粮对湖羊结肠黏膜紧密连接蛋白和上皮细胞增殖与凋亡相关基因的表达产生了显著影响。紧密连接蛋白ZO-1和Occludin的基因表达量显著降低，这表明高精料日粮破坏了结肠黏膜紧密连接蛋白的表达，导致肠道屏障功能受损。紧密连接蛋白是维持肠道屏障功能的重要组成部分，其表达的降低会使肠道上皮细胞间的紧密连接受到破坏，增加肠道的通透性，使有害物质更容易进入机体，引发炎症反应。上皮细胞增殖相关基因PCNA的表达量显著降低，表明高精料日粮抑制了湖羊结肠上皮细胞的增殖；凋亡相关基因Bax的表达量显著升高，抗凋亡基因Bcl-2的表达量显著降低，Bax/Bcl-2比值显著升高，表明高精料日粮促进了结肠上皮细胞的凋亡。细胞增殖和凋亡的失衡会影响肠道上皮的正常更新和修复，进一步损害肠道的功能。相关性分析结果表明，VFA浓度与结肠上皮细胞增殖和凋亡相关基因的表达密切相关，VFA浓度的升高会抑制结肠上皮细胞的增殖，促进其凋亡。这可能是因为过高的VFA浓度对结肠上皮细胞产生了毒性作用，导致细胞增殖受到抑制，凋亡增加。高精料日粮抑制了湖羊结肠黏膜VFA吸收与代谢相关基因MCT1和CD36的表达，从而影响了结肠上皮对VFA的吸收和代谢功能。MCT1是一种重要的VFA转运蛋白，负责将VFA从肠腔转运到上皮细胞内；CD36则参与脂肪酸的转运和代谢。它们表达量的降低会减少结肠上皮对VFA的吸收和利用，影响反刍动物的能量供应。高精料日粮对湖羊结肠发酵及上皮形态和功能产生了显著的负面影响。高精料日粮导致结肠内环境酸化，VFA和氨态氮含量升高，损害了结肠上皮的形态结构，破坏了肠道屏障功能，影响了上皮细胞的增殖与凋亡平衡，抑制了VFA吸收与代谢相关基因的表达。这些结果提示，在湖羊养殖中，应合理控制精料的饲喂水平，以维持结肠的健康，保障湖羊的生产性能和机体健康。三、湖羊结肠内容物细菌菌群对高精料日粮的适应性应答3.1材料与方法3.1.1实验动物饲养与样品采集本实验所用湖羊与第二章中高精料日粮对湖羊结肠发酵及上皮形态和功能影响实验中的湖羊为同一批，即选取18只健康状况良好、体重相近（约为[X]kg）的6月龄湖羊公羊，随机分为3组，每组6只羊，分别为对照组（CON组）、实验组1（HG1组）和实验组2（HG2组），并按照相应的日粮配方进行饲养，具体饲养管理和试验设计同2.1.1和2.1.2。在正试期结束后，每组随机选取3只羊进行屠宰，迅速采集结肠内容物样品。将采集到的约10g结肠内容物置于无菌冻存管中，立即放入液氮中速冻，随后转移至-80℃冰箱保存，用于后续的DNA提取及相关分析。3.1.2结肠内容物细菌总DNA提取使用PowerSoilDNAIsolationKit（MoBioLaboratories,Inc.,Carlsbad,CA,USA）试剂盒提取结肠内容物细菌总DNA，具体操作步骤按照试剂盒说明书进行。首先，将1g结肠内容物样品加入到装有陶瓷珠的PowerBeadTubes中，再加入600μLSolutionC1，涡旋振荡30s使样品充分混匀。接着，将PowerBeadTubes置于FastPrep-245G仪器（MPBiomedicals,LLC,Solon,OH,USA）中，以6.0m/s的速度振荡40s，以破碎细菌细胞。然后，将PowerBeadTubes在13,000×g条件下离心5min，将上清液转移至新的离心管中。向上清液中加入200μLSolutionC2，涡旋振荡10s，室温孵育5min。之后，将离心管在13,000×g条件下离心5min，将上清液转移至装有SpinFilter的CollectionTube中，在13,000×g条件下离心1min，弃去滤液。向SpinFilter中加入500μLSolutionC3，在13,000×g条件下离心1min，弃去滤液。再向SpinFilter中加入700μLSolutionC4，在13,000×g条件下离心1min，弃去滤液。重复此步骤一次。最后，将SpinFilter转移至新的离心管中，向SpinFilter中央加入50μLSolutionC6，室温孵育5min，在13,000×g条件下离心1min，收集含有DNA的滤液。提取的DNA用超微量分光光度计（NanoDrop2000，ThermoFisherScientific公司，美国）测定浓度和纯度，要求A260/A280比值在1.8-2.0之间，将合格的DNA样品保存于-20℃冰箱备用。3.1.3PCR扩增及测序以提取的结肠内容物细菌总DNA为模板，对16SrRNA基因的V3-V4可变区进行PCR扩增。PCR扩增引物为338F（5'-ACTCCTACGGGAGGCAGCAG-3'）和806R（5'-GGACTACHVGGGTWTCTAAT-3'），引物两端分别添加了Illumina平台的接头序列和Index序列。PCR反应体系为25μL，包括12.5μL2×TaqMasterMix（VazymeBiotechCo.,Ltd.,Nanjing,China）、1μL上游引物（10μM）、1μL下游引物（10μM）、1μLDNA模板（50-100ng/μL）和9.5μLddH2O。PCR反应条件为：95℃预变性3min；95℃变性30s，55℃退火30s，72℃延伸30s，共30个循环；最后72℃延伸5min。PCR扩增产物通过2%琼脂糖凝胶电泳进行检测，切取目的条带，使用AxyPrepDNAGelExtractionKit（AxygenBiosciences,UnionCity,CA,USA）试剂盒回收纯化。纯化后的PCR产物使用QuantiFluor™-ST蓝色荧光定量系统（PromegaCorporation,Madison,WI,USA）进行定量。将定量后的PCR产物按照等摩尔浓度混合，构建测序文库，使用IlluminaMiSeq测序平台（Illumina,Inc.,SanDiego,CA,USA）进行双端测序（PE300）。3.1.4测序数据分析测序得到的原始数据首先使用Fastp软件（version0.20.0）进行质量控制和过滤。去除低质量的reads（Q值低于20的碱基占比超过10%）、含有接头序列的reads以及长度小于50bp的reads。经过质量控制后的cleanreads使用FLASH软件（version1.2.11）进行拼接，将双端测序得到的reads拼接成完整的序列。拼接后的序列使用Usearch软件（version11.0.667）进行去噪和嵌合体去除，得到高质量的序列。使用Uparse算法将高质量的序列按照97%的相似性进行聚类，得到操作分类单元（OTUs）。利用RDPClassifier（version2.11）对每个OTU的代表性序列进行物种注释，注释数据库为Silva13816SrRNA基因数据库。利用Qiime2（version2020.8）软件计算样品的Alpha多样性指数（包括Chao1、Ace、Shannon和Simpson指数）和Beta多样性指数（包括UnweightedUniFrac和Bray-Curtis距离）。Alpha多样性指数用于评估样品中微生物群落的丰富度和多样性，Chao1和Ace指数反映群落的丰富度，Shannon和Simpson指数反映群落的多样性。Beta多样性指数用于评估不同样品之间微生物群落组成的差异，通过主坐标分析（PCoA）和非度量多维尺度分析（NMDS）对Beta多样性数据进行可视化。使用LEfSe（LineardiscriminantanalysisEffectSize）分析方法，寻找在不同组间具有显著差异的微生物类群（LDAscore>3.0，P<0.05）。3.1.5菌群功能预测利用PICRUSt（PhylogeneticInvestigationofCommunitiesbyReconstructionofUnobservedStates，version1.1.4）软件对结肠内容物细菌菌群的功能进行预测。首先，根据OTUs的代表性序列构建系统发育树，然后利用已知的微生物基因组信息，预测每个OTU的基因功能。将预测得到的基因功能注释到KEGG（KyotoEncyclopediaofGenesandGenomes）数据库中，分析不同组间菌群功能的差异。通过计算不同组间KEGG通路的相对丰度，使用STAMP（StatisticalAnalysisofMetagenomicProfiles，version2.1.3）软件进行统计学分析，筛选出在不同组间具有显著差异的KEGG通路（P<0.05）。3.1.6实时荧光定量PCR测定总菌、厚壁菌门和拟杆菌门基因拷贝数以提取的结肠内容物细菌总DNA为模板，采用实时荧光定量PCR（qPCR）技术测定总菌、厚壁菌门和拟杆菌门的基因拷贝数。总菌的引物为27F（5'-AGAGTTTGATCCTGGCTCAG-3'）和1492R（5'-GGTTACCTTGTTACGACTT-3'）；厚壁菌门的引物为513F（5'-GTGCCAGCMGCCGCGGTAA-3'）和907R（5'-CCGTCAATTCMTTTRAGTTT-3'）；拟杆菌门的引物为555F（5'-GTGYCAGCMGCCGCGGTAA-3'）和907R（5'-CCGTCAATTCMTTTRAGTTT-3'）。qPCR反应体系为20μL，包括10μL2×SYBRGreenPremixExTaqII（TaKaRa公司，日本）、0.8μL上游引物（10μM）、0.8μL下游引物（10μM）、2μLDNA模板（50-100ng/μL）和6.4μLddH2O。qPCR反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3个重复，以10倍梯度稀释的含有目的基因片段的质粒作为标准品，制作标准曲线，根据标准曲线计算样品中目的基因的拷贝数。3.2结果与分析3.2.1结肠内容物菌群丰富度与多样性通过对湖羊结肠内容物16SrRNA基因测序数据的分析，得到各组的Alpha多样性指数，结果如表8所示。Chao1指数和Ace指数用于衡量菌群的丰富度，Shannon指数和Simpson指数用于衡量菌群的多样性。与CON组相比，HG1组和HG2组的Chao1指数和Ace指数均显著降低（P<0.05），且HG2组低于HG1组，差异显著（P<0.05），表明高精料日粮降低了湖羊结肠内容物菌群的丰富度，且精料水平越高，菌群丰富度下降越明显。在多样性方面，HG1组和HG2组的Shannon指数显著降低（P<0.05），Simpson指数显著升高（P<0.05），且HG2组与HG1组之间差异显著（P<0.05），说明高精料日粮使湖羊结肠内容物菌群的多样性降低，菌群结构发生改变。组别Chao1指数Ace指数Shannon指数Simpson指数CON组[X]±[X][X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X][X]±[X]为了进一步分析不同组间结肠内容物菌群组成的差异，进行了Beta多样性分析，并通过主坐标分析（PCoA）进行可视化，结果如图1所示。基于UnweightedUniFrac距离的PCoA分析结果显示，CON组、HG1组和HG2组的样品点在PCoA图上明显分开，表明不同精料水平下湖羊结肠内容物菌群组成存在显著差异。且随着精料水平的升高，菌群组成的差异逐渐增大，说明高精料日粮对湖羊结肠内容物菌群组成产生了明显的影响。注：图中不同颜色的点代表不同组别的样品，PC1和PC2分别代表主坐标1和主坐标2，其贡献率表示该主坐标对数据变异的解释程度。3.2.2结肠内容物细菌菌群组成在门水平上，湖羊结肠内容物中相对丰度排名前10的菌门如图2所示。湖羊结肠内容物中主要的菌门为厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidetes）、变形菌门（Proteobacteria）、放线菌门（Actinobacteria）和疣微菌门（Verrucomicrobia）等。与CON组相比，HG1组和HG2组厚壁菌门的相对丰度显著增加（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；拟杆菌门的相对丰度显著降低（P<0.05），HG2组低于HG1组，差异显著（P<0.05）。变形菌门、放线菌门和疣微菌门的相对丰度在各组间也存在一定差异，但差异不显著（P>0.05）。注：图中不同颜色的柱子代表不同的菌门，横坐标为组别，纵坐标为相对丰度。在属水平上，湖羊结肠内容物中相对丰度排名前10的菌属如图3所示。主要的菌属包括瘤胃球菌属（Ruminococcus）、拟杆菌属（Bacteroides）、丁酸弧菌属（Butyrivibrio）、粪杆菌属（Faecalibacterium）等。与CON组相比，HG1组和HG2组瘤胃球菌属的相对丰度显著增加（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；拟杆菌属的相对丰度显著降低（P<0.05），HG2组低于HG1组，差异显著（P<0.05）。丁酸弧菌属和粪杆菌属的相对丰度在各组间也有变化，但差异不显著（P>0.05）。注：图中不同颜色的柱子代表不同的菌属，横坐标为组别，纵坐标为相对丰度。通过LEfSe分析，进一步筛选出在不同组间具有显著差异的微生物类群（LDAscore>3.0，P<0.05），结果如图4所示。在CON组中，拟杆菌门及其下属的多个菌属（如拟杆菌属、普雷沃氏菌属等）具有较高的相对丰度，是该组的标志性微生物类群；在HG1组中，瘤胃球菌属、丁酸弧菌属等厚壁菌门下属菌属的相对丰度较高；在HG2组中，厚壁菌门的相对丰度进一步增加，其中瘤胃球菌属的相对丰度显著高于其他两组，成为该组的标志性微生物类群。这些结果表明，高精料日粮改变了湖羊结肠内容物中微生物类群的组成，使不同组间具有显著差异的微生物类群发生了变化。注：图中圆圈从内到外依次代表门、纲、目、科、属水平；不同颜色的节点代表在相应组中具有显著差异的微生物类群，节点大小表示相对丰度高低；黄色节点表示在各组间无显著差异的微生物类群。3.2.3结肠内容物菌群功能预测利用PICRUSt软件对湖羊结肠内容物细菌菌群的功能进行预测，将预测得到的基因功能注释到KEGG数据库中，分析不同组间菌群功能的差异。通过计算不同组间KEGG通路的相对丰度，筛选出在不同组间具有显著差异的KEGG通路（P<0.05），结果如表9所示。与CON组相比，HG1组和HG2组中碳水化合物代谢、氨基酸代谢、能量代谢等相关通路的相对丰度发生了显著变化。在碳水化合物代谢通路中，淀粉和蔗糖代谢、果糖和甘露糖代谢等通路的相对丰度显著增加（P<0.05），表明高精料日粮促进了结肠内碳水化合物的代谢。在氨基酸代谢通路中，丙氨酸、天冬氨酸和谷氨酸代谢、甘氨酸、丝氨酸和苏氨酸代谢等通路的相对丰度显著增加（P<0.05），说明高精料日粮也影响了结肠内氨基酸的代谢。在能量代谢通路中，氧化磷酸化通路的相对丰度显著增加（P<0.05），表明高精料日粮可能提高了结肠内的能量代谢水平。KEGG通路CON组HG1组HG2组P值淀粉和蔗糖代谢[X]±[X][X]±[X][X]±[X][X]果糖和甘露糖代谢[X]±[X][X]±[X][X]±[X][X]丙氨酸、天冬氨酸和谷氨酸代谢[X]±[X][X]±[X][X]±[X][X]甘氨酸、丝氨酸和苏氨酸代谢[X]±[X][X]±[X][X]±[X][X]氧化磷酸化[X]±[X][X]±[X][X]±[X][X]3.2.4结肠内容物中总菌、厚壁菌门和拟杆菌门数量采用实时荧光定量PCR（qPCR）技术测定湖羊结肠内容物中总菌、厚壁菌门和拟杆菌门的基因拷贝数，结果如表10所示。与CON组相比，HG1组和HG2组结肠内容物中总菌和厚壁菌门的基因拷贝数显著增加（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；拟杆菌门的基因拷贝数显著降低（P<0.05），HG2组低于HG1组，差异显著（P<0.05）。这与16SrRNA基因测序分析中厚壁菌门和拟杆菌门相对丰度的变化趋势一致，进一步表明高精料日粮改变了湖羊结肠内容物中总菌以及厚壁菌门和拟杆菌门的数量。组别总菌（log10拷贝数/g）厚壁菌门（log10拷贝数/g）拟杆菌门（log10拷贝数/g）CON组[X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X]3.3讨论本研究结果表明，高精料日粮显著改变了湖羊结肠内容物细菌菌群的丰富度、多样性和群落组成。随着精料水平的升高，湖羊结肠内容物菌群的丰富度和多样性显著降低，这与之前的相关研究结果一致。高精料日粮中富含大量易发酵的碳水化合物，会导致结肠内环境发生改变，如pH值降低、VFA浓度升高等，这些变化可能不利于一些微生物的生长和生存，从而导致菌群丰富度和多样性下降。研究发现，当瘤胃内pH值低于5.5时，一些对酸性环境敏感的微生物会受到抑制，其数量和种类都会减少。本研究中高精料组湖羊结肠内容物的pH值显著降低，可能是导致菌群丰富度和多样性下降的重要原因之一。在菌群组成方面，高精料日粮使湖羊结肠内容物中厚壁菌门的相对丰度显著增加，拟杆菌门的相对丰度显著降低。厚壁菌门和拟杆菌门是反刍动物肠道内的主要菌门，它们在肠道内发挥着不同的功能。厚壁菌门具有较强的碳水化合物发酵能力，能够利用淀粉等碳水化合物产生大量的VFA。在本研究中，高精料日粮中富含淀粉等碳水化合物，为厚壁菌门的生长提供了丰富的底物，从而使其相对丰度显著增加。拟杆菌门则主要参与纤维素等多糖的降解，同时在维持肠道屏障功能和免疫调节方面发挥着重要作用。高精料日粮中纤维素含量较低，可能导致拟杆菌门的生长受到限制，其相对丰度降低。拟杆菌门相对丰度的降低可能会影响肠道屏障功能和免疫调节，增加肠道疾病的发生风险。瘤胃球菌属作为厚壁菌门下属的重要菌属，在高精料组湖羊结肠内容物中的相对丰度显著增加。瘤胃球菌属具有较强的淀粉降解能力，能够将淀粉分解为葡萄糖等小分子物质，进而发酵产生VFA。在高精料日粮条件下，瘤胃球菌属能够更好地利用丰富的淀粉底物，大量繁殖，导致其相对丰度显著升高。拟杆菌属作为拟杆菌门的代表菌属，其相对丰度在高精料组显著降低。这可能是由于高精料日粮改变了结肠内的生态环境，不利于拟杆菌属的生长和繁殖。拟杆菌属相对丰度的降低可能会削弱肠道对纤维素等多糖的降解能力，影响肠道的正常消化功能。通过LEfSe分析筛选出的不同组间具有显著差异的微生物类群，进一步表明高精料日粮改变了湖羊结肠内容物中微生物类群的组成。在CON组中，拟杆菌门及其下属的多个菌属具有较高的相对丰度，这与该组低精料水平、高纤维含量的日粮组成有关，这些微生物能够适应富含纤维素的环境，发挥对纤维素的降解作用。在HG1组和HG2组中，厚壁菌门下属的瘤胃球菌属、丁酸弧菌属等菌属的相对丰度较高，这与高精料日粮中碳水化合物含量高，促进了这些具有碳水化合物发酵能力的菌属的生长有关。瘤胃球菌属在HG2组中的相对丰度显著高于其他两组，成为该组的标志性微生物类群，这进一步说明了随着精料水平的升高，瘤胃球菌属在结肠内容物菌群中的优势地位更加明显。菌群功能预测结果显示，高精料日粮使湖羊结肠内容物细菌菌群在碳水化合物代谢、氨基酸代谢、能量代谢等相关通路的相对丰度发生了显著变化。在碳水化合物代谢通路中，淀粉和蔗糖代谢、果糖和甘露糖代谢等通路的相对丰度显著增加，这与高精料日粮中富含碳水化合物，刺激了结肠内碳水化合物代谢相关菌群的生长和代谢活动有关。这些菌群能够利用碳水化合物产生更多的VFA，为湖羊提供能量。在氨基酸代谢通路中，丙氨酸、天冬氨酸和谷氨酸代谢、甘氨酸、丝氨酸和苏氨酸代谢等通路的相对丰度显著增加，这可能是由于高精料日粮中蛋白质和氨基酸的含量较高，或者是由于碳水化合物代谢的改变影响了氨基酸的代谢。在能量代谢通路中，氧化磷酸化通路的相对丰度显著增加，表明高精料日粮可能提高了结肠内的能量代谢水平，这可能与碳水化合物和氨基酸代谢的增强，产生了更多的能量底物有关。实时荧光定量PCR结果表明，高精料日粮改变了湖羊结肠内容物中总菌以及厚壁菌门和拟杆菌门的数量，这与16SrRNA基因测序分析中菌群相对丰度的变化趋势一致。总菌和厚壁菌门的基因拷贝数显著增加，拟杆菌门的基因拷贝数显著降低，进一步证明了高精料日粮对湖羊结肠内容物细菌菌群数量的影响。厚壁菌门基因拷贝数的增加，反映了其在高精料日粮条件下的生长优势；而拟杆菌门基因拷贝数的降低，则表明其生长受到了抑制。高精料日粮对湖羊结肠内容物细菌菌群产生了显著影响，改变了菌群的丰富度、多样性、群落组成和功能。这些变化可能与高精料日粮导致的结肠内环境改变有关，如pH值降低、VFA浓度升高、底物组成改变等。菌群的变化又会进一步影响结肠的发酵功能、营养物质代谢以及肠道健康。在实际湖羊养殖中，应充分考虑高精料日粮对结肠微生物区系的影响，合理调整日粮配方，以维持结肠微生物的平衡，保障湖羊的健康和生产性能。四、湖羊结肠黏膜微生物区系和炎症因子表达对高精料日粮的适应性应答4.1材料与方法4.1.1实验动物饲养与样品采集本实验所用湖羊与前文研究为同一批，即18只健康状况良好、体重相近（约为[X]kg）的6月龄湖羊公羊，随机分为3组，每组6只羊，分别为对照组（CON组）、实验组1（HG1组）和实验组2（HG2组），并按照相应的日粮配方进行饲养，具体饲养管理和试验设计同2.1.1和2.1.2。在正试期结束后，每组随机选取3只羊进行屠宰，迅速采集结肠黏膜样品。用无菌镊子轻轻刮取结肠黏膜表面的黏液和附着的微生物，将采集到的黏膜样品置于无菌冻存管中，立即放入液氮中速冻，随后转移至-80℃冰箱保存，用于后续的DNA提取及相关分析。同时，采集部分结肠黏膜组织，用于提取RNA，检测炎症因子相关基因的表达水平。4.1.2结肠黏膜细菌总DNA提取使用QIAampDNAStoolMiniKit（Qiagen公司，德国）试剂盒提取结肠黏膜细菌总DNA，具体操作步骤按照试剂盒说明书进行。首先，取约200mg结肠黏膜样品加入到装有缓冲液ASL的离心管中，涡旋振荡30s使样品充分混匀。然后，将离心管在95℃孵育10min，以裂解细菌细胞。接着，将离心管在13,000×g条件下离心5min，将上清液转移至新的离心管中。向上清液中加入10μL蛋白酶K和200μL缓冲液AL，涡旋振荡10s，室温孵育10min。之后，加入200μL无水乙醇，涡旋振荡15s，将混合液转移至QIAampMiniSpinColumn中，在13,000×g条件下离心1min，弃去滤液。向QIAampMiniSpinColumn中加入500μL缓冲液AW1，在13,000×g条件下离心1min，弃去滤液。再向QIAampMiniSpinColumn中加入500μL缓冲液AW2，在13,000×g条件下离心3min，弃去滤液。将QIAampMiniSpinColumn转移至新的离心管中，向SpinColumn中央加入50μL缓冲液AE，室温孵育5min，在13,000×g条件下离心1min，收集含有DNA的滤液。提取的DNA用超微量分光光度计（NanoDrop2000，ThermoFisherScientific公司，美国）测定浓度和纯度，要求A260/A280比值在1.8-2.0之间，将合格的DNA样品保存于-20℃冰箱备用。4.1.3PCR扩增及测序以提取的结肠黏膜细菌总DNA为模板，对16SrRNA基因的V3-V4可变区进行PCR扩增。PCR扩增引物为338F（5'-ACTCCTACGGGAGGCAGCAG-3'）和806R（5'-GGACTACHVGGGTWTCTAAT-3'），引物两端分别添加了Illumina平台的接头序列和Index序列。PCR反应体系为25μL，包括12.5μL2×TaqMasterMix（VazymeBiotechCo.,Ltd.,Nanjing,China）、1μL上游引物（10μM）、1μL下游引物（10μM）、1μLDNA模板（50-100ng/μL）和9.5μLddH2O。PCR反应条件为：95℃预变性3min；95℃变性30s，55℃退火30s，72℃延伸30s，共30个循环；最后72℃延伸5min。PCR扩增产物通过2%琼脂糖凝胶电泳进行检测，切取目的条带，使用AxyPrepDNAGelExtractionKit（AxygenBiosciences,UnionCity,CA,USA）试剂盒回收纯化。纯化后的PCR产物使用QuantiFluor™-ST蓝色荧光定量系统（PromegaCorporation,Madison,WI,USA）进行定量。将定量后的PCR产物按照等摩尔浓度混合，构建测序文库，使用IlluminaMiSeq测序平台（Illumina,Inc.,SanDiego,CA,USA）进行双端测序（PE300）。4.1.4测序数据分析测序得到的原始数据首先使用Fastp软件（version0.20.0）进行质量控制和过滤。去除低质量的reads（Q值低于20的碱基占比超过10%）、含有接头序列的reads以及长度小于50bp的reads。经过质量控制后的cleanreads使用FLASH软件（version1.2.11）进行拼接，将双端测序得到的reads拼接成完整的序列。拼接后的序列使用Usearch软件（version11.0.667）进行去噪和嵌合体去除，得到高质量的序列。使用Uparse算法将高质量的序列按照97%的相似性进行聚类，得到操作分类单元（OTUs）。利用RDPClassifier（version2.11）对每个OTU的代表性序列进行物种注释，注释数据库为Silva13816SrRNA基因数据库。利用Qiime2（version2020.8）软件计算样品的Alpha多样性指数（包括Chao1、Ace、Shannon和Simpson指数）和Beta多样性指数（包括UnweightedUniFrac和Bray-Curtis距离）。Alpha多样性指数用于评估样品中微生物群落的丰富度和多样性，Chao1和Ace指数反映群落的丰富度，Shannon和Simpson指数反映群落的多样性。Beta多样性指数用于评估不同样品之间微生物群落组成的差异，通过主坐标分析（PCoA）和非度量多维尺度分析（NMDS）对Beta多样性数据进行可视化。使用LEfSe（LineardiscriminantanalysisEffectSize）分析方法，寻找在不同组间具有显著差异的微生物类群（LDAscore>3.0，P<0.05）。4.1.5菌群功能预测利用PICRUSt（PhylogeneticInvestigationofCommunitiesbyReconstructionofUnobservedStates，version1.1.4）软件对结肠黏膜细菌菌群的功能进行预测。首先，根据OTUs的代表性序列构建系统发育树，然后利用已知的微生物基因组信息，预测每个OTU的基因功能。将预测得到的基因功能注释到KEGG（KyotoEncyclopediaofGenesandGenomes）数据库中，分析不同组间菌群功能的差异。通过计算不同组间KEGG通路的相对丰度，使用STAMP（StatisticalAnalysisofMetagenomicProfiles，version2.1.3）软件进行统计学分析，筛选出在不同组间具有显著差异的KEGG通路（P<0.05）。4.1.6实时荧光定量PCR测定总菌、厚壁菌门和拟杆菌门基因拷贝数以提取的结肠黏膜细菌总DNA为模板，采用实时荧光定量PCR（qPCR）技术测定总菌、厚壁菌门和拟杆菌门的基因拷贝数。总菌的引物为27F（5'-AGAGTTTGATCCTGGCTCAG-3'）和1492R（5'-GGTTACCTTGTTACGACTT-3'）；厚壁菌门的引物为513F（5'-GTGCCAGCMGCCGCGGTAA-3'）和907R（5'-CCGTCAATTCMTTTRAGTTT-3'）；拟杆菌门的引物为555F（5'-GTGYCAGCMGCCGCGGTAA-3'）和907R（5'-CCGTCAATTCMTTTRAGTTT-3'）。qPCR反应体系为20μL，包括10μL2×SYBRGreenPremixExTaqII（TaKaRa公司，日本）、0.8μL上游引物（10μM）、0.8μL下游引物（10μM）、2μLDNA模板（50-100ng/μL）和6.4μLddH2O。qPCR反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3个重复，以10倍梯度稀释的含有目的基因片段的质粒作为标准品，制作标准曲线，根据标准曲线计算样品中目的基因的拷贝数。4.1.7结肠黏膜RNA提取及RT-qPCR使用Trizol试剂（Invitrogen公司，美国）提取结肠黏膜组织的总RNA，按照试剂盒说明书的操作步骤进行。提取的RNA用超微量分光光度计（NanoDrop2000，ThermoFisherScientific公司，美国）测定浓度和纯度，要求A260/A280比值在1.8-2.0之间。取1μg总RNA，利用反转录试剂盒（PrimeScriptRTreagentKitwithgDNAEraser，TaKaRa公司，日本）将其反转录为cDNA。以cDNA为模板，采用SYBRGreen荧光定量PCR试剂盒（TBGreenPremixExTaqII，TaKaRa公司，日本）在实时荧光定量PCR仪（CFX96Touch，Bio-Rad公司，美国）上进行qPCR反应。反应体系为20μL，包括10μLSYBRGreenPremixExTaqII、0.8μL上游引物、0.8μL下游引物、2μLcDNA模板和6.4μLddH2O。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3个重复，以β-actin作为内参基因，采用2-ΔΔCt法计算目的基因的相对表达量。本实验检测的炎症因子相关基因包括白细胞介素-1β（IL-1β）、白细胞介素-6（IL-6）、肿瘤坏死因子-α（TNF-α）等，引物序列如表11所示：基因名称引物序列（5'-3'）IL-1βF:[X]R:[X]IL-6F:[X]R:[X]TNF-αF:[X]R:[X]β-actinF:[X]R:[X]4.2结果与分析4.2.1结肠黏膜微生物菌群的丰度与多样性通过对湖羊结肠黏膜16SrRNA基因测序数据的分析，得到各组的Alpha多样性指数，结果如表12所示。Chao1指数和Ace指数反映菌群的丰富度，Shannon指数和Simpson指数反映菌群的多样性。与CON组相比，HG1组和HG2组的Chao1指数和Ace指数均显著降低（P<0.05），且HG2组低于HG1组，差异显著（P<0.05），表明高精料日粮降低了湖羊结肠黏膜菌群的丰富度，且精料水平越高，菌群丰富度下降越明显。在多样性方面，HG1组和HG2组的Shannon指数显著降低（P<0.05），Simpson指数显著升高（P<0.05），且HG2组与HG1组之间差异显著（P<0.05），说明高精料日粮使湖羊结肠黏膜菌群的多样性降低，菌群结构发生改变。组别Chao1指数Ace指数Shannon指数Simpson指数CON组[X]±[X][X]±[X][X]±[X][X]±[X]HG1组[X]±[X][X]±[X][X]±[X][X]±[X]HG2组[X]±[X][X]±[X][X]±[X][X]±[X]为进一步分析不同组间结肠黏膜微生物菌群组成的差异，进行了Beta多样性分析，并通过主坐标分析（PCoA）进行可视化，结果如图5所示。基于UnweightedUniFrac距离的PCoA分析结果显示，CON组、HG1组和HG2组的样品点在PCoA图上明显分开，表明不同精料水平下湖羊结肠黏膜微生物菌群组成存在显著差异。且随着精料水平的升高，菌群组成的差异逐渐增大，说明高精料日粮对湖羊结肠黏膜微生物菌群组成产生了明显的影响。注：图中不同颜色的点代表不同组别的样品，PC1和PC2分别代表主坐标1和主坐标2，其贡献率表示该主坐标对数据变异的解释程度。4.2.2结肠黏膜细菌菌群组成在门水平上，湖羊结肠黏膜中相对丰度排名前10的菌门如图6所示。湖羊结肠黏膜中主要的菌门为厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidetes）、变形菌门（Proteobacteria）、放线菌门（Actinobacteria）和疣微菌门（Verrucomicrobia）等。与CON组相比，HG1组和HG2组厚壁菌门的相对丰度显著增加（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；拟杆菌门的相对丰度显著降低（P<0.05），HG2组低于HG1组，差异显著（P<0.05）。变形菌门、放线菌门和疣微菌门的相对丰度在各组间也存在一定差异，但差异不显著（P>0.05）。注：图中不同颜色的柱子代表不同的菌门，横坐标为组别，纵坐标为相对丰度。在属水平上，湖羊结肠黏膜中相对丰度排名前10的菌属如图7所示。主要的菌属包括瘤胃球菌属（Ruminococcus）、拟杆菌属（Bacteroides）、丁酸弧菌属（Butyrivibrio）、粪杆菌属（Faecalibacterium）等。与CON组相比，HG1组和HG2组瘤胃球菌属的相对丰度显著增加（P<0.05），且HG2组高于HG1组，差异显著（P<0.05）；拟杆菌属的相对丰度显著降低（P<0.05），HG2组低于HG1组，差异显著（P<0.05）。丁酸弧菌属和粪杆菌属的相对丰度在各组间也有变化，但差异不显著（P>0.05）。注：图中不同颜色的柱子代表不同的菌属，横坐标为组别，纵坐标为相对丰度。通过