不同稻田栽培模式下土壤細菌群落多樣性分析

宋 宇,王 鵬,韋月平

(遼東學院, 遼寧 丹東 118003)

【研究意義】土壤微生物數量巨大，種類繁多，細菌是土壤中數量最為豐富、分布最為廣泛的微生物類群[1]。土壤細菌在有機物質分解、養(yǎng)分運轉、土壤結構維持和溫室氣體產生等方面發(fā)揮著重要的調節(jié)作用[2-3]。同時，土壤細菌的數量和結構受土壤類型、土壤理化性質、植被類型和農田管理制度等多方面因素的影響[4-6]。由于世界人口的增長，糧食安全正成為一個全球性的問題。大米是大約30億人日常飲食中的主要成分。因此水稻生產是全球糧食安全的重要組成部分。過去幾十年中，由于大量使用化肥和農藥，水稻產量有所增加，但化肥和農藥的過量使用會導致土壤酸化板結、溫室氣體排放增加、水體富氧化等一系列生態(tài)環(huán)境問題[7-9]。如今，以當地和可再生資源為基礎，將個別農場的作物和動物生產有效結合起來以提高養(yǎng)分循環(huán)效率為目的的生態(tài)循環(huán)農業(yè)則受到人們廣泛關注。相關研究描述了水稻-魚、稻-鴨、稻-龜和稻蛙種植生態(tài)系統并說明了這些做法如何降低化肥和農藥施用量[10]，但對水稻種植模式的研究大多集中在病蟲害、土壤理化性質的影響上[11-13]，而對不同栽培模式下稻田土壤細菌多樣性的研究則鮮見。【前人研究進展】傳統的土壤微生物分離培養(yǎng)方法難以反映土壤微生物群落的結構、豐度及其作用機制[14],20世紀末，隨著PCR、DGGE、T-RFLP等分子生物技術的發(fā)展為研究環(huán)境樣品中微生物的群落多樣性和豐度成為可能，但檢測到的細菌群落豐度相對較低，范圍較小且操作復雜成本高[15]。高通量測序技術具有靈敏度高的優(yōu)點且低錯誤率、低成本可較快速準確地獲取土壤微生物信息[16]因而被應用于細菌群落多樣性分析中。【本研究切入點】本研究使用高通量測序技術開展3種不同栽培模式下稻田表層土壤細菌多樣性研究?！緮M解決的關鍵問題】探討不同栽培模式稻田表層土壤細菌多樣性，不同分類水平的優(yōu)勢類群及各類群的相對豐度，細菌多樣性和土壤理化性質的關系。為生態(tài)循環(huán)農業(yè)提升稻田土壤質量合理選擇水稻栽培模式提供理論基礎。

1 材料與方法

1.1 試驗地概況與試驗設計

試驗于2018年9月28日水稻結穗期在丹東東港市北井子鎮(zhèn)蓋家壩農場進行,供試稻田分別是河蟹共作稻田(SDHXZ)、非洲雁共作稻田(SDFYZ)和單作稻田(SDKBZ)，分別在3個模式稻田進行5點取樣，采集0～5 cm的表層土壤充分混勻，每份土樣分成2份，1份去用于細菌多樣性分析，1份用于測土壤理化性質。

1.2 測定方法

1.2.1 土壤理化性質的測定 土壤pH值采用酸度計測定，氨氮和硝氮測定采用氯化鉀浸提-分光光度法，速效鉀測定采用氯化鈣浸提-火焰光度計法,速效磷測定采用碳酸氫鈉浸提-鉬銻抗比色法，有機質測定采用重鉻酸鉀氧化-分光光度法。

1.2.2 土壤細菌總DNA提取及高通量測序 樣品土壤細菌總DNA提取及高通量測序均在上海美吉生物醫(yī)藥科技有限公司進行。細菌多樣性鑒定區(qū)域：16SV3～V4區(qū)。通用引物為338F(5’-ACTCCTACGGGAGGGAGCA-3’)和806R(5’-GGACTCHVGGGTWTCTAAT-3’)。

1.3 數據處理與生物信息統計方法

使用Trimmomatic和Flash軟件進行原始雙端序列去雜和拼接。采用Usearch軟件進行OTU分析并采用RDP classifier 貝葉斯算法對97 %相似水平的OTU代表序列進行分類學分析。利用mothur軟件計算Chao1指數、Ace指數、shannon指數等。利用Qiime計算beta多樣性距離矩陣并利用R語言作圖。通過statistic10.0軟件進行相關性分析。

2 結果與分析

2.1 不同栽培模式稻田土壤的理化性質

由表1可以看出,3種栽培模式下稻田土壤的理化性質差異較大，與水稻單作組比，河蟹共作稻田pH值下降而非洲雁共作稻田則上升。河蟹稻田和非洲雁稻田的氨氮，硝氮、速效磷和有機質含量都顯著提高，河蟹稻田的速效鉀值較單作稻田降低但非洲雁稻田則顯著提升。

2.2 土壤樣品測序

2.2.1 OTU及分類學分析 通過Illumina MiSeq 高通量測序技術測定9個土壤樣本細菌的16SrRNA V3～V4變異區(qū)序列，共獲得240 070個有效序列數，平均測序長度為420～460 bp。通過Usearch軟件結合RDP classifier 貝葉斯算法對97 %相似水平的OTU代表序列進行分類學分析后共獲得3157個OTU,分屬于48個門，115個綱，224個目，391個科，721個屬，1343個種。

表1 不同栽培模式的稻田土壤理化性質

注:數據由平均值和標準差組成,不同的字母代表處理間存在顯著差異(P<0.05)相同的字母代表處理間不存在顯著差異(P>0.05)。

Notes:The data consists of averages and standard deviations, different letters mean that there are significant differences(P<0.05),while the same letters mean no significant differences(P>0.05).

2.2.2 物種Venn圖分析 Venn圖是利用R語言對相似水平為97 %的OTU進行統計和作圖，可用于多個樣本中所共有和獨有的物種(如OTU)數目，可以比較直觀的表現環(huán)境樣本的物種(如OTU)組成相似性及重疊情況。從圖1可以看出，3種栽培模式稻田土壤細菌共有的OTU數為1525個，SDKBZ 和SDHXZ共有的OTU數最多為2183個，SDKBZ和SDFYZ共有的OTU數最少為1610個。不同的OTU代表不同的物種，SDFYZ特有的OTU數最多為357個，說明特有的物種數最多。

2.3 土壤細菌群落多樣性分析

2.3.1 α-多樣性分析 α多樣性是指某個群落或生境內部的種的多樣性，常用的度量標準有Chao1、Shannon、Ace、Coverage等。Chao1指數用來估計群落中含OTU數目的指數。Shannon指數可用來反映樣品的多樣性程度，其值越高說明群落物種的多樣性越高。Ace指數用來估算群落物種的豐富度，其值越高表明物種群落豐度越高。Coverage指數可用來反映測序的深度。從表2可以看出，3種栽培模式的測序覆蓋度都達到了99 %，說明各樣品的取樣都能達到一定深度，能反映各樣品土壤細菌群落的真實情況。由表2可以看出，SDHXZ的Chao1、Shannon和Ace 3種多樣性指數均高于SDKBZ和SDFYZ，說明有機稻田在引入河蟹后土壤細菌的多樣性和豐富度均有提升。而SDFYZ的3種多樣性指數卻低于SDKBZ和SDHXZ說明非洲雁糞肥的稻田土壤，雖然土壤肥力提升顯著但其細菌的多樣性和豐度并未得到提升。

圖1 不同栽培模式稻田土壤細菌Venn圖分析

2.3.2 土壤細菌多樣性與理化因子的相關性分析 從表3可以看出，Ace指數和pH及氨態(tài)氮呈顯著負相關。Shannon指數和硝態(tài)氮呈正相關。Chao1和Ace指數均和6種理化因子呈負相關。

2.3.3 β多樣性分析 β多樣性表示的是微生物群落構成的比較, 可評估微生物群落間的差異。從圖2可以看出，SDFYZ和SDHXZ之間及二者與SDKBZ之間的距離都較遠，表明3種栽培模式稻田土壤細菌群落的OTU組成差異較大，細菌群落結構的差異也較大。

表2 不同栽培模式下稻田土壤細菌多樣性指數

表3 稻田土壤理化因子與細菌群落多樣性指數相關性分析

圖2 3種栽培模式稻田土壤細菌群落主成分(PCoA)分析圖

2.4 3種栽培模式稻田土壤細菌群落組成分析

2.4.1 細菌優(yōu)勢菌門分析 從圖3可知，3個樣本中含量最高的門是變形菌門(Proteobacteria)，平均豐度為35.51 %，其SDKBZ、SDHXZ和SEDFYZ的相對豐度依次為33.41 %±3.50 %、36.73 %±3.48 %、36.37 %±0.52 %；其次為擬桿菌門(Bacteroidetes)平均豐度為13.69 %，相對豐度依次為12.12 %±2.04 %、15.62 %±3.24 %、13.32 %±5.32 %；綠彎菌門(Chloroflexi)，平均豐度為13.66 %，相對豐度依次為12.22 %±2.32 %、13.64 %±2.48 %、15.12 %±2.63 %；放線菌門(Actinobacteria)，平均豐度為9.67 %，相對依次為8.23 %±1.32 %、6.87 %±1.16 %、13.92±7.28 %；酸桿菌門(Acidobacteria),平均豐度為7.90 %，相對豐度依次為10.31 %±1.05 %、6.88 %±2.37 %、6.52 %±1.34 %；藍藻菌門(Cyanobacteria)平均豐度為5.38 %，相對豐度依次為8.48 %±1.77 %、5.15 %±2.90 %、2.51 %±1.24 %。其余的硝化螺旋菌門(Nitrospirae)、厚壁菌門(Firmicutes)、疣微菌門(Verrucomicrobia)等相對豐度較低。說明在3種栽培模式下稻田土壤細菌群落的優(yōu)勢菌門為變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、綠彎菌門(Chloroflexi)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria)和藍藻菌門(Cyanobacteria)。

圖3 3種栽培模式稻田表層土壤優(yōu)勢菌門的微生物組成豐度圖

2.4.2 細菌優(yōu)勢菌屬分析 從圖4可知，平均豐度低于0.01 %的物種被合并為圖中的others，其在SDKBZ、SDHXZ和SDFYZ中分別占46.28 %、45.03 %和44.33 %。相對豐度>1 %的優(yōu)勢菌屬有18個。共有22個屬因在分類學數據庫分類學譜系的中間等級，沒有科學名稱，所以用norank作為標記。主要優(yōu)勢菌屬為norank-厭氧繩菌屬其在SDKBZ、SDHXZ和SDFYZ中的平均豐度依次為4.49 %、4.24 %和7.00 %；norank-藍細菌屬其平均豐度依次為6.31 %、3.24 %和1.38 %；norank-擬桿菌屬-VadinHA17其平均豐度依次為1.85 %、4.74 %和2.63 %；norank-酸桿菌屬其平均豐度依次為3.23 %、1.93 %和3.85 %；norank-p-Saccharibacteria其平均豐度依次為2.68 %、3.10 %和1.83 %。

3 討 論

本研究利用Illumina Miseq平臺對16srDNA基因的V3～V4變異區(qū)進行高通量測序分析，成功獲得3種栽培模式水稻表層土壤細菌群落結構多樣性的全面菌群信息。共獲得有效序列數240070個，3157個OTUs.其中SDKBZ的OTUs數為2443個、SDHXZ為2540個、SDFYZ為2168個。Venn 圖表明3種栽培模式共有的細菌 OTU 數為1525個，SDFYZ特有的 OTU 數為357個，高于SDHXZ 56.3 %,SDKBZ 50.98 %。雖然SDFYZ的 OTU 數低于SDHXZ和SDKBZ，但特有的OTU 數卻明顯高于SDHXZ和SDKBZ說明僅依靠非洲雁糞肥的稻田土壤可獲得更多的細菌特有物種。

圖4 3種栽培模式下水稻表層土壤優(yōu)勢屬的微生物組成豐度圖

本研究通過土壤理化因子測定和α-多樣性分析表明，由于河蟹的排泄作用，改變了稻田的pH值并增加了氨態(tài)氮、硝態(tài)氮及速效磷和有機質的含量，改善了稻田土壤細菌的生長條件同時河蟹通過爬行、打洞等活動有效地增強了土壤的透氣性，因此增加了河蟹稻田土壤細菌的多樣性和豐度。這和趙翔剛等[17-18]的研究結果一致。而SDFYZ通過非洲雁的排泄作用在沒有施加任何肥料的情況下雖然顯著提升了土壤各種養(yǎng)分含量但土壤細菌的多樣性和豐度并未得到提升，可能和非洲雁是站立行走而非爬行同時雁蹼平鈍而非銳利因此會增加土壤的緊密性而影響了土壤的透氣性有關。通過β多樣性分析表明3種稻田栽培模式的土壤細菌群落結構彼此差異都較大。

有研究表明土壤理化因子對土壤細菌群落有顯著的影響[19-20]，本文進行了土壤細菌多樣性與理化因子的相關性分析，Ace指數和pH和氨態(tài)氮呈顯著負相關。Shannon指數和硝態(tài)氮呈正相關。Chao1和Ace指數均和6種理化因子呈負相關。

通過樣品細菌群落組成分析表明3種栽培模式的土壤表層土壤細菌優(yōu)勢菌門為變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、綠彎菌門(Chloroflexi)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria)和藍藻菌門(Cyanobacteria)。其中以變形菌門(Proteobacteria)為最大優(yōu)勢菌門，這和以前的一些研究結果一致[21-23]。相對豐度>1 %的優(yōu)勢菌屬有18個。

4 結 論

稻田中放入一定量的河蟹既增加了土壤的養(yǎng)分含量又提升了土壤細菌的多樣性，因此具有更好的微生物生態(tài)基礎。稻田中放入非洲雁，可增加土壤肥力但并未提升土壤細菌的多樣性，若采取一定措施增加土壤透氣性這能否進一步提升土壤細菌的多樣性和豐度還有待于進一步探討。該研究可為合理選擇和改進水稻栽培模式提供一定理論支持。