江漢油田區典型農田土壤烴類降解微生物功能基因bssA的遺傳多樣性研究*

左小虎 王明霞 姚炎紅 李振輪 周志峰

（西南大學資源環境學院，重慶 400716）

江漢油田區典型農田土壤烴類降解微生物功能基因bssA的遺傳多樣性研究*

左小虎 王明霞 姚炎紅 李振輪 周志峰?

（西南大學資源環境學院，重慶 400716）

油田區土壤易受烴類物質影響并可能富集了特異的石油烴降解微生物類群。針對江漢油田區5個不同油井口附近的典型旱地農田土壤，采用石油烴（Petroleum hydrocarbons，PHs）中苯系物代謝的關鍵功能基因-bssA（苯甲基琥珀酸合成酶基因）作為分子標識物，通過克隆文庫結合末端限制性片段長度多樣性（Terminal-restriction fragment length polymorphism，T-RFLP）的方法，研究該油田區土壤含有bssA基因的烴類降解微生物群落結構，并探討其環境驅動機制。結果表明，土壤中PAHs含量在0.21～2.01 mg kg-1之間，石油烴污染程度較低。T-RFLP的分析表明不同土壤樣品中的bssA基因多樣性差異明顯，PAHs（Polycyclic Aromatic Hydrocarbons，多環芳烴）含量最高土壤中bssA基因多樣性最高，其優勢bssA基因類群與硫酸鹽還原菌或地桿菌有較近的親緣關系。冗余分析進一步表明，土壤硝態氮、有效磷、PAHs含量均是影響bssA基因多樣性的重要因子。這些結果表明：江漢油田區典型農田土壤中含有bssA基因的主要類群為β-變形菌和δ-變形菌，并與地桿菌屬（Geobacter）、索氏菌屬（Thauera）和固氮菌屬（Azoarcus）具有較近的系統發育親緣關系。這些微生物可能通過硝酸鹽、硫酸鹽及鐵還原代謝過程降解土壤PAHs。

土壤；厭氧降解；苯甲基琥珀酸合成酶基因（bssA）；群落結構

石油烴（Petroleum hydrocarbons，PHs）污染是油田開發面臨的重要問題，尤其是油田區域的農田土壤，其污染物微生物修復得到了廣泛關注。石油烴是由烷烴、烯烴、芳香烴、雜環芳烴等多組分的復雜分子組成的均質混合物，其中的多環芳烴（Polycyclic Aromatic Hydrocarbons，PAHs）和苯系物（苯、甲苯、乙苯、二甲苯等）備受關注，是土壤中常見的有機污染物［1-3］，不僅嚴重危害土壤環境質量，也會直接或間接危害人類健康［4-5］。因此，土壤中PAHs和苯系物降解微生物一直是學術界關注的熱點［6-7］。一般認為，在純培養體系有氧條件下PAHs及苯系物的微生物降解更易發生，但在復雜土壤和沉積物等自然環境中，石油烴大多閉蓄于微孔介質等厭氧或微氧生境，厭氧微生物降解可能發揮了重要的作用［3，8］。例如，在苯系物、烴類及萘等低分子量PAHs的厭氧降解過程中，由甘氨酰自由基合成酶所催化的延胡索酸加入途徑起著十分重要的紐帶作用［9］。該途徑最早發現于β-變形菌厭氧降解甲苯的過程中，該過程中的關鍵環節是由苯甲基琥珀酸合成酶（Benzylsuccinate synthase，BSS）基因所編碼的BBS催化延胡索酸加入到甲苯［10］，形成的（R）-苯甲基琥珀酸酯［11-13］可經隨后的一系列反應最終轉化為關鍵中間代謝物苯甲酰-CoA［9］。隨后的研究表明，該反應普遍存在于苯系物，烷烴及PAHs的厭氧降解過程中。如，硫酸鹽還原菌（Sulfate-reducing bacteria，SRB）和硝酸鹽還原菌（Nitrate-reducing bacteria，NRB）在厭氧條件下也可以對n-烷烴進行類似于甲苯添加延胡索酸的活化降解，該過程由烷基琥珀酸合成酶（ASS）催化完成［14-15］，而SRB則是通過萘基-2-甲基-琥珀酸合成酶（NMS）催化延胡索酸加入2-甲基萘［16-17］，來完成厭氧代謝2-甲基萘的第一步反應。因此，近年來，編碼BSS，ASS及NMS的關鍵功能基因（bssA；assA；nmsA）被嘗試作為分子標識物，研究環境中苯系物等PHs類化合物厭氧代謝的微生物類群［10，18-19］，表征環境受污染程度及微生物降解的潛在能力。如，Frederick等［20］以bssA、nmsA和assA作為特異性標識物，研究了一些石油烴污染環境（陸地及海洋）中的厭氧PHs降解菌群，并發現了一些與16S rRNA基因分類迥然不同的δ-變形菌和梭菌類群。Alejandro 等［21］比較了bssA、nmsA和assA基因的遺傳多樣性，研究了原油泄漏區域海水及沉積物中PHs厭氧降解菌的群落結構，發現這些基因的豐度及多樣性具有明顯的差異，并與受污染的程度密切相關。但已有的研究均表明，與nmsA和assA分子標靶基因相比較，以bssA為代表的編碼延胡索酸加入的關鍵基因具備較強的專一性，是未來研究石油污染環境微生物修復的重要標識物。但目前相關的bssA基因研究多見于海洋石油污染環境，油田區土壤中可能富集了大量的石油烴類降解微生物［7，22］，研究其中的關鍵功能基因多樣性對污染環境生物修復具有重要意義。

油田區域的土壤有受原油污染的風險，而土壤中極可能存在著厭氧/微氧環境［23］。據此，本研究針對具有50年油氣開發歷史的江漢油田區，在約20 km2的范圍內，選擇了5個不同油井附近的典型旱地農田土壤，以bssA基因為分子標識物，通過克隆文庫結合T-RFLP的方法研究了油田區典型農田土壤中苯系物厭氧降解微生物的群落結構，初步探討了其群落結構與環境因子之間的關系，為進一步深入研究PHs類有機污染物的土壤厭氧代謝的微生物機理提供了一定的理論依據。

1 材料與方法

1.1 土壤樣品采集

江漢油田位于湖北省潛江市境內，該油田已有50多年的油氣開發歷史，是新中國最早開發的油氣田之一。本研究的土壤樣品采集于2014年7月。在江漢油田區域共選擇5個不同的油井，編號為JH-1 至JH-5，每個油井口30～40 m半徑范圍內，采集農田0～20 cm表層新鮮土壤（潮土，黏壤）。在每個采樣點中心周圍30 m2的范圍內，隨機選擇9個點，進行多點混合采樣并初步剔除石塊和雜物，以完成單個土樣的采集，且各樣點內均按上述方法采集3個重復的土壤樣品。采集好的樣品裝于自封袋后，置于放有冰袋的箱子中運回實驗室。隨后，過2 mm篩后分裝，部分放于4 ℃用于基本性質測定，部分置于-20 ℃冰箱用于DNA的提取。

1.2 土壤性質及PAHs含量測定

土壤pH采用2.5∶1的水土比，用“PHS-3C型”pH儀測定。土壤總有機質、有效磷、銨態氮和硝態氮的測定參照《土壤農業化學分析方法》［24］進行，分別使用重鉻酸鉀氧化還原滴定法（K2Cr2O7-H2SO4法）、鉬銻抗比色法、靛酚藍分光光度法和萘乙二胺分光光度法測定。

土壤中PAHs含量的測定采用經典的索氏提取法［25］。具體步驟如下：稱取10 g經冷凍干燥的土壤樣品、10 g無水硫酸鈉及少量銅粉置于同一燒杯中混勻，無損移入濾紙桶，加入150 ml丙酮–正己烷（1∶1）混合液浸泡土壤樣品12 h，將濾紙桶裝入索氏提取器中75 ℃提取6 h；提取液于-46 kPa，45 ℃下旋轉蒸發濃縮至1～2 ml，然后用硅膠層析柱（1 cm無水硫酸鈉，6 cm活化硅膠和3 cm氧化鋁）凈化；以15 ml色譜級正己烷預洗后，用70 ml二氯甲烷和正己烷的混合液（3∶7，V/V）洗脫；收集洗脫液，洗脫液旋轉蒸發至1～2 ml，再加入25 ml甲醇進行溶劑置換，再旋轉蒸發濃縮并用甲醇定容至1 ml；最后用高效液相色譜儀（HITACHI L-7100，日本）測定各樣品中PAHs含量。

1.3 土壤DNA提取

土壤DNA提取：使用PowerSoil?DNA Isolation Kit（MO BIO，美國）試劑盒提取土壤DNA。稱取0.5 g土壤樣品，具體操作步驟按說明書進行。提取到的總DNA，用1%的瓊脂糖凝膠電泳檢測提取的DNA產物，并用DU 800 spectrophotometer （Beckman Coulter Inc，美國）測定其濃度，最后置于-20 ℃保存備用。

1.4 克隆文庫構建及T-RFLP分析

為了全面解析所采土壤樣品中bssA、nmsA和assA基因的類群，將所有5個土壤樣品的bssA、nmsA和assA基因PCR擴增產物分別混合并構建其混合克隆文庫，為后續T-RFLP分析中限制性內切酶的選擇提供依據，同時也有助于確定T-RFLP分析中所獲取的各末端片段所對應的類群。具體步驟如下：分別以5個不同土壤樣品的DNA為模板，利用表1中的引物及相應的擴增條件對上述各目標基因進行擴增，聚合酶鏈式反應主要包括：0.25 μl TaKaRa Ex Taq（5 U μl-1）、2.5 μl 10×Ex Taq Buffer、2 μl dNTP Mixture（2.5 m mol-1）、前后引物各0.5 μl（10 μm mol-1）、2 μl DNA模板，加入滅菌水補足至25 μl反應體系。將不同土壤樣品bssA、nmsA及assA基因的PCR產物分別混合，混合后的PCR產物用Wizard SV Gel and PCR Clean-Up System（Promega，美國）進行純化，純化后的PCR產物與pGEM-T Easy Vector（Promega，美國）連接，之后轉入Escherichia coli JM109 （TakaRa，日本）感受態細胞，通過藍白斑篩選，各自隨機挑選陽性克隆子送上海生物工程有限公司進行測序，從而完成各目的基因混合克隆文庫的構建。

表1 PCR擴增引物及反應程序Table 1 Primer sets and PCR procedures

使用DNAMAN（Version 6.0.3.48，Lynnon Biosoft，美國）對所得序列進行同源性分析，相似度＞9 5%的序列歸入同一分類操作單元（Operational Taxonomic Unit，OTU），每個OTU選擇一條具有代表性的序列，采用鄰位相鄰法（Neighbor-Joining），用MEGA（version 6.0）軟件構建系統發育樹。之后，使用DNAMAN軟件對各OTU的代表性序列進行酶切位點分析，結果表明限制性內切酶Alu I（Takara，日本）能夠較好地對bssA基因克隆文庫中優勢的OTU進行分型，可作為后續T-RFLP分析中的限制性內切酶。本研究所獲得的序列已提交至GenBank并獲取了相應的序列號（KX148522-KX148545）。據此，進一步分別以各土壤樣品的DNA為模板，用前引物5′端帶FAM標記的引物對（7772f/8546r）對bssA基因進行PCR擴增。獲得的PCR產物純化后，用Alu I在37 ℃酶切4 h，酶切產物用醋酸鈉和冰乙醇純化［26］，隨后用核酸測序儀ABI PRISM3700（Applied Biosystems，美國）進行基因掃描，用GeneMapper （Applied Biosystems，美國）軟件對T-RFLP圖譜進行分析。

1.5 數據處理

所有數據采用SPSS 17.0進行處理分析，樣品之間的平均值差異采用One–Way ANOVA單因素方差分析，p＜0.05表示差異顯著。用CANOCO （4.5）中的冗余分析（Redundancy analysis，RDA）對bssA基因多樣性與環境因子之間的關系進行分析。

2 結 果

2.1 土壤PAHs含量及理化性質

不同油田附近典型農田土壤PAHs含量及基本性質見表2。各土壤樣品PAHs含量具有明顯差別，其中JH-5（2.01 mg kg-1）和JH-1（1.23 mg kg-1）的含量較高，其余各樣品均低于1 mg kg-1。供試土壤均呈堿性，且各樣品pH之間的差異性不顯著。不同土壤樣品的有機質含量差異性顯著（p ＜0.05），其中樣品JH-5和JH-2的有機質含量較高，分別為40.2 g kg-1和38.3 g kg-1，而JH-9的有機質含量最低，僅為11.3 g kg-1。有效磷含量最高的樣品是JH-5（25.4 mg kg-1），其次是JH-3（15.9 mg kg-1）和JH-4（11.5 mg kg-1），其余各樣品的有效磷含量均低于10 mg kg-1。各土壤樣品的銨態氮含量差異不明顯，均在4.10～4.30 mg kg-1之間。硝態氮含量最高的是JH-5（31.1 mg kg-1），其次是JH-3和JH-4，其余樣品均低于5 mg kg-1。

表2 樣點土壤基本性質Table 2 Basic properties of the tested soil samples

2.2 混合克隆文庫中bssA序列及其末端限制性片段組成

對bssA、nmsA和assA克隆測序結果進行比對分析，結果表明只成功擴增到土壤樣品中的bssA基因。從bssA混合克隆文庫中隨機挑選220個陽性克隆子測序，并在GenBank中進行Blast比對，共得到95條有效的bssA序列。在95%的相似水平下，可分為24個OTUs。各OTU代表性序列經Alu I酶切后的末端限制性片段（Terminal-restriction fragment，T-RF）長度及其所含序列占克隆文庫總序列數的百分比如下表3所示。可見，其序列數占克隆文庫總序列數的比例超過10%的優勢OTUs為OTU-3 （10.5%）、5（11.55%）、6（17.85%）和23 （13.65%）。從表中還可看出，限制性內切酶Alu I能夠較好地將這些優勢OTUs進行分型，適于后續的T-RFLP分析。

從每個OTU中選擇一條代表性的序列來構建bssA基因的系統發育樹（圖1）。從圖1可知，該油田區典型農田土壤bssA基因可分為4個不同的簇（Clusters1-4）。其中，Cluster1含有8個OTUs，且包括占優勢的OTU-3（571bp T-RF）和OTU-6 （409bp T-RF），它們與β-變形菌有較近的親緣關系。而Cluster2和Cluster4中各含有2個和4個OTUs，其所含序列數分別占克隆文庫總序列數的4.2%和11.2%。其余的10個OTUs 均屬于Cluster3，包含占據優勢的OTU-5（388bp T-RF）和OTU-23（433bp T-RF）。此外，Cluster2、 Cluster3 和 Cluster4均與δ-變形菌的親緣關系較近。通過對bssA基因混合克隆文庫的構建及相應的系統發育分析，可為后續T-RFLP中確認不同長度末端限制性片段（T-RFs）所代表的具體bssA類群提供依據。

2.3 土壤bssA基因的遺傳多樣性及其影響因子

各樣品bssA基因的T-RFs組成如圖2所示。由圖可知，樣品JH-1和JH-2中bssA的群落組成較為相似，均含有3種T-RFs類群（96、154、409 和433bp），且409bp和433bp的T-RFs是其共有的優勢類群。結合2.2的結果可知，409bp T-RF可能屬于OTU-1、OTU-6或OTU-18，在系統發育樹中位于Cluster1，與β-變形菌門的紅環菌科（Rhodocyclaceae）下的索氏菌屬（Thauera）和固氮弓菌屬（Azoarcus）有較近的親緣關系。而433bp的T-RFs則屬于OTU-23，位于Cluster3，與δ-變形菌門的地桿菌科（Geobacteraceae）的地桿菌屬（Geobacter）有較近的親緣關系。樣品JH-3中只檢測到一個252bp（不確定類群）獨特的bssA類群。樣品JH-4中的bssA由89bp（Cluster2或4）、175bp（Cluster3）、205bp（不確定類群）和295bp（不確定類群），共4種T-RFs所組成。

而在PAHs含量最高的JH-5中，含有的T-RFs類群最多，即89bp（Cluster4）、96bp（不確定類群）、175bp（Cluster3）、205bp（不確定類群）、295bp（不確定類群）、388bp（Cluster2或3）、409bp（Cluster1）及433bp（Cluster3），共8種T-RFs。其中388bp T-RF（Cluster2或3，硫酸鹽還原菌或地桿菌）的類群是其特異性類群。值得注意的是，在PAHs含量較高的JH-1和JH-5中，盡管bssA基因的組成差異明顯，但是組成JH-1中3種bssA類群（96，409和433bp）均出現在后者中，且96bp（不確定類群）是其所獨有的T-RFs。由此可見，96bp及388bp 的bssA類群可能受到土壤PAHs的富集。將土壤基本性質及T-RFLP圖譜（去除過于單一的JH-3）分別作為環境變量和物種變量，對土壤環境因子與bssA群落結構進行冗余分析的結果（圖3）表明，四個樣品的bssA基因群落組成差異明顯，均不在一個象限。此外，除有效磷、硝態氮、銨態氮等因子外，土壤PAHs含量也是影響bssA群落組成的關鍵因子之一，且96bp（不確定類群）和388bp（硫酸鹽還原菌或地桿菌）T-RFs所代表的bssA類群與PAHs含量呈正相關。

表3 克隆文庫中的OTUs及Alu I 酶切的末端片段Table 3 OTUs in the bssA clone library and Alu I-digested TRFs

3 討 論

圖1 基于bssA克隆文庫構建的N-J（Neighbor-joining）系統發育樹Fig. 1 Phylogenetic tree based on bssA clone library

本研究中，土壤樣品的PAHs含量范圍為0.21～2.01 mg kg-1，低于中華人民共和國環境保護部2008年頒布的《土壤環境質量標準》修訂版（征求意見稿）中擬定的農業用地土壤16種優先控制PAHs總含量的篩選值（4.90～10.1 mg kg-1）［27］。據報道，長江上游重慶金佛山土壤16種優先控制PAHs的含量范圍為0.24～2.12 mg kg-1［28］。而對位于長江下游長江三角洲的蘇州市土壤PAHs含量的調查表明，該區域土壤總PAHs含量范圍為0.045～3.70 mg kg-1［29］。由此可見，本研究各樣點土壤中PAHs含量均在上述調查結果的范圍之內，雖源于油井口附近，卻未遭受嚴重的污染。同時檢測到了苯系物厭氧代謝的關鍵基因-bssA，表明這些農田土壤具有降解石油烴污染物的遺傳基礎。克隆文庫的結果表明，該油田區土壤bssA類群主要為β-變形菌（Cluster1）和δ-變形菌（Cluster2、Cluster3和Cluster4），它們所含序列數分別占克隆文庫總序列數的42.6%和57.4%。已有的研究表明，變形菌門是多數PHs污染土壤中的優勢bssA類群［20-21］，這與本研究的結論一致。其中，β-變形菌門的一些反硝化細菌可以芳香烴為唯一碳源進行生長［30-31］。而δ-變形菌中的一些硫酸鹽還原菌、鐵還原菌等也具備厭氧代謝芳香烴的能力［32］。從屬的水平來看，江漢油田區土壤中β-變形菌的bssA類群與索氏菌屬（Thauera）和固氮弓菌屬（Azoarcus）有較近的親緣關系，它們均具備在反硝化條件下厭氧代謝芳香烴化合物的能力［33］。而該區域土壤中δ-變形菌的bssA類群則與硫酸鹽還原菌TRM1（Sulfate-reducing bacterium TRM1）及地桿菌屬（Geobacter）有較近的親緣關系，上述兩種菌常見于PHs污染土壤中，它們分別能夠以硫酸根及鐵離子為電子受體在厭氧條件下降解芳香烴化合物［34］。如，Winderl等［18］報道，德國慕尼黑煤氣廠附近的水體沉積物中的主要bssA類群為地桿菌屬。Sun等［22］則在分別添加了硝酸鹽和硫酸鹽的甲苯污染土壤中，分別檢測到索氏菌DNT-1（Thauera sp. DNT-1）和硫酸鹽還原菌PRTOL1（sulfate-reducing bacterium PRTOL1）。此外，據報道，δ-變形菌N47和 NaphS2均含有nmsA基因［17，20］，而脫硫腸桿菌屬（Desulfotomaculum）的一些菌株則含有assA基因［19］。本研究中含有bssA基因的類群，目前還沒有發現它們同時含有nmsA或assA基因，這可能是后兩種基因不能成功從該區域土壤中成功擴增的原因。

圖2 土壤樣品中bssA群落結構的組成Fig. 2 Structural composition of the bssA community in the soil samples

圖3 基于bssA基因T-RFLP剖面與土壤性質的RDA分析Fig. 3 RDA analysis based on T-RFLP of bssA and soil properties

本研究中T-RFLP的結果（圖2）表明，該油田區域不同樣點土壤樣品中的bssA微生物群落結構差異明顯。其中，在PAHs含量最高JH-5中，具有最多的bssA基因類群，且388bp（Cluster2或3，硫酸鹽還原菌或地桿菌）是其特異性的優勢T-RFs，這可能表明較高含量的PHs有機污染物對土壤中硫酸鹽還原菌或鐵還原菌具有一定富集作用。Alejandro等［21］也發現，海洋沉積物中受原油污染程度越大、時間越久，bssA基因的多樣性越高。結合上述分析可知，硫酸鹽還原或鐵還原條件可能是重度PHs污染土壤中污染物去除的重要厭氧代謝途徑［31］。在PAHs含量僅次于JH-5的JH-1中，雖然僅檢測到3種bssA類群，但它們均出現在JH-5的bssA類群中，且96bp（不確定類群）的bssA類群只出現在它們中，這可能體現出PHs類污染物對土壤bssA類群的選擇性。此外，在PAHs含量較低的樣品中也出現了上述優勢bssA類群，如409bp和433bp 也是JH-2的優勢bssA類群，而JH-4的bssA類群多于JH-1，且其bssA類群均出現在JH-5中。可能原因有以下兩點：首先，所有的土壤樣品均采自同一區域，且類型相同，因此其土著微生物群落結構（包含bssA類群）也應該具有一定的同源性。其次，除PAHs外，土壤的其他理化性質，如NO3-、NH4+、有效磷等也會對微生物的群落結構產生影響［35］。JH-3中僅檢測到一個bssA類群，這可能是由于某些未知的土壤環境因子或人為干擾所導致。隨后的RDA分析（圖3）進一步表明，除能夠作為微生物生長所必須的銨態氮、硝態氮及有效磷外，土壤PAHs含量也是影響bssA群落結構的重要環境因子。與T-RFLP結果一致的是，96bp（不確定類群）和388bp（硫酸鹽還原菌或地桿菌）T-RFs所代表的bssA類群與土壤PAHs含量呈正相關，表明PHs污染可能對它們有較強的富集作用。綜上所述，江漢油田區域土壤含有PHs類污染物厭氧降解的遺傳基礎，可能通過不同的細菌在多種代謝途徑下進行，如以硝酸鹽為電子受體的反硝化作用，以硫酸鹽為電子受體的硫酸鹽還原作用及以鐵離子為電子受體的鐵還原作用等。同時，土壤環境條件及受污染程度均可能影響bssA類群及其對PHs的降解效率。

4 結 論

本研究以典型PHs類污染物-苯系物厭氧代謝的功能基因bssA作為標識，通過克隆文庫結合T-RFLP的方法，研究了典型油田區域土壤

PHs降解微生物的遺傳多樣性。該油田區主要的bssA類群為β-變形菌和δ-變形菌，包括索氏菌屬（Thauera）、固氮弓菌屬（Azoarcus）、硫酸鹽還原菌以及地桿菌屬（Geobacter），已知的這些微生物具有硝酸鹽還原、硫酸鹽還原及鐵還原代謝能力，可能是江漢油田區土壤厭氧微域中PHs降解的主要微生物類群。土壤氮、磷及PHAs含量均是影響bssA基因遺傳多樣性的關鍵因子，較高含量的PAHs可能對含有bssA基因的硫酸鹽還原菌及地桿菌屬的微生物具有一定的富集作用。

［1］ Anderson R T，Lovley D R. Ecology and biogeochemistry of in situ groundwater bioremediation. Advances in Microbial Ecology，1997，15：289—350

［2］ Lin B，van Verseveld H W，R?ling W F. Microbial aspects of anaerobic BTEX degradation. Biomedical and Environmental Sciences，2002，15（2）：130—144

［3］ Chakraborty R，Coates J D. Anaerobic degradation of monoaromatic hydrocarbons. Applied Microbiology and Biotechnology，2004，64（4）：437—446

［4］ Cai Q Y，Mo C H，Li Y H，et al. Occurrence and assessment of polycyclic aromatic hydrocarbons in soils from vegetable fields of the Pearl River Delta，South China. Chemosphere，2007，68（1）：159—168

［5］ Puglisi E，Cappa F，Fragoulis G，et al. Bioavailability and degradation of phenanthrene in compost amended soil. Chemosphere，2007，67（3）：548—556

［6］ Heider J，Spormann A M，Beller H R，et al. Anaerobic bacterial metabolism of hydrocarbons. FEMS Microbiology Reviews，1998，22（5）：459—473

［7］ Kleinsteuber S，Schleinitz K M，Vogt C. Key players and team play：Anaerobic microbial communities in hydrocarbon-contaminated aquifers. Applied Microbiology and Biotechnology，2012，94（4）：851—873

［8］ Grossi V，Cravo L C，Guyoneaud R，et al. Metabolism of n-alkanes and n-alkenes by anaerobic bacteria：A summary. Organic Geochemistry，2008，39（8）：1197—1203

［9］ Heider J. Adding handles to unhandy substrates：Anaerobic hydrocarbon activation mechanisms. Current Opinion in Chemical Biology，2007，11（2）：188—194

［10］ Beller H R，Kane S R，Legler T C，et al. A realtime polymerase chain reaction method for monitoring anaerobic，hydrocarbon-degrading bacteria based on a catabolic gene. Environmental Sciences andTechnology，2002，36（18）：3977—3984

［11］ Biegert T，Fuchs G，Heider J. Evidence that anaerobic oxidation of toluene in the denitrifying bacterium Thauera aromatica is initiated by formation of benzylsuccinate from toluene and fumarate. European Journal of Biochemisty，1996，238（3）：661—668

［12］ Beller H R，Spormann A M. Anaerobic activation of toluene and o-xylene by addition to fumarate in denitrifying strain T. Journal of Bacteriology，1997，179（3）：670—676

［13］ Leuthner B，Leutwein C，Schulz H，et al. Biochemical and genetic characterization of benzylsuccinate synthase from Thauera aromatica：A new glycyl radical enzyme catalysing the first step in anaerobic toluene metabolism. Molecular Microbiology，1998，28（3）：615—628

［14］ Rabus R，Wilkes H，Behrends A，et al. Anaerobic initial reaction of n-alkanes in a denitrifying bacterium：Evidence for（1-methylpentyl）succinate as initial product and for involvement of an organic radical in n-hexane metabolism. Journal of Bacteriology，2001，183（5）：1707—1715

［15］ Kropp K G，Davidova I A，Suflita J M. Anaerobic oxidation of n-dodecane by an addition reaction in a sulfate-reducing bacterial enrichment culture. Applied and Environmental Microbiology，2000，66（12）：5393—5398

［16］ Selesi D，Jehmlich N，Bergen M，et al. Combined genomic and proteomic approaches identify gene clusters involved in anaerobic 2-methylnaphthalene degradation in the sulfate-reducing enrichment culture N47. Journal of Bacteriology，2010，192（1）：295—306

［17］ Musat F，Galushko A，Jacob J，et al. Anaerobic degradation of naphthalene and 2-methylnaphthalene by strains of marine sulfate-reducing bacteria. Environmental Microbiology，2009，11（1）：209—219

［18］ Winderl C，Schaefer S，Lueders T. Detection of anaerobic toluene and hydrocarbon degraders in contaminated aquifers using benzylsuccinate synthase （bssA）genes as a functional marker. Environmental Microbiology，2007，9（4）：1035—1046

［19］ Callaghan A V，Davidova I A，Savage-Ashlock A K，et al. Diversity of benyzl- and alkylsuccinate synthase genes in hydrocarbon-impacted environments and enrichment cultures. Environmental Science and Technology，2010，44（19）：7287—7294

［20］ Frederick V N，Giovanni P，Sara K，et al. Enhanced gene detection assays for fumarate-adding enzymes allow uncovering of anaerobic hydrocarbon degraders in terrestrial and marine systems. Applied and Environmental Microbiology，2013，79（2）：543—552

［21］ Alejandro A G，Ramon R M，Silvia M. Diversity of benzylsuccinate synthase-like（bssA）genes in hydrocarbon-polluted marine sediments suggests substrate-dependent clustering. Applied and Environmental Microbiology，2013，79（12）：3667—3676

［22］ Sun W M，Sun X X，Alison M，et al. Presence，diversity and enumeration of functional genes（bssA and bamA）relating to toluene degradation across a range of redox conditions and inoculum sources. Biodegradation，2014，25（2）：189—203

［23］ Johnsen A R，Wick L Y，Harms H. Principles of microbial PAH-degradation in soil. Environmental Pollution，2005，133（1）：71—84

［24］ 魯如坤. 土壤農業化學分析方法. 北京：中國農業科技出版社，2000

Lu R K. Analytical methods for soil and agro-chemistry （In Chinese）. Beijing：China Agricultural Science and Technology Press，2000

［25］ 趙有政，裔兆軍. 江蘇鹽城灘涂土壤中多環芳烴的測定. 污染防治技術，2011，24（5）：41—43

Zhao Y Z，Yi Z J. Determination of polycyclic aromatic hydrocarbons（PAHs）in the Coastal soil of Yancheng in Jinagsu Province（In Chinese）. Pollution Control Technology，2011，24（5）：41—43

［26］ Chen X，Zhang L M，Shen J P，et al. Abundance and community structure of ammonia-oxidizing archaea and bacteria in an acid paddy soil. Biology and Fertility of Soils，2011，47（3）：323—331

［27］ 中華人民共和國環境保護部. 土壤環境質量標準（征求意見稿）. GB15618-2008，2008

Ministry of Environmental Protection of the People's Republic of China. Environmental quality standard for soils（exposure draft）（In Chinese）. GB15618-2008，2008

［28］ 師陽，孫玉川，梁作兵，等. 重慶金佛山土壤中PAHs含量的海拔梯度分布及來源解析. 環境科學，2015，36（4）：1417—1424

Shi Y，Sun Y C，Liang Z B，et al. Altitudinal gradient distribution and source analysis of PAHs content of topsoil in Jinfo Mountian，Chongqing（In Chinese）. Environmental Science，2015，36（4）：1417—1424

［29］ 劉增俊，騰應，黃標，等. 長江三角洲典型地區農田土壤多環芳烴分布特征與源解析. 土壤學報，2010，47 （6）：1110—1117

Liu Z J，Teng Y，Huang B，et al. Distributionand sources analysis of PAHs in farmland soils in areas typical of the Yangtze River Delta，China（In Chinese）. Acta Pedologica Sinica，2010，47（6）：1110—1117

［30］ Shinoda Y，Sakai Y，Uenishi H，et al. Aerobic and anaerobic toluene degradation by a newly isolated denitrifying bacterium，Thauera sp. strain DNT-1. Applied and Environmental Microbiology，2004，70 （3）：1385—1392

［31］ Valderrama J A，Shingler V，Carmona M，et al. AccR is a master regulator involved in carbon catabolite repression of the anaerobic catabolism of aromatic compounds in Azoarcus sp. CIB. Journal of Biological Chemistry，2014，289（4）：1892—1904

［32］ Boll M，L?ffler C，Morris B E L，et al. Anaerobic degradation of homocyclic aromatic compounds via arylcarboxyl-coenzyme A esters：Organisms，strategies and key enzymes. Environmental Microbiology，2014， 16（3）：612—627

［33］ Mechichi T，Stackebrandt E，Gad’on N，et al. Phylogenetic and metabolic diversity of bacteria degrading aromatic compounds under denitrifying conditions，and description of Thauera phenylacetica sp. nov.，Thauera aminoaromatica sp. nov.，and Azoarcus buckelii sp. nov. Archives of Microbiology，2002，178（1）：26—35

［34］ Kunapuli U，Jahn M K，Lueders T，et al. Desulfitobacterium aromaticivorans sp. nov. and Geobacter toluenoxydans sp. nov.，iron-reducing bacteria capable of anaerobic degradation of monoaromatic hydrocarbons. International Journal of Systematic and Evolutionary Microbiology，2010，60（3）：686—695

［35］ Sarkar J，Kazy S K，Gupta A，et al. Biostimulation of indigenous microbial community for bioremediation of petroleum refinery sludge. Frontiers in Microbiology，2016，7，DOI：10.3389/fmicb.2016.01407

Genetic Diversity of Hydrocarbons Degrading Microbial Functional Gene （bssA）in the Farmland Soil Typical of Jianghan Oilfield

ZUO Xiaohu WANG Mingxia YAO Yanhong LI Zhenlun ZHOU Zhifeng?
（College of Resources and Environment，Southwest University，Chongqing 400716，China）

Farmland soils in oil fields are liable to get polluted with petroleum hydrocarbons（PHs），which consequently leads to enrichment of some special PHs degrading microbial groups in the soils. In this study，soil samples were collected from typical farmland fields near five different oil wells in the Jianghan Oilfield located in Qianjiang City，Hubei Province for analysis of PHs degrading microbial groups，using bssA （benzylsuccinate synthase gene），a functional gene key to anaerobic benzenes degradation as biomarker，and community structure of the PHs degrading microbial groups was determined with T-RFLP（Terminalrestriction fragment length polymorphism）and clone library. On such a basis，environmental factors affecting community composition of the soil bssA were discussed. Results show that the content of polycyclic aromatic hydrocarbons（PAHs）in the soil varied in the range of 0.21～2.01 mg kg-1，which indicates that the soil in this region was not seriously contaminated with PHs. Meanwhile，the T-RFLP analysis indicates that bssA diversity varied significantly from soil sample to soil sample in this oilfield，being the highest in the soil sample the highest in PAHs content，and its dominant group was in a fairly close kinship with sulfatereducing bacterium or Geobacter. Furthermore，RDA（Redundancy analysis）reveals that the contents of soil available nitrogen，phosphorus and PAHs were all main factors affecting soil bssA diversity. To sum up，all the findings indicate that the bssA-bearing microbial groups in the soil of this oilfield are Beta-proteobacteria and Deta-proteobacteria，and in close kinships with Thauera，Azoarcus and Geobacter in phylogenesis. All these microbes might degrade PAHs through reducing metabolic processes of nitrate，sulfate，and iron.

Soil；Anaerobic degradation；Benzylsuccinate synthase gene（bssA）；Community structure

S154.1

A

（責任編輯：盧 萍）

* 國家自然科學基金項目（41371477）和中央高校專項基金項目（XDJK2014B047）資助 Supported by the National Natural Science Foundation of China（No. 41371477）and the Fundamental Research Funds for the Central Universities（No. XDJK2014B047）

? 通訊作者 Corresponding author，E-mail：zhouzhf@swu.edu.cn

左小虎（1990—），男，湖北漢川人，碩士研究生，主要從事土壤微生物學研究。E-mail：375332084@qq.com

2016-11-04；

2017-01-10；優先數字出版日期（www.cnki.net）：2017-02-28

10.11766/trxb201611040324