楊承劍，韋升菊，梁辛，梁賢威*，鄒彩霞

(1.中國(guó)農(nóng)業(yè)科學(xué)院廣西水牛研究所，廣西 南寧 530001；2.農(nóng)業(yè)部(廣西)水牛遺傳繁育重點(diǎn)實(shí)驗(yàn)室，廣西 南寧 530001)

德昌水牛是中國(guó)地方水牛優(yōu)良品種之一，屬于沼澤型水牛，其體格較大，體質(zhì)結(jié)實(shí)，生長(zhǎng)發(fā)育快，適應(yīng)性廣，役力強(qiáng)，主要分布于海拔為480～2 520 m的四川省涼山彝族自治州安寧河流域的德昌、西昌、冕寧、會(huì)理等縣和渡口市的米易縣[1]。深入研究德昌水牛的消化生理以及營(yíng)養(yǎng)特性有助于地方水牛品種的保護(hù)以及開發(fā)利用。

瘤胃是反芻動(dòng)物獨(dú)特的消化器官，其中棲息有數(shù)以億計(jì)的微生物，包括細(xì)菌、產(chǎn)甲烷古菌、原蟲和真菌等。這些微生物之間的互作對(duì)于飼料的消化利用具有重要作用，可以脂肪酸和微生物蛋白的形式為宿主提供能量。近年來，對(duì)于產(chǎn)甲烷古菌的研究已成為研究者們關(guān)注的重點(diǎn)。飼料被消化過程中所產(chǎn)生的中間代謝產(chǎn)物，如氫、甲酸、乙酸、甲醇等，可通過產(chǎn)甲烷菌還原二氧化碳進(jìn)而生成甲烷[2]。反芻動(dòng)物的甲烷排放意味著能量的損失以及溫室氣體的增加[3]。目前，已有研究[4–5]表明，多種動(dòng)物體內(nèi)均可分離出產(chǎn)甲烷菌。部分研究也利用非培養(yǎng)方式，如16S rRNA 基因克隆文庫(kù)分析、real–time PCR 等方法研究瘤胃內(nèi)產(chǎn)甲烷菌的多樣性及數(shù)量[6–8]。不同肥胖程度的人類或者不同品種的豬腸道中的產(chǎn)甲烷菌結(jié)構(gòu)也存在差異[9–10]。這些研究表明，胃腸道中產(chǎn)甲烷菌的結(jié)構(gòu)組成可能與物種的特異性或與特定功能緊密相關(guān)[11]。作為四川特色遺傳資源品種之一，德昌水牛瘤胃微生物可能具有其獨(dú)特的微生物區(qū)系，但目前關(guān)于德昌水牛瘤胃微生物的組成，尤其是產(chǎn)甲烷菌多樣性方面的研究仍處于空白。筆者利用16S rRNA 基因克隆文庫(kù)技術(shù)分析德昌水牛瘤胃中產(chǎn)甲烷菌的結(jié)構(gòu)組成及多樣性，旨在為進(jìn)一步探索德昌水牛瘤胃微生物生態(tài)系統(tǒng)功能以及品種資源保護(hù)提供參考依據(jù)。

1 材料與方法

1.1 供試水牛及飼料

選用3 頭體況基本一致的5 歲左右的健康成年雌性德昌水牛，飼養(yǎng)于廣西壯族自治區(qū)水牛研究所水牛場(chǎng)(國(guó)家級(jí)重點(diǎn)種畜場(chǎng))。供試水牛每天08:00和16:00 定時(shí)飼喂，自由飲水。自由采食粗料。補(bǔ)充適量精料(每天每頭5.0 kg)。

粗料由玉米青貯料、木薯渣、啤酒渣、新鮮象草組成。精料配方：玉米53.5 %、豆粕6.0%、麥麩15.0%、棉粕17.0%、磷酸氫鈣1.5%、貝殼粉1.0%、食鹽2.0%、小蘇打3.0%、預(yù)混料1.0%。精料營(yíng)養(yǎng)水平：NE 6.83 MJ/Kg，CP 16.75/%，Ca 0.80/%，P 0.94/%，CF 4.54/%。預(yù)混料組成：V–E 3 000 IU/kg，V–D 150 kIU/kg，V–A 500 kIU/kg，Cu 1.3 g/kg，F(xiàn)e 4.0 g/kg，Mn 3.0 g/kg，I 80 mg/kg，Zn 6.0 g/kg，Co 80 mg/kg，Se 50 mg/kg。

1.2 方 法

試驗(yàn)期30 d，第30 天晨飼前通過胃管式瘤胃液采樣器采集瘤胃內(nèi)容物250 mL，裝入干凈充滿二氧化碳厭氧瓶中，置于碎冰塊中冷藏，迅速帶回實(shí)驗(yàn)室，保存于–80 ℃冰箱，備用。

1.3 總DNA 的提取

用寬口吸頭從充分混勻的瘤胃內(nèi)容物樣品(約250 mL)中吸取1.5 mL(包括固液兩相)，12 000 g 離心5 min，去除上清液，收集沉淀，依據(jù)文獻(xiàn)[12]所描述的機(jī)械破壁及試劑盒提取相結(jié)合的方法從沉淀樣品中提取總DNA。

1.4 PCR 擴(kuò)增及16S rRNA 克隆文庫(kù)構(gòu)建

利用產(chǎn)甲烷菌特異性引物Met86F/Met1340R[6]擴(kuò)增16S rRNA 基因序列。引物序列如下： Met86F，5′–GCTCAGTAACACGTGG–3′；Met1340R，5′–CG GTGTGTGCAAGGAG–3′。PCR 反應(yīng)體系如下：10×Buffer 5μL；dNTP(10 mmoL/L) 1 μL；MgCL2(25 mmoL/L) 4 μL；TaqDNA聚合酶(5 U/μL) 0.8 μL；l μL總DNA(l～10 ng/mL)，Met86F 和Met1340R 引物(10 μmoL/L)各l μL；加無(wú)菌雙蒸水至50 μL。反應(yīng)條件如下：94 ℃預(yù)變性3 min；94 ℃變性30 s，58 ℃退火30 s，72 ℃延伸90 s，40 個(gè)循環(huán)；72 ℃終延伸10 min。利用瓊脂糖凝膠純化回收試劑盒(美國(guó)Promega)純化PCR 產(chǎn)物。將來自3 頭牛的等體積PCR 產(chǎn)物混合后，按照廠家說明書連接到載體pMD19–T(日本Takara)上進(jìn)行克隆，經(jīng)藍(lán)白斑篩選，隨機(jī)挑取白色克隆，用pMD19–T 通用引物M13–47 (5′–CGCCAGGGTTTTCCCAGTCACGAC–3′)和RV –M(5′–GAGCGGATAACAATTTCACACAGG–3′)對(duì)陽(yáng)性克隆進(jìn)行驗(yàn)證，將PCR 檢測(cè)為陽(yáng)性的克隆送上海生工生物技術(shù)有限公司測(cè)序。

1.5 16S rRNA 基因序列分析

對(duì)測(cè)序獲得的序列用Decipher軟件去除嵌合體序列[13]。使用Ribosomal Database Project (RDP) Release 11 中的Classifier 程序?qū)λ行蛄羞M(jìn)行分類[14]。用Mothur 軟件劃分分類操作單元(OTU)，將相似性大于97%的序列視為同一個(gè)OTU[15]。選取每個(gè)OTU代表序列，利用Blast 程序在Genbank 中搜索相似性最高序列。用Clustal X Version 1.83 軟件進(jìn)行多重比對(duì)，通過MEGA 4.0(Molecular evolutionary genetics analysis，MEGA)軟件計(jì)算出序列的系統(tǒng)進(jìn)化距離，采用Neighbor–Joining 方法構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹，1 000 次隨機(jī)抽樣，計(jì)算自舉值(Bootstrap)以評(píng)估系統(tǒng)發(fā)育樹的置信度。本研究所獲得的序列在DDBJ數(shù)據(jù)庫(kù)里的登錄號(hào)為AB905615～ AB905718。

2 結(jié)果與分析

2.1 16S rRNA 基因克隆文庫(kù)分析

從建立的德昌水牛瘤胃產(chǎn)甲烷菌16S rRNA 基因克隆文庫(kù)中，隨機(jī)挑選107 個(gè)陽(yáng)性克隆進(jìn)行測(cè)序。對(duì)測(cè)得的序列進(jìn)行處理，剔除嵌合序列后獲得99個(gè)序列，RDP 分析結(jié)果表明，99 個(gè)序列均屬于古細(xì)菌域(Archaea)序列，其中93 個(gè)序列(占總序列數(shù)的94.2%)歸屬M(fèi)ethanobrevibacter 的16S rRNA 序列、5 個(gè)(占總序列數(shù)的4.8%)歸屬M(fèi)ethanosphaera (甲烷球菌屬)的16s rRNA 序列、還有1 個(gè)(占總序列數(shù)的1.0%)歸屬unclassified_ Methanobacteriaceae (未分類甲烷桿菌)的16S rRNA 序列。以97%序列相似性為閾值，用Mothur 軟件對(duì)序列進(jìn)行分析，99條序列可分為19 個(gè)OTUs(表1)。

表1 德昌水牛瘤胃產(chǎn)甲烷古菌16S rRNA 基因序列分析Table 1 Sequence analysis of archaeal 16S rRNA gene from the rumen of Dechang buffaloes

將99 條序列與Genbank 中收錄的序列進(jìn)行Blast 分析，結(jié)果表明，99 個(gè)序列中，有96 個(gè)序列(17 個(gè)OTUs)占總序列的97.0%，與已知細(xì)菌16S rRNA 序列相似性≥97%；其中，3 個(gè)OTUs(約占總序列的57.6%)與Methanobrevibacter gottschalkii 的16S rRNA 相似；10 個(gè)OTUs(約占總序列的32.3%)與Methanobrevibacter millerae 的16S rRNA 相似；2 個(gè)OTUs(約占總序列的4.0%)與Methanosphaera stadtmanae 的16S rRNA 相似；1 個(gè)OTUs(約占總序列的1.0%)與Methanobrevibacter ruminantium 的16S rRNA 相似；1 個(gè)OTUs(約占總序列的1.0%)與Methanobrevibacter boviskoreani 的16S rRNA 相似；1 個(gè)OTUs(約占總序列的1.0%)與Methanobrevibacter smithii 的16S rRNA 相似。另外3 個(gè)序列(2 個(gè)OTUs)與已培養(yǎng)菌的16S rRNA 序列的相似性為90%～97%，其中1 個(gè)OTUs(約占總序列的2.0%)與Methanosphaera stadtmanae 的16S rRNA 相似；1 個(gè)OTUs(約占總序列的1.0%)與Methanobrevibacter olleyae 的16S rRNA 相似。

2.2 16S rRNA 基因系統(tǒng)進(jìn)化分析

以Thermosphaera aggregans(X99556)和Pyrolobus fumartus (X99555)作為外群(Outgroup)，將19 個(gè)OTUs 代表序列與其他最相似已知序列構(gòu)建系統(tǒng)發(fā)育樹。結(jié)果(圖1)表明，D51、D66 與Methanosphaera stadtmanae 聚集為一簇，屬于Methanosphaera；D11、D64、D24、D4、19、D62、32、D6、6、30 與Methanobrevibacter smithii、Methanobrevibacter gottschalkii、Methanobrevibacter millerae 聚集為一簇(SGMT 簇)，屬于Methanobrevibacter；D49 與Methanobrevibacter ruminantium、Methanobrevi- bacter olleyae 聚集為一簇(RO 簇)，但與Methano- brevibacter 中其他序列相隔較遠(yuǎn)；SGMT 簇序列和 RO 簇序列所占總序列的比例分別為75.8%、1.0%；D17、D16、3、35、7、D59 單獨(dú)聚集為一簇，與任何已知相近序列都相隔較遠(yuǎn)，它們可能代表Methanobrevibacter 中新的種。

圖1 德昌水牛瘤胃產(chǎn)甲烷古菌16S rRNA 基因序列系統(tǒng)發(fā)育分析Fig.1 Phylogenetic analysis of methanogens 16S rRNA gene sequences from the rumen of Dechang water buffaloes

3 討 論

在瘤胃厭氧條件下，多種厭氧微生物通過共同作用完成有機(jī)物的降解。被動(dòng)物消化的飼料成分，包括植物性結(jié)構(gòu)碳水化合物、蛋白質(zhì)和其他有機(jī)聚合物，通過初次厭氧發(fā)酵后都被降解成簡(jiǎn)單的物質(zhì)。這些簡(jiǎn)單的化合物通過具有水解復(fù)雜化合物功能的微生物(初次發(fā)酵微生物)和其他不具有水解復(fù)雜化合物功能的微生物(次級(jí)發(fā)酵微生物)的共同作用，被進(jìn)一步轉(zhuǎn)變成脂肪酸、二氧化碳和氫氣等。在瘤胃內(nèi)利用二氧化碳作為碳源，氫作為主要的電子供體的氫利用途徑是主要的產(chǎn)甲烷途徑[16]。甲烷的生成可避免氫分壓升高，以免抑制電子傳遞反應(yīng)中的微生物酶的正常功能，尤其是NADH 脫氫酶，它若被抑制可導(dǎo)致NADH 的積累，最終降低瘤胃發(fā)酵效率。深入了解水牛瘤胃產(chǎn)甲烷菌的結(jié)構(gòu)組成及其多樣性，有助于調(diào)控水牛瘤胃甲烷生成，減少甲烷排放量。產(chǎn)甲烷菌屬于古菌，基于已有的數(shù)據(jù)分析，許多被檢測(cè)到的瘤胃古菌都屬于產(chǎn)甲烷微生物，其中主要是Methanobrevibacter，占瘤胃古菌近2/3(61.6%)。剩余的1/3 的古細(xì)菌的系統(tǒng)發(fā)育型大致均等，屬于Methanomicrobium 的占14.9%，被標(biāo)記為不可培養(yǎng)的瘤胃C 簇或RCC 的占15.8%[17]。目前，已有許多關(guān)于反芻動(dòng)物，如綿羊[18–19]、黃牛[7,20]以及水牛[21–23]的瘤胃微生物菌群結(jié)構(gòu)組成的報(bào)道，其中水牛以印度水牛為主。關(guān)于中國(guó)地方品種水牛瘤胃微生物菌群結(jié)構(gòu)組成，尤其是在產(chǎn)甲烷菌多樣性方面的研究仍十分薄弱。據(jù)筆者所知，本研究是國(guó)內(nèi)外第一個(gè)報(bào)道中國(guó)地方品種德昌水牛瘤胃產(chǎn)甲烷菌多樣性的。

瘤胃微生物菌群結(jié)構(gòu)組成受地理位置、周邊環(huán)境、飼料營(yíng)養(yǎng)水平以及飼喂頻率、添加劑、動(dòng)物種類以及動(dòng)物遺傳背景等的影響。一般認(rèn)為，瘤胃中主要的甲烷菌為瘤胃甲烷短桿菌(Methanobrevibacter ruminantium)和巴氏甲烷八疊球菌(Methanosarcina barkerii)。本研究通過16S rRNA 基因克隆文庫(kù)分析，表明德昌水牛瘤胃內(nèi)產(chǎn)甲烷菌以Methanobrevi- bacter 為主，這與國(guó)內(nèi)外(委內(nèi)瑞拉、澳大利亞、中國(guó))許多在其他反芻動(dòng)物(綿羊、袋鼠、牦牛)上的研究結(jié)果相一致[14,18,24–26]。但是，Chaudhary 等[21]的研究表明，印度摩拉水牛瘤胃產(chǎn)甲烷菌以Methano- microbium 為主；Singh 等[27]的研究表明，蘇替水牛(Surti)瘤胃內(nèi)產(chǎn)甲烷菌以Methanobacteriales 為主；Hook 等[28]的研究表明，在澳大利亞綿羊瘤胃內(nèi)則以未知古菌為主；Shin 等[29]以及Tajima 等[30]的研究發(fā)現(xiàn)，奶牛瘤胃古菌以Methanomicrobium 為主(分別為85.6%和60.9%)，但Methanobrevibacter 數(shù)量很少。也有研究[31]表明，奶牛體內(nèi)存在大量的Methanosphaera(26.8%)和Methanimicrococcus(14.6%)，但是Methanobrevibacter 占優(yōu)勢(shì)(58.5%)。說明不同動(dòng)物種類之間的瘤胃產(chǎn)甲烷菌結(jié)構(gòu)組成也存在差異。這些異同與研究方法、動(dòng)物種類、日糧、動(dòng)物管理等之間的相互關(guān)系還需要進(jìn)一步研究。

反芻動(dòng)物瘤胃內(nèi)仍存在許多未知的還未培養(yǎng)出來的產(chǎn)甲烷菌[32]。本研究的系統(tǒng)進(jìn)化樹分析結(jié)果也表明，德昌水牛瘤胃內(nèi)存在許多未知的Methano- brevibacter 產(chǎn)甲烷菌。Paul 等[33]的研究表明，Methanobrevibacter 能夠利用二氧化碳和氫或者甲酸生成甲烷；Rea 等[34]從黃牛瘤胃中分離到的Methanobrevibacter millerae 能夠利用甲酸生長(zhǎng)；Samuel 等[35]發(fā)現(xiàn)Methanobrevibacter smithii PS 可利用二氧化碳和氫或者甲酸生成甲烷；Denman 等[7]的研究表明，反芻動(dòng)物體內(nèi)的Methanobrevibacter smithii 能夠影響日糧中的多糖成分的利用效率。以上研究表明，Methanobrevibacter 產(chǎn)甲烷菌對(duì)于反芻動(dòng)物瘤胃甲烷生成以及其他生理功能具有重要作用。關(guān)于Methanobrevibacter 產(chǎn)甲烷菌在德昌水牛瘤胃生理功能的作用還有待進(jìn)一步研究。

St–pierre 等[36]提議將Methanobrevibacter 中的產(chǎn)甲烷菌序列分為兩大類：與Methanobrevibacter smithii(S) 、 Methanobrevibacter gottschalkii(G) 、Methanobrevibacter millerae (M)和Methanobrevibacter thaurei(T)聚集在一塊的統(tǒng)稱為SGMT 簇；與M. ruminantium (R) 及Methanobrevibacter olleyae (O)聚集在一起的稱為RO 簇。不同種類動(dòng)物或不同條件下的動(dòng)物胃腸道的SGMT 和RO 簇序列分布規(guī)律可能存在一定差異。在同等日糧條件下，荷斯坦奶牛瘤胃中的SGMT 簇序列的比例要高于娟姍奶牛[20]。在本研究中SGMT 簇序列和RO 簇序列所占總序列的比例分別為75.8%、1.0%，即在德昌水牛瘤胃中SGMT 簇序列占大多數(shù)。不同動(dòng)物種類及飼養(yǎng)條件造成的SGMT 簇和RO 簇序列分布的差異，是否可以用于估測(cè)瘤胃甲烷生成，進(jìn)而用于調(diào)節(jié)瘤胃甲烷生成，還有待深入研究。

[1]《四川家畜家禽品種志》編輯委員會(huì)．四川家畜家禽品種志[M]．成都：四川科學(xué)技術(shù)出版社，1987：36–38．

[2]Samuel B S，Gordon J I．A humanized gnotobiotic mouse model of host-archaeal-bacterial mutualism[J]．Proc Natl Acad Sci U S A，2006，103(26)：10011–10016．

[3]Johnson K A，Johnson D E．Methane emissions from cattle[J]．Journal of Animal Science，1995，73(8)：2483–2492．

[4]Jarvis G N，Strompl C，Burgess D M，et al．Isolation and identification of ruminal methanogens from grazing cattle[J]．Curr Microbiol，2000，40(5)：327–332．

[5]Miller T L W M，Kusel E．Isolation and characteri- zation of methanogens from animal feces[J]． Syst Appl Microbiol，1986，8(3)：234–238．

[6]Wright A D G，Williams A J，Winder B，et al．Molecular diversity of rumen methanogens from sheep in western Australia[J]．Appl Environ Microbiol，2004，70(3)：1263–1270．

[7]Denman S E，Tomkins N W，Mcsweeney C S. Quantitation and diversity analysis of ruminal metha- nogenic populations in response to the antimethanogenic compound bromochloromethane[J]．FEMS Microbiol Ecol,2007，62(3)：313–322．

[8]Pei C X，Mao S Y，Cheng Y F，et al．Diversity，abundance and novel 16S rRNA gene sequences of methanogens in rumen liquid，solid and epithelium fractions of Jinnan cattle[J]．Animal，2010，4(1)：20–29．

[9]DiBaise J K，Zhang H，Crowell，M D，et al．Gut microbiota and its possible relationship with obesity[J]. MAYO Clinic Proceedings，2008，83(4)：460–469．

[10]Luo Y H，Su Y，Wright A D，et al．Lean breed Landrace pigs harbor fecal methanogens at higher diversity and density than obese breed Erhualian pigs[J]．Archaea，2012：http://dx.doi.org/10.1155/2012/605289．

[11]Shi P J，Meng K，Zhou Z G，et al．The host species affects the microbial community in the goat rumen[J]．Letters in Applied Microbiology，2008，46(1)：132–135．

[12]Denman S E，Mcsweeney C S．Development of a real- time PCR assay for monitoring anaerobic fungal and cellulolytic bacterial populations within the rumen[J]. FEMS Microbiol Ecol，2006，58(3)：572–582．

[13]Wright E S，Yilmaz L S，Noguera D R．DECIPHER，a search-based approach to chimera identification for 16S rRNA sequences[J]．Appl Environ Microbiol，2012，78(3)：717–725．

[14]Cole J R，Wang Q，Cardenas E，et al．The ribosomal database project：Improved alignments and new tools for rRNA analysis[J]．Nucleic Acids Res，2009，37(Sup.1)：D141–145．

[15]Schloss P D，Westcott S L，Ryabin T，et al．Introducing mothur：Open-source，platform-independent，community- supported software for describing and comparing microbial communities[J]．Appl Environ Microbiol，2009，75(23)：7537–7541．

[16]Hungate R E．Hydrogen as an intermediate in rumen fermentation[J]．Archives of Microbiology，1967，59 (1/3)：158–164．

[17]Janssen P H，Kris M．Structure of the archaeal community of the rumen[J]．Appl Environ Microbiol，2008，74(12)：3619–3625．

[18]Wright A D，Ma X，Obispo N E．Methanobrevibacter phylotypes are the dominant methanogens in sheep from Venezuela[J]．Microb Ecol，2008，56(2)：390–394．

[19]Wright A D，Toovey A F，Pimm C L．Molecular identification of methanogenic archaea from sheep in Queensland，Australia reveal more uncultured novel archaea[J]．Anaerobe，2006，12(3)：134–139．

[20]King E E，Smith R P，St–pierre B，et al．Differences in the rumen methanogen populations of lactating Jersey and Holstein dairy cows under the same diet regimen[J]. Appl Environ Microbiol，2011，77(16)：5682–5687．

[21]Chaudhary P P，Sirohi S K．Dominance of Methano- microbium phylotype in methanogen population present in Murrah buffaloes (Bubalus bubalis)[J]．Lett Appl Microbiol，2009，49(2)：274–277．

[22]Tan H Y，Sieo C C，Lee C M，et al．Diversity of bovine rumen methanogens in vitro in the presence of condensed tannins，as determined by sequence analysis of 16S rRNA gene library[J]．J Microbiol，2011，49(3)：492–498．

[23]Tatsuoka N，Mohammed N，Mitsumori M，et al. Phylogenetic analysis of methyl coenzyme-M reductase detected from the bovine rumen[J]．Lett Appl Microbiol，2004，39(3)：257–260．

[24]Hook S E，Wright A D，Mcbride B W．Methanogens：Methane producers of the rumen and mitigation strategies[J]．Archaea，2010：http://dx.doi.org/10.1155/ 2010/945785．

[25]Evans P N，Hinds L A，Sly L I，et al．Community composition and density of methanogens in the foregut of the Tammar wallaby (Macropus eugenii)[J]．Appl Environ Microbiol，2009，75(8)：2598–2602．

[26]An D，Dong X，Dong Z．Prokaryote diversity in the rumen of yak (Bos grunniens) and Jinnan cattle (Bos taurus) estimated by 16S rDNA homology analyses[J]. Anaerobe，2005，11(4)：207–215．

[27]Singh K M，Pandya P R，Parnerkar S，et al．Methanogenic diversity studies within the rumen of Surti buffaloes based on methyl coenzyme M reductase A (mcrA) genes point to methanobacteriales[J]．Polish Journal of Micro- biology，2010，59(3)：175–178．

[28]Hook S E，Northwood K S，Wright A D G，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，374–380．

[29]Shin E C，Choi B R，Lim W J，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．

[30]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．

[31]Whitford M F，Teather R M，F(xiàn)orster R J．Phylogenetic analysis of methanogens from the bovine rumen[J]．BMC Microbiol，2001：http://www.biomedcentral.com/1471- 2180/1/5．

[32]Tymensen L D，Mcallister T A．Community structure analysis of methanogens associated with rumen protozoa reveals bias in universal archaeal primers[J]．Appl Environ Microbiol，2012，78(11)：4051–4056．

[33]Paul K，Nonoh J O，MikulskiI L，et al. “Methanoplas- matales,” thermoplasmatales-related archaea in termite guts and other environments，are the seventh order of methanogens[J]．Appl Environ Microbiol，2012，78(23)：8245–8253．

[34]Rea S，Bowman J P，Popovski S，et al. Methano- brevibacter millerae sp．nov．a(chǎn)nd Methanobrevibacter olleyae sp．nov．，methanogens from the ovine and bovine rumen that can utilize formate for growth[J]．Int J Syst Evol Microbiol，2007，57(Pt 3)：450–456．

[35]Samuel B S，Hansen E E，Manchester J K，et al．Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut[J]．Proc Natl Acad Sci U S A，2007，104(25)：10643–10648．

[36]St–pierre B，Wright A D．Molecular analysis of methano- genic archaea in the forestomach of thealpaca (Vicugna pacos)[J]．BMC Microbiol，2012，12(1)：doi：10.1186/ 1471-2180-12-1．