4種網絡交易的糞肥中抗生素抗性基因多樣性分析和傳播風險評估

劉 燕,朱 冬,2，何 燕，楊宇帆，彭 超，錢永明，王 東,唐 赟,路 璐*

(1.西華師范大學環境科學與工程學院,四川 南充 637009;2.中國科學院生態環境研究中心,北京 100085;3.西華師范大學生命科學學院,四川 南充 637002;4.西華師范大學法學院,四川 南充 637002)

我國是抗生素生產和使用大國，每年抗生素生產量約21萬t，其中9.7萬t(占年總產量的46.1%)抗生素用于畜禽養殖業。在畜禽養殖業中抗生素被應用于疾病預防和治療，同時作為飼料添加劑促進動物生長。然而，長期濫用或不合理使用抗生素會誘導畜禽糞便中產生大量的抗生素耐藥微生物和抗生素抗性基因。大量研究發現，畜禽糞便中含有高度富集的抗生素抗性基因，這些糞便會對養殖區域及周圍環境造成潛在的抗性基因污染。同時，抗性基因的污染不同于傳統的化學污染物，抗性基因可通過各種可移動遺傳元件，如整合子、轉座子或質粒等，通過轉化、轉導和接合等基因水平轉移方式在同一菌種的個體之間甚至在不同菌種之間傳遞，并通過轉化、轉導和接合的途徑進入其他環境細菌和致病菌中，對人類健康和生態環境造成潛在的威脅。

目前，我國對產生的大量畜禽糞便主要是通過肥料化技術進行資源化利用。雖然堆肥技術可以消減畜禽糞便中抗生素的殘留、抗生素抗性基因的豐度以及致病微生物，但是一些研究發現堆肥后糞肥中仍有較高豐度的抗性基因存在，表示堆肥可能是良好的生物反應器，能促進其抗性微生物的增殖。此外，農田土壤中施用糞肥顯著提高了抗生素抗性基因的豐度和多樣性,可能導致耐藥性微生物和抗性基因在農田生態系統中的擴散與傳播。然而，除了批量地施用于農田土壤外，這些糞肥還可用于家庭種植花卉或者果蔬。目前隨著網絡交易平臺的成熟，糞肥也大量在網絡平臺售賣，比如淘寶平臺糞肥銷售廠家約有700家，每天交易量達數百起。這種網絡銷售途徑可能是抗生素抗性基因在更廣范圍上傳播的新途徑，進而導致被施用土壤中抗生素抗性基因的污染。此外，有研究表明肥料中的抗生素抗性基因能遷移至地上種植的植物或蔬菜根際、葉際中，進而可能通過食物鏈進入人體，對人類健康造成潛在的威脅。因此，探究網絡銷售的糞肥中抗生素抗性基因的多樣性和豐度，并評估其潛在的傳播風險顯得尤為緊迫。

本研究以在網絡平臺購買的牛糞、羊糞、雞糞和雞羊混合糞4種糞肥作為供試材料，采用高通量熒光定量PCR技術對4種網絡交易的糞肥中抗生素抗性基因的豐度和多樣性進行比較分析，評估網絡交易的糞肥是否為抗生素抗性基因的傳播途徑，進而評價其傳播風險。本文采用的296個引物涵蓋了283個抗生素抗性基因(針對9種主要的抗生素：氨基糖苷類、β-內酰胺類、氯霉素類、MLSB類、磺胺類、四環素類、萬古霉素類、多重抗性類和其他類)、12個可移動遺傳元件標記基因和細菌16S rRNA基因，以期對網絡源糞肥中抗生素抗性基因污染情況有較全面的認識，并為可能的潛在生態環境風險評價提供理論支持。

1 材料與方法

1.1 供試糞肥信息

于2018年10月8日在淘寶在線交易平臺，分別在4個月銷量大于1 300 kg的廠家購買了牛糞肥、雞糞肥、羊糞肥和雞羊混合肥4種發酵糞肥，發貨地址分別是河南鄭州、河南安陽、四川成都和安徽合肥。其中，雞糞肥、羊糞肥、雞羊混合糞肥是將其糞便經過堆肥發酵而形成的畜禽類糞肥；牛糞肥為蚯蚓堆肥，即使用發酵過的純牛糞堆成壟狀飼養蚯蚓，經蚯蚓消化后的肥料，再經過堆肥過篩再次發酵，形成松散的均勻粒狀的牛糞肥。將各種糞肥過2 mm孔徑篩并裝入塑封袋中置于-20℃冰箱，以備后續分子生物學分析。

1.2 總DNA的提取

采用Fast DNA Spin Kit for Soil試劑盒(MP Biomedicals,Cleveland,OH,USA) 提取糞肥中的微生物總DNA，每個樣品3個重復。樣品中DNA濃度和質量采用NanoDrop分光光度計 (Nanodrop Technologies Inc.,Wilmington,DE) 進行測定和評估，所提取的DNA樣品的A260/A280值在1.8～2.0之間，提取的DNA樣品保存于-20℃以備后續分析。

1.3 熒光定量PCR分析

實時熒光定量PCR擴增使用SYBRPremix ExTaq TM(TliRNaseH Plus)試劑盒在CFX96 Optical Real-Time PCR System熒光定量儀(Bio-Rad Laboratories Inc.,Hercules,CA,USA)上進行分析,對所有樣品的16S rRNA基因豐度進行絕對定量。實時熒光定量PCR標準曲線采用含有16S rRNA基因的質粒，質粒標準曲線按10倍梯度稀釋，得到7個數量級的標準曲線，其標準曲線質粒濃度的變化范圍為10～10copies/μL。實時熒光定量PCR擴增體系為10 μL：5.0μL SYBRPremixEx Taq TM(TliRNaseH Plus)，1.0 μL DNA模板，10 μmol/L正、反向引物各0.1 μL，3.8 μL滅菌雙蒸水，陰性對照采用滅菌雙蒸水代替樣品DNA，擴增效率為90.6%～101.3%，

R

在0.97～0.99的范圍內。對4種糞肥的3個重復樣品開展微生物總DNA的定量測定。

1.4 高通量熒光定量PCR測定抗生素抗性基因和可移動遺傳元件

糞肥中的抗生素抗性基因在SmartChip Real-Time PCR Systems(WaferGen Inc.,USA)高通量實時熒光定量PCR平臺測定。抗性基因引物參照文獻[5]，使用了283對抗性基因引物、8對轉座酶基因、4對整合酶-整合子基因和1對細菌通用16S rRNA基因引物。高通量熒光定量PCR擴增體系的體積為100 nL，其中各試劑終濃度為：1×LightCycler 480 SYBRGreen IMaster Mix (Roche，USA),Nuclease-free PCR-Grade water，1 μg/μL的牛血清蛋白(BSA)，5 ng/μL的DNA模板和1 μmol/L的上下游引物。高通量熒光定量PCR反應條件為：95℃預變性10 min；95℃變性30 s，60℃退火延伸30 s，總共40個循環；熱循環儀Cycler的程序自動升溫進行熔解曲線分析。高通量熒光定量PCR每個芯片都有不添加DNA模板(用滅菌超純水代替)的陰性對照。高通量熒光定量PCR反應得到的數據通過Cycler預設定的篩選條件(擴增效率介于1.8～2.2)進行導出。根據Smart Chip Real-Time PCR Systems的靈敏度和精確度，確定循環次數

C

值為31次時作為儀器的檢測限，每個樣品的2個重復均被檢測出，且擴增效率在90%～110%內被視為有效數值。

1.5 統計分析

抗生素抗性基因和可移動遺傳元件的相對拷貝數參照Zhu等的研究，采用公式(1)進行估算，并根據Stalder等的研究，采用公式(2)計算抗生素抗性基因的絕對豐度。具體計算公式如下：

C

=10[(31-)(103)]

(1)

基因的絕對拷貝數

(2)

式中:

C

為抗生素抗性基因的相對拷貝數(copies/bacterial cell)；

C

為PCR反應收集到特定熒光時的循環次數(次)；

A

為抗生素抗性基因的絕對豐度(copies/g)，16S rRNA基因的拷貝數(copies/g)。

本文利用Microsoft Excel 2016進行計算分析，使用SPSS 23.0(SPSS Inc.，USA)進行單因素方差顯著性分析和相關性統計分析，利用 Origin Pro 2018(Origin Lab，USA)制圖。

2 結果與討論

2.1 4種網絡交易的糞肥中檢測出的抗生素抗性基因和可移動遺傳元件的種類

4種網絡交易的糞肥中檢測出的抗生素抗性基因和可移動遺傳元件的種類，見圖1。

圖1 4種網絡交易的糞肥中檢測出的抗生素抗性基因和可移動遺傳元件的種類Fig.1 Diversity of antibiotic resistance genes and Mobile Genetic Elements(MGEs) in the four manure

由圖1可以看出：

(1) 4種網絡交易的糞肥中共檢測出165種不同的抗生素抗性基因亞型，其中羊糞肥中檢出的抗生素抗性基因最多，有130種，牛糞肥和雞羊混合糞肥中檢出的抗生素抗性基因最少，均為89種，雞糞肥中檢出98種抗生素抗性基因[見圖1(a)]，這一發現與施用豬糞的水稻土壤(107種)、大型養豬場糞便(136種)、雞糞肥(127種)和牛糞肥(107種)等環境中檢出的抗生素抗性基因的個數相當，表明網絡交易的糞肥中存在多樣性豐富的抗生素抗性基因。4種網絡交易的糞肥中共檢出的可移動遺傳元件的數量為10種，包括3個整合子和7個轉座子，其中羊糞肥中檢出的可移動遺傳元件的個數最多，為8種[見圖1(a)]。已有研究表明，抗生素抗性基因在環境中的水平基因轉移與可移動遺傳元件密切相關，抗生素抗性基因可以隨著可移動元件(如整合子、質粒等)進行水平基因轉移，使得抗性基因在不同種屬微生物間傳播，更是加劇了抗性基因在不同環境介質中的傳播風險。可移動遺傳元件作為基因水平轉移重要的分子傳播元件，在4種網絡交易的糞肥中均被檢出，表明在這4種糞肥中存在抗生素抗性基因水平轉移的潛在風險。

(2) 4種網絡交易的糞肥中檢出的抗生素抗性基因主要以氨基糖苷類、大環內脂類-林肯酰胺類-鏈陽性菌素B類(MLSB類)、多重耐藥類和四環素類為主，分別占整體的比例為22%、21%、17%和14%[見圖1(b)]，該發現與以往在牛糞肥、雞糞肥和豬糞肥中檢出的抗生素抗性基因的抗性類型相同。這可能與我國養殖業中使用的抗生素類型相關，因為抗生素的添加和殘留是畜禽糞便抗生素抗性基因富集的主要原因之一，即使抗生素的濃度非常低，也能對環境中的抗性微生物和抗性基因產生選擇性壓力。且也有數據統計表明，我國養殖業飼料中主要添加隸屬于氨基糖苷類的新霉素、慶大霉素，四環素類的土霉素、金霉素，林肯酰胺類的林可霉素，大環內酯類的泰樂菌素等抗生素。

(3) 4種網絡交易的糞肥中抗生素抗性基因的抗性類型組成存在一定的差異。在羊糞肥中檢出的氨基糖苷類(30±0.71)、MLSB類(29±0.71)、多重耐藥類(26±0.00)、四環素類(16±0.71)抗性基因的個數顯著高于其他類糞肥(

p

<0.05)，而在牛糞肥中檢出的氯霉素類(3±0.71)抗性基因的個數高于其他類糞肥，但4種糞肥中檢出的β-內酰胺類、萬古霉素類和其他類抗性基因的種類沒有顯著差別。4種網絡交易的糞肥中抗生素抗性基因抗性類型的不同，可能是由于其糞肥對應的養殖場抗生素使用情況不同，且與堆肥的方式不同有關。宋婷婷等研究發現，堆肥方式的不同(溫度、濕度、pH值、堆肥材料的性質、糞便的類型等)對糞肥中抗生素抗性基因的抗性類型有顯著的影響。根據對抗生素抗性機制的不同，將檢出的抗生素抗性基因分為抗生素失活、核糖體保護、外排泵作用和其他/未知4大類。4種網絡交易的糞肥中抗生素主要的抗性機制為抗生素失活(43%)，其次為外排泵作用(31%)和核糖體保護(23%)，其他/未知僅為3%[見圖1(c)]，這與Zhu等在3個大型養豬場豬糞及周邊土壤的研究結果類似，其在豬糞中也發現抗生素失活和外排泵作用為主導的抗生素抗性機制。

2.2 4種網絡交易的糞肥中抗生素抗性基因的豐度分析

4種網絡交易的糞肥中抗生素抗性基因和可移動遺傳元件的絕對豐度和相對拷貝數，見圖2。

由圖2可以看出：

(1) 4種糞肥中抗生素抗性基因的絕對豐度介于7.18×10～3.99×10copies/g(干重)之間，且4種糞肥中抗生素抗性基因的絕對豐度間有顯著差異(

p

<0.01)。其中，牛糞肥中抗生素抗性基因的絕對豐度最低，為7.18×10copies/g(干重)；羊糞肥中抗生素抗性基因的絕對豐度最高，為3.99×10copies/g(干重)，是牛糞肥的5 555.2倍；雞羊混合糞肥中抗生素抗性基因的絕對豐度也顯著高于牛糞肥和雞糞肥，為2.66×10copies/g(干重)[見圖2(a)]。羊糞肥和雞羊混合糞肥中抗生素抗性基因的絕對豐度顯著高于以往在農田土壤[10～10copies/g(干重)]、湖泊和水庫[10～10copies/g(干重)]、河口底泥[10～10copies/g(干重)]等環境樣品中所檢測出的抗生素抗性基因的豐度，表明糞肥的網絡交易可能導致抗生素抗性基因的高能擴散，其攜帶的抗性基因能通過施肥等方式擴散至不同的環境中。此外，該些網絡交易的羊糞肥、雞糞肥和雞羊混合糞肥所檢測出的抗生素抗性基因豐度也均高于Zhang等在雞糞肥、牛糞肥和豬糞肥[10～10copies/g(干重)]中所檢出的抗生素抗性基因的豐度，進一步揭示網絡交易的糞肥存在抗生素抗性基因傳播的風險，因此亟需制定網絡交易糞肥抗性基因的控制標準。

圖2 4種網絡交易的糞肥中抗生素抗性基因和可移動遺傳元件的絕對豐度和相對拷貝數Fig.2 Absolute abundance and relative abundance of antibiotic resistance genes and mobile genetic elements in the four types of manure

(2) 從各抗性基因類型的絕對豐度計算結果來看，4種網絡交易的糞肥中主導的抗生素抗性基因類型是氨基糖苷類(47.10%)、磺胺類(22.04%)和多重耐藥類(20.26%)，占總的絕對豐度的89.4%。這一結果與Zhang等對雞糞肥、牛糞肥和豬糞肥中抗性基因類型豐度的研究結果相似，其也發現氨基糖苷類(20%)為豐度最高的抗性基因。而氨基糖苷類和磺胺類抗性基因的富集也可能與養殖業中這兩類抗生素的廣泛使用相關。多重耐藥類抗性基因在4種網絡交易的糞肥中的發現需要引起更多關注，已有研究表明多重耐藥類抗性基因的傳播對人類和動物抗感染治療是極大的威脅。多重耐藥類抗性基因的檢出和較高豐度進一步表明需對網絡交易的糞肥做微生物安全的防控。此外，不同類型的抗性基因絕對豐度在4種網絡交易的糞肥中的比例也極顯著不同(

p

<0.01)，如：雞羊混合糞肥中氨基糖苷類抗性基因的絕對豐度顯著高于其他3種糞肥(

p

<0.05)，其絕對豐度分別是牛糞肥、雞糞肥和羊糞肥的10 266.1倍、110.7倍和1.2倍；磺胺類抗性基因的絕對豐度在羊糞肥中最高，是牛糞肥、雞糞肥和羊糞肥的13 354.23倍、6.87倍和1.62倍。不同糞肥中抗生素抗性基因絕對豐度的差異可能與不同養殖場抗生素使用類型、用量以及后續發酵處理的差異等因素相關。如Zhu等在分析比較北京、嘉興和莆田3個養豬場糞肥中抗生素抗性基因的絕對豐度時也發現，3個養殖場糞肥中抗生素抗性基因的種類和豐度組成不同，該結果說明不同網絡交易的糞肥中抗生素抗性基因的承載量以及類型、豐度有顯著的差異。

(3) 為評估抗生素抗性基因在整個細菌群落中的豐度及其富集程度，本文計算了4種網絡交易的糞肥中平均每個細菌細胞中抗生素抗性基因的相對拷貝數[見圖2(b)]，相對拷貝數越高，表明該樣品中抗性基因的富集密度和傳播風險越大。4種網絡交易的糞肥中抗生素抗性基因的相對拷貝數范圍在0.013～0.727 copies/bacterial cell之間，平均每個細胞有0.253個抗性基因。4種網絡交易的糞肥中抗生素抗性基因的相對拷貝數變化規律與抗生素抗性基因絕對豐度的變化規律一致，羊糞肥中抗生素抗性基因的相對拷貝數仍是最高，為0.727 copies/bacterial cell，牛糞肥中抗生素抗性基因的相對拷貝數最低，為0.013 copies/bacterial cell,僅是羊糞肥的1/56。4種網絡交易的糞肥中抗生素抗性基因的相對拷貝數介于前期研究的堆肥污泥(0.77～2.50 copies/bacterial cell)、河口底泥(0.21～1.93 copies/bacterial cell)以及雞糞肥、牛糞肥和豬糞(0.15～1.81 copies/bacterial cell)等環境樣品中抗生素抗性基因的相對拷貝數檢測值。而這4種網絡交易的糞肥中抗生素抗性基因的相對拷貝數均顯著高于何燕等在7個稻田土壤糞肥中抗生素抗性基因的檢測范圍(0.006～0.025 copies/bacterial cell)，同時也高于其他農田土壤(平均值為0.074 copies/bacterial cell)、淡水湖泊(0.01～0.38 copies/bacterial cell)等環境樣品中抗生素抗性基因在微生物群落中的富集程度。這些結果表明網絡交易的糞肥中存在抗生素抗性基因的富集，具有較高的傳播風險。

(4) 從各類型抗性基因的相對拷貝數計算結果來看，4種網絡交易的糞肥中的抗生素抗性基因種類以氨基糖苷類(39.91%)、磺胺類(24.71%)和多重耐藥類(22.60%)3類為主導。4種網絡交易的糞肥中不同抗生素抗性基因的相對拷貝數也有顯著的差異(

p

<0.01)。如在牛糞肥、羊糞肥和雞羊混合糞肥3種糞肥中，主導的抗生素抗性基因為氨基糖苷類，占總抗性基因相對拷貝數的23.64%～65.57%，其次是多重耐藥類，占總抗性基因相對拷貝數的9.55%～28.11%。然而，雞糞肥中主導的抗性基因為磺胺類(72.35%)，其次才是氨基糖苷類(9.17%)，兩者占總抗性基因相對拷貝數的81.54%。值得注意的是，在4種網絡交易的糞肥中，牛糞肥中抗生素抗性基因的種類和豐度均極顯著低于其他糞肥(

p

<0.01)，且顯著低于以往在牛糞肥等糞肥中檢出的抗生素抗性基因種類和豐度。本研究中的牛糞肥經過了蚯蚓進一步分解，其中的有機質進一步得到分解。有研究表明，在長期施用污泥或雞糞肥的土壤中發現，土壤動物蚯蚓腸道中的抗生素抗性基因豐度顯著增加。而本研究中網絡購買的牛糞肥蚯蚓堆肥中較低的抗生素抗性基因豐度可能受很多未知因素的影響，其抗生素抗性基因豐度是否與抗性基因向蚯蚓腸道中遷移有關、抑或蚯蚓降低了糞肥中抗性基因的豐度，或是否由于抗性基因本身在糞源中含量較低，仍需要未來進一步設計試驗來進行研究驗證。該結果也進一步表明網絡交易的不同糞肥中含有的抗生素抗性基因的豐度有較大的差異。

2.3 4種網絡交易的糞肥中可移動遺傳元件豐度比較及其相對拷貝數與抗生素抗性基因相對拷貝數的相關性分析

由圖2(c)、(d)可見，4種網絡交易的糞肥中可移動遺傳元件(轉座子和整合子)的絕對豐度和相對拷貝數分別介于2.61×10～2.90×10copies/g(干重)和0.007～0.088 copies/bacterial cell之間，并且不同的糞肥其絕對豐度和相對拷貝數均有顯著性的差異(

p

<0.05)。4種網絡交易的糞肥中可移動遺傳元件的絕對豐度顯著高于河口底泥[5.66×10～4.23×10copies/g(干重)]，與稻田土壤的可移動遺傳元件絕對豐度[1.14×10～2.29×10copies/g(干重)]相當。本研究中的糞肥樣品中可移動遺傳元件的相對豐度顯著高于森林土壤，且其每克糞肥樣品中可移動遺傳元件的絕對豐度比受人為干擾極少的北極圈土壤高出10～10copies/g(干重)。羊糞肥中可移動遺傳元件的絕對豐度和相對拷貝數均顯著高于其他糞肥樣品，分別高出其他3種糞肥約4～1 110倍和6～12倍。可移動遺傳元件作為基因水平轉移的重要傳播元件，表明羊糞肥中抗性基因可能有較高的傳播風險。4種網絡交易的糞肥中各抗性基因相對拷貝數與可移動遺傳元件相對拷貝數的相關性分析結果，見表1。

表1 4種網絡交易的糞肥中各抗生素抗性基因相對拷貝數與可移動遺傳元件相對拷貝數的相關性分析Table 1 Pearson correlation for the relative copy mumber of antibiotic resistance genes and mobile genetic elements in the four online-trading manure

通過分析4種網絡交易的糞肥中各抗性基因相對拷貝數與可移動遺傳元件相對拷貝數的相關性(見表1)，結果表明：4種網絡交易的糞肥中總抗性基因相對拷貝數與總可移動遺傳元件的相對拷貝數呈極顯著正相關(

p

<0.01)，這一結果與在水體、雞糞、牛糞和豬糞、施肥后的土壤等不同環境樣品中的研究結果相同，表明抗生素抗性基因可通過可移動遺傳元件進行基因的水平轉移，使得抗生素抗性基因在不同微生物間遷移傳播，進而增加其傳播風險，說明這些糞肥中抗生素抗性基因的水平轉移可能是抗性基因遷移傳播的重要途徑。本研究還發現，多重抗藥類抗性基因、四環素類抗性基因、氯霉素類抗性基因、MLSB類抗生素抗性基因的相對拷貝數與可移動遺傳元件的相對拷貝數都呈極顯著正相關(

p

<0.01)，這一結果與dela Cruz等的研究結果相似，其發現四環素類和氯霉素類抗性基因相對拷貝數與可移動遺傳元件相對拷貝數具有一定的相關性，進一步表明這4種糞肥中可移動遺傳元件可能促進了這些抗生素抗性基因的遷移、傳播和擴散。此外，所有的抗生素抗性基因類型都與整合子有顯著的相關性(

p

<0.05)。整合子被認為是導致耐藥基因在細菌間水平轉移的重要原因之一，表明抗性基因的水平轉移可能在糞肥中抗生素抗性基因的富集中起著關鍵的作用，并可能加劇糞肥中抗性基因傳播的風險。然后，這些糞肥在施入土壤后，其中的抗生素抗性基因在土壤中的歸趨行為、或傳播至地上植物的遷移行為等，仍需要進一步開展模擬試驗研究來直接揭示抗性基因的傳播行為模式和風險。

2.4 4種網絡交易的糞肥中共有、獨有和主導的抗生素抗性基因種類分析

4種網絡交易的糞肥中共有、獨有和主導的抗生素抗性基因種類分析結果，見圖3和圖4。

圖3 4種網絡交易的糞肥中檢出的抗生素抗性基因和可移動遺傳元件種類數的韋恩圖Fig.3 Venn based on the number of antibiotic tresisance genes and mobile genetic elements detected in the four online-trading manure

圖4 4種網絡交易的糞肥中抗生素抗性基因和可移動遺傳元件的相對豐度熱圖Fig.4 Heat map based on the relative abundance of antibiotic resistance genes and mobile genetic elementsin the four online-trading manure

由圖3和圖4可見，4種糞肥中共有的抗生素抗性基因亞型有35種，可移動遺傳元件有4種。4種網絡交易的糞肥中共有的主導抗性基因亞型和可移動遺傳元件為

sul

2(磺胺類)，

aadA

-02、

aadA

2-03、

aadA

-01、

aadA

1和

aadA

2-02(氨基糖胺類)，

tetL

-02(四環素類)，

intI

-1

LC

(整合子類)。其中，

sul

2抗性基因在牛糞肥、雞糞肥、羊糞肥和雞羊混合糞肥中分別占總抗性基因相對拷貝數的4.49%、58.26%、18.67%和18.56%，

sul

2抗性基因已在河流沉積物、土壤、養豬場豬糞及其周邊土壤等環境樣品中檢出并且豐度較高，這可能是由于

sul

2抗性基因通常定位于較小的非接合型質粒或較大的可移動多重抗性質粒上，而質粒作為可移動元件是抗生素抗性基因水平基因轉移的重要途徑，可能導致了

sul

2抗性基因在各環境條件下的廣泛分布；

aadA

類基因在牛糞肥、雞糞肥、羊糞肥和雞羊混合糞肥中占總抗性基因相對拷貝數的5.30%、5.36%、31.81%和61.97%，

aadA

類抗性基因是在污水、廢水和動物糞便等環境樣品中檢出最為頻繁的一類基因，且

aadA

類抗性基因通常存在于整合子類

intI

-1的基因盒中，

aadA

類抗性基因在4種網絡交易的糞肥中均被檢出進一步表明糞肥中抗性基因存在很高的水平基因遷移風險；編碼外排泵類四環素抗性基因

tetL

-02在4種網絡交易的糞肥中均被檢出，并且有較高的豐度，

tetL

-02抗性基因在豬糞、土壤、廢水等環境樣品中均被檢出，此外有研究發現，

tetL

-02抗性基因與轉座子

tnpA

抗性基因呈顯著正相關，而轉座子作為可移動遺傳元件可促進基因水平轉移，進而增加其潛在的傳播風險；

intI

-1

LC

抗性基因是對抗生素抗性基因水平遷移有關鍵作用的整合子基因，

intI

-1

LC

類整合子是判定環境是否受人為源污染的標志基因之一，此外

intI

-1

LC

抗性基因可與很多抗生素抗性基因的基因盒進行整合連接，如氨基糖苷類抗性基因

aadA

、四環素類抗性基因

tetG

等都在

intI

-1類整合子的基因盒中發現，從而使細菌產生多重耐藥性。此外，不同糞肥中的主導抗性基因也各不相同(見圖4)。如牛糞肥中主導的抗性基因亞型為

aacC

(氨基糖胺類)、

mphA

-01(

MLSB

類)、

sul

2(磺胺類)和

oprD

(多重耐藥類)，占總抗性基因相對豐度的45.19%；而羊糞肥主導的抗性基因亞型則為

sul

2(磺胺類)、

qacEdelta

1-02和

mexF

(多重耐藥類)、

aadA

9-02、

aadA

9-01和

aadA

2-01(氨基糖胺類)，占總抗性基因相對豐度的58.15%。4種網絡交易的糞肥中均檢出其獨有的抗生素抗性基因(見圖3)。如羊糞肥中檢出24種其獨有的抗性基因，其中多重耐藥類

mexF

抗性基因的相對豐度最高，占總抗性基因相對豐度的28.12%。有研究發現

mexF

抗性基因在施用了糞肥的土壤中顯著富集；黃福義等也研究發現施用了豬糞的水稻土壤中

mexF

基因的相對豐度相對于對照組土壤增加了1 791倍。牛糞肥中檢出10種獨有的抗性基因，其中萬古霉素類

vanC

-03抗性基因占總萬古霉素類抗性基因相對豐度的70.41%。萬古霉素作為人類抵抗細菌感染的“最后一道防線”之一，然而，隨著人類醫療中萬古霉素的廣泛使用，已在很多環境中富集，如河口沉積物、土壤、水體，萬古霉素類抗性基因的檢出和富集表明牛糞肥可能受到人為污染源的影響，并警示糞肥中抗性基因的污染亟需管控。此外，

qacEdelta

1-02抗性基因僅在羊糞肥和雞羊混合糞肥中檢出，且分別占兩者總抗性基因相對豐度的68.11%和22.76%，且該抗性基因與可移動遺傳元件呈顯著正相關(

R

=0

.

993，

p

<0.01)。

qacEdelta

1抗性基因一般存在于Ⅰ類整合子的3′保守區，該抗性基因在糞肥中的富集表明這些糞肥施用至土壤環境中會存在潛在的水平基因轉移。

3 結論與建議

(1) 4種網絡交易的糞肥中存在豐富多樣的抗生素抗性基因，糞肥的網絡交易可能是抗生素抗性基因潛在的遷移和傳播載體。

(2) 購買于不同商家的4種糞肥中抗生素抗性基因的豐度均有顯著的差異(

p

<0.01)，且其主導的抗性基因類型以及亞基型組成均存在差異，這可能與生產糞肥的養殖場類型、抗生素使用類型以及發酵技術相關。更大尺度和廣度的網絡源糞肥中抗性基因的分布規律仍需進一步探究。(3) 4種網絡交易的糞肥中均檢測出較高豐度的可移動遺傳元件，且抗生素抗性基因相對拷貝數與可移動遺傳元件的相對拷貝數呈極顯著正相關(

p

<0.01)，表明糞肥中的抗生素抗性基因存在轉移的風險，因此需加強對糞肥抗生素抗性基因水平轉移的風險評估，并建議國家相關農業部門加強糞肥的監查管理，制定網絡交易糞肥的微生物安全標準，并完善相關法律管控。