低溫環境下野外自然好氧堆肥理化特征及微生物群落結構變化

曾橋 施春陽 李俊 李祥 張璐璐

摘要 ［目的］研究低溫環境下野外自然好氧堆肥理化特征及微生物群落結構變化，實現低溫季節養殖廢棄物的資源化利用。［方法］以豬糞和稻殼為原料在低溫環境下進行自然好氧堆肥處理，對堆肥過程中溫度、pH、含水率、NH+4-N、NO-3-N和種子發芽指數（GI）進行分析，選取細菌16S rDNA V3～V4區和真菌ITS序列進行高通量測序，比較好氧堆肥過程中不同階段微生物群落多樣性及組成差異。［結果］當環境溫度為-6.7～6.4 ℃時，整個堆肥過程高溫階段持續11 d，經過38 d的好氧堆肥處理，堆體含水率、pH、NH+4-N、NO-3-N等指標均趨于穩定，堆肥結束時GI達到了122.12%，表明堆肥產品完全腐熟。低溫條件下堆肥不同階段樣本微生物群落分析表明，Acinetobacter是堆肥初期和升溫階段優勢細菌屬，而Pseudofulvimonas和Ensifer在堆肥高溫階段豐度最高，Parapedobacter為堆肥末期優勢細菌屬。Ascobolus是堆肥初期優勢真菌屬，而Mycothermus為好氧堆肥升溫、高溫和腐熟階段的優勢菌群。細菌群落功能預測結果表明，氨基酸轉運和代謝、碳水化合物運輸和代謝等是堆肥過程中細菌的主要代謝功能。理化指標與微生物多樣性相關性分析表明，堆肥過程中pH與屬水平上豐度較高的微生物相關性強，含水率、溫度、NH+4-N和發芽指數次之。［結論］低溫條件下野外自然好氧堆肥可順利進行并完全腐熟，這為低溫環境下野外自然好氧堆肥提供了依據。

關鍵詞 好氧堆肥；低溫環境；理化特征；高通量測序；群落結構；功能預測

中圖分類號 X713? 文獻標識碼 A? 文章編號 0517-6611（2024）06-0051-09

doi：10.3969/j.issn.0517-6611.2024.06.013

開放科學（資源服務）標識碼（OSID）：

Changes of Physico-chemical Characteristics and Microbial Community Structure During Natural Aerobic Composting Under Low Temperature Environment in the Field

ZENG Qiao1，2， SHI Chun-yang1，2， LI Jun3 et al

（1.School of Food and Biological Engineering，Shaanxi University of Science&Technology，Xian，Shaanxi 710021;2.Shaanxi Research Institute of Agricultural Products Processing Technology， Xian，Shaanxi 710021;3.Xian Agricultural Technology Extension Center，Xian，Shaanxi 710061）

Abstract ［Objective］To study the changes of physico-chemical characteristics and microbial community structure during natural aerobic composting under low temperature in the field，realize the resource utilization of breeding waste in low temperature season.［Method］Using pig manure and rice husk as raw materials for aerobic composting in the field， the temperature， pH， water content， NH+4-N， NO-3-N and germination index （GI） of seeds during composting process were analyzed， the bacterial 16S rDNA V3-V4 region and fungal ITS region were further selected for high-throughput sequencing.The diversity and composition of microbial communities in different stages of aerobic composting process were compared.［Result］When the environmental temperature was -6.7-6.4 ℃， the high-temperature stage of the entire composting process lasted for 11 days. After 38 days of aerobic composting treatment，the moisture content， pH， NH+4-N， NO-3-N and other indicators of the pile tended to be stable. At the end of composting， the GI reached 122.12%， indicating that the compost product was completely decomposed.The analysis of microbial communities in samples from different stages of composting under low temperature conditions showed that Acinetobacter was the dominant bacterial genus in the early and heating stages of composting， Pseudofulvimonas and Ensifer had the highest abundance in the high temperature stage of composting， and Parapedactor was the dominant bacterial genus in the late stage of composting. Ascobolus was the dominant fungal genus in the early stages of composting， while Mycothermus was the dominant microbial community in the heating， high temperature and maturation stages of aerobic composting.The predicted results of bacterial community function showed that amino acid transport and metabolism，carbohydrate transport and metabolism et al were the main metabolic functions of bacteria in composting process.The correlation analysis between physicochemical indexes and microbial diversity showed that the pH and genus abundance were highly correlated， followed by water content， temperature， NH+4-N and GI.［Conclusion］Natural aerobic compost can be carried out smoothly and completely decomposed in the field under low temperature conditions，it provide a basi for natural aerobic composting in the field under low temperature environment.

Key words Aerobic composting;Low temperature;Physico-chemical characteristics;High-throughput sequencing;Community structure;Function prediction

近年來，隨著我國人民生活水平的提高，對肉蛋奶等畜禽產品的需求量不斷增長，畜禽養殖業朝著集約化、工廠化、規模化方向迅速發展，從而導致畜禽糞便的排放量不斷增加。我國是世界上最大的畜禽養殖國，每年產生的畜禽糞便廢棄物高達40億t［1］，而規?；B殖場中有90%以上都缺乏必要的糞污治理措施，農村也存在大量未經處理的畜禽糞便被隨意處置和排放的現象，這不僅浪費了大量的養分資源，而且給環境帶來了極大的壓力。

好氧堆肥處理是實現畜禽糞便無害化和資源化利用的一個重要途徑，經過好氧堆肥處理，其發酵產物施用于土壤不僅可以增加土壤肥力，提高土壤保水、透氣的能力，還可以調節土壤中微生物菌群種類和數量，從而改善土壤結構，推進農業提質增效，同時避免了畜禽糞便低質量處理所導致的環境問題。野外自然好氧堆肥堆體適中，易于操作，不需要額外添加菌劑，適用于農村較為分散的農戶對在畜禽養殖和農業生產過程中產生的廢棄物進行資源化利用。環境溫度對好氧堆肥進程有較大的影響，環境溫度高，堆肥過程啟動快，堆肥效率也較高。我國地域遼闊，北方地區冬季時間長，溫度低，但是有關低溫環境下的野外自然好氧堆肥研究報道較少。

好氧堆肥的實質是利用微生物的生長和代謝活動轉化畜禽糞便中的有機物為無機物和穩定的腐殖質的過程［2］，微生物群落的代謝能力和菌群結構是影響好氧堆肥效果的關鍵因素［3］。堆肥過程中參與的微生物數量龐大，種類繁多，作用復雜，由于絕大部分微生物難以實現純培養，因此，無法采用傳統方法對其進行分離和鑒定。近年來，高通量測序技術被廣泛應用于堆肥環境中微生物多樣性分析，該技術可檢測到樣本中傳統純培養不能發現的低豐度微生物菌群，有助于更加全面準確地揭示環境中微生物群落的復雜度和多樣性［4］。該研究選取豬糞和稻殼作為堆肥原料，在冬季低溫條件下采用好氧堆肥工藝，研究野外堆肥過程中主要理化參數的動態變化，并對堆肥過程中微生物群落變化進行分析，從而深入了解低溫環境下堆肥過程中微生物群落的變化規律，為開發適宜農村分散的小堆體畜禽糞便冬季自然好氧堆肥化利用提供依據，從而促進畜禽糞便的資源化利用。

1 材料與方法

1.1 材料來源

新鮮豬糞取自興平市星光良種豬繁養殖有限公司，稻殼購自湖北京和米業有限公司。

1.2 試驗方法

1.2.1 好氧堆肥試驗設計。

試驗在冬季進行，環境溫度為-6.7～6.4 ℃，試驗地點為陜西西安。

將稻殼粉碎，過12目篩，按一定比例與新鮮豬糞混合均勻，調節碳氮比為30左右，調節堆體含水量為60%左右，堆體為圓錐體。每3 d翻堆1次，進入降溫階段后，每7 d翻堆1次，進入腐熟期停止翻堆，整個堆肥周期為38 d。

1.2.2 采樣方法。

從堆制0 d起，每2 d于堆體內部分別距離底部30 cm處和距離外部15～30 cm處的截面上多點共取混合樣品500 g（鮮重），樣品混勻后，采用四分法取樣，樣品分為兩份，一份用于測定理化指標，另一份貯存于-20 ℃冰箱備用。同時，于堆肥的0、4、10、16、38 d分別取樣，編號依次為1、2、3、4、5，貯存于-80 ℃冰箱用于高通量測序。

1.2.3 理化指標分析。

采用自動溫度記錄儀測定堆體溫度，將溫度記錄儀探頭插入堆體距離底部30 cm處和距離外部15～30 cm處的截面上隨機排布5個點，取平均溫度作為堆體實際溫度。取適量樣品于105 ℃下干燥至恒重，計算含水率。pH采用pH計測定。

硝態氮（NO-3-N）按NY/T 1116—2014規定的方法測定。

銨態氮（NH+4-N）采用靛酚藍比色法測定［5］。

種子發芽指數：取新鮮堆肥樣品3 g加入30 mL去離子水中，置于水平搖床上振蕩2 h后過濾，吸取濾液5 mL，加到鋪有濾紙的9 cm培養皿中。每個培養皿內點播20粒顆粒飽滿的小白菜種子，置于培養箱中于28 ℃培養，避光培養3 d，以去離子水作對照，按下式計算發芽指數［6］：

GI=（A×B）/（C×D）×100%（1）

式中：A為堆肥浸提液種子發芽率；B為堆肥浸提液種子根長；C為蒸餾水種子發芽率；D為蒸餾水種子根長。

1.2.4 微生物群落結構分析。

1.2.4.1 DNA的提取。

樣品DNA采用E.Z.N.ATM Mag-Bind Soil DNA Kit試劑盒（OMEGA，美國）按操作流程進行提取。

1.2.4.2 16S rDNA和ITS序列高通量測序。

采用特異性引物341F（CCTACGGGNGGCWGCAG）、805R（GACTACHVGGGTATCTAATCC）對細菌16S rDNA 基因的V3～V4可變區進行PCR擴增。采用特異性引物ITS1F（CTTGGTCATTTAGAGGAAGTAA）、ITS2R（GCTGCGTTCTTCATCGATGC）對真菌ITS rDNA基因的ITS1～ITS2區進行PCR擴增。PCR初次擴增條件：94 ℃ 3 min，94 ℃ 30 s，45 ℃ 20 s，65 ℃ 30 s，5個循環；94 ℃ 20 s，55 ℃ 20 s，72 ℃ 30 s，20個循環；72 ℃ 5 min。PCR二次擴增條件：95 ℃ 3 min，94 ℃ 20 s，55 ℃ 20 s，72 ℃ 30 s，5個循環；72 ℃ 5 min。進一步對DNA進行純化，利用Qubit3.0 DNA 檢測試劑盒對回收的DNA 精確定量，等量混合后送生工生物工程（上海）股份有限公司基于Illumina平臺進行高通量測序分析。

1.3 數據分析

理化指標測定相關試驗均重復3次，試驗結果均采用Graph Pad Prism軟件進行處理。高通量測序序列結果采用QIIME進行分析。序列質控過程將堿基平均質量低于20以及片段長度低于200 bp 的序列刪除［7］，同時，將含有模糊堿基的序列也進行刪除。采用“join paired ends.py”命令對兩端測序序列進行拼接。采用USEARCH61（Edgar 2010）比對 SILVA和UNITE數據庫進行嵌合體去除［8-9］。對去噪后序列采用QIIME進行多樣性分析，在97%的水平上進行分類操作單元（OTU）聚類［10］，并分別與SILVA和UNITE進行比對，確定微生物的分類信息?；赪eighted UniFrac采用“beta diversity through plots.py”在QIIME中進行主成分分析（PCA）。

2 結果與分析

2.1 好氧堆肥過程中溫度的變化

從圖1可以看出，整個堆肥過程環境溫度為-6.7～6.4 ℃，堆體溫度在6 d時達到50.7 ℃，在10 d時達到峰值54.3 ℃，至15 d開始，堆體溫度逐步下降，至17 d時，堆體溫度下降至50 ℃以下，且下降較快；從25 d起，堆體溫度下降緩慢，并逐步穩定。整個堆肥過程由于溫度的不同明顯分為3個階段，0～5 d為升溫階段，6～16 d為高溫階段，17～38 d為降溫腐熟階段。該試驗堆肥初期堆體升溫較慢，且堆體最高溫度較相關報道［11］的溫度低，可能的原因包括一方面此次堆肥試驗在冬季進行，環境溫度較低，進而造成堆體溫度也較低；另一方面，該試驗為自然堆肥，未添加外源菌劑，因此啟動較慢。從圖1可以發現，整個堆肥過程中，堆體溫度在50 ℃以上維持了11 d，符合國家《糞便無害化衛生要求》（GB 7959—2012）規定，達到了畜禽糞便無害化處理的要求，說明自然好氧堆肥在冬季較低環境溫度下可以順利進行。

2.2 好氧堆肥過程中含水率和pH的變化

有文獻報道，合適的堆肥初始含水率為50%～60%［12］。水分對于堆肥過程中微生物的生長繁殖過程必不可少，一方面，堆肥物料中的營養成分溶解于水中被微生物生長和代謝所利用，另一方面，堆肥過程中由于微生物的作用，堆體溫度逐步升高，水分的蒸發有利于調節堆體溫度，從而保證堆肥進程的順利進行。由圖2可知，堆肥初始含水率為58.69%，隨著堆肥過程的進行，含水率逐步降低，至20 d，含水率降至47.27%；從22 d 起，含水率下降速率大大降低，至38 d，含水率為43.90%，整個過程堆肥含水率下降了14.79百分點。該試驗含水率下降較慢，且堆肥結束時產品含水率也較高，這與張蘊琦等［13］采用牛糞和水稻秸稈在低溫下進行好氧堆肥的研究結果一致，主要是由于堆肥過程在冬季進行，環境溫度較低，因此，堆體水分蒸發較慢，從而導致產品含水率較高。好氧堆肥過程中pH從初始的7.02逐步上升至8.10，從6 d起，pH逐步下降，至16 d，pH逐步穩定，至堆肥結束時，pH為7.86。在堆肥前期，堆體pH迅速上升，這可能是由于堆體溫度上升，大量有機質被分解產生了大量NH+4-N；而隨著堆肥的進行，堆肥過程中產生的NH3大量揮發，以及微生物發酵產生了有機酸，從而使得堆體pH略下降，至堆肥結束，堆體逐步穩定，pH也趨于穩定。總體來看，冬季低溫環境下堆肥pH變化較其他季節溫和，這是由于在氣溫相對較高的季節微生物繁殖旺盛，堆肥過程相對劇烈。

2.3 好氧堆肥過程中NO-3-N、NH+4-N的變化

從圖3可以看出，NH+4-N含量從堆肥開始明顯上升，至4 d達到最大值，總體來說，NH+4-N含量呈先上升后下降的趨勢，這是由于隨著堆肥的進行，溫度逐漸升高，微生物活動加劇，同時由于礦化作用，因而產生了大量的NH+4-N。然而進入高溫階段后，由于溫度的升高以及翻堆，促使NH+4-N以NH3釋放出來，從而造成NH+4-N含量的下降。與圖1對比可以發現，NH+4-N含量的最高點在堆肥進入高溫階段前出現，這說明相關微生物在堆肥升溫階段已經開始大量生長和代謝，并通過脫氨作用將NH3釋放，從而使得NH+4-N上升。NO-3-N變化趨勢與NH+4-N完全相反，堆肥初期，NO-3-N含量先下降，而后逐漸上升，至20 d，NO-3-N含量達到峰值，此后NO-3-N含量緩慢下降并逐步保持穩定。堆肥過程中NO-3-N含量的變化與硝化細菌的活性密切相關［14］，高溫對硝化細菌的活性有一定的抑制作用，因此，NO-3-N含量的峰值較堆肥高溫階段滯后。至堆肥結束，NO-3-N和NH+4-N含量比值接近1，說明在低溫環境下進行的好氧堆肥已經達到腐熟。

2.4 好氧堆肥過程中發芽指數的變化

堆肥腐熟度是評價堆肥質量的重要指標，未腐熟的堆肥含有大量對植物種子有毒害作用的小分子有機物、NH3、揮發性脂肪酸和酚酸等［15］，從而抑制種子的生長。發芽指數是通過檢驗堆肥浸提液對植物種子發芽是否產生抑制作用來評價堆肥無害化、穩定化程度的指標，但是目前對于種子發芽指數達到多少表明堆肥已經腐熟，目前還尚需進一步研究確定。有研究認為，種子發芽指數大于50%則可以認為堆肥已達腐熟［16］，一般達到80%時，可認為堆肥已經完全腐熟［17］，也有學者認為，發芽指數達到110%是評價堆肥腐熟的較為合適的指標［18］。從圖4可以看出，在堆肥結束的38 d，種子發芽指數達到了122.12%。此外，NO-3-N/NH+4-N接近1，同時溫度、含水率、pH等指標均趨于穩定，可以判斷低溫環境下的野外自然堆肥產品完全腐熟，這為冬季低溫環境下堆肥提供了依據。

2.5 低溫環境下堆肥微生物群落組成分析

2.5.1 堆肥不同階段樣本α多樣性分析。

α多樣性可以反映微生物群落的豐度和多樣性，Chao指數和Ace指數越大，說明微生物群落豐富度越高；Shannon指數越大，說明群落多樣性越高；Simpson指數越大，說明群落多樣性越低；Coverage指數實際反映了測序結果是否代表樣本的真實情況，其數值越高，則樣本中序列沒有被測出的概率越低［19］。堆肥各階段樣品細菌和真菌群落α多樣性分析結果見圖5，從圖5可以看出，Coverage指數均在0.93以上，說明測序結果能較真實地反映樣品中微生物群落［20］；由Chao指數和Ace指數可知，細菌和真菌群落豐富度變化趨勢一致，均呈下降—上升—下降—上升。由Shannon指數可知，細菌群落多樣性在升溫和高溫階段呈逐步上升趨勢，高溫階段末期和腐熟階段呈下降趨勢；而真菌群落多樣性在升溫階段先上升，至高溫階段下降，在高溫末期則又上升，腐熟階段多樣性下降。α多樣性分析表明，豬糞堆肥樣品的微生物多樣性初期和升溫階段是最豐富的，而田偉［18］以牛糞為原料堆肥時，腐熟后的堆肥樣品中的微生物多樣性最豐富，說明不同的氮源對微生物多樣性有較大影響。

2.5.2 堆肥過程中微生物群落結構。

2.5.2.1 細菌群落結構。

從圖6可以看出，Acinetobacter是堆肥初期優勢細菌屬，其豐度為65.91%，其次豐度較高的細菌屬為Pseudomonas、Psychrobacter、Devosia、Brevundimonas、Comamonas，其豐度為1.19%～9.74%。隨著堆肥的進行，堆體溫度逐漸升高，在升溫階段，Acinetobacter豐度下降至33.2%，但仍為豐度最高的細菌屬，Devosia、Leadbetterella、Pseudomonas、Brevundimonas、Stenotrophomonas、Comamonas、Gemmobacter等為豐度次之的細菌屬，其豐度為1.65%～8.40%。隨著堆肥進入高溫階段，Acinetobacter的豐度迅速下降至1.54%，Pseudofulvimonas豐度最高，達19.95%，Ensifer豐度次之，為4.25%，Verrucomicrobium的豐度為3.76%，其余豐度較高的細菌屬為Devosia、Luteimonas、Pseudoxanthomonas、Parapedobacter、Solibacillus、Kofleria、Povalibacter等，其豐度為1.84%～2.88%。高溫階段末期，Pseudofulvimonas、Ensifer、Parapedobacter、Gp6、Luteimonas豐度均略有升高，累計豐度達37.47%；Verrucomicrobium豐度略有下降，為2.38%；Bordetella、Devosia、Povalibacter、Pseudoxanthomonas的豐度為1.76%～2.27%；Parapedobacter為堆肥末期優勢細菌屬，豐度為12.37%，Pseudofulvimonas、Pseudomonas、Luteimonas、Rhodanobacter、Phaselicystis、Arachidicoccus、Ensifer等豐度為1.69%～3.43%。

從圖6b可以發現，整個好氧堆肥過程中，總體豐度較高的細菌屬分別為Acinetobacter、Pseudofulvimonas、Parapedobacter、Pseudomonas、Devosia、Ensifer、Luteimonas、Verrucomicrobium和Gp6等，它們是堆肥過程中主要的細菌屬。Acinetobacter主導升溫階段，Pseudofulvimonas主導高溫階段和降溫階段，Parapedobacter主導腐熟階段。

2.5.2.2 真菌群落組成。

從圖7a可以看出堆肥過程中真菌群落在屬水平上的分布，Ascobolus是堆肥初期優勢真菌屬，其豐度為18.53%，其次為Verticillium，豐度為6.65%，其余豐度在1.00%以上的真菌屬包括Alternaria、Papiliotrema。在堆肥升溫階段，真菌群落結構發生了較為明顯的變化，Ascobolus豐度下降很快，下降幅度達84%，但仍然為升溫階段豐度較高的真菌屬，其余豐度較高的包括Arthrographis、Fusarium、Anthopsis、Cystofilobasidium、Cryptococcus、Alternaria等，其豐度為0.36%～2.02%。隨著堆肥過程的進行，Mycothermus成為好氧堆肥高溫階段的絕對優勢菌群，其豐度為72.75%，Anthopsis豐度略高于升溫階段，達到3.48%，其余豐度較高的分別為Candida、Aspergillus、Verticillium、Ascobolus、Mortierella等，豐度為0.32%～0.80%。高溫階段后期，Mycothermus豐度略有下降，但仍然為優勢菌群，Anthopsis豐度上升至4.94%，Verticillium、Candida、Aspergillus、Alternaria、Archaeorhizomyces、Mortierella、Ascobolus等豐度為0.46%～1.66%。在堆肥腐熟階段，Mycothermus豐度進一步下降，但仍然為優勢菌群，至好氧堆肥末期，其豐度下降至45.10%；Mortierella豐度上升較快至3.80%，Chrysosporium、Myceliophthora、Verticillium、Aspergillus等豐度為0.39%～0.72%。

從圖7b可以發現，好氧堆肥過程中，總體豐度較高的真菌屬為Mycothermus、Ascobolus、Anthopsis、Verticillium、Mortierella等，其中Ascobolus是堆肥初期和升溫階段主導真菌屬，Mycothermus主導高溫階段、降溫和腐熟階段。

2.5.3 堆肥不同階段樣本PCA和OTU分析。

基于細菌門水平相對豐度PCA分析結果（圖8a）表明，PCA1和PCA2的主要貢獻之和為90%，表明在這2個維度上能夠反映樣本的實際情況。從圖8a可以看出，在堆肥過程中細菌群落結構發生了明顯的變化，不同階段樣本發生了明顯的分離，其中樣本3和樣本4距離較近，說明菌群結構相似。這可能是由于樣本3和樣本4分別為堆肥第10天和第16天，堆體溫度差異不大，因而菌群結構也較相似，隨著堆肥逐步進入腐熟階段，堆體溫度下降較快，細菌群落結構變化也較大。利用Venn圖研究不同階段樣品中細菌分布的變化，結果如圖8b所示，5個樣本中共檢出16 511個OTUs，其中僅114個OTUs是共有的，占總OTUs數目的0.69%。其中，堆肥初期樣本中特有OTUs數目為3 103個，占總OTUs數目的18.79%；升溫階段樣本中特有OTUs數目為3 183個，占總OTUs數目的19.28%；高溫期和高溫末期特有的OTUs數目分別為2 898和2 642個，分別占總OTUs數目的17.55%和16.00%；堆肥末期樣本中特有OTUs數目為3 679個，占總OTUs數目的22.28%。堆肥各不同階段樣本中特有OTUs數目總計15 505個，共占總OTUs數目的93.90%。該結果表明，在低溫環境下好氧堆肥不同階段細菌群落結構發生了明顯的變化。

基于真菌門水平相對豐度PCA分析結果（圖9a）表明，PCA1和PCA2的主要貢獻之和為88%，表明在這2個維度上能夠反映樣本的實際情況。從圖9a可以看出，在堆肥過程中真菌群落結構發生了明顯的變化，不同階段樣本發生了明顯的分離，其中樣本3和樣本4距離較近，說明菌群結構相似。這可能是由于樣本3和樣本4分別為堆肥第10天和第16天，堆體溫度差異不大，因而菌群結構也較相似，隨著堆肥逐步進入腐熟階段，堆體溫度下降較快真菌群落結構變化也較大。利用Venn圖研究不同階段樣品中真菌分布的變化，結果如圖9b所示，5個樣本中共檢出3 498個OTUs，其中26個OTUs是共有的，占總OTUs數目的0.74%，堆肥初期、升溫階段、高溫期、高溫末期、堆肥末期樣本中特有OTUs數目分別為584、436、691、603和734個，共占總OTUs數目的87.14%。該結果表明，在低溫環境下好氧堆肥不同階段真菌群落結構發生了明顯的變化。

2.5.4 堆肥過程中細菌群落功能變化。

從圖10可以看出，在好氧堆肥初期，細菌菌群對應的主要功能為氨基酸轉運和代謝（amino acid transport and metabolism），轉錄（transcription），能量產生與轉化（energy production and conversion），翻譯、核糖體結構和合成（translation，ribosomal structure and biogenesis），細胞壁、細胞膜和細胞殼的合成（Cell wall/membrane/envelope biogenesis），無機離子轉運和代謝（inorganic ion transport and metabolism）等。而在升溫階段，主要功能為氨基酸轉運和代謝，轉錄，能量產生與轉化，碳水化合物運輸和代謝（carbohydrate transport and metabolism），細胞壁、細胞膜和細胞殼的合成，無機離子代謝，信號傳導機制（singal transduction mechanisms），翻譯、核糖體結構和合成等。高溫和降溫階段，細菌菌群主要代謝功能與升溫階段類似?？傮w來看，氨基酸轉運和代謝、碳水化合物運輸和代謝等是堆肥過程中細菌的主要代謝功能，堆肥原料中豐富的纖維素、蛋白質、脂肪等通過氨基酸轉運和代謝、碳水化合物運輸和代謝所產生的氨基酸、大分子糖類、多酚、多肽等是腐殖酸合成的前體物質，具有促進腐殖酸合成和提高腐殖酸產量及分子量的作用，從而加快堆肥進程和促進堆肥的腐熟［21-22］。

2.6 微生物多樣性與理化指標相關性分析

將堆肥過程中理化指標與豐度較高的細菌屬和真菌屬進行相關性分析，結果如表1所示。pH與Acinetobacter、Pseudomonas、Luteimonas、Pseudoxanthomonas、Ascobolus、Verticillium、Anthopsis高度相關（|r|>0.8），且與Acinetobacter、Ascobolus和Verticillium呈顯著負相關（P<0.05），與Luteimonas呈顯著正相關（P<0.05）；含水率與Acinetobacter、Pseudofulvimonas、Pseudoxanthomonas、Parapedobacter、Povalibacter、Verticillium、Anthopsis、Mycothermus高度相關（|r|>0.8），其中與Anthopsis呈顯著負相關（P<005）；溫度與Acinetobacter、Pseudofulvimonas、Pseudoxanthomonas、Parapedobacter、Povalibacter、Mycothermus高度相關（|r|>0.8），與Parapedobacter呈顯著負相關（P<0.05），與Mycothermus則呈顯著正相關（P<0.05）；NH+4-N與Pseudofulvimonas、Pseudoxanthomonas、Parapedobacter、Povalibacter、Verticillium、Anthopsis、Mycothermus高度相關（|r|>0.8），與Pseudoxanthomonas呈極顯著正相關（P<0.01），與Anthopsis呈顯著負相關（P<0.05）；NO-3-N除與Pseudoxanthomonas高度相關外，與其他屬微生物相關性較弱；發芽指數與Acinetobacter、Pseudofulvimonas、Pseudoxanthomonas、Parapedobacter、Povalibacter、Anthopsis、Mycothermus高度相關（|r|>0.8），與Anthopsis呈顯著正相關（P<0.05）。

3 結論

（1）低溫環境下野外自然好氧堆肥啟動慢，高溫階段溫度相對較低，但經過38 d好氧堆肥過程，堆體各項理化特征指標趨于穩定，種子發芽指數達到了122.12%，說明低溫環境下堆肥樣品可以達到腐熟。

（2）堆肥升溫階段細菌和真菌多樣性較為豐富，細菌和真菌多樣性變化存在一定差異。堆肥過程中微生物群落結構發生了明顯變化，細菌屬水平上優勢菌群隨著堆肥過程的進行而發生改變，但真菌屬水平上優勢菌群在進入高溫階段后至堆肥結束主要為Mycothermus。氨基酸轉運和代謝、碳水化合物運輸和代謝等是堆肥過程中細菌的主要代謝功能，它們加快了堆肥的進程和促進了堆肥的腐熟。

（3）理化指標與微生物多樣性相關性分析表明，堆體pH、含水率、溫度、NH+4-N和發芽指數等因素與多種屬水平上豐度較高的微生物存在顯著相關性（P<0.05），而NO-3-N與屬水平上主要菌群相關性均不顯著（P>0.05）。

參考文獻

［1］ 閻波杰，趙春江，潘瑜春，等.大興區農用地畜禽糞便氮負荷估算及污染風險評價［J］.環境科學，2010，31（2）：437-443.

［2］ 劉立，張朝升，趙美花.生物炭對好氧堆肥化處理影響的研究進展［J］.環境科學與技術，2018，41（9）：170-175，182.

［3］ 王秀紅，李欣欣，史向遠，等.好氧堆肥微生物代謝多樣性及其細菌群落結構［J］.環境科學研究，2018，31（8）：1457-1463.

［4］ 牛世全，龍洋，李海云，等.應用Illumina MiSeq 高通量測序技術分析河西走廊地區鹽堿土壤微生物多樣性［J］.微生物學通報，2017，44（9）：2067-2078.

［5］ 郭文淼，辛宇，張金堯，等.全波長掃描式多功能讀數儀－靛酚藍比色法測定水樣銨態氮含量［J］.中國土壤與肥料，2018（4）：166-170.

［6］ 王若斐，劉超，操一凡，等.不同碳氮比豬糞堆肥及其產品肥效［J］.中國土壤與肥料，2017（6）：127-134，154.

［7］ CAPORASO J G，KUCZYNSKI J，STOMBAUGH J，et al.QIIME allows analysis of high-throughput community sequencing data［J］.Nat Methods，2010，7（5）：335-336.

［8］ QUAST C，PRUESSE E，YILMAZ P，et al.The SILVA ribosomal RNA gene database project：Improved data processing and web-based tools［J］.Nucleic Acids Res，2013，41：D590-D596.

［9］? EDGAR R C.Search and clustering orders of magnitude faster than BLAST［J］.Bioinformatics，2010，26（19）：2460-2461.

［10］ ZGADZAJ R，GARRIDO-OTER R，JENSEN D B，et al.Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere，root，and nodule bacterial communities［J］.Proc Nat Acad Sci USA，2016，113（49）：E7996-E8005.

［11］ 曹云，黃紅英，孫金金，等.超高溫預處理對豬糞堆肥過程碳氮素轉化與損失的影響［J］.中國環境科學，2018，38（5）：1792-1800.

［12］ ABBASI T，GAIALAKSHMI S，ABBASI S A，et al.Towards modeling and design of vermicomposting systems：Mechanisms of composting/vermicomposting and their implications［J］.Indian J Biotechnol，2009，8：177-182.

［13］ 張蘊琦，徐鳳花，張晴，等.低溫環境下水稻秸稈與牛糞混合堆肥效果研究［J］.安徽農業科學，2016，44（29）：112-114，126.

［14］ 楊萍萍，尹華，彭輝，等.外接菌種對污泥堆肥效能及堆體細菌群落的影響［J］.環境科學，2017，38（8）：3536-3543.

［15］ KIRCHMANN H，WIDN P.Separately collected organic household wastes.Chemical composition and composting characteristics［J］.Swed J Agr Res，1994，24（1）：3-12.

［16］ ZUCCONI F，MONACO A，FORET M，et al.Phytotoxins during the stabilization of organic matter［M］//GASSER J K R.Compositing of agricultural and other wastes.New York：Elsevier Applied Science Publishers，1985：73-85.

［17］ 徐江兵，王艷玲，羅小三，等.納米氧化鋅對堆肥過程中細菌群落演替的影響［J］.應用與環境生物學報，2017，23（6）：1166-1171.

［18］ 田偉.牛糞高溫堆肥過程中的物質變化、微生物多樣性以及腐熟度評價研究［D］.南京：南京農業大學，2012：39.

［19］ 黃雅楠，王曉慧，曹琦，等.高通量測序技術分析豬糞堆肥過程中微生物群落結構變化［J］.微生物學雜志，2018，38（5）：21-26.

［20］ 任省濤，郭夏麗，蘆阿虔，等.林可霉素菌渣堆肥微生物群落多樣性分析［J］.環境科學，2018，39（10）：4817-4824.

［21］ ZHANG Y C，YUE D B，MA H.Darkening mechanism and kinetics of humification process in catechol-Maillard system［J］.Chemosphere，2015，130：40-45.

［22］ HARDIE A G，DYNES J J，KOZAK L M，et al.The role of glucose in abiotic humification pathways as catalyzed by birnessite［J］.Journal of molecular catalysis A：Chemical，2009，308（1/2）：114-126.