环境微生物转录组学研究进展

环境微生物转录组学研究进展

自20世纪90年代以来，生物技术的快速发展极大地促进了人们对自然微生物多样性的了解。特别是近年来随着高流量统计数据的开发，全面微生物完整序列的研究已成为一种传统的实验手段。但是基因组序列本身不能提供微生物的基因功能信息,因此开展微生物转录组学研究尤为重要。环境微生物转录组学以环境样品中微生物的全部转录本即mRNA为研究对象,从群体水平上研究环境微生物功能基因的表达水平及其在不同环境条件下的转录调控规律,是一门研究微生物与自然环境相互关系的新兴学科。然而,环境微生物转录组学的研究进展非常缓慢,主要因为:(1)mRNA占环境样品总RNA中提取量的比例较低,通常仅为1–5%(Neidhardt&Umbarger,1996);(2)mRNA非常不稳定,极易被RNases降解,半衰期在几秒到几分钟之间(Deutscher,2006);(3)原核生物mRNA通常无多聚A尾,无法像真核生物mRNA那样可通过oligo(dT)引物直接选择(Poretskyetal.,2005);(4)后续的分析中均涉及到大量的酶促反应如PCR扩增和测序,但因环境样品组成复杂,其中的重金属和腐殖质等物质都可能降低RNA纯度,严重干扰下游分析;(5)测序技术是环境转录组学发展的重要技术瓶颈,单个微生物细胞的基因通常高达上千个,复杂环境中微生物的基因数量更是难以估算,而传统的测序技术通量较低,每次仅获得几百个mRNA序列,信息量较小,很难获得令人信服的微生物生理代谢过程的原位证据;同时传统测序技术依赖于克隆文库构建,费时费力且成本高昂,无法应对转录组研究对大规模测序的需求。自2008年以来,DNA测序技术取得了重大突破,基于不同原理的新一代高通量测序技术逐渐成为目前的主流平台(Shendure&Ji,2008)。与传统的Sanger测序法相比,新一代测序技术的效率大幅提高,一次测序可以获得高达上百万的通量,单个碱基的测序成本也大幅下降。新一代测序技术的快速发展催生了RNA-Seq(High-throughputRNAsequencing)技术的出现。RNA-Seq是利用新一代高通量测序技术对转录组进行分析的技术(Wangetal.,2009b),相比其他基于杂交技术的芯片技术(microarrays)(Schenaetal.,1998)、基于序列分析的基因表达系列分析(serialanalysisofgeneexpression,SAGE)(Velculescuetal.,1995)及其改进版的大规模平行信号测序系统(massivelyparallelsignaturesequencing,MPSS)(Brenneretal.,2000),RNA-Seq具有明显的优势,如检测范围广、无需特异性探针、高灵敏度等(Qietal.,2011)。RNA-Seq使复杂环境样品中微生物转录组学研究变得简单易行,近年来在水体(Frias-Lopezetal.,2008;Gilbertetal.,2008;Hollibaughetal.,2011)、土壤(Leiningeretal.,2006;Urichetal.,2008)、污泥(Yu&Zhang,2012)、沉积物(Millsetal.,2012)、动物肠道(Poroykoetal.,2010;Turnbaughetal.,2010)等多种环境中得到广泛应用,在环境微生物学领域受到极大关注。1环境中的微生物转移技术流程环境微生物转录组学的技术流程可大致分为以下几个步骤:1.1采集和保存方法环境样品一般分为水样、土壤和寄主样品(如人体微生物)三大类,不同样品类型决定了具体的采集和保存办法。一般的原则是严格控制样品采集的时间,样品采集后要利用液氮原位快速冷冻,并置于-70℃长期保存,或者将样品浸入RNA保存液(如Ambion公司的RNAlater保存液)中,-20℃长期保存。1.2rna提取试剂盒不同环境来源的样品RNA的提取方法各不相同,但一般都包括细胞裂解、蛋白抽提、核酸沉淀、DNA酶解等步骤。另外,有多种商业化的RNA提取试剂盒可供选择,如RNeasyMiniKit(Qiagen)、mirVanamiRNAIsolationKit(Ambion)、PowerSoil™TotalRNAIsolationKit(MoBioLaboratories)等。1.3非试剂盒测定由于mRNA在总RNA中含量很低,所以通常情况下需要对mRNA进行富集。由于原核生物的mRNA不含多聚A尾,所以不能用oligo(dT)直接进行选择,只能通过消减rRNA的方法间接富集mRNA。这一步骤一般采用试剂盒进行,如Ambion公司的MICROBExpress™BacterialmRNAEnrichmentKit、EpiCentre公司的mRNA-ONLY™ProkaryoticmRNAIsolationKit等,还有一些非试剂盒的物理化学方法。通常情况下,mRNA的含量会提高数倍到数十倍,最高可达95%左右(Giannoukosetal.,2012)。此外,mRNA可以用化学水解法或酶降解法生成较小的片段以适应测序平台的需要,例如,Illumina或SOLiD测序平台适合200–250nt大小的片段(Chu&Corey,2012)。片段化步骤也可以在cDNA合成后进行。1.4c有效使用可实现cd的第二链以富集后的mRNA为模板,经反转录酶催化生成cDNA的第一链,再经DNA聚合酶等多种酶的共同作用生成cDNA的第二链。这一步骤的实现也有多种试剂盒可供选择,如Invitrogen的SuperScript®Double-StrandedcDNASynthesisKit,Promega公司的UniversalRiboClone®cDNASynthesisSystem等。引物一般选择随机六聚体引物。1.5测序步骤不同纯化后的cDNA经过一系列的前处理后,进行高通量测序。这些步骤随不同的测序平台而有差异。通常情况下,高通量测序的数据为FASTQ格式,数据分析的过程因样品类型和研究目的不同而具有明显的差异,在数据处理及分析部分会详细阐述。2目标微生物的rna含量和测序量环境微生物样品和纯培养菌株样品在许多方面存在差异,首先是环境样品微生物多样性虽然高,但目标微生物的RNA含量非常低,受到样品量的限制,大多数情况下不能满足测序量的要求;其次是环境样品中存在多种杂质如腐植酸的干扰,另外物种组成非常复杂,rRNA消减程度有限,这些都加大了环境样品研究的难度。但近年来,这些方面的研究都有不同程度的进展。2.1添加样品量、纯化程序及添加量高通量测序一般需要微克级的RNA,而RNA的产量受诸多因素的影响。以土壤样品为例,土壤中活性微生物的多少、RNA提取方法的选择以及样品的纯化步骤在很大程度上都可能影响RNA的产量。不同类型的土壤中活性微生物含量差别巨大,每克土壤RNA的提取量多则几微克,少时仅几十纳克(Wangetal.,2012),因此增加样品量是提高RNA产量的方法之一。RNA提取的方法有很多种,目前似乎没有一种方法可适合于所有的样品类型,即使同样的样品,用不同的方法提取,RNA产量也会有显著差异。因此,需要根据样品类型及研究目的选择适当的方法或改进相应的技术。环境微生物转录组学研究一般包含多次RNA纯化步骤,每次纯化过程都会造成一定量的RNA损失,因此选择质量好的纯化试剂盒能够最大程度地减少这一损失。除了在提取和纯化过程中对流程进行优化,在cDNA合成之前也可对mRNA进行扩增,Ambion公司的MessageAmp™BacteriaIIRNAAmplificationKit试剂盒可使RNA数量增加200–1,500倍,已在多个研究中得到应用(Frias-Lopezetal.,2008;Hewsonetal.,2010;deMenezesetal.,2012)。2.2rna的提取获取高质量的RNA是转录组学研究的基础。RNA纯化试剂盒并不能除掉所有的杂质。以土壤样品RNA提取为例,主要的挑战是腐植酸的污染,它会在核酸提取过程中与核酸共沉淀,并对后续的酶反应造成不同程度的抑制。近年来有多种方法被证明能够有效地去除腐植酸,并可应用到RNA提取过程中的不同阶段,如控制裂解的温度(Wangetal.,2009a)、低pH条件(Metteletal.,2010)、抽提液中加入CTAB(cetyltrimethylammoniumbromide)(Griffithsetal.,2000)或PVPP(polyvinylpolypyrrolidone)(Rajendhran&Gunasekaran,2008)、用PEG来沉淀RNA(Bürgmannetal.,2003)等。但有时仅改进提取条件仍然不够,对于腐植酸含量高的样品,需要对提取后的RNA进行进一步的纯化,最近的研究认为柱纯化方法,如SephacrylS-400凝胶过滤(Wangetal.,2009a)和Q-SepharoseFastFlow离子交换柱(Metteletal.,2010)是相对高效的方法。2.3mrna的富集和净化尽管高通量测序可以获得数量巨大的序列数据,但由于mRNA在总RNA中的比例很低,为了提高测序的有效性同时也为了避免浪费,需要去除总RNA中的rRNA以便对mRNA进行富集,这也是转录组学研究中的一个重要步骤。目前mRNA的富集方法主要包括:基于消减杂交法的Ambion公司的MICROBExpress™BacterialmRNAEnrichmentKit和EpiCentre公司的Ribo-Zero™rRNARemovalKit,以及基于核酸外切酶法的EpiCentre公司的mRNA-ONLY™ProkaryoticmRNAIsolationKit。MICROB-Express™BacterialmRNAEnrichmentKit较早进入市场,目前使用非常广泛,但是对于不同的环境样品的富集效率差异较大(Stewartetal.,2010)。mRNA-ONLY™ProkaryoticmRNAIsolationKit的富集效果一般,相关研究报道很少。Ribo-Zero™rRNARemovalKit是由EpiCentre最新开发的,Giannoukos等(2012)评价了该试剂盒的rRNA去除效果,对于3种GC含量不同的纯培养细菌的实验结果,rRNA含量均可降至1%以下。通过与其他试剂盒的比较,作者认为Ribo-Zero™rRNARemovalKit效率最高。此外,该试剂盒对人粪便样品来源的混合微生物中rRNA的去除效果也较好,rRNA也可降至5%以下,但对于其他类型的环境样品,该试剂盒的处理效果还不清楚。最近,利用“not-so-random”六聚体引物在cDNA合成过程中选择性富集mRNA的基本原理(Armouretal.,2009),NuGEN公司开发了一种Ovation®ProkaryoticRNA-SeqSystem试剂盒,该试剂盒可以单独使用或者和已有的rRNA消减方法结合使用。另外,还有一种基于双链特异性核酸酶(duplex-specificnuclease,DSN)的cDNA文库标准化方法(Zhulidovetal.,2004),DSN在cDNA双链变性和复性的过程中选择性地降解数量较多的cDNA双链,以此减少各种转录子的含量差异,有利于大量rRNA的去除和稀有转录子的检测(Yietal.,2011)。但由于该方法会改变不同转录子的相对含量,因此不适合转录子的定量研究。除了以上提到的试剂盒方法,还有一些非试剂盒方法如切胶纯化法(McGrathetal.,2008),通过将琼脂糖电泳图中rRNA之间的胶切下来回收达到去除rRNA富集mRNA的目的。但这种方法需要的RNA上样量较大,而且无法回收和rRNA迁移率一样的mRNA。另外,还可根据样品的物种组成自行设计样品特异性的反义rRNA探针用于消减rRNA(Stewartetal.,2010),但由于步骤繁琐,尚未得到广泛使用。3其他典型的测序平台2005年以来,陆续诞生了多种新一代高通量测序技术平台,如罗氏公司(Roche)的454高通量测序平台、Illumina公司的Solexa平台以及ABI公司的SOLiD平台。由于这些平台基于不同的测序原理,所以在平均读长、数据量、错误类型、运行时间以及测序费用上均有一定差异,可根据自己的研究需要进行选择(Chu&Corey,2012)。关于这些平台的测序原理及性能比较的文献已有很多(Shendure&Ji,2008),因此本文仅对这些具有代表性的平台作简要介绍。需要注意的是这些测序平台的名称在不同的文献中可能有所不同,如有些文献以其公司名指代,有些以具体的测序仪型号指代,还有的以其核心技术指代,如Riche454高通量测序平台也称焦磷酸测序(pyrosequencing)。3.1基于gsflx的新能源测序平台2007年罗氏公司在454公司原有的GenomeSequencer20(GS20)测序仪的基础上,推出第二代高通量测序仪GSFLX,系统性能有了显著提高,这就是Riche454高通量测序平台,是第一个商业化的新一代测序平台。其基础是焦磷酸测序技术。与其他测序平台相比,该技术平台最显著的优势是序列有效读长可达800bp,耗时较短,目前已发表的环境微生物转录组学相关研究文献多数采用这种测序平台。但缺点是测序通量较小、测序费用相对较高、准确率较低。3.2a测序平台的核心技术Illumina公司的新一代测序仪GenomeAnalyzer(GA,通常称为Solexa仪)最早由Solexa公司研发,后来包括Solexa公司和Illumina公司在内的4个公司合并后,研发了目前的Illumina/Solexa测序平台,其核心技术是“DNA簇”和“可逆性末端终结”(reversibleterminator)。该测序平台具有多项优越性能,如测序通量大、准确率高、成本低、所需样品量少,最近逐渐开始在环境微生物转录组学研究中应用(Giannoukosetal.,2012;Xiongetal.,2012;Yu&Zhang,2012),具有很好的前景。其缺点是读长较短、运行时间较长,但随着技术的改进,这些不足正逐步得到改善。目前已有多种型号的测序仪可供选择,如IlluminaHi-Seq和Mi-Seq高通量测序系统在序列读长方面都有一定的提高,分别达到了150×2bp和250×2bp。鉴于微生物基因的大小通常为400bp左右,因此250×2bp的测序读长可基本满足序列分析的要求。3.3solid和其它高通量测序平台SOLiD(supportedoligoligationdetetion)平台是ABI公司自主研发的高通量测序平台,与454及Solexa平台基于边合成边测序(sequencingbysynthesis,SBS)不同,SOLiD的测序原理为基于磁珠的大规模并行克隆连接DNA测序法,其优势是测序通量大、准确率高、成本低,但同样存在读长短、运行时间长的缺点。使用该平台的研究较少(Liuetal.,2011)。4数据分析工具与形式序列数据的分析是环境微生物转录组学研究的重点和难点。转录组数据分析涉及到一系列的工具和资源,形式主要有数据库、门户网站、网络服务及一些单机程序(Cardenas&Tiedje,2008)。主要的分析步骤及常用的方法如下。4.1测序的去除、自适应调整RNA-Seq数据通常以FASTQ格式输出和存储,这个格式包括每个读段(reads)的ID号、序列和质量分数。在进行正式的序列分析之前,一般需要对序列数据作一些预处理,如除掉一些短的、质量差的和低复杂度的序列(low-complexityreads,序列一般为几个碱基的简单重复),以及纠正一些测序错误。每个被测的碱基都有一个质量分数(qualityscore),代表着测序的精确度,低的质量分数预示着可能的错误。对于454和Illumina数据来说,目前已有相应的算法来完成错误的检验、去除或更正(Rougemontetal.,2008;Quinceetal.,2009)。还有一些测序前人为加上的测序接头、标签等也需要去除。另外,也可对序列的末端进行修剪,因为末端的错误率通常较高(Balzeretal.,2010)。由于测序长度通常局限在500bp以下,所以需要对这些读段进行组装以获得全长的转录子序列。对于不同类型的数据,有3种序列组装策略可供选择:有参考基因组的组装(reference-based或abinitio装配)策略、重头组装策略(denovostrategy)及上述两种组装策略的结合。每种组装策略在原理、步骤、相应的序列定位及组装工具及使用优缺点上都存在差异,关于这些组装策略及相关软件的详细介绍可参考Martin和Wang(2011)。对于微生物组成复杂的环境样品来说,通常情况下并没有可供参考的基因组,所以一般选择重头组装策略,组装的效果主要受测序深度的影响。经过序列组装可以大大减少后续的比对工作量,初步组装后获得的较长片段称为Contig或序列重叠群,将已知的Contig用代表未知序列的N连接,可以进一步组装成更长的片段,称为Scaffold,再进一步进行补洞处理,可以得到含N量最少且两端不能继续延长的片段,称为Unigene。组装后的序列即可进行后续的序列比对及各种分析。4.2rrna序列的确定4.3基于序列的功能分析进行mRNA分析的工具较多,目前环境转录组学研究中最常用到的工具有MG-RAST(MetaGenomeRapidAnnotationusingSubsystemsTechnology)(Meyeretal.,2008)和MEGAN。MG-RAST是一种在线宏基因组和宏转录组分析工具,MEGAN是单机版软件,目前MEGAN4版本的功能已有较大拓展,可同时对多个样品进行基因功能比较分析,在环境转录组研究中应用非常广泛。掌握这两种工具即可完成环境转录组学分析的大部分工作。常见的几类mRNA分析有分类学分析(taxonomicanalysis)、基因功能分析(functionalanalysis)、KEGGpathway分析和不同样品间的基因表达差异分析。分类学分析的目的是确定mRNA来源于哪些物种。用blast软件中BLASTx程序将除去rRNA后剩下的RNA序列翻译为6种可能的氨基酸序列,与蛋白数据库进行比对,常用的蛋白数据库有NCBInr(noredundant)、COG(ClustersofOrthologousGenes)、UniProtKB/Swiss-Prot、PIR(TheProteinInformationResource)、KEGG(KyotoEncyclopediaofGenesandGenomes)等。Bitscores大于设定值(如30),则认为该序列属于putativemRNA转录子,结果输出后用MEGAN等软件进行后续的分析。MEGAN软件使用LCA(LeastCommonAncestor)算法,参照NCBI分类学方法将序列进行分类,分类结果以树状图的方式呈现,树的节点代表每个分类单元,匹配到该分类单元的序列数量也会标注出来,以方便比较这些分类单元在系统发育树上的分布情况。分类结果由序列和分类单元间的匹配特异性程度决定:特异性程度越低,分类越偏向于分类树的根部。功能分析指mRNA编码基因的注释(annotation)及其功能分类。基因注释的依据是序列之间的同源性,同源性高意味着功能上的高度相似性。通过BLASTx在各个蛋白数据库中进行搜索,查找序列相似性最高的蛋白,从而获知该mRNA编码基因的功能注释信息。基因的注释是后续基因功能分类的基础。利用GO(GeneOntology)、COG(ClustersofOrthologousGroups)及SEED等数据库可以将已进行功能注释的基因分到不同的基因家族中。GO是基因本体联合会(GeneOnotologyConsortium)所建立的数据库,分为分子功能(MolecularFunction)、细胞成分(CellularComponent)和生物过程(BiologicalProcess)三大部分,是应用最广泛的一种基因功能分类体系(Ashburneretal.,2000)。使用Blast2GO软件(Conesaetal.,2005)可对功能注释后的基因进行GO功能注释,然后再用WEGO软件(Yeetal.,2006)统计GO功能分类,即可得知某个样品中全部转录子的基因功能分布情况。COG数据库即直系同源序列聚类数据库,是根据系统进化关系对细菌、藻类和真核生物的数十个完整基因组的编码蛋白分类构建而成,用于预测蛋白质的功能。利用COG下的COGNITOR程序,可以把某个蛋白质与所有COGs中的蛋白质进行比对,并把它归入适当的COG簇。基因功能分类还有一种重要的方法叫做SEEDsubsystem分类(Overbeeketal.,2005),通过MG-RAST可以对mRNA编码基因进行SEEDsubsystem分类,SEED数据库将不同功能的基因划分到不同的subsystem中。MG-RAST是以RAST服务器(Azizetal.,2008)为基础的改进版,通过一个免费注册的账号就可以自动进行功能基因的系统发育及功能分析。其优点是能够快速注释大量的短片段DNA序列。如果研究人员不想将自己未发表的数据上传到网上进行分析,还可下载单机版的MEGAN软件。比较分析表明,在进行SEEDsubsystem分类时,对于同一样品,MG-RAST和MEGAN匹配到各个subsystem的序列数量是非常相似的(Mitraetal.,2011)。KEGGpathway分析虽然本质上也是一种基因功能分析,但是相对于上述的功能分类方法更加系统化。pathway指的就是基因所参与的各种生物学过程,KEGG据此定义了自己的直系同源群,称为KEGGOrthology(KO),每个KO都有一个收录号,对应一个特定的pathway。通过KEGG可查询包括碳水化合物、各种氨基酸、核酸及其他多种有机物的生物降解的可能代谢途径,对代谢过程中相关的酶也进行了注释,是一种生物代谢过程分析的强大工具。通过将待分析序列匹配到KO收录号已知的参考序列,MG-RAST和MEGAN都可以完成未知样品序列的KEGGpathway分析。但是在给出的pathway结果中,MG-RAST仅用一种颜色标示出现在pathway中的酶,而MEGAN可以用不同的颜色梯度来表征这些酶的相对含量,这非常有助于对酶动力学调控的理解(Mitraetal.,2011)。KEGGpathway分析有助于进一步了解mRNA编码的基因产物参与的代谢通路及其行使的生物学功能(Altermann&Klaenhammer,2005;Kanehisaetal.,2008)。在环境转录组学研究中常常需要比较不同环境条件下样品中的基因表达差异,MG-RAST和MEGAN均可完成这类比较分析。此外,RPKM(readsperkilobasespermillionreads)指数(Mortazavietal.,2008)在评估基因的表达水平时使用非常广泛,计算公式为RPKM=109×(C/NL)。其中C代表匹配到某个基因的序列数,N代表匹配到所有基因的序列数,L代表该基因的核酸序列长度。RPKM方法消除了序列长度及测序深度对基因表达量评估的可能影响,可直接用于比较不同样品间的基因表达水平差异。已有公开发表的软件用于计算RPKM值,如Cufflinks(Trapnelletal.,2010)和DEGseq(Wangetal.,2010)。4.4提交序列提交5环境中微生物转移的应用和前景5.1淡水生态环境基于RNA-Seq技术的环境微生物转录组学研究在近几年发展非常迅速,目前已应用于水体、土壤、淤泥、沉积物、动物肠道等多种环境类型。水体的研究目前大多集中于海水样品,以美国Georgia大学Moran教授课题组及美国麻省理工学院Parsons实验室的DeLong教授课题组为代表,近年来对多个地区海水微生物样品开展了一系列的转录组学研究,研究区域覆盖太平洋、大西洋等多个海区,类型涉及到表层海水(Frias-Lopezetal.,2008;Hewsonetal.,2010)、河口(Hollibaughetal.,2011)、人工模拟海洋生态系统(Gilbertetal.,2008)、海洋生物(水华藻类、海绵等)的附生微生物(Hewsonetal.,2009;Radaxetal.,2012)、海洋不同水层等(Stewartetal.,2011b)。对于淡水生态环境及内陆水体的研究目前还很少见(Liuetal.,2011)。关于土壤、淤泥及沉积物的研究也有一些,涉及到的样品类型有农业土壤(Leiningeretal.,2006)、贫瘠的沙质草地(Urichetal.,2008)、温带森林表层土壤(Stewartetal.,2011a)、人工模拟土壤环境(deMenezesetal.,2012)、污水处理厂活性淤泥(Yu&Zhang,2012)、海底沉积物(Millsetal.,2012)等。此外,关于人和动物的肠道(Poroykoetal.,2010;Turnbaughetal.,2010;Gosalbesetal.,2011;Xiongetal.,2012)以及排泄物(Booijinketal.,2010)的微生物转录组学研究也有不少。5.2其他控制变量的研究RNA-Seq技术最显著的特点是高通量,借助该技术,每个样品可得到数十万甚至上百万条序列数据,通过对这些海量数据的分析可以获得传统方法无法得到的信息。第一,高通量的数据促进了大量的功能未知的新基因及小RNA的发现。例如,在对太平洋表层海水样品进行微生物转录组学分析后,发现高达50%的基因为新基因(Frias-Lopezetal.,2008);Gilbert等(2008)在挪威海岸一个人工模拟海洋生态系统(mesocosm)中进行的微生物转录组学分析显示,RNA-Seq技术在环境微生物新基因的发现上具有强大的能力,在一些高度表达的大基因家族中,约91%的基因是新成员。另外,借助RNA-Seq技术,Shi等(2009)发现在海水中存在独特的微生物小RNA群,它们和已知的小RNA不同,位于基因组上的基因间隔区,推测可能具有基因调节功能。第二,通过比较微生物在不同环境(自然环境或人工控制环境)条件下的高通量转录组数据,结合宏基因组及定量PCR等其他研究手段,有助于发现环境条件对微生物代谢活性的影响以及微生物应对环境变化而进行的转录调控,了解微生物的转录活性对时空变化的响应模式。这一类应用是多方面的,研究对象可以是样品中的整个微生物群落。如Poretsky等(2009)通过比较北太平洋副热带环流系的微生物转录组数据,发现该区域的微生物在白天的代谢活动主要集中在光合作用、氧化磷酸化及C1化合物的合成,到了夜间以细胞膜、氨基酸及维生素的合成代谢为主。Hewson等(2010)比较了全球开阔大洋中8个不同采样点的微生物转录组数据,发现不同位点的微生物基因表达差异主要集中在几个关键的基因表达路径上,包括光系统I、II和氨吸收,原因可能是海洋表层水中的优势微生物为光能自养的原绿球藻,不同采样点微生物转录组的差异主要来源于这种蓝细菌光营养代谢的差异。通过控制实验或者与定量PCR等其他研究手段结合,也可以对某类功能菌群进行有针对性的研究。例如,通过分析加入菲的土壤样品中的微生物转录组数据,发现涉及芳香族化合物代谢及胁迫应答的转录子显著增加,可以获知多环芳烃(PAH)这类有毒污染物对土壤微生物活性的影响以及土壤微生物对PAH胁迫的应答模式;另外还第一次发现了重金属P型ATP酶和硫氧还蛋白与PAH胁迫相关联(deMenezese