瘤胃微生物是指栖息在反刍动物瘤胃中的微生物，对反刍动物的生长，特别是消化代谢具有特殊的意义^{[1, 2]}。瘤胃微生物与反刍动物的共生现象由来已久，与宿主表现出共进化的特点^[3]。因此，瘤胃微生物群落结构是常被用来反映宿主遗传背景的重要指标^{[4, 5, 6]}。尤其是其中的瘤胃细菌，种类最为繁多、数量最为巨大，不仅是瘤胃微生物中最主要的功能群，而且，相比原虫而言，瘤胃细菌种群的种间差异要远远大于种内差异^[6]，更能反映宿主品种的差异性。中国是拥有牛品种最多的国家，但目前在牛方面对瘤胃微生物的研究多以荷斯坦牛等大品种为研究对象^{[7, 8]}，对小品种牛，特别是我国特有的牛品种的瘤胃微生物群落结构及多样性的了解有限。锦江牛属于我国南方小型黄牛，具有耐热性好、耐粗饲、役力较强等特点，是极其宝贵的地方品种资源，主要分布于江西省西北部、锦江中游的低丘岗地和平原地区，2006年调查存栏量超过100万头。目前对锦江牛的研究多集中在资源调查及营养调控上，对锦江牛瘤胃微生物的研究尚未见报道。本研究利用MiSeq高通量测序技术，分析锦江牛瘤胃菌群结构及多样性，旨在为我国南方黄牛种质资源的开发利用奠定基础，也为锦江牛微生态营养的深入研究提供参考。

选择3头安装有永久性瘤胃瘘管的锦江牛，体重(400±20) kg，供瘤胃液采集。试验动物单栏饲喂、自由饮水。饲粮精料为玉米-豆粕型，粗料为稻草，精粗比20 ∶ 80。每天08:00和17:00饲喂。采样日于晨喂前采集瘤胃液，4层灭菌纱布挤压过滤后，分装后迅速放入液氮罐冷冻，-80 ℃保存备用。

取0.5 mL瘤胃液样品进行瘤胃微生物基因组DNA的提取，DNA提取按照天根生化科技(北京)有限公司生产的细菌基因组DNA提取试剂盒说明书进行。采用NanoDrop-ND1000测定提取的DNA浓度，并利用1%的琼脂糖凝胶电泳检测DNA提取质量。提取质量合格的DNA样品，送上海美吉生物公司对细菌16S rDNA V4~V5区(515F~907R)进行测序，测序使用Illumina MiSeq PE250测序平台。

对测得的原始数据进行质量控制，舍弃低质量序列(read尾部碱基质量<20，质控后的read<50 bp)，获得的高质量序列用Usearch 7.1软件进行聚类，以16S rDNA序列97%的相似度作为分类操作单元(operational taxonomic unit，OTU)的划分标准。获得的OTU与RDP数据库(Release 11.1，http://rdp.cme.msu.edu/)比对，通过RDP Classifier鉴定OTU代表性序列的微生物分类地位。丰富度指数(Chao、Ace)和α-多样性指数(Shannon、Simpson)的计算利用Mothur 1.30.1软件完成。

经数据前处理，从3个样本中获得56 763条高质量16S rDNA序列，平均长度393.74 bp，根据序列相似性97%水平划分得到的OTU总数为946个。在OTU水平，3个样品的平均丰富度指数Chao指数、Ace指数分别为836、841，α-多样性指数Shannon指数、Simpson指数分别为4.96、0.021 5。

本研究中，锦江牛瘤胃细菌有15个门被鉴定，拟杆菌门(Bacteroidetes)(64.37%)、厚壁菌门(Firmicutes)(21.20%)和变形菌门(Proteobacteria)(7.59%)是相对含量最多的菌门(表1)。其他菌门，如纤维杆菌门(Fibrobacteres)、螺旋体门(Spirochaetes)、黏胶球形菌门(Lentisphaerae)、柔膜菌门(Tenericutes)、疣微杆菌门(Verrucomicrobia)、浮霉菌门(Planctomycetes)、绿弯菌门(Chloroflexi)、互养菌门(Synergistetes)、TM7、迷踪菌门(Elusimicrobia)、蓝菌门(Cyanobacteria)和装甲菌门(Armatimonadetes)也存在所有个体中，但相对含量较低。在纲水平，鉴定出23个纲，其中拟杆菌纲(Bacteroidia)和梭菌纲(Clostridia)相对含量最多，分别占总菌的52.88%和17.03%。其次是鞘脂杆菌纲(Sphingobacteria)、γ-变形菌纲(Gammaproteobacteria)和Negativicute。

进一步细化分类水平，仅有70.29%和54.26%的序列可以鉴定到科和属。26个科在所有个体中被鉴定，其中普雷沃氏菌科(Prevotellaceae)(36.26%)、瘤胃菌科(Ruminococcaceae)(9.33%)和毛螺菌科(Lachnospiraceae)(4.89%)相对含量最多，其他优势科，包括鞘脂杆菌科(Sphingobacteriaceae)、理研菌科(Rikenellaceae)、琥珀酸弧菌科(Succinivibrionaceae)、紫单胞菌科(Porphyromonadaceae)和氨基酸球菌科(Acidaminococcaceae)等。在属水平，44种属在所有样品中被鉴定。普雷沃氏菌属(Prevotella)相对含量最多，占总菌的28.97%。其次是帕拉普氏菌属(Paraprevotella)、理研菌属(Rikenella)、梭菌_Ⅳ(Clostridium_Ⅳ)、解琥珀酸菌属(Succiniclastium)和瘤胃球菌属(Ruminococcus)。

群落生态学可通过单样品的OTU和α-多样性分析来反映微生物群落的丰度和多样性。本研究通过MiSeq测序获得的锦江牛瘤胃细菌16S rDNA序列共56 763条，划分为946个OTU，低于荷斯坦牛(1 800~2 939)^{[7, 8]}及牦牛(6 642)^[9]的研究结果。Chao指数和Shannon指数也相对荷斯坦牛和牦牛较低(836 vs. 1 311~4 242和4.96 vs. 6.23~8.28)^{[8, 9, 10]}。这说明锦江牛瘤胃细菌的丰富度和多样性相对较低。瘤胃微生物的多样性是宿主和瘤胃微生物之间强烈选择和协同进化的结果，宿主动物的进化历程的差异可能决定其瘤胃的差异^[11]。锦江牛主产区在江西省西北部、锦江中游，属亚热带湿润气候区，狭窄的分布范围以及温和的气候可能导致其瘤胃菌群丰富度和多样性低。另外，饲粮结构也是影响瘤胃微生物的重要因素。Kobayashi等^[12]研究发现野生动物和半家养动物瘤胃菌群丰富度高于家养动物，认为这是由于野生动物和半家养动物饲料的复杂度相对较大。因此，我们估计饲粮结构的相对单一[以天然牧草(春夏秋)和稻草(冬)为主，冬季补饲少量玉米]可能是引起锦江牛瘤胃菌群丰富度和多样性低的又一重要原因。

表 1(Table 1)

表1 锦江牛瘤胃细菌在门、纲、科和属 4个水平的相对含量 Table 1 Relative abundance of bacteria phylum,class,family and genus in the rumen of Jinjiang cattle % 分类 Classification 相对含量 Relative abundance 拟杆菌门 Bacteroidetes 64.37 拟杆菌纲 Bacteroidia 52.88 普雷沃氏菌科 Prevotellaceae 36.26 普雷沃氏菌属 Prevotella 28.97 帕拉普氏菌属 Paraprevotella 3.08 紫单胞菌科 Porphyromonadaceae 1.89 Paludibacter 1.17 理研菌科 Rikenellaceae 2.28 理研菌属 Rikenella 2.28 鞘脂杆菌纲 Sphingobacteria 7.73 鞘脂杆菌科 Sphingobacteriaceae 3.93 厚壁菌门 Firmicutes 21.20 梭菌纲 Clostridia 17.03 瘤胃菌科 Ruminococcaceae 9.33 梭菌_Ⅳ Clostridium_Ⅳ 1.77 醋酸弧菌属 Acetivibrio 1.11 瘤胃球菌属 Ruminococcus 1.53 毛螺菌科 Lachnospiraceae 4.89 丁酸弧菌属 Butyrivibrio 0.81 Negativicutes 2.06 氨基酸球菌科 Acidaminococcaceae 1.62 解琥珀酸菌属 Succiniclasticum 1.62 变形菌门 Proteobacteria 7.59 γ-变形菌纲 Gammaproteobacteria 6.12 琥珀酸弧菌科 Succinivibrionaceae 2.03 琥珀酸弧菌属 Succinivibrio 0.92

表1 锦江牛瘤胃细菌在门、纲、科和属 4个水平的相对含量 Table 1 Relative abundance of bacteria phylum,class,family and genus in the rumen of Jinjiang cattle

在大多数研究中，牛瘤胃中最主要的优势菌群是厚壁菌门，其次是拟杆菌门^{[9, 13, 14, 15]}。本研究中，锦江牛瘤胃内最主要的优势菌是拟杆菌门(64.37%)，厚壁菌门占第2位(21.20%)。也有部分关于荷斯坦奶牛的研究结果^{[7, 16]}与我们一致，这可能是因为饲粮差异的影响。李亚丹^[17]发现持续喂养芒草18个月后，湘西黄牛瘤胃中拟杆菌门比例由16.33%上升到28.15%，厚壁菌门比例由68.88%下降到60.92%。当给安格斯和荷斯坦牛饲喂高纤维的狗牙根时，也发现瘤胃细菌的拟杆菌门增加^[10]。低厚壁菌门与拟杆菌门比例可能更常见于高纤维饲粮。与北方牧草相比，南方牧草普遍存在更加易粗老，粗纤维含量高的特点^[18]，长期高纤维饲粮可能导致锦江牛瘤胃内拟杆菌门数量增加。

瘤胃微生物在种及亚种的水平上具有较高的多样性，在属的阶元不多。Jami等^[7]采集了16头泌乳期的奶牛瘤胃内容物样品，通过焦磷酸测序分析，发现不同个体中有51%的细菌分类相似，核心菌属仅有32个。本研究中，有44个属存在所有个体瘤胃中，其中普雷沃氏菌属是数量最多的属(28.97%)。普雷沃氏菌属遗传多样性高^{[19, 20]}，它在瘤胃中发挥着重要作用，是重要的淀粉降解菌，同时还具有降解蛋白质、小肽、多糖等多种代谢产物的能力^{[21, 22]}。因此，在牦牛^[9]、奶牛^{[10, 15]}及其他反刍动物，如驯鹿^[23]与山羊^{[24, 25]}中，都发现普雷沃氏菌属是瘤胃最重要的菌属。

本研究通过MiSeq高通量测序技术分析了饲喂低精料饲粮的锦江牛瘤胃细菌区系。结果证明，相对其他大品种牛而言，锦江牛瘤胃细菌多样性较低；锦江牛瘤胃中最优势菌门为拟杆菌门，其次是厚壁菌门；普雷沃氏菌属是锦江牛瘤胃中最优势细菌属。