动物营养学报    2018, Vol. 30 Issue (10): 4071-4080    PDF    
高通量测序技术分析蓝狐肠道菌群多样性
陈双双, 司华哲, 穆琳琳, 李光玉, 刘晗璐     
中国农业科学院特产研究所, 长春 130117
摘要: 本试验旨在研究蓝狐肠道菌群组成及其多样性。选取8只健康育成期(5~6月龄)雄性舍饲蓝狐,采集新鲜粪便后采用高通量测序技术对蓝狐肠道菌群组成及多样性进行分析。结果显示:从8只健康蓝狐粪便中共获得569 930条有效序列、操作分类单元(OTU)范围为468~574,分布于16个门的209个属。在门水平上以厚壁菌门(Firmicutes)、拟杆菌门(Baeteroidetes)、放线菌门(Actinobacteria)、变形菌门(Proteobacteria)、梭杆菌门(Fusobacteria)为主,其中厚壁菌门所占比例最高,为62.97%,其余依次为22.05%、8.89%、5.15%、0.88%。在属水平上以链球菌属(Streptococcus)所占比例最高,为11.75%,其次为乳杆菌属(Lactobacillus)(9.86%)、普雷沃菌属(Prevotella)(9.28%)、巨型球菌属(Megasphaera)(8.21%)、柯林斯菌属(Collinsella)(7.27%)、Blautia(7.08%)、拟杆菌属(Bacteroides)(5.64%)。通过高通量测序技术发现蓝狐具有复杂的肠道菌群结构,且个体间存在较大差异。
关键词: 蓝狐     肠道菌群     高通量测序     多样性    
Analysis of Intestinal Microbiota Diversity of Blue Foxes by High-Throughput Sequencing Technology
CHEN Shuangshuang, SI Huazhe, MU Linlin, LI Guangyu, LIU Hanlu     
Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun 130117, China
Abstract: The aim of this experiment was to study the composition and diversity of intestinal microbiota of blue foxes (Alopex lagopus). Eight healthy growing female blue foxes (5 to 6 months of age) under the stall-feeding condition were selected to collect fresh faecal samples, and then the composition and diversity of intestinal microbiota were analyzed by the high-throughput sequencing technology. The results showed as follows:a total of 569 930 valid sequences were identified in faecal samples from 8 healthy blue fox. The numbers of operational taxonomic unit (OUT) ranged from 468 to 574, representing 16 phyla and 209 genera. The dominant bacterial phyla in turn were Firmicutes (62.97%), Bacteroidtes (22.05%), Actinobacteria (8.89%), Proteobacteri (5.15%) and Fusobacteria (0.88%). The dominant bacterial genera in turn were Streptococcus (11.75%), Lactobacillus (9.86%), Prevotella (9.28%), Megasphaera (8.21%), Collinsella (7.27%), Blautia (7.08%) and Bacteroides (5.64%). The present study represents that there is a complex intestinal microbiota structure of blue foxes by the high-throughput sequencing technology, and there were bigger differences among individuals.
Key words: blue foxes     intestinal microbiota     high-throughput sequencing     diversity    

蓝狐(Alopex lagopus),又名北极狐,属于脊索动物门,哺乳纲,食肉目,犬科,北极狐属[1]。蓝狐作为毛皮动物在世界范围内被大量饲养,我国的养殖地主要集中于胶东半岛、河北以及东北三省等地。蓝狐为适应寒冷的自然气候条件,体内常储存大量脂肪,对高蛋白质、高脂肪的动物性饲料有较高的利用率,并可较大范围的利用动物性饲料成分,其中包括各种鱼类、畜禽屠宰副产品(去骨的牛肉、羊肉、猪肉、兔肉、鸡肉等)以及畜禽的肝脏(如兔肝、牛肝、羊肝、鸡肝、鸭肝、鹅肝)等。目前,大多数研究集中于蓝狐的繁殖育种、营养需要量、毛皮发育等研究方面[2-5],鲜有关于蓝狐肠道菌群组成和多样性的研究。肠道微生物系统是动物本身最复杂的微生态系统,肠道中大量微生物参与营养物质的吸收、分布、代谢以及机体的免疫生理,肠道菌群的失调通常会导致动物机体功能紊乱[6-7]。因此,了解健康动物正常的肠道菌群分布,对维持肠道菌群平衡、保证动物健康有重要意义。作为肉食性单胃动物,狐的胃肠道微生物主要为大肠微生物[8],其肠道中具体含有哪些微生物,以及这些微生物在消化生理中发挥的作用,仍不清楚。本试验采用高通量测序技术对蓝狐肠道菌群组成情况进行分析,研究蓝狐肠道菌群多样性,为进一步揭示蓝狐肠道菌群对营养物质的代谢利用的作用提供微生物学基础。

1 材料与方法 1.1 试验材料

2017年9月17号从河北省昌黎县(东经39.72°、北纬119.15°)选取8只健康的育成期(5~6月龄)雄性舍饲蓝狐(F1~F8),其精神状态良好、采食均衡正常,生活环境干净卫生。现场采集蓝狐新鲜粪便时,等待蓝狐排出粪便后立即用无菌分装袋收集,防止落地被污染。采集蓝狐新鲜粪便15 g分装于5 mL的冻存管中(分装于3个管,每管5 g),放入液氮罐中带回实验室,备后续分析。蓝狐饲粮为养殖场中日常饲粮,其组成及营养水平如表 1所示。

表 1 蓝狐饲粮组成及营养水平(风干基础) Table 1 Composition and nutrient levels of the diet for blue foxes (air-dry basis)
1.2 样品处理 1.2.1 微生物总DNA的提取

蓝狐粪便样品中微生物基因组DNA的提取按照Fast DNA® SPIN as Kit for Feces说明书进行,具体如下:1)取0.5 g粪便样品放入Lysing Matrix E管中,加入825 μL磷酸盐缓冲液和275 μL PLS溶液,混合均匀。2)14 000×g离心5 min,弃上清。3)加入978 μL磷酸盐缓冲液和122 μL MT缓冲液,混合均匀。4)将样品放入FastPrep® 24设备中,以6.0 m/s运行40 s。5)14 000×g离心5 min,取上清放入2 mL离心管。6)取250 μL PPS溶液加入到上清中,混合均匀后放于4 ℃静止10 min,14 000×g离心2 min。7)取上清放入2 mL离心管中,加入1 mL Binding Matrix Solution放入振动器振荡3 min。8)14 000×g离心2 min,弃上清。9)加入1 mL清洗液1混合均匀,取600 μL转移到SPIN Filter管中,14 000×g离心2 min,弃掉底部液体。10)加入500 μL清洗液, 混合均匀,14 000×g离心2 min,弃掉底部液体。11)继续14 000×g离心2 min,弃掉底部液体。12)将整个SPIN Filter管转移到新的Catch管中,加入100 μL TES溶液,14 000×g离心2 min,弃掉整个SPIN Filter管,将DNA产物放于-20 ℃保存备用。

基因组DNA的鉴定:提取的基因组DNA经1%的琼脂糖凝胶电泳(电压120 V,电泳23 min,Marker为λⅢ DNA Maker)检测,DNA产物在14 000 bp左右,并且DNA条带亮度较高,无明显降解、无杂质,可进行下一步试验[9]。将目的DNA送至北京诺禾致源科技股份有限公司,基于Illumina Miseq测序平台,设备名称为HiSeq PE250,测序长度为250 bp。

1.2.2 目的基因的扩增与测序

以样品中提取的总DNA为模板,细菌选用16S rRNA V3/V4区合成特定引物,利用338F(5′-ACTCCTACGGGAGGCAGCA-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)引物进行PCR扩增。PCR扩增体系(50 μL):上游引物(338F)1 μL,下游引物(806R)1 μL,DNA样品1 μL,Taq酶25 μL,ddH2O 22 μL。PCR扩增条件:95 ℃预变性5 min;95 ℃变性30 s,55 ℃退火30 s,72 ℃复性2 min,25个循环;72 ℃延伸10 min。扩增产物4 ℃保存,经2%的琼脂糖凝胶电泳检测。

1.2.3 生物信息学分析

通过对原始数据进行拼接、过滤去杂得到有效信息,从而使信息结果更加准确、可靠。对有效信息进行聚类分析和物种分类分析,设定具有97%以上相似度的序列为1个操作分类单元(OTU),并在各分类水平统计每个样品的菌群组成。根据OTU聚类结果,进行菌群丰富度、Alpha多样性计算,其中Ace和Chao1指数用来表示菌群的丰富度,Simpson和Shannon指数用来表示菌群的多样性。在门和属2个分类水平上统计样品物种的丰度和多样性。

1.3 数据处理

结果以平均值±标准误表示,数据用Excel 2013进行整理并用SAS 9.4软件进行分析。

2 结果 2.1 多样性指数

表 2所示,8个样本共得到654 716条序列,单个样本序列范围为57 002~95 638条,其中共得到569 930条有效序列(有效序列指tags,下文同),单个样本序列范围为52 758~84 112条。对样本进行聚类分析,在97%相似度下的OTU生物信息统计显示,共得到4 071个OTU(样本范围为468~574)。丰富度指数:Ace指数平均值为559,范围为480~699;Chao1指数平均值为595,范围为465~1036。多样性指数:Shannon指数平均值为5.27,范围为4.05~5.80;Simpson指数平均值为0.936,范围为0.847~0.957。每个样本的具体参数见表 2

表 2 样本测序数据分析 Table 2 Analysis of sample sequencing data

利用mothur对97%相似度的OTU做稀释曲线(图 1-A),可知,当测序数量超过40 000时,各样本(F8除外)曲线逐渐趋于平缓,同时Shannon曲线(图 1-B)也趋于平坦。

图 1 稀释曲线(A)和Shannon曲线(B)分析 Figure 1 Analysis of rarefaction curves (A) and Shannon curves (B)
2.2 门分类水平下菌群结构

在门分类水平上分析,这些菌群在分类学上归属于16个门,其中厚壁菌门(Firmicutes)在蓝狐肠道菌群中的所占比例最多,占62.97%;拟杆菌门(Bacteroidetes)为第二大优势菌门,所占比例为22.05%;随后依次为放线菌门(Actinobacteria),占8.89%;变形菌门(Proteobacteria),占5.15%;梭杆菌门(Fusobacteria),占0.88%。以上5个菌门为蓝狐肠道菌群中的优势菌门,共占99.94%,具体门分类水平的分布见图 2。8只蓝狐中F4与其他蓝狐的肠道菌群相似性较低,F4肠道中Bacteroidetes所占比例最高,而其他蓝狐肠道中Firmicutes为最优势菌门。

Firmicutes:厚壁菌门;Bacteroidetes:拟杆菌门;Actinobaeteria:放线菌门;Proteobacteria:变形菌门;Fusobacteria:梭杆菌门;Cyanobacteria:蓝藻菌门;Verrucomicrobia:疣微菌门;Spirochaetes:螺旋体门;Tenericutes:柔壁菌门;Fibrobacteres:纤维杆菌门;Others:其他。 图 2 基于Unweighted Unifrac距离的UPGMA聚类树 Figure 2 UPGMA clustering tree based on Unweighted Unifrac distance
2.3 属分类水平下菌群结构

在属分类水平上分析,这些菌群在分类学上归属于209个属,其中有30个属的相对丰度占总体的85.95%。这30个属中链球菌属(Streptococcus)所占比例最高,为11.75%;随后比例较高的依次为乳杆菌属(Lactobacillus),占9.86%;普雷沃菌属(Prevotella),占9.28%;巨型球菌属(Megasphaera),占8.21%;柯林斯菌属(Collinsella),占7.27%;Blautia,占7.08%;拟杆菌属(Bacteroides),占5.64%;拟普雷沃菌属(Alloprevotella),占5.50%;Peptoclostridium,占5.07%;巨单胞菌属(Megamonas),占3.55%。具体属分类水平的相对丰度见图 3

Others:其他;Fusobacterium:梭杆菌属; Roseburia:罗斯氏菌属;Dialister:杆菌属;Enterococcus:肠球菌属;Ruminococcus:瘤胃球菌属;Phascolarctobacterium:考拉杆菌属;Anaerobiospirillum:厌氧螺菌属;Bifidobacterium:双歧杆菌属;Succinivibrio:琥珀酸弧菌属;Desulfovibrio:脱硫弧菌属;Sutterella:萨特氏菌属;Megamonas:巨单胞菌属;Bacteroides:拟杆菌属;Alloprevotella:拟普雷沃菌属;Collinsella:柯林斯菌属;Prevotella:普雷沃菌属;Megasphaera:巨型球菌属;Lactobacillus:乳杆菌属;Streptococcus:链球菌属。 图 3 蓝狐肠道菌群属分类的相对丰度 Figure 3 Relative abundance of intestinal bacterial communities at genus level of blue foxes

在属水平上每只蓝狐肠道菌群中排名靠前的10种菌属所占比例见表 3。F1肠道菌群中最优势菌属为Collinsella,占17.45%;F2、F3、F4肠道菌群中最优势菌属为Prevotella,分别占14.52%、13.06%、22.16%;F5、F7肠道菌群中最优势菌属为Lactobacillus,分别占17.56%、21.55%;F6、F8肠道菌群中最优势菌属为Streptococcus,分别占15.80%、25.72%。与其他蓝狐相比,F8肠道菌群中StreptococcusLactobacillusMegasphaera所占比例较高,3个菌属共占71.00%。

表 3 蓝狐肠道不同菌属所占比例 Table 3 Proportion of intestinal flora at different genus of blue foxes
3 讨论 3.1 多样性指数分析

一般通过Alpha多样性分析反映样本的丰富度和多样性。本试验中,所有样本测序覆盖率(coverage)均在0.997以上,说明测序结果对蓝狐粪便样品中的细菌覆盖率高,测序结果合适,可满足肠道菌群多样性的分析需求[10]。当测序数量超过40 000时,各样本(F8除外)曲线逐渐趋于平缓,表明该测序数据量合理,随着数据量的增加几乎不产生新的OTU[11]。同时,Shannon曲线也趋于平坦,表明测序数据量足够大,从而能够反映出样品中绝大多数菌群[12]。Shannon指数越大则表示该样本中微生物多样性越丰富,由表 2中Shannon指数数据可知,F8肠道中菌群多样性较低。由图 1-A显示,较其他个体而言,F8个体的稀释曲线存在一定的坡度,可能是其取样深度不够,从而导致F8个体肠道菌群多样性较低。从整体来看,蓝狐肠道菌群中OTU平均值为509,Chao1指数平均值为595,Shannon指数平均值为5.27。范忠原等[13]研究发现,舍饲水貂肠道菌群中OTU为294,Chao1指数为350,Shannon指数为3.74。因此,与毛皮动物水貂相比,蓝狐肠道菌群具有较高的丰度与多样性。

3.2 门分类水平下菌群结构分析

从门分类水平上看,蓝狐肠道菌群中主要优势菌门为Firmicutes和Bacteroidetes。8个样本中F4肠道中Bacteroidetes所占比例最高,而其他蓝狐肠道中Firmicutes为最优势菌门。研究表明,肠道中的Firmicutes主要是对碳水化合物及蛋白质起水解作用,Bacteroidetes则是对类固醇、胆汁酸及多糖起代谢作用,从而促进多糖的吸收以及蛋白质的合成[14-15]。肠道内Firmicutes与Bacteroidetes的比例对动物的肥胖程度有影响。Turnbaugh等[16]发现肥胖小鼠的肠道微生物中Firmicutes细菌的丰富度增加,Bacteroidetes细菌的丰富度降低,并发现饮食过程中对能量摄入量增多。饮食的改变会引起动物肠道菌群的改变,当动物饮食从低脂肪、高纤维素的饮食转变为高脂肪、高糖的饮食时,肠道菌群中Bacteroidetes细菌所占比例显著下降[17],从而导致肥胖增加。因此,F4肠道中菌群可能对高脂、高糖的代谢作用较强,但本试验未对试验动物的体重水平进行观察。除此之外,动物肠道菌群的组成和数量受动物的年龄、性别、生活环境、饲粮组成、生理状态、遗传因素[18-20]等的影响。8只蓝狐在生活环境和饲粮组成相同的情况下,F4和F8个体的肠道菌群组成与其他个体的肠道菌群组成存在差异性,其原因可能与动物个体的生理状态和遗传因素有关系。因此,相同环境下的不同个体间肠道菌群组成也存在一定的差异性。

研究发现,犬科动物肠道菌群在种类及数量方面存在较大差异。生活在野生环境下的健康狼(Canis lupus)的粪便中最丰富的菌门是Firmicutes(60.0%),其次是Bacteroidetes(16.9%)、Proteobacteria(9.2%)、Fusobacteria(9.2%)和Actinobacteria(4.6%)[21];豺(Cuon alpinus)的肠道菌群中较丰富的5个菌门依次为Firmicutes(20.97%~44.01%)、Bacteroidetes(21.63%~38.97%)、Proteobacteria(9.33%~17.60%)、Fusobacteria(9.11%~17.90%)和Actinobacteria(1.22%~2.87%)[22];家犬的肠道菌群中最丰富的菌门为Firmicutes(47.7%),其次为Proteobacteria(23.3%)、Fusobacteria(16.6%)和Bacteroidetes(12.4%)[23]。宠物犬的肠道菌群中最丰富的菌门为Firmicutes(64.17%),其次是Bacteroidetes(19.89%)、Fusobacteria(13.58%)、Actinobacteria(1.5%)和Proteobacteria(0.86%)[24]。与其他犬科动物相比,蓝狐肠道菌群中Actinobacteria所占比例较高,该菌门中大部分细菌为红蝽菌科(Coriobacteriaceae)和双歧杆菌科(Bifidobacteriaceae)(共占Actinobacteria的99.99%),而在狼和豺肠道菌群中Bifidobacteriaceae的数量很少[21-22]。蓝狐肠道菌群Bifidobacteriaceae中仅有双歧杆菌(bifidobacteria)。研究发现人类从出生就开始携带双歧杆菌,但随着年龄的增加双歧杆菌的数量逐渐减少,并且发现双歧杆菌的主要产物为乳酸、乙酸等,其可调节肠道pH,抑制有害菌增殖[25];同时,研究发现双歧杆菌在能够维持机体肠道健康和减缓肠道老化方面具有重要作用[26]。本试验中育成期蓝狐年龄(5~6月龄)较小,因此其肠道中Bifidobacteria数量较多。与其他动物相比,蓝狐肠道菌群中Actinobacteria所占比例较高。研究发现,在应激情况下Coriobacteriaceae细菌会明显增加,该科细菌的大量存在易引起动物焦虑、紧张等情绪[27]。研究表明,Coriobacteriaceae还与高胆固醇血症[28]、肥胖、2型糖尿病[29]等密切相关。与狼和蓝狐比较,犬的肠道菌群中Fusobacteria所占比例较高[21-24, 30],Bacteroidetes所占比例较低。根据《伯吉氏系统细菌学手册》,可把Bacteroidetes分为4个菌纲:拟杆菌纲(Bacteroidia)、黄杆菌纲(Flavobacteria)、鞘脂杆菌纲(Sphingobacteria)、纤维黏网菌纲(Cytophagia)[31]。蓝狐肠道菌群中Bacteroidia在Bacteroidetes中所占比例为99.77%。家犬肠道菌群中Bacteroidetes比例较低的主要原因可能在于家犬的饮食受到人为干预,饮食以狗粮为主,其饲粮中除蛋白质以外碳水化合物成分增加,与其他犬科动物相比饮食中的脂肪含量下降。将本试验结果与其他犬科动物肠道菌群结构进行比较,发现蓝狐肠道菌群多样性特征与其他犬科动物基本相同,仅在Fusobacteria、Actinobacteria的数量上存在较大差别,但蓝狐肠道菌群差异对宿主的生活习性是否存在影响还需进一步深入研究。

水貂作为我国饲养的主要毛皮动物之一,在国内被大量养殖,研究发现其肠道菌群中较丰富的5个菌门为Firmicutes(60.0%)、Bacteroidetes(16.2%)、Fusobacteria(11.5%)、Actinobacteria(5.9%)和Proteobacteria(5.3%)[13]。与蓝狐肠道菌群相比,除Fusobacteria在水貂肠道菌群中所占比例较高外,其他菌门种类与蓝狐相似。研究发现,饲粮中纤维含量的增加会导致Fusobacteria细菌数量减少[32]。水貂为肉食动物,与蓝狐饲粮相比,水貂饲粮中纤维含量较低可能导致水貂肠道细菌中Fusobacteria数量增加。研究发现,马作为草食动物,其肠道菌群中没有发现Fusobacteria细菌的存在[33]。同时,研究发现人肠道菌群中Fusobacteria所占比例较低,不足1%[34]。因此,Fusobacteria细菌可能受饲粮中植物性饲料组成的影响。

3.3 属分类水平下菌群结构分析

从属分类水平上看,蓝狐肠道菌群中最大的优势菌属为Streptococcus,与舍饲水貂的研究结果一致[13]Streptococcus在动物肠道中能够发酵简单的糖类、分解蛋白质,其与宿主体内能量代谢紧密相关[35]Lactobacillus是蓝狐肠道菌群中的第二大菌属。Lactobacillus作为益生菌,近年来受到越来越多的关注。Lactobacillus可与其他益生菌相互作用,治疗肠易激综合症(IBS)等疾病[36],其在动物肠道中能够发酵分解碳水化合物,产生乳酸、丁酸等,使肠道内pH下降,抑制有害菌的生长;Lactobacillus还可以产生过氧化氢,抑制肠道中念珠菌等真菌的毒性和生长[37];同时,Lactobacillus还可以激活免疫细胞,破坏机体肠道中的病原微生物,从而提高机体的免疫力[38]StreptococcusLactobacillus同属于Firmicutes中的乳杆菌目(Lactobacillales),该目在蓝狐肠道菌群中所占比例为22.5%。Prevotella为Bacteroidetes中的拟杆菌目(Bacteroidales),具有促进蛋白质分解的作用,其在瘤胃中还能协同其他微生物分解纤维素,提高纤维素的降解能力[39]。研究发现,非洲布基纳法索(Burkina Faso)儿童肠道中普氏菌属占53%,而同龄的欧洲儿童却很少[40],饮食习惯与肠道菌群的组成密切相关,以蛋白质和动物性脂肪为主要食物的西方人肠道菌群多以Bacteroides为主,而以碳水化合物为主要食物的非洲人,肠道菌群则主要以Prevotella为主[41]Collinsella属于Actinobacteria中的Coriobacteriales,研究发现该属中的细菌能够发酵碳水化合物产生甲酸、乳酸盐、乙醇等,并与体重、体重指数(BMI)呈负相关[42]Blautia属于Firmicutes,该菌属能将肠道中的气体转化成乙酸,从而可清除肠道中的气体[10]Blautia属于毛螺旋菌科(Lachnospiraceae),该菌科在瘤胃动物中主要具有降解半纤维素的能力[43],而这些细菌同样大量存在于蓝狐肠道中,但该类细菌在蓝狐肠道中的作用仍有待进一步研究确定。因此,蓝狐肠道菌群中包括能够促进类固醇、胆汁酸及多糖代谢的Bacteroidetes,能够促进蛋白质分解的Streptococcus,能够发酵分解碳水化合物、抑制有害菌的生长的LactobacillusCollinsella,以及能够促进纤维降解的PrevotellaBlautia等,从而组成丰富多样化的肠道菌群。

蓝狐为肉食性动物,但人们为减少蓝狐的养殖成本,在养殖过程中对蓝狐饲粮进行人为干预,使蓝狐饲粮的动物性蛋白质在饲粮中的比例降低,而植物性蛋白质在饲粮中的比例增加。目前并未发现对野生蓝狐和以动物性蛋白质为主的蓝狐肠道菌群的报道,本试验蓝狐饲粮以植物性蛋白质为主,因此其肠道菌群不能代表以肉食(野生蓝狐或者以动物性蛋白质为主的蓝狐)为主的蓝狐的肠道菌群,以植物性蛋白质为主的蓝狐肠道菌群和以肉食为主的蓝狐肠道菌群差异性无从得知。但在野生犬科动物(狼、沙狐、赤狐、乌苏里貉、豺)肠道菌群中并未发现纤维素降解菌Prevotella的存在[8],而在本试验饲粮条件下蓝狐肠道菌群中Prevotella为优势菌属,蓝狐肠道菌群中Prevotella增加的原因可能是其饲粮组成的改变导致其肠道菌群结构的改变。因此,后续试验需对野生蓝狐和以动物性蛋白质为主的蓝狐的肠道菌群组成进行研究,并对以植物性蛋白质为主和以肉食为主的蓝狐肠道菌群组成的差异性进行研究。

属水平上个体间差异性分析显示,F8肠道菌群中StreptococcusLactobacillusMegasphaera所占比例较高,其他菌属所占比例较低,这种现象可能是样本的测序取样深度不够所致。除此之外,动物肠道菌群组成和多样性还受到动物的遗传因素的影响[20],因此,相同环境下的不同个体间肠道菌群组成和多样性也存在一定的差异性。这些优势菌属与营养物质代谢存在密切的联系,在动物肠道中主要起到分解蛋白质以及提高纤维素的降解能力的作用,但对于这种联系的机理仍需继续开展大量的研究。

4 结论

本试验饲粮条件下的8只蓝狐肠道菌群分布于16个菌门,主要为Firmicutes(62.97%)、Bacteroidetes(22.05%)、Actinobacteria(8.89%)、Proteobacteria(5.15%)和Fusobacteria(0.88%),占菌门总量的99.94%;所有菌门分布于209个菌属,组成丰富多样化的肠道菌群,在属的分类水平上Streptococcus所占比例最高,为11.75%,其次为Lactobacillus,占9.86%。

参考文献
[1]
HERSTEINSSON P. Population genetics and ecology of different colour morphs of arctic foxes Alopex lagopus in Iceland[J]. Finnisch Game Reserve, 1989, 46: 64-78.
[2]
郭俊刚, 贡筱, 张铁涛, 等. 饲粮中添加益生菌对冬毛期蓝狐生长性能、营养物质消化率及毛皮品质的影响[J]. 动物营养学报, 2014, 26(8): 2232-2239. DOI:10.3969/j.issn.1006-267x.2014.08.026
[3]
李仲玉, 任海轮, 王超, 等. 脂肪来源、添加量对蓝狐生产性能及毛皮品质的影响[J]. 黑龙江畜牧兽医, 2017(13): 240-243.
[4]
郭俊刚, 张铁涛, 吴学壮, 等. 饲粮蛋氨酸水平对繁殖期蓝狐繁殖性能的影响[J]. 动物营养学报, 2015, 27(1): 165-170. DOI:10.3969/j.issn.1006-267x.2015.01.020
[5]
张海华, 李光玉, 刘佰阳, 等. 狐营养研究进展[J]. 饲料工业, 2008, 29(13): 55-58. DOI:10.3969/j.issn.1001-991X.2008.13.017
[6]
胡健微.肠道微生物相互作用的研究及体外肠道模拟系统的构建和应用[D].硕士学位论文.长春: 吉林大学, 2017.
[7]
叶雷, 闫亚丽, 陈庆森, 等. 高通量测序技术在肠道微生物宏基因组学研究中的应用[J]. 中国食品学报, 2016, 16(7): 216-223.
[8]
曹先海.犬科动物粪便菌群组成及多样性分析[D].硕士学位论文.曲阜: 曲阜师范大学, 2011.
[9]
石聪, 李世瑞, 李跑, 等. 基于高通量测序浏阳豆豉不同发酵阶段微生物多样性分析[J]. 食品与发酵工业, 2017, 44(2): 27-32, 39.
[10]
曲巍, 张智, 马建章, 等. 高通量测序研究益生菌对小鼠肠道菌群的影响[J]. 食品科学, 2017, 38(1): 214-219.
[11]
杨俊花, 赵志辉, 郭文博, 等. 应用Illumina-MiSeq高通量测序技术分析脱氧雪腐镰刀菌烯醇对小鼠肠道菌群的影响[J]. 动物营养学报, 2017, 29(1): 158-167. DOI:10.3969/j.issn.1006-267x.2017.01.018
[12]
邓晓影, 张宾, 汤贺, 等.基于高通量测序的南美白对虾中微生物群落分析[J/OL].食品科学, [2017-10-30].http://kns.cnki.net/kcms/detail/11.2206.TS.20171030.1202.002.html.
[13]
FAN Z Y, LIU H L, WANG F, et al. Determination of bacterial community composition in the distal gut of mink (Mustela vison) by high-throughput sequencing[J]. Microbiology China, 2016, 43(1): 123-130.
[14]
BÄCKHED F, LEY E, SONNENBURY J L, et al. Host-bacterial mutualism in the human intestine[J]. Science, 2005, 307(5717): 1915-1920.
[15]
XU J, BJURSELL M K, HIMROD J, et al. A genomic view of the human-Bacteroides thetaiotaomicron symbiosis[J]. Science, 2003, 299(5615): 2074-2076. DOI:10.1126/science.1080029
[16]
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.
[17]
THOMAS F, HEHEMANN J H, REBUFFET E, et al. Environmental and gut Bacteroidetes:the food connection[J]. Frontiers in Microbiology, 2011, 2: 93.
[18]
SIMPSON J M, MARTINEAU B, JONES W E, et al. Characterization of fecal bacterial populations in canines:effects of age, breed and dietary fiber[J]. Microbial Ecology, 2002, 44(2): 186-197. DOI:10.1007/s00248-002-0001-z
[19]
LU J R, IDRIS U, HARMON B, et al. Diversity and succession of the intestinal bacterial community of the maturing broiler chicken[J]. Applied and Environmental Microbiology, 2003, 69(11): 6816-6824. DOI:10.1128/AEM.69.11.6816-6824.2003
[20]
谭振, 翟丽维, 陈少康, 等. 肠道微生物与宿主遗传背景互作关系的研究进展[J]. 中国畜牧杂志, 2016, 52(5): 84-88. DOI:10.3969/j.issn.0258-7033.2016.05.018
[21]
ZHANG H H, CHEN L. Phylogenetic analysis of 16S rRNA gene sequences reveals distal gut bacterial diversity in wild wolves (Canis lupus)[J]. Molecular Biology Reports, 2010, 37(8): 4013-4022.
[22]
WU X Y, ZHANG H H, CHEN J, et al. Comparison of the fecal microbiota of dholes high-throughput Illumina sequencing of the V3-V4 region of the 16S rRNA gene[J]. Applied Microbiology and Biotechnology, 2016, 100(8): 3577-3586.
[23]
SUCHODOLSKI J S, CAMACHO J, STEINER J M. Analysis of bacterial diversity in the canine duodenum, jejunum, ileum, and colon by comparative 16S rRNA gene analysis[J]. FEMS Microbiology Ecology, 2008, 66(3): 567-578. DOI:10.1111/fem.2008.66.issue-3
[24]
KIM J, AN J U, KIM W, et al. Differences in the gut microbiota of dogs (Canis lupus familiaris) fed a natural diet or a commercial feed revealed by the Illumina MiSeq platform[J]. Gut Pathogens, 2017, 9: 68.
[25]
DI GIOIA D, ALOISIO I, MAZZOLA G, et al. Bifidobacteria:their impact on gut microbiota composition and their applications as probiotics in infants[J]. Applied Microbiology and Biotechnology, 2014, 98(2): 563-577.
[26]
KIM O S, CHO Y J, LEE K, et al. Introducing EzTaxon-e:a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species[J]. International Journal of Systematic and Evolutionary Microbiology, 2012, 62(3): 716-721.
[27]
BENDTSEN K M B, KRYCH L, SØRENSEN D B, et al. Gut microbiota composition is correlated to grid floor induced stress and behavior in the BALB/c mouse[J]. PLoS One, 2012, 7(10): e46231. DOI:10.1371/journal.pone.0046231
[28]
ISLAM K B M S, FUKIYA S, HAGIO M, et al. Bile acid is a host factor that regulates the composition of the cecal microbiota in rats[J]. Gastroenterology, 2011, 141(5): 1773-1781. DOI:10.1053/j.gastro.2011.07.046
[29]
YIN X C, PENG J H, ZHAO L P, et al. Structural changes of gut microbiota in a rat non-alcoholic fatty liver disease model treated with a Chinese herbal formula[J]. Systematic and Applied Microbiology, 2013, 36(3): 188-196.
[30]
MENTULA S, HARMOINEN J, HEIKKILÄ M, et al. Comparison between cultured small-intestinal and fecal microbiotas in Beagle dogs[J]. Applied and Environmental Microbiology, 2005, 71(8): 4169-4175. DOI:10.1128/AEM.71.8.4169-4175.2005
[31]
KRIEG N R, LUDWIG W, WHITMAN W, et al. Bergey's manual of systematic bacteriology[M]. New York: Springer, 2010: 89-100.
[32]
MIDDELBOS I S, BOLER B M V, QU A, et al. Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing[J]. PLoS One, 2010, 5(3): e9768. DOI:10.1371/journal.pone.0009768
[33]
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.
[34]
YATSUNENKO T, REY F E, MANARY M J, et al. Human gut microbiome viewed across age and geography[J]. Nature, 2012, 486(7402): 222-227. DOI:10.1038/nature11053
[35]
LEY R E, PETERSON D A, GORDON J I. Ecological and evolutionary forces shaping microbial diversity in the human intestine[J]. Cell, 2006, 124(4): 837-848. DOI:10.1016/j.cell.2006.02.017
[36]
FORD A C, QUIGLEY E M M, LACY B E, et al. Efficacy of prebiotics, probiotics, and synbiotics in irritable bowel syndrome and chronic idiopathic constipation:systematic review and meta-analysis[J]. The American Journal of Gastroenterology, 2014, 109(10): 1547-1561.
[37]
ERDOGAN A, RAO S S. Small intestinal fungal overgrowth[J]. Current Gastroenterology Reports, 2015, 17(4): 16. DOI:10.1007/s11894-015-0436-2
[38]
范忠原.家养水貂肠道微生物多样性分析[D].硕士学位论文.北京: 中国农业科学院, 2015.
[39]
杨琦玥, 黄勇, 陈亚冰, 等. 16S rRNA高通量测序技术筛选牦牛瘤胃细菌基因组DNA提取方法及菌群结构[J]. 中国农业科学, 2017, 50(5): 932-941.
[40]
DE FILIPPO C, CAVALIERI D, DI PAOLA M, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa[J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(33): 14691-14696.
[41]
WU G D, CHEN J, HOFFMANN C, et al. Linking long-term dietary patterns with gut microbial enterotypes[J]. Science, 2011, 334(6052): 105-108.
[42]
熊健会, 李琳琳, 毛新民, 等. 新疆维、哈两民族粪样中奇异菌簇和产气柯林斯菌属水平与2型糖尿病的相关性研究[J]. 中国微生态学杂志, 2015, 27(06): 633-637, 641.
[43]
肖文萍.藏猪肠道微生物多样性的研究[D].硕士学位论文.杨凌: 西北农林科技大学, 2012.