畜牧兽医学报  2023, Vol. 54 Issue (6): 2619-2630. DOI: 10.11843/j.issn.0366-6964.2023.06.038    PDF    
大熊猫肠道菌群年龄演替规律分析
李蔚1, 张强2, 瞿嘉豪1, 吴亚平2, 胡若辰1, 贾若艺2, 郭如海1, 马清义2, 潘广林2, 王兴龙1     
1. 西北农林科技大学动物医学院, 杨凌 712100;
2. 秦岭大熊猫研究中心(陕西省珍稀野生动物救护基地), 周至 710400
摘要:本研究旨在探究不同年龄段大熊猫(Ailuropoda melanoleuca)肠道菌群多样性、结构及功能的演替规律。于秦岭大熊猫研究中心采集16只不同年龄段大熊猫的新鲜粪便, 基于16S rDNA高通量测序技术及层次聚类方法分为5个年龄组进行分析。结果表明, 不同年龄段大熊猫共有优势菌门为厚壁菌门(Firmicutes, 60.44%)和变形菌门(Proteobacteria, 38.69%), 且随年龄增长, 厚壁菌门的相对丰度逐渐增加, 变形菌门的相对丰度逐渐降低; 共有优势菌属为埃希菌属-志贺菌属(Escherichia-Shigella, 57.35%)和狭义梭菌属1(Clostridium_sensu_stricto_1, 15.17%), 且性别差异对大熊猫肠道菌群结构影响不大。不同年龄段大熊猫的肠道菌群Alpha多样性存在显著差异(P < 0.05), 呈现随年龄增长先增加后降低的趋势, 在18~22岁达最高; Beta多样性结果表明, 组间差异大于组内差异, 不同分组对样品差异的解释具有可信度。13~14岁成年大熊猫特有的优势菌属——魏斯菌属(Weissella, 9.975%)为肠道益生菌, 34岁老年大熊猫在氨基酸转运和代谢、信号转导、防御机制等通路的相对丰度不同于其他年龄段, 其特有的优势菌属——八叠球菌属(Sarcina, 53.29%)为条件致病菌。综上, 该研究认为大熊猫肠道菌群多样性随年龄增大先增加后降低, 不同年龄段共有和特有优势菌丰度及菌群功能存在差异, 随年龄变化呈现一定演替规律, 提示在大熊猫不同年龄阶段(尤其是老年阶段)的饲养过程中可补充相应缺乏的营养物质, 以维持肠道内微生态平衡。
关键词高通量测序    大熊猫    肠道菌群    年龄    
Analysis on the Age Succession of Intestinal Flora of Giant Panda
LI Wei1, ZHANG Qiang2, QU Jiahao1, WU Yaping2, HU Ruochen1, JIA Ruoyi2, GUO Ruhai1, MA Qingyi2, PAN Guanglin2, WANG Xinglong1     
1. College of Veterinary Medicine, Northwest A&F University, Yangling 712100, China;
2. Research Center for Qinling Giant Panda(Shaanxi Rare Wildlife Rescue Base), Zhouzhi 710400, China
Abstract: The study aimed to explore the succession of intestinal flora diversity, structure and function of giant pandas at different ages. Fresh feces of 16 giant pandas of different ages were collected in Research Center for Qinling Giant Panda, divided into 5 groups and analyze based on 16SrDNA high-throughput sequencing technology and hierarchical clustering method. The results showed that the main bacterial phyla in giant pandas at different ages consisted of Firmicutes (60.44%) and Proteobactera (38.69%), and the relative abundance of Firmicutes increased while Proteobacteria decreased with age; at the genus level, the bacteria consisted mainly of Escherichia Shigella (57.35%) and Clostridium_sensu_stricto_1 (15.17%), while the gender difference has little effect on the intestinal flora structure of giant pandas. There was a significant difference in the Alpha diversity of intestinal flora (P < 0.05), showing a trend of first increasing and then decreasing with age, reaching the highest in 18-22 years old. Beta diversity results show that the difference between groups is greater than the difference within groups, and the interpretation of sample differences by different groups has credibility. The dominant bacteria unique to adult giant pandas aged 13-14 years-Weissella (9.975%) is intestinal probiotics. The relative abundance in amino acid transport and metabolism, signal transduction, defense mechanisms and other pathways of 34 year old giant pandas is different from other age groups. Its unique dominant bacteria-Sarcina (53.29%) is conditioned pathogen. To sum up, this study believes that the intestinal flora diversity of giant pandas first increases and then decreases with age, and the abundance and function of common and unique dominant bacteria in different age groups are different, showing a certain succession law with age, suggesting that the corresponding lack of nutrients can be supplemented in the feeding process of giant pandas at different ages to maintain the intestinal microecology balance.
Key words: high throughput sequencing    panda    intestinal flora    age    

大熊猫(Ailuropoda melanoleuca)是以竹为主食的食肉目(Carnivora)动物,是我国特有的珍稀保护动物[1],因其具有独特的饮食结构备受关注。消化系统疾病是导致大熊猫等野生动物数量减少以及种群濒危的重要因素之一[2],当大熊猫的生活方式、饮食结构等因素发生变化时,可能导致肠道菌群的失衡,进而引发肠道性疾病并使机体处于亚健康状态甚至出现死亡[3-4],因此研究大熊猫肠道菌群的结构组成,对保护大熊猫具有重要的意义。

肠道菌群在大熊猫的生长发育、营养吸收、代谢以及免疫等方面发挥重要作用[5-6],由于目前暂未发现大熊猫体内含有与纤维素降解直接相关的基因,因此推测其主要依赖肠道菌群产生纤维素酶或编码纤维素代谢相关的基因来进行降解[7]。先前的研究表明,厚壁菌门(Firmicutes)和变形菌门(Proteobacteria)是大熊猫肠道内主要优势菌门[8],且年龄对大熊猫肠道菌群结构和数量会产生一定的影响[9],但对于不同年龄段大熊猫在肠道菌群多样性、结构组成和功能上的动态变化及演替规律研究尚浅。因此本试验采集了16只不同年龄段大熊猫的新鲜粪便,基于16 S rRNA高通量测序技术[10]和层次聚类结果分为不同年龄段的5组进行分析,以探究不同年龄段大熊猫肠道菌群多样性、结构组成和功能的演替规律,为大熊猫肠道性疾病的防治提供思路。

1 材料与方法 1.1 样本采集与分组

在秦岭大熊猫研究中心采集不同年龄段大熊猫的新鲜粪便共16份(表 1),采样期间大熊猫状态良好,无疾病等异常情况,粪样采集后用冰盒保存带回实验室,-20 ℃冷冻保存。

表 1 大熊猫基本信息 Table 1 Basic information of giant panda
1.2 总DNA提取、文库构建及测序

依照DNA提取试剂盒(天根生化科技有限公司)说明书步骤进行操作,提取粪便中所有细菌DNA,利用PCR扩增试剂盒结合引物515F-806R对细菌DNA的16 S rDNA V4区域进行扩增(前端引物:5′-GTGCCAGCMGCCGCGGTAA-3′,后端引物:5′-GGACTACHVGGGTWTCTAAT-3′),并对PCR产物的片段长度和浓度进行筛选,主带长度筛选条件为290~310 bp。按照标准流程进行建库操作,检测合格的文库用Illumina Nova 6000平台进行双末端250 bp测序。

1.3 数据分析

fastp 0.23.2[11]软件对测序原始数据进行质控,usearch 10.0软件的UPARSE聚类方法对非重复序列进行聚类,以97% 相似性水平将有效序列聚类成操作分类单元(operational taxonomic units, OTU),并基于SILVA数据库进行物种注释。ggplot2 3.3.6绘制样品的丰富度稀释曲线(richness rarefaction curve)以反映样品中微生物多样性信息,使用平均连接法(average-linkage clustering)的聚类分析方法和Chisq距离算法对样品进行层次聚类(hierarchical clustering),并将聚类结果与属水平上物种相对丰度整合。QIIME 2 2 021.2软件计算样品的Chao1、Simpson和Dominance等Alpha多样性指数[12-13],并基于单因素方差分析(one-way ANOVA)比较各指数间差异,R软件的vegan包进行非参数多元方差分析(non-parametric multivariate analysis of variance,ADONIS)和非度量多维尺度分析(non-metric multidimensional scaling,NMDS),均采用bray_curtis距离算法。

ggplot2包绘制门和属水平肠道菌群共有以及特有OTUs的venn图,及肠道菌群相对丰度的柱状堆叠图,基于PICRUSt2软件分析各组间在COG(clusters of orthologous groups)的L2水平和KEGG通路(KEGG pathway database)的L2水平上的功能组成,进行通路聚类分析,其中聚类方式选择average。数据以“ $ \bar x \pm {s_{\bar x}}$”表示,P>0.05为差异不显著,P<0.05为差异显著,P<0.01为差异极显著。

2 结果 2.1 测序数据信息

对下机数据进行质控、拼接和去噪后,16份粪便样品共获得1 235 229条高质量序列,平均每个样品77 202条序列。基于97%相似性水平进行OTU聚类,共获得1 707个OTU,包括12个门,22个纲,47个目,68个科,90个属。

以重抽样抽出的序列数比例作为横坐标,OTU水平下物种丰富度(species richness)作为纵坐标,进一步构造物种稀释性曲线。由图 1可知,随着抽出序列数比例的增加,肠道菌丰富度先显著升高,后趋于平缓,这表明随着样品序列数量不断增加,检测到大熊猫肠道菌的种类数目逐渐稳定,即便样品序列数继续增加,检测到新的肠道菌种类也是有限的。因此本试验的测序数据量足够大,能够检测到样品中绝大部分的微生物信息,可用于后续分析。

图 1 97%相似度水平下各样品的肠道菌群丰富度稀释曲线 Fig. 1 Dilution curve of intestinal flora abundance of each sample at 97% similarity level
2.2 大熊猫肠道菌群层次聚类分析

层次聚类分析能反映样品的肠道菌群结构及样品间的相似度。结果显示(图 2),编号为16的老年大熊猫肠道菌群结构单独聚为一个类别,其余大熊猫聚为一个类别。编号为1、2、3、4的亚成年大熊猫位于同一树杈分支点下,编号为5、6,编号为9、10及编号为13、14、15的成年大熊猫分别聚类到不同的树杈分支点下,且相同年龄的雌性和雄性大熊猫肠道菌群没有明显差异。依据层次聚类结果,将16只大熊猫分为5组(表 2),分别为2~3岁亚成年组(G1)、8~9岁成年组(G2)、13~14岁成年组(G3)、18~22岁成年组(G4)、34岁老年组(G5)。年龄相近的大熊猫肠道菌群结构相似性较高,年龄差异越大的肠道菌群结构差异越大,从而形成鲜明的菌群差异派别。

聚类树后的数字为大熊猫编号,黑色数字表示雌性大熊猫,棕色数字表示雄性大熊猫 The number behind the cluster tree is the giant panda number, black numerals refer to female giant pandas and brown numerals refer to male giant pandas 图 2 各样品在属水平的层次聚类分析 Fig. 2 Hierarchical cluster analysis of samples at genus level
表 2 大熊猫分组信息 Table 2 Group information of Giant pandas
2.3 大熊猫肠道菌群的多样性分析

Alpha多样性指数包括Chao 1指数、Simpson指数和Dominance指数等,能够分别衡量大熊猫肠道菌群丰富度、群落多样性以及物种均匀度等,其中Simpson指数值与群落多样性呈负相关。结果可知,不同年龄组的大熊猫在Chao1指数上无显著差异(P>0.05),2~3岁亚成年组和34岁老年组的Simpson指数显著高于18~22岁成年组(P<0.05),老年组的Dominance指数显著低于18~22岁成年组(P<0.05)(表 3)。上述结果表明,不同年龄段大熊猫在菌群丰富度上差异不大,2~3岁亚成年大熊猫的肠道菌群在群落多样性上明显低于18~22岁大熊猫,34岁老年大熊猫在肠道微生物群落多样性以及均匀度上明显低于18~22岁大熊猫。综合3个指数结果表明,大熊猫肠道菌群Alpha多样性随着年龄增大呈现先增加后降低的趋势,且在18~22岁达到最高。

表 3 不同年龄组大熊猫肠道菌群Alpha多样性指数($ \bar x \pm {s_{\bar x}}$) Table 3 Alpha diversity index of intestinal microflora of giant panda in different age groups($ \bar x \pm {s_{\bar x}}$)

Beta多样性能对不同样本间的微生物群落构成进行比较。NMDS分析能在二维平面判断群落组成差异,结果如图 3所示,2~3岁亚成年大熊猫的肠道菌群较为聚集,与其余4组存在一定距离,而8~9、13~14及18~22岁成年大熊猫肠道菌群分离较明显,组间距离近。ADONIS分析能判断分组效果是否显著,结果可知,除仅有34岁大熊猫所在的G6不参与组间差异比较外,其余各比较组在OTU水平差异程度均>0.25且<0.5,说明组间存在差异,但差异不大,且除G3 vs G4比较组外,其余差异比较组的P值均<0.05,说明不同分组对样品差异的解释具有可信度(表 4)。综合以上结果表明,虽然各组内肠道菌群存在差异,但组间差异仍大于组内差异,即该分组具有可信度。

图 3 不同年龄组间肠道菌群NMDS分析 Fig. 3 NMDS analyze of intestinal flora between different age groups
表 4 不同年龄比较组间ADONIS统计结果 Table 4 Statistical results of ADONIS between different age groups
2.4 不同年龄段大熊猫肠道菌群在门和属水平的演替

首先从门水平对大熊猫共有以及特有肠道菌群进行分析,如图 4A所示,5个年龄段组包含6个共有肠道菌门,34岁组包含1个特有菌门,为蓝细菌门(Cyanobacteria)。如图 5A所示,5个年龄组中共有优势菌门为厚壁菌门(Firmicutes)、变形菌门(Proteobacteria)和Epsilonbacteraeota,平均丰度为60.44%、38.69%和5.37%;变形菌门在2~3岁组的相对丰度最高,为66.50%;厚壁菌门在34岁组的相对丰度最高,为70.73%,由此推断出,随着大熊猫年龄增长,厚壁菌门的相对丰度逐渐增加,变形菌门的相对丰度逐渐降低。

A. 在门水平的花瓣Venn图;B.在属水平的花瓣Venn图。图中不同颜色代表不同的分组,重叠区域表示5个年龄段共有的功能数目,非重叠区域表示5个年龄段各自特有的功能数目 A. Venn diagram of petals at the door level; B. Venn diagram of petals at the genus level. Different colors in the figure represent different groups. Overlapping areas represent the number of functions shared by the five age groups, while non overlapping areas represent the number of functions unique to each of the five age groups 图 4 大熊猫肠道菌群在门水平(A)和属水平(B)的Venn图 Fig. 4 Venn diagram of intestinal flora of giant pandas at phyla level (A) and genus level (B)
A. 门水平的物种相对丰度;B. 属水平的物种相对丰度。不同颜色表示不同肠道菌种类 A. Species relative abundance at the phyla level; B. Species relative abundance at the genus level. Different colors indicate different types of intestinal bacteria 图 5 大熊猫肠道菌群在门水平(A)和属水平(B)的物种相对丰度 Fig. 5 Species relative abundance of giant panda intestinal flora at phyla level (A) and genus level (B)

在属水平上,如图 4B所示,5个年龄段组包含42个共有肠道菌属,2~3岁组包含7个特有菌属,占比前两个分别为冷杆菌属(Psychrobacter)和梭状芽胞杆菌属(Shuttleworthia);8~9岁组包含7个特有菌属,占比前两个分别为Ruminiclostridium_5和帕拉苏特菌属(Parasutterella);13~14岁组包含5个特有菌属,占比前两个分别为海洋芽胞杆菌属(Oceanobacillus)和根瘤菌属(Rhizorhapis);18~22岁组包含3个特有菌属,占比前两个分别为黄杆菌属(Luteolibacter)和假黄色单胞菌属(Pseudoxanthomonas);34岁组包含1个特有菌属,为Aliterella_CENA595。如图 5B所示,5个年龄组中共有优势菌属为变形菌门的埃希菌属-志贺菌属(Escherichia/Shigella)和厚壁菌门的狭义梭菌属1(Clostridium_sensu_stricto_1),平均丰度分别为57.35%和15.17%;13~14岁组的优势菌属为魏斯菌属(Weissella),其平均丰度为9.975%,显著高于其他四组(P<0.05),34岁组的优势菌属为八叠球菌属(Sarcina),其平均丰度为53.29%,显著高于其他4组(P<0.05)。

2.5 大熊猫肠道菌群基因功能预测

通过统计COG的L2水平和KEGG通路的L2水平上的功能组成和丰度信息,并进行聚类分析,能预测出不同年龄段大熊猫的肠道菌群基因在功能上的差异。如图 6所示,在COG的L2水平上,13~14岁成年大熊猫肠道菌群基因在细胞运动(cell motility)功能的相对丰度最低,在脂质转运与代谢(lipid transport and metabolism)的相对丰度最高;34岁老年大熊猫在碳水化合物转运和代谢(carbohydrate transport and metabolism)、氨基酸转运和代谢(amino acid transport and metabolism)的相对丰度低于其他年龄段大熊猫,在转录(transcription)、防御机制(defense mechanisms)、信号转导机制(signal transduction mechanisms)的相对丰度高于其他年龄段大熊猫。如图 7所示,在KEGG通路的L2水平上,G2、G3、G4 3个成年大熊猫组别在代谢相关的途径上相对丰度相似,2~3岁亚成年大熊猫在信号传导途径(signal transduction)的相对丰度高于成年和老年大熊猫;34岁老年大熊猫在其他次级代谢物的生物合成途径(biosynthesis of other secondary metabolites)、复制和修复途径(replication and repair)、细胞生长与死亡途径(cell growth and death)等方面的相对丰度高于其他年龄段。由此可推断出,不同年龄段大熊猫肠道菌群中基因功能存在差异。

横轴为不同分组,纵轴为各个COG编号,左侧聚类树为功能聚类树,上方聚类树为样品聚类树,色块表示每个分组对应的COG相对丰度,丰度高低用颜色深浅来表征 The horizontal axis refers to different groups, and the vertical axis refers to each COG number, the left clustering tree refers to the functional clustering tree, the upper clustering tree refers to the sample clustering tree. The color block represents the relative abundance of COG corresponding to each group, and the abundance is represented by the color depth 图 6 大熊猫肠道菌群基因功能在COG的L2水平的聚类热图 Fig. 6 Cluster heat map of gene function of giant panda intestinal flora at L2 level of COG
横轴为不同分组,纵轴为各KO编号,左侧聚类树为功能聚类树,上方聚类树为样品聚类树,色块表示每个分组对应的KO相对丰度,丰度高低用颜色深浅来表征 The horizontal axis refers to different groups, the vertical axis refers to each KO number, the left clustering tree refers to the functional clustering tree, the upper clustering tree refers to the sample clustering tree, and the color block refers to the relative abundance of KO corresponding to each group. The abundance is represented by the color depth 图 7 大熊猫肠道菌群基因功能在KEGG通路的L2水平的聚类热图 Fig. 7 Cluster heat map of gene function of giant panda intestinal flora at L2 level of KEGG pathway
3 讨论

机体的消化代谢、免疫功能等与肠道菌群多样性息息相关,肠道菌群能通过刺激肠道上皮分泌黏液等多种途径加强肠道的屏障保护作用,以维持肠道结构完整性,当肠道菌群多样性降低时,肠道黏膜修复能力下降,更易引起病原体的侵入和肠道性疾病的发生[14-15]。Xue等[16]报道,雌性和雄性大熊猫肠道菌群没有显著性差异,本试验层次聚类结果也表明,相同年龄的雌性和雄性大熊猫肠道菌群结构没有明显差异,所以性别差异对本试验没有影响。

大熊猫年龄划分一直存在争议,胡锦矗曾根据大熊猫咬节值及粪团直径将大熊猫分为幼年组(0~1.5岁)、青年组(1.5~5.5岁)、成年组(5.5~20岁)和老年组(20岁以上),而胡杰等[17]基于离差平方和聚类的方法对大熊猫年龄划分进行了优化,并认为该方法避免了人为因素,值得进一步实践和推广。但上述基于咬节值的年龄划分结果中,成年组跨度均较大,直接应用于大熊猫肠道菌群结构分析难以具有说服力,因此本试验采用层次聚类方法,通过聚类图对11只成年大熊猫进一步分组,最终将16只大熊猫分为2~3岁亚成年组(4只)、8~9岁成年组(4只)、13~14岁成年组(4只)、18~22岁成年组(3只)及34岁老年组(1只),并通过Beta多样性分析得出,除仅有34岁大熊猫所在的G5不参与组间差异比较外,其余各组的组间差异均大于组内差异,即该分组具有可信度。层次聚类结果还显示,34岁老年大熊猫与亚成年、成年大熊猫肠道菌群的结构组成差异较大,罗亚等的研究结果也表明,亚成年、成年大熊猫在肠道菌群结构和优势菌群上与老年大熊猫具有显著差异[14]。这提示大熊猫肠道菌群结构随年龄增长处于动态变化过程中,饲养人员需要在大熊猫不同年龄阶段(尤其是老年阶段)补充不同的营养物质以维持肠道菌群的稳态。

崔明全等[18]的研究显示,成年圈养大熊猫的肠道菌群Alpha多样性最高,其次分别是亚成体和老年体。本试验进一步完善了崔明全等[18]的研究结论,即不同年龄段大熊猫肠道菌群Alpha多样性具有显著差异(P<0.05),其中群落多样性及均匀度呈现随年龄增长先增加后降低的趋势,在18~22岁达到最高。

大熊猫的肠道优势菌群以及与纤维素、半纤维素降解相关酶的丰度与食肉动物接近,而显著不同于食性相似的食草动物[19]。Tun等[20]证实,厚壁菌门和变形菌门是老年和成年大熊猫中最丰富的微生物群,且不同年龄段的大熊猫肠道内优势菌群既有相似性,也存在一定差异[21]。在门水平,本试验中不同年龄段的大熊猫共有优势菌门为厚壁菌门、变形菌门和Epsilonbacteraeota,且随着大熊猫年龄增长,厚壁菌门的相对丰度逐渐增加,变形菌门的相对丰度逐渐降低。目前对Epsilonbacteraeota的研究极少,据Waite等[22]报道,该门与变形菌门联系紧密,近些年才将变形菌门下的一些菌属划分到Epsilonbacteraeota,Yang等[23]的试验结果也证实,给高脂肪小鼠饮食补充香菇多糖时,变形菌门和Epsilonbacteraeota的变化趋势一致,由此推测Epsilonbacteraeota在大熊猫肠道中的功能和变形菌门相似。此外,老年大熊猫包含1个特有菌门,为蓝细菌门。蓝细菌是广泛分布于各种水体、土壤等地的光合细菌,在大熊猫采食竹子或饮水时进入体内,未被消化的蓝细菌随粪便排出体外,王立志和徐谊英[24]也在大熊猫粪便中发现过蓝细菌,推测它仅作为路过菌存在于大熊猫肠道中,与消化吸收无关。

在属水平上,本试验不同年龄段大熊猫共有优势菌属为变形菌门的埃希氏菌属-志贺氏菌属和厚壁菌门的狭义梭菌属1,梭菌属是常见的纤维素降解菌,埃希氏菌属能够维持大熊猫肠道内微生态平衡,这些菌属与大熊猫降解竹子中纤维素关系密切[25]。优势菌属间的相互作用会影响大熊猫对纤维素的消化能力,如Guo等[26]的研究结果显示,圈养大熊猫纤维素酶活性偏低,导致消化竹子的能力弱于野生大熊猫,其原因与肠道菌群结构和优势菌属的差异相关。本试验中13~14岁大熊猫特有的肠道优势菌属为魏斯菌属,34岁老年大熊猫特有的优势菌属为八叠球菌属。魏斯菌是广泛分布于发酵食物、发酵植物和土壤等处的一类乳酸菌,可以抑制人体肠道内腐败细菌的增殖,降低腐败霉菌素的产生[27-28]。魏斯菌为成年大熊猫肠道内的优势菌,说明其肠道处于较为健康的状态[29],这也提示魏斯菌可进一步分离研究以制成大熊猫益生菌菌剂,从而抑制肠道内病原菌的滋生,提高大熊猫免疫水平。八叠球菌为革兰阳性厌氧性球菌,Debey等[30]和Owens等[31]的研究结果表明,八叠球菌为致病菌,黑猩猩感染后出现的“动物流行性神经和胃肠综合征”以及人类感染后出现的胃扩张和肺气肿性胃炎与其相关,但暂无大熊猫感染的相关报道,且本次采样时老年大熊猫表观健康,因此推测八叠球菌属是大熊猫肠道内的条件致病菌,在特定条件下会引起肠道疾病,且在抵抗力较弱的老年大熊猫肠道内相对丰度最高,不过这种推测有待进一步研究确认。

对大熊猫肠道菌群功能进行预测可得,不同年龄段大熊猫肠道菌群的功能存在差异,尤其是34岁老年大熊猫在碳水化合物转运和代谢、氨基酸转运和代谢方面相对丰度最低,在转录、防御机制、信号转导机制的相对丰度高于其他年龄段。Yang等[32]的研究结果表明,圈养野猪的小肠中由于梭状芽胞杆菌等微生物的富集,致使其在小分子营养物质代谢、氨基酸的合成、碳水化合物代谢等功能上丰度较高。猪作为一种杂食动物,与大熊猫在肠道菌群组成上存在相似性,由此可推测大熊猫肠道菌群在碳水化合物转运和代谢、氨基酸的转运和代谢功能上的富集与梭菌属的相对丰度存在一定联系。此外,转录通路和信号转导通路的激活与肠道炎性反应、代谢紊乱等密切相关[33-34],防御机制通路的激活也表明肠道中可能存在应激源,这提示34岁老年大熊猫更易受到病原体的侵害,所以在饲喂老年尤其是高龄大熊猫时,需要格外注意其饮食卫生,可适量补充抗生素。

4 结论

大熊猫肠道菌群多样性、结构及功能随年龄变化呈现一定演替规律。随年龄增长,厚壁菌门的相对丰度逐渐增加,变形菌门的相对丰度逐渐降低,肠道菌群Alpha多样性先增加后降低,在18~22岁达到最高。13~14岁大熊猫特有的优势菌属——魏斯菌属为肠道益生菌,34岁老年大熊猫特有的八叠球菌属为条件致病菌。结果为大熊猫肠道菌群研究数据作了进一步完善,同时可为大熊猫的饲养管理和疾病预防提供指导。

参考文献
[1]
CHEN D Y, LI C W, FENG L, et al. Analysis of the influence of living environment and age on vaginal fungal microbiome in giant pandas (Ailuropoda melanoleuca) by high throughput sequencing[J]. Microb Pathog, 2018, 115: 280-286. DOI:10.1016/j.micpath.2017.12.067
[2]
甄军爱, 杨亚昕, 朱轶强, 等. 基于粪便组学的野生哺乳动物肠道微生态研究进展[J]. 野生动物学报, 2022, 43(3): 836-844.
ZHEN J A, YANG Y X, ZHU Y Q, et al. Advances on intestinal microecology based upon omics analyses of feces samples in wild mammals[J]. Chinese Journal of Wildlife, 2022, 43(3): 836-844. DOI:10.12375/ysdwxb.20220333 (in Chinese)
[3]
程鸣, 邹书珍, 廖梓彤, 等. 野生和圈养大熊猫肠道微生态环境特征及其差异性分析[J]. 西华师范大学学报: 自然科学版, 2020, 41(2): 117-124.
CHENG M, ZOU S Z, LIAO Z T, et al. Characteristics of intestinal micro-ecological environment of wild and captive giant pandas and their differences[J]. Journal of China West Normal University: Natural Sciences, 2020, 41(2): 117-124. (in Chinese)
[4]
晋蕾, 何永果, 杨晓军, 等. 卧龙国家级自然保护区大熊猫主食竹的营养成分与微生物群落结构[J]. 应用与环境生物学报, 2021, 27(5): 1210-1217.
JIN L, HE Y G, YANG X J, et al. Nutrient composition and microbial communities of bamboo at Wolong national nature reserve[J]. Chinese Journal of Applied and Environmental Biology, 2021, 27(5): 1210-1217. DOI:10.19675/j.cnki.1006-687x.2020.08039 (in Chinese)
[5]
WEI F W, WANG X, WU Q. The giant panda gut microbiome[J]. Trends Microbiol, 2015, 23(8): 450-452. DOI:10.1016/j.tim.2015.06.004
[6]
DAVARI S, TALAEI S A, ALAEI H, et al. Probiotics treatment improves diabetes-induced impairment of synaptic activity and cognitive function: behavioral and electrophysiological proofs for microbiome-gut-brain axis[J]. Neuroscience, 2013, 240: 287-296. DOI:10.1016/j.neuroscience.2013.02.055
[7]
FANG W, FANG Z M, ZHOU P, et al. Evidence for lignin oxidation by the giant panda fecal microbiome[J]. PLoS One, 2012, 7(11): e50312. DOI:10.1371/journal.pone.0050312
[8]
何明清. 动物微生态学[M]. 北京: 中国农业出版社, 1994: 146-172.
HE M Q. Animal microecology[M]. Beijing: China Agriculture Press, 1994: 146-172. (in Chinese)
[9]
罗亚, 唐贇, 张丁, 等. 不同年龄大熊猫肠道菌群及其酶活特征分析[J]. 畜牧兽医学报, 2020, 51(4): 763-771.
LUO Y, TANG Y, ZHANG D, et al. Research on microbiota and enzyme activity characteristics in the feces of giant pandas of different ages[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(4): 763-771. (in Chinese)
[10]
WANG A S, ZHAN M Y, PEI E L. Succession of intestinal microbial structure of giant pandas (Ailuropoda melanoleuca) during different developmental stages and its correlation with cellulase activity[J]. Animals (Basel), 2021, 11(8): 2358.
[11]
CHEN S F, ZHOU Y Q, CHEN Y R, et al. fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890. DOI:10.1093/bioinformatics/bty560
[12]
张晓文. 湖南舜皇山植物群落α多样性分析[J]. 中国农学通报, 2019, 35(19): 45-48.
ZHANG X W. Plant communities in Shunhuang Mountain, Hunan: α diversity analysis[J]. Chinese Agricultural Science Bulletin, 2019, 35(19): 45-48. DOI:10.11924/j.issn.1000-6850.casb18100114 (in Chinese)
[13]
WILLIAMS C L, DILL-MCFARLAND K A, VANDEWEGE M W, et al. Dietary shifts may trigger dysbiosis and mucous stools in giant pandas (Ailuropoda melanoleuca)[J]. Front Microbiol, 2016, 7: 661.
[14]
石满杰, 胡瑞宇, 李冰, 等. 肠道菌群失衡在常见呼吸系统疾病中的研究进展[J]. 中国呼吸与危重监护杂志, 2022, 21(3): 215-220.
SHI M J, HU R Y, LI B, et al. Research progress of gut microbiota dysbiosis in common respiratory diseases[J]. Chinese Journal of Respiratory and Critical Care Medicine, 2022, 21(3): 215-220. (in Chinese)
[15]
LEE M, CHANG E B. Inflammatory bowel diseases (IBD) and the microbiome-searching the crime scene for clues[J]. Gastroenterology, 2021, 160(2): 524-537. DOI:10.1053/j.gastro.2020.09.056
[16]
XUE Z S, ZHANG W P, WANG L H, et al. The bamboo-eating giant panda harbors a carnivore-like gut microbiota, with excessive seasonal variations[J]. mBio, 2015, 6(3): e00022-15.
[17]
胡杰, 胡锦矗, 屈植彪, 等. 黄龙大熊猫种群数量及年龄结构调查[J]. 动物学研究, 2000, 21(4): 287-290.
HU J, HU J C, QU Z B, et al. Population and age structure of giant panda in Huanglongsi nature reserve, Sichuan, China[J]. Zoological Research, 2000, 21(4): 287-290. DOI:10.3321/j.issn:0254-5853.2000.04.006 (in Chinese)
[18]
崔明全, 何廷美, 钟志军, 等. 大熊猫粪便菌群ERIC-PCR指纹图谱的分析及优势菌群的鉴定[J]. 畜牧与兽医, 2013, 45(9): 6-11.
CUI M Q, HE T M, ZHONG Z J, et al. Analysis of characteristics of faecal flora pandas by ERIC-PCR fingerprinting and identification of the dominant flora[J]. Animal Husbandry & Veterinary Medicine, 2013, 45(9): 6-11. (in Chinese)
[19]
GUO W, MISHRA S, ZHAO J C, et al. Metagenomic study suggests that the gut microbiota of the giant panda (Ailuropoda melanoleuca) may not be specialized for fiber fermentation[J]. Front Microbiol, 2018, 9: 229. DOI:10.3389/fmicb.2018.00229
[20]
TUN H M, MAUROO N F, YUEN C S, et al. Microbial diversity and evidence of novel homoacetogens in the gut of both geriatric and adult giant pandas (Ailuropoda melanoleuca)[J]. PLoS One, 2014, 9(1): e79902. DOI:10.1371/journal.pone.0079902
[21]
李果, 王鑫, 李才武, 等. 圈养老年大熊猫肠道内菌群结构研究[J]. 黑龙江畜牧兽医, 2019(16): 160-164, 185-186.
LI G, WANG X, LI C W, et al. Study on intestinal flora structure of captive aged giant panda[J]. Heilongjiang Animal Science and Veterinary Medicine, 2019(16): 160-164, 185-186. (in Chinese)
[22]
WAITE D W, VANWONTERGHEM I, RINKE C, et al. Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to epsilonbacteraeota (phyl. nov.)[J]. Front Microbiol, 2017, 8: 682. DOI:10.3389/fmicb.2017.00682
[23]
YANG X Y, ZHENG M X, ZHOU M L, et al. Lentinan supplementation protects the gut-liver axis and prevents steatohepatitis: the role of gut microbiota involved[J]. Front Nutr, 2022, 8: 803691. DOI:10.3389/fnut.2021.803691
[24]
王立志, 徐谊英. 圈养大熊猫粪便中微生物多样性的研究[J]. 四川动物, 2016, 35(1): 17-23.
WANG L Z, XU Y Y. Diversity of microorganism in the feces of captive giant pandas[J]. Sichuan Journal of Zoology, 2016, 35(1): 17-23. (in Chinese)
[25]
晋蕾, 周应敏, 李才武, 等. 野化培训与放归、野生大熊猫肠道菌群的组成和变化[J]. 应用与环境生物学报, 2019, 25(2): 344-350.
JIN L, ZHOU Y M, LI C W, et al. Composition and variation of gut microbiome of trained, preparatory reintroduced, reintroduced and wild giant pandas[J]. Chinese Journal of Applied and Environmental Biology, 2019, 25(2): 344-350. (in Chinese)
[26]
GUO W, MISHRA S, WANG C D, et al. Comparative study of gut microbiota in wild and captive giant pandas (Ailuropoda melanoleuca)[J]. Genes (Basel), 2019, 10(10): 827. DOI:10.3390/genes10100827
[27]
TANG Q S, YI H B, HONG W B, et al. Comparative effects of L. plantarum CGMCC 1258 and L. reuteri LR1 on growth performance, antioxidant function, and intestinal immunity in weaned pigs[J]. Front Vet Sci, 2021, 8: 728849.
[28]
杨旭, 张硕, 王芳, 等. 大熊猫粪便乳酸菌的分离鉴定[J]. 生物资源, 2019, 41(4): 335-341.
YANG X, ZHANG S, WANG F, et al. Isolation and identification of lactic acid bacteria from giant panda feces[J]. Biotic Resources, 2019, 41(4): 335-341. (in Chinese)
[29]
XIONG L, NI X Q, NIU L L, et al. Isolation and preliminary screening of a Weissella confusa strain from giant panda (Ailuropoda melanoleuca)[J]. Probiotics Antimicrob Proteins, 2019, 11(2): 535-544.
[30]
DEBEY B M, BLANCHARD P C, DURFEE P T. Abomasal bloat associated with Sarcina-like bacteria in goat kids[J]. J Am Vet Med Assoc, 1996, 209(8): 1468-1469.
[31]
OWENS L A, COLITTI B, HIRJI I, et al. Publisher correction: a Sarcina bacterium linked to lethal disease in sanctuary chimpanzees in Sierra Leone[J]. Nat Commun, 2021, 12(1): 2035.
[32]
YANG H, HUANG X C, FANG S M, et al. Uncovering the composition of microbial community structure and metagenomics among three gut locations in pigs with distinct fatness[J]. Sci Rep, 2016, 6: 27427.
[33]
程晓范, 朱胡太龙, 刘灵, 等. STAT3信号在炎症性肠病及结肠炎相关结直肠癌发生发展中的作用和相关药物研究进展[J]. 药学学报, 2022, 57(8): 2253-2261.
CHENG X F, ZHU H T L, LIU L, et al. The role of STAT3 in inflammatory bowel disease and colitis- associated cancer and research progress of the related drugs[J]. Acta Pharmaceutica Sinica, 2022, 57(8): 2253-2261. (in Chinese)
[34]
SCHIERING C, KRAUSGRUBER T, CHOMKA A, et al. The alarmin IL-33 promotes regulatory T-cell function in the intestine[J]. Nature, 2014, 513(7519): 564-568.

(编辑   白永平)