萊蕪豬和杜長大豬盲腸和結腸微生物菌群結構組成和功能分析

李 超,趙雪艷,王永軍,王彥平,任一帆,李菁璇,王懷中,王繼英*,宋勤葉

(1.河北農業大學動物醫學院,保定 071000;2.山東省農業科學院畜牧獸醫研究所山東省畜禽疫病防治與繁育重點實驗室,濟南 250100;3.農業農村部畜禽生物組學重點實驗室,濟南 250100;4.濟南市動物疫病預防與控制中心,濟南 250000)

動物腸道微生物是一個非常復雜的生態系統,影響許多重要的宿主生理功能,如食物攝入、新陳代謝調節、免疫系統激活、上皮細胞增殖和抗感染[1]。大量研究表明,不同品種豬具有不同的腸道菌群結構,與引進豬種相比,地方豬種腸道微生物組成更豐富[2]。近年來,研究發現腸道微生物菌群對豬脂肪沉積、飼料利用率、免疫抗性等有重要影響。He等[3]、Yang等[4]和Chen等[5]通過分析腸道微生物菌群與背膘厚、脂肪率等性狀的關系,鑒定出普雷沃氏菌屬(Prevotella)、大腸桿菌(Escherichiacoli)等影響脂肪沉積的微生物類群。Bergamaschi 等[6]研究發現,直腸微生物菌群中的厭氧弧菌屬(Anaerovibrio)、擬桿菌屬(Bacteroides)、布氏桿菌屬(Blautia)等影響杜洛克、大約克和長白等瘦肉型豬的飼料利用率。Zhao等[7]通過研究宿主-微生物群相互作用,揭示乳桿菌屬(Lactobacillus)、梭菌屬(Clostridia)和真桿菌屬(Eubacterium)等潛在有益菌介導了民豬對炎癥性腸病的抗性。

萊蕪豬是中國著名的地方豬品種,具有肉質好、抗逆性強、耐粗飼等優良種質特性[8],尤其是肌內脂肪沉積能力特別強,含量高達10%[9],是我國地方豬種資源的寶貴基因庫。系統研究萊蕪豬和杜長大(DLY)豬盲腸和結腸微生物菌群結構組成和功能,鑒定影響萊蕪豬特性的微生物菌屬,將為從腸道菌群層面深入理解萊蕪豬特性奠定基礎。本研究以DLY豬為對比,系統研究自然生產狀態下萊蕪豬盲腸和結腸腸道食糜微生物群落的結構組成和功能,分析腸道菌群與背膘厚、血清生化指標的相關性,研究腸道中優勢菌種與其表型性狀的潛在聯系,為從腸道菌群層面深入理解萊蕪豬高脂肪沉積、耐粗飼和抗逆等特性作用機理奠定理論基礎。

1 材料與方法

1.1 試驗動物和樣品采集

本研究中的12頭萊蕪豬去勢公豬來自萊蕪豬種豬繁育有限公司,6頭DLY豬去勢公豬來自萊蕪豬種豬繁育有限公司臨近的商品豬場。所有試驗豬健康無腹瀉,且屠宰前3個月內均未使用過抗生素。萊蕪豬和DLY豬分別參照國家標準《豬的營養需要量》(GB/T 39235—2020)的脂肪型豬和瘦肉型豬的標準配制。飼糧組成及營養水平詳見表1。

表1 飼糧組成及營養水平(風干基礎)Table 1 Diet composition and nutrient level in this study(air-dry basis)

萊蕪豬日齡303 d,宰前體重為(95±5.8)kg,杜長大豬日齡192 d,宰前體重為(120±4.8)kg,均處于肥育后期。12頭試驗動物禁食24 h,電擊致暈后按照標準流程屠宰。屠宰時采血,3 000 r·min-1離心5 min,吸取血清于離心管中,置于-20 ℃冰箱中保存。收集盲腸和結腸食糜于無菌離心管中,液氮中速凍,置于-80 ℃冰箱中保存。

1.2 背膘厚和血清生化指標測定

根據《瘦肉型豬胴體性狀測定技術規范》( NY/T 825—2004)[10]測量胴體背中線肩部最厚處、最后肋、腰薦結合處三點的皮下脂肪厚度,取3點平均值作為平均背膘厚。采用日立7100自動生化分析儀測定血清中的總蛋白(TP)、甘油三酯(TG)、葡萄糖(GLU)、高密度脂蛋白膽固醇(HDL-C)、低密度脂蛋白膽固醇(LDL-C)、尿素氮(BUN)、總膽固醇(TCH)和二胺氧化酶(DAO)。

1.3 DNA 提取、PCR擴增和16S rDNA高通量測序

使用土壤和糞便基因組DNA提取試劑盒(Qiagen, Valencia, California, USA)提取盲腸和結腸食糜中總基因DNA。使用NanoDrop ND-1000分光光度計(NanoDrop,德國)測量提取的DNA的濃度和純度,并在1% (w/v)瓊脂糖凝膠上評估完整性。質控合格后,使用特異性引物341F(CCTAYGGGRBGCASCAG)和806R(GGACTACNNGGGTATCTAAT)擴增16S rRNA基因的V3-V4區域[11]。根據PCR產物濃度等量混樣,充分混勻后使用2%濃度的瓊脂糖凝膠電泳分離PCR產物,并使用Universal DNA Gel Extraction Kit (TianGen, 北京, 中國)純化目的條帶。使用NEB-Next?Ultra DNA Library Prep Kit (Illumina, San Diego, CA, USA)構建測序文庫。文庫用Qubit 2.0熒光計(Thermo Scientific)量化后,由諾禾致源(中國北京)利用Illumina NovaSeq 6000進行雙末端測序,測序讀段長度為250 bp。

1.4 16S rDNA測序數據處理

根據Barcode和PCR擴增引物序列從下機數據中拆分出各樣本數據,截去Barcode和引物序列后使用FLASH(V1.2.7)[12]對reads進行拼接,得到原始Tags數據(Raw Tags);原始Tag使用QIIME(V1.9.1)[13]進行如下質量控制得到高質量的Tags(Clean tags):a)Tag截取:將原始Tags從連續低質量值(默認質量閾值為≤19)堿基數達到設定長度(默認長度值為3)的第一個低質量堿基位點截斷;b)Tag長度過濾:經過截取后得到的Tags數據集,進一步過濾掉其中連續高質量堿基長度小于Tags長度75%的Tags。經過以上處理后得到的高質量Tag序列通過VSEARCH[14]與物種注釋數據庫比對檢測和去除嵌合體序列,并得到最終的有效Tag數據(Effective Tags)。

利用Uparse算法(Uparse v7.0.1001)[15]對樣本的全部有效Tags進行聚類,以默認的97%序列一致性(Identity)將有效Tag序列聚類成為分類操作單元(OTUs),同時選取出現頻數最高的序列作為OTUs的代表序列。使用Mothur方法與SILVA138[16]的SSUrRNA數據庫[17]對OTUs代表序列進行物種注釋分析(設定閾值為0.8～1),獲得分類學信息并分別在各個分類水平(門、綱、目、科、屬、種)統計各樣本的微生物菌群組成。使用MUSCLE[18](Version 3.8.31)軟件進行快速多序列比對,得到所有OTUs代表序列的系統發生關系。最后以樣本中數據量最少的為標準進行均一化處理,后續的Alpha多樣性分析和Beta多樣性分析都是基于均一化處理后的數據。

1.5 腸道菌群多樣性及品種差異分析

使用QIIME(V 1.9.1)計算樣品的Simpson和Shannon等Alpha多樣性指標,并使用Wilcox檢驗分析Alpha多樣性指數組間差異;使用QIIME軟件(V 1.9.1)計算 Bray-Curtis相異度以估計群落結構的差異,使用主坐標分析(PCoA)對樣本的相似性或差異性進行可視化呈現。

使用LEfSe[19]開展品種間物種差異判別分析以鑒定顯著性影響萊蕪豬和DLY豬的腸道微生物菌群,所用分析參數為:使用Kruskal-Wallis檢驗比較OTU的相對豐度,P水平為0.05;將所得顯著不同的OTU用作線性判別分析(LDA)的輸入,以計算它們的效應大小,選擇對數LDA得分閾值>3.5的OTU作為組間顯著不同的OTU。使用Metastats[20]對門、屬、種水平的各個分類單元在萊蕪豬和DLY豬之間的豐度差異進行比較檢驗。

1.6 腸道菌群與背膘厚、血清生化指標的相關性分析

使用R語言計算背膘厚和血清生化指標與豐度前40的腸道微生物菌屬豐度之間的Spearman相關系數,將獲得的相關系數數值矩陣通過R語言pheatmap程序包繪制相關性熱圖。

1.7 腸道菌群功能預測

使用Tax4Fun[21]對16S rDNA基因測序數據進行功能預測,并使用熱圖可視化前35個KEGG通路的統計差異。

2 結 果

2.1 背膘厚和血清生化指標檢測

萊蕪豬和DLY豬的背膘厚和血清生化指標檢測結果見表2。可以看出,萊蕪豬和DLY豬具有不同的品種特點,二者的背膘厚、部分血清生化指標存在顯著差異(P<0.05)。萊蕪豬的背膘厚和TP顯著高于DLY豬(P<0.05),GLU、TG和HDL-C均顯著低于DLY豬(P<0.05)。

表2 萊蕪豬和杜長大豬背膘厚和血清生化指標比較Table 2 Comparison of backfat thickness and serum biochemical parameters between Laiwu and DLY pigs

2.2 測序數據概述

萊蕪豬和DLY豬的盲腸、結腸腸道微生物菌群16S rRNA 測序數據如表3所示。本研究共得到2 105 633個高質量的有效tag序列,平均每個樣本為 58 490個,有效tag序列的平均長度為409 bp,Q20和Q30的平均值分別為98.39% 和94.79%,表明本研究測序數據質量良好。

表3 16S rDNA測序數據統計性描述Table 3 Summary statistics for the 16S rDNA sequencing data

2.3 腸道微生物菌群多樣性分析

基于97%的序列一致性將有效Tag序列聚類,共獲得了6 642個OTUs。萊蕪豬和DLY豬盲腸和結腸共有和特異的OTU韋恩圖(圖1)顯示,萊蕪豬、DLY豬盲腸和結腸4個腸段中共有的OTU數量為1 130個,占所有OTU總數的17.01%;萊蕪豬盲腸和結腸共有OTU數 1 722個,DLY豬盲腸和結腸的共有OTU 1 251個;萊蕪豬和DLY豬盲腸和結腸共有OTU 分別為1 454和 1 372個。特異OTU分析顯示,萊蕪豬盲腸和結腸特異OTU分別為708和2 825個,DLY豬盲腸和結腸特異OTU分別為690和283個。萊蕪豬盲腸和DLY豬盲腸的特異OTU數量基本一致,而萊蕪豬結腸的特異OTU數量遠高于DLY豬結腸的特異OTU數量。

DLYCA代表杜長大豬盲腸,DLYCO代表杜長大豬結腸,LWCA代表萊蕪豬盲腸,LWCO代表萊蕪豬結腸,下同DLYCA means caecum of DLY pig, DLYCO means colon of DLY pig, LWCA means caecum of Laiwu pig, LWCO means colon of Laiwu pig, the same as below圖1 萊蕪豬和杜長大豬盲腸和結腸共有和特異分類操作單元(OTUs)韋恩圖Fig.1 Venn diagram of common and unique operational taxonomic units (OTUs) of caecum and colon in Laiwu and DLY pigs

萊蕪豬和DLY豬盲腸和結腸微生物菌群的Alpha多樣性見圖2所示,可以看出利用Simpson和Shannon指數分析的Alpha多樣性結果基本一致。萊蕪豬盲腸Simpson指數高于DLY豬,且萊蕪豬與DLY豬在盲腸的Simpson指數上的差異達到了顯著水平(P<0.05),表明萊蕪豬比DLY豬具有更豐富的盲腸微生物多樣性。

左為Simpson指數,右為Shannon指數,*表示P<0.05Left. Simpson index; Right. Shannon index, *. P<0.05圖2 萊蕪豬和杜長大豬盲腸和結腸微生物的Alpha多樣性比較Fig.2 Comparison of microbial alpha diversity of caecum and colon in Laiwu and DLY pigs

為衡量樣本或組間的微生物組差異,根據Bray-Curtis相異度繪制PCoA圖(圖3)。可以看到,萊蕪豬盲腸(LWCA)、萊蕪豬結腸(LWCO)、DLY豬盲腸(DLYCA)和DLY豬結腸(DLYCO)4組樣本相對聚集在一起,組間有一定程度重疊。DLY豬盲腸樣本分布最為集中,表明DLY豬樣本盲腸微生物一致性最高。DLY豬盲腸與結腸樣本只有少許重疊,而萊蕪豬盲腸與結腸樣本重疊較多,表明萊蕪豬盲腸與結腸微生物菌群的相似度更高。萊蕪豬盲腸樣本與DLY豬盲腸樣本、萊蕪豬結腸樣本與DLY豬結腸樣本均有部分重疊,表明不同豬種同一腸段的菌群之間有一定的相似性。

圖3 萊蕪豬和杜長大豬盲腸和結腸微生物菌群的主坐標分析Fig.3 Principal coordinate analysis by the bacterial microbiota of caecum and colon in Laiwu and DLY pigs

2.4 腸道微生物菌群分類學統計

本研究檢測的6 642個OTUs中,4 750個(占總數的 71.51%)可在Silva數據庫中注釋,門、綱、目、科、屬和種水平的注釋率分別為64.50%、62.15%、56.55%、46.69%、31.11%、6.68%。

在門水平上(圖 4A),萊蕪豬和DLY豬均以厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidota)為優勢菌門,厚壁菌門(Firmicutes)和擬桿菌門(Bacteroidota)豐度之和分別占萊蕪豬盲腸微生物的92.86%和結腸微生物的91.15%,占DLY豬盲腸微生物的96.13%和結腸微生物的96.51%。品種間比較,萊蕪豬盲腸和結腸中的厚壁菌門(Firmicutes)的相對豐度占比均低于DLY豬,其中盲腸的差別達到顯著水平(69.16%vs. 89.59%,P<0.05,圖5A)。此外,萊蕪豬結腸和盲腸中的廣古菌門(Euryarchaeota)和梭桿菌門(Fusobacteriota)相對豐度分別為2.47%和2.54%,遠高于DLY豬結腸和盲腸中的這2個門的相對豐度(0.06%和0.01%),但沒有達到顯著差異(P>0.05)。

在屬水平上(圖4B),萊蕪豬和DLY豬盲腸和結腸微生物豐度最高的菌屬均為土孢桿菌屬(Terrisporobacter),占比為15.64%～32.50%。在萊蕪豬和DLY豬盲腸中微生物豐度排名第2、3位的屬分別為蘇黎世桿菌屬(Turicibacter)和狹義梭菌屬(Clostridium_sensu_stricto_1),在萊蕪豬結腸中微生物豐度排名第2、3位的屬分別為狹義梭菌屬(Clostridium_sensu_stricto_1)和蘇黎世桿菌屬(Turicibacter),而在DLY豬結腸中微生物豐度排名第2、3位的屬分別為鏈球菌屬(Streptococcus)和狹義梭菌屬(Clostridium_sensu_stricto_1)。品種間比較,萊蕪豬盲腸中(Oscillospiraceae_UCG-007)和候選單胞生糖菌屬(Candidatus_Saccharimonas)的相對豐度顯著低于DLY豬(0.004 347 % 和0.019 69%vs.0.013 30% 和0.158%,P<0.05,圖5B、5C)。另外,萊蕪豬結腸中甲烷短桿菌(Methanobrevibacter)、活潑瘤胃球菌([Ruminococcus]_gnavus_group)和盲腸中的梭桿菌屬(Fusobacterium)的相對豐度遠高于DLY豬對應菌屬的相對豐度(2.47%、3.11%和2.54%vs.0.03%、0.22%和0.01%),萊蕪豬結腸中的鏈球菌屬(Streptococcus)和乳酸菌屬(Lactobacillus)相對豐度遠小于DLY豬對應菌屬的相對豐度(0.7%和2.21%vs.15.50%和15.37%),但差異均沒有達到顯著水平(P>0.05)。

圖4 萊蕪豬和杜長大豬盲腸和結腸微生物菌群前10個門(A)和前10個屬(B)相對豐度組成Fig.4 Relative abundance of the top 10 bacterial microbiota at phylum (A) and genus (B) levels of caecum and colon in Laiwu and DLY pigs

A. 厚壁菌門;B. Oscillospiraceae_UCG-007菌屬;C. 候選單胞生糖菌屬。* 表示差異顯著(P<0.05)A. Firmicutes; B. Oscillospiraceae_UCG-007; C. Candidatus_Saccharimonas. Significance is presented as * (P<0.05).圖5 萊蕪豬和杜長大豬盲腸微生物顯著差異菌門和屬Fig.5 Significantly different microbiota at phylum and genus levels of caecum in Laiwu and DLY pigs

此外,為了鑒定顯著影響萊蕪豬和DLY豬盲腸和結腸的特定微生物菌群,本研究進一步使用LEfSe開展組間物種差異判別分析(圖6)。從圖6中可以看出,萊蕪豬盲腸中檢測到的特異菌屬最多,包括變形菌門(Proteobacteria)中的大腸志賀桿菌屬(Escherichia_Shigella)、厚壁菌門(Firmicutes)中的克里斯滕森菌科_R-7菌屬(Christensenellaceae_R_7_group)和顫螺旋菌科(UCG-005屬)、梭桿菌門(Fusobacteria)中的梭菌屬(Fusobacterium)等特異菌屬。在DLY豬盲腸中檢測到的特異菌屬包括厚壁菌門(Firmicutes)中的蘇黎世桿菌屬(Turicibacter)和消化鏈球菌科(Peptostreptococcaceae)中的Intestinibacter、真桿菌屬成員(unidentified_Eubacterium_coprostanoligenes_group)等。在萊蕪豬結腸檢測到的特異菌屬包括厚壁菌門(Firmicutes)中的狹義梭菌屬(Clostridium_sensu_stricto_1)、放線菌門(Actinobacteriota)中的棒狀桿菌屬(Corynebacterium)等。在DLY豬結腸中檢測到的特異菌屬包括厚壁菌門(Firmicutes)中的乳酸菌屬(Lactbbacillus)、鏈球菌屬(Streptococcus)和(Subdoigranulum)等。

圖6 萊蕪豬和杜長大豬盲腸和結腸微生物群落LEfSe多級物種差異判別分析結果Fig.6 Bacterial taxa differentially represented in the corresponding segment of caecum and colon of Laiwu and DLY pigs identified by LEfSe

2.5 腸道菌群豐度與背膘厚和血清生化指標的相關性分析

為探討腸道菌群豐度與背膘厚和血清生化指標間的關系,本研究計算了腸道微生物菌群前40個菌屬的豐度與背膘厚、TP、GLU、TG和HDL-C在兩品種顯著差異指標間的Spearman相關系數,詳見圖7。由圖7可知,腸道微生物菌屬豐度與背膘厚和部分血清生化指標具有的顯著相關性(P<0.05),且盲腸和結腸的微生物菌屬與各性狀的相關模式不同。例如盲腸中背膘厚與擬桿菌屬(Bacteroides)、顫螺菌科成員(Oscillospiraceae-UCG-005)、擬梭菌屬(Clostridioides)、歐陸森氏菌屬(Olsenella)和小桿菌屬(Dialister)呈顯著正相關,與蘇黎世桿菌屬(Turicibacter)、狹義梭菌屬(Clostridium_sensu_stricto_1)和腸桿菌屬(Intestinibacter)呈顯著負相關;結腸中背膘厚與狹義梭菌屬(Clostridium_sensu_stricto_1)、甲烷短桿菌屬(Methanobrevibacte)、蘇黎世桿菌屬(Turicibacter)、土壤芽孢桿菌屬(Solibacillus)、擬梭菌屬(Clostridioides)、羅氏菌屬(Rothia)和不動桿菌屬(Acinetobacter)與呈顯著正相關,與鏈球菌屬(Streptococcus)、罕見小球菌屬(Subdoligranulum)和陰溝桿菌屬(Agathobacter)呈顯著負相關。

圖7 豬盲腸(左)和結腸(右)微生物菌群與背膘厚和血清生化指標相關性熱圖Fig.7 Correlation heatmap between caecum(left) and colon(right) bacterial microbiota and backfat thickness and serum biochemical parameters in Laiwu and DLY pig

2.6 腸道細菌群落功能預測分析

為探討腸道菌群影響宿主表型的機理,利用Tax4fun對萊蕪豬和DLY豬盲腸和結腸微生物菌群進行了功能預測分析(圖8)。可以看出,盲腸和結腸微生物菌群的功能主要富集在代謝、免疫、信息處理等通路。與結腸相比,盲腸富集了更多的通路。品種間比較,萊蕪豬和DLY豬微生物菌群功能富集的通路有很大不同。萊蕪豬盲腸和結腸富集了異種生物降解和代謝、輔助因子和維生素代謝、能量代謝、其他次生代謝產物的生物合成、氨基酸代謝、多酮類和萜類化合物的代謝等代謝相關通路;細胞過程和信號、細胞群落-原核生物、細胞生長和死亡、細胞運動等細胞過程相關通路;衰老、癌癥、內分泌和代謝疾病、免疫系統等生物體系統和疾病相關通路。DLY豬盲腸和結腸富集了其他氨基酸代謝、核苷酸代謝、碳水化合物代謝、脂類物質代謝、多糖生物合成與代謝等代謝相關通路;轉錄、復制和修復、翻譯、折疊、分類和降解等遺傳信息處理相關通路;抗藥性、傳染病等疾病相關通路。

圖8 基于Tax4fun在第二層級上的功能預測的熱圖聚類Fig.8 The heatmap of microbial functions at the second level predicted by Tax4fun

3 討 論

3.1 盲腸與結腸微生物菌群的多樣性和組成

目前,豬腸道微生物研究主要集中在對消化道的糞便或直腸微生物菌群研究[22-24]。然而,由于解剖結構和生理特性的差異,豬的腸道各部位菌群在種類、豐度上也明顯不同[25-26]。糞便只能定性地反映腸道的菌群種類,卻不能定量反映不同腸段的菌群豐度[27]。盲腸的結構有利于食糜停留,結腸中食糜的流通速度緩慢,豬的盲腸和結腸給予細菌大量增殖的機會,是豬主要的發酵器官。本研究中萊蕪豬盲腸和結腸的Shannon指數為5.52和5.31,低于譚占坤等[28]報道的放牧藏豬、舍飼藏豬和DLY豬盲腸的Shannon指數8.68、7.48和7.21以及Shang等[29]報道的藏豬和約克夏豬結腸的Shannon指數5.68和6.22。且本研究中萊蕪豬盲腸和結腸的微生物多樣性高于Yang等[4]報道的萊蕪豬的空腸和回腸的Shannon指數(1.8～2.4)以及Holman等[26]薈萃分析的十二指腸、空腸和回腸的Shannon指數(1.87、3.06和2.72),表明盲腸和結腸的微生物多樣性高于小腸腸段。

本研究中,厚壁菌門(Firmicutes)和擬桿菌門(Bacteroidota)是盲腸和結腸腸道微生物中豐度占比最高的菌門(圖4),厚壁菌門和擬桿菌門豐度之和分別占萊蕪豬盲腸和結腸微生物的92.86%和91.15%,占DLY豬盲腸和結腸微生物的96.13%和96.51%。這與Isaacson和Kim[30]研究中豬盲腸和結腸有超過90%的細菌歸屬于厚壁菌門和擬桿菌門的結果相一致。而且,本研究發現,與DLY豬相比較,萊蕪豬具有獨特腸道菌群組成。在門水平上,萊蕪豬盲腸和結腸中的厚壁菌門(Firmicutes)的相對豐度占比均低于DLY豬,其中盲腸的差別達到顯著水平(69.16%vs. 89.59%)。對人類和小鼠的研究表明,在肥胖個體糞便微生物菌群研究中,厚壁菌門(Firmicutes)的豐度占比例較大[31],而擬桿菌門(Bacteroidota)的比例較低[32],與本文中的研究不一致。這種差異可能是由于物種以及糞便微生物組和盲腸腸道微生物組成的差異所致。萊蕪豬盲腸和結腸微生物中占比例較高菌屬為土孢桿菌屬(Terrisporobacter)、蘇黎世桿菌屬(Turicibacter)和狹義梭菌屬(Clostridium_sensu_stricto_1),這與石寶明等[2]報道的金華豬和藏豬腸道微生物豐度較高的菌屬種類有著較高的一致性。

3.2 盲腸與結腸微生物菌群對宿主背膘厚、血清生化指標表型的影響

萊蕪豬是我國華北黑豬的典型代表,具有很強的脂肪沉積、抗病和耐粗飼能力。DLY豬是我國養豬生產中應用最為廣泛的商品豬種,以生長速度快、瘦肉率高著稱。本研究發現,除了背膘厚外,萊蕪豬與DLY豬的部分血清生化指標(TP、GLU、TG和HDL-C)也存在著顯著差異性(P<0.05,表1)。血清生化指標是反映機體健康狀況、內環境平衡及養分消化代謝的重要參數,這些指標可能通過影響消化代謝進而影響生產性能。本研究開展了腸道微生物菌屬豐度與背膘厚、兩品種間顯著差異血清生化指標的相關性分析,發現了一些與這些性狀顯著相關的腸道微生物類群(P<0.05),這些微生物類群可能對萊蕪豬和DLY豬的重要性狀有重要影響。

結腸中的甲烷短桿菌屬(Methanobrevibacter)與背膘厚呈顯著正相關。甲烷短桿菌屬(Methanobrevibacter)屬于廣古菌門(Euryarchaeota),可能在宿主的能量代謝和脂肪沉積中起著至關重要的作用,將人類甲烷短桿菌接種到無菌小鼠的模型試驗證明了這一點[33]。而且,人類肥胖癥研究中發現,肥胖個體利用h2產甲烷菌數量明顯高于瘦人個體[34]。梭桿菌屬(Fusobacterium)和大腸桿菌-志賀菌屬(Escherichia-Shigella)與TG呈顯著負相關。人類肥胖研究證明梭桿菌屬(Fusobacterium)和大腸桿菌-志賀菌屬(Escherichia-Shigella)的豐度占比在肥胖人群中同樣大幅度增加[35]。而且梭桿菌屬(Fusobacterium)和大腸桿菌-志賀菌屬(Escherichia-Shigella)均為萊蕪豬LEfSe差異判別分析結果中的特異菌,說明梭桿菌屬(Fusobacterium)和大腸桿菌-志賀菌屬(Escherichia-Shigella)可能通過影響萊蕪豬血清中TG從而影響了萊蕪豬的脂肪沉積。萊蕪豬結腸中的活潑瘤胃球菌([Ruminococcus]_gnavus_group)相對豐度遠高于DLY豬(3.11%vs. 0.22%),擬桿菌屬(Bacteroides)為LEfSe物種差異判別分析中鑒定出的萊蕪豬盲腸微生物菌屬。以往研究表明,擬桿菌屬(Bacteroides)在飼喂高纖維日齡的豬腸道中豐度顯著提高[36],活潑瘤胃球菌([Ruminococcus]_gnavus_group)是豬腸道微生物中纖維素和半纖維素分解的重要細菌菌屬[37],這兩個菌屬,通過分解宿主消化系統的纖維素來獲取營養,可能是導致萊蕪豬高纖維消化能力的主要因素之一,當然這也可能是由于萊蕪豬飼喂高纖維傳統日糧造成的菌群差異。羅伊氏乳桿菌(Lactobacillus_reuteri)是國際上公認的益生乳酸菌,可改善腸道菌群分布,拮抗有害菌定植,避免罹患腸道疾病。前期研究表明,羅伊氏乳桿菌Lactobacillus_reuteri與背膘厚呈負相關[38],具有降低膽固醇[39]和預防肥胖[40]的作用。本研究中羅伊氏乳桿菌(Lactobacillus_reuteri)在DLY豬結腸中豐度占比遠高于萊蕪豬(2.34%vs. 0.33%),為DLY豬LEfSe差異判別分析結果中的特異菌,故可能是影響DLY豬低背膘厚、高生長速度的菌屬之一。

動物腸道微生物與宿主之間相互選擇與適應,構成一個相對穩定的超有機體。宿主遺傳特征、飲食、環境和抗生素使用等因素均可影響腸道微生物的多樣性和功能[41]。本研究初步探索自然生產狀態下萊蕪豬腸道微生物的特點,以后研究中需要在控制飼養環境、飼糧和體重(或日齡)一致的情況下進一步研究分析,以驗證本研究中結果。

4 結 論

本研究分析了自然生產狀態下萊蕪豬和DLY豬育肥階段腸道微生物菌群組成、結構和功能,初步探索了萊蕪豬腸道微生物的特點。萊蕪豬具有與DLY豬不同的盲腸和結腸微生物菌群多樣性和相對豐度,且萊蕪豬盲腸微生物菌群多樣性顯著高于DLY豬;萊蕪豬腸道內的一些微生物菌群與脂肪沉積、耐粗飼和抗逆性強等有關,可能是影響萊蕪豬特性的原因之一。