毒死蜱降解菌的馴化篩選及高活性毒死蜱降解菌群的構建

孫孝文 劉誠 王慧敏

摘要 [目的]通過對3種活性污泥的長期馴化、微生物多樣性分析及優(yōu)勢菌篩選與復配，構建一種對毒死蜱具有高降解活性的菌群。 [方法]以農藥公司的3種不同廢水處理工序的污泥為材料，利用毒死蜱長期脅迫馴化以獲得高效降解毒死蜱的活性污泥;隨后，在馴化后的3種活性污泥中篩選優(yōu)勢菌株，并以篩選的菌株構建復合菌群用于毒死蜱的降解。[結果]以毒死蜱為碳源，在馴化后的活性污泥中篩選到3株具有毒死蜱降解活性的菌株（Paracoccus MLCa-1、Klebsiella MLCp-2、Serratia MLCs-1），并對3種菌株進行復配，構建了復合菌群MLCF7。復合菌群MLCF7對100 mg/L毒死蜱降解率達82.84%。 [結論]復合菌群對毒死蜱具有高降解活性，表明其在農藥污染生物修復中具有應用潛力。

關鍵詞 毒死蜱;活性污泥;脅迫馴化;復合菌群;生物降解

中圖分類號 X-172? 文獻標識碼 A? 文章編號 0517-6611（2021）13-0064-05

doi：10.3969/j.issn.0517-6611.2021.13.017

開放科學（資源服務）標識碼（OSID）：

Domestication and Screening of Chlorpyrifos degrading Bacteria and Construction of Composited Flora with High Chlorpyrifos degrading Capacity

SUN Xiao wen1，2，LIU Cheng2，WANG Hui min1

（1. Hubei Key Laboratory of Environmental and Health Effects of Persistent Toxic Substances， Institute of Environment and Health， Jianghan University， Wuhan，Hubei 430056;2. School of Life Sciences， Hubei University，Wuhan，Hubei 430062）

Abstract [Objective] Based on three activated sludges， through long term domestication， analysis of microbial diversity， and dominant bacteria screening and composition to construct a composited flora with high chlorpyrifos degrading capacity. [Method]Using sludge from 3 different wastewater treatment processes of pesticide companies as materials，the chlorpyrifos long term stress domestication was used to obtain activated sludge that efficiently degrades chlorpyrifos;subsequently， the dominant strains were screened from the three domesticated activated sludges， and the screened strains were used to construct a complex flora for the degradation of chlorpyrifos.[Result] Using chlorpyrifos as a carbon source， three dominant bacterial strains （Paracoccus MLCa 1， Klebsiella MLCp 2， and Serratia MLCs 1） were screened from these domesticated activated sludges when chlorpyrifos as the carbon source. A chlorpyrifos degrading flora was constructed by composition of these three dominant bacterial strains. The degradation capacity of the composited flora MLCF7 was reached 82.84% towards 100 mg/L chlorpyrifos. [Conclusion] The easily constructed composited flora MLCF7 has high degradation activities to chlorpyrifos， indicating it has the application potential for bioremediation of pesticide contamination.

Key words Chlorpyrifos;Activated sludge;Stress domestication;Composited flora;Biodegradation

毒死蜱是一種有機磷廣譜殺蟲劑，廣泛應用于經(jīng)濟作物的病蟲害防治[1-2]。作為一種乙酰膽堿酯酶（AChE）的抑制劑，毒死蜱可以結合AChE并使其失去分解乙酰膽堿的能力，導致大量乙酰膽堿積累，從而使中樞神經(jīng)系統(tǒng)功能喪失，最終導致昆蟲和哺乳動物死亡[3-4]。毒死蜱在農業(yè)上的大量持續(xù)使用，導致其在土壤和水中過量累積，從而對環(huán)境造成了嚴重的污染[1，5]。目前，已報道多種處理毒死蜱污染的方法，相較于傳統(tǒng)的物理化學方法（紫外光解、fenton 氧化法、臭氧氧化法和多相光催化降解等[6-9]），微生物修復法具有高效、廉價、無二次污染等特點，被認為是安全有效的毒死蜱環(huán)境污染修復方法[10]?，F(xiàn)已發(fā)現(xiàn)了一些來源于毒死蜱污染的土壤或者水體的毒死蜱降解菌株，包括Enterobacter、Pseudomonas、Serratia、Klebsiella、Bacillus、Paracoccus、Sphingopyxis、Ochrobactrum、Shewanella和Tistrella等[11-17]。相比單一微生物，復合菌群由于包含更健全的毒死蜱降解代謝基因，可有效避免降解途徑中中間產(chǎn)物積累且對環(huán)境具有更強的適應性，更適合毒死蜱環(huán)境污染的修復。Vidya等[18]報道了來源于農藥污染的土壤菌群（包含P.aeruginosa、B.cereus、Klebsiella sp.、S.marscecens）可以在30 d完全降解50 mg/L。Barathidasan等[19]報道了一種由Cellulomonas fimi 和Phanerochaete chrysosporium組成的菌群可以完全礦化毒死蜱，且其效果強于單種菌株。因此，開發(fā)新的毒死蜱降解菌群，將對毒死蜱農藥污染修復有著積極意義。該研究使用來源于農藥公司廢水處理工序的3種污泥，利用毒死蜱長期脅迫馴化以獲得高效降解毒死蜱的活性污泥;隨后，在馴化后的3種活性污泥中篩選優(yōu)勢菌株，并以篩選的菌株構建了復合菌群用于毒死蜱的降解。

1 材料與方法

1.1 試驗材料

1.1.1 活性污泥。3種活性污泥分別采集于湖北省農藥公司中農藥制備污水處理的不同工序。活性污泥MLCa采集于好氧池I、Ⅱ、Ⅲ、Ⅳ和V（1∶1∶1∶1∶1比例混合）;活性污泥MLCp采集于推流曝氣反應池;活性污泥MLCs采集于SBR反應器。

1.1.2 主要試劑。毒死蜱標準品（純度≥99.0%）、3，5，6-三氯-2-吡啶醇（TCP）（純度≥99.0%）購于Aladdin（中國，上海）。色譜級乙腈、二氯甲烷（純度≥99.8%）購于MREDA（中國，北京）。

1.1.3 培養(yǎng)基。LB培養(yǎng)基：胰蛋白胨10.0 g、酵母提取物5.0 g、NaCl 10.0 g，pH 7.0～7.2。MSM培養(yǎng)基：Na2HPO4 5.8 g、KH2PO4 3.0 g、NaCl 0.5 g、NH4Cl 1.0 g、MgSO4·7H2O 0.25 g，pH 7.0～7.2。馴化補料培養(yǎng)基：葡萄糖10.0 g、（NH4）2SO4 0.5 g、MgSO4·7H2O 0.3 g、NaCl 0.3 g、KCl 0.3 g、FeSO4·7H2O 0.03 g、MnSO4·4H2O 0.03 g、酵母粉0.4 g，pH 7.0～7.2。

1.2 毒死蜱脅迫下活性污泥的馴化

將新鮮采集的3種活性污泥分別按30%的沉降比接種到含有90 mL無菌水的500 mL錐形瓶中，并補加50 mL滅菌的馴化補料培養(yǎng)基，使錐形瓶中液體的總體積為200 mL。28 ℃、150 r/min條件下進行馴化7 d后，靜置過夜，倒出50 mL馴化體系上清，補加50 mL 新鮮的馴化補料培養(yǎng)基，并添加5 mg毒死蜱脅迫進行下一輪馴化。此后，后一輪馴化體系比前一輪毒死蜱添加量遞增5 mg。循環(huán)馴化20次后，維持每輪100 mg毒死蜱添加量不變，繼續(xù)馴化4次，結束馴化。

1.3 活性污泥微生物多樣性檢測

1.3.1 活性污泥取樣。分別取3種活性污泥不同馴化時期（初始樣品，第10輪馴化樣品和最終馴化樣品）的活性污泥樣品各5 mL，8 000 r/min離心10 min去上清，活性污泥沉淀送至派森諾生物科技股份有限公司（中國，上海）進行微生物多樣性分析。每個樣品取3組平行樣。

1.3.2

活性污泥微生物多樣性檢測及分析。微生物多樣性測序平臺為Illumina MiSeq。使用QIIME v1.8.0[20]對獲得的高通量測序數(shù)據(jù)剔除不可靠的序列。對得到的序列進行合并，按97%的序列相似度對操作分類單元（OUT）進行劃分;選取每個OTU中最豐富的序列作為OTU的代表序列。然后利用每個OUT的代表序列進行分類狀態(tài)識別和系統(tǒng)發(fā)育分析。使用Chao1和Shannon指數(shù)考察微生物群落的Alpha多樣性。使用R軟件，對豐度前50位的屬進行聚類分析并繪制熱圖。使用Metastats[21]分析樣品之間在屬水平上豐度的差異（P值<0.05具有統(tǒng)計學差異）。采用基于加權UniFrac距離的主坐標分析方法（PCoA）對群落相似性進行多維分析[22]。

1.4 毒死蜱檢測

1.4.1 HPLC法毒死蜱濃度標準曲線繪制。色譜級乙腈配制0、50、100、150、200和250 mg/L的毒死蜱溶液，使用高效液相色譜儀LC-20A（SHIMADZU，日本）測定毒死蜱的濃度。色譜柱為Hypersil ODS-C18 反向色譜柱（4.6 mm×150 mm×5 μm）。流動相為乙腈∶水=70∶30，流速1.0 mL/min，柱溫40 ℃，進樣20 μL。檢測波長290 nm，毒死蜱的保留時間為7.8 min，峰面積和毒死蜱濃度呈正相關，回歸方程為峰面積= 9 951×毒死蜱濃度-22 414（R2=0.998）。

1.4.2 樣品中毒死蜱含量的檢測。取含毒死蜱溶液5.0 mL，加入2 倍體積的二氯甲烷萃取2次，回收下層有機相溶液，經(jīng)無水Na2SO4 脫水后，在烘箱中吹干，并重溶于5.0 mL的乙腈中。預試驗顯示，萃取的回收率達90%。使用上述液相條件檢測，并根據(jù)毒死蜱標準曲線計算樣品中毒死蜱濃度。

1.5 毒死蜱降解菌的篩選與鑒定

1.5.1 毒死蜱降解菌的初篩。 將馴化后的活性污泥在初篩培養(yǎng)基I（MSM培養(yǎng)基添加100 mg/L毒死蜱）上進行劃線分離;挑取生長良好的單菌落過夜活化，然后再在初篩培養(yǎng)基Ⅱ（MSM培養(yǎng)基添加200 mg/L毒死蜱）上進行劃線分離;再將生長良好的單菌落過夜活化，最后在初篩培養(yǎng)基Ⅲ（MSM培養(yǎng)基添加500 mg/L毒死蜱）上進行劃線分離，最終得到能夠以500 mg/L毒死蜱為碳源生長的菌株。

1.5.2 毒死蜱降解菌的復篩。將上述分離的菌株過夜活化，按1%接種量接種于含100 mg/L 毒死蜱的MSM 液體培養(yǎng)基中，28 ℃、150 r/min培養(yǎng)7 d，8 000 r/min離心10 min檢測上清液殘留毒死蜱濃度，并計算毒死蜱的降解率：

降解率=（1-Ct/C0）×100%，

式中，Ct為農藥殘留濃度，C0為農藥初始濃度。

1.5.3 毒死蜱降解菌的鑒定。挑取毒死蜱降解菌的單菌落，使用通用引物27F（5′-AGAGTTTGATCCTGGCTCAG-3′）和1492R（5′-TACGGCTACCTTGTTACGACTT-3′），PCR擴增16S rRNA基因序列。PCR 產(chǎn)物回收后與載體pClone007 Blunt Vector 連接，轉化感受態(tài)細胞E.coli DH5α后送至武漢擎科生物技術公司進行測序。將測序得到的結果在EzTaxon server 2.1 數(shù)據(jù)庫[23]中進行比對后，采用CLUSTAL W對16S rRNA 基因序列進行多序列比對，并利用MEGA 6.0 軟件[24]neighbor-joining method構建系統(tǒng)發(fā)育樹分析，完成菌株鑒定。

1.6 復合菌群的構建及其降解毒死蜱的效率測定

1.6.1 復合菌群的構建。復合菌群由3 株經(jīng)過交叉平板劃線培養(yǎng)后沒有拮抗活性的毒死蜱降解菌株組成。將過夜活化的3 株單菌（MLCa-1、MLCp-2 和MLCs-1），按照1∶1∶1 的體積比比例混合，并按1%的接種量將混合菌液接種到新鮮的LB 培養(yǎng)基中。28 ℃、150 r/min培養(yǎng)24 h，即得到復合菌群（MLCF7）。

1.6.2 復合菌群降解毒死蜱的效率測定。將OD600 nm=1.0 的MLCF7按2%的接種量接種到含有100 mg/L毒死蜱MSM 培養(yǎng)基中，28 ℃、150 r/min培養(yǎng)7 d。采用上述HPLC法檢測毒死蜱的降解效率。

2 結果與分析

2.1 馴化后的活性污泥對毒死蜱的降解效果

3種來源于農藥公司污水處理系統(tǒng)不同處理工序的活性污泥，在毒死蜱脅迫下，經(jīng)過6個月的馴化后，回收了降解活性穩(wěn)定的活性污泥。馴化后的活性污泥，以1%的添加量添加至含100 mg/L毒死蜱的MSM培養(yǎng)基中，在28 ℃、150 r/min的條件下，培養(yǎng)5 d，檢測毒死蜱降解效率。由圖1可見，隨著培養(yǎng)時間的增加，3種活性污泥對毒死蜱的降解效率也不斷增強。3種活性污泥對毒死蜱的第5天降解效率排序為活性污泥MLCp（100.00%）> 活性污泥MLCs（87.59%）> 活性污泥MLCa（82.34%）。3種馴化后的活性污泥對毒死蜱的降解率均達到了80%以上，說明經(jīng)過毒死蜱的長期脅迫，活性污泥被成功馴化。

2.2 3種活性污泥在不同馴化時期的微生物多樣性分析

2.2.1 樣本微生物群落的Alpha多樣性。對3種活性污泥在不同馴化時期進行了微生物多樣性分析。樣品包含初始樣品（MLCa-0、MLCp-0、MLCs-0）、第10輪馴化樣品（MLCa-2、MLCp-2、MLCs-2）和最終馴化樣品（MLCa-4、MLCp-4、MLCs-4）。每個樣本重復3次，使用A、B和C區(qū)別。

使用Chao1和Shannon指數(shù)評價樣本微生物群落的Alpha多樣性。Chao1指數(shù)值越大，表明群落的豐富度越高。Shannon指數(shù)值越高，表明群落的多樣性越高。由表1可知，3種活性污泥在毒死蜱的長期脅迫下的馴化過程中，微生物群落豐富度和多樣性呈現(xiàn)了不同的變化趨勢。來源于好氧池的活性污泥MLCa在馴化過程中微生物群落豐富度和多樣性均不斷下降;來源于推流曝氣反應池的活性污泥MLCp在馴化過程中微生物群落豐富度不斷下降，多樣性先上升后下降;來源于SBR反應器的活性污泥MLCs在馴化過程中微生物群落豐富度和多樣性均不斷上升。這可能與3種活性污泥原始微生物菌群組成差異有關。

2.2.2 基于UniFrac距離的PCoA主坐標分析。?? 由圖2可知，每個點代表一個樣本，不同顏色的點屬于不同樣本（組），兩點之間的距離越近，表明2個樣本之間的微生物群落結構相似度越高，差異越小。坐標軸括號中的百分比代表了對應的主坐標所能解釋的原始數(shù)據(jù)中差異的比例。結果表明3種活性污泥在馴化的過程中，微生物群落結構都發(fā)生了明顯的變化。菌群MLCp和MLCs馴化后微生物群落結構較為相近，與MLCa馴化后微生物群落結構差異較大。

2.2.3

結合聚類分析的群落組成熱圖。使用R 軟件，對樣本中豐度前50 位的屬進行聚類分析并繪制熱圖。由圖3可知，馴化結束后，活性污泥MLCa 中主要優(yōu)勢菌群有Paracoccus sp.、Enterobacter sp.和Methylobacillus sp.等;活性污泥MLCp 中主要優(yōu)勢菌群有Klebsiella sp.、Pseudoxanthomonas sp.和Achromobacter sp.等;活性污泥MLCs 中主要優(yōu)勢菌群有Serratia sp.、Bacillus sp.和Pseudomonas sp.等。

2.3 毒死蜱降解菌的篩選及鑒定

2.3.1 毒死蜱降解菌的篩選。將3種馴化后的活性污泥在含毒死蜱的MSM固體培養(yǎng)基中經(jīng)過多次劃線分離篩選，共得到7株可以以500 mg/L毒死蜱為碳源生長的菌株，分別為來源于活性污泥MLCa的MLCa-1和MLCa-2，來源于活性污泥MLCp的MLCp-1、MLCp-2和MLCp-3以及來源于活性污泥MLCs的MLCs-1和MLCs-2。所有毒死蜱降解菌株經(jīng)過過夜活化后，按2%的接種量接種于含有100 mg/L 毒死蜱的MSM 培養(yǎng)基中，28 ℃、150 r/min培養(yǎng)7 d后，測定了菌株對毒死蜱的降解率。由圖4可知，所有毒死蜱降解菌株對毒死蜱均有降解作用，其中，降解率前三的菌株為MLCa-1（降解率為53.71%）、MLCp-2（降解率為51.85%）和MLCs-1（降解率為65.37%）。

2.3.2 毒死蜱降解菌的鑒定。篩選的毒死蜱降解菌株的16S rRNA基因，經(jīng)過CLUSTAL W多序列比對，MEGA 6.0構建系統(tǒng)發(fā)育樹對菌株進行鑒定。由圖5可知，MLCa-1與Paracoccus bengalensis JJJT的序列相似性為99.40%，MLCp-2與Klebsiella singaporensis LX3T的序列相似性為99.37%，MLCs-1與Serratia nematodiphila DSM 21420T的序列相似性為99.58%。因此將MLCa-1鑒定為Paracoccus sp.，MLCp-2鑒定為Klebsiella sp.，MLCs-1鑒定為Serratia sp.。

2.4 復合菌群MLCF7的構建及其對毒死蜱的降解

2.4.1 MLCF7降解毒死蜱的生長曲線和降解曲線測定。復合菌群MLCF7按1%的接種量接種于含100 mg/L毒死蜱的MSM培養(yǎng)基中，150 r/min培養(yǎng)7 d，監(jiān)測毒死蜱的降解率。由圖6可知，在7 d的降解過程中，降解體系的菌液OD600 nm不斷上升，第5天后達到穩(wěn)定期。相應地，在7 d的降解過程中，降解體系降解曲線也呈現(xiàn)不斷的上升趨勢。這表明復合菌群MLCF7可以降解毒死蜱，并利用毒死蜱為碳源生長。在第7天時，復合菌群MLCF7對毒死蜱的降解率達82.84%，由圖4可知，3株單菌對毒死蜱的降解率都在50%～60%。復合菌群對毒死蜱的降解顯著性高于3株單菌對毒死蜱的降解率。結果表明，3株單菌存在協(xié)同作用，提高了對毒死蜱的降解能力。

2.4.2 復合菌群MLCF7降解毒死蜱中間代謝產(chǎn)物研究。使用HPLC法檢測第1天、第3天、第7天復合菌群MLCF7降解毒死蜱體系代謝產(chǎn)物，由圖7可知，隨著降解時間增長，毒死蜱濃度不斷下降，中間產(chǎn)物TCP的濃度先上升后下降，直至消失。因此，復合菌群MLCF7 首先將毒死蜱降解成TCP，再進一步地TCP 被降解成其他的小分子物質。

3 結論與討論

對3種來源于農藥公司廢水處理系統(tǒng)的污泥進行毒死蜱長期脅迫下的馴化，得到了3種對100 mg/L毒死蜱降解率達80%以上的活性污泥。對3種活性污泥馴化的過程進行微生物多樣性分析，發(fā)現(xiàn)毒死蜱脅迫下3種活性污泥的菌群發(fā)生了明顯的變化。其中，副球菌屬、克雷伯氏菌屬和沙雷式菌屬等菌屬在馴化的過程中被富集。通過以毒死蜱為碳源，在馴化后的活性污泥種篩選了3株毒死蜱降解細菌Paracoccus MLCa-1、Klebsiella MLCp-2和Serratia MLCs-1。以3種毒死蜱降解細菌構建了復合菌群MLCF7，對100 mg/L毒死蜱降解率達82.84%。

3種毒死蜱降解菌對毒死蜱的降解活性較低，但復合后降解活性明顯增加，這種協(xié)同機制需進一步研究。復合菌群MLCF7降解毒死蜱會形成中間產(chǎn)物TCP，然而TCP如何繼續(xù)被降解，也需進一步研究。復合菌群MLCF7構建簡單，易循環(huán)構建使用，可避免活性污泥降解易退化的問題，在農藥污染生物修復中具有應用潛力。

參考文獻

[1] JOHN E M，SHAIKE J M.Chlorpyrifos：Pollution and remediation[J].Environ Chem Lett，2015，13（3）：269-291.

[2] 李亞楠，劉晶晶，陳少華，等.毒死蜱的應用現(xiàn)狀及降解研究進展[J].廣東農業(yè)科學，2011，38（6）：92-96.

[3] 張葉翠，李翎，胡晨陽，等.毒死蜱生殖毒性與神經(jīng)毒性研究進展[J].中國職業(yè)醫(yī)學，2019，46（5）：628-632.

[4] EATON D L，DAROFF R B，AUTRUP H，et al.Review of the toxicology of chlorpyrifos with an emphasis on human exposure and neurodevelopment[J].Crit Rev Toxicol，2008，38（S2）：1-125.

[5] 雷文娟，周向陽.生物炭對農藥降解產(chǎn)物三氯吡啶醇在土壤中遷移的影響研究[J].農業(yè)工程學報，2019，35（10）：173-180.

[6] MERCI S，SALJOOQI A，SHAMSPUR T，et al.Investigation of photocatalytic chlorpyrifos degradation by a new silica mesoporous material immobilized by WS2 and Fe3O4 nanoparticles：Application of response surface methodology[J].Appl Organomet Chem，2020，34（3）：e5343.

[7] MU OZ A，VERA T，SIDEBOTTOM H，et al.Studies on the atmospheric degradation of chlorpyrifos methyl[J].Environ Sci Technol，2011，45（5）：1880-1886.

[8] RIBEIRO A R，NUNES O C，PEREIRA M F，et al.An overview on the advanced oxidation processes applied for the treatment of water pollutants defined in the recently launched Directive 2013/39/EU[J].Environ Int，2015，75：33-51.

[9] CAO X，YAN C，YANG X，et al.Photolysis induced neurotoxicity enhancement of chlorpyrifos in aquatic system：A case investigation on Caenorhabditis elegans[J].J Agric Food Chem，2020，68（2）：461-470.

[10] 黃強，吳祥為，花日茂，等.農藥毒死蜱環(huán)境污染微生物修復研究進展[J].安徽農業(yè)科學，2008，36（22）：9682-9683.

[11] 陳琳，祁靜，李祖明，等.白菜葉際細菌多樣性與毒死蜱降解菌篩選及分離鑒定[J].食品工業(yè)科技，2018，39（22）：107-112，120.

[12] 胡嘉博，何繼烈.毒死蜱降解菌DH6的分離鑒定及降解機理研究[J].科學與信息化，2018（35）：188-191.

[13] NAYAK T，PANDA A N，KUMARI K，et al.Comparative genomics of a paddy field bacterial isolate Ochrobactrum sp.Cpd 03：Analysis of chlorpyrifos degradation potential[J].Indian J Microbiol，2020，60（3）：325-333.

[14] AHIR U N，VYAS T K，GANDHI K D，et al.In vitro efficacy for chlorpyrifos degradation by novel isolate Tistrella sp.AUC10 isolated from chlorpyrifos contaminated field[J].Curr Microbiol，2020，77（9）：2226-2232.

[15] GHANEM I，ORFI M，SHAMMA M.Biodegradation of chlorphyrifos by Klebsiella sp.isolated from an activated sludge sample of waste water treatment plant in Damascus[J].Folia Microbiol，2007，52（4）：423-427.

[16] GOVARTHANAN M，AMEEN F，KAMALA KANNAN S，et al.Rapid biodegradation of chlorpyrifos by plant growth promoting psychrophilic Shewanella sp.BT05：An eco friendly approach to clean up pesticide-contaminated environment[J].Chemosphere，2020，247：1-7.

[17] ZHU J W，ZHAO Y，RUAN H H.Comparative study on the biodegradation of chlorpyrifos methyl by Bacillus megaterium CM Z19 and Pseudomonas syringae CM Z6[J].An Acad Bras Cienc，2019，91（3）：e20180694.

[18] VIDYA LAKSHMI C，KUMAR M，KHANNA S.Biodegradation of chlorpyrifos in soil by enriched cultures[J].Curr Microbiol，2009，58（1）：35-38.

[19] BARATHIDASAN K，REETHA D，JOHN MILTON D，et al.Biodegradation of chlorpyrifos by co culture of Cellulomonas fimi and Phanerochaete chrysosporium[J].Afr J Microbiol Res，2014，8（9）：961-966.

[20] CAPORASO J G，KUCZYNSKI J，STOMBAUGH J，et al.QIIME allows analysis of high throughput community sequencing data[J].Nat Methods，2010，7（5）：335-336.

[21] WHITE J R，NAGARAJAN N，POP M.Statistical methods for detecting differentially abundant features in clinical metagenomic samples[J].PLoS Comput Biol，2009，5（4）：1-11.

[22] RAMETTE A.Multivariate analyses in microbial ecology[J].FEMS Microbiol Ecol，2007，62（2）：142-160.

[23] YOON S H，HA S M，KWON S，et al.Introducing EzBioCloud：A taxonomically united database of 16S rRNA gene sequences and whole genome assemblies[J].Int J Syst Evol Microbiol，2017，67（5）：1613-1617.

[24] TAMURA K，STECHER G，PETERSON D，et al.MEGA6：Molecular evolutionary genetics analysis version 6.0[J].Mol Biol Evol，2013，30（12）：2725-2729.