南京农业大学学报2021，44(3):548-553 Journal of Nanjing Agricultural University
nauxb.njau.edu D01:10.7685/jnau.202003019
王智博，姚晓磊，张红琳，等.不同产羔数和发情时期湖羊代谢产物及肠道菌结构的差异分析研究[j ].南京农业大学学报，2021，44 (3): WANG Zhibo ，YA0 Xiaolei，ZUANG Uonglin，et al. Changes of rectal contents microbiota and their metabolites associated with prolificacy and estrus cycle in Uu sheep[J]. Journal of Nanjing Agricultural Lniversitv ，2021，44(3) ：548-553.
(1.南京农业大学江苏省家畜胚胎工程实验室，江苏南京210095;2.南京晓庄学院食品科学学院，江苏南京211171)
摘要：[目的]本试验旨在探究不同繁殖力湖羊直肠内容物代谢产物及肠道菌结构的变化特点，分析
湖羊繁殖性能与肠道 微生物间的相关性。[方法]选取体型均一、生理状态相似的连续三胎产三羔（高繁组，H P 组）和产单羔（低繁组，L P 组）空 怀湖羊各6只，每组3只，对其进行同期发情处理，分别在撤栓时（发情前，L P 1组、H P 1组）和发情时（L P 2组、H P 2组）采 集直肠内容物，测量其挥发性脂肪酸的含量，并采用高通量测序技术对其直肠微生物进行分析。[结果]发情前，L P 组湖羊 乙酸含量显著低于H P 组（P<0.05) ;L P 组湖羊在发情前直肠内容物的丙酸和丁酸含量显著低于发情时(P <0.05)，H P 组湖 羊在发情时的乙酸、丙酸和丁酸含量显著高于发情前(P <0.05)。高通量测序分析结果显示：在发情时，L P 组与H P 组湖羊 菌结构差异显著(P <0.05)舆情控制
多样性的箱型图分析显示L P 组和H P 组湖羊菌结构差异在发情时与发情前均极显著
(P <0.01)。[结论]不同繁殖力和发情时期湖羊的肠道微生物及代谢产物呈现较明显的差异性，为湖羊繁殖力与肠道菌
关系的研究提供了理论基础。
关键词：湖羊；繁殖力；肠道微生物；挥发性脂肪酸（VFA)中图分类号:S826.8
文献标志码：A
文章编号：1000-2030(2021) 03-0548-06
Changes of rectal contents microbiota and their metabolites
associated with prolificacy and estrus cycle in Hu sheep
W A N G Z h ib o 1 ,Y A 0 X ia o le i1 ,Z H A N G H o n g lin 2 ,Z H A N G Y a n li1*
A b s tra c t ： [ Objectives]This study aimed to explore the characteristics of Hu sheep fecal metabolites and intestinal microflora struc- ture，and to analyze the correlation betw-een the reproductive perfor^nance and intestinal microflora. [ Methods ] A total of ten Hu sheeps with identical lambing records(3 records)were selected for estrus treatment at the same tim e，and divided into a high prolificacy group( H P ，litter size = 3) and a low prolificacy group ( L P ，litter size - 1). The feces samples were collected before estrus or during estrus ，respectively. Then the volatile fatty acid( V F A ) content was measured，and the intestinal microorganism was analyzed by high-throughput sequencing. [Results] Before estrus，the content of acetic acid in LP group was significantly lower than that in HP group( P <0.05). The content of propionic acid and butyric acid in LP group was significantly low-er than that of HP group before estrus ( P<0.05) ；the contents of acetic acid，propionic acid and butyric acid during estrus in HP group were significantly higher than those before estrus(P<0.05). The results of high-throughput sequencing showed that there was a sig
nificant difference in the structure of the flora between the LP group and the HP group( P <0.05)，and the difference between the LP group and the HP group was significant(P<0.01)at the peak of estrus and before estrus period. [ Conclusions]There are obvious differences in intestinal microorganisms and metabolites of Hu sheep in different fecundity and estrous periods，which indicate that there may be a certain correlation betw，een reproductive performance and intestinal microbial changes of Hu sheep.Keywords ：Hu sheep ；prolificacy；intestinal microflora；volatile fatty acid(V F A )
动物胃肠道微生物区系平衡直接关系宿主（动物）的健康状况和生产性能，当外界环境、日粮或应激 状态改变时，微生物区系平衡会受到冲击，进而影响机体生理状态[1-3]。有报道，肠道微生物与机体内生
收稿日期：2020-03-09
舞弊三角理论基金项目：江苏现代农业（肉羊）产业技术体系建设专项（JATS 〔 2020〕427)
作者简介：王智博，博士研究生。*通信作者：张艳丽，教授，研究方向为肉羊遗传育种与繁殖,E-mail:**********************.cii 。
548-553.
不同产羔数和发情时期湖羊代谢产物及
肠道菌结构的差异分析研究
王智博1，姚晓磊1，张红琳2，张艳丽1
(1. Jiangsu Livestock Engineering Laboratory，Nanjing Agricultural University，Nanjing 210095，China;
2.School of Food Science，Nanjing Xiaozhuang University，Nanjing 211171，
China)
第3期王智博，等：不同产羔数和发情时期湖羊代谢产物及肠道菌结构的差异分析549
殖激素的浓度有关[4_5]，其中肠道菌结构的变化是调节机体雌激素水平的重要方式[6]。人类医学研究也表明，肠道微生物会参与调节机体雌激素的代谢和卵巢发育。K rita S等[7]报道在日粮中添加益生菌后可以缩短母猪的繁殖周期。唐义梅[8]通过对梅山猪进行粪菌移植，发现肠道微生物可以提早母猪初情日龄，影响发情周期激素水平，促进卵泡发育。这些已有的研究提示：肠道微生物结构与动物繁殖力可能存在一定的相关性。繁殖性能作为动物重要的生产性能指标，也可能与肠道微生物息息相关[
7]。然而，目前关于肠道微生物的研究主要集中在动物营养、免疫和生理等方面，而有关肠道菌与繁殖性能之间的相关研究很少。
湖羊作为我国特有的高繁殖力绵羊品种，产羔率达220%以上。排卵率和产羔数是高繁殖性能绵羊的重要经济性状，而繁殖作为一个复杂的生理过程，其是否与肠道微生物存在相关性尚未见报道。此外，乙酸、丙酸、丁酸和异戊酸等挥发性脂肪酸（v o la tile fa tty a d d,V F A)是消化道内生物代谢的产物，其浓度对动物的肠道健康具有影响[9]。因此，本试验首次尝试以单、多羔湖羊为研究对象，分别在发情前（撤栓时）和发情时2个阶段采集粪样，测定其V F A浓度，并进行高通量测序分析，探究不同繁殖力、不同发情时期湖羊的肠道菌及V F A变化规律，为利用肠道微生物调控母羊繁殖性能提供理论依据和新的思路。
1材料与方法
1.1试验动物及处理
试验在泰州海伦国家级湖羊种羊场进行，依据生产记录选取体重40〜50 k g，年龄2~3岁，连续三胎产三羔（高繁组，h igh Pr o lif k a c y,H P组）和产单羔（低繁组，lo w pr〇l i f k a c y,L P组）的健康空怀母羊12只，进行检疫消毒驱虫防疫。试验羊舍通风良好，自然光照，试验前用消毒剂对羊舍地面、墙壁和羊栏等进行彻底消毒，预饲期间对试验羊进行健康检查、驱虫和健胃。试验羊采用分栏饲养，
试验开始第1天对试验羊进行放栓处理，同时注射氯前列烯醇,11d后撤栓，注射氯前列烯醇，次日用试情公羊开始试情（早晨、中午、傍晚各1次），母羊接受公羊爬跨视为发情。
每只羊采集5粒直肠内容物（在落地污染前用表1试验分组
20 m L固定管直接从羊肛门中取出），分别在撤栓时Table 1 Experim ental group
和发情时采集，采集后样品于-80丈保存。具体试，发情期低繁组.高繁组
Estrus stage Low prolificacy ( LP) High prolificacy ( HP)
验设计见表1。发情前Before estrus LP1 HP1
1.2 试验方法 发情时 During estrus_________LP2________________HP2________ 1.
2.1 V F A的测定参照刘文涛等[l0]的方法，其中气相谱仪测定改为采用毛细吸管柱，柱温120丈，汽化温度180丈，采用氢离子火焰检测器，检测温度180丈，载气为氦气，压力60 k P a，氢气压力50 k P a。1.2.2生物信息学分析将样品送至百迈客生物科技有限公司的Illu m in a M iS e q平台进行高通量测序。利用双末端测序的方法，构建小片段文库进行测序。通过对r e a d s拼接过滤，进行操作分类单元(operational taxonom ic u n it s,0T U)聚类和物种分类分析，并将序列相似度>97%聚类同一' 0T
U。同时使用M o th u r( v丄30)软件，对样品a多样性指数进行评估，得到稀释度曲线，香农指数（S ha nn on)、C H A O、A C E 和S im p o n指数用来评估样品的测序深度、菌多样性。使用Q I I M E软件进行0多样性分析，得到0距离箱线图，并进行U P G M A( unw eighted p a ir-g ro u p m ethod w ith a rith m e tic m e a n)分析和主坐标分析（p rin c ip a l coordinates a n a ly s is,P C o A)。
1.3数据处理
性人君子V F A含量及测序数据用E xce l 2019进行初步整理，采用SPSS 25.0软件进行统计分析，结果数据均以“平均数±标准差（尤±切）”表示，差异显著水平设为0.05。
2结果与分析
2.1单羔和多羔湖羊直肠内容物中V F A含量的差异分析
各组直肠内容物中V F A水平如表2所示。发情前，各组直肠内容物中丙酸、丁酸、异戊酸含量均无显著差异(P>0.05)，但L P组直肠内容物中乙酸含量显著低于H P组（P<0.05)。发情时各组直肠内容物中乙酸、丙酸、丁酸及异戊酸含量均无显著差异(P>0.05)。对于L P组湖羊，发情时直肠内容物中丙酸及丁酸含量显著高于发情前(P<0.05)。发情时,H P组乙酸、丙酸及丁酸含量较发情前显著提高（户<0.05)，但
550南京农业大学学报第44卷
异戊酸含量在各组间差异不显著(P >0.05)。
表2
繁殖力和不同发情时期对湖羊直肠内容物中挥发性脂肪酸的影响
Table 2 Effect of prolificacy and estrus cycle on volatile fa tty acid( V F A ) in H u sheep rectal contents
mm 〇l-L _1
组别
Group 乙酸
Acetic Acid 丙酸
Propionic acid 丁酸
Butyric acid 异戊酸
Isovaleric acid LP1 6.08±0.51a 4.02±0.22a 1.27±0.23a 0.86±0.15HP18.81±1.11b A 5.56±1.08a 1.81±0.61a 0.90±0.24LP214.51±2.8010.92±3.70b 2.91±0.27b 0.84±0.14HP2
工业设备
18.78±4.74b
9.94±2.72b
2.50±0.15b
0.84±0.24
注：不同字母表示差异显著(P<0.05)，小写字母为相同时期不同繁殖力间的比较，大写字母为相同繁殖力不同发情时期的比较。
Note：Different letters indicate significant difference(P <0.05) ^lowercase letters are the comparison of different prolificacy in the sam e period^and
uppercase letters are the caparison of the sam e prolificacy in the different eatrus periods.
2.2单羔和多羔湖羊肠道微生物的a 多样性分析
如图1所示：测序共获得O U T 928条，通过微生物a 多样性分析发现各组等级丰度曲线、稀释性曲线 和香农指数均逐渐趋于平缓，说明样品的测序深度已满足试验要求。通过计算S ha nn on 、C H A O 、Simpson 和A C E 指数，不同发情时期及不同产羔数对肠道菌构成有一定的影响。
800
200
Cl C2 C3 D1 D2 D3 A1 A2 A3 B1 B2 B3 Total
SamplelD
20 000 40 000样本序列数
Number of sequences sampled 60 00020 000 40 000样本序列数
Number of sequences sampled
60 000
图1不同繁殖力及发情期对湖羊肠道微生物菌多样性指数的影响
Fig. 1 Effect of protificacy and estrus stage on intestinal m icroflora diversity index of H u sheep
* P<0.05, ** P<0.01.
2.3发情期湖羊肠道微生物物种注释及差异分析
测序结果表明，在门水平共测得15种菌，属水平上发现172个。湖羊直肠内容物中主要菌构成 为厚壁菌门（F irm icu te s )、拟杆菌门（Bacteroidetes )、疣微菌门（V e rru c o m ic ro b ia )和变形菌门
(P ro te o b a c te ria )。微生物菌门水平分析表明（图2-A ，表3)，H P 组和L P 组湖羊发情时，直肠内容物中
厚壁菌门和蓝藻门菌的比例均低于发情前，而拟杆菌门占比则均有所上升，且厚壁菌门和拟杆菌门与未 发情组存在显著差异。为进一分析样品菌的差异，本研究又进行了属水平的分析（图2-B ，表4)。结果 显示，湖羊发情后，L P 2组和H P 2组的[E ^a c ie r ^m J —coprostanoligenes—g ro u p 和A M e r m a W a 丰度均显著下 降（P <0.05)，R ikene llacea e_R C 9_ _gut_group 和拟杆菌属（B a c ie ro id e s )丰度显著提高（P <0.05)。
2.4不同发情期及产羔数对母羊肠道微生物0多样性的影响
基于W eig hted U n iF ra c 距离的组间0多样性分析表明：与发情前相比，H P 组和L P 组湖羊在发情时微 生物落构成差异均极显著（P <0.01)；与L P 2组相比，H P 2组的菌结构也有显著差异（P <0.05)
(图3-A )。基于W eighted U n iF r a c 的U P G M A 聚类树显示：样品菌基于是否发情大致聚集成2类，LP1
丁酮肟
组和H P 1组间距离较近，而L P 2组和H P 2组则距离较远（图3-B )，证明发情期对肠道菌具有较大的影 响。用P C o A 分析肠道微生物菌结构变化，发现H P 2组和L P 2组呈现区域性聚集，H P 1组和L P 1组所 处区域距离较远（图3-C )，与U P G M A 结果一致。证明繁殖力和发情期对肠道菌均具有较大的影响
。
第3期王智博，等：不同产羔数和发情时期湖羊代谢产物及肠道菌结构的差异分析
551
Genus:I Unclassified Others
unculturedbacteriumfLachnospiraceae Ruminococcaceae_NK4A214_group Prevotellaceae_UCG-003Phocaeicola Alistipes
unculturedbacteriumfRuminococcaceae Ruminococcaceae_UCG-013 Akkermansia
[Eubacterium]_coiprostano\igQnQS_group
Bacteroides
Rikenellaceae_RC9_gut_group Ruminococcaceae_UCG-010 Ruminococcaceae_UCG-005 Christensenellaceae_R-7-group
LP2
H P1
分组Group LP1
LP2 H P1
分组Group
图2不同发情期和繁殖力湖羊肠道微生物区系中门（A )和属（B )水平的相对丰度差异
Fig. 2 Relative abundance o f phylum (A ) and genus (B ) of gut m icroflora of H u sheep w ith different
estrou periods and fecundity
表3不同发情期湖羊肠道微生物区系中门水平的相对丰度差异
Table 3 The relative abundance difference o f phylum in the gut m icroflora of
H u sheep in different estrous period
丰度/% Abundance
尸值尸-value
Phylum
LP1LP2HP1HP2LP HP 发情前
发情时
Firmicutes 62.4556.1965.9758.340.100.010.210.41Bacteroidetes 22.3928.8522.1527.950.160.000.790.79V e rrucomicrobia 7.46 5.57 6.26 5.880.340.810.350.88Proteobacteria 1.53 2.66 1.52 2.420.140.460.990.83Lentisphaerae 1.23 2.63 1.01 1.740.090.300.370.36Cyanobacteria 0.990.590.890.490.320.840.670.77Spirochaetae 1.04 1.370.650.940.670.130.200.60Planctomycetes 0.850.600.290.180.540.300.070.33Actinobacteria 0.500.280.560.550.070.940.490.19Tenericutes
0.21
0.60
0.20
0.21
0.34
0.94
0.51
0.36
表4不同发情期湖羊肠道微生物区系中属水平的相对丰度差异
Table 4 The relative abundance difference o f genus in the gut m icroflora of
H u sheep in different estrous period
丰度/% Abundance
尸值尸-value
G enus
LP1LP2HP1HP2LP HP 发情前
发情时
Christensenellaceae_R-7_group 12.0912.1212.0411.700.990.860.980.72Ruminococcaceae_UCG-0059.798.379.889.580.460.710.910.53Ruminococcaceae_UCG-010 6.147.62 6.427.880.480.430.670.92Rikenellaceae_RC9_gut_group 4.55 5.76 4.83 6.830.060.030.450.19[
]_coprostanoligenes_group
6.02 3.83 6.40 4.560.050.010.520.35Bacteroides 4.47 6.13 4.36
7.610.010.020.820.97Akkermansia
5.85 2.03 3.71 2.060.010.020.150.39Ruminococcaceae_UCG-013
3.21 2.39 3.03 1.670.250.020.800.04uncultured_bacterium_f_Ruminococcaceae
2.73 2.29 2.58 2.140.220.070.620.54Alistipes
1.85
2.35 1.70 2.560.420.010.620.71Prevotellaceae_UCG-003
2.58 1.01 1.91 1.800.080.800.370.13Phocaeicola
1.72
2.47 1.36 2.020.410.170.390.62uncultured_bacterium_f_Lachnospiraceae
1.42
1.63
1.49
2.10
0.75
0.25
0.71
0.55
00806040200
1
1A u o L m p u n q B
9>;BP^
%/鹋
卅友
r
552南京农业大学学报第44卷
^ntrol2 HP1 HP2 Contrail LP1美丽加减法
分组Groups
P C oA -P C l v s PC2
-4
4
8
P C I, 32.08%
12
图3
发情周期对胡羊母羊肠道微生物多样性组成的影响
Fig. 3 Effect o f estrous cycle on m icrobial diversity in the intestine o f H u sheep
A .基于0 多样性箱型图 0 diversity box plot;B. UPGMA 聚类树 UPGMA clustering tree;C.多样本主坐标分析（n = 3) Multisample principal coordinate analysis(PCoA).
3
讨论
动物生殖与肠道微生物关系密切，越来越多的研究发现，肠道中栖息着庞杂繁多的微生物菌，发挥
着动物第二“大脑”的作用，不仅调控着宿主营养吸收、能量供应和疾病的发生与发展，而且影响着动物的 情绪、行为、生理状态等，其肠道菌与动物生殖生理密切相关[11-14]，可能会影响到动物的繁殖性能。饮
食诱导的代谢紊乱绵羊模型中，肠道菌失调与胆汁酸水平降低影响脂溶性维生素A 在肠-睾轴的运输, 睾丸中维生素A 缺乏导致精子发生异常[15]。肠道微生物可以通过代谢生物体内雌激素以调控雌激素功 能活性[4]。湖羊因其多胎性能而倍受关注，在繁殖、营养、生理等方面已有较多研究[16-17]，但有关其繁殖 性能与肠道菌结构的相关性研究未见报道。本研究结果显示，不同繁殖力湖羊的肠道微生物存在差异， 提示肠道微生物对动物的繁殖力有一定影响。研究表明机体雌激素的代谢与肠道菌有关[6，18-19]。据此 推测微生物可能通过影响激素代谢进而影响动物繁殖力，但是具体机制尚需进一步研究。
本试验结果显示，在门水平上，湖羊发情时，厚壁菌门和蓝藻菌门丰度会显著降低；在属水平上，湖羊 发情后，高繁组和低繁组中R u m in o c o c c a c e a e -U G G -013丰度均显著下降，R ikene llacea e—R C 9—g u t_g ro u p 和拟 杆菌属B ac 如以如
丰度显著提高。另夕卜，雌激素也会影响肠道菌中R u m in o co cca ce a e 科的丰度[2Q ] ;K w a
研究发现,
和
R u m ino cocca cea e 参与人体雌激素代谢[5]。而雌激素与动物的生殖系统发育以及
发情周期的调控密切相关[21-22]。本试验基于U P G M A 聚类树和P C o A 结果可以看出，发情时，高繁组和低 繁组湖羊直肠菌结构差异明显，提示湖羊繁殖力的高低与菌结构可能存在一定关联。
V F A 作为肠道微生物的代谢产物，对动物发情周期有影响。田鼠发情会改变其胃肠道中V F A 含 量[23]；动物肠道中V F A 含量与直肠内容物总菌数、肠杆菌数、乳酸杆菌数有关，妊娠动物直肠内容物中
V F A
与对照组差异显著时，肠道菌也差异显著[24]。本研究中，在低繁组和高繁组中，未发情时丙酸及丁
酸含量均显著低于发情组；在高繁组中，未发情时直肠内容物中乙酸含量显著低于发情组。由于乙酸及丙 酸是肠道菌主要的代谢产物，这同样说明不同发情时期湖羊的肠道菌结构呈现不同变化。总之，本试验通过V F A 含量的测定和对菌结构分析，发现不同繁殖力和发情时期的湖羊在肠道菌 结构上均存在差异。表明湖羊繁殖性能与肠道微生物变化可能存在一定的相关性，但具体的调控机制 还有待进一步研究。本试验结果为利用肠道微生物调控母羊繁殖性提供了理论依据和新思路。
参考文献References:
[1]
Andreas K .影响猪的肠道菌的日粮营养策略[J ].中国畜牧兽医,2006,33(4):63-65.
AndreasK. Dietary nutrition strategies affecting pig intestinal flora[J]. China Animal Husbandry & Veterinary Medicine,2006,33(4) ：63-65 (in Chinese).[2]
赵桂英，杨亮宇，段纲，等.断奶仔猪胃肠道正常菌的数量和分区[J ].黑龙江畜牧兽医，2003(9):26-27.
Zhao G Y,Yang L Y, Duan G, et al. Number and division of normal bacterial flora in weaning piglets [ J ]. Heilongjiang Journal of Animal Science and Veterinary Medicine,2003(9) ：26-27(in Chinese).[3]
Hang S Q,Zhu W Y. Gut bacterial and lactobacilli communities of weaning piglets in response to mannan oligosaccharide and sugar beet pulp
c
4
0-4承
c n l n l t o d
A C A C C A D D
B B D B
5 0
5
5 3 2 2 1 1 0 0.
0.0.0.0.0.
90
u a
.s 2p Q o B
.g }u f l P 05.2P Q M
A 诞
s t u l u i s -s A
V
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议