未培養(yǎng)技術(shù)在瘤胃產(chǎn)甲烷菌群研究中的應(yīng)用

王興文王加啟趙圣國(guó)李發(fā)弟卜登攀

（1. 中國(guó)農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所 動(dòng)物營(yíng)養(yǎng)學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100193；2. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，蘭州 730070）

未培養(yǎng)技術(shù)在瘤胃產(chǎn)甲烷菌群研究中的應(yīng)用

王興文1，2王加啟1趙圣國(guó)1李發(fā)弟2卜登攀1

（1. 中國(guó)農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所 動(dòng)物營(yíng)養(yǎng)學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100193；2. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，蘭州 730070）

瘤胃中的產(chǎn)甲烷菌為嚴(yán)格厭氧微生物，很難通過(guò)分離培養(yǎng)的方法進(jìn)行研究。現(xiàn)已經(jīng)培養(yǎng)并公布的瘤胃產(chǎn)甲烷菌只有4種，成千上萬(wàn)種類的瘤胃產(chǎn)甲烷菌不能培養(yǎng)出來(lái)。未培養(yǎng)技術(shù)的發(fā)展與應(yīng)用突破了傳統(tǒng)分離培養(yǎng)檢測(cè)方法的限制，使瘤胃產(chǎn)甲烷菌的研究有了較大的進(jìn)展。綜述了實(shí)時(shí)定量PCR、16S rDNA文庫(kù)、DGGE和高通量測(cè)序技術(shù)等6項(xiàng)具有代表性的未培養(yǎng)技術(shù)及其在瘤胃產(chǎn)甲烷菌研究中的應(yīng)用，為以后瘤胃產(chǎn)甲烷菌的研究提供方法上的參考。

瘤胃 產(chǎn)甲烷菌 未培養(yǎng)技術(shù) 全基因組測(cè)序 高通量測(cè)序

反芻動(dòng)物瘤胃中甲烷的生成會(huì)造成飼料能量的損失，同時(shí)甲烷也是造成溫室效應(yīng)的主要?dú)怏w之一。它對(duì)溫室效應(yīng)的作用是CO2的25倍，全球每年由反芻動(dòng)物排放的甲烷達(dá)8.1×107-9.2×107t，這相當(dāng)于人類制造的甲烷的23%-27%［1］。產(chǎn)甲烷菌不足瘤胃微生物總量的4%，但幾乎產(chǎn)生了反芻動(dòng)物釋放的全部甲烷［2］。將產(chǎn)甲烷菌的活性和數(shù)量控制在最適水平，達(dá)到反芻動(dòng)物最佳的生產(chǎn)性能而對(duì)環(huán)境不造成影響是目前研究的重要方向。產(chǎn)甲烷菌是嚴(yán)格厭氧的古菌，一些在體內(nèi)可以發(fā)現(xiàn)的產(chǎn)甲烷菌在體外很難分離培養(yǎng)。所以到目前為止，只有很少種類的產(chǎn)甲烷菌從瘤胃中分離出來(lái)，已經(jīng)分離培養(yǎng)出來(lái)的瘤胃產(chǎn)甲烷菌主要為甲酸甲烷桿菌（Methanobacterium formicicum）、布氏甲烷桿菌（Methanobacterium byantii）、甲烷短桿菌（Methanobrevibacter ruminantium）和可移動(dòng)的甲烷微菌（Methanomicrobium mobile）［3］。隨著分子生物學(xué)的發(fā)展，未培養(yǎng)技術(shù)已越來(lái)越多地應(yīng)用到瘤胃產(chǎn)甲烷菌的研究中，通過(guò)未培養(yǎng)技術(shù)可以不需要產(chǎn)甲烷菌的分離培養(yǎng)就可以準(zhǔn)確地研究瘤胃中的產(chǎn)甲烷菌性質(zhì)。本文主要介紹幾種未培養(yǎng)技術(shù)以及其在產(chǎn)甲烷菌研究上的應(yīng)用。

1 產(chǎn)甲烷菌多樣性

1.1 產(chǎn)甲烷菌的引物

產(chǎn)甲烷菌研究中目的基因主要為16S rDNA和mcrA基因，16S rDNA是編碼核糖體RNA（rRNA）分子的對(duì)應(yīng)的DNA序列，既具有保守性，又具有高變性。保守性能夠反應(yīng)出生物物種的親緣關(guān)系，為系統(tǒng)發(fā)育重建提供線索；高變性則能揭示出生物物種的特征核酸序列。mcrA基因編碼甲基輔酶M還原酶（MCR）的α亞基，甲基輔酶M還原酶是甲烷產(chǎn)生過(guò)程中的一種關(guān)鍵酶，這種酶將與它的連接的甲基還原為CH4［4］，它是產(chǎn)甲烷菌所獨(dú)有的目的基因。16S rDNA的常見(jiàn)引物已在表1中列出，mcrA基因的常見(jiàn)引物已在表2中列出。

1.2 變性梯度凝膠電泳（Denaturing gradient gel electrophoresis，DGGE）

DGGE技術(shù)以DNA指紋圖譜為基礎(chǔ)，利用了DNA分子從雙螺旋型到局部變性時(shí)電泳遷移率會(huì)急劇下降的原理。在得到的圖譜中每一個(gè)條帶都代表著某個(gè)微生物優(yōu)勢(shì)菌群。

Zhou等［26］以16S rDNA為目的基因利用DGGE研究?jī)煞N不同日糧：生長(zhǎng)期的日糧（高能量）和育肥期的日糧（低能量）對(duì)瘤胃產(chǎn)甲烷菌的影響。從DGGE的圖譜中發(fā)現(xiàn)，兩種日糧的優(yōu)勢(shì)條帶不同，測(cè)序發(fā)現(xiàn)低能量日糧主要以Methanobrevibacter ruminantiumNT7為主，而高精料日料含有Methanobrevibacter smithii、Methanobrevibactersp. AbM4和M.ruminantiumNT7。Popova等［34］用高精料日糧和低精料日糧飼喂奶牛，以mcrA基因作為目的基因用DGGE的方法研究產(chǎn)甲烷菌的多樣性，發(fā)現(xiàn)高精料日糧飼喂的奶牛產(chǎn)甲烷菌的多樣性明顯低于高粗料日糧組，日糧中的細(xì)胞壁成分的含量與瘤胃內(nèi)的產(chǎn)甲烷菌的多樣性成正比。Monica等［35］利用DGGE觀察馴鹿兩個(gè)時(shí)期（即馴鹿采食量的高峰期和低谷期）瘤胃產(chǎn)甲烷菌多樣性發(fā)現(xiàn)，兩組沒(méi)有較大差別，說(shuō)明采食量不會(huì)引起馴鹿瘤胃微生物的多樣性的變化。

Zhou等［36］在奶牛飼糧中添加外源纖維素酶，甲烷生成量增加，為了解外源纖維分解酶對(duì)瘤胃產(chǎn)甲烷菌多樣性和數(shù)量的影響以及甲烷產(chǎn)量和產(chǎn)甲烷菌群落結(jié)構(gòu)的關(guān)系，應(yīng)用實(shí)時(shí)定量PCR和DGGE技術(shù)，發(fā)現(xiàn)瘤胃產(chǎn)甲烷菌的數(shù)量沒(méi)有改變，但是瘤胃產(chǎn)甲烷菌的多樣性發(fā)生了改變。Cheng等［37］用3對(duì)引物（344fGC/915r、519f/915rGC和1106fGC/1378r）作DGGE研究中國(guó)西部綿羊產(chǎn)甲烷菌的多樣性發(fā)現(xiàn)，引物519f/915rGC產(chǎn)生的條帶分辨率最高。

DGGE直接利用微生物樣品中的DNA或RNA對(duì)微生物群體加以區(qū)別，能解決微生物無(wú)法進(jìn)行體外培養(yǎng)的難題，可判斷出微生物間的種屬關(guān)系，甚至還可以發(fā)現(xiàn)新的微生物。但是DGGE技術(shù)也存在一些缺點(diǎn)，如DNA的分離技術(shù)不同，細(xì)胞壁的破裂程度不同，DNA的回收率不同以及PCR不均等擴(kuò)增都會(huì)造成試驗(yàn)結(jié)果的不同。此外，當(dāng)微生物的DNA量不到樣品DNA量的1%時(shí)，很難檢測(cè)其條帶。有一些產(chǎn)甲烷菌的基因組中含有多個(gè)目的基因的拷貝，就會(huì)過(guò)高的估計(jì)產(chǎn)甲烷菌的數(shù)量。

1.3 16S rDNA/mcrA基因建庫(kù)

構(gòu)建16S rDNA/mcrA文庫(kù)可以對(duì)某一特定環(huán)境中甲烷菌群落結(jié)構(gòu)和多樣性進(jìn)行分析。用PCR技術(shù)擴(kuò)增特異的產(chǎn)甲烷菌的基因，然后將這種特異的基因與載體連接，轉(zhuǎn)入大腸桿菌形成文庫(kù)，選取陽(yáng)性克隆測(cè)序，在測(cè)序的結(jié)果中去除載體和嵌合體序列，得到的序列與GenBank/RDP數(shù)據(jù)庫(kù)比對(duì)，就可以知道每個(gè)克隆中的片段屬于那一種微生物。

Franzolin等［38］研究不同粗飼料類型（牧草、玉米秸稈和甘蔗）對(duì)瘤胃產(chǎn)甲烷菌的影響，經(jīng)16S rDNA建庫(kù)后，得到467個(gè)克隆，包含19個(gè)操作分類單元：其中有4個(gè)在3個(gè)文庫(kù)中都存在，8個(gè)是特異的，6個(gè)存在于玉米秸稈組和牧草組，1個(gè)存在于甘蔗組和牧草組，由此發(fā)現(xiàn)不同粗飼料類型會(huì)改變瘤胃產(chǎn)甲烷菌的多樣性。Huang等［39］為了解牦牛甲烷產(chǎn)量少于青藏高原牛的原因，構(gòu)建牦牛和青藏高原牛瘤胃產(chǎn)甲烷菌16S rDNA文庫(kù)，在牦牛的文庫(kù)中發(fā)現(xiàn)了209個(gè)克隆，青藏高原牛文庫(kù)中205個(gè)克隆，將相似度在98%及以上歸為一個(gè)操作分類單元，共發(fā)現(xiàn)了95個(gè)操作分類單元：牦牛單獨(dú)含有46個(gè)操作分類單元，青藏高原牛含有34個(gè)操作分類單元，在兩個(gè)文庫(kù)中都存在的操作分類單元有15個(gè)。牦牛瘤胃的產(chǎn)甲烷多樣性要高于青藏高原牛。

表1 產(chǎn)甲烷菌16S rDNA引物［5］

表2 產(chǎn)甲烷菌mcrA基因引物［5］

裴彩霞等［40］采用3對(duì)特異性引物擴(kuò)增瘤胃古菌16S rRNA基因分別建立克隆文庫(kù)研究晉南牛瘤胃產(chǎn)甲烷菌的多樣性，引物Archf364/r1386建立的克隆庫(kù)中，克隆分為4類，分別與Methanobrebacter. ap.1Y（61%）、SM9（23%）、NT7（14%）和AK-87（2%）相似。引物1Af/1100Ar建立的克隆庫(kù)中，克隆分為兩類，分別與Methanobacterium aarhusense（72%）和Methanosphaera stadtmanaeDSM3091（28%）相似。引物Met86F/Met1340R建立的克隆文庫(kù)反映的古菌種類較為全面，除4種甲烷短桿菌（所占比例分別為47%、26%、11%和3%）外，還有Methanomicrobium mobile（2%）、以及類似Methanobacterium aarhusense（1%）和Methanosphaera stadtmanae（3%）的序列，還有7%的未匹配序列。Tymensen等［41］分別利用16S rDNA和mcrA建立文庫(kù)研究瘤胃液和附著在原蟲(chóng)上的產(chǎn)甲烷菌的多樣性，結(jié)果表明瘤胃液中產(chǎn)甲烷菌的多樣性要顯著高于附著在原蟲(chóng)上的產(chǎn)甲烷菌。Hook等［42］通過(guò)16S rDNA建庫(kù)和DGGE的方法研究莫能菌素對(duì)瘤胃產(chǎn)甲烷菌多樣性的影響發(fā)現(xiàn)，莫能菌素不會(huì)影響產(chǎn)甲烷菌的多樣性。

建立基因文庫(kù)較為準(zhǔn)確，可以精確地了解瘤胃內(nèi)產(chǎn)甲烷菌的系統(tǒng)分類。在實(shí)際研究中常常與DGGE相結(jié)合。但是要得到全面準(zhǔn)確的數(shù)據(jù)，就要得到足夠的克隆，挑取克隆的過(guò)程和步驟較為繁瑣，測(cè)序的費(fèi)用也相當(dāng)高。而且即使有相當(dāng)量的克隆也有可能遺漏一些豐度較低的基因。

1.4 限制性片段長(zhǎng)度多態(tài)性（Restriction fragment length polymorphism，RFLP）技術(shù)

RFLP（T-RFLP）技術(shù)通過(guò)檢測(cè)DNA在限制性內(nèi)切酶酶切后形成的特定DNA片段的大小以識(shí)別限制性片段長(zhǎng)度的多態(tài)性?，F(xiàn)在已經(jīng)發(fā)展出末端限制性片段長(zhǎng)度多態(tài)性技術(shù)等，T-RFLP在技術(shù)上與RFLP相似，只是其中的一個(gè)引物的5'端用熒光物質(zhì)標(biāo)記，將PCR產(chǎn)物進(jìn)行特異的酶切后就產(chǎn)生了許多不同長(zhǎng)度的限制性片段，其中帶有熒光標(biāo)記的目標(biāo)片段就可以被DNA自動(dòng)測(cè)序儀檢測(cè)到。

Frey等［43］為了了解泌乳牛瘤胃、十二指腸、回腸和糞便中產(chǎn)甲烷菌的多樣性，應(yīng)用RFLP技術(shù)研究發(fā)現(xiàn)在4個(gè)部位的產(chǎn)甲烷菌多樣性差異不顯著。Danielsson等［44］應(yīng)用T-RFLP技術(shù)和16S rDNA建庫(kù)結(jié)合研究不同日糧精粗比例對(duì)瘤胃產(chǎn)甲烷菌的影響，發(fā)現(xiàn)日糧精粗比與瘤胃產(chǎn)甲烷菌的多樣性成反比。

1.5 高通量測(cè)序

高通量測(cè)序是通過(guò)對(duì)小亞基核糖體（如16S rRNA）基因高變區(qū)進(jìn)行測(cè)序快速的檢測(cè)樣本中產(chǎn)甲烷菌的多樣性。在轉(zhuǎn)錄組學(xué)研究中，RNA-seq直接定量檢測(cè)瘤胃產(chǎn)甲烷菌的轉(zhuǎn)錄出來(lái)的mRNA的數(shù)量，從而了解產(chǎn)甲烷菌中代謝所需酶的活性［45，46］。

Fout等［47］利用高通量測(cè)序來(lái)檢測(cè)瘤胃產(chǎn)甲烷菌多樣性，得到138個(gè)序列，97%相似度序列歸類分析，發(fā)現(xiàn)了20個(gè)OTU。其中18個(gè)OTU為數(shù)據(jù)庫(kù)中已存在的序列，2個(gè)OTU為新菌，分別與Methanobrevibacter和Thermogymnomonas最為相似。

Hristov等［48］利用高通量測(cè)序檢測(cè)日糧中添加不飽和脂肪酸對(duì)瘤胃中產(chǎn)甲烷菌的影響，發(fā)現(xiàn)添加不飽和脂肪酸可以明顯增加Methanosphaera的數(shù)量而減少M(fèi)ethanobrevibacter的數(shù)量。

Lee等［49］研究韓國(guó)本地1種羊和3種牛的瘤胃產(chǎn)甲烷菌群落結(jié)構(gòu)，高通量測(cè)序發(fā)現(xiàn)在4種動(dòng)物的體內(nèi)都不存在未知的產(chǎn)甲烷菌，都可以在數(shù)據(jù)庫(kù)中找到。

2 產(chǎn)甲烷菌的定量分析

實(shí)時(shí)定量PCR（Real-time quantitative PCR，RTPCR）是在PCR的體系中加入熒光標(biāo)記物后，模板擴(kuò)增產(chǎn)物的量與熒光信號(hào)的強(qiáng)弱成正比，即每個(gè)模板的CT值與該模板的起始拷貝數(shù)的對(duì)數(shù)存在線性關(guān)系。利用已知拷貝數(shù)的標(biāo)準(zhǔn)品可做出標(biāo)準(zhǔn)曲線，只要得到未知樣品的CT值，就能算出起始拷貝數(shù)。

Zhou等［25］以16S rDNA為目的基因用3種不同的特異性引物AbM4-F/AbM4-R，Stad-F/Stad-R和uniMet1-F/uniMet1-R分別定量測(cè)定Methanobrevibactersp. strain AbM4，M. stadtmanae和Total methanogens來(lái)研究日糧消化效率對(duì)牛瘤胃產(chǎn)甲烷菌數(shù)量的影響，發(fā)現(xiàn)消化率對(duì)總產(chǎn)甲烷菌的數(shù)量沒(méi)有影響，但是高消化率組的Methanosphaera stadtmanae和Methanobrevibactersp. strain AbM4的數(shù)量是低消化組的1.92和2.26倍。Tymensen等［50］利用引物NestmetF/ NestmetR，Nest-MbbF/NestMbbR，NestRCCF/NestRCCR和NestMmF/NestMmR分別擴(kuò)增總甲烷菌Methanobrevibacterspp.，RCC和Methanomicrobiumspp.來(lái)研究與原蟲(chóng)共生的產(chǎn)甲烷菌，發(fā)現(xiàn)Methanobrevibacterspp.和RCC是與原蟲(chóng)連接的主要產(chǎn)甲烷菌。Hook等［51］以mcrA基因?yàn)槟康幕驊?yīng)用實(shí)時(shí)定量PCR和建庫(kù)研究能量高低不同的日糧對(duì)奶牛瘤胃內(nèi)的產(chǎn)甲烷菌的影響發(fā)現(xiàn)，日糧能量的高低對(duì)瘤胃產(chǎn)甲烷菌的影響不大。Denman等［30］用mcrA基因作實(shí)時(shí)定量PCR研究溴氯甲烷對(duì)瘤胃產(chǎn)甲烷菌的影響發(fā)現(xiàn)，日糧中添加溴氯甲烷可顯著減少瘤胃中產(chǎn)甲烷菌的數(shù)量。David等［52］利用mcrA基因定量總產(chǎn)甲烷菌發(fā)現(xiàn)飼喂干草的綿羊瘤胃內(nèi)的產(chǎn)甲烷菌要多于飼喂混合日糧的產(chǎn)甲烷菌。

實(shí)時(shí)定量PCR技術(shù)能夠準(zhǔn)確地定量研究瘤胃中產(chǎn)甲烷菌，但是要研究每一種產(chǎn)甲烷菌的數(shù)量必須得到相應(yīng)的引物，但因一些產(chǎn)甲烷菌的引物還沒(méi)有得到或在PCR過(guò)程中效率太低或特異性不好，因此不能定量研究。

3 產(chǎn)甲烷菌基因組

全基因組測(cè)序是對(duì)未知基因組序列的物種進(jìn)行個(gè)體的基因組測(cè)序，從整體上研究微生物的功能和組成。在新西蘭、澳大利亞和韓國(guó)正在進(jìn)行著瘤胃產(chǎn)甲烷菌的基因組測(cè)序（表3）。

到目前為止，已有2種產(chǎn)甲烷菌測(cè)序完成并公布，第一個(gè)為Methanobrevibacter ruminantiumM1［53］，它是由Bryant在1965年分離得到的，全長(zhǎng)2.93 Mb。通過(guò)全基因組測(cè)序在M1的基因組中發(fā)現(xiàn)前噬菌體的基因組、分解酶基因以及內(nèi)異肽酶PeiR基因，這些基因與M1的代謝息息相關(guān)。此外，還在M1的基因組中還發(fā)現(xiàn)了非核糖體肽合成酶基因，這是此基因在古菌基因組中第一次發(fā)現(xiàn)，從而揭示一些未知的代謝機(jī)理。另一個(gè)測(cè)序完成的產(chǎn)甲烷菌為Methanobrevibactersp. JHI，它是在韓國(guó)土牛瘤胃中分離得到的［54］，通過(guò)基因組測(cè)序發(fā)現(xiàn)JHI的基因組與Methanobrevibactersp. AbM4（從新西蘭羊皺胃分離，并測(cè)序完成，未發(fā)表）相似，說(shuō)明它們?cè)谙到y(tǒng)發(fā)育學(xué)上屬于同一種類。但是在JHI的基因組中發(fā)現(xiàn)了前噬菌體基因，而AbM4沒(méi)有，說(shuō)明兩者在進(jìn)化上出現(xiàn)了分支。全基因組測(cè)序在JHI的基因組中也發(fā)現(xiàn)了與甲烷生成利用的H2、CO2和甲酸鹽相關(guān)的酶。

通過(guò)全基因組測(cè)序可以清楚地了解目的產(chǎn)甲烷菌的代謝過(guò)程和系統(tǒng)發(fā)育的位置。在此基礎(chǔ)之上，應(yīng)用免疫方法和代謝阻遏物可以控制代謝過(guò)程中酶的活性來(lái)減少甲烷的排放［55］。但是全基因組測(cè)序是一項(xiàng)很大的工程，而且只有分離培養(yǎng)出來(lái)的產(chǎn)甲烷菌才能進(jìn)行全基因組測(cè)序，成功分離培養(yǎng)的產(chǎn)甲烷菌只占瘤胃產(chǎn)甲烷菌很小的一部分。

表3 瘤胃產(chǎn)甲烷菌基因組測(cè)序計(jì)劃［56］

在實(shí)際的研究過(guò)程中，DGGE技術(shù)、16S rDNA/

mcrA基因建庫(kù)、RFLP技術(shù)以及高通量測(cè)序主要用于研究產(chǎn)甲烷菌群落結(jié)構(gòu)，包括菌落的多樣性，豐富度。實(shí)時(shí)定量PCR主要用于定量產(chǎn)甲烷菌。若要了解產(chǎn)甲烷菌的功能和性質(zhì)則需要全基因組測(cè)序。有時(shí)，為了更加準(zhǔn)確具體地研究產(chǎn)甲烷菌，幾項(xiàng)技術(shù)可以同時(shí)運(yùn)用，克服彼此的缺點(diǎn)。如在一些研究中DGGE與16S rDNA/mcrA基因建庫(kù)聯(lián)合使用，DGGE可以發(fā)現(xiàn)主要的產(chǎn)甲烷菌，而16S rDNA/mcrA基因建庫(kù)可以提高覆蓋率全面地研究菌種的多樣性，兩者相互補(bǔ)充，互相驗(yàn)證。16S rDNA/mcrA基因建庫(kù)與定量PCR的聯(lián)合使用常常用于研究影響產(chǎn)甲烷菌的因素。

4 展望

未培養(yǎng)技術(shù)已經(jīng)應(yīng)用到了很多瘤胃產(chǎn)甲烷菌的研究中。隨著高通量測(cè)序的普及，基于高通量測(cè)序的多樣性研究是未來(lái)的發(fā)展趨勢(shì)，全基因組測(cè)序也是未來(lái)研究產(chǎn)甲烷菌代謝途徑和系統(tǒng)發(fā)育的重要內(nèi)容。利用這些技術(shù)可以了解瘤胃中產(chǎn)甲烷菌的群落結(jié)構(gòu)和生理狀態(tài)，為合理有效控制甲烷的生成提供理論依據(jù)。

［1］ IPCC（Intergovernmental Panel on Climate Change）, Climate Change 2007［M］. The Scientific Basis. Cambridge, UK：Cambridge University Press, 2007.

［2］ Lin C, Raskin L, Stahl DA. Microbial community structure in gastrointestinal tracts of domestic animals：comparative analyses using rRNA-targeted oligonucleotide probes［J］. FEMS Microbiol Ecol, 1997, 22：281-294.

［3］ Janssen PH, Kirs M. Structure of the archaeal community of the rumen［J］. Appl Environ Microbiol, 2008, 74（12）：3619-3625.

［4］ Ellermann J, Hedderich R, B?cher R, et al. The final step in methane formation. Investigations with highly purified methyl-CoM reductase（component C）fromMethanobacterium thermoautotrophicum（strain Marburg）［J］. Eur J Biochem, 1988, 172（3）：669-677.

［5］ Zhou M, McAllister TA, Guan LL. Molecular identification of rumen methanogens：Technologies, advances and prospects［J］. Anim Feed Sci Tech, 2011, 166-167：76-86.

［6］ Yanagita K, Kamagata Y, Kawaharasaki M, et al. Phylogenetic analysis of methanogens in sheep rumen ecosystem and detection ofMethanomicrobiummobile By fluorescencein situhybridization［J］. Biosci Biotechnol Biochem, 2000, 64（8）：1737-1742.

［7］ Whitford MF, Teather RM, Forster RJ. Phylogenetic analysis of methanogens from the bovine rumen［J］. BMC Microbiol, 2001, 1：5.

［8］ Skillman LC, Evans PN, Naylor GE, et al. 16S ribosomal DNA-directed PCR primers for ruminal methanogens and identification of methanogens colonising young lambs［J］. Anaerobe, 2004, 10（5）：277-285.

［9］ Tajima K, Nagamine T, Matsui H, et al. Phylogenetic analysis of archaeal 16S rRNA libraries from the rumen suggests the existence of a novel group of archaea not associated with known methanogens［J］. FEMS Microbiol Lett, 2001, 200（1）：67-72.

［10］ Wright AD, Williams AJ, Winder B, et al. Molecular diversity of rumen methanogens from sheep in Western Australia［J］. Appl Environ Microbiol, 2004, 70（3）：1263-1270.

［11］ Shin EC, Choi BR, Lim WJ, et al. Phylogenetic analysis of archaea in three fractions of cow rumen based on the 16S rDNA sequence［J］. Anaerobe, 2004, 10（6）：313-319.

［12］ Whitehead TR, Cotta MA. Phylogenetic diversity of methanogenic archaea in swine waste storage pits［J］. FEMS Microbiol Lett, 1999, 179（2）：223-226.

［13］ DeLong EF. Archaea in coastal marine environments［J］. Proc Natl Acad Sci USA, 1992, 89：5685-5689.

［14］ Tokura M, Chagan I, Ushida K, et al. Phylogenetic study of methanogens associated with rumen ciliates［J］. Curr Microbiol, 1999, 39（3）：123-128.

［15］ Marchesi JR, Weightman AJ, Cragg BA, et al. Methanogen and bacterial diversity and distribution in deep gas hydrate sediments from the Cascadia Margin as revealed by 16S rRNA molecular analysis［J］. FEMS Microbiol Ecol, 2001, 34（3）：221-228.

［16］ Embley TM, Finlay BJ, Thomas RH, et al. The use of rRNA sequences and fluorescent probes to investigate the phylogenetic positions of the anaerobic ciliateMetopus palaeformisand its archaeobacterial endosymbiont［J］. J Gen Microbiol, 1992, 138（7）：1479-1487.

［17］ Barns SM, Fundyga RE, Jeffries MW, et al. Remarkable archaeal diversity detected in a Yellowstone National Park hot spring environment［J］. Proc Natl Acad Sci USA, 1994, 91（5）：1609-1613.

［18］ Baker GC, Smith JJ, Cowan DA. Review and re-analysis of domainspecific 16S primers［J］. J Microbiol Methods, 2003, 55（3）：541-555.

［19］ Yu Z, García-González R, Schanbacher FL, et al. Evaluations of different hypervariable regions of archaeal 16S rRNA genes in profiling of methanogens by Archaea-specific PCR and denaturing gradient gel electrophoresis［J］. Appl Environ Microbiol, 2008, 74（3）：889-893.

［20］ Bano N, Ruffin S, Ransom B, et al. Phylogenetic composition of Arctic Ocean archaeal assemblages and comparison with Antarctic assemblages［J］. Appl Environ Microbiol, 2004, 70（2）：781-789.

［21］ Pinar G, Saiz-Jimenez C, Schabereiter-Gurtner C, et al. Archaeal communities in two disparate deteriorated ancient wall paintings：detection, identification and temporal monitoring by denaturing gradient gel electrophoresis［J］. FEMS Microbiol Ecol, 2001, 37：45-54.

［22］ Kleikemper J, Pombo SA, Schroth MH, et al. Activity and diversity of methanogens in a petroleum hydrocarbon-contaminated aquifer［J］. Appl Environ Microbiol, 2005, 71（1）：149-158.

［23］ Wright AD, Pimm C. Improved strategy for presumptive identification of methanogens using 16S riboprinting［J］. J Microbiol Methods, 2003, 55（2）：337-349.

［24］ van Hoek AH, van Alen TA, Sprakel VS, et al. Multiple acquisition of methanogenic archaeal symbionts by anaerobic ciliates［J］. Mol Biol Evol, 2000, 17（2）：251-258.

［25］ Amann RI, Ludwig W, Schleifer KH. Phylogenetic identification andin situdetection of individual microbial cells without cultivation［J］. Microbiol Rev, 1995, 59（1）：143-169.

［26］ Zhou M, Hernandez-Sanabria E, Guan LL. Assessment of the microbial ecology of ruminal methanogens in cattle with different feed efficiencies［J］. Appl Environ Microbiol, 2009, 75（20）：6524-6533.

［27］ Hales BA, Edwards C, Ritchie DA, et al. Isolation and identification of methanogen-specific DNA from blanket bog peat by PCR amplification and sequence analysis［J］. Appl Environ Microbiol, 1996, 62（2）：668-675.

［28］ Luton PE, Wayne JM, Sharp RJ, et al. The mcrA gene as an alternative to 16S rRNA in the phylogenetic analysis of methanogen populations in landfill［J］. Microbiology, 2002, 148（11）：3521-3530.

［29］ Hallam SJ, Girguis PR, Preston CM, et al. Identification of methyl coenzyme M reductase A（mcrA）genes associated with methaneoxidizing archaea［J］. Appl Environ Microbiol, 2003, 69（9）：5483-5491.

［30］ Denman SE, Tomkins NW, McSweeney CS. Quantitation and diversity analysis of ruminal methanogenic populations in response to the antimethanogenic compound bromochloromethane［J］. FEMS Microbiol Ecol, 2007, 62（3）：313-322.

［31］ Scanlan PD, Shanahan F, Marchesi JR. Human methanogen diversity and incidence in healthy and diseased colonic groups using mcrA gene analysis［J］. BMC Microbiol, 2008, 8：79.

［32］ Steinberg LM, Regan JM. Phylogenetic comparison of the methanogenic communities from an acidic, oligotrophic fen and an anaerobic digester treating municipal wastewater sludge［J］. Appl Environ Microbiol, 2008, 74（21）：6663-6671.

［33］ Juottonen H, Galand PE, Yrj?l? K. Detection of methanogenic Archaea in peat：comparison of PCR primers targeting the mcrA gene［J］. Res Microbiol, 2006, 157（10）：914-921.

［34］ Popova M, Martin C, Eugène M, et al. Effect of fibre and starchrich finishing diets on methanogenicArchaeadiversity and activity in the rumen of feedlot bull［J］. Anim Feed Sci Tech, 2011, 166-167：113-121.

［35］ Sundset MA, Edwards JE, Cheng YF, et al. Rumen microbial diversity in Svalbard reindeer, with particular emphasis on methanogenic archaea［J］. FEMS Microbiol Ecol, 2009, 70（3）：553-562.

［36］ Zhou M, Chung YH, Beauchemin KA, et al. Relationship between rumen methanogens and methane production in dairy cows fed diets supplemented with a feed enzyme additive［J］. J Appl Microbiol, 2011, 111（5）：1148-1158.

［37］ Cheng YF, Mao SY, Liu JX, et al. Molecular diversity analysis of rumen methanogenic Archaea from goat in eastern China by DGGE methods using different primer pairs［J］. Lett Appl Microbiol, 2009, 48（5）：585-592.

［38］ Franzolin R, St-Pierre B, Northwood K, et al. Analysis of rumen methanogen diversity in water buffaloes（Bubalus bubalis）under three different diets［J］. Microb Ecol, 2012, 64（1）：131-139.

［39］ Huang XD, Tan HY, Long R, et al. Comparison of methanogen diversity of yak（Bos grunniens）and cattle（Bos taurus）from

the Qinghai-Tibetan plateau, China［J］. BMC Microbiol, 2012, 12：237.

［40］ 裴彩霞, 毛勝勇, 朱偉云. 晉南牛瘤胃中古菌分子多樣性的研究［J］. 微生物學(xué)報(bào), 2008, 48（1）：8-14.

［41］ Tymensen LD, Beauchemin KA, McAllister TA. Structures of freeliving and protozoa-associated methanogen communities in the bovine rumen differ according to comparative analysis of 16S rRNA and mcrA genes［J］. Microbiology, 2012, 158（7）：1808-1817.［42］ Hook SE, Northwood KS, Wright AD, et al. Long-term monensin supplementation does not significantly affect the quantity or diversity of methanogens in the rumen of the lactating dairy cow［J］. Appl Environ Microbiol, 2009, 75（2）：374-380.

［43］ Frey JC, Pell AN, Berthiaume R, et al. Comparative studies of microbial populations in the rumen, duodenum, ileum and faeces of lactating dairy cows［J］. J Appl Microbiol, 2010, 108（6）：1982-1993.

［44］ Danielsson R, Schnürer A, Arthurson V, et al. Methanogenic population and CH4production in swedish dairy cows fed different levels of forage［J］. Appl Environ Microbiol, 2012, 78（17）：6172-6179.

［45］ Zened A, Combes S, Cauquil L, et al. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets［J］. FEMS Microbiol Ecol, 2013, 83（2）：504-514.

［46］ Thoetkiattikul H, Mhuantong W, Laothanachareon T, et al. Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing［J］. Curr Microbiol, 2013, 67（2）：130-137.

［47］ Fouts DE, Szpakowski S, Purushe J, et al. Next generation sequencing to define prokaryotic and fungal diversity in the bovine rumen［J］. PLoS One, 2012, 7（11）：48289.

［48］ Hristov AN, Callaway TR, Lee C, et al. Rumen bacterial, archaeal, and fungal diversity of dairy cows in response to ingestion of lauric or myristic acid［J］. J Anim Sci, 2012, 90（12）：4449-4457.

［49］ Lee HJ, Jung JY, Oh YK, et al. Comparative survey of rumen microbial communities and metabolites across one caprine and three bovine groups, using bar-coded pyrosequencing and1H nuclear magnetic resonance spectroscopy［J］. Appl Environ Microbiol, 2012, 78（17）：5983-5993.

［50］ Tymensen L, Barkley C, McAllister TA. Relative diversity and community structure analysis of rumen protozoa according to T-RFLP and microscopic methods［J］. J Microbiol Methods, 2012, 88（1）：1-6.

［51］ Hook SE, Steele MA, Northwood KS, et al. Impact of high-concentrate feeding and low ruminal pH on methanogens and protozoa in the rumen of dairy cows［J］. Microb Ecol, 2011, 62（1）：94-105.

［52］ Yá?ez-Ruiz DR, Macías B, Pinloche E, et al. The persistence of bacterial and methanogenic archaeal communities residing in the rumen of young lambs［J］. FEMS Microbiol Ecol, 2010, 72（2）：272-278.

［53］ Leahy SC, Kelly WJ, Altermann E, et al. The genome sequence of the rumen methanogenMethanobrevibacter ruminantiumreveals new possibilities for controlling ruminant methane emissions［J］. PLoS One, 2010, 5（1）：8926.

［54］ Lee JH, Rhee MS, Kumar S, et al. Genome sequence ofMethanobrevibactersp. strain jh1, isolated from rumen of Korean native cattle［J］. Genome Announc, 2013, 1（1）：2-13.

［55］ Wedlock DN, Janssen PH, Leahy SC, et al. Progress in the development of vaccines against rumen methanogens［J］. Animal, 2013, 7（2）：244-252.

［56］ Morgavi DP, Kelly WJ, Janssen PH, et al. Rumen microbial（meta）genomics and its application to ruminant production［J］. Animal, 2013, 7（1）：184-201.

（責(zé)任編輯 狄艷紅）

The Application of Uncultured Methods in the Study of Ruminal Methanogen Population

Wang Xingwen1，2Wang Jiaqi1Zhao Shengguo1Li Fadi2Bu Dengpan1
（1. State Key Laboratory of Animal Science，Institute of Animal Science，Chinese Academy of Agricultural Science，Beijing 100193；2. College of Animal Science and Technology，Gansu Agricultural University，Lanzhou 730070）

The ruminal methanogen are strictly anaerobic, which are hard to culture.The ruminal methanogen which have been cultured successfully are only 4 species.while hundreds of species of methanogen could not be cultured.The application and development of uncultured technology promotes the research of ruminal methanogen, which overcome the limitation of the tranditional technology of isolation and culturing. Several dominant culture-independent technologies used to study ruminal methanogens were introduced in this paper to provide

on methods for future ruminal methanogen study.

Rumen Methanogen Uncultured methods Whole genome sequencing High throughput DNA sequencing

2013-09-03

國(guó)家自然科學(xué)基金項(xiàng)目（31261140365），科技支撐計(jì)劃項(xiàng)目（2012BAD12B02），中國(guó)科學(xué)院知識(shí)創(chuàng)新工程重要方向項(xiàng)目（XDA01010301）

王興文，男，碩士研究生，研究方向：反芻動(dòng)物營(yíng)養(yǎng)；E-mail：xing.wen111@163.com

卜登攀，副研究員，碩士生導(dǎo)師，研究方向：反芻動(dòng)物營(yíng)養(yǎng)與代謝調(diào)控；E-mail：burdengpan@gmail.com