网络出版地址:http://www.cnki.net/kcms/detail/11.2031.Q.20160413.1055.010.html
第36卷第24期2016年12月
ACTAECOLOGICASINICA
生态学报
Vol.36,No.24Dec.,2016
DOI:10.5846/stxb201506041132
姬洪飞,王颖.分子生物学方法在环境微生物生态学中的应用研究进展.生态学报,2016,36(24): ⁃ .
JiHF,WangY.Advancesinmolecularapproachapplicationsinmicrobialecologystudies.ActaEcologicaSinica,2016,36(24): ⁃ .
分子生物学方法在环境微生物生态学中的应用研究进展
姬洪飞1,2,王 颖1,2,∗
1西北农林科技大学水土保持研究所黄土高原土壤侵蚀与旱地农业国家重点实验室,杨凌 7121002中科院水利部水土保持研究所,杨凌 712100
摘要:随着分子生物学方法的不断发展和改进,微生物在生态系统中的作用被更好的挖掘出来。目前快速发展的先进的分子生物学技术,已经开始应用于分析环境微生物的多样性、微生物的生物地理学及微生物对气候变化的响应等。一般环境微生物的研究目标主要有3个,即确定微生物的种类和多样性、微生物的功能或潜在作用及在特定时间点活跃的微生物等。然而,现有微生物的研究方法复杂多样,容易给研究者在方法的选择上带来困惑。本文将从微生物的多样性和功能研究两个方面介绍和分析相应的分子生物学方法,尤其是近年来快速发展的高通量测序、宏组学和单细胞水平研究方法(如纳米二次离子质谱与荧光原位杂交相结合的方法)等新技术及其应用情况,以期为研究者选择合适的研究方法进行环境微生物的研究提供依据。关键词:微生物多样性;微生物功能;同位素标记;高通量测序;宏基因组/宏转录组;单细胞水平研究方法Advancesinmolecularapproachapplicationsinmicrobialecologystudies
1StateKeyLaboratoryofSoilErosionandDrylandFarmingonLoessPlateau,InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling2InstituteofSoilandWaterConservation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling712100,China712100China
JIHongfei1,2,WANGYing1,2,∗
Abstract:Microorganismsinhabitalmosteveryimaginableenvironmentandareimportantplayersinbiogeochemicalcycles.Withthedevelopmentandimprovementofmolecularbiology⁃basedapproaches,thefunctionsandcommunitycompositionofbeenappliedtoquestionsregardingmicrobialdiversity,biogeography,andresponsestoenvironmentalchanges.Studiesofmicroorganismsintheenvironmentgenerallyfocusonthreeobjectives—determiningwhichmicroorganismsarepresent,whatmicroorganismsinecosystemscanbeexploredinmoredetailthaneverbefore.Rapidlyadvancingmoleculartechniqueshave
theirfunctionsare,andwhichareactiveatagiventime.Comprehendingtherangeoftechniquescurrentlyavailablecanbedaunting.Tofacilitatetheselectionoftheappropriateapproachbyresearcherstostudymicrobialcommunitiesintheenvironment,weintroducemolecularmethodsaccordingtomicrobialdiversityandfunction,suchastherapidlydevelopingspectrometry⁃fluorescenceinsituhybridization,NanoSIMS⁃FISH),andtheirapplicationsinmicrobialecologystudies.high⁃throughputsequencing,meta⁃omics,andmicrobialsingle⁃cellapproaches(e.g.,nano⁃scalesecondaryionmass
KeyWords:microbialdiversity;microbialfunction;isotopelabeling;high⁃throughputsequencing;metagenomics/metatranscriptomics;single⁃cellapproaches
基金项目:国家自然科学基金项目(41301322);黄土高原土壤侵蚀与旱地农业国家重点实验室基金项目(K318009902⁃1303);陕西省自然科学基金(2015JQ4111)
收稿日期:2015⁃06⁃04; 网络出版日期:2015⁃00⁃00∗通讯作者Correspondingauthor.E⁃mail:yingwang@nwsuaf.edu.cn
http://www.ecologica.cn
2 生 态 学 报 36卷
微生物包括细菌、古菌、真菌等,是生态系统的重要组成部分,几乎存在于所有已知的环境中,在生态系统中起着非常重要的作用,如维持生态系统功能、直接参与碳、氮、硫、磷和一些金属的生物地球化学循环[1]、加快污染环境修复进程等[2]。此外,微生物在生态系统对气候变化的响应方面也可能起着重要作用[3]。基于渐浓。
目前环境微生物生态学的研究主要集中在确定所研究的环境中有哪些微生物存在和分析微生物在所研究的环境中起什么样的作用或执行哪个生态功能[4]两个方面,此外,微生物之间以及微生物与环境因子之间的相互作用受到越来越多的关注。由于微生物个体小(平均直径约1μm),肉眼不可见且未知,微生物的研究方法的发展加速了研究者对微生物的认知。在过去的十多年中,微生物的研究方法得到了快速发展,且商业化迅速,特别是高通量测序、宏组学和单细胞水平研究方法(如纳米二次离子质谱与荧光原位杂交相结合的方法)等(图1)技术手段使微生物研究的难度和成本都得到了显著改善。随着研究方法的发展,微生物在生态系统中的作用可以被更好的挖掘出来。然而,快速发展的微生物研究方法,以及之前已有的研究方法,给不熟悉这些技术的研究者在方法的选择上带来了困惑。因此,本文将从微生物的多样性和功能分析两个方面对相应的研究方法进行阐述,并以1—2篇文献为例具体介绍相关方法在微生物生态学研究中的应用,以期为研究者选择合适的研究方法进行环境微生物的研究提供参考。1 环境微生物的多样性分析
在很多微生物生态学的研究中,首先是分析所研究的环境中有哪些微生物,以及微生物的多样性和分布情况。微生物研究的传统方法首先是获取环境样品,然后对其进行分离培养得到微生物的单菌株。然而,很多微生物难以通过现有的培养方法得到[5],使微生物的多样性分析受到很大。与传统分离培养方法不同,现代分子生物学技术直接以从环境样品中提取的核酸(DNA和RNA)为分析对象,以特定的核酸片段为生物标志物,理解该环境中微生物的组成和多样性,即将所提取的DNA进行目标片段的PCR扩增,得到的PCR产物通过多种技术来鉴定其组成的多态性,如群落结构分析的指纹图谱、微列阵(Microarray)和高通量测序等。微生物多样性分析的研究中,PCR扩增的目标片段可以选择微生物的一段相对保守而又足以区分不同类群的基因区域,如细菌和古菌常用的16SrRNA[6]、真菌的18SrRNA或ITS区[7]等;同时,一些微生物的功能基因也可用于分析功能微生物的多样性,如氨氧化细菌的氨单加氧酶amoA基因(氨单加氧酶,催化硝化作用的第一步)。得到目标基因的PCR产物后,可根据需求选择下面的方法进行微生物的多样性分析1.1 变性梯度凝胶电泳和末端性片段长度多态性变性梯度凝胶电泳(DGGE)和末端性片段长度多态性(T⁃RFLP)均是通过电泳对PCR产物进行分型,从而分析微生物的多样性和相对丰富度,在分析环境中优势微生物群落(相对丰富度>1%的微生物种群)[8]的变化方面具有较多的应用(表1)。DGGE和T⁃RFLP是应用于环境微生物多样性分析中较早的分子生物学方法,不仅应用方便,分辨率高,且价格低廉。DGGE主要用于分析低于500bp的片段;T⁃RFLP对低于500bp的酶切片段分型精度较高,而高于500bp的酶切片段分型精度则不足。DGGE条带可通过切胶回收、构建克隆文库和测序来了解感兴趣的条带的系统发育信息或物种信息,而对于复杂的微生物群落组成(如土壤细菌的群落组成)的分析,T⁃RFLP的分辨率高于DGGE。1.2 微列阵(Microarray)
基因微列阵用来分析样品中微生物的组成和多样性,已经用了将近20年,经历了数代的改进,常用的微列阵有系统发育芯片(PhyloChips,用于识别微生物以及微生物之间的系统发育联系,分析微生物的多样性)和功能基因芯片(GeoChips,用于研究功能基因的多样性和功能微生物的活性)。微列阵即为排满探针(用于和样品DNA进行杂交的已知短序列)的芯片,通过探针提供生物的系统发育信息或功能性质信息、或者二者
http://www.ecologica.cn
(图1)。
微生物在生态系统中所发挥的重要作用,环境和生态学研究的学者对理解生态系统中微生物的动态变化兴趣
24期 姬洪飞 等:分子生物学方法在环境微生物生态学中的应用研究进展 3
图1 微生物生态学研究中分子生物学方法的应用选择图[4]
Fig.1 Decisiondiagramforchoosingamolecularapproachforuseinmicrobiologystudies
均包含。当样品中的序列(纯菌、环境样品DNA或RNA的目标片段的荧光标记PCR产物或随机引物的荧光标记PCR产物)与探针进行杂交后,即可计算与探针匹配序列的相对荧光比例,从而获得微生物的多样性和相对丰富度信息。
基因微列阵的方法尤其适合识别不同时间、地点和处理之间的代表微生物或微生物群落的差别,此外,功能基因芯片还可以定量C、N、S和P循环、有机污染物降解和胁迫响应等相关的功能基因的变化。Aronson等[19]利用3代PhyloChip(16SrRNA基因微列阵,可以提供60000个左右不同OTU的信息),研究不同实验条件下松林土壤甲烷循环的相关微生物,分析了土壤中甲烷氧化菌(Methanotrophs)和产甲烷菌(Methanogens)的多样性及其与土壤中甲烷通量变化之间的联系。该方法不仅检测到了目标群落在不同时间、地点和氮水平处理之间的差异,同时发现了群落中占少数的微生物的变化,而在克隆筛选的方法或深度不够的高通量测序中,这些少数微生物的差异均可能被忽略。Cong等[20]利用GeoChip5.0(覆盖了393个功能基因家族)分析了
http://www.ecologica.cn
4
segatnav点da缺siDydgeozlyolaibnoarsc况ime情ml标tif用atn指o应e析sep的mn分yT中ori学vn态e生ni物se生hc微ao境rp环pa在ras法lue方cel类pyt学o种el物mf品p生os样maS子no分itdaecs iu1lp表pAnoi st1e术eaunilb技qbiaT用nmhocc联eTnino率itu辨los分eR法方d用oht应eM生 态 学 报
36卷
兴⁃T,在要但偏次测;]由已得量高要镜好发采卤改胞感;种,]需多在31完,[与获求是需微更自过用来细对9物PL[,列的]的序差成列列难需不,来的通采]71物,个F0序来种证一R⁃1T[条带基测偏组,序序很)显记Sm合品可物验表,移万火碱an生短条的,或[)生结样,]6微等个来代估飘百退一im产更整到时标μ子IM相H号1合[一序定低的至段同ul者得l易得对素S0电o5n等SI信止测一的定甚片的I;容或获,比位a—描N1扫F号结响]81不隆不一千中中]2基楚行同;(、法光。信相影[C]1时]能克段性现几阵11列[碱使清进谱的41小)方]荧5光谱会21[1[可片样多出定测列确序量域个即;不区图移%较ME镜[的荧质能的带行FR落会确检及准条定区百难还能确偏m域T(o微标针强子可离新条进-能来时以不一确列几困能功准谱ta区显目探增也次入程引个要T群可合结同计果对准序为定功在的光0测镜1微线测盖检掩来二过如一需个物中相以设结,到的多一的潜能到定或H米备,S的带每生中做高列有因或功测到检显射锁响IF纳制成的微程库可的使序难量序因基因到检达且,子X并影⁃D与的组E条过隆阵针可GP对测很含的基多基列n需术电或,技射)择易R(品素G的LF成泳克列探也4复CG到的很的序am高量透ME选光A针C探样位;同D趣R造电与微是差54重定得整于知由aR较通与S(的荧用素善的如化功菌细群的N微细,9-1微菌氧及化化的因、成真化的跃C组、境,C,的能烷以氧氧菌基、氧等活环像2-1功甲菌的菌氨菌中量菌氨烷含等成(N关和真因细如化境数真烷菌甲的及、包息光像5相1C环菌、甲基的,信和菌和量能产度应以菌的跃菌氧环菌化质能荧成21循,成细菌数功和富响细应活能烷或古硝性功测素,N-的化和及菌丰迫的响中功甲和量落、反组的和检等位N和跃氧成组以化对胁度菌、落跃迫境落类号菌4活氨落菌氧相和富真菌群活环及菌数群胁度或以化和个群种信解同1CC、2中如群古烷和环丰、化的、成循和中和富中菌硝成整物的n降和1菌细菌氧等境菌的菌甲、组N成境环丰境古反组中、、落境生物am奈素,-,等)物细氨菌环能等细菌落、C组环循和环菌菌群环微生a,R菌元C3-1F生品品品品品品品微养样样样样样样样的培境境境境境境境中体环环环环环环环物水等等等等等等等养体体体体体体体培或水水水水水水水体体、、、、、、、物AA物A物AA物A物A水水胞、细积N物N或物物R积R积NR积N物NR积R积NR积NR沉体积生//A沉A沉/A沉//A沉A沉/A沉/A水胞沉微、、、、、、、细、的壤N、ND壤D壤ND壤ND壤ND壤ND壤ND菌物菌中土的土的土的土的土的土的土的纯生纯物定定,4,库,,记54HS隆记记记记IH记克标库标标标/序标,标,素隆R素素素素H素pC位克,ihp位位位测序位SI位P同RCi同F⁃量定CPoh同R同a测nlC同yo定D同R定稳量he定定CP定i酸m磷PG稳稳量稳ullI焦稳A定C稳富群丰物对相生(微平平平的平平水水水%平水水胞胞落1水>)落落落细细群度落群群群单单性P杂制LF位限R⁃g原端Tn/i末EcGn交光e/Gu杂荧泳Dq组e位⁃电性s录sci原谱Ht胶态yaurp转m/光质SI凝多rao序hg宏ostcpii荧子F⁃H⁃离S度度rci测u/ror组mcos谱SIF次IMS梯长M量ht光⁃n二ona性段列通⁃因nnearhg基gattat曼am米N变片微高iH宏eMem拉aR纳交http://www.ecologica.cn
24期 姬洪飞 等:分子生物学方法在环境微生物生态学中的应用研究进展 5
热带雨林土壤微生物的功能基因的组成和多样性,发现所检测到的基因类型与生物地球化学过程有关,如碳、氮和磷循环相关的功能基因,其中,碳循环相关基因包括易降解和难降解底物的相关微生物基因,在三个取样地点之间显著不同,典范对应分析(canonicalcorrespondenceanalysis)和多元回归树分析(multivariate1.3 高通量测序
regressiontreeanalysis)均显示土壤有效氮含量与土壤微生物功能基因的组成和代谢潜力显著相关。
Sanger测序已经用了几十年,测序质量高,序列长度达750—1000bp。该方法无法测混合序列,需首先建
立克隆文库,即将目标序列转移到寄主细胞(一般为Escherichiacoli)中进行培养,形成单菌株,然后对单菌株中的目标序列进行测序,在分析微生物多样性较高的环境样品时,费时且费用高。
在过去十年中,测序方法得到了很大改进,在Sanger测序方法的基础上,用不同颜色的荧光标记四种不同的dNTP,当DNA聚合酶合成互补链时,每添加一种dNTP就会释放出相应的荧光,捕捉到的荧光信号经过特定的计算机软件处理,从而获得待测DNA的序列信息,即边合成边测序,也称为二代测序。该方法以微生物目标基因的PCR产物为样本进行测序,一个反应得到几万至几百万条序列,测序的广度和深度大大提高。目前用于环境样品中微生物多样性分析的常用方法有2种,即RocheInc公司的454焦磷酸测序和IlluminaInc公司的MiSeq/HiSeq测序,这两种高通量测序方法通过在每个样品的引物上加标签来识别不同样品的序FluidigmCorporation,SouthSanFrancisco,CA,USA)通过使用阵列芯片,可同时进行2304个PCR的扩增,将样品定量PCR和序列文库PCR在一个反应中同时完成,该PCR产物可直接用于二代测序分析,仅需5h,与二代测序相结合,可同时得到微生物的群落组成和丰富度信息,简化二代测序样本的制备程序[21]。目前,很多研究中应用的测序方法依然是454焦磷酸测序。然而,Roche于2013年10月声明454焦磷酸测序将在2015年底停止运营,但现有的454测序仪还将继续运行。与454焦磷酸测序相比,Illumina平台测序的通量大大提高,序列长度不断改进,费用也有所降低。因此,Illumina平台测序以及其他一些测序方法(如Helicons、PacificBiotechnology和OxfordNanoporeSystem的单分子测序)在今后的环境微生物研究中将会有更多的应用。
二代测序的测序广度和深度较克隆文库的方法大大提高,尤其是提高了对数量上占少数的微生物群落的覆盖,不仅能从群落水平上揭示微生物群落组成的变化,还能从更细的微生物分类水平上显示微生物群落的具体变化,在环境样品的16SrRNA、18SrRNA、真菌的ITS区和功能基因的分析中均有应用[6⁃7,22⁃23]。寒带针叶林在全球的占地面积约为11%,占土壤有机碳储量的16%,一般认为地上部的植物残体是土壤有机碳累积的主要来源。而Clemmensen等[24]等通过分析北方森林土壤年代序列的有机碳组成发现,50%—70%的土壤有机碳储量源于根系和根系伴生微生物,且稳定同位素标记与真菌ITS区的454焦磷酸测序相结合的研究结果表明,在不同地点之间,有机碳累积差异最大的土壤剖面中,根系伴生真菌均占优势,说明根系伴生真菌在生态系统的碳动态中起重要的调节作用。Hermann⁃bank[21]将高通量列阵与二代测序的方法相结合,分析了肠道微生物的丰富度和组成,结果表明,不同肠道位置或腹泻情况下的肠道微生物组成和丰富度不同,腹泻与链球菌属(Streptococcus)成员的数量显著降低相关,尤其是非解乳链球菌(S.alactolyticus)。尽管目前高通量列阵与二代测序相结合的方法的应用较少,该研究依然显示了高通量列阵在定量不同样品中微生物的丰富度的同时,通过与二代测序相结合具有快速得到微生物系统发育信息的潜力。2 环境微生物的功能性质和基因表达分析
除了理解微生物的多样性和群落组成之外,很多时候还需要分析这些微生物成员的功能,如微生物的多样性与其功能性质有什么样的联系,以及微生物与其所在生态系统中的能量和养分流动是如何相互作用的等。微生物功能分析常用的方法包括稳定同位素标记与微生物DNA或RNA分析相结合的方法、功能基因微列阵(见1.2)、微生物宏基因组和宏转录组分析及单细胞水平研究方法来识别目标微生物的种类和功能的
http://www.ecologica.cn
列,均可以同时分析多个环境样品,并得到大量序列。此外,高通量列阵(TheAccessArray48.48,AA48.48,
6 生 态 学 报 36卷
方法。
2.1 稳定同位素标记与DNA或RNA相结合
将水体、土壤、沉积物以及其他生态系统中的微生物的种类和功能联系起来一直是微生物生态学研究的主要目标之一。过去二十几年的研究中,环境微生物的种类和功能研究已经得到了大量结果,而二者的直接联系,仍然是个薄弱环节[25]。稳定同位素和放射性同位素技术将不可培养的微生物种类和代谢功能相联系起来,是微生物生态学发展中的重要进步。稳定同位素探针技术(Stableisotopeprobing,SIP))采用稳定同位素(如13C,15N和18O等)标记的底物[26⁃28],使利用这些底物的微生物的标志物(DNA、RNA、蛋白质或磷脂脂肪酸(PLFA))被标记(表2),从而从群落水平识别环境中活跃的微生物以及微生物之间的相互作用。
表2 DNA⁃SIP、RNA⁃SIP、PLFA⁃SIP和protein⁃SIP方法方面的比较[2]
Table2 ComparisonofmethodologicalconsiderationsforDNA⁃SIP,RNA⁃SIP,PLFA⁃SIPandprotein⁃SIP
性质Trait
敏感性Sensitivity培养时间
Incubationtime分类学分辨率
Taxonomicresolution代谢活性预测
Indicationofmetabolicability提取容易程度Easeofextraction稳定性Stability
应用比较Comparisonof
biomarkers′applicabilityDNA<RNA<PLFA<蛋白质DNA>RNA>PLFA>蛋白质DNA≈RNA>蛋白质>PLFADNA<RNA<蛋白质DNA≈PLFA>RNA>蛋白质DNA≈PLFA>蛋白质>RNA解析
Explanation
DNA⁃SIP需要同位素的标记丰度较高,13C达到15—20%,15N达到30%左右,而蛋白质标记只需1%的标记丰度;RNA的标记速度是DNA的6.5倍
培养时间与敏感性直接相关,DNA⁃SIP需要的培养时间最长,也更容易带来微生物交叉取食的偏差PLFA⁃SIP只区分较大的分类群,DNA或RNA⁃SIP可提供属水平或者更低水平的分类信息。与16SrRNA相比,蛋白序列的数据库现有信息受到很大程度的蛋白质是代谢活性最明确的指示指标,DNA只显示代谢潜力
DNA和PLFA在不同介质中均可以按照步骤提取获得,而RNA和蛋白质在环境样品中的提取有一定困难DNA和PLFA相对稳定,蛋白质可能变性,mRNA非常容易降解
唯一可以将微生物种类和代谢功能直接联系起来的方法,还可用于分析植物和微生物之间的相互作用。其中,DNA⁃SIP和RNA⁃SIP均可将微生物的系统发育信息或功能基因信息与其代谢过程相联系起来,操作较容易,且分辨率高,在环境样品,尤其是较复杂的环境样品的分析中(如土壤)应用广泛,并常与其他DNA或RNA分析的方法相结合来分析环境样品中活跃代谢标记底物的微生物群落组成,如荧光实时定量PCR、DGGE[22]、T⁃RFLP[23]、微列阵[29]和高通量测序(目标基因或宏组学)[22⁃23,30]等(表1)。缺氧的水稻土壤微生物对甲烷排放的贡献占全球甲烷排放的10%—25%,而甲烷产生的主要碳源是植物残体和根系分泌物,识别活跃代谢植物残体和根系分泌物且产生甲烷的微生物种群具有重要的环境意义。Lu等[31]通过采用13CO2代替空气中的CO2对水稻进行脉冲标记,分析根际土壤活跃代谢根系分泌物的微生物,发现根际土壤中古菌RiceClusterI显著同化了13C,说明该类古菌在植物分泌物或残体分解产生甲烷的过程中起着重要的作用。DNA)相比,重DNA组分(理论上为13C⁃DNA)中酶切片段为485bp的峰占优势,系统发育上属于硫单胞杆菌属(Desulfuromonas);通过宏基因组学的方法,发现该菌除了具有严格厌氧菌的II型苯甲酰辅酶A还原酶之外,还具有兼性厌氧菌的I⁃型苯甲酰辅酶A还原酶,另外,还识别出了各种厌氧电子受体和编码末端氧化酶的相关基因,从而推断出该属微生物代谢甲苯的潜力。2.2 宏基因组学和宏转录组学
微生物基因组和转录组,即提取一种已知生物的DNA或RNA,不进行目标片段的PCR扩增,随机打断成几百个碱基的片段进行测序,得到全面的遗传信息,也称为鸟测序法。与基因组学和转录组学不同的是,宏
http://www.ecologica.cn
SIP实验虽然只能应用于富集培养物或环境群落的密闭培养实验(bottleexperiment),但是该方法是目前
Kim等[23]采用13C⁃甲苯标记海洋沉积物中微生物,通过DNA⁃T⁃RFLP分析发现,与轻DNA组分(理论上为12C⁃
24期 姬洪飞 等:分子生物学方法在环境微生物生态学中的应用研究进展 7
基因组学和宏转录组学反映的是多种生物或环境中整个群落的特征,包含环境微生物的全部遗传信息,相比于16SrRNA,除了群落中各种微生物群落的分类信息外,宏组学更包含了所有微生物的基因信息,有助于对微生物群落的潜在功能进行深入分析[32]。宏基因组学的挑战有2个方面,一是从随机测序的短序列中识别和定量基因功能,二是将这些功能基因片段与其他可以作为分类信息参考的基因联系起来,一般功能基因主要通过将测得的短序列与已经发表的或已注释的基因组进行比对来识别[4],与已有的参考数据库信息直接相关。宏转录组学和宏基因组学的区别在于,宏转录组学可以实时反映微生物群落的基因表达情况,得到的信息包括样品中哪些生物在当前是活跃的,以及它们在做什么。微生物细胞和环境样品中(如土壤)RNA的滞留时间显著低于DNA,因此,环境样品中提取的RNA理论上是该样品当前生物基因表达的一个快照。深度宏组学研究常用来理解复杂微生物群落的动态,然而,这类项目(尤其是在土壤中)所需的计算资源依然受到,如计算机硬件和数据处理方面。目前,虽然计算所需硬件的仍然存在,但是数据处理方面发展迅速,如WebCarma[33],MG⁃RAST[34],Galaxy[35]和GenePattern[36]等工具已经得到了广泛的应用。息;与稳定同位素标记相结合,较其他微生物多样性分析的方法,宏组学可以更全面的反应代谢特定底物的微生物的信息。随着测序技术的发展,宏组学的方法在环境微生物的研究中已有较多的应用,包括土壤、湖泊和海洋等生态系统[23,30,37]。Xia等[38]采用宏基因组学和宏转录组学的方法分析了高温下纤维素厌氧降解产生甲烷的过程的微生物群落,结果表明,在纤维素高温降解过程中,热袍菌(Thermotogales)对β⁃糖的消耗是纤维素水解过程中的关键步骤;尽管纤维素发酵过程产生大量的乙酸,嗜乙酸甲烷八叠球菌(Methanosarcinales)的丰富度较嗜氢甲烷杆菌(Methanobacteriales)低60%,而关键基因的转录活性结果表明,乙酸降解产甲烷途径的活跃程度显著高于嗜氢产甲烷途径,说明嗜乙酸产甲烷菌在数量上占少数并不代表其代谢活性低。
2.3 单细胞水平研究方法识别目标微生物的种类和功能稳定同位素与微生物的DNA或RNA(DNA/RNA⁃SIP)相结合,可以从群落水平将微生物的种类、组成与其功能相联系起来,但是对微生物单细胞水平的代谢活性以及不同种类微生物细胞的代谢活性的差异等信息依然难以揭示。随着单细胞成像方法的改进,结合稳定同位素或放射性同位素标记和荧光原位杂交(Fluorescenceinsituhybridization,FISH),可将微生物的数量、种类和功能以及一些微生物之间的相互作用从视觉上展现在图片中。放射性自显影技术(Microautoradiography,MAR)是最早的微生物单细胞活性的测定方法之一,灵敏度高,放射性底物的检测量最低可达10-10—10-15mg[39],然而,该方法的应用受到标记元素的限Raman)和纳米二次离子质谱(Nano⁃scalesecondaryionmassspectrometry,nanoSIMS)技术的发展,使单细胞成像的方法得到补充。
拉曼光谱技术是根据原子质量高的重同位素(如13C和15N)产生的拉曼光谱波长较轻同位素短,来分析制,标记元素需要同时具有放射性和合适的半衰期两个条件。拉曼光谱学显微技术(Ramanspectroscopy,
宏基因组学和宏转录组学呈现的是多种生物或环境中整个群落的特征,包含环境微生物的全部遗传信微生物细胞的代谢功能。例如,当微生物利用13C标记的底物以后,一些生物标志物的拉曼波长变短,如核酸、蛋白质、碳水化合物和脂类等[14]。此外,尽管每个物种均可以检测到核酸、蛋白质、碳水化合物和脂类等生物标志物,但是不同物种所产生图谱的峰高有一定差异,因此,每个物种产生的拉曼图谱不同,说明拉曼光谱技术还具有生成全细胞指纹图谱(whole⁃organismfingerprints)的潜力[40]。
纳米二次离子质谱可以检测任何稳定性同位素和具有合适半衰期的任何放射性同位素,常用于分析生物
样品中的C、N、P、S、O、H和卤族元素等[41],得到复杂微生物群落中微生物细胞代谢活性的成像。此外,测定同位素的比值信息(例如13C/12C),可以用于计算单细胞的元素吸收速率,分析元素的代谢途径[25],为环境中微生物的单细胞生态学和代谢潜力的研究打开了一个前所未有的窗口。该方法对同位素比值的分析具有非常高的灵敏度和准确度(理论上为1‰)[25,42],与其他单细胞的研究方法相比,具有非常高的优势。另外,纳米二次离子质谱与荧光原位杂交的方法相结合时,探针中的荧光标记由卤素代替,即可在同一张图像中同时
http://www.ecologica.cn
8 生 态 学 报 36卷
显示微生物的种类和功能[17],可以避免将荧光原位杂交的图像与NanoSIMS分析的图像难以正确联系起来的问题。
拉曼光谱技术和纳米二次离子质谱的分辨率可达微米甚至纳米级别,通过放射性或稳定性同位素底物标记分析纯培养或环境微生物细胞的代谢活性,结合原位杂交来确定目标微生物的系统发育信息[43],尽管这些方法目前主要应用于相对简单的环境,如水体、纯培养微生物或富集培养物等[43⁃44],在土壤微生物的研究中的应用还需要更多的探索和优化,这些方法的应用增强了对环境微生物种类和功能之间联系的理解,同时为将微生物如动植物一样在单个生物水平上进行研究提供了契机。Raman⁃FISH作为一种单细胞水平非破坏性的方法,在微生物生态学研究中的应用还较少。到目前为止,只有13C标记的研究结果。通过13C10⁃萘标记结果表明,地表水中绿脓杆菌(Pseudomonasspp.)可以显著吸收同化13C10⁃萘,且单个绿脓杆菌细胞之间13C的含NanoSIMS⁃FISH在微生物生态学中的应用较Raman⁃FISH广泛,在微生物碳、氮和硫循环的研究中均有应用[43,46⁃47]。Musat等[47]研究寡营养湖水中3种厌氧光合细菌,闪囊菌(Lamprocystispurpurea)、奥氏着色菌(Chromatiumokenii)和泥生绿菌(Chlorobiumclathratiforme)对无机碳(H13CO⁃3)和氮(15NH+4)的吸收特征,结果表明,仅占细胞总数0.3%的奥氏着色菌贡献了光合细菌碳吸收总量的70%和氮吸收总量的40%,说明丰富度较低的奥氏着色菌在该湖水的碳氮循环过程起着主要作用。Dekas等[48]采用15N2标记培养实验与NanoSIMS⁃FISH分析相结合,发现厌氧甲烷氧化团聚颗粒可以进行生物固氮,且生物固氮由厌氧甲烷氧化古菌完成,此外,部分固定的氮可以转运至与其共生的非固氮的硫还原细菌中,该研究不仅揭示了全球碳、氮和硫循环之间的联系,还显示出将同位素标记技术与NanoSIMS⁃FISH联合使用,使微生物细胞的生理组成元素及其相关代谢和转运过程的研究可视化的潜力。3 总结与展望
微生物在全球尺度上着生物地球化学循环过程,影响生态系统的功能。理解生态系统功能,从群落和单细胞水平揭示自然环境中微生物的种类和功能、微生物之间以及微生物与环境因子之间的相互作用是必不可少的。分子生物学方法的发展和进步为环境微生物的研究提供了方便,但是不同方法有其内在的局限性(表1),在方法的选择和应用上需根据研究具体情况进行选择。例如,高通量测序、宏组学和单细胞水平研究的方法均是近些年快速发展的方法,与DGGE和T⁃RFLP等传统方法相比,微生物的多样性或功能基因的分析的广度、深度和分辨率均有显著提高,显示出了其前所未有的优势和应用前景。然而,这些方法本身仍存在一些问题:(1)高通量测序得到的多为几百个碱基或者更短的序列,组成完整的基因有一定困难;即使获得整条序列,由于很多基因的功能还不清楚,很难获得由序列到功能的准确图谱;(2)单细胞水平研究方法中,Raman检测到光谱偏移的同位素标记量需达到10atom%[14];FISH样品的自发荧光易影响或掩盖探针的荧光信号;样品的制备过程中,固定和脱水等步骤也可能会影响微生物细胞的同位素组成[18]等。此外,基于PCR的方法,由于引物的选择以及PCR本身的问题易对微生物的理解带来偏差,如产生假序列(chimeras)(可通过生物信息学分析移除,如Bellerophon[49]、ChimeraSlayer[50]等软件),另外,得到的信息主要以易于被PCR扩增出的优势种群为主,而对环境中丰度较低的稀有微生物种群认识不足。未来的工作中,宏组学和单细胞水平研究的方法将有助于揭示自然环境中微生物整个群落的分类、功能信息以及丰度较低但发挥重要功能作用的微生物种群。由于目前用于阐明微生物群落组成和功能的方法均基于参考数据库,这些参考数据库主要是由克隆、Sanger测序或基因组等组成,即使新测序技术和单细胞水平技术有巨大的应用潜力,基于培养方法的研究依然要继续,如果没有这些基础数据,不仅很多序列的功能都难以辨认,用于识别环境微生物种类或功能的核酸探针的准确性也受到。随着微生物分子生物学方法(测序技术、单细胞水平研究技术等)的不断改进和数据分析工具的不断发展,对环境中微生物群落组成、微生物之间的相互作用的认识以及对生态系统的功能的理解也在不断增加,为检验和构思生态学理论提供新机遇[4],也为提出环境修复、自然生态系统和
http://www.ecologica.cn
量不同[45],说明这种方法除了同时检测微生物的种类和功能外,还具有识别微生物细胞代谢异质性的潜力。
24期 姬洪飞 等:分子生物学方法在环境微生物生态学中的应用研究进展 9
农业生态系统的可持续管理的新方法奠定基础。
参考文献(References):
[1] BarnardR,LeadleyPW,HungateBA.Globalchange,nitrification,anddenitrification:areview.GlobalBiogeochemicalCycles,2005,19[2] UhlikO,LeewisMC,StrejcekM,MusilovaL,MackovaM,LeighMB,MacekT.Stableisotopeprobinginthemetagenomicsera:abridge[3] MackelprangR,WaldropMP,DeAngelisKM,DavidMM,ChavarriaKL,BlazewiczSJ,RubinEM,JanssonJK.Metagenomicanalysisofa[4] ZimmermanN,IzardJ,KlattC,ZhouJZ,AronsonE.Theunseenworld:environmentalmicrobialsequencingandidentificationmethodsfor[5] EpsteinSS.Thephenomenonofmicrobialuncultivability.CurrentOpinioninMicrobiology,2013,16(5):636⁃642.
ecologists.FrontiersinEcologyandtheEnvironment,2014,12(4):224⁃231.
permafrostmicrobialcommunityrevealsarapidresponsetothaw.Nature,2011,480(7377):368⁃371.towardsimprovedbioremediation.BiotechnologyAdvances,2013,31(2):154⁃165.(1):GB1007.
[6] TringeSG,HugenholtzP.Arenaissanceforthepioneering16SrRNAgene.CurrentOpinioninMicrobiology,2008,11(5):442⁃446.
emergingsequencingtechnologies.FEMSMicrobiologyLetters,2009,296(1):97⁃101.MicrobiologicalMethods,2004,58(2):169⁃188.
[7] NilssonRH,RybergM,AbarenkovK,SjökvistE,KristianssonE.TheITSregionasatargetforcharacterizationoffungalcommunitiesusing[8] KirkJL,BeaudetteLA,HartM,MoutoglisP,KlironomosJN,LeeH,TrevorsJT.Methodsofstudyingsoilmicrobialdiversity.Journalof[9] 邢德峰,任南琪.应用DGGE研究微生物群落时的常见问题分析.微生物学报,2006,46(2):331⁃335.[10] 余素林,吴晓磊,钱易.环境微生物群落分析的T⁃RFLP技术及其优化措施.应用与环境生物学报,2006,12(6):861⁃868.
159⁃175.
[11] GentryTJ,WickhamGS,SchadtCW,HeZ,ZhouJ.Microarrayapplicationsinmicrobialecologyresearch.MicrobialEcology,2006,52(2):[12] LomanNJ,MisraRV,DallmanTJ,ConstantinidouC,GharbiaSE,WainJ,PallenMJ.Performancecomparisonofbenchtophigh⁃throughput[13] MinocheAE,DohmJC,HimmelbauerH.Evaluationofgenomichigh⁃throughputsequencingdatageneratedonIlluminaHiSeqandGenome[14] HuangWE,StoeckerK,GriffithsR,NewboldL,DaimsH,WhiteleyAS,WagnerM.Raman⁃FISH:combiningstable⁃isotopeRaman
1878⁃1889.
Analyzersystems.GenomeBiology,2011,12:R112⁃R112.sequencingplatforms.NatureBiotechnology,2012,30(5):434⁃439.spectroscopyandfluorescenceinsituhybridizationforthesinglecellanalysisofidentityandfunction.EnvironmentalMicrobiology,2007,9(8):
[15] Pett⁃RidgeJ,WeberPK.NanoSIP:NanoSIMSapplicationsformicrobialbiology.MethodsinMolecularBiology,2012,881:375⁃408.
deposition(CARD⁃FISH).SoilBiologyandBiochemistry,2008,40(7):1883⁃1891.NanoSIMS.EnvironmentalMicrobiology,2008,10(3):580⁃588.Spectrometry,2004,15(4):478⁃485.
[16] EickhorstT,TippkotterR.Improveddetectionofsoilmicroorganismsusingfluorescenceinsituhybridization(FISH)andcatalyzedreporter[17] LiTL,WuTD,MazéasL,ToffinL,Guerquin⁃KernJL,LeblonG,BouchezT.Simultaneousanalysisofmicrobialidentityandfunctionusing[18] PeteranderlR,LecheneC.Measureofcarbonandnitrogenstableisotoperatiosinculturedcells.JournaloftheAmericanSocietyforMass[19] AronsonEL,DubinskyEA,HellikerBR.Effectsofnitrogenadditiononsoilmicrobialdiversityandmethanecyclingcapacitydependondrainage[20] CongJ,LiuXD,LuH,XuH,LiYD,DengY,LiDQ,ZhangYG.Availablenitrogenisthekeyfactorinfluencingsoilmicrobialfunctional[21] Hermann⁃BankML,SkovgaardK,StockmarrA,LarsenN,MølbakL.Thegutmicrobiotassay:ahigh⁃throughputqPCRapproachcombinablewith[22] LuL,JiaZJ.Ureasegene⁃containingArchaeadominateautotrophicammoniaoxidationintwoacidsoils.EnvironmentalMicrobiology,2013,15[23] KimSJ,ParkSJ,ChaIT,MinD,KimJS,ChungWH,ChaeJC,JeonCO,RheeSK.Metabolicversatilityoftoluene⁃degrading,iron⁃
reducingbacteriaintidalflatsediment,characterizedbystableisotopeprobing⁃basedmetagenomicanalysis.EnvironmentalMicrobiology,2014,16(1):189⁃204.(6):1795⁃1809.
nextgenerationsequencingtostudygutmicrobialdiversity.BMCGenomics,2013,14:788⁃801.genediversityintropicalrainforest.BMCMicrobiology,2015,15:167⁃167.conditionsinapineforestsoil.SoilBiologyandBiochemistry,2013,62:119⁃128.
[24] ClemmensenKE,BahrA,OvaskainenO,DahlbergA,EkbladA,WallanderH,StenlidJ,FinlayRD,WardleDA,LindahlBD.Rootsand
associatedfungidrivelong⁃termcarbonsequestrationinborealforest.Science,2013,339(6127):1615⁃1618.
http://www.ecologica.cn
10 生 态 学 报 36卷
[25] KuypersMMM,JørgensenBB.Thefutureofsingle⁃cellenvironmentalmicrobiology.EnvironmentalMicrobiology,2007,9(1):6⁃7.
[26] RadajewskiS,InesonP,ParekhNR,MurrellJC.Stable⁃isotopeprobingasatoolinmicrobialecology.Nature,2000,403(6770):646⁃649.[27] CadischG,EspanaM,CauseyR,RichterM,ShawE,MorganJAW,RahnC,BendingGD.Technicalconsiderationsfortheuseof
15
[28] SchwartzE.CharacterizationofgrowingmicroorganismsinsoilbystableisotopeprobingwithH218O.AppliedandEnvironmentalMicrobiology,[29] PutkinenA,LarmolaT,TuomivirtaT,SiljanenHMP,BodrossyL,TuittilaES,FritzeH.Peatlandsuccessioninducesashiftinthecommunity[30] DumontMG,PommerenkeB,CasperP.UsingstableisotopeprobingtoobtainatargetedmetatranscriptomeofaerobicmethanotrophsinlakecompositionofSphagnum⁃associatedactivemethanotrophs.FEMSMicrobiologyEcology,2014,88(3):596⁃611.2007,73(8):2541⁃2546.
stable⁃isotopeprobingforfunctionalmicrobialactivityinsoils.RapidCommunicationsinMassSpectrometry,2005,19(11):1424⁃1428.
N⁃DNA
[31] sediment.EnvironmentalMicrobiologyReports,2013,5(5):757⁃764.
[32] Lu[33] 刘莉扬YH,Gerlach,W,崔鸿飞ConradJünemann,R.田埂Insitu.S,高通量测序技术在宏基因组学中的应用stableisotopeprobingofmethanogenicarchaea.中国医药生物技术inthericerhizosphere.,2013,Science,8(3):196⁃2005,200.
309(5737):1088⁃unassembledmetagenomicreads.TilleBMCF,GoesmannBioinformatics,A,Stoye2009,J.WebCARMA:10:430⁃430.
awebapplicationforthefunctionalandtaxonomicclassification[34] MeyerThemetagenomicsF,PaarmannRASTD,D′Souzaserver⁃aM,publicOlsonresourceR,GlassforEtheM,automaticKubalM,phylogeneticPaczianT,RodriguezandfunctionalA,StevensanalysisR,ofWilkemetagenomes.A,WilkeningBMCJ,Bioinformatics,
Edwards[35] 2008,Goecks9:computationalJ,386⁃Nekrutenko386.
researchA,intheTaylorlifeJ,sciences.TheGalaxyGenomeTeam.Biology,Galaxy:2010,acomprehensive11(8):R86⁃R86.approachforsupportingaccessible,reproducible,and[36] [37] ReichFiererM,N,LiefeldLauberT,CL,GouldRamirezJ,LernerKS,J,ZaneveldTamayoJ,P,BradfordMesirovMJA,P.KnightGenePatternR.Comparative2.0.Naturemetagenomic,Genetics,2006,phylogenetic38(5):and500⁃physiological500.
[38] soilXiamicrobialrecharacterizedY,WangcommunitiesYB,FangacrossHHP,nitrogenJinT,gradients.ZhongHTheZ,ISMEZhangJournal,T.Thermophilic2012,6:microbial1007⁃1017.cellulosedecompositionandmethanogenesis[39] SorokinYI.RadioisotopicbymetatranscriptomicMethodsinandHydrobiology.metagenomicHeidelberg:analysis.ScienceJ,WhiteleySpringer,Reports,AS.Raman1999.2014,4:6708⁃6708.
[40] 2004,Huang76(15):WE,Griffiths4452⁃R4458.I,ThompsonIP,BaileyMmicroscopicanalysisofsinglemicrobialcells.Analytical[41] 2009,BoxerS38:G,53⁃Kraft74.
ML,WeberPK.Advancesinimagingsecondaryionmassspectrometryforbiologicalsamples.AnnualReviewof[42] OrphanbacteriaVwithinJ,TurkstructuredKA,GreensyntrophicAM,consortiaHouseCrevealedH.PatternsbyFISH⁃SIMS.of
15
Nassimilation[43] MusatN,FosterR,VagnerT,AdamB,KuypersMMM.DetectingEnvironmentalandgrowthmetabolicactivitiesMicrobiology,ofmethanotrophicinsingle2009,cells,11(7):ANME⁃2withemphasis1777⁃archaea1791.
andsulfate⁃reducing
MicrobiologyReviews,2012,36(2):486⁃511.onnanoSIMS.[44] Abraham2014,98(11):WR.Applications4817⁃4828.
andimpactsofstableisotopeprobingforanalysisofmicrobialinteractions.AppliedMicrobiologyandBiotechnology,[45] HuangwithinWE,FergusonA,SingerAC,LawsonK,ThompsonIP,KalinRM,LarkinMJ,BaileyMJ,Whiteley[46] hybridization.microbial[47] 胡行伟Musat,张丽梅Appliedpopulations:,贺纪正and.Environmentalinsituanalyses纳米二次离子质谱技术Microbiology,usingrRNA(NanoSIMS)2009,and75(1):mRNAstable在微生物生态学研究中的应用234⁃241.isotopeprobingcoupledwithsingle⁃cellAS.Resolvingraman⁃fluorescencegenetic.生态学报,2013,33(2):ontheN,ecophysiologyHalmH,WinterhollerB,HoppeP,bacteria.PeduzziS,ProceedingsHillionF,ofHorreardtheNationalF,AmannAcademyR,JorgensenofSciencesBofB,theKuypersUnitedMStatesMM.ofA348⁃America,single⁃cell[48] 105(46):DekasAE,17861⁃Poretsky17866.ofanaerobicphototrophic(5951):422⁃426.
RS,OrphanVJ.Deep⁃seaarchaeafixandsharenitrogeninmethane⁃consumingmicrobialconsortia.Science,[49] Huber20(14):T,2317⁃Faulkner2319.
G,HugenholtzP.Bellerophon:aprogramtodetectchimericsequencesinmultiplesequencealignments.Bioinformatics,[50] HaasDeSantisBJ,GeversD,EarlAM,FeldgardenM,WardDV,GiannoukosG,CiullaD,TabbaaD,HighlanderSK,SodergrenE,SangerandTZ,454⁃pyrosequencedTheHumanMicrobiomePCRamplicons.Consortium,GenomePetrosinoResearch,JF,Knight2011,R,21(3):Birren494⁃BW504..Chimeric16SrRNAsequenceformationandhttp://www.ecologica.cn
ofRA.
ofpathwaysFEMSfunctions
insitu2008,
view
3262004,B,
in1090.transparentanalysesChemistry,Biophysics,357.
2009,Methédetection
因篇幅问题不能全部显示,请点此查看更多更全内容