吳波, 馮凱, 職曉陽(yáng), 何強(qiáng),5, 許玫英, 鄧曄, 肖凡書, 汪善全, 于玲,魯祺鴻,連英麗,羅麗娟,原珂,陳保衛(wèi),顏慶云, 仇榮亮, 欒天罡, 賀志理,9
(1.中山大學(xué)環(huán)境科學(xué)與工程學(xué)院環(huán)境微生物組研究中心,廣東 廣州 510006;2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心,北京100085;3.中國(guó)科學(xué)院大學(xué),北京100049;4.云南大學(xué)生命科學(xué)學(xué)院微生物研究所, 云南 昆明 650091;5.田納西大學(xué)城市與環(huán)境工程系及環(huán)境生物技術(shù)研究中心,田納西 諾克斯維爾 37996;6.廣東省微生物研究所,廣東 廣州 510070;7.中山大學(xué)生命科學(xué)學(xué)院,廣東 廣州 510275;8.中山大學(xué)海洋科學(xué)學(xué)院,廣東 廣州 510275;9.俄克拉荷馬大學(xué)環(huán)境基因組學(xué)研究所, 俄克拉荷馬 諾曼 73019)
環(huán)境微生物組多樣性及功能研究進(jìn)展
吳波1, 馮凱2,3, 職曉陽(yáng)4, 何強(qiáng)1,5, 許玫英6, 鄧曄2, 肖凡書1, 汪善全1, 于玲1,魯祺鴻1,連英麗1,羅麗娟7,原珂8,陳保衛(wèi)8,顏慶云1, 仇榮亮1, 欒天罡7, 賀志理1,9
(1.中山大學(xué)環(huán)境科學(xué)與工程學(xué)院環(huán)境微生物組研究中心,廣東 廣州 510006;2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心,北京100085;3.中國(guó)科學(xué)院大學(xué),北京100049;4.云南大學(xué)生命科學(xué)學(xué)院微生物研究所, 云南 昆明 650091;5.田納西大學(xué)城市與環(huán)境工程系及環(huán)境生物技術(shù)研究中心,田納西 諾克斯維爾 37996;6.廣東省微生物研究所,廣東 廣州 510070;7.中山大學(xué)生命科學(xué)學(xué)院,廣東 廣州 510275;8.中山大學(xué)海洋科學(xué)學(xué)院,廣東 廣州 510275;9.俄克拉荷馬大學(xué)環(huán)境基因組學(xué)研究所, 俄克拉荷馬 諾曼 73019)
快速發(fā)展的測(cè)序技術(shù)和各種高通量手段使我們能夠深入探索環(huán)境微生物的多樣性及其功能。本文中提及的環(huán)境微生物組(environmental microbiome)界定為特定環(huán)境中所有微生物基因組的集合。環(huán)境微生物組研究的主要目的是解讀特定環(huán)境中微生物的多樣性、結(jié)構(gòu)、功能、相互作用和進(jìn)化,以及它們與環(huán)境因子和生態(tài)系統(tǒng)功能之間的關(guān)系等。本文重點(diǎn)關(guān)注土壤、水體及空氣中的微生物組研究,回顧了目前的高通量測(cè)序技術(shù)和主要的環(huán)境微生物組學(xué)數(shù)據(jù)分析工具,然后綜述了不同的生態(tài)環(huán)境中微生物組研究的主要的進(jìn)展,最后,提出了環(huán)境微生物組研究中的關(guān)鍵科學(xué)問(wèn)題,并突出了在未來(lái)環(huán)境微生物組研究中微生物功能和生態(tài)機(jī)制的重要性。
環(huán)境微生物組;微生物功能;宏基因組;生物信息學(xué)工具
微生物廣泛分布于我們生活的環(huán)境甚至是極端環(huán)境中,且在全球生產(chǎn)力、人類健康、生物地球化學(xué)循環(huán)(如碳、氮、磷、硫、金屬等)、食品和生物資源、以及環(huán)境污染物修復(fù)等方面發(fā)揮著重要的作用。研究表明,土壤微生物能幫助植物抵抗疾病、促進(jìn)生長(zhǎng)、增強(qiáng)抗逆性和影響作物的產(chǎn)量及質(zhì)量等[1-5]。
微生物通常分為3大類:細(xì)菌,古菌和真核生物,其個(gè)體大小差異非常大,從0.2 μm到200 μm不等,是公認(rèn)的多樣性最高的生物類群。但是因?yàn)闄z測(cè)限的制約,使得稀有微生物多樣性的檢測(cè)受到了極大的挑戰(zhàn)[6]。最近,基于35 000個(gè)樣地的高通量分子數(shù)據(jù),通過(guò)相似律結(jié)合對(duì)數(shù)正態(tài)模型預(yù)測(cè)全球微生物大概有1012種,但目前只有0.001%的種類被鑒定[7]。在此,“微生物組”指的是特定環(huán)境中的微生物基因組的集合,包括主體相關(guān)的生態(tài)系統(tǒng),而“微生物群”是指微生物的集合。環(huán)境微生物組學(xué)是指分析各種環(huán)境中的微生物基因組信息,包括自然環(huán)境(如土壤、沉積物、水和空氣)和人工環(huán)境(如生物反應(yīng)器、建筑)。環(huán)境微生物組研究的目的是理解微生物群落的多樣性、結(jié)構(gòu)、功能、相互作用、動(dòng)態(tài)演化等及其與環(huán)境因素和生態(tài)系統(tǒng)功能的關(guān)系。
相對(duì)于人類、動(dòng)物和植物微生物組,環(huán)境微生物組研究面臨著更多的挑戰(zhàn)。 第一,環(huán)境微生物的多樣性和豐度特別高,它們交織在一起形成極為復(fù)雜的關(guān)系網(wǎng)。1 g土壤可能含有數(shù)以百萬(wàn)計(jì)的微生物種類[8-9];全球海底沉積微生物豐度約為2.9×1029個(gè)細(xì)胞[10]。因此,表征多樣性極高的微生物群體仍然是一個(gè)非常艱巨的任務(wù)。其次,大部分(>99%)的微生物還無(wú)法培養(yǎng)[7, 11],所以研究它們?cè)诃h(huán)境中的功能更加困難。第三,由于環(huán)境變化,微生物會(huì)快速的適應(yīng)和進(jìn)化,產(chǎn)生了時(shí)間變化和生物地理學(xué)的動(dòng)態(tài)模式。此外,只有少部分的環(huán)境微生物被測(cè)序作為參考基因組,這限制了對(duì)微生物在環(huán)境中的功能作用和潛在機(jī)制的了解。
本文重點(diǎn)關(guān)注土壤、水體和空氣環(huán)境中的微生物組,并回顧了目前環(huán)境微生物組學(xué)技術(shù)的發(fā)展。我們也提出了要解決的關(guān)鍵科學(xué)問(wèn)題,突出了在未來(lái)環(huán)境微生物組研究中對(duì)微生物功能、相互作用和潛在機(jī)制的理解。
高通量測(cè)序技術(shù)的發(fā)展對(duì)我們探索微生物多樣性、組成、功能、作用和演變等基本問(wèn)題至關(guān)重要。同時(shí),生物信息學(xué)工具的開(kāi)發(fā)在解析和理解這些由高通量測(cè)序產(chǎn)生的海量數(shù)據(jù)中具有不可或缺的作用。因此,我們首先回顧了當(dāng)前的測(cè)序技術(shù)和與環(huán)境微生物組學(xué)分析相關(guān)的生物信息學(xué)工具。
測(cè)序技術(shù)已發(fā)展到第3代:包括第1代傳統(tǒng)的Sanger測(cè)序、第2代高通量測(cè)序、第3代單分子測(cè)序。目前,第2代測(cè)序技術(shù),尤其是Illumina 平臺(tái)(例如,Miseq, Hiseq, Hiseq X, Moleculo, NovaSeq等)對(duì)于環(huán)境微生物組學(xué)研究有許多優(yōu)點(diǎn),具有通量高、錯(cuò)誤率和成本低等特征[12-13]。Illumina 測(cè)序平臺(tái)具有能獲得不同測(cè)序讀長(zhǎng)和通量的測(cè)序儀,例如,Moleculo能夠產(chǎn)生更長(zhǎng)的讀長(zhǎng)(>8 kb)[14]。
單分子測(cè)序技術(shù)的開(kāi)發(fā)為微生物的研究提供了更便利的條件。由于不需進(jìn)行DNA擴(kuò)增,因此極大地提高了測(cè)序精度,降低了錯(cuò)誤率[12]。目前,最廣泛使用的是太平洋生物科學(xué)公司(PacBio)基于單分子實(shí)時(shí)(SMRT)測(cè)序技術(shù)的平臺(tái)[15],這種測(cè)序方法取決于DNA聚合酶和零模波導(dǎo)孔所指示的光能。PacBio測(cè)序儀能在很短的時(shí)間內(nèi)產(chǎn)生長(zhǎng)達(dá)10~25 kb的讀長(zhǎng),并且能檢測(cè)DNA的甲基化修飾[16]。另外一種有前景的單分子測(cè)序由牛津納米科技公司( Oxford Nanopore Technologies )開(kāi)發(fā)的納米孔測(cè)序,其中單個(gè)DNA分子在通過(guò)納米孔時(shí)被直接測(cè)序[17]。雖然PacBio和納米孔測(cè)序技術(shù)都可以檢測(cè)DNA修飾堿基[17],但是都存在測(cè)序通量低,測(cè)序錯(cuò)誤率高的問(wèn)題[13],從而限制了它們?cè)诃h(huán)境微生物組研究中的應(yīng)用。而當(dāng)單分子測(cè)序技術(shù)和第二代測(cè)序技術(shù)結(jié)合使用時(shí),Tsai等[18]能夠檢測(cè)微生物的單核苷酸多態(tài)性水平的變化,表明這種結(jié)合使用的方法將成為環(huán)境微生物研究中進(jìn)行宏基因組定量、重建以及功能注釋的強(qiáng)有力手段。
通常環(huán)境微生物研究可以產(chǎn)生兩種類型的數(shù)據(jù):擴(kuò)增子及鳥槍法宏基因組/宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)。由于它們的分析方法不同,因此相應(yīng)地開(kāi)發(fā)了多種數(shù)據(jù)分析工具。
1.2.1 擴(kuò)增子測(cè)序數(shù)據(jù) 系統(tǒng)發(fā)育和功能性標(biāo)記基因?qū)τ谔剿鳝h(huán)境微生物的分類學(xué)、系統(tǒng)發(fā)育以及種群多樣性及其功能具有重要意義。系統(tǒng)進(jìn)化分子標(biāo)記基因中,用于標(biāo)記細(xì)菌和古細(xì)菌的16S rRNA基因[19]及用于真菌的核糖體DNA內(nèi)轉(zhuǎn)錄間隔區(qū)(ITS)[20]最為常見(jiàn),因此針對(duì)系統(tǒng)發(fā)育標(biāo)記基因擴(kuò)增產(chǎn)物的分析開(kāi)發(fā)出了許多生物信息學(xué)工具(例如,QIIME、Mothur)。
功能基因擴(kuò)增子測(cè)序可以提供與微生物群落功能直接相關(guān)的信息,例如nifH基因與固氮菌[21],及amoA基因與硝化細(xì)菌[22]功能的相關(guān)性,而這些在系統(tǒng)發(fā)育基因中的研究較少。功能基因擴(kuò)增子分析與基于系統(tǒng)發(fā)育標(biāo)記的分析類似,但功能基因擴(kuò)增子分析需要額外的步驟進(jìn)行移碼校正來(lái)預(yù)測(cè)蛋白質(zhì)編碼區(qū),并與相應(yīng)的功能基因數(shù)據(jù)庫(kù)比對(duì)以獲得相應(yīng)的功能信息。但如果沒(méi)有對(duì)應(yīng)的功能基因數(shù)據(jù)庫(kù),就需要為每個(gè)功能基因家族構(gòu)建一個(gè)數(shù)據(jù)庫(kù)。值得關(guān)注的是,F(xiàn)unGene提供多種功能基因數(shù)據(jù)庫(kù),涉及多種生態(tài)途徑(如C、N循環(huán))以及集成功能基因分析工具,如mcClust, RDPTools及FrameBot[23]。值得一提的是,其中的一些具體算法和工具,如PICRUSt,Tax4Fun,F(xiàn)APROTAX,F(xiàn)UNGuild及PAPRICA可通過(guò)系統(tǒng)發(fā)育標(biāo)記基因序列推斷微生物功能特征。
1.2.2 鳥槍法宏基因組及宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù) 鳥槍法宏基因組及宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)確保了信息從遺傳物質(zhì)到功能角度的傳遞,并為探索所有遺傳要素及其在環(huán)境中的動(dòng)態(tài)提供了可能[24]。通常宏基因組測(cè)序數(shù)據(jù)分析包括基于讀長(zhǎng)的分析和基于組裝的分析?;谧x長(zhǎng)的分析包括以下5個(gè)步驟:① 對(duì)原始讀長(zhǎng)進(jìn)行預(yù)處理;② 將序列比對(duì)到分類群;③ 基因預(yù)測(cè);④ 功能基因注釋;⑤ 豐度計(jì)算及統(tǒng)計(jì)分析?;诮M裝的分析與基于讀長(zhǎng)的分析類似,除了在基因預(yù)測(cè)之前需要使用多種組裝方法將短的原始讀長(zhǎng)組裝成更長(zhǎng)的片段,組裝方法包括基于參考序列的組裝和從頭組裝。
由于對(duì)計(jì)算資源的高度要求,基于宏基因組測(cè)序數(shù)據(jù)的組裝分析一般是以網(wǎng)絡(luò)為基礎(chǔ)的工作流程[25]。目前,常用的宏基因組分析工具/流程是MEGAN[26],MG-RAST[27],IMG/M[28]和CAMERA[29]。最近開(kāi)發(fā)的MEGAN CE可用于大規(guī)模的微生物測(cè)序數(shù)據(jù)分析,是MEGAN的更新版本[30]。MEGAN CE結(jié)合了DIAMOND[31]與MganServer,具有大規(guī)模鳥槍法宏基因組測(cè)序數(shù)據(jù)分析的能力[30]。MG-RAST[27]使用SEED子系統(tǒng)[32]對(duì)宏基因組測(cè)序數(shù)據(jù)分析,包含了公共宏基因組數(shù)據(jù)的注釋。IMG/M[28]使用多個(gè)能在微生物宏基因組深入分析時(shí)使用的參考數(shù)據(jù)庫(kù)(如IGM系統(tǒng)的基因組)以進(jìn)行宏基因組測(cè)序數(shù)據(jù)的注釋。CAMERA[29]包含大規(guī)模數(shù)據(jù)庫(kù)及多種注釋工具,最近已被納入QIIME[33]。
由于宏基因組測(cè)序數(shù)據(jù)分析需要大量的計(jì)算機(jī)資源并且非常耗時(shí),多種通過(guò)宏基因組測(cè)序數(shù)據(jù)進(jìn)行微生物種群檢測(cè)或鑒定的算法或方法已被開(kāi)發(fā)。GSMer是一個(gè)新的基于k-mer的方法,其利用現(xiàn)有微生物基因組特異標(biāo)記(genome-specific markers,GSMS)從菌株/物種水平對(duì)原始宏基因組測(cè)序數(shù)據(jù)進(jìn)行標(biāo)識(shí)[34]。潛在菌株分析(latent strain analysis,LSA)可以檢測(cè)宏基因組測(cè)序數(shù)據(jù)集中的菌株并識(shí)別低豐度的菌株(相對(duì)豐度小于0.000 01%)[35]。ConStrains采用單核苷酸多態(tài)性模式,從宏基因組測(cè)序數(shù)據(jù)中識(shí)別同種細(xì)菌菌株,并重建這些物種的親緣關(guān)系[36]。這些方法的發(fā)展提供了對(duì)宏基因組測(cè)序數(shù)據(jù)的快速、高分辨分析的替代手段,進(jìn)一步提高我們對(duì)微生物群落及其在環(huán)境中功能的認(rèn)識(shí)。
宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)分析也有與宏基因組測(cè)序數(shù)據(jù)分析類似的步驟和工具,但仍需要特殊的組裝及rRNA去除步驟及工具。宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)的組裝工具是基于de Bruijn圖形開(kāi)發(fā)用于從頭組裝,例如TAG[37],IDBA-MT[38],Trinity[39]及Oases[40]。由于參照宏基因組碎片化,TAG有更多的轉(zhuǎn)錄組丟失風(fēng)險(xiǎn)[37]。IDBA-MT[38]是一種從頭組裝方法,又進(jìn)一步發(fā)展為IDBA-MTP[41],通過(guò)利用已知的蛋白數(shù)據(jù)庫(kù)及相應(yīng)的核甘酸序列,來(lái)進(jìn)行宏轉(zhuǎn)錄組的組裝,從而有效的組裝mRNAs。從鳥槍法宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)去除rRNA的工具在有或沒(méi)有任何參考數(shù)據(jù)庫(kù)的情況下開(kāi)發(fā)。riboPicker[42]及rRNASelector[43]都可以從宏轉(zhuǎn)錄組序列中選擇并去除rRNA,而rRNAFilter[44]可用在沒(méi)有任何參考序列的情況下去除rRNA讀長(zhǎng)。此外,SortMeRNA可以識(shí)別mRNA和rRNA,從宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)集中快速過(guò)濾、去除rRNA基因片段[45]。對(duì)于其他類型的RNA,例如非編碼RNA以及小分子RNA(sRNA)[46],INFERNAL[47]能夠依靠Rfam數(shù)據(jù)庫(kù)搜索RNA結(jié)構(gòu)[48],對(duì)已知的sRNA進(jìn)行分配和進(jìn)一步功能注釋,這與通過(guò)對(duì)多種數(shù)據(jù)庫(kù)進(jìn)行比對(duì)的宏基因組測(cè)序數(shù)據(jù)分析類似。BioBakery是其中一種宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)分析綜合流程,SAMSA[49]是近年開(kāi)發(fā)的宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)綜合分析工具。
用于擴(kuò)增子、宏基因組及宏轉(zhuǎn)錄組測(cè)序數(shù)據(jù)分析的生物信息學(xué)工具的發(fā)展,使我們能夠表征環(huán)境微生物,并探究它們的多樣性及其在環(huán)境中的功能作用。這樣的數(shù)據(jù)集在大小和樣本數(shù)量上不斷增加,因此需要更加快速、高效、自動(dòng)、準(zhǔn)確及強(qiáng)大的生物信息工具與流程,這將帶領(lǐng)環(huán)境微生物的研究進(jìn)入一個(gè)新紀(jì)元。
測(cè)序技術(shù)、組學(xué)技術(shù)以及數(shù)據(jù)處理技術(shù)的發(fā)展,逐漸打破了環(huán)境微生物組研究中的技術(shù)瓶頸。在此,我們對(duì)具有代表性的生境中的微生物組研究進(jìn)行了歸納總結(jié),包括土壤(農(nóng)田、草地、森林、苔原)、水體(海洋和淡水)以及空氣微生物組。這些研究在環(huán)境微生物的多樣性、組成、功能、相互作用、進(jìn)化、演替以及環(huán)境生態(tài)功能等方面均有新的突破和見(jiàn)解。
土壤中的微生物多樣性極高,據(jù)估計(jì)每克土壤有109個(gè)微生物[9],隸屬104個(gè)不同的物種[50]。土壤微生物在地球化學(xué)循環(huán)過(guò)程中扮演著重要的角色[11],同時(shí)也是生態(tài)系統(tǒng)的維護(hù)者,如維持植物生產(chǎn)力[1]、實(shí)現(xiàn)污染物的生物修復(fù)[51]等。已有研究證明生物多樣性是影響生態(tài)系統(tǒng)功能的主要因素[52-54]。但是,以往大部分的研究主要通過(guò)分析16S rRNA 和 ITS基因來(lái)研究微生物的群落結(jié)構(gòu),而無(wú)法提供重要的功能信息。
宏基因組學(xué)方法(shotgun sequencing)則彌補(bǔ)上述缺陷,為我們提供微生物的功能信息。最近有研究利用宏基因組學(xué)的方法揭示了氮代謝途徑與生境類型、土壤碳、氮含量顯著相關(guān),同時(shí),該研究有效的將物種與功能聯(lián)系起來(lái),指出在土壤中,一些δ-變形菌(δ-Proteobacteria)具有完整的氮代謝途徑,而參與氮循環(huán)的藍(lán)細(xì)菌(Cyanobacteria)則只含有一些特定的氮代謝的途徑[55]。本文將重點(diǎn)介紹近年來(lái)自農(nóng)田、草原、森林、苔原等生境中土壤微生物組研究進(jìn)展,尤其是應(yīng)用了宏基因組和宏轉(zhuǎn)錄組學(xué)方法的研究。
2.1.1 農(nóng)田 作物相關(guān)微生物被認(rèn)為是作物的第二基因組,其對(duì)農(nóng)作物的健康和產(chǎn)量有著至關(guān)重要的作用。研究表明,使用PhyloChip研究抑病土壤的微生物群落,發(fā)現(xiàn)土壤微生物中γ-變形菌(γ-Proteobacteria)通過(guò)非核糖體肽合酶介導(dǎo)的抑制真菌的代謝途徑與農(nóng)作物病害抑制率密切相關(guān)[56]。另外,有關(guān)大豆根際微生物群落分類和功能的研究表明,生境選擇(niche-based selection)主導(dǎo)了微生物群落的組成[57]。同時(shí),通過(guò)宏轉(zhuǎn)錄組的方法研究大豆根系微生物對(duì)高濃度的CO2和O3的響應(yīng),發(fā)現(xiàn)碳、氮循環(huán)中的關(guān)鍵基因的表達(dá)被高含量的CO2刺激上調(diào),而被O3刺激下調(diào),從而揭示了土壤微生物應(yīng)對(duì)全球氣候變化的潛在機(jī)制[58]。這些研究都表明組學(xué)方法在揭示土壤微生物在維持植物健康和應(yīng)對(duì)環(huán)境變化的研究中是一個(gè)有利的工具。
2.1.2 草地 全球約有40%的土地面積被草地所覆蓋[59]。草地在生物地球化學(xué)循環(huán)中起著重要的作用,草地土壤微生物對(duì)環(huán)境及氣候變化的響應(yīng)成為研究熱點(diǎn)[60-63]。英國(guó)洛桑研究所通過(guò)宏基因組分析公園草地中的土壤微生物發(fā)現(xiàn),其功能變化受季節(jié)及垂直距離的影響異常地低。同時(shí)也發(fā)現(xiàn),碳代謝相關(guān)基因的豐度最高,cAMP信號(hào)因子、Ton和Tol轉(zhuǎn)運(yùn)系統(tǒng)等的表達(dá)量高[64]。
為了洞悉草地對(duì)全球氣候變化的相應(yīng),研究者在加溫2 ℃長(zhǎng)達(dá)10 年的草地土壤中,分析了微生物的宏基因組數(shù)據(jù),發(fā)現(xiàn)碳、氮循環(huán)的關(guān)鍵代謝途徑(如纖維所降解、反硝化等)在加溫條件下被顯著性富集[65]。另一組學(xué)研究調(diào)查了長(zhǎng)期高含量CO2對(duì)草地土壤微生物組成及其功能的影響,研究表明高含量的CO2提高了有機(jī)物降解途徑、異化硝酸鹽還原、氮固定中關(guān)鍵基因的豐度,而谷氨酰胺合成和厭氧氨氧化途徑的關(guān)鍵基因則降低了[66]??傊?,草地土壤微生物的功能特性在研究草地生態(tài)功能及草地應(yīng)對(duì)環(huán)境變化過(guò)程中有著十分重要的指示作用。
2.1.3 森林 森林生態(tài)系統(tǒng)是復(fù)雜且差異性大的陸地生境之一,為地球碳循環(huán)和生物多樣性的維持提供了重要的生態(tài)作用[67-68]。森林土壤微生物對(duì)森林生態(tài)系統(tǒng)的地球元素循環(huán)有著重要作用[69-70]。有研究利用宏基因組學(xué)分析方法發(fā)現(xiàn)人為擾動(dòng)對(duì)森林土壤微生物功能特性有著顯著性的影響,降低了木質(zhì)素、纖維素、半纖維素和果膠降解相關(guān)功能基因的豐度[71]。而利用宏轉(zhuǎn)錄組方法分析加氮肥對(duì)森林土壤微生物影響發(fā)現(xiàn),氮的添加會(huì)引起碳代謝途徑中酶的過(guò)表達(dá)但不會(huì)影響群落結(jié)構(gòu)和生物量的變化[72]。此外,通過(guò)系統(tǒng)發(fā)育和功能基因的研究發(fā)現(xiàn),重金屬污染會(huì)顯著影響森林土壤微生物的群落結(jié)構(gòu),但該研究也證實(shí)了森林土壤微生物在重金屬污染壓力下的恢復(fù)性,驗(yàn)證了金屬濃度與金屬抗性基因之間的關(guān)系,以此闡明了微生物功能分析的重要性[73]。利用宏基因組學(xué)研究森林土壤微生物中編碼鹵化和脫鹵酶這類特異性的功能基因表明,森林土壤微生物中涉及鹵素循環(huán)基因多樣性高,且發(fā)現(xiàn)了多種非特異性脫鹵以及協(xié)同脫鹵基因,表明了自然脫鹵的穩(wěn)健性[70]。以上研究表明,森林土壤微生物有多種多樣的已知和未知新型的功能,包括營(yíng)養(yǎng)循環(huán)、生物修復(fù)和其他生態(tài)系統(tǒng)功能。
2.1.4 苔原 苔原土壤主要分布在在北極圈和高山(高山林線以上),是一類沒(méi)有樹木生長(zhǎng)的生境[74]。其土壤層作為永久凍土帶,儲(chǔ)存了含量巨大的碳元素[75]。由于氣候變暖以及凍土解凍,逐漸升高的地表溫度影響著微生物的活性和相關(guān)的地球元素循環(huán)[76-78]。最近,越來(lái)越多的宏基因組學(xué)研究手段應(yīng)用于苔原微生物組學(xué)的研究。例如,對(duì)短時(shí)(1周)凍土解凍的苔原土壤微生物進(jìn)行的深度測(cè)序,研究表明碳氮循環(huán)相關(guān)基因,尤其是涉及甲烷代謝的相關(guān)基因發(fā)生了快速的響應(yīng)[79]。而且,多組學(xué)研究手段包含16S rRNA擴(kuò)增子測(cè)序,宏基因組以及宏轉(zhuǎn)錄組測(cè)序被用于分析不同凍土(完整的凍土,季節(jié)融化的活性層,和熱融沼澤)的微生物組,檢測(cè)其分類組成、功能活性的變化,結(jié)果展示了不同代謝(在沼澤中甲烷合成)速率和代謝相關(guān)的組學(xué)數(shù)據(jù)良好的相關(guān)性,而在凍融活性層則發(fā)現(xiàn)了大量的冷激蛋白以及趨化、運(yùn)動(dòng)相關(guān)的蛋白的基因[80]。這些有關(guān)苔原的研究闡明了凍土微生物組對(duì)全球氣候變化的快速響應(yīng),但也表明不同的功能微生物菌群在不同的時(shí)間和空間上響應(yīng)均有差異。因此,以上的研究有助于我們預(yù)測(cè)未來(lái)全球變暖對(duì)苔原生態(tài)功能的影響。
微生物廣泛分布于各類水環(huán)境,并在生物地球化學(xué)過(guò)程中扮演著重要的角色,包括水生生態(tài)系統(tǒng)的營(yíng)養(yǎng)循環(huán)和能量流動(dòng)[81]。水體生態(tài)系統(tǒng)中的浮游細(xì)菌是地球上最廣泛分布的一類微生物[82],如1mL表層海水中包含了104~106個(gè)細(xì)菌[83],全球海洋中的浮游細(xì)菌總數(shù)> 3.1×1028[82]。它們不僅在表層水中比較常見(jiàn),而且在深海海底也有廣泛的分布[84],被認(rèn)為是海洋C、N和S循環(huán)的引擎[85]。越來(lái)越多的研究發(fā)現(xiàn)這些生物地球化學(xué)過(guò)程是由環(huán)境中所有的成員(即群落)發(fā)生復(fù)雜交互作用引起的,而非由特定的某些物種引起[86]。因此,在群落水平上研究水體微生物群落是非常必要的,這樣我們才可以對(duì)生態(tài)系統(tǒng)功能有一個(gè)比較全面的認(rèn)識(shí)。
2.2.1 海洋 有研究提出海洋微生物的種類是無(wú)處不在的[87-88],但是近來(lái)越來(lái)越多的研究認(rèn)為浮游微生物群落的分布主要基于生態(tài)位過(guò)程的影響,如發(fā)現(xiàn)某些物種只分布在一些特定的區(qū)域[89-90]。此外,海洋中的細(xì)菌群落可反復(fù)出現(xiàn)特定的格局表明微生物群落和組成能夠通過(guò)環(huán)境參數(shù)來(lái)進(jìn)行預(yù)測(cè)[91]。有些環(huán)境因子,如鹽度[92]、溫度[93]能夠顯著影響浮游細(xì)菌的組成。值得一提的是,最近發(fā)現(xiàn)海洋生態(tài)系統(tǒng)中浮游細(xì)菌在長(zhǎng)時(shí)間尺度會(huì)隨環(huán)境的周期性改變而重復(fù)演替[94-95]。
雖然仍有不足,但是近期已經(jīng)使用宏基因組學(xué)、宏轉(zhuǎn)錄組學(xué)和單細(xì)胞測(cè)序等技術(shù)證明海洋微生物群落的功能特征。例如,基于單細(xì)胞基因組學(xué)對(duì)來(lái)自海洋“微生物暗物質(zhì)”的201株野生古細(xì)和細(xì)菌進(jìn)行的單細(xì)胞測(cè)序,揭露出野生古菌和細(xì)菌細(xì)胞令人意想不到的代謝特征[96-97]。利用宏基因組技術(shù)檢測(cè)紅海微生物群落的過(guò)量表達(dá)(over-representation)基因表明海洋微生物群落能獲得特定的環(huán)境適應(yīng)性[98]。對(duì)采自Tara海243個(gè)樣品的宏基因組測(cè)序數(shù)據(jù)分析顯示, 73%以上的海洋微生物與人類消化道菌群具有同樣的功能[99]。一項(xiàng)最近針對(duì)紅海微生物的宏基因組研究表明,只有功能協(xié)變模式不會(huì)隨著環(huán)境梯度的改變而改變,但是大多數(shù)的變化都是由溫度的改變導(dǎo)致的[100]。
2.2.2 淡水 淡水通常被定義為溶解鹽小于1 000 mg/L的天然水體[101],包括湖水、河水和地下水等。單細(xì)胞基因組學(xué)和宏基因組學(xué)測(cè)序分析顯示,淡水湖泊中好氧的化能自養(yǎng)菌比以前的預(yù)測(cè)更加普遍[102]。這些研究也為種群間的基因遷移和視紫紅質(zhì)基因重組提供了證據(jù)。宏基因組測(cè)序發(fā)現(xiàn)已知同源物微囊藻毒素-降解基因并沒(méi)有過(guò)多表達(dá),但是生物異源物質(zhì)代謝的基因在微囊藻毒素-修正微生態(tài)中過(guò)量表達(dá)[103],表明異生基因在微囊藻毒素降解中非常重要。
宏基因組測(cè)序和基于功能推理的方法發(fā)現(xiàn)“核心功能特性(core functional traits)”在密西西比河中是保守的[104]。利用宏基因組和轉(zhuǎn)錄組對(duì)哥倫比亞河不同鹽度的水樣品進(jìn)行測(cè)序發(fā)現(xiàn),宏基因組中的功能基因表達(dá)譜是非常類似的,但是與轉(zhuǎn)錄組相關(guān)的基因表達(dá)譜隨著河水鹽度梯度的變化而高度變異[105]。在時(shí)間尺度上,通過(guò)對(duì)河水微生物群落的研究,發(fā)現(xiàn)在水質(zhì)不好的環(huán)境中營(yíng)養(yǎng)代謝菌和噬菌體大量存在,主要是因?yàn)檗r(nóng)田徑流增加[106]。在空間上,涉及N循環(huán)的大量功能基因被檢測(cè)到沿著河流有明顯的變化[107]。
地下水是另一種重要的淡水資源,受到了人類活動(dòng)的威脅,包括工業(yè)副產(chǎn)物、農(nóng)業(yè)徑流和人類/動(dòng)物糞便的污染[108]。近年來(lái)淡水微生物群落受環(huán)境擾動(dòng)的響應(yīng)逐漸成為研究的熱點(diǎn)[108-110]宏基因組的比較分析表明受到污染的淡水微生物群落比未受污染的淡水微生物群落更脆弱,功能也衰退,從而營(yíng)養(yǎng)循環(huán)和適應(yīng)性減弱[108]。Hemme等[111]還發(fā)現(xiàn)淡水微生物群落在進(jìn)化過(guò)程中受到重金屬污染時(shí)會(huì)發(fā)生基因水平轉(zhuǎn)移。
大氣環(huán)境中的微生物在物種多樣性及豐度(如25.6~124.4 cfu /m3)[112]上都遠(yuǎn)低于土壤及水環(huán)境。但大氣微生物組不僅對(duì)大氣過(guò)程有顯著影響,而且對(duì)大氣物理、化學(xué)及氣候變化等都具有重要作用[113]。此外,大氣微生物組對(duì)人類健康、農(nóng)業(yè)生產(chǎn)力和生態(tài)系統(tǒng)穩(wěn)定性也極其重要[114]。然而,目前對(duì)于大氣微生物組的研究還處于起步階段,主要集中在對(duì)于室內(nèi)空氣微生物的研究方面[115]。
截至目前,關(guān)于大氣微生物的研究主要是針對(duì)住所、辦公區(qū)、教室、醫(yī)院、公共休息室、博物館和地鐵等環(huán)境的調(diào)查。通常,建筑內(nèi)部的空氣微生物主要來(lái)自人類、寵物、植物、城市管道系統(tǒng)、供暖/通風(fēng)/空調(diào)系統(tǒng)、霉菌和灰塵重懸浮以及室外空氣[116]。盡管微生物組成隨著地理位置、季節(jié)、用途和樣品基質(zhì)(是否接近表面)發(fā)生改變,但是這些建筑環(huán)境中豐度最高的通常是變形菌門、厚壁菌門和放線菌門[115]。將建筑周圍空氣中的微生物組與居住者攜帶的微生物組進(jìn)行對(duì)比分析,發(fā)現(xiàn)建筑周圍的空氣是建筑環(huán)境中微生物組的主要來(lái)源,而居住者皮膚攜帶的微生物組也是重要來(lái)源但其貢獻(xiàn)相對(duì)較小。
然而空氣流通、潮濕以及其他因素也可能對(duì)大氣微生物組產(chǎn)生重要影響,但很少有研究系統(tǒng)探討這些機(jī)制對(duì)微生物多樣性、組成及功能的影響。因此,需要開(kāi)展相關(guān)工作,以評(píng)價(jià)這些機(jī)制在大氣環(huán)境中微生物組形成中的重要性。而大氣環(huán)境微生物組如何影響人類健康則是另一個(gè)值得關(guān)注的科學(xué)問(wèn)題。盡管有許多研究對(duì)空氣環(huán)境中微生物的16S rRNA基因進(jìn)行了測(cè)序分析,但關(guān)于大氣環(huán)境中微生物功能和活性的信息卻很少,導(dǎo)致我們很難評(píng)估建筑周圍空氣中微生物對(duì)人類健康的潛在影響。近期,一種針對(duì)宏基因組測(cè)序的大氣微生物DNA提取方法被開(kāi)發(fā)出來(lái),在一定程度上緩解了大氣環(huán)境微生物的生物量過(guò)低導(dǎo)致研究難以開(kāi)展的問(wèn)題[117]。
多組學(xué)技術(shù)的發(fā)展,如DNA和RNA測(cè)序技術(shù)、蛋白和代謝譜分析、功能基因矩陣等提高了我們分析環(huán)境微生物的能力。隨著計(jì)算工具、統(tǒng)計(jì)和建模手段的發(fā)展,研究人員在過(guò)去10年取得了前文所述各項(xiàng)重要進(jìn)展。然而,大部分研究都集中在不同地點(diǎn)環(huán)境微生物組的生物多樣性,及其對(duì)環(huán)境波動(dòng)(特別是全球變化和環(huán)境污染)和時(shí)空變化的響應(yīng)方面,但很少有針對(duì)環(huán)境微生物組的功能、代謝活性、相互作用、進(jìn)化和潛在機(jī)制的研究。未來(lái)我們?cè)诮鉀Q微生物生態(tài)學(xué)和環(huán)境科學(xué)方面的實(shí)際問(wèn)題時(shí)必須重視以下6個(gè)方面的研究?jī)?nèi)容:① 微生物組在環(huán)境中的功能是什么?② 環(huán)境微生物多樣性及功能的形成和維持的機(jī)制是什么?③ 微生物群落中的各微生物如何相互作用?④ 關(guān)鍵微生物、優(yōu)勢(shì)和稀有物種在生態(tài)系統(tǒng)穩(wěn)定性和生產(chǎn)力中扮演怎樣的角色?⑤ 微生物組、環(huán)境和生態(tài)系統(tǒng)之間如何相互關(guān)聯(lián)?⑥ 如何通過(guò)微生物數(shù)據(jù)預(yù)測(cè)生態(tài)系統(tǒng)功能?盡管當(dāng)代的測(cè)序技術(shù)和相關(guān)的生物信息工具已極大地提高了我們對(duì)環(huán)境微生物組的認(rèn)識(shí)能力,但相關(guān)的技術(shù)挑戰(zhàn)仍然存在,包括測(cè)序手段的進(jìn)一步發(fā)展、更高效的計(jì)算和分析軟件、更全面的環(huán)境基因和基因組數(shù)據(jù)庫(kù)以及環(huán)境微生物研究者更多的跨學(xué)科訓(xùn)練。
高通量測(cè)序和其他組學(xué)技術(shù)不斷革新了我們認(rèn)識(shí)各種生境如土壤、沉積物、水及空氣中微生物的多樣性、組成、功能、相互作用、進(jìn)化和演替的能力。在過(guò)去10年間,這些研究為我們了解環(huán)境微生物組提供了新的見(jiàn)解。首先,我們對(duì)微生物多樣性的認(rèn)識(shí)得到了加強(qiáng);其次,新的微生物基因、基因組、微生物群體和通路的發(fā)現(xiàn)拓展了我們對(duì)環(huán)境中微生物功能的認(rèn)識(shí)。此外,科學(xué)家也開(kāi)始通過(guò)宏轉(zhuǎn)錄組和宏蛋白組分析研究原位微生物群落的活性。不僅如此,也有研究在微生物生態(tài)系統(tǒng)中發(fā)現(xiàn)了宏觀生態(tài)學(xué)中的一些時(shí)空變化規(guī)律,這能促使我們?nèi)ヌ剿魃锵到y(tǒng)中的普遍規(guī)律。通過(guò)上述的研究手段,我們將宏觀生態(tài)學(xué)及環(huán)境科學(xué)的實(shí)驗(yàn)策略、理論和假說(shuō)應(yīng)用于環(huán)境微生物組研究。最后,研究發(fā)現(xiàn)環(huán)境微生物組對(duì)導(dǎo)致微生物群落形成和調(diào)控其功能的全球環(huán)境變化十分敏感,相應(yīng)地,環(huán)境微生物也可能對(duì)生態(tài)系統(tǒng)功能產(chǎn)生反饋。
然而,我們?cè)趯?duì)環(huán)境微生物組的功能和機(jī)制,以及它們與環(huán)境驅(qū)動(dòng)力和生態(tài)系統(tǒng)功能的相關(guān)性的理解上仍然面臨許多挑戰(zhàn)。因此,未來(lái)環(huán)境微生物組研究的方向?qū)⒓杏冢孩?開(kāi)發(fā)新的測(cè)序方法和策略以及更高效的微生物組信息分析工具;② 探討不同生境中微生物功能和活性的機(jī)制;③ 了解關(guān)鍵微生物、優(yōu)勢(shì)和稀有物種在生態(tài)系統(tǒng)穩(wěn)定性和生產(chǎn)力中的作用;④ 理解微生物之間以及微生物和宿主間的相互作用;⑤ 通過(guò)新的分離方法和單細(xì)胞基因組學(xué)技術(shù)獲得更多的參考基因組用于以基因組為核心的研究;⑥ 利用微生物數(shù)據(jù)增強(qiáng)對(duì)生態(tài)系統(tǒng)功能的預(yù)測(cè)能力;⑦ 通過(guò)合成生物學(xué)和系統(tǒng)微生物學(xué)方法構(gòu)建理想的功能微生物組;⑧ 利用基因工程微生物提高環(huán)境健康。
[1] van der HEIJDEN M G A, BARDGETT R D, van STRAALEN N M. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems[J]. Ecology Letters, 2008,11(3): 296-310.
[2] BERENDSEN R L, PIETERSE C M, BAKKER P A. The rhizosphere microbiome and plant health[J]. Trends in Plant Science, 2012,17(8): 478-486.
[3] MENDES R, GARBEVA P, RAAIJMAKERS J M.The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms[J]. FEMS Microbiology Reviews, 2013,37(5): 634-663.
[4] BERG G, GRUBE M, SCHLOTER M, et al.Unraveling the plant microbiome: looking back and future perspectives[J]. Frontiers in Microbiology, 2014(5):148.
[5] HUNTER P. Plant microbiomes and sustainable agriculture[J]. EMBO Reports, 2016: e201643476-n/a.
[6] HAEGEMAN B, HAMELIN J, MORIARTY J,et al.Robust estimation of microbial diversity in theory and in practice[J]. The ISME Journal, 2013, 7(6): 1092-1101.
[7] LOCEY K J, LENNON J T. Scaling laws predict global microbial diversity[J]. Proceedings of the National Academy of Sciences, 2016,113(21): 5970-5975.
[8] TORSVIK V, OVREAS L.Microbial diversity and function in soil: from genes to ecosystems[J]. Current Opinion in Microbiology, 2002,5(3): 240-245.
[9] GANS J, WOLINSKY M, DUNBAR J. Computational improvements reveal great bacterial diversity and high metal toxicity in soil[J]. Science, 2005, 309(5739): 1387-1390.
[10] KALLMEYER J, POCKALNY R, ADHIKARI R R, et al. Global distribution of microbial abundance and biomass in subseafloor sediment[J]. Proceedings of the National Academy of Sciences, 2012, 109(40): 16213-16216.
[11] WHITMAN W B, COLEMAN D C, WIEBE W J. Prokaryotes: The unseen majority[J]. Proceedings of the National Academy of Sciences, 1998, 95(12): 6578-6583.
[12] HEATHER J M, CHAIN B. The sequence of sequencers: The history of sequencing DNA[J]. Genomics, 2016, 107(1): 1-8.
[13] WHITE III R A, CALLISTER S J, MOORE R J, et al.The past, present and future of microbiome analyses[J]. Nature Protocols, 2016,11(11): 2049-2053.
[14] VOSKOBOYNIK A, NEFF N F, SAHOO D, et al. The genome sequence of the colonial chordate, Botryllus schlosseri[J]. Elife, 2013, 2:e00569.
[15] van DIJK E L, AUGER H, JASZCZYSZYN Y, et al. Ten years of next-generation sequencing technology[J]. Trends in Genetics, 2014, 30(9): 418-426.
[16] MURRAY I A, CLARK T A, MORGAN R D, et al. The methylomes of six bacteria[J]. Nucleic Acids Research, 2012,40(22): 11450-11462.
[17] LAVER T, HARRISON J, O’NEILL P A, et al. Assessing the performance of the oxford nanopore technologies minion[J]. Biomolecular Detection and Quantification, 2015, 3: 1-8.
[18] TSAI Y C, CONLAN S, DEMING C, et al. Resolving the complexity of human skin metagenomes using single-molecule sequencing[J]. Mbio, 2016, 7(1): e01948-15.
[19] MUKHOPADHYAY A,HE Z, ALM E J, et al. Salt stress inDesulfovibriovulgarisHildenborough: an integrated genomics approach[J]. Journal of Bacteriology, 2006,188(11): 4068-4078.
[20] SCHOCH C L, SEIFERT K A, HUHNDORF S, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi[J]. Proceedings of the National Academy of Sciences, 2012, 109(16): 6241-6246.
[21] ZHOU J, DENG Y, SHEN L, et al. Temperature mediates continental-scale diversity of microbes in forest soils[J]. Nature Communications, 2016, 7: 12083.
[22] SANCHEZ O, FERRERA I, GONZALEZ J M, et al. Assessing bacterial diversity in a seawater-processing wastewater treatment plant by 454-pyrosequencing of the 16S rRNA and amoA genes[J]. Microbial Biotechnology, 2013, 6(4): 435-442.
[23] FISH J A, CHAI B, WANG Q, et al. FunGene: the functional gene pipeline and repository[J]. Frontiers in microbiology, 2013, 4: 291.
[24] SINGH A H, DOERKS T, LETUNIC I, et al. Discovering functional novelty in metagenomes: Examples from light-mediated processes[J]. Journal of Bacteriology, 2009, 191(1): 32-41.
[25] DUDHAGARA P, BHAVSAR S, BHAGAT C, et al. Web resources for metagenomics studies[J]. Genomics, Proteomics & Bioinformatics, 2015, 13(5): 296-303.
[26] HUSON D H, MITRA S.Introduction to the analysis of environmental sequences: metagenomics with MEGAN[J]. Methods in Molecular Biology, 2012, 856: 415-29.
[27] MEYER F, PAARMANN D, D’SOUZA M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes[J]. Bmc Bioinformatics, 2008, 9: 386.
[28] MARKOWITZ V M, IVANOVA N N, SZETO E, et al.IMG/M: a data management and analysis system for metagenomes[J]. Nucleic Acids Research, 2008,36: D534-D538.
[29] SUN S L, CHEN J, LI W Z, et al. Community cyberinfrastructure for advanced microbial ecology research and analysis: the CAMERA resource[J]. Nucleic Acids Research, 2011, 39: D546-D551.
[30] HUSON D H, BEIER S, FLADE I, et al. MEGAN community edition-interactive exploration and analysis of large-scale microbiome sequencing data[J]. Plos Computational Biology, 2016,12(6): e1004957.
[31] BUCHFINK B, XIE C, HUSON D H. Fast and sensitive protein alignment using DIAMOND[J]. Nat Meth, 2015, 12(1): 59-60.
[32] OVERBEEK R, BEGLEY T, Butler R M, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes[J]. Nucleic Acids Research, 2005, 33(17): 5691-5702.
[33] CAPORASO J G, BITTINGER K, BUSHMAN F D, et al. PyNAST: a flexible tool for aligning sequences to a template alignment[J]. Bioinformatics, 2010,26(2): 266-267.
[34] TU Q C, HE Z L, ZHOU J Z.Strain/species identification in metagenomes using genome-specific markers[J]. Nucleic Acids Research, 2014,42(8).
[35] CLEARY B, BRITO I L, HUANG K, et al.Detection of low-abundance bacterial strains in metagenomic datasets by eigengenome partitioning[J]. Nature Biotechnology, 2015, 33(10): 1053-1060.
[36] LUO C, KNIGHT R, SILJANDER H, et al. ConStrains identifies microbial strains in metagenomic datasets[J]. Nature Biotechnology, 2015, 33(10): 1045-1052.
[37] YE Y Z, TANG H X.Utilizing de Bruijn graph of metagenome assembly for metatranscriptome analysis[J]. Bioinformatics, 2016, 32(7): 1001-1008.
[38] LEUNG H C M, YIU S M, PARKINSON J, et al. IDBA-MT: De novo assembler for metatranscriptomic data generated from next-generation sequencing technology[J]. Journal of Computational Biology, 2013,20(7): 540-550.
[39] GRABHERR M G, HAAS B J, YASSOUR M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology, 2011, 29(7): 644-U130.
[40] SCHULZ M H, ZERBINO D R, VINGRON M, et al. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels[J]. Bioinformatics, 2012, 28(8): 1086-1092.
[41] LEUNG H C M, YIU S M, CHIN F Y L. IDBA-MTP: A hybrid metatranscriptomic assembler based on protein information[J]. Journal of Computational Biology, 2015, 22(5): 367-376.
[42] SCHMIEDER R, LIM Y W, EDWARDS R. Identification and removal of ribosomal RNA sequences from metatranscriptomes[J]. Bioinformatics, 2012,28(3): 433-435.
[43] LEE J H, YI H, CHUN J.rRNASelector: A computer program for selecting ribosomal RNA encoding sequences from metagenomic and metatranscriptomic shotgun libraries[J]. Journal of Microbiology, 2011, 49(4): 689-691.
[44] WANG Y,HU H, LI X. rRNAFilter: a fast approach for ribosomal RNA read removal without a reference database[J]. Journal of Computational Biology, 2017,24(4):368-375.
[45] KOPYLOVA E, NOE L, TOUZET H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data[J]. Bioinformatics, 2012,28(24): 3211-3217.
[46] SHI Y, TYSON G W, DELONG E F. Metatranscriptomics reveals unique microbial small RNAs in the ocean’s water column[J]. Nature, 2009, 459(7244): 266-269.
[47] NAWROCKI E P, EDDY S R. Infernal 1.1: 100-fold faster RNA homology searches[J]. Bioinformatics, 2013,29(22): 2933-2935.
[48] GRIFFITHS-JONES S, MOXON S, Marshall M, et al. Rfam: annotating non-coding RNAs in complete genomes[J]. Nucleic Acids Research, 2005,33: D121-D124.
[49] WESTREICH S T, KORF I, MILLS D A, et al. SAMSA: A comprehensive metatranscriptome analysis pipeline[J]. BMC Bioinformatics,2016,17:399.
[50] CURTIS T P, SLOAN W T, SCANNELL J W. Estimating prokaryotic diversity and its limits[J]. Proceedings of the National Academy of Sciences, 2002,99(16): 10494-10499.
[51] PAUL D, PANDEY G, PANDEY J, et al. Accessing microbial diversity for bioremediation and environmental restoration[J]. TRENDS in Biotechnology, 2005, 23(3): 135-142.
[52] BALVANERA P, PFISTERER A B, BUCHMANN N, et al. Quantifying the evidence for biodiversity effects on ecosystem functioning and services[J]. Ecology Letters, 2006, 9(10): 1146-1156.
[53] WAGG C, BENDER S F, WIDMER F, et al. Soil biodiversity and soil community composition determine ecosystem multifunctionality[J]. Proceedings of the National Academy of Sciences, 2014,111(14): 5266-5270.
[54] MENDES L W, TSAI S M, NAVARRETE A A, et al. Soil-borne microbiome: linking diversity to function[J]. Microbial Ecology, 2015,70(1): 255-265.
[55] NELSON M B, MARTINY A C, MARTINY J B. Global biogeography of microbial nitrogen-cycling traits in soil[J]. Proceedings of the National Academy of Sciences of the United States of America, 2016, 113(29): 8033-8040.
[56] MENDES R, KRUIJT M, de BRUIJN I, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria[J]. Science, 2011, 332(6033): 1097-1100.
[57] MENDES L W, KURAMAE E E, NAVARRETE A A, et al. Taxonomical and functional microbial community selection in soybean rhizosphere[J]. The ISME Journal, 2014, 8(8): 1577-1587.
[58] HE Z L, XIONG J B, KENT A D, et al. Distinct responses of soil microbial communities to elevated CO2and O3in a soybean agro-ecosystem[J]. The ISME Journal, 2014, 8(3): 714-726.
[59] REVENGA C, BRUNNER J, HENNINGER N, et al. Pilot analysis of global ecosystem:Grassland ecosystems[R]. World Resources Institute & International Food Policy Research Institute, 2000:103-118.
[60] GRAYSTON S J, GRIFFITH G S, MAWDSLEY J L, et al. Accounting for variability in soil microbial communities of temperate upland grassland ecosystems[J]. Soil Biology & Biochemistry, 2001, 33(4/5): 533-551.
[61] HU S, FIRESTONE M K, FIELD C B, et al. Nitrogen limitation of microbial decomposition in a grassland under elevated CO2[J]. Nature, 2001,409(6817): 188-191.
[62] XUE K, XIE J P, ZHOU A F, et al. Warming Alters Expressions of Microbial Functional Genes Important to Ecosystem Functioning[J]. Frontiers in Microbiology, 2016, 7:668.
[63] ZHOU J, XUE K, XIE J, et al. Microbial mediation of carbon-cycle feedbacks to climate warming[J]. Nature Climate Change, 2012,2(2): 106-110.
[64] DELMONT T O, PRESTAT E, KEEGAN K P, et al. Structure, fluctuation and magnitude of a natural grassland soil metagenome[J]. The ISME Journal, 2012, 6(9): 1677-1687.
[65] LUO C, RODRIGUEZ R L, JOHNSTON E R, et al. Soil microbial community responses to a decade of warming as revealed by comparative metagenomics[J]. Applied and Environmental Microbiology, 2014,80(5): 1777-86.
[66] TU Q, HE Z, WU L, et al. Metagenomic reconstruction of nitrogen cycling pathways in a CO2-enriched grassland ecosystem[J]. Soil Biology & Biochemistry, 2017,106: 99-108.
[67] UROZ S, BUEE M, DEVEAU A, et al. Ecology of the forest microbiome: Highlights of temperate and boreal ecosystems[J]. Soil Biology & Biochemistry, 2016, 103: 471-488.
[68] PAN Y D, BIRDSEY R A, FANG J Y, et al. A large and persistent carbon sink in the world’s forests[J]. Science, 2011, 333(6045): 988-993.
[69] ROSCH C, MERGEL A, BOTHE H. Biodiversity of denitrifying and dinitrogen-fixing bacteria in an acid forest soil[J]. Applied and Environmental Microbiology, 2002, 68(8): 3818-3829.
[70] WEIGOLD P, EL-HADIDI M, RUECKER A, et al. A metagenomic-based survey of microbial (de)halogenation potential in a German forest soil[J]. Scientific Reports, 2016,6:28958.
[71] CARDENAS E, KRANABETTER J M, HOPE G, et al. Forest harvesting reduces the soil metagenomic potential for biomass decomposition[J]. The ISME Journal, 2015,9(11): 2465-2476.
[72] HESSE C N, MUELLER R C, VUYISICH M, et al. Forest floor community metatranscriptomes identify fungal and bacterial responses to N deposition in two maple forests[J]. Frontiers in Microbiology, 2015: 6337.
[74] BONAN G B,III F S C, THOMPSON S L. Boreal forest and tundra ecosystems as components of the climate system[J]. Climatic Change, 1995,29(2): 145-167.
[75] TARNOCAI C, CANADELL J G, SCHUUR E A G, et al. Soil organic carbon pools in the northern circumpolar permafrost region[J]. Global Biogeochemical Cycles, 2009,23(2): 2607-2617.
[76] HEIMANN M, REICHSTEIN M.Terrestrial ecosystem carbon dynamics and climate feedbacks[J]. Nature, 2008, 451(7176): 289-292.
[77] McCALLEY C K, WOODCROFT B J, HODGKINS S B, et al. Methane dynamics regulated by microbial community response to permafrost thaw[J]. Nature, 2014,514(7523): 478-481.
[78] COOLEN M J L, ORSI W D. The transcriptional response of microbial communities in thawing Alaskan permafrost soils[J]. Frontiers in Microbiology, 2015, 6:197.
[79] MACKELPRANG R, WALDROP M P, DEANGELIS K M, et al. Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw[J]. Nature, 2011, 480(7377): 368-U120.
[80] HULTMAN J, WALDROP M P, MACKELPRANG R, et al. Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes[J]. Nature, 2015, 521(7551): 208-212.
[81] AZAM F. Microbial control of oceanic carbon flux: the plot thickens[J]. Science, 1998, 280(5364): 694-696.
[82] KARNER M B, DELONG E F, KARL D M. Archaeal dominance in the mesopelagic zone of the Pacific Ocean[J]. Nature, 2001,409(6819): 507-510.
[83] DUCKLOW H. Bacterial production and biomass in the oceans[M]// KIRCHMAN D L, ed. Microbial ecology of the oceans.New York: Wiley,2000:85-120.
[84] DAS S, LYLA P S, KHAN S A. Marine microbial diversity and ecology: importance and future perspectives[J]. Current Science, 2006,90(10): 1325-1335.
[85] FUHRMAN J A. Microbial community structure and its functional implications[J]. Nature, 2009, 459(7244): 193-199.
[86] STROM S L. Microbial ecology of ocean biogeochemistry: a community perspective[J]. Science, 2008,320(5879): 1043-1045.
[87] FINLAY B J, CLARKE K J.Ubiquitous dispersal of microbial species[J]. Nature, 1999, 400(6747): 828-828.
[88] FINLAY B J. Global dispersal of free-living microbial eukaryote species[J]. Science, 2002, 296(5570): 1061-1063.
[89] EILER A, HAYAKAWA D H, RAPPE M S. Non-random assembly of bacterioplankton communities in the subtropical North Pacific Ocean[J]. Frontiers in Microbiology, 2011, 2: 140.
[90] LANGENHEDER S, BERGA M, OSTMAN O, et al. Temporal variation of β-diversity and assembly mechanisms in a bacterial metacommunity[J]. The ISME Journal, 2012, 6(6): 1107-1114.
[91] FUHRMAN J A, HEWSON I, SCHWALBACH M S, et al. Annually reoccurring bacterial communities are predictable from ocean conditions[J]. Proceedings of the National Academy of Sciences of the United States of America, 2006,103(35): 13104-13109.
[92] LOGARES R, LINDSTROM E S, LANGENHEDER S, et al. Biogeography of bacterial communities exposed to progressive long-term environmental change[J]. The ISME Journal, 2013,7(5): 937-948.
[93] FUHRMAN J A, STEELE J A, HEWSON I, et al.A latitudinal diversity gradient in planktonic marine bacteria[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008,105(22): 7774-7778.
[94] CHOW C-E T, SACHDEVA R, CRAM J A, et al. Temporal variability and coherence of euphotic zone bacterial communities over a decade in the Southern California Bight[J]. The ISME Journal, 2013, 7(12): 2259-2273.
[95] CRAM J A, CHOW C-E T, SACHDEVA R, et al. Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years[J]. The ISME Journal, 2015, 9(3): 563-580.
[96] RINKE C, SCHWIENTEK P, SCZYRBA A, et al.Insights into the phylogeny and coding potential of microbial dark matter[J]. Nature, 2013, 499(7459): 431-437.
[97] GAWAD C, KOH W, QUAKE S R.Single-cell genome sequencing: current state of the science[J]. Nature Reviews Genetics, 2016, 17(3): 175-188.
[98] THOMPSON L R, FIELD C, ROMANUK T, et al. Patterns of ecological specialization among microbial populations in the Red Sea and diverse oligotrophic marine environments[J]. Ecology and Evolution, 2013, 3(6): 1780-1797.
[99] SUNAGAWA S, COELHO L P, CHAFFRON S, et al.Structure and function of the global ocean microbiome[J]. Science, 2015,348(6237): 1261359.
[100] THOMPSON L R, WILLIAMS G J, HAROON M F, et al. Metagenomic covariation along densely sampled environmental gradients in the Red Sea[J]. The ISME Journal, 2017,11(1): 138-151.
[101] OKAFOR N. Ecology of microorganisms in freshwater[J]. Environmental Microbiology of Aquatic and Waste Systems, 2011: 111-122.
[102] MARTINEZ-GARCIA M, SWAN B K, POULTON N J, et al. High-throughput single-cell sequencing identifies photoheterotrophs and chemoautotrophs in freshwater bacterioplankton[J]. The ISME Journal, 2012, 6(1): 113-123.
[103] MOU X, LU X, JACOB J, et al. Metagenomic identification of bacterioplankton taxa and pathways involved in microcystin degradation in Lake Erie[J]. PLoS ONE, 2013, 8(4): e61890.
[104] STALEY C, GOULD T J, WANG P, et al. Core functional traits of bacterial communities in the Upper Mississippi River show limited variation in response to land cover[J]. Frontiers in Microbiology, 2014,5: 414.
[105] FORTUNATO C S, CRUMP B C. Microbial gene abundance and expression patterns across a river to ocean salinity gradient[J]. PLoS One, 2015, 10(11): e0140578.
[106] van ROSSUM T, PEABODY M A, UYAGUARI-DIAZ M I, et al. Year-long metagenomic study of river microbiomes across land use and water quality[J]. Frontiers in Microbiology, 2015. 6:1405.
[107] YAN Q, BI Y, DENG Y, et al. Impacts of the Three Gorges Dam on microbial structure and potential function[J]. Scientific Reports, 2015, 5:8605.
[108] HEMME C L, TU Q, SHI Z, et al.Comparative metagenomics reveals impact of contaminants on groundwater microbiomes[J]. Frontiers in Microbiology, 2015,6:1205.
[109] HEMME C L, DENG Y, GENTRY T J, et al. Metagenomic insights into evolution of a heavy metal-contaminated groundwater microbial community[J]. The ISME Journal, 2010,4(5): 660-672.
[110] LENSKI R E. Experimental evolution and the dynamics of adaptation and genome evolution in microbial populations[J]. The ISME Journal, 2017.DOI:10.1038/ismej.2017.69.
[111] HEMME C L, GREEN S J, RISHISHWAR L, et al. Lateral gene transfer in a heavy metal-contaminated-groundwater microbial community[J]. MBio, 2016,7(2):e02234-15.
[112] ORTIZ G, YAGUEE G, SEGOVIA M, et al. A study of air microbe levels in different areas of a hospital[J]. Current Microbiology, 2009, 59(1): 53-58.
[113] YOO K, LEE T K, CHOI E J, et al. Molecular approaches for the detection and monitoring of microbial communities in bioaerosols: a review[J]. Journal of Environmental Sciences, 2016,51:234-247.
[114] BRODIE E L, DESANTIST Z, PARKER J P M, et al. Urban aerosols harbor diverse and dynamic bacterial populations[J]. Proceedings of the National Academy of Sciences of the United States of America, 2007,104(1): 299-304.
[115] ADAMS R I, BATEMAN A C, Bik H M, et al. Microbiota of the indoor environment: a meta-analysis[J]. Microbiome, 2015, 3(1): 1-18.
[116] PRUSSIN A J, MARR L C. Sources of airborne microorganisms in the built environment[J]. Microbiome, 2015, 3: 78.
[117] JIANG W J, LIANG P, WANG B Y, et al. Optimized DNA extraction and metagenomic sequencing of airborne microbial communities[J]. Nature Protocols, 2015,10(5): 768-779.
Progressesinenvironmentalmicrobiomediversityandfunctionresearch
WUBo1,FENGKai2,3,ZHIXiaoyang4,HEQiang1,5,XUMeiying6,DENGYe2,XIAOFanshu1,WANGShanquan1,YULing1,LUQihong1,LIANYingli1,LUOLijuan7,YUANKe8,CHENBaowei8,YANQingyun1,QIURongliang1,LUANTiangang7,HEZhili1,9
(1. Environmental Microbiome Research Center and the School of Environmental Science and Engineering, Sun Yat-sen University, Guangzhou 510006, China; 2. Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China; 3.University of Chinese Academy of Sciences, Beijing 100049, China; 4.Yunnan Institute of Microbiology, School of Life Science, Yunnan University, Kunming 650091, China; 5.Department of Civil and Environmental Engineering, The University of Tennessee, Knoxville, TN, 37996 ,USA; 6.Guangdong Institute of Microbiology, Guangzhou 510070, China;7.School of Marine Science, Sun Yat-sen University, Guangzhou 510275, China;8.School of Life Science, Sun Yat-sen University, Guangzhou 510275, China;9. Institute for Environmental Genomics, University of Oklahoma, Norman, OK, 73019, USA)
Currently developed sequencing and other high throughput technologies enable us to explore the role of microorganisms in the environment. Here environmental microbiome is defined as the collection of microbial genomes in a particular environment, and the environmental microbiome research (or environmental microbiomics) aims to understand the biodiversity, composition, structure, function, interaction, evolution and dynamics of microbial communities as well as their linkages with environmental drivers and ecosystem functioning. This review is focused on environmental microbiomes from soil, water and air. Firstly, the current advanced sequencing technologies and associated tools for environmental microbiome data analysis are reviewed, then the current research progress in environmental microbiomes were summarized with foci on addressing key scientific questions in microbial ecology and environmental science. This review proposes key scientific questions, and highlights the importance of our understanding of microbial functions and underlying mechanisms in future environmental microbiome studies.
environmental microbiomes; microbial function; metagenome sequencing; bioinformatics tools
X172
A
0529-6579(2017)05-0001-11
10.13471/j.cnki.acta.snus.2017.05.001
2017-07-18
中山大學(xué)千人計(jì)劃項(xiàng)目(38000-18821105);中山大學(xué)百人計(jì)劃項(xiàng)目(38000-18821107)
吳波(1986年生),女;研究方向微生物生態(tài);E-mail: wubo28@mail.sysu.edu.cn 馮凱(1992年生),男;研究方向微生物生態(tài)和生物信息;E-mail: kaifeng_st@rcees.ac.cn(以上2人并列第1作者)
賀志理(1965年生),男;研究方向:微生物組學(xué);E-mail: hezhili@mail.sysu.edu.cn