劉金雨，張瑞，羅亮，王世會(huì)，郭坤，白慶利，徐偉，趙志剛

（1.中國(guó)水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所，黑龍江省冷水性魚類種質(zhì)資源及增養(yǎng)殖重點(diǎn)開放實(shí)驗(yàn)室，黑龍江 哈爾濱 150070；2.上海海洋大學(xué)水產(chǎn)與生命學(xué)院，上海 201306）

稻漁綜合種養(yǎng)是將水稻種植和水產(chǎn)養(yǎng)殖相結(jié)合的復(fù)合農(nóng)業(yè)生產(chǎn)方式[1]，利用水稻與水生生物之間的互利共生關(guān)系，提升種養(yǎng)系統(tǒng)物質(zhì)循壞和能量流動(dòng)的效率[2]。目前，稻漁綜合種養(yǎng)已形成稻-魚、稻-蝦、稻-蟹、稻-鱉、稻-鰍等5 大類19 種典型模式[3]，具有穩(wěn)糧、促漁、提質(zhì)、增效和可持續(xù)發(fā)展的特點(diǎn)[4,5]。全國(guó)稻蟹綜合種養(yǎng)面積占稻漁綜合種養(yǎng)總面積的4.97%，主要分布在北方地區(qū)[6]。黑龍江省作為我國(guó)水稻種植面積最大的省區(qū)之一，種植面積達(dá)400 多萬(wàn)公頃，發(fā)展稻漁綜合種養(yǎng)潛力巨大，稻蟹綜合種養(yǎng)也是目前黑龍江省發(fā)展稻漁綜合種養(yǎng)的主要模式之一。

關(guān)于寒區(qū)稻蟹綜合種養(yǎng)目前多集中于通過(guò)機(jī)體生長(zhǎng)、水質(zhì)理化因子、浮游生物和底棲生物等參數(shù)來(lái)評(píng)價(jià)種植、施肥和投喂等技術(shù)層面的研究[7-11]，缺乏種養(yǎng)系統(tǒng)土壤微生物評(píng)價(jià)等的基礎(chǔ)研究[12,13]，尤其對(duì)寒區(qū)不同稻蟹共作模式下土壤微生物群落結(jié)構(gòu)和多樣性研究尚未見報(bào)道。稻田土壤系統(tǒng)中的微生物是土壤物質(zhì)循環(huán)的重要組成部分，在稻田系統(tǒng)中發(fā)揮著至關(guān)重要的作用[14,15]。水產(chǎn)動(dòng)物的攝食、排泄等活動(dòng)通過(guò)食物網(wǎng)影響稻漁綜合種養(yǎng)系統(tǒng)的物質(zhì)循環(huán)和能量流動(dòng)[16]。研究稻田土壤微生物群落組成和多樣性有助于明晰稻漁綜合種養(yǎng)的共作機(jī)制。細(xì)菌是微生物的主要組成部分，研究土壤細(xì)菌的群落結(jié)構(gòu)和多樣性可以在一定程度上掌握土壤微生物的結(jié)構(gòu)特征。土壤中大部分細(xì)菌因較難分離培養(yǎng)，傳統(tǒng)方法無(wú)法滿足土壤微生物群落結(jié)構(gòu)的研究需要[17,18]。Illumina MiSeq 測(cè)序平臺(tái)高效準(zhǔn)確的特點(diǎn)更適用于土壤微生物群落結(jié)構(gòu)和多樣性的研究[19,20]。本研究采用高通量測(cè)序技術(shù)比較分析了寒區(qū)稻-蟹共作（RC）、稻-蟹-鰍共作（RCL）和水稻單作（CK）三種模式下表層土壤細(xì)菌（16S rDNA）的多樣性，探討稻-蟹共作和稻-蟹-鰍共作模式對(duì)水稻表層土壤微生物群落結(jié)構(gòu)和多樣性的影響，為寒區(qū)稻蟹綜合種養(yǎng)模式的建立提供理論依據(jù)。

1 材料與方法

1.1 實(shí)驗(yàn)設(shè)計(jì)

實(shí)驗(yàn)在黑龍江省綏化市國(guó)家級(jí)稻漁綜合種養(yǎng)示范區(qū)進(jìn)行。實(shí)驗(yàn)田采用：稻田單作模式（CK）、稻蟹共作模式（RC）和稻-蟹-鰍共作模式（RCL），分別設(shè)3 個(gè)重復(fù)，共9 個(gè)單元，每個(gè)單元面積為300 m2，栽培管理方式一致。各單元內(nèi)，在距離田埂0.6 m 處挖口寬0.5 m，底寬0.3 m，深0.4 m 的環(huán)溝，溝田比為10%。在各單元的外圍利用塑料布設(shè)高于地面0.4 m 的防逃墻。共作模式的水稻均采用大壟雙行，邊行加密的栽培模式，水稻單作組仍采用等間距的常規(guī)稻田種植方式，水稻品種為“綏粳4 號(hào)”。稻蟹模式放養(yǎng)河蟹（Eriocheir sinensis）體質(zhì)量為（21.45±1.45）g，放養(yǎng)密度為500 ind./667 m2；稻-蟹-鰍模式，河蟹體質(zhì)量為（21.45±1.45）g/ind.，放養(yǎng)密度為350 ind./667 m2，體質(zhì)量為（0.47±0.07）g 泥鰍（Misgurnus anguillicaudatus）的放養(yǎng)密度為2 223 inds/667 m2。河蟹和泥鰍的放養(yǎng)時(shí)間分別為2020 年6 月25日和7 月5 日。養(yǎng)殖期間投喂配合飼料，每天投喂一次，投喂量為河蟹體體質(zhì)量的2%～3%。河蟹和泥鰍于2020 年9 月10 日收獲，水稻于9 月20 日收割。施有機(jī)肥為基肥，化肥為提苗肥，實(shí)驗(yàn)期間均不使用農(nóng)藥。

1.2 土壤樣品采集和處理

2020 年9 月20 日分別在各單元的稻田中隨機(jī)選5 個(gè)點(diǎn)采集表層土壤樣品，混合為1 個(gè)樣品。樣品編號(hào)分別為CK_1、CK_2、CK_3、RC_1、RC_2、RC_3、RCL_1、RCL_2 和RCL_3。鮮土去除雜質(zhì)混勻裝入50 mL 離心管后于-80 ℃冷凍保存，用于微生物多樣性分析。

1.3 土壤微生物DNA 提取與檢測(cè)

用試劑盒提取微生物基因組DNA，NanoDrop 2000 紫外-可見分光光度計(jì)檢測(cè)DNA 純度和濃度?；蚪MDNA 抽提后，用1%瓊脂糖凝膠電泳檢測(cè)DNA 完整性，PCR 預(yù)擴(kuò)增檢測(cè)樣品。

1.4 PCR 擴(kuò)增及高通量測(cè)序

采用通用引物338F（5'ACTCCTACGGGAGGGA GCA-3'）[21]和806R（5'-GGACTCHV GGGTWTCTAAT-3'）[22]擴(kuò)增細(xì)菌16S rRNA 基因V3-V4 區(qū)，PCR 采 用TransGen AP221-02：TransStart Fastpfu DNA Pliymerase，PCR 儀為ABI GeneAmp9700 型。PCR 反應(yīng)使用以下程序進(jìn)行：95 ℃預(yù)變性120 s，95℃變性30 s，55 ℃退火30 s，72 ℃延伸30 s，持續(xù)25個(gè)循環(huán)周期，最后72 ℃延伸5 min，每個(gè)樣品3 個(gè)重復(fù)。使用AxyPrepDNA 凝膠回收試劑盒切膠回收PCR 產(chǎn)物，Tris_HCL 洗脫；利用2%瓊脂糖凝膠電泳檢測(cè)PCR 擴(kuò)增效果，借助QuantiFluorTM-ST 藍(lán)色熒光定量系統(tǒng)對(duì)PCR 產(chǎn)物進(jìn)行檢測(cè)定量。高通量測(cè)序采用Illumina 公司的Miseq 平臺(tái)，土壤微生物總DNA 提取和測(cè)序服務(wù)委托上海美吉生物醫(yī)藥科技公司完成。

1.5 生物信息分析

利用QIIME1.9.1 對(duì)原始DNA 序列進(jìn)行過(guò)濾處理，去除嵌合體，嵌合體由uchime 確定，得到優(yōu)化序列[23]；利用UPARSE 7.1 軟件按照97%相似度將非重復(fù)序列進(jìn)行聚類，得到OTU（Operational Taxonomic Unit）的代表序列，在此基礎(chǔ)上進(jìn)行分類學(xué)分析，統(tǒng)計(jì)各個(gè)分類水平上樣品群落組成。

1.6 數(shù)據(jù)分析

實(shí)驗(yàn)結(jié)果采用平均值±標(biāo)準(zhǔn)誤差（Mean±SE）表示，數(shù)據(jù)進(jìn)行單因素方差（ANOVA）分析，并進(jìn)行Duncan 氏多重比較，以P<0.05 作為差異顯著性標(biāo)志。用Venn 圖統(tǒng)計(jì)多個(gè)樣本中所共有和獨(dú)有的物種數(shù)目。Alpha 多樣性分析包括Chao 指數(shù)、ACE 指數(shù)、Shannon 指數(shù)和Coverage 指數(shù)分析。Beta 多樣性分析通過(guò)未加權(quán)（Unweighted UniFrac）計(jì)算距離矩陣，根據(jù)得出的距離矩陣進(jìn)行主坐標(biāo)分析（PCoA）。通過(guò)土壤門和屬水平的細(xì)菌菌群結(jié)構(gòu)，分析各模式的微生物組成。

2 結(jié)果與分析

2.1 DNA 測(cè)序數(shù)據(jù)和OTUs

由表1 可知，RC、RCL 和CK 組的9 個(gè)樣品中土壤細(xì)菌V3-V4 區(qū)測(cè)序獲得的平均序列數(shù)分別為69 561.33、65 007.00 和65 872.00；平均堿基數(shù)分別為28 977 495.33、27 084 291.33 和27 493 837.00；序列平均長(zhǎng)度分別為416.57 bp、416.61 bp 和417.35 bp。以97%相似度劃分，共得到6 404 個(gè)OTUs，平均OTU 數(shù)分別為5 026、5 305 和4 969。用Coverage 反映測(cè)序深度，各組覆蓋率超過(guò)97.0%（表1）。用稀釋曲線來(lái)反映測(cè)序深度，稀釋曲線達(dá)到飽和平臺(tái)表明測(cè)序數(shù)據(jù)量合理，圖1 表示測(cè)序深度足以描述微生物多樣性。

表1 不同模式下細(xì)菌高通量測(cè)序結(jié)果統(tǒng)計(jì)Tab.1 Illumina MiSeq sequencing results of the bacterial community in different modes

圖1 三種模式稻田土壤微生物的稀釋曲線Fig.1 Rarefaction curve analysis of microbial diversity in the paddy soils in three modes

Venn 圖評(píng)估不同組OTU 的相似性和重疊情況見圖2。圖2 中數(shù)字表示所在區(qū)域的OTU 數(shù)目。結(jié)果顯示：三種模式共有3 679 個(gè)OTUs，占總OTU 的57.45%，說(shuō)明三種模式中的微生物種類較為相似。兩兩共享的OTU 中，RCL 和RC 共享數(shù)量最多，共有623 個(gè)，其次為RCL 和CK，共有550 個(gè)OTUs，最低的為CK 和RC，共有365 個(gè)OTUs。不同的OTU代表了不同的物種，特有的OTU 中，RCL 最高，為453 個(gè)，表示RCL 的特有物種最多；其次為CK，375個(gè)OTUs；RC 特有物種最少，為359 個(gè)OTUs。

圖2 三種模式稻田土壤微生物的韋恩圖Fig.2 Venn diagram of microbial diversity in the paddy soils in three modes

2.2 多樣性指數(shù)

α 多樣性是指群落內(nèi)部種的多樣性，常用指標(biāo)包 括ACE 指數(shù)、Chao 指數(shù)、Shannon 指數(shù)以及Simpson 指數(shù)等。ACE 指數(shù)用來(lái)估算群落物種豐富度，Chao 指數(shù)用來(lái)估計(jì)群落中含OTU 數(shù)目的指數(shù)，ACE 指數(shù)和Chao 指數(shù)越大說(shuō)明細(xì)菌群落豐度越高。Shannon 指數(shù)和Simpson 指數(shù)用來(lái)反映群落物種多樣性，Shannon 指數(shù)值越大或Simpson 指數(shù)越小，說(shuō)明細(xì)菌群落多樣性越高。與CK 相比，RC 和RCL 模式土壤細(xì)菌的ACE 指數(shù)增加了0.13%和1.58%，Chao 指數(shù)增加了0.88%和2.79%，Shannon指數(shù)增加了1.54%和2.71%（表2）。方差分析結(jié)果表明，2 種共作模式的稻田土壤細(xì)菌群落豐度和多樣性均高于水稻單作，但差異并不顯著（P>0.05）。

表2 不同模式下細(xì)菌群落多樣性指數(shù)Tab.2 Bacterial diversity index of microbial diversity in the paddy soils in three modes

β 多樣性表示的是微生物群落之間的差異。對(duì)三種模式下稻田土壤的OTU 進(jìn)行PCoA 分析，第一主軸和第二主軸對(duì)細(xì)菌結(jié)構(gòu)變異的貢獻(xiàn)率分別為56.97%和12.60%（圖3）。圖中每個(gè)點(diǎn)表示一個(gè)樣品，點(diǎn)之間的距離表示樣本間微生物群落結(jié)構(gòu)的相似程度，距離越近越相似。PCoA 圖顯示RC 和RCL距離較近，同時(shí)二者與CK 的距離較遠(yuǎn)，且出現(xiàn)顯著性差異（R=0.4239，P=0.045）。

圖3 不同模式微生物群落的主坐標(biāo)分析（PCoA）Fig.3 Principal coordinate analysis（PCoA）of microbial diversity in the paddy soils in three modes

2.3 細(xì)菌門水平分類

本次測(cè)序中，得到的所有序列均為細(xì)菌序列。9個(gè)稻田土壤樣品的6 404 個(gè)OTUs 分別屬于59 個(gè)門、181 個(gè)綱、397 個(gè)目、1 093 個(gè)屬、2 245 個(gè)種。圖4為門水平上的細(xì)菌分類，平均豐度低于1%的物種被合并為圖中的others。相對(duì)豐度大于5%的優(yōu)勢(shì)菌門均為相同的7 個(gè)菌門，分別為放線菌門（Actinobacteria）、綠彎菌門（Chloroflexi）、變形菌門（Proteobacteria）、擬桿菌門（Bacteroidota）、酸桿菌門（Acidobacteria）、厚壁菌門（Firmicutes）以及脫硫桿菌門（Desulfobacterota）。上述菌門在三種模式中豐度不同，但均未出現(xiàn)顯著性差異（P>0.05）。CK 中相對(duì)豐度從高到低的菌門依次為：放線菌門（18.68%）、變形菌門（15.39%）、綠彎菌門（14.20%）、擬桿菌門（13.34%）、厚壁菌門（10.81%）、酸桿菌門（8.05%）以及脫硫桿菌門（7.68%）。RCL 相對(duì)豐度從高到低的菌門依次為：綠彎菌門（19.08%）、放線菌門（16.07%）、變形菌門（12.44%）、酸桿菌門（11.46%）、擬桿菌門（10.82%）、厚壁菌門（9.89%）和脫硫桿菌門（6.36%）。RC 相對(duì)豐度從高到低的菌門依次為：放線菌門（20.31%）、綠彎菌門（17.76%）、酸桿菌門（12.54%）、變形菌門（11.80%）、擬桿菌門（10.42%）、厚壁菌門（8.92%）和脫硫桿菌門（5.34%）。

圖4 三種模式稻田土壤微生物在門水平的細(xì)菌物種組成Fig.4 Bacterial community composition in the paddy soils in three modes at the phylum level

2.4 細(xì)菌屬水平分類

如圖5 所示，熱圖顯示了三種模式在屬水平上微生物群落的相對(duì)豐度的差異。平均豐度低于1%的物種被合并為others，其在CK、RC 和RCL 中分別占59.17%、55.74%和57.62%。相對(duì)豐度>1%的優(yōu)勢(shì)菌屬有22 個(gè)，其中一些屬?zèng)]有科學(xué)名稱，以norank 作為標(biāo)記。CK 中相對(duì)豐度前五的優(yōu)勢(shì)菌屬為稻田土壤菌屬（Oryzihumus，4.32%）、擬桿菌屬（norank_f_Bacteroidetes_vadinHA17，3.71%）、假節(jié)桿屬（Pseudarthrobacter，3.85%）、（norank_f_norank_o_SJA-15，3.57 %）和（norank_f_Lentimicrobiaceae,2.84%）。RC 中相對(duì)豐度前五的菌屬為稻田土壤菌屬（4.86%）、假節(jié)桿屬（4.08%）、（norank_f_ norank_o_norank_c_KD4-96，4.01%）、擬桿菌屬（3.86%）以及（norank_f_norank_o_RBG-13-54-9，2.44%）。RCL中相對(duì)豐度前五的菌屬分別為擬桿菌屬（4.02%）、（norank_f_norank_o_norank_c_KD4-96，3.97%）、（no rank_f_norank_o_RBG-13-54-9，2.96%）、稻田土壤菌屬（2.87%）以及假節(jié)桿屬（2.81%）。

圖5 三種模式稻田土壤微生物在屬水平的豐度熱圖Fig.5 Heatmap of microbial diversity in the paddy soils in three modes at the phylum level

3 討論

3.1 不同稻蟹共作模式土壤的細(xì)菌多樣性

微生物多樣性指數(shù)是評(píng)價(jià)微生物群落豐富度和多樣性的綜合指標(biāo)，較高的生物多樣性更能保證生態(tài)系統(tǒng)的穩(wěn)定[24]。微生物的多樣性變化可直接或間接影響土壤肥力[25]。在稻漁綜合種養(yǎng)系統(tǒng)中，灌溉、施肥、栽培和共作方式均會(huì)影響土壤微生物[26,27]，稻漁共作過(guò)程中養(yǎng)殖動(dòng)物和飼料投喂水平也會(huì)造成土壤微生物多樣性的變化[28,29]。研究表明，水產(chǎn)動(dòng)物與水稻共作可提升稻田土壤微生物的群落結(jié)構(gòu)和多樣性[30-33]。本研究結(jié)果顯示，在寒區(qū)進(jìn)行RC 和RCL 后的稻田土壤微生物多樣性指數(shù)均高于水稻單作，與上述結(jié)果一致，表明這兩種寒區(qū)稻蟹共作模式均提升了稻田土壤細(xì)菌群落的多樣性。同時(shí)，RCL 的Shannon 指數(shù)、ACE 指數(shù)和Chao 指數(shù)均高于RC，表明稻-蟹-鰍共作模式較稻-蟹模式更能提升稻田土壤微生物多樣性，在一定程度上揭示了稻-蟹-鰍共作模式下，土壤生態(tài)系統(tǒng)相對(duì)更為穩(wěn)定。宋宇等[33]研究結(jié)果中的PCoA 圖顯示，稻蝦共作和稻蟹共作均與水稻單作的距離較遠(yuǎn)，而兩種共作模式之間的距離較近，這與本實(shí)驗(yàn)PCoA 分析圖中RC 和RCL 與CK 在不同區(qū)域聚類，RC 和RCL距離較近的結(jié)果相似。本實(shí)驗(yàn)土壤微生物不同區(qū)域聚類表明稻蟹共作模式與水稻單作模式的稻田土壤微生物群落結(jié)構(gòu)差異較大，而兩種稻蟹共作模式稻田土壤微生物的群落結(jié)構(gòu)較為相近。稻蟹共作模式在種養(yǎng)過(guò)程中，其殘餌、糞便可為土壤微生物提供營(yíng)養(yǎng)；河蟹和泥鰍的活動(dòng)也對(duì)土壤產(chǎn)生擾動(dòng)，增加土壤的通透性，這些因素均可能顯著影響稻蟹共作系統(tǒng)中稻田土壤微生物的組成和多樣性。

3.2 不同稻蟹共作模式土壤的細(xì)菌群落結(jié)構(gòu)組成

高通量測(cè)序技術(shù)可檢測(cè)到更多豐度較低的微生物類群[34]。本研究中，3 種模式在97%的相似性下共獲得6 404 個(gè)OTUs，檢測(cè)到59 個(gè)門1 093 個(gè)屬的微生物種群。在門的分類水平上，三種模式的優(yōu)勢(shì)菌門均為放線菌門、綠彎菌門、變形菌門、擬桿菌門、酸桿菌門、厚壁菌門及脫硫桿菌門。放線菌可加速土壤中動(dòng)植物殘?bào)w的腐爛分解[35]。宋宇等[36]研究認(rèn)為，稻蟹共作土壤放線菌的豐度較水稻單作增加了1.36%。本研究中，稻蟹模式土壤放線菌的豐度增加了1.63%，土壤放線菌豐富度的增加有助于加速稻田土壤系統(tǒng)的物質(zhì)循環(huán)進(jìn)程。變形菌門大多為革蘭氏陰性菌，其中大部分細(xì)菌均具有固氮能力[37]。成永旭等[38]研究表明，稻蟹共作模式土壤細(xì)菌相對(duì)豐度最高的菌門為變形菌門。羅衡等[39]對(duì)稻鱉共作模式土壤微生物的研究也認(rèn)為，土壤微生物以變形菌門為主導(dǎo)。本研究在寒區(qū)稻蟹共作模式中也得到了相似的結(jié)果。稻魚共作土壤系統(tǒng)中綠彎菌門豐度大幅上升，成為優(yōu)勢(shì)菌門[30]。本實(shí)驗(yàn)稻蟹共作模式也有效提升了土壤綠彎菌門的相對(duì)豐度，可能因共作模式稻田水淹程度較大，土壤中形成的微厭氧環(huán)境適宜該菌的生長(zhǎng)。擬桿菌門、酸桿菌門、厚壁菌門及脫硫桿菌門均為稻田常見的細(xì)菌類群[40]。本研究在寒區(qū)兩種稻蟹共作模式和水稻單作模式中也分別監(jiān)測(cè)到了這些菌群。稻蟹共作模式中，宋宇等[33]檢測(cè)的土壤細(xì)菌優(yōu)勢(shì)菌門與本實(shí)驗(yàn)相似，而在屬分類水平上差異較大，前者優(yōu)勢(shì)屬菌為藍(lán)藻屬（norank_c_Cyarbacteria）、硫曲菌屬（Sulfuricurvum）和厭氧繩菌屬（norank_f_Anaerolineareae），而本研究中的優(yōu)勢(shì)菌屬分別為稻田土壤菌屬、假節(jié)桿屬、norank_f_norank_o_norank_c_KD4-96、擬桿菌屬和norank_f_norank_o_RBG-13-54-9，主要原因可能是不同實(shí)驗(yàn)地點(diǎn)的氣候環(huán)境差異所致。本研究在土壤主要優(yōu)勢(shì)菌屬中，對(duì)比水稻單作，稻-蟹共作和稻-蟹-鰍共作土壤中的norank_f_norank_o_no－rank_c_KD4-96、擬桿菌屬和norank_f_norank_o_RBG-13-54-9 的相對(duì)豐度均有所增加，可能因兩種稻蟹共作模式中殘餌、糞便提升了稻田土壤系統(tǒng)的肥力，在一定程度上為上述細(xì)菌類群提供了更為適宜的生長(zhǎng)環(huán)境。

3.3 結(jié)論

本研究中，寒區(qū)“稻-蟹”和“稻-蟹-鰍”兩種共作模式下稻田土壤細(xì)菌群落結(jié)構(gòu)和多樣性結(jié)果表明，與水稻單作相比，稻蟹共作模式和稻-蟹-鰍共作模式均能改變稻田土壤的細(xì)菌群落結(jié)構(gòu)，增加稻田土壤細(xì)菌群落的豐富度和多樣性。