油菜是我国的主要油料作物之一，种植面积居世界前列[1]。目前，油菜秸秆的处置方式多为露天焚烧和自然降解，会造成资源浪费和环境污染，因此探究优质高效的秸秆资源综合利用方法成为热点话题[2]。但油菜秸秆木质化程度高、生物质结构复杂，造成其利用难、利用率低。高效且低成本的降解技术是实现资源化利用油菜秸秆的亟须解决的难题。目前，油菜秸秆的降解方法主要有物理处理法、化学处理法、生物处理法，其中生物降解技术操作简单、生产费用较低且效用高，不会污染环境，在降解纤维素类物质方面具有巨大的潜力[3]。秸秆的主要成分包括木质素、纤维素和半纤维素，木质素和半纤维素以共价键的形式形成一种天然的屏障包裹纤维素，防止纤维素被酶解[4]。木质素结构复杂，对纤维素和半纤维素的酶解产生可发酵糖及高附加值产品的开发形成了一定阻碍，因此能否有效去除木质素也会影响纤维素和半纤维素的有效利用[5]。生物降解木质素主要借助于自然界中真菌、放线菌、细菌等，与细菌相比，真菌降解木质素的能力较强，且具有产生多种氧化还原酶的能力[6-7]，在木质素解聚方面起重要作用，其中关于白腐真菌的研究最多。但白腐真菌存在效率低、碳水化合物损失大、覆盖面积大、发酵周期长等缺点，因此有需要筛选新菌株、探索新方法降解木质素[8]。本研究从腐烂木头及土壤中筛选、分离、鉴定具有木质素降解能力的真菌，进一步探究了该菌株以油菜秸秆为底物的固态发酵产酶优化以及降解效果，为后续油菜秸秆生物降解提供高效微生物资源，为秸秆的资源化利用提供一定参考。1材料与方法1.1试验材料1.1.1样品来源土壤和腐烂木头采自四川轻化工大学汇南校区树林里含腐烂木头的土壤地表下2~5 cm处，使用无菌袋包装，冰袋保存带回实验室，4 °C冰箱保存。1.1.2培养基富集培养基：葡萄糖5 g、蛋白胨10 g、牛肉膏5 g、NaCl 10 g、蒸馏水1 000 mL，pH值自然，121 ℃高压灭菌20 min。初筛培养基：木质素2 g、硫酸铵2 g、磷酸二氢钾1 g、磷酸氢二钠0.2 g、琼脂20 g，pH值7.0，蒸馏水1 000 mL，121 ℃高压灭菌20 min，倒平板备用[9]。复筛培养基：愈创木酚-PDA：马铃薯200 g、葡萄糖20 g、琼脂20 g、蒸馏水定容至1 000 mL，pH值6.0~6.5，121 ℃高压灭菌20 min，冷却至不烫手，使用0.22 μm滤膜过滤，除菌，加入0.1%愈创木酚[9]。苯胺蓝-PDA：马铃薯200 g、酵母膏20 g、蛋白胨5 g、葡萄糖20 g、NaCl 5 g、蒸馏水定容至1 000 mL，pH值5.5左右，琼脂20 g；灭菌，冷却至不烫手，用0.22 μm滤膜过滤除菌加入0.1 g苯胺蓝，倒平板备用[10]。木质素降解培养基[10]：碱性木质素2 g、KH2PO4 1 g、MgSO4·7H2O 1 g、蒸馏水1 000 mL，pH值7.0。初始固态发酵产酶培养基[10]：预处理过的油菜秸秆3.0 g、麸皮2 g装入150 mL锥形瓶中，以料液比1∶2加入营养液摇匀。营养液为KH2PO4 1.0 g/L、MgSO4·7H2O 1.0 g/L、(NH4)2SO4 10 g/L、蛋白胨1.0 g/L。所有培养基均需1×105 Pa灭菌30 min后冷却。1.1.3试剂与仪器1.1.3.1试验试剂木质素、碱性木质素、铁氰化钾、乳酸钠购自上海麦克林生化科技有限公司；葡萄糖、氯化铁购自成都市科隆化学品有限公司；2,2-二氮-双（3-乙基苯并噻唑-6-磺酸）购自合肥巴斯夫生物科技有限公司；苯胺蓝、愈创木酚购自上海泰坦科技股份有限公司；藜芦醇购自上海阿拉丁生化科技股份有限公司。1.1.3.2试验仪器HH-4恒温水浴锅购自上海力辰邦西仪器科技有限公司，LRH-300生化培养箱购自常州诺基仪器有限公司，UV-1800紫外/可见分光光度计购自上海美谱达仪器有限公司，SW-CJ-2D型超净工作台购自上海苏净实业有限公司，PHS-3C酸度计购自成都世纪方舟科技有限公司，LS-50HJ立式压力蒸汽灭菌器购自江阴滨江医疗设备有限公司，F2000全自动纤维分析仪购自济南海能仪器股份有限公司）。1.2试验方法1.2.1油菜秸秆的预处理参考王辉等[11]预处理秸秆方法，称取20.00 g粉碎至过20目筛的油菜秸秆粉，置于250 mL锥形瓶中，加入200 mL 1.1% NaOH（液料比为10 mL/g），密封，置于45 ℃环境中浸泡36 h，置于滤袋中洗至中性，烘干至恒重，称量。1.2.2木质素降解菌株的分离筛选取10 g土壤样品于已灭菌的50 mL富集培养基中，放入摇床28 ℃、180 r/min富集7 d。取出在超净工作台内取上清液1 mL梯度稀释至10-1~10-8，吸取10-5、10-6、10-7、10-8各200 μL涂布于初筛培养基上，每个梯度设置3个重复。28 ℃恒温培养箱中培养3 d，选择生长旺盛的菌株转接至初筛分离培养基上进行划线分离纯化，直至得到纯菌，并对筛选得到菌株进行编号。挑取初筛培养基上的单菌落，使用牙签点种至愈创木酚-PDA、苯胺蓝-PDA培养基上的中心，28~30 ℃恒温培养5 d，分别挑取周边有红色显色圈、脱色圈明显的菌株，反复划线纯化，保藏。1.2.3木质素降解菌的鉴定将筛选出的菌株M20根据菌落培养特征和菌体结构特征，参考《真菌鉴定手册》[12]进行初步鉴定，将菌种送至上海生工进行内转录间隔区（ITS）测序，测序结果使用BLAST程序在NCBI GenBank数据库进行同源性分析，利用MEGA 5.0软件计算序列的相似度，使用Neighbor-joining法构建菌株系统发育树，确定物种的系统发育地位[13]。1.2.4木质素降解菌株降解率测定采用普鲁士法[14]测木质素浓度，将菌株涂布于初筛培养基平板上，均匀铺满平板后，使用打孔器做成直径8 mm菌塞，每个木质素降解培养基中无菌操作接种5个菌塞，28 ℃、150 r/min降解5 d，取出，摇匀，过滤掉菌的孢子，滤液使用0.22 μm滤膜抽滤，所得滤液测木质素含量。木质素含量标准曲线操作步骤参考文献[13]，得到回归方程y=0.028 2x+0.030 1（R2=0.992 2）。1.2.5木质素降解菌株酶活测定木质素胞外漆酶（Lac）、木质素过氧化物酶（Lip）和锰过氧化物酶（MnP）活性的测定具体操作步骤参考文献[15-17]。1.2.6固态发酵产酶单因素优化碳、氮源优化：在初始固态发酵产酶培养基的基础上，分别选取秸秆与麸皮比为4∶1、3∶1、2∶1、3∶2、1∶1，选取蛋白胨、硫酸铵、硝酸铵作为氮源，初始添加量为3%。选出最优氮源后，分别选取1%、2%、3%、4%添加量，在发酵条件为接种量10%、28 ℃、5 d进行试验。采用最优碳氮源，即秸秆、麸皮比3∶1，硫酸铵2%，分别选取液料比1、2、3、4 L/g，发酵温度24、28、32、36、40 ℃，发酵时间为2、3、4、5、6、7 d，接种量分别取5%、10%、20%、25%、30%进行单因素优化。1.2.7固态发酵产酶优化正交试验设计根据单因素试验结果，选择对固态发酵产酶影响较大的因素设计正交试验，进一步获得最佳工艺条件，以秸秆与麸皮比（A）、时间（B）和氮源添加量（C）作为3个因素，选取3个水平进行试验。采用L9（34）正交表进行正交试验设计，确定固态发酵产酶的最佳工艺。正交试验因素及水平设计见表1。10.13557/j.cnki.issn1002-2813.2023.09.018.T001表1正交试验因素水平设计水平秸秆与麸皮比A时间B/d氮源添加量C/%空列-14∶132103∶143212∶15431.2.8验证试验以最优发酵工艺检验正交法所获得的结果的可靠性，通过比较初筛发酵条件和优化后的酶活，评价优化前后发酵条件对固态发酵产酶的影响。1.2.9油菜秸秆降解率的测定方法参照GB/T 20805—2006、GB/T 6434—2006，采用滤袋法分别测定饲料中的酸性洗涤木质素含量和粗纤维含量。将发酵油菜秸秆粉碎过40目筛，称取约0.5 g左右发酵油菜秸秆于滤袋封口，将滤袋放入烧杯中用石油醚浸没10 min脱脂，在通风橱晾干。放入F2000全自动纤维分析仪中，设置不同程序进行消煮和洗涤，取出滤袋，沥干水分，放入丙酮中浸没5 min，通风橱晾干，放入烘箱（100±2）℃烘干至恒重，马弗炉灰化至恒重，计算酸性洗涤木质素含量、粗纤维含量。2结果与分析2.1木质素降解菌的筛选从初筛培养基上共挑出56株菌，共有20株菌出现褪色圈或/和显色圈。其中14株菌株同时在愈创木酚-PDA培养基、苯胺蓝-PDA培养基上显色/褪色圈出现时间较快且直径圈较大，表明其产漆酶和过氧化物酶能力较强。但仅通过褪色圈/显色圈只能简单判断产酶的特性，无法具体了解酶的活性，因此还需后续试验加以验证。将14株菌株分别接入木质素降解培养基中通过发酵后测木质素降解率进行下一步复筛，结果显示，降解率大于30%的有10株菌，再分别接入以油菜秸秆为唯一碳源的初始固态产酶培养基中，以筛选出高效降解油菜秸秆的菌株。10株木质素降解菌的木质素酶活见表2。10.13557/j.cnki.issn1002-2813.2023.09.018.T002表210株木质素降解菌的木质素酶活菌株漆酶锰过氧化物酶木质素过氧化物酶M83.27±0.110431.18±5.02M203.04±0.3395.89±12.691 129.03±19.99M152.87±0.9794.10±6.57834.41±10.34M241.48±0.53142.02±9.46912.90±64.32M8-22.29±0.5796.94±8.91670.97±14.90M4-23.02±0.42158.33±2.74997.85±55.96M41.97±0.3567.72±2.11778.49±53.53M151.91±0.0888.76±1.23960.22±16.08M411.29±0.2637.80±5.13883.87±18.06M562.77±0.9927.71±1.23855.91±18.86U/g由表2可知，菌株M8的漆酶最高，但其余两种酶活很低；M20、M4-2的3种酶活均较高，M20的均略高于M4-2，且M20的苯胺蓝褪色效果以及木质素降解效果比M4-2好，因此最终确定M20为目的菌株进行后续研究。2.2木质素降解菌的鉴定2.2.1菌株M20的形态特征及显微观察结果（见图1）由图1可知，菌株M20在苯胺蓝-PDA、愈创木酚-PDA平板上均呈黄绿色，丝绒质地，菌丝为致密短绒状，向四周扩展，边缘不规整，在平板上呈中央凸起生长。显微镜观察菌体形态，菌株具有分生孢子梗，孢梗茎壁薄且长，壁平滑，顶端有球形顶囊，由顶囊的小梗产生孢子，分生孢子头辐射状。结合邓维琴等[18]对米曲霉的鉴定研究结果，初步判断筛选出的M20菌株为米曲霉。10.13557/j.cnki.issn1002-2813.2023.09.018.F001图1菌株M20的形态特征及显微观察结果2.2.2菌株M20的ITS序列分析及系统发育树双向测序获得菌株M20的ITS基因序列，经拼接后在NCBI数据库中进行序列同源性比对，选取BLAST中相似度较高的真菌序列，利用MEGA5.0计算真菌序列相似度，采用其中的Neighbor Joining法构建系统发育树，结果见图2。由图2可知，M20与Aspergillus oryzae BS A2在同一个分支上，序列相似性为100%。因此，将该菌株初步命名为A. oryzae M20。10.13557/j.cnki.issn1002-2813.2023.09.018.F002图2菌株M20系统发育树2.3固态发酵产酶优化单因素试验结果（见图3~图9）由图3可知，当秸秆与麸皮比例为3∶1时，木质素过氧化物酶和锰过氧化物酶均达到峰值，漆酶的酶活略低于秸秆与麸皮比为2∶1和1∶1；考虑秸秆的利用率最大化，选择秸秆与麸皮比例4∶1~2∶1进行后续试验。10.13557/j.cnki.issn1002-2813.2023.09.018.F003图3秸秆与麸皮比对固态发酵产酶的影响由图4、图5可知，以硫酸铵为氮源时，3种酶活性均最高。随着硫酸铵含量的增加，3种酶的活性均上升再下降，其中漆酶和木质素过氧化物酶的在硫酸铵含量为2%时达到峰值，锰过氧化物酶在硫酸铵含量为3%时达到峰值；当硫酸铵含量为4%时，3种酶活性均下降。因此，选择2%~4%硫酸铵含量进行后续试验。10.13557/j.cnki.issn1002-2813.2023.09.018.F004图4不同氮源对固态发酵产酶的影响10.13557/j.cnki.issn1002-2813.2023.09.018.F005图5氮源添加量对固态发酵产酶的影响由图6可知，当接种量在5%~20%范围内时，3种酶活波动幅度不大。因此，接种量对产酶结果影响不大，分析原因是接种量过多，微生物生长旺盛，使环境竞争加剧，微生物群提前衰亡期，且在进入主要代谢阶段之前底物已经大量被消耗，使最终产物减少。10.13557/j.cnki.issn1002-2813.2023.09.018.F006图6接种量对固态发酵产酶的影响由图7可知，随着发酵时间延长，3种酶活均逐渐升高，到6 d时，3种酶活性均下降。10.13557/j.cnki.issn1002-2813.2023.09.018.F007图7发酵时间对固态发酵产酶的影响因此，本试验选择发酵时间为3~5 d进行后续试验。由图8可知，24~32 ℃范围内，3种酶活性随着温度的升高而逐渐增大，在米曲霉最适生长温度32~36 ℃范围内[19]，3种酶活性均逐渐下降。因此，选取32 ℃作为后续试验。10.13557/j.cnki.issn1002-2813.2023.09.018.F008图8温度对固态发酵产酶的影响由图9可知，液料比为1~2时，漆酶与木质素过氧化物酶活性均随着液料比增加而增加，而锰过氧化物酶活性呈下降趋势；液料比大于2时，随时液料比增加，3种酶活性均下降，可能是固态发酵适宜霉菌生长，水分含量过大会导致培养基透气性不足，影响菌株生长，产酶降低[20]。10.13557/j.cnki.issn1002-2813.2023.09.018.F009图9液料比对固态发酵产酶的影响2.4固态发酵产酶正交试验优化结果为考察各因素之间的关系及影响，进一步优化各因素的最佳量，应用SPSS 24软件，设计秸秆与麸皮比、发酵时间、氮源添加量3因素3水平的正交试验，正交试验结果见表3，方差分析见表4~表6。10.13557/j.cnki.issn1002-2813.2023.09.018.T003表3正交试验结果及极差分析项目ABC空列木质素过氧化物酶/(U/g)锰过氧化物酶/(U/g)漆酶/(U/g)111111 141.9436.931.85212221 818.2883.592.03313331 042.0123.080.92421231 797.85125.133.15522311 541.94102.560.74623121 046.2487.692.85731321 339.7877.441.97832131 434.41143.593.52933211 661.29152.823.97k11 334.0771 426.5231 207.5301 448.390k21 462.0101 598.2101 759.1401 401.433k31 478.4931 249.8471 307.9101 424.757R木144.416348.363551.61046.957k'147.86779.83389.40397.437k'2105.127109.913120.51382.907k'3124.61787.86367.69397.267R锰76.75030 .08052.82014.530k''11.6002.3232.7402.187k''22.2472.0973.0502.283k''33.1532.5801.2102.530R漆1.5530.4831.8400.34310.13557/j.cnki.issn1002-2813.2023.09.018.T004表4锰过氧化物酶方差分析项目平方和自由度均方F值P值修正模型15 234.939a62 539.15612.1710.078截距77 065.461177 065.461369.3870.003A9 550.02424 775.01222.8870.042B1 455.7442727.8723.4890.223C4 229.17122 114.58610.1360.090误差417.2612208.631总计92 717.6629修正后总计15 652.2008注：R2=0.973，调整后R2=0.893。10.13557/j.cnki.issn1002-2813.2023.09.018.T005表5漆酶方差分析项目平方和自由度均方F值P值修正模型9.82761.63817.4170.055截距49.000149.000521.0920.002A3.65321.82719.4240.049B0.35120.1751.8660.349C5.82322.91130.9600.031误差0.18820.094总计59.0159修正后总计10.0158注：1.R2=0.981，调整后R2=0.925。2.P0.05表示影响显著，P0.01表示影响极显著；下表同。10.13557/j.cnki.issn1002-2813.2023.09.018.T006表6木质素过氧化物酶方差分析来源平方和自由度均方F值P值修正模型737 500.264a6122 916.7174.3270.013截距18 272 037.03118 272 037.0311 049.0490.000A37 494.63218 747.3211.3360.081B182 047.62291 023.8155.0420.180C517 958.022258 979.01156.6040.006误差3 307.4421 653.72总计19 012 844.739修正后总计740 807.7048注：R2=0.996，调整后R2=0.982。由表3~表6可知，氮源添加量对3种酶活力影响均显著（P0.05），发酵时间仅对过氧化物酶活力影响显著（P0.05），对其余两种酶活力影响不显著（P0.05），而秸秆与麸皮比例对锰过氧化物酶和漆酶活力影响显著（P0.05），对木质素过氧化物酶活性影响不显著（P0.05），因此最佳方案为A3B2C2，即产木质素过氧化物酶最佳的发酵条件为秸秆与麸皮比为2∶1、时间4 d、(NH4)2SO4添加量为3%；产锰过氧化物酶最佳的发酵条件为秸秆与麸皮比为2∶1、时间4 d、(NH4)2SO4添加量为3%；产漆酶最佳的发酵条件为秸秆与麸皮比为2∶1、时间5 d、(NH4)2SO4添加量为3%。根据极差分析结果，产漆酶的最佳方案为A3B3C2，即产漆酶的最佳发酵条件中仅发酵时间与另两种酶发酵时间不同，但发酵时间仅对过氧化物酶活力影响显著，因此以过氧化物酶的最佳发酵时间为准，最终最佳方案为A3B2C2。2.5最优发酵条件的验证以正交试验得到的最优发酵培养条件组合为秸秆与麸皮比为2∶1、(NH4)2SO4添加量为3%、发酵4 d进行试验验证。优化后木质素过氧化酶从1 121.03 U/g上升到1 796.74 U/g，提高了60.28%；锰过氧化酶从101.13 U/g上升至146.78 U/g，提高了45.14%；漆酶从3.135 U/g上升至3.979 U/g，提高了27.00%。2.6油菜秸秆中木质素、纤维素的降解率在最优发酵培养基条件下，未接菌的秸秆固体培养基粗纤维含量47.76%、木质素含量为13.34%；接入菌株M20经过发酵后的秸秆固体培养基粗纤维33.09%、木质素含量为6.71%。菌株M20对油菜秸秆中木质素降解率达49.70%，粗纤维降解率达30.72%，均高于Ćilerdžić等[21]利用Ganoderma applanatum对秸秆的木质素、纤维素的降解率。将未发酵的油菜秸秆固体培养基发酵基质与发酵过后的秸秆培养基基质分别烘干、粉碎、过筛，用导电双面胶带固定在样品台上，在真空环境下进行镀金处理，物料表面形成一层导电膜后，用扫描电子显微镜进行观察秸秆表面形态变化[22]，结果见图10。10.13557/j.cnki.issn1002-2813.2023.09.018.F010图10M20固态发酵4 d前后的油菜秸秆扫描电镜结果由图10可知，未发酵的秸秆表面光滑，结构紧密且完整，试验组秸秆经菌株发酵处理后可见，菌体孢子覆盖、夹杂在秸秆细胞中，秸秆致密完整的结构被破坏。3讨论李佳少[23]将米曲霉CGMCC5992进行固态发酵后，发现锰过氧化物酶达到2.08 U/g，木质素过氧化物酶达到1.79 U/g，而通过单因素和响应面试验优化了米曲霉CGMCC5992固态发酵培养的工艺条件，得到木质素过氧化物酶最大酶活值可达到6.62 U/g，秸秆木质素降解率达40.25%。本研究筛选的米曲霉M20产的锰过氧化物酶、木质素过氧化物酶活均远高于上述研究结果，同时还能够产漆酶。张芳芳等[4]利用白腐真菌接种到玉米秸秆固体培养基中，发酵过程中漆酶活性最高值为21.36 U/g，锰过氧化物酶活性峰值0.32 U/g，木质素过氧化物酶峰值为15 U/g。与张芳芳等[4]的研究结果相比，本研究中的菌株M20的漆酶活性较低，但木质素过氧化物酶和锰过氧化物酶分别是其119倍和456倍。由M20的固态发酵产酶结果发现，米曲霉M20在油菜秸秆为底物的固态发酵中，木质素过氧化物酶在木质素降解中起主要作用，与Pinto等[24]报道的木质素过氧化物酶在固态发酵占优势的结果一致，与Fackler等[25]用白腐菌降解木质素时，锰过氧化物酶起主要作用结果不同，可能是不同菌种代谢差异导致[26]。此外，与已报道其他种属的木质素降解菌相比，M20在短时间、低能耗的发酵条件下对水溶碱性木质素和油菜秸秆（天然木质素）均具有较强的降解能力。李灵灵等[10]以白腐菌BYL-7发酵5 d对碱性木质素降解率为36.5%，培养8 d对秸秆木质素降解率为32.8%，而M20在液体培养基中发酵5 d对碱性木质素降解率约为40%，固态发酵5 d后对秸秆木质素降解率达49.70%，分别提高了4.5%和16.9%。Zhu等[27]以变色栓菌为木质素降解试验菌株，发现秸秆木质素降解率为34.8%，本试验中M20的木质素降解率提高了14.9%，且明显缩短了发酵周期。范寰等[28]利用复合木质素菌系发酵30 d时，秸秆的木质素降解率达到（57.0±0.2）％，虽M20秸秆木质素降解率比其低（7.3±0.2）%，但发酵周期缩短了7.5倍。4结论本研究分离筛选的菌株米曲霉M20具有较强的产木质素酶活和木质素降解能力，对其利用油菜秸秆固态发酵培养基进行发酵优化，在初始发酵培养基的基础上改变秸秆与麸皮比为2∶1、(NH4)2SO4添加量为3%、发酵4 d可明显提高其产木质素酶活能力，且油菜秸秆中木质素降解率达49.70%，粗纤维降解率达30.72%。该菌株在油菜秸秆生物降解方面具有较大的开发利用前景，一定程度上为秸秆资源化利用提供了新的微生物资源。后续还可利用此菌株与其他菌株结合处理油菜秸秆，多菌协同产酶，使酶系统更加完整，提高油菜秸秆的降解率。