摘 要：利用MiSeq高通量测序技术分析了延长某页岩气井田场地老化石油污染土壤的微生物群落多样性，采用冗余分析探究了微生物群落结构与土壤环境因子的关系。结果表明，石油污染可改变土壤微生物群落结构及多样性，石油污染土壤微生物的优势菌门为变形菌门(Proteobacteria)、放线菌门(Actinobacteria);鞘氨醇单胞菌属(Sphingomonas)、分枝杆菌属(Mycobacterium)、Saccharibacteria＿genera＿incertae＿sedis、芽孢杆菌属(Bacillus)、节杆菌属(Arthrobacter)、交替红色杆菌属(Altererythrobacter)和诺卡氏菌属(Nocardioides)等7种菌属具有石油烃类污染物降解功能，分枝杆菌属和短波单胞杆菌属(Brevundimonas)具有重金属降解功能；影响微生物群落多样性的环境主控因子为pH值、石油烃含量和重金属含量，石油烃的摄入会导致土壤中重金属含量升高。
关键词:老化土壤;石油污染;微生物群落;多样性;环境主控因子;
Characteristics of Microbial Community Diversity in Aged Petroleum-contaminated Soil
GUAN Shengdi SU Yingxuan LI Mengli
LUO Haimei YU Yan
School of Chemical & Environmental Engineering,China University of Mining and
Technology(Beijing)
Abstract：We analyzed the microbial community diversity in the aged petroleum-contaminated soil of the well field of Yanchang Shale Gas by MiSeq high-throughput sequencing technology, and explored the relationship between microbial community structure and soil environmental factors through redundancy analyses.The results show that petroleum pollution can change the structure and diversity of microbial community in soil.Proteobacteria and Actinobacteria are the most abundant at phylum level.Sphingomonas,Mycobacterium,Saccharibacteria＿genera＿incertae＿sedis,Bacillus,Arthrobacter,Altererythrobacter,and Nocardioides have the function of degrading petroleum hydrocarbon pollutants, while Mycobacterium and Brevundimonas have the function of degrading heavy metals.The main-control environmental factors of microbial community diversity are pH value, petroleum hydrocarbon content, and heavy metal content, and the intake of petroleum hydrocarbon can lead to the increase of heavy metal content in soil.
Keyword：aged soil; petroleum contamination; microbial community; diversity; main-control environmental factor;
我国经济高速发展，能源需求迅速增长，能源供需矛盾日益突出，非常规能源越来越受到人们的重视[1]。但是页岩气开采造成土壤中石油含量增加，显著影响微生物群落，导致土壤微生物数量减少及生物多样性降低[2]。土壤微生物是土壤的重要组成部分，其种类较多，在土壤的形成发育、物质循环和肥力演变等过程中发挥关键作用。页岩气开采造成的土壤污染具有隐蔽性、滞后性、累积性、不可逆转性等特点[3]。因此，切实有效地解决石油污染土壤的修复问题，成为国内外研究的热点[4]。石油污染土壤的修复技术主要有物理修复、化学修复和生物修复[5]。其中，生物修复技术因具有成本低、几乎不产生二次污染和可进行原位修复等优点而受到广泛关注[6,7,8]。
延长油气田位于鄂尔多斯盆地，鄂尔多斯盆地延长组长 7 段烃源岩分布广泛，处于生油窗后期——生湿气窗的高峰期[9],油气共生，部分成熟度较高的区域产气强度较大，具备页岩气成藏的物质基础，因此，页岩气开采发展迅速。但在页岩气开采过程中产生的油基岩屑会污染土壤[10]。2015年，我国出台了禁止油基岩屑落地的政策，因此后期没有新增石油污染土壤，但是了解石油污染土壤的污染现状对于指导土壤生物修复具有重要意义。 基于此，作者采集延长某页岩气井田场地老化石油污染土壤，以周围未受石油污染土壤作为对照，利用 MiSeq 高通量测序技术分析土壤微生物群落多样性，采用冗余分析(RDA)探究微生物群落结构与土壤环境因子的关系，为探究土壤微生物对石油的耐性机制、开发土壤修复技术提供理论支撑。
1 实验
1.1 土壤样品采集
2020年10月在陕西省延长某页岩气井田场地按照棋盘布点法设置6个点位，每个点位采集50 cm×50 cm小样方的四角及中心0～20 cm深的土壤样品，混合得到该点位的土壤样品，依次编号为YC1～YC6;同时在井田场地外采集未污染土壤样品作为对照，编号为YC。所有土壤样品剔除草根和石块后储存于采样袋内，置于装有冰袋的恒温箱中，立即带回实验室4 ℃保存，用于微生物高通量测序的土壤样品于-80 ℃下保存。
1.2 土壤样品的理化性质分析
将采集的土壤样品自然风干，研磨过筛，采用火焰原子吸收分光光度法测定总铬，采用石墨炉原子吸收分光光度法测定铅，采用碳酸氢钠浸提-钼锑抗分光光度法测定有效磷，采用重铬酸钾氧化-分光光度法测定有机碳，采用凯氏定氮法测定总氮，采用电位法测定pH值，采用重量法测定含水率和含盐量，采用气相色谱法测定石油烃。采用spss 22软件对土壤样品的各项理化指标的变化量进行相关性分析。
1.3 微生物DNA的提取及高通量测序
对土壤样品进行预处理，用DNA提取试剂盒(E.Z.N.ATM Mag-Bind Soil DNA Kit, Omega)提取DNA。对样本进行质检，样本浓度在6.04～24.60 ng·μL-1之间。采用16S rRNA通用引物341F(5′-CCTACGGGNGGCWGCAG-3′)、805R(5′-GACTACHVGGGTATCTAATCC-3′)对V3-V4区目的基因进行PCR扩增，在第二轮扩增中引入Illumina桥式PCR兼容引物。最后，在MiSeq平台对16S rRNA的V3-V4高变区进行高通量测序。
1.4 测序数据分析
对MiSeq测序所得原始数据进行拼接，区分样本后对序列质量进行质控和过滤，然后对所得操作单元(operational taxonomic unit, OTU)进行聚类分析和物种分类学分析。基于OTU聚类分析结果，对OTU进行多种多样性指数分析，计算alpha多样性指数(Chao1指数、Shannon指数、Simpson指数等),以表征微生物群落多样性，鉴别采样点的石油降解优势菌属。利用R语言绘制OTU聚类韦恩图、聚类堆积柱状图，并对优势菌门、菌属的丰度与理化性质进行冗余分析。
2 结果与讨论
2.1 土壤样品的理化性质(表1)
表1 土壤样品的理化性质
采样点地处西北，干旱少雨，盐碱积累，土壤的pH值和含盐量高。由表1可知，污染土壤样品的pH值范围在8.45～9.98之间，均高于未污染土壤样品；污染土壤样品的有机碳含量在0.81～5.18 g·kg-1之间，均高于未污染土壤样品；未污染土壤样品的含盐量较低，远远低于污染土壤样品，但含水率远远高于污染土壤样品，主要是因为石油具有疏水性，石油污染土壤的疏水性也随之升高，导致土壤储水能力下降，含水率除低[11];污染土壤样品的铅含量在3.03～17.10 mg·kg-1之间；污染土壤样品的石油烃含量在38.81～320.69 mg·kg-1之间，均高于未污染土壤样品，其中YC4样品的石油烃含量最高，是未污染土壤样品的40倍左右。石油烃含量较高的土壤样品(YC1、YC3、YC4、YC6)的碳氮比失衡，均高于100/10,主要是由于陕北地区黄绵土氮含量低，碳含量高，导致碳氮比较高[12]。
将污染土壤样品与未污染土壤样品的理化性质进行Person相关性分析，结果见表2。
表2 土壤理化性质的Person相关性分析
由表2可知，pH值与有机碳含量呈极显著负相关性，总氮含量与有效磷含量、铅含量与石油烃含量均呈极显著正相关性，有效磷含量与含水率呈显著正相关性，其它指标间的相关性不显著。
2.2 微生物群落多样性分析
2.2.1 OTU聚类分析(图1)
由图1可知，7个土壤样品共同含有532个 OTU。OTU数量顺序为：YC>YC6>YC3>YC1>YC4>YC5≈YC2,其中YC样品含有的OTU数量最多，YC5和YC2样品含有的OTU数量最少，说明石油污染土壤中微生物数量减少，主要是受石油污染物的胁迫和环境因素的影响[13]。
2.2.2 多样性指数分析
在97%相似度水平下，计算alpha多样性指数，结果见表3。
表3 alpha多样性指数
Chao1指数和 Shannon指数反映了微生物群落的丰富度和多样性，数值越大说明群落丰富度和多样性越高； Simpson 指数反映了微生物群落的均匀性，数值越小，说明群落分布越均匀。由表3可知，YC和YC6样品的Chao1指数和Shannon指数较大，说明微生物群落丰富度和多样性较高；YC样品的Simpson 指数最小，说明未污染土壤样品的微生物群落均匀性高于污染土壤样品。对比发现，石油烃的摄入降低了土壤微生物丰富度，石油烃含量越高，多样性越低，群落分布越不均匀。
2.2.3 石油降解微生物优势菌落分析
门水平物种相对丰度如图 2 所示。
由图2可知，变形菌门(Proteobacteria)的相对丰度最高，为优势菌门，广泛存在于石油污染土壤中，说明其对老化石油污染土壤环境具有较强的适应性；放线菌门(Actinobacteria)在7个样品中同样表现出了丰度优势；此外，酸杆菌门(Acidobacteria)、拟杆菌门(Bacteroidetes)、芽单胞菌门(Gemmatimonadetes)在7个样品中均有一定丰度。污染土壤样品中， YC3样品与其它样品在门水平上有较大的差异，其余样品微生物组成略有差异，YC2与YC6样品在门水平上一致。相较于YC样品，由于石油烃类污染物的摄入，变形菌门、绿弯菌门(Chloroflexi)的相对丰度明显上升，其中石油烃含量较高的YC4样品中变形菌门微生物的增加尤为明显，而放线菌门的相对丰度有所降低。表明在石油污染土壤的老化过程中，变形菌门微生物增多，而放线菌门微生物减少。
当土壤被石油烃污染时，相对丰度显著升高或新出现的细菌可以被视为石油烃降解菌(嗜油性细菌),而相对丰度显著降低的细菌可以被视为疏油性细菌[14]。因此推断延长地区嗜油性菌门为变形菌门和绿弯菌门，这与辽河石油污染土壤微生物群落结构多样性的研究[15]及陕北石油污染土壤微生物多样性的研究[16]结果一致。
将相对丰度前15的物种进行属水平分析，结果如图 3 所示。
由图3可知，未污染土壤样品与污染土壤样品中优势菌群组间差异较大。不同采样点属水平上群落结构有差异，其中YC2与YC6样品在属水平上相似性最高。根据各样品理化性质，可以看出属水平相似性较高的样品，其pH值、有机碳含量相似性最高。污染土壤样品中主要菌属有鞘氨醇单胞菌属(Sphingomonas)、分枝杆菌属(Mycobacterium)、Saccharibacteria_genera_incertae_sedis、芽孢杆菌属(Bacillus)、节杆菌属(Arthrobacter)、交替红色杆菌属(Altererythrobacter)、诺卡氏菌属(Nocardioides)。
对属水平占比较高的优势菌群进行功能性分析，结果见表4。
表4 占比较高的优势菌群的功能性分析
石油烃降解菌能将碳氢化合物转化为能源、生物质及生物废料副产品[26]。石油污染土壤样品中，共有7个属具有石油烃类污染物降解功能，其中诺卡氏菌属和芽孢杆菌属可产表面活性剂，降低石油烃类污染物表面张力，实现其乳化和增溶；鞘氨醇单胞菌属含有一组独特的芳香化合物降解基因，可通过调节细胞表面疏水性来适应石油污染环境。
假单胞菌属(Pseudomonas)、红球菌属(Rhodococcus)、微杆菌属(Microbacterium)、芽孢杆菌属、产碱杆菌属等是常见的石油烃降解菌[22]。王耀民[27]从大港油田长期被原油污染的土壤中分离、筛选出4 株能以原油为唯一碳源的菌株，其中2株原油降解率高的菌株分别为不动杆菌属和芽孢杆菌属。
2.2.4 微生物群落结构与土壤环境因子的关系
微生物群落结构是土壤微生态环境的响应，要结合理化指标进行分析。对7个土壤样品微生物门和属水平与理化性质进行冗余分析，结果如图4所示。
由图4a可知，第一主轴和第二主轴分别解释了总方差变量的45.38%和30.21%,第一主轴上总氮含量是主要的影响因子，第二主轴上pH值和总铬含量是主要的影响因子。其中pH值、总氮含量、总铬含量是微生物门水平的主要影响因子，石油烃含量也对微生物门水平有影响，尤其是与变形菌门具有显著的正相关性。土壤中微生物数量最多的变形菌门受pH值的影响较大。铅含量与石油烃含量存在显著的相关性，说明石油污染会导致土壤中重金属含量升高，进而加重土壤污染，与Person相关性分析(表2)结果一致。受石油烃的影响，微生物菌门在结构上会发生变化， 其中YC2与YC6样品的菌门结构相似性高，其微生物菌门主要受石油烃含量、铅含量、pH值的影响；YC4样品主要受总氮含量的影响；YC1与YC5样品的菌门结构相似性较高，与门水平聚类堆积柱状图(图2)结果一致。
由图4b可知，第一主轴和第二主轴分别解释了总方差变量的45.88%和27.69%。第一主轴上pH值、铅含量、石油烃含量是主要的影响因子，第二主轴上石油烃含量和总铬含量是主要的影响因子，因此，与微生物群落结构相关性较高的土壤理化指标是石油烃含量、pH值和重金属含量。Saccharibacteria_genera_incertae_sedis、分枝杆菌属与石油烃含量存在显著正相关性，推测其对石油烃具有良好的降解潜能；分枝杆菌属与铅含量具有显著正相关性，其对铅具有良好的降解性能；鞘氨醇单胞菌属与pH值具有显著正相关性，与石油烃含量具有正相关性，与荆佳维等[13]的研究结果一致。Gp3、短波单胞杆菌属(Brevundimonas)与总铬含量存在显著正相关性，推测其对总铬有降解能力。YC1、YC2、YC5、YC6样品的菌属相似性较高，YC3、YC4样品与其它样品不同，未污染土壤与污染土壤样品的菌属差异性大。
3 结论
利用 MiSeq 高通量测序技术分析了延长某页岩气井田场地老化石油污染土壤的微生物群落多样性，采用冗余分析探究了微生物群落结构与土壤环境因子的关系，结论如下：
(1)石油污染土壤偏碱性，土壤污染不均匀，石油烃含量差异较大，由于石油烃污染造成碳氮比失衡，氮含量不足。
(2)石油污染土壤微生物的优势菌门为变形菌门和放线菌门；鞘氨醇单胞菌属、分枝杆菌属、Saccharibacteria_genera_incertae_sedis、芽孢杆菌属、节杆菌属、交替红色杆菌属和诺卡氏菌属等 7 种菌属具有石油烃类污染物降解功能，分枝杆菌属和短波单胞杆菌属具有重金属降解功能。
(3)石油污染土壤微生物的多样性和均匀度均有所下降，石油污染会导致鞘氨醇单胞菌属、分枝杆菌属、芽孢杆菌属等菌属的相对丰度有所上升。
(4)石油污染土壤微生物群落多样性主要受pH值、石油烃含量和重金属含量的影响，铅含量与石油烃含量存在显著正相关性，说明石油烃的摄入会导致土壤中重金属含量上升、营养失衡，进而加重土壤污染。
[1] 程鹏立.页岩气勘查开发的环境风险及其应对机制研究[J].南京工业大学学报(社会科学版),2016,15(2):61-67.CHENG P L.Environmental risks and prevention mechanism of shale gas exploration[J].Journal of Nanjing Tech University(Social science Edition),2016,15(2):61-67.
[2] 吴盼云，晁群芳，赵亚光，等.克拉玛依石油污染土壤微生物群落结构及其代谢特征[J].基因组学与应用生物学，2019,38(5):2062-2069.WU P Y,CHAO Q F,ZHAO Y G,et al.Microbial community structure and metabolic characteristics of oil-contaminated soil in Karamay[J].Genomics and Applied Biology,2019,38(5):2062-2069.
[3] 陶贵荣，王莉衡.浅析土壤污染与修复技术[J].化学与生物工程，2010,27(6):4-6,12.TAO G R,WANG L H.Disscussion on soil pollution and remediation technology[J].Chemistry & Bioengineering,2010,27(6):4-6,12.
[4] OSSAI I C,AHMED A,HASSAN A,et al.Remediation of soil and water contaminated with petroleum hydrocarbon:a review[J].Environmental Technology & Innovation,2020,17:100526.
[5] 毛健，杨代凤，刘腾飞，等.多环芳烃污染土壤的菌群-植物联合修复效应研究[J].环境与可持续发展，2017,42(4):108-110.MAO J,YANG D F,LIU T F,et al.Interactive effect of microbial consortium and three plants in remediation of PAHs contaminated soil[J].Environment and Sustainable Development,2017,42(4):108-110.
[6] 张晓庆，徐丽，齐悦，等.紫松果菊对多环芳烃重污染土壤修复效能[J].生态学杂志，2018,37(2):492-497.ZHANG X Q,XU L,QI Y,et al.Remediation efficiency of Echinacea purpurea for heavy PAHs contaminated soils[J].Chinese Journal of Ecology,2018,37(2):492-497.
[7] SARMA H,NAVA A R,PRASAD M N V.Mechanistic understanding and future prospect of microbe-enhanced phytoremediation of polycyclic aromatic hydrocarbons in soil[J].Environmental Technology & Innovation,2019,13:318-330.
[8] CZARNY J,STANINSKA-PIETA J,PIOTROWSKA-CYPLIK A,et al.Acinetobacter sp.as the key player in diesel oil degrading community exposed to PAHs and heavy metals[J].Journal of Hazardous Materials,2020,383:121168.
[9] 张洪洁，王凤琴，刘航，等.鄂尔多斯盆地东南部长7段页岩气含气量及影响因素[J].新疆石油地质，2019,40(6):666-672.ZHANG H J,WANG F Q,LIU H,et al.Shale gas content and its influencing factors in Chang 7 Member of southeastern Ordos Basin[J].Xinjiang Petroleum Geology,2019,40(6):666-672.
[10] 李学林，杨诗雨.页岩气开发中的环境问题及对策——以长宁页岩气开发区为例[J].环境保护与循环经济，2020,40(6):78-84.LI X L,YANG S Y.Environmental problems and countermeasures in shale gas development:a case study of Changning shale gas development area[J].Liaoning Urban and Rural Environmental Science & Technology,2020,40(6):78-84.
[11] 贾建丽，刘莹，李广贺，等.油田区土壤石油污染特性及理化性质关系[J].化工学报，2009,60(3):726-732.JIA J L,LIU Y,LI G H,et al.Contamination characteristics and its relationship with physicochemical properties of oil polluted soil in oilfields of China[J].CIESC Journal,2009,60(3):726-732.
[12] 许殷瑞，吴蔓莉，王丽，等.陕北石油污染土壤微生物种群变化及影响因素[J].中国环境科学，2021,41(9):4349-4359.XU Y R,WU M L,WANG L,et al.The influences of petroleum pollution on the microbial population distribution in Northern Shaanxi Province of China[J].China Environmental Science,2021,41(9):4349-4359.
[13] 荆佳维，王卅，郭书海.典型油田区油污土壤微生物群落区域性分布研究[J].环境科学学报，2021,41(11):4660-4675.JING J W,WANG S,GUO S H.Regional distribution of crude oil-contaminated soil microbial communities in typical oilfields[J].Acta Scientiae Circumstantiae,2021,41(11):4660-4675.
[14] 王宝山，温成成，孙秦川，等.石油烃类污染对青藏高原北麓河地区冻区土壤微生物多样性的影响[J].环境工程学报，2018,12(10):2917-2928.WANG B S,WEN C C,SUN Q C,et al.Effects of petroleum hydrocarbon pollution on soil microbial diversity in permafrost region of Beiluhe area of Qinghai-Tibet plateau[J].Chinese Journal of Environmental Engineering,2018,12(10):2917-2928.
[15] 赵明阳，王卅，李凤梅，等.辽河油田土壤石油污染及其微生物群落特征[J].应用生态学报，2020,31(12):4215-4224.ZHAO M Y,WANG S,LI F M,et al.Petroleum pollution and microbial community structure in the soil of Liaohe Oilfield[J].Chinese Journal of Applied Ecology,2020,31(12):4215-4224.
[16] 张宝宝.陕北石油污染土壤微生物多样性及降解菌特性研究[D].延安：延安大学，2021.ZHANG B B.Microbial diversity of petroleum contaminated soil in Northern Shaanxi and characteristics of isolated petroleum[D].Yan′an:Yan′an University,2021.
[17] ZHOU L S,LI H,ZHANG Y,et al.Sphingomonas from petroleum-contaminated soils in Shenfu,China and their PAHs degradation abilities[J].Brazilian Journal of Microbiology,2016,47(2):271-278.
[18] CHEN L,CHEN W F,XU Z L,et al.Sphingomonas oleivorans sp.nov.,isolated from oil-contaminated soil[J].International Journal of Systematic & Evolutionary Microbiology,2018,68(12):3720-3725.
[19] DAI Y Y,LIU R,ZHOU Y M,et al.Fire phoenix facilitates phytoremediation of PAH-Cd co-contaminated soil through promotion of beneficial rhizosphere bacterial communities[J].Environment International,2020,136:105421.
[20] WU M L,MA C,WANG D,et al.Nutrient drip irrigation for refractory hydrocarbon removal and microbial community shift in a historically petroleum-contaminated soil[J].Science of the Total Environment,2020,713:136331.
[21] 王娟丽.给水厂污泥改良生物滞留系统对氮磷去除的优化探究[D].北京：北京建筑大学，2019.WANG J L.Study on optimization of nitrogen and phosphorus removal in bioretention systems modified with water treatment residuals[D].Beijing:Beijing University of Civil Engineering and Architecture,2019.
[22] 沈聪，刘爽，王春霞，等.盐池采油区污染土壤石油降解菌的筛选鉴定及其降解特性[J].生物技术通报，2021,37(6):127-135.SHEN C,LIU S,WANG C X,et al.Screening,identification and characteristics of petroleum degrading bacteria from the contaminated soil of Yanchi aera[J].Biotechnology Bulletin,2021,37(6):127-135.
[23] 李娟，UWAREMWE C,冷艳，等.节杆菌属细菌处理有机物和重金属污染物的研究进展[J].环境科学与技术，2017,40(10):89-97.LI J,UWAREMWE C,LENG Y,et al.Progress on the study of biodegradation of organic pollutants and adsorption of heavy metals with Arthrobacter strains[J].Environmental Science & Technology,2017,40(10):89-97.
[24] 马聪，秦丹，孙倩，等.耐盐类固醇激素降解菌交替赤杆菌MH-B5的降解特性、降解途径及其固定化[J].环境科学，2015,36(12):4651-4658.MA C,QIN D,SUN Q,et al.Degradation of steroidal hormones by salt-tolerant Altererythrobacter strain MH-B5:degradation characteristics,metabolites and its immobilization[J].Environmental Science,2015,36(12):4651-4658.
[25] 李文均，张永光.拟诺卡氏菌属放线菌研究进展[J].微生物学通报，2016,43(5):1123-1135.LI W J,ZHANG Y G.Advances in studies on the genus Nocardiopsis[J].Microbiology China,2016,43(5):1123-1135.
[26] VARJANI S J,GNANSOUNOU E,PANDEY A.Comprehensive review on toxicity of persistent organic pollutants from petroleum refinery waste and their degradation by microorganisms[J].Chemosphere,2017,188:280-291.
[27] 王耀民大港油田污染土壤烷烃高效降解菌的筛选及特性研究[D].石家庄:河北地质大学，2020.WANG Y M.Screening and characterization of alkane degrading bacteria in polluted soil from Dagang Oilfield [D] .Shijiazhuang:Hebei GEO University,2020.