飼料磷真消化率不同的山羊空腸菌群結(jié)構(gòu)差異研究

金 磊，王立志，王之盛，薛 白，彭全輝

(四川農(nóng)業(yè)大學 動物營養(yǎng)研究所，四川 成都 611130)

磷是動物體內(nèi)一種重要的常量礦物元素[1-3]，飼料中未被動物胃腸道消化吸收的磷，通過糞便排放到體外后，會造成嚴重的環(huán)境污染[4-5]。因此，提高動物對磷的消化利用效率對于降低動物養(yǎng)殖成本和減少環(huán)境污染均具有極其重要的意義。反芻動物飼料中的磷主要在小腸中進行消化吸收[6-7]，小腸微生物可能是影響其對飼料磷消化吸收的一個重要因素[8-10]。已有關(guān)于山羊胃腸道微生物多樣性的報道[11-12]，但磷消化率不同的山羊小腸微生物的結(jié)構(gòu)與組成是否具有差異仍未知。明晰飼料中磷消化率是否受到小腸微生物多樣性的影響，將有助于提高飼料磷利用率相關(guān)技術(shù)的研發(fā)。另外，由于消化道內(nèi)源磷的存在，采用表觀消化率并不能準確衡量動物對磷的消化吸收情況[13-14]。與表觀消化率相比，真消化率才能真實且準確地反映出宿主對磷的生物學利用效率[15-17]。本研究采用16S rRNA高通量測序技術(shù)比較飼料磷真消化率不同的山羊空腸微生物組成的差異性，旨在促進人們更清晰地認識空腸微生物與宿主飼料磷真消化率之間的關(guān)系，增進對空腸微生物與宿主磷消化利用關(guān)系的理解。

1 材料與方法

1.1 試驗動物與飼養(yǎng)管理

選取24只10月齡，雌性，平均體質(zhì)量(24.25±2.47)kg的努比亞黑山羊，試驗在四川農(nóng)業(yè)大學動物營養(yǎng)研究所實驗基地進行，整個試驗期試羊單籠飼養(yǎng)，每日飼料(干物質(zhì))按體質(zhì)量3.5%供給，自由飲水。參照我國《肉羊飼養(yǎng)標準》(NY/T 816—2004)，以每天每頭增重0.1 kg為標準配制日糧，日糧配方和營養(yǎng)水平詳見表1。

1.2 試驗設(shè)計與樣品采集

整個試驗由2期(Ⅰ期和Ⅱ期)構(gòu)成，每期包括14 d預飼期和6 d代謝試驗期，共40 d。代謝試驗期內(nèi)，每天記錄每只山羊的飼料供給量和殘留量，并采集當天飼喂的飼料及剩料樣品，于-20 ℃保存?zhèn)溆谩２捎萌占S法收集代謝試驗期內(nèi)山羊糞便，取糞總量的10%，用10%鹽酸固氮后于-20 ℃貯存?zhèn)溆谩＂衿诤廷蚱诜謩e飼喂對應日糧，試驗第41天晨飼前，屠宰所有山羊，迅速分割出空腸并進行結(jié)扎，以防止相鄰腸斷的內(nèi)容物流入空腸。無菌采集空腸內(nèi)容物樣品于凍存管中，-80 ℃保存?zhèn)溆谩?/p>

表1 日糧組成及營養(yǎng)水平(干物質(zhì)基礎(chǔ))

Table 1 The composition and nutrient level of the diet (dry matter basis) %

1) 預混料是由同一類的多種添加劑或不同類的多種添加劑按一定配比制作而成的勻質(zhì)混和物，預混料為每kg日糧提供：Fe(FeSO4)30 mg，Cu(CuSO4)10 mg，Zn(ZnSO4)50 mg，Mn(MnSO4)60 mg，VA2 937 IU，VD343 IU，VE30 IU； 2) 營養(yǎng)水平是指日糧中營養(yǎng)成分的含量，代謝能為計算值，其余為實測值。

1) Premix was a homogeneous mixture made from a variety of additives of the same class or a variety of additives of different kinds according to a certain proportion，premix provided the following per kg of the diet: Fe(FeSO4)30 mg， Cu(CuSO4)10 mg， Zn(ZnSO4)50 mg， Mn(MnSO4)60 mg， VA2 937 IU， VD343 IU， VE30 IU; 2) Nutrient level referred to the content of nutrients in the diet， metabolic energy was a calculated value and others were measured values.

根據(jù)山羊飼料采食量，飼料和糞便樣品中的磷含量等數(shù)據(jù)，計算出24個山羊個體有機磷的真消化率：

磷真消化率(%)=[(Ⅱ期食入磷-Ⅰ期食入磷)-(Ⅱ期排泄磷-Ⅰ期排泄磷)]÷(Ⅱ期食入磷-Ⅰ期食入磷)×100%[18]。

統(tǒng)計24只山羊?qū)︼暳狭渍嫦实钠骄导皹藴什睿瑢⒘渍嫦蚀笥?群體平均值+0.5倍標準差)的山羊作為高磷真消化率組(high true digestibility of phosphorus，HP)，磷真消化率小于(群體平均值-0.5倍標準差)的山羊作為低磷真消化率組(low true digestibility of phosphorus，LP)。

1.3 DNA提取和PCR擴增

提取HP組和LP組空腸內(nèi)容物微生物總DNA，DNA的提取和純化參照文獻[19]方法進行[19]。用細菌通用引物對：515F(5′-GTGCCAGCMGCCGCGGTAA-3′)和806R(5′-GGACTACVSGGGTATCTAAT-3′)，以空腸內(nèi)容物菌群總DNA為模板，針對細菌16S rRNA基因的V4區(qū)用PCR儀進行擴增。PCR反應體系為：5 U·μL-1Taq酶0.25 μL，10×Buffer 5.0 μL，10 mmol·L-1dNTPs 1.0 μL，10 μmol·L-1引物各1.25 μL，50 ng·μL-1DNA模板1.0 μL，補ddH2O至50 μL。擴增反應條件為：95 ℃ 2 min；95 ℃ 30 s，55 ℃ 30 s，72 ℃ 30 s，循環(huán)30次；72 ℃ 5 min。PCR產(chǎn)物用含溴化乙錠的2%瓊脂糖凝膠電泳鑒定并用PCR純化試劑盒純化和回收。回收的PCR產(chǎn)物用QuantiFluorTM-ST fluorometer定量，送北京諾禾致源生物技術(shù)有限公司構(gòu)建文庫，利用Illumina MiSeq測序平臺進行高通量測序。

1.4 生物信息學分析

測序平臺得到的原始數(shù)據(jù)經(jīng)過質(zhì)控、拼接、Tags過濾和去除剔除嵌合體[20-23]后，使用Uclust[24]聚類算法將序列聚類為OTU(operational taxonomic units)，并使用RDF分類器[25]對OTU代表序列從門到屬進行物種注釋。基于OTU table和rep_set. tree文件及抽樣的最大深度，繪制OTU稀釋曲線并計算α多樣性指數(shù)。結(jié)合各個樣品的OTU種類及其豐度進行計算，獲得樣品間Unweighted unifrac距離矩陣，并利用加權(quán)組平均法(UPGMA)進行聚類分析，然后繪制PCoA聚類圖[26-27]。根據(jù)門和屬水平上物種的結(jié)構(gòu)和組成，用OriginPro 9.0和R x64 3.0.2軟件繪制優(yōu)勢菌門柱形圖和共享屬聚類熱圖。基于16s rRNA擴增子測序結(jié)果，用PICRUSt(phylogenetic investigation of communities by reconstruction of unobserved states)，通過比對KEGG(kyoto encyclopedia of genes and genomes， http://www.genome.jpkegg)數(shù)據(jù)庫，對空腸微生物的基因進行功能預測，并挑選主要功能通路進行組間差異性分析[28-29]。

1.5 數(shù)據(jù)分析

采用SPSS 18.0統(tǒng)計軟件，對LP組和HP組之間細菌的相對豐度進行非參數(shù)檢驗，將P<0.05視為差異顯著，P<0.01為差異極顯著。對組間差異顯著的微生物，用其相對豐度與宿主對飼料磷的真消化率進行相關(guān)性(Spearman)分析。

2 結(jié)果與分析

2.1 山羊有機磷真消化率

如圖1所示，24只山羊?qū)θ占Z有機磷真消化率的平均值為(78.61±7.65)%，變異系數(shù)為9.73%。HP和LP組各挑選6只山羊，飼料有機磷的真消化率平均值分別為(89.20±0.01)%和(68.95±0.05)%，HP組山羊飼料有機磷真消化率極顯著高于LP組(P<0.01)。

2.2 空腸內(nèi)容物16S rRNA基因的高通量測序和OUT統(tǒng)計

對12個樣本進行測序，共獲得900 688條有效序列。基于相似性大于97%的原則，將獲得的有效序列進行聚類，共獲得5 611個OTU。其中HP組有3 830個OTU，平均每個樣品含(2 313±179)個OTU；LP組有3 977個OTU，平均每個樣品含(2 085±271)個OTU；兩組間共享的OTU有2 196個(圖2)。

圖1 磷真消化率Fig.1 True digestibility of phosphorus

各樣品稀釋曲線如圖3所示。在本試驗的測序深度下，各條曲線最終均趨于平緩，說明試驗的測序量足以覆蓋各樣品的絕大多數(shù)微生物。

2.3 山羊空腸微生物的α多樣性分析

α多樣性指數(shù)主要用于評價樣品中微生物的豐富性和均勻性。表2列出了HP組和LP組樣品的Shannon、Chao1、Goods coverage和PD_whole_tree指數(shù)。在同一測序深度下，HP組的α多樣性指數(shù)均大于LP組，但僅Chao1指數(shù)在組間具有顯著差異(P<0.05)。

圖2 OTU維恩圖Fig.2 OTU venn diagram

圖3 各樣品的稀釋曲線Fig.3 Rarefaction curves of each sample

2.4 山羊空腸微生物門水平和屬水平上的微生物組成

將所得有效序列在不同分類水平上進行物種注釋和統(tǒng)計，將相對豐度大于0.1%的微生物視為優(yōu)勢菌群。HP和LP組均有8個優(yōu)勢菌門，2個組相對豐度排名前10的門均依次分別為厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidetes)、變形菌門(Proteobacteria)、黏膠球形菌門(Lentisphaerae)、螺旋體門(Spirochaetes)、軟壁菌門(Tenericutes)、放線菌門(Actinobacteria)、螺旋體菌門(Saccharibacteria)、疣微菌門(Verrucomicrobia)和綠彎菌門(Chloroflexi)(圖4)。厚壁菌門(Firmicutes)在HP和LP組的相對豐度分別達到了59.75%和70.69%，處于絕對優(yōu)勢地位。對HP和LP組的優(yōu)勢菌門進行組間差異分析(表3)，HP組厚壁菌門(Firmicutes)的相對豐度極顯著低于LP組(P<0.01)，擬桿菌門(Bacteroidetes)極顯著高于LP組(P<0.01)。

在屬水平上，2個組相對豐度大于0.1%的微生物屬共有47個(圖5)。對HP和LP組的優(yōu)勢菌屬進行組間差異性分析，其中有顯著性差異的菌屬如表3所示，HP組Prevotella、Rikenellaceae_RC9_gut_group、Fibrobacter、Ruminococcus_2和Clostridium_sensu_stricto_1的相對豐度顯著高于LP組(P<0.05)，而Desulfovibrio和Victivallis的相對豐度顯著低于LP組(P<0.05)。

如表3所示，在門水平，宿主的磷真消化率與擬桿菌門(Bacteroidetes)的相對豐度極顯著正相關(guān)(P<0.01)，與厚壁菌門(Firmicutes)的相對豐度極顯著負相關(guān)(P<0.01)。在屬水平，山羊磷真消化率與Prevotella、Fibrobacter和Victivallis的相對豐度顯著正相關(guān)(P<0.05)，與Desulfovibrio和Ruminococcus_2的相對豐度極顯著正相關(guān)(P<0.01)，與Clostridium_sensu_stricto_1的相對豐度極顯著負相關(guān)(P<0.01)。

表2 組間α多樣性指數(shù)的差異性比較

Table 2 Comparison of α diversity index between the two groups

組別GroupShannonChao1Goods coveragePD_whole_treeHP7.38±0.852170.24±202.840.99±0.01107.26±7.30LP6.60±1.521779.45±304.720.99±0.0193.31±15.13P值P value0.3020.0260.1610.07

圖4 門水平上的菌群組成Fig.4 Composition of microorganism at phylum level

圖5 共享菌群在屬水平上的聚類熱圖Fig.5 Log-scaled percentage heatmap of bacteria community structure at the shared-genus level

表3 HP和LP組在門和屬水平上差異的菌群及其與磷真消化率的相關(guān)性

Table 3 The difference between HP and LP groups at the level of portal and genera and their correlation with true phosphorus digestibility

分類單元Taxa相關(guān)豐度 Relative abundanceHPLPP值 P value相關(guān)系數(shù)rCorrelation coefficient門Phylum厚壁菌門 Firmicutes59.753±4.76770.692±6.3270.007-0.813**擬桿菌門 Bacteroidetes10.032±2.9165.704±1.4310.0080.718**屬GenusPrevotella5.341±2.6631.922±0.9110.0140.686*Rikenellaceae_RC9_gut_group 2.640±0.2961.580±0.2500.0010.401Fibrobacter0.825±0.2600.376±0.1730.0170.623*Desulfovibrio0.134±0.0190.358±0.0320.0030.845**Victivallis3.555±2.3330.724±0.1990.0040.695*Ruminococcus_23.393±0.2762.468±0.4310.0060.843**Clostridium_sensu_stricto_1 0.851±0.1195.083±0.5590.002-0.995**

**表示極顯著相關(guān)(P<0.01)，*表示顯著相關(guān)(P<0.05)。

**Represented extremely significant correlation(P<0.01)， * Represented significant correlation(P<0.05).

2.5 PICRUSt基因預測

基于本試驗16s rRNA擴增子的測序結(jié)果，用PICRUSt對空腸微生物的基因進行功能預測，2個組共預測到了35個功能基因。圖6展示了相對豐度處于前20位的功能基因，這20種功能基因分別占HP組和LP組總基因豐度的93.49%和94.08%。圖6由下到上基因的相對豐度逐步降低，相對豐度最高的基因其功能依次為膜轉(zhuǎn)運、碳水化合物代謝、氨基酸代謝、能量代謝、輔助因子和維生素代謝、脂質(zhì)代謝等。這20種功能基因，其相對豐度在組間的差異均不顯著(P>0.05)。

3 討論

本試驗采用限飼，磷的供給量低于肉羊飼養(yǎng)標準，滿足采用差量法測定飼料磷真消化率的條件。本研究測定的山羊?qū)︼暳狭椎恼嫦势骄禐?8.61%，略高于單胃動物豬[30]和雞[31]。分析原因，可能是由于反芻動物胃腸道具有更豐富的胃腸道微生物，微生物可能會顯著影響宿主對飼料中磷的消化吸收。動物對日糧中磷消化的影響因素主要包括日糧磷的來源、磷的水平、鈣磷比和日糧植酸的水平等[32-35]。一般日糧中磷含量越高，無機磷的含量越低，動物對磷的消化利用效率越低[36]。但本試驗結(jié)果發(fā)現(xiàn)，飼喂相同飼糧且遺傳背景相同的山羊，飼料磷的真消化率也存在顯著的個體差異。

圖6 20種優(yōu)勢功能基因的相對豐度組成Fig.6 The composition of the relative abundance contributed by the top 20 genes

學者們普遍認為小腸是磷消化吸收的主要場所，但有關(guān)宿主對飼料磷的消化吸收與小腸中微生物的相關(guān)性研究還未見報道。本研究采用高通量測序技術(shù)探索宿主磷真消化率與空腸微生物組成關(guān)系，較全面地覆蓋了研究對象的空腸微生物，能真實全面地反映其微生物組成情況。研究結(jié)果發(fā)現(xiàn)，HP組α多樣性指數(shù)中的Chao1指數(shù)顯著高于LP組，暗示HP組微生物的豐富度或者是均勻度高于LP組；但是HP組和LP組之間OTU數(shù)并無顯著差異，表明2組的微生物豐富度并無顯著性差異。綜合OTU數(shù)和Chao1指數(shù)的分析結(jié)果可以得出，HP組在空腸微生物的均勻度上要顯著高于LP組，PCoA聚類圖也直觀地證明了這一點。

在本研究中，HP組和LP組空腸微生物相對豐度最高的門均是厚壁菌門(Firmicutes)，肉牛和奶牛上的研究都表明[37-38]，厚壁菌門(Firmicutes)在空腸微生物中處于絕對優(yōu)勢地位，這可能與反芻動物日糧中粗飼料含量較高有關(guān)，有研究表明，胃腸道中厚壁菌門(Firmicutes)的比例與日糧中粗飼料的含量呈顯著正相關(guān)[39]。在門和屬水平上，HP和LP組之間均存在差異顯著的菌群，并且這些差異菌群與宿主對飼料磷的真消化率之間存在顯著的相關(guān)性，說明空腸微生物的結(jié)構(gòu)會顯著影響山羊?qū)︼暳狭椎南铡Ｔ诒狙芯恐校琀P組Bacteroidetes屬的相對豐度極顯著高于LP組，同時隸屬于Bacteroidetes門的Prevotella屬和Rikenellaceae_RC9_gut_group屬的相對豐度也顯著高于LP組，并且三者的相對豐度均與磷真消化率之間呈顯著的正相關(guān)。這些結(jié)果表明，空腸中的Bacteroidetes門菌群可能有助于宿主對飼料磷的消化吸收。尤其值得注意的是Bacteroidetes門的Prevotella屬，其所能消化分解的底物非常廣泛，例如隸屬于Prevotella屬的短普雷沃氏菌，不但具有很高的蛋白酶活性，還能發(fā)酵葡萄糖、乳糖和纖維二糖；棲瘤胃普雷沃氏菌和布氏普雷沃氏菌，能分泌木聚糖酶和羧甲基纖維素酶以降解淀粉、木聚糖等營養(yǎng)物質(zhì)[40]。因此，普雷沃氏菌屬幾乎可以單獨完成蛋白質(zhì)、糖類和碳水化合物等營養(yǎng)物質(zhì)的降解任務(wù)。本研究中，Prevotella屬菌群可能也參與了飼料中植酸磷的分解，使得植物性飼料中所含的原本不能被宿主利用的磷被釋放出來，從而被山羊消化吸收。已有研究表明，Prevotella屬菌群能分泌微生物植酸酶催化植酸的水解，釋放飼料中的植酸磷，提高宿主對飼料植物磷的消化率[41]。

以有研究表明，F(xiàn)ibrobacter和Ruminococcus_2均是動物胃腸道中重要的纖維分解菌。Fibrobacter能產(chǎn)生多糖酶，這些酶能夠降解日糧中的多種纖維類物質(zhì)[42-44]。隸屬于Ruminococcus屬的白色瘤胃球菌和黃色瘤胃球菌都是主要的纖維分解菌種，在纖維素和半纖維素的降解中發(fā)揮著重要作用。HP組的Fibrobacter和Ruminococcus_2的相對豐度顯著高于LP組，且相對豐度都與山羊?qū)︼暳狭椎南曙@著正相關(guān)。產(chǎn)生這種現(xiàn)象的原因可能是飼料纖維結(jié)構(gòu)的破壞促進了飼料中被束縛的磷的釋放，纖維物質(zhì)的降解發(fā)酵可以提高飼料磷的消化吸收率[45]。

目前為止，人們對腸道微生物在宿主營養(yǎng)物質(zhì)消化代謝的過程中發(fā)揮的作用認知十分有限。通過PICRUSt基因預測，空腸微生物最主要的功能類別為代謝功能類，如能量代謝、糖酵解途徑、碳水化合物代謝、脂肪代謝和氨基酸代謝等，這與前人的研究結(jié)果一致[46-47]。必須指出的是，雖然本研究表明，空腸微生物與宿主對飼料磷的真消化率存在顯著的相關(guān)性，但目前尚缺乏直接的證據(jù)證明微生物的差異直接影響了山羊?qū)︼暳狭椎南铡５豢煞裾J的是，本研究首次發(fā)現(xiàn)了空腸微生物的組成與山羊磷的真消化率存在關(guān)聯(lián)性，這為今后研究調(diào)控微生物以促進磷的消化利用提供了嶄新的思路。