劉蘭英，黃 薇，呂 新，何肖云，陳麗華，宋永康

（福建省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所/福建省農(nóng)產(chǎn)品質(zhì)量安全重點(diǎn)實(shí)驗(yàn)室，福建 福州 350003）

0 引言

【研究意義】沼液中含有豐富的營(yíng)養(yǎng)元素、氨基酸、B 族維生素和植物生長(zhǎng)素等生理活性物質(zhì)[1]，將沼液作為有機(jī)肥料施用于農(nóng)業(yè)土壤，可以為作物生長(zhǎng)提供充足的養(yǎng)分需求，同時(shí)也是處理畜禽養(yǎng)殖廢棄物的主要途徑之一。但畜禽養(yǎng)殖過程投入的大量抗生素只有少部分被動(dòng)物吸收，大部分最終會(huì)殘留在沼液中，而隨著沼液的還田利用，部分抗生素也隨之進(jìn)入土壤環(huán)境中[2]。土壤中抗生素的積累可能會(huì)誘導(dǎo)微生物產(chǎn)生抗生素抗性基因（Antibiotic resistance genes,ARGs），造成ARGs 的傳播擴(kuò)散[3]。研究表明，ARGs 能在環(huán)境中持久性殘留，其污染風(fēng)險(xiǎn)可能要比抗生素本身的風(fēng)險(xiǎn)更大[4]。因此，研究沼液施用到土壤后對(duì)土壤ARGs 的影響，對(duì)于進(jìn)一步評(píng)估沼液還田的安全性具有很重要的意義?！厩叭搜芯窟M(jìn)展】目前，關(guān)于糞肥施用對(duì)農(nóng)田土壤中ARGs 影響的報(bào)道較多。Peng 等[5]發(fā)現(xiàn)長(zhǎng)期的豬糞施用顯著增加了土壤中ARGs 的多樣性，豬糞處理的土壤比施用牛糞土壤含有更豐富的ARGs。Li 等[6]發(fā)現(xiàn)長(zhǎng)期的牛糞和雞糞施用均顯著增加了溫室土壤中抗生素和ARGs 的殘留水平，抗生素濃度與ARGs 豐度顯著正相關(guān)。Zhao 等[7]發(fā)現(xiàn)長(zhǎng)期的雞糞施用顯著提高了農(nóng)田土壤中ARGs 的豐度水平。McKinney 等[8]發(fā)現(xiàn)牛糞的施用顯著增加了土壤中ARGs 的豐度，且ARGs 豐度會(huì)隨土壤施肥量增加而增加?！颈狙芯壳腥朦c(diǎn)】然而，在眾多研究中，對(duì)于沼液施用對(duì)土壤中ARGs的影響及其季節(jié)變化方面的研究較少。唐鵬[2]報(bào)道了沼液施用對(duì)橡膠人工林土壤抗生素殘留的影響，但并未涉及ARGs 方面的研究。【擬解決的關(guān)鍵問題】本研究選取養(yǎng)豬場(chǎng)周邊長(zhǎng)期施用沼液的土壤為研究對(duì)象，同時(shí)采集當(dāng)?shù)匚词┯眠^沼液的土壤作為對(duì)照，采用熒光定量PCR 方法，對(duì)土壤中常見的14 種四環(huán)素類ARGs（tetA、tetB、tetC、tetG、tetH、tetK、tetM、tetO、tetQ、tetS、tetT、tetW、tetY、tetZ）和3 種磺胺類ARGs（sul1、sul2、sul3）進(jìn)行檢測(cè)，分析長(zhǎng)期施用沼液對(duì)土壤中ARGs 多樣性和豐度的影響，并探討其季節(jié)分布特征，以期為沼液的安全農(nóng)用提供理論依據(jù)。

1 材料與方法

1.1 采樣點(diǎn)基本情況

土壤樣品采集于建甌市某養(yǎng)豬場(chǎng)糞污處理沼氣工程周邊的農(nóng)田（27°13.30′N，118°12.80′E），已持續(xù)施用該養(yǎng)豬場(chǎng)沼液6 年。該養(yǎng)豬場(chǎng)目前生豬存欄數(shù)2 500 頭，生豬生產(chǎn)區(qū)及污水處理區(qū)占地面積約6 hm2，豬場(chǎng)污水日排放量約30 m3·d-1。施肥用沼液是由沼氣工程充分發(fā)酵后經(jīng)氧化塘處理的液體（沼液與氧化塘的水配比為1∶2）。供試沼液基本理化性質(zhì)為：pH 7.99，化學(xué)需氧量（COD）285 mg·L-1，五日生化需氧量（BOD5）263 mg·L-1，氨氮475 mg·L-1，總氮978 mg·L-1，總磷22.50 mg·L-1。采樣點(diǎn)土壤質(zhì)地為砂壤土，主要種植作物為檳榔芋，每年只種1 季，3 月中下旬種植，11 月初收獲，冬季閑田。沼液灌溉的田塊按照基肥和追肥方式灌入，在檳榔芋種植前，每667 m2灌入沼液量約15 t，與田土混勻作基肥；隨后每隔15 d 進(jìn)行1 次追肥，每次追肥灌入沼液量約6 t，大概8 月下旬進(jìn)行最后1 次追肥，整個(gè)生長(zhǎng)期都不施用其他肥料。其他農(nóng)藝管理措施均按照當(dāng)?shù)爻Ｒ?guī)習(xí)慣統(tǒng)一實(shí)施。

1.2 樣品采集與制備

分別于2019 年5 月（春季）、7 月（夏季）、10月（秋季）和2020 年1 月（冬季），在檳榔芋種植地塊上進(jìn)行根際土壤樣品的采集，其中春季和夏季的取樣時(shí)間均為下一次沼液施入的前兩天。采樣時(shí)除去表層枯枝落葉等雜物，用取土器沿著植株的根部采集0～20 cm 土層的土壤樣品，標(biāo)記為BS+，同時(shí)采集周邊一處地勢(shì)較高未施用過沼液的檳榔芋土壤作為對(duì)照樣品，標(biāo)記為BS-，這兩塊檳榔芋田同屬于一個(gè)農(nóng)戶，種植方式和田間管理方式相同。每塊樣地設(shè)置3 份重復(fù)，每份重復(fù)樣品均由五點(diǎn)取樣法獲得5 個(gè)小樣組成，現(xiàn)場(chǎng)混合均勻后裝入無菌塑封袋中，放入帶有冰塊的保溫箱中運(yùn)回實(shí)驗(yàn)室，-80 ℃冷凍保存待測(cè)。

1.3 土壤基因組DNA 提取

采 用TIANamp Soil DNA Kit（TIANGEN，北京）試劑盒，參照說明書，對(duì)土壤微生物總DNA 進(jìn)行提取，將提取的DNA 用1% 瓊脂糖凝膠進(jìn)行電泳，使用超微量分光光度計(jì)（Quawell Q3000，美國(guó)）檢測(cè)其濃度和純度。

1.4 土壤ARGs 測(cè)定

采用StepOnePlus? Real-Time PCR System（ABI，美國(guó)）實(shí)時(shí)熒光定量反應(yīng)平臺(tái)，對(duì)土壤中14 種四環(huán)素類ARGs（tetA、tetB、tetC、tetG、tetH、tetK、tetM、tetO、tetQ、tetS、tetT、tetW、tetY、tetZ）、3 種磺胺類ARGs（sul1、sul2、sul3）及16S rRNA 基因進(jìn)行qPCR 檢測(cè)和定量分析，其中四環(huán)素類ARGs 和16S rRNA 基因所用引物參考Zhu 等[9]，磺胺類sul1和sul2 基因所用引物參考Lu 等[10]，sul3 基因參考Xiong 等[11]。qPCR 擴(kuò)增體系參照TB Green? Premix Ex Taq? II（Takara，日本）試劑盒說明書，10 μL 的qPCR體系包括5.0 μL 的TB Green Premix ExTaqII（Tli RNaseH Plus）（2×），0.4 μL 的上下游引物（10 μmol·L-1），0.2 μL 的ROX Reference Dye（50×），1.0 μL 的DNA 模板，ddH2O 補(bǔ)充體積至10 μL。擴(kuò)增條件為：95 ℃預(yù)變性30 s；95 ℃變性5 s，60 ℃退火30 s，擴(kuò)增40 個(gè)循環(huán)。每個(gè)樣品設(shè)3 個(gè)平行，以無菌水作為陰性對(duì)照，3 個(gè)平行均被擴(kuò)增出來即為有效擴(kuò)增。

1.5 土壤ARGs 相對(duì)豐度計(jì)算

選取16S rRNA 基因?yàn)閰⒈然?，采用△Ct 方法計(jì)算土壤樣品中ARGs 的相對(duì)豐度，具體計(jì)算公式如下[12]：

式中，Ct（ARG）為目標(biāo)基因的熒光信號(hào)達(dá)到設(shè)定閾值時(shí)所需的循環(huán)數(shù)；Ct（16S rRNA）為參比基因的熒光信號(hào)達(dá)到設(shè)定閾值時(shí)所需的循環(huán)數(shù)；F 值為相對(duì)豐度。

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

采用Microsoft Excel 2007 和SPSS Statistics 13.0軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，并采用Duncan 檢驗(yàn)對(duì)不同樣品進(jìn)行差異性分析（P＜0.05），柱形圖通過OrginPro 2015 軟件繪制。

2 結(jié)果與分析

2.1 不同季節(jié)土壤中ARGs 的多樣性

不同季節(jié)土壤樣品中14 種四環(huán)素類ARGs（tet）和3 種磺胺類ARGs（sul）的檢測(cè)結(jié)果如表1 所示。從表1 可以看出，在17 種目標(biāo)ARGs 中，長(zhǎng)期施用沼液的土壤樣品中可以檢出11～13 種ARGs，其中磺胺類ARGs 的檢出率達(dá)到100%，檢測(cè)到的四環(huán)素類ARGs 主要包括5 種編碼外排泵的基因（tetB、tetG、tetH、tetY、tetZ）和6 種編碼核糖體保護(hù)蛋白的基因（tetM、tetO、tetQ、tetS、tetT、tetW）。在未施用沼液的土壤樣品中檢出6～7 種ARGs，磺胺類ARGs 的檢出率也達(dá)到了100%，但檢測(cè)到的四環(huán)素類ARGs 僅有3～4 種（tetB、tetG、tetY、tetZ）。

表1 不同季節(jié)土壤中檢測(cè)到的抗性基因（n=3）Table 1 Resistance genes detected in the soils in different seasons

長(zhǎng)期施用沼液土壤樣品中的ARGs 多樣性在季節(jié)之間波動(dòng)較明顯，4 個(gè)季節(jié)土壤樣品的ARGs 多樣性表現(xiàn)為春季＞冬季＞夏季＝秋季，春季明顯高于夏季和秋季。未施用沼液土壤樣品中一年四季的ARGs多樣性變化不大，僅在秋季略有下降?？傮w而言，長(zhǎng)期施用沼液比未施用沼液土壤在各季節(jié)均有更高的ARGs 多樣性，說明長(zhǎng)期施用沼液明顯增加了土壤中ARGs 的多樣性。

2.2 不同季節(jié)土壤中四環(huán)素類ARGs 的相對(duì)豐度

通過計(jì)算不同季節(jié)土壤中四環(huán)素類ARGs 總和（sumtet）的相對(duì)豐度，比較了土壤中四環(huán)素類ARGs 的存在水平，結(jié)果如圖1 所示。

從圖1 可以看出，長(zhǎng)期施用沼液的土壤中sumtet的相對(duì)豐度均在10-4數(shù)量級(jí)，未施用沼液的土壤中sumtet的相對(duì)豐度在10-5～10-4數(shù)量級(jí)。長(zhǎng)期施用沼液土壤中四環(huán)素類ARGs 總和的相對(duì)豐度，均顯著高于未施用沼液土壤中相應(yīng)ARGs 總和的相對(duì)豐度（P＜0.05），表明長(zhǎng)期施用沼液會(huì)增加土壤中四環(huán)素類ARGs 的存在水平。從季節(jié)上來看，在長(zhǎng)期施用沼液的土壤中，春季四環(huán)素類ARGs 總和的相對(duì)豐度顯著高于其他3 個(gè)季節(jié)（P＜0.05）。對(duì)于未施用沼液的土壤，四環(huán)素類ARGs 在春季也呈現(xiàn)出較高的相對(duì)豐度。

圖1 不同季節(jié)土壤中四環(huán)素類ARGs 相對(duì)豐度總和Fig.1 Sum of relative abundance of tetracycline ARGs in soil specimens collected in different seasons

進(jìn)一步統(tǒng)計(jì)了長(zhǎng)期施用沼液和未施用沼液土壤中同時(shí)檢測(cè)到的tetG、tetY 和tetZ 等基因的相對(duì)豐度，結(jié)果如圖2 所示。從圖2 可以看出，施用過沼液和未施用過沼液土壤中tetG、tetY 和tetZ 基因的相對(duì)豐度分別在10-6～10-4、10-7～10-6和10-6～10-5數(shù)量級(jí)，表明采集的所有土壤中都有著豐富的tetG和tetZ 基因。與未施用沼液土壤相比，長(zhǎng)期施用沼液土壤中tetG 基因的相對(duì)豐度基本顯著增加，而tetY 基因的相對(duì)豐度則顯著減少（秋季除外），說明長(zhǎng)期施用沼液能夠提高土壤中微生物攜帶tetG 基因的水平。經(jīng)過連續(xù)6 年的沼液施肥，土壤中tetZ 基因的相對(duì)豐度與未施用沼液土壤相比無顯著差異（P＞0.05），表明該基因的相對(duì)豐度沒有受到沼液施用的顯著影響。從季節(jié)上來看，在長(zhǎng)期施用沼液的土壤中，tetG 基因在春季土壤中的相對(duì)豐度要顯著高于其他3 個(gè)季節(jié)（P＜0.05），而tetY 和tetZ 基因的相對(duì)豐度在不同季節(jié)之間均無顯著差異。對(duì)于未施用沼液的土壤，tetG 和tetY 基因在春季也都呈現(xiàn)出較高的相對(duì)豐度，而tetZ 基因的相對(duì)豐度在不同季節(jié)之間無顯著差異，表明該基因的相對(duì)豐度沒有受到季節(jié)的顯著影響。

圖2 不同季節(jié)土壤中四環(huán)素類ARGs 的相對(duì)豐度Fig.2 Relative abundance of tetracycline ARGs in soil specimens collected in different seasons

在長(zhǎng)期施用沼液土壤中還檢測(cè)到tetM、tetO、tetQ、tetS、tetT、tetW 等編碼核糖體保護(hù)蛋白的基因（圖3），其相對(duì)豐度在10-7～10-4數(shù)量級(jí)，這些基因在未施用沼液的土壤中未能檢出。從圖3 可以看出，在春季和秋季，各基因的相對(duì)豐度排序?yàn)閠etM＞tetW＞tetO＞tetQ＞tetT，tetS 基因未檢出；在夏季，各基因的相對(duì)豐度排序?yàn)閠etM＞tetO＞tetW＞tetQ＞tetT，tetS 基因未檢出；在冬季，各基因的相對(duì)豐度排序?yàn)閠etM＞tetO＞tetW＞tetQ＞tetS，tetT 基因未檢出?？梢?，tetM、tetO 和tetW 基因在不同季節(jié)長(zhǎng)期施用沼液土壤中都呈現(xiàn)出較高的相對(duì)豐度，其中tetM 基因的相對(duì)豐度在10-5～10-4數(shù)量級(jí)，tetO 和tetW 基因的相對(duì)豐度在10-5數(shù)量級(jí)，說明這3 種核糖體保護(hù)蛋白基因在長(zhǎng)期施用沼液土壤中的存在水平比其他基因更高。方差分析結(jié)果顯示，tetM和tetW 基因在春季土壤中的相對(duì)豐度要顯著高于其他3 個(gè)季節(jié)，tetO 基因在夏季土壤中的相對(duì)豐度要顯著高于其他3 個(gè)季節(jié)（P＜0.05）。

圖3 長(zhǎng)期施用沼液的土壤中四環(huán)素類ARGs 的相對(duì)豐度Fig.3 Relative abundance of tetracycline ARGs in soils with long-term application of biogas slurry

2.3 不同季節(jié)土壤中磺胺類ARGs 的相對(duì)豐度

通過計(jì)算不同季節(jié)土壤中磺胺類ARGs 總和（sumsul）的相對(duì)豐度，比較了土壤中磺胺類ARGs 的存在水平，結(jié)果如圖4 所示。

圖4 不同季節(jié)土壤中磺胺類ARGs 相對(duì)豐度總和Fig.4 Sum of relative abundance of sulfonamide ARGs in soil specimens collected in different seasons

從圖4 可以看出，長(zhǎng)期施用沼液的土壤中sumsul的相對(duì)豐度在10-4～10-3數(shù)量級(jí)，未施用沼液的土壤中sumsul的相對(duì)豐度在10-5～10-4數(shù)量級(jí)。長(zhǎng)期施用沼液土壤中磺胺類ARGs總和的相對(duì)豐度，均顯著高于未施用沼液土壤中相應(yīng)ARGs 總和的相對(duì)豐度（P＜0.05），表明長(zhǎng)期施用沼液會(huì)增加土壤中磺胺類ARGs 的存在水平。從季節(jié)上來看，在長(zhǎng)期施用沼液的土壤中，磺胺類ARGs 總和的相對(duì)豐度呈現(xiàn)春季＞夏季＞秋季＞冬季的趨勢(shì)，春季的相對(duì)豐度顯著高于其他季節(jié)（P＜0.05）。對(duì)于未施用沼液的土壤，磺胺類ARGs 的相對(duì)豐度在季節(jié)間的波動(dòng)沒有明顯的規(guī)律，但在春季也呈現(xiàn)出較高的相對(duì)豐度。

進(jìn)一步統(tǒng)計(jì)了各土壤樣品中不同種磺胺類ARGs的相對(duì)豐度，其季節(jié)分布特征如圖5 所示，可以看出，在分析的3 種磺胺類ARGs 中，sul1 和sul2 基因在各季節(jié)土壤中的相對(duì)豐度普遍高于sul3 基因，表明采集的所有土壤中都有著豐富的sul1 和sul2 基因。與未施用沼液土壤相比，長(zhǎng)期施用沼液土壤中sul1 和sul2 基因的相對(duì)豐度都顯著增加，這2 種基因在施用沼液土壤中的相對(duì)豐度均達(dá)到了10-4數(shù)量級(jí)水平，說明長(zhǎng)期施用沼液能夠提高土壤中微生物攜帶sul1 和sul2 基因的水平。從季節(jié)上來看，在長(zhǎng)期施用沼液的土壤中，sul2 基因在春季土壤中的相對(duì)豐度比在其他3 個(gè)季節(jié)高出1 個(gè)數(shù)量級(jí)，差異達(dá)到顯著性水平（P＜0.05）；夏季sul3 基因的相對(duì)豐度顯著高于其他3 個(gè)季節(jié)。對(duì)于未施用沼液的土壤，sul1、sul2 和sul3 基因均在春季土壤中具有更高的相對(duì)豐度。

圖5 不同季節(jié)土壤中磺胺類ARGs 的相對(duì)豐度Fig.5 Relative abundance of sulfonamide ARGs in soil specimens collected in different seasons

3 討論

在ARGs 多樣性分析中，發(fā)現(xiàn)在未施用沼液的土壤樣品中同樣能檢測(cè)到一些ARGs，包括tetB、tetG、tetY、tetZ、sul1、sul2 和sul3，這些ARGs 都存在于長(zhǎng)期施用沼液的土壤樣品中，說明這些ARGs 在土壤中廣泛存在，可能是土壤中的背景ARGs，也可能是污染殘留造成的。Knapp 等[13]研究發(fā)現(xiàn)在無抗生素污染歷史的土壤中就已經(jīng)存在四環(huán)素類、β-內(nèi)酰胺類和大環(huán)內(nèi)酯類ARGs；Brion 等[14]調(diào)查發(fā)現(xiàn)，在使用抗生素前的土壤樣本中同樣能檢測(cè)到tetC 等四環(huán)素類ARGs。Tang 等[15]對(duì)稻田土壤的研究也發(fā)現(xiàn)不施糞肥土壤中ARGs 的存在。以上研究均表明土壤環(huán)境中本身就存在一定濃度的ARGs，這些ARGs可能是由土壤中固有的抗生素抗性微生物所攜帶的。相關(guān)研究表明，土壤中的一些土著微生物為了能夠更好地適應(yīng)環(huán)境，自身會(huì)產(chǎn)生低濃度的抗生素以抵制競(jìng)爭(zhēng)者的生存[16]。

分析沼液施用后土壤中ARGs 的變化情況，發(fā)現(xiàn)與未施用沼液的土壤相比，長(zhǎng)期施用沼液的土壤中具有更高的ARGs 多樣性，檢測(cè)出了7 種后者不具有的ARGs，包括tetH、tetM、tetO、tetQ、tetS、tetT 和tetW，這些ARGs 可能是由沼液帶入土壤的外源抗性基因。有研究指出，tetM、tetO、tetQ、tetS、tetT、tetW 等四環(huán)素類ARGs 在畜禽糞便中較為常見[17-18]。根據(jù)定量分析結(jié)果可知，與未施用沼液的土壤相比，長(zhǎng)期施用沼液后土壤中四環(huán)素類和磺胺類ARGs的總豐度呈顯著增加的趨勢(shì)（P＜0.05），表明長(zhǎng)期沼液的施用顯著提高了土壤中ARGs 的整體存在水平。這種增加可能是由于沼液提供的養(yǎng)分引起土壤中某些攜帶四環(huán)素類或磺胺類ARGs 微生物的大量繁殖，也有可能是沼液中攜帶的ARGs 及抗性微生物在土壤中的存活和繁殖引起的。Han 等[19]研究表明，施用動(dòng)物糞肥通過影響土壤的理化特性，改變土壤中的細(xì)菌數(shù)量，進(jìn)而影響土壤中ARGs 的存在水平。大量研究也表明，有機(jī)肥的施用會(huì)增加農(nóng)田土壤中ARGs 的多樣性和豐度[9,20-21]。然而，本研究發(fā)現(xiàn)長(zhǎng)期施用沼液土壤中tetY 和sul3 基因在個(gè)別季節(jié)的相對(duì)豐度并沒有高于未施用沼液土壤，tetZ 基因受到沼液施用的影響相對(duì)較小。這可能與土壤中宿主微生物數(shù)量的消長(zhǎng)以及基因發(fā)生水平轉(zhuǎn)移有關(guān)；此外，不同ARGs 之間或許存在互補(bǔ)作用，某一類ARGs 豐度的增加也可能伴隨著另一類ARGs 豐度的減少同時(shí)發(fā)生[22]。

分析不同季節(jié)土壤中ARGs 的分布情況，發(fā)現(xiàn)在長(zhǎng)期施用沼液的土壤中，tetG、tetM、tetQ、tetW 和sul2 等基因的相對(duì)豐度具有明顯的季節(jié)特點(diǎn)，表現(xiàn)出春季的相對(duì)豐度要顯著高于其他3 個(gè)季節(jié)（P＜0.05）的趨勢(shì)，尤其sul2 基因在春季土壤中的相對(duì)豐度比在其他3 個(gè)季節(jié)高出1 個(gè)數(shù)量級(jí)。在未施用沼液的土壤中，檢測(cè)到的tetG、tetY、sul1、sul2 和sul3 等基因的相對(duì)豐度在季節(jié)之間波動(dòng)更明顯，但季節(jié)之間所呈現(xiàn)出的規(guī)律與長(zhǎng)期施用沼液土壤相類似，都在春季土壤中具有更高的相對(duì)豐度?？傮w而言，春季各土壤樣品中四環(huán)素類ARGs 總和與磺胺類ARGs總和的相對(duì)豐度均要高于其他季節(jié)土壤。這可能是由于春季的氣候條件更適合土壤中tetG、sul2 等基因宿主細(xì)菌的生存，進(jìn)而導(dǎo)致四環(huán)素類和磺胺類ARGs的總豐度在春季更高。張俊等[23]在研究長(zhǎng)期施用四環(huán)素殘留豬糞的土壤樣本中發(fā)現(xiàn)，適宜的氣候條件對(duì)耐藥菌及ARGs 的生成具有很好的促進(jìn)作用，但該研究指出秋季涼爽的氣候更適合土壤中ARGs的生成。也有相關(guān)研究表明，冬季較低的氣溫減緩了抗生素在土壤環(huán)境中的降解，加上冬季養(yǎng)殖場(chǎng)為預(yù)防疫病而加大了抗生素的使用量，導(dǎo)致冬季土壤中ARGs 豐度更高[24]。而本研究發(fā)現(xiàn)，除了春季，不同ARGs 在其他季節(jié)之間并沒有一致的變化規(guī)律，關(guān)于季節(jié)對(duì)土壤中不同ARGs 形成的影響機(jī)制仍需進(jìn)一步長(zhǎng)期的跟蹤調(diào)查。

分析土壤中各種目標(biāo)ARGs 的分布情況，發(fā)現(xiàn)四環(huán)素類的2 種ARGs（tetG 和tetZ）與磺胺類的2 種ARGs（sul1 和sul2），在不同季節(jié)長(zhǎng)期施用沼液和未施用沼液土壤樣品中都占據(jù)較高的比例，屬于優(yōu)勢(shì)ARGs。此外，tetM、tetO 和tetW 基因也在長(zhǎng)期施用沼液的土壤中占據(jù)較高的比例，其中tetM 基因達(dá)到10-5～10-4數(shù)量級(jí)水平，tetO 和tetW 基因在10-5數(shù)量級(jí)水平。樓晨露[24]在研究長(zhǎng)期施用豬糞的稻田土壤中發(fā)現(xiàn)，磺胺類ARGs 中sul1 和sul2 基因的相對(duì)豐度較高，介于10-4～10-3數(shù)量級(jí)水平，與本研究長(zhǎng)期施用沼液土壤中相應(yīng)基因的豐度水平相當(dāng)；四環(huán)素類ARGs 中tetG 基因的豐度最高，介于10-4～10-3數(shù)量級(jí)水平，高于本研究（10-5～10-4數(shù)量級(jí)水平）。彭雙等[22]在研究連續(xù)施用豬糞的稻麥輪作土壤中也發(fā)現(xiàn)，在施用豬糞土壤中tetG 基因的相對(duì)豐度最高（接近10-2數(shù)量級(jí)水平），其次是tetZ 基因（10-3數(shù)量級(jí)水平），均高于本研究（圖2）。相關(guān)研究也表明，四環(huán)素類ARGs 中tetM、tetO、tetW 基因在養(yǎng)豬場(chǎng)排泄物及周邊土壤中大量存在[25]。以上研究表明，四環(huán)素類ARGs 中tetG、tetZ、tetM、tetO 和tetW 基因與磺胺類ARGs 中sul1 和sul2 基因在土壤中的傳播風(fēng)險(xiǎn)高于其他ARGs，因此養(yǎng)殖場(chǎng)仍需優(yōu)化發(fā)酵工藝以降低糞肥農(nóng)用帶來的ARGs 污染風(fēng)險(xiǎn)。

4 結(jié)論

（1）長(zhǎng)期施用沼液使土壤中ARGs 的多樣性增加，在所分析的17 種目標(biāo)ARGs 中，tetG、sul1 和sul2 基因的相對(duì)豐度基本顯著增加，土壤中四環(huán)素類和磺胺類ARGs 的整體相對(duì)豐度均明顯提高，因此沼液農(nóng)用前的ARGs 污染水平評(píng)價(jià)顯得尤為重要。

（2）四環(huán)素類和磺胺類ARGs 在土壤中的存在水平具有明顯的季節(jié)特點(diǎn)，兩類ARGs 總豐度在春季土壤中更高，春季長(zhǎng)期施用沼液和未施用沼液土壤中都有著豐富的tetG 和sul2 基因，關(guān)于季節(jié)對(duì)土壤中ARGs 形成的影響機(jī)制仍需進(jìn)一步探究。

（3）在所分析的17 種目標(biāo)ARGs 中，tetG、tetZ、tetM、tetO、tetW、sul1 和sul2 基因在長(zhǎng)期施用沼液的不同季節(jié)土壤中均大量存在，為土壤中的優(yōu)勢(shì)ARGs，未來需加強(qiáng)優(yōu)化畜禽廢棄物發(fā)酵工藝，以削減ARGs 帶來的污染風(fēng)險(xiǎn)。