李曉斌 趙國棟 劉 振 吳婷婷 何周瑞 鄧海峰 楊開倫(新疆農(nóng)業(yè)大學，新疆肉乳用草食動物營養(yǎng)實驗室，烏魯木齊830052)

3～6月齡伊犁馬腸道微生物群落多樣性的研究

李曉斌 趙國棟 劉 振 吳婷婷 何周瑞 鄧海峰 楊開倫*
(新疆農(nóng)業(yè)大學，新疆肉乳用草食動物營養(yǎng)實驗室，烏魯木齊830052)

本試驗旨在研究3～6月齡伊犁馬腸道微生物群落的多樣性，揭示斷奶前馬駒腸道菌群演替的規(guī)律，從微生物學角度為此階段馬駒生長發(fā)育及腸道健康研究提供理論依據(jù)。試驗選取平均體重為(89.75±8.81) kg、出生日期相同的3月齡伊犁馬5匹，進行為期90 d的飼養(yǎng)試驗。分別在試驗的第0天(即試驗開始的前1天)、第30天、第60天、第90天采集馬駒糞便樣品，從每份樣品中提取微生物基因組總DNA，采用Illumina HiSeq測序技術檢測樣品微生物群落多樣性。結(jié)果顯示：1)對5匹馬駒的20份糞樣進行測序，共獲得有效序列數(shù)157 665條，獲得平均操作分類單元(OTUs)1 117個。2)馬駒糞便中微生物群落α多樣性指數(shù)(ACE、Chao1、Shannon和Simpson指數(shù))隨著馬駒月齡的增加呈波動性變化，但各時間點之間均無顯著性變化(P>0.05)。3)在門水平上，馬駒糞便中10大優(yōu)勢菌分別為厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidetes)、變形菌門(Proteobacteria)、疣微菌門(Verrucomicrobia)、螺旋體菌門(Spirochaetes)、梭桿菌門(Fusobacteria)、無壁菌門(Tenericutes)、放線菌門(Actinobacteria)、TM7、廣古菌門(Euryarchaeota)，其中以厚壁菌門、擬桿菌門、變形菌門及疣微菌門的豐度較高；在科水平上，馬駒糞便中10大優(yōu)勢菌分別為芽孢桿菌科(Bacillaceae)、莫拉氏菌科(Moraxellaceae)、動球菌科(Planococcaceae)、肉桿菌科(Carnobacteriaceae)、BS11、RFP12、乳桿菌科(Lactobacillaceae)、瘤胃球菌科(Ruminococcaceae)、毛螺旋菌科(Lachnospiraceae)、紫單胞菌科(Porphyromonadaceae)；在屬水平上，馬駒糞便中10大優(yōu)勢菌分別為不動桿菌屬(Acinetobacter)、德庫菌屬(Desemzia)、乳桿菌屬(Lactobacillus)、Ureibacillus、Paludibacter、芽孢桿菌屬(Bacillus)、埃希式桿菌屬(Escherichia)、肉食桿菌屬(Carnobacterium)、密螺旋體屬(Treponema)、艱難桿菌屬(Mogibacterium)。由此得出，采用Illumina HiSeq測序技術能夠準確地對3～6月齡伊犁馬腸道微生物群落進行分類研究；厚壁菌門、擬桿菌門、變形菌門、疣微菌門是3～6月齡伊犁馬腸道內(nèi)的優(yōu)勢菌門。

伊犁馬；Illumina HiSeq測序技術；腸道微生物群落；多樣性

馬是單胃動物，其盲腸相當于反芻動物的瘤胃，可容納28～36 L的消化物和消化液。通過盲腸內(nèi)容物或糞便培養(yǎng)檢測發(fā)現(xiàn)，其微生物種類有400～500種，數(shù)量可達1010～1012個/g[1-2]。定植于腸道內(nèi)的微生物參與腸道營養(yǎng)物質(zhì)的消化、吸收及合成，是動物腸道重要的組成部分。因此，促進馬匹腸道微生物區(qū)系的建立，維持腸道正常菌群的平衡對馬匹消化吸收、腸道健康和免疫調(diào)控意義重大[3]。馬駒出生時腸道是無菌的，在分娩、哺乳及與周圍環(huán)境接觸過程逐漸感染細菌。隨著馬駒年齡、生活環(huán)境、飼糧結(jié)構(gòu)、健康狀況等的改變，腸道微生物的種類、結(jié)構(gòu)、豐度隨之發(fā)生變化[4]。在應激條件下這些因素的改變會導致馬駒腸道微生物區(qū)系紊亂，并引起腹瀉，嚴重時甚至死亡。腸道菌群的平衡在馬的生長發(fā)育、健康狀況方面扮演著重要角色[5]。早期補飼對馬駒生長發(fā)育意義重大，過早的攝入補喂料對馬駒腸道微生物結(jié)構(gòu)改變產(chǎn)生一定的影響。3月齡馬駒隨著母乳減少逐漸開始增加補喂料，飼糧由液態(tài)逐漸轉(zhuǎn)為固態(tài)，其腸道微生物的種類和豐度也隨之發(fā)生變化，從而影響馬駒腸道健康和生長發(fā)育。鑒于此，本試驗以3～6月齡的伊犁馬為研究對象，通過檢測糞便中微生物群落的多樣性，研究馬駒從3月齡到6月齡內(nèi)腸道微生物定植、演替的規(guī)律，為研究開發(fā)適宜馬駒健康的微生物制劑或腸道微生物疾病防治等提供參考依據(jù)。

1 材料與方法

1.1 試驗時間及地點

試驗時間：2014年7月至2014年10月；試驗地點：新疆伊犁哈薩克自治州昭蘇馬場。

1.2 試驗設計

本試驗選取平均體重為(89.75±8.81) kg、出生日期相同的3月齡伊犁馬哺乳馬駒5匹。5匹馬駒與其母馬來自同一放牧草場，且母馬年齡(11歲)、胎次(7胎)均相同。在哺乳條件下，每天每匹馬駒按照體重的0.6%補喂精料補充料(每30 d調(diào)整1次補喂量，具體見表1)，并投喂等量的優(yōu)質(zhì)苜蓿干草，進行為期90 d的飼養(yǎng)試驗。分別在試驗的第0天(即試驗開始的前1天)、第30天、第60天、第90天12:00前采集馬駒糞便樣品。

1.3 飼養(yǎng)管理及飼糧組成

試驗馬駒與母馬均飼養(yǎng)于室外馬舍(長75 m×寬40 m)。馬駒08:00—18:00拴系飼養(yǎng)，其余時間自由活動。將馬駒全天的精料補充料分成3等份，分別在08:00、13:00、17:00使用料兜(料兜由本試驗設計，其大小和深淺符合馬駒采食精料，確保不影響馬駒呼吸及防止馬駒仰頭將精料撒出)補喂。所有試驗馬駒每天提供等量[(1.00±0.25) kg]的優(yōu)質(zhì)苜蓿干草，自由飲水。樣品采集當天馬駒室內(nèi)圈舍飼養(yǎng)，08:00將馬駒牽入室內(nèi)圈舍，單槽單欄位飼養(yǎng)，防止馬駒糞便互相混合。飼喂后打掃圈舍，為12:00前采集糞樣做準備。精料補充料組成及營養(yǎng)水平見表2。

表1 精料補充料補喂量

1.4 馬駒糞便樣品采集

于試驗的第0天、第30天、第60天、第90天的08:00將馬駒牽入室內(nèi)圈舍，單欄位拴系并套帶收糞袋。將每匹馬駒自進圈后至12:00排出的全部糞便收集于集糞桶內(nèi)，混合均勻后，采用5點取樣法采集馬駒糞樣，每個點采集5 g，再將5個點共25 g糞樣混合均勻后立即裝入密封袋，-20 ℃冷凍保存，待測。

1.5 微生物測定樣品分組編號

在后續(xù)樣品測定及分析中分別將試驗的第0天、第30天、第60天、第90天的樣品記為4個組，每個時間點5匹馬駒的樣品視為1個組。

1.6 糞樣高通量測序分析

總DNA提取、PCR擴增及Illumina HiSeq測序及結(jié)果分析均由北京諾禾致源生物信息科技有限公司協(xié)助完成。

1.7 數(shù)據(jù)處理

Metastats分析使用R軟件在各分類水平下做組間的置換檢驗(permutation test)，得到P值，然后利用Benjamini and Hochberg False Discovery Rate方法對于P值進行修正。組間差異顯著的物種分析利用R軟件做組間t檢驗(t-test)。

表2 精料補充料組成及營養(yǎng)水平(干物質(zhì)基礎)

1)預混料為每千克精料補充料提供The premix provided the following per kg of the concentrate supplement：VA 3 000 IU，VB120 mg，VB220 mg，VB66 mg，VC 20 mg，VD 1 000 IU，VE 500 IU，泛酸 pantothenic acid 10 mg，煙酰胺 nicotinamide 100 mg，Cu 25 mg，F(xiàn)e 107 mg，Mn 81 mg，Zn 74 mg，I 6 mg，Se 14 mg，Co 3 mg，氯化膽堿 choline chloride 120 mg。

2)營養(yǎng)水平為實測值。Nutrient levels were measured values.

2 結(jié)果與分析

2.1 測序結(jié)果及α多樣性分析

測序結(jié)果及α多樣性數(shù)據(jù)見表3。對5匹馬駒的20份糞樣進行測序，共獲得有效序列數(shù)157 665條，以97%的序列一致性(identity)，將序列共聚類成4 471個操作分類單元(OTUs)。4個時間點的測序覆蓋率均達到96%以上，說明本試驗中獲得的序列代表了微生物界96%以上的細菌種系型，測序覆蓋率好，且各時間點均觀察到較豐富的物種。馬駒糞便中微生物群落α多樣性指數(shù)(ACE、Chao1、Shannon和Simpson指數(shù))隨著馬駒月齡的增加呈波動性變化，但各時間點之間均無顯著性差異(P>0.05)。

2.2 β多樣性分析

基于OTUs水平的主成分分析(PCA)結(jié)果見圖1。圖中主成分1(PC1)為第1主坐標，對檢測到的總微生物的代表性貢獻率為17.28%，主成分2(PC2)為第2主坐標，貢獻率為9.68%；4個組20個樣品能夠較明顯的區(qū)分開，第0天、第30天能夠明顯與第60天、第90天分開，第60天與第90天能夠明顯分開，說明同一時間點微生物群落的組成相似。

圖2左側(cè)顯示的是非加權組平均法(UPGMA)聚類樹結(jié)構(gòu)，右側(cè)顯示的是各組在門水平上的物種相對豐度分布。可以看出，第0天、第30天聚為一類，第60天、第90天聚為一類，說明隨著馬駒月齡的增加，微生物群落趨于復雜化和多樣化。

表3 測序結(jié)果及α多樣性分析

同列數(shù)據(jù)肩標無字母或相同字母表示差異不顯著(P>0.05)，不同小寫字母表示差異顯著(P<0.05)，不同大寫字母表示差異極顯著(P<0.01)。

In the same column, values with no letter or the same letter superscripts mean no significant difference (P>0.05), while with different small letter superscripts mean significant difference (P<0.05), and with different capital letter superscripts mean significant difference (P<0.01).

圖1 主成分分析結(jié)果

2.3 馬駒腸道微生物構(gòu)成及分布豐度

綜合分析測序樣品的序列數(shù)，本試驗研究了4個時間點馬駒糞便中微生物群落結(jié)構(gòu)，同時各時間點之間進行了比較。根據(jù)物種注釋結(jié)果，選取每個樣品在各分類水平上最大豐度排名前10的物種，生成物種相對豐度柱形累加圖，以便直觀查看各樣品在不同分類水平上相對豐度較高的物種及其比例。

2.3.1 3～6月齡馬駒門水平上腸道微生物豐度

在門水平上馬駒腸道微生物豐度數(shù)據(jù)見表4。由表中數(shù)據(jù)可知，馬駒糞便中10大優(yōu)勢菌門分別為厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidetes)、變形菌門(Proteobacteria)、疣微菌門(Verrucomicrobia)、螺旋體菌門(Spirochaetes)、無壁菌門(Tenericutes)、梭桿菌門(Fusobacteria)、放線菌門(Actinobacteria)、TM7、廣古菌門(Euryarchaeota)，占菌群分類的98%以上，1%左右的為其他未被分類的菌門，門水平上腸道微生物豐度柱狀圖如圖3所示。其中厚壁菌門、擬桿菌門、變形菌門及疣微菌門為優(yōu)勢菌門。隨著月齡的增加，馬駒糞便中各微生物的豐度發(fā)生了不同程度的變化。厚壁菌門的豐度在試驗第0天、第30天、第60天和第90天時分別為33.19%、32.90%、73.72%和67.89%，第60天時極顯著高于第0天和第30天時(P<0.01)，第90天時顯著高于第0天和第30天時(P<0.05)。放線菌門的豐度在試驗第60天時達到1.34%，顯著高于試驗第0天和第30天時(P<0.05)。無壁菌門的豐度隨著馬駒月齡的增加而增加，在試驗第90天時為3.05%，顯著高于第0天和第30天時(P<0.05)。擬桿菌門的豐度在試驗第0天、第30天、第60天和第90天時分別為40.97%、43.44%、0.49%和2.75%，第60天和第90天時極顯著低于第0天和第30天時(P<0.01)。螺旋體菌門的豐度隨著馬駒月齡的增加先降低后趨于平穩(wěn)，試驗第0天時為3.13%，與之相比，在試驗第30天、第60天和第90天時顯著降低(P<0.05)。變形菌門、疣微菌門、梭桿菌門、TM7、廣古菌門的豐度在整個試驗期無顯著性變化(P>0.05)，其中變形菌門的豐度在試驗第30天時降低，第60天和第90天時又有所升高，與第0天時接近。

Firmicutes：厚壁菌門；Bacteroidetes：擬桿菌門；Proteobacteria：變形菌門；Verrucomicrobia：疣微菌門；Spirochaetes：螺旋體菌門； Tenericutes：無壁菌門；Fusobacteria：梭桿菌門；Actinobacteria：放線菌門；Euryarchaeota：廣古菌門；Others：其他。圖3同 The same as Fig.3.

圖2 基于Weighted Unifrac距離的UPGMA聚類樹

同行數(shù)據(jù)肩標無字母或相同字母表示差異不顯著(P>0.05)，不同小寫字母表示差異顯著(P<0.05)，不同大寫字母表示差異極顯著(P<0.01)。下表同。

In the same row, values with no letter or the same letter superscripts mean no significant difference (P>0.05), while with different small letter superscripts mean significant difference (P<0.05), and with different capital letter superscripts mean significant difference (P<0.01). The same as below.

圖3 3～6月齡馬駒門水平上腸道微生物豐度柱狀圖

2.3.2 3～6月齡馬駒科水平上腸道微生物豐度

在科水平上馬駒腸道微生物豐度數(shù)據(jù)見表5。由表中數(shù)據(jù)可知，馬駒糞便中10大優(yōu)勢菌科分別為芽孢桿菌科(Bacillaceae)、莫拉氏菌科(Moraxellaceae)、動球菌科(Planococcaceae)、肉桿菌科(Carnobacteriaceae)、BS11、RFP12、乳桿菌科(Lactobacillaceae)、瘤胃球菌科(Ruminococcaceae)、毛螺旋菌科(Lachnospiraceae)、紫單胞菌科(Porphyromonadaceae)，科水平上腸道微生物豐度柱狀圖如圖4所示。芽孢桿菌科的豐度在試驗第60天時為22.31%，極顯著高于第0天和第30天時(P<0.01)，顯著高于第90天時(P<0.05)。乳桿菌科的豐度在試驗第90天時為5.77%，極顯著高于第0天和第30天時(P<0.01)。瘤胃球菌科的豐度在試驗第0天、第30天、第60天和第90天時分別為11.93%、10.19%、3.44%和7.66%，第60天時極顯著低于第0天和第30天時(P<0.01)，顯著低于第90天時(P<0.05)。毛螺旋菌科的豐度在試驗第0天、第30天、第60天和第90天時分別為6.98%、6.45%、2.96%和6.88%，第60天時顯著低于第0天、第30天和第90天時(P<0.05)。紫單胞菌科的豐度在試驗第0天、第30天、第60天和第90天時分別為2.58%、0.40%、0.02%和0.14%，第60天時極顯著低于第0天和第30天時(P<0.01)。BS11的豐度在試驗第0天、第30天、第60天和第90天時分別為8.95%、15.98%、0.06%和0.53%，第60天和第90天時顯著低于第0天和第30天時(P<0.05)。莫拉氏菌科、動球菌科、肉桿菌科、RFP12的豐度在整個試驗期無顯著性變化(P>0.05)

表5 3～6月齡馬駒科水平上腸道微生物豐度

Others：其他；Porphyromonadaceae：紫單胞菌科；Lachnospiraceae：毛螺旋菌科；Ruminococcaceae：瘤胃球菌科；Lactobacillaceae：乳桿菌科；Carnobacteriaceae：肉桿菌科；Planococcaceae：動球菌科；Moraxellaceae：莫拉氏菌科；Bacillaceae：芽孢桿菌科。

圖4 3～6月齡馬駒科水平上腸道微生物豐度柱狀圖

Fig.4 The histograms of abundances of micrbiota in intestine of 3- to 6-month-old colts at family level

2.3.3 3～6月齡馬駒屬水平上腸道微生物豐度

在屬水平上馬駒腸道微生物豐度數(shù)據(jù)見表6。由表中數(shù)據(jù)可知，馬駒糞便中10大優(yōu)勢菌屬分別是不動桿菌屬(Acinetobacter)、德庫菌屬(Desemzia)、乳桿菌屬(Lactobacillus)、Ureibacillus、Paludibacter、芽孢桿菌屬(Bacillus)、埃希式桿菌屬(Escherichia)、肉食桿菌屬(Carnobacterium)、密螺旋體屬(Treponema)、艱難桿菌屬(Mogibacterium)，屬水平上腸道微生物豐度柱狀圖如圖5。乳桿菌屬的豐度在試驗第90天時為5.77%，極顯著高于第0天、第30天和第60天時(P<0.01)。芽孢桿菌屬的豐度在試驗第60天時為4.11%，極顯著高于第0天和第30天時(P<0.01)。

3 討 論

3.1 3～6月齡馬駒腸道微生物群落多樣性

近年，Illumina高通量測序技術在微生物研究領域廣泛使用。利用Illumina高通量測序技術能夠更深入、直觀、準確地通過腸道微生物群落的穩(wěn)定性和多樣性的變化反映動物的胃腸道健康和生長發(fā)育情況[3,6-8]。Proughbred等[9]對8匹去勢純血馬的14份糞樣進行測序，共獲得有效序列488 213條，OTUs的平均范圍為1 200～3 000個。Yatsunenk等[10]研究顯示，應用高通量測序技術在人糞便中共檢測出97%以上覆蓋率的微生物細菌，OTUs超過2 000個。Costa等[11]在11匹馬糞便中共獲得有效序列數(shù)6 536 523條。本試驗結(jié)果顯示，通過Illumina HiSeq測序獲得5匹馬駒4個階段20個糞樣細菌16S rDNA有效序列數(shù)157 665條，且隨著馬駒月齡的增加獲得有效序列數(shù)在增加，從試驗第0天的23 549條增加到第90天的53 617條；平均獲得1 117個OTUs，且OTUs數(shù)也隨著馬駒月齡的增加在增加，從試驗初的1 005個增加到試驗結(jié)束時的1 347個。通過α多樣性分析，本試驗結(jié)果顯示，隨著馬駒月齡的增加，糞便樣品的α多樣性指數(shù)中ACE、Chao1、Shannon指數(shù)逐漸升高，Simpson指數(shù)則隨著馬駒月齡的增加逐漸降低，說明隨著馬駒月齡的增長、飼糧的變更，馬駒腸道微生物種類越來越豐富，與Sakaitani等[4]報道的結(jié)果一致。隨著馬駒月齡的增加，馬駒采食的固體飼草料成為腸道系統(tǒng)營養(yǎng)物質(zhì)的主要來源，飼糧類型的變更是影響腸道微生物結(jié)構(gòu)變化的主要因素[12]。本試驗中，馬駒逐漸用精料補充料及苜蓿干草替代母乳，隨著母馬泌乳量的減少，馬駒在補喂定量精料補充料的基礎上，采食的粗飼料(苜蓿干草)增加，飼糧類型逐漸改變。因此，隨著馬駒的生長，馬駒糞便中檢測出微生物的種類和數(shù)量更豐富、更多樣、更復雜。該結(jié)果與Combes等[13]在家兔上的研究結(jié)果一致，即日齡對腸道菌群結(jié)構(gòu)、多樣性和復雜性有顯著影響，其原因與飼糧的改變有直接的關系。

表6 3～6月齡馬駒屬水平上腸道微生物豐度

Others：其他；Mogibacterium：艱難桿菌屬；Treponema：密螺旋體屬；Carnobacterium：肉食桿菌屬；Escherichia：埃希式桿菌屬；Bacillus：芽孢桿菌屬；Lactobacillus：乳桿菌屬；Desemzia：德庫菌屬；Acinetobacter：不動桿菌屬。

Others：其他；Mogibacterium：艱難桿菌屬；Treponema：密螺旋體屬；Carnobacterium：肉食桿菌屬；Escherichia：埃希式桿菌屬；Bacillus：芽孢桿菌屬；Lactobacillus：乳桿菌屬；Desemzia：德庫菌屬；Acinetobacter：不動桿菌屬。

圖5 3～6月齡馬駒屬水平上微生物豐度柱狀圖

Fig.5 The histograms of abundances of micrbiota in intestine of 3- to 6-month-old colts by genera level

3.2 3～6月齡馬駒腸道微生物群落結(jié)構(gòu)的變化

馬腸道微生物定植的研究相對于反芻動物較少。馬駒胎糞中沒有細菌，與馬駒腸道無菌結(jié)果相一致[4]。乳酸桿菌、梭狀芽胞桿菌、腸桿菌在馬駒出生后第3天從糞便中分離出來，12周齡的馬駒腸道纖維水解菌和乳酸菌的水平和成年馬相似[14]。當馬駒到42日齡時，在其糞便中檢測到比母馬更豐富的微生物種群，且與母馬的相似度達到60%。在其他一些研究中，馬駒在42～80日齡與母馬腸道微生物種群相似度達到50%。因此，研究馬駒腸道微生物區(qū)系早期的建立對馬駒腸道健康及生長發(fā)育意義重大。

許宇靜等[15]采用Illumina高通量測序技術研究斷奶幼兔盲腸微生物區(qū)系，發(fā)現(xiàn)其在門水平上主要由厚壁菌門、擬桿菌門、變形菌門和疣微菌門組成；在科水平上主要由瘤胃球菌科、毛螺旋菌科、擬桿菌科和紫單胞桿菌科組成。白秀娟等[16]利用聚合酶鏈式反應-變性梯度凝膠電泳(PCR-DGGE)技術分析發(fā)現(xiàn)，斷奶仔兔腸道主要優(yōu)勢菌為厚壁菌門、毛螺旋菌科等。谷莉[17]采用454高通量測序技術分析發(fā)現(xiàn)，不同飲食結(jié)構(gòu)條件下小鼠腸道厚壁菌門和擬桿菌門為主要優(yōu)勢菌群。本試驗結(jié)果顯示，采用Illumina HiSeq測序技術在馬駒糞便中檢測到豐度占99%以上門為厚壁菌門、擬桿菌門、變形菌門、疣微菌門、螺旋體菌門、無壁菌門、放線菌門、TM7、廣古菌門，從豐度來看厚壁菌門、擬桿菌門、變形菌門、螺旋體菌門為主要的優(yōu)勢菌門。與以上結(jié)果對比發(fā)現(xiàn)，厚壁菌門和擬桿菌門是單胃大型動物和小型動物腸道主要的優(yōu)勢菌。本試驗中，馬駒在門水平上發(fā)現(xiàn)較多的菌群，這與種間差異有關。Dougal等[6]研究結(jié)果顯示，在5匹純血馬和5匹矮種馬的大腸內(nèi)容物中檢測到的主要優(yōu)勢菌為厚壁菌門、擬桿菌門、螺旋體菌門、變形菌門和放線菌門。由此說明在相同的種屬間腸道微生物菌群組成更接近。

此外，在鼠、豬、馬、牛等的研究中發(fā)現(xiàn)，厚壁菌門和擬桿菌門是其胃腸道主要的優(yōu)勢菌群[18-20]，且厚壁菌門是攝食纖維飼料動物胃腸道促進纖維分解的主要菌群[21]。厚壁菌門能從纖維飼料中獲取更多能量滿足動物生長需要[22]。胃腸道菌群的豐度與消化能力有直接的關系[23-24]。厚壁菌門為馬駒優(yōu)勢菌群，這與盲腸、結(jié)腸為馬屬動物主要的消化場所有關。擬桿菌門雖然隨著馬駒月齡的增加豐度有所降低，但仍然是馬駒腸道的優(yōu)勢菌群。擬桿菌門是草食動物利用碳水化合物的主要菌群[25]，這與馬駒攝入纖維飼料有直接的關系。在單胃動物豬的研究中，飼糧中添加抗生素會減少擬桿菌門的數(shù)量[26]。本試驗精料補充料中并未添加有抗生素成分的藥品，但在試驗前母馬使用過驅(qū)蟲藥(主要成分為碘醚柳胺)，是否因為這個原因?qū)е埋R駒腸道擬桿菌門豐度降低有待于進一步的研究。

本試驗中，隸屬于厚壁菌門的芽孢桿菌科、乳桿菌科和瘤胃球菌科，隸屬于擬桿菌門的紫單胞菌科及螺旋體菌門的毛螺旋菌科隨著馬駒月齡的增加呈顯著性變化。研究顯示，草食動物反芻或非反芻類腸道內(nèi)包括毛螺旋科、纖維桿菌屬等在內(nèi)的細菌與其食草特性相關[27-28]。12周齡(3月齡左右)的馬駒腸道纖維降解菌數(shù)量和母馬相似[4]。在科水平上這些細菌豐度的變化與飼糧纖維含量的增加有直接的關系。毛螺旋菌科的細菌參與纖維素的降解[29]，其主要作用是降解植物飼料細胞壁中的果膠、果糖、纖維二糖等，產(chǎn)生揮發(fā)性脂肪酸(乙酸、丙酸、丁酸等)，為動物及腸道微生物提供能量。研究表明，飼糧變更是導致腸道微生物改變的首要原因[12]。在本試驗中，隨著母馬泌乳量減少，馬駒采食固體飼糧的量逐漸增加，故而纖維攝入量增加，這可能是促進與纖維降解相關菌種豐度增加的原因。

此外，變形菌門是動物消化道的優(yōu)勢菌門，其菌群龐大、種類多樣，在臨床研究中對胃腸道疾病診斷意義重大。很多致病菌包括大腸桿菌、沙門氏桿菌、霍亂弧菌、幽門螺桿菌都屬于變形菌門的γ-變形綱[30]。本試驗結(jié)果顯示，變形菌門在整個試驗期間豐度沒有發(fā)生顯著性變化，在科、屬水平上排在前10的優(yōu)勢菌群中未發(fā)現(xiàn)屬于變形菌門的有害菌。由此說明，隨著馬駒月齡的增加，厚壁菌門、擬桿菌門等有益菌豐度的增加抑制或穩(wěn)定了屬于變形菌門的有害菌豐度的增加，從而維持馬駒腸道微生物細菌的穩(wěn)定性，保證馬駒腸道的健康。

4 結(jié) 論

① 本試驗采用Illumina HiSeq測序技術獲得獲得5匹馬駒4個時間點的20個糞樣的有效序列數(shù)為157 665條，平均OTUs 1 117個，滿足測序條件。

② 厚壁菌門、擬桿菌門、變形菌門、疣微菌門為3～6月齡伊犁馬腸道主要的優(yōu)勢菌群。

[1] PIETZAK M.Bacterial colonization of the neonatal gut[J].Journal of Pediatric Gastroenterology and Nutrition,2004,38(4):389-391.

[2] GORBACK S L,PLAUT A G,HAHAS L,et al.Studies of intestinal microflora.Ⅱ.Microorganisms of the small intestine and their relations to oral and fecal flora[J].Gastroenterology,1967,53(6):856-867.

[3] DALY K,STEWART C S,FLINT H J,et al.Bacterial diversity within the equine large intestine as revealed by molecular analysis of cloned 16S rRNA genes[J].FEMS Microbiology Ecology,2001,38(2/3):141-151.

[4] SAKAITANI Y,YUKI N,NAKAJIMA F,et al.Colonization of intestinal microflora in newborn foals[J].Journal of Intestinal Microbiology,1999,13(1):9-14.

[5] BLASER M J,FALKOW S.What are the consequences of the disappearing human microbiota[J].Nature Reviews Microbiology,2009,7(12):887-894.

[6] DOUGAL K,DE AL FUENTE G,HARRIS P A,et al.Identification of a core bacterial community within the large intestine of the horse[J].PLoS One,2013,8(10):e77660.

[7] O’DONNELL M M,HARRIS H M B,JEFFERY I B,et al.The core faecal bacterial microbiome of Irish Thoroughbred racehorses[J].Letters in Applied Microbiology,2013,57(6):429-501.

[8] STEWART C J,MARRS E C L,MAGORRIAN S,et al.The preterm gut microbiota:changes associated with necrotizing enterocolitis and infection[J].Acta Paediatrica,2012,101(11):1121-1127.

[9] PROUGHBRED C J,HUNTER J O,DARBY A C,et al.Characterisation of the faecal metabolome and microbiome of thoroughbred racehorses[J].Equine Veterinary Journal,2015,47(5):580-586.

[10] YATSUNENK O T,REY F E,MANARY M J,et al.Human gut microbiome viewed across age and geography[J].Nature,2012,486(7402):222-227.

[11] COSTA M C,STMPFLI H R,ALLEN-VERCOE E,et al.Changes in the equine fecal microbiota associated with the use of systemic antimicrobial drugs[J].BMC Veterinary Research,2015,11:19.

[12] FONTY G,CHAUCHEYRAS-DURAND F.Lescosystèmes digestifs[M].Pairs:Tech & Doc Publishing,2007.

[13] COMBES S,MICHELLAND R J,MONTEILS V,et al.Postnatal development of the rabbit caecal microbiota composition and activity[J].FEMS Microbiology Ecology,2011,77(3):680-689.

[14] JULLIAND V,DE VAUX A,VILLARD L,et al.Preliminary studies on the bacterial flora of faeces taken from foals,from birth to twelve weeks.Effect of the oral administration of a commercial colostrum replace[J].Pferdeheilkunde,1996,12(3):209-212.

[15] 許宇靜,張煜坤,沈雪梅,等.采用Illumina MiSeq測序技術分析斷奶幼兔盲腸微生物群落的多樣性[J].動物營養(yǎng)學報,2015,27(9):2793-2802.

[16] 白秀娟,劉誠剛,杜智恒,等.PCR-DGGE技術分析斷奶仔兔腸道微生物菌群結(jié)構(gòu)及多樣性[J].東北農(nóng)業(yè)大學學報,2012,43(9):64-69.

[17] 谷莉.不同鼠種的腸道菌群在不同飲食結(jié)構(gòu)干預中的組成改變[D].博士學位論文.長沙:中南大學,2014.

[18] LEY R E,LOZUPONE C A,HAMADY M,et al.Worlds within worlds:evolution of the vertebrate gut microbiota[J].Nature Reviews Microbiology,2008,6(10):776-788.

[19] CALLAWAY T R,DOWD S E,EDRINGTON T S,et al.Evaluation of bacterial diversity in the rumen and feces of cattle fed different levels of dried distillers grains plus solubles using bacterial tag-encoded FLX amplicon pyrosequencing[J].Journal of Animal Science,2010,88(12):3977-3983.

[20] DE MENEZES A B,LEWIS E,O’DONOVAN M,et al.Microbiome analysis of dairy cows fed pasture or total mixed ration diets[J].FEMS Microbiology Ecology,2011,78(2):256-265.

[21] BRULC J M,ANTONOPOULOS D A,MILLER M E B,et al.Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases[J].Proceedings of the National Academy of Sciences of the United States of America,2009,106(6):1948-1953.

[22] RAY A K,BAIRAGI A,SARKAR G K,et al.Optimization of fermentation conditions for cellulase production byBacillussubtilisCY5 andBacilluscirculansTP3 isolated from fish gut[J].Acta Ichthyologica et Piscatoria,2007,37(1):47-53.

[23] JIAO J Z,WANG P P,HE Z X,et al.Invitroevaluation on neutral detergent fiber and cellulose digestion by post-ruminal microorganisms in goats[J].Journal of the Science of Food and Agriculture,2014,94(9):1745-1752.

[24] LOOFT T,ALLEN H K,CANTAREL B L,et al.Bacteria,phages and pigs:the effects of in-feed antibiotics on the microbiome at different gut locations[J].The ISME Journal,2014,8(8):1566-1576.

[25] SPENCE C,WELLS W G,SMITH C J.Characterization of the primary starch utilization operon in the obligate anaerobeBacteroidesfragilis:regulation by carbon source and oxygen[J].Journal of Bacteriology,2006,188(13):4663-4672.

[26] LOOFT T,JOHNSON T A,ALLEN H K,et al.In-feed antibiotic effects on the swine intestinal microbiome[J].Proceedings of the National Academy of Sciences of the United States of America,2012,109(5):1691-1696.

[27] BIAN G R,MA L,SU Y,et al.The microbial community in the feces of the white rhinoceros (Ceratotheriumsimum) as determined by barcoded pyrosequencing analysis[J].PLoS One,2013,8(7):e70103.

[28] DALY K,PROUDMAN C J,DUNCAN S H,et al.Alterations in microbiota and fermentation products in equine large intestine in response to dietary variation and intestinal disease[J].British Journal of Nutrition,2012,107(7):989-995.

[29] SHINKAI T,UEKI T,KOBAYASHI Y.Detection and identification of rumen bacteria constituting a fibrolytic consortium dominated byFibrobactersuccinogenes[J].Animal Science Journal,2010,81(1):72-79.

[30] EVANS N J,BROWN J M,MURRAY R D,et al.Characterization of novel bovine gastrointestinal tractTreponemaisolatesand comparison with bovine digital dermatitis treponemes[J].Applied and Environmental Microbiology,2011,77(1):138-147.

*Corresponding author, professor, E-mail： yangkailun2002@aliyun.com

(責任編輯 菅景穎)

A Study on Intestinal Microbiota Diversity of 3- to 6-Month-OldYiliHorses

LI Xiaobin ZHAO Guodong LIU Zhen WU Tingting HE Zhourui DENG Haifeng YANG Kailun*
(XinjiangKeyLaboratoryofHerbivoreNutritionforMeat&MilkProduction,XingjiangAgriculturalUniversity,Urumqi830052,China)

The current study was conducted to evaluate the intestinal microbiota diversity of 3- to 6-month-oldYilihorses, to reveal the regulation of succession of intestinal micrbiota before weaning, and to provide a theoretical basis for the researches of growth and intestinal health of young horse. Five 3-month-oldYilihorses with the average body weight of (89.75±8.81) kg and the same birth day were selected to conduct a 90-day feeding trial. The feces samples were collected on days 0 (1 day before the experiment), 30, 60 and 90 of the experiment, respectively. The microbial genomic total DNA was extracted from each feces sample. Microbiota diversity of feces was evaluated using Illumina HiSeq sequencing technology. The results showed as follows: 1) twenty feces samples for 5 colts were sequenced, a total of 157 665 effective sequences and 1 117 average operational taxonomic units (OTUs) were acquired. 2) The alpha diversity indexes (ACE, Chao1, Shannon and Simpson indices) of feces microbiota showed fluctuated changes with the age increasing, but the changes were not significant among different time points (P>0.05). 3) At the phylum level, ten dominant bacteria in feces of colts were Firmicutes, Bacteroidetes, Proteobacteria, Verrucomicrobia, Spirochaetes, Fusobacteria, Tenericutes, Actinobacteria, TM7 and Euryarchaeota, and the abundances of Firmicutes, Bacteroidetes, Proteobacteria and Verrucomicrobia were higher. At the family level, ten dominant bacteria in feces of colts were Bacillaceae, Moraxellaceae, Planococcaceae, Carnobacteriaceae, BS11, RFP12, Lactobacillaceae, Ruminococcaceae, Lachnospiraceae and Porphyromonadaceae. At the genera level, ten dominant bacteria in feces of colts wereAcinetobacter,Desemzia,Lactobacillus,Ureibacillus,Paludibacter,Bacillus,Escherichia,Carnobacterium,TreponemaandMogibacterium. Base on the results, it is concluded that the intestinal microbiota of 3- to 6-month-oldYilihorses can be accurately classified using Illumina HiSeq sequencing technology. The Firmicutes, Bacteroidetes, Proteobacteria, Verrucomicrobia were the dominant bacteria in intestine of 3- to 6-month-oldYilihorses.[ChineseJournalofAnimalNutrition, 2017, 29(5)：1535-1544]

Yilihorses; Illumina HiSeq sequencing technology; intestinal microbiota; diversity

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

2016-11-11

“十二五”農(nóng)村領域國家科技計劃課題(2012BAD45B01)

李曉斌(1988—)，男，甘肅天水人，博士研究生，研究方向為草食動物營養(yǎng)代謝。E-mail： 172387243@qq.com

*通信作者:楊開倫，教授，博士生導師，E-mail： yangkailun2002@aliyun.com

S852.6

A

1006-267X(2017)05-1535-10