李文慧 王繼濤 安明遠(yuǎn) 徐廣亞 尹翠 游宏建 曹云娥，*

（1寧夏大學(xué)葡萄酒與園藝學(xué)院，寧夏 銀川 750021；2寧夏回族自治區(qū)園藝技術(shù)推廣站，寧夏 銀川 750001；3寧夏共享人力集團(tuán)有限公司，寧夏 銀川 750027）

番茄（SolanumlycopersicumL.）在我國種植廣泛且經(jīng)濟(jì)利益較高，近年來，其生產(chǎn)方式基本由傳統(tǒng)的露地栽培轉(zhuǎn)變?yōu)樵O(shè)施栽培［1］。然而，設(shè)施農(nóng)業(yè)高集約化和高復(fù)種指數(shù)導(dǎo)致土壤微生物群落失衡，土傳病害加重及土壤性質(zhì)破壞等一系列土壤退化問題［2-3］。采用基質(zhì)栽培是緩解此類問題的有效措施之一［4］。前人研究以菇渣、牛糞、蛭石為原料，進(jìn)行了番茄基質(zhì)配方的篩選，發(fā)現(xiàn)當(dāng)牛糞∶菇渣∶蛭石體積比為6∶3∶1時，番茄的產(chǎn)量和品質(zhì)顯著提高［5］。

然而，目前無土栽培相關(guān)的國標(biāo)、行標(biāo)等都較少，市面銷售的栽培基質(zhì)物料配比、養(yǎng)分構(gòu)成也不盡相同，難以開展標(biāo)準(zhǔn)化無土栽培體系研究［6］。蚯蚓（Pheretima）作為生態(tài)系統(tǒng)中常見的分解者，可以將有機(jī)固體廢棄物分解轉(zhuǎn)化為腐殖質(zhì)，故蚯蚓堆肥作為一種富含有機(jī)質(zhì)、養(yǎng)分和腐殖質(zhì)的新型基質(zhì)而備受關(guān)注［7-8］。相關(guān)研究表明，施用蚯蚓堆肥可使黑籽（NigellasativaL.）籽粒氮、磷、鉀含量提高13.50%、10.80%和14.10%，黑籽油中的亞油酸、油酸和亞麻酸含量提高53.53%、25.36%和0.60%［9］。為使蚯蚓堆肥長期穩(wěn)定地向栽培體系輸送養(yǎng)分，保證番茄的穩(wěn)產(chǎn)高質(zhì)，本研究特引入“蚯蚓原位堆肥”（in-situ vermicomposting）——一種改良的蚯蚓堆肥技術(shù)。該技術(shù)將蚯蚓及共生菌群引入有機(jī)廢棄物（畜禽糞便、餐廚垃圾、農(nóng)田果園殘留物等）中，其分解、轉(zhuǎn)化有機(jī)廢棄物的全過程不僅直接在待處理地塊自發(fā)完成，亦結(jié)合當(dāng)?shù)馗髦贫葘?shí)時補(bǔ)充有機(jī)廢棄物物料，以保證養(yǎng)分長期供應(yīng)充足。此外，生物炭作為一種物理結(jié)構(gòu)穩(wěn)定、化學(xué)性質(zhì)呈現(xiàn)惰性的富碳材料，被廣泛應(yīng)用于農(nóng)業(yè)研究［10］。生物炭可影響土壤微生物群落結(jié)構(gòu)及多樣性，調(diào)控環(huán)境因子和微生物之間的互作關(guān)系，促進(jìn)并改良微生物生態(tài)系統(tǒng)，從而對受微生物驅(qū)動的氮、磷、鉀營養(yǎng)元素循環(huán)轉(zhuǎn)化產(chǎn)生影響，最終改變土壤質(zhì)量和功能［11］。Huang等［12］發(fā)現(xiàn)玉米芯生物炭既促進(jìn)了細(xì)菌的數(shù)量和多樣性，又能進(jìn)一步加速脫水污泥蚯蚓堆肥過程中有機(jī)物的礦化和腐殖化。同時，生物炭可有效減少NH3的揮發(fā)，進(jìn)而將NH4+轉(zhuǎn)化為NO3-，達(dá)到氮肥的緩釋效果，為作物和根系生長提供充足的養(yǎng)分，提高氮素利用效率［13-14］。

目前針對蚯蚓（原位）堆肥和生物炭作為栽培基質(zhì)的應(yīng)用均較成熟，但將蚯蚓原位堆肥與生物炭配合施用，探究其對養(yǎng)分、菌群及農(nóng)產(chǎn)品產(chǎn)量與品質(zhì)構(gòu)成的影響仍鮮有報道。因此，本研究通過向蚯蚓原位堆肥中添加不同比例的生物炭，探究比較不同處理的理化性質(zhì)、菌群組成及番茄代謝發(fā)育特征，旨在為此后相關(guān)研究提供基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)于2021年2月24日—7月2日在寧夏賀蘭縣園藝產(chǎn)業(yè)園區(qū)的日光溫室（106.33°E，38.58°N）進(jìn)行。該基地屬于國家級農(nóng)業(yè)示范園區(qū)，年降雨量180～200 mm，80%的降雨量集中在6—8月，年平均氣溫8.8 ℃，無霜期160～170 d，屬典型的中溫帶干旱氣候區(qū)，具有大陸性氣候特征。

1.2 試驗(yàn)材料

試驗(yàn)番茄品種為粉得力（武漢楚為生物科技股份有限公司）；試驗(yàn)所需蚯蚓和牛糞購自寧夏萬輝生物環(huán)?？萍加邢薰?；生物炭購自上海海諾碳工業(yè)有限公司。生物炭和蚯蚓原位堆肥基本特性見表1。

表1 基礎(chǔ)材料理化性質(zhì)Table 1 Physicochemical properties of basic materials

1.3 試驗(yàn)方法

采用單因素隨機(jī)區(qū)組設(shè)計。共設(shè)6個處理，分別為蚯蚓原位堆肥（CK）、蚯蚓原位堆肥+1%生物炭（C1）、蚯蚓原位堆肥+2%生物炭（C2）、蚯蚓原位堆肥+3%生物炭（C3）、蚯蚓原位堆肥+5%生物炭（C4）、蚯蚓原位堆肥+10%生物炭（C5），每個處理重復(fù)5次，隨機(jī)區(qū)組排列。

試驗(yàn)將基質(zhì)槽一分為二，種植壟和養(yǎng)殖壟尺寸均為4.1 m（長）×0.35 m（寬）×0.2 m（深），對照壟每壟鋪設(shè)蚯蚓原位堆肥226.73 kg，鋪設(shè)牛糞0.38 m3。根據(jù)等碳原則，種植壟添加生物炭替代蚯蚓原位堆肥中部分總碳，養(yǎng)殖壟投放一定密度蚯蚓處理牛糞產(chǎn)生蚯蚓堆肥。待種植壟一茬作物收獲，養(yǎng)殖壟中牛糞全部處理為蚯蚓堆肥后，養(yǎng)殖壟與種植壟模式調(diào)換，原養(yǎng)殖壟原位處理生成的蚯蚓堆肥中定植作物，原種植壟中重新鋪設(shè)牛糞進(jìn)行蚯蚓原位處理，為下一茬作物定植做準(zhǔn)備，最終實(shí)現(xiàn)作物原位輪作。番茄單行栽植于種植壟，株距15 cm，行距140 cm，所有處理采用水肥一體化滴管裝置進(jìn)行水肥管理，番茄整個生育期灌水量約1 400 t·hm-2，肥料選用寧夏萬輝生物環(huán)?？萍加邢薰旧a(chǎn)的茄果類專用滴灌配方肥（24.0% N+3.0% P2O5+21.0% K2O），全生育期共投入氮肥378 kg·hm-2、磷肥47.25 kg·hm-2、鉀肥330.75 kg·hm-2。

1.4 測定項目與方法

1.4.1 基質(zhì)理化性質(zhì)測定 于番茄盛果期（2021年5月18日），采用5點(diǎn)取樣法從栽培基質(zhì)中部取基質(zhì)500 g，每個處理5次重復(fù)，風(fēng)干后過1 mm篩。根據(jù)馬全會等［15］的方法測定基質(zhì)容重、孔隙度、含水量。使用pH儀和DDS-307A EC儀（上海精密科學(xué)儀器有限公司）測定基質(zhì)pH值和電導(dǎo)率（electric conductivity，EC）（土∶水=1∶10）。有機(jī)質(zhì)含量測定采用重鉻酸鉀氧化法；全氮（total N，TN）含量測定采用硫酸消煮-凱氏定氮-硼酸指示劑-0.01 mol·L-1硫酸滴定法；全磷（total P，TP）含量測定采用HCIO4-H2SO4消化-鉬銻抗比色法［16］。分別采用0.5 mol·L-1K2SO4、0.5 mol·L-1NaHCO3萃取-鉬銻法和NH4OAc萃取-火焰光度法測定速效氮（available N，AN）、速效磷（available P，AP）和速效鉀（available K，AK）含量［16］。

土壤脲酶活性采用苯酚鈉法測定；蔗糖酶活性采用3，5-二硝基水楊酸法測定；過氧化氫酶活性采用高錳酸鉀法測定；堿性磷酸酶活性采用磷酸苯二鈉法測定［16］。

1.4.2 微生物豐度指標(biāo)測定 于番茄盛果期（2021年5月18日）采集新鮮基質(zhì)樣品進(jìn)行微生物群落多樣性分析。根據(jù)QIAamp DNA MiniKit試劑盒（Qiagenn，德國）說明書進(jìn)行微生物群落總DNA抽提，使用1%的瓊脂糖凝膠電泳檢測DNA的提取質(zhì)量，使用NanoDrop 2000（美國Thermo Fisher Scientific公司）測定其DNA純度和濃度。采用引物338F（5'-ACTCCTACGGGAG GCAGCAG-3'）和806R（5'-GGACTACHVGGGTWTCT AAT-3'）對細(xì)菌16S rRNA基因V3-V4的可變區(qū)進(jìn)行PCR擴(kuò)增及高通量測序；引物ITS2R（5'-GCTGCGTTC TTCATCGATGC-3'）和ITS1F（5'-CTTGGTCATTTAGA GGAAGTAA-3'）對真菌進(jìn)行PCR擴(kuò)增及高通量測序。細(xì)菌PCR反應(yīng)體系：5×FastPfu Buffer緩沖液4 μL，2.5 mmol·L-1dNTPs 2 μL，5 μmol·L-1上、下游引物各0.8 μL，F(xiàn)astPfu Polymerase DNA聚合酶0.4 μL，模板DNA 10 ng，ddH2O補(bǔ)足至20 μL。真菌PCR反應(yīng)體系：10×Buffer緩沖液2 μL，2.5 mmol·L-1dNTPs 2 μL，5 μmol·L-1上、下游引物各0.8 μL，rTAP Polymerase DNA聚合酶0.2 μL，模板DNA 10 ng，ddH2O補(bǔ)足至20 μL。程序?yàn)椋?5 ℃預(yù)變性3 min；95 ℃變性30 s，55 ℃退火30 s，72 ℃延伸30 s，細(xì)菌27個循環(huán)，真菌35個循環(huán)；72 ℃終延伸10 min。測序數(shù)據(jù)優(yōu)化處理、操作分類單元（operational taxonomic unit，OTU）聚類分析及注釋、多樣性及群落結(jié)構(gòu)分析等均在上海美吉生物科技有限公司的I-Sanger平臺進(jìn)行，用QIIME 1.9.1軟件進(jìn)行稀疏分析，以區(qū)分OTU豐富度及α多樣性指數(shù)［覆蓋度（Coverage）、Shannon、Simpson、Chao1、Ace指數(shù)］［17］；使用CANOCO 5.0軟件對基質(zhì)微生物門、屬水平群落結(jié)構(gòu)、相對豐度等與基質(zhì)理化性質(zhì)進(jìn)行冗余分析（redundancy analysis，RDA）。

1.4.3 番茄生長量、品質(zhì)和產(chǎn)量測定 于盛果期，每個處理隨機(jī)選取植株6株，取下大小均勻果實(shí)10個，測定果實(shí)質(zhì)量?？扇苄怨绦挝锖坎捎肨D-45數(shù)字折光儀（中國金科利達(dá)公司）進(jìn)行測定。有機(jī)酸和維生素C（vitamin C，Vc）含量分別通過NaOH滴定法和2，6 -二氯苯酚直接滴定法進(jìn)行測定。此外，分別通過二硝基水楊酸和蒽酮比色法測定硝酸鹽和可溶性糖含量［18］。全生育期記錄每個處理產(chǎn)量，按照小區(qū)面積產(chǎn)量折合成番茄產(chǎn)量。

測定果實(shí)品質(zhì)的同時，在每個處理中隨機(jī)選取植株6株。用天平測定地上部與地下部鮮重后，每個植株裝進(jìn)信封并標(biāo)記，將樣品置于烘箱105 ℃下殺青30 min，80 ℃烘干至恒重。用電子天平稱量地上部與地下部干重。

1.4.4 高效液相色譜-質(zhì)譜聯(lián)用（chighperfomanceliquid chromatography-tandem mass spectrometry, LC-MS）檢測 在1 000 mg番茄鮮樣中加入1.3 mL提取液（甲醇∶水=4∶1，V∶V），利用Wonbio-96c型冷凍組織研磨儀（上海萬柏生物科技有限公司）破碎，低溫渦旋混勻處理30 min。4 ℃、13 000 r·min-1離心15 min，取上清后重復(fù)上述萃取操作3次，將上清液在氮?dú)庀禄旌蠞饪s。最終，樣品用乙腈∶水=1∶100 μL（V∶V）復(fù)溶后用于液相色譜-質(zhì)譜分析。上機(jī)完成之后，在美吉生物云平臺（https：//cloud.majorbio.com）進(jìn)行差異代謝物分析，篩選出差異代謝物，并通過KEGG（kyoto encyclopedia of genes and genomes）數(shù)據(jù)庫進(jìn)行代謝通路注釋，獲得差異代謝物參與的通路［19］。

1.5 數(shù)據(jù)統(tǒng)計與分析

采用Excel 2021對數(shù)據(jù)進(jìn)行整理與統(tǒng)計，采用SPSS 26.0軟件對數(shù)據(jù)進(jìn)行單因素方差分析（ANOVA）及相關(guān)性分析，采用Origin 2021軟件作圖。

2 結(jié)果與分析

2.1 不同比例生物炭對蚯蚓原位堆肥理化性質(zhì)的影響

由表2可知，添加生物炭增大了蚯蚓原位堆肥的總孔隙度，減小了容重，但各處理間無顯著差異（P＞0.05）。

表2 生物炭對蚯蚓原位堆肥物理性質(zhì)的影響Table 2 Effect of biochar on the physical properties of in-situ vermicomposting

由表3可知，與CK相比，添加生物炭能顯著提高蚯蚓原位堆肥的pH和EC值（P＜0.05），并且隨著生物炭添加比例的增大，pH和EC值呈持續(xù)上升趨勢。有機(jī)質(zhì)含量也表現(xiàn)出相同趨勢，從121.50 g·kg-1增加到310.50 g·kg-1。添加生物炭提高了蚯蚓原位堆肥的養(yǎng)分含量，并且隨著生物炭添加比例的增大，蚯蚓原位堆肥中全氮、速效氮、速效磷、全磷和速效鉀含量均呈現(xiàn)先增加后降低趨勢。其中，全磷、速效磷和速效鉀含量以C3處理最高，分別較CK顯著提高了46.27%、10.10%和68.43%；全氮和速效氮含量以C4處理最高，分別較CK顯著提高了15.96%和92.20%。

表3 生物炭對蚯蚓原位堆肥化學(xué)性質(zhì)的影響Table 3 Effect of biochar on the chemical properties of in-situ vermicomposting

2.2 不同比例生物炭對蚯蚓原位堆肥酶活性的影響

由圖1可知，添加生物炭能提高蚯蚓原位堆肥中的酶活性，并且隨著生物炭添加比例的增大，過氧化氫酶活性呈持續(xù)增大趨勢；蔗糖酶、脲酶、堿性磷酸酶活性均呈先增大后減小趨勢。其中，蔗糖酶和堿性磷酸酶活性以C3處理最高，分別較CK顯著提高了52.09%和18.66%；脲酶活性以C4處理最高，較CK顯著提高了86.49%。

圖1 生物炭對蚯蚓原位堆肥酶活性的影響Fig.1 Effect of biochar on enzyme activity of in-situ vermicomposting

2.3 不同比例生物炭對蚯蚓原位堆肥微生物多樣性的影響

由圖2可知，細(xì)菌群落中各處理共產(chǎn)生5 567個OTU，其中共有OTU數(shù)量為2 422，占總數(shù)的43.51%。CK～C5處理特有OTU數(shù)分別為41、46、44、65、51、68。相較CK，C1～C5處理特有OTU數(shù)量均有不同程度增加，以C5處理特有OTU數(shù)量最多。這可能是因?yàn)樯锾慷嗫椎慕Y(jié)構(gòu)有利于微生物的生長發(fā)育。由此表明，生物炭能增加蚯蚓原位堆肥中的特有細(xì)菌數(shù)量，對細(xì)菌的產(chǎn)生可能具有一定的促進(jìn)作用。

圖2 生物炭對蚯蚓原位堆肥細(xì)菌（A）和真菌（B）OTU水平的影響Fig.2 Effect of biochar on the OTU levels bacteria （A） and fungi （B） of in-situ vermicomposting

真菌群落中各處理共產(chǎn)生1 671個OTU，其中共有OTU數(shù)量為400，占總數(shù)的23.94%。CK～C5處理特有OTU數(shù)分別為70、91、67、83、125、89，除C2處理外，其他處理的特有OTU數(shù)量均高于CK，變化趨勢與細(xì)菌基本一致，以C4處理特有OTU最多。C5處理特有OTU小于C4處理的原因可能是過量生物炭對真菌的產(chǎn)生具有抑制作用。

為了進(jìn)一步分析生物炭對微生物群落多樣性的影響，采用五種指標(biāo)（覆蓋度、Ace指數(shù)、Chao1指數(shù)、Simpson指數(shù)以及Shannon指數(shù)）計算蚯蚓原位堆肥中的細(xì)菌和真菌豐度。由表4可知，細(xì)菌和真菌覆蓋度（Coverage）分別達(dá)96%和99%以上，說明蚯蚓原位堆肥中所測出的細(xì)菌、真菌序列數(shù)足以反映其真實(shí)情況。不同處理間細(xì)菌Chao1指數(shù)和Ace指數(shù)無顯著差異；隨著生物炭添加比例的增大，各處理細(xì)菌Shannon指數(shù)呈先降低后升高趨勢，Simpson指數(shù)呈先升高后降低趨勢，其中，C3處理的Shannon指數(shù)和Simpson指數(shù)相較CK差異顯著，細(xì)菌多樣性最低。真菌Chao1指數(shù)、Ace指數(shù)、Shannon指數(shù)和Simpson指數(shù)在不同處理間均無顯著差異。

表4 生物炭對蚯蚓原位堆肥細(xì)菌和真菌群落Alpha多樣性指數(shù)的影響Table 4 Effect of biochar on Alpha diversity index of biochar on Alpha diversity index of in-situ vermicomposting with bacterial and fungal communities

2.4 不同比例生物炭對蚯蚓原位堆肥微生物群落組成的影響

由圖3可知，在細(xì)菌門水平上，不同處理得到12個相對豐度＞1%的細(xì)菌菌門，其中優(yōu)勢菌門為放線菌門（Actinobacteria）、變形菌門（Proteobacteria）、綠彎菌門（Chloroflexi）和厚壁菌門（Firmicutes），共占到全部菌門相對豐度的87.44%～89.51%。CK～C5處理的綠彎菌門（Chloroflexi）相對豐度分別為12.69%、15.87%、13.93%、17.85%、16.39%和26.87%，變形菌門（Proteobacteria）和厚壁菌門（Firmicutes）相對豐度較CK均有所降低，以C3處理最低，分別為17.56%和11.76%。（圖3-A）。在屬水平上，共有11個優(yōu)勢細(xì)菌屬（相對豐度＞2%），其中相對豐度排名前三的細(xì)菌屬分別為norank_f__norank_o__Actinomarinales、unclassified_f__Intrasporangiaceae和沉積巖桿菌屬（Ilumatobacte）。相較于CK，生物炭處理中的norank_f__norank_o__Actinomarinale相對豐度下降了0.03～0.58個百分點(diǎn)，unclassified_f__Intrasporangiaceae相對豐度增加了0.29～0.59個百分點(diǎn)，以C4處理相對豐度最大（圖3-B）。

圖3 生物炭對蚯蚓原位堆肥細(xì)菌和真菌群落組成的影響Fig.3 Effect of biochar on the relative abundance of microorganisms of in-situ vermicomposting at the phylum and genus level

在真菌菌落中，不同處理檢測出7個相對豐度＞1%的真菌菌門，其中優(yōu)勢菌門為子囊菌門（Ascomycota）、被孢霉菌門（Mortierellomycota）和羅茲菌門（Rozellomycota）。各處理中子囊菌門（Ascomycota）相對豐度為48.05%～56.43%。相較CK，生物炭的添加使被孢霉菌門（Mortierellomycota）相對豐度減少。除C4處理外，其他處理的羅茲菌門（Rozellomycota）相對豐度均高于CK。此外，隨著生物炭添加比例的增加，油壺菌門（Olpidiomycota）作為潛在的優(yōu)勢菌門，其相對豐度在C4和C5處理中均高于CK（圖3-C）。在屬水平上，共有18個優(yōu)勢真菌屬（相對豐度＞2%），其中相對豐度排名前三的真菌屬為未分類菌屬（unclassified_k__Fung）、被孢霉屬（Mortierella）和羅茲菌屬（Rozellomycota）。相較CK，生物炭處理降低了被孢霉屬（Mortierella）相對豐度，降幅為1.4～4.09個百分點(diǎn)。除C4處理外，其他處理的未分類真菌屬（unclassified_k__Fung）和羅茲菌屬（Rozellomycota）相對豐度較均高于CK（圖3-D）。

2.5 相關(guān)指標(biāo)的冗余分析

將細(xì)菌和真菌門、屬水平相對豐度前10的種群與基質(zhì)理化性質(zhì)進(jìn)行RDA分析。結(jié)果表明（圖4），基質(zhì)理化性質(zhì)對細(xì)菌群落豐度解釋度為40.23%（門）和32.58%（屬），真菌群落豐度解釋度為60.54%（門）和40.39%（屬）。

在門水平上，影響細(xì)菌群落相對豐度的主要環(huán)境因子是全磷（TP）、速效鉀（AK）和速效磷（AP）。其中，放線菌門（Actinobacteria）、綠彎菌門（Chloroflexi）相對豐度與TP、AP、AK含量呈正相關(guān)，變形菌門（Proteobacteria）相對豐度與TP、AP、AK含量呈負(fù)相關(guān)（圖4-A）。影響屬水平細(xì)菌群落相對豐度的主要環(huán)境因子是全氮（TN）、速效鉀（AK）和全磷（TP）。沉積巖桿菌屬（Ilumatobacte）、unclassified_f__Intrasporangiace相對豐度與AK含量呈正相關(guān)，norank_f__norank_o__Vicinamibacterales相對豐度與AK含量呈負(fù)相關(guān)（圖4-B）。

對門水平的真菌群落變化產(chǎn)生關(guān)鍵作用的環(huán)境因子為全氮（TN）、速效氮（AN）和速效鉀（AK）。其中，子囊菌門（Ascomycota）相對豐度與TN、AN和AK含量呈正相關(guān)，被孢霉菌門（Mortierellomycota）、羅茲菌門（Rozellomycota）相對豐度與AK、AP含量呈負(fù)相關(guān)（圖4-C）。影響屬水平真菌群落相對豐度的主要環(huán)境因子是全氮（TN）、速效氮（AN）和速效鉀（AK）。枝孢菌屬（Cladosporiu）相對豐度與TN含量呈正相關(guān)，瓶毛殼屬（Lophotrichus）相對豐度與TN含量呈負(fù)相關(guān)，unclassified_p__Rozellomycota相對豐度與AK含量呈負(fù)相關(guān)（圖4-D）。

2.6 蚯蚓原位堆肥與不同比例生物炭對番茄品質(zhì)和產(chǎn)量的影響

由表5可知，生物炭處理整體顯著增加了番茄植株地下部生物量（P＜0.05）。與CK相比，C1～C5處理植株地下部生物量整體顯著增加16.84%～112.19%。隨著生物炭添加比例的增大，番茄植株地下部生物量呈先增大后減小趨勢。其中，C2處理番茄植株地下部生物量達(dá)到最大值。同時，生物炭處理整體顯著提高了番茄的產(chǎn)量（P＜0.05）。與CK相比，C1～C5處理番茄產(chǎn)量增加11.89%～36.36%，以C3處理產(chǎn)量最高（9 411.72 kg·667m-2）。

表5 蚯蚓原位堆肥與不同比例生物炭對番茄品質(zhì)及產(chǎn)量的影響Table 5 Effect of in-situ vermicomposting with different ratio of biochar on quality and yield of tomatoes

添加生物炭均可促進(jìn)果實(shí)品質(zhì)相關(guān)參數(shù)的增加。與CK相比，C1～C5處理Vc、可溶性固形物和可溶性糖含量分別增加6.74%～39.39%、0.10～0.40個百分點(diǎn)和0.72～3.92個百分點(diǎn)。隨著生物炭添加比例的增大，番茄Vc含量呈持續(xù)上升趨勢；可溶性固形物和可溶性糖含量呈先增大后減小趨勢，二者含量均以C3處理最高，分別為5.87%和14.92%。

2.7 蚯蚓原位堆肥與不同比例生物炭對番茄果實(shí)代謝的影響

為了更好地了解蚯蚓原位堆肥與生物炭對番茄果實(shí)代謝物的影響，在CK和C3處理之間進(jìn)行了代謝物分析，通過偏最小二乘法判別分析（partial least squares discriminant analysis，PLS-DA）圖來評估植物組織樣品之間的差異，結(jié)果顯示（圖5），PC1的貢獻(xiàn)率為30.70%，PC2的貢獻(xiàn)率為29.50%。從圖中易看出C3與CK間分離趨勢明顯，且因變量擬合指數(shù)（R2Y）＞0.9，模型預(yù)測指數(shù)（Q2）＞0.5，表明處理之間存在顯著差異，模型具備較好的穩(wěn)定性，獲得的數(shù)據(jù)可繼續(xù)進(jìn)行后續(xù)分析。此外，潛在結(jié)構(gòu)判別分析（orthogonal partial least squares discriminant analysis，OPLS-DA）建模的監(jiān)督正交投影也揭示了兩樣本的顯著差異。

圖5 番茄果實(shí)的LC-MS評價（n=3）Fig.5 Evaluation of LC-MS of tomato fruits （n=3）

由圖6可知，通過變量重要性投影（variable importance in projection，VIP）≥1、P值＜0.05、差異倍數(shù)（fold change，F(xiàn)C）＞1或＜1的篩選標(biāo)準(zhǔn)，在C3 vs CK中顯著鑒定了16種差異代謝物，其中有14種差異代謝物上調(diào)（fold change＞1），2種差異代謝物下調(diào)（fold change≤1）。由表6可知，相較于CK，C3處理中上調(diào)物質(zhì)包括1-己醇阿拉伯糖基葡糖苷、L-4-羥基谷氨酸半醛、4-羥基雙苯喹啉、5-甲氧基水楊酸等，下調(diào)物質(zhì)包括反式祖馬林和L-谷氨酰胺。

圖6 番茄果實(shí)差異代謝物火山圖（A）和熱圖（B）代謝途徑分析Fig.6 Volcano map （A） and heat map （B） of tomato fruit differential metabolites

表6 番茄果實(shí)差異代謝物Table 6 Differential metabolites of tomato fruit

為進(jìn)一步了解代謝途徑中差異代謝物的變化，將篩選出的差異代謝物映射到KEGG通路上進(jìn)行注釋和通路富集分析。由圖7可知，生物炭的施用對氨酰tRNA生物合成、ABC轉(zhuǎn)運(yùn)體蛋白（ATP-binding cassette）、乙醛酸和二羧酸代謝、嘌呤代謝和嘧啶代謝等代謝途徑有極顯著影響（P＜0.01）。由表7可知，參與乙醛酸和二羧酸代謝、氨酰tRNA生物合成途徑的差異代謝物為L-絲氨酸和L-谷氨酰胺；參與ABC轉(zhuǎn)運(yùn)蛋白途徑的差異代謝物為L-絲氨酸、L-谷氨酰胺和尿苷；參與嘌呤代謝途徑的差異代謝物為L-谷氨酰胺和黃嘌呤；參與嘧啶代謝途徑的差異代謝物為尿苷。

圖7 番茄果實(shí)差異代謝物途徑Fig.7 Differential metabolite pathways in tomato fruits

表7 差異代謝物的代謝途徑歸屬Table 7 Metabolic pathway attribution of differential metabolites

3 討論

研究表明，容重在0.1～0.8 g·cm-3、總孔隙度在54%～96%、氣水比在0.25～0.50范圍內(nèi)的基質(zhì)最適宜植物生長［20］。本研究所有處理的容重和總孔隙度均滿足上述要求。生物炭的施入整體顯著提高了蚯蚓原位堆肥中的有機(jī)質(zhì)和養(yǎng)分含量，原因可能是生物炭較高的比表面積和吸附能力能固持蚯蚓原位堆肥中的氮、磷、鉀離子，減少養(yǎng)分的流失，增加其可利用養(yǎng)分的總量，且生物炭能夠延緩養(yǎng)分的釋放，使肥效變緩，補(bǔ)充作物后期生長發(fā)育所需的養(yǎng)分［21-23］。但隨著生物炭施用量的增加，蚯蚓原位堆肥養(yǎng)分含量呈先上升后下降的趨勢，并且酶活性的變化趨勢與之相似，其潛在原因是過量生物炭會吸附更多酶分子，抑制酶促反應(yīng)進(jìn)行，導(dǎo)致酶活性降低［24-25］。

適量生物炭對蚯蚓原位堆肥中的細(xì)菌豐富度有顯著影響，且細(xì)菌α多樣性指數(shù)大于真菌，這歸因于生物炭結(jié)構(gòu)疏松多孔，自身呈弱堿性，配施蚯蚓原位堆肥能夠刺激微生物的代謝，促進(jìn)細(xì)菌群落多樣性的提升，這與程揚(yáng)等［26］的研究結(jié)果相吻合。赤子愛勝蚓腸道內(nèi)的細(xì)菌類群主要為放線菌門、綠彎菌門和厚壁菌門等，上述菌群也被認(rèn)為是蚯蚓堆肥中常見的優(yōu)勢微生物［27-28］。這與本研究的結(jié)論相似，但各處理間菌群的相對豐度存在差異。一方面，綠彎菌門是一種兼性厭氧菌，可通過反向三羧酸循環(huán)來固定碳、氧化硫化物和Fe2+等物質(zhì)，促進(jìn)土壤養(yǎng)分循環(huán)和碳固定［29］，添加生物炭提高了綠彎菌門的相對豐度，這表明生物炭能增加蚯蚓原位堆肥的碳匯，提升有機(jī)碳的固持能力。另一方面，生物炭的添加降低了變形菌門和厚壁菌門的相對豐度。變形菌門作為細(xì)菌中最大的門類，包括如T.Escherich 1885、沙門氏菌等在內(nèi)的諸多病原菌，施加生物炭會抑制變形菌門活性，進(jìn)而降低相關(guān)病原菌的相對豐度［30］。此外，Cleveland等［31］發(fā)現(xiàn)土壤穩(wěn)定碳含量的增加會導(dǎo)致γ-變形菌門和厚壁菌門相對豐度的減弱。真菌類群中的優(yōu)勢菌門為子囊菌門，其次是被孢霉菌門和羅茲菌門。子囊菌門的主要功能是參與土壤有機(jī)質(zhì)分解，尤其具有分解木質(zhì)纖維素的能力，因此子囊菌門相對豐度的提高有利于有機(jī)質(zhì)的礦化和養(yǎng)分釋放［32］。被孢霉菌門中含有一些具有溶磷作用的真菌種群，可以溶解土壤中難以利用的磷元素［33］。此外，隨著生物炭施用量的增加，油壺菌門的相對豐度提高，這可能是由于生物炭中某種微量物質(zhì)的積累刺激了油壺菌門的繁殖，但具體原因有待進(jìn)一步研究。

蚯蚓原位堆肥中添加生物炭提高了番茄的產(chǎn)量和品質(zhì)。生物炭可加速蚯蚓原位堆肥的腐殖化過程，生物炭的表面也可以被堆肥中的微生物群落和腐殖質(zhì)氧化，兩者的協(xié)同作用增加了基質(zhì)中可用營養(yǎng)物質(zhì)的數(shù)量，促進(jìn)了番茄植株的生長，改善了番茄的產(chǎn)量和品質(zhì)［34］。然而，隨著生物炭添加比例的增大，番茄產(chǎn)量、可溶性固形物和可溶性糖含量呈先增大后減小趨勢，以C3處理最高。因此，生物炭對作物增產(chǎn)效應(yīng)存在峰值［35］。此外。本研究通過LC-MS測定番茄果實(shí)代謝物，結(jié)果表明，當(dāng)生物炭添加比例為3%時，KEGG通路富集產(chǎn)生了5條差異極顯著代謝途徑（P＜0.01），以氨酰tRNA生物合成代謝富集率最高。L-絲氨酸和L-谷氨酰胺參與了此過程，絲氨酸作為生糖氨基酸可參與糖異生途徑轉(zhuǎn)變?yōu)槠咸烟?，而L-絲氨酸又是由糖酵解中間體3-磷酸甘油酸酯通過三個催化步驟從頭合成的；谷氨酰胺可經(jīng)過代謝轉(zhuǎn)化為三羧酸循環(huán)的中間產(chǎn)物α-酮戊二酸，進(jìn)而實(shí)現(xiàn)氨基酸與糖類之間的轉(zhuǎn)化。此外，本研究還確定了尿苷-5'-二磷酸葡萄糖醛酸、1-己醇阿拉伯葡萄糖苷等差異代謝物。尿苷-5'-二磷酸葡萄糖醛酸可在UGlc AE酶的作用下轉(zhuǎn)化生成尿苷二磷酸半乳糖醛酸，進(jìn)入高爾基體內(nèi)參與果膠的生物合成。同時，它還可生成尿苷二磷酸木糖，間接形成木酮糖［36］。這些物質(zhì)的轉(zhuǎn)化可能是C3處理糖類物質(zhì)含量顯著高于CK的主要原因（表5）。綜上，生物炭的施用可能促進(jìn)了果實(shí)中氨基酸與糖類之間的轉(zhuǎn)化，從而提高了果實(shí)品質(zhì)。

4 結(jié)論

本研究結(jié)果表明，生物炭可提高蚯蚓原位堆肥中的養(yǎng)分含量和酶活性，優(yōu)化蚯蚓堆肥菌群結(jié)構(gòu)，進(jìn)而達(dá)到改善番茄品質(zhì)、提高番茄產(chǎn)量的目標(biāo)。其中，當(dāng)生物炭替代比例為3%時，蚯蚓原位堆肥中的全磷、速效磷、速效鉀含量及蔗糖酶和堿性磷酸酶活性最高。上述營養(yǎng)狀況的改變增加了細(xì)菌群落的豐富度和多樣性，促使綠彎菌門（Chloroflexi）相對豐度提高，變形菌門（Proteobacteria）、厚壁菌門（Firmicutes）和被孢霉門（Mortierellomycota）相對豐度降低。此外，LC-MS分析發(fā)現(xiàn)，生物炭通過改變番茄果實(shí)尿苷-5'-二磷酸葡萄糖醛酸、L-絲氨酸、L-谷氨酰胺等的代謝水平而影響果實(shí)中糖類物質(zhì)含量，從而保證其品質(zhì)與產(chǎn)量保持在高水平。