肖英平 王軍軍 李天天 唐 標(biāo) 何祥祥 李開鋒 楊 華*

(1.浙江農(nóng)業(yè)科學(xué)院農(nóng)產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所，杭州310021；2.北京食品營養(yǎng)與人類健康高精尖創(chuàng)新中心，北京100193；3. 中國農(nóng)業(yè)大學(xué)動物科技學(xué)院，北京100193)

金華豬和長白豬糞便微生物移植小鼠的腸道古菌結(jié)構(gòu)分析

肖英平1王軍軍2,3李天天2,3唐 標(biāo)1何祥祥1李開鋒1楊 華1*

(1.浙江農(nóng)業(yè)科學(xué)院農(nóng)產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所，杭州310021；2.北京食品營養(yǎng)與人類健康高精尖創(chuàng)新中心，北京100193；3. 中國農(nóng)業(yè)大學(xué)動物科技學(xué)院，北京100193)

本試驗旨在比較肥胖型的金華豬和瘦肉型的長白豬腸道古菌結(jié)構(gòu)的差異，分析菌群移植小鼠腸道古菌結(jié)構(gòu)。采集金華豬和長白豬的新鮮糞便，通過灌胃移植到經(jīng)廣譜抗生素處理過的小鼠腸道中，提取2種豬糞便及其移植小鼠的盲腸內(nèi)容物總基因組DNA，使用古菌通用引物對古菌16S rRNA的V4區(qū)進(jìn)行PCR擴(kuò)增，通過Illumina Hiseq測序平臺對擴(kuò)增產(chǎn)物進(jìn)行高通量測序，使用QIIME等軟件對測序序列進(jìn)行分析統(tǒng)計。結(jié)果表明：金華豬和長白豬糞便及其移植小鼠盲腸內(nèi)容物中的古菌主要包含3個門，分別為奇古菌門(Thaumarchaeota)、廣古菌門(Euryarchaeota)和泉古菌門(Crenarchaeota)，其中金華豬糞便中廣古菌門的豐度高于長白豬；在屬水平上主要包含4個已知屬，優(yōu)勢屬為餐古菌屬(Cenarchaeum)。二者移植小鼠的盲腸內(nèi)容物古菌結(jié)構(gòu)，在門水平上，移植小鼠的奇古菌門比例占70%以上；奇古菌門的古菌幾乎全部來源于餐古菌屬；在門和屬的水平上，金華豬和長白豬糞便移植小鼠的盲腸古菌結(jié)構(gòu)無顯著差異(P>0.05)，但基于操作分類單元(OTU)水平的主成分分析顯示出很好的組內(nèi)聚集。由此可見，奇古菌門的餐古菌屬是金華豬和長白豬糞便和其移植小鼠盲腸中的優(yōu)勢古菌；金華豬糞便中廣古菌門的豐度高于長白豬糞便，且二者移植小鼠盲腸內(nèi)容物之間古菌結(jié)構(gòu)在OTU水平上存在一定的差異。

古菌；金華豬；長白豬；糞菌移植；高通量測序

豬肉歷來是我國居民肉類消費品的主要來源，在國計民生中占有重要地位。金華豬是我國著名的地方品種，有“中國熊貓豬”的美譽(yù)，其頭部和尾部為黑色，故又名兩頭烏。具有皮薄骨細(xì)、肉質(zhì)鮮美、性成熟早、繁殖能力強(qiáng)、雜交優(yōu)勢顯著的特點，但其生長速度慢、瘦肉率低[1]。而長白豬(Landrace pig)是優(yōu)良的外來品種，原產(chǎn)于丹麥，體軀瘦長，皮膚、被毛呈白色，生長速度快，瘦肉率高，但存在體質(zhì)較弱、抗逆性較差、對飼養(yǎng)條件要求較高、繁殖能力低下等缺點[2]。

腸道微生物區(qū)系是哺乳動物體內(nèi)最復(fù)雜和最大的微生態(tài)系統(tǒng)，對宿主的健康、進(jìn)化及行為產(chǎn)生明顯的影響[3-4]。古菌(Archaea)是腸道微生物的重要組成部分之一。古菌最初發(fā)現(xiàn)于高熱、高鹽度、缺氧的極端環(huán)境中(如火山口、鹽湖)，后通過16S rRNA的比較，發(fā)現(xiàn)其明顯不同于細(xì)菌(Bacteria)和真核生物(Eukarya)，與后兩者并列稱為三域系統(tǒng)[5-6]。古菌廣泛分布于各種自然環(huán)境中，如海洋、熱泉、鹽堿湖、土壤、沼澤以及動物的消化道。目前廣古菌門(Euryarchaeota)和泉古菌門(Crenarchaeota)是研究較多的古菌，尤其是前者，包含了一大類產(chǎn)甲烷菌，廣泛存在于各類動物的消化道，是動物營養(yǎng)研究的熱點之一[7-9]。腸道古菌也有助于維持消化道中良好的微生態(tài)系統(tǒng)。而奇古菌門(Thaumarchaeota)是近年來發(fā)現(xiàn)的一個新的古菌類群，廣泛存在于中溫環(huán)境中，可能在自然界重要元素(氮、碳)的生物地球化學(xué)循環(huán)中發(fā)揮著重要作用[10-12]，但其在動物消化道中的分布及功能報道很少[13-14]。

近年來，糞菌移植(fecal microbiota transplantation，F(xiàn)MT)作為一種歷史悠久且可以重建腸道菌群的療法重新被臨床所關(guān)注，為治療腸道內(nèi)及腸道外疾病提供幫助[4,15]。與此同時，作為一種功能獲得性(gain-of-function)試驗手段，F(xiàn)MT也越來越多被應(yīng)用于動物的微生物結(jié)構(gòu)功能研究[16]，小鼠是其中最常用的動物模型之一[17-18]。研究發(fā)現(xiàn)，將豬和斑馬魚的腸道微生物移植到小鼠腸道中，移植供體的菌群結(jié)構(gòu)在一定程度被宿主影響和重塑[19-21]。但目前在菌群移植過程中，普遍關(guān)注的是細(xì)菌結(jié)構(gòu)，而對于古菌在移植狀態(tài)下的情況研究較少。本試驗通過糞菌移植，分別構(gòu)建了模擬金華豬和長白豬2種豬腸道微生物環(huán)境的小鼠模型，采用高通量測序技術(shù)分析了2種豬糞便及其移植小鼠盲腸內(nèi)容物的古菌結(jié)構(gòu)，為研究跨物種糞便微生物群移植過程中古菌結(jié)構(gòu)的變化提供了數(shù)據(jù)資料。

1 材料與方法

1.1 試驗動物

動物試驗在北京斯貝福實驗動物科技公司進(jìn)行。選取28日齡斷奶的無特定病原菌(SPF)C57BL/6J小鼠24只，雌雄各12只，飼養(yǎng)在無菌隔離包中，飼喂無菌鼠糧。室內(nèi)溫度保持恒溫23 ℃，相對濕度為63%，晝夜光照交替時間為12 h/12 h。

1.2 小鼠抗生素處理

試驗采取雌雄分養(yǎng)，每個籠子6只小鼠，參照Wang等[22]的方法，采用廣譜抗生素處理小鼠以構(gòu)建無菌小鼠模型。具體方法為：在小鼠飲用水中加入0.5 g/L萬古霉素、1 g/L硫酸新霉素、1 g/L甲硝唑和1 g/L氨芐西林，自由口服28 d。

1.3 豬糞便處理

分別采集1頭產(chǎn)2胎的經(jīng)產(chǎn)金華豬母豬和長白豬母豬新鮮糞便，參照Pang等[23]的方法，制備移植供體的腸道微生物。具體方法為：將新鮮糞便迅速懸浮于預(yù)還原處理過的無菌生理鹽水中(每克糞便加入4 mL生理鹽水，即稀釋5倍)，漩渦混勻1 min，靜置 1 min，快速吸取上清，盛于1.5 mL滅菌處理過的離心管中于-20 ℃保存用于菌群移植；同時另取一部分上清用于高通量測序分析古菌結(jié)構(gòu)。

1.4 試驗鼠接種和取樣

取上述抗生素處理的小鼠，分成2組，每組雌雄各6只，采用灌胃針分別用金華豬和長白豬糞便菌懸液進(jìn)行灌胃，每次0.2 mL，連續(xù)灌注7 d。接種后連續(xù)飼養(yǎng)28 d(即腸上皮及其上覆黏液層經(jīng)過幾輪更新后)，屠宰小鼠取盲腸內(nèi)容物用于高通量測序分析古菌結(jié)構(gòu)。選取盲腸是由于其管腔內(nèi)容物可以簡便且可靠地獲得，且其中的微生物種群豐富[20]。

1.5 DNA提取和古菌16S rRNA基因擴(kuò)增

用QIAamp DNA Stool Mini Kit(QIAGEN公司)試劑盒抽提金華豬和長白豬糞便菌懸液及其灌胃小鼠盲腸內(nèi)容物總基因組DNA，采用古菌通用引物519F(5′-CAGCCGCCGCGGTAA-3′)和Arch915R(5′-GTGCTCCCCCGCCAATTCCT-3′)對古菌16S rRNA基因V4區(qū)進(jìn)行擴(kuò)增。

1.6 高通量測序分析

測序由北京諾禾致源生物信息科技有限公司完成。采用Illumina Hiseq高通量測序平臺對金華豬和長白豬糞便菌懸液及其灌胃小鼠盲腸內(nèi)容物古菌16S rRNA基因V4區(qū)進(jìn)行測序。

1.7 數(shù)據(jù)處理與分析

根據(jù)條碼序列(barcode)將雙端測序序列(read)與樣品進(jìn)行匹配，而后切除條碼序列及引物序列。使用Flash 1.2.7軟件(http://ccb.jhu.edu/software/FLASH/)將雙端序列進(jìn)行融合(merge)，獲得未加工序列。隨后使用QIIME 1.7.0軟件(http://qiime.org/scripts/split_libraries_fastq.html)對序列進(jìn)行質(zhì)控和過濾，獲得優(yōu)質(zhì)序列[24-25]。再通過UCHIME算法將其與參考數(shù)據(jù)庫(Gold database，http://drive5.com/uchime/uchime_download.html)進(jìn)行比對，發(fā)現(xiàn)并去除嵌合體序列，從而獲得有效序列[26-27]。

使用Uparse 7.0.1001軟件(http://drive5.com/uparse/)根據(jù)相似性≥97%原則將通過質(zhì)控的有效序列聚類成為操作分類單元(operational taxonomic unit，OTU)[28]。采用QIIME 1.7.0默認(rèn)參數(shù)計算各樣品的alpha多樣性指數(shù)和物種分布。為了比較2種豬糞便移植的小鼠盲腸內(nèi)容物的古菌構(gòu)成的差異，使用R 2.15.3軟件中的FactoMineR軟件包對樣品之間的古菌結(jié)構(gòu)相似度進(jìn)行主成分分析(PCoA)。使用ggplot 2軟件包繪制散點圖。

2 結(jié) 果

2.1 測序數(shù)據(jù)及OTU序列信息統(tǒng)計

由表1可見，經(jīng)過質(zhì)控、過濾、拼接、去除嵌合體序列等處理，本試驗最終獲得1 428 899條古菌有效序列，其中金華豬糞便樣品有60 979條有效序列，長白豬糞便樣品有61 369條有效序列，24例移植小鼠盲腸樣品中的古菌有效序列最少為46 362條，最多為61 524條，個體間有效序列數(shù)差異性不大。所獲得的古菌有效序列根據(jù)97%相似性水平進(jìn)行OTU聚類分析，金華豬糞便樣品最終獲得307個OTU，長白豬糞便樣品最終獲得202個OTU，24例移植小鼠盲腸樣品最終獲得104～133個OTU。各樣品的覆蓋度指數(shù)為0.999～1.000，表明樣本中序列沒有被測出的概率極低，因此本次測序結(jié)果能夠代表樣本的真實多樣性組成。

2.2 alpha多樣性分析

根據(jù)97%相似性水平下OTU的信息，使用Chao Ⅰ指數(shù)、ACE指數(shù)、香濃(Shannon)指數(shù)及辛普森(Simpson)指數(shù)對樣本內(nèi)古菌群落的豐富度和多樣性進(jìn)行評估，結(jié)果見表1。OTU數(shù)量(123±2 vs. 117±2，P=0.040 5)、Chao Ⅰ指數(shù)(130.180±11.084 vs. 122.663±8.538，P=0.076 2)和ACE指數(shù)(130.940±10.915 vs. 123.683±7.905，P=0.076 8)比較一致地顯示出金華豬糞便移植小鼠盲腸內(nèi)容物(JM1～12)的古菌群落豐富度略高于長白豬糞便移植小鼠(LM1～12)，但在Chao Ⅰ指數(shù)和ACE指數(shù)上差異不顯著(P>0.05)。Shannon指數(shù)(3.811±0.065 vs. 3.898±0.063，P=0.345 4)和Simpson指數(shù)(0.882±0.007 vs. 0.889±0.006，P=0.507 6)顯示，金華豬糞便移植小鼠盲腸內(nèi)容物(JM1～12)的古菌群落多樣性略高于長白豬糞便移植小鼠(LM1～12)，但差異不顯著(P>0.05)。

表1 樣品測序結(jié)果Table 1 Sequencing results of samples

續(xù)表1樣品Sample序列數(shù)量SequenceNo.OTU數(shù)量OTUNo.香濃指數(shù)Shannonindex辛普森指數(shù)SimpsonindexChaoⅠ指數(shù)ChaoⅠindexACE指數(shù)ACEindex覆蓋度指數(shù)CoverageindexLM3580751133.8350.900114.750116.4811.000LM4525121133.9020.900117.714120.8681.000LM5535791284.0880.896140.000137.3091.000LM6545161204.0190.882127.333127.3001.000LM7580311113.9990.905116.500117.4481.000LM8591971123.6390.881116.231119.0191.000LM9513281123.6430.879117.000118.3711.000LM10574611224.1950.917128.600128.6551.000LM11559141153.6480.866121.000123.4091.000LM12546971103.6710.842112.143112.3291.000

JD：金華豬糞便樣品feces sample ofJinhuapig；LD：長白豬糞便樣品 feces sample of Landrace pig；JM：金華豬糞便移植小鼠盲腸內(nèi)容物cecal contents of mice transplanted withJinhuapig feces；LM：長白豬糞便移植小鼠盲腸內(nèi)容物cecal contents of mice transplanted with Landrace pig feces。下表同 The same as below。

2.3 移植供體古菌結(jié)構(gòu)

將作為菌群移植供體的金華豬和長白豬糞便樣品中獲得的OTU分別在門和屬水平上進(jìn)行物種注釋，具體結(jié)果見表2和圖1?？梢钥闯觯驹囼炈玫降墓wOTU近1/2都屬于奇古菌門(Thaumarchaeota)，而廣古菌門(Euryarchaeota)和泉古菌門(Crenarchaeota)豐度相對較低。金華豬糞便樣品中廣古菌門豐度高于長白豬糞便樣品。在屬水平上，奇古菌門的古菌幾乎全部來源于餐古菌屬(Cenarchaeum)，長白豬糞便樣品中還有少量來源于亞硝化球菌屬(Nitrososphaera)；廣古菌門中，金華豬糞便樣品中甲烷短桿菌屬(Methanobrevibacter)豐度相對較高，此外還有極少量的待定屬CandidatusMethanomethylophilus，而長白豬糞便樣品中甲烷球菌屬(Methanosphaera)豐度較高。

表2 移植供體糞便在不同分類水平上的古菌組成及其豐度Table 2 Archaeal composition and their abundance at different classification levels in feces of transplant donors

A：門水平 phylum level；B：屬水平 genus level。

圖1 金華豬供體(JD)和長白豬供體(LD)糞便古菌結(jié)構(gòu)

Fig.1 Fecal archaea community structures ofJinhua(JD) and Landrace pig donors (LD)

2.4 移植小鼠盲腸內(nèi)容物古菌結(jié)構(gòu)

通過灌胃的方式，將金華豬和長白豬糞便微生物分別移植到抗生素處理的無菌小鼠體內(nèi)，接種后28 d取小鼠盲腸內(nèi)容物用于古菌結(jié)構(gòu)測序分析。所獲得的OTU分別在門和屬水平上進(jìn)行物種注釋，具體結(jié)果見表3和圖2?？梢钥闯?，在門水平，移植小鼠的古菌OTU中奇古菌門(Thaumarchaeota)的比例增加到70%以上，而廣古細(xì)菌門(Euryarchaeota)和泉古細(xì)菌門(Crenarchaeota)幾乎未檢出。而在屬水平上，奇古菌門的古菌幾乎全部來源于餐古菌屬(Cenarchaeum)。2組移植小鼠之間古菌結(jié)構(gòu)在門和屬水平上無顯著差異(P>0.05)，但與供體的古菌結(jié)構(gòu)差別較大。

表3 移植小鼠盲腸內(nèi)容物在不同分類水平上的古菌組成及其豐度Table 3 Archaeal composition and their abundance at different classification levels in cecal contents of transplanted mice

2.5 移植小鼠樣品的beta多樣性分析

根據(jù)各個樣品OTU的種類及其豐度計算各個樣品間的加權(quán)UniFrac距離(weighted UniFrac distance)，并基于此，對樣品進(jìn)行了主成分分析(圖3)?？梢钥闯?，除了JM12和LM12距離較遠(yuǎn)外，JM組和LM組內(nèi)其他樣品聚集度都較高，且2組在PC2維度上依據(jù)其組別顯示出明顯的組內(nèi)聚集。

A：移植小鼠門水平的古菌結(jié)構(gòu) archaea community structure of transplanted mice at phylum level；B：24只移植小鼠門水平的古菌平均值 mean of archaea of 24 transplanted mice at phylum level；C：移植小鼠屬水平的古菌結(jié)構(gòu) archaea community structure of transplanted mice at genus level；D：24只移植小鼠屬水平的古菌平均值 mean of archaea of 24 transplanted mice at genus level。

圖2 金華豬糞便移植小鼠(JM)和長白豬糞便移植小鼠(LM)盲腸內(nèi)容物古菌結(jié)構(gòu)

Fig.2 Archaea community structures of cecal contents of transplanted mice using feces formJinhua(JM) and Landrace pigs (LM)

圖3 金華豬糞便移植小鼠(JM)和長白豬糞便移植小鼠(LM)古菌結(jié)構(gòu)的主成分分析Fig.3 Principal coordinate analysis of archaea community structures between transplanted mice using feces form Jinhua (JM) and Landrace pigs (LM)

3 討 論

3.1 金華豬和長白豬糞便古菌結(jié)構(gòu)的異同

金華豬是我國著名的地方品種，生長速度慢，是脂肪型豬，而長白豬是優(yōu)良的外來品種，生長速度快，是瘦肉型豬[1]。對不同品種豬腸道微生物進(jìn)行跨物種移植以分析其功能具有重要的研究意義[19]，但腸道古菌在此過程中的作用目前研究還較少。本試驗通過糞菌移植，分別構(gòu)建了模擬2種豬腸道微生物環(huán)境的小鼠模型，采用高通量測序技術(shù)，分析金華豬和長白豬糞便古菌結(jié)構(gòu)差異，并探討了其糞菌移植小鼠的古菌結(jié)構(gòu)異同，為進(jìn)一步的功能及機(jī)制研究奠定了基礎(chǔ)。

以往的研究顯示，對動物和人腸道內(nèi)古菌的研究主要集中在產(chǎn)甲烷菌上。研究表明，動物腸道產(chǎn)甲烷菌的多樣性與其生存環(huán)境、發(fā)育情況及飼料等因素密切相關(guān)，因而，不同品種動物腸道內(nèi)產(chǎn)甲烷菌的結(jié)構(gòu)和組成存在一定差異[7-9]。在脂肪代謝方面，有報道指出，瘦肉型的長白豬糞便中產(chǎn)甲烷古菌的密度和多樣性均高于肥胖型二花臉豬，其原因在于甲烷的生成會影響能量獲取及體脂肪形成[7]；而Zhang等[29]的研究則顯示，肥胖個體腸道中利用氫氣(H2)的產(chǎn)甲烷古菌含量顯著高于瘦的或接受了胃腸旁路手術(shù)的個體，同時其產(chǎn)H2細(xì)菌含量也較高，細(xì)菌和古菌之間跨物種的H2轉(zhuǎn)移被認(rèn)為是肥胖個體大腸能量攝取增加的重要機(jī)制。由此可見，由于研究對象和研究方法的不同，對于動物和人腸道古菌的研究結(jié)果還存在一定的差異。而本研究顯示，廣古菌門雖然在2種豬糞便中含量均不高，但其在金華豬糞便中的豐度明顯高于長白豬糞便。其對于2種豬的生長和代謝的影響，還有待進(jìn)一步的研究。

奇古菌除了在某些白蟻和蠶中[13-14]以外，在人和動物消化道中的分布僅有少量報道，且含量極低[30]。但本研究卻發(fā)現(xiàn)，在金華豬及長白豬糞便的古菌中，奇古菌門是優(yōu)勢菌，占所有古菌的將近1/2，其中餐古菌為優(yōu)勢菌屬。這一結(jié)果與以往的研究不一致，其原因及意義還需要進(jìn)一步試驗驗證和探討。

3.2 糞菌移植過程中古菌結(jié)構(gòu)的變化

2006年John Rawls等研究人員將斑馬魚的腸道微生物群到無菌的小鼠腸道，結(jié)果發(fā)現(xiàn)移植后菌群的構(gòu)成發(fā)生了明顯變化，反之亦然。這表明，在跨物種動物間糞便微生物群移植過程中，供體微生物群會被宿主影響和重塑[20]。本試驗結(jié)果證實，在從豬到小鼠的糞菌移植過程中，古菌結(jié)構(gòu)同樣會受到宿主的影響。在本次試驗所選取的2個豬品種，金華豬糞便中廣古菌門的豐度明顯高于長白豬糞便，但移植到小鼠腸道內(nèi)后這種差異不再明顯，且從門到屬的水平上，古菌結(jié)構(gòu)均發(fā)生了較大的變化，奇古菌門的餐古菌屬優(yōu)勢更加明顯。雖然在在門和屬水平上，金華豬和長白豬糞便移植小鼠腸道古菌結(jié)構(gòu)無明顯變化，但在基于OTU水平的主成分分析中，二者有顯著不同的聚類，說明古菌結(jié)構(gòu)在OTU水平上存在一定的差異。

4 結(jié) 論

① 奇古菌門為金華豬和長白豬糞便古菌組成中的優(yōu)勢菌。

② 廣古菌門在金華豬糞便中的豐度明顯高于長白豬糞便。

③ 將2種豬的糞菌移植到廣譜抗生素處理過的小鼠腸道后，移植小鼠的盲腸內(nèi)容物古菌結(jié)構(gòu)與供體相比差異明顯，同時2種豬糞便微生物移植小鼠的古菌結(jié)構(gòu)在OTU水平上存在一定的差異。

[1] MIAO Z G,WANG L J,XU Z R,et al.Developmental changes of carcass composition,meat quality and organs in theJinhuapig and Landrace[J].Animal,2009,3(3):468-473.

[2] GUO J,SHAN T,WU T,et al.Comparisons of different muscle metabolic enzymes and muscle fiber types inJinhuaand Landrace pigs[J].Journal of Animal Science,2011,89(1):185-191.

[3] HARTSTRA A V,BOUTER K E,BCKHED F,et al.Insights into the role of the microbiome in obesity and type 2 diabetes[J].Diabetes Care,2015,38(1):159-165.

[4] KONTUREK P C,HAZIRI D,BRZOZOWSKI T,et al.Emerging role of fecal microbiota therapy in the treatment of gastrointestinal and extra-gastrointestinal diseases[J].Journal of Physiology and Pharmacology,2015,66(4):483-491.

[5] EME L,DOOLITTLE W F.Archaea[J].Current Biology,2015,25(19):R851-R855.

[6] GRIBALDO S,BROCHIER-ARMANET C.The origin and evolution of Archaea:a state of the art[J].Philosophical Transactions of the Royal Society of London.Series B:Biological Sciences,2006,361(1470):1007-1022.

[7] 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,2012:605289.

[8] SAMUEL B S,GORDON J I.A humanized gnotobiotic mouse model of host-archaeal-bacterial mutualism[J].Proceedings of the National Academy of Sciences of the United States of America,2006,103(26):10011-10016.

[9] SAMUEL B S,HANSEN E E,MANCHESTER J K,et al.Genomic and metabolic adaptations ofMethanobrevibactersmithiito the human gut[J].Proceedings of the National Academy of Sciences of the United States of America,2007,104(25):10643-10648.

[10] BROCHIER-ARMANET C,BOUSSAU B,GRIBALDO S,et al.Mesophilic Crenarchaeota:proposal for a third archaeal phylum,the Thaumarchaeota[J].Nature Reviews Microbiology,2008,6(3):245-252.

[11] YOU J,DAS A,DOLAN E M,et al.Ammonia-oxidizing archaea involved in nitrogen removal[J].Water Research,2009,43(7):1801-1809.

[12] PESTER M,SCHLEPER C,WAGNER M.The Thaumarchaeota:an emerging view of their phylogeny and ecophysiology[J].Current Opinion in Microbiology,2011,14(3):300-306.

[13] SHI Y,HUANG Z,HAN S,et al.Phylogenetic diversity of Archaea in the intestinal tract of termites from different lineages[J].Journal of Basic Microbiology,2015,55(8):1021-1028.

[14] LI G N,XIA X J,TANG W C,et al.Intestinal microecology associated with fluoride resistance capability of the silkworm (BombyxmoriL.)[J].Applied Microbiology and Biotechnology,2016,100(15):6715-6724.

[15] LESZCZYSZYN J J,RADOMSKI M,LESZCZYSZYN A M.Intestinal microbiota transplant-current state of knowledge[J].Reumatologia,2016,54(1):24-28.

[16] BOJANOVA D P,BORDENSTEIN S R.Fecal transplants:what is being transferred?[J].PLoS Biology,2016,14(7):e1002503.

[17] TURNBAUGH P J,LEY R E,MAHOWALD M A,et al.An obesity-associated gut microbiome with increased capacity for energy harvest[J].Nature,2006,444(7122):1027-1031.

[18] KULECKA M,PAZIEWSKA A,ZEBER-LUBECKA N,et al.Prolonged transfer of feces from the lean mice modulates gut microbiota in obese mice[J].Nutrition & Metabolism,2016,13:57.

[19] DIAO H,YAN H L,XIAO Y,et al.Intestinal microbiota could transfer host Gut characteristics from pigs to mice[J].BMC Microbiology,2016,16:238.

[20] RAWLS J F,MAHOWALD M A,LEY R E,et al.Reciprocal gut microbiota transplants from zebrafish and mice to germ-free recipients reveal host habitat selection[J].Cell,2006,127(2):423-433.

[21] MCFALL-NGAI M.Love the one you’re with:vertebrate guts shape their microbiota[J].Cell,2006,127(2):247-249.

[22] WANG Z N,KLIPFELL E,BENNETT B J,et al.Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease[J].Nature,2011,472(7341):57-63.

[23] PANG X Y,HUA X G,YANG Q,et al.Inter-species transplantation of gut microbiota from human to pigs[J].The ISME Journal,2007,1(2):156-162.

[24] CAPORASO J G,KUCZYNSKI J,STOMBAUGH J,et al.QIIME allows analysis of high-throughput community sequencing data[J].Nature Methods,2010,7(5):335-336.

[25] BOKULICH N A,SUBRAMANIAN S,FAITH J J,et al.Quality-filtering vastly improves diversity estimates fromIlluminaampliconsequencing[J].Nature Methods,2013,10(1):57-59.

[26] EDGAR R C,HAAS B J,CLEMENTE J C,et al.UCHIME improves sensitivity and speed of chimera detection[J].Bioinformatics,2011,27(16):2194-2200.

[27] HAAS B J,GEVERS D,EARL A M,et al.Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons[J].Genome Research,2011,21(3):494-504.

[28] EDGAR R C.UPARSE:highly accurate OTU sequences from microbial amplicon reads[J].Nature Methods,2013,10(10):996-998.

[29] ZHANG H,DIBAISE J K,ZUCCOLO A,et al.Human gut microbiota in obesity and after gastric bypass[J].Proceedings of the National Academy of Sciences of the United States of America,2009,106(7):2365-2370.

[30] GACI N,BORREL G,TOTTEY W,et al.Archaea and the human gut:new beginning of an old story[J].World Journal of Gastroenterology,2014,20(43):16062-16078.

*Corresponding author, senior livestock engineer, E-mail: yanghua806@hotmail.com

(責(zé)任編輯 王智航)

Intestinal Archaea Community Structure Analysis of Mice Transplanted withJinhuaand Landrace Pig Feces

XIAO Yingping1WANG Junjun2,3LI Tiantian2,3TANG Biao1HE Xiangxiang1LI Kaifeng1YANG Hua1*

(1.InstituteofQualityandStandardforAgro-Products,ZhejiangAcademyofAgriculturalSciences,Hangzhou310021，China; 2.BeijingAdvancedInnovationCenterforFoodNutritionandHumanHealth,Beijing100193,China; 3.CollegeofAnimalScienceandTechnology,ChinaAgriculturalUniversity,Beijing100193,China)

This experiment was conducted to explore the differences of the intestinal archaea community structures between the obese breedJinhuapig and the lean breed Landrace pig, and to investigate the changes of the archaea community structure during the inter-species fecal microbiota transplantation (FMT). The fresh feces ofJinhuaand Landrace pigs were collected and introduced into broad spectrum antibiotics treated germ-free mice guts by gavage, respectively. The total genomic DNA was extracted from the pig fecal samples and the transplanted mice cecal contents, and was used as templates to amplify the V4 region of the 16S rRNA gene by PCR with the universal primers of archaea. The PCR products were then subjected to high-throughput sequencing on an Illumina HiSeq sequencing platform． The obtained sequences were subsequently processed and analyzed with a series of software such as QIIME. The results showed as follows: at the phylum level, three main phyla of archaea, Thaumarchaeota, Euryarchaeota, and Crenarchaeota, were found in pig fecal samples and transplanted mice cecal contents, with Thaumarchaeota as the dominant phylum; while at the genus level, five known genera of archaea were observed, withCenarchaeumas the dominant genus. The abundance of phylum Euryarchaeota inJinhuapig’s feces was higher than that in Landrace pig’s feces. Thaumarchaeota dominate the archaea community structure of transplanted mice cecal contents, with abundance of more than 70%, and Thaumarchaeota mainly consisted ofCenarchaeum. It was observed that no significant difference was exhibited between the two groups of mice at the phylum and genus levels (P>0.05), but the operational taxonomic unit (OTU) based principal coordinate analysis showed a good aggregation within each group. It is concluded that genusCenarchaeumof phylum Thaumarchaeota is the dominant archaea in feces ofJinhuaand Landrace pigs and cecal contents of mice recipients. Phylum Euryarchaeota is more abundant inJinhuapig’s feces than in Landrace pig’s feces， and there is difference in archaea community structures between their respective mice recipients at the OTU level.[ChineseJournalofAnimalNutrition, 2017, 29(6)：1895-1903]

archaea;Jinhuapig; Landrace pig; fecal microbiota transplantation; high-throughput sequencing

10.3969/j.issn.1006-267x.2017.06.010

2016-12-01

浙江省自然科學(xué)基金重點項目(LZ15C170001)；北京食品營養(yǎng)與人類健康高精尖創(chuàng)新中心開放基金；國家自然科學(xué)基金(31402083)

肖英平(1984—)，男，江西興國人，助理研究員，博士，研究方向為畜禽腸道微生物。E-mail: ypxiaozju@126.com

*通信作者:楊 華，高級畜牧師，碩士生導(dǎo)師，E-mail: yanghua806@hotmail.com

S828

A

1006-267X(2017)06-1895-09