网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

冻融对燕麦青贮有氧暴露阶段细菌和真菌演替的影响  PDF

  • 李海萍 1,2
  • 贾志锋 3
  • 刘文辉 3
  • 马祥 3
  • 周青平 2
  • 关皓 2
1. 青海师范大学 生命科学学院,青海 西宁; 2. 西南民族大学,四川若尔盖高寒湿地生态系统国家野外科学观测研究站,四川 成都; 3. 青海大学,青海省青藏高原优良牧草种质资源利用重点实验室,青海 西宁

最近更新:2025-03-07

DOI: 10.13343/j.cnki.wsxb.20240523

CSTR: 32112.14.j.AMS.20240523

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

目的

为探索冻融条件对燕麦青贮在有氧暴露阶段的影响。

方法

将燕麦在冻融条件(20 ℃和-5 ℃每12 h交替1次,S组)下青贮60 d后进行有氧暴露监测,并以恒温20 ℃ (20组)作为对照。分别在青贮60 d和有氧暴露1、3、5 d取样进行发酵品质和营养品质的测定,同时进行16S rRNA基因和ITS测序。

结果

随着有氧暴露时间的延长pH值升高。与原冻融温度下的有氧暴露相比,室温下的有氧暴露导致pH值和氨态氮含量上升更快,乳酸和乙酸含量下降也更快。S组,特别是室温下的有氧暴露条件下,丙酸和丁酸的生成速度更快。随着有氧暴露的推进,肠杆菌和酵母菌数量不断增加,而乳酸菌数量不断减少(P<0.05)。细菌的Shannon指数和Simpson指数随着有氧暴露时间的延长而增加,乳酸菌的丰度不断降低。有氧暴露第3天时,细菌群落结构出现明显差异。与原温度下的有氧暴露相比,室温下的有氧暴露加速了微生物的更替。S组中指征腐败的酵母菌和霉菌的种类及数量更多。

结论

冻融条件加速了燕麦青贮在有氧暴露阶段的变败过程,尤其是在室温下进行有氧暴露时。本研究为高寒区燕麦青贮的高品质调制及存放提供理论指导。

发酵后的出料标志着青贮过程的结束。除了良好的发酵品质,有氧稳定性也是确保饲喂效果的关键因素。通过研究冻融对燕麦青贮发酵过程的影响,发现冻融不仅作用于发酵阶段的青贮品质,而且在开窖后对青贮的有氧稳定性也起着至关重要的作[

1]。在青贮发酵过程中,产生的乳酸为微生物提供了生长所需的底物,因此,当青贮窖(或青贮包)被打开时,发酵好的青贮料再次与氧气接[2],导致一些好氧腐败微生物(如酵母等)变得活跃并大量繁[3],从而引发二次发酵。这不仅会造成干物质的大量损[4]和营养物质的严重流[5],还会导致青贮料温度升高,进一步加剧青贮变质。此外,有氧腐败还会增加潜在致病性微生物或其他不良微生物(如霉菌、肠杆菌和单核增生李斯特菌)的增殖风险,对青贮饲料的卫生质量造成严重的负面影[6],进而影响饲喂效果。

冻融研究是发酵过程研究的重要组成部分。冻融引发的应激是一个复杂的过程,可能包括脱水作用、物理应激、低温应激和氧化应激[

7]。对土壤微生物的冻融影响研究发现,冻融不会导致土壤微生物α多样性的减少,但会显著影响微生物的生物量和群落结构,导致微生物C和N代谢功能下[8]。Dong[9]通过研究红三叶经受冻融后再进行青贮的情况,发现原料的冻融促进了青贮微生物区系中乳酸菌的生长,降低了发酵期间泛菌和肠杆菌的相对丰度,对发酵品质有一定的促进作用;然而,不动杆菌和假单胞菌在青贮末期再次出现,影响了青贮品质和有氧稳定性。随着分子技术的快速发展,高通量测序技术为深入了解青贮微生物群落变化提供了技术支撑。然而,近年来,研究多集中于对青贮发酵过程的分析,而对于有氧暴露阶段的关注则相对较少。

为全面了解冻融对燕麦青贮有氧稳定性的影响,本研究在燕麦青贮60 d后进行有氧暴露试验,旨在探究其在有氧暴露阶段的品质变化及微生物演替情况。

1 材料与方法

1.1 青贮调制

试验材料选用乳熟期的燕麦,留茬高度为10 cm。刈割后的燕麦使用揉丝机切割并揉丝成1-2 cm长度,经短暂晾晒后充分混匀。在所有材料充分混合后,随机取约700 g的材料装入无菌抽真空袋(直径×高度:20 cm×30 cm),并进行抽真空处理。试验设计采用双因素完全随机方式,温度分别是恒温20 ℃ (20组)和冻融条件(S组,在20 ℃和-5 ℃每12 h交替1次,即20 ℃/-5 ℃,温度切换时以每分钟0.3 ℃的速率变化直到稳定),取样时间点设定为青贮后的60 d,以及有氧暴露后的1、3、5 d,每个处理设置3个重复。

1.2 有氧稳定性监测

在青贮60 d后,从每个处理中均匀取出1 kg的样品放入2 L的灭菌烧杯中压实,并用两层灭菌棉纱布覆盖。将样品分别在其原始处理温度和室温(25±2) ℃条件下放置5 d。使用实时温度记录仪(神华科技股份有限公司)每5 min测量1次青贮核心区域(深度为10 cm)的温度,持续检测。分别在有氧暴露后的1、3、5 d分别取样,进行微生物计数和发酵品质分析。有氧稳定性的评估基于暴露于空气中的青贮饲料的温度超过基准环境温度2 °C所持续的时[

10]

1.3 青贮发酵品质和营养成分分析

样品于65 ℃烘箱中烘干至恒重,以测定其干物质含量(dry matter, DM)。粉碎后过0.6 mm筛用于测定其他营养成分。粗蛋白含量(crude protein, CP)采用杜马斯燃烧定氮法测[

11]。可溶性糖含量使用植物可溶性糖(water-soluble carbohydrates, WSC)含量试剂盒(苏州科铭生物技术有限公司)进行测定。中性洗涤纤维(neutral detergent fiber, NDF)和酸性洗涤纤维(acid detergent fiber, ADF)依据Van Soest[12]方法进行测定。取20 g样品置于180 mL无菌蒸馏水中浸提24 h后过滤,一部分滤液用pH计(上海仪电分析仪器有限公司)测量pH值,剩余滤液分别用于氨态氮和有机酸的测定。氨态氮采用苯酚-次氯酸钠比色法进行测[13]。有机酸则使用超高效液相色谱仪进行测定,滤液通过0.22 μm滤头过滤至进样瓶中,采用赛默飞UltiMate 3000超高效液相色谱仪测定乳酸、乙酸、丙酸和丁酸。色谱条件为:RSpak KC-811色谱柱,流动相0.1% H3PO4,流速0.5 mL/min,柱温55 ℃。

1.4 微生物计数

青贮微生物计数采用平板计数[

14]。取20 g样品浸没于180 mL灭菌生理盐水(0.85%)中,在4 ℃、120×g振荡1 h,然后用灭菌蒸馏水依次稀释至10-5,从中取原液、10-3、10-5这3个梯度各0.1 mL进行涂布。分别使用MRS培养基37 ℃厌氧培养48 h对乳酸菌进行计数,使用马铃薯葡萄糖琼脂(PDA)培养基30 ℃有氧培养72 h对酵母菌和霉菌进行计数,使用结晶紫中性红胆盐琼脂(VRBA)培养基在37 ℃有氧培养24 h对肠杆菌进行计数,每个梯度每种培养基设置3个重复。

1.5 微生物测序

将20 g的青贮样品中加入180 mL灭菌蒸馏水,在4 ℃振荡30 min,通过两层灭菌纱布过滤后4 ℃、10 000×g离心10 min,倒掉上清液,加入3 mL灭菌蒸馏水重制成混合菌液。采用CTAB方法对样品的基因组DNA进行提取,之后使用琼脂糖凝胶电泳检测提取出的DNA纯度和浓度。随后取适量DNA于离心管中,并用无菌水稀释至1 ng/μL。以稀释后的基因组DNA为模板,根据扩增区域选择使用带barcode的特异性引物、New England Biolabs公司的Phusion® High-Fidelity PCR Master Mix with GC Buffer和高效高保真酶进行PCR,以确保扩增效率和准确性。细菌序列使用特异性引物515F (5′-GTTTCGGTGCCAGCMGCCGCGGTAA-3′)和806R (5′-GCCAATGGACTACHVGGGTWTCTAAT- 3′)扩增16S rRNA基因的V4区;真菌序列使用特异性引物ITS5-1737F (5′-GGAAGTAAAAGT CGTAACAAGG-3′)和ITS2-2043R (5′-GCTGC GTTCTTCATCGATGC-3′)对样品真菌基因组DNA进行扩增。PCR反应体系:PCR混合液加入15 µL Phusion® High-Fidelity PCR Master Mix(New England Biolabs公司)、0.2 µmol/L引物和10 ng基因组DNA模板。反应条件:98 ℃预变性1 min;98 ℃变性10 s,50 ℃退火30 s,72 ℃延伸30 s,共30次循环;72 ℃保持5 min。PCR产物使用2%浓度的琼脂糖凝胶进行电泳检测;对检测合格的PCR产物进行磁珠纯化,并采用酶标定量。根据PCR产物浓度进行等量混样,充分混匀后使用2%的琼脂糖凝胶电泳检测PCR产物。对目的条带使用胶回收试剂盒(Qiagen公司)回收产物。使用NEBNext® Ultra™ II DNA Library Prep Kit进行文库构建,将构建好的文库进行Qubit和qPCR定量。待文库合格后,使用NovaSeq 6000进行上机测序,并进行后续的数据分析。使用FLASH (v1.2.7)对每个样本的reads进行拼接。参照QIIME (v1.9.1)的Tags质量控制流程,得到最终的有效序列。利用Uparse算法(v7.0.1001)对所有样本的全部有效序列进行聚类,默认以97%的一致性将序列聚类成为OTUs。使用MUSCLE (3.8.31)软件进行快速多序列比对,得到所有OTUs代表序列的系统发生关系。使用QIIME软件(v1.9.1)计算α和β多样性指数,并用R软件(v2.15.3)进行可视化展示。微生物测序数据已上传至NCBI数据库(https://www.ncbi.nlm.nih.gov/),登录号为PRJNA1184690。

采用SPSS (IBM SPSS Statistics 26)软件进行双因素方差分析,旨在探讨温度和青贮阶段及其交互作用对有氧暴露阶段所有指标的影响,并应用Tukey’s HSD检验方法确定差异显著性水平(P<0.05)。

2 结果与分析

2.1 发酵品质分析

青贮60 d及开袋后有氧暴露的发酵特性如表1所示,随着有氧暴露时间的延长pH值升高。在室温下有氧暴露1 d后,S组的pH值升高至7.10,5 d后最终超过了8.00。在整个有氧暴露过程中,S组在室温下的pH值始终显著高于同一时段在原温度下的pH值(P<0.05)。20组除第1天外,第3天和第5天的pH值也呈现出同样的趋势。NH3-N含量不断上升,而乳酸含量快速下降,特别是S组。乙酸含量先上升后下降,而丙酸和丁酸在S组有氧暴露后被检测到,并且随着有氧暴露时间的延长其含量逐渐增加。与原温度下的有氧暴露相比,室温下的有氧暴露导致pH值和氨态氮含量上升更快,乳酸和乙酸含量下降更快。在S组,尤特别是在室温下有氧暴露时,丙酸和丁酸的生成速度更快。乳酸和乙酸与时间段、温度以及两者的交互作用均存在极显著的相关性(P<0.01)。pH值和丁酸与时间段、温度也呈现极显著的相关性(P<0.01),但与两者的交互作用无显著相关性。

表1  燕麦青贮60 d和有氧暴露过程中的发酵特性
Table 1  Fermentation characteristics of oat silage ensiled for 60 days and during aerobic exposure
时间Time温度T/pH

氨态氮

NH3-N (%)

乳酸

Lactic acid

(mg/g)

乙酸

Acetic acid

(mg/g)

丙酸

Propionic acid

(mg/g)

丁酸

Butyric acid

(mg/g)

Day 60 20 4.62k 3.50f 48.03a 11.87c ND ND
20/-5 4.92h 4.24de 28.06c 12.74c ND ND
AE Day 1 AE 20 4.67hi 3.88ef 41.36b 18.26a ND ND
AE 20/-5 5.24g 4.19de 19.95d 15.62b ND ND
RT AE 20 4.73hi 3.86ef 18.44e 18.64a ND ND
RT AE 20/-5 7.10e 4.25de 5.75h 14.71b 0.44d 1.31d
AE Day 3 AE 20 6.71f 4.25de 13.90f 13.08c ND ND
AE 20/-5 7.76d 4.43de 3.05jk 9.84d ND ND
RT AE 20 8.03c 4.61d 3.48i 12.19c ND ND
RT AE 20/-5 8.46b 4.68d 3.99i 7.18e 2.44b 3.11b
AE Day 5 AE 20 7.71d 6.29c 7.68g 9.84d ND ND
AE 20/-5 8.09c 6.30c 2.17kl 6.44ef ND 1.15e
RT AE 20 8.45b 7.61b 1.89l 9.22d 1.35c 2.13c
RT AE 20/-5 8.81a 8.58a 1.65l 5.48f 2.71a 4.40a
SEM 0.091 0.103 1.003 0.166 0.133 0.175
Period ** * ** ** * **
Temperature ** ** ** ** NS **
Period×Temperature NS NS ** ** NS NS

同列不同小写字母表示在0.05水平上差异显著(*:P<0.05;**:P<0.01);NS:不显著;SEM:标准误;DM:干物质;AE:有氧暴露;RT:室温;ND:未检测到。

The same column followed by different lowercase letters are significantly different at the 0.05 level (*: P<0.05; **: P<0.01); NS: Not significant; SEM: Standard error of mean; DM: Dry matter; AE: Aerobic exposure; RT: Room temperature; ND: Not detected.

2.2 营养成分分析

表2所示,有氧暴露5 d后,20组和S组的灰分均显著升高(P<0.05),同时中性洗涤纤维的含量也显著增加。然而,在有氧暴露5 d后,20组和S组之间的中性洗涤纤维和酸性洗涤纤维含量无著差异。粗蛋白的含量显著下降(P<0.05),但不同处理间无显著差异。S组的可溶性糖含量在有氧暴露5 d后显著下降(P<0.05),并且在室温下进行有氧暴露的样品中,可溶性糖含量的下降速度更快(P<0.05)。

表2  燕麦青贮60 d和有氧暴露5 d的营养品质
Table 2  Nutritional components in silage ensiled for 60 days and aerobic exposure for 5 days
时间Time温度T/灰分Ash(%)中性洗涤纤维Neutral detergent fiber (%)酸性洗涤纤维Acid detergent fiber (%)

粗蛋白

Crude protein (%)

可溶性糖Water-soluble carbohydrates (%)
Day 60 20 10.14d 56.98b 35.97b 8.67a 4.96a
20/-5 11.26c 59.03b 40.81ab 8.93a 4.84ab
AE Day 5 AE 20 11.91b 69.64a 43.85a 7.54b 4.33ab
AE 20/-5 11.45c 67.83a 41.95ab 7.97b 4.24ab
RT AE 20 12.79a 70.25a 41.01ab 7.82b 4.11b
RT AE 20/-5 12.47a 67.22a 37.26ab 7.62b 3.33c
SEM 0.105 1.939 2.157 0.062 0.096

同列不同小写字母表示在0.05水平上差异显著;SEM:标准误;DM:干物质;AE:有氧暴露;RT:室温。

The same column followed by different lowercase letters are significantly different at the 0.05 level. SEM: Standard error of mean; DM: Dry matter; AE: Aerobic exposure; RT: Room temperature.

2.3 微生物计数

表3所示,在青贮60 d开袋时,S组中检测到了肠杆菌,而20组和S组中均未检测到酵母菌。随着有氧暴露时间的推移,肠杆菌和酵母菌的数量不断增加,乳酸菌的数量则不断减少(P<0.05)。在S组有氧暴露3 d后检测到霉菌,而20组在有氧暴露5 d后才检测到霉菌,S组在室温下进行有氧暴露时,霉菌的数量显著多于其他处理(P<0.05)。有氧暴露5 d后,各处理下肠杆菌的计数差异显著(P<0.05)。肠杆菌的数量与时间阶段和温度均呈极显著相关性,但与两者的交互作用无相关性。乳酸菌的计数与时间段、时间段与温度的交互作用极显著相关(P<0.01),与温度也呈显著相关性(P<0.05)。酵母菌的计数与时间阶段和时间阶段与温度的交互作用极显著相关(P<0.01),但与温度无显著相关性。霉菌的计数与时间段、温度及两者的交互作用均呈极显著相关性(P<0.01)。

表3  燕麦青贮60 d和有氧暴露过程中的微生物计数
Table 3  Microorganism counts of oat silage ensiled for 60 days and during aerobic exposure (lg CFU/g)
时间Time

温度

T/

肠杆菌

Coliform bacteria

乳酸菌

Lactic acid bacteria

酵母菌

Yeasts

霉菌

Molds

Day 60 20 ND 5.04de ND ND
20/-5 4.42ef 5.37c ND ND
AE Day 1 AE 20 ND 5.41c 3.82f ND
AE 20/-5 4.25f 6.01a 4.18ef ND
RT AE 20 ND 4.32g 3.00g ND
RT AE 20/-5 4.3ef 5.27cd 4.42de ND
AE Day 3 AE 20 ND 5.39c 5.31ab ND
AE 20/-5 6.21a 5.77ab 4.82c 2.62d
RT AE 20 4.53d 5.66b 5.46ab ND
RT AE 20/-5 5.12c 5.76ab 4.68cd 3.93b
AE Day 5 AE 20 2.15g 4.96ef 5.06bc 2.54de
AE 20/-5 6.24a 4.79ef 4.89c 2.97c
RT AE 20 4.97c 4.71f 5.69a 2.43e
RT AE 20/-5 5.61b 4.35g 5.45ab 4.71a
SEM 0.187 0.053 0.053 0.075
Period ** ** ** **
Temperature ** * NS **
Period×Temperature NS ** ** **

同列不同小写字母表示在0.05水平上差异显著(*:P<0.05;**:P<0.01);NS:不显著;SEM:标准误;AE:有氧暴露;RT:室温;ND:未检测到。

The same column followed by different lowercase letters are significantly different (*: P<0.05; **: P<0.01). NS: Not significant; SEM: Standard error of mean; AE: Aerobic exposure; RT: Room temperature; ND: Not detected.

2.4 微生物演替分析

通过α多样性分析,发现细菌的ACE和Chao1这2个指数无显著差异,而Shannon指数和Simpson指数差异显著,20组和S组在室温下进行有氧暴露时,Shannon指数和Simpson指数不断增加,5 d后达到最高,且显著高于其他处理(P<0.05);而在原温度下进行有氧暴露时,Shannon指数和Simpson指数则先减小后增加。Shannon指数和Simpson指数与时间段和温度均呈极显著相关性(P<0.01),而时间段与温度的交互作用与Shannon指数无显著相关性,但与Simpson指数呈显著相关性(P<0.05)。真菌的α多样性指数在有氧暴露阶段均降低,且各处理间的差异显著(P<0.05)。ACE指数和Chao1指数与时间段、温度以及两者的交互作用均呈极显著相关性(P<0.01),Shannon指数和Simpson指数与时间段和温度极显著相关,但与两者的交互作用无显著相关性(表4)。

表4  燕麦青贮60 d和有氧暴露过程中的微生物多样性指数
Table 4  Alpha diversity of oat silage ensiled for 60 days and during aerobic exposure

时间

Time

温度

T/

细菌Bacteria真菌Fungi
ShannonSimpsonChao1ACEShannonSimpsonChao1ACE
Day 60 20 3.37d 0.81de 283.48 281.74 4.26a 0.87a 549.32a 562.44a
20/-5 2.97e 0.73g 252.28 261.59 3.29b 0.80a 297.60bc 195.57d
AE Day 1 AE 20 2.96e 0.79e 227.27 233.76 1.28def 0.41cd 172.93def 190.97d
AE 20/-5 2.53f 0.69h 236.86 250.93 1.44cdef 0.38cde 206.11cde 212.14d
RT AE 20 3.50cd 0.84cd 280.91 285.38 1.81cd 0.48bc 274.41bcd 289.26c
RT AE 20/-5 3.00e 0.74fg 253.33 265.44 1.11ef 0.27de 175.26def 195.83d
AE Day 3 AE 20 3.02e 0.79e 226.15 241.89 1.72cde 0.62b 94.95efg 137.5def
AE 20/-5 2.98e 0.78ef 252.55 258.01 0.92f 0.22e 43.05g 45.69g
RT AE 20 3.75bc 0.87bc 242.27 250.92 2.05c 0.61b 356.63b 359.74b
RT AE 20/-5 3.95b 0.88b 248.86 266.45 1.15def 0.32cde 77.80fg 85.70efg
AE Day 5 AE 20 3.65cd 0.87bc 230.30 242.80 1.74cde 0.63b 55.79fg 55.41g
AE 20/-5 3.46d 0.84cd 243.80 259.22 1.16def 0.28de 69.30fg 72.53fg
RT AE 20 4.44a 0.91ab 285.34 282.53 1.21def 0.42cd 84.38efg 85.29efg
RT AE 20/-5 4.70a 0.94a 286.79 286.45 1.50cdef 0.40cde 148.38efg 157.55de
SEM 0.061 0.007 5.037 4.896 0.061 0.018 13.865 10.661
Period ** ** NS NS ** ** ** **
Temperature ** 0.003** NS NS ** ** 0.001** **
Period×Temperature 0.054 0.031* NS NS NS 0.05 0.005** **

同列不同小写字母表示在0.05水平上差异显著(*:P<0.05;**:P<0.01);NS:不显著;SEM:标准误;AE:有氧暴露;RT:室温。

The same column followed by different lowercase letters are significantly different (*: P<0.05; **: P<0.01). NS: Not significant; SEM: Standard error of mean; AE: Aerobic exposure; RT: Room temperature.

S组中促生乳杆菌属(Levilactobacillus) 的相对丰度较高,而20组中广布乳杆菌属(Latilactobacillus)和芽孢杆菌属(Bacillus)的相对丰度则较为显著(图1A)。观察各个时间段的细菌相对丰度图(图1B)可以发现,20组中的广布乳杆菌属(Latilactobacillus)在有氧暴露后先增加,5 d后减少,而S组中逐渐减少。芽孢杆菌属(Bacillus)在2个处理中均是先增加后减少,且室温下有氧暴露时其增加的趋势出现得更早。明串珠菌属(Leuconostoc)随着有氧暴露时间的延长而逐渐减少,赖氨酸芽孢杆菌属(Lysinibacillus)在室温下有氧暴露5 d时大幅增加,尤其是20组。Hafnia-Obesumbacterium在S组原温度有氧暴露的第3天和第5天时具有较高的丰度。类芽孢杆菌属(Paenibacillus)随着有氧暴露时间的延长而逐渐增加。细菌韦恩图(图1E)显示,有氧暴露5 d后,20组和S组共有OTUs为167个,20组在原温度和室温下有氧暴露处理的特有OTUs分别为42个和47个,而S组特有OTUs分别为51个和71个。根据真菌相对丰度图(图1C)可知,与S组相比,20组中的青霉菌属(Penicillium)和Wickerhamomyces的相对丰度较高,而毕赤酵母(Pichia)的相对丰度较低。S组中的Issatchenkia、镰刀菌属(Fusarium)和毛霉菌属(Mucor)的相对丰度高于20组。在整个有氧暴露过程中(图1D) 20组和S组的真菌组成结构不同,变化也不同,20组的Wickerhamomyces先增加后减少,而室温下有氧暴露的S组中的Issatchenkia的丰度逐渐增加。根据真菌韦恩图(图1F)发现,20组在原温度和室温下有氧暴露处理的特有OTUs分别为7个和14个,而S组原温度和室温下暴露处理的特有OTUs分别为29个和222个。

fig

图1  有氧暴露过程中细菌和真菌群落分析。A:20组与S组在有氧暴露阶段的细菌相对丰度比较;B:20组与S组在有氧暴露各个阶段的细菌相对丰度;C:20组与S组在有氧暴露阶段的真菌相对丰度比较;D:20组与S组在有氧暴露各个阶段的真菌相对丰度;E:有氧暴露5 d后各个处理细菌韦恩图;F:有氧暴露5 d后各个处理真菌韦恩图。A表示原温度有氧暴露;RA表示室温下有氧暴露。

Figure 1  Bacterial and fungal community during aerobic exposure. A: Comparison of the relative abundance of bacteria during the aerobic exposure phase between group 20 and group S; B: The relative abundance of bacteria at various days of aerobic exposure between group 20 and group S; C: Comparison of relative abundance of fungi during the aerobic exposure phase between group 20 and group S; D: The relative abundance of fungi at various days of aerobic exposure between group 20 and group S; E: Venn plots of bacteria after 5 days aerobic exposure; F: Venn plots of fungi after 5 days aerobic exposure. A: Aerobic exposure under freeze-thaw; RA: Aerobic exposure at room temperature.

细菌LEfSe分析(图2)进一步揭示了不同处理组在不同时间段的生物标志物。青贮60 d后,20组主要以明串珠菌属(Leuconostoc)为生物标志物(biomaker) (图2A)。有氧暴露1 d后(图2B),20组在室温下有氧暴露的biomakers为广布乳杆菌属(Latilactobacillus)和芽孢杆菌属(Bacillus)。相比之下,S组在原处理温度下有氧暴露的biomakers为促生乳杆菌属(Levilactobacillus)、L.brevisHafnia-Obesumbacterium;在室温下有氧暴露的biomakers为Enterobacteriacese。有氧暴露3 d后(图2C),20组在原温度下有氧暴露的biomakers为广布乳杆菌属(Latilactobacillus),室温下有氧暴露的biomakers为芽孢杆菌属(Bacillus)和赖氨酸芽孢杆菌属(Lysinibacillus);S组在原温度下有氧暴露的biomakers为促生乳杆菌属(Levilactobacillus)、L. brevisHafnia-Obesumbacterium,室温下有氧暴露的biomakers为肠球菌属(Enterococcus)和类芽胞杆菌(Paenibacillus)。有氧暴露5 d后(图2D),20组在原温度下有氧暴露的biomakers为芽孢杆菌属(Bacillus)和广布乳杆菌属(Latilactobacillus),在室温下有氧暴露的biomakers为赖氨酸芽孢杆菌属(Lysinibacillus)、普罗威登斯菌属(Providencia)和类芽孢杆菌属(Paenibacillus);S组在原温度下有氧暴露的biomakers为促生乳杆菌属(Levilactobacillus)、L. brevisHafnia-Obesumbacterium和肠球菌属(Enterococcus);在室温下有氧暴露的biomakers为节杆菌属(Arthrobacter)、布鲁氏菌属(Brucella)、纤维微菌属(Cellulosimicrobium)和Brodetella

fig

图2  有氧暴露阶段细菌LEfSe分析。A:60 d开袋;B:有氧暴露1 d;C:有氧暴露3 d;D:有氧暴露5 d。

Figure 2  Bacterial variations in samples during aerobic exposure using LEfSe analysis. A: 60 days silage; B: Aerobic exposure for 1 day; C: Aerobic exposure for 3 days; D: Aerobic exposure for 5 days.

根据真菌LEfSe分析(图3),青贮60 d后20组中的biomakers有青霉菌属(Penicillium)、VishniacozymaWickerhamomyces;S组中的biomakers为镰刀菌属(Fusarium)、毛霉菌属(Mucor)和Gibberella (图3A)。有氧暴露1 d后(图3B),20组在原温度下有氧暴露的biomakers为青霉菌属(Penicillium);在室温下有氧暴露的biomakers为AcremoniumWickerhamomycesSaccharomycetaceae_sp.和Blumeria。S组在原温度下有氧暴露的biomakers为Issatchenkia;在室温下有氧暴露的biomakers为毕赤酵母(Pichia)。有氧暴露3 d后(图3C),20组在原温度下有氧暴露的biomakers为青霉菌属(Penicillium)和Wickerhamomyces;在室温下有氧暴露的biomakers为Saccharomycetaceae_sp.和Alternaria。S组在原温度下有氧暴露的biomakers为毕赤酵母(Pichia),在室温下有氧暴露的biomakers为Issatchenkia。有氧暴露5 d后(图3D),20组在原温度下有氧暴露的biomakers为青霉菌属(Penicillium)、WickerhamomycesSaccharomycetaceae_sp.;S组在原温度下有氧暴露的biomaker为Fusarium_sp.;在室温下有氧暴露的biomaker为Issatchenkia

fig

图3  有氧暴露阶段真菌LEfSe分析。A:60 d开袋;B:有氧暴露1 d;C:有氧暴露3 d;D:有氧暴露5 d。

Figure 3  Fungal variations in samples during aerobic exposure using LEfSe analysis. A: 60 days silage; B: Aerobic exposure for 1 day; C: Aerobic exposure for 3 days; D: Aerobic exposure for 5 days.

通过matestate分析(图4),20组和S组之间的差异细菌(图4A)主要有片球菌属(Pediococcus)、广布乳杆菌属(Latilactobacillus)、乳酪杆菌属(Lacticaseibacillus)和肠球菌属(Enterococcus)等,其中前4个菌属在20组中的数量显著多于S组,而后4个菌属在S组中的数量显著多于20组(P<0.05)。20组和S组之间的差异真菌(图4B)主要有青霉菌属(Penicillium)、毕赤酵母(Pichia)、IssatchenkiaWickerhamiellaWickerhamomyces等,其中除青霉菌属(Penicillium)和Wickerhamomyces在20组中的数量显著多于S组外,其余均是S组显著多于20组(P<0.05)。

fig

图4  细菌(A)和真菌(B)群落的matestate分析

Figure 4  Matestate analysis of bacterial (A) and fungal (B) communities. *: P<0.05.

通过细菌与发酵指标的相关性分析(图5),在20组(图5A)中pH值与明串珠菌属(Leuconostoc)、促生乳杆菌属(Levilactobacillus)和广布乳杆菌属(Latilactobacillus)呈极显著负相关,与芽孢杆菌属(Bacillus)、赖氨酸芽孢杆菌属(Lysinibacillus)、类芽胞杆菌(Paenibacillus)、腐败乳杆菌属(Loigolactobacillus)和肠球菌属(Enterococcus)呈极显著正相关(P<0.01)。LA与明串珠菌属(Leuconostoc)、促生乳杆菌属(Levilactobacillus)、Hafnia-ObesumbacteriumSphingomonas呈极显著正相关,与芽孢杆菌属(Bacillus)、赖氨酸芽孢杆菌属(Lysinibacillus)、类芽胞杆菌(Paenibacillus)、腐败乳杆菌属(Loigolactobacillus)和肠球菌属(Enterococcus)呈极显著负相关(P<0.01)。AA与明串珠菌属(Leuconostoc)、促生乳杆菌属(Levilactobacillus)、广布乳杆菌属(Latilactobacillus)呈极显著正相关,与芽孢杆菌属(Bacillus)和类芽胞杆菌(Paenibacillus)呈极显著负相关(P<0.01)。在S组(图5B)中,促生乳杆菌属(Levilactobacillus)和广布乳杆菌属(Latilactobacillus)与pH和氨态氮呈极显著负相关,与LA和AA呈极显著正相关(P<0.01)。芽孢杆菌属(Bacillus)、赖氨酸芽孢杆菌属(Lysinibacillus)、类芽胞杆菌(Paenibacillus)和肠球菌属(Enterococcus)与pH值呈极显著正相关(P<0.01)。类芽胞杆菌(Paenibacillus)和肠球菌属(Enterococcus)与LA和AA呈极显著负相关,芽孢杆菌属(Bacillus)与PA和BA呈极显著正相关。

fig

图5  有氧暴露阶段20(A)S(B)的细菌群落与发酵指标的相关性分析

Figure 5  The correlation analysis of silage quality and differential bacteria communities at 20 ℃ (A) and freeze-thaw condition (B) during aerobic exposure. *: P<0.05; **: P<0.01.

通过真菌与发酵指标的相关性分析(图6)可知,在20组中(图6A),青霉菌属(Penicillium)与pH值呈极显著负相关,与LA和AA呈显著正相关(P<0.05),而毕赤酵母(Pichia)和Wickerhamomyces与pH和氨态氮呈极显著正相关(P<0.01),与LA呈极显著负相关(P<0.01)。在S组中(图6B),Candida与pH值极显著正相关,与LA极显著负相关。BlumeriaVishniacozyma与pH和氨态氮呈极显著负相关(P<0.01),与LA和AA呈极显著正相关(P<0.01)。

fig

图6  有氧暴露阶段20(A)S(B)的真菌群落与发酵指标的相关性分析

Figure 6  The correlation analysis of silage quality and differential fungal communities at 20 ℃ (A) and freeze-thaw condition (B) during aerobic exposure. *: P<0.05, **: P<0.01.

3 讨论

3.1 有氧暴露后青贮发酵品质和营养成分的变化

青贮饲料的有氧稳定性是确保其营养价值不受损失、有效抑制霉菌孢子及霉菌毒素污染,进而为动物提供优质饲料的关键。青贮60 d后,S组的pH值显著高于20组,且2个处理下的pH值均高于优质青贮通常要求的pH 4.2,这可能受到贮存温度和原料附生乳酸菌数量的共同影响。研究表明低温(<20 ℃)会抑制青贮发酵,导致pH值偏高、pH下降速率减[

15]以及有机酸含量降低,而冻融条件会进一步加剧这种抑制作用。无论是20组还是S组,在室温条件下进行有氧暴露都会加速青贮饲料的腐败过程。这可能是由于在高温高湿的有氧环境中,酵母菌、肠杆菌以及霉菌得以迅速增殖。这些微生物会消耗发酵过程中产生的乳酸及其他可利用底物,特别是在S组中,肠杆菌、酵母菌及霉菌的增长更为显著。研究发现,有氧变质过程中主要涉及的生化和微生物因素包括DM、乙酸、丁酸以及酵母菌和霉菌,其中由于霉菌具有孢子且繁殖需要时间,因此在有氧暴露过程中的微生物演替中较晚出现,而有氧变质初期主要是酵母菌起主导作[16]。酵母菌兼性厌氧且耐酸,因此是有氧暴露阶段最先出现的微生[17],这与本研究结果一致,即霉菌较酵母菌出现得晚。值得注意的是,S组在青贮60 d开袋时即已检测到肠杆菌的存在,这合理解释了为何20组的氨态氮含量显著低于S组。肠杆菌和梭菌作为关键的蛋白水解微生物,随着其数量的增加,蛋白质不断水解,导致氨态氮含量持续上升。此外,20组和S组在有氧暴露1 d后的肠杆菌数量和氨态氮含量与青贮60 d时的差异并不显著,显示出肠杆菌与氨态氮之间的密切相关性。

3.2 有氧暴露后细菌多样性的变化及其影响

青贮60 d后,20组和S组中的乳酸菌相关OTUs占主导地位,但其丰度相对较低,并非绝对优势菌,这与前人的研究结果不一[

18]。这可能是由于本研究采用的青贮原料附生的乳酸菌数量较少或竞争力不足,难以维持稳定的高乳酸菌丰度。如表4所示,有氧暴露后细菌的ACE指数和Chao1指数并未发生显著变化,表明这一过程中细菌的种类数与丰富度变化不显[1],Shannon指数和Simpson指数随着有氧暴露时间的延长而逐渐增加,反映出乳酸菌丰度的不断降低,特别是在有氧暴露第3天时,细菌群落结构出现明显不同,与原始温度下的有氧暴露相比,室温下的有氧暴露加速了微生物的更替,这与韦恩图的分析结果一致,说明室温下有氧暴露导致杂菌增多且出现得更早,尤其是经过冻融处理的样品。片球菌属(Pediococcus)、广布乳杆菌属(Latilactobacillus)、乳酪杆菌属(Lacticaseibacillus)等对发酵有益的益生菌在20组中的数量显著多于S组,而变形菌属(Proteus)、肠球菌属(Enterococcus)、FusobacteriumHafnia-Obesumbacterium在S组中的数量则显著多于20组,变形菌属(Proteus)属于肠杆菌科(Enterobacteriaceae),不发酵乳糖,会产生苯丙氨酸和赖氨酸,具有脱氨作[19],对热敏感且不耐酸,能在-20 ℃存活2-3[20],这可能是其在S组中显著增多的原因。肠球菌属(Enterococcus)作为乳酸菌能够产生胞外多糖(exopolysaccharides, EPS),具有抗氧化、抗菌、抗生物膜等活[21],但是它并非绝对的益生菌,在某些条件下与生物胺的产生和食物变质有关,且产酸效率较[22],在本研究中,肠球菌属与发酵特性呈负相关,尤其是在S组中,可能是在冻融条件下菌株产生了脱氨和脱羧作用,促进了蛋白质水解,从而增加了NH3-N含[22]。在青贮60 d后及整个有氧暴露期间,20组的芽孢杆菌属(Bacillus)数量相较于S组更为丰富。芽孢杆菌属(Bacillus)种类繁多,功能各异:部分种类能够促进发酵进程,有效抑制有害菌群的生长,并有助于降低蛋白水解程度和纤维含[23];然而,也有部分种类会引发罐头食品腐败及食物中毒问[24-25]。通过相关性分析,本研究中发现芽孢杆菌属(Bacillus)与LA和AA呈负相关,与pH值和NH3-N呈正相关,表明其抑制了发酵过程,且S组中的芽孢杆菌属(Bacillus)还与PA和BA呈极显著正相关,可能是在冻融条件下芽孢杆菌属(Bacillus)引起了丁酸堆积和蛋白水[26]。S组中的Hafnia-Obesumbacterium数量多于20组。Wang[27]认为Hafnia-Obesumbacterium在青贮中与更高的乙酸、丙酸、乙醇和氨态氮含量有关,但本研究的结果与此不完全一致,在冻融条件下Hafnia-Obesumbacterium显著多于20组,但是在冻融条件下该菌与发酵无相关性而在20组中却表现出促进发酵的作用。通过细菌LEfSe分析,有氧暴露1 d时,室温下暴露的S组的biomarker是肠杆菌科(Enterobacteriaceae),它将竞争性地消耗底物,发生脱氨和脱羧作[28],加速腐败进程。到第3天时,室温下暴露的S组的biomarker是类芽胞杆菌(Paenibacillus)和肠球菌属(Enterococcus),类芽胞杆菌(Paenibacillus)是从芽孢杆菌属后期划分出来[29],在有氧变质的玉米青贮中有报道,且饲料中出现的类芽胞杆菌(Paenibacillus)是乳制品生产中非常受关注的问[30-31]。在厌氧产孢菌中,类芽胞杆菌(Paenibacillus)在pH值大于4.5,温度高于30 ℃时占主导地位,有研究发现该菌属与玉米青贮的有氧稳定性密切相关,可通过多种手段来抑制类芽胞杆菌(Paenibacillus)的产生以提高玉米青贮的有氧稳定[30],而在本研究的20组中室温下有氧暴露5 d才成为biomarker,这表明20组的有氧稳定性更好。

3.3 有氧暴露后真菌多样性的变化及其影响

青贮开袋前后的真菌群落结构存在显著差异。与发酵指标呈正相关的BlumeriaVishniacizyma等真菌迅速减少,毕赤酵母(Pichia)、WickerhamomycesIssatchenkia等成为优势菌种。在大多数青贮饲料有氧暴露后,酵母菌是导致其腐败的主要原因,其中具有消耗乳酸作用的酵母,如IssatchenkiaSaccharomycesCandida和毕赤酵母(Pichia),通常是青贮饲料有氧腐败的主要引发[

2]。S组中的毕赤酵母(Pichia)、IssatchenkiaCandida的数量显著多于20组,尤其是S组中的毕赤酵母(Pichia)在有氧暴露第1天时其丰度就超过了75%,这一现象在室温下有氧暴露时尤为明显。Wickerhamomyces在20组中的数量显著多于S组,有研究发现该属中的某些种会破坏敏感酵母的细胞[32],这可能也是20组中酵母菌数量相对较少的原因之一。除了酵母菌,霉菌也是有氧暴露阶段需要重点关注的真菌,常见的主要有镰刀菌(Fusarium)、青霉菌(Penicillium)、曲霉菌(Aspergillus)[33]。在20组中,青霉菌属(Penicillium)的数量显著多于S组,在有氧暴露第1天时,其丰度高达73.6%,但在第3天时却骤然下降,可能是因为毕赤酵母(Pichia)在有氧暴露后迅速增加,而研究表明毕赤酵母(Pichia)会抑制青霉菌(Penicillium)的生[34],这也说明S组中青霉菌属(Penicillium)丰度较低的原因。另外,通过相关性分析发现青霉菌(Penicillium)与pH值呈极显著负相关,与LA和AA呈显著正相关,即青霉菌有利于发酵,有研究还发现该菌能够通过果胶酶、淀粉酶、蛋白酶、脂肪酶等分解食品中的复杂结构元素,释放生长所需的基质,因此广泛应用于食品发酵[35],在青贮过程中,青霉菌可能也发挥了类似的功能,且受到处理温度的影响,因为在S组中并未观察到此相关性。青贮60 d后,S组中的镰刀菌属(Fusarium)、毛霉菌属(Mucor)和Gibberella等是biomaker,在有氧暴露过程中酵母菌迅速繁殖成为主要的biomaker,但是有氧暴露5 d后,S组中的Fusarium_sp.和A. restrictus为主要的biomarker。镰刀菌属(Fusarium)是食品和饲料中最重要的产真菌毒素的真菌属之一,其几乎所有种类均具备产生真菌毒素的能力,而常见的镰刀菌真菌毒素主要有脱氧雪腐镰刀菌烯醇、玉米赤霉烯酮和伏马菌素B1[36]。曲霉菌属(Aspergillus)和毛霉菌属(Mucor)是引起青贮腐败发霉的主要真[37]。有氧暴露过程中,好氧微生物通过消耗青贮饲料中的营养物质而产生热量,酵母和霉菌通过消耗有机酸(如乳酸和乙酸)提高了青贮饲料的pH值和温度,刺激了好氧和厌氧孢子形成菌的生[30],这就解释了类芽胞杆菌(Paenibacillus)在有氧暴露后期丰度增加的原因,它与酵母菌和霉菌之间存在相互联系。此外,研究还发现,反复的冻融过程会对Bacillus的形态产生显著影响,并削弱其抑菌活[13]。因此,在S组中,芽孢杆菌属(Bacillus)的相对丰度明显低于20组。这一现象充分说明,在有氧暴露阶段,细菌与真菌之间存在着复杂的相互作用。

4 结论

有氧暴露后,酵母和霉菌迅速代谢乳酸,导致pH值不断升高,乳酸和乙酸含量降低。同时,细菌和真菌的种群结构也发生了变化,S组的群落更替较快,尤其是在室温下有氧暴露时,S组中的好氧腐败菌更多。这表明冻融加速了有氧变败过程,尤其是在室温下有氧暴露。

作者贡献声明

李海萍:研究构思和设计,论文撰写和修改;贾志锋:数据收集和处理;刘文辉:数据收集和处理;马祥:数据收集与处理;周青平:论文修改;关皓:研究构思和设计。

利益冲突

作者声明不存在任何可能会影响本文所报告工作的已知经济利益或个人关系。

参考文献

1

LI HP, GUAN H, JIA ZF, LIU WH, MA X, LIU Y, WANG H, ZHOU QP. Freeze-thaw condition limits the fermentation process and accelerates the aerobic deterioration of oat (Avena sativa) silage in the Qinghai-Xizang Plateau[J]. Frontiers in Microbiology, 2022, 13: 944945. [百度学术] 

2

PAHLOW G, MUCK R, DRIEHUIS F, OUDE ELFERINK S, SPOELSTRA SF. Microbiology of Ensiling[M]. Agronomy Publication No. 42. Madison, WI: American Society of Agronomy, 2003: 31-93. [百度学术] 

3

WOOLFORD MK. The detrimental effects of air on silage[J]. The Journal of Applied Bacteriology, 1990, 68(2): 101-116. [百度学术] 

4

BOLSEN KK, DICKERSON JT, BRENT BE, SONON RN, DALKE BS, LIN C, BOYER JE. Rate and extent of top spoilage losses in horizontal silos[J]. Journal of Dairy Science, 1993, 76(10): 2940-2962. [百度学术] 

5

KUNG L, SHEPERD AC, SMAGALA AM, ENDRES KM, BESSETT CA, RANJIT NK, GLANCEY JL. The effect of preservatives based on propionic acid on the fermentation and aerobic stability of corn silage and a total mixed ration[J]. Journal of Dairy Science, 1998, 81(5): 1322-1330. [百度学术] 

6

LINDGREN S, PAHLOW G, OLDENBURG E. Influence of microbes and their metabolites on feed and food quality[c]. In Proceedings of the 19th General Meeting of the EGF, La Rochelle, France, 2002, 503-511. [百度学术] 

7

JUN SM, TAKAGI H. Stress-tolerance of baker’s-yeast (Saccharomyces cerevisiae) cells: stress-protective molecules and genes involved in stress tolerance[J]. Biotechnology and Applied Biochemistry, 2009, 53(Pt 3): 155-164. [百度学术] 

8

JI XM, LIU MH, YANG JL, FENG FJ. Meta-analysis of the impact of freeze-thaw cycles on soil microbial diversity and C and N dynamics[J]. Soil Biology and Biochemistry, 2022, 168: 108608. [百度学术] 

9

DONG ZH, LI JF, CHEN L, WANG SR, SHAO T. Effects of freeze-thaw event on microbial community dynamics during red clover ensiling[J]. Frontiers in Microbiology, 2019, 10: 1559. [百度学术] 

10

RANJIT NK, KUNG L. The effect of Lactobacillus buchneri, Lactobacillus plantarum, or a chemical preservative on the fermentation and aerobic stability of corn silage[J]. Journal of Dairy Science, 2000, 83(3): 526-535. [百度学术] 

11

LEE MH. Official methods of analysis of AOAC International (16th edn)[J]. Trends in Food Science and Technology, 1995, 6(11):382. [百度学术] 

12

Van SOEST PJ, ROBERTSON JB, LEWIS BA. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition[J]. Journal of Dairy Science, 1991, 74(10): 3583-3597. [百度学术] 

13

BRODERICK GA, KANG JH. Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media[J]. Journal of Dairy Science, 1980, 63(1): 64-75. [百度学术] 

14

BEUCHAT LR, COPELAND F, CURIALE MS, DANISAVICH T, GANGAR V, KING BW, LAWLIS TL, LIKIN RO, OKWUSOA J, SMITH CF, TOWNSEND DE. Comparison of the SimPlateTM total plate count method with PetrifilmTM, RedigelTM, and conventional pour-plate methods for enumerating aerobic microorganisms in foods[J]. Journal of Food Protection, 1998, 61(1): 14-18. [百度学术] 

15

ALI M, CONE JW, KHAN NA, HENDRIKS WH, STRUIK PC. Effect of temperature and duration of ensiling on in vitro degradation of maize silages in rumen fluid[J]. Journal of Animal Physiology and Animal Nutrition, 2015, 99(2): 251-257. [百度学术] 

16

OHYAMA Y, HARA S, MASAKI S. Analysis of factors affecting aerobic deterioration of grass silages[C]. Forage Conservation in the 80s BGS Occasional Symposium No 11, United Kingdom, 1980: 257-261. [百度学术] 

17

ÁVILA CLS, CARVALHO BF. Silage fermentation-updates focusing on the performance of micro-organisms[J]. Journal of Applied Microbiology, 2020, 128(4): 966-984. [百度学术] 

18

DUNIERE L, XU SW, LONG J, ELEKWACHI C, WANG YX, TURKINGTON K, FORSTER R, McALLISTER TA. Bacterial and fungal core microbiomes associated with small grain silages during ensiling and aerobic spoilage[J]. BMC Microbiology, 2017, 17(1): 50. [百度学术] 

19

DAVID B. Feigin and cherry’s textbook of pediatric infectious diseases[J]. Journal of Paediatrics and Child Health, 2011, 47(5): 319-320. [百度学术] 

20

KUSHWAHA K, BABU D, JUNEJA VK. Proteus[M]//Encyclopedia of Food Microbiology. Oxford: Academic Press, 2014: 238-243. [百度学术] 

21

KAVITAKE D, DEVI PB, DELATTRE C, REDDY GB, SHETTY PH. Exopolysaccharides produced by Enterococcus genus: an overview[J]. International Journal of Biological Macromolecules, 2023, 226: 111-120. [百度学术] 

22

GIRAFFA G. Lactic acid bacteria: enterococcus in milk and dairy products[M]//Encyclopedia of Dairy Sciences (Third Edition). Oxford: Academic Press, 2022: 151-159. [百度学术] 

23

BAI J, XU DM, XIE DM, WANG MS, LI ZQ, GUO XS. Effects of antibacterial peptide-producing Bacillus subtilis and Lactobacillus buchneri on fermentation, aerobic stability, and microbial community of alfalfa silage[J]. Bioresource Technology, 2020, 315: 123881. [百度学术] 

24

刘巧云. 干香菇中热噬淀粉芽孢杆菌耐热性的研究[J]. 食品安全导刊, 2024, (13): 47-50. [百度学术] 

25

DRIEHUIS F. Silage and the safety and quality of dairy foods: a review[J]. Agricultural and Food Science, 2013, 22: 16-34. [百度学术] 

26

SILVA VP, PEREIRA OG, LEANDRO ES, da SILVA TC, RIBEIRO KG, MANTOVANI HC, SANTOS SA. Effects of lactic acid bacteria with bacteriocinogenic potential on the fermentation profile and chemical composition of alfalfa silage in tropical conditions[J]. Journal of Dairy Science, 2016, 99(3): 1895-1902. [百度学术] 

27

WANG SR, LI JF, ZHAO J, DONG ZH, DONG D, SHAO T. Silage fermentation characteristics and microbial diversity of alfalfa (Medicago sativa L.) in response to exogenous microbiota from temperate grasses[J]. World Journal of Microbiology & Biotechnology, 2021, 37(12): 204. [百度学术] 

28

ROOKE J, HATFIELD R. Biochemistry of ensiling[M]//Silage Science and Technology. vol. 42. Madison: Amer Society of Agronomy, 2003: 95-139. [百度学术] 

29

宋丽丽, 张雪, 张更新, 魏鑫丽, 王先平. 类芽孢杆菌(Paenibacillus sp.)蛋白质组学研究中三种蛋白提取方法的比较分析[J]. 微生物学报, 2024, 64(2): 623-632. [百度学术] 

SONG LL, ZHANG X, ZHANG GX, WEI XL, WANG XP. Comparison of three protein extraction methods for the proteomics research on Paenibacillus sp.[J]. Acta Microbiologica Sinica, 2024, 64(2): 623-632 (in Chinese). [百度学术] 

30

BORREANI G, DOLCI P, TABACCO E, COCOLIN L. Aerobic deterioration stimulates outgrowth of spore-forming Paenibacillus in corn silage stored under oxygen-barrier or polyethylene films[J]. Journal of Dairy Science, 2013, 96(8): 5206-5216. [百度学术] 

31

BORREANI G, FERRERO F, NUCERA D, CASALE M, PIANO S, TABACCO E. Dairy farm management practices and the risk of contamination of tank milk from Clostridium spp. and Paenibacillus spp. spores in silage, total mixed ration, dairy cow feces, and raw milk[J]. Journal of Dairy Science, 2019, 102(9): 8273-8289. [百度学术] 

32

GUO FJ, MA Y, XU HM, WANG XH, CHI ZM. A novel killer toxin produced by the marine-derived yeast Wickerhamomyces anomalus YF07b[J]. Antonie Van Leeuwenhoek, 2013, 103(4): 737-746. [百度学术] 

33

ORSI RB, CORRÊA B, POSSI CR, SCHAMMASS EA, NOGUEIRA JR, DIAS SMC, MALOZZI MAB. Mycoflora and occurrence of fumonisins in freshly harvested and stored hybrid maize[J]. Journal of Stored Products Research, 2000, 36(1): 75-87. [百度学术] 

34

BOYSEN ME, BJÖRNEHOLM S, SCHNÜRER J. Effect of the biocontrol yeast Pichia Anomala on interactions between Penicillium roqueforti, Penicillium carneum, and Penicillium paneum in moist grain under restricted air supply[J]. Postharvest Biology and Technology, 2000, 19(2): 173-179. [百度学术] 

35

MULLAN WMA. STARTER CULTURES. | Importance of Selected Genera[M]//Encyclopedia of Food Microbiology (Second Edition). Oxford: Academic Press, 2014: 515-521. [百度学术] 

36

De NIJS M. Mycotoxins[M]//Encyclopedia of Food Microbiology. Oxford: Academic Press, 1999: 1520-1526. [百度学术] 

37

ALONSO VA, PEREYRA CM, KELLER LAM, DALCERO AM, ROSA CAR, CHIACCHIERA SM, CAVAGLIERI LR. Fungi and mycotoxins in silage: an overview[J]. Journal of Applied Microbiology, 2013, 115(3): 637-643. [百度学术]