氨氮污染源区包气带剖面微生物生态分布特征的研究.docx

1、Vol.32,No.2296305文章编号：0254-5357(2013)02-0296-10氨氮污染源区包气带剖面微生物生态分布特征的研究程雅楠口，关翔宇'，曲文龙3,陈鸿汉，刘菲），谢宇轩'，朱玲玲'（1.水资源与环境工程北京市重点实验室，中国地质大学（北京），北京100083；2.中交铁道勘察设计院有限公司，北京100088；3.沈阳水务集团，辽宁沈阳110003）摘要：包气带土壤中氨氮污染迁移规律已成为近年来国内外学者研究的重点，目前在氨氮污染迁移规律研完中采用柱实验模拟以及软件模拟的方法较多，土壤生物技术则广泛应用于污染物降解领域，而在污染物迁移规律研完中应用

2、较少。本文采用变性梯度凝胶电泳（DGGE）技术、16SrRNA序列分析技术以及典范时应分析相结合,对华北平原3个典型氨氮污染区土壤表层至包气带剖面微环境的细菌垂直分布特征及群落结构进行研究。结合污染区土壤理化性质分析，认为包气带上壤剖面的细菌群落中存在着与氮循环、硫酸盐代谢等过程偶联的优势细菌类群，说明土壤微环境中细菌群落分布明显受氨态氮、硝态氮、亚硝态氮和硫酸盐的分布影响，进一步表明污染土壤优势茵群的群落结构信息是描述包气带土壤环境氨氮污染物迁移规律的重要参数。关键词：包气带；氨氮污染；微生物群落；变性梯度凝胶电泳；l6SrR、A中图分类号：S151.93；X820.4文献标识码：A包气带作

3、为全球生物地球化学循环的重要地质介质，对地下水生态系统影响显著，其理化因素的空间特异性能够使相应上壤环境中微生物群落发生改变，进而影响氨氮污染物在包气带土壤生态系统中的运移，以及氨氮污染降解菌的分布。近年来，现有的物理化学技术对处理大面积的土体和地下水污染已不具优势，污染土壤生态环境的研究重点已开始转向分子生物修复领域.采用PVA1799冷冻改良技术固定的优势菌群，能够使浓度为565mg/L的苯酚去除率达到94%以上。在硝化性能明显的氮污染土壤中分离得到的一种微生物Bacillussp.LY能够在有机物去除的同时,将氨氮转化为氮气，并能以多种有机污染物包括壬基酚聚氧乙烯醮（NPEOs）、邻苯二

4、酚及苯酚等为唯一碳源生长，是新型高效生物脱氮联合去除有机污染物技术的重要突破电。土壤环境中分离出具有特殊作用的微生物，一株鉴定为OchrobactrumintermediumDN2的烟碱降解细菌经证实，能够去除或减少烟草废弃物中烟碱”。在苯乙酸废水污染的高盐土壤分离出的中度嗜盐菌Halomotuissp.BYS-1中分离提取的耐盐基因片段，为耐盐遗传机制的研究提供了高质虽实验材料，在盐污染土壤治理中得到有效利用们。收稿日期：2012-03-07；接受日期：2012-07-24基金项目：国家自然科学基金项目面上基金（51078338）；中央高校基本科研业务费自由探索项目（2010ZY05）作者简

5、介：程雅楠.硕士研究生.从事地卜水及土壤污染微生物修复技术研究。E-mail：ptinny_5696163.conu,通讯作者：陈鸿汉，教授，主要从事地下水科学与环境工程教学与科研工作。E-mail：ehenhonghaneugb.«lu.en0反渗透等物理方法以及化学还原脱氮法是去除包气带土壤氨氮污染的传统方法。但是，这些方法并非完全将硝酸盐降解为无毒的气态貌化合物和氮气，仅为硝酸盐的转移，还会产生二次污染问题。近年来，生物脱觎技术尤其是分子生物技术在土壤氨氮污染去除中逐渐得到关注。变性梯度凝胶电泳（简称DGGE）技术能够直接利用DNA及DNA对微生物遗传特性进行表征，不但避免r耗

6、时的传统菌种分离，更可鉴定出无法利用传统方法分离的萌种，它的分辨精度相比琼脂糖凝胶电泳和聚丙烯酰胺凝胶电泳更加精确，可以检测到一个核甘酸水平的差异。DGGE胶体图谱中出现的可见条带，分别代表土壤样品中优势细菌16SrRNA的目的片段，条带的数鼠反映出细菌群落组成中的优势群体的数量，条带的强度代表细菌种群中细菌的含量，与模板25.StefanJG,OrnP,ThomasMG.Denitrifyingbacteriaisolatedfromterrestrialsubsurfacesedimentsexposedtomixed-wastecontamination_J_.AppliedatulEn

7、viron-menlalMicrobiology,2010(76):3244-3254.26ChristianeW,AndreaT,AntjeW.Horizon-specificbacterialcommunitycompositionofgennangrasslandsoils,asrevealedbypyrose(,uencing-basedanalysisofI6SrRNAgenesJ.AppliedandEnvironmentalMicrobiology,2010(76):6751-6759.MicrobialEcologicalDiversityCharacteristicsofth

8、eSoilProfileintheVadoseZonePollutedbyAmmoniaNitrogenCHENGYa-nan1'2,GUANXiang-yu1,QUWen-long3,CHENIlong-han*,LIUFei,XIEYu-xuan,ZHULing-ling(1.BeijingKeyLaboratoryofWaterResourcesandEnvironmentalEngineering,ChinaUniversityofGeosciences(Beijing),Beijing100083,China；2. RailwaySurveyandDesignInstitut

9、eofChinaCommunicationsCo.,LTD,Beijing100088,China；3. ShenyangWaterGroupCo.LTD,Shenyang110003,China)Abstract:Scholarsathomeandabroadinrecentyearshavefocusedonresearchingammonianitrogenpollutionmigrationrulesinthevadosezonesoil.Columnsimulationandsoftwaresimulationwereprimarilyconductedforthestudyofam

10、monianitrogenpollutionmigration.Biotechnologywaswidelyusedindegradationofpollutantinsoil,butlessappliedtothestudyofpollutionmigration.ThisstudyappliedDenaturingGradientGelElectrophoresis(DGGE)andthesequenceanalysisoftheV3RegionofI6SrRNAcombinedwithcanonicalcorrespondenceanalysistocharacterizethebact

11、eriaverticaldistributioncharaclerislicsandbacterialcommunitystructureinthesoilfromthreetypicalcontaminatedzonesintheNorthChinaPlain.Accordingtotheanalysisofphysicalandchemicalsoilpropertiesinpollutedareas,thereweresomedominantindividualbacteriainthekeypathwaysofthenitrogencycleandsulfatemetabolism.I

12、tissuggestedthatthebacterialcommunitiesareaffectedbythedistributionsofammonia,nitrateandnitritenitrogen,showingthatthecommunitystructureinfonnationofthedominantpopulationincontaminatedsoilisanimportantparametertostudytheammonianitrogenpollutionmigiationrules.Keywords:vadosezone；ammoniapollution；micr

13、obialcomniunities;DenaturingGradientGelElectrophoresis(DGGE)；16SrKNA27NicoleD,BrunoG,SteffenK.MethanolrophiccommunitiesinBrazilianFerralsolsfromnaturallyforested,afforested,andagriculturalsitesJ'.AppliedandEnvironinenial.Microbiology,2010(76)：1307-1310.DNA片段含量成正比。然而该技术的缺点在于,胶体条带的分辨率影响部分条带的生物信息，不

14、能获取到全部的样品微生物信息。此外16SrRNA克隆文库分析也是有效鉴定样品微生物分予生物技术,该方法能够详细地揭示样本中微生物的多样性,能够将样品微生物分类鉴定至种的水平。具体过程是将总DNA山样本中提取出，应用通用引物对其进行16SrRNA基因的扩增.进而建立基因克隆文库.再以限制性内切牺对其进行酶切自由，通过对其I6SrR、A基肉前切图潜的分析鉴定微生物的具体分类，得到微生物群落多样性的信息。针对DGGE技术有限分辨率造成的信息缺失以及克隆文库分析的覆盖范围局限，本文将DGGE技术与克隆序列结合，对微生物群落丰富度的差异进行深层研究，同时通过典范对应分析对污染源区微生物群落多样性及生态分

15、布特征进行评估，以期揭示氨氮污染源区包气带土壤理化性质对其细菌群落结构的塑造作用，为通过细菌生态系统变化反映污染物在包气带中的分布和迁移规律提供理论依据。1实验部分11仪器和装冒PCR扩增仪（CFX9，美国Bio-rad公司），凝胶成像分析系统（GelDocXR,美国Bio-rad公司），变性梯度凝胶电泳仪（DCode,美国Bio-rad公司），琼脂糖凝胶电泳仪（DYCP-44P,北京六一仪器生产公司），小型高速台式离心机（5418型，德国Eppdorf公司），高速冷冻离心机（3K3O,美国Sigma公司）,-80Y冰箱（美国ThennoElectronCorporation公司），高温高压灭

16、菌锅（LDZM-60KCS,上海医疗器械厂），旋涡混合器（Vortex-Genie2T,美国ScientificIndustries公司），洁净工作台（SW-CJ-1FD,苏州安泰），恒温摇床（THZ-82A型，苏州威尔），电热恒温水槽（DK8I）,上海风棱），电热恒温培养箱（D、P-9052,苏州威尔），电子分析天平（AI204,美国梅特勒公司），哈纳多参数测定仪（HI832O型，意大利Hanna公司）。12主要试剂去离子水（Millipore纯化柱），丙烯酰胺（纯度>99.0%,美国Sigma公司），双丙烯酰胺（纯度>99.9%,美国Promega公司），尿素（纯度N99.5%

17、,美国Sigma公司）,Tris碱（纯度'99.9%,美国Promega公司），EDTA二钠（EDTA-Na2-2H20,纯度"99.9%,美国Promega公司），去离子甲酰胺（纯度99.5%,美国Amresco公司），甲酰胺（纯度99.5%,美国Amresco公司），TEMED（纯度N99%，美国Sigma公司），过硫酸被（纯度N99.9%,美国Promega公司）,CTAB（生化级，纯度N99.0%,美国Sigma公司）。乙酸（分析纯，纯度N99.5%）,乙醇（分析纯,纯度N99.7%）,NaOH（分析纯，纯度N99.5%）,甲醛（分析纯，纯度299.5%）,NaCl（

18、分析纯，纯度99.5%）,Na2HPO4（分析纯，纯度N99.5%）,KH2PO4（分析纯，纯度N99.5%）:这些试剂均购自北京化学试剂公司。EB（10mg/mL母液，上海生工生物有限公司）,FastDNASpinKitforSoil（美国公司）,DNA胶网收试剂盒（北京天根生化科技有限公司）oTris-盐酸溶液（1mol/L,pH8.0）：称取Tris碱6.06g,加超纯水40mL溶解，滴加浓盐酸约2.1mL调pH至8.0,定容至50mLoEDTA（0.5mol/L,pH8.0）：称取9.306g的EDTA二钠，加超纯水35mL,剧烈搅拌，用约1gNaOH颗粒调pH至8.0,定容至50mL

19、。TAE（50x）缓冲液：Tris-HC12mol/L,CTAB50minol/L,EDTA（0.5mol/L,pH=8.0）20mmol/L。PBS缓冲液（0.01mol/L,pH=7.4）配方:NaCl8g/L,KC10.2g/L,Na2HPO41.44g/L,KH,P040.24g/L0TE（1x）缓冲液：1mol/LTris-HCl（pH8.0）1mL,0.5mmol/LEDTA（0.5mol/L,pH=8.0）0.2mL。超纯水定容至100mL013场地概况与样品采集采样区位于华北平原温榆河畔与北京城区之间，区域平均海拔约20m,上覆厚度为1015m黏土层，含水层为410m,具体由埋

20、深02m的杂性土壤层、埋深为210m的粉质黏土夹粉土层以及埋深为10-25m的粉砂-细沙层组成。取样区域单层含水层厚度一般小于5m,部分地区可大于5m,累计厚度在20m左右，局部地段大于30m,整体渗透系数在30-50m/d之间。土壤样品的采集点（T1、T2和T3）地处东经116°33北纬39。59,范围内,】'】位点位于靠近城区的生活潜在污染区,T2位点位于市郊果蔬种植场地,T3位点为靠近温榆河岸及垃圾填埋厂的污染区域。各站点土壤的野外采集完成后，迅速将样品统一置于无菌盒中，贴好标签密封，返回实验室后置于-70无超低温冰箱内长期保存，用于后续微生物计数和核酸提取。14样品环

21、境因子测定土壤样品pH值及电导率等参数利用哈纳多H18320参数测定仪进行现场测定，颗粒组成、土壤含水率、三氮含量及其他化学成分由后期实验测定,其他具体理化参数测定遵循表1所列的标准方法。表1上壤样品环境因子测定方法Table1Determinationofenvironmentalfactorsinsoilsamples测试项目测试方法测试项目测试方法土填含水址质量法态氮、总貌紫外分光光度法pH值电化学分析法亚硝态氮N-（l-恭基）-乙二胺光度法缜化还原电位电位法131囱字电感耦合等离子体颗粒组成比收计法I4J罔丁发射光谐法铉态氮纳氏比色法其他阴离r离子色谙法1.5荧光染料（DAPI）i，t

22、数取0.5g冻存土样和3.5mL无菌水，用匀浆器充分混匀，再加1.5mL含有0.05%TritonX-KX）的3xPBS缓冲液，750g离心1min（约2500r/min）,取1mL上清液，与3mL4%多聚甲醛均匀混合,4弋放置过夜以固定细胞。1mL细胞固定液中加染料DAPI至终浓度为200ng/mL,染色10min后用细幽滤器过滤，使菌体集中在0.22jun聚碳酸滤膜（Walemum滤膜）上,取下滤膜置于载玻片上迅速进行镜检;随机计数10个视野中的菌体个数取平均值，计算出每克土样中萌体细胞总数5。16样品中总DNA提取及细菌16SrRNA-V3叮变区片段的扩增取5g冻存样品装入10mL离心管

23、中，加入4mL裂解缓冲液3%CTABJOOmmol/LTris-盐酸（pH=8.0）,100mmol/LEDTA（pH8.0）,1.4mol/LNaCl,2%R筑基乙醇（V7U）,1%PVP,充分研磨,旋涡振荡仪剧烈振荡，尽虽使沉积物悬浮于裂解缓冲液中。加入蛋白丽K至终浓度1mg/mL,55勾温浴1h,其间每15min补充一次蛋白酶K。将裂解后的上清液分装到4个1.5mL离心管中，加入等体积的氯仿-异戊醇（体积比24：1）,以12000g转速4Y离心10min,抽提两次。取上清液加1/10体积3mol/L醋酸钠和预冷的等体积异丙醇，于-20幻放置1h后以12000g的转速4无离心29815mi

24、n,弃上清液，剩余沉淀用1mL70%乙醇洗涤一次，自然干燥，用50皿IxTE缓冲液（pH=8.0）溶解沉淀"门。利用16SrRNA-V3区引物341F/534R扩增样品DNA序列，上游引物341F包括一段5黏性末端连接GC夹子结构的前导区，具体序列成分见表2i。PCR反应体系为：DNA模板I皿；引物10pmol；TaqDNA聚合的1U；2.5mol/LdNTP,4p.L；5xPCR缓冲液，10piL；Mg2+溶液1.6皿；去离了水补至50puLoPCR反应条件为：95丁预变性4min，然后94丁变性30s,55幻退火30$和72X延伸1min,共35个循环，最后72T延伸10mino

25、PCR反应产物用1.2%琼脂糖凝胶电泳检测，得到的PCR产物利用产物纯化试剂盒（北京天根生化科技有限公司）进行纯化°表2引物序列Table2Primersequence引物序列34IF5'-CCTACCCCACGCAGCAG-3，534R5'-ATIACCGCGGCTCCTCC-3'34IF-GCS'一CCCCCGCCGCGCGCCGCGGGCGGGGCCCGGGCACGCCCCGCCTACCGCAGGCAGCAG-3'MI3-475'-TAATACGACTCACTATAGGGC-3，RV-M5*-ATTTAGGTGACACTATAGA

26、ATACTC-3'1.7细菌16SrRNAV3可变区片段的DGGE分肉与聚类分析将40jiL的16SrRNA-V3区PCR胶P）收产物，用变性梯度凝胶电泳仪（D-Codesystem,BIO-RAD,Hercules,CA,USA）进行DGGE分离,PAGE胶浓度为8%,变性梯度30%60%.60Y,80V电压，1xTAE（20mmol/LTris基质，10mmol/L醋酸钠,0.5mmol/LEDTA）缓冲液中运行13ho电泳后EB染色30min,置于凝胶成像仪紫外光下检测并照相。利用NTsys软件对DCGE凝胶电泳胶体上条带位置及条带亮度进行测定'，对DGGE条带分离出的优

27、势菌群进行LPGMA聚类及丰度分析，得到以皮尔森相似度为基础的聚类分析结果，并通过离差平方和法绘制聚类树状图。聚类树状图中，纵轴站点两个对象延长线交线垂直于横轴相似系数，对应横轴上的数据即为两者之间的相似度。同样,4个层位延长线与另外4个层位延长线的垂直交线对应数据.即为两站点间的相似程度。1.8细御16SrRNA-V3区序列的克隆和测序紫外灯下对成型的DGGE凝胶进行检测，将紫外光激发下可见的凝胶条带切下，用去离子水冲洗3次后，浸入50似，去离子水中,4丁过夜沉淀。取各样品浸出液2jxL为模板再次进行PCR扩增，二次PCR产物经胶回收试剂盒纯化后，连接到PMD18-T载体（TAKARA，大连

28、）上,经热激法转化到大肠杆菌TOP10细胞中，涂布含有Amp/X-Gal/IPTG的LB平板,37勾静置培养16ho然后4幻静置Ih,使蓝色充分显现后，随机调取白色函落到新的培养平板内过夜培养。随机挑取单个白色菌落，加入含100jiLddH2O的PCR管内.置于PCR仪99丁加热裂解细胞c取2皿裂解液作为模板，利用M13-47和RV-M为引物（具体序列列于表2）,进行后续的PCR扩增，验证插入片断的大小，筛选阳性克隆。将筛选好的阳性克隆保存于1mL含有20%H油的液体LB培养星中，每个平板随机选取10株已得到验证的阳性克隆进行测序（BGI,北京）o1.9典范对应分析（CCA）将三个站点土壤样品

29、的理化因子数据，包括pH值、Eh值、含水率、硝态氮含筮、亚硝态氮含虽、氨态敏含量、纽离子及氯离含量，与克隆文库测得的上壤细菌群落组成进行了典范对应分析,分为两个数据矩阵,即土壤理化因子数据矩阵以及土壤细菌排落组成数据矩阵。通过Canoe。软件计算分析出生态重要理化因了与土壤细菌群落分布之问的相丘影响关系,从而分析出单个环境因子对细弟群落组成的具体影响。1.10核酸序列世录号通过随机测序共获得107个16SrRNA-V3区序列，利用Check-Chimera检测（http:/cgis/llimera.cgi?su=SSL），将其中的嵌合体等不正常的序列剔除,81个

30、属于正常序列。本实验中采用的全部序列已提交到GenBank,序列号（AccessionNumber）为HQ916750HQ916804。111系统发育与统计分析获得的16SrRNA-V3区序列通过BLAST工具从GenBank获得与本实验所得到序列亲缘关系最近的类群序列。使用Thompson等"刁的方法对这些16SrRNA序列进行比对，然后利用MEGA3.0中的Neighbor-Joining模型进行系统发育与进化分析。并利用DOTUR软件（以DNA序列间距划分微生物分类及多样性的软件）版本1.53比较本实验三个站点共12个土层样品的细菌多样性,将相似性298%的细菌I6SrRNA基

31、因序列归为同一可操作分类单元（OUT）。2结果与讨论2.1士壤样品细曲总数统计及环境理化因素测定站点细电总数/酸碱度pH含水率/%代化还原电位Eli/mV干重mgkg'1)土峰岩性(CUFL)深度/mNH4-NNO,-Nno2-nSO?T1-1粉质黏土5.37x1091.948.7320.53247.002.240.460.0856.22Tl-2黏土夹砂7.27x10s3.278.7316.75201.000.520.470.1236.61T1-3粉细砂8.17xlO67.458.6318.34133.100.420.690.139.82Tl-4细砂2.64xI058.218.6817

32、.26107.200.531.960.126.62T2-1黏土8.87xlO90.958.1417.77566.0021.0422.2116.8819.04T2-2粉土1.53xl0R2.218.7117.75187.604.885.090.3524.38T2-3粉细砂8.31xlO64.008.333.83284.009.972.040.4349.83T2-4粉土9.53xIO68.158.569.83226.0066.5911.610.8133.2913-1黏上夹砂3.22xlO70.278.550.17491.000.464.561.34448.0513-2福土含砂1.54xlO81.2

33、38.6319.17372.00l.ll4.062.4729.10T3-3中细砂8.00xlO63.258.6218.0290.702.483.210.5945.75T3-4粗砂1.60xlO54.218.7316.8786.005.163.680.5228.00表3各点位土壤样品理化参数及微生物数量统计Table3physicochemicalfactorsandmicroorganismsofsoilsamples实验区三个站点土壤样品的pH值均为8.0以上;含水量在各点的变化较小，均在15%以上;氨氮含量及硝基氮含在T2站点积累较为明显,在站点T3积累其次，而T1站点三氮含量:检出较站点

34、T2与T3明显减少，亚硝氮在站点T3的表层位点含量较高，在T1站点及T2站点中没有明显积累;三个站点中硫酸盐含量不均，在6.62448.05mg/kg之间（表3）。各站点12个土壤样品中，土壤岩性细菌总数在1.6x1058.87xIO。范围之内。研究区剖面上壤样品各环境因素的测定结果表明，总氮及各类污染氮素的分布及含童在各采样站点差异明显，最高站点约是最低站点的60倍。取样位点的土壤环境中NO2-N、NO-N、Nf-N、so：-等阴离了的含鼠过高，表明研究区土壤处于氮素及无机盐类污染状态，地下水环境极有可能因此受到威胁。其中，研究区T1位点土壤的总微生物量随地层深度加深明显减少，说明细菌数量随

35、地层增加、有机质含量减少而减少。同时，表3数据还反映出，研究区内部水平方向微生物筋群分布基本均匀.但影响细菌总数的因素并不单-,研究区T2站点中深层位点（地下8m左右）与中层位点（地下3m左右）的细曲总量同属一个数量级，说明其他土壤理化因素的改变，如污染氮素的异常积累、土壤含水率及溶氧技等同样是影响土壤微生物总虽的关键因素。研究区站点同样存在特例,表层的细菌总数略低地下浅层的细菌总数，说明土壤样品中细菌的含虽并非随地层加深而增多，站点中土壤颗粒组成变化仍是影响细前总量的主要因素。2.2基于DGGF：的微生物多样性分析对研究区T1、T2及T3共12个土层进行变性梯度凝胶电泳（DGGE）的测定，结

36、果表明三个站点的上壤细菌群落具有较为明显的差异（如图1所示）。12个泳道中产生的条带均聚集在胶体的中间部位（该区域的变性浓度属于高低混合）。DGGE胶体图谱中共出现85条可见条带，并标注出35条亮度明显、位置清晰的DGGE胶体条带，代表土壤样品中优势细菌16SrRNA目的片段，进行后续序列分析。DGGE胶体图谱中，站点T1中表层土壤及浅层土壤的DGGE条带分布较为相似，各层位条带随深度增加亮度降低、种类也相应减少，部分菌群（条带1、3、4、6、8、9）表现尤为明显。结合表3中各站点土壤理化参数分析，包气带土壤剖面表层至底层之间微生物多样性发生的明显改变，与层位加深后上壤颗粒结构改变造成的有机质

37、积累变化相关。同时表现在DGGE图谱中，站点T2各层位的条带随层位加深变化并不明显，除部分菌群（条带17、18、19、21、24）为个别层位的特别优势菌群外，其他菌群（条带12、13、14、15、20）均贯穿于所有层位土壤之中，这与该位点各层位相似的土壤颗粒组成（表层至底层多以黏土及粉质黏土为主），及相似的污染物积累状态（氮污染积累较为明显）关系密切。站点T3在DGGE图谱中的条带分布表明,表300图1三个站点12个层面上壤样品16SrR、A-V3可变区的DGGE图潜Fig.ID（；（；Eprofilesof（he16SrDNA-V3regionof（hrtwelvelayersalongth

38、esoilsamplethreestations层及浅层I.壤的细菌群落分布相似性很高，条带2734所对应的细菌菌群在两层位中的含量存在差异，而细菌种类儿乎没有改变，但该层位见水较浅，且土壤颗粒变化较大，因此表层及浅层中大量存在的优势细菌群落在位处浅层他水带的13-3及T3-4层位中逐渐消亡，仅有少量优势菌群存在（条带29、30、31）并新增一种优势菌（条带35）,说明该位点土壤颗粒吸附细菌的能力随土壤剖面加深而减弱，同时较深层位污染物浓度的增加,也是造成土壤剖而细菌多样性分布差异明显的主要原因。使用软件NTsysVersion2.0131根据DGGE图谱中每条带所建立的矩阵进行聚类分析，通过

39、非加权组平均法（UPGMA）分析表明三个站点之间的群落结构相似程度均存在变化，聚类图如图2所示，横线链接的末端垂线所对应的数值即为两者的相似程度，如图中站点T1表层Tl-1与T1-3的微生物群落结构最相近，相似度为90%左右。说明T1-3中部分微生物种群在较深地层仍然存在。地下浅层T1-2与表层T1-1相比，微生物群落结构变化稍明显，相似度为86%。而随着地层的加深，底层T1-4的微生物群落分布明显不同于另外三层，相似度约只有75%。同样站点T2各层位的微生物群落结构也存在差异,T2-2层与其他三层差别最明显，相似度只有76%左右，而地层12-4与T2-I及T2-3两层的群落结构相似度则保持在

40、80%以上，其中表层T2-1与地下深层T2-3的微生物群落结构相似度最高，达到88%。丫-变形箫酸杆菌7%放线菌3%拟杆菌4%厚壁南4%芽布杆菌5%产黄杆陶2%a-变形菌7%&变形菌7%绿弯菌3%站点层位|T1-1'TI-3TI-2ITI-4|T3-1T3-2一|T3-3T3-4T2-1IT2-3T2-4T2-20.570.650.730.800.880.96相似系数图2DGGE分离所得细常I6SrRNA-V3|乂序列的NTSYS聚类图Fig.2NTSYSdendrogrambasedon16SrRNA-V3regionsequencesfromDGGE站点T3的各层位群落结构

41、差异较为明显，反映在图谱中的成层分布差异较大。表层T3-1与T3-2之间的相似度达到96%,而深层T3-3与T3-4之间的群落结构相似度同样为96%。但较浅两层与较深两层之间的差异明显，只达到80%,说明浅层与深层土壤剖而的微生物群落结构变化较大。反映在三个站点之间的群落结构相似性，站点T1与站点T3的群落结构较为接近，相似度达到75%,而站点T2与其他两个站点T1、T3差异较大，相似度只有57%。说明微生物菌群在站点T1及T2之间变化不大，但站点T3与站点T1及T2的细曲.群落结构分布差异均比较明显，原因在于站点T3紧靠温榆河畔，包气带地层见水层位较浅，土质成多为砂质土壤，造成站点T3的微生

42、物群落结构变化明显。2.3细菌群落结构及系统发育分析将DGGE分离出的优势细两条带（分别标记为条带135）在NCBI基因库（NCBI,/）中进行序列比对发现，实验区12个土壤样品中提取出的81个正常细菌16SrRNA序列共54个不同亲缘类型（OUT）,分属于11个系统发育门类。如图3所示，各分属具体细菌门类及所占细菌总亲缘类型比例分别为：a-变形细茜（7%）,8-变形细菌（10%）"-变形细菌（25%）,S-变形细菌（7%）壁前（4%），酸杆菌（7%）,放线菌（3%），绿弯菌（3%）,拟杆菌（4%）,产黄杆曲（2%），芽单胞幽（5%

43、）和未分类细菌（26%基于16SrDNA基因序列相似性298%的准则，0变形菌未分类细菌26%a-变形菌绿弯菌。酸杆荫8变形菌厚壁菌拟杆荫Y-变形菌产黄杆菌芽弛杆菌&变形菌放线维未分类细茜图3研究区三站点各层位总微生物胖落组成Fig.3Microbialcommunitycompositionofsoilhorizonsfromthreestationsinresearcharea选取DGGE图谱（图1）中得到的部分代表条带进行序列系统发育分析，结果表明细菌类群在不同站点和层位中的分布存在显著差异。如图4所示,T1站点4个层位样品中的优势细菌类群是y-变形细萌（代表条带3、7）,3-变

44、形细菌（代表条带2、8）-变形细菌（代表条带5）,绿弯菌（代表条带1），厚壁菌（代表条带4、9）和芽抱杆菌（代表条带6）,另外在该站点深层区域T1-4上壤样品中，检出少量产黄杆萌（代表条带10）和放线菌（代表条带11）。根据测序结果，发现该站点土壤总细弟中y-变形细菌在所占比例最大为36.4%,其次是厚壁菌和3-变形细菌（占到18.2%），绿弯菌占到13.6%，芽弛杆菌所占比例为9%,产黄杆笛及放线菌所占总菌类的比例仅为4.5%，其余菌种属于未分类细菌。随着深度的增加.底层T1-4中的优势细菌主要为产黄杆前。该站点由于表层覆盖建筑垃圾杂填土，浅层土壤颗粒为粉质黏土，深层主要由粉细砂、细砂组成。

45、地卜,污染氮素及以硫酸根高子为代表的无机盐在粉质黏土及黏土层积累明显，在以细砂为主的土层中含最迅速降低,y-变形细曲占绝对优势。另外厚壁街、5-变形细菌、芽弛杆菌和绿弯菌同样是表层的优势细萌群，其中厚壁菌类细南与来自植物根系土壤的细菌亲缘性较高），尤其发现一类异养脱氮细菌Bacillussp,菌属与氮循环密切相0000.«.、/、Z7ZVzz/-L80206040i11.z/A1蓍1A!11z.，/'/ii£|1£a1sl1§811*未分类细的芽抱杆情拟杆摘酸杆俩口放线菌H产黄杆菌厚壁的绿弯的方变形菌大变形前a缶变形菌甲变形菌T2.22-32-.

46、4T2T33T3-4n图4三个站点各层位细偶组成统计Fig.4Bacterialcom|w）sitionfromthreestationsinresearcharea关|叫。另有分属于5-变形菌的除硫单胞菌属（Desulfuromonas）被认为在硫循环中具有重要的作用“，同样仅存在于土质为粉质黏土和黏土夹细砂的Tl-1及T1-2层位中，说明相对砂质土壤硫酸根离子-蓄积量较高的黏七质土壤为是该类细菌.生存的首选环境,T2站点4个层位样品中的优势细菌类群是变形细菌（代表条带I4、23）,y-变形细|将（代表条带12、19、24、25）,a-变形细菌（代表条带20,21）.拟杆菌（代表条带15.2

47、2）,5-变形细菌（代表条带I6J7、】8）,产黄杆菌（代表条带26）和酸杆菌（代表条带13）o序列测试结果显示，各类细菌所占该站点总菌比例分别为/?-变形细菌（21.6%）、y-变形细菌（19.6%）、a-变形细菌（17.6%）、拟杆菌（17.6%）、5-变形细菌（5.9%）、酸杆菌（7.8%）和产黄杆菌（2%）.其余属于未分类细菌。该站点位于长期受到含氮类农药污染的城郊果蔬种植园,分析其土壤环境发现，该区域氨态氮积累极为严重，尤其在地卜浅层与饱和带之间积累更加显著。硝态氮及亚硝态氮含量随层位加深呈递减趋势,说明在这些层位的土壤微环境中发生了较为明显的反硝化作用，其中占绝对优势的细菌均属于反

48、硝化细菌类别，如假单胞菌和反硝化|W（ParacoccusfMintotrophusATCC35512）O此类细菌能够利用硝酸根、亚硝酸根、氮钗化物作为末端电子受体，以还原性硫化物为能源生长F-22.进而解释了站点T2中亚硝态瓶随层位增加的现象。另外，分属变形菌y亚纲的高度耐盐W盐单胞菌（Hulonunwssp.）生存于硫酸盐含原较高的T2站点粉质土壤中，与其亲缘关系最近的物种来自于苯乙酸废水的高盐污泥'，这类细菌的存在多数受污染盐类含量的影响。T2站点的各层位中还存在着各类有机污染的降解菌，如分属于3-变形细菌的地杆菌（Geobaeteraceae）、放线菌、变形菌以及拟杆菌等，这些

49、门类的细菌多来自混合污染土壤及淤泥中：24_251o但草地土壤街、农田）等环境的常见细菌绿弯菌、芽抱杆曲及厚壁菌未在T2站点中检出。T3站点4个层位样品中y-变形细菌（代表条带30、31、32、35）占较大优势，是总菌量的17.6%,6-变形细幽（代表条带33）愣-变形细曲.（代表条带28）、芽抱杆菌（代表条带29）、拟杆菌（代表条带34）和酸杆菌（代表条带27）所占比例均在8.0%左右，厚壁两及产黄杆菌所占比例仅为2.1%,其余曲种属于未分类细菌。毗邻河道的站点T3,与站点T1及T2相比，饱水带层位较浅，土壤密度较大,砂含量较高。由于主要颗粒幻成为粉土夹粉细砂的表层土壤T1，相比有黏土夹层的

50、浅层土壤12积累污染物质及吸附微生物的能力略差，剖面表层T3-1、T3-2层位的优势细菌群组成属于群落多样性较为丰富的一层。站位T3氨态氮含量由浅层到深层饱和层逐渐增加，黏质土壤层位T3-1中硝态知及亚硝态氮也有明显积累,微生物数量和种类较为丰富。但随着层位的加深,土壤受污染程度加大，包气带砂质土壤中微生物吸附能力衰减，使微生物尤其是优势硝化菌群减种类大幅下降，硝化反应程度减弱，氨态氮转化为硝态氮及亚硝态氮的程度也连锁性降低，造成细菌的群落多样性随着取样地层的加深而明显减少。2.4细慎类群的分布与环境因*的典范对应分析（CCA）各站点三氮虽、含水率、pH值、总氮、阴离子含珀等环境因子与克隆文库

51、细菌群落的典范对应分析结果显示（图5）,实验区各采样站点微生物组成在二维排序图上的分布与各种理化因子具有一定的联系。位于二维排序图左上方的物种有变形菌、放线萌、芽单胞菌和拟杆菌等,在实验区土壤中的分布与总氮含成呈明显正相关。放线菌及a-变形幽的分布与阴离子因素有一定的相关性，同时。-变形菌门与硝基氮含最显著正相关，而亚硝基氮含量及氨敏含量是影响排序图中右下方变形萌分布的关键参数。5-变形细菌和含水最具有显著的正相关性，与其他阴高于因素，如氟离子等，也有一'定的相关性，表明这类细菌主要分布在盐类含量较高的地区。图5左下方的绿弯菌、壁厚菌受含水量及pH值影响较大,这两种参数值的平稳与否决定

52、了两类萌群生长的状态。实验区存在的未分类菌群分布受pH值及含水量影响较大，同时与亚硝基氮及氨短含址呈负相关。图5实验区三位点中各细儡类群与环境关系的CCA二维排序图Fig.5CCAbiplotofsamplesandenvironmentalfactorsfromthreestationsofresearcharea图中各箭头标注代表的环境因子为:TN一总氮；Eh轼化还原电位;NO；-N暗基氮；NO；-N-亚硝基氮；NH；-N-氨氮；WC含水缺；pH-pH值；FCl"-MT；SOf-硫酸根离;NO；-无机盐硝酸根。三角图拆代表文库所包含的所有微生物种群分属的门类:Act放线菌；Un未

53、分类序列；gay-变形Bal拟杆菌；aia-变形菌；be尸-变形州deS-变形菌；Fir-壁厚菌；Chi一绿弯菌；Gem芽单胞菌；Fla-产黄杆爵3结语本研究样品所处的地下水系统中具有显著的空间特异性，地下水系统中微生物群落沿不同地层的成层分布特征与地球化学特征的成层分布相吻合。地下微环境中存在的微生物群落多为与氮循环有关的优势细菌个体，而总体群落结构分析表明该区域存在的微生物及其群落还参与着生态系统中硫酸盐等无机离子代谢的过程。研究结果表明研究区各站点土壤颗粒结构是影响包气带环境中多种污染物迁移转化及微生物群落多样性的关键，污染物在土层中的分布变化是细前群落结构更替的主要驱动力，而土壤中的微

54、生物也是污染物迁移过程中不可或缺的作用因素。通过对细菌群落结构分布与包y带土壤环境的地球化学性质分析，本研究揭示了地下水微环境中典型污染物迁移与细慎生态系统具有密切相关性，将为进一步研究包气带土壤理化参数对地下水环境的影响奠定基础。4参考文献:IWilliamJII.Vadosczonemierobialbioliarriersremovenitratefrompercolatinggroundwater'J.CurrentMicrobiology,2009(58)：622-627.'2MamieNI.KateMS,DennisER.Reductionofperrhloralt

55、*andnitratebymicrobialcommunitiesinvadosesoil:J.AppliedandEnvironmentalMicrobiology,2005,71(7)；3928-3934.3 卞华松，张仲燕.冷冻固定化优势萌群处理含甲醛苯酚废水J.环境科学,1998.19(2)：39-42.4 赵娟.吕剑，何义亮，靳强.张文英，何霞.异养脱敏菌株Bacillus$p.LY降解有毒有机污染物的研究J.环境科学.2007,28(12):2838-2842.5 袁勇军.陆兆新，黄丽金.吕凤霞.别小妹.烟碱降解细菌的分离鉴定及其降解性能的初步研究J.微生物学报,2005.45(2

56、)：181-184.6 洪青.张忠辉.张晓舟.徐剑宏.李顺胳.中度嗜盐偏Halomonassp.BYS21启动子的克隆和测序J.应用与环境生物学报,2005,11(6)：729-732.7 MyersBM,FischerSG,ManiatisT.M(HlificationofIhemeltingpropertiesofduplexDNAbyattachmentofaGC-richDNAsequenceasdeterminedbydenaturinggradientgelelectrophoresis_JJ.NucleicAcidsResearch,1985,13：3111-3129.I8Muy

57、zerG,BrinkhoffT,NubelU.Denaturinggradientgelelectrophoresis(DGGE)inmicrobialecologyJ.MolecularMicrobialEcologyManual,1998,3(44)：I-27.91ZweifelUL,HagstrornA.Totalcountsofmarinebaderiaincludealargefra(*tionofnon-nucleoid-containingIwiclrria(ghosts)J1.AppliedandEnviron-mentalMicrobiology,1995.61(6)：2180-2185.10SchallenbeiM,Kalfl*J,RasmussenJB.Solutionstoproblemsinenumeratingsedimentbacteriabydirectcount.J.Applie