馮丹妮 伍鈞 楊虎德

摘要：采用PCR-DGGE克隆測(cè)序技術(shù)，分析沼液連續(xù)施用水稻油菜輪作耕層土壤細(xì)菌群落結(jié)構(gòu)、豐富度指數(shù)和多樣性指數(shù)等。通過(guò)DGGE圖譜分析得出，沼液連年施用不利于豐富土壤細(xì)菌群落多樣性，當(dāng)連續(xù)3 a沼液施用總量控制為339.31～396.81 t/hm2時(shí)，細(xì)菌群落多樣性增加，種群功能豐富、信息復(fù)雜、種群穩(wěn)定，過(guò)高或過(guò)低施用沼液都會(huì)顯著降低細(xì)菌群落多樣性;通過(guò)UPGMA分析表明，長(zhǎng)期連續(xù)施用沼液后土壤細(xì)菌遺傳相似性降低，細(xì)菌種群結(jié)構(gòu)變化越來(lái)越大。

關(guān)鍵詞：水稻;油菜;輪作;沼液施用量;土壤細(xì)菌多樣性;PCR-DGGE

中圖分類號(hào)：S154.32? ? ? 文獻(xiàn)標(biāo)志碼：A? ? ? 文章編號(hào)：1001-1463（2019）08-0041-09

Abstract：PCR-DGGE cloning sequencing technology was used to analyze and compare the bacterial community structure，richness index and diversity index of soil bacteria in rice-rape rotation cultivation layer. Through the analysis of the DGGE profiles concluded that biogas slurry in successive years application does not favor the rich soil bacterial community diversity， when biogas slurry application of total amount controled in 339.31～396.81 t/hm2， bacterial community diversity increases， bacterial population feature rich， complicated information and population stable， too high or too low application of biogas slurry can significantly reduce the bacterial community diversity; According to UPGMA analysis， the genetic similarity of soil bacteria decreased after long-term continuous application of biogas slurry， and the bacterial population structure changed bigger and bigger.

Key words：Rice;Rape;Rotation;Amount of biogas slurry;Bacterial community diversity;PCR-DGGE

隨著規(guī)模化養(yǎng)殖業(yè)迅猛發(fā)展，畜禽糞便等廢棄物排放成為突出的環(huán)境問(wèn)題，沼氣工程是處理畜禽糞便等廢棄物最有效的途徑，為此國(guó)家增加投入建設(shè)沼氣工程，隨之產(chǎn)生大量的沼液。沼液中含有氮、磷、鉀等作物生長(zhǎng)必須的營(yíng)養(yǎng)元素和鐵、鋅、錳等多種微量元素[1 ]，同時(shí)富含各類氨基酸、維生素、蛋白質(zhì)、生長(zhǎng)素等有益的有機(jī)成分，對(duì)提高農(nóng)產(chǎn)品品質(zhì)、改良土壤有積極作用[2 - 4 ]。但現(xiàn)代規(guī)?；B(yǎng)殖為加快畜禽生長(zhǎng)速度及提高飼料利用率和預(yù)防疾病，常在飼料中添加激素類、抗生素等物質(zhì)，它們大部分殘留于畜禽糞便中，對(duì)沼液的資源化利用提出了新的挑戰(zhàn)，不合理地施用沼液可能提高農(nóng)業(yè)面源污染風(fēng)險(xiǎn)。

微生物作為土壤生態(tài)系統(tǒng)中最活躍的部分，受土壤環(huán)境條件影響較為敏感。越來(lái)越多的研究表明，土壤微生物多樣性在評(píng)價(jià)土壤生態(tài)系統(tǒng)功能、 維護(hù)土壤生態(tài)系統(tǒng)平衡及促進(jìn)土壤可持續(xù)利用等方面發(fā)揮著重要作用[5 - 6 ]。近年來(lái)，分子生物學(xué)分析方法不斷進(jìn)步，并被廣泛應(yīng)用于微生物生態(tài)學(xué)研究，使土壤微生物多樣性研究有了突破式的進(jìn)展。多聚酶擴(kuò)增-變性梯度電泳法（PCR-DGGE）在土壤微生物多樣性分析中應(yīng)用的較為成熟，可以較好的反映土壤微生物的結(jié)構(gòu)變化[7 - 8 ]。我們于2009 — 2012年以資源化、無(wú)害化處理沼液，防止其造成面源污染為目的，在四川南方山丘區(qū)水稻油菜輪作土壤中進(jìn)行連續(xù)3 a施用沼液試驗(yàn)，通過(guò)監(jiān)測(cè)0～20 cm耕層土壤細(xì)菌群落多樣性，評(píng)價(jià)長(zhǎng)期施用沼液對(duì)土壤細(xì)菌群落多樣性的影響，以期為沼液的農(nóng)田消解方式提供更全面的參考和理論依據(jù)。

1? ?材料與方法

1.1? ?試驗(yàn)地概況

試驗(yàn)設(shè)在四川省邛崍市固驛鎮(zhèn)黑石村，地勢(shì)平坦向陽(yáng)，排灌方便。供試土壤為黃壤性水稻土。原始耕層土壤肥力中等，含全氮2.11 g/kg、堿解氮173.4 mg/kg、有機(jī)質(zhì)35.65 g/kg、速效磷57.32 g/kg、速效鉀31.48 mg/kg，pH為4.8。

1.2? ?供試材料

指示水稻品種為秈型水稻，油菜品種為宜油15。供試沼液為四川金利有限公司生豬黑石養(yǎng)殖場(chǎng)以通威飼料養(yǎng)殖的豬場(chǎng)糞尿經(jīng)厭氧發(fā)酵完全的沼液，其成分含量為全氮876.23 mg/L、銨態(tài)氮 634.7 mg/L、硝態(tài)氮 7.27 mg/L、全磷83.19 mg/L、全鉀378.5 mg/L。供試化肥為尿素（含N≥46%）、含有機(jī)質(zhì)磷肥（速效磷≥12%，有機(jī)質(zhì)≥3.0%）、氯化鉀（K2O≥60%）、硼肥（硼酸鈉鹽含量99%，純硼含量15%）、有機(jī)復(fù)合肥（N- P2O5-K2O為20-20-10）。

1.3? ?試驗(yàn)設(shè)計(jì)

試驗(yàn)采用水稻油菜輪作。自2009年播種水稻至2012年收割油菜，連續(xù)3 a種植6季作物（水稻油菜輪作），每季作物根據(jù)其生長(zhǎng)所需養(yǎng)分及以往施肥經(jīng)驗(yàn)確定不同的沼液施用量。設(shè)置12個(gè)處理（表1），包括1個(gè)清水對(duì)照（CK1）、1個(gè)常規(guī)施肥對(duì)照（CK2）和10個(gè)純沼液處理（BS1～BS10）。隨機(jī)區(qū)組排列，3次重復(fù)，小區(qū)面積20 m2（5 m×4 m），處理間間隔40 cm，重復(fù)間間隔50 cm，四周設(shè)保護(hù)行。水稻、油菜種植均按當(dāng)?shù)爻Ｒ?guī)方法育苗，水稻4月播種，5月移栽，按寬窄行移栽，每小區(qū)定植288穴;油菜9月份播種，10月份移栽，株距35 cm，每小區(qū)種植168株。

1.4? ?樣品采集

連續(xù)3 a 6季水稻油菜輪作，油菜收割后分小區(qū)多點(diǎn)采集0～20 cm土樣，混勻，四分至1 kg左右，用滅菌鑷子去除植物根系、石礫等雜物，將土樣收集于無(wú)菌封口袋中，密封帶回實(shí)驗(yàn)室儲(chǔ)存于-70 ℃冰箱備用。

1.5? ?分析方法

1.5.1? ? 土壤總DNA提取? ? 采用美國(guó)MOBIO公司的PowerSoilTM DNA Isolation Kit（DNA 強(qiáng)力提取試劑盒）提取，依照試劑盒設(shè)計(jì)，土壤加樣量設(shè)定為0.25 g，提取完畢用1%的瓊脂糖凝膠電泳進(jìn)行質(zhì)量檢測(cè)。

1.5.2? ? PCR-DGGE分析? ? PCR擴(kuò)增：選用細(xì)菌16S rDNA V6-V8區(qū)通用引物（上海生工生物工程技術(shù)服務(wù)有限公司合成），具體引物名稱與序列見(jiàn)表2[9 - 10 ]。PCR擴(kuò)增采用Nested-PCR（巢式PCR）程序。為減少擴(kuò)增過(guò)程中產(chǎn)生非特異性產(chǎn)物，試驗(yàn)采用Touch- down反應(yīng)程序。

第1次PCR反應(yīng)體系：PCR Master Mix （TIANGEN BIOTECH， BEIJING） 25 μL，引物BSF8/20和BSR1541/20各0.5 μL，DNA模板8 ng，加ddH2O至50 μL。反應(yīng)程序：94 ℃變性10 min，隨后94 ℃變性1 min、65～55 ℃退火50 s（每次循環(huán)降低0.5 ℃）、72 ℃延伸90 s，循環(huán)20次，而后94 ℃變性1 min、55 ℃退火50 s、72 ℃延伸90 s，15個(gè)循環(huán)，最后在72 ℃條件下再延伸10 min，產(chǎn)物保存于4 ℃。

第2次PCR擴(kuò)增：將第1輪PCR擴(kuò)增產(chǎn)物1∶10稀釋液作為DNA模板。體系為：25 μL PCR Master Mix，引物GC-968F和1401R各0.5 μL，DNA模板1 μL，ddH2O 23 μL（總體系為50 μL）。反應(yīng)體系在第1次PCR程序基礎(chǔ)上，退火溫度20個(gè)循環(huán)改為63～53 ℃、15個(gè)循環(huán)改為53 ℃，其他條件不變。

1.5.3? ? 變性凝膠梯度電泳（DGGE）? ? 凝膠制備：按DGGE濃度8%、變性劑高低濃度40%～55%，并分別加入25 μL TEMED、50 μL 10%過(guò)硫酸銨溶液，搖勻，迅速手動(dòng)灌膠，插入梳子，置于50 ℃烘箱30 min，而后置于水平桌面過(guò)夜。

1.5.4? ? 電泳? ? 輕輕拔出梳子并用緩沖液清洗點(diǎn)樣孔，安裝膠板并置于裝有預(yù)熱（60 ℃）的1×TAE緩沖液的電泳槽中，按電壓50 V電泳30 min，再按電壓180 V電泳7 h。

1.5.5? ? 染色與成像? ? 電泳結(jié)束之后，采用銀染法對(duì) DGGE 凝膠進(jìn)行染色。其主要流程如下：①固定：將凝膠置于1 L固定液中（10%無(wú)水乙醇，0.5%冰醋酸），水平振蕩15 min;②染色：取出膠片用去離子水沖洗凝膠2～3次，然后加入1 L銀染液（0.2% AgNO3，200 μL 甲醛），避光水平振蕩20 min;③顯色：倒棄銀染液用去離子水沖洗凝膠2～3次，加入1 mL顯色液（1.5% NaOH，2 mL甲醛），振蕩至出現(xiàn)較清楚條帶;④再固定：倒棄顯色液放入1 L固定液（可重復(fù)使用）固定10 s;⑤清洗膠片。

成像：將膠片平整放在膠片觀察燈上，用相機(jī)清晰拍下膠片。

1.5.6? ? DGGE圖譜分析? ? DGGE凝膠上的每個(gè)條帶都代表一個(gè)單獨(dú)的序列類型或系統(tǒng)發(fā)育類型。采用所得圖像用Bio-Rad Quantuty One 4.62 對(duì)DGGE圖譜進(jìn)行定量分析，如條帶位置、相對(duì)亮度等。采用Excel 2007計(jì)算 Simpson 指數(shù)、 Shannon-Wiener 指數(shù)和均勻度指數(shù)（EH）。

式中Ni為第i種物種的光密度值，N為所有條帶的總光密度值。

根據(jù)條帶的出現(xiàn)或消失來(lái)建立二元矩陣，經(jīng) NTSYS 軟件計(jì)算得出 Jaccard 系數(shù)，并用 UPGMA（Unweighted pair group method using arithmetic averages）進(jìn)行聚類分析。

1.6? ?數(shù)據(jù)處理

將Quantity One定量分析數(shù)據(jù)拷貝到Excel 2007中進(jìn)行轉(zhuǎn)換運(yùn)算，再由SPSS18.0進(jìn)行統(tǒng)計(jì)分析，用鄧肯（Duncan’s test）檢測(cè)法檢驗(yàn)其差異顯著性水平（P < 0.01）。

2? ?結(jié)果與分析

2.1? ?土壤微生物總DNA提取與PCR擴(kuò)增

采用DNA強(qiáng)力提取試劑盒提取土壤微生物總DNA，所提取DNA溶液濃度為6.61～14.11 ng/μL（表3）。第1輪PCR擴(kuò)增應(yīng)用細(xì)菌通用引物BSF8/20和BSR1541/20獲取DNA模板，將其稀釋10倍，取1 μL稀釋液于50 μL體系，采用GC-968F和1401R引物對(duì)進(jìn)行第2輪擴(kuò)增得到目標(biāo)片段，取3 μL第二輪PCR產(chǎn)物，經(jīng)1%瓊脂糖凝膠電泳檢測(cè)， PCR產(chǎn)物片段大小均勻，特異性高（圖1），無(wú)需純化適合進(jìn)行DGGE操作。

2.2? ?DGGE圖譜分析

根據(jù)DGGE原理，DGGE圖譜中每個(gè)條帶大致代表一種微生物優(yōu)勢(shì)菌群，條帶染色后，條帶強(qiáng)度代表該菌種相對(duì)豐度，條帶數(shù)多寡則代表細(xì)菌種群多樣性高低。由圖2可以看出，不同處理樣品泳道間的公共條帶代表供試樣品中存在共有細(xì)菌，說(shuō)明該種細(xì)菌類群穩(wěn)定，不易受土壤環(huán)境變化影響;不同處理樣品泳道DGGE條帶數(shù)量、遷移位置、亮度和粗細(xì)均存在差異，說(shuō)明不同處理樣品中細(xì)菌類群種類和數(shù)量差異較大。與CK1相比，CK2樣品中部分條帶亮度減弱或消失;沼液施用處理樣品中雖有部分條帶消失，同時(shí)也有新的條帶出現(xiàn)，隨沼液施用量增加條帶位置和亮度逐步發(fā)生明顯變化。

2.3? ?細(xì)菌群落結(jié)構(gòu)分析

應(yīng)用 Quantity One 4.62對(duì) DGGE 圖譜進(jìn)行數(shù)字化處理，計(jì)算不同處理樣品泳道的豐富度指數(shù)、Shannon-Wiener 指數(shù)、Simpson 指數(shù)和均勻度指數(shù)等。

從表4可以看出，清水處理（CK1）土壤細(xì)菌豐富度為30、Shannon-Wiener 指數(shù)3.390，極顯著高于化肥處理和多數(shù)BS處理，說(shuō)明清水處理土壤中細(xì)菌類群繁多、功能性強(qiáng)。長(zhǎng)期施用化肥的CK2處理土壤中細(xì)菌類群豐富度和多樣性指數(shù)極顯著小于清水對(duì)照處理。長(zhǎng)期施用沼液處理土壤細(xì)菌類群也發(fā)生明顯變化，隨沼液施用量增加，土壤細(xì)菌種群豐富度、多樣性指數(shù)大體呈先升高后降低的趨勢(shì)，并以BS4處理（沼液施用總量396.81 t/hm2）達(dá)最大值，其中豐富度為32、多樣性指數(shù)為3.470，極顯著大于CK1、CK2和其他BS處理（BS3除外）。當(dāng)沼液施用量低于396.81 t/hm2時(shí)（即BS1、BS2處理），土壤細(xì)菌豐富度和多樣性指數(shù)極顯著小于CK1、大于CK2，高于該施用量（即BS5、BS6、BS8、BS9、BS10）的處理細(xì)菌種群豐富度和多樣性逐漸降低且極顯著小于CK1、CK2，說(shuō)明長(zhǎng)期不合理的施用沼液會(huì)導(dǎo)致土壤細(xì)菌豐富度和多樣性指數(shù)減小。

各處理土壤 Simpson 指數(shù)呈完全相反變化規(guī)律。當(dāng)土壤細(xì)菌 Shannon-Wiener指數(shù)較高時(shí)，Simpson 指數(shù)較低;Shannon- Wiener 指數(shù)降低時(shí) Simpson 指數(shù)升高。說(shuō)明土壤細(xì)菌類群繁多時(shí)優(yōu)勢(shì)類群不突出，土壤細(xì)菌類群數(shù)量減少時(shí)優(yōu)勢(shì)類群明顯。

不同處理?xiàng)l件下細(xì)菌種群均勻度差異不顯著，說(shuō)明連續(xù)3 a施肥處理存活的土壤細(xì)菌類群通過(guò)群落結(jié)構(gòu)的不斷變化來(lái)優(yōu)化群落的整體功能，從而適應(yīng)環(huán)境條件，降低物種奇異度。

2.4? ?細(xì)菌群落相似性分析

通過(guò) UPGMA 方法對(duì)不同施肥處理土壤細(xì)菌 16S rRNA DGGE 圖譜條帶進(jìn)行群落結(jié)構(gòu)相似性聚類分析。Lane1—Lane36每3條為同一處理土壤樣品，依次代表BS10～BS1、CK2、CK1處理土壤樣品。由圖3可知，不同處理土壤細(xì)菌群落結(jié)構(gòu)有明顯變化，CK1與CK2處理土壤細(xì)菌群落相似度大于70%;CK1與BS1、BS2沼液處理相似度高于67%，與BS3、BS4、BS5、BS6相似度僅37%，與BS7、BS8、BS9、BS10相似度僅34%。說(shuō)明經(jīng)長(zhǎng)期不同施肥處理，土壤細(xì)菌類群繁育會(huì)出現(xiàn)差異，以適應(yīng)不同的土壤環(huán)境。

隨沼液施用量增加，土壤細(xì)菌群落結(jié)構(gòu)逐步發(fā)生變化。根據(jù)相似度大小將BS處理樣品分組，表明沼液施用量相近的土壤細(xì)菌種群相似度較高，其中BS1與BS2（Lane7- 12）聚為一類，為Ⅰ組，BS3、BS4、BS5、BS6（Lane13-24）聚類為Ⅱ組，BS7、BS8、BS9、BS10（Lane25-36）聚類為Ⅲ組。不同組別之間相似度范圍在34%～37%，其中BS1、BS2、BS3、BS4、BS5、BS6（Ⅰ組、Ⅱ組）與BS7、BS8、BS9、BS10（Ⅲ組）之間相似性最小，僅34%。說(shuō)明不同沼液施用量對(duì)土壤細(xì)菌群落相似性影響很大。

3? ?小結(jié)與討論

通過(guò)試驗(yàn)得出，長(zhǎng)期施用化肥導(dǎo)致土壤中某些細(xì)菌類群消失，但遺傳相似性較高，說(shuō)明細(xì)菌群落穩(wěn)定性較好;長(zhǎng)期施用沼液土壤中不僅有某些細(xì)菌類群消失，同時(shí)伴有新的類群出現(xiàn)，遺傳相似性降低。適量施用沼液能夠提高土壤細(xì)菌豐度和Shannon-Wiener指數(shù)，過(guò)多則效果降低，在本研究條件下，連續(xù)3 a沼液施用總量為339.31～396.81 t/hm2較適宜。

通過(guò)DGGE圖譜分析得出，不同處理泳道中多數(shù)條帶相同，說(shuō)明這些細(xì)菌種群穩(wěn)定，不易受到土壤環(huán)境變化的影響，但同時(shí)也有部分條帶缺失或特異性條帶出現(xiàn)。Marschner P等[11 ]研究指出，不同肥料種類、用量及配施方式均可改變土壤養(yǎng)分含量及組成等，進(jìn)而對(duì)土壤微生物產(chǎn)生影響。通過(guò)DGGE圖譜的數(shù)字化分析看出，不同施肥處理耕層土壤細(xì)菌豐富度、Shannon-Wiener指數(shù)、Simpson指數(shù)差異極顯著，均勻度指數(shù)差異不顯著。本研究中，經(jīng)化肥長(zhǎng)期施用土壤中細(xì)菌豐富度、Shannon-Wiener指數(shù)極顯著小于清水處理，與李英秀等[12 ]、徐永剛? ? ?等[13 ]研究結(jié)論一致;沼液處理土壤細(xì)菌豐富度和多樣性指數(shù)多小于清水處理，與朱金山等[14 ]的研究結(jié)果相似。化肥和沼液的大量施用有效增加了土壤有機(jī)質(zhì)、堿解氮、有效磷、速效鉀等養(yǎng)分含量，刺激喜富營(yíng)養(yǎng)微生物生長(zhǎng)繁育，但土壤自然條件下該類微生物種類較少，從而導(dǎo)致群落多樣性降低;相反，貧營(yíng)養(yǎng)微生物和土著微生物種類繁多，在清水處理和低沼液處理的養(yǎng)分較低的土壤中迅速繁衍，使得土壤細(xì)菌種群多樣性提高。同時(shí)，土壤養(yǎng)分本底值含量也是影響細(xì)菌種群豐富度的重要因素，如在土壤有機(jī)質(zhì)本底值含量較高時(shí)施入有機(jī)肥，其土壤微生物細(xì)菌多樣性低于在低有機(jī)質(zhì)含量土壤中施入有機(jī)肥[15 - 16 ]，連續(xù)的施肥處理不斷向較高養(yǎng)分含量土壤中施入更多養(yǎng)分，會(huì)導(dǎo)致土壤細(xì)菌種群多樣性降低。

隨沼液施用量增加，土壤細(xì)菌多樣性差異極顯著。當(dāng)沼液施用量低于396.81 t/hm2時(shí)，土壤細(xì)菌群落豐度、多樣性指數(shù)極顯著高于清水處理和常規(guī)施肥處理，超過(guò)此沼液施用量則土壤細(xì)菌群落多樣性減小，甚至低于化肥處理，這與黃繼川等[17 ]研究結(jié)果相同。土壤微生物群落結(jié)構(gòu)與土壤pH、有機(jī)質(zhì)、土壤養(yǎng)分及土壤酶活性密切相關(guān)，這些因素的變化均能影響土壤微生物群落結(jié)構(gòu)，改變微生物群落組成和數(shù)量[18 - 19 ]。李忠佩等[20 ]通過(guò)擬合稻田土壤中細(xì)菌Shannon-Wiener指數(shù)與速效磷的一元二次方程發(fā)現(xiàn)，土壤速效磷對(duì)土壤細(xì)菌多樣性影響存在明顯轉(zhuǎn)折點(diǎn)，即在速效磷含量超過(guò)37.8～49.3 mg/kg時(shí)土壤細(xì)菌多樣性指數(shù)下降。

相應(yīng)地，細(xì)菌Simpson指數(shù)與豐富度和Shannon-Wiener指數(shù)呈反比，當(dāng)豐富度與多樣性指數(shù)增大時(shí)細(xì)菌種群功能多樣、信息復(fù)雜，優(yōu)勢(shì)類群不突出;當(dāng)豐富度和多樣性指數(shù)降低時(shí)一些細(xì)菌種群消失，使適應(yīng)環(huán)境的細(xì)菌種群有更充足的養(yǎng)分、生長(zhǎng)繁殖旺盛，成為優(yōu)勢(shì)種，優(yōu)勢(shì)度較高。

為進(jìn)一步了解各處理土壤細(xì)菌群落結(jié)構(gòu)的變化，本研究通過(guò)UPGMA方法對(duì)樣品DGGE圖譜進(jìn)行聚類分析。結(jié)果表明，長(zhǎng)期施用化肥土壤中細(xì)菌種群遺傳相似性高達(dá)67%，這與孟慶杰等[21 ]、Zhong W H等[22 ]的結(jié)果一致。隨沼液施用量增加，土壤細(xì)菌遺傳相似性呈階段性變化，聚為一類的組內(nèi)處理細(xì)菌群落相似度均高于50%，組間相似性逐級(jí)降低。不同處理土壤微生物群落結(jié)構(gòu)受土壤C/N比影響深刻，不同施肥方式下的不同微域生態(tài)環(huán)境也會(huì)導(dǎo)致微生物群落結(jié)構(gòu)變化。在研究中，隨沼液施用量增加，土壤養(yǎng)分含量、水熱狀況、微生物區(qū)系組成等存在差異，土壤中存活的優(yōu)勢(shì)細(xì)菌種群通過(guò)改變?nèi)郝浣Y(jié)構(gòu)、優(yōu)化群落功能來(lái)適應(yīng)新的環(huán)境條件，從而導(dǎo)致代謝功能變化、遺傳相似性降低。

參考文獻(xiàn)：

[1] WENTZEL S，SCHMIDT R，PIEPHO H P， et al.? Response of soil fertility indices to longterm application of biogas and raw slurry under organic farming[J].? Applied Soil Ecology，2015，96（4）：99-107.

[2] 王遠(yuǎn)遠(yuǎn)，劉榮厚.? 沼液綜合利用研究進(jìn)展[J].? 安徽農(nóng)業(yè)科學(xué)，2007，35（4）：1089-1091.

[3] 張亞莉，董仁杰，劉玉青.? 沼肥在農(nóng)業(yè)生產(chǎn)中的應(yīng)用[J].? 安徽農(nóng)業(yè)科學(xué)，2007，35（35）：11549-11550.

[4] 黃世文，廖西元.? 沼肥用于水稻的現(xiàn)狀及展望[J].? 中國(guó)沼氣，2005，23（2）：23-26.

[5] CABALLEROMELLADO J.? Agriculture microbiology and microbe interaction with plants[J].? Revista Latinoamericana De Microbiología， 2006，48（2）：154-161.

[6] 譚益民，何苑皞，郭文平.? 基于分子技術(shù)的土壤微生物多樣性研究進(jìn)展[J].? 中南林業(yè)科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2014（10）：1-9.

[7] 邢華銘，杜海濤，張黎黎，等.? PCR-DGGE與Biolog 技術(shù)在土壤微生物多樣性研究中的比較[J].? 農(nóng)業(yè)開(kāi)發(fā)與裝備，2013（10）：48;68.

[8] CAN S，LI PL，YAN Q D，et al.? Culture- independent methods for studying environmental microorganisms：methods，application， and perspective[J].? Applied Microbiology & Biotechnology，2011，93（3）：993-1003.

[9] GU Y F，ZHANG X P，TU S H，et al.? Soil microbial biomass，crop yields，and bacterial community structure as affected by long-term fertilizer treatments under wheat-rice cropping[J].? European Journal of Soil Biology，2009，45（3）：239-246.

[10] HEUER H，KRSEK M，BAKER P，et al. Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients[J].? Applied and environmental microbiology，1997，63（8）：3233-3241.

[11] MARSCHNER P，CROWLEY D，YANG C H.? Development of specific rhizosphere bacterial communities in relation to plant species， nutrition and soil type[J].? Plant and Soil，? 2004，261（1-2）：199-208.

[12] 李秀英，趙秉強(qiáng)，李絮花，等.? 不同施肥制度對(duì)土壤微生物的影響及其與土壤肥力的關(guān)系[J].? 中國(guó)農(nóng)業(yè)科學(xué)，2005，38（8）：1591-1599.

[13] 徐永剛，宇萬(wàn)太，馬? ?強(qiáng)，等.? 長(zhǎng)期不同施肥制度對(duì)潮棕壤微生物生物量碳、氮及細(xì)菌群落結(jié)構(gòu)的影響[J].? 應(yīng)用生態(tài)學(xué)報(bào)，2010，10（8）：2078-2085.

[14] 朱金山，張? ?慧，馬連杰，等.? 不同沼灌年限稻田土壤微生物群落分析[J].? 環(huán)境科學(xué)，2018，39（5）：2400-2411.

[15] KONG W D，LIU K X， LIAO Z W.? Effects of organic matters on metabolic functional diversity of soil microbial community under pot incubation conditions[J]. Acta Ecologica Sinica，2005（9）：2291-2296.

[16] BENDING G D，TURNER M K，JONES J E. Interactions between crop residue and soil organic matter quality and the functional diversity of soil microbial communities[J].? Soil Biology and Biochemistry，2002，34（8）：1073-1082.

[17] 黃繼川，徐培智，彭智平，等.? 基于稻田土壤肥力及生物學(xué)活性的沼液適宜用量研究[J].? 植物營(yíng)養(yǎng)與肥料學(xué)報(bào)，2016，22（2）：362-371.

[18] WALSH J J，ROUSK J，EDWARDS-JONES G，et al.? Fungal and bacterial growth following the application of slurry and anaerobie digestate of livestock manure to temperate pasture soils[J].? Biology and Fertility of Soils， 2012，48：889-897.

[19] MANGALASSERY S，MOONEY S J，SPA-RKES D L，et al.? Impacts of zero tillage on soil enzyme activities， microbial characteristics and organic matter functional chemistry in temperate soils[J].? European Journal of Soil Biology， 2015，68：9-17.

[20] 李忠佩，吳曉晨，陳碧云.? 不同利用方式下土壤有機(jī)碳轉(zhuǎn)化及微生物群落功能多樣性變化[J].? 中國(guó)農(nóng)業(yè)科學(xué)，2007，40（8）：1712-1721.

[21] 孟慶杰，許艷麗，李春杰，等.? 不同施肥、土地利用方式對(duì)黑土細(xì)菌多樣性的影響[J].? 大豆科學(xué)，2008（3）：480-486.

[22] ZHONG W H，CAI Z C，YIN L C，et al. Effects of the long-term application of inorganic fertilizers on microbial community diversity in rice-planting red soil as studied by using PCR-DGGE[J].? Acta Ecologica Sinica，2007，27（10）：4011-4018.

（本文責(zé)編：陳? ? 偉）