-
化肥一直以来作为农业增产的重要手段，我国化肥消费量达5.56×1010 kg，占世界化肥总消费量的34%[1]，我国大部分耕地土壤贫瘠[2]，加上长期施用化肥、除草剂和杀虫剂及无节制的掠夺式生产模式，破坏了土壤结构，造成土壤营养元素严重缺失，土壤酸化，土地生产能力显著下降[3-4]，基于目前化肥不当使用带来的危害以及国家对于减少化肥使用的政策驱动，生物有机肥以具有改良土壤肥力、促进作物营养吸收、提高农作物产量及改善品质等独特的优势成为化肥的替代品，对农业的发展具有重大意义[5-7]。
-
海藻肥是以天然海藻为原料，通过物理和生化工艺提取、精制或再混配以氮磷钾和微量元素制成的生物肥，含有海藻特有的海藻酸钠、矿物质、氨基酸和植物自然生长调节物质等[8]。很多研究表明，施用海藻肥后可促进植物根系发育和主茎的伸长，提高农作物产量，改善果实品质，增强对干旱环境的抵抗力[9-10]，还具有抗病抑菌作用和防止病虫害的发生[11]。葡萄籽肥是以天然葡萄籽制成的肥料，含有多酚化合物、氨基酸、矿物质、脂肪等有益于植物生长的物质[12-14]，多用于医学，可以温和抑制口腔病原菌微生物，是一种非常有前景的天然抗菌物质[15]。目前，对于生物有机肥的研究主要集中在不同类型功能微生物的筛选、生产工艺和技术以及对农作物产量的影响[16]，而关于生物有机肥对土壤微生物的研究却少见报道。微生物是土壤中重要的组成部分，它们在改善土壤质量、提高土地生产力和土壤健康可持续发展方面发挥着重要作用[17]，不同的施肥处理很大程度上影响土壤微生物的群落结构与数量[18-19]。因此，对施加生物有机肥后土壤微生物群落结构变化的研究具有重要意义。本研究以7 年生烟富3 苹果为试材，对施加45%有机肥、葡萄籽肥和海藻肥以后树体生物量、产量及果实品质进行了测定，同时，采用高通量测序技术，对土壤微生物群落结构进行研究，为今后生物有机肥在果园生产上的有效应用提供理论依据。
-
1 材料与方法
-
1.1 试验园概况
-
试验在文登市泽头镇峰山村果园进行，该园地处东经121°43′～122°19′，北纬36°52′～ 37°23′，年均降水量762.2 mm，降水分布不均， 夏季集中，6 ～ 9 月份降水量约占全年70%，春秋季降水偏少，常发生干旱；年平均气温11.5℃，极端最高气温36.4℃，极端最低气温-25.5℃，年日照时数2 540.7 h，无霜期194 d。试验地土壤类型为棕壤土，建园前为杨树林，土壤肥力水平较差。
-
1.2 试验材料
-
供试品种为7 年生‘烟富3’，砧木为平邑甜茶（Malus hupenensis Rehd.），2010 年春季定植建园，株行距3.0 m×5.0 m，南北行向，树形为小冠疏层形，行间自然生草，树盘起垄栽培，有滴灌系统，管理水平偏上。
-
1.3 试验设计
-
试验设3 个有机肥处理（表1），分别为海藻有机肥11 250 kg/hm2、葡萄籽有机肥4 500 kg/hm2、 45%有机肥11 250 kg/hm2（对照CK），其中，海藻有机肥有机质含量≥ 45%，N+P+K ≥ 5%；葡萄籽有机肥有机质含量≥ 90%，N+P2O5+K2O ≥ 12%； 45%有机肥有机质含量≥ 45%，有效活菌数≥ 0.2 亿/g，N+P+K ≥ 5%。选择树势均匀一致的苹果树为样本树，10 株树为一小区，3 次重复，每处理共计30 株树。2016 年10 月，在样本树四周距树干30 cm左右挖深度约30 cm的放射沟，将有机肥进行沟施，然后覆土、浇水。生长季除6 月份每株树追施500 g左右果树专用复合肥（N ≥ 16%， P2O5 ≥ 6%，K2O ≥ 32%）外，不再施其他肥料， 其它管理按常规进行。
-
1.4 样品采集
-
1.4.1 土壤样品采集
-
于果实成熟期进行土壤样品采集。采用五点取样法，分别在处理的30 株样本树东、南、西、北4 个方位距树干60 cm处用直径4 cm土壤取样器采集0 ～ 20、20 ～ 40 cm土壤，去除杂质后进行混合，装进封口袋放进干冰中带回实验室，-80℃冰箱保存用于微生物多样性分析。
-
1.4.2 果实采集
-
果实成熟后在样本树树冠外围1.3 ～ 1.5 m高度随机采果50 个，运回实验室测定单果重、纵横径、着色指数、光洁度指数、果肉硬度、可溶性固形物含量等品质指标。
-
1.5 测定项目及方法
-
1.5.1 高通量测序方法
-
采用MIO-BIO PowerSoil DNA Isolation Kit试剂盒对样本的基因组DNA进行提取，之后取3 μL进行1.2%琼脂糖凝胶电泳，检测DNA的纯度和浓度。之后使用ABI9700 智能梯度PCR仪进行PCR扩增，第一次PCR扩增体系包括：10 μL 5xBuffe， 1 μL dNTP（10 mmol/L），1U Phusion超保真DNA聚合酶，F/R内侧引物（10 μmol/L）各1 μL，20 ng模板，ddH2O补至50 μL；第二次PCR扩增体系包括8 μL 5xBuffe，1 μL dNTP（10mM），0.8 U Phusion超保真DNA聚合酶，1 μL F/R外侧引物（10 μmol/L），5 μL模板，ddH2O补至40 μL。引物对应区域：16S V4-V5 区引物（515F和926R）， 鉴定细菌多样性；ITS1 区引物（ITS1F和ITS1R）， 鉴定真菌多样性。
-
采用AXYGEN公司的AxyPrepDNA凝胶回收试剂盒回收PCR扩增产物，采用FTC-3000TMrealtime PCR仪对凝胶回收产物进行定量，使用Tru Seq® DNAPCR-Free Sample Preparation Kit建库试剂盒进行文库构建，构建好的文库用FTC-3000TM real-time PCR仪进行实时荧光定量，文库合格后， 使用Hi Seq2500 PE250 进行上机测序。
-
1.5.2 树体结构及枝类组成的调查测定
-
于落叶后用卷尺及游标卡尺测量树高、干径、 冠径、外围新梢长度，并统计叶丛枝、短枝、中枝、长枝、徒长枝数量，计算枝类比例。
-
1.5.3 产量调查
-
果实成熟后单株数果，单位枝采果、称重，计算平均单果重和折合每公顷产量。
-
1.5.4 果实品质的测定
-
单果重用1/100 电子天平称量，果实纵、横径用数显游标卡尺测量，果实去皮硬度用GY-1 型果实硬度计测量，可溶性固形物含量用数显糖量计测定，果面色泽用日本产CI-410 色差计测定。
-
果面着色指数=∑（各级果数 × 代表级值）/（总果数 × 最高级值）×100%
-
光洁度指数=∑（各级果数 × 代表级值）/（总果数 × 最高级值）×100%
-
分级标准如表2 所示。
-
1.6 数据分析
-
测序数据采用双端fastq格式保存（fq1 和fq2）， 并根据Barcode信息对原始序列进行拆分，确定各样本对应的有效序列数。用Trimmomatic软件对有效序列进行质量过滤，用FLASH软件对每个样品的reads进行拼接[20]；再使用Mothur 1.33.3 软件去除（screen），得到优化序列（clean Tags）[21]。利用Uparse在97%相似度下进行优化序列聚类，得到OTU的代表序列，再用Uchime去除嵌合体[22]，其中16S嵌合体数据库为gold database（v20110519）， ITS嵌合体数据库为UNITE（v20140703）；然后用usearch_global筛选OTV的代表序列，通过Mothur（classify.seqs）软件将OTU代表序列与数据库比对进行物种注释，置信度阈值为0.8。
-
使用Mothur软件分析Alpha多样性[23]，包括Chao1、Ace及Simpson指数；使用R软件绘制稀释曲线、聚类图等，使用Excel 2010 和SPSS 23.0 进行数据处理。
-
2 结果与分析
-
2.1 施用不同肥料对红富士树体生长的影响
-
由表3 可以看出，与45%有机肥处理（CK） 相比，经过葡萄籽肥G和海藻肥F处理后红富士的株高和地径均有所提高，差异表现不显著；冠径在东西方向表现为F>G>CK，差异显著，而南北方向则相反，其中CK与F之间差异显著；外围新梢表现为葡萄籽肥G显著高于海藻肥F和CK；枝类比例中海藻肥处理叶丛和短枝比例最高，占总枝量的59.53%，而长枝和徒长枝比例最低，为33.23%， 其次是葡萄籽肥，分别为54.42%和37.22%，而CK比例最低，分别为53.80%和37.62%。
-
注：同列数据后不同小写字母表示处理间差异显著（P<0.05）。下同。
-
2.2 施用不同肥料对红富士产量的影响
-
由表4 可以看出， 与有机肥处理（CK） 相比，经过葡萄籽肥G处理后红富士的平均单株果数、平均单果质量、单株产量和每公顷产量均显著增加，分别增加了13.55%、4.31%、18.44%和18.44%；经过海藻肥F处理后红富士的平均单果质量、单株产量和每公顷产量均显著增加，分别增加了5.18%、1.11%和1.11%；
-
2.3 施用不同肥料对红富士果实品质的影响
-
果实品质测定结果（表5）显示，红富士着色指数和光洁度指数均表现为海藻肥F> 葡萄籽肥G> 有机肥（CK），表现出显著性差异；果形指数以及可溶性固形物含量在处理间均无明显差异；色泽L、 a和b均表现为海藻肥F最大，其中色泽L葡萄籽肥G与CK差异不显著；阳面硬度均高于阴面硬度， 但在处理间差异均不显著，综合而言，生物有机肥处理对果实外观及内在品质优于CK，其中，海藻肥处理果实品质最好，果实着色面积最大（着色指数88.89%），颜色最红（a值39.33），果面光洁度最高（光洁度指数85.48%、L值49.80），果肉阳面硬度最高（9.18 kg/cm2），可溶性固形物含量较高。
-
2.4 样品测序结果评估
-
通过数据统计软件处理各阶段样品序列数， 评估数据质量。各样品测序数据评估结果如表6 所示。
-
稀释性曲线可以用来比较测序数据量不同的样本中物种的丰富度[24]，也可以用来说明样本的测序数据量是否合理。当曲线趋向平坦时，说明测序数据量合理，由图1 可以看出，在一定范围内，随着测序数的增多，若曲线表现为急剧上升，则表示群落中有大量物种被发现；当曲线趋于平缓时， 则表示此环境中的物种并不会随着测序数量的增加而明显增多[25]。样品稀释曲线分别在测序数为20 000 条处逐渐趋向平坦，而表6 显示，细菌的优化测序数量超过25 000 条，真菌的优化序列超过40 000 条，说明样品序列充分，测序数量合理。
-
注：E_1、E_2，F_1、F_2，G_1、G_2 分别代表45%有机肥（CK）、海藻有机肥、葡萄籽有机肥的0 ～ 20、20 ～ 40 cm土样，下同。
-
图1 土壤样本细菌和真菌稀释曲线
-
2.5 施用不同肥料对土壤微生物群落组成的影响
-
从图2 可以看出，不同生物有机肥处理细菌和真菌门水平相对丰度和群落组成发生一定改变，土壤中的优势细菌群落包括酸杆菌门（Acidobacteria）、变形菌门（Proteobacteria）、芽单胞菌门（Gemmatimonadetes）、 拟杆菌门（Bacteroidetes）、 放线菌门（Actinobacteria）、 绿弯菌门（Chloroflexi），非优势细菌包括衣原体门（Gracilibacteria、 Chlamydiae）、梭杆菌门（Omitrophica、Fusobacteria）等（图2 左），同时0 ～ 20 和20 ～ 40 cm土层优势群落的相对丰度会发生调整。与有机肥处理（CK）相比，海藻肥F处理后0 ～ 20 和20 ～ 40 cm土层酸杆菌门的相对丰度均有所提高，较对照分别提高了13.18%和191.42%；降低了变形菌门和拟杆菌门的相对丰度，较对照分别降低了14.39%、13.27%和64.39%、29.67%；降低了0 ～ 20 cm土层芽单胞菌门的相对丰度，提高了20 ～ 40 cm土层芽单胞菌门的相对丰度。葡萄籽肥处理后0 ～ 20 和20 ～ 40 cm土层酸杆菌门的相对丰度均有所降低，较对照分别降低了43.52%和16.67%；提高了0 ～ 20 cm土层变形菌门的相对丰度，降低了20 ～ 40 cm土层变形菌门的相对丰度；降低了0 ～ 20 cm土层芽单胞菌门和拟杆菌门的相对丰度，较对照分别降低了18.15%和10.49%，提高了20 ～ 40 cm土层芽单胞菌门和拟杆菌门的相对丰度。值得注意的是，海藻肥处理后20 ～ 40 cm的第一大优势门是酸杆菌门，其他处理的第一大优势菌门是变形菌门。子囊菌门（Ascomycota）和接合菌门（Zygomycota）是土壤中含量最丰富的两大真菌群落（图2 右），对0 ～ 20 cm土层来说，海藻肥和葡萄籽肥处理均提高了子囊菌门的相对丰度， 其相对丰度表现为G_1>F_1>E_1，降低了接合菌门的相对丰度，其相对丰度表现为E_1>F_1>G_1； 对20 ～ 40 cm土层来说，海藻肥和葡萄籽肥处理均提高了接合菌门的相对丰度，其相对丰度表现G_1>F_1>E_1，其中海藻肥处理提高了子囊菌门的相对丰度，较对照提高了56.88%，葡萄籽肥处理降低了子囊菌门的相对丰度，较对照降低了47.04%。不同生物有机肥处理均不同程度上改变了不同土层细菌和真菌的群落组成。
-
图2 门分类水平土壤样本细菌和真菌群落组成
-
2.6 施用不同肥料对土壤微生物Alpha多样性的影响
-
Alpha多样性（Alpha diversity） 包括Chao1 指数、Ace指数、Shannon指数以及Simpson指数等。 前面3 个指数越大，最后1 个指数越小，说明样品中的物种越丰富[26-27]。由表7 可以看出，对于0 ～ 20 cm土层，海藻肥和葡萄籽肥处理后土壤细菌Ace指数和Chao1 指数变化较小；海藻肥处理后降低了土壤细菌Simpson优势度指数，提高了土壤细菌Shannon指数，葡萄籽肥处理正好与此相反，3 种生物有机肥处理土壤细菌Simpson优势度指数以G最高，达到0.055 3。海藻肥和葡萄籽肥处理均提高了土壤真菌Ace、Chao1 丰富度指数和Simpson优势度指数，显著降低Shannon指数。对于20 ～ 40 cm土层，海藻肥处理后降低了土壤细菌和真菌Ace、 Chao1 丰富度指数和Shannon指数，显著提高Simpson优势度指数；葡萄籽肥处理则提高了土壤细菌和真菌Ace、Chao1 丰富度指数和Shannon指数，显著降低了Simpson优势度指数。由此可见，不同生物有机肥处理后在不同程度上改变了土壤细菌和真菌的群落多样性，海藻肥和葡萄籽肥处理后表层土壤细菌多样性增加，真菌多样性降低。
-
注：显著性检验为同一土层处理间比较。
-
2.7 土壤物种丰度聚类热图
-
Heatmap是以颜色梯度来代表数据矩阵中数值的大小并根据物种或样品丰度相似性进行聚类的1 种图形展示方式。将高丰度和低丰度的物种分块聚集，通过颜色梯度及相似程度来反映多个样品群落组成的相似性和差异性。纵向聚类表示不同物种在各样品间丰度的相似情况，2 个物种间距离越近，枝长越短，说明这2 个物种在各样品间的丰度越相似；横向聚类表示不同样品的各物种丰度的相似情况，与纵向聚类一样[16]。对土壤中相对丰度前50 的细菌和真菌属进行聚类分析（图3）得知，葡萄籽肥处理后0 ～ 20 和20 ～ 40 cm土层细菌和真菌群落均聚为1 个分支，表明其组内差异较小，E_1、F_1、G_1、G_2 和E_2 样品最终聚为一大类，最后与F_2 的样品聚为一类，表明E_1、 F_1、G_1、G_2 和E_2 的组间差异相对较小，而与F_2 在属水平的物种组成上有一定差异。在属水平上的真菌多样性丰度图（图4 左）显示，海藻肥和葡萄籽肥处理后0 ～ 20 和20 ～ 40 cm土层真菌优势属的丰度与有机肥处理（CK）差异较小， Geratobasidlum 优势属的丰度在土壤真菌群落中显著降低，其中F_2 处理 Clonostachys（螺旋聚孢霉属）、Lecanicillium（蚧霉属）、Exophiala（外瓶霉属）、Mrakia（木拉克酵母属）、Gliomastix（粘鞭霉属）、Coprinus（鬼伞属）等14 个优势属在土壤真菌群落中显著降低；在属水平上的细菌多样性丰度图（图4 右）显示，葡萄籽肥处理后G_1（0～ 20 cm）和G_2（20 ～ 40 cm）优势属的丰度差异较小， 海藻肥处理后F_1（0 ～ 20 cm）和F_2（20 ～ 40 cm）优势属的丰度差异显著，F_1（0 ～ 20 cm） 和CK的E_2（20 ～ 40 cm） 优势属的丰度差异较小。
-
图3 不同生物有机肥处理在属水平上的细菌和真菌群落结构聚类图
-
图4 属水平土壤样品细菌和真菌丰度聚类热图
-
3 讨论
-
3.1 生物有机肥对红富士树体生物量的分析
-
中国虽是世界第一大苹果生产国，但并不是苹果生产强国[28]。由于过度追求苹果大小和产量， 长期化肥单施和过量施用导致土壤板结、酸化、盐渍化，果树根系生长环境恶化，使得苹果产量和品质下降[29]。宋以玲等[30]研究发现用不同比例和种类的生物有机肥代替化学氮肥的投入在一定程度上可以提高作物产量及品质，这与本试验研究结果一致，3 种有机肥处理均在不同程度上提高了红富士树体的生长。有研究发现，生物有机肥是采用完全腐熟的加工工艺，虫卵死亡率能够达到95%以上[31]，还具有使用方便，营养元素丰富，改善土壤条件，提升作物生产质量，改善作物根系微生物群等特点[32-34]，为树体的生长提供了良好的土壤环境，从而促进了果树的株高和地径，张奇等[35] 研究发现生物有机肥在施用量为6 750 kg/hm2 时土壤养分含量和生物功能改善效果明显，作物产量在短期内可以显著得到提升，本试验研究也发现海藻有机肥在施用量为11 250 kg/hm2、葡萄籽有机肥在施用量为4 500 kg/hm2 时可以明显提高红富士果树的产量，且果实品质优良，这可能与生物有机肥处理可以改善土壤性质，提高土壤养分含量有关[36]，这与姜洁等[37]研究结果一致。
-
3.2 生物有机肥对红富士土壤微生物群落的分
-
微生物群落在土壤有机质分解和营养元素循环中有着至关重要的作用，因此，土壤微生物群落的数量、结构和活性与生态系统功能密切相关，只有健康、良性循环的土壤生态环境才能保证植物的健康生长[38]，曲成闯等[39]、袁英英等[40]研究表明，施用生物有机肥可以提高黄瓜和西红柿土壤的酶活性，调控土壤微生物群落结构，增强土壤微生物生态系统的稳定性。随着分子生态学研究技术的发展，高通量测序技术成为了研究土壤微生物多样性和群落结构特征的重要途径[41]，本研究利用高通量测序技术对施加生物有机肥后不同土层的土壤进行研究发现，施加生物有机肥之后增加了0 ～ 20 cm土层细菌的多样性和丰富度，减少了真菌的丰富度和多样性，对于20 ～ 40 cm土层，海藻肥处理后减少了土壤细菌和真菌的丰富度和多样性；葡萄籽肥处理则提高了土壤细菌和真菌的丰富度和多样性，这刘振香等[42]的研究结果一致；贺文员等[43]研究发现施加生物有机肥（QD）之后变形菌门（Proteobacteria）、子囊菌门（Ascomycota）分别是细菌和真菌中的优势物种，α-多样性分析显示QD使用前后真菌的多样性显著减少，这与本试验研究结果一致，施用海藻肥处理后不同土层子囊菌门的相对丰度上升，施用葡萄籽肥处理后0 ～ 20 cm土层子囊菌门的相对丰度上升，20 ～ 40 cm土层子囊菌门的相对丰度下降；子囊菌门大多数为腐生菌，可以分解木质素、角质素等难降解的有机质， 在养分循环中担任着重要角色[44-45]；酸杆菌门能够降解复杂的木质素和纤维素，为土壤提供养分， 还会分泌植物激素IAA从而促进作物生长[46-47]，变形菌包括了能与植物共生的固氮细菌，对土壤氮、磷、硫和有机质循环有重要作用[48]，本试验研究结果发现，海藻肥处理0 ～ 20 和20 ～ 40 cm土层酸杆菌门（Acidobacteria）的相对丰度提高， 变形菌门（Proteobacteria）的相对丰度降低；葡萄籽肥处理后降低了0 ～ 20 和20 ～ 40 cm土层酸杆菌门的相对丰度，提高了0 ～ 20 cm土层变形菌门的相对丰度，降低了20 ～ 40 cm土层变形菌门的相对丰度。土壤物种丰度聚类热图结果显示，海藻肥和葡萄籽肥处理后在属水平上细菌差异较小，真菌中被孢霉属（Mortierella）和酵母属（Guehomyces）明显增加，李发虎等[49]研究发现施入生物炭后增加了被孢霉科真菌比例，可促进土壤有机质和养分含量的提高，从而利于黄瓜根系生长并提高根系体积。被孢霉属是一类可以产生多烯不饱和脂肪酸，如花生四烯酸的真菌，被广泛应用于食品、饲料和医用领域[50]，还具有溶解土壤中磷的功能[51]，酵母属（Guehomyces）也属于有益菌， 可以有效抑制致病菌对于植物的伤害，此外，酵母属菌种还有一定的重金属吸附作用，可以有效减少重金属对植物的伤害[52-53]。综上所述，生物有机肥的使用一定程度上对果园土壤的改良产生了积极的影响。
-
4 结论
-
生物有机肥（海藻肥和葡萄籽肥）处理后提高了树体的株高、冠径和产量，且果实品质优良；海藻肥处理后树体长势、叶丛枝和短果枝比例、平均单果质量以及果实品质均优于葡萄籽肥，但其公顷产量低于葡萄籽肥。
-
海藻肥处理后0 ～ 20 和20 ～ 40 cm土层酸杆菌门（Acidobacteria）的相对丰度提高，变形菌门（Proteobacteria）和拟杆菌门（Bacteroidetes）的相对丰度降低；葡萄籽肥处理后降低了0 ～ 20 和20 ～ 40 cm土层酸杆菌门的相对丰度，提高了0 ～ 20 cm土层变形菌门的相对丰度，降低了20 ～ 40 cm土层变形菌门的相对丰度，拟杆菌门相对丰度的变化与变形菌门相反。值得注意的是，海藻肥处理后20 ～ 40 cm的第一大优势菌门是酸杆菌门，其他处理的第一大优势菌门是变形菌门。
-
子囊菌门（Ascomycota）和接合菌门（Zygomycota） 是土壤中含量最丰富的两大真菌群落，与45%有机肥处理相比较，对0 ～ 20 cm土层来说，海藻肥和葡萄籽肥处理均提高了子囊菌门的相对丰度，其相对丰度表现为G_1>F_1>E_1，降低了接合菌门的相对丰度，其相对丰度表现为E_1>F_1>G_1；对20 ～ 40 cm土层来说，海藻肥和葡萄籽肥处理均提高了接合菌门的相对丰度，其相对丰度表现为G_1>F_1>E_1，其中海藻肥处理提高了子囊菌门的相对丰度，较CK提高了56.88%，葡萄籽肥处理降低了子囊菌门的相对丰度。
-
海藻肥和葡萄籽肥处理后增加了0 ～ 20 cm土层细菌的多样性和丰富度，减少了真菌的丰富度和多样性，3 种生物有机肥处理土壤细菌Simpson优势度指数以G最高，达到0.055 3；对于20 ～ 40 cm土层，海藻肥处理后减少了土壤细菌和真菌的丰富度和多样性；葡萄籽肥处理则提高了土壤细菌和真菌的丰富度和多样性。
-
参考文献
-
[1] 朱兆良，金继运．保障我国粮食安全的肥料问题［J］．植物营养与肥料学报，2013，19（2）：259-273．
-
[2] 杨卿，郎南军，苏志豪，等．土壤退化研究综述［J］．林业调查规划，2009，34（1）：20-24．
-
[3] 薛仁风，丰明，赵阳，等．不同生物有机肥对绿豆生长与生理特性的影响［J］．东北农业科学，2019（4）：9-12，71．
-
[4] 刘戈，易玉林．微生物肥料的发展现状与前景展望［J］．安徽农业科学，2007（11）：3318，3332．
-
[5] 周玉飞，康专苗，彭竹晶．珍稀濒危铁皮石斛的研究进展 ［J］．基因组学与应用生物学，2018，37（4）：1629-1635．
-
[6] 任海虹，王景雪，聂菁．铁皮石斛原球茎高效再生体系的研究［J］．中草药，2017，48（19）：4057-4061．
-
[7] 付小猛，毛加梅，沈正松，等．中国生物有机肥的发展现状与趋势［J］．湖北农业科学，2017，56（3）：401-404．
-
[8] 袁蕊，王学江，李峰．海藻肥中海藻酸钠的提取及分子量检测［J］．浙江农业科学，2019（8）：1275-1277，1279．
-
[9] 王健健，蒲静，刘筱，等．不同浓度海藻肥对白三叶生长的影响［J］．贵州农业科学，2018，46（3）：50-52．
-
[10] 董华芳，李松，许延波，等．海藻有机肥对草莓叶绿素和果实可溶性固形物的影响［J］．现代农业科技，2018（12）：52-54．
-
[11] 涂海华，周坚，毛宇，等．天然海藻肥对“夏黑”葡萄植株生长及果实品质的影响［J］．中国土壤与肥料，2019（4）：213-217．
-
[12] 张国治，韩宝丽，王伟玲，等．葡萄籽油中脂肪酸成分分析 ［J］．粮食科技与经济，2011，36（1）：49-50．
-
[13] 陈仁．营养保健食品［M］．北京：中国轻工业出版社，2001．75-82．
-
[14] 周建设，王建新．原花青素抗氧化的临床研究［J］．中国医院药学杂志，2007，27（12）：1717-1718．
-
[15] Savickiene N，Jekabsone A，Raudone L，et al．Efficacy of proanthocyanidins from root extract in reducing viability while preserving oral commensal［J］．Materials（Basel），2018，11（9）：1499-1513．
-
[16] 付小猛，毛加梅，沈正松，等．中国生物有机肥的发展现状与趋势［J］．湖北农业科学，2017，56（3）：401-404．
-
[17] 李世贵，吕天晓，顾金刚，等．土壤微生物分子生态学研究方法［J］．中国土壤与肥料，2008（6）：1-4，14．
-
[18] 朱英，郭永婷，田兴武，等．施肥对设施番茄生长及土壤化学性质的影响［J］．农业工程技术，2018，38（34）：71-74．
-
[19] Celik I，Gunal H，Budak M，et al．Effects of long-term organic and mineral fertilizers on bulk density and penetration resist ance in semi-arid Mediterranean soil conditions［J］．Geoderma，2010，160（2）：236-243.
-
[20] Magoč T，Salzberg S L．FLASH：fast length adjustment of short reads to improve genome assemblies［J］．Bioinformatics，2011，27（21）：2957-2963．
-
[21] Bokulich N A，Subramanian S，Faith J J，et al．Qualityfiltering vastly improves diversity estimates from Illumina amplicon sequencing［J］．Nature Methods，2013，10（1）：57-59．
-
[22] Haas B J，Gevers D，Earl A M，et al．Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons［J］．Genome Research，2011，21（3）：494-504．
-
[23] Edgar R C．UPARSE：highly accurate OTU sequences from microbial amplicon reads［J］．Nature Methods，2013，10（10）：996-998．
-
[24] Wang Y，Sheng H F，He Y，et al．Comparison of the levels of bacterial diversity in freshwater，intertidal wetland，and marine sediments by using millions of illumina tags［J］．Applied and Environmental Microbiology，2012，78（23）：8264-8271．
-
[25] 王宁，刘铜，靳亚忠，等．木霉菌对土壤微生物多样性及草莓生长和发病的影响［J］．江苏农业科学，2018，46（18）：108-112．
-
[26] Pitta D W，Parmar N，Patel A K，et al．Bacterial diversity dynamics associated with different diets and different primer pairs in the rumen of Kankrej cattle［J］．Plos One，2014，9（11）：e111710．
-
[27] Sun R B，Zhang X X，Guo X S，et al．Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw［J］． Soil Biology & Biochemistry，2015，88（4）：9-18．
-
[28] 林惠娇，张海娥，蒋湘，等．中国苹果产业进出口及其植物检疫情况［J］．落叶果树，2019，51（4）：22-26．
-
[29] 刘宝芝．烟台苹果产业现状、存在问题及发展建议［J］．农家参谋，2019（7）：59．
-
[30] 宋以玲，于建，陈士更，等．化肥减量配施生物有机肥对油菜生长及土壤微生物和酶活性影响［J］．水土保持学报，2018，32（1）：352-360．
-
[31] 詹永海．黄皮果树配合使用生物有机肥有利其优质丰产稳产［J］．粮食科技与经济，2019（8）：75-76，96．
-
[32] 张丽英，张正斌，徐萍，等．黄淮小麦农艺性状进化及对产量性状调控机理的分析［J］．中国农业科学，2014，47（5）：1013-1028．
-
[33] 和凤美，朱芮，朱永平，等．甜玉米自交系性状相关分析和主成分分析［J］．作物杂志，2014（3）：32-35．
-
[34] 薛云云，白冬梅，田跃霞，等．24 份山西花生资源农艺性状的相关性和主成分分析［J］．山西农业科学，2017（10）：1587-1590．
-
[35] 张奇，张振华，刘丽珠，等．增施有机肥对黄泛冲积区贫瘠土壤养分和玉米产量的影响［J］．江苏农业科学，2019（17）：271-276．
-
[36] 刘斌祥，王兴龙，周芳，等．减氮配施不同种类有机肥对玉米物质分配、转运与产量的影响［J］．生态学杂志，2020（1）：130-138．
-
[37] 姜洁，龚一富，郭蓉，等．海藻生物肥对草莓产量和品质的影响［J］．核农学报，2019，33（5）：1032-1037．
-
[38] Zeller V，Bardgett R D，Tappeiner U．Site and management effects on soil microbial properties of subalpine meadows：a study of land abandonment along anorth-south gradient in the European Alps［J］．Soil Biology and Biochemistry，2001，33（4）：639-649．
-
[39] 曲成闯，陈效民，韩召强，等．施用生物有机肥对黄瓜不同生育期土壤肥力特征及酶活性的影响［J］．水土保持学报，2017，31（6）：279-284．
-
[40] 袁英英，李敏清，胡伟，等．生物有机肥对番茄青枯病的防效及对土壤微生物的影响［J］．农业环境科学学报，2011，30（7）：1344-1350．
-
[41] 张风革，霍云倩，孙艺，等．连续施用生物有机肥对草地生物量及土壤微生物区系的影响［J］．南京农业大学学报，2018，41（2）：382-388．
-
[42] 刘振香，刘鹏，贾绪存，等．不同水肥处理对夏玉米田土壤微生物特性的影响［J］．应用生态学报，2015，26（1）：113-121．
-
[43] 贺文员，宋清晖，杨尚霖，等．生物有机肥对水稻土壤酶活性及微生物群落结构的影响［J］．中国农学通报，2019，35（27）：106-113．
-
[44] Beimforde C，Feldberg K，Nylinder S，et al．Estimating the phanerozoic history of the ascomycota lineages：Combining fossil and molecular data［J］．Molecular Phylogenetics and Evolution，2014，78：386-398．
-
[45] 赵帆，赵密珍，王钰，等．基于高通量测序研究草莓根际微生物群落结构和多样性［J］．土壤，2019，51（1）：51-60．
-
[46] Kielak A M，Cipriano M A，Kuramae E E．Acidobacteria strains from subdivision 1 act as plant growth-promoting bacteria［J］． Archives of Microbiology，2016，198（10）：987-993．
-
[47] Lopez-Mondejar R，Voriskova J，Vetrovsky T，et al．The bacterial community inhabiting temperate deciduous forests is vertically stratified and undergoes seasonal dynamics［J］．Soil Biology and Biochemistry，2015，87：43-50．
-
[48] 王鹏，陈波，张华．基于高通量测序的鄱阳湖典型湿地土壤细菌群落特征分析［J］．生态学报，2017，37（5）：1650-1658．
-
[49] 李发虎，李明，刘金泉，等．生物炭对温室黄瓜根际土壤真菌丰度和根系生长的影响［J］．农业机械学报，2017，48（4）：265-270，341．
-
[50] 顾晓颖，王文思，陈家，等．高产花生四烯酸的高山被孢霉菌株发酵条件优化［J］．西南农业学报，2018，31（11）：2300-2306．
-
[51] 段倩倩，杨晓红，黄先智．石漠化地区生态桑林根际土壤微生物资源多样性研究［J］．西南大学学报（自然科学版），2017，39（7）：25-29．
-
[52] 王祥红，贾仁洁，李静，等．嗜杀酵母的生物学功能及其应用［J］．生物技术通报，2013（4）：33-38．
-
[53] 李明春，姜恒，侯文强，等．酵母菌对重金属离子吸附的研究［J］．菌物系统，1998（4）：80-86．
-
摘要
以 7 年生烟富 3 苹果为试材，探讨了 45% 有机肥和生物有机肥（葡萄籽肥和海藻肥）对果树生物量及土壤微生物环境的影响，为果树优质高效生产和减量施用化肥提供科学依据。试验共设置 3 个处理：45 % 有机肥 11250 kg/hm2 （E）、海藻有机肥 11250 kg/hm2 （F）、葡萄籽有机肥 4500 kg/hm2 （G），于果实成熟期采集土壤样品及果实，测定其生物量及果实品质等相关指标，应用高通量测序技术研究不同处理后不同土层微生物多样性的变化。与 45% 有机肥处理对照相比较，生物有机肥处理后提高了树体的株高、冠径、产量和品质，海藻肥处理后树体叶丛枝和短果枝比例、平均单果质量以及果实品质均优于葡萄籽肥，外围新梢表现为葡萄籽肥最高。对于 0 ～ 20 cm 土层，海藻肥处理后显著降低了土壤细菌 Simpson 指数，提高了土壤细菌 Shannon 指数，葡萄籽肥处理结果与此相反；海藻肥和葡萄籽肥处理均提高了土壤 0 ～ 20 cm 土层真菌 Ace、Chao1 和 Simpson 指数，降低了 Shannon 指数。对于 20 ～ 40 cm 土层，海藻肥处理后降低了土壤细菌和真菌 Ace、Chao1 和 Shannon 指数，显著提高了 Simpson 指数，葡萄籽肥处理结果与此相反。海藻肥处理后提高了 0 ～ 20 和 20 ～ 40 cm 土层酸杆菌门（Acidobacteria）的相对丰度，降低了变形菌门（Proteobacteria）和拟杆菌门（Bacteroidetes）的相对丰度；葡萄籽肥处理后降低了 0 ～ 20 和 20 ～ 40 cm 土层酸杆菌门的相对丰度，提高了 0 ～ 20 cm 土层变形菌门的相对丰度，降低了 20 ～ 40 cm 土层变形菌门的相对丰度，降低了 0 ～ 20 cm 土层拟杆菌门的相对丰度，提高了 20 ～ 40 cm 土层拟杆菌门的相对丰度。子囊菌门（Ascomycota）和接合菌门（Zygomycota）是土壤中含量最丰富的两大真菌，对 0 ～ 20 cm 土层来说，海藻肥和葡萄籽肥处理均提高了子囊菌门的相对丰度，降低了接合菌门的相对丰度， 对 20 ～ 40 cm 土层来说，海藻肥和葡萄籽肥处理均提高了接合菌门的相对丰度，海藻肥处理提高了子囊菌门的相对丰度，较对照提高了 56.88%，葡萄籽肥处理降低了子囊菌门的相对丰度。土壤物种丰度聚类热图结果显示不同土层之间存在差异，葡萄籽肥处理与对照处理之间差异较大。葡萄籽肥和海藻肥处理均有助于红富士树体树高、 干径、东西冠径以及外围新梢的增大，且枝类结构更为合理，不同生物有机肥处理后改变了土壤细菌和真菌的群落多样性及结构，土壤细菌多样性增加，真菌多样性降低，改善了土壤微生物环境，为生物有机肥的推广奠定了基础。
Abstract
The effects of biological organic fertilizers（grape seed fertilizer and seaweed fertilizer）on fruit tree biomass and soil microbial environment were studied by using 7-year-old Yanfu 3 apple as test material and 45% organic fertilizer as control （CK），which provides the scientific basis for high quality and efficiency production of fruit trees and fertilizer reduction. The experiment consisted of three treatments：45% organic fertilizer with an application rate of 11250 kg/hm2 （E），seaweed organic fertilizer with an application rate of 11250 kg/hm2 （F）and grape seed organic fertilizer with an application rate of 4500 kg/hm2 （G）. Soil samples and fruit were collected at maturity. The biomass and fruit quality were measured，and the changes of soil microbial diversity after different soil treatments were studied by high throughput sequencing technique. Compared with CK，the plant height，ground diameter，yield and quality of the tree were improved after the treatment of biological organic fertilizer. The proportion of leafy and short fruit branches，average single fruit quality and fruit quality of seaweed fertilizer treatment were better than those of grape seed fertilizer treatment. However，the outer tips of branch slightly showed the highest with the grape seed fertilizer treatment. For 0 ～ 20 cm soil layer，the bacterial Simpson index was significantly decreased and the bacterial Shannon index was increased after seaweed fertilizer treatment，but the result of grape seed fertilizer treatment was the opposite. The treatments of the seaweed fertilizer and the grape seed increased the Ace，Chao1 and Simpson index of the soil 0 ～ 20 cm，and the Shannon index was reduced. For 20 ～ 40 cm soil layer， seaweed fertilizer treatment decreased the Ace，Chao1 and Shannon index of soil bacteria and fungi，and significantly increased the Simpson index，but the result of grape seed fertilizer treatment was the opposite. After seaweed fertilizer treatment，the relative abundance of Acidobacteria in 0 ～ 20 and 20 ～ 40 cm soil layer was increased，which decreased the relative abundance of Proteobacteria and Bacteroidetes. With grape seed fertilizer treatment，the relative abundance of Acidobacteria in 0 ～ 20 and 20 ～ 40 cm soil layer decreased，The relative abundance of Proteobacteria in 0 ～ 20 cm soil layer was improved，the relative abundance of Bacteroidetes in 20 ～ 40 cm soil layer was reduced，the relative abundance of Proteobacteria in the 0 ～ 20 cm soil layer was reduced，and the relative abundance of Bacteroidetes in 20 ～ 40 cm soil layer was improved. Ascomyces and Zygomycota were the two most abundant fungi in soil. For 0 ～ 20 cm soil layer， both seaweed manure and grape seed fertilizer treatment increased the relative abundance of Ascomyces and decreased the relative abundance of Zygomycota. For 20 ～ 40 cm soil layer，both seaweed fertilizer and grape seed fertilizer treatment increased the relative abundance of Ascomyces，seaweed fertilizer treatment increased the relative abundance of Ascomyces， which was 56.88% higher than that of the control. And grape seed fertilizer treatment decreased the relative abundance of Ascomyces. The results of cluster heat map of soil species abundance showed that there were differences between different soil layers and between grape seed fertilizer treatment and control treatment. Grape seed fertilizer and seaweed fertilizer treatment were helpful to the increase of height，dry diameter，east-west crown diameter and peripheral fresh of red Fuji trees，and the branch structure was more reasonable. The community diversity and structure of soil bacteria and fungi were changed after different bioorganic fertilizer treatments，the diversity of soil bacteria increased，the diversity of fungi decreased，and the microbial environment of soil was improved，which laid a foundation for the promotion of high quality and high yield and biological organic fertilizer.
Keywords
grape seed fertilizer ； seaweed fertilizer ； apple ； high throughput sequencing ； bacteria ； fungi